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

摘要：使用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