CAP方法反演2016年唐山ML4.3地震震源機制解

朱音杰，劉 檀，趙英萍，劉 新，李冬圣，丁 成

(1.河北省地震局石家莊中心臺，石家莊 050021；2.河北省地震局，石家莊 050021；3.河北省地震局豐寧地震臺，河北 承德 067000)

CAP方法反演2016年唐山ML4.3地震震源機制解

朱音杰1，劉 檀1，趙英萍2，劉 新2，李冬圣2，丁 成3

(1.河北省地震局石家莊中心臺，石家莊 050021；2.河北省地震局，石家莊 050021；3.河北省地震局豐寧地震臺，河北 承德 067000)

利用河北臺網(wǎng)及臨近區(qū)域臺網(wǎng)的地震波形資料，采用CAP方法反演了2016年9月10 日河北唐山ML4.3地震的震源機制解和震源深度并將反演結(jié)果作為已知輸入，利用F-K方法計算理論地震波形，與觀測記錄進行對比，以驗證結(jié)果的可靠性。結(jié)果顯示，此次唐山地震為右旋—正傾斷層，矩震級為4.4，最佳雙力偶解為節(jié)面I：走向238.0°，傾角44.0°，滑動角-143.0°；節(jié)面Ⅱ：走向119.5°，傾角65.3°，滑動角-52.4°，最佳震源深度為6.5 km。

唐山地震；CAP方法；震源機制解；余震

0 引言

2016年9月10日河北唐山古冶區(qū)及開平區(qū)交界處發(fā)生ML4.3地震，市區(qū)有明顯震感。本文利用河北臺網(wǎng)及臨近區(qū)域臺網(wǎng)寬頻帶數(shù)字地震儀的地震波形記錄，對唐山地震精確定位并反演其震源機制解。根據(jù)地震活動性和地質(zhì)構(gòu)造特征探討唐山ML4.3地震的發(fā)震構(gòu)造，并對其是否屬于唐山大地震余震進行討論分析。

在震源機制的研究中，學(xué)者們在這方面積累了大量寶貴的經(jīng)驗。Zhao和Helmberger等[1]提出的CAP(Cut and Paste)波形反演求取震源機制解的方法，采用近震寬頻帶波形資料，利用分段的波形擬合，通過濾波去除長周期的地脈動及漂移，充分利用波形中所包含的震源信息，彌補了P波初動法在臺站分布不均、P波極性判斷及計算地震深度的問題，使得反演結(jié)果更加穩(wěn)定、可靠。CAP方法獨特的波形分離、浮動擬合技術(shù)，即使在速度模型不佳、格林函數(shù)不理想、甚至是臺站分布不均勻的情況下，仍能獲得不錯的波形擬合效果[2]。

1 發(fā)震背景

1.1 地質(zhì)構(gòu)造背景

唐山地區(qū)位于華北板塊的北邊界陰山EW向構(gòu)造帶南緣—燕山褶斷帶和冀中塊體東邊界帶—滄東斷裂的交匯部位[3]。唐山極震區(qū)構(gòu)造顯著的特點是：由寧河—昌黎斷裂、豐臺—野雞坨大斷裂(現(xiàn)榛子鎮(zhèn)斷裂)、灤縣—樂亭斷裂和薊運河斷裂的切割形成一個NEE向的菱形塊體。

唐山斷裂帶發(fā)育在菱形塊體內(nèi)部的蓋層中，并由3條近平行的斷層組成(圖1)：

1)陡河斷層。全長約20 km，緊靠地壘構(gòu)造的西側(cè)展布，東北段為傾向NW的正斷層，西南段由2條平行的小斷層組成，斷面都傾向NW，西邊一條為正斷層，東邊一條為逆斷層。

2)唐山—巍山—長山南坡斷層。全長約21 km，由一些斷續(xù)出露的斷層組成，走向NE，傾向NW，傾角達80°，以擠壓逆沖性質(zhì)為主，斷層多沿地層層面分布。

3)唐山—古冶斷層。長約32 km，總體走向NE至NEE，其中西南段走向NE 30°，東北段走向NE50°，傾向SE，傾角50°～60°，屬于走滑性質(zhì)。

由震中位置可以推測，此次唐山ML4.3地震的發(fā)震斷層可能為唐山—古冶斷層，與唐山MS7.8地震發(fā)震斷層一致。

1.2 地震活動性背景

一次較大地震發(fā)生后，余震序列總的趨勢是隨時間有起伏地逐漸減弱[4]。余震是在大地震，即中強度以上地震之后在地球內(nèi)部發(fā)生主震的同一地方發(fā)生。因發(fā)震位置的不同，也能造成余震延續(xù)時間的長短不一，例如板緣地震衰減比較快，可能只需幾天或十幾天時間，板內(nèi)地震衰減會比較慢，可能延續(xù)幾十年乃至上百年，唐山地震余震序列便如上所述。所以判斷一次地震是否屬于余震可以從空間、時間以及破裂機制等方面進行討論。

圖1 唐山ML4.3級地震區(qū)域構(gòu)造圖

唐山7.8地震的余震序列總體上呈衰減特征，但是衰減過程中有階段性起伏，這個唐山ML4.3地震是近期一次新的起伏。唐山7.8地震已經(jīng)過去40年，但是其余震序列仍在持續(xù)。根據(jù)仲秋等[5]的研究，該地震的余震序列還將持續(xù)100年左右。

2004年以來，唐山老震區(qū)的ML4.0以上地震呈現(xiàn)出2年左右的準(zhǔn)周期特征，準(zhǔn)周期約為22～26個月，同時兼有雙震特征[6](圖2)。

圖2 研究區(qū)域內(nèi)4.0級及以上地震M-t圖

從圖2中可以看出，從2004年到2014年中基本呈現(xiàn)出這個特征，但在2015年09月發(fā)生ML4.7地震，打亂了這個規(guī)律。并且從地震活動性來看，2015年之后3.5級及以上地震明顯增加，地震活動性增強，具體原因有待更深一步的研究。

2 資料與求解

2.1 資料選取及數(shù)據(jù)處理

通常情況下，區(qū)域速度結(jié)構(gòu)的橫向不均勻性隨著震中距增加而增大，區(qū)域地震臺站記錄到的震相復(fù)雜性也隨著震中距增大而增大。因此，在利用區(qū)域地震波形反演中強地震的震源機制解時，較多地采用震中距250 km范圍內(nèi)的臺站記錄[7]。

CAP方法反演震源機制解主要通過對比給出震源機制解的理論波形與實際觀測波形，擬合殘差最小的機制解為最終反演結(jié)果。在研究中將原始速度記錄扣除儀器響應(yīng)，并旋轉(zhuǎn)成Z-R-T分量，并將其分成Pnl和Snl兩部分。將Pnl部分經(jīng)帶寬為0.05～0.2 Hz、Snl部分經(jīng)帶寬為0.05～0.1 Hz的4階Butterworth帶通濾波器進行濾波，提高信噪比。根據(jù)Pnl和Snl部分的波形特點，設(shè)置2者相對權(quán)重為2∶1，更好地約束2者的振幅比對震源深度及機制解的影響。在計算過程中，采用距離影響因子、頻率—波數(shù)法(F-K方法)等解決了震中距過小及給定速度模型下各個臺站間格林函數(shù)的計算結(jié)果問題的影響，得到各種頻率下的體波和面波波形，由震源函數(shù)及格林函數(shù)合成理論波形。經(jīng)上述處理，得到每個臺Pnl部分的垂直分量、徑向分量和Snl部分的3個分量。

由于CAP方法對速度結(jié)構(gòu)依賴相對較小，綜合研究區(qū)域內(nèi)地形因素，采用與河北臺網(wǎng)MSDP軟件對應(yīng)一致的速度模型(表1)。

表1 研究地區(qū)地殼速度模型

2.2 震源機制解結(jié)果分析

根據(jù)近震臺站記錄和華北地區(qū)的速度模型，利用CAP方法反演得到的最佳震源機制解、理論波形與實際觀測波形的對比(圖3)。本文采用的是震中距200 km范圍內(nèi)的臺站。在圖3顯示的理論合成波形與實測波形的比較中，9個臺站Pnl部分和面波部分共42個震相，相關(guān)系數(shù)大于0.8的有31個，達到了70%以上，屬于強度相關(guān)。其中，相關(guān)系數(shù)大于0.9的有23個，在強度相關(guān)中超過了70%，反演方差為9.547×10-4，反演結(jié)果理論地震波形圖與實測地震波形圖吻合得較好。該深度對應(yīng)的雙力偶解為最佳震源機制解(表2)，其結(jié)果為矩震級4.4。

圖3 CAP方法反演的震源機制解理論波形圖(紅色)與實測波形圖(黑色)對比

表2 唐山ML4.3級地震反演結(jié)果

根據(jù)反演結(jié)果與研究區(qū)域地質(zhì)構(gòu)造圖推測，節(jié)斷面Ⅰ為發(fā)震斷層面；由震中位置圖可以推斷此次地震發(fā)生在唐山—古冶斷裂帶。由于唐山—古冶斷裂帶呈現(xiàn)NE走向，唐山—古冶斷裂帶的東北段走向為50°，與震源機制解結(jié)果一致，因此得出此次地震發(fā)生在唐山斷裂帶中NE向唐山—古冶構(gòu)造帶東北段上。反演結(jié)果中壓應(yīng)力軸方向為近WEE向，張應(yīng)力軸近NS向，根據(jù)斷層走向分析得到，發(fā)震斷層為右旋—正傾斷層。

由最終選取的臺站分布圖看(圖4)， 臺站間最大空隙角小于180°,臺站個數(shù)為9個，臺站分布參數(shù)deltaU=0.30，小于0.35，該地震利用CAP方法反演得到的機制解可以采用。

圖4 參與計算臺站分布圖

為了驗證結(jié)果的可靠性，同時利用初動及振幅比方法[8]計算該地震的震源機制解。初動及振幅比方法主要以計算直達波PG、SG振幅比為主，初至波初動方向為輔進行震源機制解計算。只要有準(zhǔn)確的直達波觀測振幅比及觀測值歸算到震源球面上的位置正確，振幅比觀測對震源機制參數(shù)有較強的約束能力[9]。通過初動及振幅比方法得到的震源機制解為節(jié)面Ⅰ走向為229.0°，傾角為36.0°，滑動角為-173.0°；節(jié)面Ⅱ走向為134.0°，傾角為86.0°，滑動角為-54.0°(圖5)。

a 初動及振幅比方法 b CAP方法圖5 分別利用初動及振幅比和CAP方法計算震源機制解沙灘球?qū)Ρ葓D

表3 2種方法計算震源機制解結(jié)果對比

從圖5及表3可以看出，利用初動—振幅比方法計算震源機制解與CAP方法計算結(jié)果基本一致，得到的斷層類型均為右旋—正傾滑動。2種方法略有差異，在初動—振幅比方法中對振幅選取精確度要求較高。

劉桂萍等在計算唐山大地震應(yīng)力場變化一文中，對唐山7.8級主震斷層位錯模型參數(shù)中總結(jié)到主震斷層走向為229°，傾角80°[10]，與此次地震震源機制節(jié)面Ⅰ走向基本一致，因此也可以推斷2016年9月10日唐山ML4.3級地震為唐山大地震的余震序列。

2.3 震源深度

一般來說，利用CAP方法可以獲得較準(zhǔn)確的震源深度。然而震源深度可能受到反演過程中多種因素的影響，比如面波和體波的相對權(quán)重，以及不同的濾波頻段等因素[11]。圖6給出了唐山ML4.3地震震源機制解隨不同深度變化的取值情況，可以看出，震源深度在6.5 km處震源機制解的反演方差相對較小，且每個深度反演得到的震源機制節(jié)面參數(shù)基本相同，反演結(jié)果比較穩(wěn)定。

圖6 不同深度誤差和震源機制隨不同震源深度的變化

在震源深度方面，如果某一臺站離震中較近或者就在震中區(qū)，則該臺震源距與震源深度相近[12]。根據(jù)震源距公式：

D=VΦ×Δt

(1)

VΦ=(VP·VS)/(VP-VS)

(2)

式中：VΦ為虛波速度，由速度模型中的第2層計算，取6.01 km/s；Δt為某一臺站記錄的S波與P波的到時差，由最近的臺站陡河臺(39.75°N，118.33°E)的觀測記錄得到，取1.27 s。據(jù)此推算，震源深度約為7.6 km，而CAP反演該地震的最佳擬合深度為6.5 km，反演結(jié)果略小于理論推算值。一方面可能是濾波頻段的選擇及人工標(biāo)注震相走時差造成的誤差；另一方面，CAP反演的是矩心深度，反映的是能量釋放最大處的深度，物理意義上有一定區(qū)別，因此會造成誤差。

3 結(jié)論

本文利用河北臺網(wǎng)及臨近臺網(wǎng)寬頻帶地震儀的近震波形數(shù)據(jù)，采用CAP方法對2016年9月10日唐山ML4.3地震進行反演，得出以下結(jié)論。

1)此次地震的震源機制解節(jié)面走向與唐山大地震主震斷層基本一致，結(jié)合地震活動性背景來看，此次地震應(yīng)屬于唐山大地震的余震序列。發(fā)震斷層可能是唐山斷裂帶中NE向唐山—古冶斷裂。

2)該反演地震在矩心深度為6.5 km時震源機制解反演方差達到較小值，波形相關(guān)系數(shù)絕大多數(shù)都大于0.7，表明反演結(jié)果可靠。

3)CAP方法反演得到的ML4.3地震最佳震源機制解結(jié)果為：節(jié)面Ⅰ，走向為238.0°，傾角為44.0°，滑動角為-143.0°；節(jié)面Ⅱ，走向為119.5°，傾角為65.3°，滑動角為-52.4°。節(jié)面Ⅰ與唐山—古冶斷裂帶東北段走向基本一致，故推斷此斷裂帶為唐山ML4.3地震的發(fā)震構(gòu)造，震源斷層面取向NEE，NEE向區(qū)域構(gòu)造應(yīng)力為孕震和發(fā)震的主要力源。發(fā)震斷層為右旋—正傾斷層。

致謝：感謝審稿老師給予的寶貴意見。

[1] Zhao L S, Helmberger D V. Source estimation from broadband regional seismograms[J]. Bulletin of the Seismological Society of America, 1994, 84(1): 91-104.

[2] 段剛. CAP方法反演福建仙游ML5.0、ML4.5級地震震源機制解[J]. 地球物理學(xué)報, 2016, 31(4): 1432-1437.

[3] 國家地震局《一九七六年唐山地震》編輯組. 一九七六年唐山地震[M]. 北京: 地震出版社, 1982.

[4] 曲均浩, 蔣海昆. 余震活動機理研究綜述[J]. 中國地震, 2012, 28(2): 109-120.

[5] 仲秋, 史保平. 1976年唐山MS7.8地震余震序列持續(xù)時間及對地震危險性分析的意義[J]. 地震學(xué)報, 2012, 34(4): 494-508.

[6] 張躍剛, 王玉珍, 尹寶軍, 等. 2012年5月28日河北省唐山4.8級地震[J]. 中國地震, 2013, 29(2): 219-229.

[7] 呂堅, 鄭秀芬, 肖健, 等. 2012年9月7日云南彝良MS5.7、MS5.6地震震源破裂特征與發(fā)震構(gòu)造研究[J]. 地球物理學(xué)報, 2013, 56(8): 2645-2654.

[8] 梁尚鴻, 李幼銘, 束沛鎰, 等. 利用區(qū)域地震臺網(wǎng)、振幅比資料測定小震震源參數(shù)[J]. 地球物理學(xué)報, 1984, 27(3): 249-257.

[9] 李冬圣, 蔡玲玲, 李小軍, 等. 河北贊皇小震群震源斷層分析[J]. 地震地磁觀測與研究, 2015, 36(6): 19-25.

[10] 劉桂萍, 傅征祥. 1976年7月28日唐山7.8級地震觸發(fā)的區(qū)域地震活動和靜應(yīng)力場變化[J]. 地震學(xué)報, 2000, 22(1): 17-26.

[11] 謝祖軍, 鄭勇, 倪四道, 等. 2011年1月19日安慶ML4.8地震的震源機制解和深度研究[J]. 地球物理學(xué)報, 2012, 55(5): 1624-1634.

[12] 曲均浩, 劉瑞峰, 李金, 等. CAP方法反演2014年山東乳山M4.2、M4.0地震震源機制解[J]. 地震工程學(xué)報, 2014, 36(4): 1076-1080.

Focal Mechanism Inversion of the 2016 TangshanML4.3 Earthquake Using the CAP Method

ZHU Yin-jie1, LIU Tan1, ZHAO Ying-ping2, LIU Xin2, LI Dong-sheng2, DING Cheng3

(1. Huangbizhuang Seismic Station, Earthquake Administration of Hebei Province, Shijiazhuang 050021, China;2. Earthquake Administration of Hebei Province, Shijiazhuang 050021, China; 3. Fengning Seismic Station, Earthquake Administration of Hebei Province, Chengde 067000, China)

On September 10, 2016, aML4.3 earthquake occurred in Tangshan city which belongs to the Tangshan earthquake sequence. In this study, on basis of waveform data from digital seismic network of Hebei province and its adjacent regions, the focal mechanism solution and focal depth of the earthquake is inverted using the CAP method. The results of focal mechanism and focal depth are confirmed by comparing the observed records with the synthetic waveform computed using the F-K method. The results show that the earthquake mechanism is a strike-slip fault with normal thrust component. The moment magnitude is 4.4. The strike, dip and slip angles of two nodal planes of the best double couple solution is 238.0°, 44.0° , -143° and 119.5° , 65.3° , -52.4° , the estimated focal depth is 6.5 km.

Tangshan earthquake; CAP method; focal mechanism; aftershock

朱音杰，劉檀，趙英萍，等. CAP方法反演2016年唐山ML4.3地震震源機制解[J]．華北地震科學(xué)，2017，35(2):50-55.

2016-12-02

測震臺網(wǎng)青年骨干培養(yǎng)專項(20150404)

朱音杰(1986—)，男，山東文登人，助理工程師，現(xiàn)主要從事地震監(jiān)測工作． E-mail：zhuyinjie_361@163.com

P315.33

A

1003-1375(2017)01-0050-06

10.3969/j.issn.1003-1375.2017.01.008