地磁加卸載響應比方法最佳閾值研究1

賈立峰 喬子云 張國苓 丁志華 張建國

?

地磁加卸載響應比方法最佳閾值研究1

賈立峰1）喬子云1）張國苓1）丁志華1）張建國2）

1）河北省地震局，石家莊 050021?2）河北省地震局邯鄲中心臺，邯鄲 056001

地磁加卸載響應比方法（LURR）是從測震學中移植過來的、用于地震中期分析預測的一種地磁學方法。筆者應用紅山地震臺地磁數字化觀測數據進行了地磁加卸載響應比的計算分析，經過統計得出在閾值選擇為3.0時對地震進行預測的效果比較好，其與地震的對應概率為42%。同時，筆者還認為可能受限于方法本身及外界因素，地磁加卸載響應比方法在進行地震預測時，其虛報率及漏報率較高。

地磁 加卸載響應比 閾值 地震預報

引言

地磁場是地球固有的地球物理場，能夠深刻反應地下物質信息的變化（徐文耀，2003；2009），同時地磁場也是一種用來研究地震、預測地震的前兆手段（丁鑒海等，1994）。地磁加卸載響應比方法（LURR）便是其中的方法之一，目前這一方法在全國地震系統應用較為廣泛，已經成為一種用來進行中短期地震分析預測的地磁手段。這一方法的提出源于測震學中的加卸載響應比理論，即在非線性系統失穩的過程中伴隨著外部應力加載與卸載的不平衡作用過程（尹祥礎等，1991）。隨后，加卸載響應比理論被應用于地震活動性預測當中并取得了一定的效果，同時，這一方法在地形變、地下水位、重力等地震前兆資料的分析中也得到了較好的驗證（宋治平等，1996）。例如：曾小蘋等（1996）將這一方法運用于地磁分析中，將磁暴發生過程中地磁場的變化視為磁場對太陽風的加卸載響應，分析了磁暴場暴時擾日變化()的響應比值()隨時間的變化，發現()的高值與地磁臺站周邊一定范圍內S≥5.5級的中強主震有較好的對應結果，隨后，作者使用這一方法對我國華北和西南兩個地震活動帶進行了分析檢驗，同樣發現在地震發生前的不同時間段內()值有出現高值的現象發生；馮志生等（2000）使用這一方法研究了江蘇地磁臺網1983—1997年間的加卸載響應比異常，以及這一期間臺同內和網緣附近S4.6級以上地震之間的關系，并對其內檢預報效能進行了評價，同時提出了適用于江蘇地磁臺分析預測的標志體系；李偉等（2013）的研究發現，同一臺站分別使用57型模擬磁力儀和GM4型數字化磁力儀分別觀測得到的()值存在不同程度的差異，在()>2.5的情況下，兩者的差異比重范圍為-19.51%—0.39%，作者利用()>2.5的頻次分析地震時，由GM4型數字化磁力儀觀測結果計算得到的()值在反映地震前兆異常的能力上較強。

自2002年起，河北省紅山地震臺在繼續進行模擬地磁觀測的基礎上，增加了一套FHD-1型分量質子磁力儀，數據產品為地磁分鐘值相對記錄數據。在“十五”時期，根據中國地震前兆數字化觀測網絡的要求，該分量質子磁力儀升級為FHD-2B型并增加了網絡傳輸協議，同時也提高了儀器的觀測精度。本文對紅山臺周邊一定范圍內的地震數據進行了分析研究，計算并提取了適用于該地區()的閾值，可為紅山臺所在地區震情的判定提供標志依據。

1 紅山臺背景

1966年3月在邢臺地區發生了S7.2級地震，同年4月由中國科學院地球物理研究所在邢臺地區負責籌建紅山地震臺，并于當年7月開始投入使用，之后，在1971年10月移交給河北省地震局負責。

紅山臺地處華北平原南部邢臺市山口鎮，其地理位置如圖1所示。紅山臺位于太行山山前斷裂帶分支斷裂即隆堯束鹿斷裂的南段，海拔高度37m，基巖裸露，周圍是大面積沖積平原，地表出露地層為震旦紀石英砂巖，層位較低，基巖出露面積小。同時，紅山臺遠離城鎮和村莊，干擾源較少，是比較理想的地震及地磁觀測場所。

紅山臺開展地磁觀測較早，自建臺以來已積累了30多年的豐富地磁觀測資料，數據質量穩定可靠。同時作為全國地磁觀測的基準臺，其地磁觀測數據還承擔著對外交換的任務。紅山臺數字化地磁觀測從2002年1月開始，現有的絕對觀測儀器為G856型質子核旋儀及Mingeo-DI型儀，相對記錄儀為GM4型磁通門磁力儀及FHD-2B型分量質子磁力儀。

2 原理與方法

太陽風主要以紫外輻射和粒子流輻射兩種輻射形式影響地球磁場，從而使地磁場出現變化磁場。變化磁場主要有兩種：一是地磁場變化即太陽靜日變化q；二是磁暴場。其中，磁暴場的構成為：=st+s+p()，式中st為暴時變化，s為暴時擾日變化，p()為極區亞暴。

按照曾小蘋等（1996）的地磁加卸載響應比理論，地球自轉使太陽的紫外輻射對地球每天加卸載一次，其地磁效應即是q；太陽暴風粒子流輻射的地磁效應就是全球同時發生的磁暴現象，由統計可知，平均大約5—6天地球便會被太陽高速暴風粒子流加卸載一次。另外，考慮到地震是一種非線性失穩現象，在孕震區地下介質系統由穩態變為非穩態的過程中，介質的物理性質，其中包括電導率將會發生相應的反應。因而可以認為，不同的地磁測點，穩定地區和非穩定地區暴時擾日變化不同；同一地磁測點非穩定時期暴時擾日變化與正常時期的變化又不同。因此，可通過計算地磁加卸載響應比()值來得到其異常信息。

地磁垂直分量與地下介質及其變化較其它地磁分量的關系更密切，因此可以取分量作為計算s()的加卸載響應參量：

()=s()+/s()-（1）

式中，s()為分量地磁擾動場的日變幅；標志“+”表示加載；標志“-”表示卸載。

在實際運用中可將地磁場垂直分量日變幅+和-看作是太陽風對地磁場的加載和卸載，此時計算公式可簡化為：

()=Δ+/Δ-（2）

3 計算結果及分析

目前，紅山臺所使用的地磁相對記錄儀器有兩種：一種為FHD型分量質子磁力儀；另一種為GM4型磁通門磁力儀。本文中采用的數據為FHD型分量質子磁力儀所產出的分鐘值相對記錄數據。該儀器總強度觀測精度為0.3nT，水平分量觀測精度為0.6nT，磁偏角觀測精度為0.1'。

根據紅山臺的地理位置及周邊歷史地震記錄，筆者以距離臺站300km范圍為限，選取了自紅山臺2002年1月數字化地磁觀測以來，臺站周邊發生的S4.0級以上地震作為研究對象，地震目錄及參數如表1所示，地震空間分布如圖1所示。

表1 紅山臺周邊300km以內MS4.0級以上地震統計

紅山臺2012年1月至2013年12月的地磁加卸載響應比計算結果如圖2所示，并在計算前進行了細致校對，剔除了因缺數造成的錯誤日變幅數據。

圖2中出現的高值突跳即為地磁加卸載響應比異常，而且該類異常通常具有成叢性。圖中所標注的地震為臺站周邊300km以內發生的所有S4級以上地震，從圖中可以看出，這些地震的發生和該地磁加卸載響應比的高值分布有一定的對應關系。

同時，對加卸載響應比計算結果進行閾值的選取是確定異常的關鍵。為了確定響應比是否為異常的映震效果，筆者嘗試規定了以下準則：①確定閾值后，高于閾值的值定為異常，異常間隔小于3個月定為一組異常，每組異常有效期為12個月，有效期開始時間以最后高值出現的日值為準；②地震可以在異常結束后發生，也可以發生在異常期間內（震前或震后異常）；③若地震發生在若干組異常的共同有效期內，則判定各組異常均為有效異常，但對應的地震次數只記1次。

按照上述準則，筆者對紅山臺的()計算結果嘗試進行了不同閾值的選擇和統計，對異常閾值的起始值劃分為5個級別即從閾值2.6劃分至閾值3.4，將高于閾值的加卸載高值定為異常，同時根據規定的準則統計各次異常分組，其結果如表2所示。

表2 紅山臺地磁加卸載響應比映震效能檢驗統計

例如從表2可以看出，當()的閾值選擇為3.0時，共統計出異常12組，對應地震次數為5次，對應概率為42%；其對應的地震分別為2006年4月9日河南范縣S4.1級地震、2006年7月4日河北文安S5.1級地震、2008年3月10日河南封丘S4.4級地震、2010年4月4日山西大同S4.2級地震、2010年6月5日山西陽曲S4.6級地震（見表3）。同時，表3還給出了其他閾值異常與地震的對應情況。

表3 紅山臺地磁加卸載響應比（P(Z)=3.0）對應地震統計表

續表

事件序號異常組號異常日期（年-月-日）異常值對應地震發震間隔（天）震中距（km） 1392010-12-163.84虛報 2011-01-133.17 2011-03-043.83 14102012-10-113.31虛報 15112013-10-163.67虛報 2013-12-233.08 2014-01-113.09 16122014-08-304.09虛報

另外，從表3還可以看出，上述5次對應發生的地震與響應比異常的最短間隔為30天，最長間隔為331天。所以從統計結果看，閾值選擇為3.0時，對地震的預測效果較其他閾值更為理想。盡管閾值3.0低于閾值3.2和閾值3.4對應的地震發生概率，但閾值3.0時的地震虛報率和漏報率較低，顯示出了明顯的統計優勢。

4 討論

筆者將紅山臺數字化地磁數據進行了加卸載響應比計算及分析，確定了可以利用的合理閾值起始值。同時還得出了以下幾點初步認識：

（1）地磁加卸載響應比方法在統計學意義上對地震有一定的預測性。紅山臺的計算結果顯示，選擇合適的異常閾值后可以在未來1年的時間尺度內，對S4.0級及以上地震作出40%以上的預測，此時的虛報率為50%。

（2）加卸載響應比方法有時會出現連續性的漏報及虛報。從表3的統計結果可以看出，紅山臺對最早的3次震例均為漏報，造成漏報率高的直接原因可能是由于缺乏2001年的響應比數據，因此，滿足一定時間長度的地磁數據樣本是對該方法進行計算、驗證及研究的前提條件。

另外，該方法虛報率較高也是一個無法回避的問題。筆者在對紅山臺進行過的所有閾值劃定結果中，其虛報率始終不低于50%。造成這種結果的原因除了該方法本身外，對地震目錄選擇限制過多也可能是另一種原因，假如對300km以內及S4級以上兩個地震目錄的限制條件有所改變，其響應比計算結果可能會有所不同。這一點筆者將在以后的研究中對其進行進一步的計算驗證。

總之，地磁加卸載響應比方法是從測震學中移植借鑒的一種方法，其使用的是能較好地反應地下介質導電性質的響應比()作為響應函數，并以它在正常情況下的時間變化和空間分布作為背景，研究震前由于電導率變化引起的地磁響應比異常來預測地震。在我國目前的地震科研及地震預測預報體系中，對預測方法的研究及建立一定的映震指標體系是一項重要的工作，以此為基礎可以對地磁前兆數據進行合理的異常量化及判定，這對地震分析預報工作有著積極的意義。

丁鑒海，盧振業，黃雪香，1994．地震地磁學．北京：地震出版社．

馮志生，林蕓芳，王建宇等，2000．江蘇地磁加卸載響應比的異常標志體系．地震，20（2）：61—68．

李偉，龔耀，2013．上海地磁加卸載響應比資料分析與地震研究．華南地震，33（4）：31—37．

宋治平，尹祥礎，陳學中，1996．加卸載響應比的時空演變特征及其對地震三要素的預測意義．地震學報，18（2）：179—186．

徐文耀，2003．地磁學．北京：地震出版社．

徐文耀，2009．地球電磁現象物理學．北京：中國科學技術大學出版社．

尹祥礎，尹燦，1991．非線形系統的失穩前兆與地震預測．中國科學（B輯），5（5）：512—518．

曾小蘋，續春榮，趙明等，1996．地球磁場對太陽風的加載響應與地震．地震地磁觀測與研究，17（1）：49—53．

The Study of Optimal Threshold of Geomagnetic Load-Unloading Response Ratio

Jia Lifeng1), Qiao Ziyun1), Zhang Gouling1), Ding Zhihua1)and Zhang Jianguo2)

1) Earthquake Administration of Hebei Province, Shijiazhuang 050021, China?2) Handan Seismic Station, Handan 066001, China

As a middle-term earthquake prediction method, the geomagnetic load-unload response ratio (LURR) was transplanted from the seismological method. Based on digital observation data, we calculated the LURR of Hongshan station. Our results show that when the threshold is 3.0, the accurate rate of the earthquake prediction probability is about 42%. Meantime, find that the false and missing ratio is relatively high when we applied the LURR to predict earthquake.

Geomagnetic; Load-unload response ratio; Threshold; Earthquake prediction

國家自然科學基金（41274079）；地震行業科研專項（201508009-09）

2014-04-09

賈立峰，男，生于1981年。碩士。主要從事地震地磁分析工作。E-mail：729102262@qq.com