馮麗麗 李霞 趙玉紅 張朋濤

1)中國地震局青海格爾木青藏高原內(nèi)部地球動力學(xué)野外科學(xué)觀測研究站，西寧 810001 2)青海省地震局，西寧 810001

0 引言

極化法(Polarization)由Hayakawa等(1996)提出，因其能夠反映地震地磁擾動異常的產(chǎn)生、傳播和分布特征，被廣泛應(yīng)用于國內(nèi)外地震電磁擾動異常提取(Hayakawa et al，2000；Molchanov et al，2003；Prattes et al，2008； 馮志生等，2010； 李琪等，2015； 廖曉峰等，2019； 艾薩·伊斯馬伊力等，2020； 樊文杰等，2021)。但以往研究主要以單臺資料分析為主，對震前電磁輻射異常的空間特征研究較少。2019年，基于2015—2018年中國地磁臺網(wǎng)磁通門磁力儀產(chǎn)出的秒采樣觀測數(shù)據(jù)，在對分析結(jié)果進(jìn)行了適當(dāng)?shù)奶幚砗?，總結(jié)了地磁垂直強(qiáng)度極化高值異常與后續(xù)地震的時空關(guān)系，建立了完整的地震預(yù)報指標(biāo)體系(中國地震局監(jiān)測預(yù)報司，2020)，使得該方法得到了較大發(fā)展。賀曼秋等(2019)在2017年九寨溝7.0級和精河6.6級地震前，在2個震中附近同時提取到了地磁垂直強(qiáng)度極化異常； 李霞等(2021)在2017年西藏米林6.9級地震前10天提取到了青藏高原出現(xiàn)的大范圍地磁垂直強(qiáng)度極化高值異常； 管貽亮等(2021)對2015年尼泊爾地震前的地磁垂直強(qiáng)度極化異常進(jìn)行了分析并從構(gòu)造動力演化角度對地磁極化異常的孕育和發(fā)展進(jìn)行了解釋； 馮麗麗等(2021)總結(jié)了2015—2020年中國大陸西部地磁垂直強(qiáng)度極化異常和地震的關(guān)系，認(rèn)為地磁垂直強(qiáng)度極化異常出現(xiàn)后半年內(nèi)，空間分布圖中0.2閾值線附近發(fā)生6級以上強(qiáng)震的可能性較大，且高值異常區(qū)面積與震級存在正相關(guān)關(guān)系。

地震前也常觀測到地下流體的異常變化。在2021年青海瑪多7.4級(鐘駿等，2021)、2017年新疆精河6.6級(鐘駿等，2018)、2010年青海玉樹7.1級地震(李滔等，2012)等中國大陸西部強(qiáng)震前，均在震中周邊多個臺站觀測到了地下流體化學(xué)組分的濃度變化、溫度變化或水位變化等。

1 地磁垂直強(qiáng)度極化法

Hayakawa等(1996)提出的地磁垂直強(qiáng)度極化法能夠突出源于震源區(qū)的異常信號，抑制來源于外源場的信號。具體數(shù)學(xué)表述如下

(1)

式中，Z(ω)為地磁垂直分量的譜幅度值，H(ω)為地磁水平分量的譜幅度值；HX(ω)為地磁南北分量譜幅度值，HY(ω)為地磁東西分量譜幅度值。地磁分量單位為nT，ω為圓頻率。

1.1 數(shù)據(jù)處理

以往震例表明，震前的極化異常在不同臺站出現(xiàn)時間可能略有差異，在研究異?？臻g特征時，可能因為前后出現(xiàn)時間的差異導(dǎo)致異常的空間范圍與實際不相符的情況，因此對逐日值進(jìn)行了5日滑動平均處理，得到Y(jié)ZH1。

1.2 歸一化

因觀測環(huán)境及儀器的差異，各臺YZH1波動范圍不盡相同。在提取YZH1高值進(jìn)行空間特征研究時，可能出現(xiàn)因臺站差異導(dǎo)致的失真現(xiàn)象。為了消除臺站影響，需要對各臺YZH1進(jìn)行歸一化。為了使結(jié)果更直觀，對歸一化后的結(jié)果進(jìn)行了置零處理。具體算法為

(2)

我們稱YZH2為歸一置零極化值，當(dāng)其高于0時，等同于YZH1高于2倍方差。在空間上對現(xiàn)有臺站的歸一置零極化值進(jìn)行克里金法數(shù)學(xué)插值，即可得到大陸西部地區(qū)每日的地磁垂直強(qiáng)度極化空間分布。

2 異常分析

2.1 地磁垂直強(qiáng)度極化異常

2021年10月27日前后，全國范圍內(nèi)77個參與計算臺站中共有17個臺站出現(xiàn)高值異常。從青海地區(qū)的YZH1時序曲線來看(圖1)，從2021年10月20日起青海地區(qū)多個臺站的YZH1出現(xiàn)持續(xù)升高，至10月27日前后達(dá)到極大值，后持續(xù)回落。這是2021年6月(瑪多7.4級地震)以來青海地區(qū)最顯著、幅度最大的一次地磁垂直強(qiáng)度極化高值異常。

圖1 2021年1月以來青海地區(qū)各臺地磁垂直強(qiáng)度極化時序曲線

圖2 2021年10月27日地磁垂直強(qiáng)度極化YZH2異?？臻g分布

經(jīng)過后續(xù)處理和分析，發(fā)現(xiàn)在異常幅度最大、異常臺站最多的時間，即2021年10月27日，在青海西北部形成1個面積大于6×104km2的高值異常區(qū)(圖2)，異常面積為6.6×104km2。大于閾值0.2的范圍內(nèi)有3個異常臺站，分別為德令哈臺、貴德臺和金銀灘臺，詳細(xì)信息見表1。其中歸一置零極化值最高的為金銀灘臺，達(dá)到0.9。

表1 2021年10月27日青海地區(qū)垂直強(qiáng)度極化信息

2.2 地下流體準(zhǔn)同步短期異常

2021年7月至11月，青海地區(qū)陸續(xù)出現(xiàn)多項流體短臨異常。其中佐署臺動水位的下降異常及淺層水溫的上升異常出現(xiàn)時間最早，2個測項出現(xiàn)異常的時間一致，均為7月初。而平安靜水位破年變異常隨后在7月底至8月初出現(xiàn)。門源靜水位突降異常、樂都逸出氣氡突升異常以及西寧逸出氣氡濃度突降異常均出現(xiàn)在2021年9月。上述異常大多僅持續(xù)較短的時間，而西寧逸出氣氡濃度突降異常一直持續(xù)到11月初，幾乎與此同時，共和的2套中層水溫也同時出現(xiàn)了突降異常，并于12月初再次出現(xiàn)同步變化。德令哈水溫在12月上、中旬出現(xiàn)突升—轉(zhuǎn)平—突降的“倒U”形變化，整個變化過程持續(xù)約半個月。至此，青海地區(qū)地下流體異常數(shù)量達(dá)到8項(共和2套水溫觀測儀器相同且裝置于同一井孔中，因此計為1項)，空間分布見圖3，時序曲線見圖4。

圖3 門源地震前青海地區(qū)地下流體異常臺站分布

圖4 2021年6月以來青海多個臺站地下流體準(zhǔn)同步異常

從震中距及異常出現(xiàn)時間來看，最遠(yuǎn)的德令哈臺和相對稍遠(yuǎn)的共和臺出現(xiàn)異常的時間最晚。而最近的門源臺于震前4個月出現(xiàn)異常，在此次準(zhǔn)同步變化過程中處于中段。西寧臺、佐署臺、平安臺、樂都臺及共和臺震中距相差不大，但異常出現(xiàn)的時間最大相差了5個月，無明顯的時空演化特征。

從異常類別來看，共3個臺出現(xiàn)水位變化，3個臺的4套儀器觀測顯示水溫變化，2個臺出現(xiàn)逸出氣氡濃度變化。其中平安臺、佐署臺和門源臺出現(xiàn)了水位變化，平安臺打破了原本的年變化上升趨勢，水位轉(zhuǎn)平下降； 門源臺出現(xiàn)了水位上升變化(平安臺與門源臺為靜水位觀測，數(shù)值下降表明水位上升)； 而佐署臺則呈現(xiàn)了持續(xù)性的下降變化。水溫變化在德令哈臺和佐署臺主要表現(xiàn)為升高變化，而共和臺2套儀器的觀測結(jié)果均顯示階降變化。而從逸出氣氡濃度來看，樂都臺出現(xiàn)稍早，主要為上升變化，在樂都臺開始下降階段西寧臺出現(xiàn)了上升變化。

青海地區(qū)地下流體測項在2021年5月22日瑪多7.4級地震前1～3個月，曾提出多項準(zhǔn)同步異常變化，包括玉樹水溫(蘇維剛等，2021)、佐署動水位、共和逸出氣氡及共和中層水溫等(劉磊等，2021)。2016年門源6.4級地震前1～3個月，也曾出現(xiàn)過3項準(zhǔn)同步變化，包括玉樹水溫“V”型變化、樂都?xì)怆笨焖偕弋惓?劉磊等，2020)、平安靜水位破年變變化等(楊曉霞等，2016)，部分異常臺站和測項與此次異常特征相似。且對應(yīng)地震的時間均在3個月內(nèi)。

3 “時、空、強(qiáng)”綜合分析

3.1 時間判定

表2 列出了青海地區(qū)地磁垂直強(qiáng)度極化異常信息。自2015年以來，青海地區(qū)共出現(xiàn)8次地磁垂直強(qiáng)度極化異常，其中6次在異常出現(xiàn)后8個月內(nèi)異常區(qū)內(nèi)及附近發(fā)生了5.6級以上地震，對應(yīng)率為75%，具有較高的信度。而其中4次地震發(fā)生在異常出現(xiàn)后的60天內(nèi)，表明該異常具有一定的短期指示意義。從異常出現(xiàn)時間推測優(yōu)勢發(fā)震時間在2021年12月27日前。但考慮到目前震例較少，且部分震例對應(yīng)時間較長，因此在發(fā)震時間判定中同時考慮了地下流體準(zhǔn)同步變化的時間。2021年7月以來，青海地區(qū)陸續(xù)出現(xiàn)多項流體異常。2021年11月再次出現(xiàn)多項準(zhǔn)同步變化，這種現(xiàn)象對于后續(xù)3個月內(nèi)青海地區(qū)的強(qiáng)震有較好的指示意義。綜合上述分析，判斷發(fā)震時間為2021年12月至2022年1月。

表2 青海及周邊地區(qū)極化異常信息

圖5 青海地區(qū)地磁垂直強(qiáng)度極化YZH2異?？臻g分布

3.2 地點判定

圖5給出了青海及周邊地區(qū)垂直強(qiáng)度極化異常空間分布與后續(xù)地震的震中位置，震中全部位于異常區(qū)內(nèi)或閾值線(歸一置零極化值為0.2的等值線)附近。因此，預(yù)測發(fā)震地點在異常區(qū)內(nèi)或閾值線附近。

3.3 震級判定

由表2 可知，青海地區(qū)震例震級在5.6～7.4級之間。根據(jù)以往對中國大陸西部的震例研究(馮麗麗等，2021)，后續(xù)地震震級與高值異常區(qū)面積有近似正相關(guān)關(guān)系，并存在擬合關(guān)系。此次異常區(qū)面積與夏河5.7級地震前異常區(qū)面積相當(dāng)。綜合分析認(rèn)為，目標(biāo)地震可能在6級左右，因此最終給出預(yù)測震級為5.6～6.4級。

3.4 綜合預(yù)測結(jié)果

基于青海及周邊地區(qū)地磁垂直強(qiáng)度極化異常，并綜合2021年7—11月以來青海地區(qū)地下流體準(zhǔn)同步異常變化，基于以往震例分析認(rèn)為2021年11月23日至2022年1月23日，青海西北部地磁垂直強(qiáng)度極化高值區(qū)內(nèi)可能發(fā)生5.6～6.4級左右地震，并填報地震短臨預(yù)測卡片。此次門源6.9級地震發(fā)生于異常后的73天，震中位于預(yù)測區(qū)邊界，發(fā)震時間和震中位置均符合預(yù)期，但震級比預(yù)期高0.5級。

4 結(jié)論與討論

(1)此次門源6.9級地震前73天在青海西北部出現(xiàn)了地磁垂直強(qiáng)度極化異常，高值區(qū)面積為6.6×104km2，空間上位于此次地震震中西南。震前，通過對該異常在青海地區(qū)的預(yù)測指標(biāo)進(jìn)一步分析認(rèn)為異常出現(xiàn)后60天至稍長時間，高值區(qū)附近可能發(fā)生6級左右地震。

(2)此次門源6.9級地震前1～6個月，青海地區(qū)出現(xiàn)了多項地下流體異常，在時間上表現(xiàn)為準(zhǔn)同步變化。根據(jù)異常指標(biāo)，認(rèn)為發(fā)震時間可能在2021年12月至2022年1月。

(3)此次門源6.9級地震的發(fā)生表明震前的分析和判斷基本正確。地磁垂直強(qiáng)度極化異?？臻g閾值線是地震地點預(yù)測的良好指標(biāo)。該異常在青海地區(qū)可作為短臨異常，指示未來2～3個月內(nèi)的中強(qiáng)以上地震，地下流體的準(zhǔn)同步異常具有較強(qiáng)的時間指示意義。

優(yōu)質(zhì)的時間、空間預(yù)測指標(biāo)是本次地震短臨跟蹤分析成功的關(guān)鍵。但同時應(yīng)指出，此次預(yù)測的震級仍存在0.5級的偏差，主要原因是目前在青海地區(qū)尚未梳理出高質(zhì)量的震級預(yù)測指標(biāo)，需要在下一步工作中繼續(xù)加強(qiáng)此方面工作。