張琳琳 敖雪明 聶曉紅

新疆維吾爾自治區(qū)地震局，烏魯木齊市新市區(qū)科學(xué)二街338號(hào) 830011

0 引言

2017年8月9日精河發(fā)生6.6級(jí)地震后，9月16日相距240km的庫(kù)車縣又發(fā)生了5.7級(jí)地震，這2次中強(qiáng)地震均發(fā)生在天山中部地區(qū)，該區(qū)域2011年11月1日尼勒克、鞏留和2012年6月30日新源、和靜相繼發(fā)生了6.0級(jí)、6.6級(jí)地震，但是，近年來在此范圍內(nèi)5.5級(jí)以上中強(qiáng)地震（震群型地震除外）在1.5個(gè)月的時(shí)間內(nèi)接連發(fā)生的情況實(shí)屬罕見。因此，加強(qiáng)對(duì)該區(qū)域地震活動(dòng)性的研究有著重要的意義。

“地震窗”的地震活動(dòng)可用來表征區(qū)域構(gòu)造應(yīng)力場(chǎng)，即主要通過弱震活動(dòng)形式來反映區(qū)域應(yīng)力場(chǎng)變化的特殊構(gòu)造部位?！暗卣鸫啊睘樾陆粘５卣鸨O(jiān)測(cè)預(yù)測(cè)的重要手段之一，20世紀(jì)80年代“七五”科技攻關(guān)前后，新疆地震局開展了“地震窗”預(yù)測(cè)方法研究與應(yīng)用，并取得了不同程度的效果。利用“地震窗”預(yù)測(cè)地震主要基于“地震窗”內(nèi)周、月頻次變化及其與中強(qiáng)地震間的對(duì)應(yīng)關(guān)系來進(jìn)行分析研究的（敖雪明等，1998；鄭兆苾等，1999；陳運(yùn)平等，2001；劉小鳳等，2003；王培玲等，2004；姜秀娥等，1989；薛丁等，2012）。李瑩甄等（2000）對(duì)伽師強(qiáng)震群活動(dòng)期間“喀什地震窗”和“阿圖什余震窗”異常特征進(jìn)行研究發(fā)現(xiàn)，1997年幾組后續(xù)6級(jí)地震前小震頻次明顯增加。王桂嶺等（1991）對(duì)1990年新疆烏恰6.4級(jí)地震前“地震窗”的異常特征進(jìn)行分析，并討論了所得結(jié)果可能的物理機(jī)制。但是，隨著新疆區(qū)域地震臺(tái)的地震數(shù)據(jù)由模擬式記錄改換為數(shù)字化記錄，導(dǎo)致了原有“地震窗”資料的中斷。因此，為滿足短期地震預(yù)測(cè)的需要，以便于更加深入和系統(tǒng)地研究新疆各地震臺(tái)數(shù)字記錄數(shù)據(jù)，有必要重新普查并構(gòu)建新的“地震窗”?！皫?kù)米什地震窗”就是在這項(xiàng)研究中篩選出來的預(yù)測(cè)效能相對(duì)較高的窗口。

以庫(kù)米什臺(tái)為中心，選取2008年以來S-P≤10s范圍內(nèi)ML≥1.0地震的月頻次進(jìn)行分析，開窗指標(biāo)主要以頻次異常為主，2017年8月以來，天山中部地區(qū)相繼發(fā)生了精河6.6級(jí)、庫(kù)車5.7級(jí)地震，2次中強(qiáng)地震前“庫(kù)米什地震窗”均表現(xiàn)出不同程度的異常，通過對(duì)“庫(kù)米什地震窗”異常與北天山中強(qiáng)地震之間的關(guān)系進(jìn)行分析，獲取“庫(kù)米什地震窗”的異常指標(biāo)，同時(shí)對(duì)其預(yù)測(cè)效能進(jìn)行評(píng)價(jià)，以期為北天山中強(qiáng)地震活動(dòng)趨勢(shì)研判提供參考依據(jù)。

1 “庫(kù)米什地震窗”概況

“庫(kù)米什地震窗”位于北天山地震帶，歷史地震活動(dòng)水平較高，1600年以來，史料記載該區(qū)曾經(jīng)發(fā)生過32次6級(jí)以上地震，其中，6級(jí)地震25次；7級(jí)地震 6次；8級(jí)地震1次，即1812年3月8日尼勒克8級(jí)地震。2008年以來，發(fā)現(xiàn)該區(qū)域微震、小震活動(dòng)頻繁，并且這些小震活動(dòng)與北天山地區(qū)中強(qiáng)地震間具有一定相關(guān)性。

北天山地震帶的強(qiáng)震活動(dòng)以近EW向逆斷裂為主，與逆沖-褶皺活動(dòng)構(gòu)造帶關(guān)系密切?！皫?kù)米什地震窗”構(gòu)造區(qū)主要位于博羅科努-阿其克庫(kù)都克斷裂、伊犁盆地北緣、包爾圖斷裂和焉耆斷裂等交匯區(qū)，這些斷層性質(zhì)主要以右旋走滑和逆沖為主（圖1），周圍地下構(gòu)造背景復(fù)雜。

該窗口的選取是以庫(kù)米什地震臺(tái)為中心，以該臺(tái)記錄的S-P≤10s范圍內(nèi)小地震為資料作為常震的“庫(kù)米什地震窗”。庫(kù)米什臺(tái)自2008年以來使用數(shù)字化地震儀記錄數(shù)據(jù)，能記錄到 100km范圍內(nèi)ML≥1.0地震，該臺(tái)地理坐標(biāo)為 42°18′N、88°12′E。

2 資料和分析方法

2.1 研究數(shù)據(jù)

本文所采用的數(shù)據(jù)均取自新疆庫(kù)米什地震臺(tái)數(shù)字地震儀觀測(cè)記錄的微震及小震資料。研究時(shí)段為2008年1月～2017年9月，新疆天山中部地區(qū)為選取的研究區(qū)，“庫(kù)米什地震窗”所對(duì)應(yīng)中強(qiáng)地震的震源參數(shù)來自中國(guó)地震臺(tái)網(wǎng)統(tǒng)一地震目錄。

對(duì)于“庫(kù)米什地震窗”所對(duì)應(yīng)中強(qiáng)地震的統(tǒng)計(jì)方法作如下規(guī)定：若是主-余型地震，則僅取主震；若是震群型、雙震型地震，則取其最大地震或是第1個(gè)地震作為1次大震事件。

2.2 研究方法

Molchan圖表法是20世紀(jì)80～90年代開始研究者對(duì)已開展的經(jīng)驗(yàn)性地震預(yù)測(cè)進(jìn)行科學(xué)總結(jié)以解決固定研究區(qū)強(qiáng)震時(shí)間預(yù)測(cè)問題，并試圖給出概率解釋而逐漸發(fā)展起來的統(tǒng)計(jì)檢驗(yàn)方法（孫麗娜等，2012；蔣長(zhǎng)勝等，2011；Jackson，2004）。該方法涉及的主要變量為：

υ（漏報(bào)率）——預(yù)測(cè)“無震”而實(shí)際發(fā)震的地震數(shù)/空間網(wǎng)格數(shù)與總的實(shí)發(fā)地震數(shù)/所占空間網(wǎng)格數(shù)之比；

τ（異常的時(shí)空占有率）——發(fā)出預(yù)測(cè)“警報(bào)”的時(shí)空范圍與總的時(shí)空范圍之比；

h（擊中數(shù)）——預(yù)測(cè)“有震”而實(shí)際發(fā)震的地震數(shù)/空間網(wǎng)格數(shù)；

H（擊中率）——預(yù)測(cè)“有震”而實(shí)際發(fā)震的地震數(shù)/空間網(wǎng)格數(shù)與總的實(shí)發(fā)地震數(shù)/所占空間網(wǎng)格數(shù)之比。

Molchan圖表法主要使用τ和υ來進(jìn)行統(tǒng)計(jì)評(píng)分，最佳預(yù)測(cè)效能對(duì)應(yīng)在最大預(yù)測(cè)成功（υ→0）下付出最小的代價(jià)（τ→0）。τ在計(jì)算中常根據(jù)預(yù)測(cè)檢驗(yàn)策略或要求的不同進(jìn)行加權(quán)處理。

Molchan圖表法的概率增益定義為

概率增益 Gain在圖表中即為（τ，υ）至（0，1）連線的斜率，而當(dāng)連線為（0，1）至（1，0）的直線時(shí)，Gain=1，則表示預(yù)測(cè)方法無統(tǒng)計(jì)顯著性。在 Molchan圖表中繪制τ～υ曲線，其中，τ∈［0，1］，υ∈［0，1］。地震預(yù)測(cè)效果的評(píng)估主要是在直觀上比較（τ，υ）曲線與 Molchan圖表的邊界線τ=0、υ=0所包圍的面積，面積越小，則預(yù)測(cè)效果越好。同時(shí)，還需考慮擊中數(shù)h所對(duì)應(yīng)的顯著性水平。

圖1 北天山地震帶及震情窗口位置分布

3 “庫(kù)米什地震窗”的映震分析

3.1 “庫(kù)米什地震窗”預(yù)測(cè)指標(biāo)的確定

以庫(kù)米什地震臺(tái)為中心，分析S-P≤10s范圍內(nèi)ML≥1.0地震數(shù)據(jù)資料，研究該地區(qū)的小震活動(dòng)，即為“庫(kù)米什地震窗”（或“庫(kù)米什窗”）的小震活動(dòng)。根據(jù)歷史震例和異常的最佳對(duì)應(yīng)原則的分析，定義“庫(kù)米什地震窗”小震月頻次≥54（即54次/月）為“庫(kù)米什地震窗”的異常標(biāo)準(zhǔn)。

圖2為“庫(kù)米什窗”2008年1月～2017年9月小震月頻次的變化及其與附近中強(qiáng)地震間的關(guān)系。由圖2可見，小震月頻次在小于54次的范圍內(nèi)隨機(jī)波動(dòng)，但是，其中有7個(gè)點(diǎn)呈現(xiàn)出大于54次/月的異常高值，而其后6個(gè)月內(nèi)有6次MS≥5.5中強(qiáng)地震與之對(duì)應(yīng)。這6次地震震中均位于北天山，其中，4次對(duì)應(yīng)的中強(qiáng)地震震中位于北天山西段。小震月頻次異常特征量以及所對(duì)應(yīng)的中強(qiáng)地震參數(shù)見表1。

圖2 “庫(kù)米什地震窗”小震月頻次與中強(qiáng)地震（MS≥5.5）間的對(duì)應(yīng)關(guān)系

表1 “庫(kù)米什地震窗”小震月頻次異常與地震間的對(duì)應(yīng)關(guān)系

“庫(kù)米什地震窗”所對(duì)應(yīng)的6次中強(qiáng)地震震中距為60～490km，其中，小于400km的有4次。分析“庫(kù)米什地震窗”的實(shí)際資料異常與地震對(duì)應(yīng)情況發(fā)現(xiàn)，7組異常對(duì)應(yīng)了6次中強(qiáng)地震。依據(jù)許紹燮（1989）的地震預(yù)報(bào)能力評(píng)分方法，對(duì)“庫(kù)米什地震窗”窗口異常進(jìn)行R值評(píng)分檢驗(yàn)，所得R值評(píng)分為0.45，高于具有97.5%置信度水平的R值（0.356），表明該方法的短期預(yù)測(cè)效能較高。

2017年以來，天山中部地區(qū)5級(jí)以上中強(qiáng)地震活動(dòng)水平較高，尤其自8月以來，8月9日精河、9月16日庫(kù)車相繼發(fā)生了MS6.6、MS5.7地震，2次中強(qiáng)地震前“庫(kù)米什地震窗”均出現(xiàn)了異常，其中，精河6.6級(jí)地震前3.3個(gè)月“庫(kù)米什地震窗”出現(xiàn)異常，震情窗小震月頻次4月份出現(xiàn)55次/月的異常，“庫(kù)米什地震窗”8月份出現(xiàn)月頻次54次的異常，異常結(jié)束后15天庫(kù)車發(fā)生了5.7級(jí)地震。由表1可見，“庫(kù)米什地震窗”對(duì)應(yīng)中強(qiáng)地震的震中距有4次為300～490km，而且這4次中強(qiáng)地震均位于天山中部地區(qū)。距離最近的僅為60km，即在“庫(kù)米什地震窗”附近區(qū)域發(fā)震。

3.2 基于Molchan模型的“庫(kù)米什地震窗”預(yù)測(cè)效能評(píng)價(jià)

首先，將2008年以來“庫(kù)米什地震窗”的小震月頻次數(shù)據(jù)以及所對(duì)應(yīng)的中強(qiáng)地震進(jìn)行整理；然后，基于Matlab程序?qū)π≌鹪骂l次數(shù)據(jù)和所對(duì)應(yīng)地震之間的關(guān)系進(jìn)行分析；再后，在統(tǒng)計(jì)檢驗(yàn)時(shí)，可以根據(jù)映震時(shí)間的長(zhǎng)短進(jìn)行不同開窗指標(biāo)的驗(yàn)證；最后，獲取用于繪制Molchan檢驗(yàn)圖的數(shù)據(jù)，同時(shí)，依據(jù)地震漏報(bào)率與預(yù)測(cè)地震時(shí)空占有率間的比率進(jìn)行效能評(píng)價(jià)。

當(dāng)Gain=1時(shí)，則表示預(yù)測(cè)方法無統(tǒng)計(jì)顯著性。預(yù)測(cè)效果以（τ，υ）曲線與Molchan圖表的邊界線τ=0、υ=0所包圍的面積來表示，面積越小，則預(yù)測(cè)效果越好（圖3）。由圖3可見，（τ，υ）曲線與Molchan圖表的邊界線τ=0、υ=0所包圍的面積多數(shù)位于 Gain=1以左的范圍內(nèi)，預(yù)測(cè)效果相對(duì)較好。同時(shí)，結(jié)合上述R值評(píng)分檢驗(yàn)結(jié)果分析認(rèn)為，“庫(kù)米什地震窗”預(yù)測(cè)效能較好，可將其作為中短期異常指標(biāo)進(jìn)行日常跟蹤。

圖3 “庫(kù)米什地震窗”Molchan檢驗(yàn)效果

4 結(jié)論與討論

（1）以庫(kù)米什地震臺(tái)為中心，選取S-P≤10s范圍內(nèi)ML≥1.0小地震的月頻次，當(dāng)小震月頻次≥54次時(shí)即為異常。2008年以來共出現(xiàn)異常7次，6次異常結(jié)束后均有5.5級(jí)以上地震發(fā)生，異常結(jié)束后中強(qiáng)地震發(fā)震的優(yōu)勢(shì)時(shí)段為6個(gè)月，優(yōu)勢(shì)發(fā)震區(qū)域?yàn)樘焐街胁康貐^(qū)。2017年8月份以來，精河6.6級(jí)、庫(kù)車5.7級(jí)2次地震前“庫(kù)米什地震窗”均出現(xiàn)異常。

（2）當(dāng)區(qū)域應(yīng)力場(chǎng)增強(qiáng)時(shí)，潛在地震區(qū)均為應(yīng)力集中點(diǎn)，而“庫(kù)米什地震窗”介質(zhì)的抗剪切強(qiáng)度因降低而出現(xiàn)小震頻次增高異常，且可能反復(fù)出現(xiàn)，即多點(diǎn)的月頻次異常對(duì)應(yīng)地震孕育過程。通過“庫(kù)米什地震窗”的異常特征可以判斷北天山地震發(fā)震時(shí)間，但對(duì)具體發(fā)震地點(diǎn)難以判定。引起區(qū)域應(yīng)力場(chǎng)增強(qiáng)的因素較多，如大震應(yīng)力觸發(fā)、地震構(gòu)造活動(dòng)等（張彬等，2008）。

（3）基于Molchan模型的“庫(kù)米什地震窗”預(yù)測(cè)效能評(píng)價(jià)結(jié)果表明，“庫(kù)米什地震窗”預(yù)測(cè)效能相對(duì)較好，可將其作為中短期異常指標(biāo)進(jìn)行日常跟蹤?；贛olchan模型對(duì)“庫(kù)米什地震窗”的檢驗(yàn)屬于初步嘗試，一些參數(shù)的設(shè)置需要經(jīng)過不斷的調(diào)試進(jìn)行完善，如開窗指標(biāo)的條件和時(shí)間，都需要經(jīng)過對(duì)不同“地震窗”預(yù)測(cè)效能的檢驗(yàn)評(píng)價(jià)進(jìn)行總結(jié)。