利用兩次新疆于田MS7.3地震對(duì)圖像信息方法可預(yù)測(cè)性的檢驗(yàn)

夏彩韻 張永仙 張小濤 吳永加

1）中國(guó)北京100036中國(guó)地震局地震預(yù)測(cè)研究所

2）中國(guó)北京100045中國(guó)地震臺(tái)網(wǎng)中心

引言

圖像信息方法（pattern informatics method，簡(jiǎn)寫為PI）是由Rundle等（2000a，b）提出的用于地震活動(dòng)性研究的地震統(tǒng)計(jì)物理方法，并一直應(yīng)用于地震活動(dòng)的中長(zhǎng)期預(yù)測(cè)研究中，具有比較高的空間分辨能力和較好的預(yù)測(cè)效果（蔣長(zhǎng)勝等，2009；Zhang et al，2013）.最初該方法被稱為“相位動(dòng)力學(xué)概率變化”（phase dynamical probability change，簡(jiǎn)寫為PDPC）.他們把地震活動(dòng)看作是自組織閾值系統(tǒng)的一個(gè)特例，并利用PDPC方法計(jì)算出2000年1月—2009年12月美國(guó)南加州有可能發(fā)生M≥5.0地震的地區(qū).事后檢驗(yàn)大約90%的目標(biāo)地震發(fā)生在所預(yù)測(cè)的危險(xiǎn)區(qū)內(nèi)（Rundle，2000a，2002）.研究結(jié)果表明，該方法在地震中長(zhǎng)期危險(xiǎn)區(qū)預(yù)測(cè)中具有較高的空間分辨率.Holliday等（2005）利用PI方法對(duì)全球2000年1月1日—2009年12月31日發(fā)生M≥7.0地震的可能性進(jìn)行了預(yù)測(cè)，回溯性檢驗(yàn)結(jié)果顯示，2004年12月23日麥闊里島M8.1地震和2004年12月26日蘇門達(dá)臘MS9.0地震的震中均位于“地震熱點(diǎn)”上或其附近.Nanjo等（2006）利用PI方法對(duì)日本（33°—38°N，136°—142°E）2000年1月—2009年12月發(fā)生 M≥5.0地震的可能性進(jìn)行預(yù)測(cè)，回溯性檢驗(yàn)結(jié)果顯示，2000年6—8月發(fā)生在Miyake島的MS7.0地震群位于“地震熱點(diǎn)”上，2004年10月23日新潟（Niigata）ML6.8地震及其余震的震中均位于“地震熱點(diǎn)”附近.Chen等（2005）將PI方法應(yīng)用到我國(guó)臺(tái)灣地區(qū)（21°—26°N，119°—123°E）地震危險(xiǎn)性的研究中，回溯性檢驗(yàn)結(jié)果顯示1999年9月21日臺(tái)灣集集MS7.6地震的震中位于計(jì)算所得到的“地震熱點(diǎn)”內(nèi).蔣長(zhǎng)勝和吳忠良（2008）通過回溯性研究認(rèn)為，PI方法對(duì)川滇地區(qū)的中長(zhǎng)期地震活動(dòng)具有一定的預(yù)測(cè)能力.張小濤（2009）用PI方法對(duì)我國(guó)云南、大華北等地區(qū)進(jìn)行中長(zhǎng)期地震危險(xiǎn)性的預(yù)測(cè)研究，通過多次回溯性檢驗(yàn)篩選出適合于不同研究區(qū)域的計(jì)算參數(shù)，并利用這些參數(shù)對(duì)研究區(qū)域的地震危險(xiǎn)性進(jìn)行預(yù)測(cè).Zhang等（2013）以2008年于田MS7.3地震和汶川MS8.0地震為例，對(duì)我國(guó)西部地區(qū)（20°—50°N，70°—110°E）分別進(jìn)行了回溯性預(yù)測(cè)檢驗(yàn)研究，通過調(diào)整PI方法中的計(jì)算參數(shù)范圍，利用R值評(píng)分方法（Xu，1989；石耀霖等，2000）和ROC方法（Swets，1973；Molchan，1997）評(píng)價(jià)兩次地震的預(yù)測(cè)效能，給出預(yù)測(cè)效能最高的計(jì)算參數(shù)選取范圍，為后續(xù)PI方法在中國(guó)西部的應(yīng)用提供了參考.

新疆是中國(guó)強(qiáng)震多發(fā)區(qū)之一，昆侖山、天山及阿爾泰山均為地震活動(dòng)較為頻繁的區(qū)域（尹光華等，2008）.2014年2月12日新疆于田（36.1°N，82.5°E）發(fā)生的 MS7.3地震，是繼2008年3月21日新疆于田（35.6°N，81.6°E）發(fā)生MS7.3地震6年后的又一次大地震.2008年新疆于田MS7.3地震震中位于西昆侖地震帶與阿爾金斷裂帶的交匯區(qū)（尹光華等，2008），而2014年于田MS7.3地震震中與其僅相距100km，也在阿爾金斷裂帶上（王曉欣等，2014）.從大的活動(dòng)地塊來說，這一地區(qū)屬于青藏和西域兩大活動(dòng)地塊邊界區(qū)（張培震，1999）.

Zhang等（2013）的研究結(jié)果表明，2008年于田MS7.3地震震中附近出現(xiàn)了PI熱點(diǎn)圖像，由此引發(fā)的問題是，本次于田MS7.3地震前震中附近是否存在PI熱點(diǎn)圖像.此外，在相隔100km的于田地區(qū)不到6年的時(shí)間內(nèi)發(fā)生了兩次MS7.3地震，這為檢驗(yàn)PI方法的可預(yù)測(cè)性提供了很好的檢驗(yàn)條件.由于兩次于田MS7.3地震均位于青藏構(gòu)造塊體西北緣，本文選取青藏構(gòu)造塊體為研究區(qū)域，采用以往預(yù)測(cè)M7.0地震的較好計(jì)算參數(shù)（網(wǎng)格尺度為1平方度，預(yù)測(cè)時(shí)間窗口為8年）（Zhang et al，2013），系統(tǒng)研究青藏構(gòu)造塊體1993年以來每8年時(shí)間窗的“地震熱點(diǎn)”圖像，特別研究?jī)纱斡谔颩S7.3地震前后的圖像信息演化過程，并通過ROC方法檢驗(yàn)和R值評(píng)分方法，對(duì)PI方法在該地區(qū)的預(yù)測(cè)效能進(jìn)行評(píng)價(jià)；在此基礎(chǔ)上嘗試對(duì)青藏構(gòu)造塊體潛在的MS7.0地震危險(xiǎn)區(qū)域進(jìn)行預(yù)測(cè).

1 PI方法

Rundle（2000b）提出的PI方法是通過將研究區(qū)域進(jìn)行時(shí)空網(wǎng)格劃分，分別對(duì)落入每個(gè)網(wǎng)格內(nèi)的地震活動(dòng)構(gòu)建其頻度時(shí)間序列，通過統(tǒng)計(jì)方法計(jì)算每個(gè)網(wǎng)格中的地震活動(dòng)強(qiáng)度函數(shù)對(duì)背景的偏離程度，并歸算到顯著地震事件發(fā)生的概率，再扣除背景概率后檢測(cè)發(fā)震概率高的網(wǎng)格，即“地震熱點(diǎn)”.本文將參考Holliday等（2005）對(duì)PI方法的實(shí)現(xiàn)過程，系統(tǒng)計(jì)算青藏構(gòu)造塊體1993年以來的“地震熱點(diǎn)”圖像.

2 研究區(qū)域及計(jì)算參數(shù)選取

2.1 研究區(qū)域

本文選取青藏構(gòu)造塊體所包含的矩形區(qū)域（21.0°—41.0°N，74.0°—106.0°E）為研究對(duì)象.該區(qū)域受印度－歐亞板塊強(qiáng)烈碰撞作用的影響，新生代以來先后形成了青藏活動(dòng)地塊區(qū)的多個(gè)活動(dòng)構(gòu)造地塊（尹光華等，2008），使得該區(qū)域地質(zhì)構(gòu)造復(fù)雜，是現(xiàn)今地殼運(yùn)動(dòng)最為活躍的地區(qū)之一，也是強(qiáng)震活動(dòng)比較復(fù)雜的地區(qū)之一（徐錫偉等，2005）.據(jù)中國(guó)地震臺(tái)網(wǎng)中心①http：∥10.5.202.22／bianmu／validate.jsp.統(tǒng)計(jì)自1970年1月1日以來，該研究區(qū)域記錄到M6.0—6.9地震135次，M7.0—7.9地震17次，M≥8.0地震2次.本文研究區(qū)域中，1993—2014年共發(fā)生9次M≥7.0地震，具體震例信息見表1.

表1 本文研究震例Table 1 Earthquake examples used in this study

2.2 PI方法計(jì)算參數(shù)的設(shè)定

本文使用的地震目錄引自中國(guó)地震臺(tái)網(wǎng)中心①http：∥10.5.202.22／bianmu／validate.jsp..在整個(gè)青藏構(gòu)造塊體內(nèi)，西藏地區(qū)的監(jiān)測(cè)能力較低，目前只能監(jiān)測(cè)到M≥3.0以上地震②黃志斌在中國(guó)地震臺(tái)網(wǎng)內(nèi)部交流圖..由于本文研究方法需使用1970年以來的地震目錄，因此對(duì)西藏地區(qū)（26.5°—36.5°N，77°—100°E）1970年以來的地震目錄做G-R關(guān)系圖（圖1），顯示該區(qū)的完備震級(jí)下限為ML4.5.在PI方法中，選取的地震目錄下限M0與預(yù)測(cè)的目標(biāo)地震震級(jí)Mf有關(guān)（Mf＝M0＋2）（Tiampo et al，2002；Rundle et al，2002；Holliday et al，2005），因此采用Mc＝ML5.0為研究區(qū)域的震級(jí)下限，不僅可以滿足選取地震目錄完備性的要求，還可以滿足研究區(qū)域M≥7.0地震的預(yù)測(cè)需求.

本文在采用PI方法研究?jī)纱斡谔颩S7.3地震時(shí)，參照了Zhang等（2013）以大陸西部為研究區(qū)域，分析汶川MS8.0、于田MS7.3地震前的PI異常圖像演化過程的研究結(jié)果，即計(jì)算參數(shù)閾值取0.6、網(wǎng)格大小為1°×1°和預(yù)測(cè)時(shí)間段為8年.為了檢驗(yàn)兩次于田地震熱點(diǎn)演化過程的穩(wěn)定性，我們把預(yù)測(cè)時(shí)間窗以1年為步長(zhǎng)逐年向前滑動(dòng)，給出1993年以來的地震熱點(diǎn)圖像，以獲取該研究地區(qū)較長(zhǎng)時(shí)間的PI圖像演化過程，探索PI圖像演化與M≥7.0地震的關(guān)系，特別是兩次于田MS7.3地震前后的PI演化圖像.

圖1 青藏高原地區(qū)地震目錄G-R關(guān)系圖Fig.1 Completeness of catalogue test by G-R relation in Qinghai-Xizang（Tibet）region

3 計(jì)算結(jié)果分析

依據(jù)2.2節(jié)的計(jì)算參數(shù)，共獲得22個(gè)預(yù)測(cè)時(shí)間段的PI地震熱點(diǎn)分布，如圖2所示.在這些預(yù)測(cè)時(shí)段內(nèi)，研究區(qū)域共發(fā)生9次M≥7.0地震（表1），每個(gè)預(yù)測(cè)時(shí)段內(nèi)所發(fā)生的M≥7.0地震均標(biāo)注于預(yù)測(cè)圖中（圖2中空心圓圈）.

分析圖2中地震熱點(diǎn)的演化圖像與預(yù)測(cè)時(shí)段所發(fā)生的M≥7.0地震，可得到以下認(rèn)識(shí)：

1）2008年和2014年兩次于田MS7.3地震前震中所在網(wǎng)格及鄰近網(wǎng)格均出現(xiàn)地震熱點(diǎn).從圖2i開始，PI異常圖像顯示出西昆侖地震帶與阿爾金斷裂帶的交匯地區(qū)存在地震熱點(diǎn)；將預(yù)測(cè)時(shí)間窗向后推一年（圖2j），該處地震熱點(diǎn)分布范圍未改變，但發(fā)震概率增高.將預(yù)測(cè)時(shí)間窗繼續(xù)向后逐年類推（圖2k-n），該處地震熱點(diǎn)收縮至于田附近，可以明顯看到2008年于田MS7.3地震發(fā)生前的PI熱點(diǎn)演變過程.隨著預(yù)測(cè)時(shí)間窗繼續(xù)逐年向后類推，該地區(qū)熱點(diǎn)依然繼續(xù)存在，對(duì)2014年于田MS7.3地震的發(fā)生具有指示意義.

圖2 不同預(yù)測(cè)時(shí)間窗口的PI地震熱點(diǎn)演化圖 （a-f）小色塊為計(jì)算所得熱點(diǎn)，閾值lg（ΔP／ΔPmax）?。?.6，空心圈表示在相應(yīng)預(yù)測(cè)時(shí)段內(nèi)該地區(qū)發(fā)生的M＞7.0地震Fig.2 The PI anomalies in different periods（a-f）The colored boxes represent the calculated hotspots under the threshold possibility lg（ΔP／ΔPmax）is taken as－0.6.Circles represent the M＞7.0earthquakes in different prediction period

圖2 不同預(yù)測(cè)時(shí)間窗口的PI地震熱點(diǎn)演化圖 （g－－p）小色塊為計(jì)算所得熱點(diǎn)，閾值lg（ΔP／ΔPmax）?。?.6，空心圈表示在相應(yīng)預(yù)測(cè)時(shí)段內(nèi)該地區(qū)發(fā)生的M＞7.0地震Fig.2 The PI anomalies in different periods（g－－p）The colored boxes represent the calculated hotspots under the threshold possibility lg（ΔP／ΔPmax）is taken as－0.6.Circles represent the M＞7.0earthquakes in different prediction period

圖2 不同預(yù)測(cè)時(shí)間窗口的PI地震熱點(diǎn)演化圖 （q-v）小色塊為計(jì)算所得熱點(diǎn)，閾值lg（ΔP／ΔPmax）?。?.6，空心圈表示在相應(yīng)預(yù)測(cè)時(shí)段內(nèi)該地區(qū)發(fā)生的M＞7.0地震Fig.2 The PI anomalies in different periods（q-v）The colored boxes represent the calculated hotspots under the threshold possibility lg（ΔP／ΔPmax）is taken as－0.6.Circles represent the M＞7.0earthquakes in different prediction period

從圖2o開始，2014年于田MS7.3地震落入預(yù)測(cè)時(shí)間窗內(nèi)的地震熱點(diǎn)上，但預(yù)測(cè)時(shí)間窗繼續(xù)向后類推時(shí)，該處的地震熱點(diǎn)變?nèi)酰▓D2q-t），直到預(yù)測(cè)時(shí)間窗為2013年1月1日—2020年12月31日（圖2u-v）開始，該處的地震熱點(diǎn)開始變強(qiáng)，2014年于田MS7.3地震重新落入地震熱點(diǎn)內(nèi)，這可能預(yù)示著該處在未來的時(shí)間仍有發(fā)生強(qiáng)震的可能，值得我們關(guān)注.

2）研究區(qū)域發(fā)生的9次M≥7.0地震中，有6次地震前震中所在網(wǎng)格和鄰近網(wǎng)格持續(xù)出現(xiàn)穩(wěn)定的地震熱點(diǎn)，2次地震前在個(gè)別預(yù)測(cè)時(shí)間段內(nèi)出現(xiàn)地震熱點(diǎn)，1次地震前震中所在網(wǎng)格和鄰近網(wǎng)格無地震熱點(diǎn)出現(xiàn).圖2b-d顯示出1997年西藏瑪尼MS7.5地震的周圍存在地震熱點(diǎn)，并逐年增強(qiáng)；2001年昆侖山口西MS8.1地震也在地震熱點(diǎn)上.將預(yù)測(cè)時(shí)間窗繼續(xù)向后類推（圖2e-i），可以看到昆侖山口西MS8.1地震所在的地震熱點(diǎn)呈現(xiàn)出從弱到強(qiáng)，再變?nèi)醯囊粋€(gè)演化過程.圖2l顯示出2010年青海玉樹MS7.1地震在地震熱點(diǎn)上，隨著預(yù)測(cè)時(shí)間窗向后類推（圖2o），該處的地震熱點(diǎn)經(jīng)歷了逐漸變?nèi)踔敝料У难莼^程.對(duì)于2013年蘆山MS7.0地震而言，雖然其僅穩(wěn)定出現(xiàn)在由2010年1月1日類推到2013年1月1日的3年預(yù)測(cè)時(shí)間窗（圖2r-u）的計(jì)算熱點(diǎn)上，但本文與張小濤等（2014）以川滇地區(qū)作為研究區(qū)域的計(jì)算結(jié)果進(jìn)行對(duì)比，發(fā)現(xiàn)其異常演化過程基本一致.而喀喇昆侖山口MS7.1及汶川MS8.0地震均在個(gè)別預(yù)測(cè)時(shí)間窗內(nèi)出現(xiàn)地震熱點(diǎn)（圖2d，p），麗江MS7.0地震則沒有地震熱點(diǎn)出現(xiàn).通過前人利用PI方法對(duì)川滇等地區(qū)進(jìn)行回溯性檢驗(yàn)，這些地區(qū)在發(fā)震前后均存在地震熱點(diǎn)（蔣長(zhǎng)勝等，2009；蔣卉等；2013；Zhang et al，2013）.對(duì)于造成這一差異的原因，可能是與研究區(qū)域的選取及計(jì)算參數(shù)的設(shè)定有很大的關(guān)系，尚需進(jìn)一步深入討論.

3）兩次于田MS7.3地震均發(fā)生在震中附近網(wǎng)格發(fā)震概率值升高到極大值的年份.由于每個(gè)預(yù)測(cè)時(shí)間窗覆蓋了8年，為了定量刻畫兩次于田MS7.3地震震中附近網(wǎng)格發(fā)震概率的時(shí)間變化過程與地震之間的關(guān)系，本文進(jìn)行了如下處理：分別以兩次于田地震震中所在網(wǎng)格為中心，為避免地震目錄定位帶來的誤差問題，地震發(fā)生在閾值之上的“熱點(diǎn)”及周圍8個(gè)網(wǎng)格內(nèi)的均被認(rèn)為是“命中”，鄰近的8個(gè)方框?qū)狳c(diǎn)框而言稱之為“摩爾鄰近”（Wolfram，2002），所以我們?cè)诿恳粋€(gè)預(yù)測(cè)時(shí)間窗內(nèi)將包括震中在內(nèi)的9個(gè)網(wǎng)格的發(fā)震概率取平均值，然后平均到8年窗口的每一個(gè)年份，再把每一個(gè)年份的8個(gè)概率值累加，獲得每一個(gè)年份的平均發(fā)震概率值.對(duì)于2008年于田MS7.3地震，網(wǎng)格選取范圍為（34.0°—37.0°N，80.0°—83.0°E）；對(duì)于2014年于田MS7.3地震，網(wǎng)格選取范圍為（35.0°—38.0°N，81.0°—84.0°E）.兩次于田地震震中所在網(wǎng)格和摩爾鄰近網(wǎng)格的平均發(fā)震概率隨時(shí)間變化如圖3所示.可以看出：對(duì)于2008年于田地震而言，發(fā)震概率逐年增高，到2008年達(dá)到一個(gè)峰值，隨后開始逐年下降，于田地震發(fā)生在概率峰值的2008年；對(duì)于2014年于田地震而言，該地區(qū)發(fā)震概率逐年上升，2004年達(dá)到高值，之后持續(xù)波動(dòng)，在2014年時(shí)達(dá)到最高，這次于田地震也發(fā)生在概率峰值的2014年.值得注意的是：2008年于田MS7.3震中附近的概率值也出現(xiàn)升高現(xiàn)象，可能是由于其網(wǎng)格范圍與2014年于田MS7.3地震網(wǎng)格范圍重疊，包含了2014年地震孕育信息的原因.此外，2014年于田MS7.3地震震中附近的發(fā)震概率一直高于2008年于田MS7.3地震震中附近的發(fā)震概率，且時(shí)間變化過程具有一定的相關(guān)性.這值得深入研究.

圖3 2008年于田MS7.3地震（實(shí)心圓連線）和2014年2月12日于田MS7.3地震（三角形連線）的平均發(fā)震概率隨時(shí)間的變化Fig.3 Variation of the average seismogenic probability for the 2008Yutian MS7.3（solid circle curve）and 2014Yutian MS7.3（solid triangle curve）earthquakes with time

4 PI預(yù)測(cè)效能檢驗(yàn)

為了定量評(píng)價(jià)PI方法在不同預(yù)測(cè)時(shí)間窗的預(yù)測(cè)效能，本文采用ROC曲線方法和R值評(píng)分方法進(jìn)行統(tǒng)計(jì)檢驗(yàn).

1）ROC方法是用在所設(shè)定的危險(xiǎn)性概率閾值下預(yù)測(cè)結(jié)果的“命中率”和“虛報(bào)率”來評(píng)估計(jì)算方法的預(yù)測(cè)效能.“命中率”是由預(yù)測(cè)“有震”而實(shí)際發(fā)震的空間網(wǎng)格數(shù)與總體實(shí)際發(fā)生地震所占的空間網(wǎng)格數(shù)之比來表示；“虛報(bào)率”是由預(yù)測(cè)“有震”而沒有發(fā)生地震的空間網(wǎng)格數(shù)與實(shí)際未發(fā)生地震的空間網(wǎng)格數(shù)之比來表示（Holliday et al，2005）.因此，本文參照Zhang等（2013）的研究方法，將ROC檢驗(yàn)圖中PI方法預(yù)測(cè)實(shí)線與隨機(jī)預(yù)測(cè)虛線所包絡(luò)的面積定義為有效預(yù)測(cè)系數(shù)Ef，即Ef值越大，則表明預(yù)測(cè)效能越好.由于篇幅有限，本文僅給出ROC曲線評(píng)價(jià)最高和最低的兩張圖（圖4）.可以看出，PI方法的預(yù)測(cè)效能顯著優(yōu)于隨機(jī)預(yù)測(cè).表2給出了PI方法在不同預(yù)測(cè)時(shí)間窗的有效預(yù)測(cè)系數(shù)Ef值.可以看出，PI方法的Ef值均比較穩(wěn)定，是可以應(yīng)用于M7.0地震的中長(zhǎng)期預(yù)測(cè)的較好方法.

圖4 PI算法回溯性研究的ROC曲線檢驗(yàn)圖Fig.4 ROC test for the retrospective study using PI method

表2 不同預(yù)測(cè)時(shí)間窗下PI方法的R值評(píng)分和Ef值Table 2 Rscore and Efvalues for different predict windows by PI method

2）R值評(píng)分方法是由Xu（1989）提出用來評(píng)價(jià)地震預(yù)測(cè)效能的方法，目前常用的計(jì)算方式是將研究區(qū)域進(jìn)行1°×1°網(wǎng)格化，R值即為報(bào)震的成功率減去危險(xiǎn)區(qū)占用的預(yù)報(bào)網(wǎng)格數(shù)與總的網(wǎng)格數(shù)之比，R值越大，表明預(yù)測(cè)效果越好（石耀霖等，2000）.本文將研究區(qū)域中的中國(guó)大陸按1°×1°劃分網(wǎng)格，共占390個(gè)網(wǎng)格，以圖2h為例，此預(yù)測(cè)時(shí)間段窗口中只發(fā)生1次地震并命中在地震熱點(diǎn)上，其報(bào)震成功率為1，地震熱點(diǎn)占54個(gè)網(wǎng)格，其R值為

按此方法依次計(jì)算不同預(yù)測(cè)時(shí)間窗的R值，其評(píng)價(jià)結(jié)果列于表2.可以看出，R值評(píng)分僅有一個(gè)時(shí)間窗低于0，其余12個(gè)時(shí)間窗均大于0，一些預(yù)測(cè)時(shí)段的R值評(píng)分高達(dá)0.8，表明PI方法對(duì)M7.0地震的可預(yù)測(cè)性.

5 討論與結(jié)論

1）本文利用PI方法，參照前人研究工作的結(jié)果選取預(yù)測(cè)效能較高的計(jì)算參數(shù)，研究了于田兩次MS7.3地震前后青藏高原塊區(qū)的異常圖像信息演化過程.結(jié)果表明，1993—2014年，研究區(qū)域內(nèi)共發(fā)生9次M≥7.0地震，其中包括于田兩次MS7.3地震在內(nèi)的6次地震震前均出現(xiàn)持續(xù)穩(wěn)定的地震熱點(diǎn)，兩次地震在個(gè)別預(yù)測(cè)時(shí)間段里出現(xiàn)地震熱點(diǎn)，1次地震沒有“命中”在地震熱點(diǎn)及其摩爾鄰近區(qū)內(nèi).本文進(jìn)一步采用ROC曲線方法和R值評(píng)分方法進(jìn)行檢驗(yàn)，其結(jié)果均表明PI方法在該地區(qū)具有較高的預(yù)測(cè)效能.于田兩次地震均發(fā)生在平均發(fā)震概率的峰值點(diǎn)年份，表明對(duì)PI熱點(diǎn)的發(fā)震概率時(shí)間過程追蹤可能是對(duì)M7.0地震時(shí)間預(yù)測(cè)逐步逼近的途徑.

2）PI方法的地震熱點(diǎn)的檢驗(yàn)是將區(qū)域網(wǎng)格化，進(jìn)而檢驗(yàn)各網(wǎng)格地震活動(dòng)狀態(tài)偏離度高的區(qū)域.在本文的研究中，1996年2月3日麗江MS7.0及2008年5月12日汶川MS8.0等3次地震震前并沒有出現(xiàn)持續(xù)穩(wěn)定的地震熱點(diǎn)演化過程，該結(jié)果與前人的研究結(jié)果有所不同.其原因可能是由于地理范圍的選取及計(jì)算參數(shù)的差異性造成的，這與強(qiáng)震附近的地震目錄是否記錄了較為完整的地震活動(dòng)性資料，以及強(qiáng)震前若干年沒有中小地震活動(dòng)可能都有關(guān)系.對(duì)于如何選取合適的地理范圍作為研究區(qū)域及合理的參數(shù)設(shè)定，這一問題仍值得深入研究與探討.

3）覆蓋2014年以后的預(yù)測(cè)時(shí)間窗內(nèi)還存在3個(gè)地震熱點(diǎn)顯著的區(qū)域（圖2p-v），即川滇藏交界（28.0°—32.0°N，97.0°—103.0°E）、甘青川交界（30.0°—34.0°N，100.0°—104.0°E）、滇西南－滇西地區(qū)（21.0°—26.0°N，101.0°—104.0°E）.這3個(gè)區(qū)域有熱點(diǎn)持續(xù)穩(wěn)定的異常演化過程，可能是未來發(fā)生M7.0地震概率較高的地區(qū).

4）本文中計(jì)算所得到的“地震熱點(diǎn)”是描述與自己背景歷史相比及計(jì)算區(qū)域所有網(wǎng)格相比地震活動(dòng)狀態(tài)偏離度較大的地區(qū).從圖2中可看出“熱點(diǎn)”還具有一定的不穩(wěn)定性，這可能與滑動(dòng)時(shí)間窗所包含的地震活動(dòng)信息不同有關(guān).根據(jù)這些“熱點(diǎn)”的時(shí)空演化過程可以考察出地震長(zhǎng)期孕育的過程，對(duì)地震的發(fā)生具有一定的指示意義.

蔣卉，吳忠良，馬騰飛，蔣長(zhǎng)勝.2013.對(duì)圖像信息學(xué)（PI）算法的一個(gè)回溯性預(yù)測(cè)檢驗(yàn)：四川蘆山7.0級(jí)地震［J］.物理，42（5）：334-340.

Jiang H，Wu Z L，Ma T F，Jiang C S.2013.Retrospective test of the PI forecast：Case study of the April 20，2013，Lushan，Sichuan，China，MS7.0earthquake［J］.Physics，42（5）：334-340（in Chinese）.

蔣長(zhǎng)勝，吳忠良.2008.對(duì)地震預(yù)測(cè)的一個(gè)統(tǒng)計(jì)物理算法在川滇地區(qū)的回溯性預(yù)測(cè)檢驗(yàn)［J］.中國(guó)科學(xué)：D輯，38（7）：852-861.

Jiang C S，Wu Z L.2008.Sichuan-Yunnan：A approach based on statistical physics and retrospective forecast test［J］.Science in China：Series D，38（7）：852-861（in Chinese）.

蔣長(zhǎng)勝，吳忠良，馬宏生，周龍泉.2009.PI算法用于川滇—安達(dá)曼—蘇門答臘地區(qū)7.0級(jí)以上強(qiáng)震危險(xiǎn)性預(yù)測(cè)的回溯性檢驗(yàn)［J］.地震學(xué)報(bào)，31（3）：307-318.

Jiang C S，Wu Z L，Ma H S，Zhou L Q.2009.Sichuan-Yunnan versus Andaman-Sumatra：PI approach and retrospective forecast test［J］.Acta Seismologica Sinica，31（3）：307-318（in Chinese）.

石耀霖，劉杰，張國(guó)民.2000.對(duì)我國(guó)90年代年度地震預(yù)報(bào)的評(píng)估［J］.中國(guó)科學(xué)院研究生院學(xué)報(bào)，17（1）：63-69.

Shi Y L，Liu J，Zhang G M.2000.The evaluation of Chinese annual earthquake prediction in the 90s［J］.Journal of the Graduate School Academia Sinica，17（1）：63-69（in Chinese）.

王曉欣，劉瑞豐，鄒立曄，袁菲，陳宏峰，陳運(yùn)泰，張立文，梁建宏.2014.2014年2月12日新疆于田MS7.3地震參數(shù)測(cè)定［J］.地震學(xué)報(bào)，36（3）：522-528.

Wang X X，Liu R F，Zou L Y，Yuan F，Chen H F，Chen Y T，Zhang L W，Liang J H.2014.Determination of the parameters for the 12February 2014Yutian MS7.3earthquake in Xinjiang Uygur Autonomous Region［J］.Acta Seismologica Sinica，36（3）：522-528（in Chinese）.

徐錫偉，張培震，聞學(xué)澤，秦尊麗，陳桂華，朱艾斕.2005.川西及其鄰近地區(qū)活動(dòng)構(gòu)造基本特征與強(qiáng)震復(fù)發(fā)模型［J］.地震地質(zhì)，27（3）：446-461.

Xu X W，Zhang P Z，Wen X Z，Qin Z L，Chen G H，Zhu A L.2005.Features of active tectonics and recurrence behaviours of strong earthquakes in the western Sichuan Province and its adjacent regions［J］.Seismology and Geology，27（3）：446-461（in Chinese）.

尹光華，蔣靖祥，吳國(guó)棟.2008.2008年3月21日于田7.4級(jí)地震的構(gòu)造背景［J］.干旱區(qū)地理，31（4）：543-549.

Yin G H，Jiang J X，Wu G D.2008.Tectonic background of the MS7.4earthquake at Yutian on March 21，2008［J］.Arid Land Geography，31（4）：543-549（in Chinese）.

張培震.1999.中國(guó)大陸巖石圈最新構(gòu)造變動(dòng)與地震災(zāi)害［J］.第四紀(jì)研究，19（5）：404-413.

Zhang P Z.1999.Late Quaternary tectonic deformation and earthquake hazard in continental China［J］.Quaternary Sciences，19（5）：404-413（in Chinese）.

張小濤.2009.圖像信息（PI）方法在我國(guó)大陸部分地區(qū)的應(yīng)用研究［D］.北京：中國(guó)地震局地震預(yù)測(cè)研究所：1-64.

Zhang X T.2009.Study on Pattern Informatics Method to Apply on Some Regions of China［D］.Beijing：Institute of Earthquake Science，China Earthquake Administration：1-64（in Chinese）.

張小濤，張永仙，夏彩韻，吳永加，余懷忠.2014.利用圖像信息方法研究蘆山MS7.0地震前川滇及附近地區(qū)的圖像異常［J］.地震學(xué)報(bào)，36（5）：780-789.

Zhang X T，Zhang Y X，Xia C Y，Wu Y J，Yu H Z.2014.Anomalous seismic activities in the Sichuan-Yunnan region and its adjacent areas before the Lushan MS7.0earthquake by the pattern informatics method［J］.Acta Seismologica Sinica，36（5）：780-789（in Chinese）.

Chen C C，Rundle J B，Holliday J R，Nanjo K Z，Turcotte D L，Li S C，Tiampo K F.2005.The 1999Chi-Chi，Taiwan，earthquake as a typical example of seismic activation and quiescence［J］.Geophys Res Lett，32（22）：L22315.doi：10.1029／2005GL023991.

Holliday J R，Nanjo K Z，Tiampo K F，Rundle J B，Turcotte D L.2005.Earthquake forecasting and its verification［J］.Nonlin Processes Geophys，12（6）：965-977.

Molchan G M.1997.Earthquake prediction as a decision-making problem［J］.Pure Appl Geophys，149（1）：233-247.

Nanjo K Z，Rundle J B，Holliday J R，Turcotte D L.2006.Pattern informatics and its application for optimal forecasting of large earthquakes in Japan［J］.Pure Appl Geophys，163（11／12）：2417-2432.

Rundle J B，Klein W，Gross S J，Tiampo K F.2000a.Dynamics of seismicity patterns in systems of earthquake faults［C］∥Geocomplexity and the Physics of Earthquakes：Vol.120 of Geophys Monogr Ser.Washington D C：AGU：127-146.

Rundle J B，Klein W，Turcotte D L，Malamud B D.2000b.Precursory seismic activation and critical-point phenomena［J］.Pure Appl Geophys，157（11／12）：2165-2182.

Rundle J B，Tiampo K F，Klein W，Martins J S S.2002.Self-organization in leaky threshold systems：The influence of near-mean field dynamics and its implications for earthquakes，neuro-biology，and forecasting［J］.Proc Natl Acad Sci USA，99：2514-2521.

Swets J A.1973.The relative operating characteristic in psychology：A technique for isolating effects of response bias finds wide use in the study of perception and cognition［J］.Science，182（4116）：990-1000.

Tiampo K F，Rundle J B，McGinnis S A，Klein W.2002.Pattern dynamics and forecast methods in seismically active regions［J］.Pure Appl Geophys，159（10）：2429-2467.

Wolfram S.2002.A New Kind of Science［M］.Champaign，IL：Wolfram Media，Inc，1197pp.

Xu S X.1989.Mark evaluation for earthquake prediction efficacy［G］∥Collected Papers of Research on Practical Methods of Earthquake Prediction：Volume of Seismology.Beijing：Academic Books and Periodical Press：586-590（in Chinese）.

Zhang Y X，Zhang X T，Wu Y J，Yin X C.2013.Retrospective study on the predictability of pattern informatics to the Wenchuan M8.0and Yutian M7.3earthquakes［J］.Pure Appl Geophys，170（1／2）：197-208.