鐘駿 王博 周志華 晏銳

中國地震臺網中心，北京 100045

0 引言

據中國地震臺網測定，北京時間2021年5月22日2時4分青海果洛藏族自治州瑪多縣(34.59°N，98.34°E)發(fā)生MS7.4地震，震源深度17km，震中最高烈度Ⅹ度。此次地震造成數人受傷，大量房屋倒塌，部分道路、橋梁等基礎設施被破壞或受損(李智敏等，2021)。地震地下流體是地震前兆觀測的重要手段，此次震前是否觀測到表征地震即將發(fā)生的地下流體異常信號備受關注。

地下流體是地震孕育與發(fā)生的主要影響因素，在地震預測研究和震情跟蹤工作中發(fā)揮著重要作用(王鐵城等，1994；Roeloffs，1998；車用太等，2006；孫小龍等，2020)。據不完全統(tǒng)計，中國歷史地震資料中與地下流體有關的地震前兆異常占全部異常的50%以上(汪成民，1990)。Montgomery等(2003)通過研究不同空間尺度地下水對地震的響應特征發(fā)現，地下水變化的響應范圍可達數千千米。李軍等(2005)對理塘毛埡溫泉三十余年的水溫觀測資料研究發(fā)現，該溫泉水溫異常與川滇地區(qū)的6級以上地震存在較好的對應關系。Ghosh等(2009)指出，土壤氣濃度異常變化與地震震級呈正相關，即土壤氣異常幅度隨著地震震級的上升而增大。劉耀煒等(2015)總結魯甸MS6.5地震震前地下流體典型異常特征，認為區(qū)域應力加載作用可能與斷層裂隙的開啟和閉合有關，從而引起水位、水溫、深部氣體以及水化學離子組分發(fā)生顯著變化。此外，眾多學者基于《中國震例》研究發(fā)現地下流體前兆異常在時空上具有階段性和不均勻性，演化過程主要表現為“向震中收縮”“構造控制”“相對集中”等3種典型特征，特別是進入短臨階段后，異常向震源區(qū)收縮的趨勢更為明顯；而震前異常在數量上主要表現為“持續(xù)增長”型和“先增后減”型2類，且以“先增后減”型居多，趨勢性異常轉折結束、新突變異常增多往往是地震孕育進入短臨階段的標志(張肇誠，1988；晏銳等，2004；付虹等，2008；陳棋福，2008；蔣海昆，2018、2019；蔣海昆等，2009；孫小龍等，2016)。

本文通過系統(tǒng)梳理瑪多MS7.4地震前地下流體前兆異常，研究這些異常的震前動態(tài)演化特征并進行強震異常特征對比，進而積累震例資料，為提升地下流體前兆異常認識水平及未來震情趨勢判定提供借鑒和參考。

1 觀測資料概述

青藏高原北部地下流體觀測點分布不均，主要集中在青藏高原東北緣，其他地區(qū)監(jiān)測能力較弱?，敹郙S7.4地震震中200km范圍內僅有2個地下流體觀測點，500km范圍內存在21個地下流體觀測點，共38個測項(圖1)，其中水溫18項，水位8項，氣氡7項，水氡5項，以水溫觀測為主。按照測項預報效能分類，可分為A類1項，B類和C類各12項，D類4項，未評估9項，以B類和C類為主。

圖1 瑪多MS7.4地震震中距500km范圍內前兆異?？臻g分布震源機制來自全球質心矩張量數據庫(GCMT)；藍色圓圈為MS3.0余震；時間截至2021年8月19日

瑪多MS7.4地震前共出現5項地下流體異常，分別為玉樹水溫、共和水溫(2套探頭)、共和氣氡、佐署水位、門源水位，基本特征見表1。

表1 瑪多MS7.4地震震中距500km范圍內地下流體異?；拘畔?/p>

1.1 玉樹水溫

玉樹臺位于巴顏喀拉山南麓的甘孜-玉樹斷裂附近，水溫觀測井深105m，觀測儀器為SZW-Ⅱ型數字水溫儀。含水層巖性為中生代侏羅紀花崗巖(楊曉霞等，2016)。自2007年8月開始觀測以來，整體呈多年穩(wěn)定的趨勢上升變化，歷史資料分析顯示，玉樹水溫異常一般出現在震前1～3個月，異常持續(xù)時間在12～41天，基本上遵循“下降—上升”的“V”字形變化過程，重復性較好(王博等，2016)。該異常曾在2008年3月21日新疆于田MS7.3、2008年5月12日四川汶川MS8.0、2010年4月14日青海玉樹MS7.1、2015年4月25日尼泊爾MS8.1、2016年1月21日青海門源MS6.4和2016年10月17日青海雜多MS6.2等地震前出現。

2021年3月15日玉樹水溫在正常背景下出現快速下降，幅度約0.0027℃，16日開始轉折上升，過程中發(fā)生了3月19日西藏比如MS6.1地震，3月21日受異常核實影響再次出現下降，幅度約0.0023℃，3月22日再次轉折上升，至4月23日恢復至正常背景，29天后發(fā)生瑪多MS7.4地震，異常過程整體表現出不規(guī)則的“V”字形變化特征(圖2(a))。

圖2 瑪多MS7.4地震地下流體異常時序圖

1.2 共和氣氡、水溫

共和臺地處昆侖-秦嶺緯向構造帶與河西系構造復合部位的沉降帶，位于共和盆地北緣斷裂與共和盆地南緣斷裂之間，觀測井深174.35m，氣氡和水溫(2套探頭)觀測分別使用BG2015R型數字氣氡儀、SZW-I和SZW-Ⅱ型數字水溫儀。含水層巖性為砂質黏土、礫石(邱鵬成等，2010)。自2011年開始對共和氣氡進行觀測，觀測數據較穩(wěn)定，呈現冬高夏低的年變特征，2019年5月因儀器老化對其進行更換，更換后觀測時間較短，無正常年變特征；分別于2011年和2015年開始對共和水溫2套儀器進行觀測，觀測數據同步性較好，無年變特征。共和氣氡和水溫(2套探頭)觀測數據于2021年2月28日出現同步上升變化，幅度分別為17.4Bq/L、0.047℃和0.07℃。氣氡濃度于3月1—10日持續(xù)呈小幅波動上升變化，累計上升幅度為29Bq/L，之后開始緩慢下降，73天后發(fā)生瑪多MS7.4地震，震后氣氡濃度維持下降狀態(tài)(圖2(b))；水溫于3月1—2日同步下降恢復至變化前的趨勢水平，幅度分別為0.072℃和0.21℃，81天后發(fā)生瑪多MS7.4地震(圖2(c)～(d))。

1.3 佐署水位

佐署臺位于拉脊山構造帶的北緣、日月山斷裂和拉脊山斷裂交匯區(qū)域，觀測井深107m，觀測儀器為SWY-Ⅱ型數字水位儀。含水層巖性為砂礫巖(張昱等，2008)。自1986年5月開始對佐署水位進行觀測，期間受數字化改造、臺站改造、儀器故障以及修路施工等影響，水位觀測極不穩(wěn)定，數據質量較差，于2019年10月施工結束后才開始正常觀測。自2019年10月觀測以來水位呈波動趨勢下降變化，直至2021年2月25日出現突降變化，幅度為0.268m，26日開始呈緩慢上升趨勢，85天后發(fā)生瑪多MS7.4地震，并出現幅度為0.097m的同震階升變化(圖2(e))。

1.4 門源水位

門源水位地處北祁連新元古代-早古生代縫合帶，達坂山斷裂北約15m，觀測井深105.4m，觀測儀器為SWY-Ⅱ型數字水位儀。含水層巖性為砂巖、間夾泥巖(邱鵬成等，2010；汪發(fā)耀等，2020)。自2017年11月開始對門源水位進行觀測，觀測資料連續(xù)穩(wěn)定，無正常年變特征。2019年9月24日—12月4日門源水位出現大幅上升，幅度約6.85m，12月5日開始轉折快速下降，至2020年5月1日恢復至以往平穩(wěn)變化趨勢，累計下降4.14m，異常結束386天后發(fā)生2021年5月22日瑪多MS7.4地震(圖2(f))。

2 異常特征分析

2.1 數量特征

瑪多MS7.4地震震中500km范圍內共5項異常，約占測項總數的13%。其中，水溫異常2項，水位異常2項，氣氡異常1項，分別占各觀測手段的11%、25%和14%，異常測項以水位和水溫為主。震中100km范圍內無流體觀測；100～200km范圍內流體觀測2項，無異常；200～300km范圍內流體觀測5項，異常3項，異常測項比為60%；300～400km范圍內流體觀測16項，異常1項，異常測項比約為6%；400～500km范圍內流體觀測15項，異常1項，異常測項比約為7%。綜上，異常以震中距200～300km為主(圖3)。

圖3 不同震中距范圍內的地下流體異常測項數、背景測項數及異常測項比

2.2 時間特征

從時間演化上看(圖4)，震前地下流體異常主要為3個月內的短期異常，有4項異常在2021年2月之后出現，占異?？倲档?0%。異常起始和結束時間的階段性特征明顯，僅門源水位1項趨勢異常出現在震前19個月，持續(xù)213天后，在震前12個月結束；其余4項短期異常均出現在震前3個月內并在震前結束，異常持續(xù)時間較短(2～40天)。整個演化過程呈現出異常數量由少增多、在震前3個月內顯著加速的特點。

圖4 瑪多MS7.4地震前地下流體異常時間進程

2.3 時空特征

從時空演化上看，震中距500km范圍內的38項地下流體觀測主要集中在震中東北方向，震中距200km內僅有的2個地下流體觀測點，未出現異常，5項異常均位于震中200km以外，主要集中在200～300km范圍內(3項，占異?？倲档?0%)，表現出流體觀測和異常測項存在空間分布上的不均勻性。異常最早出現在距震中較遠的位置(門源水位，震中距421km)，隨著時間的推移，在震前3個月至震前2個月，震中400km內的異常開始集中出現并向震中位置遷移，其中距震中最近的異常為玉樹水溫，震中距為213km(圖5)。

圖5 瑪多MS7.4地震前地下流體異常時空演化特征

2.4 形態(tài)特征

前人通過震例研究發(fā)現，部分流體異常形態(tài)存在明顯的重現性(王博等，2016；Sunetal，2018；蘇維剛等，2021)。玉樹水溫異常在周邊MS6.0以上地震以及巴顏喀拉塊體邊界附近發(fā)生的多次MS7.0以上地震前重復出現，異常形態(tài)主要表現為幅度不同的“下降—上升”的“V”字形變化(圖6)，該測項在瑪多MS7.4地震前出現同樣的異常形態(tài)，對于地震的時、空、強判定有一定的指示意義(表2)。而共和氣氡、共和水溫、佐署水位和門源水位異常屬于觀測以來首次出現的前兆異常變化，相關異常特征及震例信息有待積累和進一步驗證。

圖6 玉樹水溫2007—2021年觀測曲線

表2 玉樹水溫異常變化與對應地震的關系

3 討論

截至2021年8月，瑪多MS7.4地震是繼汶川MS8.0地震之后發(fā)生在中國大陸的最大地震，震中位于青藏高原北部巴顏喀拉塊體內部的昆侖山口-江錯斷裂(王未來等，2021)。該斷裂整體走向N230°W～270°，傾向NE，傾角80°，運動方式以左旋走滑為主，兼有少量正斷分量(華俊等，2021)。歷史上，巴顏喀拉塊體周邊活動斷裂發(fā)育，構造運動劇烈，強震頻繁發(fā)生，尤其是巴顏喀拉塊體東邊界附近曾多次發(fā)生MS7.0以上地震，如1973年四川爐霍MS7.6、1976年四川松潘MS7.2、1990年青海共和MS7.0、2008年四川汶川MS8.0、2010年青海玉樹MS7.1、2013年四川蘆山MS7.0和2017年四川九寨溝MS7.0地震等(朱亞戈等，2021)。根據震例總結統(tǒng)計(張肇誠，1988；蔣海昆，2018、2019；晏銳等，2018；中國地震臺網中心，2019)，這些地震震前500km內均存在一定的地下流體異常變化(表3)。

表3 巴顏喀拉塊體東邊界MS7.0以上地震震中距500km范圍內地下流體異常情況

其中，1973年爐霍MS7.6地震前流體異常僅有1項，為姑咱水氡，震中距為220km；1976年松潘MS7.2地震前流體異常共4項，分別為江油水位、理塘水溫、松潘水氡和姑咱水氡，震中距范圍50～460km；1990年共和MS7.0地震前流體異常共13項，分別為佐署水位、湟源流量、湟源水氡、西寧水氡、長寧水氡、樂都水氡、民和水氡、張掖水氡、山丹水氡、五泉山水氡、格爾木水氡、松潘拱北水氡和山丹氣汞，震中距范圍110～495km；2008年汶川MS8.0地震前地下流體異常共24項，分別為水富水位、蒲江水位、德陽水位、小金水位、瀘州水位、邛崍水位、二道橋水位、北川水位、巴塘水溫、理塘水溫、鄉(xiāng)城水溫、龍頭溝水溫、道孚水溫、二道橋水溫、大足水溫、武都殿溝水氡、武山1號泉水氡、武山22號井水氡、清水溫泉水氡、天水花牛水氡、昭覺CO2、普格CO2、布拖CO2、龍頭溝氣體總量，震中距范圍71～492km；2010年玉樹MS7.1地震前流體異常共2項，為玉樹水溫和德令哈水溫，震中距分別為22～480km；2013年蘆山MS7.0地震前流體異常共10項，分別為邛崍水位、瀘州水位、理塘水溫、二道橋水溫、龍頭溝水溫、瀘定水溫、理縣水氡、普格CO2、布拖CO2和鹽源CO2，震中距范圍27～350km；2017年九寨溝MS7.0地震前流體異常共14項，分別為干鹽池水位、大足水位、北碚水位、清水李溝流量、二道橋水溫、干鹽池水溫、瀘定水溫、巴塘水溫、理塘水溫、武都殿溝水氡、武山1號泉水氡、武山22號井水氡、勉縣水氡和洋縣水氡，震中距范圍108～493km。

用濃度為0.1%的MKnO4對厚萼凌霄種子消毒10 min，然后用流水沖洗。將海水配制成1%、5%、10%、20%、30%、40%的6組濃度梯度，供實驗備用。在每個干燥潔凈的培養(yǎng)皿中放入濾紙2張和經滅菌處理的厚萼凌霄種子100粒。其中，實驗組培養(yǎng)皿分別加入6個不同濃度梯度的海水，對照組培養(yǎng)皿加入蒸餾水。將實驗組、對照組培養(yǎng)皿置于光照培養(yǎng)箱中培養(yǎng)種子，培養(yǎng)條件均為25℃、12 h光照（1 000 lx）和12 h黑暗，并在實驗組培養(yǎng)皿中定期添加適量相應濃度的海水，以保證種子被充分浸潤。

數量特征上(圖7)，震前異常數量小于5個的地震有3次，分別為爐霍地震、松潘地震和玉樹地震，其中，松潘地震與瑪多地震異常數量相近，觀測手段相同，相似性較好；時間特征上(表3，圖8)，除瑪多地震外，震前短期異常占比超過80%的僅有爐霍地震和玉樹地震(均為100%)，但這兩次地震前流體異常數較少，時間特征統(tǒng)計意義較弱。其余地震短期異常占比均小于30%，震前3個月內短期異常加速特征不顯著。時空特征上(圖8)，共和地震和汶川地震前6個月，異常存在由外圍向震中收縮的遷移特征。進入震前6個月內，遷移特征消失，異常在距震中不同范圍內均有出現，與瑪多地震由外圍向震中收縮的遷移特征差異較大。其余地震異常時空演化特征不明顯。形態(tài)特征上，如前面所述，玉樹水溫在汶川地震和玉樹地震均出現“V”字形變化，與瑪多地震前形態(tài)十分相似(圖6)。因此，從整體異常特征上看，瑪多地震MS7.4震前地下流體異常和松潘MS7.2、玉樹MS7.1地震相似性較高。

圖7 MS7.0以上地震震中距500km范圍內地下流體異常數量對比

圖8 MS7.0以上地震500km范圍內地下流體異常時空對比

研究認為，流體前兆響應和地震孕育過程中的區(qū)域應力加載作用密切關系(車用太等，1999)。一般情況下，7級左右地震在地下流體方面的異常持續(xù)時間一般在1年左右，甚至更短(蔣海昆等，2009)。松潘地震地下流體異常的持續(xù)時間基本在1年左右，與前人經驗一致?，敹嗟卣鹎?年尺度的前兆異常僅有門源水位1項，玉樹地震前無1年尺度的異常信息出現，可能反映了此時段處于應力加載初期緩慢的彈性積累階段，異?，F象并不突出。隨著應力的不斷積累，進入長期應力加載后的巖石膨脹階段，短期前兆異常開始集中顯現直至地震發(fā)生(張學民等，2018)。而瑪多地震和玉樹地震震中附近臺站稀少，可能是臨震破裂階段震源區(qū)附近未觀測到異?，F象的主要原因。整體上看，瑪多地震和玉樹地震震中空間位置較為接近、構造環(huán)境類似，震前3個月內地下流體短期異常在時空分布上也具有較好的相似性，可能具有類似的應力加載過程。而松潘地震相對較遠且震前無短期異常出現，說明其孕育過程可能存在一定的差別，需要進一步深入研究。

瑪多MS7.4地震前，中國地震臺網中心預報部地下流體研究室依據玉樹水溫異常的高重現性、高映震率以及相對固定的震例發(fā)生地點等優(yōu)勢給出了明確預報意見，證明了地下流體異常在短期地震預測預報中發(fā)揮著不可替代的作用。雖然此次玉樹水溫異常幅度較以往小一個數量級，導致預報震級偏低，但共和氣氡、共和水溫以及佐署水位等3項短期異常的集中出現以及瑪多地震的發(fā)生，很可能說明震前出現的這些地下流體短期異常與地震孕育過程存在某些有待探索的物理聯系。結合歷史震例對比結果，認為即使地下流體異常個數較少或異常測項比較低，巴顏喀拉塊體發(fā)生MS7.0以上地震的可能性仍然較高，特別是出現類似玉樹水溫這種關鍵性短期異常，以及出現群體異常集中分布或由外圍向震中收縮等，這些特征對于未來地震趨勢研判均具有重要的參考意義。

在過去的半個多世紀里，我國地震地下流體觀測網從無到有、由點及面，觀測臺站數量逐年增多，已發(fā)展成為世界上規(guī)模最大、地震監(jiān)測能力最好的現代化地震地下流體觀測網(車用太等，2017)。然而，研究表明，絕大多數地震前的地下流體異常數量占比較低，基本不超過20%(張立，2013；孫小龍等，2016；楊曉霞等，2016；鐘駿等，2018；王俊等，2020)，根據震前捕捉到的流體異常提出預測意見的地震次數更是不到10%(王廣才等，2003)。這一事實固然與地震的孕震環(huán)境不同以及觀測井孔所在構造條件不同有關，但更突出反映了現有地下流體觀測網空間分布明顯存在不均勻的弊端。同時，本次地震震中500km范圍內測項的預報效能主要為B類和C類，其次為D類，而A類僅有1項，未評估甚至高達9項，說明該區(qū)域內地下流體觀測資料的預報能力整體偏弱，導致監(jiān)測預報水平十分有限，遠不能滿足防震減災的需求。因此，有必要進一步優(yōu)化與升級，尤其是在臺網布局上，亟需增強地震弱監(jiān)視區(qū)的監(jiān)測水平，并從基礎資料、觀測質量、影響因素和震例信息等方面對每個測項的預報效能進行常態(tài)化動態(tài)評估，從而不斷提升我國地下流體臺網的監(jiān)測預報水平以及地震分析預報能力。

4 結論

通過對瑪多MS7.4地震前地下流體觀測資料及異常特征分析，結合比巴顏喀拉塊體多次MS7.0以上地震前的地下流體異常特征，得出以下結論：

(1)瑪多MS7.4地震前地下流體異常的時間階段性特點明顯，整個演化過程呈現出由外向震中收縮的遷移特征，短期階段異常數量呈現加速增多的特征。

(2)玉樹水溫異常具有明顯的重現性。異常形態(tài)多表現為“下降—上升”的“V”字形變化，對于地震的時、空、強判定有一定的指示意義。

(3)瑪多MS7.4地震震前地下流體異常和松潘MS7.2地震、玉樹MS7.1地震具有一定的相似性；在區(qū)域應力加載作用方面，瑪多MS7.4地震與玉樹MS7.1地震可能具有類似的應力加載過程。