汪 建 向 涯 楊亞運(yùn) 李翠平

1）中國(guó)重慶 400020 重慶市地震局

2）中國(guó)武漢 430071 中國(guó)地震局地震研究所

3）中國(guó)武漢 430071 中國(guó)地震局地震大地測(cè)量重點(diǎn)實(shí)驗(yàn)室

0 引言

地震波在地球介質(zhì)中經(jīng)多次散射形成尾波。尾波在地球介質(zhì)中傳播時(shí)間長(zhǎng)，由于重復(fù)采樣、疊加和放大作用，其對(duì)介質(zhì)的微小變化更加敏感，能識(shí)別直達(dá)波不能識(shí)別的介質(zhì)的細(xì)微變化（張金川等，2014）?；谠撎匦?，可利用尾波干涉方法識(shí)別震源、散射體位置、地震波速度與介質(zhì)的細(xì)微變化，并廣泛應(yīng)用于斷層活動(dòng)性、工程材料檢測(cè)、冰川活動(dòng)、火山監(jiān)測(cè)等領(lǐng)域（肖卓等，2015）。近年來(lái)，諸多地震學(xué)者利用尾波干涉方法開(kāi)展了一系列研究，如：宋麗莉等（2012）通過(guò)室內(nèi)試驗(yàn)，利用尾波干涉方法觀測(cè)到直達(dá)波難以分辨的巖石物性變化；張金川（2014）利用尾波干涉技術(shù)研究廣西龍灘水庫(kù)重復(fù)地震，發(fā)現(xiàn)散射體運(yùn)移可能與庫(kù)水向地下裂隙中的運(yùn)移有關(guān)；肖卓等（2016）運(yùn)用尾波干涉技術(shù)分析2014 年云南盈江雙震期間地殼介質(zhì)狀態(tài)的變化；李樂(lè)等（2017）利用重復(fù)地震尾波研究汶川強(qiáng)震破裂區(qū)震后波速的微變化；鮑子文（2019）利用重復(fù)地震尾波研究安徽霍山地區(qū)介質(zhì)變化，發(fā)現(xiàn)地震波速度存在季節(jié)性變化，與大氣降水有較好的相關(guān)性。

2017 年11 月23 日17 時(shí)43 分，重慶武隆區(qū)發(fā)生MS5.0 地震。該地震事件是有地震儀器記錄以來(lái)重慶渝東南地區(qū)發(fā)生的最大地震。文中利用尾波干涉原理，研究武隆MS5.0地震發(fā)生后，余震期間震源區(qū)地下介質(zhì)的變化，有助于了解強(qiáng)震發(fā)生時(shí)地殼介質(zhì)性質(zhì)和應(yīng)力狀態(tài)的變化。

1 資料選取

武隆區(qū)地處重慶市東南部，位于烏江下游峽谷區(qū)，巖石成分以碳酸鹽、砂巖和頁(yè)巖為主，巖溶地貌發(fā)育，形成大量溶洞群、天坑群、地縫等自然景觀。該地區(qū)附近分布方斗山基底斷裂、七曜山—金佛山基底斷裂及彭水基底斷裂，中小斷裂有三會(huì)沖斷裂、芙蓉江斷裂、文復(fù)斷裂、馬武斷裂、火石埡斷裂、郁山斷裂、和尚巖斷裂（戴應(yīng)洪等，2020）（圖1）。李翠平等（2019）分析認(rèn)為，震中附近流體的深部滲透作用導(dǎo)致局部斷層活化，誘發(fā)此次武隆MS5.0 地震，發(fā)震斷層與文復(fù)斷裂關(guān)系較為密切。

圖1 研究區(qū)地震臺(tái)站、震中及區(qū)域構(gòu)造Fig.1 Seismic stations,epicenters,and regional structure in the study area

武隆MS5.0 地震發(fā)生后，在震中附近架設(shè)3 套流動(dòng)臺(tái)（L5541、L5550、L5551），均配備短周期地震計(jì)（采樣率100 Hz），進(jìn)行地震監(jiān)測(cè)。各流動(dòng)臺(tái)臺(tái)基噪聲測(cè)試參數(shù)見(jiàn)表1，可見(jiàn)環(huán)境地噪聲水平均在Ⅱ級(jí)或以上。截至12 月23 日，重慶數(shù)字地震臺(tái)網(wǎng)共記錄到武隆及鄰區(qū)（29°00′—29°42′N(xiāo)，107°12′—108°24′E）發(fā)生ML≥1.0 地震108 次（含MS5.0 主震），最大余震為2017 年11 月29 日10 時(shí)34 分武隆ML3.1 地震。選取108 次ML≥1.0 地震波形數(shù)據(jù)和觀測(cè)報(bào)告，分析地殼介質(zhì)變化。

表1 流動(dòng)臺(tái)站臺(tái)基噪聲參數(shù)Table 1 Noise parameters of three mobile stations

2 重復(fù)地震挑選

重復(fù)地震概念于20 世紀(jì)60 年代被提出，目前尚無(wú)統(tǒng)一定義（李宇彤等，2011），Nadeau 等（1995）將發(fā)生在同一斷層位置，震級(jí)和復(fù)發(fā)間隔相近，震源機(jī)制與波形高度相似的一組地震稱(chēng)為重復(fù)地震。Schaff等（2004）則將重復(fù)地震定義為被至少一個(gè)臺(tái)站記錄到，且波形相關(guān)系數(shù)不小于0.8的地震對(duì)事件。文中通過(guò)計(jì)算波形互相關(guān)系數(shù)（cross-correlation簡(jiǎn)稱(chēng)cc）和約束震源位置的方法，來(lái)挑選武隆MS5.0 地震發(fā)生后的重復(fù)地震。

2.1 地震序列重定位

汪建等（2020）利用武隆地區(qū)2010—2017 年重復(fù)地震和觀測(cè)報(bào)告資料，發(fā)現(xiàn)該區(qū)地震定位和震相拾取精度較高。為得到更精確的震源位置，利用雙差定位法，采用武隆地區(qū)精細(xì)地殼速度結(jié)構(gòu)模型（表2）（王小龍等，2013），對(duì)所選108 次地震序列進(jìn)行重定位。雙差定位結(jié)果顯示，MS5.0 主震發(fā)生在三會(huì)沖、芙蓉江、文復(fù)斷裂與馬武斷裂交匯部位，震中位于烏江上游彭水電站與下游銀盤(pán)電站水位段之間，與烏江水系最近距離約2 km，震源深度約10 km，余震序列呈條帶分布，與附近斷裂走向平行，地震序列破裂面呈NE—SW向單側(cè)破裂（圖1）。

表2 武隆地區(qū)速度結(jié)構(gòu)模型（據(jù)王小龍等，2013）Table 2 The seismic velocity model used in the location (According to Wang Xiaolong et al,2013)

2.2 重復(fù)地震挑選

將地震波形轉(zhuǎn)換為SAC 格式，選取垂直分量，進(jìn)行去傾斜、去線(xiàn)性趨勢(shì)及0.5—10 Hz 帶通濾波預(yù)處理（汪建等，2018），挑選重復(fù)地震，步驟如下：①計(jì)算同一臺(tái)站中任意2個(gè)地震波形的互相關(guān)系數(shù)，計(jì)算起點(diǎn)為P 波到時(shí)前1 s，計(jì)算長(zhǎng)度選擇4 倍Sg、Pg 震相走時(shí)差；②找出同時(shí)被3 個(gè)臺(tái)站（L5541、L5550、L5551）記錄且各臺(tái)互相關(guān)系數(shù)（cc）≥0.9的地震對(duì)事件作為相似地震；③結(jié)合雙差定位結(jié)果，尋找間距（不考慮震源深度的影響）小于1 km 的相似地震。

經(jīng)挑選，得到3 組滿(mǎn)足條件的地震對(duì)，均包含2 個(gè)地震事件。重復(fù)地震以D1、D2、D3編號(hào)，主震設(shè)為D0。重復(fù)地震與流動(dòng)臺(tái)站分布見(jiàn)圖2，可見(jiàn)D1、D2、D3均分布在主震D0南側(cè)，由北向南依次擴(kuò)展，結(jié)合發(fā)震時(shí)刻可知，重復(fù)地震對(duì)均發(fā)生在余震活動(dòng)相對(duì)減少階段。

圖2 臺(tái)站與重復(fù)地震分布Fig.2 Distribution of seismic stations and repeating earthquakes

對(duì)武隆MS5.0 地震序列重復(fù)地震對(duì)參數(shù)信息進(jìn)行統(tǒng)計(jì)，結(jié)果見(jiàn)表3，可知：D1、D2、D3的震級(jí)均小于ML2.0，震級(jí)差值最大ML0.5，最小ML0.3；震源深度差值最大2 km，最小1 km；地震對(duì)重復(fù)周期具有明顯變化，時(shí)間間隔從十幾小時(shí)至幾天不等，最大間隔近4 天，D1、D2時(shí)間間隔相近。

表3 武隆MS 5.0 地震及重復(fù)地震對(duì)的參數(shù)統(tǒng)計(jì)Table 3 Parameter statistics of the Wulong MS 5.0 earthquake and repeating earthquakes

3 尾波干涉分析

3.1 尾波干涉原理

Snieder（2002）基于前人研究提出尾波干涉方法，并對(duì)其原理進(jìn)行了系統(tǒng)闡述。尾波干涉法基于路徑疊加理論，認(rèn)為波場(chǎng)是沿可能散射路徑傳播的子波的疊加，重復(fù)地震與臺(tái)站之間的傳播路徑大致相同，波形差異由介質(zhì)變化所引起，重復(fù)地震尾波的變化反映了以臺(tái)站和重復(fù)地震為焦點(diǎn)的橢球體上介質(zhì)的變化，可以利用重復(fù)地震尾波研究地殼介質(zhì)的波速變化，進(jìn)而分析引起波速變化的原因（肖卓等，2016）。

波形數(shù)據(jù)原始采樣率為100 Hz，即兩列波形走時(shí)延遲最小為0.01 s。為獲得更小的走時(shí)延遲，在尾波干涉測(cè)量之前，采用插值法將原始信號(hào)的采樣率提升至10 000 Hz（Peng et al，2006），即兩列波形的走時(shí)延遲最小為0.000 1 s。當(dāng)波速在介質(zhì)中發(fā)生變化時(shí)，重復(fù)地震的尾波穿過(guò)該介質(zhì)時(shí)也將發(fā)生走時(shí)擾動(dòng)，在接收點(diǎn)表現(xiàn)為波形的走時(shí)偏移。文中采用移動(dòng)窗口互相關(guān)法（Niu et al，2003；肖卓等，2016）計(jì)算重復(fù)地震之間的走時(shí)偏移，即將發(fā)震時(shí)刻較早的重復(fù)地震當(dāng)作參考地震，從P 波初至到時(shí)起算，對(duì)重復(fù)地震進(jìn)行尾波干涉分析。通過(guò)移動(dòng)時(shí)間窗，計(jì)算每一個(gè)時(shí)間窗內(nèi)波形的走時(shí)延遲，得到走時(shí)延遲及不相關(guān)系數(shù)隨流逝時(shí)間的分布。結(jié)合前人經(jīng)驗(yàn)（張金川，2014；肖卓等，2016），選取P 波初至到時(shí)至S 波到時(shí)的2 倍加4 s 作為尾波干涉測(cè)量窗口，移動(dòng)窗口長(zhǎng)度為1 s，移動(dòng)步長(zhǎng)為0.05 s，即每1 s 內(nèi)有20 個(gè)干涉測(cè)量點(diǎn)。

假設(shè)介質(zhì)中波速變化是均勻的，若波速下降（增大），則兩列波形的走時(shí)偏移將隨著流逝時(shí)間線(xiàn)性上升（減?。⊿nieder，2006）。設(shè)原始波速為v，介質(zhì)整體發(fā)生大小為δv的波速變化，則走時(shí)偏移τ與波速變化δv的對(duì)應(yīng)關(guān)系（Snieder，2006；Payan et al，2011）為

式中，走時(shí)偏移τ與流逝時(shí)間t經(jīng)線(xiàn)性擬合得到的直線(xiàn)斜率即為介質(zhì)波速的相對(duì)變化。在計(jì)算每小時(shí)時(shí)間窗j內(nèi)波形的走時(shí)偏移τj時(shí)，其測(cè)量誤差估計(jì)（Liu et al，2010）為

式中，tj為小時(shí)間窗j的中心位置，M為時(shí)間窗的數(shù)量，ητ為擬合殘差，表示為

測(cè)量誤差下限στ需滿(mǎn)足Cramer-Rao Lower Bound 法則（Quazi，1981；Carter，1987），公式如下

式中，f0為信號(hào)主頻，B為信號(hào)的頻寬與主頻之比，T為窗口長(zhǎng)度，ρ為波形相關(guān)系數(shù)，SNR 為信噪比。f0和B由波形數(shù)據(jù)控制，測(cè)量誤差主要取決于互相關(guān)系數(shù)與信噪比，選取互相關(guān)系數(shù)及信噪比較高的尾波段進(jìn)行尾波干涉測(cè)量較為合理（肖卓等，2016）。

3.2 尾波干涉結(jié)果

利用上述原理，采用尾波干涉技術(shù)，分別對(duì)重復(fù)地震對(duì)D1、D2、D3在3 個(gè)流動(dòng)臺(tái)的波形記錄進(jìn)行尾波干涉處理。

3.2.1 L5541 臺(tái)站記錄尾波干涉處理。分別對(duì)重復(fù)地震對(duì)D1、D2、D3在L5541 臺(tái)的記錄進(jìn)行尾波干涉處理，結(jié)果見(jiàn)圖3、圖4、圖5（0 s 設(shè)為P 波初至到時(shí)）。

圖3 L5541 臺(tái)站D1 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.3 The seismic waveforms,uncorrelation coefficient and delay time of D1 repeating earthquakes with time at the L5541 station

圖4 L5541 臺(tái)站D2 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.4 The seismic waveforms,uncorrelation coefficient and delay time of D2 repeating earthquakes with time at the L5541 station

圖5 L5541 臺(tái)站D3 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.5 The seismic waveforms,uncorrelation coefficient and delay time of D3 repeating earthquake with time at the L5541 station

（1）重復(fù)地震對(duì)D1。D1的2 次地震分別發(fā)生在主震D0后第7 和第8 天，分析L5541 臺(tái)站的2 次地震預(yù)處理波形及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間（從P 波初至到時(shí)起算）的變化，結(jié)果見(jiàn)圖3。重復(fù)地震對(duì)D1的2 次地震事件震級(jí)相差不大（0.5），但其相關(guān)系數(shù)卻達(dá)0.937。由圖3 可見(jiàn)：①在P 波后續(xù)部分（1.00—1.95 s），走時(shí)延遲基本保持不變，后迅速下降至零值附近波動(dòng)；②在S 波振幅較大部分（1.95—3.20 s），走時(shí)延遲緩慢下降，S 波后續(xù)部分（3.20—4.70 s）的走時(shí)延遲呈先增加后減小的變化趨勢(shì)；③隨著流逝時(shí)間的增加，S 波早期尾波段（5.25—6.75 s）走時(shí)延遲明顯增加，意味著波速的減小，S 波尾波后續(xù)部分（6.75—7.75 s）的走時(shí)延遲明顯減小，最大走時(shí)延遲約65×10-3s；④結(jié)合移動(dòng)窗口內(nèi)波形的信噪比、不相關(guān)系數(shù)和走時(shí)延遲的線(xiàn)性變化規(guī)律，選取S 波早期尾波段（5.25—6.75 s），由式（1）計(jì)算得到，其相對(duì)波速變化約為2.1‰。

（2）重復(fù)地震對(duì)D2。D2的2 次地震分別發(fā)生在主震D0后第8 和第9 天，圖4 給出重復(fù)地震對(duì)D2在L5541 臺(tái)站的尾波干涉結(jié)果。由圖4 可見(jiàn)：①在P 波后續(xù)部分（1.00—1.55 s），走時(shí)延遲迅速下降至零值附近波動(dòng)；②在S 波振幅最大部分（1.55—2.40 s），走時(shí)延遲基本保持不變，隨著流逝時(shí)間的增加，S 波后續(xù)部分（2.50—3.75 s）的走時(shí)延遲有緩慢增加的趨勢(shì)；③在S 波早期尾波段（3.80—4.90 s），走時(shí)延遲明顯增加，變化較為線(xiàn)性，S 波尾波后續(xù)部分（4.90—6.35 s）的走時(shí)延遲明顯減小，最大走時(shí)延遲約19.1×10-3s；④選取S 波早期尾波段（3.80—4.90 s）計(jì)算其相對(duì)波速變化，結(jié)果約為13.4‰，波速變化較大。

重復(fù)地震對(duì)D2距L5541 臺(tái)站最近，傳播路徑不穿過(guò)烏江流域，可能受到斷層擠壓作用，巖石裂隙增多，導(dǎo)致地下介質(zhì)劇烈變化，波速下降值較大。若選取S 波后續(xù)部分（2.50—3.75 s）計(jì)算其相對(duì)波速變化，結(jié)果約1.4‰，約為S 波早期尾波段計(jì)算結(jié)果的1/10，正是由于尾波的多次散射、疊加、放大，以及比S 波更敏感的特點(diǎn)，尾波識(shí)別的波速變化比S 波放大近10 倍。

（3）重復(fù)地震對(duì)D3。D3的2 次地震分別發(fā)生在主震D0后第9 和第14 天，圖5 給出重復(fù)地震對(duì)D3在L5541 臺(tái)站的尾波干涉結(jié)果。由圖5 可見(jiàn)：①S 波衰減較快，在S 波部分（2.35—3.25 s）走時(shí)延遲緩慢增加，在S波早期尾波段（3.25—4.60 s）走時(shí)延遲增加較為明顯，S 波后續(xù)尾波部分（4.60—7.35 s）走時(shí)延遲變化復(fù)雜，整體呈下降趨勢(shì)。選取S 波早期尾波段（3.25—4.60 s），計(jì)算得到其相對(duì)波速變化約為2.2‰；②在P波尾波部分（1.25—1.45 s）發(fā)現(xiàn)走時(shí)延遲和不相關(guān)系數(shù)均存在1 個(gè)短周期“尖峰”（spike）變化，且變化比較尖銳，相比于P 波尾波，S 波及其尾波變化更為平滑。

Niu 等（2003）在美國(guó)圣安德列斯斷層帕克菲爾德區(qū)域，對(duì)重復(fù)微震進(jìn)行數(shù)值模擬，發(fā)現(xiàn)若只是單一散射體的位置發(fā)生變化，走時(shí)擾動(dòng)與不相關(guān)系數(shù)均將出現(xiàn)一個(gè)短周期的“尖峰”變化，若震源位置或介質(zhì)波速發(fā)生變化，隨著時(shí)間的流逝，走時(shí)擾動(dòng)與不相關(guān)系數(shù)將呈現(xiàn)長(zhǎng)周期變化趨勢(shì)，根據(jù)實(shí)驗(yàn)結(jié)果，P 波尾波的“尖峰”變化可能與地下介質(zhì)中散射體的運(yùn)移有關(guān)。由圖2 可知，重復(fù)地震對(duì)D3到L5541 臺(tái)站的傳播路徑穿過(guò)烏江流域。該區(qū)巖溶發(fā)育，遍布溶洞和地下暗河，因S 波不能在液體中傳播，P 波尾波出現(xiàn)的“尖峰”變化可能與震后庫(kù)水在地下裂隙中的運(yùn)移有關(guān)。

3.2.2 L5550 臺(tái)站記錄尾波干涉處理。分別對(duì)重復(fù)地震對(duì)D1、D2、D3在L5541 臺(tái)的記錄進(jìn)行尾波干涉處理，可知D1在該臺(tái)站尾波干涉結(jié)果未見(jiàn)明顯變化，D2、D3的尾波干涉結(jié)果見(jiàn)圖6、圖7。

圖6 L5550 臺(tái)站D2 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.6 The seismic waveforms,uncorrelation coefficient and delay time of D2 repeating earthquakes with time at the L5550 station

圖7 L5550 臺(tái)站D3 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.7 The seismic waveforms,uncorrelation coefficient and delay time of D3 repeating earthquakes with time at the L5550 station

（1）重復(fù)地震對(duì)D2。圖6 給出重復(fù)地震對(duì)D2在L5550 臺(tái)站的尾波干涉結(jié)果。由圖6可見(jiàn)：①在P 波尾波部分（1.00—2.55 s）走時(shí)延遲無(wú)明顯變化；②走時(shí)延遲在S 波早期部分（2.55—4.00 s）減小，后續(xù)部分（4.00—6.05 s）增加，在S 波早期尾波段（6.05—7.25 s）明顯減小，意味著波速的增加，在S 波后續(xù)尾波部分（7.25—9.40 s）保持在零值附近波動(dòng)，最大走時(shí)延遲約為6.2×10-3s；③選取S波早期尾波段（6.05—7.25 s）計(jì)算其相對(duì)波速變化，結(jié)果約為6.7‰，可能是由于此段波形相關(guān)系數(shù)不高，由式（2）—（3）計(jì)算，測(cè)量誤差與波速變化在同一個(gè)數(shù)量級(jí)。

（2）重復(fù)地震對(duì)D3。圖7給出重復(fù)地震對(duì)D3在L5550臺(tái)站的尾波干涉結(jié)果。由圖7可見(jiàn)：①P 波尾波部分（1.00—2.70 s）的相關(guān)系數(shù)較低，走時(shí)延遲先增加后減??；②在S 波部分（2.70—5.65 s），走時(shí)延遲變化復(fù)雜，呈減小—增加—減小的變化趨勢(shì)，S 波早期尾波段（5.75—6.85 s）走時(shí)延遲明顯增加，變化較為線(xiàn)性，S 波后續(xù)尾波部分（6.85—8.25 s）走時(shí)延遲下降至零值保持不變；③在不同波段（P、P 尾波、S、S 尾波），走時(shí)延遲的變化趨勢(shì)不同且復(fù)雜，但在S 波早期尾波段（5.65—6.85 s），走時(shí)延遲有明顯增加，且較為線(xiàn)性，最大走時(shí)延遲約為9.5×10-3s，其相對(duì)波速變化約為6.1‰。

3.2.3 L5551 臺(tái)站記錄尾波干涉處理。分別對(duì)重復(fù)地震對(duì)D1、D2、D3在L5551 臺(tái)的記錄進(jìn)行尾波干涉處理，可知D3在該臺(tái)尾波干涉結(jié)果未見(jiàn)明顯變化，D1、D2的尾波干涉結(jié)果見(jiàn)圖8、圖9。

（1）重復(fù)地震對(duì)D1。圖8 給出重復(fù)地震對(duì)D1在L5551 臺(tái)站的尾波干涉結(jié)果。由圖8可見(jiàn)：①在P 波尾波部分（1.35—1.75 s），發(fā)現(xiàn)走時(shí)延遲和不相關(guān)系數(shù)均存在2 個(gè)短周期“尖峰”變化，S 波及其尾波的走時(shí)延遲也有數(shù)個(gè)“尖峰”變化，但其最大變化幅度約為P 波的1/3，根據(jù)Niu 等（2003）的實(shí)驗(yàn)結(jié)果，此現(xiàn)象可能與震源位置或介質(zhì)波速變化有關(guān)，若震源位置發(fā)生變化，直達(dá)波與尾波變化的量級(jí)應(yīng)基本一致（Grêt et al，2006），假設(shè)由于介質(zhì)波速變化導(dǎo)致P 波尾波、S 波及其尾波出現(xiàn)數(shù)個(gè)“尖峰”變化，則可能與地下介質(zhì)的不均勻性或復(fù)雜性有關(guān)；②在S 波及其尾波段（2.10—7.40 s），走時(shí)延遲整體呈上升趨勢(shì)；③選取S 波早期尾波段（4.00—6.05 s）計(jì)算其相對(duì)波速變化，結(jié)果約為2.1‰。

圖8 L5551 臺(tái)站D1 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.8 The seismic waveforms,uncorrelation coefficient,and delay time of D1 repeating earthquakes with time at the L5551 station

（2）重復(fù)地震對(duì)D2。圖9 給出重復(fù)地震對(duì)D2在L5551 臺(tái)站的尾波干涉結(jié)果。由圖9 可見(jiàn)：①P 波尾波部分（1.00—2.10 s）相關(guān)系數(shù)不高，走時(shí)延遲變化較為復(fù)雜；②S 波振幅較大部分（2.10—4.65 s）走時(shí)延遲基本保持不變，在S 波早期尾波段（4.65—5.95 s）走時(shí)延遲明顯減小，意味著波速的增加，在S 波后續(xù)尾波部分（5.95—7.50 s）走時(shí)延遲保持在零值附近波動(dòng)；③選取S 波早期尾波段（4.65—5.95 s）計(jì)算其相對(duì)波速變化，結(jié)果約為4.5‰。

圖9 L5551 臺(tái)站D2 波形記錄及其不相關(guān)系數(shù)和走時(shí)延遲隨流逝時(shí)間的變化Fig.9 The seismic waveforms,uncorrelation coefficient,and delay time of D2 repeating earthquakes with time at the L5551 station

4 結(jié)論與討論

綜上，絕大多數(shù)的測(cè)量誤差比波速變化測(cè)量值小1 個(gè)數(shù)量級(jí)，所得波速變化結(jié)果較為可靠。圖2 可見(jiàn)，在L5541 臺(tái)站（距震源區(qū)最近），3 組重復(fù)地震記錄的S 波尾波部分走時(shí)延遲均有明顯增加現(xiàn)象，可能意味著地震波速的持續(xù)減小，其中D1和D3的波速相對(duì)變化小于3‰，D2的波速相對(duì)變化最大，約為13.4‰，意味著震源區(qū)附近巖石遭受破壞，結(jié)構(gòu)弱化，應(yīng)力得到較大釋放；在D3的P 波尾波部分發(fā)現(xiàn)了走時(shí)延遲和不相關(guān)系數(shù)均存在1個(gè)短周期“尖峰”變化，可能與庫(kù)水在地下裂隙中的運(yùn)移有關(guān)。在L5550 臺(tái)站，D2的S 波尾波部分走時(shí)延遲明顯減小，波速上升約6.7‰；D3的S 波尾波部分走時(shí)延遲明顯增加，波速下降約6.1‰，二者得到的波速變化值接近，可能表明震后局部地區(qū)地震波速度存在上升—下降的恢復(fù)過(guò)程，該階段應(yīng)力狀態(tài)可能處于由累積到釋放的快速調(diào)整過(guò)程。在L5551 臺(tái)站，D1的S 波尾波部分走時(shí)延遲明顯增加，波速下降約2.1‰；D2的S 波尾波部分走時(shí)延遲明顯減小，波速上升約4.5‰。

選取布設(shè)在武隆MS5.0 地震震中附近3 個(gè)流動(dòng)臺(tái)的地震事件波形和觀測(cè)報(bào)告，利用波形互相關(guān)技術(shù)，結(jié)合雙差定位結(jié)果，選出3 組符合條件的重復(fù)地震。利用尾波干涉原理，通過(guò)移動(dòng)窗口互相關(guān)方法，計(jì)算了每組重復(fù)地震之間的走時(shí)延遲變化。尾波干涉結(jié)果顯示，在P 波的尾波部分，走時(shí)延遲變化較為復(fù)雜，可能與P 波和S 波不同的變化規(guī)律和PS 轉(zhuǎn)換波復(fù)雜的震相有關(guān)。在S 波振幅較大部分，相關(guān)系數(shù)較高，走時(shí)延遲變化不大，可能是由于S 波在重復(fù)地震識(shí)別過(guò)程中參與互相關(guān)計(jì)算的權(quán)重較高。在S 波早期的尾波部分，走時(shí)延遲變化最為明顯且較為線(xiàn)性，可能是由地殼介質(zhì)的地震波速變化引起的。尾波干涉結(jié)果反映的是整個(gè)傳播路徑上介質(zhì)變化的平均結(jié)果，距離震源區(qū)較遠(yuǎn)的XNS、LAX、XIT 等臺(tái)站，走時(shí)延遲變化并不明顯，說(shuō)明震源區(qū)外的地殼介質(zhì)變化微弱。

選擇流動(dòng)臺(tái)站的數(shù)據(jù)進(jìn)行尾波干涉分析，優(yōu)勢(shì)在于流動(dòng)臺(tái)站距震中位置足夠近，布局更加合理，且識(shí)別出的重復(fù)地震數(shù)量?jī)?yōu)于固定臺(tái)站，更能捕捉到地下介質(zhì)的細(xì)微變化。因武隆MS5.0 地震前近450 天研究區(qū)內(nèi)無(wú)地震記錄，不能通過(guò)尾波干涉方法得到震前地下介質(zhì)變化信息。

感謝江西省地震局查小惠工程師為本研究提供波形互相關(guān)和尾波干涉程序。