周 依 王 想 孫麗娜

1 河北省地震局，石家莊市槐中路262號，050021

據(jù)中國地震臺網(wǎng)中心測定，北京時間2020-07-12 06:38河北省唐山市古冶區(qū)發(fā)生MS5.1地震，震源深度13 km，震中與1976年唐山MS7.8地震相距約28 km，唐山及周邊地區(qū)震感明顯。本次地震發(fā)生在唐山老震區(qū)范圍內(nèi)，是1991年唐山MS5.1地震后該區(qū)發(fā)生的最大地震事件，與區(qū)域地震序列強(qiáng)度衰減趨勢相比，本次地震震級明顯偏高，表明老震區(qū)中等地震活動出現(xiàn)活躍態(tài)勢，開展地震的發(fā)震構(gòu)造分析有助于了解唐山老震區(qū)的地震活動特征。本文利用河北區(qū)域數(shù)字地震臺網(wǎng)提供的地震資料，采用CAP波形反演方法求取主震震源機(jī)制解，同時利用近震sPL深度震相波形擬合方法，得到精確的主震震源深度，進(jìn)而探討唐山MS5.1地震的發(fā)震構(gòu)造。

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

唐山老震區(qū)位于陰山-燕山褶皺帶與冀魯斷塊拗陷的交匯處，被NE向的寧河-昌黎斷裂、豐臺-野雞坨斷裂及NW向的薊運(yùn)河斷裂、灤縣-樂亭斷裂圍限，中間發(fā)育唐山斷裂帶(圖1)。研究區(qū)基底為下古生界變質(zhì)巖系，蓋層為古生界沉積巖系，地表發(fā)育第四系松散層[1-2]。地震活動圖像顯示，唐山老震區(qū)的主要地震構(gòu)造為唐山斷裂帶、薊運(yùn)河斷裂、灤縣-樂亭斷裂和盧龍斷裂[3]。2000年以來的小震精定位結(jié)果顯示，在唐山斷裂帶北側(cè)，沿徐家樓斷裂出現(xiàn)了NNE向密集的小震活動，被認(rèn)為是由徐家樓斷裂的深部活動引起，且推測該斷裂可能為近年來該區(qū)多次MS4.0以上中等強(qiáng)度地震的發(fā)震斷層[4]。張素欣等[5]根據(jù)1996年以來唐山老震區(qū)ML4.0以上地震的空間分布推測，在唐山斷裂帶北段存在NW向斷裂活動，可能為未來中等強(qiáng)度以上地震的發(fā)震場所。本次唐山MS5.1地震發(fā)生在豐臺-野雞坨斷裂、唐山斷裂帶和徐家樓斷裂的圍空區(qū)域，該區(qū)域可能存在未探明的隱伏斷裂，地質(zhì)結(jié)構(gòu)復(fù)雜。

圖1 唐山MS5.1地震區(qū)域構(gòu)造背景Fig.1 Regional tectonic settings of the Tangshan MS5.1 earthquake

2 主震震源參數(shù)測定

2.1 CAP波形反演方法測定震源機(jī)制

本文采用CAP波形反演方法[6]，利用近震波形數(shù)據(jù)計算唐山MS5.1地震的震源機(jī)制解?；诠俚萚7]對唐山地區(qū)進(jìn)行速度結(jié)構(gòu)分析時使用的一維速度模型(表1)，選取震中距250 km范圍內(nèi)有較好方位角覆蓋的12個寬頻帶數(shù)字地震臺站數(shù)據(jù)。首先將去除儀器響應(yīng)的三分量地震波形旋轉(zhuǎn)至大圓路徑，得到徑向、切向和垂向的位移記錄；然后將實(shí)際觀測數(shù)據(jù)截斷為Pnl波和S波2個部分，并分別在0.05～0.17 Hz和0.05～0.1 Hz頻率范圍對其進(jìn)行帶通濾波；最后計算各臺站的格林函數(shù)，以5°的網(wǎng)格步長搜索不同震源深度的最佳雙力偶解。通過不同深度下的擬合誤差分布(圖2)可以看出，擬合深度在14 km處的誤差最小，該深度即為反演得到的矩心深度。對比最佳擬合深度14 km處的理論地震波形和實(shí)際觀測波形(圖3)發(fā)現(xiàn)，超過80%分量的擬合相關(guān)系數(shù)大于80%，其中接近70%分量的擬合相關(guān)系數(shù)大于90%，表明擬合程度較高。此次地震的最佳雙力偶震源機(jī)制解的節(jié)面Ⅰ走向143°、傾角83°、滑動角-14°，節(jié)面Ⅱ走向235°、傾角76°、滑動角-173°。為檢驗(yàn)結(jié)果的準(zhǔn)確性，本文同時采用P波初動方法計算震源機(jī)制解，得到節(jié)面Ⅰ走向245°、傾角79°、滑動角-158°，節(jié)面Ⅱ走向151°、傾角69°、滑動角-11°。P波初動方法與CAP波形反演方法計算得到的震源機(jī)制解參數(shù)大致相同，表明結(jié)果較為可靠。

表1 一維速度模型

圖2 唐山MS5.1地震震源機(jī)制反演誤差隨深度的變化Fig.2 Variation of fitting error with depth during the focal mechanism inversion for the TangshanMS5.1 earthquake

紅線為理論地震波形，黑線為實(shí)際觀測波形；波形左側(cè)上方為臺站名，下方數(shù)字表示震中距(單位km)和相對偏移時間(單位s)；波形下方的兩行數(shù)字為理論地震波形相對實(shí)際觀測波形的位移時間(單位s)及相關(guān)系數(shù)(單位%)圖3 唐山MS5.1地震震源機(jī)制解Fig.3 Focal mechanism solution of the Tangshan MS5.1 earthquake

從震源機(jī)制解結(jié)果來看，節(jié)面Ⅰ走向NW，呈左旋走滑運(yùn)動性質(zhì)，節(jié)面Ⅱ走向NE，為右旋走滑，2個節(jié)面均表現(xiàn)為高傾角，表明唐山MS5.1地震發(fā)生在具有高傾角并以走滑為主的斷層上。

2.2 sPL深度震相方法測定震源深度

sPL深度震相對震源深度變化較為敏感，因此利用其求取震源深度可顯著提高結(jié)果精度[8]。當(dāng)從震源出發(fā)的SV波入射到自由表面下方時，有一部分能量會轉(zhuǎn)換為P波；而當(dāng)其臨界入射時，轉(zhuǎn)換P波將沿地表傳播。崇加軍等[9]將該波列及其經(jīng)過多次反射折射形成的波列定義為sPL震相，該震相的徑向分量及垂向分量的振幅較強(qiáng)，切向分量振幅很弱，且與直達(dá)P波(Pg)的到時差和震源深度呈線性關(guān)系[10]。因此，若將不同深度下理論地震波的Pg震相和sPL震相的到時差與實(shí)際觀測波形進(jìn)行對比，挑選出到時差吻合程度最高時的深度，該深度即為擬合的震源深度。

在唐山MS5.1地震記錄中，QIX臺和CLI臺能夠清晰地觀測到sPL深度震相，因此本文基于由表1一維速度模型及CAP波形反演方法計算得到的震源機(jī)制解結(jié)果，采用頻率-波數(shù)(F-K)法模擬了這2個臺站在不同深度下的理論地震波形，經(jīng)過1 Hz低通濾波后，比較理論地震波形與實(shí)際觀測波形的Pg震相和sPL震相到時差的擬合情況，以此確定震源深度。從圖4可以看出，2個臺站觀測到的sPL震相清晰，處于直達(dá)P波和S波之間，其能量主要集中在徑向和垂向分量，切向分量振幅不明顯。在徑向和垂向分量的理論地震波形(圖4(a)、4(b)、4(d)及4(e))中，sPL震相和Pg震相之間的到時差與震源深度呈較好的線性關(guān)系，表現(xiàn)出sPL震相對震源深度變化敏感的特征。同時，當(dāng)震源深度為14 km時，在徑向和垂向分量的理論地震波形中，Pg震相和sPL震相之間的到時差與實(shí)際觀測波形的到時差基本一致，由此判斷唐山MS5.1地震的震源深度約為14 km。該結(jié)果與利用CAP波形反演方法得到的震源深度結(jié)果基本一致。

黑線為理論地震波形，紅線為實(shí)際觀測波形圖4 sPL震相確定唐山MS5.1地震震源深度Fig.4 Focal depth determination for Tangshan MS5.1 earthquake with sPL wave

3 發(fā)震構(gòu)造討論

除震源機(jī)制外，地震序列的空間分布特征也是確定地震發(fā)震構(gòu)造的重要依據(jù)[11]。唐山MS5.1地震發(fā)生后，截至2020-07-27共記錄到ML0.0以上余震107次，利用雙差定位方法[12]對余震序列進(jìn)行重定位(圖5)，發(fā)現(xiàn)唐山MS5.1地震序列分布較為集中，優(yōu)勢展布方向?yàn)镹E向。參考地震烈度等震線分布(圖1)發(fā)現(xiàn)，唐山MS5.1地震等震線呈近橢圓分布，長軸為NE向，與震源機(jī)制解節(jié)面Ⅱ的走向大致相符。因此判定該地震的發(fā)震斷層面為節(jié)面Ⅱ，發(fā)震斷層可能為走向NE、傾向NW的陡傾走滑斷裂。

圖5 唐山MS5.1地震序列重定位震中分布Fig.5 Relocated epicenters of the Tangshan MS5.1 earthquake sequence

重定位后發(fā)現(xiàn)，唐山MS5.1地震震中位于唐山斷裂帶北段東側(cè)和徐家樓斷裂西側(cè)，北側(cè)為豐臺-野雞坨斷裂東段。唐山斷裂帶總體走向NE，徐家樓斷裂走向近NS，豐臺-野雞坨斷裂東段沿NE向逐漸彎折呈近EW向，唐山斷裂帶和豐臺-野雞坨斷裂走向均與節(jié)面Ⅱ相符。研究顯示，唐山斷裂帶為褶皺斷裂帶[13]，其北段斷裂主要包括陡河斷裂、巍山-長山南坡斷裂和唐山-古冶斷裂等，其中陡河斷裂為傾向SE的逆斷層，巍山-長山南坡斷裂向NW陡傾并具有逆沖性質(zhì)，唐山-古冶斷裂為傾向NW的正斷層[14]，而豐臺-野雞坨斷裂東段傾向SE[1]。比對發(fā)現(xiàn)，唐山-古冶斷裂的產(chǎn)狀與唐山MS5.1地震的發(fā)震斷層參數(shù)較為吻合。根據(jù)唐山地區(qū)深地震反射剖面可知，唐山-古冶斷裂向W陡傾并延伸至18～20 km的上、下地殼分界面附近[15]，CAP波形反演和sPL深度震相方法確定的唐山MS5.1地震的震源深度為14 km，可能位于唐山-古冶斷裂的中下段。此外，根據(jù)地質(zhì)調(diào)查和小震活動分析結(jié)果可知，1976年唐山地震以來，陡河斷裂和巍山-長山南坡斷裂不活動或活動幅度較小，區(qū)域地震主要沿唐山-古冶斷裂分布[16]。由此判斷，本次唐山MS5.1地震的發(fā)震斷層可能為唐山-古冶斷裂。

4 結(jié) 語

為研究唐山MS5.1地震的發(fā)震構(gòu)造，本文利用河北區(qū)域數(shù)字地震臺網(wǎng)提供的地震資料，采用CAP波形反演和sPL深度震相方法得到了唐山MS5.1地震的震源參數(shù)。結(jié)果表明，唐山MS5.1地震的發(fā)震斷層走向?yàn)?35°，傾角為76°，滑動角為-173°，震源深度為14 km。綜合地震序列分布特征、地震烈度等震線分布、區(qū)域地質(zhì)構(gòu)造、深地震反射剖面及區(qū)域地震活動特征等資料，判斷此次地震的發(fā)震斷層為唐山-古冶斷裂。

致謝：朱露培教授提供CAP和F-K程序及Waldhauser提供雙差定位程序，文中圖件使用GMT軟件繪制，一并表示感謝。