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

    利用背景噪聲互相關研究安寧河—則木河斷裂帶周邊地區(qū)波速變化

    2024-01-01 00:00:00張佩周聰王文青柴旭超劉庚史克旭
    地震研究 2024年4期
    關鍵詞:安寧河背景噪聲

    摘要:使用2013—2018年安寧河—則木河斷裂帶周邊地區(qū)23個寬頻帶流動臺站的連續(xù)波形記錄,通過背景噪聲干涉測量,計算獲取了臺站對之間的相對地震波速變化;計算了5個不同周期的相對波速變化,重點分析了2018年石棉MS4.5地震前后的波速變化。結(jié)果顯示:在1~3 s周期,相對波速變化顯示出明顯的季節(jié)性變化特征,即每年的6—10月,相對波速變化顯示出緩慢的降低趨勢,誤差也較小,而每年的1—5月的相對波速變化呈現(xiàn)出隨機變化特征,且誤差較大。結(jié)合當?shù)貧庀髷?shù)據(jù)推測,地震波速變化可能與季節(jié)性降水引起的孔隙壓力變化有關。3~8 s周期的波速變化呈離散型分布,更長周期的6~15 s和8~20 s波速變化較弱,在零值線附近波動,表明中下地殼介質(zhì)處于相對穩(wěn)定狀態(tài)。高頻段1~4 Hz結(jié)果顯示,石棉地震后近震源區(qū)的部分臺站出現(xiàn)波速降低,而未見震源區(qū)地表破裂,推測波速變化可能是由地震動引起淺部裂隙的張開或閉合所致。

    關鍵詞:背景噪聲;西昌臺陣;安寧河—則木河斷裂帶;波速變化

    中圖分類號:P315.2文獻標識碼:A文章編號:1000-0666(2024)04-0493-10

    doi:10.20015/j.cnki.ISSN1000-0666.2024.0045

    0引言

    在地殼演化過程中,介質(zhì)的物理性質(zhì)將發(fā)生一系列改變,如出現(xiàn)微破裂、擴容、塑性硬化等,因此當?shù)卣鸩ㄍㄟ^地殼傳播時,波速也將發(fā)生相應的變化,所以通過波速的變化研究地下介質(zhì)的物性變化是可行的(馮德益,1981)。早在20世紀六七十年代,研究人員就關注到地震前波速會有微小變化,而預期會存在很大變化的結(jié)果都爭議較大。地質(zhì)環(huán)境中影響波速走時改變的因素很多,包括潮汐應力(De Fazio et al,1973)、地下水的季節(jié)性改變(Sens-Schnfelder,Wegler,2006)、地震或火山活動導致的應力改變(Rubinstein,Beroza,2005;Peng,Ben-Zion,2006)、震前的流體遷移和擴容(Aggarwal et al,1973)等。因此研究地震波速變化對地震監(jiān)測預報和火山預警等都有重要的意義。

    主動源探測技術(Vidale,Li,2003;Wang et al,2008)和重復地震的相似波形(Peng,Ben-Zion,2006;Schaff,Beroza,2004)實現(xiàn)了對地下時變的監(jiān)測,但它們在時空上的不均勻性使得連續(xù)長期的地震波速監(jiān)測依舊困難重重。近年來,得益于地震背景噪聲技術的發(fā)展,使用背景噪聲互相關得到的相似波形推動了連續(xù)監(jiān)測地震波速變化技術的快速發(fā)展。已有很多研究利用這一技術來監(jiān)測地下介質(zhì)的時間變化,包括地下波速的季節(jié)性變化(Sens-Schnfelder,Wegler et al,2006;Wang et al,2017)、與火山噴發(fā)相關的波速變化以及與強震相關的波速變化等(Brenguier et al,2008a,b;Chen et al,2010;Hobiger et al,2016;Liu et al,2014;Nimiya et al,2017;Wu et al,2016;楊建文等,2021)。

    本文收集了2013—2018年安寧河—則木河斷裂帶周邊地區(qū)23個寬頻帶地震儀記錄的連續(xù)波形數(shù)據(jù),利用背景噪聲互相關技術計算了各臺站對之間的波速變化及時空分布特征,討論了波速變化與信號頻帶之間的聯(lián)系、波速變化的可能機理和環(huán)境氣候因素等對介質(zhì)變化的影響。

    1數(shù)據(jù)和方法

    西昌臺陣主要布設在鮮水河斷裂帶南段、安寧河—則木河斷裂帶北段,包含36個寬頻帶地震儀,儀器記錄起始于2013年1月??紤]到數(shù)據(jù)缺失會影響波速變化計算的精度,本文選用了2013—2018年西昌臺陣中23個連續(xù)性較好的寬頻帶地震儀(圖1)記錄的連續(xù)波形數(shù)據(jù)。數(shù)據(jù)預處理參照Liu等(2014)的方式進行:對每天的三分量波形進行10 Hz的重采樣;為了降低地震信號、儀器周邊干擾及其它非穩(wěn)態(tài)噪聲源的影響,得到高質(zhì)量的單臺背景噪聲記錄,對重采樣波形記錄進行了0.02~1.0 Hz帶通濾波、時域歸一化處理及0.05~1 Hz頻帶譜白化處理。

    當臺間距較大時,很難在短時間內(nèi)穩(wěn)定地重構(gòu)出經(jīng)驗格林函數(shù),因此本文僅對臺間距小于200 km的臺站對做進一步分析。我們計算了4個不同周期(1~3 s,3~8 s,6~15 s及8~20 s)的互相關函數(shù)。為了提高信噪比,對于不同周期數(shù)據(jù),將各臺站對當天及其前、后15天(共31天)的互相關函數(shù)疊加來代表該天的經(jīng)驗格林函數(shù),而參考經(jīng)驗格林函數(shù)則由各臺站對2013—2018年整個研究時段的互相關函數(shù)進行疊加得到。

    Froment等(2010)研究發(fā)現(xiàn)噪聲源的季節(jié)性變化會嚴重影響噪聲互相關的面波信號,而其對尾波的影響則很小,所以本文選用噪聲互相關函數(shù)的尾波部分估算地震波速變化。參照Liu等(2018)的測試結(jié)果,在1~3 s周期,選取瑞利波到達后30~90 s的時間窗來計算dt;對于3~8 s周期,dt的時間窗為瑞利波到達后的40~130 s;6~15 s及8~20 s周期計算的dt時間窗則分別為瑞利波到達之后的40~160 s和50~200 s。不同周期的dt計算窗口的選擇如圖2所示,圖中白色粗實線為各周期的平均瑞利波走時曲線,兩條白色虛線之間的部分即為該周期計算dt所選擇的信號時窗。

    假設地殼介質(zhì)相對地震波速在空間上是均勻變化的,相對波速變化可以通過測量每天的經(jīng)驗格林函數(shù)與參考經(jīng)驗格林函數(shù)之間的相對走時偏移來獲得,即dv/v=-dt/t(Sens-Schnfelder,Wegler,2006;Snieder et al,2002)。獲取走時偏移的方法主要有兩種:一是頻率域的移動窗互譜方法(Poupinet et al,1984;Brenguier et al,2008a,b;Clarke et al,2011)。它將用于計算的經(jīng)驗格林函數(shù)尾波部分分成若干個相互重疊的窗口,再將每個窗口信號變換到頻率域,通過互相關譜的計算得到走時偏移,進而獲得波速變化;另一種是時間域的壓縮—拉伸方法(Lobkis,Weaver,2003;Sens-Schnfelder,Wegler,2006)。該方法通過適當拉伸或壓縮每天的經(jīng)驗格林函數(shù),以使經(jīng)驗格林函數(shù)與參考經(jīng)驗格林函數(shù)的相關性達到最大,進而獲得波速相對變化。噪聲的季節(jié)性變化會干擾壓縮-拉伸方法的測量結(jié)果(劉志坤,黃金莉,2010;Zhan et al,2013)。因此,本文采用移動窗互譜方法計算相對走時偏移dt/t。

    2結(jié)果分析

    2.1不同周期的臺站對波速變化

    圖3展示了采用移動窗互譜方法計算出的2013—2018年部分典型臺站對1~3 s周期的相對波速變化。圖3a中各臺站對相對波速變化呈現(xiàn)出明顯的周期性變化:每年6—10月,相對波速變化顯示出緩慢的降低趨勢,計算誤差也較小,這種降低趨勢在圖3b中更加明顯,如紅色矩形框所示。每年的1—5月,相對波速呈現(xiàn)出無規(guī)律的隨機變化特征,并且在隨機變化的時段計算誤差相對也較大,如圖3中暖色所示。部分臺站對相對波速出現(xiàn)較大誤差,如圖3中暗紅色標示,經(jīng)異常核實發(fā)現(xiàn),較大誤差是由于該臺站對在此時間段波形數(shù)據(jù)缺失,用于當天經(jīng)驗格林函數(shù)疊加計算的天數(shù)不足31天,從而導致信噪比較低,計算誤差偏大。

    典型臺站對3~8 s周期的相對波速變化如圖4a所示。各臺站對的相對波速變化總體呈現(xiàn)出無規(guī)律的隨機性:不同位置的臺站對在不同時間段內(nèi),相對波速變化在(-0.4%,+0.4%)范圍內(nèi)呈彌散性分布。3~8 s周期相對波速變化結(jié)果中XC04-XC06部分時間段計算誤差較大,主要是由于XC06臺站缺少波形數(shù)據(jù)。其它臺站對的相對波速變化的計算誤差較低,表明了該周期相對波速變化計算結(jié)果的穩(wěn)定性。圖4b展示了6~15 s周期的相對波速變化,相比于3~8 s周期的結(jié)果,該周期相對波速變化趨于平緩,基本圍繞在零值線附近波動;8~20 s周期的相對波速變化呈現(xiàn)出與6~15 s周期相同的特征。而XC06-XC09臺站對之間出現(xiàn)了誤差較大的相對波速變化異常,經(jīng)檢驗發(fā)現(xiàn)是由XC09臺站部分時段數(shù)據(jù)缺失所致。

    2.2波速變化的空間分布

    對于各臺站對的相對波速變化時間序列,我們僅在1~3 s周期觀測到季節(jié)性波動,其它3個周期的相對波速變化均未呈現(xiàn)出明顯的規(guī)律性分布。為了更好地確認1~3 s周期相對波速變化的季節(jié)性變化特征,將2016—2018年每年1—5月和6—10月的各臺站相對波速度變化分別疊加平均,得到各臺站對之間平均波速變化的空間分布(圖5)。圖5a中干季(1—5月)的臺站對平均波速變化顯示整體波速變化的增大或減小都較大,且有一定的空間分界特征:29°N以北的臺站對的平均波速變化以增大為主,29°N以南的臺站對的平均波速變化以減小為主。相比于干季的平均波速變化,濕季各臺站對的平均波速變化明顯減弱,僅有部分臺站對出現(xiàn)微弱的平均波速度增大或減小(圖5b),說明該地區(qū)的平均波速變化受氣候影響顯著。

    2018年5月16日四川省石棉地區(qū)接連發(fā)生了3次地震(MS3.3、MS4.5、MS4.3)。為了更好地探究地震對震源區(qū)周圍介質(zhì)的影響范圍,本文以石棉MS4.5地震為例,分別計算了地震前2個月(2018年3月15日—2018年5月15日)和地震后2個月(2018年5月17日—2018年7月17日)各臺站對的平均波速變化(圖6)。圖6a顯示,地震前XC06-XC34、XC02-XC05與XC05-XC08臺站對平均波速明顯增大,而XC04-XC06、XC04-XC07、XC07-XC09與XC07-XC34臺站對波速顯著降低。地震后XC04-XC06與XC07-XC08臺站對出現(xiàn)明顯的波速降低現(xiàn)象,XC07-XC34臺站對的波速繼續(xù)降低(圖6b),XC06-XC34臺站對的波速增大現(xiàn)象消失,XC08-XC09臺站對則由震前的波速降低轉(zhuǎn)變?yōu)檎鸷蟮牟ㄋ僭黾印?/p>

    總體來說,波速變化主要局限在震源區(qū)較近范圍內(nèi),且隨著臺站對之間距離的增大,波速變化的平均效應增強,即臺間距較大的臺站對呈現(xiàn)相對弱的波速變化。

    3討論

    3.1降雨對波速變化的影響

    前人研究結(jié)果表明波速變化可能受季節(jié)性變化影響較大(Sens-Schnfelderand,Wegler,2006;Hobiger et al,2016;王俊等,2020)。圖3也顯示1~3 s周期的相對波速變化具有明顯的季節(jié)性變化特征。為了更好地研究波速變化與環(huán)境因素之間的聯(lián)系,將石棉地震震中周邊地區(qū)的相對波速變化與漢源氣象站的日降水量數(shù)據(jù)進行比較,降雨量使用中國氣象局數(shù)據(jù)中心提供的距離石棉地震區(qū)最近的漢源氣象臺數(shù)據(jù)。對比不同位置的臺站對相對波速變化與降雨量(圖7)可以發(fā)現(xiàn):每年的濕季(6—10月)波速變化平緩,干季(1—5月)波速變化劇烈且隨機。尤其是2018年濕季,降雨量明顯增多,波速降低明顯,濕季過后地震波速開始緩慢回升。

    水力作用對地震波速的影響引起了學者的廣泛關注,地下水和降雨會導致地震波速度降低0.01%~0.1%(Sens-Schnfelder,Wegler,2006;Meier et al,2010;Wang et al,2017;Lecocq et al,2017;胡淑芳等,2021)。隨著降雨量的增加,雨水滲透到地下致使地殼最上部數(shù)百米深度范圍內(nèi)孔隙壓力增大,進而導致了剪切模量的降低與地震波速的減小。Wang等(2017)在對日本波速變化的研究中發(fā)現(xiàn)了幅度高達10-4的季節(jié)性變化,并將這種夏季波速的降低歸因于強降雨導致的地殼流體孔隙壓力的增大。謝凡等(2020)在滑坡監(jiān)測中也發(fā)現(xiàn)降雨導致介質(zhì)內(nèi)部的流變特性的改變,特別是其內(nèi)部顆粒體間含水量的增大,介質(zhì)內(nèi)部有效應力的降低,從而導致波速隨彈性模量下降而下降。本文采用Wang等(2017)研究公式計算了石棉地區(qū)由于降雨滲流導致的孔隙壓力變化,結(jié)果表明孔隙壓力的變化相對降雨量有一定的時間延遲,這是由雨水從地表緩慢滲透到地下一定深度范圍內(nèi)所導致的。圖7中各臺站對波速降低的極值與降雨量的最大值也有一定的時間延遲,據(jù)此我們推測在1~3 s周期觀測到的地震波速度的季節(jié)性變化主要受到區(qū)域季節(jié)性降水影響。

    3.2地震對波速變化的影響

    大地震的發(fā)生往往伴隨震源區(qū)的破壞和地表破碎帶的顯現(xiàn)?;诒尘霸肼暤恼鹪磪^(qū)介質(zhì)變化研究為探測地殼深部應力的積累與釋放提供了獨特的視角,發(fā)現(xiàn)了大量同震波速降低和震后快速恢復的證據(jù)(Wegler et al,2009;Brenguier et al,2008a;Chen et al,2010;Cheng et al,2010;趙盼盼等,2012;Froment et al,2013;Liu et al,2014,2018)。研究認為地震導致的波速變化的主要機理為淺部的動態(tài)震動(Sleep,2015)和深部介質(zhì)應力狀態(tài)的改變。

    為了避免直達面波受噪聲源強烈的季節(jié)性變化的影響,本文采用噪聲互相關波形的尾波計算波速變化。噪聲互相關波形尾波的早期部分主要是面波成分,因此可以通過測量不同周期的波速變化來得到不同深度的波速變化(Liu et al,2014;Wu et al,2016)。Froment等(2013)和Liu等(2014)計算了汶川震區(qū)內(nèi)瑞利波對剪切波速度擾動的靈敏度曲線,發(fā)現(xiàn)1~3 s周期面波對3 km以上的淺層地殼更為敏感,而3~8 s和6~15 s面波周期則分別對2~10 km 和5~20 km的地殼結(jié)構(gòu)更為敏感。Liu等(2018)采用10年數(shù)據(jù)對龍門山斷裂帶的波速變化進行了研究,發(fā)現(xiàn)1~3 s和3~8 s面波周期均存在明顯的同震波速降低,而在6~15 s和8~20 s周期未觀測到與汶川地震有關的明顯信號。說明同震波速變化在短周期比較明顯,反映地殼淺部物質(zhì)的狀態(tài)變化。本文對安寧河—則木河斷裂帶周邊地區(qū)的地震波速變化研究發(fā)現(xiàn),波速變化在1~3 s周期具有明顯的季節(jié)性變化特征,在3~8 s周期未發(fā)現(xiàn)與地震相關的信號。

    使用重復地震和噪聲互相關技術的研究表明,波速瞬時變化主要限制在地殼最上部的幾百米范圍內(nèi)(Rubinstein,Beroza,2005;Peng,Ben-Zion,2006;Wegler et al,2009;Takagi et al,2012)。地殼最淺部的這種波速變化主要與強震動下裂隙的張開有關,是一種非線性的場地響應。謝凡等(2020)在2~25 Hz的較高頻帶觀測到了四川長寧MS6.0地震后波速的短時下降(約0.5%)以及恢復過程,認為地球表層介質(zhì)受到地震影響導致顆粒間隙松動,從而使整體彈性模量下降。隨后在重力以及溫度場的共同作用下,介質(zhì)內(nèi)部顆粒間隙重新減小,導致波速變化在同震后出現(xiàn)快速下降,而震后介質(zhì)彈性模量的恢復則促使了地震波速度復原的進程(Niu et al,2008)。為探究地震與震源區(qū)周圍介質(zhì)波速變化之間的聯(lián)系,本文將背景噪聲互相關的頻帶提高到1~4 Hz,進一步研究近地表物質(zhì)的變化情況。結(jié)果顯示,距離2018年石棉MS4.5地震震中較近的幾個臺站對均出現(xiàn)了明顯的同震波速降低,隨后迅速恢復的現(xiàn)象(圖8)。隨著震中距的增大,其它較遠的臺站對并未發(fā)現(xiàn)同震信號。這說明作為一個中等強度地震,石棉地震的影響范圍有限,一般不會出現(xiàn)明顯的地表破裂帶,地震波速的變化更可能是由地震引起淺部裂隙的張開或閉合所致。

    3.3波速變化與測量周期關系

    利用背景噪聲互相關技術計算獲得了2013—2018年安寧河—則木河斷裂帶周邊地區(qū)相對波速的時空變化。結(jié)果表明,僅在1~3 s周期觀測到比較顯著的季節(jié)性變化,其它3個周期(3~8 s、6~15 s及8~20 s)的波速變化無明顯信號特征。已有研究表明,1~3 s周期面波反映地殼之上3 km以內(nèi)的介質(zhì)特性,相應地,3~8 s和6~15 s周期則分別對地下2~10 km和5~20 km的結(jié)構(gòu)更為敏感(Froment et al,2013;Liu et al,2014)。已有研究成果也證實,地震引起的波速變化主要局限于3~4 km的淺層地殼(Peng,Ben-Zion,2006;Wegler et al,2009),甚至只發(fā)生在地表以下數(shù)百米的深度(Rubinstein,Beroza,2005),其機制為強震引起地殼淺部裂隙產(chǎn)生及發(fā)展而導致介質(zhì)的物理破壞。研究區(qū)近年來中小地震頻發(fā),但未發(fā)生造成淺地表的介質(zhì)嚴重破壞的強震。結(jié)合前文降雨量數(shù)據(jù)的分析,區(qū)域季節(jié)性降水是在1~3 s周期觀測到相對波速變化的主要原因。

    對于更長周期的波速變化測量,F(xiàn)roment等(2013)和Obermann等(2014)利用川西流動地震臺陣數(shù)據(jù)檢測到12~20 s周期的波速延遲變化,但Liu等(2018)在對龍門山斷裂帶10年尺度的波速變化研究中并未觀測到較長周期(8~20 s)、與地震有關的、可靠的地震波速變化。安寧河—則木河斷裂帶周邊地區(qū)小震活動較強,但研究時段內(nèi)未發(fā)生大地震,因而不會造成斷層內(nèi)部及其周邊巖石的嚴重破壞,也不會導致斷層內(nèi)部介質(zhì)結(jié)構(gòu)和應力狀態(tài)的強烈變化;另一方面,氣候的季節(jié)性變化主要影響淺地表物質(zhì),降雨導致的地下孔隙壓力變化也主要作用于地殼淺部。因此,我們僅在反映淺部結(jié)構(gòu)特征的1~3 s周期觀測到相對波速變化的季節(jié)性特征,而未能在較長周期的觀測中發(fā)現(xiàn)明顯的相對波速變化規(guī)律。

    4結(jié)論

    本文利用2013—2018年西昌臺陣記錄的連續(xù)波形資料,通過背景噪聲互相關計算了安寧河—則木河斷裂帶及周邊地區(qū)23個臺站之間的相對波速變化,得到如下主要結(jié)論:

    (1)地震波速變化在短周期(1~3 s)呈現(xiàn)出明顯的季節(jié)性變化特征,短周期結(jié)果主要反映地殼最上部的介質(zhì)變化,與降雨主要影響地殼最淺部具有較強的一致性,因此認為在1~3 s周期觀測到的相對波速的季節(jié)性變化主要受到區(qū)域季節(jié)性降水影響。其它周期(3~8 s,6~15 s及8~20 s)沒有觀測到明顯的波速變化特征。

    (2)2018年5月16日石棉地震發(fā)生后,靠近震源區(qū)的部分臺站在高頻段(1~4 Hz)出現(xiàn)震后波速降低,而未見震源區(qū)地表破裂,推測地震波速的變化可能是由震動引起淺部裂隙的張開或閉合所致。但受觀測數(shù)據(jù)和數(shù)據(jù)處理技術等的限制,目前尚未有與石棉MS4.5地震相當震級地震相關的可靠波速變化的報道。

    (3)基于背景噪聲互相關方法計算出的地震波速變化是臺站對之間一定范圍內(nèi)地殼介質(zhì)狀態(tài)的一種平均效應,介質(zhì)變化的精確時間和空間位置及物理機制都尚未有明確的論斷。本文僅觀測到了反映地殼最上部幾千米的短周期波速變化和地震后靠近震源區(qū)部分臺站對的震后波速降低,未見明顯的前震信號或類似火山噴發(fā)前的前兆信號,關于介質(zhì)速度能否對中小地震孕育發(fā)生過程產(chǎn)生明顯的響應還需要進一步研究。中國地震局地球物理研究所“中國地震科學探測臺陣數(shù)據(jù)中心”為本研究提供地震波形數(shù)據(jù),漢源氣象臺降雨數(shù)據(jù)來自于國家氣象科學數(shù)據(jù)中心(http://data.cma.cn/site/index.html),審稿專家對本文提出了寶貴的修改意見,在此一并表示感謝。

    參考文獻:

    馮德益.1981.地震波速異常[M].北京:地震出版社.

    Feng D Y.1981.Velocity anomalies of seismic waves[M].Beijing:Siesmological Press.(in Chinese)

    胡淑芳,李軍,戴宗輝,等.2021.利用背景噪聲研究福建金鐘庫區(qū)地殼介質(zhì)波速變化[J].地震研究,44(1):49-55.

    Hu S F,Li J,Dai Z H,et al.2021.Crustal media seismic velocity changes in Jinzhong reservoir area in Fujian by using ambient noise[J].Journal of Seismological Research,44(1):49-55.(in Chinese)

    劉志坤,黃金莉.2010.利用背景噪聲互相關研究汶川地震震源區(qū)地震波速度變化[J].地球物理學報,53(4):853-863.

    Liu Z K,Huang J L.2010.Temporal changes of seismic velocity around the Wenchuan earthquake fault zone from ambient seismic noise correlation[J].Chinese Journal of Geophysics,53(4):853-863.(in Chinese)

    王俊,鄭定昌,張金川,等.2020.2013年蘆山地震震源區(qū)地殼介質(zhì)地震波速變化的特征分析[J].地球物理學報,63(2):517-531.

    Wang J,Zheng D C,Zhang J C,et al.2020.Seismic velocity changes in the epicentral region of the 2013 Lushan earthquake measured from ambient seismic noise[J].Chinese Journal of Geophysics,63(2):517-531.(in Chinese)

    謝凡,夏開文,黃會寶,等.2020.基于多重散射波波速變化的滑坡實時監(jiān)測方法與應用研究[J].巖石力學與工程學報,39(11):2274-2282.

    Xie F,Xia K W,Huang H B,et al.2020.A new rock mass wave velocity monitoring method of slopes based on multiple scattering waves and its field application[J].Chinese Journal of Rock Mechanics and Engineering,39(11):2274-2282.(in Chinese)

    楊建文,金明培,葉泵,等.2021.利用背景噪聲研究漾濞MS6.4地震前區(qū)域波速變化[J].地震地質(zhì),43(5):1171-1187.

    Yang J W,Jin M P,Ye B,et al.2021.A study on the regional wave velocity changes before the MS6.4 Yangbi earthquake using ambient noise[J].Seismology and Geology,43(5):1171-1187.(in Chinese)

    趙盼盼,陳九輝,Campillo M,等.2012.汶川地震區(qū)地殼速度相對變化的環(huán)境噪聲自相關研究[J].地球物理學報,55(1):137-145.

    Zhao P P,Chen J H,Campillo M,et al.2012.Crustal velocity changes associated with the Wenchuan M8.0 earthquake by auto-correlation function analysis of seismic ambient noise[J].Chinese Journal of Geophysics,55(1):137-145.(in Chinese)

    Aggarwal Y P,Sykes L R,Armbruster J,et al.1973.Premonitory changes in seismic velocities and prediction of earthquakes[J].Nature,241(5385):101-104.

    Brenguier F,Campillo M,Hadziioannou C,et al.2008b.Postseismic relaxation along the San Andreas fault at Parkfield from continuous seismological observations[J].Science,321(5895):1478-1481.

    Brenguier F,Shapiro N M,Campillo M,et al.2008a.Towards forecasting volcanic eruptions using seismic noise[J].Nature Geoscience,1(2):126-130.

    Chen J H,F(xiàn)roment B,Liu Q Y,et al.2010.Distribution of seismic wave speed changes associated with the 12 May 2008 MW7.9 Wenchuan earthquake[J].Geophysical Research Letters,37(18):L18302.

    Cheng X,Niu F L,Wang B.2010.Coseismic velocity change in the rupture zone of the 2008 MW7.9 Wenchuan earthquake observed from ambient seismic noise[J].Bulletin of the Seismological Society of America,100(5B):2539-2550.

    Clarke D,Zaccarelli L,Shapiro N M,et al.2011.Assessment of resolution and accuracy of the Moving Window Cross Spectral technique for monitoring crustal temporal variations using ambient seismic noise[J].Geophysical Journal International,186(2):867-882.

    De Fazio T L,Aki K,Alba J.1973.Solid earth tide and observed change in the in situ seismic velocity[J].Journal of Geophysical Research,78(8):1319-1322.

    Froment B,Campillo M,Chen J H,et al.2013.Deformation at depth associated with the 12 May 2008 MW7.9 Wenchuan earthquake from seismic ambient noise monitoring[J].Geophysical Research Letters,40(1):78-82.

    Froment B,Campillo M,Roux P,et al.2010.Estimation of the effect of nonisotropically distributed energy on the apparent arrival time in correlations[J].Geophysics,75(5):SA85-SA93.

    Hobiger M,Wegler U,Shiomi K,et al.2016.Coseismic and post-seismic velocity changes detected by passive image interferometry:comparison of one great and five strong earthquakes in Japan[J].Geophysical Journal International,205(2):1053-1073.

    Lecocq T,Longuevergne L,Pedersen H A,et al.2017.Monitoring ground water storage at mesoscale using seismic noise:30 years of continuous observation and thermoelastic and hydrological modeling[J].Scientific Reports,7(1):1-16.

    Liu Z K,Huang J L,He P,et al.2018.Ambient noise monitoring of seismic velocity around the Longmenshan fault zone from 10 years of continuous observation[J].Journal of Geophysical Research:Solid Earth,123(10):8979-8994.

    Liu Z K,Huang J L,Peng Z G,et al.2014.Seismic velocity changes in the epicentral region of the 2008 Wenchuan earthquake measured from threecomponent ambient noise correlation techniques[J].Geophysical Research Letters,41(1):37-42.

    Lobkis O I,Weaver R L.2003.Coda-wave interferometry in finite solids:Recovery of P-to-S conversion rates in an elasto-dynamic billiard[J].Physical Review Letters,90(25):254302.

    Meier U,Shapiro N M,Brenguier F.2010.Detecting seasonal variations in seismic velocities within Los Angeles basin from correlations of ambient seismic noise[J].Geophysical Journal International,181(2):985-996.

    Nimiya H,Ikeda T,Tsuji T.2017.Spatial and temporal seismic velocity changes on Kyushu Island during the 2016 Kumamoto earthquake[J].Science Advances,3(11):e1700813.

    Niu F L,Silver P G,Daley T M,et al.2008.Preseismic velocity changes observed from active source monitoring at the Parkfield SAFOD drill site[J].Nature,454(7201):204-208.

    Obermann A,F(xiàn)roment B,Campillo M,et al.2014.Seismic noise correlations to image structural and mechanical changes associated with the MW7.9 2008 Wenchuan earthquake[J].Journal of Geophysical Research:Solid Earth,119(4):3155-3168.

    Peng Z G,Ben-Zion Y.2006.Temporal changes of shallow seismic velocity around the Karadere-Duzce branch of the north Anatolian fault and strong ground motion[J].Pure and Applied Geophysics,163(2):567-600.

    Rubinstein J L,Beroza G C.2005.Depth constraints on nonlinear strong ground motion from the 2004 Parkfield earthquake[J].Geophysical Research Letters,32(14):1-5.

    Schaff D P,Beroza G C.2004.Coseismic and postseismic velocity changes measured by repeating earthquakes[J].Journal of Geophysical Research:Solid Earth,109(B10):B10302.

    Sens-Schnfelder C,Wegler U.2006.Passive Image Interferemetry and seasonal variations of seismic velocities at Merapi Volcano,Indonesia[J].Geophysical Research Letters,33(21):229-236.

    Sleep N H.2015.Shallow S-wave well logs as an indicator of past strong shaking from earthquakes on the Newport-Inglewood fault[J].Bulletin of the Seismological Society of America,105(5):2696-2703.

    Snieder R,Grêt A,Douma H,et al.2002.Coda wave interferometry for estimating nonlinear behavior in seismic velocity[J].Science,295(5563):2253-2255.

    Takagi R,Okada T,Nakahara H,et al.2012.Coseismic velocity change in and around the focal region of the 2008 Iwate-Miyagi Nairiku earthquake[J].Journal of Geophysical Research:Solid Earth,117(B6):B06315.

    Vidale J E,Li Y G.2003.Damage to the shallow Landers fault from the nearby Hector Mine earthquake[J].Nature,421(6922):524-526.

    Wang B S,Zhu P,Chen Y,et al.2008.Continuous subsurface velocity measurement with coda wave interferometry[J].Journal of Geophysical Research:Solid Earth,113(B12):B12313.

    Wang Q Y,Brenguier F,Campillo M,et al.2017.Seasonal crustal seismic velocity changes throughout Japan[J].Journal of Geophysical Research:Solid Earth,122(10):7897-8002.

    Wegler U,Nakahara H,Sens-Schnfelder C,et al.2009.Sudden drop of seismic velocity after the 2004 MW6.6 mid-Niigata earthquake,Japan,observed with Passive Image Interferometry[J].Journal of Geophysical Research:Solid Earth,114(B6):B06305.

    Wu C,Delorey A,Brenguier F,et al.2016.Constraining depth range of S wave velocity decrease after large earthquakes near Parkfield,California[J].Geophysical Research Letters,43(12):6129-6136.

    Zhan Z W,Tsai V C,Clayton R W.2013.Spurious velocity changes caused by temporal variations in ambient noise frequency content[J].Geophysical Journal International,194(3):1574-1581.

    Study on the Seismic Wave Velocity Variation around the AnningheZemuhe

    Fault Zone by Ambient Noise Cross-correlation

    ZHANG Pei1,2,ZHOU Cong1,Wang Wenqing1,Chai Xuchao1,LIU Geng1,SHI Kexu1

    (1.The Second Monitoring and Application Center,China Earthquake Administration,Xian 710054,Shaanxi,China)

    (2.Institute of Geophysics,China Earthquake Administration,Beijing 100081,China)

    Abstract

    We presented the 6-year continuous changes of the seismic wave velocity from 2013 to 2018 around the Anninghe-Zemuhe fault zone.The continuous waveforms recorded by the Xichang array consisting of 23 broadband stations,and ambient noise correlation techniques were used to calculate the wave velocity changes in five different period bands.A significant characteristic of the seasonal variation of the wave velocity changes was observed in the period band of 1-3 s:from June to October each year,the relative wave velocity changes(dv/v)showed a trend of slow decreasing,and the error of the speed change was small.The relative speed change from January to May showed the characteristics of random change with larger errors.Referring to the local meteorological data,we speculated that the wave velocity changes may be related to the pore pressure changes caused by seasonal precipitation.The velocity changes in the 3-8 s period band were discretely distributed,while velocity changes in the period bands 6-15 s and 8-20 s were weaker and fluctuated near the zero line,this indicated that the middle and lower crustal media were in a relatively stable state.In the focal area of the Shimian MS4.5 earthquake,we had also detected wave velocity reduction at some pairedstation in the 1-4 Hz frequency band.Although there were frequent small earthquakes in Shimian region,no surface rupture had been spotted yet.It is more likely that the opening or closing of shallow cracks induced by ground motions contributes to the changes in wave velocity.

    Keywords:ambient noise;the Xichang array;the AnningheZemuhe fault zone;wavevelocity variation

    猜你喜歡
    安寧河背景噪聲
    窄帶電力線通信信道背景噪聲抑制方法
    環(huán)境背景噪聲對飛機噪聲監(jiān)測結(jié)果的影響
    利用背景噪聲研究福建金鐘庫區(qū)地殼介質(zhì)波速變化
    地震研究(2021年1期)2021-04-13 01:05:04
    突然安寧
    揚子江(2020年4期)2020-08-04 11:59:33
    安寧河東支斷裂主要控制點分析
    利用GPS資料研究安寧河
    地震研究(2018年3期)2018-11-26 09:32:10
    八月三號, 晴欲雨(外一首)
    中國詩歌(2018年5期)2018-11-14 20:52:30
    安寧河涼山段重金屬時空分布及主要因子分析
    應用背景噪聲成像研究祁連山地區(qū)地殼S波速度結(jié)構(gòu)
    地震研究(2017年3期)2017-11-06 23:38:05
    海上單道地震勘探中船舶等背景噪聲的影響分析及壓制
    男人舔女人下体高潮全视频| 波多野结衣高清作品| 少妇的逼水好多| 亚洲第一电影网av| 1024手机看黄色片| 搡老岳熟女国产| 午夜福利18| 又粗又爽又猛毛片免费看| 日韩中字成人| 成人国产一区最新在线观看| 午夜福利在线在线| 亚洲片人在线观看| 亚洲欧美日韩无卡精品| 久久亚洲精品不卡| 亚洲精品一区av在线观看| 精品久久久久久久久亚洲 | 亚洲色图av天堂| 免费看光身美女| 欧美精品啪啪一区二区三区| 国产极品精品免费视频能看的| 亚洲 国产 在线| 美女免费视频网站| 欧美bdsm另类| 久久久久久久久大av| 国产高清激情床上av| 99热只有精品国产| 他把我摸到了高潮在线观看| 亚洲成人精品中文字幕电影| 国产一区二区激情短视频| 色吧在线观看| 五月玫瑰六月丁香| 国产真实伦视频高清在线观看 | 亚洲成人中文字幕在线播放| 国产精品影院久久| www.色视频.com| 亚洲精品粉嫩美女一区| 我要搜黄色片| 狠狠狠狠99中文字幕| a在线观看视频网站| 麻豆国产97在线/欧美| 综合色av麻豆| 精品久久久久久久久av| 欧美一区二区亚洲| 色播亚洲综合网| 香蕉av资源在线| 在线观看一区二区三区| 国产精品伦人一区二区| 最新在线观看一区二区三区| 男人的好看免费观看在线视频| 欧美国产日韩亚洲一区| 97碰自拍视频| a级毛片a级免费在线| 久久国产精品人妻蜜桃| 欧美国产日韩亚洲一区| 国产成人影院久久av| 亚洲国产精品成人综合色| 岛国在线免费视频观看| 欧美性猛交黑人性爽| 少妇的逼好多水| av专区在线播放| 国产精品,欧美在线| 日韩欧美三级三区| 亚洲欧美日韩无卡精品| 国产老妇女一区| 亚洲一区二区三区不卡视频| 老司机深夜福利视频在线观看| 999久久久精品免费观看国产| 99国产极品粉嫩在线观看| 搡女人真爽免费视频火全软件 | 俺也久久电影网| 狠狠狠狠99中文字幕| 亚洲最大成人手机在线| 又爽又黄无遮挡网站| 高清在线国产一区| 午夜视频国产福利| 黄色日韩在线| 精品人妻偷拍中文字幕| 亚洲av免费在线观看| 中文字幕熟女人妻在线| 最近在线观看免费完整版| 夜夜夜夜夜久久久久| 十八禁国产超污无遮挡网站| 午夜视频国产福利| 不卡一级毛片| 国产三级在线视频| 国产精品一区二区三区四区久久| 少妇裸体淫交视频免费看高清| 一级毛片久久久久久久久女| 成年人黄色毛片网站| 亚洲国产精品久久男人天堂| 99久国产av精品| 窝窝影院91人妻| 免费看光身美女| 99久久精品一区二区三区| 91九色精品人成在线观看| 亚洲在线自拍视频| 国产色婷婷99| 国内精品久久久久精免费| 欧美绝顶高潮抽搐喷水| 国产探花极品一区二区| 亚洲色图av天堂| 激情在线观看视频在线高清| 免费看光身美女| 一区二区三区高清视频在线| 久久热精品热| 久久人人精品亚洲av| 无遮挡黄片免费观看| 中文字幕人成人乱码亚洲影| 99在线视频只有这里精品首页| 日本一二三区视频观看| 老司机午夜十八禁免费视频| 99热这里只有精品一区| 人妻丰满熟妇av一区二区三区| 内射极品少妇av片p| 国产极品精品免费视频能看的| 亚洲成人精品中文字幕电影| 美女黄网站色视频| 69人妻影院| 国产欧美日韩一区二区精品| 淫秽高清视频在线观看| 欧美黄色片欧美黄色片| 亚洲18禁久久av| 少妇裸体淫交视频免费看高清| 国产亚洲精品久久久com| 18禁黄网站禁片午夜丰满| 麻豆国产av国片精品| 人妻久久中文字幕网| 欧美另类亚洲清纯唯美| 两人在一起打扑克的视频| 色视频www国产| 亚洲人成伊人成综合网2020| 国产精品98久久久久久宅男小说| 欧美激情久久久久久爽电影| 国产高清视频在线播放一区| 亚洲精品亚洲一区二区| 五月伊人婷婷丁香| 在线播放国产精品三级| 午夜精品一区二区三区免费看| 欧美日韩瑟瑟在线播放| 国产美女午夜福利| 欧美激情久久久久久爽电影| 国产大屁股一区二区在线视频| 日日摸夜夜添夜夜添小说| 三级男女做爰猛烈吃奶摸视频| 欧美日本亚洲视频在线播放| 99在线人妻在线中文字幕| 99riav亚洲国产免费| 成年女人毛片免费观看观看9| 麻豆一二三区av精品| 久久精品91蜜桃| 久久亚洲精品不卡| 国产精品女同一区二区软件 | 国产精品女同一区二区软件 | 日韩欧美免费精品| 国产伦一二天堂av在线观看| 观看免费一级毛片| 国产成人a区在线观看| 久久久久国产精品人妻aⅴ院| 窝窝影院91人妻| 欧美最新免费一区二区三区 | 动漫黄色视频在线观看| 日韩中字成人| or卡值多少钱| 久99久视频精品免费| 午夜久久久久精精品| 伦理电影大哥的女人| 欧美日本亚洲视频在线播放| 亚洲无线在线观看| 男人和女人高潮做爰伦理| 一卡2卡三卡四卡精品乱码亚洲| 国产在线男女| 免费在线观看影片大全网站| 日韩欧美三级三区| 一区二区三区免费毛片| www.熟女人妻精品国产| 精品一区二区三区人妻视频| 亚洲av日韩精品久久久久久密| 日韩有码中文字幕| 精品久久久久久成人av| 国产视频一区二区在线看| 国产av不卡久久| 夜夜夜夜夜久久久久| 又紧又爽又黄一区二区| 18禁在线播放成人免费| 国产av一区在线观看免费| 一区二区三区激情视频| 国产精品精品国产色婷婷| 校园春色视频在线观看| 国产三级中文精品| www.www免费av| 波多野结衣巨乳人妻| 3wmmmm亚洲av在线观看| 精品熟女少妇八av免费久了| av天堂中文字幕网| 久久久精品欧美日韩精品| 中文字幕免费在线视频6| 成人av一区二区三区在线看| 观看美女的网站| 亚洲男人的天堂狠狠| 精品国产亚洲在线| 俄罗斯特黄特色一大片| 亚洲欧美日韩高清在线视频| 最新中文字幕久久久久| 精品熟女少妇八av免费久了| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 国产黄a三级三级三级人| 国产成人欧美在线观看| 99久久久亚洲精品蜜臀av| 日日干狠狠操夜夜爽| 亚洲av成人精品一区久久| 午夜免费成人在线视频| 国产黄色小视频在线观看| 亚洲av成人精品一区久久| 少妇人妻一区二区三区视频| 欧美黑人巨大hd| 麻豆一二三区av精品| 久久婷婷人人爽人人干人人爱| 欧美色视频一区免费| 欧美激情国产日韩精品一区| 有码 亚洲区| 亚洲美女搞黄在线观看 | 色视频www国产| 熟女电影av网| 超碰av人人做人人爽久久| 老女人水多毛片| 国产三级在线视频| 国产美女午夜福利| 亚洲中文字幕一区二区三区有码在线看| 美女xxoo啪啪120秒动态图 | 日本a在线网址| 亚洲在线观看片| 深夜a级毛片| а√天堂www在线а√下载| 国产久久久一区二区三区| 久久99热6这里只有精品| 人人妻人人澡欧美一区二区| 欧美黑人欧美精品刺激| 色精品久久人妻99蜜桃| 国产aⅴ精品一区二区三区波| 午夜福利在线观看免费完整高清在 | 婷婷丁香在线五月| 白带黄色成豆腐渣| 日韩欧美一区二区三区在线观看| 欧美xxxx黑人xx丫x性爽| 直男gayav资源| 亚洲国产高清在线一区二区三| 久久久久久久久久黄片| 熟女人妻精品中文字幕| 中文字幕久久专区| ponron亚洲| 综合色av麻豆| 在线看三级毛片| 欧美乱色亚洲激情| 精品久久久久久久久亚洲 | 少妇被粗大猛烈的视频| 亚洲,欧美,日韩| 欧美日韩综合久久久久久 | 日本一本二区三区精品| 精品国产亚洲在线| 天堂动漫精品| 美女高潮喷水抽搐中文字幕| 国产一区二区三区在线臀色熟女| 久久久久九九精品影院| 国产精品伦人一区二区| 在线十欧美十亚洲十日本专区| 丰满的人妻完整版| 美女cb高潮喷水在线观看| 国产亚洲精品久久久com| 淫秽高清视频在线观看| 欧美丝袜亚洲另类 | 99精品在免费线老司机午夜| 中文字幕人成人乱码亚洲影| 麻豆成人午夜福利视频| 日日夜夜操网爽| 九九在线视频观看精品| 日本 av在线| 色综合站精品国产| 国产亚洲精品久久久com| 色尼玛亚洲综合影院| 人妻夜夜爽99麻豆av| 免费电影在线观看免费观看| 国产精品嫩草影院av在线观看 | 蜜桃久久精品国产亚洲av| 欧美日韩综合久久久久久 | 亚洲成人中文字幕在线播放| 嫩草影院入口| 午夜精品一区二区三区免费看| 亚洲精品一区av在线观看| 国产精品久久久久久人妻精品电影| 中文字幕免费在线视频6| 永久网站在线| 国产私拍福利视频在线观看| 高清毛片免费观看视频网站| 观看美女的网站| 在线免费观看不下载黄p国产 | 欧美高清成人免费视频www| 欧美日韩黄片免| 极品教师在线视频| 91狼人影院| 亚洲国产精品999在线| 成人午夜高清在线视频| 高清毛片免费观看视频网站| 欧美潮喷喷水| 亚洲国产精品sss在线观看| 在线观看午夜福利视频| av在线老鸭窝| 亚洲中文字幕日韩| 国产精品98久久久久久宅男小说| 天堂√8在线中文| 亚洲国产精品sss在线观看| 热99在线观看视频| 国产老妇女一区| 国产成年人精品一区二区| 麻豆国产av国片精品| 看片在线看免费视频| 日本免费一区二区三区高清不卡| 噜噜噜噜噜久久久久久91| 亚洲av电影不卡..在线观看| 欧美日韩综合久久久久久 | 亚洲国产精品久久男人天堂| 少妇高潮的动态图| 国产成年人精品一区二区| 成年女人毛片免费观看观看9| 欧美日韩黄片免| а√天堂www在线а√下载| 久久精品国产亚洲av天美| 国产精品98久久久久久宅男小说| 国产精品三级大全| 特大巨黑吊av在线直播| 欧美区成人在线视频| 午夜福利免费观看在线| 国产麻豆成人av免费视频| 麻豆一二三区av精品| 2021天堂中文幕一二区在线观| 日韩欧美精品免费久久 | 亚洲,欧美精品.| 可以在线观看的亚洲视频| 内射极品少妇av片p| 国产精品一及| 最近视频中文字幕2019在线8| 91在线精品国自产拍蜜月| 少妇的逼好多水| 男人的好看免费观看在线视频| 十八禁国产超污无遮挡网站| 久久久久久九九精品二区国产| 日日干狠狠操夜夜爽| 免费人成在线观看视频色| 午夜福利在线在线| 日本 欧美在线| 在线观看舔阴道视频| 久久人人爽人人爽人人片va | 男女床上黄色一级片免费看| 麻豆国产97在线/欧美| 亚洲人成伊人成综合网2020| 搡老熟女国产l中国老女人| 久久精品久久久久久噜噜老黄 | av在线观看视频网站免费| 毛片一级片免费看久久久久 | 最近最新中文字幕大全电影3| 免费在线观看影片大全网站| 一个人免费在线观看电影| 天堂影院成人在线观看| 999久久久精品免费观看国产| 久久精品夜夜夜夜夜久久蜜豆| 日韩亚洲欧美综合| 人妻制服诱惑在线中文字幕| 亚洲熟妇中文字幕五十中出| aaaaa片日本免费| 精品午夜福利视频在线观看一区| 午夜福利在线观看免费完整高清在 | 窝窝影院91人妻| 在线a可以看的网站| 国产成人福利小说| 能在线免费观看的黄片| 在线十欧美十亚洲十日本专区| 97超视频在线观看视频| xxxwww97欧美| 网址你懂的国产日韩在线| 最新中文字幕久久久久| 国产真实乱freesex| 伊人久久精品亚洲午夜| 欧美黄色片欧美黄色片| 精品一区二区三区视频在线| 国产精品电影一区二区三区| 最近最新中文字幕大全电影3| 午夜福利欧美成人| 51国产日韩欧美| 成年女人看的毛片在线观看| 国产一区二区三区在线臀色熟女| 老女人水多毛片| av国产免费在线观看| 国产美女午夜福利| 又紧又爽又黄一区二区| 少妇裸体淫交视频免费看高清| 免费在线观看日本一区| av欧美777| 精品无人区乱码1区二区| 成年女人看的毛片在线观看| 久久久久久久精品吃奶| 国产亚洲精品av在线| 十八禁网站免费在线| 国产69精品久久久久777片| av在线蜜桃| 欧美成人一区二区免费高清观看| 特级一级黄色大片| 国产一区二区亚洲精品在线观看| 午夜福利成人在线免费观看| 国产爱豆传媒在线观看| 午夜福利在线观看吧| 国产伦精品一区二区三区四那| 久久精品国产清高在天天线| 久久久国产成人精品二区| 日韩人妻高清精品专区| 最近最新免费中文字幕在线| 欧美绝顶高潮抽搐喷水| 午夜福利免费观看在线| 高潮久久久久久久久久久不卡| avwww免费| 97超级碰碰碰精品色视频在线观看| 国产免费一级a男人的天堂| 欧美成人一区二区免费高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久国产a免费观看| 内地一区二区视频在线| 久久国产乱子伦精品免费另类| 别揉我奶头~嗯~啊~动态视频| 亚洲av五月六月丁香网| 中文字幕熟女人妻在线| 99精品在免费线老司机午夜| 国产精品免费一区二区三区在线| 成年免费大片在线观看| 国产精品一区二区性色av| 啪啪无遮挡十八禁网站| 欧美+日韩+精品| a级毛片a级免费在线| 亚洲在线观看片| 香蕉av资源在线| 成熟少妇高潮喷水视频| 无遮挡黄片免费观看| 桃色一区二区三区在线观看| 91午夜精品亚洲一区二区三区 | 久久精品综合一区二区三区| 精品人妻一区二区三区麻豆 | 成人无遮挡网站| 熟妇人妻久久中文字幕3abv| 女生性感内裤真人,穿戴方法视频| 日韩欧美国产在线观看| 赤兔流量卡办理| 国产免费av片在线观看野外av| 又爽又黄无遮挡网站| 91午夜精品亚洲一区二区三区 | 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| 国产av不卡久久| 真人一进一出gif抽搐免费| 国产黄片美女视频| 日本 欧美在线| 久久久久免费精品人妻一区二区| 国产欧美日韩精品亚洲av| 日韩欧美在线二视频| 久久精品人妻少妇| 真人一进一出gif抽搐免费| 色播亚洲综合网| 久久久久久久久久成人| 哪里可以看免费的av片| 久久人妻av系列| 精品久久久久久久久亚洲 | 欧美乱妇无乱码| 三级男女做爰猛烈吃奶摸视频| 岛国在线免费视频观看| 国产午夜精品论理片| 99在线人妻在线中文字幕| 国产高清三级在线| 高清日韩中文字幕在线| 国产又黄又爽又无遮挡在线| 男女下面进入的视频免费午夜| 日本黄色片子视频| 身体一侧抽搐| 精品久久久久久久末码| 伦理电影大哥的女人| 少妇裸体淫交视频免费看高清| 在线播放无遮挡| 天堂动漫精品| 日本一本二区三区精品| 亚洲激情在线av| 欧美高清性xxxxhd video| 日韩中字成人| 免费看日本二区| 人人妻,人人澡人人爽秒播| 欧美日韩福利视频一区二区| 人妻夜夜爽99麻豆av| 日本一二三区视频观看| 欧美一区二区精品小视频在线| 男人狂女人下面高潮的视频| 亚洲综合色惰| 国产色婷婷99| 国产精品日韩av在线免费观看| 国内精品久久久久精免费| 亚洲无线观看免费| 午夜影院日韩av| 99久久精品国产亚洲精品| 精品一区二区三区视频在线观看免费| 成人无遮挡网站| 中文亚洲av片在线观看爽| 97热精品久久久久久| 黄色日韩在线| 天美传媒精品一区二区| 赤兔流量卡办理| 综合色av麻豆| 亚洲成av人片在线播放无| 国产国拍精品亚洲av在线观看| 亚洲精品乱码久久久v下载方式| 婷婷六月久久综合丁香| 欧美日韩亚洲国产一区二区在线观看| 自拍偷自拍亚洲精品老妇| 色视频www国产| 午夜福利高清视频| 国产私拍福利视频在线观看| 乱码一卡2卡4卡精品| 男女做爰动态图高潮gif福利片| 亚洲不卡免费看| 18禁黄网站禁片免费观看直播| 美女高潮的动态| 97人妻精品一区二区三区麻豆| 在线免费观看不下载黄p国产 | 脱女人内裤的视频| 国产高清视频在线观看网站| 校园春色视频在线观看| 国产真实伦视频高清在线观看 | 日韩高清综合在线| av天堂中文字幕网| 亚洲av第一区精品v没综合| 极品教师在线视频| 国产三级中文精品| 一级黄片播放器| 精品人妻1区二区| aaaaa片日本免费| 久久欧美精品欧美久久欧美| 中文字幕精品亚洲无线码一区| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 国产av在哪里看| 欧美黄色淫秽网站| 91av网一区二区| 校园春色视频在线观看| 日本 av在线| 精品久久久久久久久久久久久| 国内揄拍国产精品人妻在线| 欧美国产日韩亚洲一区| 国产欧美日韩一区二区三| 九色国产91popny在线| 一进一出好大好爽视频| 久久国产精品人妻蜜桃| 天堂动漫精品| 看片在线看免费视频| 国产一区二区三区在线臀色熟女| 国产亚洲精品av在线| 一个人免费在线观看的高清视频| 国产伦一二天堂av在线观看| 日日干狠狠操夜夜爽| 又紧又爽又黄一区二区| 99热这里只有是精品在线观看 | 欧美高清成人免费视频www| 亚洲精品久久国产高清桃花| 午夜福利18| 国产高清三级在线| 亚洲av熟女| 欧美+亚洲+日韩+国产| 午夜久久久久精精品| 午夜亚洲福利在线播放| 国产一区二区三区视频了| 成熟少妇高潮喷水视频| a级毛片免费高清观看在线播放| 精品不卡国产一区二区三区| 国产精品久久久久久人妻精品电影| 成人国产一区最新在线观看| 精品日产1卡2卡| 国产成年人精品一区二区| 又黄又爽又免费观看的视频| 美女 人体艺术 gogo| 免费av毛片视频| 国产免费男女视频| 久99久视频精品免费| 亚洲综合色惰| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 91午夜精品亚洲一区二区三区 | 国产一区二区三区视频了| 色av中文字幕| 亚洲狠狠婷婷综合久久图片| 午夜两性在线视频| 色综合欧美亚洲国产小说| 色哟哟·www| 亚洲av二区三区四区| 很黄的视频免费| 亚洲激情在线av| 亚洲欧美日韩东京热| 三级毛片av免费| 亚洲精品色激情综合| 国产午夜精品论理片| 国产真实乱freesex| 久久精品国产清高在天天线| 日本熟妇午夜| 亚洲成人久久爱视频| 丰满人妻一区二区三区视频av| 又黄又爽又免费观看的视频| 无人区码免费观看不卡| 丰满人妻一区二区三区视频av| 国产毛片a区久久久久| 国产日本99.免费观看| 久久精品影院6| 国产老妇女一区| 亚洲av不卡在线观看| 乱人视频在线观看|