2013年吉林前郭強震群序列的譜振幅相關系數變化特征

康建紅，鄭國棟，張 宇，張洪艷

（1.吉林省地震局，吉林 長春 130117；2.長春地震臺，吉林 長春 130117）

2013年吉林前郭強震群序列的譜振幅相關系數變化特征

康建紅1，鄭國棟2，張 宇2，張洪艷1

（1.吉林省地震局，吉林 長春 130117；2.長春地震臺，吉林 長春 130117）

前郭強震群；譜振幅相關系數；震源機制一致性；區域應力場

0 引言

2013年10月31日11時3分在吉林省松原市前郭縣發生MS5.5地震，震源深度8km，7分鐘之后，同一震中區再次發生MS5.0地震。之后地震活動進入衰減階段，余震序列地震數量偏少，較大震級余震成對發生；11月22日16時18分震中區再次發生MS5.3地震，次日在6時4分、6時32分連續發生MS5.8、MS5.0地震。地震序列十分復雜。如何根據當前序列發展，判斷未來震情發展趨勢，一直是科研人員致力解決的難題。

對地震序列的研究方法，尤其是基于地震目錄的分析方法很多，如吳開統的地震序列概論[1]、蔣海昆的中國大陸地震序列研究[2]等。近年來，隨著數字地震學的飛速發展，利用波形資料提取震源信息，并將其應用于地震序列的判斷，逐漸成為地震序列原有研究方法的重要補充。

震源機制解描述震源的幾何特性，根據震源機制解可以判斷發震斷層的錯動方式，推斷應力場的方向。陳颙[3]通過研究邢臺、海城等強震，提出使用震源機制一致性作為判別是否為前震的新指標；刁桂苓等[4-6]通過研究國內外多個震例，得出主震前中小地震震源機制解趨于一致的結論；王俊國等[7]在研究千島弧地區震源機制一致性特征時，提出了震源機制和構造應力場的一致性參數概念，并利用這一參數進行強震預測。

實踐中，由于受臺網布局及密度的影響，能記錄到小地震的臺站數量有限，因此，對于小地震，求解震源機制存在較大的不確定性，難以精確求解。為此，Lund等[8]曾提出小震體波譜振幅相關分析方法，很多學者利用這一方法計算地震序列中小地震體波譜振幅相關系數，對震源機制一致性進行歸一化量化，直接反映地震序列的震源機制相對變化過程。朱航等[9]利用該方法對2001—2003年川滇地區若干次地震序列的震源機制變化過程進行研究，結果表明，強震序列與小震序列的譜振幅相關系數曲線形態存在顯著差異，強震發生之前，震源機制有從散亂到一致的變化過程。崔子健等[10]利用該方法對云南西北部的五個震群的譜振幅相關系數進行計算，結果顯示，強震前的小震序列譜振幅相關系數接近1.0，比一般性小震序列明顯偏高，說明其震源機制一致性較高。

本文采用譜振幅相關分析法，計算了2013年前郭強震群序列中小地震體波譜振幅相關系數，總結其變化特征。

1 前郭強震群概況

截至2016年12月31日，前郭震群共記錄地震1356次，其中：MS5.0以上地震5次，MS4.0～4.9地震10次，MS2.0～3.9地震相對較少，大多數地震震級集中在MS0.0～1.0，4.0級以上地震成對發生（圖1、表1），最大地震為2013年11月23日MS5.8地震，是自1960年榆樹土橋5.8級地震以來吉林省遭遇的最強地震活動。依據中國地震局《地震現場工作大綱和技術指南》[11]，前郭地震序列為震群型，10月31日MS5.5地震為首發地震[12]。在G-R關系中，其序列b值為0.4，顯著低于東北地區平均水平（0.6～0.8）[13]。

圖1 前郭強震群序列M-T、頻度圖Fig.1 M-T and frequency diagrams of Qianguo strong earthquake swarm sequence

表1 前郭震群地震序列目錄（MS≥4.0）

2 資料選取

自2013年10月31日首震發生以來至2016年12月31日，前郭強震群共記錄地震1356次，其中，ML≥2.5地震66次。選取松原臺（SYT）、長春臺（CN2）、乾安臺（QAT）、白城臺（BCT）、長嶺臺（CLT）5個對震群包絡較好且觀測質量較高的臺站，選取的地震事件要求5個臺站同時記錄清晰且信噪比高。本研究側重探討中小地震震源機制一致性與中強地震的關系，5個ML＞5.0地震不計算在內。綜合以上原則，選擇序列中45次地震（2.5≤ML≤5.0）作為研究對象，分析前郭強震群的譜振幅相關系數變化特征。圖2為吉林及鄰省數字化地震臺網、主要斷裂及ML≥2.5地震震中分布圖。

圖2 前郭強震群及周圍臺站分布圖Fig.2 Distribution of Qianguo strong swarm and the surrounding stations

3 方法原理

在地震學中，臺站記錄126oE 地震波位移觀測譜振幅的表達式為：

（1）式中，Uij(f)為第j個臺站記錄的第i個地震的位移譜；f為頻率；r為震源距；Mi(f)為地震震源譜；φ(θ,?)為地震波輻射花樣因子， 、θ為地震到臺站的離源角、方位角；Pij(r)為幾何擴散系數；Fij(f,r)為非彈性衰減系數；Gj(f)為臺站場地響應系數，Ij(f)為儀器響應。

由式（1）可知：當多個地震位于同一個位置，或者地震之間的距離遠小于地震到臺站的距離時，由于地震波的傳播路徑以及臺站場地完全相同，因此，這些地震到同一個臺站的幾何擴散系數、非彈性衰減系數、儀器響應系數和場地響應系數是完全相同的，同一臺站記錄到的多個地震之間的零頻觀測譜的差異主要來源于地震之間的震源機制差異。由以上分析可知：通過對同一震源區、不同地震事件在同一個臺站得到的零頻觀測譜值進行對比，可以得到不同地震之間的震源機制相似程度。

在計算過程中，對同一震源區的兩次地震事件在同一臺站記錄的波形數據，使用直達P、S波的零頻譜振幅值來計算相關系數。每個臺站的3分向波形記錄可以轉換成5個零頻譜值，即垂向P波、徑向P波、垂向S波、徑向S波和切向S波5個譜值分量。計算地震x、y的譜振幅相關系數：

式（2）中，xij、yij為地震x、y在第i個臺站上的第j個譜值分量的對數。

譜振幅相關系數越接近于1，則說明臺站記錄的兩次地震波形越相似，其震源機制解的一致性越高；反之，相關系數較低，則兩次地震的波形越不相似，其震源機制解也不相同。由此可以認為，體波譜振幅相關系數是描述震源機制一致性程度的數值。

4 計算結果

利用前郭強震群周邊松原臺（SYT）、長春臺（CN2）、乾安臺（QAT）、白城臺（BCT）、長嶺臺（CLT）5個地震臺記錄到的45次2.5≤ML≤5.0的余震觀測資料，計算前郭強震群序列的譜振幅相關系數，圖3上圖為譜振幅相關系數隨時間變化曲線，由圖看出，其數值在0.85～0.95之間，平均值為0.90，該序列具有較高的震源機制一致性。

從該序列的譜振幅相關系數變化曲線圖上看，相關系數變化大致分為5個階段：（1）2013年10月31日首震發生后至12月13日，譜振幅相關系數最高達到0.95，地震活動密集，期間共發生11次MS≥4.0地震；（2）2013年12月13日至2014年8月5日，譜振幅相關系數低于0.92，與上一階段對比，相對較低，本階段發生1次MS≥4.0地震；（3）2014年8月5日至2015年2月17日，譜振幅相關系數為0.92～0.95，處于相對高值水平，此階段發生了2次MS≥4.0地 震；（4）2015年2月17日至同年10月21日，譜振幅相關系數低于0.90，此階段未發生MS≥4.0地震；（5）2015年10月21日以后，譜振幅相關系數逐漸升高，總體高于0.90，至2016年12月31日，發生1次MS≥4.0地震，震級為MS4.5。目前震群序列譜振幅相關系數值較高，接近0.94，說明震區受區域應力場控制仍較強。

圖3 2013年前郭強震群序列體波譜振幅相關系數隨時間變化曲線（上）及震級-時間圖（下）Fig.3 Curves of the correlation coefficient of spectral amplitude varied with time and M-T plot of Qianguo strong earthquake swarm

強余震發生之前，譜振幅相關系數達到較高水平（圖4）。2013年11月22日MS5.3地震發生之前，譜振幅相關系數出現上升，達到較高水平，最高為0.94。

在較大余震發生后小震機制有一定的發散[9]，即出現震源機制不相似現象，滑動平均譜振幅相關系數時序曲線下降。如圖4所示，2013年11月22日MS5.3地震之后，譜振幅相關系數降低，11月23日MS5.8地震后，譜振幅相關系數仍較低，但很快譜振幅相關系數又恢復至較高值。說明區域應力場強度較高，2次5級以上地震的發生并未改變區域應力場強度。

圖4 前郭強震群序列體波譜振幅相關系數變化曲線Fig.4 The correlation coefficient curve of spectral amplitude of Qianguo strong earthquake swarm

5 討論與結論

通過對比同一臺站記錄的來自同一震源區的不同地震的直達波零頻觀測譜值，可以得到不同地震事件之間的震源機制一致性程度，由此得到整個序列的震源機制一致性變化過程。本文經計算得到2013年前郭強震群序列的譜振幅相關系數變化曲線，總體上，曲線高值時段較多，說明該震群序列的震源機制一致性較好，受區域應力場控制較強，這可能也是該震群序列一直持續發生的原因。

強余震發生在譜振幅相關系數較高的背景下，即震源機制一致性較高的背景下，與刁桂苓[14]、趙英萍[15]研究結果一致。11月22日MS5.3地震、11月23日MS5.8、5.0余震發生前，譜振幅相關系數曾達到較高值。

強余震發生后，出現譜振幅相關系數下降現象，即震源機制發散現象，與朱航[9]研究結果一致。5次MS≥5.0地震發生后，均出現了譜振幅相關系數明顯下降的現象。

[1] 吳開統， 焦遠碧， 呂培苓， 等. 地震序列概論[M]. 北京:北京大學出版社，1990，97-118.

[2] 蔣海昆， 傅征祥， 劉杰. 中國大陸地震序列研究[M]. 北京 : 地震出版社， 2007， 38-51.

[3] 陳颙. 用震源機制一致性作為描述地震活動性的新參數[J]. 地震物理學報， 1978， 21（2）: 142-159.

[4] 刁桂苓， 于新昌. 唐山地震前后京、津、唐、張地區的綜合斷層面解[J]. 西北地震學報， 1980， 2（3）: 39-47.[5] 刁桂苓， 于利民， 李欽祖. 強震前后震源區應力場變化一例 [J]. 地震學報， 1994， 16（1）: 64-69.

[6] 刁桂苓， 徐錫偉， 陳于高， 等. 汶川MW7.9和集集MW7.6地震前應力場轉換現象及其可能的前兆意義[J].地球物理學報， 2011， 54（1）: 128-136.

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

[8] Lund B， Boevarsson R. Correlation of microearthquake body-wave spectral amplitudes[J]. Bull. Seism. Soc. Am.，2002， 92（6）: 2419-2433.

[9] 朱航， 劉杰， 陳天長. 采用體波譜振幅相關系數方法研究地震序列的震源機制變化過程[J]. 地震， 2006， 26（2）: 1-11.

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

[11] 中國地震局. 地震現場工作大綱和技術指南[M]. 北京:地震出版社， 1998， 15-64.

[12] 蔣海昆， 曲延軍， 李永莉， 等. 中國大陸中強地震余震序列的部分統計特征[J]. 地球物理學報， 2006， 49（4）: 110-117.

[13] 劉正榮. b值特征的研究[J]. 地震研究， 1995， 18（2）:168-173.

[14] 刁桂苓， 趙英萍. 大同晚期強余震前震源機制的一致性特征 [J]. 內陸地震， 2004， 18（3）: 202-206.

[15] 趙英萍， 刁桂苓. 張北強余震前震源機制的一致性特征[J]. 華北地震科學， 2004， 22（1）: 1-4.

Abstract:Using digital seismic waves recorded by Jilin digital seismic network, we calculated the correlation coefficient of spectral amplitude ofML≥2.5 earthquakes of Qianguo strong earthquake swarm occurred in 2013. The results show that the correlation coefficients of spectral amplitude of Qianguo earthquake swarm are generally high values, which are between 0.85 and 0.95. It means more consistency of focal mechanisms of Qianguo earthquake swarm and high level of the stress field in the source region.Perhaps it is the reason why Qianguo earthquake swarm has been continuous occurred. We also have found that the correlation coefficient of spectral amplitude climbed up beforeMS5.8 earthquake, and went down after it. It means that focal mechanisms become more consistency before strong earthquake and less consistency after it. The result agrees with those of previous researches.

Key words:Qianguo strong earthquake swarm；correlation coefficient of spectral amplitude；consistency of focal mechanism；Local stress field

Characteristics of the Correlation Coefficient of Spectral Amplitude of Qianguo Strong Earthquake Swarm in 2013

KANG Jian-hong1，ZHENG Guo-dong2，ZHANG Yu2，ZHANG Hong-yan1
( 1. Earthquake Administration of Jilin Province, Jilin Changchun 130117, China;2. Changchun Seismic Station, Jilin Changchun 130117, China )

P315.6

A

10.13693/j.cnki.cn21-1573.2017.03.017

1674-8565（2017）03-0092-05

吉林省地震局合同制科研課題資助（2016024）

2017-05-09

2017-07-03

康建紅（1979-），女，河北省昌黎縣人，2005年畢業于燕山大學，碩士，工程師，現主要從地震分析預報研究方面的工作。E-mail：13169035@qq.com