華北地區(qū)地磁加卸載響應(yīng)比方法應(yīng)用研究①

賈立峰,張國(guó)苓,喬子云,梁紅杰,張　波

(1.河北省地震局,河北 石家莊 050021; 2.河北省地震局柏舍地震臺(tái),河北 邢臺(tái) 056001)

?

華北地區(qū)地磁加卸載響應(yīng)比方法應(yīng)用研究①

賈立峰1,張國(guó)苓1,喬子云1,梁紅杰1,張波2

(1.河北省地震局,河北 石家莊 050021; 2.河北省地震局柏舍地震臺(tái),河北 邢臺(tái) 056001)

運(yùn)用地磁加卸載響應(yīng)比方法(LURR)對(duì)華北地區(qū)23個(gè)地磁臺(tái)站2008—2014年的數(shù)字化地磁相對(duì)記錄數(shù)據(jù)進(jìn)行計(jì)算分析,結(jié)果顯示華北地區(qū)7年內(nèi)發(fā)生了38次共計(jì)10組加卸載響應(yīng)比異常,其中7組有震例對(duì)應(yīng),在選定的24次震例中16次有震前異常。經(jīng)統(tǒng)計(jì)發(fā)現(xiàn)各次異常發(fā)生時(shí)其異常面積與最大異常強(qiáng)度成正比關(guān)系,各組異常內(nèi)累計(jì)異常臺(tái)站數(shù)量與組內(nèi)地震輻射能量有一定相關(guān)性。對(duì)有震例的7組異常分析后認(rèn)為多數(shù)情況下組內(nèi)各次異常之間有一定聯(lián)系,主要表現(xiàn)為位置的連續(xù)性及繼承性,并且發(fā)現(xiàn)震中位置和異常位置分布有一定關(guān)系。

華北地區(qū); 地磁加卸載響應(yīng)比; 前兆異常; 地震預(yù)報(bào)

0　引言

地磁場(chǎng)是地球本身固有的一種重要物理場(chǎng),它能夠深刻反映地下物質(zhì)信息的變化情況[1-2],同時(shí)也是一種可以用來(lái)研究地震、預(yù)測(cè)地震的前兆手段[3-5]。目前地震系統(tǒng)內(nèi)部使用地磁數(shù)據(jù)進(jìn)行地球內(nèi)部孕震環(huán)境研究及地震預(yù)測(cè)的方法有很多,使用較為廣泛的有地磁日變低點(diǎn)位移法、流動(dòng)地磁分析方法及地磁諧波振幅比方法等[6-8]。

加卸載響應(yīng)比方法(LURR)的理論基礎(chǔ)為非線性系統(tǒng)失穩(wěn)的過(guò)程中伴隨著外部應(yīng)力加載與卸載的不平衡作用過(guò)程[9]。該理論被應(yīng)用于地震活動(dòng)性預(yù)測(cè)當(dāng)中并取得了一定的成果,隨后在地形變、地下水位、重力等前兆資料分析方面也有較好的內(nèi)檢結(jié)果[10]。在地震地磁學(xué)領(lǐng)域,該方法也逐漸成為一種主要的分析方法,常用于地震的中短期預(yù)測(cè)分析。本文使用華北地區(qū)23個(gè)地磁臺(tái)站2008—2014年的地磁Z分量數(shù)字化日變幅數(shù)據(jù),對(duì)該地區(qū)進(jìn)行了時(shí)間跨度較大的地磁場(chǎng)加卸載響應(yīng)比計(jì)算,對(duì)計(jì)算結(jié)果進(jìn)行了異常分析確定,統(tǒng)計(jì)了異常面積與異常強(qiáng)度、異常數(shù)量與地震輻射能量之間的關(guān)系,最后對(duì)同組內(nèi)的異常進(jìn)行了連續(xù)性及地震對(duì)應(yīng)關(guān)系研究。

1　方法原理及Z分量日變幅計(jì)算

1.1地磁加卸載響應(yīng)比P(Z)值計(jì)算原理

太陽(yáng)風(fēng)主要以紫外輻射和粒子流輻射兩種形式影響地球磁場(chǎng),從而使地磁場(chǎng)出現(xiàn)變化磁場(chǎng)。變化磁場(chǎng)主要由太陽(yáng)靜日變化Sq和磁暴場(chǎng)D引起,磁暴場(chǎng)D的構(gòu)成為:D=Dst+Ds+Dp(B),其中Dst為暴時(shí)變化,Ds為暴時(shí)擾日變化,Dp(B)為極區(qū)亞暴。

地球自轉(zhuǎn)使太陽(yáng)的紫外輻射對(duì)地球每天加卸載一次,其地磁效應(yīng)即是Sq;太陽(yáng)暴風(fēng)粒子流輻射的地磁效應(yīng)就是全球同時(shí)發(fā)生的磁暴現(xiàn)象,由統(tǒng)計(jì)可知平均大約5～6天地球便會(huì)被太陽(yáng)高速暴風(fēng)粒子流加卸載一次。另一方面,地震是一種非線性失穩(wěn)現(xiàn)象,孕震區(qū)地下介質(zhì)系統(tǒng)由穩(wěn)態(tài)變?yōu)榉欠€(wěn)態(tài)的過(guò)程中,介質(zhì)的物理性質(zhì)(包括電導(dǎo)率)將會(huì)發(fā)生相應(yīng)的變化。因而可以認(rèn)為,不同的地磁測(cè)點(diǎn),穩(wěn)定地區(qū)和非穩(wěn)定地區(qū)暴時(shí)擾日變化不同;同一地磁測(cè)點(diǎn)非穩(wěn)定時(shí)期暴時(shí)擾日變化與正常時(shí)期的變化又不同。因此可通過(guò)計(jì)算地磁加卸載響應(yīng)比P(Z)值來(lái)得到其異常信息[11]。

地磁垂直分量Z與地下介質(zhì)及其變化較其他地磁分量的關(guān)系更密切,可以取Z分量作為計(jì)算Ds(Z)的加卸載響應(yīng)參量:

P(Z)=Ds(Z)+/Ds(Z)-

(1)

其中:Ds(Z)為Z分量地磁擾動(dòng)場(chǎng)的日變幅,標(biāo)志“+”表示加載,“-”表示卸載。在實(shí)際運(yùn)用中可將地磁場(chǎng)垂直分量日變幅ΔZ+和ΔZ-看作是太陽(yáng)風(fēng)對(duì)地磁場(chǎng)的加載和卸載,此時(shí)計(jì)算公式可簡(jiǎn)化為:

P(Z)=ΔZ+/ΔZ-

(2)

1.2地磁Z分量日變幅計(jì)算

根據(jù)以上定義可知,P(Z)值的計(jì)算需要完整、精確的地磁Z分量日變幅數(shù)據(jù)做基礎(chǔ),目前中國(guó)地磁臺(tái)網(wǎng)提供了數(shù)據(jù)質(zhì)量較好的預(yù)處理地磁數(shù)據(jù),但為了最大限度地降低日變幅計(jì)算結(jié)果的失真度,對(duì)數(shù)據(jù)再做如下處理:

(1)城市軌道交通噪聲影響

我國(guó)主要大城市建設(shè)的直流驅(qū)動(dòng)城市軌道交通在運(yùn)行過(guò)程中對(duì)地漏電造成地下電性結(jié)構(gòu)改變,產(chǎn)生的干擾磁場(chǎng)嚴(yán)重影響了地磁臺(tái)站觀測(cè)數(shù)據(jù)質(zhì)量[12]。其影響方式主要是對(duì)觀測(cè)數(shù)據(jù)造成某一時(shí)段的大動(dòng)態(tài)噪聲污染,進(jìn)而造成地磁Z分量日變幅數(shù)據(jù)失真,因此在計(jì)算日變幅數(shù)據(jù)前需使用一定的降噪方法恢復(fù)其正常日變幅度。

本次研究中受城市軌道交通噪聲影響的臺(tái)站主要是北京地震臺(tái)及靜海地震臺(tái),筆者使用小波變換對(duì)二者觀測(cè)數(shù)據(jù)進(jìn)行了消噪處理。以北京地震臺(tái)2014年1月27日地磁觀測(cè)數(shù)據(jù)為例,由于干擾嚴(yán)重,臺(tái)站當(dāng)天Z分量觀測(cè)數(shù)據(jù)日變幅達(dá)28.8 nT[圖1(a)],筆者選用db2小波對(duì)其進(jìn)行降噪處理后,日變幅降為16.6nT,降幅達(dá)42.4%[圖1(b)],且降噪后的數(shù)據(jù)形態(tài)光滑、平穩(wěn),接近于不受此類(lèi)電磁影響的臺(tái)站數(shù)據(jù)。對(duì)上述兩個(gè)臺(tái)站經(jīng)過(guò)該方法處理后,地磁Z分量觀測(cè)數(shù)據(jù)失真度大大降低,已經(jīng)符合計(jì)算日變幅的要求。

圖1　城市軌道交通電磁干擾及小波消噪結(jié)果Fig.1　Urban railtransit Electromagnetic interference and wanelet de-noising result

(2)磁暴導(dǎo)致日變幅失真影響

地磁Z分量日變幅數(shù)據(jù)平均5～6天會(huì)出現(xiàn)一個(gè)高值,即太陽(yáng)暴風(fēng)粒子流平均約5～6天對(duì)地磁場(chǎng)產(chǎn)生一次明顯的加卸載作用,這是計(jì)算地磁加卸載響應(yīng)比的物質(zhì)基礎(chǔ)。但偶爾發(fā)生的大磁暴會(huì)造成地磁場(chǎng)的猛烈變化,導(dǎo)致正常的地磁場(chǎng)日變幅高值被淹沒(méi)從而使日變幅計(jì)算結(jié)果失真[圖2(a)]。經(jīng)筆者統(tǒng)計(jì),本次研究中共有超過(guò)40臺(tái)次磁暴現(xiàn)象造成響應(yīng)比計(jì)算結(jié)果嚴(yán)重失真。針對(duì)此情況,筆者手工剔除了因磁暴引起的地磁瞬時(shí)巨幅突跳數(shù)據(jù),使日變幅計(jì)算結(jié)果恢復(fù)到正常的太陽(yáng)暴風(fēng)粒子加載水平[圖2(b)]。

圖2　磁暴導(dǎo)致日變幅遠(yuǎn)高于正常暴風(fēng)加載Fig.2　Magnetic storm causes daily ranges much larger than normal

2　研究區(qū)域及地震目錄

本次華北地磁加卸載響應(yīng)比研究區(qū)范圍及所用地磁臺(tái)站分布見(jiàn)圖3,涉及地磁臺(tái)站共計(jì)23個(gè),主要觀測(cè)儀器為磁通門(mén)磁力儀、FHD分量質(zhì)子磁力儀及FHDZ-M15綜合磁力儀。研究過(guò)程優(yōu)先使用觀測(cè)連續(xù)率較好的儀器資料,當(dāng)某臺(tái)儀器出現(xiàn)長(zhǎng)時(shí)間缺數(shù)時(shí),使用該臺(tái)站其他儀器的數(shù)據(jù)進(jìn)行數(shù)據(jù)補(bǔ)缺。臺(tái)站名稱(chēng)、儀器情況見(jiàn)表1。地震目錄選取2008年1月1日—2015年3月31日研究區(qū)內(nèi)發(fā)生的M4及以上地震,震中分布及地震基本參數(shù)見(jiàn)圖3、表2。

3　地磁加卸載響應(yīng)比分析

3.1P(Z)異常判定標(biāo)準(zhǔn)確定及映震統(tǒng)計(jì)

對(duì)上述再處理的地磁場(chǎng)Z分量數(shù)據(jù)進(jìn)行日變幅計(jì)算,得出各臺(tái)站2008年1月—2014年12月的地磁加卸載響應(yīng)比P(Z)值(圖4)。結(jié)果顯示華北地區(qū)各地磁臺(tái)站加卸載響應(yīng)比計(jì)算結(jié)果在整體上有一定的相似性,數(shù)據(jù)曲線起伏及成簇性較為一致,體現(xiàn)了地磁場(chǎng)的整體性及統(tǒng)一性。從計(jì)算數(shù)值分布來(lái)看,華北地區(qū)數(shù)字化地磁P(Z)值多集中在1.0～1.5,部分地區(qū)最高值能達(dá)到5和6左右,顯示出一定的地區(qū)差異。

表1　地磁臺(tái)站及觀測(cè)儀器統(tǒng)計(jì)

表2　M4.0以上地震目錄

目前將地磁加卸載響應(yīng)比作為一種地震前兆手段來(lái)研究地震時(shí),使用較為廣泛的做法是對(duì)P(Z)值設(shè)置一定的閾值指標(biāo)P0(Z),即認(rèn)為當(dāng)P(Z)≥P0(Z)時(shí)地磁加卸載響應(yīng)比出現(xiàn)異常變化,由此來(lái)分析在該閾值指標(biāo)下P(Z)值的變化與后續(xù)地震之間的對(duì)應(yīng)關(guān)系[13-14]。

根據(jù)各省局地磁分析經(jīng)驗(yàn)及筆者反復(fù)研究,認(rèn)為當(dāng)P0(Z)=3.0時(shí)且在一定的規(guī)則前提下,研究區(qū)內(nèi)MS4.0以上地震和P(Z)高值異常之間有一定的對(duì)應(yīng)關(guān)系。具體的規(guī)則為:

(1)取P0(Z)=3.0(圖4中的橫虛線),即認(rèn)定當(dāng)P(Z)≥3.0時(shí)出現(xiàn)地磁加卸載響應(yīng)比異常;

(2)必須至少有2個(gè)或2個(gè)以上臺(tái)站同時(shí)滿(mǎn)足P(Z)≥P0(Z)時(shí)該異常才有效;

(3)各異常的日期間隔小于3個(gè)月的均歸為一組異常,每組異常的有效期為3個(gè)月,即地震發(fā)生在組內(nèi)或組后3個(gè)月均認(rèn)為與該組異常對(duì)應(yīng);

(4)取異常臺(tái)站有效映震半徑為200 km。

根據(jù)上述規(guī)則對(duì)各臺(tái)站加卸載響應(yīng)比異常進(jìn)行統(tǒng)計(jì),結(jié)果顯示,2008年1月—2014年12月華北地區(qū)出現(xiàn)有效異常38次,共計(jì)10組,具體結(jié)果見(jiàn)表3。將此結(jié)果與表1所列地震進(jìn)行對(duì)應(yīng)后認(rèn)為有7組異常和后續(xù)地震存在時(shí)間及空間上的同時(shí)對(duì)應(yīng)關(guān)系,而另外3組異常無(wú)地震對(duì)應(yīng),為虛報(bào),虛報(bào)率為30%。另外,24次地震中有8次在震前未出現(xiàn)時(shí)間及空間同時(shí)滿(mǎn)足規(guī)則的地磁加卸載響應(yīng)比異常,為漏報(bào),漏報(bào)率為33%,其余未漏報(bào)的16次地震震前均出現(xiàn)異常。

3.2異常強(qiáng)度及地震能量分析

本次研究中華北地區(qū)共出現(xiàn)異常38次,異常強(qiáng)度最大的一次出現(xiàn)在2010年12月16日,當(dāng)天的P(Z)峰值出現(xiàn)在天津靜海臺(tái)(地區(qū)),達(dá)到了6.1;而異常范圍最大的一次出現(xiàn)在2011年3月4日,共有19個(gè)臺(tái)站出現(xiàn)了高P(Z)值?？紤]到地磁場(chǎng)是一個(gè)整體性的地球物理場(chǎng),當(dāng)震前出現(xiàn)地磁異常時(shí)其異常地域面積與異常強(qiáng)度在一定程度上應(yīng)存在一定聯(lián)系,為此統(tǒng)計(jì)了38次異常發(fā)生時(shí)異常臺(tái)站數(shù)量(異常區(qū)域面積)與最大異常強(qiáng)度之間的相關(guān)性(圖5),發(fā)現(xiàn)兩者之間基本存在一種正相關(guān)性,即相比于地域面積覆蓋較小的地磁加卸載響應(yīng)比異常,較大地域面積的異常往往會(huì)出現(xiàn)更大的P(Z)峰值。

表3　華北地磁Z分量響應(yīng)比異常及地震對(duì)應(yīng)情況統(tǒng)計(jì)

圖5　38次加卸載響應(yīng)比異常面積與異常強(qiáng)度相關(guān)性Fig.5　Correlativity between area and strength of 38 LURR anomalyies

此外,在上節(jié)中提到的10組加卸載響應(yīng)比異常中有7組存在地震對(duì)應(yīng)關(guān)系,將各組的異常臺(tái)站數(shù)分別統(tǒng)計(jì),并對(duì)各組內(nèi)的對(duì)應(yīng)地震進(jìn)行地震輻射能量計(jì)算[15],最后對(duì)10組異常臺(tái)站數(shù)量與地震能量做了相關(guān)性分析(圖6),結(jié)果發(fā)現(xiàn),組內(nèi)異常臺(tái)站越多,該組可能釋放的地震能量也就越大。例如,除了第1組、第4組及第6組無(wú)地震對(duì)應(yīng)外(圖中能量值記為零),其余各組中地震能量最小的為第10組,該組累積異常臺(tái)站數(shù)為32個(gè),地震能量為0.24×1012J,而能量最大的為第8組,組內(nèi)累積異常臺(tái)站36個(gè),發(fā)生過(guò)一次遼寧燈塔MS5.1地震,地震能量為2.82×1012J。

3.3異常演化及震例分析

利用各臺(tái)站P(Z)值插值得到有震例對(duì)應(yīng)的7組加卸載響應(yīng)比異常平面圖,并標(biāo)注所有對(duì)應(yīng)震例(圖7)。結(jié)果顯示華北地區(qū)地磁加卸載響應(yīng)比異常區(qū)域展布明顯,P(Z)=3.0分界線清晰,部分組中的異常有明顯的連續(xù)變化特征,且和震例有較好的空間或時(shí)間上的對(duì)應(yīng)關(guān)系。下面對(duì)7組異常分別進(jìn)行描述:

第二組[圖7(a)～(f)]:該組異常共計(jì)6次,開(kāi)始于2008年11月27日。該組異常有明顯的位置遷移現(xiàn)象,開(kāi)始異常多集中在華北東部山東及遼南部分區(qū)域,之后逐漸西移至山西及河南一帶,隨后發(fā)生山西原平MS4.3地震,震后異常強(qiáng)度減弱,但仍集中在華北西部。最后在黃海北部發(fā)生MS4.2地震,與之前該區(qū)域附近遼寧大連臺(tái)的異常有對(duì)應(yīng)關(guān)系。

圖6　各組異常中累積異常臺(tái)站數(shù)與地震能量相關(guān)性Fig.6　Correlativity between number of stations and seismic energy in each set of LURR anomalies

第三組[圖7(g)～(l)]:該組同樣出現(xiàn)6次異常,且各次異常在位置分布上有明顯的連續(xù)性,例如第1次異常與第2次、第6次有明顯的位置變化上的聯(lián)系,異常遷移現(xiàn)象明顯,而第3次與第5次也屬于類(lèi)似情況,并且本組5次震例的震中位置與各次異常有明顯的位置對(duì)應(yīng)關(guān)系。

第五組[圖7(m)～(p)]:該組有4次異常,各次異常的中心位置較分散,位置變化的連續(xù)性不很明顯,2011年3月4日最后一次異常的面積最大,之后接連發(fā)生三次MS4以上地震,震中位置也較分散,認(rèn)為這種情況和該組異常分散的特性有關(guān)。

第七組[圖7(q)～(t)]:共有4次異常,各次異常的位置也較為分散,但其中兩次對(duì)應(yīng)震例的震中位置和異常中心有一定對(duì)應(yīng)關(guān)系。

第八組[圖7(u)～(y)]:與第三組特征相似,異常區(qū)域在不同時(shí)間有明顯的連續(xù)性,該組第1次、第2次及第4次異常有明顯的繼承性,而第3次及第5次異常有明顯的聯(lián)系。該組異常只有一次震例對(duì)應(yīng),發(fā)生在第2次及第5次異常范圍內(nèi)部。

第九組[圖7(z)～(dd)]:與第二組較為類(lèi)似,前后各次異常有一定的聯(lián)系,異常區(qū)位置有明顯的遷移現(xiàn)象,異常強(qiáng)度方面有弱-強(qiáng)-弱的變化趨勢(shì),且組內(nèi)唯一一次地震發(fā)生在異常強(qiáng)度最大的第3次異常之后。

第十組[圖7(ee)～(hh)]:共發(fā)生了4次異常,各次異常面積、強(qiáng)度及位置都不盡相同。和第五組相似的是該組異常同樣是在出現(xiàn)一次較大異常之后發(fā)生了一次顯著地震,即河北涿鹿MS4.3地震。

4　討論

本文嘗試將地磁加卸載響應(yīng)比方法應(yīng)用于大范圍的地震分析研究,通過(guò)數(shù)據(jù)處理、異常確定與映震情況分析等得到了如下認(rèn)識(shí):

(1)計(jì)算地磁加卸載響應(yīng)比時(shí)需要質(zhì)量良好的地磁連續(xù)記錄數(shù)據(jù),筆者在研究過(guò)程中發(fā)現(xiàn),持續(xù)缺數(shù)、單點(diǎn)突跳都會(huì)造成響應(yīng)比結(jié)果出現(xiàn)嚴(yán)重偏差,往往會(huì)對(duì)后期的平面異常中心位置造成嚴(yán)重失真,因此對(duì)其進(jìn)行人工校核是十分必要的。其次,需要處理電磁環(huán)境影響造成的數(shù)據(jù)噪聲,本次研究中發(fā)現(xiàn)對(duì)受此類(lèi)影響的臺(tái)站數(shù)據(jù),是否采取消噪處理會(huì)直接影響后期響應(yīng)比計(jì)算結(jié)果的本底值大小,經(jīng)過(guò)相關(guān)處理,計(jì)算結(jié)果會(huì)恢復(fù)到一個(gè)合理區(qū)間范圍。

(2)研究結(jié)果顯示各臺(tái)站地磁加卸載響應(yīng)比結(jié)果在整體上有一定的相似性,但各地區(qū)由于地下電性結(jié)構(gòu)等因素的不同造成了彼此之間的差異,這是該方法用于地震研究分析的物理基礎(chǔ)。本次研究制定了嚴(yán)格的異常確立規(guī)則及地震對(duì)應(yīng)準(zhǔn)則,統(tǒng)計(jì)得出該方法在華北地區(qū)大約有30%的虛報(bào)率,地震漏報(bào)率為33%,且認(rèn)為計(jì)算結(jié)果較為合理。

(3)嘗試分析的38次地磁加卸載響應(yīng)比異常發(fā)生時(shí),各次異常的展布面積與異常強(qiáng)度之間基本存在一種正相關(guān)。另外該10組異常中各組內(nèi)的累計(jì)異常臺(tái)站數(shù)與地震輻射能量之間也基本存在一定的正相關(guān)性,即組內(nèi)異常臺(tái)站越多,該組可能釋放的地震能量也就越大。

(4)最后對(duì)有震例對(duì)應(yīng)的7組異常進(jìn)行了平面展布分析,發(fā)現(xiàn)這7組異常中組內(nèi)各次異常之間有聯(lián)系性的占了較大比例,其主要表現(xiàn)為各次異常在位置關(guān)系上基本存在一種遷移、繼承現(xiàn)象。另外,地震震中位置與各組內(nèi)的各次異常有明顯的聯(lián)系性,大部分地震發(fā)生在曾經(jīng)發(fā)生過(guò)異常的地區(qū),并且震中位置會(huì)隨著異常的位置變化而發(fā)生一定的變化。

References)

[1]徐文耀.地磁學(xué)[M].北京:地震出版社,2003:221-262.

XU Wen-yao.Geomagnetism[M].Beijing:Seismological Press,2003:221-262.(in Chinese)

[2]徐文耀.地球電磁現(xiàn)象物理學(xué)[M].北京:中國(guó)科學(xué)技術(shù)大學(xué)出版社,2009:253-254.

XU Wen-yao.Physics of Electromagnetic Phenomena of the Earth[M].Beijing:Press of University of Science and Technology of China,2009:253-254.(in Chinese)

[3]丁鑒海,盧振業(yè),黃雪香.地震地磁學(xué)[M].北京:地震出版社,1994.

DING Jian-hai,LU Zhen-ye,HUANG Xue-xiang.Seismogeomagnetism[M].Beijing:Seismological Press,1994.(in Chinese)

[4]郭增建,郭安寧,張煒超,等.甘肅岷縣漳縣6.6級(jí)地震發(fā)生日期的觸發(fā)因素分析及對(duì)短臨預(yù)測(cè)方法的思考[J].地震工程學(xué)報(bào),2013,35(3):413-418.

GUO Zeng-jian,GUO An-ning,ZHANG Wei-chao,et al.Analysis of Trigger Factors of the Minxian—Zhangxian MS6.6 Earthquake Occurrence Date and Short-impending Prediction Methods[J].China Earthquake Engineer Journal,2013,35(3):413-418.(in Chinese)

[5]許康生,辛長(zhǎng)江,李英.汶川地震前后地磁Z分量的譜質(zhì)心變化[J].地震工程學(xué)報(bào),2014,36(2):393-397.

XU Kang-sheng,XIN Chang-jiang,LI Ying,et al.The Spectrum Centroid Variation in Geomagnetic Z Component Data before and after the Wenchuan Earthquake[J].China Earthquake Engineer Journal,2014,36(2):393-397.(in Chinese)

[6]馮志生,李琪,李鴻宇,等.地磁低點(diǎn)位移線兩側(cè)異常變化的反相位現(xiàn)象及其解釋[J].中國(guó)地震,2009,25(2):206-213.

FENG Zhi-sheng,LI Qi,LI Hong-yu,et al.Anti-phase Phenomena on Two Sides of Geomagnetic Low-point Displacement Line and Its Interpretation[J].Earthquake Research in China,2009,25(2):206-213.(in Chinese)

[7]許儀西,方韡,陳新澤,等.2008年7月5日長(zhǎng)泰ML4.7地震前流動(dòng)地磁觀測(cè)結(jié)果的分析[J].地震地磁觀測(cè)與研究,2010,31(1):39-42.

XU Yi-xi,FANG Wei,CHEN Xin-ze,et al.Analysis on the Results in Mobile Geomagnetic Survey before Changtai ML4.7 Earthquake in Fujian Province on July 5,2008[J].Seismological and Geomagnetic Observation and Research,2010,31(1):39-42.(in Chinese)

[8]馮志生,梅衛(wèi)萍,張秀霞,等.中強(qiáng)震前地磁諧波振幅比的趨勢(shì)性變化特征初步研究[J].西北地震學(xué)報(bào),2004,26(1):50-56.

FENG Zhi-sheng,MEI Wei-ping,ZHANG Xiu-xia,et al.Preliminary Study on the Characteristics of Long and Medium Term Variation for the Amplitude Radio of Geomagnetic Humorous Wave before Moderate-strong Earthquakes[J].Northwestern Seismological Journal,2004,26(1):50-56.(in Chinese)

[9]尹祥礎(chǔ),尹燦.非線形系統(tǒng)的失穩(wěn)前兆與地震預(yù)測(cè)[J].中國(guó)科學(xué):B輯,1991(5):512-518.

YIN Xiang-chu,YI Can.Unstable Precursor of Non-linear System and Earthquake Forecasting[J].China Science:Series B,1991,(5):512-518.(in Chinese)

[10]宋治平,尹祥礎(chǔ),陳學(xué)中.加卸載響應(yīng)比的時(shí)空演變特征及其對(duì)地震三要素的預(yù)測(cè)意義[J].地震學(xué)報(bào),1996,18(2):179-186.

SONG Zhi-ping,YI Xiang-chu,CHEN Xue-zhong,et al.Time and Spatial Variation Characteristic of Loading-unloading Response Ratio and Three Essential Factors of Earthquake Forecast Significance[J].Acta Seismologica Sinica,1996,18(2):179-186.(in Chinese)

[11]曾小蘋(píng),續(xù)春榮,趙明,等.地球磁場(chǎng)對(duì)太陽(yáng)風(fēng)的加卸載響應(yīng)與地震[J].地震地磁觀測(cè)與研究,1996,17(1):49-53.

ZENG Xiao-ping,XU Chun-rong,ZHAO Ming,et al.The Load-unload Response of Geomagnetic Field to Solar Wind and Earthquake Prediction[J].Seismological and Geomagnetic Observation and Research,1996,17(1):49-53.(in Chinese)

[12]謝凡,滕云田,胡星星,等.地磁臺(tái)站的城市軌道交通干擾的小波抑制方法研究——以天津軌道交通干擾為例[J].地球物理學(xué)報(bào),2011,54(10):2698-2707.

XIE Fan,TENG Yun-tian,HU Xing-xing,et al.Suppression of Magnetic Perturbation Caused by Urban Railway to Geomagnetic Observations Using the Wavelet Method——A Case Study of Tianjin Subway[J].Chinese Journal of Geophysics,2011,54(10):2698-2707.(in Chinese)

[13]張建國(guó),閆俊崗,王靜,等.地磁加卸載響應(yīng)比方法在地震預(yù)報(bào)中的應(yīng)用[J].大地測(cè)量與地球動(dòng)力學(xué),2008,28(2):45-50.

ZHANG Jian-guo,YAN Jun-gang,WANG Jing,et al.Application of Geomagnetic Loading-unloading Response Ratio Method to Earthquake Prediction[J].Journal of Geodesy and Geodynamics,2008,28(2):45-50.(in Chinese)

[14]李偉,龔耀,趙文舟,等.地磁加卸載響應(yīng)比方法在上海及其鄰區(qū)地震研究中的應(yīng)用[J].地震,2014,34(1):125-133.

LI Wei,GONG Yao,ZHAO Wen-zhou,et al.Application of Geomagnetic LURR Method in Seismic Studies in Shanghai and Its Adjacent Areas[J].Earthquake,2014,34(1):125-133.(in Chinese)

[15]中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司.地震參數(shù)——數(shù)字地震學(xué)在地震預(yù)報(bào)中的應(yīng)用[M].北京:地震出版社,2003:73-74.

Department of Monitoring and Forecasting,China Earthquake Administration.Seismic Parameters:Application of Digital Seismology in Earthquake Prediction[M].Beijing:Seismological Press,2003:73-74.(in Chinese)

Application of Geomagnetic Load/Unload Response Ratio Method in North China

JIA Li-feng1,ZHANG Gou-ling1,QIAO Zi-yun1,LIANG Hong-jie1,ZHANG Bo2

(1.Earthquake Administration of Hebei Province,Shijiazhuang 050021,Hebei,China;2.Baishe Seismic Station,Earthquake Administration of Hebei Province,Xingtai 056001,Hebei,China)

The load/unload response ratio (LURR)method has been applied for many years in seismicity forecasting and has achieved good results.Over time,it has become the main analysis method in the field of seismomagnetism and is often used in short-term earthquake analysis or prediction.In this study,we used LURR to extract anomalies from digital geomagnetic data recorded at 23 geomagnetic stations during 2008 to 2014 in North China and investigated the correspondence between earthquake occurrence and the anomalies.First,we carefully processed the geomagnetic data to eliminate any distortions caused by missing data and magnetic storms.For example,we examined a few instances of electromagnetic interference around the stations and chose the wavelet method to process the resulting noise.The research region,North China,has numerous evenly distributed geomagnetic stations,and as our research object,we selected 24 earthquakes stronger than MS4 that occurred in this region.The results show that there were ten groups of LURR anomalies in the seven-year study period in North China,of which seven were thought to be related to earthquakes and the other three were not.In addition,16 of the selected 24 earthquakes had precursory anomalies.We found the maximum anomaly intensity to be proportional to the anomaly area in each group of anomalies.Therefore,the larger the anomaly area,the greater the maximum intensity P(Z)is likely to be.Moreover,we found that seismic energy is associated with the number of stations recording anomalies in one group.Therefore,the greater the number of stations recording anomalies in one group,the larger is the earthquake radiation energy.Lastly,in the groups of anomalies related to earthquakes,we found a relationship between each anomaly,mainly reflected in the position continuity and inheritance.For example,an anomaly was more likely to occur near the position of the preceding anomaly.We also found a relationship between the epicenter of the earthquake and the anomaly location.

North China; geomagnetic load-unload response ratio; precursory anomaly; earthquake prediction

2015-07-30

國(guó)家自然科學(xué)基金項(xiàng)目(41274079);地震行業(yè)科研專(zhuān)項(xiàng)經(jīng)費(fèi)項(xiàng)目(201508013)

賈立峰(1981-)，男，河北石家莊人，工程師，主要從事地震電磁分析工作。E-mail:729102262@qq.com。

P315.73

A

1000-0844(2016)04-0588-10

10.3969/j.issn.1000-0844.2016.04.0588