云南魯甸6.5級地震序列的譜振幅相關系數變化特征*

黃 浩，付 虹，王培玲，馬玉虎

(1.青海省地震局，青海西寧810001;2.云南省地震局，云南昆明650224)

0 引言

據中國地震臺網中心測定，北京時間2014年8月3日16時3分，在云南省昭通市魯甸縣龍頭山鎮發生 MS6.5地震，震中位置 (27.1°N，103.3°E)，震源深度12 km。截至9月9日，震區共記錄到68次ML≥3.0余震，其中ML3.0～3.9地震63次，ML4.0～4.9地震5次，沒有發生ML≥5余震。

魯甸6.5級地震極震區烈度達到Ⅸ度，地震序列為主震—余震型，主震之后沒有發生更大地震，最大余震為MS4.4。余震震級偏低、中強余震次數少，是否意味著震源區仍存在發生強余震的可能?在余震活動相對減弱的情況下，利用震源機制一致性來研究這一問題也許可以有些幫助。

已有震例研究表明 (陳颙，1978;Wiemer，Wyss，2002;趙英萍等，2004;王俊國，刁桂苓，2005)，震源機制一致性可用來判斷未來是否會發生更大地震。然而，由于受到臺站方位分布和密度的限制，很難精確求解中小地震的震源機制解，因此中小地震震源機制一致性方面的研究還比較少。

Lund和B?evarsson(2002)提出微震體波譜振幅相關分析法，該方法可以不直接求解震源機制，而利用譜振幅的相關程度來描述震源機制的一致性。崔子健等 (2012)利用譜振幅相關分析法研究了發生在云南的幾個小震群序列，結果表明強震前的小震群序列的譜振幅相關系數較大，而一般性小震群序列的譜振幅相關系數較小，獲得了與理論推測相符合的結果。黃浩和付虹(2014)應用譜振幅相關分析法研究了6個發生在滇西地區不同類型的地震序列，結果表明前震序列的譜振幅相關系數較高，但譜振幅相關系數較高的不一定是前震序列。雖然譜振幅相關分析法在云南地區初步應用的結果有所不同，但還應積累更多的震例進行深入研究。本文采用譜振幅相關分析法，計算了魯甸地震序列的譜振幅相關系數，并對其變化特征進行了分析。

1 譜振幅相關分析法

譜振幅相關分析法要求參與計算的地震擁有近似的射線路徑以扣除傳播路徑效應，因此適用于時間和空間上叢集發生的小震群。崔子健等(2012)已對譜振幅相關分析法進行過詳細的描述，本文只做簡要介紹。臺站記錄的位移譜Uij(f)可表示為

式中，Uij(f)為臺站j記錄的地震的位移譜;f為頻率;Pij(f)為臺站j至地震i之間的傳播路徑效應;Si(f)為地震i包含輻射圖型因子的震源譜;L'j(f)為臺站j的局部場地效應;Nj(f)為臺站j的噪聲;Ij(f)為臺站j的儀器響應;Wj(f)為臺站j的自由表面效應。

在去除噪聲、儀器響應、臺站場地效應和傳播路徑效應后，若震源機制相同，則兩次地震的位移譜的相對大小只與譜振幅的相對大小有關。因此，可計算兩次地震的譜振幅相關系數來描述震源機制的相似程度。

譜 振 幅 相 關 系 數 rxy(Lund，B?evarsson，2002)可表示為

式中，xjl、yjl分別為地震x和y的分量的譜振幅的對數值;分別為xjl與yjl的平均值;l取值1～5，分別表示P波的徑向、垂向分量與S波的徑向、切向、垂向分量;n為參與計算的臺站數量。

計算每次地震和地震序列中發生在它前面的所有地震的相關系數，取平均后，代表這次地震與前面所有地震譜振幅的相關程度，這樣就可得到整個序列的譜振幅相關系數。

計算過程按以下步驟進行:①收集魯甸地震序列的數字波形;② 把波形的水平分量旋轉為徑向和切向分量，并采用延遲時間窗方法 (Chael，1987;黃玉龍等，2003)計算其位移譜;③利用Atkinson和Mereu(1992)方法反演魯甸震區的Q值，QP=35.5f0.831、QS=53.2f0.872;④ 采用多臺多事件聯合反演臺站的場地響應 (Moya，Jorge，2000);⑤利用經過噪聲、儀器響應、傳播路徑及場地響應校正的譜振幅，計算魯甸地震序列的譜振幅相關系數。圖1為本文使用的昭通臺記錄的一次地震波形及其位移譜，圖中，Z、R、T分別表示垂向、徑向、切向分量;PZ、PR分別表示P波垂向和徑向分量的位移譜;SZ、SR、ST分別表示S波垂向、徑向和切向分量的位移譜。

2 地震序列資料及計算結果

2.1 地震序列資料

魯甸地震發生后，截至9月9日，震源區共發生ML≥3.0地震64次。剔除部分信噪比較差或受主震影響較大的余震，本文選取了48次ML≥3.0地震計算了魯甸序列的譜振幅相關系數 (圖2)。

圖1 昭通臺記錄的一次地震的波形 (a)及其位移譜 (b)Fig.1 Waveform of a earthquake recorded by Zhaotong Station(a)，and its displacement spectra(b)

2.2 結果分析

采用昭通臺 (ZAT)、巧家臺 (QIJ)、普格臺(PGE)、雷波臺 (LBO)、東川臺 (DOC)5個地震臺記錄的48個ML≥3.0的波形資料計算了魯甸地震序列的譜振幅相關系數，其相關系數值在0.77～0.83之間變化，平均值0.82。

雖然譜振幅相關系數的變化范圍不大，但具有明顯的起伏過程 (圖3)。從相關系數隨時間變化的特點來看，大致可以以8月11日為界將譜振幅相關系數的變化過程劃分為兩個階段，第一階段相關系數波動較大 (圖3a)，第二階段相關系數保持平穩 (圖3b)。雖然兩個階段的相關系數均值同為0.82，但第一階段相關系數變化更明顯，且多數MS≥4.0余震發生在第一階段，因此本文重點分析第一階段的譜振幅相關系數變化特征。

第一階段譜振幅相關系數隨時間變化的形態比較豐富 (圖3a)，可以劃分為4個變化過程，即:①主震后3小時內，相關系數自0.81下降至0.77;②8月3日19時至8月4日5時，相關系數自整個序列的最小值0.77上升至最大值0.83;③8月4日5時至8月5日8時，相關系數總體呈下降過程;④8月5日8時至8月11日，相關系數緩慢上升至平均值并保持穩定。

譜振幅相關系數在主震后3小時內快速下降至最小值0.77，表明主震后震源機制一致性在短時間內有所降低。在第二個變化過程中，譜振幅相關系數持續升高，從序列最小值0.77上升為序列最大值0.83。在此期間，震源區相繼發生了4次4級余震，表明這段時間震源機制表現出趨于一致的現象。在第三個變化過程中，譜振幅相關系數既有下降也有上升，但總體表現為下降，說明這一時段的震源機制一致性有所降低。在第四個變化過程中，譜振幅相關系數在緩慢上升后表現平穩，表明震源機制一致性程度趨于穩定。

圖2 魯甸M L≥3.0地震序列和周圍臺站分布Fig.2 Distribution of Ludian M L≥3.0 earthquake sequence and their surrounding stations

圖3 魯甸地震序列的相關系數及震級—時間圖Fig.3 The correlation coefficient of spectral amplitude and M-t plot of Ludian earthquake sequence

3 結論與討論

魯甸地震序列的譜振幅相關系數可以大致以8月11日為界劃分為兩個階段，第一階段相關系數變化多樣，第二階段相關系數保持穩定。第一階段可根據變化特征不同，劃分為4個變化過程，其中4次4級余震發生在第二個變化過程。在第二個變化過程中，相關系數持續升高、上升幅度最大，說明這段時間震源機制表現出趨于一致的現象。

魯甸地震序列的譜振幅相關系數較高，但至今沒有發生更大地震。崔子健等 (2012)的研究結果表明，強震前的小震群序列和一般性小震群序列的譜振幅相關系數差異較大;黃浩和付虹(2014)對滇西地區的6個地震序列進行譜振幅相關分析后認為，前震序列的譜振幅相關系數較高，但譜振幅相關系數較高的序列不一定是前震序列。因此據以往震例，魯甸地震序列的譜振幅相關系數較高，只能說明區域應力水平較高，不能判定該序列不會發生較大地震，也不能說明其后還有較大地震發生。目前，譜振幅相關分析法的應用還不多，需要積累更多的震例進行驗證。

當譜振幅相關系數較高時，區域應力水平較高，震源區及其附近地區繼續發生MS≥5.0地震的可能性相對較大。8月11日以后，相關系數仍穩定在0.82左右，雖然在震源區沒有發生更大地震，但8月17日在距主震震中110 km的永善縣發生了MS5.0地震。這與2008年3月21日盈江MS5.2和2011年3月10日盈江MS5.8地震序列譜振幅相關系數較高，其后在附近地區發生MS≥5.0地震的現象類似。

云南省數字地震臺網為本文提供了相關資料，云南省地震預報中心劉自鳳和孫楠在收集資料過程中提供了許多幫助，審者提出了寶貴意見，在此表示衷心的感謝!

陳颙.1978.用震源機制一致性作為描述地震活動性的新參數［J］.地球物理學報，21(2):142-159.

崔子健，李志雄，陳章立，等.2012.判別小震群序列類型的新方法研究——譜振幅相關分析法［J］.地球物理學報，55(5):1718-1724.

黃浩，付虹.2014.2008年以來滇西地區地震序列的譜振幅相關系數變化特征［J］.地震學報，36(4):631-639.

王俊國，刁桂苓.2005.千島島弧大震前哈佛大學矩心矩張量(CMT)解一致性的預測意義［J］.地震學報，27(2):178-183.

趙英萍，刁桂苓，高景春，等.2004.張北強余震前震源機制解的一致性特征［J］.華北地震科學，22(1):1-4.

Atkinson G.M.，Mereu R.F..1992.The shape of ground motion attenuation curves in southeastern Canada［J］.Bull.Seismol.Soc.Am.，82(5):2014-2031.

Chael E..1987.Spectral scaling of earthquakes in the Miramichi region of New Brunswick［J］.Bull Seismol Soc Am，77(2):347-365.

Lund B.，B?evarsson R..2002.Correlation of microearthquake bodywave spectral amplitudes［J］.Bull.Seismol.Soc.Am.，92(6):2419-2433.

Moya A.，Jorge A..2000.Inversion of source parameters and site effects from strong ground motion records using genetic algorithms［J］.Bull.Seismol.Soc.Am.，90(4):977-992.

Wiemer S.，Wyss M..2002.Mapping spatial variability of the frequencymagnitude distribution of earthquakes［J］.Advances in Geophysics，45:259-302.