焦明若，夏彩韻，王喜龍，王 亮，孔祥瑞，張志宏

（遼寧省地震局，遼寧 沈陽(yáng) 110034）

遼寧地區(qū)中強(qiáng)震預(yù)測(cè)指標(biāo)體系研究

焦明若，夏彩韻，王喜龍，王 亮，孔祥瑞，張志宏

（遼寧省地震局，遼寧 沈陽(yáng) 110034）

分別提取地震三要素的異常指標(biāo)，建立了遼寧地區(qū)異常指標(biāo)體系。通過(guò)梳理遼寧及周邊中強(qiáng)地震前測(cè)震學(xué)科異常變化特征，提取異常指標(biāo)，研究各指標(biāo)與地震三要素的關(guān)系，給出異常指標(biāo)的權(quán)重；根據(jù)歷史震例研究前兆學(xué)科典型趨勢(shì)和短臨異常變化特征，提取異常變化幅度、持續(xù)時(shí)間和形態(tài)變化特征等定量指標(biāo)，同時(shí)給出各異常指標(biāo)與地震三要素的關(guān)系，依據(jù)地震對(duì)應(yīng)情況給出各異常指標(biāo)的權(quán)重；最后結(jié)合遼寧地區(qū)測(cè)震學(xué)前及前兆學(xué)指標(biāo)，以及各異常指標(biāo)權(quán)重信息集成，建立遼寧地區(qū)中強(qiáng)震預(yù)測(cè)指標(biāo)體系，為遼寧地區(qū)危險(xiǎn)區(qū)判定和跟蹤提供依據(jù)與技術(shù)支持。

遼寧地區(qū)；異常提?。划惓Ｖ笜?biāo)；信息集成

0 引言

自2008年以來(lái)，遼寧地區(qū)多次發(fā)生4級(jí)地震或震群活動(dòng)，特別是2013年以來(lái)遼寧地區(qū)相繼發(fā)生了燈塔M5.1和科爾沁M(jìn)5.3地震，隨后，遼寧地區(qū)又相繼發(fā)生一系列5級(jí)地震，形成了一條近南北向的5級(jí)地震條帶，該地區(qū)應(yīng)力水平明顯增強(qiáng)。遼寧地區(qū)一直是華北地區(qū)打破6級(jí)地震平靜的首發(fā)地區(qū)之一，多年來(lái)也一直是華北地區(qū)劃定的地震重點(diǎn)危險(xiǎn)區(qū)之一。因此未來(lái)一段時(shí)間遼寧地區(qū)6級(jí)左右地震的危險(xiǎn)性進(jìn)一步增強(qiáng)，震情跟蹤的任務(wù)更加艱巨?？紤]到遼寧地區(qū)連續(xù)多年被劃為6級(jí)地震重點(diǎn)危險(xiǎn)區(qū)以及當(dāng)前嚴(yán)峻的震情形勢(shì)，建立遼寧地區(qū)中強(qiáng)震的前兆指標(biāo)體系并指導(dǎo)未來(lái)地震預(yù)報(bào)實(shí)踐非常緊迫且必要。

考慮到遼寧地區(qū)中強(qiáng)地震活動(dòng)的規(guī)律性，本文采用遼寧地區(qū)1970年以來(lái)發(fā)生的10次M5.0以上地震進(jìn)行總結(jié)。遼寧地區(qū)指標(biāo)體系的建立可以通過(guò)長(zhǎng)中短臨的思路建立，盡量考慮到不同構(gòu)造分區(qū)存在指標(biāo)不同的特點(diǎn)。文中分別給出遼寧地區(qū)M5.0以上地震的不同學(xué)科異常特征和三要素預(yù)測(cè)指標(biāo)。最后對(duì)測(cè)震學(xué)和部分前兆學(xué)異常指標(biāo)進(jìn)行了信息集成。

1 選取資料和分析方法

1.1震例選取

對(duì)于遼寧地區(qū)中強(qiáng)震前異常資料的收集與分析主要基于中國(guó)震例[1-2]（1966—2006年）與前人公開(kāi)發(fā)表的相關(guān)文獻(xiàn)、著作[3-5]以及年度趨勢(shì)報(bào)告①遼寧省地震局. 遼寧省地震趨勢(shì)年度會(huì)商報(bào)告.沈陽(yáng)：遼寧省地震局1990-2016.。選取的歷史中強(qiáng)地震為1970年以來(lái)的M5.0以上地震，對(duì)所收集到的遼寧及鄰區(qū)發(fā)生的歷史中強(qiáng)地震震例進(jìn)行分析，挑選出有遼寧地區(qū)前兆異常對(duì)應(yīng)的地震（表1），而后依據(jù)長(zhǎng)、中、短臨時(shí)間尺度對(duì)所挑選的地震進(jìn)行前兆異常特征提取及分析研究。

表1 遼寧及鄰區(qū)1970年以來(lái)中強(qiáng)地震震例

1.2測(cè)震學(xué)指標(biāo)提取方法

測(cè)震學(xué)前兆異常提取的技術(shù)途徑主要采用日常會(huì)商中比較成熟的異常指標(biāo)提取方法?？煞譃榈卣鸹顒?dòng)圖像異常與地震波異常，主要提取方法包括地震空區(qū)、條帶；前兆震群、前震的識(shí)別；相關(guān)地震、余震窗、信號(hào)震、逼近地震，平靜或活躍的判別及測(cè)震學(xué)的各項(xiàng)參數(shù)（如b值、Q值、D值（危險(xiǎn)度）、C值（集中度）、U值（能量均勻度）、K值（歸一化熵值）、Rm值（小震調(diào)制比）、波速比、P波S波振幅比、應(yīng)力降、視應(yīng)力等等的異常提取。

1.3前兆學(xué)指標(biāo)提取方法

前兆異常提取的技術(shù)途徑主要采用日常會(huì)商中比較成熟的異常指標(biāo)提取方法進(jìn)行。流體學(xué)科一般采用差分、相關(guān)系數(shù)、小波分析、平滑濾波以及閾值等分析方法；形變學(xué)科一般采取平滑濾波、多項(xiàng)式擬合、卡爾曼濾波法、速度分析、周期濾波等方法；電磁學(xué)科一般采用滑動(dòng)傅氏去年變周期、濾波（滑動(dòng)、數(shù)字濾波）波譜分析以及加卸載響應(yīng)比等分析方法。

1.4地震前兆異常特征

1.4.1 M≥7地震前兆異常特征

通過(guò)對(duì)遼寧及鄰區(qū)歷史中強(qiáng)地震前兆異常特征進(jìn)行分析，發(fā)現(xiàn)M≥7地震具有如下特征：遼寧地區(qū)7級(jí)地震前兆異常震中距范圍大，可達(dá)500km以上（圖1a）；異常持續(xù)時(shí)間長(zhǎng)，最長(zhǎng)可達(dá)5年以上（圖1b）；異常數(shù)量多，定點(diǎn)前兆異常大于40項(xiàng)；異常變化幅度大，最大可達(dá)正常值幾十倍以上；異常具同步性變化特征；異常發(fā)震特征流體主要以突跳、上升發(fā)震為主，下降發(fā)震和異常恢復(fù)后發(fā)震次之，形變主要以趨勢(shì)轉(zhuǎn)向和破年變?yōu)橹?，破年變又體現(xiàn)為短期加速、年變變幅增大或較小等動(dòng)態(tài)特征。地震一般多發(fā)生在異常頻次開(kāi)始轉(zhuǎn)折降低階段，即距最高頻次1個(gè)月內(nèi)發(fā)生地震（圖2a）。宏觀異常分布范圍廣，可達(dá)上百至近千公里；異常數(shù)量巨大，可達(dá)千起（圖3）；異常早期少，臨震階段劇增，且隨時(shí)間推移由外圍向震中區(qū)遷移。

1.4.2 M 5～6地震前兆異常特征

通過(guò)對(duì)遼寧及鄰區(qū)歷史中強(qiáng)地震前兆異常特征進(jìn)行分析，發(fā)現(xiàn)M5～6地震具有如下特征：遼南地區(qū)M5～6地震流體前兆異常多集中在距震中200km范圍內(nèi)，遼西地區(qū)前兆異常震中距可達(dá)350km左右（圖1a） ，異常主要以水氡、傾斜、地磁異常為主，持續(xù)時(shí)間較長(zhǎng)，多為趨勢(shì)性異常，短臨異常具準(zhǔn)同步性特征（圖1b）；異常數(shù)量較M7明顯減少，M5～6地震前定點(diǎn)前兆異常數(shù)量?jī)H有十幾項(xiàng)左右；異常發(fā)震特征與7級(jí)地震類似，主要以突跳、上升和破年變發(fā)震為主，下降發(fā)震和異?；謴?fù)后發(fā)震次之，地震一般在異常最高頻詞開(kāi)始轉(zhuǎn)折降低后1—2個(gè)月內(nèi)發(fā)生（圖2b）。宏觀異常相對(duì)較少（20～30項(xiàng)）或無(wú)宏觀異常，空間分布范圍也相對(duì)較小。

2 三要素指標(biāo)提取及信息集成

2.1測(cè)震學(xué)指標(biāo)

通過(guò)對(duì)以上10次地震進(jìn)行震前指標(biāo)的提取，歸納出遼寧地區(qū)中強(qiáng)地震前存在的測(cè)震學(xué)異常及對(duì)應(yīng)地震情況，如表2所示。

根據(jù)遼寧地區(qū)測(cè)震異常指標(biāo)對(duì)地震三要素的不同指示意義，可將以上測(cè)震學(xué)指標(biāo)分為時(shí)間指標(biāo)（表3）、強(qiáng)度指標(biāo)（表4）、地點(diǎn)指標(biāo)（表5）三類。

表3 遼寧地區(qū)地震活動(dòng)預(yù)測(cè)時(shí)間指標(biāo)

表4 遼寧地區(qū)地震活動(dòng)預(yù)測(cè)地點(diǎn)指標(biāo)

表5 遼寧地區(qū)地震活動(dòng)預(yù)測(cè)強(qiáng)度指標(biāo)

2.2前兆學(xué)三要素指標(biāo)

通過(guò)對(duì)以上10次地震進(jìn)行震前指標(biāo)的提取，歸納出遼寧地區(qū)中強(qiáng)地震前存在的前兆學(xué)科異常及對(duì)應(yīng)地震情況，具體見(jiàn)表6—8。

2.3指標(biāo)信息集成

（1）“權(quán)重集成”介紹

前述的指標(biāo)中部分指標(biāo)在震前有的出現(xiàn)異常，有的沒(méi)有出現(xiàn)異常，同時(shí)有的指標(biāo)出現(xiàn)異常并未發(fā)震。我們要判斷未來(lái)遼寧地區(qū)發(fā)生ML5.0地震的危險(xiǎn)性時(shí)，將異常指標(biāo)用權(quán)重集成的方法進(jìn)行綜合，當(dāng)權(quán)重集成值P達(dá)到某一閾值后，就有地震危險(xiǎn)。

為了更好的衡量每個(gè)異常指標(biāo)的預(yù)測(cè)效能，我們計(jì)算了每個(gè)指標(biāo)的報(bào)準(zhǔn)率、虛報(bào)率及漏報(bào)率具體見(jiàn)表9。其中報(bào)準(zhǔn)率為預(yù)報(bào)有震次數(shù)與實(shí)際發(fā)震次數(shù)的比值、虛報(bào)率為出現(xiàn)異常并未發(fā)生地震的次數(shù)與出現(xiàn)異?？偞螖?shù)的比值、漏報(bào)率為未出現(xiàn)異常并發(fā)生地震的次數(shù)與統(tǒng)計(jì)總發(fā)震次數(shù)的比值。

該方法是對(duì)一組預(yù)報(bào)指標(biāo)yi，令當(dāng)yi出現(xiàn)異常為1，未出現(xiàn)異常為0， 以權(quán)重iw的大小表征該項(xiàng)指標(biāo)在過(guò)去地震預(yù)報(bào)中的優(yōu)劣，建立權(quán)重方程。權(quán)重集成值P由下式求出[7]：

因此各項(xiàng)指標(biāo)的權(quán)重iw，由下式求出：

式中，Ei表示第i個(gè)預(yù)報(bào)指標(biāo)過(guò)去在預(yù)報(bào)中的錯(cuò)報(bào)率。本文中Ei用地震虛報(bào)率和地震漏報(bào)率之和表示。為更好地提高該方法的準(zhǔn)確性，從表9地震異常參數(shù)中挑選出16項(xiàng)預(yù)報(bào)效能較好的指標(biāo)，將已選用的16項(xiàng)指標(biāo)列于表10中，同時(shí)根據(jù)上述公式算出每項(xiàng)指標(biāo)的權(quán)重wi，并以年為單位求出的Pi值及與地震對(duì)應(yīng)情況見(jiàn)表11和圖4。

（2）預(yù)報(bào)未來(lái)一年有無(wú)ML≥5.0地震的綜合預(yù)報(bào)指標(biāo)—權(quán)重集成值Pi

為判斷未來(lái)一年遼寧地區(qū)ML≥5.0地震危險(xiǎn)性，選用16項(xiàng)有年度預(yù)報(bào)能力的指標(biāo)，先根據(jù)它們過(guò)去預(yù)報(bào)地震的優(yōu)劣求出每個(gè)指標(biāo)的權(quán)重wi，再以年為單位統(tǒng)計(jì)出每年異常出現(xiàn)與否，若出現(xiàn)為1，未出現(xiàn)為0，最后根據(jù)求出權(quán)重集成值Pi。

選用的16項(xiàng)指標(biāo)列于表10中。表10中還列出每個(gè)指標(biāo)的權(quán)重wi，以年為單位求出的Pi值及與地震對(duì)應(yīng)情況見(jiàn)表11和圖4。從圖4和表11中可以看出，當(dāng)Pi值取0.174時(shí)，地震對(duì)應(yīng)率為52%，概括率為100%；就發(fā)生在遼寧省內(nèi)ML≥5.0地震而言，地區(qū)當(dāng)Pi值取0.265時(shí)，地震對(duì)應(yīng)率為50%，概括率為100%，而當(dāng)Pi值取0.437以上時(shí)，地震對(duì)應(yīng)率為100%，概括率為100%，對(duì)應(yīng)率非常好，可以看出當(dāng)遼寧地區(qū)Pi值為0.437以上時(shí)，遼寧地區(qū)發(fā)生中強(qiáng)地震的危險(xiǎn)性很高。

表6 遼寧及鄰區(qū)流體學(xué)科預(yù)測(cè)指標(biāo)

表7 遼寧及鄰區(qū)形變學(xué)科預(yù)測(cè)指標(biāo)

表7 遼寧及鄰區(qū)形變學(xué)科預(yù)測(cè)指標(biāo)

表8 遼寧及鄰區(qū)電磁學(xué)科預(yù)測(cè)指標(biāo)

表9 遼寧地區(qū)異常及地震對(duì)應(yīng)表

表10 各項(xiàng)指標(biāo)的權(quán)重值

表11 Pi與ML≥5.0地震關(guān)系表

（3）預(yù)報(bào)地點(diǎn)的綜合指標(biāo)—權(quán)重集成值Pp

若前述判斷未來(lái)一年ML≥5.0地震的異常指標(biāo)已出現(xiàn)，那么就需要判斷發(fā)震的危險(xiǎn)區(qū)。選取5項(xiàng)有預(yù)報(bào)發(fā)震地點(diǎn)能力的指標(biāo)，求出它們的權(quán)重，在平面上“集成”，找出綜合預(yù)報(bào)指標(biāo)Pp。選用的指標(biāo)列于表12中。將求出的權(quán)重iw按0.5×0.5平方度分配在有異常的地區(qū)，然后用公式求出每一小塊內(nèi)的權(quán)重集成值Pp，地震就將發(fā)生在Pp值最大的地區(qū)，也就是多種異常重疊的地區(qū)。由于地震空區(qū)、地震條帶由于資料的完備度所限，其從1990年開(kāi)始計(jì)算，故圖5為1999年岫巖M5.4和2013年燈塔M5.1地震的Pp值的空間分布，可見(jiàn)兩次地震均發(fā)生Pp值最大的地區(qū)。

表12 求Pp值的預(yù)報(bào)指標(biāo)

（4）R值評(píng)分

R值評(píng)分方法是由許紹燮[8]提出用來(lái)評(píng)價(jià)地震預(yù)測(cè)效能的方法。根據(jù)R值評(píng)分的定義其為有震報(bào)準(zhǔn)率與預(yù)報(bào)占時(shí)率（或預(yù)報(bào)占時(shí)空域率）之差。其需根據(jù)不同的前兆類型選用，即

以前兆震群指標(biāo)為例，在1970年至2013年期間，其預(yù)報(bào)研究總時(shí)間為43年，預(yù)報(bào)占用時(shí)間為3年，則其預(yù)報(bào)占時(shí)率為0.07，報(bào)準(zhǔn)率為0.3。根據(jù)R值評(píng)分公式計(jì)算得出R值評(píng)分為0.23。若R=1，表示全部報(bào)對(duì)；R=-1，表示全部報(bào)反；R=0，則表示預(yù)報(bào)沒(méi)有起作用。因此，按此方法依次計(jì)算不同地震指標(biāo)的R值評(píng)分，結(jié)果如表13。

3 主要結(jié)論

（1）通過(guò)研究傳統(tǒng)地震學(xué)、數(shù)字地震學(xué)以及定點(diǎn)前兆資料，對(duì)遼寧地區(qū)1970年以來(lái)10次ML≥5.0地震進(jìn)行震例回溯性檢驗(yàn)研究，提出了16項(xiàng)單項(xiàng)預(yù)報(bào)指標(biāo)，并分別以時(shí)間、地點(diǎn)、強(qiáng)度預(yù)測(cè)進(jìn)行分裂，給出其異常判定、震前特征、預(yù)測(cè)實(shí)效以及各項(xiàng)所占的權(quán)重，并利用R值評(píng)分對(duì)單項(xiàng)指標(biāo)的預(yù)測(cè)效能給予評(píng)價(jià)。

表13 不同地震預(yù)測(cè)指標(biāo)的R值評(píng)分

（2）利用權(quán)重集成法將單項(xiàng)指標(biāo)進(jìn)行綜合，提出了預(yù)報(bào)未來(lái)一年有無(wú)ML≥5.0地震的綜合預(yù)報(bào)指標(biāo)Pi，當(dāng)Pi值取0.174時(shí)，地震對(duì)應(yīng)率為52%；就發(fā)生在遼寧省內(nèi)ML≥5.0地震而言，當(dāng)Pi值 取0.265時(shí)，地震對(duì)應(yīng)率為50%，概括率為100%，而當(dāng)Pi值取0.437以上時(shí)，地震對(duì)應(yīng)率為100% ；對(duì)應(yīng)率非常好，可以看出當(dāng)遼寧地區(qū)Pi值為0.437以上時(shí)，遼寧地區(qū)發(fā)生中強(qiáng)地震危險(xiǎn)性很高。

（3）以定性和定量指標(biāo)建立遼寧地區(qū)近期中強(qiáng)震危險(xiǎn)區(qū)的定量指標(biāo)體系，是依據(jù)震例總結(jié)，并從實(shí)際應(yīng)用出發(fā)，在預(yù)報(bào)未能實(shí)現(xiàn)根本解決時(shí)所用的一種經(jīng)驗(yàn)總結(jié)性指標(biāo)體系，具有一定的地區(qū)參考意義。

（4）前述指標(biāo)是在研究震例中挑選提取出來(lái)的。用統(tǒng)計(jì)學(xué)方法進(jìn)行綜合，但尚存一定的局限性，與主觀臆斷性，且任何指標(biāo)對(duì)應(yīng)地震時(shí)都存在一定的漏報(bào)和虛報(bào)，部分指標(biāo)虛報(bào)率還比較高。

[1] 陳棋福，鄭大林，高榮林，等. 中國(guó)震例（1997-1999）[M]. 北京：地震出版社，2003：1-466.

[2] 車時(shí)，蔣海昆，付虹，等. 中國(guó)震例（2003-2006）[M]. 北京：地震出版社，2014：1-774.

[3] 王喜龍，焦明若，王海燕，等.中強(qiáng)震前遼寧地區(qū)流體異常特征與地震預(yù)測(cè)研究[J].地震，2016，36（4）：131-143.

[4] 曹鳳娟，焦明若，王海燕，等. 遼寧地區(qū)前兆異?？臻g分布的集中性與非均勻性研究[J]. 地震，2013，（01）:127-140.

[5] 焦明若，王海燕，曹鳳娟，等. 遼寧地區(qū)地震前兆指標(biāo)與機(jī)理研究[M]. 北京：地震出版社，2011.

[6] 朱鳳鳴，吳戈.一九七五海城地震[M].北京，地震出版社，1982：1-220.

[7] 羅平，范楊，馬宏生，等. 云南地區(qū)近期強(qiáng)震危險(xiǎn)區(qū)定量指標(biāo)體系的研究[J].地震學(xué)報(bào)， 1994，16(增刊):168-174.

[8] 許紹燮. 地震預(yù)報(bào)能力評(píng)分[A].國(guó)家地震局科技監(jiān)測(cè)司.地震預(yù)報(bào)方法實(shí)用化研究文集地震學(xué)專輯[C]. 北京：地震出版社，1989：586-589．

RESEARCH ON THE PREDICTIVE INDEX SYSTEM OF MEDIUM-STRONG EARTHQUAKES IN LIAONING PROVINCE.

JIAO Ming-ruo，XIA Cai-yun， WANG Xi-long，WANG Liang，KONG Xiang-rui，ZHANG Zhi-hong

（Earthquake Administration of Liaoning Province，Liaoning Shenyang 110034，China）

ract: This paper respectively extracted the three elements of earthquake for the anomaly index, which established the anomaly index system in Liaoning province. Through features of the seismic anomaly variations before moderate-strong earthquakes in Liaoning and around areas .We extract the anomaly index ,study the relations between every index and the three elements of earthquake and give the weight of anomaly index. To study typical tendency anomalies and short-term anomaly features in every subject and extract several quantitative indicators for amplitude of the anomaly changes , duration and the features of form changes ,we obtain the relations between anomalies index and three elements of earthquakes .At the same time ,the anomaly index weights are given . Finally ,combined with the geological structures in Liaoning province and the anomaly index weight information integration , we established the index system to predict the strong earthquakes in Liaoning province. It will provide the technical support and basis for judging and tracking the danger area in Liaoning province.

ds: Liaoning area; extracting the abnormal information; abnormal indicators; information integration

P315.3

A

10.13693/j.cnki.cn21-1573.2016.04.001

1674-8565（2016）04-0001-11

2016-08-06

2016-09-28

焦明若(1963-)，男，遼寧省法庫(kù)縣人，博士研究生，研究員，現(xiàn)主要從事地震預(yù)報(bào)、前兆機(jī)理方面的研究工作。E-mail: 359203477@qq.com