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

關(guān)鍵詞：背景噪聲；西昌臺(tái)陣；安寧河—?jiǎng)t木河斷裂帶；波速變化

中圖分類號(hào)：P315.2文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2024）04-0493-10

doi：10.20015/j.cnki.ISSN1000-0666.2024.0045

0引言

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

主動(dòng)源探測(cè)技術(shù)（Vidale，Li，2003；Wang et al，2008）和重復(fù)地震的相似波形（Peng，Ben-Zion，2006；Schaff，Beroza，2004）實(shí)現(xiàn)了對(duì)地下時(shí)變的監(jiān)測(cè)，但它們?cè)跁r(shí)空上的不均勻性使得連續(xù)長(zhǎng)期的地震波速監(jiān)測(cè)依舊困難重重。近年來(lái)，得益于地震背景噪聲技術(shù)的發(fā)展，使用背景噪聲互相關(guān)得到的相似波形推動(dòng)了連續(xù)監(jiān)測(cè)地震波速變化技術(shù)的快速發(fā)展。已有很多研究利用這一技術(shù)來(lái)監(jiān)測(cè)地下介質(zhì)的時(shí)間變化，包括地下波速的季節(jié)性變化（Sens-Schnfelder，Wegler et al，2006；Wang et al，2017）、與火山噴發(fā)相關(guān)的波速變化以及與強(qiáng)震相關(guān)的波速變化等（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年安寧河—?jiǎng)t木河斷裂帶周邊地區(qū)23個(gè)寬頻帶地震儀記錄的連續(xù)波形數(shù)據(jù)，利用背景噪聲互相關(guān)技術(shù)計(jì)算了各臺(tái)站對(duì)之間的波速變化及時(shí)空分布特征，討論了波速變化與信號(hào)頻帶之間的聯(lián)系、波速變化的可能機(jī)理和環(huán)境氣候因素等對(duì)介質(zhì)變化的影響。

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

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

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

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

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

2結(jié)果分析

2.1不同周期的臺(tái)站對(duì)波速變化

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

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

2.2波速變化的空間分布

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

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

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

3討論

3.1降雨對(duì)波速變化的影響

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

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

3.2地震對(duì)波速變化的影響

大地震的發(fā)生往往伴隨震源區(qū)的破壞和地表破碎帶的顯現(xiàn)?；诒尘霸肼暤恼鹪磪^(qū)介質(zhì)變化研究為探測(cè)地殼深部應(yīng)力的積累與釋放提供了獨(dú)特的視角，發(fā)現(xiàn)了大量同震波速降低和震后快速恢復(fù)的證據(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）。研究認(rèn)為地震導(dǎo)致的波速變化的主要機(jī)理為淺部的動(dòng)態(tài)震動(dòng)（Sleep，2015）和深部介質(zhì)應(yīng)力狀態(tài)的改變。

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

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

3.3波速變化與測(cè)量周期關(guān)系

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

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

4結(jié)論

本文利用2013—2018年西昌臺(tái)陣記錄的連續(xù)波形資料，通過(guò)背景噪聲互相關(guān)計(jì)算了安寧河—?jiǎng)t木河斷裂帶及周邊地區(qū)23個(gè)臺(tái)站之間的相對(duì)波速變化，得到如下主要結(jié)論：

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

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

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

參考文獻(xiàn)：

馮德益.1981.地震波速異常［M］.北京：地震出版社.

Feng D Y.1981.Velocity anomalies of seismic waves［M］.Beijing：Siesmological Press.（in Chinese）

胡淑芳，李軍，戴宗輝，等.2021.利用背景噪聲研究福建金鐘庫(kù)區(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.利用背景噪聲互相關(guān)研究汶川地震震源區(qū)地震波速度變化［J］.地球物理學(xué)報(bào)，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］.地球物理學(xué)報(bào)，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）

謝凡，夏開(kāi)文，黃會(huì)寶，等.2020.基于多重散射波波速變化的滑坡實(shí)時(shí)監(jiān)測(cè)方法與應(yīng)用研究［J］.巖石力學(xué)與工程學(xué)報(bào)，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ū)地殼速度相對(duì)變化的環(huán)境噪聲自相關(guān)研究［J］.地球物理學(xué)報(bào)，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