葉向頂，顧申宜，張 慧

(海南省地震局，海南 海口 570203）

0 前言

近年來，利用地磁垂直分量日變化出現(xiàn)的極小值的時間以及其在空間分布上產(chǎn)生的異常變化，進(jìn)行短期的或臨震的地震預(yù)報倍受關(guān)注。經(jīng)過多年的震例和經(jīng)驗積累，分析和利用低點異常或雙低點異常預(yù)測地震亦取得了一定的進(jìn)展[1～4]。盡管其機(jī)理成因等還需要進(jìn)一步研究，但通過大量實例的統(tǒng)計分析證明，其在地震預(yù)測上有一定的實效性。瓊中地磁臺(地理坐標(biāo)為北緯19°51′，東經(jīng)109°01′，海拔高度為230 m)地處中國的最南端的海南省，是離赤道最近、緯度最低的地磁臺，其地磁緯度為7.8°，在全國地磁測點中具有地緣性的優(yōu)勢。為此，詳細(xì)地分析其垂直分量的低點時間變化特征，有助于我們利用低點異常現(xiàn)象做為震磁前兆異常，提高地震預(yù)測的成功率。

1 數(shù)據(jù)與資料處理

地磁日變化中，Z分量日變極小值出現(xiàn)的時間為低點時間，其在空間的異常特征稱之為 “低點位移”現(xiàn)象，是少見的地磁異?，F(xiàn)象。在地球北半球，地磁Z分量的日變化正常形態(tài)是：總體上呈現(xiàn)周期性變化；每日中午出現(xiàn)極小值，上午、下午各出現(xiàn)一個極大值。

本文采用了2002至2012瓊中地磁臺GM3磁通門磁力儀Z分量數(shù)據(jù)，首先統(tǒng)計了每一年中各個月的低點時間的平均值，然后再對各年的某個月數(shù)據(jù)作平均，最后得到一年中12個月低點時間的分布特征。同時分別對各個季節(jié)低點時間的分布和標(biāo)準(zhǔn)偏差進(jìn)行分析，確定其平均值分布區(qū)間。

2 地磁Z分量低點時間的特征分析

2.1 Z值低點時間的平均分布狀態(tài)

地磁場“日變化”是地球變化磁場的重要組成部分，是地球物理學(xué)地磁觀測中的基本規(guī)律之一。它起源于太陽輻射在電離層中產(chǎn)生的電流體系的影響，包含著空間和地球內(nèi)部各種電磁過程豐富的信息。因此其變化也是比較復(fù)雜的，但單臺低點時間的平均分布狀態(tài)又是很有規(guī)律的，它集中分布在固定的時間段內(nèi)，在極大值兩側(cè)出現(xiàn)頻次迅速平緩，符合正態(tài)分布：樣本(低點時間)x落入任意區(qū)間(a，b)的概率，記作

落在(μ-σ， μ+σ)的概率是68.3% 左右(如圖1)。

圖1 地磁Z值低點時間正態(tài)分布統(tǒng)計圖 （UTC）Fig.1 The normal distribution chart of the low point time of geomagnetic Z-component

2.2 Z值低點時間月均值和年均值

為了能夠準(zhǔn)確地統(tǒng)計出瓊中臺Z分量低點時間，我們設(shè)計了一些小程序，統(tǒng)計了2002年至2012年11年的數(shù)據(jù)，繪制地磁各個季低點時間分布狀況圖（圖1），并算出一年中各個月的月均值和年均值，如表1。

由表1可知，瓊中臺地磁Z分量全部資料年平均低點時間為地方時12.87時即12點52分，由式（1）可計算出其主要的低點時間分布區(qū)域在地方時12點10分至13點34分（UTC 04：10～05：34）；各個月平均低點時間分布如圖2，一年中各個月的低點時間近似于一個字母“W”形狀，均值線恰好地把一年切成上下六個月，其中首尾均處在最高值，顯示一年中有半年的時間低點出現(xiàn)時間是滯后的，分別為1、2、3月和10、11、12月，而另半年則為4月至9月份，即Z值低點時間的出現(xiàn)具有半年為周期的分布特征?！癢”形的兩個低值時間分別是5月份和九月份，最早的低點時間（趨前）出現(xiàn)在5月份，最晚出現(xiàn)則是11月份，正好也相隔半年時間，也是具有半年周期特征。5月份的平均低點時間比年平均低點時間提前了差不多一個小時（地方時12點02分），11月份最晚則延遲至地方時13點45分，年平均低點時間也延時了幾乎1個小時，即最早和最晚之間兩者相差幾乎有兩個小時。然而，應(yīng)該指出的是，這僅僅是平均低點時間分布狀況，且它也只是與概率上的68.3%，少量資料即使連續(xù)的靜日低點時間的出現(xiàn)也可能有更大的差異，可能更為分散，相差有可能達(dá)好幾個小時。因此，我們在統(tǒng)計低點時間，判別低點位移異常時，應(yīng)該充分考慮不同月份的低點出現(xiàn)情況，同時也更要考慮個別情況的特殊性。因此，區(qū)域性的低點位移是必須要考慮的，以便更準(zhǔn)確提出異常地震前兆信息。

表1 瓊中臺歷年地磁Z分量低點時間月均值 （北京時）Table 1 The monthly average of the low point time of geomagnetic Z-component over the years in Qiongzhong station

圖2 瓊中臺地磁Z分量低點時間月均值分布圖Fig.2 The monthly average of the low point time of the geomagnetic Z-component of Qiongzhong station

2.3 低點時間與低點位移異常的季節(jié)性分布

由圖2也可看出，低點時間的出現(xiàn)是具有季節(jié)性的：在1、2月和11、12月在均值之上，按地磁季節(jié)性來分，正好處在冬季，即冬季低點出現(xiàn)的時間相對延遲（滯后）；春秋季（3、4、9、10月）則處在均值點，最接近年平均低點時間；夏季5～8月則處在均值線之下，即低點時間趨前。各個季節(jié)低點時間出現(xiàn)的頻次亦有所不同。表2資料表明，瓊中臺低點平均時間分布的離差σ大概在半小時內(nèi)，相對比較集中。但μ值相差較大，以夏季(12:19)和冬季(13:31)最為明顯，差值達(dá)一個多小時。對應(yīng)的表3白家疃地磁臺[1]，其季節(jié)性的分布則明顯不同，白家疃臺整個數(shù)據(jù)的平均低點時間分布離差σ較大，近1 h，顯然其分布時間上較為離散；與瓊中臺不同的是，μ值的差別最大的是春秋季和夏季，也達(dá)到1 h，夏季與冬季之差只在半小時以內(nèi)。相對瓊中臺，其各個季節(jié)之間分布的差別很小。日變低點時間顯然主要是受經(jīng)度影響，經(jīng)度每變化15°，則低點時間（北京時）相差1 h。出現(xiàn)這種季節(jié)上明顯不同的原因可能是緯度的差異和地理位置不同所造成。

表2 瓊中臺正態(tài)分布的均值和標(biāo)準(zhǔn)偏差 （北京時）Table 2 The average and standard deviation of the normal distribution of Qiongzhong station(cst)

表3 白家疃臺正態(tài)分布的均值和標(biāo)準(zhǔn)偏差 （北京時）Table 3 The average and standard deviation of the normal distribution of Baijiatong station(cst)

以上所示是為Z值低點時間的平均分布狀況，以此為基礎(chǔ)，我們分別又統(tǒng)計了全年的低點位移異常情況。表4是在充分考慮到不同月份、季節(jié)之后，統(tǒng)計出的2002～2012年11年間的低點位移異常。從表4和圖3上看，低點位移異常的出現(xiàn)以2月和11月最多，而9月份異常的出現(xiàn)為零；即，冬季低點異常頻度最大、春秋季最小，11年間的數(shù)據(jù)季節(jié)性變化十分明顯，其中春秋季與夏季兩個季節(jié)加起來還沒有冬季異常頻度高，充分顯示了低點異常在冬季高度集中，占全年的60%。同時我們也發(fā)現(xiàn)低點位移異常主要出現(xiàn)在低點時間的滯后區(qū)，即1、2、3月和10、11、12月，占全年的74%。因此，冬季出現(xiàn)的低點位移異常，虛報的概率相對較高。

3 瓊中臺地磁低點位移異常震例

海南島陸及鄰區(qū)中強(qiáng)地震較少，據(jù)歷史記載最大地震為1605年7月13日海南省瓊山市71/2級地震?，F(xiàn)就近十幾年來其所發(fā)生的少量的中強(qiáng)地震中，就其地磁低點位移狀況作一些簡單的探討。

表4 瓊中臺低點異常季節(jié)分布Fig 4 The seasonal distribution of abnormal low point of Qiongzhong station

圖3 瓊中地磁低點位移異常季節(jié)分布曲線Fig.3 The seasonal distribution of geomagnetic abnormal low point displacement of Qiongzhong station

表5 海南及其鄰近地區(qū)地震目錄表Table 5 The earthquake catalog table of Hainan and its adjacent areas

海南及其鄰區(qū)十幾年間所發(fā)生的MS4～5級的地震也有好幾次，對本島近海地震如1992～1995年東方近海震群 （最大震級MS4.8）、1999年萬寧近海MS4.8和2004年廣東陽江MS4.9級地震，其地磁低點無異常反應(yīng)，即漏報。然而對最近廣東東源縣發(fā)生的三次地震，震前瓊中臺地磁低點似乎有異常與之對應(yīng)，且低點位移異常的出現(xiàn)大多也是連續(xù)性出現(xiàn)。對此，我們有所疑惑：地磁低點位移現(xiàn)象是否存在著“盲點”或“敏感”區(qū)域？因資料有限，這里也只是針對瓊中臺資料簡要提點，以資討論。

（1）海南島陸及鄰區(qū)地震在以上11年間6級以上的中強(qiáng)地震僅有兩次，即1994年12月31日和1995年1月10日北部灣6.1、6.2級雙震，震前1994年11月18、20、25日(圖4（a）、（b）、（c）)連續(xù)三天均有低點位移異常出現(xiàn)，低點出現(xiàn)時間均在UTC時間08點左右（北京時16點）。

圖4 北部灣6.1、6.2級地震瓊中低點位移異常時均值圖Fig.4 The average curve of geomagnetic abnormal low point displacement of Qiongzhong station in M6.1、M6.2 earthquake in the North Bay

（2）見圖5廣東河源三次地震瓊中地磁低點位移異常情況：2012年2月16日MS4.8級地震，震前瓊中臺2011年11月22和2012年1月18日均出現(xiàn)低點位移異常（圖5a、5b）；2012年8月31日MS4.2級，震前8月10日瓊中地磁也出現(xiàn)了低點位移異常（圖5c）；2013年2月22日MS4.8級，1月25日、2月9日瓊中地磁連續(xù)出現(xiàn)低點位移異常（圖5d、5e）。

從以上的震例中我們也看到，其大部分地震多發(fā)生在冬季。而上面統(tǒng)計數(shù)據(jù)顯示，冬季地磁低點位移異常最為頻繁和集中、虛報率最高。因此，這些震例中所出現(xiàn)的低點位移異常也極有可能是氣候、季節(jié)原因所產(chǎn)生的地磁異常而非地震地磁異常。因本區(qū)域地震震例相當(dāng)少，這里只能以之作為地磁低點位移現(xiàn)象來探討，是否為地震地磁前兆異常還有待進(jìn)一步研究和更多的震例來說明。

4 結(jié)語

圖5 海南及其鄰區(qū)MS4.0級以上地震低點位移異常圖Fig.5 The abnormal low point displacement of earthquakes with MS≥4 in Hainan and its adjacent areas

4.1 主要認(rèn)識

（1）瓊中臺地磁Z分量平均低點時間出現(xiàn)在北京時12點52分，其平均分布狀態(tài)符合正態(tài)分布特征，主要分布區(qū)域為北京時12:10～13:34。

（2）瓊中臺地磁Z分量平均低點時間具有明顯的季節(jié)性：一般春秋季出現(xiàn)時間與年均值點最接近，夏季則處在均值點之下，即低點時間相對提前，冬季平均低點時間處在年均值點之上，即低點時間相對滯后。

（3）相對表3白家疃等臺站，瓊中臺地磁Z分量平均低點時間分布特征較為明顯，區(qū)間分布基本在一個小時內(nèi)，顯然更為集中。

（4）瓊中臺地磁Z分量低點位移異常的出現(xiàn)也具有季節(jié)性，冬季出現(xiàn)的頻率最高，占全年的60%，異常的出現(xiàn)主要在平均低點時間的滯后區(qū)，因此，虛報率也是最高的。

4.2 討論

地磁低點位移現(xiàn)象是否存在“盲點”或“敏感”區(qū)域：瓊中臺地磁Z分量低點位移異常與本區(qū)域中強(qiáng)地震（MS6.0以上）有一定的對應(yīng)性，但因震例較少，只能作為參考；而對本島近海MS4.5級左右的地震無異常反應(yīng)和對鄰區(qū)MS4.2級以上的地震如廣東河源地震似乎有所對應(yīng)。對此，我們有所疑惑：這是否存在所謂的“盲點”或“敏感”現(xiàn)象？因本區(qū)域地震震例相對較少，這里也只是針對瓊中臺單臺資料簡要提點，是否為地震地磁前兆異常還有待進(jìn)一步研究和更多的震例來說明。

[1]丁鑒海，黃雪香.地磁垂直分量日變“低點位移”預(yù)報地震方法實用化研究[C]//國家地震局科技監(jiān)測司編.地震預(yù)報方法實用化研究文集地磁地電專輯.北京：學(xué)術(shù)書刊出版社，1990.

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

[3]陳紹明，謝美娟.地磁垂直日變低點位移預(yù)報方法研究[C]//國家地震局科技監(jiān)測司編.地震預(yù)報方法實用化研究文集地磁地電專輯.北京：學(xué)術(shù)書刊出版社，1990.

[4]李明光，梁秀蘭.利用地磁Z分量日變低點預(yù)報地震的新探索[C]//國家地震局科技監(jiān)測司編.地震預(yù)報方法實用化研究文集地磁地電專輯.北京：學(xué)術(shù)書刊出版社，1990.