裴東洋 戴勇 王磊 格根

（中國(guó)呼和浩特 010010 內(nèi)蒙古自治區(qū)地震局）

0 引言

地震與地磁場(chǎng)的關(guān)系是地震電磁學(xué)的主要研究方向之一（馮志生等，2019）。地震地磁場(chǎng)異常主要研究方法有地磁日變化、磁擾動(dòng)、基于磁測(cè)深的異常分析方法等等，這些方法已在國(guó)內(nèi)外得到廣泛應(yīng)用（丁鑒海等，2009；艾薩· 伊斯馬伊力等，2020）。其中，地磁低點(diǎn)位移法具有較顯著的地震短期預(yù)測(cè)意義?！暗卮诺忘c(diǎn)位移”現(xiàn)象是指一個(gè)大區(qū)域的地磁垂直分量日變化極小值出現(xiàn)時(shí)間（簡(jiǎn)稱低點(diǎn)時(shí)間）明顯地與另一個(gè)大區(qū)域的低點(diǎn)時(shí)間不同，而每個(gè)大區(qū)域內(nèi)部低點(diǎn)時(shí)間又基本一致，2 個(gè)區(qū)域之間有明顯的突變分界線（簡(jiǎn)稱低點(diǎn)位移線）,并且2 個(gè)區(qū)域之間的低點(diǎn)時(shí)間相差2 h 以上（馮志生等，2009）。國(guó)內(nèi)在地磁低點(diǎn)位移現(xiàn)象形成機(jī)理及與地震間的關(guān)系的研究中已取得一些進(jìn)展（李軍輝等，2013）。趙潔等（2005）認(rèn)為，地下深部高溫物質(zhì)的對(duì)流施加于地殼時(shí)，會(huì)引起地殼構(gòu)造變動(dòng)而形成地磁低點(diǎn)位移異常。張學(xué)民等（2008）研究了地磁低點(diǎn)位移的頻譜特征及機(jī)理，認(rèn)為低點(diǎn)位移發(fā)生時(shí)，8.5 h、13.7 h 周期的頻譜值增加，地磁低點(diǎn)位移的發(fā)生與空間電子密度的增加間有較好的對(duì)應(yīng)性。馮志生等（2009）依據(jù)分布于低點(diǎn)位移線兩側(cè)的臺(tái)站觀測(cè)到的地磁垂直分量日變化反相位畸變變化，推測(cè)在低點(diǎn)位移線下方有外空變化磁場(chǎng)感應(yīng)電流聚集。丁鑒海等（1988，2008）研究了1966—2008 年低點(diǎn)位移異常與地震間的相關(guān)性，發(fā)現(xiàn)低點(diǎn)位移法可預(yù)測(cè)強(qiáng)震，特別是對(duì)6 級(jí)以上地震效果較好。

據(jù)中國(guó)地震臺(tái)網(wǎng)正式測(cè)定，2019 年9 月16 日20 時(shí)48 分甘肅張掖市甘州區(qū)（38.6°N，100.35°E）發(fā)生5.0 級(jí)地震，震源深度11 km。此次地震前29 天，即8 月18 日在甘肅中部地區(qū)有明顯的地磁低點(diǎn)位移現(xiàn)象。10 月28 日1 時(shí)56 分，甘肅甘南州夏河縣（35.1°N，102.69°E）發(fā)生5.7 級(jí)地震，震源深度12 km。此次地震前42 天，即9 月16 日在甘肅中東部有明顯的地磁低點(diǎn)位移現(xiàn)象。2 次地震震中位于地磁低點(diǎn)位移線明顯轉(zhuǎn)折處附近。本文擬在地磁低點(diǎn)位移法的基礎(chǔ)上，研究低點(diǎn)位移異常日地磁垂直分量低點(diǎn)時(shí)間的梯度分布特征，并分析地磁低點(diǎn)位移異常與中強(qiáng)地震間的關(guān)系。

1 原理與資料

1.1 原理

丁鑒海等（1988，2003，2004，2008）通過(guò)統(tǒng)計(jì)單臺(tái)地磁垂直分量低點(diǎn)時(shí)間的平均分布和多臺(tái)低點(diǎn)時(shí)間的空間變化規(guī)律給出了低點(diǎn)時(shí)間的分布規(guī)律及低點(diǎn)位移的判定標(biāo)準(zhǔn)、預(yù)測(cè)規(guī)則，即單臺(tái)低點(diǎn)時(shí)間集中分布在固定時(shí)間段內(nèi)，一般為地方時(shí)12 h 左右?？臻g上低點(diǎn)時(shí)間主要分布特征為隨經(jīng)度而變化，經(jīng)度相差15°時(shí)，低點(diǎn)時(shí)間相差1 h。如果偏離此時(shí)間過(guò)多，則認(rèn)為單臺(tái)低點(diǎn)發(fā)生位移。

隨著地震地磁臺(tái)網(wǎng)密度增加、震例積累、低點(diǎn)位移形成機(jī)理研究的深入，一些研究者對(duì)現(xiàn)有異常判定方法和預(yù)測(cè)規(guī)則作了進(jìn)一步完善。戴勇等（2019）提出低點(diǎn)時(shí)間變化的梯度值可以大幅縮小預(yù)測(cè)區(qū)域，梯度值1.5 h/（°）可以作為判定異常的閾值。本文將基于上述低點(diǎn)位移梯度方法，分析2019 年甘肅張掖5.0 級(jí)、夏河5.7 級(jí)地震前低點(diǎn)時(shí)間變化的梯度分布特征，研究低點(diǎn)時(shí)間梯度帶與中強(qiáng)地震震中分布間的關(guān)系，以期給出新的地磁低點(diǎn)位移判定指標(biāo)。

1.2 資料選取及處理過(guò)程

選擇中國(guó)大陸（15°—55°N，70°—140°E）地磁臺(tái)站垂直分量低點(diǎn)時(shí)間數(shù)據(jù)作為研究對(duì)象，數(shù)據(jù)來(lái)源于中國(guó)地震臺(tái)網(wǎng)中心國(guó)家地震科學(xué)數(shù)據(jù)中心（http://data.earthquake.cn）。在數(shù)據(jù)處理過(guò)程中對(duì)于同一臺(tái)站有多套儀器的情況，若儀器記錄數(shù)據(jù)真實(shí)有效且無(wú)明顯干擾，取各儀器記錄的低點(diǎn)時(shí)間的平均值為該臺(tái)站的低點(diǎn)時(shí)間。對(duì)于存在輕微干擾的臺(tái)站，采用低通濾波濾除干擾，若濾波后已無(wú)明顯干擾則納入研究范圍；否則，不采用此臺(tái)站數(shù)據(jù)。若因儀器故障導(dǎo)致當(dāng)日缺數(shù)情況嚴(yán)重或干擾嚴(yán)重，則不采用此臺(tái)站數(shù)據(jù)。

全國(guó)地震地磁臺(tái)站分布極不均勻，其中，華北、東南、南北地震帶地區(qū)臺(tái)間距多為幾十千米，而新疆維吾爾自治區(qū)、西藏自治區(qū)、青海省、內(nèi)蒙古自治區(qū)臺(tái)間距多為幾百千米（戴勇等，2019），故無(wú)法直接采用梯度法進(jìn)行計(jì)算。綜合考慮不同插值方法的效果，在計(jì)算過(guò)程中采用克里格方法（吳長(zhǎng)祥等，2008；姜金征等，2015）對(duì)中國(guó)大陸低點(diǎn)時(shí)間數(shù)據(jù)進(jìn)行空間插值處理。對(duì)于位于地震臺(tái)網(wǎng)內(nèi)部的區(qū)域，采用內(nèi)插算法；對(duì)于位于臺(tái)網(wǎng)外部的區(qū)域，則采用外推算法。在插值處理過(guò)程中，對(duì)網(wǎng)格數(shù)據(jù)文件間隔調(diào)試結(jié)果顯示，間隔大小對(duì)整體特征影響不大（戴勇等，2019；格根等，2020），最終確定經(jīng)向點(diǎn)數(shù)為1 401 個(gè)，緯向點(diǎn)數(shù)為801 個(gè)，經(jīng)向、緯向點(diǎn)間距均為0.05°的密度較合適。

低點(diǎn)時(shí)間隨經(jīng)度變化的特征稱為經(jīng)度效應(yīng)（丁鑒海等，1988；馮志生等，2009）。梯度計(jì)算過(guò)程中，首先將120°E 經(jīng)線作為基準(zhǔn)線以剔除經(jīng)度效應(yīng)；再對(duì)研究區(qū)內(nèi)每個(gè)網(wǎng)格點(diǎn)計(jì)算經(jīng)度向、緯度向梯度值的矢量和，將其作為該網(wǎng)格點(diǎn)的梯度值；最后繪制中國(guó)大陸地磁低點(diǎn)時(shí)間梯度空間分布圖。

2 2019 年甘肅2 次5 級(jí)以上中強(qiáng)地震震例

2.1 張掖5.0 級(jí)地震前地磁低點(diǎn)位移特征

2019 年9 月16 日張掖5.0 級(jí)地震前29 天，即8 月18 日中國(guó)大陸出現(xiàn)明顯的地磁低點(diǎn)位移現(xiàn)象。低點(diǎn)位移線經(jīng)過(guò)內(nèi)蒙古西部、甘肅中部、寧夏南部、山西南部、河北南部，將中國(guó)大陸明顯分為2 個(gè)區(qū)域，2 個(gè)區(qū)域低點(diǎn)時(shí)間相差約4 h。通過(guò)分析世界地磁中心網(wǎng)站（http://wdc.kugi.kyoto-u.ac.jp/dstdir/index.html）提供的磁暴環(huán)電流Dst指數(shù)和中國(guó)科學(xué)院國(guó)家空間科學(xué)中心網(wǎng)站（http://www.sepc.ac.cn/）提供的Kp指數(shù)發(fā)現(xiàn)，8 月18 日前后幾天Dst指數(shù)絕對(duì)值均小于20 nT（圖1），Kp指數(shù)均小于5，中國(guó)大陸地磁臺(tái)站沒(méi)有受到明顯磁暴、亞暴現(xiàn)象的影響。因此認(rèn)為，此次低點(diǎn)位移異常真實(shí)存在，符合地磁低點(diǎn)位移法的判別方法。

圖1 2019 年8 月磁暴環(huán)電流Dst 指數(shù)變化Fig.1 Magnetic storm ring current index variation curve in August 2019

對(duì)2019 年8 月的低點(diǎn)時(shí)間進(jìn)行梯度法計(jì)算。掃描結(jié)果顯示，無(wú)地磁低點(diǎn)位移異常時(shí)中國(guó)大陸大部分區(qū)域梯度值小于1.0 h/（°）。如圖2 所示，以8 月10 日為例，中國(guó)大陸地磁臺(tái)站低點(diǎn)時(shí)間空間分布沒(méi)有出現(xiàn)明顯的分區(qū)特征，同時(shí)也沒(méi)有出現(xiàn)明顯的高梯度帶。8 月18 日出現(xiàn)地磁低點(diǎn)位移異常當(dāng)天，同步出現(xiàn)了大于1.0 h/（°）的高梯度帶，其局部地區(qū)梯度值大于2 h/（°），高梯度帶主要分布在甘肅中部地區(qū)（圖3）。9 月16 日張掖5.0 級(jí)地震震中在高梯度帶邊緣1.0 h/（°）等值線附近，該處也是低點(diǎn)位移線出現(xiàn)明顯轉(zhuǎn)折的位置。

圖2 正常日中國(guó)大陸地磁低點(diǎn)時(shí)間、梯度分布（以2019 年8 月10 日為例）Fig.2 Distribution of low-point time and gradient of geomagnetism in Chinese mainland at nonanomalous date (on 10 August 2019)

圖3 2019 年8 月18 日中國(guó)大陸地磁低點(diǎn)位移線分布及低點(diǎn)時(shí)間梯度分布特征Fig.3 Distribution characteristics of low-point displacement line and low-point time’s gradient in Chinese mainland on 18 August 2019

2.2 夏河5.7 級(jí)地震前地磁低點(diǎn)位移特征

2019 年10 月28 日甘肅夏河5.7 級(jí)地震前42 天，即9 月16 日中國(guó)大陸出現(xiàn)明顯的地磁低點(diǎn)位移現(xiàn)象，低點(diǎn)位移線經(jīng)過(guò)西藏自治區(qū)、青海省、穿過(guò)甘肅省東南、經(jīng)過(guò)寧夏回族自治區(qū)、陜西省，在內(nèi)蒙古自治區(qū)中部結(jié)束，將中國(guó)大陸明顯分為2 個(gè)區(qū)域，2 個(gè)區(qū)域低點(diǎn)時(shí)間相差約2 h。通過(guò)分析磁暴環(huán)電流Dst指數(shù)和Kp指數(shù)發(fā)現(xiàn)，9 月16 日前后幾天，Dst指數(shù)絕對(duì)值均小于20 nT（圖4），Kp指數(shù)均小于5，中國(guó)大陸地磁臺(tái)站沒(méi)有受到明顯磁暴、亞暴現(xiàn)象的影響。因此認(rèn)為，此次低點(diǎn)位移異常真實(shí)存在，符合地磁低點(diǎn)位移法的判別方法。

圖4 2019 年9 月1—29 日磁暴環(huán)電流Dst 指數(shù)變化Fig.4 Magnetic storm ring current index variation curve in September 2019

對(duì)2019 年9 月的低點(diǎn)時(shí)間進(jìn)行梯度法計(jì)算。掃描結(jié)果顯示，無(wú)地磁低點(diǎn)位移異常時(shí)中國(guó)大陸沒(méi)有出現(xiàn)梯度值大于1.0 h/（°）的區(qū)域。如圖5 所示，以9 月20 日為例，中國(guó)大陸地磁臺(tái)站低點(diǎn)時(shí)間空間分布沒(méi)有出現(xiàn)明顯的分區(qū)特征，同時(shí)也沒(méi)有出現(xiàn)明顯的高梯度帶。9 月16 日出現(xiàn)地磁低點(diǎn)位移異常當(dāng)天，同步出現(xiàn)了梯度值大于1.0 h/（°）的高梯度帶，其主要分布在甘肅中東部地區(qū)（圖6）。10 月28 日夏河5.7 級(jí)地震震中在高梯度帶邊緣1.0 h/（°）等值線附近，該處也是低點(diǎn)位移線出現(xiàn)明顯轉(zhuǎn)折的位置。

圖5 正常日中國(guó)大陸地磁低點(diǎn)時(shí)間分布和梯度分布特征（以2019 年9 月20 日為例）Fig.5 Distribution of low-point time and gradient of geomagnetism in Chinese mainland on nonanomalous date (20 September 2019)

圖6 2019 年9 月16 日中國(guó)大陸地磁低點(diǎn)位移線分布及低點(diǎn)時(shí)間梯度分布Fig.6 Distribution of low-point displacement line and low-point time’s gradient in Chinese mainland on 16 September 2019

2.3 異常日地磁垂直分量日變化特征

根據(jù)低點(diǎn)位移形成機(jī)理認(rèn)為，異常日臺(tái)站地磁垂直分量日變化形態(tài)畸變是由于感應(yīng)電流重新分布而造成的（馮志生等，2009；戴勇等，2019）。異常日在低點(diǎn)位移線下方短時(shí)間形成高導(dǎo)通道，感應(yīng)電流在高導(dǎo)通道內(nèi)聚集產(chǎn)生的磁場(chǎng)造成低點(diǎn)位移線兩側(cè)臺(tái)站地磁垂直分量日變化形態(tài)畸變和低點(diǎn)時(shí)間偏移。

以2019 年8 月18 日低點(diǎn)位移異常為例，低點(diǎn)位移線兩側(cè)200 km 范圍內(nèi)，山丹地震臺(tái)、湟源地震臺(tái)、固原地震臺(tái)、中衛(wèi)黑山嘴地震臺(tái)等的垂直分量明顯出現(xiàn)日變化形態(tài)畸變和低點(diǎn)時(shí)間偏移的現(xiàn)象，具體結(jié)果如圖7 所示。由圖7 可見(jiàn)，低點(diǎn)位移線兩側(cè)臺(tái)站地磁垂直分量日變化形態(tài)均由無(wú)異常日的“V”字型轉(zhuǎn)為“W”型，低點(diǎn)位移線上方固原地震臺(tái)、中衛(wèi)地震臺(tái)、大同地震臺(tái)等低點(diǎn)時(shí)間出現(xiàn)明顯滯后（偏移），低點(diǎn)位移線下方低點(diǎn)時(shí)間基本在地方時(shí)12 時(shí)左右，未出現(xiàn)明顯偏移。

圖7 2019 年8 月18 日地磁低點(diǎn)位移線兩側(cè)200 km 范圍內(nèi)臺(tái)站地磁垂直分量日變化Fig.7 Daily variation pattern of geomagnetic vertical component of stations within 200 km on both sides of low-point displacement line on 18 August 2019

綜上所述，低點(diǎn)位移異常日在低點(diǎn)位移線附近臺(tái)站地磁垂直分量日變化形態(tài)出現(xiàn)明顯畸變，低點(diǎn)時(shí)間相較無(wú)異常日明顯偏移。低點(diǎn)位移線下方感應(yīng)電流形成的磁場(chǎng)在臺(tái)站本身地磁垂直分量上疊加，造成垂直分量日變化形態(tài)由“V”字型轉(zhuǎn)變?yōu)椤癢” 字型。低點(diǎn)位移線轉(zhuǎn)折較明顯的位置有利于感應(yīng)電流形成的磁場(chǎng)互相疊加，容易形成高梯度帶，該處也是中強(qiáng)地震容易發(fā)生的位置。

3 結(jié)論與討論

（1）2019 年8 月18 日中國(guó)大陸地區(qū)出現(xiàn)明顯地磁低點(diǎn)位移現(xiàn)象，在甘肅中部地區(qū)存在明顯地磁低點(diǎn)時(shí)間高梯度帶。此次低點(diǎn)位移發(fā)生后29 天，即9 月16 日發(fā)生張掖5.0 級(jí)地震，地震震中在低點(diǎn)位移線發(fā)生明顯轉(zhuǎn)折處及低點(diǎn)時(shí)間高梯度帶1.0 h/（°）等值線邊緣；9 月16 日中國(guó)大陸地區(qū)再次出現(xiàn)明顯地磁低點(diǎn)位移現(xiàn)象，在甘肅東南部地區(qū)存在明顯地磁低點(diǎn)時(shí)間高梯度帶。此次低點(diǎn)位移發(fā)生后42 天，即10 月18 日甘肅夏河發(fā)生5.7 級(jí)地震，地震震中同樣在低點(diǎn)位移線發(fā)生明顯轉(zhuǎn)折處及低點(diǎn)時(shí)間高梯度帶1.0 h/（°）等值線邊緣。通過(guò)以上2 次震例的分析認(rèn)為，在甘肅中東部地區(qū)，出現(xiàn)低點(diǎn)位移的情況下對(duì)低點(diǎn)時(shí)間進(jìn)行梯度分析能夠有效縮小發(fā)震地點(diǎn)預(yù)測(cè)范圍，高梯度帶邊緣可能是易發(fā)生中強(qiáng)地震的位置。

（2）低點(diǎn)位移異常日低點(diǎn)位移線200 km 范圍內(nèi)臺(tái)站地磁垂直分量日變化形態(tài)、日變化幅度、低點(diǎn)時(shí)間等相較無(wú)異常日都有明顯變化。低點(diǎn)位移線下方感應(yīng)電流形成的磁場(chǎng)在臺(tái)站地磁垂直分量上疊加，當(dāng)二者方向相同時(shí)，造成垂直分量日變化形態(tài)“V”字型幅度變大，當(dāng)二者方向相反時(shí)，造成垂直分量日變化形態(tài)由“V”字型轉(zhuǎn)變?yōu)椤癢”字型；日變化幅度由于感應(yīng)電流產(chǎn)生的磁場(chǎng)疊加相較無(wú)異常日呈現(xiàn)明顯的偏大或偏小；低點(diǎn)位移線兩側(cè)低點(diǎn)時(shí)間變化不一致，其中，某一側(cè)相較無(wú)異常日出現(xiàn)明顯偏移，而另外一側(cè)低點(diǎn)時(shí)間未發(fā)生明顯變化。2 個(gè)震例中低點(diǎn)位移線上方的臺(tái)站低點(diǎn)時(shí)間明顯滯后于當(dāng)?shù)貢r(shí)間12 h，而低點(diǎn)位移線下方臺(tái)站低點(diǎn)時(shí)間變化不大。低點(diǎn)位移線轉(zhuǎn)折較明顯的位置有利于感應(yīng)電流形成的磁場(chǎng)互相疊加，容易形成高梯度帶，該處也是中強(qiáng)地震容易發(fā)生的位置。

（3）每年8—9 月中國(guó)大陸地磁低點(diǎn)時(shí)間分布掃描結(jié)果顯示，存在多次低點(diǎn)位移線經(jīng)過(guò)甘肅中東部地區(qū)的異常。其中，僅有2019 年8 月18 日、9 月16 日的低點(diǎn)位移異常同時(shí)伴有低點(diǎn)時(shí)間高梯度帶存在，而眾多低點(diǎn)位移異常中僅有這2 次異常分別對(duì)應(yīng)了張掖5.0 級(jí)、夏河5.7 級(jí)地震。因此，對(duì)低點(diǎn)時(shí)間進(jìn)行梯度分析能夠改進(jìn)地磁低點(diǎn)位移方法在地震預(yù)測(cè)中的應(yīng)用，有利于減少虛報(bào)率。