利用近震源區(qū)強(qiáng)震記錄識(shí)別魯?shù)?.5級(jí)地震后疊加的余震事件

崔慶谷 林國(guó)良 李 倩

1 云南省地震局，昆明市北辰大道148，650224

序列分析是震后地震趨勢(shì)判斷的重要手段，一個(gè)中強(qiáng)地震發(fā)生后，后續(xù)小震活動(dòng)是判斷地震類(lèi)型（主余型、震群型、前震－主震型）的重要依據(jù)。對(duì)于有前震的強(qiáng)震，前震序列對(duì)后續(xù)地震具有極其重要的指示意義。有學(xué)者認(rèn)為，在所有的短臨前兆中，前震是公認(rèn)的預(yù)報(bào)強(qiáng)震最有效的指標(biāo)之一［1］；對(duì)于主－余型或震群型地震，第1次強(qiáng)震后的小震活動(dòng)是判斷后續(xù)地震危險(xiǎn)性的重要依據(jù)［2］。

實(shí)際觀測(cè)中，緊隨主震后的余震往往被主震淹沒(méi)并導(dǎo)致地震目錄不完整，最終影響后續(xù)地震危險(xiǎn)性的判定。為從疊加的記錄數(shù)據(jù)中辨識(shí)出后續(xù)小震事件，前人做了大量工作。劉希強(qiáng)等利用小波變換將疊加起始時(shí)刻精確地估計(jì)出來(lái)，從而實(shí)現(xiàn)疊加小震事件的檢測(cè)［3］。Schaff等利用記錄波形與已知地震事件的互相關(guān)，將疊加的地震事件識(shí)別出來(lái)［4］。譚培毅等［5］利用匹配濾波技術(shù)檢測(cè)出地震目錄中被遺漏的小震事件。但上述方法的使用有先決條件，即疊加小震與之前地震的尾波在頻率、相位、振幅等某一方面存在顯著差異［6］。可見(jiàn)，基于數(shù)據(jù)處理的震相識(shí)別和疊加事件分離方法受諸多客觀條件的限制，如果數(shù)據(jù)本身不具備上述條件，純粹依靠數(shù)據(jù)處理并不能將疊加事件分離，此時(shí)必須尋找新的數(shù)據(jù)源。

微震觀測(cè)數(shù)據(jù)是地震編目中的傳統(tǒng)數(shù)據(jù)源，但是，對(duì)于中強(qiáng)以上地震，使用常規(guī)的數(shù)據(jù)處理方法難以有效分離疊加事件。而近震源區(qū)的強(qiáng)震記錄波列長(zhǎng)度短，余震疊加的幾率小且不易限幅，對(duì)識(shí)別疊加的強(qiáng)余震具有特別的優(yōu)勢(shì)，可以作為新的數(shù)據(jù)源用于識(shí)別疊加地震。

本文通過(guò)對(duì)比不同震中距記錄的魯?shù)?.5級(jí)地震，分析了震中距對(duì)地震波波列長(zhǎng)度的影響，并利用震源區(qū)強(qiáng)震記錄識(shí)別提取出兩次被淹沒(méi)的余震事件，其中一次事件在地震目錄中沒(méi)有記載。

1 魯?shù)?.5級(jí)地震波波列隨震中距的變化

在散射很弱的情況下，近源區(qū)臺(tái)站記錄到的地震波表現(xiàn)為突然急始后快速消失的孤立子波［7］。1993年印度Killari－Latur地震的余震記錄（震中距為5.3km）表明，震源區(qū)的地震記錄只有持續(xù)時(shí)間很短、清晰急始的P 波和S 波。但是，在地震臺(tái)站密度較低的情況下，要記錄到上述短時(shí)子波很難，因?yàn)槭孪群茈y把地震臺(tái)站準(zhǔn)確架設(shè)在震源區(qū)。

隨著震中距的增大，由于傳輸路徑上介質(zhì)不均勻性增加，P波和S波開(kāi)始出現(xiàn)強(qiáng)散射，表現(xiàn)為初至子波之后持續(xù)一段時(shí)間的混合信號(hào)，這一混合信號(hào)由噪聲和尾波組成。震中距越大，源于介質(zhì)不均勻性的地震波波列也越長(zhǎng)，地震事件疊加的可能性也就越大。當(dāng)震中距增加到一定距離后，面波的頻散形成更長(zhǎng)的波列，此時(shí)地震波波列長(zhǎng)度急劇增加，短時(shí)間內(nèi)的幾次地震事件可能因疊加而無(wú)法分離。因此，利用震源區(qū)的強(qiáng)震記錄來(lái)識(shí)別余震，理論上是一種可行的方法。因?yàn)檎鹬芯嘣叫。瑐鞑ヂ窂缴系慕橘|(zhì)越單一，記錄波列也就越短，后續(xù)事件出現(xiàn)重疊的概率也就越小。這一觀點(diǎn)在2014云南魯?shù)?.5級(jí)地震的記錄數(shù)據(jù)中得到了很好的印證。

魯?shù)?.5 級(jí)地震的震中位于魯?shù)榭h龍頭山鎮(zhèn)，地震發(fā)生后，分布于震中龍頭山鎮(zhèn)及附近區(qū)域的臺(tái)站完整記錄到了該次地震。這些臺(tái)站的記錄情況見(jiàn)表1。

表1 魯?shù)镸6.5級(jí)地震近震源區(qū)（震中距Δ＜80km）臺(tái)站一覽表Tab.1 Stations which near the epicentral area of Ludian M6.5（Δ＜80km）

魯?shù)辇堫^山臺(tái)布設(shè)在主震的震中區(qū)，主震記錄數(shù)據(jù)中的S－P 約為1s，對(duì)應(yīng)的震中距離約為8.33km，比中國(guó)地震臺(tái)網(wǎng)中心基于微震臺(tái)網(wǎng)觀測(cè)數(shù)據(jù)測(cè)定的震源深度（12km）略小。從震中的龍頭山臺(tái)到鉛廠、茨院，震源距逐漸增加，記錄波形也逐漸變長(zhǎng)（圖1）。

圖1 龍頭山、鉛廠、茨院強(qiáng)震記錄對(duì)比（橫坐標(biāo)起始零時(shí)刻為16：29：52）Fig.1 Comparing strong acceleration recording（zero point of time in abscissa is 2014－11－29 16：29：52）

圖1中，龍頭山、鉛廠、茨院的震源距分別為8.33、18.61、32.60km，隨著震源距的增加，地震波波列逐漸變長(zhǎng)。由于龍頭山觀測(cè)點(diǎn)離震源最近，其記錄數(shù)據(jù)中隱約可以識(shí)別出3次較小的地震事件（圖1（a）），時(shí)間坐標(biāo)分別為29～31s、47s、66.5s。在鉛廠的記錄中，時(shí)間坐標(biāo)33s處也隱約可見(jiàn)一次疊加小震的痕跡（圖1（b）），但其波形疊加到主震尾波中難以識(shí)別。當(dāng)震中距增加到32.60km 時(shí)，小地震事件已經(jīng)完全淹沒(méi)在主震數(shù)據(jù)中（圖1（c））。

在震源區(qū)附近80km 以內(nèi)的臺(tái)站中，巧家、昭通兩個(gè)臺(tái)站同時(shí)架設(shè)有強(qiáng)震和微震兩種儀器。進(jìn)一步對(duì)比微震觀測(cè)數(shù)據(jù)發(fā)現(xiàn)，震源距47.35km的昭通地震臺(tái)微震記錄已經(jīng)嚴(yán)重限幅，地震波波列也比龍頭山臺(tái)的強(qiáng)震記錄長(zhǎng)，以至于所有的余震事件完全淹沒(méi)在主震尾波中難以識(shí)別。兩個(gè)臺(tái)站NS分向記錄數(shù)據(jù)見(jiàn)圖2。

圖2 龍頭山強(qiáng)震、昭通臺(tái)微震記錄對(duì)比（橫坐標(biāo)起始零時(shí)刻為16：29：52）Fig.2 Comparingstrong acceleration recording at Longtoushan site（zero point of time in abscissa is 2014－11－29 16：29：52）

可見(jiàn)，近震源區(qū)的記錄受傳輸路徑的干擾較小，包含了更多的震源信息，因此可以用龍頭山臺(tái)的記錄數(shù)據(jù)提取疊加的余震事件。

2 利用震源區(qū)強(qiáng)震記錄檢測(cè)疊加的地震事件

圖2（a）中，在時(shí)間坐標(biāo)29、47、66.5、103s左右可能存在余震事件，但經(jīng)核實(shí)，具有完整P 波和S 波特征的只有29s、103s處的兩個(gè)事件：29s處為第一個(gè)余震事件，103s處為第二個(gè)余震事件。近震源區(qū)龍頭山臺(tái)的強(qiáng)震儀記錄的魯?shù)?.5級(jí)主震波形如圖3。

對(duì)比圖3中垂直分向與水平分向的記錄可以看出，垂直分向的高頻成分遠(yuǎn)遠(yuǎn)大于兩個(gè)水平分向，垂直向的P 波比水平向的P 波更發(fā)育，而S波則正好相反。如果用垂直向的振幅譜與水平向的振幅譜進(jìn)行對(duì)比，每一次新的地震事件，譜比值都會(huì)經(jīng)歷一次從大到小的變化，由此可以識(shí)別出疊加地震事件的到時(shí)。第1次余震事件與主震尾波疊加的細(xì)節(jié)見(jiàn)圖4。

圖3 龍頭山強(qiáng)震記錄波形圖（橫坐標(biāo)起始零時(shí)刻為16：29：52）Fig.3 Strong acceleration recording at Longtoushan site（zero point of time in abscissa is 2014－11－29 16：29：52）

圖4 與主震尾波疊加的第1次余震事件Fig.4 The first aftershock overlap with mainshock’s coda wave in strong acceleration recording at Longtoushan site

圖4中，自時(shí)間坐標(biāo)29.3s起（虛線所示），UD 向記錄數(shù)據(jù)中的高頻成分與兩個(gè)水平分向相比突然增強(qiáng)，代表一個(gè)新的地震事件的到來(lái)。用等長(zhǎng)度的滑動(dòng)窗沿時(shí)間軸滑動(dòng)，并在窗內(nèi)求UD分向與兩個(gè)水平分向功率譜的比值，可以得到3維坐標(biāo)（時(shí)間－頻率－譜比）中譜比值變化的圖像。由于余震事件波形的頻率區(qū)域?yàn)?～30 Hz，從3維坐標(biāo)中截取5～30Hz頻段的譜比值正視圖，得到疊加前后譜比值隨時(shí)間變化的圖像（圖5）。

圖5中的譜比值自時(shí)間坐標(biāo)29.3s起開(kāi)始增大，標(biāo)志著第1次余震事件的到來(lái)；30.5s后P波結(jié)束，譜比值恢復(fù)正常。這一結(jié)果與圖4中的疊加情況相符合。

自29.3s開(kāi)始截取第1 次余震事件的強(qiáng)震記錄，經(jīng)數(shù)值積分后可轉(zhuǎn)換為速度記錄［8］，計(jì)算結(jié)果見(jiàn)圖6。

圖5 垂直分向與水平分向功率譜比值隨時(shí)間變化圖Fig.5 The spectrum ratio of UD／EW and UD／NS in strong acceleration recording at Longtoushan site

圖6 第1次余震事件強(qiáng)震記錄數(shù)據(jù)計(jì)算得到的速度記錄波形Fig.6 The first aftershock’s velocity recording get from strong acceleration recording at Longtoushan site

與此相似，由強(qiáng)震記錄計(jì)算第2次事件的速度記錄波形，結(jié)果見(jiàn)圖7。

圖7 由第2次余震事件強(qiáng)震記錄數(shù)據(jù)計(jì)算得到的速度記錄波形Fig.7 The second aftershock’s velocity recording get from strong acceleration recording at Longtoushan site

利用圖6、圖7 中的波形數(shù)據(jù)可以計(jì)算出兩次余震的里式震級(jí)。里式震級(jí)的計(jì)算公式為［9］：

式中，Vmax為最大速度振幅，單位為μm，實(shí)際計(jì)算中取兩個(gè)水平分向最大振幅的平均值；T為最大振幅對(duì)應(yīng)的周期，單位為s；R（Δ）為量歸函數(shù)；Δ為震中距；s為臺(tái)站校正值。

s為未知，可以利用已知的主震震級(jí)求解，求解過(guò)程如下。

將M6.5主震的強(qiáng)震記錄數(shù)據(jù)積分，得到其速度記錄；將速度記錄中的最大振幅值代入式（1），計(jì)算得到主震震級(jí)為：ML＝lg（0.8×55×104／2π）＋1.8＋s＝6.65＋s；根據(jù)中國(guó)地震臺(tái)網(wǎng)中心的地震目錄，魯?shù)榈卣鹱詈蟠_定的震級(jí)為M6.6，根據(jù)震級(jí)轉(zhuǎn)換公式M＝1.13×ML－1.08＝6.434 8＋1.13s＝6.6，得到s≈0.1。

將圖6中的速度記錄峰值數(shù)據(jù)代入式（1），得到第1次強(qiáng)余震的震級(jí)為：ML＝lg（0.25×3.15×104／2π）＋1.8＋0.1＝5.0，轉(zhuǎn)換為M 震級(jí)：M＝1.13×5.0－1.08＝4.6。

第2次余震的震級(jí)為：ML＝lg（0.2×5×103／2π）＋1.8＋s＝4.2，轉(zhuǎn)換為M 震級(jí)：M＝1.13×4.2－1.08＝3.6。

進(jìn)一步利用龍頭山臺(tái)強(qiáng)震數(shù)據(jù)計(jì)算較強(qiáng)余震的震級(jí)并將計(jì)算結(jié)果與地震目錄對(duì)比，結(jié)果見(jiàn)表2。

表2 地震目錄中的魯?shù)?.5級(jí)系列震級(jí)（M≥3.5）與強(qiáng)震記錄計(jì)算結(jié)果對(duì)比Tab.1 Comparion of magnitude get from strong accelerator records with existing catalog earthquake swarm of Ludian M6.5（M≥3.5）

從表2中可以看出，用強(qiáng)震數(shù)據(jù)計(jì)算得到的震級(jí)比地震目錄中的震級(jí)略大，但誤差在可以接受的范圍，說(shuō)明計(jì)算結(jié)果是可靠的。

3 結(jié) 語(yǔ)

近震源區(qū)的強(qiáng)震記錄對(duì)于疊加強(qiáng)余震的識(shí)別具有獨(dú)特的作用。

2014年魯?shù)?.5級(jí)地震的實(shí)例表明，近震源區(qū)強(qiáng)震記錄的波列長(zhǎng)度短，余震疊加的幾率小，從該記錄數(shù)據(jù)中有效識(shí)別并提取常規(guī)記錄中難以識(shí)別的兩次強(qiáng)余震事件，震級(jí)分別為M4.6、M3.6。

此外，在近震源區(qū)的強(qiáng)震記錄中，龍頭山臺(tái)記錄的6.5級(jí)主震從P波初至到S波最大振幅處只需約3s且不限幅。即使是在震源距18km 以外的鉛廠，從P 波初至到S 波最大振幅處也只需8s。可見(jiàn)，近震源區(qū)的強(qiáng)震加速度記錄在地震預(yù)警中的優(yōu)勢(shì)是明顯的，這與當(dāng)前國(guó)外地震預(yù)警的實(shí)踐經(jīng)驗(yàn)相吻合［10］。

［1］李麗，陳颙.2009年4月6日意大利拉奎拉地震的前震及其預(yù)測(cè)意義［J］中國(guó)地震，2009，25（2）：151－158（Li Li，Chen Yong.Light of Earthquake Prediction：Shocks before the L’Aquila Earthquake of April 6，2009［J］.2009，25（2）：151－158）

［2］馮建剛，蔣長(zhǎng)勝，韓立波，等 甘肅測(cè)震臺(tái)網(wǎng)監(jiān)測(cè)能力及地震目錄完整性分析［J］.地震學(xué)報(bào)，2012，34（5）646－658（Feng Jiangang，Jiang Changsheng，Han Libo，et al.Analysis on the Monitoring Capability of Seismic Networks and Completeness of Earthquake Catalogues in Gansu Region［J］.Acta Seismologica Sinica，2012，34（5）：646－658）

［3］劉希強(qiáng)，周蕙蘭，鄭治真，等.基于小波包變換的弱震相識(shí)別方 法［J］.地震學(xué)報(bào)，1998，20（4）：373－380（Liu Xiqiang，Zhou Huilan，Zheng Zhizhen，et al.Identification Method of Weak Seismic Phases on the Basis of Wavelet Packet Transform［J］.Acta Seismologica Sinica，1998，11（4）：431－439）

［4］Schaff D P，Waldhauser F.Waveform Cross－crorrelationbassed Differential Travel－Time Measurements at the Northern California Seismic Network［J］.Bull Seism Soc Am，2005，95（6）：2 446－2 661

［5］譚毅培，曹井泉，劉文兵，等，2013年3月涿鹿微震群遺漏地震事件檢測(cè)和發(fā)震構(gòu)造分析［J］.地球物理學(xué)報(bào)，2014，57（6）：1 847－1 856（Tan Yipei，Cao Jingquan，Liu Wenbing，et al.Missing Earthquakes Detection and Seismogenic Structure Analysis of the Zhuolu Micro－earthsquake Swarm in March 2013［J］.Chinese J Geophys，2014，57（6）：1 847－1 856）

［6］金龍，陳小宏，李景葉.基于誤差準(zhǔn)則和循環(huán)迭代的時(shí)移地震匹配濾波方法［J］.地球物理學(xué)報(bào)，2005，48（3）：698－703（Jin Long，Chen Xiaohong，Li Jingye.A New Method for Time－Lapse Seismic Matching Filter Based on Error Criteria and Cyclic Iteration［J］.Chinese J Geophys，2005，48（3）：698－703）

［7］Borman P.新地震觀測(cè)實(shí)踐手冊(cè)［M］.北京：地震出版社，2006（Borman P.New Manual of Seismological Observatory Practice［M］.Beijing：Seismological Press，2006）

［8］Yorihiko O.地震動(dòng)的譜分析入門(mén)［M］.北京：地震出版社，2008（Yorihiko O.Shin Jishindo No Spectre Kaiseki Nyumon［M］.Beijing：Seismological Press，2008）

［9］陳運(yùn)泰，劉瑞豐.地震的震級(jí)［J］.地震地磁觀測(cè)與研究，2004，25（6）：1－12（Chen Yuntai，Liu Ruifeng.Earthquake Magnitude［J］.Chinese J Seismological and Geomagnetic Observation and Research，2004，25（6）：1－12）

［10］Gasparini P，Manfredi G，Zschau J.地震預(yù)警系統(tǒng)［M］.北京：地震出版社，2014（Gasparini P，Manfredi G，Zschau J.Seismic Early Warning System［J］.Beijing：Seismological Press，2014）