孟令媛 史海霞 高金哲

1)中國(guó)地震臺(tái)網(wǎng)中心，北京市西城區(qū)三里河南橫街5號(hào) 100045

2)吉林省地震局，長(zhǎng)春 132000

0 引言

2013年10月31日11時(shí)03分吉林省松原市前郭縣(44.6°N，124.2°E)發(fā)生MS5.5地震，震源深度8km，當(dāng)天11時(shí)10分再次發(fā)生MS5.0地震。同年11月22日16時(shí)18分又發(fā)生MS5.3地震，次日06時(shí)04分及06時(shí)32分再次發(fā)生MS5.8、MS5.0地震。11月22、23日發(fā)生的地震3次MS≥5.0地震均位于10月31日MS5.5、MS5.0地震的震源區(qū)，這一系列地震屬于震群活動(dòng)。截至2014年5月31日前郭震群共記錄到MS≥0.0地震872次，其中MS≥3.0的25次，包含MS3.0～3.9地震13次，MS4.0～4.9地震7次，MS5.0～5.9地震5次，最大地震為2013年11月23日MS5.8地震(圖1;表1、2)。

圖1 截至2014年5月31日前郭震群M S≥5.0地震序列及此序列中M S≥0.0的余震分布(a)、前郭震群和CN2臺(tái)的位置示意圖(b)

表1 2013年前郭M S 5.8震群震源機(jī)制參數(shù)

表2 2013年前郭震群M S≥4.0地震目錄*

2013年前郭MS5.8震群的震中距離長(zhǎng)春市平均約120km，距哈爾濱市約230km，前郭震群5次MS≥5.0地震事件發(fā)生時(shí)，長(zhǎng)春市和哈爾濱市均有強(qiáng)烈震感。前郭地區(qū)在24天內(nèi)連續(xù)發(fā)生5次MS≥5.0地震，已構(gòu)成典型的爆發(fā)式強(qiáng)震群活動(dòng)。前郭MS5.8震群是2006年乾安-前郭MS5.0地震后吉林松原地區(qū)發(fā)生的又一次強(qiáng)震群活動(dòng)。歷史上公元1119年吉林省松原地區(qū)發(fā)生過一次級(jí)地震，震中(44.9°N，124.8°E)位于前郭縣套浩太附近，距2013年前郭震群約60km，1119年松原級(jí)地震是東北地區(qū)迄今為止有記載的震級(jí)最大的淺源地震，極震區(qū)烈度達(dá)IX度，現(xiàn)有研究成果顯示其發(fā)震構(gòu)造位于前郭-扶余一帶(李傳友等，1999)。

2013年前郭震群臨近NE向的扶余-肇東斷裂附近的深部高速體邊界帶，扶余-肇東斷裂控制了松遼盆地內(nèi)部的次級(jí)塊體并被多條NW向的次級(jí)斷裂切割，構(gòu)造相對(duì)較為破碎。2006年3月31日，距本次前郭震群以西約10km處曾發(fā)生乾安-前郭MS5.0地震(44.7°N，124.0°E)，該地震也曾經(jīng)造成震區(qū)部分房屋的倒塌損毀和較大的經(jīng)濟(jì)損失，其余震在空間上呈現(xiàn)NW、NNW或近NS向分布，主要受NE向的大安-德都斷裂和NW向的通榆-長(zhǎng)春斷裂的控制及影響(盤曉東等，2007)。

本文針對(duì)2013年前郭MS5.8震群的震源參數(shù)特征展開相關(guān)研究，選取震群中5次MS≥5.0主要地震事件，基于國(guó)家級(jí)臺(tái)站長(zhǎng)春臺(tái)(簡(jiǎn)稱CN2)的地震動(dòng)觀測(cè)記錄，計(jì)算5次主要地震事件的單臺(tái)地震動(dòng)譜;進(jìn)一步結(jié)合震群型地震震源譜的特點(diǎn)，對(duì)比前郭震群中5次主要地震事件的遠(yuǎn)場(chǎng)地震波輻射能的大小并分析CN2單臺(tái)視應(yīng)力的特征;最后，針對(duì)前郭震群的起伏特征、地震波能量釋放等給出基于震源譜特征的相關(guān)解釋。

1 震群的震源譜參數(shù)特征

在地震學(xué)中，地震記錄在時(shí)間域的一般表達(dá)式為

式中，Ω為震源過程的函數(shù);P是幾何擴(kuò)散函數(shù);F是非彈性衰減函數(shù);R是由臺(tái)站下方的淺層結(jié)構(gòu)所決定的場(chǎng)地放大函數(shù);I是地震儀的響應(yīng)函數(shù)。式(1)中地震動(dòng)觀測(cè)記錄在時(shí)間域可表達(dá)為震源過程、地震波傳播路徑及臺(tái)站影響等褶積的結(jié)果。進(jìn)一步，將式(1)轉(zhuǎn)換為頻率域，則有

式(2)實(shí)現(xiàn)了由式(1)的時(shí)間域褶積計(jì)算到頻率域的點(diǎn)乘計(jì)算，簡(jiǎn)化了計(jì)算過程。盡管如此，對(duì)于單一地震事件，由已知地震動(dòng)譜求解震源譜的過程仍然十分復(fù)雜，需要消除地震波傳播路徑及臺(tái)站的影響才能直接得到震源譜參數(shù)(陳運(yùn)泰等，2003)。對(duì)于震群型地震事件，可將式(2)轉(zhuǎn)換為

式中，i=1，2，…，n，i為震群事件中第i個(gè)主要地震事件;n為震群中主要地震事件的總數(shù)。將式(3)展開并簡(jiǎn)化則有)

對(duì)于單一地震臺(tái)站，若其距離震群型地震的源震區(qū)相對(duì)較遠(yuǎn)，則可忽略震群中主要地震事件震中位置的相對(duì)距離，即可將較遠(yuǎn)的單一地震臺(tái)站相對(duì)于震群近似地認(rèn)為是同一點(diǎn)，由式(2)～(4)可得

對(duì)于距離震群相對(duì)較遠(yuǎn)的臺(tái)站，震群中地震的幾何擴(kuò)散函數(shù)(P)、非彈性衰減函數(shù)(F)、場(chǎng)地放大函數(shù)(R)和地震儀的響應(yīng)函數(shù)(I)近似相同(式(4))。即消除了地震波傳播路徑和臺(tái)站的影響因素之后，震群中主要地震事件的震源譜特征的研究可以通過分析臺(tái)站記錄到的地震動(dòng)譜參數(shù)特征來實(shí)現(xiàn)。

2 單臺(tái)地震動(dòng)譜參數(shù)

2013年前郭MS5.8震群發(fā)生后，CN2臺(tái)記錄到了該震群的地震動(dòng)觀測(cè)數(shù)據(jù)。CN2臺(tái)(125.26°E，43.48°N)為國(guó)家基準(zhǔn)臺(tái)、大震速報(bào)臺(tái)和國(guó)際資料交換臺(tái)，其觀測(cè)手段主要以測(cè)震綜合觀測(cè)為主，該臺(tái)地處吉林省中部丘陵與平原區(qū)的結(jié)合部位(圖1(b))①http://data.earthquake.cn/datashare/datashare_network_national_basic.jsp。CN2臺(tái)距離前郭震群震源區(qū)約160km，用式(5)研究前郭震群中5次主要地震事件的震源譜，則有

式中第1～5個(gè)地震事件依次表示前郭震群依據(jù)發(fā)震時(shí)間順序排列的5個(gè)主要地震事件(表1)。

圖2給出了CN2臺(tái)記錄到的前郭震群5次主要地震事件垂向(UD)、東西向(EW)和南北向(NS)三分量在時(shí)間域的地震波形。參照表1、2可以初步將前郭震群的5次主要地震事件分為3組:第1組為2013年10月31日先后發(fā)生的MS5.5、MS5.0地震，后文中分別簡(jiǎn)稱No.1、No.2地震事件;第2組為2013年11月22日發(fā)生的MS5.3地震，簡(jiǎn)稱No.3地震事件;第3組為2013年11月23日先后發(fā)生的MS5.8、MS5.0地震，分別簡(jiǎn)稱No.4、No.5地震事件，這里分組后給出的震群中主要地震事件編號(hào)與表1、式(6)中的一致。

對(duì)比CN2臺(tái)記錄到的前郭震群5次主要地震事件時(shí)間域的速度波形可以發(fā)現(xiàn)如下特點(diǎn):①時(shí)間域內(nèi)速度波形幅值最大的是No.4事件，這與其震級(jí)具有一定的一致性;②第1組中的No.2事件相對(duì)于No.1事件的衰減特征并不明顯，其中NS向分量相對(duì)偏大，這一點(diǎn)與二者的震級(jí)具有明顯的非一致性;③第2組中的No.3事件相對(duì)于第1組中的No.1、No.2事件的衰減特征也不明顯，且No.3、No.1事件的波形相似程度較高;④第3組中的No.5事件相對(duì)于No.4事件的衰減特征十分明顯，這一點(diǎn)明顯有別于第1組中的2次地震事件;⑤第1組中的No.2、第3組中的No.5兩次地震事件震級(jí)相當(dāng)，但No.5事件的波形、幅值相對(duì)于No.2事件明顯偏小。

圖2 前郭震群5次M S≥5.0地震事件速度波形三分量(垂向、東西、南北)實(shí)測(cè)記錄

前文通過對(duì)比圖2中給出的CN2臺(tái)記錄到的前郭震群5次主要地震事件的地震波形，明確指出了該震群中5次主要地震事件速度三分量與其震級(jí)存在較為明顯的非一致性特征，這種非一致性的特征表現(xiàn)在時(shí)間域中是由速度波形的形態(tài)、幅值、持續(xù)時(shí)間等呈現(xiàn)出來的。本研究進(jìn)一步對(duì)比了CN2臺(tái)觀測(cè)數(shù)據(jù)在頻率域的特征(圖3)。圖3中給出了CN2臺(tái)記錄到的前郭震群5次主要地震事件三分量實(shí)測(cè)記錄的傅里葉頻譜特征。第1組中No.2事件相較于No.1事件在大于1.0Hz的相對(duì)高頻范圍內(nèi)衰減特征不明顯，水平向2個(gè)分量(EW和NS)的高頻成份相對(duì)偏少，垂向(UD)分量與No.1事件相當(dāng);第2組中No.3事件的UD向分量與第1組中No.1事件的一致性程度較高，EW和NS向在小于1Hz的低頻范圍內(nèi)與No.1事件一致性程度較高，但在大于1Hz的相對(duì)高頻范圍內(nèi)則更接近于No.2事件的特征，尤其在0.5～1.0Hz的范圍內(nèi)，No.3、No.2事件均出現(xiàn)了1個(gè)相對(duì)平緩的區(qū)域，這一點(diǎn)與前郭震群中其它地震事件均有所不同;第3組中No.5事件相較于No.4事件的衰減特征在頻率域仍較為顯著，且相較于震級(jí)相當(dāng)?shù)腘o.2事件也明顯幅值偏低，這一點(diǎn)與二者在時(shí)間域的特征是趨于一致的。

前郭震群的5次主要地震事件隨時(shí)間在發(fā)震過程中的起伏與其震級(jí)存在的非一致性特征，能夠從圖2、3中的時(shí)間域和頻率域特征方面給出一定的描述。但前文中的描述多為定性描述而并未定量的展開，針對(duì)震群型地震事件，進(jìn)一步分析這種非一致性特征產(chǎn)生的原因，我們認(rèn)為可能是由5次主要地震事件各自的震源特征所致。本研究在式(1)～(6)的推導(dǎo)過程中已指出，對(duì)于距離震群相對(duì)較遠(yuǎn)的單臺(tái)，在消除了地震波傳播路徑(P和F)和臺(tái)站(R和I)的影響因素之后，可以通過研究單臺(tái)記錄到的地震動(dòng)譜參數(shù)來對(duì)比分析震群中主要地震事件的震源譜特征。以2013年前郭震群為例，進(jìn)一步可以基于CN2臺(tái)的地震動(dòng)觀測(cè)記錄對(duì)比其震源譜參數(shù)特征，進(jìn)而針對(duì)該震群中5次主要地震事件相對(duì)于震級(jí)的非一致性起伏特征給出基于震源特征的相關(guān)解釋。

圖3 前郭震群5次M S≥5.0地震事件速度波形三分量實(shí)測(cè)記錄的傅里葉頻譜特征

3 基于單臺(tái)地震動(dòng)譜的震源譜及遠(yuǎn)場(chǎng)地震波輻射能比較

對(duì)于單一地震事件而言，地震矩率函數(shù)在時(shí)間域中表示為標(biāo)量地震矩乘以震源時(shí)間函數(shù)，在頻率域中則表示為震源譜(Brune，1970、1971;陳運(yùn)泰，2003)

式中，f為頻率;f0為拐角頻率;γ為高頻衰減常數(shù)。通常把振幅譜高頻漸近趨勢(shì)(包絡(luò)線)和低頻趨勢(shì)(零頻水平)的交點(diǎn)稱作拐角，與拐角相對(duì)應(yīng)的頻率稱作拐角頻率f0。通常情況下，f0很難在頻譜圖上直接準(zhǔn)確測(cè)得，這里我們采用Brune模型(1970、1971)，即令γ=2。

當(dāng)震源距r很小、非彈性衰減的影響也很小且可以忽略時(shí)，Andrews(1986)提出可以通過積分的方式求解速度平方譜及位移平方譜等震源譜參數(shù)，即對(duì)于近震觀測(cè)的情況，強(qiáng)地面運(yùn)動(dòng)位移的頻譜可以表示為

式中，D(f)為位移的頻譜;D2(f)為位移的平方譜(D2(f)=V2(f)·(2πf)-2，V2(f)為速度的平方譜);f為頻率，則SD2為位移平方譜的積分。同理有

且有拐角頻率f0可通過式(8)、(9)的積分求得，即

式中，f0的計(jì)算方式可以應(yīng)用于任何頻譜，因此可將其作為拐角頻率的定義式。Andrews(1986)同時(shí)提出，對(duì)于點(diǎn)源或圓盤斷層模型，其地震波輻射能Erad、震源譜平方及地震矩M0可以由式(11)～(13)近似求得

式中，ρ為密度;β為S波的波速。值得注意的是，Andrews(1986)提出的基于地震動(dòng)譜求解震源譜參數(shù)方法的前提條件是震源距r相對(duì)較小，這一點(diǎn)與震源譜求解的歸一化條件也是一致的。

本研究基于單臺(tái)地震動(dòng)譜，在比較分析震群地震事件震源譜參數(shù)特征的過程中，進(jìn)行如下設(shè)定:假設(shè)單臺(tái)S觀測(cè)到的相對(duì)較遠(yuǎn)的某一震群O的n次地震事件來自于震源距r相對(duì)較小的某一震群Oapp的napp次地震事件(napp=n)。如此，一方面可以通過式(5)基于S臺(tái)的地震動(dòng)觀測(cè)譜實(shí)現(xiàn)震群O的n次地震的比較研究;另一方面，基于前文中式(8)～(13)，可以計(jì)算得到的震群Oapp“近源區(qū)”的近場(chǎng)地震波輻射能及地震矩等震源參數(shù)。其中相對(duì)于震群Oapp得到的可以近似地認(rèn)為是由相對(duì)較遠(yuǎn)的震群O傳播至S臺(tái)的遠(yuǎn)場(chǎng)地震波輻射能Efar-field，則有

進(jìn)一步結(jié)合視應(yīng)力σapp的概念(Wyss et al，1968)，定義單臺(tái)視應(yīng)力

以前郭震群為例，針對(duì)該震群中5次主要地震事件 No.1～No.5，結(jié)合式(6)、(11)和(14)，則有

式(16)實(shí)現(xiàn)了基于單臺(tái)速度的平方譜對(duì)比震群中不同地震事件遠(yuǎn)場(chǎng)地震波輻射能相對(duì)大小的研究目的。進(jìn)一步對(duì)于普遍的震群型地震事件，可將式(16)轉(zhuǎn)換為

式中，i=1，2，…，n，i為震群中第i個(gè)主要地震事件;n為震群中主要地震事件的總數(shù)。

4 計(jì)算結(jié)果

圖4是基于式(11)～(15)的遠(yuǎn)場(chǎng)地震波輻射能及單臺(tái)視應(yīng)力計(jì)算結(jié)果的對(duì)比圖，其參照了前郭震群5次主要地震事件震級(jí)的相對(duì)變化。其中，圖4(a)、4(b)分別為前郭震群傳播至CN2臺(tái)的遠(yuǎn)場(chǎng)地震波輻射能和單臺(tái)視應(yīng)力隨著5次主要事件的變化曲線，因比較研究的需要，進(jìn)行了歸一化處理;圖4(c)為CN2臺(tái)接收到的前郭震群5次主要地震事件的遠(yuǎn)場(chǎng)地震波輻射能Efar-field(式(15))，單位為J;圖4(d)為CN2臺(tái)的單臺(tái)視應(yīng)力(式(16))，單位為 MPa。

圖4 CN2臺(tái)前郭震群5次M S≥5.0地震事件遠(yuǎn)場(chǎng)地震波輻射能(a)、(c)及單臺(tái)視應(yīng)力的對(duì)比(b)、(d)

圖4(a)中前郭震群5次主要地震事件歸一化的遠(yuǎn)場(chǎng)地震波輻射能Efar-field與其震級(jí)的非一致性特征十分明顯，且相較于圖2、3中更為直觀。首先，No.4事件釋放的地震波輻射能為5次事件中最高的，這一點(diǎn)與其為前郭震群中震級(jí)最大的事件相一致;其次，No.5事件所釋放的地震波輻射能相對(duì)最低，這一點(diǎn)與其在時(shí)間順序上為整個(gè)震群中最后1個(gè)主要事件并不矛盾;最后，No.2、No.5事件的震級(jí)均為震群中最小，但No.2事件的遠(yuǎn)場(chǎng)地震波輻射能Efar-field為5次主要事件中第2高值，這一點(diǎn)為前郭震群發(fā)展過程中在時(shí)間域和頻率域表現(xiàn)出來的與震級(jí)非一致性特征的關(guān)鍵，即No.2事件雖然在震級(jí)上表現(xiàn)出了相對(duì)于No.1事件的衰減，但地震波輻射能量并未衰減，反而有所增加，某種程度上加強(qiáng)了前郭震群序列整體能量的釋放進(jìn)程;同時(shí)No.1事件所釋放的能量與No.3事件的相當(dāng)，略高于No.5事件，且No.1、No.3和No.5事件各自釋放的能量大體相當(dāng)于No.2事件的50%。

進(jìn)一步，圖4(c)、4(d)分別給出了前郭震群5次地震事件的遠(yuǎn)場(chǎng)地震波輻射能及單臺(tái)視應(yīng)力的定量化特征。圖4(c)中5次地震事件的遠(yuǎn)場(chǎng)地震波輻射能Efar-field為0.4×1011～11×1011J，圖4(d)中單臺(tái)視應(yīng)力為0.03～0.31MPa。參照地震波能量與震級(jí)的經(jīng)驗(yàn)關(guān)系lgES=1.5MS+4.8，可以計(jì)算出前郭震群的地震波輻射能為0.2×1013～3.2×1013J(Choy et al，1995;陳運(yùn)泰等，2003)。因此，本文基于CN2臺(tái)計(jì)算得到的遠(yuǎn)場(chǎng)地震波輻射能Efar-field相當(dāng)于整體地震波能量ES的約1.9% ～3.4%，即前郭震群5次地震事件所釋放的地震波能量經(jīng)過約160km的傳播至CN2臺(tái)時(shí)絕大部分已經(jīng)衰減，CN2臺(tái)記錄到的地震波觀測(cè)記錄，僅僅反映了該震群對(duì)應(yīng)地震事件地震波輻射能的很小一部分。

5 結(jié)論和討論

震群型地震的主要特點(diǎn)之一即為其序列能量是通過多次震級(jí)相近的地震釋放出來的，沒有突出的主震。最大地震在全序列中所占能量比例一般均小于80%，能量的釋放有明顯的起伏，衰減速度慢，活動(dòng)的持續(xù)時(shí)間長(zhǎng)，且震源通常較淺。由圖4(a)、4(c)可見，前郭震群5次主要地震事件的地震能量釋放的起伏特征明顯，且震級(jí)最大的No.4事件的地震波輻射能量大致相當(dāng)于其他4次主要地震事件釋放能量之和，滿足震群型地震的主要特征。

2013年吉林松原前郭震群的特點(diǎn)在于5次MS≥5.0主要地震事件短時(shí)間內(nèi)爆發(fā)式活動(dòng)，且前2次地震事件(No.1、No.2)又具有“主-余型”地震的假象特征，從震級(jí)發(fā)展水平來看(表2)，震群中震級(jí)最大的No.4地震事件發(fā)生前，No.3事件的發(fā)生為整個(gè)序列升級(jí)過程的轉(zhuǎn)折點(diǎn)，基于圖2、3的定性分析，似乎也可以支持No.3事件為轉(zhuǎn)折事件的推測(cè)。然而，通過前文進(jìn)一步計(jì)算5次主要地震事件的遠(yuǎn)場(chǎng)地震波輻射能和單臺(tái)視應(yīng)力發(fā)現(xiàn)，No.2事件才是前郭震群真正的轉(zhuǎn)折點(diǎn)。實(shí)際上，前郭震群中No.1為震群的起點(diǎn)事件，No.2為轉(zhuǎn)折事件，No.3為過渡事件，No.4是整個(gè)震群震級(jí)最大、地震波能量釋放最多的事件，而No.5則為前郭震群爆發(fā)式活動(dòng)趨于結(jié)束的事件。

此外，本研究主要是基于震群地震活動(dòng)震源區(qū)相對(duì)集中、單一接收單臺(tái)路徑一致的特點(diǎn)，再結(jié)合地震動(dòng)譜與震源譜在頻率域的對(duì)應(yīng)關(guān)系，推導(dǎo)出了式(6)，從而消除了地震波傳播路徑和臺(tái)站的影響，實(shí)現(xiàn)了基于單臺(tái)記錄到的地震動(dòng)譜參數(shù)特征對(duì)比分析震群型地震事件震源譜特征的研究目的。在該轉(zhuǎn)換過程中，始于震源譜計(jì)算的基本公式，針對(duì)震群型地震的特殊性，選擇單一臺(tái)站的地震動(dòng)觀測(cè)數(shù)據(jù)，通過時(shí)-頻域的轉(zhuǎn)換、“近源區(qū)”震群的設(shè)定，、遠(yuǎn)場(chǎng)地震波輻射能的近似計(jì)算及單臺(tái)視應(yīng)力的求取，實(shí)現(xiàn)了針對(duì)震群型地震震源譜參數(shù)進(jìn)行比較研究的目的，并對(duì)震群的起伏特征給出了基于震源特征的解釋。

由于本文主要是針對(duì)前郭震群展開的研究，觀測(cè)臺(tái)站數(shù)據(jù)也僅取自CN2臺(tái)，因此尚未實(shí)現(xiàn)針對(duì)多臺(tái)觀測(cè)數(shù)據(jù)計(jì)算結(jié)果的比較和對(duì)比分析。實(shí)際上，本文的研究結(jié)果主要是初步驗(yàn)證了式(6)、(17)的合理性和有效性，初步實(shí)現(xiàn)了基于地震動(dòng)譜比較震群中主要地震事件震源譜參數(shù)特征的研究目的。在未來的研究工作中，如何選取合理的單臺(tái)震源距、單臺(tái)個(gè)數(shù)及如何設(shè)定震群型地震主要地震事件的震級(jí)下限等都是值得進(jìn)一步深入的研究方向。

致謝:本文在成稿過程中，劉桂萍研究員、劉杰研究員、蔣海昆研究員、張永仙研究員和周龍泉研究員給予了有益的指導(dǎo)與幫助，作者謹(jǐn)表謝意;感謝吉林省地震局提供的2013年前郭震群的實(shí)測(cè)數(shù)據(jù)。