薛 艷， 陳琳榮， 黎明曉， 劉雙慶

(1．中國地震臺網(wǎng)中心，北京 100045；2．吉林省地震局工程地震研究所，長春 130022;3. 天津市地震局，天津 300201)

華北地區(qū)震群活動與強震

薛 艷1， 陳琳榮2， 黎明曉1， 劉雙慶3

(1．中國地震臺網(wǎng)中心，北京 100045；2．吉林省地震局工程地震研究所，長春 130022;3. 天津市地震局，天津 300201)

通過系統(tǒng)研究華北地區(qū)震群活動與強震的關(guān)系，發(fā)現(xiàn)：①震群活動與區(qū)域強震關(guān)系密切，小震群頻度的起伏可以反映區(qū)域地震活動水平的強弱；②華北地區(qū)ML≥2.0震群半年累積頻度N≥8可以作為區(qū)域MS≥5.7地震的預(yù)測指標(biāo)，預(yù)測時間為異常出現(xiàn)后的10個月內(nèi)，預(yù)測效能評分R=0.615，大于具有97.5%置信水平的R0(0.437)；③震群頻度指標(biāo)可以通過Molchan模型檢驗，對于不同強度的目標(biāo)地震(MS≥5.0、MS≥5.5與MS≥5.7地震)，其預(yù)測效能均好于Poisson模型，其中60%目標(biāo)地震的預(yù)測概率增益可以達到2～4水平；④區(qū)域震群對強震地點的預(yù)測意義有限：在張北6.2級、瑞昌5.7級與前郭5.8級地震前2年，與強震相關(guān)的構(gòu)造帶震群活躍；震前中短期階段(2～10個月)，震群向震中及周邊遷移；而對其余4次強震，震群對發(fā)震地點沒有指示意義。

華北地區(qū)；震群；強震；統(tǒng)計檢驗；預(yù)測指標(biāo)

0 引言

震群與強震的關(guān)系很早就引起了地震學(xué)家的關(guān)注。Kanamori[1]發(fā)現(xiàn)，在日本和歌山附近震群活動非常突出，且有感地震頻次與菲律賓海板塊北部邊界帶的8級以上地震密切相關(guān)，并用數(shù)值模擬方法予以解釋。在我國淺源大地震前若干年，在未來震中區(qū)及其周圍大范圍地區(qū)，震群活動顯著增多[2]。

有關(guān)震群和強震的關(guān)系研究大致可以分為2個方面，一是根據(jù)震群序列本身，分別從序列在時、空、強方面表現(xiàn)出來的某些特征，以及在序列發(fā)展的物理背景上(如波形變化、震源參數(shù)等)尋找表征前兆性震群的某些特征指標(biāo)[3-7]。近年來，隨著數(shù)字化波形資料的使用，學(xué)者們通過高精度定位、震源機制解等手段確定某個震群的發(fā)震斷層及破裂機制，探討震群的成因[8-10]。另一方面的研究涉及一個較大區(qū)域范圍內(nèi)震群時空演化的群體特征，尋找與強震孕育和發(fā)生的關(guān)系。姜秀娥等[11]研究了1967—1979年華北區(qū)域震群的時空活動特征與唐山地震的關(guān)系；薛艷等[12]研究了唐山地震前后華北區(qū)域震群及前兆震群的時空演化特征；敖雪明等[13]系統(tǒng)研究了天山地震帶7級左右強震前的區(qū)域震群活動特征。研究表明，在強震前，區(qū)域震群活動增強可能是一個共性特征。

本文作者在日常震情跟蹤工作中積累了1995年以來華北及附近地區(qū)(30°～43°N，108°～124°E)的震群資料，并一直應(yīng)用震群頻次的起伏變化預(yù)測該區(qū)中強以上地震1,2,3。本文將系統(tǒng)研究區(qū)域震群活動與中強以上地震的關(guān)系，同時提煉震群預(yù)測指標(biāo)，為震情跟蹤提供判定依據(jù)。

1 華北地區(qū)1995年以來地震目錄最小完整性震級

地震目錄是進行地震活動性研究的基礎(chǔ)資料，完整的地震目錄是地震活動性研究結(jié)果可靠性的保證。黃瑋瓊等[14]研究結(jié)果顯示，北京、天津、河北、山西和山東地區(qū)從1970年開始ML3.0地震基本完整，安徽中北部及附近地區(qū)1971年開始ML3.0地震基本完整;焦遠(yuǎn)碧等[15]研究結(jié)果顯示，1979年以后華北地區(qū)(32°～42°N，110°～122°E)的測震臺網(wǎng)可以控制ML≥2.5地震；魏富勝等[16]研究了1990—2005年中國大陸地區(qū)地震觀測的完整震級下限，結(jié)果顯示，京津唐地區(qū)ML1.0以上地震基本完備,華北北部(渤海除外)、鄂爾多斯北緣、東緣、南緣、遼寧中西部、安徽和江蘇南部及湘、贛、皖交界地區(qū)ML≥1.5地震基本完整,華北其他地區(qū)ML2.0以上地震基本完整。2007年完成的“十五”計劃，使得地震觀測技術(shù)有了重要進步，大大加強了地震觀測的靈活性和可靠性。

2 震群資料的處理與震群頻度隨時間的變化

2.1 震群的選取原則

本文確定震群的原則主要參照姜秀娥等[11]的約定，具體判定原則如下：

1)同一地點(半徑20 km范圍內(nèi))的一群地震中，在半個月內(nèi)至少有3個主要地震，滿足M1>M2>M3，且ΔM1,2≤1，ΔM2,3≤1(ΔM1,2=M1-M2，ΔM2,3=M2-M3)，但M1<6.0。

2)震群活動峰值日(連續(xù)24小時)的地震頻度N滿足以下條件之一：ML≥1.0，N≥7；ML≥1.5，N≥5；ML≥2.0，N≥4；ML≥2.5，N≥3(當(dāng)這3次地震有2次3級以上時，M3≥2.0即可)。

3)對于持續(xù)活動數(shù)日至數(shù)年的震群點，當(dāng)發(fā)生新的顯著地震且伴隨小震的爆發(fā)活動，同時這些地震又滿足1)、2)兩個條件者，可作為新的一次震群事件。

4)對于強震的余震活動區(qū)，則以正常的衰減曲線為基線，凡高于基線，且滿足1)、2)兩個條件者，算作一次震群事件。

依據(jù)以上原則，華北地區(qū)(30°～43°N，108°～124°E)1995—2015年間共發(fā)生188次2級以上震群，其中2.0～2.9級79次，3.0～3.9級87次，4.0～4.9級20次，5.0～5.9級2次。平均每年大約發(fā)生9次2級以上震群。這些震群主要分布在張渤地震帶、汾渭地震帶的中北段、河北平原帶南段及遼寧、山東半島、江蘇及安徽地區(qū)(圖1a)。1996—2005年，華北及鄰區(qū)(圖1b)共發(fā)生21次MS≥5.0地震(不含余震)，其中10次MS≥5.5地震，MS≥6.0地震3次，最大為1996年5月3日內(nèi)蒙古包頭6.4級地震。

圖1 華北及鄰區(qū)震群(a: 1995—2015年)與MS≥5.0地震空間分布(b：1996—2005年)

考慮大華北地區(qū)的地震監(jiān)測能力，認(rèn)為1995年以來，除了海域，其他地區(qū)2級以上震群大體完整。圖2a～2c分別是ML≥2.0、ML≥3.0和ML≥4.0震群年頻度圖，圖2d是華北及鄰區(qū)MS≥5.5地震活動M-t圖?？梢钥闯?，1995年以來2級以上震群年頻度介于4～16次之間，沒有出現(xiàn)早期震群年頻度明顯少于近年來的情況(圖2a)。這表明，1995年以來華北地區(qū)2.0級以上震群基本完備。

2.0級以上震群頻度與3.0級以上震群頻度起伏變化大體一致。在1995—1999年、2002—2005年和2012—2015年3個時段，震群年頻度較高。與此相對應(yīng)的是，在此3個階段，華北及鄰區(qū)均發(fā)生了5.5級以上地震。而4.0級以上震群由于數(shù)量少，年頻度起伏與5.5級以上地震活動的相關(guān)性不好(圖2)。

a ML≥2.0； b ML≥3.0； c ML≥4.0； d M-t圖圖2 區(qū)域震群年頻度及大華北及鄰區(qū)MS≥5.5地震M-t圖

2.2 區(qū)域震群頻度隨時間的變化與中強以上地震

圖3是華北地區(qū)(30°～43°N，108°～124°E)2級以上震群頻度曲線(窗長6個月，滑動步長1個月)，在其中標(biāo)注了華北及鄰區(qū)(圖1b)1996年以來發(fā)生的所有5.0級以上地震(沒有標(biāo)注黃海海域發(fā)生的小于6.0級的地震)。

由圖3可以看出，震群高頻度與5.0級以上地震具有較好的相關(guān)性。值得注意的是，2008—2011年震群頻度也比較高，這可能與2008年5月12日四川汶川8.1級地震與2011年3月11日日本MW9.1地震有關(guān)。

表1列出了震群頻次N≥8開始出現(xiàn)的時間及后續(xù)5.7級以上地震，這些地震包括了華北地塊及其邊緣在此期間發(fā)生的所有MS≥5.7地震(圖1b)。可以看出，1995年以來震群頻度共出現(xiàn)8組N≥8的高值，其中7組高值后華北及鄰區(qū)發(fā)生了5.7級以上地震，1組高值后沒有發(fā)生5.7級以上地震。強震發(fā)生在震群頻度高值出現(xiàn)后的3～13個月，其中86%的地震發(fā)生在高值后的10個月內(nèi)。

圖3 華北地區(qū)ML≥2.0震群頻度曲線(6個月窗長，1個月滑動)與5.0級以上地震

序號震群高頻度(N≥8)出現(xiàn)的時間高頻度后MS≥5.7地震高頻度出現(xiàn)的時間距離中強震發(fā)生的時間差/月11996-01-251996-05-03包頭MS6.41996-11-09南黃海MS6.13.29.521997-05-191998-01-10張北MS6.27.732002-07-222003-08-16巴林左旗MS5.912.842005-08-052005-11-26瑞昌MS5.73.752012-07-29無62013-04-252013-10-31前郭MS5.52013-11-23前郭MS5.8772014-03-212015-04-15阿拉善左旗MS5.81382015-01-152015-04-15阿拉善左旗MS5.83

如果將震群頻度N≥8作為異常，則共出現(xiàn)8組高值異常。如果以10個月為預(yù)報時間，則虛報3次，報準(zhǔn)地震6次，漏報1次，異常出現(xiàn)后發(fā)震的條件概率為5/8=0.625。如果以13個月為預(yù)測時間，則虛報1次，漏報0次。1996—2015年預(yù)測區(qū)域內(nèi)5.7級以上地震在10個月內(nèi)的自然概率(平均發(fā)震率，不考慮余震)為0.38。顯然，當(dāng)出現(xiàn)震群高值異常后發(fā)生5.7級以上地震的條件概率明顯高于自然概率。

為了評價震群高頻次的地震預(yù)測效能，本文利用許紹燮[17]提出的R值評分(公式1)來評估。

(1)

將異常出現(xiàn)后10個月作為預(yù)測時間，則預(yù)測總時間為61個月(9.5+7.7+3.7+7+3+3×10=61)，預(yù)報研究的總時間為 (2015-1995+1)×12=252個月，則根據(jù)公式(1)，預(yù)測效能評分R=6/7-61/252=0.615，具有97.5%置信水平的R0=0.437，顯然R>R0，即震群半年累積頻次N≥8可以作為華北地區(qū)5.7級以上地震的預(yù)測指標(biāo)。

2.3 應(yīng)用Molchan圖表法對震群頻度的地震預(yù)測效能進行統(tǒng)計檢驗

Molchan圖表法(MED)可以對地震前兆的預(yù)報效能進行評估，并得到概率增益曲線[18-19]。定義時空占有率為τ，漏報率為υ。當(dāng)前兆異常閾值由最高變化到最低值時，時空占有率τ從0變到1，漏報率υ則從1變化到0。圖4給出了當(dāng)前兆指標(biāo)函數(shù)作為報警閾值變化時τ隨υ的變化曲線。很容易看到，在完全隨機猜測的情況下，υ=1-τ，是一條連接(0,1)和(1,0)的直線(圖4中Gain=1的實線，即Poisson模型)。一個比完全隨機的Poisson分布好的預(yù)警指標(biāo)，υ-τ曲線應(yīng)該在Gain=1實線的下方；而在Gain=1實線的上方的υ-τ曲線的預(yù)警指標(biāo)的預(yù)測效能比Poisson模型差(圖4)。

在統(tǒng)計學(xué)中，貝葉斯公式可表示為

P(E|A)=P(A|E)×P(E)/P(A) 。

(2)

式中：E表示地震；A表示前兆異常；P(E|A)表示有異常前兆下的地震發(fā)生條件概率；P(A|E)表示地震發(fā)生前，前兆異常出現(xiàn)的概率；P(E)為地震發(fā)生概率；P(A)為前兆異常出現(xiàn)的概率。概率增益可定義為：

Gain=P(E|A)/P(E) =P(A|E)/P(A) 。

(3)

根據(jù)式(3)可以得到概率增益Gain=1、Gain=2、Gain=3、Gain=4的直線，通過Molchan圖表計算可以得到如圖4所示關(guān)系圖。

圖4 Molchan圖表模板

為了評估震群頻度的預(yù)測效能，首先構(gòu)建震群頻度函數(shù)，確定6個月窗長，1個月滑動，對震群頻度進行采樣，這樣得到一個頻度序列，1個月1個值。

首先按照固定時間間隔劃定N個單元，依次為每個單元賦值。如果第i個單元的異常值，為第i+1個單元地震的前兆，則可以對每個地震都可以得到一個特征值。以該特征值為閾值，分別計算時空占有率和漏報率；依次下調(diào)閾值，計算時空占有率和漏報率，直到漏報率為零，并將結(jié)果投影到Molchan圖表。取1個月、2個月、3個月等不同的時間間隔來劃分單元，取均值作為單元的異常值。通過反復(fù)實驗，當(dāng)單元時間長度為3個月時預(yù)測效果最好。圖5a～5c分別是MS≥5.0、MS≥5.5與MS≥5.7地震的檢驗結(jié)果，可以看出，對于不同強度的目標(biāo)地震，利用震群頻度進行預(yù)測，其預(yù)測效能均好于Poisson模型。在18次5.0級以上地震中，有11次地震的預(yù)測概率增益為2～4；在10次5.5級以上地震中，有6次地震的預(yù)測概率增益可以達到2～4；在7次5.7級以上地震中，有4次地震的預(yù)測概率增益達到3～4。

圖5 震群頻度的預(yù)測檢驗結(jié)果(a：MS≥5.0；b：MS≥5.5；c：MS≥5.7；單元時間長度均為3個月)

2.4 7次MS≥5.7地震前2年區(qū)域震群時空演化特征

圖6給出了1996年5月3日包頭6.4級、1996年11月9日長江口6.1級、1998年1月10日張北6.2級、2003年8月16日巴林左5.9級、2005年11月26日瑞昌5.7級與2013年10月31日前郭5.8級地震前區(qū)域震群的空間分布。由圖可見，1995年1月至1996年包頭地震前，區(qū)域震群比較活躍，主要分布在華北北部和蘇皖浙及附近地區(qū)，其中4.0級以上震群1次，為1995年7月20日河北沙城4.1級震群(圖6a中紅色五角星)，該震群一直持續(xù)活動到1997年7月，特別是開始活動后至1995年底，共發(fā)生1.0級以上地震153次，其中2.0級以上地震48次，表現(xiàn)出很強的爆發(fā)性，該震群距離包頭地震震中約500 km。包頭6.4級地震后至南黃海6.1級地震前震群集中在張渤帶的東段—渤海附近地區(qū)，而長江口南黃海6.1級地震附近沒有發(fā)生震群(圖6b)。

1996年震群主要沿張渤帶的中東段分布，其中3.0級以上震群有3月28日遼寧新金3.0級、6月17日渤海3.4級、10月15日河北唐山3.6級及12月16日北京順義4.5級震群(圖6c中藍色圓)。1997年震群主要分布在2個區(qū)域，一個是長江口及其附近，發(fā)生了8月3日長江口4.4級和10月20日江蘇常熟3.2級地震，這2個震群發(fā)生在1996年11月9日長江口MS6.1地震后，可能是一次強震發(fā)生后應(yīng)力場調(diào)整的結(jié)果[20]；另一個震群相對集中的區(qū)域為張渤帶的西段，其中3.0級以上震群有1月10日山西大同3.3級、11月2日內(nèi)蒙古涼城3.3級和12月9日山西渾源3.6級震群(圖6c中紅色圓)。顯然，張北地震前2年區(qū)域震群主要沿張渤帶分布，特別是震前2個月華北地區(qū)僅發(fā)生3次2.0級以上小震群，均集中在張北地震震中附近(圖6d)。在這些震群中尤以北京順義4.5級震群最為突出，該震群從1996年12月16日開始一直持續(xù)至1998年3月初，共發(fā)生1.0級以上地震158次。

2002年區(qū)域震群數(shù)量顯著增多，統(tǒng)計區(qū)內(nèi)共發(fā)生8次3.0級以上震群，其中4.0級以上震群5次，分別為6月17日海城4.0級、7月14日霍山4.0級、8月3日懷來4.4級、10月12日大同4.0級和12月25日邢臺4.2級震群(圖6e藍色圓)。2003年1月至8月16日巴林左旗5.9級地震前震群數(shù)量有所減少，共發(fā)生3.0級以上震群4次，分別為1月19日邢臺3.9級、3月28日大同3.1級、5月9日安徽來安與江蘇六合交界3.1級及6月5日山東青島嶗山4.2級震群(圖6e紅色圓)。突出的震群是2003年青島嶗山4.2級震群，該震群從2003年5月21日開始活動一直持續(xù)到11月，發(fā)生30次1.0級以上地震，其中3.0級以上7次，4.0級以上地震2次。2004年11月1日該區(qū)再次發(fā)生4.1級地震，余震非常豐富。此外，值得注意的是，在統(tǒng)計區(qū)外圍，距離巴林左旗5.9級地震震中北西約220 km處，2002年10月20日發(fā)生了西烏珠穆沁旗5.2級震群。

2004年統(tǒng)計區(qū)內(nèi)震群數(shù)量不多，2.0級以上震群主要分布在遼寧、唐山及天津地區(qū)，其中3.0級以上震群2次，分別為3月7日海城3.5級和3月13日唐山3.9級震群(圖7f藍色圓)。2005年1月至瑞昌地震前震群數(shù)量顯著增多，特別是華北南部地區(qū)，震群數(shù)量增多明顯，發(fā)生了5月20日霍山3.2級、7月10日安徽金寨3.3級和10月2日山東聊城3.4級震群。其中金寨和霍山震群距離瑞昌5.7級地震震中約200 km(圖6f紅色圓)。

2012年大華北地區(qū)震群活動顯著增強，共發(fā)生3級以上震群10次，最大為2012年7月20日江蘇寶應(yīng)5.1級震群。這些震群主要分布在郯廬斷裂帶東側(cè)，即遼寧南部至江蘇地區(qū)(圖6g藍色圓)。2013年震群持續(xù)活躍，至2013年10月31日前郭地震前共發(fā)生3.0級以上地震11次，其中9次集中分布在遼蒙交界至山東半島，最大為2013年1月21日內(nèi)蒙古科爾沁左翼后旗4.3級震群(圖6g紅色圓)。前郭地震發(fā)生在震群集中區(qū)的東北部，距離最近的科爾沁4.3級震群約200 km(圖6g)。

2012年以來華北地區(qū)最突出的是遼寧蓋州震群與山東乳山震群，其中蓋州震群始于2012年2月，至2016年5月12日仍在持續(xù)，共發(fā)生1級以上地震1 587次，其中4級地震8次，最大震級ML4.7。乳山震群始于2013年10月，至2016年5月12日仍在活動，共發(fā)生ML≥1.0地震1 596次，其中4.0級以上地震4次。前郭地震發(fā)生在這2個震群活動持續(xù)當(dāng)中。

2014年至2015年4月15日阿拉善左旗5.8級地震前，區(qū)域震群仍非?；钴S，并于2015年1月15日出現(xiàn)高值異常，2015年1月18日在距離5.8級地震約170 km處的內(nèi)蒙古臨河發(fā)生了2.8級小震群，3月8日臨河再次發(fā)生3.8級地震。從空間分布看，震群主要分布在華北東部，而鄂爾多斯西緣很少(圖6h)。

以上分析表明，區(qū)域震群對強震地點的預(yù)測意義有限。在張北6.2級、瑞昌5.7級與前郭5.8級地震前2年，與主震相關(guān)的構(gòu)造帶震群活躍；震前中短期階段(2—10個月)震群向震中及周邊收縮。而在其余4次強震前震群對發(fā)震地點沒有指示意義。

a:包頭6.4級地震前(藍色與紅色圓分別表示1995年1月—1995年12月與1996年1月1日—1996年5月2日期間發(fā)生的震群)；b：包頭地震后至南黃海地震前(1996年5月3日—1996年11月8日)；c: 張北地震前 (藍色與紅色圓分別表示1996年1月—1996年12月與1997年1月—1997年12月期間發(fā)生的震群); d: 張北地震 (1997年11月1日—1998月1月9日); e: 巴林左地震前(藍色與紅色圓分別表示2002年1月—2002年12月與2003年1月1日—2003年8月23日期間發(fā)生的震群); f: 江西瑞昌地震前(藍色與紅色圓分別表示2004年1月—2004年12月與2005年1月1日—2005年11月25日期間發(fā)生的震群); g: 前郭地震前(藍色與紅色圓分別表示2012年1月—2012年12月與2013年1月1日—2013年10月30日期間發(fā)生的震群)；h：阿拉善左旗地震前(藍色圓與紅色圓分別表示2014年1月—2014年12月與2015年1月1日—2015年4月14日期間發(fā)生的震群)圖6 7次強震前2年左右區(qū)域震群空間分布

3 結(jié)論與討論

綜上所述，得到以下幾點結(jié)論：

1)華北地區(qū)震群活動與區(qū)域強震關(guān)系密切，小震群頻度的起伏可以反映區(qū)域地震活動水平的強弱。

2)華北地區(qū)ML≥2.0震群半年累積頻度N≥8可以作為區(qū)域5.7級以上地震的預(yù)測指標(biāo)，R=0.615，大于具有97.5%置信水平的R0=0.437。預(yù)測時間為異常出現(xiàn)后的10個月內(nèi)。

3)取半年震群累積頻度，滑動步長1個月為原始數(shù)據(jù)序列，應(yīng)用Molchan模型對震群頻度的檢驗。結(jié)果顯示，對于不同強度的目標(biāo)地震(MS≥5.0、MS≥5.5與MS≥5.7地震)，其預(yù)測效果均好于完全隨機狀態(tài)下的Possion分布。具體看，在18次5.0級以上目標(biāo)地震中有11次地震的預(yù)測概率增益為2～4；在10次5.5級以上目標(biāo)地震中有6次地震的預(yù)測概率增益可以達到2～4；在7次5.7級以上目標(biāo)地震中有4次地震的預(yù)測概率增益達到3～4。

4)1998年張北6.2級地震前2年，區(qū)域震群主要沿張渤地震帶分布，震前2個月震群向震中附近遷移；2005年瑞昌5.7級地震前11個月華北南部震群活躍；2013年前郭5.8級地震前2年，震群主要沿郯廬地震帶附近分布，震前10個月震群主要分布在郯廬帶的渤海至遼寧段附近。而1996年包頭6.4級、南黃海6.1級、2003年巴林左5.9級及2015年阿拉善左旗5.8級地震前2年，震群分布在整個華北地區(qū)，與主震地點的關(guān)系不明顯。

如何理解華北地區(qū)強震前大范圍的小震群活動增強呢？華北地區(qū)的小震群主要發(fā)生在強震破裂區(qū)、水庫庫區(qū)、斷裂端部與交匯部位等，多為弱化區(qū)，這些特殊區(qū)域?qū)?yīng)力場擾動比較敏感，容易形成應(yīng)力集中點，并以頻繁的小震活動來釋放應(yīng)變能。當(dāng)華北地區(qū)孕育強震(或一組強震)時，大區(qū)域應(yīng)力場擾動增強，多個弱化區(qū)應(yīng)力集中，從而形成多點震群活動。大范圍的小震群活動是強震(一組強震)的一種“場兆”。

[1] Kanamori H. Relation between tectonic stress, great earthquakes and earthquake swarms[J]. Tectonophysics, 1972, 14(1): 1-12.

[2] 梅世蓉, 薛艷, 宋治平. 華北地區(qū)強震前地震活動長期演變過程的共性[C]//國家地震局預(yù)測預(yù)防司. 地震短臨預(yù)報的理論與方法——“八五”攻關(guān)三級課題論文集. 北京: 地震出版社, 1997: 3-12.

[3] 宋俊高, 陸遠(yuǎn)忠, 施玉芳, 等. 震群用于地震預(yù)報的實用程式研究[C]//國家地震局科技監(jiān)測司. 地震預(yù)報方法實用化研究文集-地震學(xué)專集. 北京: 學(xué)術(shù)書刊出版社, 1989: 204-228.

[4] 王煒, 楊德志. 利用Weibull分布研究華北地區(qū)前兆震群的特征[J]. 中國地震, 1987, 3(4): 13-21.

[5] 陸遠(yuǎn)忠, 宋俊高, 戴維樂. 一個判斷震情的指標(biāo)—震群U值[J]. 地震學(xué)報, 1984, 6(S): 495-508.

[6] 顧浩鼎, 曹天青. 前兆震群和S波偏振[J]. 地震學(xué)報, 1980, 2(4): 343-355.

[7] 崔子健, 李志雄, 陳章立, 等. 判別小震群序列類型的新方法研究—譜振幅相關(guān)分析法[J]. 地球物理學(xué)報, 2012, 55(5):1718-1724, doi: 10.6038/j.issn.0001-5733.2012.05.028.

[8] 吳建平, 明躍紅, 張恒榮, 等. 長白山天池火山區(qū)的震群活動研究[J]. 地球物理學(xué)報, 2007, 50(4): 1089-1096.

[9] 曹鳳娟, 王亮, 焦明若, 等. 2012年2月2日蓋州ML4.7震群序列特征及發(fā)震構(gòu)造探討[J]. 地震地質(zhì), 2013, 35(4): 842-852.

[10] 鄭建常, 曲利, 曲均浩, 等. 乳山震群大小臺網(wǎng)地震定位的對比研究[J]. 地震地質(zhì), 2015, 37(4): 953-965.

[11] 姜秀娥, 陳非比. 區(qū)域震群與唐山大震[J]. 地震學(xué)報, 1983, 5(2): 145-157.

[12] 薛艷, 梅世蓉. 唐山大震前后一般震群與前兆震群時空演化特征研究[J]. 地震, 1999, 19(2): 183-189.

[13] 敖雪明, 龍海英, 王桂嶺. 新疆天山近7級地震前震群活動的時空分布演化特征[J]. 地震, 2010, 30(4): 124-132.

[14] 黃瑋瓊, 李文香, 曹學(xué)鋒. 中國大陸地震資料完整性研究之二——分區(qū)地震資料基本完整的起始年分布圖像[J]. 地震學(xué)報, 1994, 16(4): 423-432.

[15] 焦遠(yuǎn)碧, 吳開統(tǒng), 楊滿棟. 我國地震臺網(wǎng)監(jiān)測能力及臺網(wǎng)觀測條件質(zhì)量評定[J]. 中國地震, 1990, 6(4): 1-7.

[16] 魏富勝, 許忠淮, 王淑輝, 等. 中國大陸地區(qū)的地震檢測能力[J]. 地震地磁觀測與研究, 2009, 30(5): 26-33.

[17] 許紹燮. 地震預(yù)報能力評分[C]//國家地震局科技監(jiān)測司. 地震預(yù)報方法實用化研究文集—地震學(xué)專輯. 北京: 學(xué)術(shù)書刊出版社, 1989: 586-590.

[18] Chen S, Jiang C S, Zhuang J C. Statistical evaluation of efficiency and possibility of earthquake predictions with gravity field variation and its analytic signal in western China[J]. Pure and Applied Geophysics, 2016, 173(1): 305-319, doi: 10.1007/s00024-015-1114-x.

[19] 陳石, 徐偉民, 蔣長勝. 中國大陸西部重力場變化與強震危險性關(guān)系[J]. 地震學(xué)報, 2015, 37(4): 562-574, doi: 10.11939/jass.2015.04.005.

[20] 李興才，馮樹文．震群活動的一種可能的觸發(fā)機制[J]. 地震學(xué)報，1983，5(3)：297-303．

Relationship between Earthquake Swarms and the Strong Earthquake in North China

XUE Yan1, CHEN Lin-rong2, LI Ming-xiao1, LIU Shuang-qing3

(1.China Earthquake Networks Center, Beijing100045, China; 2. Institute of Engineering Mechanics, Earthquake Administration of Jilin Province, 130022, China; 3. Earthquake Administration of Tianjin Municipality, Tianjin 300201, China)

The relationship between earthquake swarms and the strong earthquake in North China has been studied, the results show: ①Earthquake swarm activity is closely related to regional strong earthquakes, and the frequency of small swarms can reflect the level of regional seismic activity. ②The half year cumulative frequency ofML≥2.0swarms great than or equal to 8 can be used as a predictive index ofMS≥5.7earthquakes. The prediction time is 10 months after the appearance of the anomaly. The value of the predictive efficacy score R is 0.615, greater than that of R0, 0.437, with a 97.5% confidence level. ③Frequency index of earthquake swarms can be tested by Molchan model. For target earthquakes withMS≥5.0,MS≥5.5 andMS≥5.7 respectively, the predicted results are all better than the Poisson model. The prediction probability gain can be achieved 2～4 for the 60% of the target earthquakes especially. ④The place prediction of earthquake swarms to the strong earthquake is limited. Two years before the ZhangbeiMS6.2, RuichangMS5.7 and QianguoMS5.8 earthquake, swarms are mainly active in the vicinity of the tectonic zones associated with the strong earthquake. While they migrate to the epicenter and the surrounding area 2 months to 10 months before the strong earthquake. But for the rest strong earthquakes, the place prediction of swarms is not significant.

North China; earthquake swarm; strong earthquake; statistical test; prediction index

2017-02-08

地震科技星火計劃項目“全球巨大地震前地震活動及序列演化特征研究”(XH17048)；震情跟蹤項目“華北地區(qū)6級以上地震危險性研究”(2017020504)

薛艷(1969—)，女，北京人，高級工程師，現(xiàn)主要從事地震預(yù)測研究．E-mail：xueyan5619@seis.ac.cn

P315．2

A

1003-1375(2017)02-0001-09

10.3969/j.issn.1003-1375.2017.02.001

薛艷，陳琳榮，黎明曉，等. 華北地區(qū)震群活動與強震[J]. 華北地震科學(xué)，2017，35(2):1-9.

1 中國地震臺網(wǎng)中心. 2012. 2013年度中國地震趨勢研究報告. 2 中國地震臺網(wǎng)中心. 2013. 2014年度中國地震趨勢研究報告. 3 中國地震臺網(wǎng)中心. 2014. 2015年度中國地震趨勢研究報告.