• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于臺陣記錄的局部場地條件地震動效應(yīng)分析

    2016-12-02 05:24:14劉必?zé)?/span>劉培玄王振宇
    地震學(xué)報(bào) 2016年2期
    關(guān)鍵詞:傅里葉山體基巖

    王 偉 劉必?zé)?劉培玄 王振宇 劉 欣

    1) 中國河北三河065201防災(zāi)科技學(xué)院 2) 中國北京100036地殼運(yùn)動監(jiān)測工程研究中心

    ?

    基于臺陣記錄的局部場地條件地震動效應(yīng)分析

    王 偉1),劉必?zé)?)劉培玄2)王振宇1)劉 欣1)

    1) 中國河北三河065201防災(zāi)科技學(xué)院 2) 中國北京100036地殼運(yùn)動監(jiān)測工程研究中心

    基于自貢西山公園山體地形影響臺陣和唐山響堂土層場地臺陣的地震動記錄, 采用互相關(guān)分析、 相干系數(shù)分析、 傅里葉譜比法和考慮相干系數(shù)的傅里葉譜比法, 對比研究了山體地形和土層場地兩種局部場地條件對地震動的影響. 結(jié)果表明, 局部場地條件下的地震動效應(yīng)與地震波的頻率成分相關(guān). 低頻地震波由于其波長較長, 易于穿過土層場地和山體地形發(fā)生衍射, 從而導(dǎo)致二者地震動顯著相關(guān), 但二者在該頻段的地震動場地效應(yīng)差異不大; 較高頻率的地震波由于其波長較短, 更容易受土層場地條件和山體地形條件的影響, 造成二者地震動場地效應(yīng)差異顯著. 此外, 由于山體頂部的尺寸變小, 低頻長波長地震波容易在該位置發(fā)生衍射, 導(dǎo)致山頂周邊測點(diǎn)的地震動相關(guān)性增加, 而較高頻率的地震波則容易引起山頂局部場地發(fā)生共振效應(yīng). 山體地形上相鄰位置的地震動差異一般較大, 相關(guān)性較低, 這與山體地形不同位置的幾何形狀對較高頻率的地震波影響有關(guān); 山底基巖的地震動受山體地形的影響較小, 該位置的地震動效應(yīng)機(jī)制尚需作進(jìn)一步研究.

    地形臺陣 土層場地臺陣 場地效應(yīng) 互相關(guān)分析 相干分析 傅里葉譜比

    引言

    局部場地條件對地震動的影響一直是地震工程領(lǐng)域的研究熱點(diǎn). 場地條件對地表地震動會產(chǎn)生明顯的影響, 是確定設(shè)計(jì)地震作用的關(guān)鍵因素. 對于土層場地, 研究人員一般從土層厚度、 土層結(jié)構(gòu)、 土的動力學(xué)參數(shù)等方面研究場地覆蓋土層對地震動的影響; 對于局部地形, 一般從地形的幾何尺寸、 物質(zhì)組成、 結(jié)構(gòu)構(gòu)造、 入射地震動特性等方面展開數(shù)值模擬分析. 強(qiáng)震動觀測的對象也多以土層場地為主, 專門針對山體地形的觀測臺陣較少, 一般很少將二者對地震動的影響加以對比分析, 而在某一局部范圍內(nèi)對土層場地和山體地形同時(shí)進(jìn)行強(qiáng)震動觀測的研究更是鮮有開展. 大型工程, 例如大跨度橋梁選址時(shí), 橋墩往往位于不同的地貌單元, 在抗震設(shè)計(jì)過程中就需要考慮地震動的多點(diǎn)輸入問題; 如果橋墩分別位于土層場地和山體上, 則兩種局部場地條件的地震動效應(yīng)差異勢必會影響橋梁的響應(yīng). Yu等(2014)對蘆山地震動記錄進(jìn)行分析的結(jié)果表明, 隨著斷層距的增加, 山區(qū)的地震動衰減速度明顯大于平原地區(qū). 在土層場地和相鄰山體布設(shè)強(qiáng)震動觀測儀組成觀測臺陣, 就能夠?qū)@兩種局部場地的地震動效應(yīng)差異進(jìn)行很好的觀測, 因此基于強(qiáng)震動觀測記錄, 研究局部范圍內(nèi)土層場地與山體地形之間的地震動效應(yīng)差異是有意義的.

    傅里葉譜比法首先由Borcherdt(1970)提出, 采用場地及其附近參考場地之間地震動記錄進(jìn)行傅里葉譜幅值比分析. 該方法原理簡單, 處理數(shù)據(jù)方便, 廣泛應(yīng)用于具有參考場地臺陣的地震動場地效應(yīng)研究(Davis, West, 1973; ?elebi, 1987; Gelietal, 1988; Bouchonetal, 1996). Geli等(1988)通過試驗(yàn)和理論分析, 認(rèn)為之前關(guān)于地形對地震動影響的研究結(jié)果只是在定性方面具有一致性, 而數(shù)值模擬結(jié)果總是低估實(shí)際觀測值. Davis和West(1973)利用傅里葉譜比法對山體的地形效應(yīng)地震動記錄進(jìn)行分析, 結(jié)果顯示時(shí)域內(nèi)的峰值加速度比通常遠(yuǎn)小于頻域內(nèi)的傅里葉譜幅值比. Sánchez-Sesma等(1986)采用均方根傅里葉譜比法衡量山體地形對地震動的影響, Paolucci(1999)則對地震動記錄中不同分量交叉偶聯(lián)效應(yīng)對地震動記錄傅里葉譜幅值比的影響進(jìn)行了研究. 基于自貢西山公園地形影響臺陣獲得的2008年汶川主震記錄, 王海云和謝禮立(2010), 楊宇等(2011)以及唐暉等(2012)分別利用傅里葉譜比法和H/V譜比法對自貢西山公園的山體地震動放大系數(shù)進(jìn)行了分析.

    本文擬基于自貢西山公園地形影響臺陣在汶川地震中獲得的主震記錄, 以及唐山響堂土層場地臺陣獲得的地震動記錄, 采用互相關(guān)分析、 相干分析、 傳統(tǒng)傅里葉譜比法和考慮相干系數(shù)的傅里葉譜比法, 對地形臺陣中不同高度的臺站和場地臺陣中不同深度的觀測點(diǎn)之間的地震動記錄分別進(jìn)行分析, 以對比研究山體地形與土層場地兩種局部場地條件對地震動場地效應(yīng)的影響問題.

    1 強(qiáng)震動觀測臺陣簡介

    自貢西山公園地形影響臺陣(以下簡稱為自貢地形臺陣), 位于四川盆地南部的自貢市自流井區(qū)西山公園內(nèi), 為“十五”國家重大項(xiàng)目“中國數(shù)字地震觀測網(wǎng)絡(luò)”中的固定臺陣. 為獲取反映地表地形對地震動影響的多點(diǎn)實(shí)測數(shù)據(jù), 臺陣在水平自由地表場地、 山腳、 山底及山脊等位置布設(shè)了8個(gè)觀測點(diǎn), 其中0號測點(diǎn)為地表土層測點(diǎn), 1號測點(diǎn)為山底基巖測點(diǎn), 其余測點(diǎn)為山體基巖測點(diǎn), 各測點(diǎn)位置如圖1所示(王偉等, 2015).

    唐山響堂土層場地臺陣始建于1992年, 位于河北省唐山市灤縣響堂鄉(xiāng)杜峪村, 是我國用于研究局部場地條件影響的第一個(gè)專用井下臺陣, 由1個(gè)基巖測點(diǎn)、 2個(gè)地表土層測點(diǎn)和3個(gè)井下測點(diǎn)(深度分別為16.0 m, 32.0 m, 47.0 m)組成. 圖2為響堂土層場地臺陣測點(diǎn)分布剖面圖, 其中1號和2號測井的局部場地條件列于表1.3號測井始建于1997年,

    圖1 自貢地形臺陣測點(diǎn)位置示意圖

    圖2 響堂土層場地臺陣測點(diǎn)分布剖面圖

    2000年12月初投入試運(yùn)行, 并于2001年6月安裝普通市話線路, 實(shí)現(xiàn)了北京強(qiáng)震觀測中心與響堂井下臺陣之間的有線遙測功能, 使該臺陣成為我國第一個(gè)具有遙測功能的局部場地條件影響臺陣(盧滔等, 2006).

    表1 響堂土層場地臺陣中1號和2號測井的場地條件

    本文分析使用的地震動記錄, 主要采用自貢地形臺陣獲得的2008年汶川MS8.0地震主震記錄1組, 以及響堂土層場地臺陣1號和2號測井獲得的1994年9月3日MS1.6地震的地震動記錄1組.

    2 臺陣記錄的互相關(guān)分析和相干系數(shù)分析

    互相關(guān)函數(shù)是描述信號x(t)某一時(shí)刻的取值與信號y(t)另一個(gè)時(shí)刻的取值之間的相關(guān)程度, 表示為

    (1)

    式中,rxy(τ)為互相關(guān)函數(shù),T為信號總持時(shí),τ為時(shí)延.

    對于兩個(gè)能量有限的離散信號x(n)和y(n), 經(jīng)過互相關(guān)分析得到的rxy(τ)即為二者之間的互相關(guān)系數(shù). 互相關(guān)分析可以研究系統(tǒng)的時(shí)間滯后性質(zhì)、 系統(tǒng)輸入信號與輸出信號的互相關(guān)函數(shù)以及臺陣中兩個(gè)臺站測點(diǎn)之間地震波傳播的時(shí)間延遲. 對自貢地形臺陣和響堂土層場地臺陣的地震動記錄進(jìn)行基線校正和濾波處理, 基線校正采用Boore(2001)的一次多項(xiàng)式方法, 濾波采用巴特沃斯(Butterworth)帶通濾波器, 頻帶范圍為0.1—25.0 Hz, 校正后的記錄基線漂移現(xiàn)象消除, 滿足使用要求.

    對自貢地形臺陣0—7號測點(diǎn)的地震動加速度依次進(jìn)行互相關(guān)分析, 考慮到EW, NS和UD方向的互相關(guān)分析結(jié)果類似, 僅將EW向的互相關(guān)系數(shù)結(jié)果示于圖3. 可以看出, 自貢地形臺陣的0—5號測點(diǎn), 其地震動加速度的互相關(guān)系數(shù)極低, 表明地震波在0—5號中的相鄰測點(diǎn)之間的地震動特性變化顯著, 差異性很大. 由互相關(guān)函數(shù)的定義可知, 對于3個(gè)變量x(t),y(t)和z(t), 如果x(t)與y(t),y(t)與z(t)的互相關(guān)系數(shù)很大, 則x(t)與z(t)的互相關(guān)系數(shù)也一定很大. 例如, 圖3中5號與6號測點(diǎn)的互相關(guān)系數(shù)r5-6和6號與7號測點(diǎn)的互相關(guān)系數(shù)r6-7的最大值分別為0.84和0.82, 5號與7號測點(diǎn)的互相關(guān)系數(shù)r5-7的最大值則為0.76, 表明5號、 6號與7號測點(diǎn)之間地震動相關(guān)性顯著.

    自貢地形臺陣中地表土層0號測點(diǎn)和山底基巖1號測點(diǎn)分別與其它山體測點(diǎn)地震動的互相關(guān)分析結(jié)果分別如圖4和圖5所示. 可以看出, 雖然0號測點(diǎn)位于地表土層場地, 但其與山體測點(diǎn)之間的互相關(guān)系數(shù)明顯大于1號山底基巖測點(diǎn)與其它山體測點(diǎn)之間的互相關(guān)系數(shù). 這是由于, 與1號山底基巖測點(diǎn)相比, 0號地表土層測點(diǎn)與其余山體基巖測點(diǎn)的位置處于相似的地震波場環(huán)境. 如果山體地形的存在對山腳基巖位置處的入射地震波產(chǎn)生反射、 折射及衍射等影響, 那么1號測點(diǎn)與其余山體基巖測點(diǎn)地震動的相關(guān)性應(yīng)該是顯著的; 而圖5中該相關(guān)性極低, 這表明山體地形的存在對山底基巖位置處的地震動影響很小, 并且進(jìn)一步證明了Blume(1969, 1972)通過對地下人工核爆破條件下與建筑物不同距離的場點(diǎn)地震動反應(yīng)的對比分析, 認(rèn)為建筑物的存在對其周圍場地地震動反應(yīng)影響很小的結(jié)論.

    圖3 自貢地形臺陣1—7號測點(diǎn)EW向加速度互相關(guān)系數(shù)r0-1表示0號與1號測點(diǎn)的互相關(guān)系數(shù), 以此類推, 下同

    相干函數(shù), 又稱為凝聚函數(shù), 其表達(dá)式為

    (2)

    式中,cxy(f)為相干函數(shù),Pxy(f)為互功率譜,Pxx(f)和Pyy(f)分別為兩個(gè)信號的自功率譜. 相干函數(shù)反映了兩個(gè)信號進(jìn)行互功率譜計(jì)算中外來不相干影響的大小, 相干系數(shù)越大表示受外來影響越?。?/p>

    對自貢地形臺陣測點(diǎn)的EW向地震動加速度進(jìn)行相干系數(shù)分析的結(jié)果如圖6所示. 可以看出, 自貢地形臺陣中0—5號測點(diǎn)的相干系數(shù)很小, 表明各測點(diǎn)之間受外來信號影響較大. 由于山體地形的起伏變化, 這些測點(diǎn)的局部地形差異很大, 因此地震動特性變化顯著, 結(jié)合圖3中1—5號測點(diǎn)之間的互相關(guān)系數(shù)較小的特點(diǎn), 進(jìn)一步說明這些測點(diǎn)之間的地震動相對獨(dú)立. 就測點(diǎn)之間的相干系數(shù)所對應(yīng)的頻段范圍來看,c0-1和c1-2在全頻段范圍內(nèi)均小于0.2;c2-3,c3-4和c4-5在5.0 Hz以下頻段總體稍高于5.0 Hz以上頻段; 5—7號測點(diǎn)在5.0 Hz以下頻段的相干系數(shù)突然增大, 整體明顯高于5.0 Hz以上頻段, 表明這3個(gè)測點(diǎn)的地震動在該頻段范圍內(nèi)相干性很強(qiáng). 另外, 考慮到5—7號這3個(gè)測點(diǎn)之間的地震動互相關(guān)系數(shù)很大, 即地震波在山頂測點(diǎn)之間發(fā)生傳播, 因此將5號和7號測點(diǎn)作為地震動輸入, 最高的6號測點(diǎn)作為輸出, 得到相應(yīng)的傅里葉譜幅值比S6/S5和S6/S7, 如圖7所示.

    圖4 自貢地形臺陣0號測點(diǎn)與其它測點(diǎn)的地震動加速度互相關(guān)系數(shù)

    由圖7可知, 山頂測點(diǎn)之間的傅里葉譜幅值比在5.0 Hz以下的放大系數(shù)約為1.0, 即沒有明顯的放大效應(yīng). 導(dǎo)致這種現(xiàn)象的原因, 一方面是在山體基巖中傳播的地震波中5.0 Hz以下的波長較長; 另一方面, 5—7號測點(diǎn)位于山頂, 山頂?shù)某叽缦鄬Φ撞匡@著變小(圖1), 因此該頻段范圍對應(yīng)的較長波長的地震波很容易在山頂附近位置發(fā)生衍射, 從而導(dǎo)致山頂測點(diǎn)的地震動在該頻段范圍內(nèi)的相干系數(shù)很大, 而傅里葉譜幅值沒有明顯放大.在5.0—6.0 Hz的頻段范圍內(nèi), 放大系數(shù)有明顯提高, 最大值可達(dá)1.70, 其原因是地震波在該頻段范圍的頻率成分與山頂?shù)淖哉耦l率相當(dāng), 從而發(fā)生共振效應(yīng); 6.0—25.0 Hz范圍內(nèi)的傅里葉譜幅值比顯著小于1, 說明山頂對該頻段范圍內(nèi)的地震動產(chǎn)生了明顯的壓制作用.

    圖6 自貢地形臺陣測點(diǎn)EW向加速度相干系數(shù)c0-1表示0號與1號測點(diǎn)之間的加速度相干系數(shù), 以此類推, 下同

    圖7 自貢地形臺陣山頂測點(diǎn)傅里葉譜幅值比

    對自貢地形臺陣的地表土層0號測點(diǎn)與山體基巖測點(diǎn)在EW向地震動加速度進(jìn)行相干系數(shù)分析, 結(jié)果如圖8所示. 可以看出, 自貢地形臺陣中0號與1號測點(diǎn)的相干系數(shù)在所有頻段范圍均極??; 0號測點(diǎn)與其余山體基巖測點(diǎn)之間的相干系數(shù)在0—2.0 Hz頻段范圍內(nèi)明顯大于2.0 Hz以上頻段. 這正好可以解釋圖4中0號測點(diǎn)與這些測點(diǎn)相關(guān)性顯著的現(xiàn)象, 即入射地震波中2.0 Hz以下頻段范圍對應(yīng)的波長較長, 其很容易穿過土層場地和山體地形發(fā)生衍射, 導(dǎo)致這兩種局部場地條件下的地震動顯著相關(guān). 在2.0 Hz以上的頻段范圍, 0號測點(diǎn)與其余基巖測點(diǎn)之間的相干系數(shù)顯著降低, 表明土層場地與山體地形這兩種不同的局部場地條件對較高頻率的地震波會產(chǎn)生明顯不同的場地效應(yīng).

    圖8 自貢地形臺陣0號測點(diǎn)與其它測點(diǎn)的EW向加速度相干系數(shù)

    3 山體地形臺陣記錄的傅里葉譜分析

    王偉等(2015)關(guān)于自貢地形臺陣地表土層測點(diǎn)與山體測點(diǎn)EW, NS和UD向3個(gè)分量的加速度傅里葉譜進(jìn)行研究的結(jié)果顯示, 對于相同分量, 各測點(diǎn)在2.0 Hz以下的低頻段傅里葉譜幅值有較好的吻合. 該現(xiàn)象與地表土層測點(diǎn)與山體基巖測點(diǎn)以及山體基巖測點(diǎn)之間地震動在低頻段范圍互相關(guān)系數(shù)、 相干系數(shù)較大的分析是一致的, 即低頻段地震波增加了土層場地和山體地形地震動的相關(guān)性和傅里葉譜幅值的一致性. 隨著測點(diǎn)高程的增加, 在2.0—5.0 Hz范圍內(nèi), 傅里葉譜幅值明顯變大, 且傅里葉譜幅值差異性顯著增大. 根據(jù)地震動互相關(guān)分析和相干系數(shù)分析可知, 這是由于土層場地和山體地形不同位置對地震動中高頻成分的場地效應(yīng)差異性所致. 地表土層測點(diǎn)相對于其它山體基巖測點(diǎn), 在5.0—9.0 Hz的頻段范圍內(nèi)出現(xiàn)明顯峰值, 下面分析出現(xiàn)這種現(xiàn)象的原因.

    自貢地形臺陣的地表土層場地地震動幅值較小, 處于線彈性狀態(tài), 且其底部無豎井臺陣, 考慮到場地條件的相似性, 決定采用響堂土層場地臺陣1994年9月3日MS1.6地震的地震動記錄分析場地不同深度地震動的傅里葉譜幅值比. 為了避免臺陣中不同深度土層下行反射波波場與上行波場發(fā)生相消干涉, 導(dǎo)致參考場地的地震動傅里葉譜在特定頻段的幅值大大降低, 產(chǎn)生譜洞, 并在其中出現(xiàn)假共振現(xiàn)象(王海云, 2014), 本文采用Steidl等(1996)提出的考慮相干系數(shù)的傅里葉譜比法對響堂土層場地臺陣地震動記錄進(jìn)行數(shù)據(jù)處理, 分析土層場地對地震動的放大效應(yīng)機(jī)制. EW, NS和UD向不同深度的加速度傅里葉譜幅值比如圖9所示.

    圖9 響堂土層場地臺陣測點(diǎn)的加速度傅里葉譜幅值比下標(biāo)1, 2和3分別代表了場地豎井臺陣中深度為0, 16 m和32 m的3個(gè)測點(diǎn)

    由圖9分析可知, 對于豎井臺陣3個(gè)分量不同深度的加速度傅里葉譜幅值比, 絕大多數(shù)在0—5.0 Hz頻段內(nèi)較低, 而在5.0 Hz以上的頻段內(nèi)顯著增大; 這與臺陣所處的土層場地條件有關(guān), 即覆蓋層厚度較薄且剪切波速較高的場地, 更有利于放大較高頻率的地震波, 而低頻段地震波的波長較長, 更容易穿過覆蓋土層發(fā)生衍射現(xiàn)象, 這就解釋了自貢地形臺陣地表土層0號測點(diǎn)地震動傅里葉譜在5.0—9.0 Hz頻段內(nèi)出現(xiàn)的峰值現(xiàn)象.

    4 討論與結(jié)論

    基于自貢西山公園地形臺陣和唐山響堂土層場地臺陣的地震動記錄, 進(jìn)行地震動的互相關(guān)分析、 相干系數(shù)分析、 傅里葉譜比法和考慮相干系數(shù)的傅里葉譜比分析. 由于自貢西山公園中地表土層測點(diǎn)與山體基巖測點(diǎn)距離較近, 其震源效應(yīng)和傳播路徑效應(yīng)基本一致, 對入射地震波產(chǎn)生的影響主要是由于這兩種局部場地條件自身差異所致. 對于山體測點(diǎn), 其位置分布在基巖上, 因此對地震動的山體地形效應(yīng)主要受山體局部幾何形狀、 尺寸和地震波入射角度等因素的影響. 如何考慮山體基巖上覆土層與山體基巖地形對地震動的耦合效應(yīng)尚需進(jìn)一步研究. 對山體地形臺陣地震動記錄分析顯示, 其地震動的互相關(guān)性和傅里葉譜形狀的特點(diǎn)均與不同頻段范圍內(nèi)的相干系數(shù)大小有較好的一致性. 通過上述分析, 得出以下結(jié)論和建議:

    1) 土層場地和山體地形的地震動效應(yīng)受地震波頻率成分的影響顯著. 對于自貢西山公園臺陣中的這兩種局部場地條件, 由于低頻地震波的波長較長, 因此容易穿過場地土層和周邊的山體地形發(fā)生衍射現(xiàn)象, 并導(dǎo)致兩種局部場地條件的地震動顯著相關(guān), 該地震動特點(diǎn)在跨越土層場地和山體地形的大跨度橋梁抗震設(shè)計(jì)中必須考慮. 另外, 這種低頻地震波的類型和成因需要進(jìn)一步研究. 這兩種局部場地對地震波中較高頻率成分的地震動影響機(jī)制不同, 導(dǎo)致相應(yīng)頻段范圍內(nèi)地震動相干系數(shù)較低, 即土層場地與山體地形在較高頻段的地震動效應(yīng)差異顯著.

    2) 山體的頂部由于局部尺寸變小, 地震波中的低頻長波長地震波更容易發(fā)生衍射, 導(dǎo)致山頂周邊測點(diǎn)的地震動在低頻段的相關(guān)系數(shù)和相干系數(shù)增加; 較高頻段的地震波則容易在山頂位置引起共振效應(yīng).

    3) 山體地形上相鄰位置的地震動差異性一般較大, 且相關(guān)性低, 這與山體地形不同位置的幾何形狀差異有關(guān); 山底基巖與其余山體基巖測點(diǎn)地震動的互相關(guān)分析和相干系數(shù)分析表明, 山底的地震動受山體地形的影響較小, 該位置的地震動效應(yīng)機(jī)制尚需進(jìn)一步研究.

    盧滔, 周正華, 周雍年, 仲維照. 2006. 關(guān)于Nakamura方法有效性的討論[J]. 地震工程與工程振動, 26(1): 43--48.

    Lu T, Zhou Z H, Zhou Y N, Zhong W Z. 2006. Discussion on validation of Nakamura’s technique[J].EarthquakeEngineeringandEngineeringVibration, 26(1): 43--48 (in Chinese).

    唐暉, 李小軍, 李亞琦. 2012. 自貢西山公園山脊地形場地效應(yīng)分析[J]. 振動與沖擊, 31(8): 74--79.

    Tang H, Li X J, Li Y Q. 2012. Site effect of topography on ground motions of Xishan Park of Zigong city[J].JournalofVibrationandShock, 31(8): 74--79 (in Chinese).

    王海云, 謝禮立. 2010. 自貢市西山公園地形對地震動的影響[J]. 地球物理學(xué)報(bào), 53(7): 1631--1638.

    Wang H Y, Xie L L. 2010. Effects of topography on ground motion in the Xishan Park, Zigong city[J].ChineseJournalofGeophysics, 53(7): 1631--1638 (in Chinese).

    王海云. 2014. 土層場地的放大作用隨深度的變化規(guī)律研究: 以金銀島巖土臺陣為例[J]. 地球物理學(xué)報(bào), 57(5): 1498--1509.

    Wang H Y. 2014. Study on variation of soil site amplification with depth: A case at Treasure Island geotechnical array, San Francisco bay[J].ChineseJournalofGeophysics, 57(5): 1498--1509 (in Chinese).

    王偉, 劉必?zé)簦?劉欣, 楊明亮, 周正華. 2015. 基于汶川MS8.0地震強(qiáng)震動記錄的山體地形效應(yīng)分析[J]. 地震學(xué)報(bào), 37(3): 452--462.

    Wang W, Liu B D, Liu X, Yang M L, Zhou Z H. 2015. Analysis on the hill topography effect based on the strong ground motion records of WenchuanMS8.0 earthquake[J].ActaSeismologicaSinica, 37(3): 452--462 (in Chinese).

    楊宇, 李小軍, 賀秋梅. 2011. 自貢西山公園山脊場地地形和土層效應(yīng)數(shù)值模擬[J]. 震災(zāi)防御技術(shù), 6(4): 436--447.

    Yang Y, Li X J, He Q M. 2011. Numerical simulation for site effect of ridge terrain and overlaying soil in Zigong Xishan Park[J].TechnologyforEarthquakeDisasterPrevention, 6(4): 436--447 (in Chinese).

    Blume J A. 1969. Response of highrise buildings to ground motion from underground nuclear detonations[J].BullSeismolSocAm, 59(6): 2343--2370.

    Blume J A. 1972. Highrise building characteristics and responses determined from nuclear seismology[J].BullSeismolSocAm, 62(2): 519--540.

    Boore D M. 2001. Effect of baseline corrections on displacements and response spectra for several recordings of the 1999 Chi-Chi, Taiwan, earthquake[J].BullSeismolSocAm, 91(5): 1199--1211.

    Borcherdt R D. 1970. Effects of local geology on ground motion near San Francisco bay[J].BullSeismolSocAm, 60(1): 29--61.

    Bouchon M, Schultz C A, Toks?z M N. 1996. Effect of three-dimensional topography on seismic motion[J].JGeophysRes, 101(B3): 5835--5846.

    ?elebi M. 1987. Topographical and geological amplifications determined from strong-motion and aftershock records of the 3 March 1985 Chile earthquake[J].BullSeismolSocAm, 77(4): 1147--1157.

    Davis L L, West L R. 1973. Observed effects of topography on ground motion[J].BullSeismolSocAm, 63(1): 283--298.

    Geli L, Bard P Y, Jullien B. 1988. The effect of topography on earthquake ground motion: A review and new results[J].BullSeismolSocAm, 78(1): 42--63.

    Paolucci R. 1999. Numerical evaluation of the effect of cross-coupling of different components of ground motion in site response analyses[J].BullSeismolSocAm, 89(4): 877--887.

    Sánchez-Sesma F J, Faccioli E, Fregonese R. 1986. An index for measuring the effects of topography on seismic ground motion intensity[J].EarthqEngStructDynam, 14(5): 719--731.

    Steidl J H, Tumarkin A G, Archuleta R J. 1996. What is a reference site?[J].BullSeismolSocAm, 86(6): 1733--1748.

    Yu T, Lu M, Li X J. 2014. Analysis of strong ground motion characterization in 2013 Lushan earthquake[J].ApplMechMater, 580/581/582/583: 1499--1505.

    Analyses on the effect of the local site conditions on the strong motion based on the array records

    Wang Wei1),Liu Bideng1)Liu Peixuan2)Wang Zhenyu1)Liu Xin1)

    1)InstituteofDisasterPrevention,HebeiSanhe065201,China2)NationalEarthquakeInfrastructureService,Beijing100036,China

    Based on the strong motion records from the topography array of Xishan Park in Zigong and the soil site array of Xiangtang in Tangshan, this paper analyzes the site effects of soil site and hill topography on the ground motion by using the cross-correlation function, coherence function, the Fourier spectral ratio methods with and without consideration of coherence function. The results show that the strong motion effect of the local site conditions is relevant to the frequency content of the input seismic wave. For the low frequency component, its wavelength is so long that it is easy to spread across the soil site and the hill topography, and the diffraction makes the strong motion of the two different local site conditions more related; however, the difference of the site effects in the low frequency band for the two local site conditions is weak. For the seismic wave with high frequency and short wavelength, the strong motion is more likely to be affected by the soil site and the hill topography, therefore the difference of the site effects on the strong motion with high frequency is obvious. In addition, with the decrease of the mountaintop dimension, the diffraction is easy to happen when the seismic wave with low frequency and long wavelength spreads to the mountaintop, leading to increase of the cross-correlation of the ground motions at the observing points surrounding the mountaintop. The resonance effect likely happens when the seismic waves with higher frequency spread to the mountaintop. The strong motion of the adjacent points on the hill topography changes greatly, and their correlation is low, which is caused by the effect of the geometrical shape on the strong motion with higher frequency on the different locations of the hill topography. The strong motion of the bedrock at the bottom of hill is weakly affected by the existence of the hill topography, and its mechanism of the strong motion needs to be further studied.

    topography array; soil site array; site effect; cross-correlation analysis; coherence analysis; Fourier spectral ratio

    王偉, 劉必?zé)簦?劉培玄, 王振宇, 劉欣. 2016. 基于臺陣記錄的局部場地條件地震動效應(yīng)分析. 地震學(xué)報(bào), 38(2): 307--317. doi:10.11939/jass.2016.02.014.

    Wang W, Liu B D, Liu P X, Wang Z Y, Liu X. 2016. Analyses on the effect of the local site conditions on the strong motion based on the array records.ActaSeismologicaSinica, 38(2): 307--317. doi:10.11939/jass.2016.02.014.

    中國地震局教師科研基金項(xiàng)目(20150107)、 國家自然科學(xué)基金(51308118)和中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(ZY20120103)聯(lián)合資助.

    2015-11-20收到初稿, 2016-01-28決定采用修改稿.

    e-mail: wwwiem@163.com

    10.11939/jass.2016.02.014

    P315.9

    A

    猜你喜歡
    傅里葉山體基巖
    濟(jì)南市山體修復(fù)中的植物應(yīng)用與技術(shù)——以濟(jì)南市臥虎山山體公園為例
    雙線性傅里葉乘子算子的量化加權(quán)估計(jì)
    基于小波降噪的稀疏傅里葉變換時(shí)延估計(jì)
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    基于改進(jìn)物元的大壩基巖安全評價(jià)
    河北省基巖熱儲開發(fā)利用前景
    基于傅里葉變換的快速TAMVDR算法
    快速離散傅里葉變換算法研究與FPGA實(shí)現(xiàn)
    電測與儀表(2015年5期)2015-04-09 11:30:44
    山體別墅設(shè)計(jì)分析
    人間(2015年21期)2015-03-11 15:23:42
    城市山體保護(hù)規(guī)劃的實(shí)踐與探索——以武漢市新洲區(qū)山體保護(hù)規(guī)劃為例
    成年女人毛片免费观看观看9| 欧美成人一区二区免费高清观看 | 亚洲国产日韩欧美精品在线观看 | 黄色视频,在线免费观看| 亚洲av成人一区二区三| 久久久久九九精品影院| 亚洲精品粉嫩美女一区| 久久久久久免费高清国产稀缺| 国产av一区在线观看免费| 免费看十八禁软件| 亚洲第一av免费看| 亚洲国产精品合色在线| 欧美一级毛片孕妇| 国产精品久久久人人做人人爽| 99国产精品一区二区三区| 精品欧美一区二区三区在线| 一级片免费观看大全| 19禁男女啪啪无遮挡网站| 午夜日韩欧美国产| 三级毛片av免费| 欧美性猛交╳xxx乱大交人| 午夜福利在线在线| 国产精品98久久久久久宅男小说| 级片在线观看| 欧美又色又爽又黄视频| 色综合欧美亚洲国产小说| 成人午夜高清在线视频 | 在线天堂中文资源库| 天天躁夜夜躁狠狠躁躁| 欧美不卡视频在线免费观看 | 国产精品,欧美在线| 麻豆成人av在线观看| 一级毛片女人18水好多| 中国美女看黄片| 88av欧美| 国产精品久久久av美女十八| 久久久久久九九精品二区国产 | av电影中文网址| 日日摸夜夜添夜夜添小说| 天堂影院成人在线观看| 色综合欧美亚洲国产小说| 首页视频小说图片口味搜索| 国产成人欧美在线观看| 麻豆av在线久日| 欧美性猛交╳xxx乱大交人| 日本五十路高清| 国产精品乱码一区二三区的特点| 99精品久久久久人妻精品| 男女之事视频高清在线观看| 国产精品九九99| 神马国产精品三级电影在线观看 | 久久精品国产清高在天天线| 国产一级毛片七仙女欲春2 | 婷婷精品国产亚洲av| 国产av一区在线观看免费| 人人妻人人澡人人看| 丝袜美腿诱惑在线| 波多野结衣高清作品| 91字幕亚洲| 欧美精品亚洲一区二区| 精品久久久久久久久久免费视频| xxxwww97欧美| 国产精品 国内视频| 免费高清在线观看日韩| 亚洲自拍偷在线| 白带黄色成豆腐渣| 天堂影院成人在线观看| 亚洲人成伊人成综合网2020| or卡值多少钱| 亚洲熟妇熟女久久| 久久精品国产亚洲av香蕉五月| 成年女人毛片免费观看观看9| 大香蕉久久成人网| 听说在线观看完整版免费高清| 久久中文看片网| 日韩欧美国产一区二区入口| 国产一级毛片七仙女欲春2 | 天天躁夜夜躁狠狠躁躁| 一本大道久久a久久精品| 中文字幕久久专区| 欧美中文综合在线视频| 非洲黑人性xxxx精品又粗又长| 国产精品香港三级国产av潘金莲| 不卡一级毛片| 国产乱人伦免费视频| 国产伦一二天堂av在线观看| 50天的宝宝边吃奶边哭怎么回事| 熟女少妇亚洲综合色aaa.| 在线看三级毛片| 男女视频在线观看网站免费 | 黄片播放在线免费| 色av中文字幕| av福利片在线| 亚洲aⅴ乱码一区二区在线播放 | 国产激情久久老熟女| 久9热在线精品视频| 免费搜索国产男女视频| 国产蜜桃级精品一区二区三区| 男人舔奶头视频| 一级黄色大片毛片| 国产精品电影一区二区三区| 午夜久久久在线观看| 中文字幕精品亚洲无线码一区 | 熟女少妇亚洲综合色aaa.| 又大又爽又粗| 黑丝袜美女国产一区| 18禁国产床啪视频网站| 在线视频色国产色| 亚洲精品在线观看二区| 999精品在线视频| 一级黄色大片毛片| 亚洲va日本ⅴa欧美va伊人久久| 高潮久久久久久久久久久不卡| 欧美另类亚洲清纯唯美| 色av中文字幕| 最新在线观看一区二区三区| 精品熟女少妇八av免费久了| 青草久久国产| 欧美性猛交黑人性爽| 丰满的人妻完整版| 美女免费视频网站| 亚洲成人久久性| 老汉色∧v一级毛片| 一区二区三区激情视频| 国产精品98久久久久久宅男小说| 国产成+人综合+亚洲专区| 精品国内亚洲2022精品成人| 国内久久婷婷六月综合欲色啪| 精品一区二区三区四区五区乱码| 午夜福利视频1000在线观看| 两个人视频免费观看高清| 久久中文字幕人妻熟女| 深夜精品福利| 99re在线观看精品视频| 女人高潮潮喷娇喘18禁视频| 美女大奶头视频| 18禁黄网站禁片午夜丰满| 亚洲激情在线av| 在线免费观看的www视频| videosex国产| 久久精品国产亚洲av高清一级| 欧美zozozo另类| 香蕉国产在线看| 国产成人欧美在线观看| 欧洲精品卡2卡3卡4卡5卡区| 黄色丝袜av网址大全| 九色国产91popny在线| 在线国产一区二区在线| 欧美成人午夜精品| 国产一级毛片七仙女欲春2 | 国产精品久久久久久人妻精品电影| 中亚洲国语对白在线视频| 在线观看日韩欧美| 可以免费在线观看a视频的电影网站| 最新美女视频免费是黄的| 真人做人爱边吃奶动态| 欧美国产精品va在线观看不卡| 亚洲五月色婷婷综合| 日韩精品青青久久久久久| 丁香欧美五月| 亚洲,欧美精品.| 亚洲av电影在线进入| av视频在线观看入口| 日韩高清综合在线| 青草久久国产| 免费在线观看亚洲国产| 国产单亲对白刺激| 一卡2卡三卡四卡精品乱码亚洲| 成人三级做爰电影| 人人澡人人妻人| 韩国av一区二区三区四区| 欧美不卡视频在线免费观看 | 国产精品国产高清国产av| 黑人巨大精品欧美一区二区mp4| 久久欧美精品欧美久久欧美| 日韩精品免费视频一区二区三区| 欧美黑人精品巨大| 19禁男女啪啪无遮挡网站| 亚洲欧美一区二区三区黑人| 欧美成狂野欧美在线观看| 日韩一卡2卡3卡4卡2021年| 国产成人啪精品午夜网站| 男女那种视频在线观看| 国产熟女xx| 久久精品91蜜桃| 自线自在国产av| 国产精品香港三级国产av潘金莲| 亚洲成av片中文字幕在线观看| 精品免费久久久久久久清纯| 久久久国产成人免费| 好男人在线观看高清免费视频 | 亚洲五月婷婷丁香| 亚洲av电影不卡..在线观看| 久久国产精品男人的天堂亚洲| 精品第一国产精品| 国产视频内射| 两人在一起打扑克的视频| 91大片在线观看| 99国产精品一区二区三区| 国产精品综合久久久久久久免费| 国产一区在线观看成人免费| 亚洲精品粉嫩美女一区| 大型av网站在线播放| 在线观看免费午夜福利视频| 亚洲欧美精品综合久久99| 高潮久久久久久久久久久不卡| 欧美性长视频在线观看| 国产亚洲av高清不卡| 午夜免费成人在线视频| 午夜a级毛片| 啦啦啦 在线观看视频| 69av精品久久久久久| 三级毛片av免费| 日韩欧美国产在线观看| 久久久久免费精品人妻一区二区 | 国产野战对白在线观看| 免费电影在线观看免费观看| 免费看美女性在线毛片视频| 久久久国产成人免费| 热99re8久久精品国产| 淫妇啪啪啪对白视频| 久久精品亚洲精品国产色婷小说| 黑人巨大精品欧美一区二区mp4| 久久欧美精品欧美久久欧美| 国产野战对白在线观看| 日韩三级视频一区二区三区| 99久久无色码亚洲精品果冻| 一级作爱视频免费观看| 亚洲男人的天堂狠狠| 黄色毛片三级朝国网站| 白带黄色成豆腐渣| 免费在线观看影片大全网站| 国内揄拍国产精品人妻在线 | 日本a在线网址| 久久九九热精品免费| 国产亚洲欧美精品永久| 亚洲av中文字字幕乱码综合 | 亚洲av电影在线进入| www.www免费av| 桃色一区二区三区在线观看| 国产激情欧美一区二区| 国产精品亚洲一级av第二区| 波多野结衣av一区二区av| 国产精品九九99| 91国产中文字幕| 国产成人欧美在线观看| 热re99久久国产66热| 亚洲精品在线美女| 91成人精品电影| 亚洲无线在线观看| 中国美女看黄片| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久久久久久久 | 动漫黄色视频在线观看| 无遮挡黄片免费观看| 99在线人妻在线中文字幕| 亚洲国产欧美一区二区综合| 99热6这里只有精品| 中文在线观看免费www的网站 | 日本熟妇午夜| 99国产极品粉嫩在线观看| 国产一区二区激情短视频| 日韩大码丰满熟妇| 男女那种视频在线观看| 在线观看www视频免费| 亚洲男人天堂网一区| 亚洲精品久久国产高清桃花| 91麻豆av在线| 久久 成人 亚洲| 欧美精品亚洲一区二区| 亚洲欧美日韩无卡精品| 免费观看人在逋| 久久草成人影院| 好看av亚洲va欧美ⅴa在| 97人妻精品一区二区三区麻豆 | xxxwww97欧美| 久久青草综合色| 久久香蕉精品热| 在线观看一区二区三区| 欧美zozozo另类| 午夜成年电影在线免费观看| 搞女人的毛片| 97人妻精品一区二区三区麻豆 | 日本撒尿小便嘘嘘汇集6| 国产在线精品亚洲第一网站| 国产精品免费视频内射| 国产精品香港三级国产av潘金莲| 狠狠狠狠99中文字幕| 国产成人av教育| 亚洲精品av麻豆狂野| 欧美不卡视频在线免费观看 | 亚洲黑人精品在线| 免费看十八禁软件| 久久性视频一级片| 波多野结衣高清作品| 黄色视频,在线免费观看| 哪里可以看免费的av片| 美女高潮到喷水免费观看| cao死你这个sao货| 一级作爱视频免费观看| 精品午夜福利视频在线观看一区| 人人澡人人妻人| 老司机在亚洲福利影院| 日韩中文字幕欧美一区二区| 亚洲专区中文字幕在线| 欧美久久黑人一区二区| 久久久久久久久中文| 日本成人三级电影网站| 黄色 视频免费看| 午夜福利高清视频| 国产av在哪里看| 免费在线观看完整版高清| 黄频高清免费视频| 日韩大码丰满熟妇| 十分钟在线观看高清视频www| 黑人欧美特级aaaaaa片| 成人永久免费在线观看视频| 麻豆久久精品国产亚洲av| 久久国产精品男人的天堂亚洲| 亚洲av日韩精品久久久久久密| 中国美女看黄片| 性欧美人与动物交配| 99精品在免费线老司机午夜| 亚洲熟女毛片儿| av欧美777| av在线天堂中文字幕| 亚洲精品久久国产高清桃花| 99久久精品国产亚洲精品| 亚洲成国产人片在线观看| 午夜免费鲁丝| 老司机深夜福利视频在线观看| 老司机靠b影院| 日本 av在线| 国产精品av久久久久免费| 18美女黄网站色大片免费观看| 成年女人毛片免费观看观看9| 中文字幕av电影在线播放| 一二三四在线观看免费中文在| 精品人妻1区二区| 国产精品九九99| 国产色视频综合| 精品久久久久久成人av| 色综合婷婷激情| 国产成人欧美| 男女下面进入的视频免费午夜 | 国产欧美日韩一区二区三| 国产高清有码在线观看视频 | 我的亚洲天堂| 欧美精品亚洲一区二区| 国产真人三级小视频在线观看| 俄罗斯特黄特色一大片| 日韩有码中文字幕| 中文在线观看免费www的网站 | av福利片在线| 黄片大片在线免费观看| 久久久精品国产亚洲av高清涩受| 一本精品99久久精品77| 国产成人啪精品午夜网站| 免费一级毛片在线播放高清视频| 在线十欧美十亚洲十日本专区| 成人av一区二区三区在线看| 欧美av亚洲av综合av国产av| 黑丝袜美女国产一区| 亚洲av第一区精品v没综合| 国产av又大| 人妻久久中文字幕网| 人人澡人人妻人| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利高清视频| 18禁国产床啪视频网站| 亚洲精品美女久久av网站| 91av网站免费观看| 俺也久久电影网| 久久久久国产一级毛片高清牌| 制服人妻中文乱码| 中文资源天堂在线| 久热爱精品视频在线9| 久久精品人妻少妇| 草草在线视频免费看| 琪琪午夜伦伦电影理论片6080| www.www免费av| 最新在线观看一区二区三区| 免费在线观看日本一区| 香蕉久久夜色| 亚洲国产高清在线一区二区三 | 国产精品98久久久久久宅男小说| 黄网站色视频无遮挡免费观看| 2021天堂中文幕一二区在线观 | 久久久国产欧美日韩av| 日韩精品免费视频一区二区三区| 极品教师在线免费播放| 亚洲成av人片免费观看| 淫妇啪啪啪对白视频| 男人舔奶头视频| 淫秽高清视频在线观看| 中文字幕人妻熟女乱码| 成人18禁高潮啪啪吃奶动态图| 国产黄片美女视频| 久久狼人影院| 国产伦人伦偷精品视频| 人人妻人人澡人人看| 看片在线看免费视频| 99久久久亚洲精品蜜臀av| 欧美黑人欧美精品刺激| 色播亚洲综合网| 可以免费在线观看a视频的电影网站| 老司机深夜福利视频在线观看| 亚洲精品一区av在线观看| 91成年电影在线观看| tocl精华| 窝窝影院91人妻| 97碰自拍视频| 99国产极品粉嫩在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产毛片av蜜桃av| 叶爱在线成人免费视频播放| 99精品久久久久人妻精品| 真人做人爱边吃奶动态| 99久久久亚洲精品蜜臀av| 免费高清视频大片| 亚洲一卡2卡3卡4卡5卡精品中文| 我的亚洲天堂| 日韩欧美 国产精品| 国产精品日韩av在线免费观看| 可以在线观看的亚洲视频| 99久久99久久久精品蜜桃| 91成人精品电影| 久久久久久久午夜电影| 在线观看免费午夜福利视频| 久久久久久久久中文| 欧美午夜高清在线| 亚洲成a人片在线一区二区| 成人特级黄色片久久久久久久| 听说在线观看完整版免费高清| 18禁黄网站禁片午夜丰满| 久久 成人 亚洲| 亚洲第一电影网av| 婷婷精品国产亚洲av| 国产蜜桃级精品一区二区三区| 免费看十八禁软件| 精品国产乱子伦一区二区三区| 丝袜在线中文字幕| 久久99热这里只有精品18| 精品国产乱码久久久久久男人| 日本撒尿小便嘘嘘汇集6| 好看av亚洲va欧美ⅴa在| 18禁观看日本| 久久中文字幕人妻熟女| 成人亚洲精品av一区二区| 国产精品免费视频内射| 亚洲av美国av| 久久香蕉激情| 久久精品国产综合久久久| 白带黄色成豆腐渣| 国产精品香港三级国产av潘金莲| 一本精品99久久精品77| 法律面前人人平等表现在哪些方面| 曰老女人黄片| 婷婷六月久久综合丁香| 久久久久国内视频| 午夜福利成人在线免费观看| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 黑丝袜美女国产一区| 亚洲成国产人片在线观看| 岛国在线观看网站| 久久狼人影院| 男女视频在线观看网站免费 | 亚洲国产精品sss在线观看| 成在线人永久免费视频| 91麻豆av在线| 国产视频一区二区在线看| 88av欧美| 美女 人体艺术 gogo| 在线十欧美十亚洲十日本专区| 少妇被粗大的猛进出69影院| 男人的好看免费观看在线视频 | 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| 白带黄色成豆腐渣| 香蕉av资源在线| 香蕉国产在线看| 91在线观看av| 亚洲国产看品久久| 18美女黄网站色大片免费观看| 怎么达到女性高潮| 97碰自拍视频| 亚洲黑人精品在线| 高潮久久久久久久久久久不卡| 人人妻,人人澡人人爽秒播| 久久人妻福利社区极品人妻图片| 国产片内射在线| 好男人电影高清在线观看| 亚洲七黄色美女视频| xxx96com| 天天一区二区日本电影三级| 久久99热这里只有精品18| 久久天躁狠狠躁夜夜2o2o| 中文资源天堂在线| 国内揄拍国产精品人妻在线 | 久久中文字幕人妻熟女| 特大巨黑吊av在线直播 | 亚洲成人国产一区在线观看| 嫁个100分男人电影在线观看| 黄色毛片三级朝国网站| 久久人人精品亚洲av| 欧美黑人巨大hd| 久久精品国产综合久久久| 国产亚洲精品第一综合不卡| 久久亚洲真实| 中文字幕高清在线视频| av欧美777| 最近在线观看免费完整版| 中文字幕精品免费在线观看视频| 国产又色又爽无遮挡免费看| 久久午夜亚洲精品久久| 韩国av一区二区三区四区| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲综合一区二区三区_| 老司机午夜十八禁免费视频| 美女国产高潮福利片在线看| 国产精品一区二区三区四区久久 | 亚洲国产欧美一区二区综合| 最新美女视频免费是黄的| 久久精品人妻少妇| 日韩av在线大香蕉| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区三区在线臀色熟女| 日韩欧美在线二视频| 巨乳人妻的诱惑在线观看| 狂野欧美激情性xxxx| 色在线成人网| 亚洲免费av在线视频| 亚洲自拍偷在线| www国产在线视频色| 丁香欧美五月| 久久草成人影院| 波多野结衣av一区二区av| 成年免费大片在线观看| 999久久久精品免费观看国产| 俺也久久电影网| 精品少妇一区二区三区视频日本电影| 久久天躁狠狠躁夜夜2o2o| 一级a爱视频在线免费观看| 色精品久久人妻99蜜桃| 侵犯人妻中文字幕一二三四区| 又大又爽又粗| 日本成人三级电影网站| 国产精华一区二区三区| 在线观看66精品国产| 这个男人来自地球电影免费观看| 午夜激情av网站| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播放欧美日韩| 18禁裸乳无遮挡免费网站照片 | 老熟妇仑乱视频hdxx| av福利片在线| 久久久久久九九精品二区国产 | 亚洲国产欧美一区二区综合| 亚洲色图av天堂| 日本免费a在线| 亚洲专区字幕在线| 搡老妇女老女人老熟妇| 亚洲va日本ⅴa欧美va伊人久久| 国产成人精品无人区| 午夜激情av网站| 又大又爽又粗| 亚洲精品中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 久久国产精品影院| 午夜免费激情av| 精品国产亚洲在线| 久久九九热精品免费| 少妇粗大呻吟视频| 国产国语露脸激情在线看| 婷婷精品国产亚洲av在线| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 在线观看日韩欧美| 人人妻,人人澡人人爽秒播| 欧美中文日本在线观看视频| 国产一区在线观看成人免费| 午夜久久久在线观看| 免费av毛片视频| 中文在线观看免费www的网站 | 中文亚洲av片在线观看爽| а√天堂www在线а√下载| 脱女人内裤的视频| 啦啦啦 在线观看视频| 午夜日韩欧美国产| 欧美日韩乱码在线| 手机成人av网站| 性色av乱码一区二区三区2| 手机成人av网站| a级毛片在线看网站| 国产精品亚洲一级av第二区| 久久性视频一级片| 一进一出抽搐gif免费好疼| 国产成人一区二区三区免费视频网站| 午夜福利在线在线| 在线av久久热| 在线观看一区二区三区| 国产高清videossex| 亚洲av成人av| 精品国产国语对白av| 国产亚洲欧美在线一区二区| 黑人操中国人逼视频| 性色av乱码一区二区三区2| 国产黄片美女视频| a在线观看视频网站| 一区二区三区高清视频在线| 久久久精品国产亚洲av高清涩受| 村上凉子中文字幕在线| 99热这里只有精品一区 | 精品熟女少妇八av免费久了|