一階差分法在地震地磁前兆信號(hào)提取中的應(yīng)用

袁桂平 馮志生 李鴻宇 姜楚峰 潘穎 錢庚 姚休義

摘要：在探索震前地磁每日一值空間相關(guān)異常機(jī)理的基礎(chǔ)上，利用2008—2018年130個(gè)臺(tái)站地磁場(chǎng)數(shù)據(jù)計(jì)算地磁場(chǎng)總強(qiáng)度F每日一值一階差分在空間上變化的不同步現(xiàn)象，并提取震前地磁異常。結(jié)果表明：當(dāng)?shù)卮艌?chǎng)總強(qiáng)度F逐漸增大這一趨勢(shì)變化出現(xiàn)不同步現(xiàn)象時(shí)，研究區(qū)內(nèi)一階差分值全正或者全負(fù)的正常狀態(tài)會(huì)被打破，從而導(dǎo)致區(qū)域內(nèi)出現(xiàn)各臺(tái)站F每日一值一階差分正、負(fù)值臺(tái)站分布較為集中的現(xiàn)象，且正、負(fù)值臺(tái)站之間出現(xiàn)明顯的異常分界線（0值線），通過(guò)異常與對(duì)應(yīng)震例統(tǒng)計(jì)，發(fā)現(xiàn)該異常發(fā)生后8個(gè)月內(nèi)在0值線300 km附近有發(fā)生5級(jí)以上地震的可能性。

關(guān)鍵詞：地磁異常;地磁場(chǎng)總強(qiáng)度;一階差分;前兆異常

中圖分類號(hào)：P315.721文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2020）04-0738-07

0引言

地震的孕育和發(fā)生是一個(gè)復(fù)雜的過(guò)程，伴隨這一過(guò)程產(chǎn)生的電磁現(xiàn)象也極其復(fù)雜。國(guó)內(nèi)外學(xué)者研究發(fā)現(xiàn)許多關(guān)于地震引起地下電導(dǎo)率變化的實(shí)驗(yàn)和野外觀測(cè)結(jié)果（Nagata，1972;Brace，Orange，1968;Barsukov，1972），這些結(jié)果引導(dǎo)人們研究地震的感應(yīng)磁效應(yīng)現(xiàn)象（中國(guó)科學(xué)院地球物理研究所第十研究室二組，1977;祁貴仲等，1981）。在相關(guān)系數(shù)異常提取方面，前人利用華北地區(qū)若干臺(tái)站資料、以局部地區(qū)地磁場(chǎng)空間相關(guān)性的分析為基礎(chǔ)，重點(diǎn)討論了“空間相關(guān)法”，并通過(guò)該方法及相關(guān)震例討論了震磁效應(yīng)的有效性（中國(guó)科學(xué)院地球物理研究所第十研究室一組，1977）;林美和沈斌（1982）認(rèn)為相關(guān)系數(shù)法是日變形態(tài)的一種發(fā)展，在地磁日變曲線產(chǎn)生畸變時(shí)，相關(guān)系數(shù)法更顯優(yōu)越性，且驗(yàn)證了該方法應(yīng)用于地震預(yù)報(bào)有一定的實(shí)際意義。在異常的定量研究方面，杜安娜和宋若薇（1982）發(fā)現(xiàn)1976年唐山7.8級(jí)地震前出現(xiàn)了相關(guān)系數(shù)低值異常，該相關(guān)系數(shù)異常量超過(guò)2.5倍均方差;1991年大同5.8級(jí)地震發(fā)生后，曾小蘋等（1992）發(fā)現(xiàn)在震中約220 km內(nèi)的地磁垂直分量Z在震前3個(gè)月和5天有明顯異常。后續(xù)其他學(xué)者在利用空間相關(guān)系數(shù)提取地震地磁異常時(shí)，也進(jìn)行了許多的研究和驗(yàn)證，針對(duì)不同地區(qū)提出了對(duì)應(yīng)的指標(biāo)，并累積了許多的震例（馮志生，王建宇，1998;賈立峰等，2017;李鴻宇等，2017，2018）。在異常機(jī)理方面，杜安娜（1998）應(yīng)用空間線性相關(guān)法對(duì)武定、麗江地震進(jìn)行了預(yù)報(bào)分析，認(rèn)為在這2次地震孕育過(guò)程中，在地應(yīng)力的作用下，通海、西昌、成都等地磁臺(tái)附近巖石電導(dǎo)率發(fā)生變化;趙凌強(qiáng)等（2015）對(duì)2013年甘肅岷縣漳縣MS6.6地震孕震環(huán)境探討中認(rèn)為松潘—甘孜地塊和隴西盆地對(duì)西秦嶺造山帶形成擠壓、阻擋作用，其能量可能通過(guò)西秦嶺造山帶中上地殼的高電阻構(gòu)造傳遞，被臨潭—宕昌斷裂帶附近低阻帶所吸收，當(dāng)這種從南向北擠壓和阻擋持續(xù)作用超過(guò)該地區(qū)介質(zhì)的應(yīng)力臨界值時(shí)，該區(qū)臨潭—宕昌低阻破碎帶顯示出了不穩(wěn)定性而發(fā)生形變，導(dǎo)致該區(qū)附近的高電阻特性的巖石產(chǎn)生破裂或?qū)娱g滑動(dòng)，進(jìn)而發(fā)生了岷縣漳縣地震。而高電阻特性的巖石產(chǎn)生破裂或?qū)娱g滑動(dòng)必然造成巖石電導(dǎo)率變化，從而導(dǎo)致感應(yīng)磁場(chǎng)發(fā)生變化。因此，地磁場(chǎng)記錄到觀測(cè)值除正常磁場(chǎng)外還包含了這一部分異常變化，這也是震前地磁場(chǎng)異常信號(hào)產(chǎn)生的根本原因。

眾多研究成果證明每日一值空間相關(guān)系數(shù)確實(shí)能提取震前地磁異常，但缺少絕對(duì)標(biāo)準(zhǔn)的臺(tái)站作為參考臺(tái)，計(jì)算結(jié)果會(huì)有一定的誤差，導(dǎo)致在異常分析判別中對(duì)異常臺(tái)站的選取有較大的困難和人為主觀性。追溯源頭，每日一值空間相關(guān)系數(shù)的降低其實(shí)是參考臺(tái)站和計(jì)算臺(tái)站的同步趨勢(shì)變化被打破所導(dǎo)致，即參考臺(tái)不變的情況下，計(jì)算臺(tái)站數(shù)據(jù)在正?；A(chǔ)上有所上升或者下降，這種變化可通過(guò)一階差分表現(xiàn)出來(lái)，且不用涉及到參考臺(tái)的選取。張素琴等（2015）利用2010—2014年我國(guó)地磁臺(tái)網(wǎng)100多個(gè)臺(tái)站的地磁場(chǎng)總強(qiáng)度F子夜均值逐日差值，研究其空間變化趨勢(shì)時(shí)發(fā)現(xiàn)，當(dāng)逐日差值在小區(qū)域內(nèi)異于周圍大趨勢(shì)增減變化的18次變化中，共有16次異常與其后幾天到半年內(nèi)該區(qū)域發(fā)生的地震有很好的對(duì)應(yīng)關(guān)系。這也說(shuō)明逐日差值法在震前地磁信號(hào)處理中有一定的效果。

地震研究43卷第4期袁桂平等：一階差分法在地震地磁前兆信號(hào)提取中的應(yīng)用本文擬通過(guò)每日一值一階差分法計(jì)算2008—2018年130個(gè)臺(tái)站記錄的世界時(shí)16～20時(shí)地磁場(chǎng)總強(qiáng)度F整點(diǎn)值的均值數(shù)據(jù)的逐日差分值，尋找異?，F(xiàn)象，并分析其和地震之間的對(duì)應(yīng)關(guān)系。1資料選取及研究方法

1.1資料選取

本文研究區(qū)為中國(guó)大陸區(qū)域（75°～129°E，19°～50°N），以各臺(tái)站儀器記錄數(shù)據(jù)的連續(xù)性和穩(wěn)定性較好為條件，從地磁臺(tái)網(wǎng)范圍內(nèi)，選取130個(gè)臺(tái)站（圖1），其中配置FHD質(zhì)子矢量磁力儀的臺(tái)站共80個(gè)，配置GSM19FD的臺(tái)站共3個(gè)，配置GSM90F1的臺(tái)站共22個(gè)，配置M15的臺(tái)站共25個(gè)。

本文采用的數(shù)據(jù)為2008—2018年130個(gè)臺(tái)站記錄的世界時(shí)16～20時(shí)（地方時(shí)0～4時(shí)）地磁場(chǎng)總強(qiáng)度F整點(diǎn)值的均值數(shù)據(jù)。由于夜間地磁場(chǎng)受空間磁場(chǎng)影響和環(huán)境干擾相對(duì)較小，其變化比其它時(shí)段更為穩(wěn)定，每個(gè)整點(diǎn)值均為該整點(diǎn)值前后5 min的分鐘值均值，再取16～20時(shí)的整點(diǎn)值的均值能夠保證數(shù)據(jù)的穩(wěn)定性，且可避免干擾和計(jì)算量過(guò)于冗繁。若有缺數(shù)，則將該天標(biāo)為缺數(shù)，該天的實(shí)際參與計(jì)算的臺(tái)站是130減去缺數(shù)臺(tái)站數(shù)，若當(dāng)天缺數(shù)臺(tái)站過(guò)多，其結(jié)果視為無(wú)效。

1.2F每日一值一階差分計(jì)算

地磁場(chǎng)由穩(wěn)定磁場(chǎng)和變化磁場(chǎng)組成。F夜間每日一值的逐日變化基本可以反映出地磁穩(wěn)定磁場(chǎng)數(shù)十年甚至上萬(wàn)年的長(zhǎng)期變化，且在一定的空間范圍內(nèi)有很好的相關(guān)性，即同步變化現(xiàn)象。變化磁場(chǎng)反映的是地磁場(chǎng)的短期變化現(xiàn)象。地震地磁異常可導(dǎo)致2個(gè)臺(tái)站之間的F每日一值空間相關(guān)系數(shù)下降，這種現(xiàn)象即為地磁短期變化所引起。相關(guān)系數(shù)的降低是由各個(gè)臺(tái)站記錄到的地磁數(shù)據(jù)不同步變化引起的，該現(xiàn)象可從F每日一值一階差分上進(jìn)行提取，從而提取地震地磁異常信號(hào)。即正常情況下區(qū)域地磁場(chǎng)是同步變化的，F(xiàn)每日一值一階差分值應(yīng)為同正或同負(fù)，當(dāng)出現(xiàn)磁異常時(shí)，F(xiàn)每日一值一階差分值正、負(fù)可能會(huì)同時(shí)出現(xiàn)，且正、負(fù)值區(qū)域間有明顯的現(xiàn)象分界線。

本文利用2008—2018年F的每日一值（世界時(shí)16，17，18，19，20時(shí)的整點(diǎn)值數(shù)據(jù)的均值）進(jìn)行一階差分計(jì)算，獲得前后兩日之間F的差值為：dF=F（t2）-F（t1）（1）式中：F（t2）為觀測(cè)日期t2的F（世界時(shí)16，17，18，19，20時(shí)的整點(diǎn)值數(shù)據(jù)的均值）的值;t2=t1+1;異常日期以F（t2）日期為準(zhǔn)。

2計(jì)算結(jié)果及震例分析

2.1研究機(jī)理

地磁主磁場(chǎng)長(zhǎng)期變化的時(shí)間特征可由長(zhǎng)期變時(shí)間譜看出，其時(shí)間尺度為數(shù)十年甚至上萬(wàn)年（徐文耀，2003）。F的長(zhǎng)期變化是相對(duì)較穩(wěn)定且緩慢的，王振東等（2017，2019）利用CHAOS-6和IGRF-12模型對(duì)2010—2016年14個(gè)地磁臺(tái)站的觀測(cè)數(shù)據(jù)描述的F都呈現(xiàn)逐漸增大趨勢(shì)。而F夜間每日一值的逐日變化基本可以將這個(gè)近似于線性增大的長(zhǎng)期變化反映出來(lái)。因此，正常情況下，研究區(qū)內(nèi)磁場(chǎng)趨勢(shì)上升/下降時(shí)，各臺(tái)站F每日一值一階差分為正/負(fù)值。只有當(dāng)研究區(qū)內(nèi)局部干擾源對(duì)磁場(chǎng)產(chǎn)生影響、導(dǎo)致區(qū)域內(nèi)部分臺(tái)站磁場(chǎng)變化不同步時(shí)，各臺(tái)站一階差分值才會(huì)出現(xiàn)正、負(fù)值同時(shí)出現(xiàn)。當(dāng)正、負(fù)臺(tái)站相對(duì)集中，將研究區(qū)在空間上劃分為2～3個(gè)區(qū)域時(shí)，說(shuō)明有較大電流通過(guò)該區(qū)域，導(dǎo)致磁場(chǎng)發(fā)生有規(guī)律的變化，即為F每日一值一階差分異常。本文以2013年5月19日各臺(tái)站F每日一值一階差分dF為示意做等值線圖，從圖2可見，dF在35°N左右以北為正、以南為負(fù)，正、負(fù)值臺(tái)站間有明顯的分界線（0值線），且將空間區(qū)域劃分為2個(gè)區(qū)域，就是本文所研究的地磁場(chǎng)總強(qiáng)度F每日一值一階差分異?，F(xiàn)象。該異常出現(xiàn)后的64天和211天，在dF的0值線附近分別發(fā)生了2013年7月22日岷縣6.6級(jí)地震和2013年12月16日湖北巴東5.1級(jí)地震。

dF的0值線出現(xiàn)是由于0值線兩側(cè)臺(tái)站dF變化不同步所致。為看清這種變化細(xì)節(jié)，筆者取0值線兩側(cè)的烏魯木齊—拉薩、銀川—邵陽(yáng)、錫林浩特—泉州等3組臺(tái)站（圖2中黑線所指臺(tái)站）的dF數(shù)據(jù)進(jìn)行分析研究，其2013年5月7—29日的逐日變化曲線如圖3所示。從圖中可以看到，各臺(tái)站dF在0值附近有輕微波動(dòng)，且基本同步。5月19日，烏魯木齊、銀川和錫林浩特臺(tái)dF同為正值（紅點(diǎn)），均位于dF正值區(qū)，而拉薩、邵陽(yáng)和泉州臺(tái)dF同為負(fù)值（藍(lán)點(diǎn)），均位于dF負(fù)值區(qū)。但如果不進(jìn)行正負(fù)標(biāo)識(shí)，僅從圖3很難看出這種差異，這也是繪制空間等值線的原因。

研究區(qū)內(nèi)各臺(tái)站dF出現(xiàn)變化不同步的現(xiàn)象，那么各臺(tái)站在日變化上必然是有所差別的。因此，取2013年5月17—20日烏魯木齊—拉薩、錫林浩特—泉州2個(gè)臺(tái)站對(duì)的地磁場(chǎng)總強(qiáng)度F日變化分鐘值（去均值）曲線（圖4）進(jìn)行詳細(xì)分析，烏魯木齊、錫林浩特臺(tái)位于0值線北側(cè)，拉薩、泉州要位于0值線南側(cè)。從圖4可以看到，5月17日和20日世界時(shí)16～20時(shí)（圖中黑框）2個(gè)臺(tái)站對(duì)數(shù)據(jù)基本重疊，而18，19日數(shù)據(jù)重疊效果較差，變化不同步，導(dǎo)致19日dF在0值線兩側(cè)出現(xiàn)一正一負(fù)的現(xiàn)象。

表1為2013年5月18—20日所選烏魯木齊—拉薩、銀川—邵陽(yáng)、錫林浩特—泉州3個(gè)臺(tái)站對(duì)的dF數(shù)據(jù)，18日6個(gè)臺(tái)站dF均為負(fù)值，20日均為正值，而19日位于0值線北側(cè)臺(tái)站dF為正，其南側(cè)臺(tái)站dF為負(fù)，與圖4中展示的結(jié)果一致，但更為清晰。19日的dF是利用19日數(shù)據(jù)減去18日數(shù)據(jù)獲取的，且18，19日磁場(chǎng)數(shù)據(jù)小于17和20日同時(shí)段數(shù)據(jù)（圖4），因此綜合分析認(rèn)為18，19日有附加磁信號(hào)，對(duì)0值線南側(cè)磁場(chǎng)的削弱效果大于其北側(cè)，造成區(qū)域出現(xiàn)一階差分異?，F(xiàn)象。此外，涉縣—信陽(yáng)、蘭州—成都、湟源—道孚、涇陽(yáng)—鐘祥等臺(tái)站數(shù)據(jù)亦出現(xiàn)同樣的現(xiàn)象，但由于這幾個(gè)臺(tái)站對(duì)距離0值線較近，其一階差分異常現(xiàn)象沒(méi)有圖4所展示得明顯。

2.2異常與地震對(duì)應(yīng)結(jié)果

異常和震例統(tǒng)計(jì)中，時(shí)間上，以8個(gè)月為期，異常出現(xiàn)后最短發(fā)震時(shí)間為0天，最長(zhǎng)為240天?？紤]到8個(gè)月時(shí)間覆蓋期較長(zhǎng)，且根據(jù)統(tǒng)計(jì)結(jié)果發(fā)現(xiàn)有震異常中6個(gè)月內(nèi)發(fā)震數(shù)占總發(fā)震數(shù)目的73%，是一個(gè)較為合理的占比，因此認(rèn)為發(fā)震優(yōu)勢(shì)時(shí)段為異常出現(xiàn)后6個(gè)月內(nèi)。地點(diǎn)上，結(jié)合磁異常震例分析認(rèn)為距0值線300 km較為合適。震級(jí)上，由于地磁場(chǎng)變化具有區(qū)域效應(yīng)，小震對(duì)地磁場(chǎng)影響較小且磁場(chǎng)變化混亂，因此對(duì)震級(jí)選擇較高，以105°E為界限，不同地區(qū)震級(jí)范圍分別為：大陸東部MS≥5.0地震、大陸西部及國(guó)界周邊MS≥6.0地震、臺(tái)灣地區(qū)MS≥6.5地震。幅度上（正值臺(tái)站差分均值減去負(fù)值臺(tái)站差分均值），對(duì)當(dāng)日1～4 nT dF幅度值進(jìn)行詳細(xì)分析發(fā)現(xiàn)，小于3 nT時(shí)一階差分等值線紊亂，大于3 nT空間上無(wú)變化，但提取異常天數(shù)急劇減小，因此幅度以3 nT為標(biāo)準(zhǔn)。

按照以上判據(jù)標(biāo)準(zhǔn)，統(tǒng)計(jì)研究區(qū)2008—2018年共4 017天內(nèi)，滿足判據(jù)的dF異常共有174組，占總天數(shù)的4%，排除57組負(fù)值臺(tái)站分布散亂未能找到明顯分界線的天數(shù)，共有117組滿足時(shí)間判據(jù)標(biāo)準(zhǔn)，僅占總天數(shù)的3%，說(shuō)明此現(xiàn)象不是頻繁出現(xiàn)的。117組異常中，在異常出現(xiàn)后8個(gè)月內(nèi)有滿足震級(jí)和距離判據(jù)標(biāo)準(zhǔn)的地震發(fā)生有96組異常，異常和地震的對(duì)應(yīng)率為82%。滿足判據(jù)標(biāo)準(zhǔn)的地震共90次，55次地震前有dF異?，F(xiàn)象出現(xiàn)，漏報(bào)地震35次，漏報(bào)率為39%;大陸區(qū)域滿足判據(jù)標(biāo)準(zhǔn)的地震共51次，漏報(bào)地震12次，漏報(bào)率為24%;有異常對(duì)應(yīng)的55次地震中，僅有6次地震為海域地震，占報(bào)對(duì)地震總數(shù)的10%，其余異常對(duì)應(yīng)地震發(fā)生于中國(guó)大陸區(qū)域或者國(guó)界附近。這是由于大陸區(qū)域臺(tái)站分布相對(duì)密集，對(duì)此信號(hào)捕捉能力較強(qiáng)，而海域僅沿海臺(tái)站對(duì)其有指示意義，對(duì)信號(hào)的捕捉能力較弱。

2.3震例分析

挑選具有代表性的幾個(gè)異常震例進(jìn)行詳細(xì)分析，包括2008年5月12日汶川8.0級(jí)地震、2016年7月31日廣西蒼梧5.4級(jí)地震、2018年9月12日陜西寧強(qiáng)5.3級(jí)地震，對(duì)震前F每日一值一階差分異常現(xiàn)象進(jìn)行分析。

（1）2008年4月4日中國(guó)大陸地區(qū)出現(xiàn)F每日一值一階差分正、負(fù)值臺(tái)站分布較為集中的現(xiàn)象，正、負(fù)臺(tái)站分界線（0值線）大致為NW—ES走向，如圖5a所示。2008年5月12日及2008年11月10日在分界線邊緣分別發(fā)生汶川8.0級(jí)地震和青海海西6.3級(jí)地震，2個(gè)地震震中均位于分界線附近300 km的范圍內(nèi)，發(fā)震日期距離異常日的天數(shù)分別38天和220天。

（2）2016年4月3日中國(guó)大陸地區(qū)出現(xiàn)F每日一值一階差分正、負(fù)值臺(tái)站分布較為集中的現(xiàn)象，正、負(fù)值臺(tái)站分界線大約為南北轉(zhuǎn)東走勢(shì)，東南部地區(qū)一階差分為負(fù)值，其余地區(qū)為正值，如圖5b所示。2016年7月31日廣西蒼梧5.4級(jí)地震震中位于正、負(fù)分界線上，且分界線兩側(cè)一階差分正、負(fù)值臺(tái)站均勻分布，發(fā)震日期和異常日間隔為94天。

（3）2018年7月14日中國(guó)大陸地區(qū)出現(xiàn)F每日一值一階差分正、負(fù)值臺(tái)站分布較為集中的現(xiàn)象，正、負(fù)臺(tái)站分界線（0值線）在中緯度地區(qū)呈W型走勢(shì)，且0值線兩側(cè)一階差分正、負(fù)值臺(tái)站均勻分布，如圖5c所示。2018年9月12日陜西寧強(qiáng)5.3級(jí)地震震中位于0值線邊緣，發(fā)震日期距離異常日天數(shù)為60天。

3討論

地磁場(chǎng)總強(qiáng)度F每日一值一階差分異?，F(xiàn)象并不是頻發(fā)現(xiàn)象，僅占研究總天數(shù)的3%;異常現(xiàn)象發(fā)生之后0值線附近確實(shí)有滿足判據(jù)標(biāo)準(zhǔn)的地震發(fā)生，報(bào)對(duì)率為82%;陸地漏報(bào)率遠(yuǎn)遠(yuǎn)小于全區(qū)漏報(bào)率，主要是因?yàn)樗x地磁臺(tái)站均在陸地，對(duì)陸地指示意義更明確。因此分析認(rèn)為震前的地震地磁異常信號(hào)確實(shí)能通過(guò)F每日一值一階差分表現(xiàn)出來(lái)。

F每日一值一階差分法是一種中短期預(yù)測(cè)方法，時(shí)間上從以8個(gè)月為異常有效期的發(fā)震結(jié)果統(tǒng)計(jì)發(fā)現(xiàn)，異常出現(xiàn)后6個(gè)月內(nèi)發(fā)震居多，占報(bào)對(duì)地震總數(shù)的73%，但無(wú)法做更精確的時(shí)間判斷。地點(diǎn)上，0值線雖然對(duì)未來(lái)的發(fā)震位置具有較強(qiáng)的指示意義，地震一般發(fā)生在0值線附近300 km左右范圍內(nèi)，但由于0值線較長(zhǎng)，因此需要結(jié)合其他前兆手段，才可詳細(xì)分析地震發(fā)生在0值線的具體位置。震級(jí)上，由于磁場(chǎng)本身就是存在較小波動(dòng)，小震產(chǎn)生的磁異?，F(xiàn)象在這種波動(dòng)下會(huì)被淹沒(méi)，難以辨別。因此，此方法所指示地震的震級(jí)相對(duì)較高，且地震震級(jí)與0值線走向和長(zhǎng)短沒(méi)有明顯的對(duì)應(yīng)關(guān)系。

在震前確實(shí)存在磁場(chǎng)在不同區(qū)域的不同步變化或變化量不一致的現(xiàn)象，說(shuō)明同一天在外空磁場(chǎng)差別不大的情況下，不同區(qū)域磁場(chǎng)的不同步變化極有可能是受到地下磁場(chǎng)局部變化所導(dǎo)致。由于地下電性結(jié)構(gòu)變化導(dǎo)致產(chǎn)生的感應(yīng)磁場(chǎng)有所變化，從而疊加在地磁場(chǎng)上，導(dǎo)致F每日一值一階差分出現(xiàn)異常。也可能是異常當(dāng)日0值線附近有感應(yīng)電流通過(guò)，該電流產(chǎn)生的磁場(chǎng)對(duì)0值線兩側(cè)磁場(chǎng)產(chǎn)生的影響不一致，從而導(dǎo)致了0值線兩側(cè)異?，F(xiàn)象出現(xiàn)，而電流的出現(xiàn)亦由地下電性結(jié)構(gòu)變化所致。因此，該異常現(xiàn)象是震前地下電性結(jié)構(gòu)不同步變化的一種反應(yīng)。

由于地磁場(chǎng)受空間天氣的干擾，如磁暴會(huì)導(dǎo)致地磁場(chǎng)發(fā)生劇烈的擾動(dòng)現(xiàn)象。雖然磁暴對(duì)地磁場(chǎng)的擾動(dòng)現(xiàn)象在地磁水平分量上表現(xiàn)更為明顯，但地磁場(chǎng)每一個(gè)分量之間都有著必然的聯(lián)系。因此，我們還對(duì)異常前后兩日磁暴指數(shù)Dst進(jìn)行了統(tǒng)計(jì)分析。結(jié)果表明，F(xiàn)每日一值一階差分異常現(xiàn)象和磁暴并沒(méi)有直接的聯(lián)系，對(duì)研究結(jié)果無(wú)影響，因篇幅較長(zhǎng)，文中沒(méi)有再進(jìn)行詳細(xì)列舉。

4結(jié)論

本文通過(guò)計(jì)算2008—2018年中國(guó)大陸區(qū)域內(nèi)130個(gè)臺(tái)站地磁場(chǎng)總強(qiáng)度F每日一值一階差分并分析其與地震之間的對(duì)應(yīng)關(guān)系，得到以下結(jié)論：

（1）當(dāng)?shù)卮艌?chǎng)總強(qiáng)度F逐漸增大這一趨勢(shì)變化出現(xiàn)不同步現(xiàn)象時(shí)，研究區(qū)內(nèi)一階差分值全正或者全負(fù)的正常狀態(tài)會(huì)被打破。

（2）區(qū)域內(nèi)出現(xiàn)各臺(tái)站F每日一值一階差分正、負(fù)值臺(tái)站分布較為集中的現(xiàn)象，且正、負(fù)區(qū)域間有明顯的分界線（0值線），為地震前兆異常。

（3）異常發(fā)生后8個(gè)月內(nèi)，0值線300 km附近有發(fā)生5級(jí)以上地震的可能性。

參考文獻(xiàn)：

杜安娜，宋若薇.1982.總磁場(chǎng)強(qiáng)度的空間相關(guān)性及唐山地震的異常[J].地震科學(xué)研究，（3）：50-52.

杜安娜.1998.地磁場(chǎng)總強(qiáng)度空間線性相關(guān)在預(yù)報(bào)武定、麗江強(qiáng)震中的應(yīng)用及機(jī)理[J].華南地震，18（3）：308-312.

馮志生，王建宇.2000.江蘇地區(qū)地磁Z21測(cè)值空間相關(guān)異常及其標(biāo)志體系[J].地震地磁觀測(cè)與研究，21（1）：44-49.

賈立峰，喬子云，張國(guó)苓，等.2017.2013年遼寧燈塔M5.1地震地磁異常變化特征[J].地震研究，40（3）：437-443.

李鴻宇，袁桂平，王俊菲，等.2017. 2013年7月22日甘肅岷縣—漳縣6.6級(jí)地震地磁總場(chǎng)F空間相關(guān)異常分析[J].地震工程學(xué)報(bào)，39（3） ：552-556.

李鴻宇，朱培育，王維，等.2018.2013年前郭5.8級(jí)震群的地磁多方法異常分析[J].地震研究，41（1）：111-117.

林美，沈斌.1982.地磁場(chǎng)垂直分量相關(guān)分析與地震的對(duì)應(yīng)關(guān)系[J].地震研究，5（2）：212-219.

祁貴仲，侯作中，范國(guó)華，等.1981.地震的感應(yīng)磁效應(yīng)（二）[J].地球物理學(xué)報(bào)，24（3）：296-309.

王振東，顧左文，陳斌，等.2017.CHAOS-6模型描述的中國(guó)地區(qū)地磁長(zhǎng)期變化及誤差分析[J].地震研究，40（3）：404-410.

王振東，王粲，袁潔浩，等.2019.中國(guó)及鄰近地區(qū)地磁長(zhǎng)期變化分析[J].地震研究，42（1）：103-111.

徐文耀.2003.地磁學(xué)[M].北京：地震出版社.

曾小蘋，林云芳，續(xù)春榮，等.1992.1991年3月26日大同5.8級(jí)地震的磁效應(yīng)初探[J].地震地磁觀測(cè)與研究，13（2）：44-52.

張素琴，胡秀娟，何宇飛，等.2015. F子夜均值逐日差空間異常變化與地震的關(guān)系研究[J].地震研究，38（1）：98-104.

趙凌強(qiáng)，詹 艷，趙國(guó)澤，等.2015.基于深部電性結(jié)構(gòu)特征的2013年甘肅岷縣漳縣MS6.6地震孕震環(huán)境探討[J].地震地質(zhì)，37（2） ：541-554.

中國(guó)科學(xué)院地球物理研究所第十研究室二組.1977.地震的感應(yīng)磁效應(yīng)（一）—三維電磁感應(yīng)的數(shù)值理論[J].地球物理學(xué)報(bào)，20（1）：70-80.

中國(guó)科學(xué)院地球物理研究所第十研究室一組.1977.地磁場(chǎng)的空間相關(guān)性及其在地震預(yù)報(bào)中的應(yīng)用[J].地球物理學(xué)報(bào)，20（3）：169-184.

Barsukov O M.1972.Variations of electric resistivity of mountain rocks connected with tectonic causes[J].Tectonophysics，14（3-4）：273-277.

Brace W F，Orange A S.1968.Electrical resistivity changes in saturated rocks during fracture and frictional sliding[J].Journal of Geophysical Research，73（4）：1433-1445.

Nagata T.1972.Application of tectonomagnetism to earthquake phenomena[J].Tectonophysics，14（3-4）：263-271.

Application of First-Order Difference Method in the Extraction

of Seismic Geomagnetism Precursor Signal

YUAN Guiping1，F(xiàn)ENG Zhisheng1，LI Hongyu1，JIANG Chufeng2，

PAN Ying1，QIAN Geng2，YAO Xiuyi3

（1.Jiangsu Earthquake Agency，Nanjing 210014，Jiangsu，China）

（2.Hubei Earthquake Agency，Wuhan 430071，Hubei，China）

（3.Catastrophic Risk Management Research Center，Yunnan University

of Financeand Economics，Kunming 650221，Yunnan，China）

Abstract

On the basis of exploring the spatial correlation anomaly mechanism of geomagnetism daily one-value before earthquake，using the total geomagnetic field data of 130 stations from 2008 to 2018，we calculated the spatial variation of the one-value first-order difference of the total geomagnetic field F and extracted the geomagnetic anomalies before earthquakes.The results show that when the non-synchronous phenomena on the gradual increase of? F was appeared，the normal state of all positive or all negative difference values in the study area will be broken.As a result，there appreaed the distribution of positive and negative value stations whose the first order difference isoline of? F daily value are more concentrated in the region，and there is an obvious abnormal boundary（zero value line）between positive and negative value stations.According to the statistics of anomalies and corresponding earthquake，it is found that there is the possibility of M≥5 earthquakes occurring within 8 months after the occurrence of the anomaly near 300 km of the zero line.

Keywords：geomagnetism anomaly;total geomagnetic intensity;first-order difference;precursor anomaly