宋程 張永仙 周少輝 畢金孟 徐小遠

摘要：應(yīng)用圖像信息（PI）方法研究2021年瑪多MS7.4地震前后的熱點特征變化。以8 a預(yù)測窗長、0.8°×0.8°網(wǎng)格尺度為主要參數(shù)，獲取2009—2028年逐年滑動的預(yù)測窗熱點分布圖像。結(jié)果顯示：① 在2009—2020年的5個連續(xù)回溯性預(yù)測窗口中，瑪多地震的震中所在網(wǎng)格及其摩爾鄰近網(wǎng)格持續(xù)存在密集熱點，對發(fā)震地點指示作用較強。② 2014—2028年，多數(shù)預(yù)測窗口內(nèi)瑪多地震震中或摩爾鄰近網(wǎng)格出現(xiàn)過PI熱點；在最后3個連續(xù)預(yù)測窗口內(nèi)（2019—2026年、2020—2027年、2021—2028年），震中附近出現(xiàn)密集深色熱點且呈逐漸收縮的趨勢，對發(fā)震緊迫性具有一定的指示意義。③ 綜合2009年以來的熱點演化圖像可知，巴顏喀拉地塊、柴達木地塊、祁連地塊和羌塘地塊內(nèi)都出現(xiàn)過熱點，出現(xiàn)在巴顏喀拉地塊與柴達木地塊交界瑪多縣附近的熱點頻次高、顏色深、覆蓋面較大，表明青藏塊體內(nèi)部巴顏喀拉地塊中北緣地震活動性較強、相對發(fā)震風(fēng)險較高。

關(guān)鍵詞：瑪多MS7.4地震；PI方法；熱點；回溯性預(yù)測

中圖分類號：P315.7?? 文獻標(biāo)識碼：A?? 文章編號：1000-0666（2023）02-0226-11

doi：10.20015/j.cnki.ISSN1000-0666.2023.0027

0 引言

2021年5月22日在青海省東部果洛州瑪多縣（34.59°N，98.34°E）發(fā)生MS7.4地震，震源深度17 km。瑪多MS7.4地震震中位于青藏塊體內(nèi)次級塊體巴顏喀拉地塊內(nèi)部，距離塊體北邊界東昆侖斷裂帶以南約70 km（岳沖等，2021；詹艷等，2021）。20世紀(jì)末以來，中國大陸7級以上地震主要分布在巴顏喀拉地塊周緣，如1997年西藏瑪尼MS7.9地震、2001年青海昆侖山口西MS8.1地震、2008年四川汶川MS8.0地震、2010年青海玉樹MS7.1地震、2013年四川蘆山MS7.0地震和2017年四川九寨溝MS7.0地震。自印度—歐亞板塊新生代碰撞以來，青藏高原及其周緣山地不斷隆升，構(gòu)造變形及應(yīng)力傳遞對整個東亞地區(qū)的構(gòu)造也都帶來了深遠的影響（Yin，Harrison，2000）。其中青藏高原中北緣巴顏喀拉地塊及周緣地區(qū)構(gòu)造變形強烈（李海兵等，2021），除其邊界斷裂的持續(xù)變形外，其內(nèi)部斷裂系統(tǒng)也較為復(fù)雜，形成了一個天然的構(gòu)造地震研究試驗場。地震目錄是跟蹤統(tǒng)計及研究天然地震的基礎(chǔ)數(shù)據(jù)資料，前人曾結(jié)合地震目錄將圖像信息方法（Pattern informatics method，簡稱PI方法）應(yīng)用于青海地區(qū)及青藏塊體的震例回溯研究（Zhang et al，2013，2017；張小濤等，2014；夏彩韻等，2015；袁伏全等，2016；余娜等，2021）。

PI方法從時、空、強三要素上將一定區(qū)域內(nèi)地震目錄建模計算后以“熱點”的形式呈現(xiàn)，是基于統(tǒng)計物理的地震預(yù)測方法（Holliday et al，2007；蔣長勝等，2009；蔣卉等，2013；Zhang et al，2013；Kawamura et al，2014；張小濤等，2014）。結(jié)合地質(zhì)構(gòu)造及斷裂分布情況，分析熱點分布、演化、遷移趨勢等，可綜合判斷未來一定時段內(nèi)可能的發(fā)震位置（Kawamura et al，2013，2014；Zhang et al，2013，2017）。近20年來，PI方法經(jīng)過國內(nèi)外眾多學(xué)者在全球不同地域的應(yīng)用和檢驗，充分驗證了其良好的預(yù)測效果，如日本局部地區(qū)（Nanjo et al，2006a，b；Kawamura et al，2013，2014；宋程等，2017，2018），印度尼西亞爪哇島地區(qū)（Wibowo et al，2017），中國川滇地區(qū)（Jiang，Wu，2010；蔣卉等，2013；Zhang et al，2013；張小濤等，2014）、華北地區(qū)（張小濤，2009；孫麗娜等，2012）、青藏地區(qū)（Zhang et al，2013，2017；余娜等，2021）以及中國臺灣地區(qū)（Chen et al，2005；Chang et al，2016）。絕大多數(shù)回溯性震例檢驗結(jié)果顯示，大地震前震中及其附近存在PI熱點。無論是板間地震還是板內(nèi)地震，在前人的研究中都驗證了PI方法熱點預(yù)測的效果明顯優(yōu)于隨機預(yù)測（張小濤等，2014；夏彩韻等，2015；尼魯帕爾· 買買吐孫等，2018；宋程等，2018；余娜等，2021）。

針對不同的目標(biāo)地震進行研究時，考慮到孕震環(huán)境的復(fù)雜性和地震活動性起伏等不同原因，應(yīng)選擇相對合理的模型參數(shù)，才能體現(xiàn)出更好的熱點效果（Zhang et al，2013；Kawamura et al，2013，2014；張小濤等，2014；宋程等，2018）。在“大震”未發(fā)生之前，充分利用研究區(qū)內(nèi)相對較小的地震“預(yù)判”后續(xù)中強震可能的發(fā)震情況，不失為一種震前提取異常的中長期方法。在前人的研究中，大面積的熱點很可能覆蓋了主震破裂區(qū)（蔣長勝，2010），這對于震前應(yīng)急防御、震后應(yīng)急救援具有一定的輔助作用。本文以2021年5月22日青?，敹郙S7.4地震為主要研究對象，使用PI方法對瑪多地震前后進行多窗口PI熱點特征研究，通過不同預(yù)測窗時段內(nèi)地震熱點的變化特征，探索識別并捕捉未來可能發(fā)生強震的異常信息。

1 方法簡介及參數(shù)設(shè)置

1.1 計算方法簡介

前人的研究中對PI 方法已有詳細介紹，本文只簡述其計算方法。主要分為4個步驟：①首先對選定的研究區(qū)域進行時空網(wǎng)格化，參與計算的目錄數(shù)據(jù)其截止震級MC在不低于最小完備性震級的同時至少要小于“目標(biāo)地震震級MT”兩個震級單位（Holliday et al，2005；Zhang et al，2013）。②對落入相應(yīng)網(wǎng)格內(nèi)的地震活動頻度構(gòu)建其時間序列，并定義為3個時段（圖1）。通常設(shè)定研究資料起始時刻為t0；滑動變化的時間序列起始時刻為tb（tb≥t0），tb從t0時刻向［t1-（t2-t1）］時刻滑動（tb

1.2 研究區(qū)及模型參數(shù)選取

針對2021年青?，敹郙S7.4地震，選?。?0°～40°N，94°～106°E）為研究區(qū)，并從中國地震臺網(wǎng)中心獲取1980年1月1日至2021年5月22日瑪多地震前研究區(qū)的地震目錄。近20年來中國大陸7級以上地震主要發(fā)生在巴顏喀拉地塊周緣，但瑪多MS7.4地震發(fā)生于巴顏喀拉地塊內(nèi)部，本文主要針對研究區(qū)內(nèi)4個二級塊體（圖2）進行分析討論。由于PI方法計算中的各項參數(shù)對于熱點效果都存在一定的影響，需要根據(jù)震例實際情況對參與計算的震級下限、震源深度、網(wǎng)格尺度、時間窗長等進行選擇設(shè)置（Zhang et al，2013；張小濤等，2014；宋程等，2018；余娜等；2021）。

在袁伏全等（2016）對青海地區(qū)MS≥5.0目標(biāo)地震和余娜等（2021）對青藏塊體（30.0°～41.0°N，88.0°～105.0°E）MS≥6.0目標(biāo)地震的回溯研究中，認(rèn)為在網(wǎng)格尺度0.2°×0.2°、預(yù)測時間窗長3 a、截止震級MC≥ML2.8（1980年以來，滿足不低于最小完備性震級）的參數(shù)組合下熱點效能較好。最小完備性震級可以有效地反映地震臺網(wǎng)的監(jiān)測能力（Schaff，2008；馮建剛等，2012），而輸入目錄數(shù)據(jù)的MC是PI方法計算中直接影響熱點效果的重要參數(shù)。MC偏高會導(dǎo)致輸入數(shù)據(jù)量不足，偏低則會導(dǎo)致輸入量過多降低熱點對于目標(biāo)地震的敏銳性。蔣長勝（2010）研究表明，對于板內(nèi)地震研究的模型參數(shù)選取中設(shè)定MC低于“目標(biāo)地震震級MT”兩個震級單位的常用標(biāo)準(zhǔn)可能會出現(xiàn)問題。而本文針對瑪多MS7.4地震預(yù)設(shè)目標(biāo)地震下限為MT=MS7.0，明顯高于5～6級震例，在保證MC+2≤MT的情況下理論上最高可取MS5.0（Rundle et al，2002；Tiampo et al，2002a；Holliday et al，2005；Nanjo et al，2006a）。為進一步選定合適的MC，分時段畫出1981年至2021年5月21日研究區(qū)的G-R關(guān)系圖（圖3），從圖3中能夠看到4個時段研究區(qū)內(nèi)整體數(shù)據(jù)擬合較好。但圖3b、d中在ML4.8處局部擬合線與整體擬合線的斜率出現(xiàn)偏轉(zhuǎn)，故綜合考慮后本文選擇輸入震級下限為ML4.8，亦可保證數(shù)據(jù)量充足。

前人應(yīng)用PI方法研究中國西部強震用到的網(wǎng)格尺度參數(shù)各有不同，包括青藏地區(qū)1.0°×1.0°（夏彩韻等，2015；Zhang et al，2017）和0.2°×0.2°（袁伏全等，2016；余娜等，2021）、新疆局部地區(qū)0.4°×0.4°（李志海，黃瑜，2010；尼魯帕爾·買買吐孫等，2018）和0.5°×0.5°（尼魯帕爾·買買吐孫，2020）、川滇地區(qū)1.0°×1.0°（蔣長勝等，2009；Zhang et al，2013）和2.0°×2.0°（張小濤等，2014）。參考前人在應(yīng)用PI方法研究中國大陸西部及國內(nèi)外7級以上震例時的參數(shù)設(shè)置（蔣長勝，2010；Zhang et al，2013，2017；蔣卉等，2013；張小濤等，2014；夏彩韻等，2015；宋程等，2018），通常網(wǎng)格尺度較大且預(yù)測窗長相對較長時，熱點效果更為集中和顯著（Zhang et al，2013，2017；張小濤等，2014；宋程等，2018）。蔣長勝（2010）對于2008年5月12日汶川MS8.0地震的PI算法回溯性研究結(jié)果表明，震前PI熱點會因為網(wǎng)格參數(shù)的變化而覆蓋到主震破裂區(qū)，說明網(wǎng)格尺度要因地、因震制宜。在確定目標(biāo)地震震級下限和截止震級后，考慮到瑪多MS7.4地震的位置特殊性、網(wǎng)格內(nèi)計算數(shù)據(jù)量充足性以及熱點遷移性，本文設(shè)定網(wǎng)格尺度為0.8°×0.8°。

瑪多地震距中國大陸上一次7級以上強震（2017年8月8日四川九寨溝MS8.0地震）發(fā)震時間間隔3.4 a，距青海省內(nèi)上一次7級以上地震（2010年4月14日青海玉樹MS7.1地震）發(fā)震時間間隔11.2 a。7級地震孕震機制較為復(fù)雜，結(jié)合之前不同學(xué)者對于國內(nèi)中強震例的研究經(jīng)驗（Zhang et al，2013；張小濤等，2014；夏彩韻等，2015；Zhang et al，2017；宋程等，2018），本文設(shè)定時間窗長為8 a，閾值選擇前人慣用的lg（ΔPi/ΔPmax）=-0.6（Zhang et al，2013；張小濤等，2014；夏彩韻等，2015；宋程等，2018）。

通常選取與研究的目標(biāo)地震震源深度相接近的范圍內(nèi)的事件參與計算，更能符合其構(gòu)造環(huán)境、孕震過程、熱點追蹤的目的和效果?？紤]到本次瑪多地震為淺緣地震（震源深度17 km）以及保證數(shù)據(jù)量的充足，本文選擇60 km以內(nèi)的淺源地震事件參與計算。

考慮到余震也是應(yīng)力釋放的一種表現(xiàn)，刪除余震可能會造成地震目錄一定程度上的畸變（張小濤等，2014），故本文并未將其刪除。本次針對的目標(biāo)地震震級為MT≥MS7.0，參與計算地震的最小震級為ML4.8，刪除目錄中的較小余震對于計算結(jié)果無實質(zhì)性影響，刪除ML>4.8的所有余震則可能無法恰當(dāng)?shù)姆磻?yīng)局部區(qū)域的熱點特征。Tiampo 等（2002b）認(rèn)為余震反映了局部區(qū)域應(yīng)力的高水平釋放；在物理上，較小的地震事件仍能攜帶較為豐富的強震危險性增強的信息（蔣長勝，2010）。

2 PI熱點圖多窗口連續(xù)演化特征分析

按照設(shè)定的時空參數(shù)，從2009年（即汶川MS8.0大地震后）開始以年尺度向未來滑動，得到圖4、5中連續(xù)的PI熱點分布圖，其中包含瑪多地震發(fā)震年份的時間窗有8個（圖5）。摩爾鄰近網(wǎng)格，指緊密包圍在某一網(wǎng)格周圍的8個網(wǎng)格（Moore，1962；Wolfram，2002）。當(dāng)目標(biāo)地震震中所在網(wǎng)格及其摩爾鄰近網(wǎng)格存在熱點時，認(rèn)為目標(biāo)地震可以被PI熱點有效預(yù)測，則針對一個目標(biāo)地震，有效熱點數(shù)最多為9個。

2.1 汶川地震后至瑪多地震前年尺度熱點演化特征分析

圖4顯示2009—2020年研究區(qū)內(nèi)5個連續(xù)預(yù)測窗口中的熱點變化情形，青?，敹嗫h附近持續(xù)出現(xiàn)成片密集熱點。雖然瑪多地震發(fā)生在2021年，但在5個預(yù)測窗口中，瑪多地震震中所在網(wǎng)格及其3個摩爾鄰近網(wǎng)格持續(xù)存在熱點，即這5個回溯性預(yù)測窗口中顯示了明確的高發(fā)震風(fēng)險覆蓋范圍，具有強震發(fā)震危險性地點指示意義。

圖4中分布于巴顏喀拉塊體南邊界的2010年青海玉樹MS7.1地震、東邊界的2013年四川蘆山MS7.0地震和2017年四川九寨溝MS7.0地震并未出現(xiàn)過有效PI熱點。甚至整個研究區(qū)內(nèi)34°N以南地區(qū)幾乎無熱點，只在2011—2018年預(yù)測窗口內(nèi)青海中南部出現(xiàn)過2個熱點。Zhang等（2013）對中國西部大陸（20°～50°N，70°～110°E）進行PI 方法回溯性預(yù)測研究時，發(fā)現(xiàn)2008 年于田 MS7.3地震和汶川MS8.0地震震中存在熱點。張小濤等（2014）對川滇地區(qū)（21°～35°N，97°～109°E）進行 PI 方法熱點回溯性預(yù)測研究得出，蘆山地震震中落入熱點網(wǎng)格內(nèi)。而圖4中蘆山地震和九寨溝地震未出現(xiàn)有效熱點的原因，除了與所選擇的研究區(qū)主體范圍不同（本文選擇研究區(qū)相對較小）外，還和PI方法對于模型參數(shù)的依賴性有關(guān)。張小濤等（2014）對蘆山MS7.0地震前的PI熱點特征研究中設(shè)置的主要參數(shù)為2°×2°/1°×1°網(wǎng)格尺度、9 a預(yù)測時間窗長。Zhang等（2013）針對2008年汶川MS8.0和于田MS7.3地震的PI熱點回溯性研究中，選取的主要參數(shù)為2°×2°/1°×1°網(wǎng)格尺度，預(yù)測窗長為4～8 a的ROC檢驗效果較好。以上研究中截止震級皆為ML4.5，不低于區(qū)域最小完備性震級，但低于本文所選的最小輸入震級ML4.8，數(shù)據(jù)量可能更大。本文所選的預(yù)測窗長為8 a，跨度時長接近或重疊于上述時間參數(shù)。差異最大的參數(shù)為網(wǎng)格尺度，本文設(shè)定為0.8°×0.8°。模型參數(shù)的差異性，可能是本文研究中其他7級以上地震缺少有效PI 熱點的主要原因。

2.2 2014—2028年瑪多MS7.4地震相關(guān)熱點演化特征分析

按照相同的時空參數(shù)，從2014年開始向未來滑動，則包含瑪多地震發(fā)震時間（2021年）的窗口有8個。從時間先后順序來看，8個預(yù)測窗口中瑪多地震的相關(guān)PI熱點呈現(xiàn)為“出現(xiàn)—減弱—消失—強密集—收縮”的趨勢變化。

在前3個窗口中（圖5a～c）熱點相對分散，且主要分布于研究區(qū)的西北部，但瑪多地震震中或摩爾鄰近網(wǎng)格出現(xiàn)了熱點。相對于圖4預(yù)測窗口中部相對穩(wěn)定連續(xù)的密集PI熱點，圖5中瑪多地震附近的熱點分布呈不連續(xù)、不穩(wěn)定、叢集性不明顯特征。在2017—2024年、2018—2025年2個預(yù)測窗口中（圖5d、e），預(yù)測窗口起始時間距離發(fā)震年份為4～5 a，即發(fā)震年份處于窗口近乎中間的年份，瑪多地震及其摩爾鄰近網(wǎng)格未出現(xiàn)相關(guān)熱點，研究區(qū)內(nèi)熱點主要分布在青海最南部，玉樹藏族自治州囊謙縣及玉樹縣（青川藏交界附近）。在 2019—2026 年預(yù)測窗口中（圖5f），瑪多地震震中及其摩爾鄰近網(wǎng)格出現(xiàn)明顯密集深色熱點，且只有2個零散熱點，熱點整體叢集性較強。至最后2個窗口時段（5g、h），瑪多地震摩爾鄰近網(wǎng)格密集熱點顏色幾乎無變化但覆蓋面積有所收縮，彰顯了發(fā)震緊迫性的增強，對瑪多地震的發(fā)震時間及位置具有一定的指示意義。

在研究區(qū)內(nèi)13個連續(xù)滑動的預(yù)測窗熱點圖像中，只有2013—2020年預(yù)測窗口出現(xiàn)了19個熱點，其余12個預(yù)測窗口中熱點總數(shù)不超過15個。在圖5中，只有2個窗口在瑪多地震震中或摩爾鄰近網(wǎng)格皆無熱點，有效熱點占比數(shù)最高的時段為2019—2026年預(yù)測窗。

2.3 熱點遷移分析

從2009年以來研究區(qū)內(nèi)的13個連續(xù)滑動預(yù)測窗熱點圖像（圖4、5）中可見，巴顏喀拉地塊內(nèi)部始終存在PI 熱點，說明其地震活躍程度相對于其他次級地塊更為活躍。其中11個預(yù)測窗中，巴顏喀拉地塊和柴達木地塊交界處存在熱點；在羌塘地塊和祁連地塊內(nèi)，各有8個預(yù)測窗口出現(xiàn)過至少1個熱點。在2009—2020年預(yù)測窗口中，5個滑動窗內(nèi)巴顏喀拉地塊和柴達木地塊交界處青?，敹嗫h附近持續(xù)存在密集熱點（圖4）。在2012—2019年、2013—2020年預(yù)測窗口中，祁連地塊出現(xiàn)了明顯熱點，主要分布在青海海西州德令哈市附近。在2014—2021年、2015—2022年、2016—2023年3個預(yù)測窗口時段內(nèi)（圖5a～c），熱點主要在南北向上跨巴顏喀拉地塊、柴達木地塊和祁連地塊3個二級塊體。

在2014—2021年預(yù)測窗口中（圖5a）無零散熱點，巴顏喀拉地塊內(nèi)瑪多地震震中所在網(wǎng)格熱點顏色最深，發(fā)震危險性最高。從2015—2022年預(yù)測窗口開始至2021—2028年窗口（圖5b～h），在羌塘地塊內(nèi)始終存在熱點。在2016—2023年預(yù)測窗口內(nèi)（圖5c），熱點總數(shù)變少，羌塘地塊以北的3個地塊內(nèi)熱點數(shù)皆有減少?，敹嗟卣饢|側(cè)巴顏喀拉地塊和柴達木地塊交界處熱點變少但顏色加深，危險性增強。至2017—2024年、2018—2025年預(yù)測窗口（圖5d、e），瑪多地震震中及其摩爾鄰近網(wǎng)格不存在有效熱點。在2017—2024年預(yù)測窗口中，柴達木地塊和祁連地塊的最北端熱點延伸到了甘肅境內(nèi)。至2018—2025年預(yù)測窗口中（圖5e），巴顏喀拉地塊以北的熱點全部消失，淺色熱點集中于巴顏喀拉地塊和羌塘地塊交匯處（青海南部），且在巴顏喀拉地塊東北緣（四川北部甘川交界附近）出現(xiàn)了一個深色熱點。

在最后3個連續(xù)的預(yù)測窗口中（圖5f～h），瑪多地震東側(cè)巴顏喀拉地塊（北部）和柴達木地塊的交匯處存在密集深色熱點，且熱點個數(shù)逐漸減少。2019—2026年預(yù)測窗口（圖5f），相對于前2個滑動窗，瑪多地震震中及摩爾鄰近網(wǎng)格出現(xiàn)深色密集熱點，發(fā)震危險性顯著增強。2020—2027年預(yù)測窗口（圖5g），相對于上一個窗口熱點總個數(shù)不變，但零散熱點危險性降低。巴顏喀拉地塊和柴達木地塊的交匯處的密集熱點危險性增強，瑪多地震震中熱點消失。至2021—2028年預(yù)測窗口（圖5h），巴顏喀拉地塊和柴達木地塊的交匯處密集熱點進一步收縮，體現(xiàn)了后續(xù)發(fā)生強震的危險性在向瑪多地震的東側(cè)遷移，并且可能和震源區(qū)及其附近仍然主要受到NEE向近水平擠壓應(yīng)力場的控制有關(guān)（呂苗苗等，2022）。

青藏高原具有持續(xù)變形的特點（許志琴等，2006），熱點的遷移也是相對變化的，本文研究中沒有任意2個預(yù)測窗口內(nèi)熱點的分布及發(fā)震概率是完全相同的。從2014—2021年預(yù)測窗口開始至2021—2028年預(yù)測窗口，瑪多地震的有效熱點的個數(shù)、分布位置和發(fā)震危險性并不具有穩(wěn)定性、持續(xù)性，這體現(xiàn)了7級以上強震孕震物理過程的復(fù)雜性。研究區(qū)內(nèi)所有熱點在不同次級塊體內(nèi)的遷移也不具有明確的規(guī)律性，其內(nèi)部地震的發(fā)生和區(qū)域構(gòu)造應(yīng)力背景與次級活動地塊息息相關(guān)。徐志國等（2021）通過應(yīng)力場反演結(jié)果得出，研究區(qū)整體的背景應(yīng)力場具有近EW向擠壓應(yīng)力場特征，與巴顏喀拉塊體的水平最大主應(yīng)力方向總體上為EW向（范桃園等，2013）具有較好的一致性。區(qū)域構(gòu)造應(yīng)力場能夠影響震源區(qū)的變形和破裂，本次瑪多MS7.4地震為雙側(cè)破裂（李智敏等，2021；王未來等，2021），但圖5不同預(yù)測窗口中瑪多地震的西側(cè)并沒有顯著熱點分布。中國大陸強震孕育發(fā)生的動力學(xué)環(huán)境整體呈“縱向分層、橫向分塊”的特征（邵志剛等，2022），塊體之間也存在一定的相互作用，而塊體運動受深部物質(zhì)作用驅(qū)動。大型走滑型地震發(fā)生時，復(fù)雜的平移錯動引起的淺部破裂面和深部破裂面是存在差異的（徐志國等，2021），其深部物質(zhì)運移和深部地下結(jié)構(gòu)異常會對后續(xù)地震產(chǎn)生影響（許志琴等，2006；王椿鏞等，2008；嘉世旭等，2017；宋向輝等，2021）。用PI方法計算得到的熱點對局部區(qū)域某個時段的發(fā)震范圍具有參考作用，但不能精準(zhǔn)預(yù)測發(fā)震地點和震源深度。研究區(qū)內(nèi)的熱點分布與遷移和淺層破裂展布并不完全對應(yīng)，可能和強震孕育過程、深部動力、地下深部結(jié)構(gòu)等因素存在一定關(guān)系。

3 討論

在本次研究中，我們注意到3個值得思考的問題。第一，雖然瑪多地震震級高達MS7.4，但在2014年以來向未來滑動的連續(xù)預(yù)測窗內(nèi)，并不是每個預(yù)測窗圖像都能提取到與瑪多地震相關(guān)的有效PI熱點信息。這不僅和PI方法對于數(shù)據(jù)的依賴性以及模型參數(shù)的選取有關(guān)，更與青藏塊體構(gòu)造環(huán)境特點和7級強震的孕震復(fù)雜性有關(guān)。第二，在回溯性預(yù)測PI熱點圖像中，雖然瑪多地震附近存在密集PI熱點，但瑪多地震的震中并不位于成片密集熱點的“中心”，這個結(jié)果顯示預(yù)測目標(biāo)地震發(fā)震地點及破裂區(qū)范圍存在困難。多數(shù)學(xué)者認(rèn)為左旋走滑的昆侖山口—江措斷裂為瑪多MS7.4地震的發(fā)震斷裂（李經(jīng)緯等，2021；李智敏等，2021；潘家偉等，2021；宋向輝等，2021；談洪波等，2021；王未來等，2021；姚生海等，2021；尹欣欣等，2021；岳沖等，2021；趙韜等，2021），瑪多MS7.4地震為雙側(cè)破裂（李智敏等，2021），但其相關(guān)熱點主要分布于震中及東側(cè)，未能完全覆蓋住破裂區(qū)及余震區(qū)。第三，在不包含瑪多地震發(fā)震年份的回溯性預(yù)測窗口內(nèi)，瑪多地震震中及鄰近網(wǎng)格存在熱點，但在包含目標(biāo)地震年份的預(yù)測窗口中，有效熱點是間歇性存在的，這進一步增加了判定發(fā)震年份的難度。

印度板塊與歐亞板塊碰撞后持續(xù)向北推擠，在青藏塊體內(nèi)部形成了若干構(gòu)造塊體分別沿大型走滑斷裂帶向東和南東方向“逃逸”的構(gòu)造環(huán)境（陳兵等，2003）。作為晚新生代以來從青藏高原朝東和東南方向“逃逸”的活動地塊之一（聞學(xué)澤，2018），巴顏喀拉塊體東邊界帶及其鄰近地區(qū)的活動構(gòu)造作用以及地震動力源，均來自該塊體的“逃逸”運動（Zhang et al，2004；聞學(xué)澤等，2011；Zhang，2013）。2021年瑪多MS7.4地震是巴顏喀拉地塊最新活動的體現(xiàn)，表明中國大陸7級地震主體活躍區(qū)域依然是巴顏喀拉塊體（姚生海等，2021），亦表明巴顏喀拉塊體具有向東走滑的特性且東向運動仍在持續(xù)中（Xu et al，2013；聞學(xué)澤，2018；尹海權(quán)等，2020），這可能正是熱點在后續(xù)的“未來窗口”中向東收縮的內(nèi)在原因。

4 結(jié)論

本文利用PI方法，參考前人對中強震震例的研究經(jīng)驗，對2021年瑪多MS7.4地震進行回溯性預(yù)測研究，選取計算參數(shù)為空間網(wǎng)格0.8°×0.8°，預(yù)測窗長8 a，震源深度60 km以內(nèi)ML≥4.8的地震事件參與PI方法計算，主要得到以下結(jié)論：

（1）在瑪多地震發(fā)生前的5個預(yù)測時間窗中，瑪多地震震中所在網(wǎng)格及摩爾鄰近網(wǎng)格持續(xù)存在熱點，且在巴顏喀拉地塊和柴達木地塊交界處始終存在成片密集熱點，這對發(fā)震地點指示作用較強。

（2）在包含2021年的8個預(yù)測窗熱點圖像中，瑪多地震震中及摩爾鄰近網(wǎng)格熱點形態(tài)變化并不穩(wěn)定，出現(xiàn)了斷續(xù)起伏的情形。

（3）在2018年之后的3個預(yù)測窗口中，瑪多地震震中附近熱點顏色逐漸加深且出現(xiàn)了向東收縮的趨勢，這對瑪多縣附近的發(fā)震危險性具有一定的指示意義。

（4）在2014—2021年、2015—2022年、2016—2023年3個窗口中，熱點主要分布在巴顏喀拉地塊、柴達木地塊和祁連地塊內(nèi)；在2019—2026年、2020—2027年和2021—2028年預(yù)測窗口中，深色密集熱點主要分布在巴顏喀拉地塊和柴達木地塊的交界處。

中國地震臺網(wǎng)中心提供地震目錄數(shù)據(jù)，本文圖件使用Matlab 和GMT軟件繪制，在此一并表示感謝！

參考文獻：

陳兵，江在森，車時，等.2003.瑪尼7.9級地震對昆侖山口西8.1級地震的觸發(fā)作用及動力背景初探［J］.中國地震，19（1）：1-7.

范桃園，陳群策，吳中海，等.2013.青藏高原東緣活動構(gòu)造與現(xiàn)今地應(yīng)力場三維粘彈性模擬研究［J］.地球物理學(xué)進展，28（3）：1140-1149.

馮建剛，蔣長勝，韓立波，等.2012.甘肅測震臺網(wǎng)監(jiān)測能力及地震目錄完整性分析［J］.地震學(xué)報，34（5）：646-658.

嘉世旭，林吉焱，郭文斌，等.2017.巴顏喀拉塊體地殼結(jié)構(gòu)多樣性探測［J］.地球物理學(xué)報，60（6）：2226-2238.

蔣長勝，吳忠良，馬宏生，等.2009.PI 算法用于川滇—安達曼—蘇門答臘地區(qū)7.0級以上強震危險性預(yù)測的回溯性檢驗［J］.地震學(xué)報，31（3）：307-318.

蔣長勝.2010.川滇地區(qū)地震活動的統(tǒng)計物理特征［J］.國際地震動態(tài)，（2）：44-45.

蔣卉，吳忠良，馬騰飛，等.2013.對圖像信息學(xué)（PI）算法的一個回溯性預(yù)測檢驗：四川蘆山7.0級地震［J］.物理，42（5）：334-340.

李海兵，潘家偉，孫知明，等.2021.大陸構(gòu)造變形與地震活動——以青藏高原為例［J］.地質(zhì)學(xué)報，95（1）：194-213.

李經(jīng)緯，陳長云，占偉，等.2021.青?，敹?.4 級地震 GNSS 同震水平位移的快速獲?。跩］.地震地質(zhì)，43（5）：1073-1084.

李志海，黃瑜.2010.興都庫什—帕米爾—新疆中南部地區(qū)PI算法的回溯性研究［J］.內(nèi)陸地震，24（3）：208-214.

李智敏，李文巧，李濤，等.2021.2021年5月22日青海瑪多MS7.4地震的發(fā)震構(gòu)造和地表破裂初步調(diào)查［J］.地震地質(zhì)，43（3）：722-737.

呂苗苗，常利軍，魯來玉，等.2022.2021年青海瑪多MS7.4地震余震序列震源機制解及其發(fā)震構(gòu)造特征［J］.地球物理學(xué)報，65（6）：1991-2005.

尼魯帕爾·買買吐孫，余懷忠.2020.利用圖像信息方法研究庫車MS5.6、于田MS6.4地震前異常演化［J］.內(nèi)陸地震，135（3）：215-222.

尼魯帕爾·買買吐孫，張小濤，張永仙，等.2018.PI算法用在天山地震帶的回溯性研究［J］.地震工程學(xué)報，40（S1）：45-52.

潘家偉，白明坤，李超，等.2021.2021年5月22日青?，敹郙S7.4地震地表破裂帶及發(fā)震構(gòu)造［J］.地質(zhì)學(xué)報，95（6）：1655-1670.

邵志剛，王武星，劉琦，等.2022.活動地塊理論框架下的地震物理預(yù)報展望［J］.科學(xué)通報，67（13）：1362-1377.

宋程，張永仙，夏彩韻，等.2017.基于圖像信息方法的日本東北MW9.0大地震回溯性預(yù)測研究［J］.地震，37（2）：47-56.

宋程，張永仙，夏彩韻，等.2018.日本東北MW9.0 地震的PI模型參數(shù)設(shè)置與預(yù)測效能回溯性檢驗［J］.地震學(xué)報，40（4）：491-505.

宋向輝，王帥軍，潘素珍，等.2021.2021年瑪多MS7.4地震的深部構(gòu)造背景［J］.地震地質(zhì)，43（4）：757-770.

孫麗娜，齊玉妍，溫超，等.2012.中長期地震預(yù)測中的PI算法改進研究及應(yīng)用［J］.地震，32（4）：44-52.

談洪波，王嘉沛，楊光亮，等.2021.2021年瑪多MS7.4地震的震后效應(yīng)模擬［J］.地震地質(zhì)，43（4）：936-957.

王椿鏞，樓海，呂智勇，等.2008.青藏高原東部地殼上地幔S波速速度結(jié)構(gòu)———下地殼流的深部環(huán)境［J］.中國科學(xué)：地球科學(xué)，38（1）：22-32.

王未來，房立華，吳建平，等.2021.2021年青?，敹郙S7.4地震序列精定位研究［J］.中國科學(xué)：地球科學(xué)，51（7）：1193-120.

聞學(xué)澤，杜方，張培震，等.2011.巴顏喀拉塊體北和東邊界大地震序列的關(guān)聯(lián)性與 2008 年汶川地震［J］.地球物理學(xué)報，54（3）：706-716.

聞學(xué)澤.2018.巴顏喀拉塊體東邊界千年破裂歷史與2008年汶川、2013年蘆山和2017年九寨溝地震［J］.地震學(xué)報，40（3）：255-267.

夏彩韻，張永仙，張小濤，等.2015.利用兩次新疆于田7.3級地震對圖像信息（PI）方法可預(yù)測性的檢驗［J］.地震學(xué)報，37（2）：312-322.

徐志國，梁姍姍，張廣偉，等.2021.2021年5月22日青海瑪多MS7.4地震發(fā)震構(gòu)造分析［J］.地球物報學(xué)報，64（8）：2657-2670.

許志琴，楊經(jīng)綏，李海兵，等.2006.青藏高原與大陸動力學(xué)：地體拼合、碰撞造山及高原隆升的深部驅(qū)動力［J］.中國地質(zhì)，33（2）：221-238.

姚生海，蓋海龍，殷翔，等.2021.青海瑪多 MS7.4 地震地表破裂帶的基本特征和典型現(xiàn)象［J］.地震地質(zhì)，43（5）：1060-1072.

尹海權(quán)，郭祥云，常明，等.2020.跨斷層與 GPS 地殼形變數(shù)據(jù)聯(lián)合反演鮮水河斷裂地震危險性［J］.地質(zhì)學(xué)報，94（8）：2487-2499.

尹欣欣，王維歡，蔡潤，等.2021.2021年青?，敹郙S7.4地震精定位和發(fā)震構(gòu)造初探［J］.地震工程學(xué)報，43（4）：834-839.

余娜，張曉清，袁伏全，等.2021.基于PI算法的青藏塊體強震回溯性研究［J］.地震地磁觀測與研究，42（1）：35-41.

袁伏全，孫世瑞，王小玲.2016.PI算法用于青海地區(qū)中強震危險性預(yù)測的回溯性檢驗研究［J］.地震研究，39（S1）：76-82.

岳沖，屈春燕，牛安福，等.2021.瑪多MS7.4地震對周邊斷層的應(yīng)力影響分析［J］.地震地質(zhì)，43（5）：1041-1059.

詹艷，梁明劍，孫翔宇，等.2021.2021年5月22日青?，敹郙S7.4地震深部環(huán)境及發(fā)震構(gòu)造模式［J］.地球物理學(xué)報，64（7）：2232-2252.

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

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

趙韜，王瑩，馬冀，等.2021.2021年青海瑪多7.4級地震序列重定位和震源機制特征［J］.地震地質(zhì)，43（4）：790-805.

Chang L Y，Chen C C，Wu Y H，et al.2016.A strategy for a routine pattern informatics operation applied to Taiwan［J］.Pure and Applied Geophysics，173（1）：235-244.

Chen C C，Rundle J B，Holliday J R，et al.2005.The 1999 Chi-Chi，Taiwan，earthquake as a typical example of seismic activation and quiescence［J］.Geophysical Research Letters，32（22）：L22315.

Holliday J R，Chen C C，Tiampo K F，et al.2007.A RELM earthquake forecast based on pattern informatics［J］.Seismological Research Letters，78（1）：87-93.

Holliday J R，Nanjo K Z，Tiampo K F，et al.2005.Earthquake forecasting and its verification［J］.Nonlinear Processes in Geophysics，12（6）：965-977.

Jiang C S，Wu Z L.2010.PI forecast for the Sichuan-Yunnan region：Retrospective test after the May 12，2008，Wenchuan earthquake［J］.Pure and Applied Geophysics，167（6/7）：751-761.

Kawamura M，Wu Y H，Kudo T，et al.2013.Precursory migration of anomalous seismic activity revealed by the pattern informatics method：A case study of the 2011 Tohoku earthquake，Japan［J］.Bulletin of the Seismological Society of America，103（2B）：1171-1180.

Kawamura M，Wu Y H，Kudo T，et al.2014.A statistical feature of anomalous seismic activity prior to large shallow earthquakes in Japan revealed by the pattern informatics method［J］.Natural Hazards and Earth System Sciences，14（4）：849.

Mohanty W K，Mohapatra A K，Verma A K，et al.2016.Earthquake forecasting and its verification in northeast India［J］.Geomatics，Natural Hazards and Risk，7（1）：194-214.

Moore E F.1962.Machine models of self-reproduction［C］//Proceedings of the Fourteenth Symposium on Applied Mathematics，American Mathematical Society，14：17-33.

Nanjo K Z，Holliday J R，Chen C C，et al.2006b.Application of a modified pattern informatics method to forecasting the locations of future large earthquakes in the central Japan［J］.Tectonophysics，424（3）：351-366.

Nanjo K Z，Rundle J B，Holliday J R，et al.2006a.Pattern informatics and its application for optimal forecasting of large earthquakes in Japan［J］.Pure and Applied Geophysics，163（11-12）：2417-2432.

Radan M Y，Hamzehloo H，Peresan A，et al.2013.Assessing performances of pattern informatics method：A retrospective analysis for Iran and Italy［J］.Natural Hazards，68（2）：855-881.

Rundle J B，Tiampo K F，Klein W，et al.2002.Self-organization in leaky threshold systems：The influence of near-mean field dynamics and its implications for earthquakes，neurobiology，and forecasting［C］.Proceedings of the National Academy of Sciences，2514-2521.

Schaff D P.2008.Semiempirical statistics of correlation-detector performance［J］.Bullentin of the Seismological Society of America，98（3）：1495-1507.

Tiampo K F，Rundle J B，McGinnis S A，et al.2002a.Pattern dynamics and forecast methods in seismically active regions［J］.Pure and Applied Geophysics，159（10）：2429-2467.

Tiampo K F，Rundle J B，McGinnis S A，et al.2002b.Mean-field threshold systems and phase dynamics：An application to earthquake fault systems［J］.Europhysics Letters，60（3）：481.

Wibowo A，Insani A，Nurdiyanto S.2017.Modifikasi pattern informatics untuk prediksi hotspot aktivitas seismik pada gempa di pulau jawa［J］.Jurnal Nasional Teknik Elektro dan Teknologi Informasi（JNTETI），6（2）：121-126.

Wolfram S.2002.A new kind of science［M］.Champaign IL：Wolfram Media，Inc，1197.

Xu X W，Tan X B，Yu G H，et al.2013.Normal- and oblique-slip of the 2008 Yutian earthquake：Evidence for eastward block motion，northern Tibetan Plateau［J］.Tectonophysics，584：152-165.

Yin A，Harrison T M.2000.Geologic evolution of the Himalayan-Tibetan Orogen［J］.Annual Review of Earth and Planetary Sciences，28：211-280.

Zhang P Z，Shen Z K，Wang M，et al.2004.Continuous deformation of the Tibetan Plateau from global positioning system data［J］.Geology，32（9）：809-812.

Zhang P Z.2013.A review on active tectonics and deep crustal processes of the Western Sichuan Region，eastern margin of the Tibetan Plateau［J］.Tectonophysics，584：7-22.

Zhang Y X，Xia C Y，Song C，et al.2017.Test of the predictability of the PI method for recent Large earthquakes in and near the Tibetan Plateau［J］.Pure and Applied Geophysics，174：2411-2426.

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

Retrospective Study on the Forecast of the 2021 MaduoMS7.4 Earthquake by PI Method

SONG Cheng1，ZHANG Yongxian2，ZHOU Shaohui3，BI Jinmeng1，XU Xiaoyuan1

（1.Tianjin Earthquake Agency，Tianjin 300201，China）

（2.Institute of Earthquake Forecasting，China Earthquake Administration，Beijing 100036，China）

（3.Shandong Earthquake Agency，Jinan 250014，Shandong，China）

Abstract

In this paper，the Pattern Informatics（PI）method was applied to the retrospective study on the forecast of the 2021 Maduo MS7.4 earthquake.The local area（30.0°N-40.0°N，94.0°E-106.0°E）of Qinghai-Tibet tectonic block was chosen as the study region.With parameters of the 8-year-long forecasting window and the 0.8°× 0.8°-sized grid，successive forecasting hotspot diagrams were obtained in the period from 2009 to 2028.The results showed that：①In five successive forecasting windows covering the period from 2009 to 2020，continuous dense hotspots appeared in the grid in which the Maduo epicenter located and in its neighboring Moore grids.These hotspots strongly indicated the location of the Maduo earthquake.②In the period from 2014 to 2028，PI hotspots appeared in the grid in which the Maduo epicenter located or in its neighboring Moore grids in most forecasting windows.In the three continuous forecasting windows after 2018（2019-2026，2020-2027 and 2021-2028），dark hotspots appeared near the epicenter of the Maduo MS7.4 earthquake and showed a gradual tendency of contraction，which indicated the impending occurrence of the Maduo MS7.4 earthquake.③According to the continuous hotspot diagrams since 2009，hotspots appeared in the Bayankala block，the Qaidam block，the Qilian block and the Qiangtang block.Relatively，hotspots appearing at the junction of the Bayankala block and the Qaidam block near Maduo county had higher frequency，darker color，larger coverage.The results indicates that the middle part of the northern margin of the Bayankala block in the Qinghai-Tibet tectonic block has stronger seismicity and higher potential earthquake risk.

Keywords：the Maduo MS7.4 earthquake；the Pattern Informatics（PI） method；hotspot；retrospective forecast

收稿日期：2022-04-07.

基金項目：震情跟蹤定向任務(wù)（2021010104，2021010128）；國家重點研發(fā)計劃“中國地震科學(xué)實驗場的地震可預(yù)測性國際合作研究（2018YFE0109700）；國家自然基金委地震聯(lián)合基金“基于數(shù)值模擬的確定性——概率地震危險性分析方法研究（U2039207）；天津市地震局局內(nèi)重點科研項目（Zd202202，Zd202304）.

第一作者簡介：宋 程（1992-），工程師，主要從事地震預(yù)報分析相關(guān)研究.E-mail：nibushiwode798@163.com.

通訊作者簡介：張永仙（1965-），研究員，主要從事地震預(yù)測理論與方法研究.E-mail：yxzhseis@sina.com.

宋程，張永仙，周少輝，等.2023.2021年瑪多MS7.4地震的PI熱點特征回溯性預(yù)測研究［J］.地震研究，46（2）：226-236，doi：10.20015/j.cnki.ISSN1000-0666.2023.0027.