付虹 劉自鳳 胡小靜 李智蓉 洪敏 彭關(guān)靈 王光明 倪喆

摘要：對云南漾濞MS6.4地震前出現(xiàn)的前震序列、地震學(xué)、定點形變和流體觀測異常進(jìn)行系統(tǒng)梳理和回顧，結(jié)果表明：前震序列較為復(fù)雜，能找到一些地震學(xué)的異常特征，但難以通過這些特征對震后趨勢進(jìn)行正確判斷;定點形變觀測異常出現(xiàn)時間早、異常幅度大，可為震級預(yù)測提供依據(jù);流體觀測異常出現(xiàn)時間短、異常相對集中，為時間預(yù)測提供了更好的支撐;形變和流體異常與亞失穩(wěn)后期階段具有較為相似的時空特征，可能與亞失穩(wěn)過程有關(guān)。前震序列的判斷結(jié)合區(qū)域地震危險性，可能是提高判斷準(zhǔn)確性更好的方法，合理利用地震孕育不同階段的異常特征，是做好綜合預(yù)測需要提倡的方法。

關(guān)鍵詞：漾濞MS6.4地震;前震序列;定點形變;流體異常;綜合預(yù)測

中圖分類號：P315.72?? 文獻(xiàn)標(biāo)識碼：A? ?文章編號：1000-0666（2021）03-0301-08

0 引言

2021年5月21日21時48分34秒，在云南省大理白族自治州漾濞縣發(fā)生MS6.4地震，震中位于（25.67°N，99.87°E），震源深度8 km。這次地震屬于1次典型的前震-主震-余震型地震事件，2021年5月18日至21日 MS6.4主震前發(fā)生了10次MS≥3.0地震，最大前震為主震前27 min發(fā)生的MS5.6地震。

地震短臨預(yù)測至今仍是世界性科學(xué)難題，探索地震預(yù)報的可能性一直是地震工作者追求的目標(biāo)。即使是發(fā)生在同一地點的2次地震，震前的異常也幾乎沒有一模一樣的，由此可見地震前兆異常的復(fù)雜性，這可能也是地震預(yù)測困難的原因之一。對每一次發(fā)生的破壞性強(qiáng)震震前異常和預(yù)測過程進(jìn)行記錄、分析、研究，汲取經(jīng)驗教訓(xùn)，是現(xiàn)階段提高地震預(yù)報水平的重要途徑。云南是中國大陸地震最頻繁的區(qū)域之一，前人對發(fā)生在云南的多個破壞性地震進(jìn)行過綜合研究（石紹先等，1991;陳立德，羅平，1992，1997;付虹等，2007，2015），但每次地震都具有不同的異常特點，為了豐富對已經(jīng)發(fā)生地震孕育過程的認(rèn)識，記錄不同地震的異常特征仍然是必要的。漾濞MS6.4地震前，云南地區(qū)地震學(xué)、流動重力、GNSS、定點形變和地下流體等觀測均有異常顯示，2021年5月18日開始還出現(xiàn)了直接前震，但跟蹤預(yù)測過程中未能正確判斷前震的發(fā)展趨勢。本文真實記錄分析地震預(yù)測過程中的異常和認(rèn)識，總結(jié)本次漾濞地震震前異常特征和跟蹤預(yù)測過程中值得汲取的經(jīng)驗教訓(xùn)，為今后的強(qiáng)震分析預(yù)測提供借鑒。

1 漾濞MS6.4地震前震序列跟蹤分析和地震學(xué)顯著異常

1.1 前震序列跟蹤分析

表1給出了云南地震臺網(wǎng)記錄的漾濞MS6.4地震主震和MS≥3.0前震的主要參數(shù)。從表中可見，主震前有直接前震發(fā)生，從2021年5月18日18時第1個MS≥3.0地震發(fā)生，到21日21時的MS6.4主震發(fā)生，中間間隔39 h。其中，主震前1 h內(nèi)發(fā)生的MS4.2、5.6、4.5前震，因與主震時間相隔太近，在跟蹤過程中沒有來得及分析、判斷，但對18和19日發(fā)生的MS4.2、4.4地震我們都做過研究。

（1）序列參數(shù)：選取1990年以來這2次4級地震周圍50 km內(nèi)的全部4級震群，用測震學(xué)地震

預(yù)報方法（國家地震局監(jiān)測預(yù)報司，1997）推薦的序列判定參數(shù)（U、k、ρ）進(jìn)行計算，將序列中達(dá)到前兆標(biāo)準(zhǔn)的參數(shù)列于表2。由表可見：前9組震群，有1組（3號）因為地震數(shù)目太少，不滿足計算條件，有1組（6號）所有的參數(shù)都未達(dá)到前震震群的判別標(biāo)準(zhǔn)，其余7組序列的U值或ρ值，或U、ρ值達(dá)到了前震的判別標(biāo)準(zhǔn)，而這7組序列中只有2016年2月8日的洱源MS4.5地震后，于2016年5月18日基本在同一震源區(qū)發(fā)生了云龍MS5.0地震，其它序列原震區(qū)都沒有發(fā)生MS≥5地震。但5月18和19日的2次MS4.2、4.4地震，21日在原震區(qū)震級升級到了MS6.4。從表2還可見，該區(qū)其它4級地震符合前震的序列里，除2002年和2007年漾濞MS4.1和巍山MS4.4地震外，其它序列的b值也都相對都較低，低于0.6。因此用b值作為該區(qū)的前震序列判別參數(shù)作預(yù)測也是困難的，虛報較高。

（2）震源機(jī)制一致性：2021年5月18—20日的前震序列中，震源機(jī)制一致性非常好（圖1a），按照前人的研究結(jié)果（陳鶞等，1978，2015），震源機(jī)制一致性好是前震序列的重要特征之一。云南滇西譜振幅相關(guān)系數(shù)的研究結(jié)果顯示（黃浩，付虹，2014），震源機(jī)制一致性好，表明區(qū)域應(yīng)力水平較高，有利于后續(xù)地震的發(fā)生，也有部分主震—余震型地震譜振幅相關(guān)系數(shù)高，震源機(jī)制一致性好，后面沒有發(fā)生更大地震的例子。例如，2017年漾濞MS5.1地震震源機(jī)制一致性非常好，其后沒有發(fā)生更大地震（圖1b）。因此，僅有震源機(jī)制一致性作為前震的判別指標(biāo)，在實際應(yīng)用中也是困難的。

（3）視應(yīng)力：陳學(xué)忠等（2003）給出了視應(yīng)力的計算方法和對多次地震序列的跟蹤分析結(jié)果，認(rèn)為低視應(yīng)力的強(qiáng)震發(fā)生后，震區(qū)發(fā)生較強(qiáng)后續(xù)地震的可能性小，高視應(yīng)力對后續(xù)地震有預(yù)報意義。2013年四川蘆山MS7.0地震前，觀測到在龍門山南段有高視應(yīng)力的地震發(fā)生（易桂喜等，2013），李艷娥等（2012，2017）對其它一些震例的研究結(jié)果也顯示主震前中小地震的視應(yīng)力是逐漸增高的。

5月20日對漾濞MS6.4地震幾次前震的視應(yīng)力計算結(jié)果顯示（圖2），2021年5月18、19日的漾濞多次3、4級地震，視應(yīng)力大部分都位于區(qū)域視應(yīng)力與震級相關(guān)的1倍均方差內(nèi)，沒有出現(xiàn)蘆山地震前視應(yīng)力高、偏離均值的狀態(tài)，因此在這次前震活動中，使用視應(yīng)力無法判別該序列震級會升級。

1.2 地震學(xué)顯著異常

漾濞MS6.4地震發(fā)生在2021年度地震重點危險區(qū)內(nèi)，年度異常的主要判據(jù)有云南省地震局.2020.2021年度云南地震趨勢研究報告.：滇西北地區(qū)的中等地震進(jìn)入活躍時段但缺6級地震、滇西北大理一帶震源機(jī)制一致性好達(dá)到預(yù)報閾值、b值在大理州及附近區(qū)域出現(xiàn)異常等，均顯示該區(qū)存在發(fā)生6級以上地震的危險，但這些依據(jù)在發(fā)震時間上不是太確定。

對時間預(yù)測貢獻(xiàn)最大的地震學(xué)異常，2021年以來滇西北地區(qū)ML≥4.0地震平靜（圖3a）。從2019年11月開始云南地區(qū)ML≥4.0地震，在滇西北的德欽至南華地區(qū)平靜457 d后，被2021年3月1日的洱源ML4.3地震打破。該區(qū)1987年有地震速報目錄以來的統(tǒng)計結(jié)果顯示（圖3b），ML≥4.0地震平靜達(dá)到235 d后，3個月內(nèi)對應(yīng)滇西北地區(qū)發(fā)生ML≥5.5地震的概率可達(dá)9/17，因此3月1日洱源ML4.3地震的發(fā)生，標(biāo)志著該區(qū)發(fā)生MS≥5.5地震的短期危險性增加。同時在漾濞MS6.4地震震中的西南地區(qū)——小滇西一帶出現(xiàn)了ML≥4.0地震叢集活動（圖3a），統(tǒng)計該區(qū)的ML≥4.0地震頻度，震前同樣達(dá)到了預(yù)報閾值（圖3c），以往滇西北地區(qū)達(dá)到閾值也對應(yīng)過發(fā)生在滇西到滇西北的MS≥5.5地震。漾濞MS6.4地震發(fā)生在滇西北地區(qū)的ML≥4.0地震平靜和小滇西地區(qū)的ML≥4.0地震活躍相交匯的??? 地區(qū)，短期階段地震學(xué)異常既滿足平靜指標(biāo)，又滿足活躍指標(biāo)，在平靜和活躍的交界地區(qū)發(fā)生強(qiáng)震的震例在云南并不多見。

2 地球物理定點形變和流體觀測異常

漾濞MS6.4地震前，云南地區(qū)的地球物理定點觀測異常中，流體異常從小滇西到西北地區(qū)最為集中，形變異常主要分布在永勝、彌渡、楚雄到通海、石屏等川滇菱形塊體的西邊界一帶，震前這些異常大部分都得到了顯著的關(guān)注。

2.1 定點形變觀測異常

漾濞MS6.4地震前定點地球物理場觀測異常中，出現(xiàn)時間最早的是形變。石屏水平擺NS向、通海垂直擺NS向、楚雄垂直擺NS向和彌渡水管NS向觀測都于震前1年多開始出現(xiàn)異常，2020年10—11月又分別出現(xiàn)楚雄、云龍鉆孔應(yīng)變和永勝水管觀測異常。從表3可見，出現(xiàn)異常最早的是震中距400 km附近的石屏水平擺NS向和通海垂直擺NS向，200 km內(nèi)的異常除彌渡水管NS向外，其余臺異常出現(xiàn)相對較晚，但所有的形變異常都開始于震前半年。漾濞MS6.4地震前的定點形變觀測異常出現(xiàn)時間早、異常幅度大，大幅度的趨勢異常出現(xiàn)時間還具有準(zhǔn)同步性（圖4），與云南地區(qū)多次6級以上地震前異常相似（付虹，劉麗芳，2003）。

2.2 流體觀測異常

漾濞MS6.4地震前的流體異常相對形變出現(xiàn)時間晚，震前最早出現(xiàn)的是龍陵測點的水氡降低、水溫和流量升高達(dá)到經(jīng)驗預(yù)報閾值，這3項異常開始于2020年8月，其后是洱源、麗江和劍川水溫，于2021年先后出現(xiàn)異常（圖5、表4）。這些異常中，麗江水溫異常是震前發(fā)現(xiàn)的，但當(dāng)時認(rèn)為近幾年每年該井水位斷流后，水溫都會出現(xiàn)下降，是干擾;震后再分析，發(fā)現(xiàn)2021年的水位斷流和水溫下降，較前幾年早3個月，且之前水位斷流引起的水溫下降會很快回升，而2021年3月的水溫下降后，一直維持在低值，與往年有區(qū)別，因此筆者認(rèn)為這是漾濞MS6.4地震的前兆異常。另外有爭議的異常還有劍川水溫，震前懷疑是抽水干擾，震后驗證發(fā)現(xiàn)抽水不會對水溫產(chǎn)生影響趙剛.2021.劍川地?zé)嵊^測情況說明.，因此認(rèn)為該項也可以作為震前異常，并將其收錄。流體異常中水溫異常最多，而且大部分都在3個月內(nèi)，說明水溫異常是短期預(yù)報較好的參量，這在云南2007年寧洱MS6.4地震前也曾出現(xiàn)過（付虹，2007）。

綜上所述，云南漾濞MS6.4地震前，定點形變觀測異常呈現(xiàn)時間早、幅度大的特點，流體異常出現(xiàn)時間短，呈現(xiàn)水溫異常數(shù)量最多（這些水溫異常點也都有水位觀測，但未見水位觀測異常）的特點。

3 結(jié)論和討論

本文對2021年漾濞MS6.4地震前云南地區(qū)地震學(xué)、定點形變和地下流體等觀測異常進(jìn)行了系統(tǒng)分析，獲得如下認(rèn)識：

（1）漾濞MS6.4地震屬于前震-主震-余震型地震。MS6.4地震前發(fā)生了10次MS≥3.0地震，按照地震活動參數(shù)計算，前震18、19日地震序列的b值、U值異常，顯示為前兆震群;震源機(jī)制一

致性也顯示該震群為前震，但視應(yīng)力和該區(qū)背景地震的水平相當(dāng)。震源區(qū)附近區(qū)域曾發(fā)生過多次4級震群序列，這些震群序列中也有地震活動參數(shù)和震源機(jī)制一致性顯示為前震序列，但其后均沒有發(fā)生更大地震，這是在跟蹤過程中，沒有把18、19日的4級地震活動判定為直接前震的主要原因。也說明該區(qū)的4級震群活動具有復(fù)雜性，很難用常規(guī)的地震學(xué)參數(shù)對序列進(jìn)行正確判斷。

（2）之前對地震視應(yīng)力的認(rèn)識是視應(yīng)力高說明區(qū)域應(yīng)力水平高，可以作為后續(xù)有破壞性地震預(yù)測的主要依據(jù)。隨著序列研究樣本的增多，目前已經(jīng)觀測到了部分序列直接前震視應(yīng)力偏低的現(xiàn)象（Chen，Shearer，2013），但震前大部分分析預(yù)報人員未能及時關(guān)注到類似震例，同時也說明視應(yīng)力在地震序列中的變化是復(fù)雜的，還需要不斷認(rèn)識和深入研究。

（3）云南地區(qū)大部分6級地震前的異常特征都是地震活動增強(qiáng)（付虹等，2015），只有滇西北地區(qū)以平靜異常為特征，最顯著的震例是1996年麗江MS7.0地震前4級地震平靜后直接發(fā)生強(qiáng)震（圖3b），因此震前特別注意到了滇西北地區(qū)4級地震平靜現(xiàn)象。震前3個月，小滇西地區(qū)4級地震活動增強(qiáng)非常顯著，4級地震頻度達(dá)到預(yù)報閾值，漾濞MS6.4地震發(fā)生在4級地震活動增強(qiáng)區(qū)與平靜區(qū)的交界區(qū)。因此可以分別用2個區(qū)域的異常特征，尋找交集進(jìn)行地點預(yù)測，這在云南的6級震例中并不多見，豐富了對該區(qū)6級地震活動異常的認(rèn)識。

（4）漾濞MS6.4地震前的定點形變觀測異常出現(xiàn)時間早、異常幅度大，大幅度的趨勢異常在震前1～2年的時間尺度出現(xiàn)。與云南以往6級地震前的異常進(jìn)行對比，2014年魯?shù)镸S6.5、景谷MS6.6地震前的形變異?；旧隙汲霈F(xiàn)在1年時間尺度內(nèi)，而以往1～2年時間尺度的異常，大部分都出現(xiàn)在7級地震前，因此漾濞MS6.4地震前的異常是否得到交代，其后是否還有地震對應(yīng)，值得進(jìn)一步深入研究。

（5）漾濞MS6.4地震前的流體異常，開始時間大部分都集中在半年時間尺度內(nèi)，表明在這個區(qū)域，流體異常作為中短期的判別指標(biāo)，對時間的預(yù)測更有意義。震前云南地區(qū)的流體異常都出現(xiàn)在震中200 km范圍內(nèi)，與定點形變的異常范圍相比，在空間上更為集中。定點形變和流體異常由外圍逐漸收縮至震中區(qū)，可能與亞失穩(wěn)有關(guān)，因為在亞失穩(wěn)態(tài)的后期階段具有相似的收縮特征何昌榮等.2020.地震亞失穩(wěn)階段識別的實驗，理論與野外觀測研究中期報告.

（6）不同學(xué)科異常在地震三要素預(yù)測中可能具有不同作用。漾濞MS6.4地震前，滇西北地區(qū)的4級地震平靜和小滇西地區(qū)的4級地震活躍，可為地點預(yù)測提供判據(jù);形變異常顯示該區(qū)未來1～2年可能發(fā)生6級以上地震，可為震級預(yù)測提供判據(jù);流體異常出現(xiàn)時間短，特別是水溫異常集中在半年時間尺度內(nèi)，為預(yù)測時間的逼近提供了依據(jù)。以上表明在地震孕育不同階段，不同學(xué)科的異常特征是有差異的，把這些差異有機(jī)的結(jié)合，是做好綜合預(yù)報的基礎(chǔ)。

（7）地震學(xué)前震判別參數(shù)用于中、小地震序列判斷時，對后續(xù)是否有更大地震發(fā)生進(jìn)行直接判斷是困難的，但漾濞MS6.4地震的跟蹤經(jīng)驗表明，如果序列發(fā)生的區(qū)域有具備發(fā)生更大地震的地震學(xué)和地球物理觀測異常時，滿足前震序列特征的中、小地震序列，將其判定后續(xù)有更大地震的概率加大，前震序列特征結(jié)合區(qū)域地震學(xué)和地球物理背景異常，可能是目前提高強(qiáng)震前震識別率更有效的方法。

參考文獻(xiàn)：

陳立德，羅平.1992.瀾滄7.6 級地震前兆異常的綜合特征（一）中期異常及其綜合特征[J].地震，（4）：18-27.

陳立德，羅平.1997.1995 年7 月12 日孟連中緬邊界7.3 級地震短臨預(yù)報及前兆異常特征（1）[J].地震研究，20（2）：151-156.

陳學(xué)忠，王小平，王琳英，等.2003.地震視應(yīng)力用于震后趨勢快速判定的可能性[J].國際地震動態(tài)，（7）：1-3.

陳鶞.1978.用震源機(jī)制一致性作為描述地震活動性的新參數(shù)[J].地球物理學(xué)報，21（2）：142-159.

陳鶞，劉杰，楊文.2015.前震序列的圖像特征研究[J].中國地震，31（2）：177-187.

付虹，劉麗芳.2003.云南地區(qū)中短期前兆場異常與強(qiáng)震關(guān)系研究[J].地震研究，26（S1）：95-100.

付虹，錢曉東，毛玉平，等.2015.2014年云南魯?shù)镸S6.5 地震異常及預(yù)測[J].地震研究，36（2）：141-147.

付虹，王世芹，秦嘉政，等.2007.2007年5—6月滇西南地區(qū)M≥5成組地震活動中、短臨異常特征及預(yù)測[J].地震研究，30（4）：303-310.

國家地震局預(yù)測預(yù)防司.1997.測震學(xué)分析預(yù)報方法[M].北京：地震出版社：105-108.

黃浩，付虹.2014.2008年以來滇西地區(qū)地震序列的譜振幅相關(guān)系數(shù)變化特征[J].地震學(xué)報，36（4）：631-639.

李艷娥，陳學(xué)忠，付虹.2012.2007年云南寧洱MS6.4地震前滇西南地區(qū)震源動力學(xué)參數(shù)時空變化特征[J].地震，32（1）：28-39.

李艷娥，陳學(xué)忠.2017.2011年日本MW9.1地震前后破裂區(qū)內(nèi)視應(yīng)力時空變化[J].地震，37（4）：10-21.

石紹先，付虹，沈斌，等.1991.1988 年瀾滄耿馬大震的預(yù)報和再認(rèn)識[J].西北地震學(xué)報，13（1）：35-40.

易桂喜，聞學(xué)澤，辛華，等.2013.龍門山斷裂帶南段應(yīng)力狀態(tài)與強(qiáng)震危險性研究[J].地球物理學(xué)報，56（4）：1112-1120.

Chen X，Shearer P M.2013.California foreshock sequences suggest aseismic triggering process[J].Geophysical Research Letters，40（11）：2602-2607.

Tracking Analysis of Seismic Anomalies of the 2021Yangbi，Yunnan MS6.4 Earthquake

FU Hong，LIU Zifeng，HU Xiaojing，LI Zhirong，HONG Min，PENG Guanling，WANG Guangming，NI Zhe

（Yunnan Earthquake Agency，Kunming 650224，Yunnan，China）

Abstract

In this study，we made a series of researches and retrospections of the foreshock sequence，as well as the anomalies of seismological，deformation and fluid observations at fixed observational sites before the occurrence of the Yangbi MS6.4 earthquake.The results are as follows：seismicity of the foreshock sequence is relatively complicated;some seismological anomalies can be extracted from the foreshock sequence.However，its difficult to make a short-term forecast based on the characteristics of the foreshock sequence.The anomalies of deformation at fixed observational sites arose early and substantially，which could shed light on the prediction of earthquake magnitude.The anomalies of fluid observation prior to the MS6.4 earthquake were transient and intensive，which could provide better evidence for forecasting of the timing of the mainshock.Features of the anomalies of deformation and underground fluid are much similar in time and space with the features of the faults later-stage meta-instability，which may be because of the faults meta-instability evolution.One could improve the accuracy of the earthquake prediction by combining the characteristics of foreshock sequence and regional seismic risk.We suggested that it was a promising method for comprehensive earthquake prediction by reasonably using the characteristics of the anomalies emerged in different seismogenic stages.

Keywords：the Yangbi MS6.4 earthquake;foreshock sequence;fixed survey of the crust deformation;underground-fluid anomalies;comprehensive prediction