王亞玲，蔡玲玲，李冬圣，王莉嬋，趙英萍，王曉山

（河北省地震局，石家莊 050021）

0 引言

強(qiáng)震發(fā)生后或者震群活動(dòng)時(shí)，在震中附近或小區(qū)域范圍內(nèi)往往會(huì)有大量余震發(fā)生，其中震級(jí)較小的地震震相往往會(huì)淹沒(méi)在背景噪聲中難以識(shí)別，成為遺漏地震。地震目錄是區(qū)域地震活動(dòng)性研究、地震危險(xiǎn)性分析的重要數(shù)據(jù)資料[1]，其完整準(zhǔn)確性直接影響區(qū)域地震活動(dòng)、地震危險(xiǎn)性的分析，并對(duì)后續(xù)科學(xué)研究的可信性和科學(xué)性產(chǎn)生影響[2-3]，尤其對(duì)發(fā)震構(gòu)造及震后趨勢(shì)分析等研究帶來(lái)直接影響。已有學(xué)者應(yīng)用匹配濾波技術(shù)對(duì)中強(qiáng)地震余震序列完整性和余震分布特征[4-7]等方面進(jìn)行研究，譚毅培等分別對(duì)2013 年3 月涿鹿微震群、2013 年8 月蔚縣震群進(jìn)行遺漏地震檢測(cè)分析其發(fā)震構(gòu)造[8-9]；單連君等基于匹配濾波技術(shù)分析河北懷來(lái)ML3.4震群余震序列[10]；王寧等將遺漏震相檢測(cè)技術(shù)應(yīng)用在2016年9 月唐山震群中[11]。本文使用匹配濾波技術(shù)識(shí)別唐山MS5.1地震的遺漏地震進(jìn)行其發(fā)震構(gòu)造的確定。

唐山斷裂帶處于近EW向陰山-燕山南緣活動(dòng)構(gòu)造帶和NNE向冀魯斷塊拗陷相匯的地方[12]。1976 年唐山7.8 級(jí)地震后，沿唐山斷裂帶余震十分活躍，并有大量關(guān)于該地震的發(fā)震構(gòu)造、余震序列研究成果產(chǎn)出。韋士忠等通過(guò)計(jì)算選頻譜，研究了唐山地震余震序列中1979 年5 月17 日5.1級(jí)地震事件，發(fā)現(xiàn)震源等效半徑較小，應(yīng)力降較高[13]，地震序列的頻譜規(guī)律為認(rèn)識(shí)孕育較大地震的過(guò)程提供信息線索；吳培稚等采用不相等間距的統(tǒng)計(jì)方法，研究唐山的余震序列展布[14]，從中得到一些預(yù)報(bào)的指標(biāo)；萬(wàn)永革等通過(guò)唐山地震序列的應(yīng)力變化研究，表明唐山地震序列中的前面大震會(huì)對(duì)后續(xù)小震的發(fā)生起到了調(diào)制作用[15]；仲秋等通過(guò)研究1976年唐山7.8級(jí)大地震的余震序列，得到相應(yīng)的余震序列持續(xù)時(shí)間公式[16]，研究學(xué)者對(duì)唐山地震余震序列的研究領(lǐng)域及方法是多方面的，并且得到相應(yīng)的數(shù)據(jù)指標(biāo)，為認(rèn)識(shí)主震與余震之間的相互關(guān)系，及對(duì)地震危險(xiǎn)性估計(jì)有一定的意義；另外楊雅瓊等基于255 個(gè)震源機(jī)制解對(duì)唐山地震序列的發(fā)震斷層進(jìn)行分段[17]。而余震展布特征不僅與斷層走向有關(guān)，也可能與斷層深部的幾何展布有關(guān)，因此研究余震幾何分布和構(gòu)造特征具有重要的科學(xué)意義。

本文基于2020 年7 月12 日6 時(shí)38 分唐山MS5.1地震，對(duì)7 月12—16 日的連續(xù)波形通過(guò)匹配濾波技術(shù)進(jìn)行遺漏地震檢測(cè)，估計(jì)遺漏地震的震級(jí)，完善地震序列目錄；并對(duì)相似波形地震進(jìn)行震相校正提高定位精度，結(jié)合震源機(jī)制解推測(cè)發(fā)震構(gòu)造；再利用 SKUA-GOCAD 軟件給出發(fā)震構(gòu)造的初步三維模型，通過(guò)二維分析與三維模型相結(jié)合達(dá)到分析此次地震構(gòu)造活動(dòng)特征的目的，為進(jìn)一步研究該地震及鄰區(qū)地震活動(dòng)性提供基礎(chǔ)數(shù)據(jù)支持。

1 構(gòu)造背景

唐山地區(qū)構(gòu)造復(fù)雜如圖1，受寧河-昌黎斷裂（F6）、榛子鎮(zhèn)斷裂（F5）、灤縣-樂(lè)亭斷裂（F1）和薊運(yùn)河斷裂（F7）的切割，唐山地區(qū)形成了一個(gè) NE 向的菱形地塊[18]。其中央是一條NE向的由陡河斷裂（F2）、巍山-長(zhǎng)山南坡斷裂（F3）和唐山-古冶斷裂（F4）組成的唐山斷裂帶，總體走向NE-NNE，傾向NW、SE，傾角70°～80°，為逆走滑斷層，長(zhǎng)約59 km。斷裂詳情見(jiàn)表1。

圖1 唐山地區(qū)構(gòu)造及MS5.1地震序列分布圖

表1 研究區(qū)域周邊斷裂參數(shù)表[19～24]

2 方法

2.1 模板地震選取

使用模版匹配濾波技術(shù)（matched filter technique，MFT）進(jìn)行遺漏地震檢測(cè)，可實(shí)現(xiàn)在低信噪比情況

下提取地震信息。依據(jù)震相報(bào)告到時(shí)數(shù)據(jù)，挑選余震序列中震級(jí)較大的、信噪比好的地震事件，以S波到時(shí)為中心，從連續(xù)波形數(shù)據(jù)中截取S波到時(shí)前2 s至后2 s波形，篩選波形三分量信噪比的平均值大于3的地震事件作為模板。本文選取7 月12—16 日ML1.0以上震中距最近3 個(gè)臺(tái)站震相清楚的余震（表2）為模板，臺(tái)站選擇遵循標(biāo)準(zhǔn)：①寬頻帶地震計(jì)的臺(tái)站，該類型臺(tái)站記錄波形質(zhì)量較好；②震中距較近的臺(tái)站，信噪比較高。因此，選取了DOH、LUX、QIX 3 個(gè)臺(tái)站參與互相關(guān)計(jì)算。

表2 模板地震選取

每個(gè)地震模板均對(duì)7 月12—16 日的連續(xù)波形進(jìn)行掃描，計(jì)算互相關(guān)系數(shù)，將每個(gè)模板三分向數(shù)據(jù)的互相關(guān)結(jié)果進(jìn)行組合，通過(guò)計(jì)算序列的絕對(duì)離差中位數(shù)（MAD）檢測(cè)遺漏地震，其表達(dá)式為：

式中：Xi代表第i個(gè)互相關(guān)系數(shù)序列；為其平均值。本次運(yùn)算中MAD的倍數(shù)閾值是9[5]。為了減少地震序列因僅有一兩個(gè)臺(tái)站距離震中較近而造成檢測(cè)結(jié)果失誤，還需設(shè)置單臺(tái)互相關(guān)系數(shù)閾值，默認(rèn)值為0.75[10]。

以Eq712064434地震模板的灤縣臺(tái)（LUX）掃描連續(xù)波形為例（圖2），圖中為掃描7 月12 日連續(xù)波形結(jié)果，掃描模板所在的連續(xù)波形互相關(guān)系數(shù)為1。連續(xù)波形通過(guò)閾值檢測(cè)得到的結(jié)果為遺漏地震文件，對(duì)照原地震目錄，序列未記錄的事件作為疑似地震事件，最后拾取遺漏地震事件S波形的最大振幅，并利用其與模板地震的振幅比來(lái)估計(jì)遺漏事件震級(jí)。

圖2 Eq712064434 模板灤縣臺(tái)的互相關(guān)掃描序列

2.2 震相校正

區(qū)域測(cè)震臺(tái)網(wǎng)震相報(bào)告給出的震相到時(shí)存在一定的人工誤差[10]，這種誤差對(duì)于地震序列時(shí)空演化的精細(xì)分析有較大的影響，因此本文利用少數(shù)幾個(gè)震級(jí)較大的地震事件，通過(guò)波形互相關(guān)的方法對(duì)波形相似的地震事件進(jìn)行震相校正，有效地減少誤差，提高定位精度。

分別采用風(fēng)速為0.2，0.6，1.0m/s3種不同風(fēng)速，在物料薄層厚度為5cm、熱風(fēng)溫度為50℃下進(jìn)行高水分小麥的干燥，探討物料含水率變化及風(fēng)速對(duì)干燥速度的影響。

以Eq712070225模板的灤縣臺(tái)（LUX）為例（圖3），對(duì)相似波形地震事件進(jìn)行震相校正。本程序得到的到時(shí)信息與臺(tái)網(wǎng)人工讀取震相無(wú)關(guān)，在一定程度上消除了地震事件在人工定位時(shí)產(chǎn)生的相對(duì)誤差，從而提高了定位結(jié)果的精度。

圖3 灤縣臺(tái)（LUX）震相標(biāo)定過(guò)程示意圖

3 遺漏地震事件檢測(cè)及發(fā)震構(gòu)造分析

3.1 遺漏地震發(fā)震時(shí)刻和震級(jí)

基于每個(gè)地震模板所識(shí)別出的地震事件的發(fā)震時(shí)刻，綜合所有模板檢測(cè)出來(lái)的地震事件，對(duì)比目錄結(jié)果便可以找出遺漏地震。由于河北測(cè)震臺(tái)網(wǎng)工作方式的改革提升，現(xiàn)采取截取實(shí)時(shí)波形的方式查看地震余震序列，一般肉眼可識(shí)別的地震已實(shí)時(shí)快報(bào)，并編寫目錄。針對(duì)此次地震來(lái)說(shuō)余震序列并不發(fā)育，因此利用遺漏檢測(cè)技術(shù)掃描7 月12—16日連續(xù)波形共檢測(cè)到的地震事件共8 個(gè)，遺漏地震震級(jí)、發(fā)震時(shí)刻、經(jīng)緯度如表3所示。其中：ML0.0～0.5地震共有7 個(gè)，ML1.4的1 個(gè)。ML0.0～0.5地震編目條目由36增至42 條，表明ML0.0～0.5之間的地震目錄的完善性得到明顯的提升。雖然大部分檢測(cè)出的遺漏地震較小，但是對(duì)分析此次MS5.1地震震中展布的幾何特征具有一定意義，也為構(gòu)建三維構(gòu)造模型提供更為完善的數(shù)據(jù)支持。

表3 遺漏地震發(fā)震時(shí)刻和震級(jí)

3.2 精定位結(jié)果及發(fā)震構(gòu)造分析

精確定位后的主震及其余震序列空間展布和震源機(jī)制解，對(duì)于確定發(fā)震構(gòu)造、判定未來(lái)地震趨勢(shì)具有重要意義[25-26]。利用CAP[27]方法反演主震的震源機(jī)制解，該方法是將近震觀測(cè)波形分解成Pnl和面波，分別進(jìn)行擬合，對(duì)3 個(gè)斷層面參數(shù)的各自取值區(qū)間進(jìn)行格點(diǎn)的取值，將得到的理論地震圖與實(shí)際觀測(cè)波形進(jìn)行擬合，波形擬合最好的一組參數(shù)即為求得的震源斷層面解。反演本次唐山MS5.1地震的震源機(jī)制解采用的速度模型表4，計(jì)算結(jié)果如圖4b。共選取10 個(gè)臺(tái)站參與，其中體波和面波各部分波形擬合互相關(guān)系數(shù)大于80%的占70%（圖5a），顯示此次計(jì)算結(jié)果擬合較好，結(jié)果可靠，深度擬合殘差顯示（圖5b），擬合深度為14 km時(shí)殘差最小，此深度即為最佳震源深度。

圖5 唐山MS5.1地震參與計(jì)算臺(tái)站的波形擬合結(jié)果及深度擬合殘差分布圖

本文對(duì)余震序列采用雙差定位方法[28]進(jìn)行精定位，速度模型采用的是郭蕾等[29]使用的一維P波速度模型（表4）。為了便于對(duì)比精定位前后的震中分布趨勢(shì)，刪除陡河臺(tái)、灤縣臺(tái)單臺(tái)定位目錄，圖4a給出了精定位前地震目錄（包括遺漏地震）的91（含主震）個(gè)地震事件的震中分布圖。對(duì)經(jīng)過(guò)波形互相關(guān)震相校正的地震事件進(jìn)行重新定位，精定位結(jié)果包含84 個(gè)地震事件。對(duì)未能標(biāo)定3 個(gè)以上臺(tái)站震相到時(shí)的地震事件，其震中的位置可以定在與其互相系數(shù)最大的模板地震震中處，最終震中分布顯示如圖4b。由圖中精定位結(jié)果可知，唐山余震序列分布趨勢(shì)更為明顯，主要集中在唐山-古冶斷裂的東北部；震中NE向的優(yōu)勢(shì)展布也更為明顯，與主震震源機(jī)制解的節(jié)面Ⅰ走向基本一致（表5）。推測(cè)節(jié)面Ⅰ為本次唐山MS5.1地震的發(fā)震斷層面，該地震是發(fā)生在走向NE、傾向NW、高傾角的右旋兼走滑斷層上，與唐山-古冶斷裂（F4）的走向、傾向、傾角基本吻合，由此推測(cè)唐山-古冶斷裂可能為其發(fā)震斷裂。

表4 唐山地區(qū)速度模型

圖4 序列精定位前后震中分布及震源機(jī)制解結(jié)果圖

表5 MS5.1主震震源機(jī)制解結(jié)果

3.3 基于GOCAD發(fā)震斷層初步三維模型

3.3.1 余震序列三維空間展布模型

由上文可知，唐山MS5.1地震發(fā)生后截止至2020 年7 月16 日，河北區(qū)域測(cè)震臺(tái)網(wǎng)記錄可定位余震82 條，經(jīng)過(guò)遺漏檢測(cè)之后增至90 條。結(jié)合圖4b所呈現(xiàn)的斷層地表幾何特征與余震展布關(guān)系，對(duì)余震空間分布的最大優(yōu)勢(shì)分布進(jìn)行分析和預(yù)處理[30-31]，基于GOCAD軟件初步建立斷層在三維空間的地表展布模型，實(shí)現(xiàn)地震三維空間展布（圖6）。

圖6 清晰展示了余震序列在二維地表展布、三維空間分布特征，更加直觀地顯示余震空間上的優(yōu)勢(shì)展布方向。

3.3.2 三維斷層幾何模型

建立震區(qū)的三維構(gòu)造模型，需要的數(shù)據(jù)資料包括以下幾個(gè)方面：斷層地表展露的幾何特征、斷層走向、傾角、傾向以及地下延展情況的詳細(xì)描述、相關(guān)地震震中展布及震源機(jī)制解約束條件等。本文根據(jù)已有的斷裂資料（圖1），初步建立研究震區(qū)灤縣-樂(lè)亭斷裂（F1）、陡河斷裂（F2）、巍山-長(zhǎng)山南坡斷裂（F3）、唐山-古冶斷裂（F4）以及榛子鎮(zhèn)斷裂（F5）、寧河-昌黎斷裂（F6）6條斷裂的三維模型（圖6）。

基于已知的斷裂數(shù)據(jù)資料，在圖6的基礎(chǔ)上通過(guò)GOCAD軟件的“make suface”功能造面成圖（圖7）。

圖6 發(fā)震區(qū)域斷層及余震三維空間展布

由圖7所示：唐山-古冶斷裂（F4）、陡河斷裂（F2）、巍山-長(zhǎng)山南坡斷裂（F3）的空間走向較為準(zhǔn)確地表達(dá)在三維空間上，實(shí)現(xiàn)對(duì)發(fā)震斷裂的三維模擬與刻畫。結(jié)合地震震源機(jī)制解結(jié)果及余震序列分布推測(cè)節(jié)面Ⅰ為同震斷裂面，唐山-古冶斷裂（F4）與節(jié)面Ⅰ的走向、傾角基本一致，結(jié)合地震的震中位置推測(cè)此次地震可能發(fā)生在唐山-古冶斷裂（F4）上。在此基礎(chǔ)上結(jié)合圖7模擬的三維空間展示，可以看出余震空間展布趨勢(shì)與唐山-古冶斷裂（F4）NE至NEE走向延展具有一致性。據(jù)此，可以更直觀、更有力地推測(cè)唐山-古冶斷裂（F4）即為河北唐山MS5.1地震的發(fā)震斷裂。

圖7 研究區(qū)主要斷裂帶的三維幾何模型

4 討論與結(jié)論

本文通過(guò)匹配濾波技術(shù)，對(duì)2020 年7 月12—16 日唐山MS5.1 地震震中附近臺(tái)站記錄的連續(xù)波形進(jìn)行遺漏地震檢測(cè)、震相校正來(lái)完善地震序列目錄，依據(jù)序列精定位、主震震源機(jī)制解分析發(fā)震構(gòu)造，并基于GOCAD軟件初步給出了發(fā)震區(qū)斷層的三維構(gòu)造模型。

1）運(yùn)用匹配濾波技術(shù)對(duì)余震進(jìn)行遺漏檢測(cè)，檢測(cè)結(jié)果發(fā)現(xiàn)了遺漏的地震事件，因此測(cè)震臺(tái)網(wǎng)人員可以通過(guò)此方法拾取震相，提高地震目錄的完整性。

2）通過(guò)波形互相關(guān)方法對(duì)震相校正，不僅可以估計(jì)遺漏地震的震級(jí)，還可以利用少數(shù)幾個(gè)震級(jí)較大的地震事件對(duì)于波形相似的地震進(jìn)行震相校正，有效地減少人工拾取震相誤差，大大提高定位精度。

3）精定位后的余震序列空間分布趨勢(shì)更為明顯，地震序列主要分布在唐山斷裂帶中NE向唐山-古冶斷裂附近，序列優(yōu)勢(shì)分布方向?yàn)镹E向，分析節(jié)面Ⅰ為同震斷層面，與唐山-古冶斷裂東北段走向、傾向、傾角大體一致，初步推測(cè)MS5.1地震發(fā)震構(gòu)造為唐山-古冶斷裂。

4）基于GOCAD軟件，結(jié)合已知的地質(zhì)背景、震源機(jī)制解等資料，初步建立發(fā)震區(qū)的三維構(gòu)造模型；在結(jié)論3的基礎(chǔ)上，進(jìn)行三維空間模擬更為直觀展示了余震序列的空間分布特征；結(jié)合唐山-古冶斷裂空間延展分布趨勢(shì)，展現(xiàn)出一致性。作為輔助工具為此次地震的發(fā)震構(gòu)造判定提供更為直觀的論據(jù)支撐，為驗(yàn)證發(fā)震構(gòu)造推斷提供一種可視化的有效途徑。

基于GOCAD建立的三維構(gòu)造模型為判斷發(fā)震機(jī)理、震源斷層的空間描述提供了一條可借鑒的有效途徑。此外，由于唐山MS5.1地震余震序列不發(fā)育，樣本數(shù)據(jù)量較小，無(wú)法有效、合理地建立三維斷層面模型，這也表明初始的三維構(gòu)造模型具有一定的局限性。未來(lái)成熟的三維模型，需要大量數(shù)據(jù)樣本，不斷提高各類數(shù)據(jù)質(zhì)量，逐漸增加各種成熟的地震活動(dòng)解譯資料及探測(cè)資料，發(fā)展構(gòu)造定量解析方法等多元數(shù)據(jù)條件約束的基礎(chǔ)上進(jìn)行豐富完善。