泉州臺地電場數據分析

倪曉寅

（福建省地震局，福建 福州 350003）

泉州臺地電場數據分析

倪曉寅

（福建省地震局，福建 福州 350003）

通過對泉州地磁臺地電場資料進行處理和分析，認識地電場無震時的變化特征，分析觀測系統的工作狀態，判別各類干擾因素，以此作為識別地震前兆異常的基礎認識，為進一步提高地電場觀測數據質量和開展地電場數據應用研究提供依據。

地電場；干擾；磁暴；小波分析

0 引言

地球表面的天然變化電場是由地球外部的各種電流體系在地球表面感應產生的，分布于整個地表的廣大地區，這種天然的的全球性或區域性的變化電場，稱為大地電場。天然的穩定電場主要由礦體、地下水和各種水系產生，分布于局部地區，這種地方性的天然穩定電場，稱為自然電場。大地電場和自然電場統稱為地電場[1]。自然電場是局部地區由各種物理化學作用產生的地表電場，成因復雜，規律性差。大地電場起源于地球外部空間電流體系，與地磁場的變化有密切的聯系，因而和地磁場一樣具有較強的規律性和周期性[2]。大地電場的觀測與研究已由100多年的歷史，在我國自行研制了ZD9A大地電場儀以后，近幾年大地電場的觀測和地震預報應用研究才較為普遍的開展起來。但由于開展觀測的時間較短，積累的經驗不多，對大地電場的的數據分析和處理研究仍處于初步階段。

本文試通過對泉州臺大地電場分鐘值觀測數據進行處理和分析，認識在無震時大地電場各參量的變化特征，分析觀測系統的工作狀態，判別各類干擾因素，以此作為識別地震前兆異常的基礎認識，為進一步提高大地電場觀測數據質量和開展大地電場數據應用研究提供依據。

1 泉州臺大地電場觀測概況

泉州南安地磁臺地電場觀測場地位于泉州南安市洪籟鎮厝斗村泉州地磁臺內及周邊水田中。2007年9月安裝ZD9A-Ⅱ型地電場儀進行大地電場觀測。電場觀測采用三角形共用中心的布極方式（圖1，O為中心NS向和EW向長短極距共用電極），共有NS向、EW向和NE向三個方向，每個方向又分長、短二種極距。其中長極距（圖1中AO和BO）277m，短極距185m（A′O和B′O；電極埋深3.5m。觀測系統的建設和布極區的環境狀況均符合觀測規范要求，場地及附近沒有高壓線盒任何大型電磁干擾源，坡度小于5%，地形較為平坦。

2 資料處理

為了認識地電場的變化特征，我們以泉州臺地電場分鐘值數據為基礎，進行分析。通過繪制每道測線的每日分鐘值曲線，可以了解地電場在各種情況下（磁靜日、磁暴、雷擊等）的變化特征，據此辨別地電場地震前兆變化。同時地電場是一個矢量，該矢量的方位角可以由電場的南北向分量ESN和EWE東西向分量計算求得α=tg-1EWE45°/ESN45°，并繪制方位角變化曲線分析其動態變化。理論上ESN和EWE在北東方向的投影量值應該與北東向觀測值ENE相同，由于實際觀測系統中的各種因素會引起誤差，該誤差Err=ESNcos45°+EWEcos135°-ENE。同時由于大地電場與起源于天電，則同一方向的長短極距日變化應該是相關的，計算每日長短極距分鐘值的相關系數，可以增加了解觀測系統的觀測質量。最后利用小波分析手段了解地電場的周期特征。

3 地電場變化特征

3.1 磁靜日時的地電場

圖2是泉州臺2008年3月12日地電場NS向長極距、NS向短極距、EW向長極距、EW向短極距、NE向長極距和NE向短極距分鐘值曲線圖。3月12日是磁靜日。六道測線的日變幅均小于39.4mV/km，有明顯的波峰波谷形態；NS向、EW向和NW向的相關系數分別是0.92、0.94和0.91，長短極距觀測到的變化一致性較好，但各測項長短極距的絕對值差異較大；白天的地電場變化幅度明顯大于夜晚；11時-18時的變化幅度最大；極值出現在11時和15時左右。地電場的日變化主要受太陽輻射影響，所以太陽輻射強的白天變化較大，同時，地電場每天的極值也出現在太陽輻射最強的午間，夜晚太陽輻射較弱的時間段，地電場則處于平穩狀態。

圖3 是同一天泉州臺長短極距的地電場方位角和誤差曲線圖。方位角α的變化幅度不大（不到6°），其日變形態與電場日變一致。誤差曲線顯示，長極距誤差幅度主要在11.8mV/km上下波動，誤差最大達到21.5mV/km。短極距誤差主要在21.1mV/km上下波動，誤差最大達到29.8mV/km。誤差在正午時分達到最大。

日變形態清晰和長短極距日變幅相當是泉州臺地電場的主要特征。

3.2 磁暴時的地電場

由于地電場中的一部分場的起源是由外空電流在地球內部產生感應電場引起的，因此磁暴等外空磁場環境的改變會引起大地電場的改變。磁暴發生時，地電場同樣會發生劇烈的變化，這種在磁暴期間記錄到的地電場劇烈變化稱為地電暴。圖4分別是泉州臺2010年4月5日地電場六測道分鐘值曲線（方框為磁暴發生時的地電場）。2010年4月5日北京時間下午16時開始發生了一次中等強度的磁暴（國家空間天氣監測預警中心數據）。從圖中可以看出，磁暴發生后地電場6測道同步記錄到持續的高頻大幅脈沖信號，磁暴發生后NS向長極距最大變化幅度達到了72.2mV/km，是磁靜日時30mV/km日變幅的兩倍以上。

3.3 雷擊

一般情況下地電場高頻變化的場源是對流層中的雷電，主要是雷電改變了區域的電場分布環境[3]。特別是在出現雷雨天氣的情況下，雷電通常離得較近，在放電的瞬間，引起數據變化幅度也較大，持續時間也相對較長。圖5是泉州臺地電場2010年5月19日分鐘值曲線，經落實21時左右開始有雷雨。從圖中可以看出，雷電期間地電場的變化特征是：地電場觀測數據出現快速的大幅度的脈沖變化，數據離散度很大，出現很明顯的畸變，最大變化幅度可達正常地電場日變幅度的數倍。

4 受干擾下的地電場變化

4.1 電極極化

電極是地電場影響觀測質量穩定的重要因素之一。目前地電場觀測使用的固體不極化電極在埋設不好或者老化的情況下會出現極化現象，具體表現主要為數據出現長期的漂移現象。如圖6所示，泉州臺的中心電極（NS向和EW向測道公用電極）極化后，短極距EW數據出現長時間的漂移現象。

4.2 電極滲水

圖7為泉州臺中心電極滲水后的長極距NS向和EW向地電場觀測數據曲線圖。由圖中可以看出，電極滲水后地電場觀測數據出現一個突降臺階，下降幅度接近200mV/km，經過一段時間后有緩慢回升至正常水平。以上現象的具體原因是水滲到電極處，使電極產生一個極化電位差導致的。而隨著時間的推移，該極化電位差慢慢消失，數據恢復到原始水平。

5 地電場的周期成分

通過小波趨勢分析發現（圖8），1階和2階小波分析下的半日波非常明顯；3階和4階能發現比較明顯全日波；5階和6階出現較長周期的半月波。

根據0.Gish編制的1936年世界時18h全球大地電流分布圖，北半球有8個電流渦旋場，中低緯度的4個渦旋中心在30附近，高緯度的4個中心接近70°附近。這些渦旋電流場基本上按地理經度等間隔分布。渦旋電流系的位置固定，白天電流強、夜間電流弱，地球自轉一周，各渦旋電流場的電流強、弱交替一次，白天兩個強電流渦旋場、夜間兩個弱電流渦旋場引起地電場經歷兩次起伏，所以產生大地電場日變化的顯著半日波周期成分[1]。同時大地電場起源于地球外部的電流體系，這些電流體系的產生則源于太陽活動。太陽風到達地球外部空間時，壓縮向陽的地球磁層，而另一面的磁層則向背陽的方向拓展，伴隨著地球的24小時自轉周期，引起了大地電場的24小時周期變化，即全日波變化。

大量的觀測表明，月球的潮汐引力變化同樣會引起大地電場的變化?？拷＿吋昂吹呐_站，受固體潮的影響較為明顯。泉州臺屬于沿海臺站，其記錄到的地電場半月周期成分應源于月球引力引起的固體潮。

6 結論

地電場記錄的資料極為豐富，通過原始曲線和小波分析可以看出，臺站不只記錄到外空場的日變特征，其成分中還包含有固體潮引發的成分。除了外空場和固體潮引起的正常變化，還受到了電極極化和滲水、磁暴和雷電等其他因素的干擾，這些干擾信息形態各異，但都具備各自獨有的形態特征：如磁暴多為持續的大幅度高頻脈沖信號，雷電則表現為瞬時的大幅變化。這些干擾信息的特征在觀測曲線上可以清晰的識別。由于臺站使用的是多極距觀測，所以在識別異常和干擾方面帶來極大的方便，也為分析可能與地震有關的異常信號打下了比較好的基礎。同時我們也發現，共用電極的布極方式雖然節約了成本，但一旦公用電極出現問題，則會造成多測道數據的同時干擾變化。為了更好的記錄到地電場變化，應在單臺增加更多極距觀測。

對泉州臺處理分析，我們發現，泉州臺的地電場觀測數據可以清晰的記錄到真實大地電場的正常日變化，周邊環境基本無干擾，但中心電極在使用一段時間后出現極化和滲水，其NS向和EW向數據受到影響失真，給前兆分析造成了困難。

本文對泉州臺的地電場數據進行了分析計算，使我們臺站的地電場變化特征、觀測狀況和數據可信度有了一個基礎的認識，為識別未來可能出現的地震前兆異常提供的依據。同時本文對大地電場特征認識的研究還不夠深入，有待未來開展進一步的研究工作。

［1］孫正江，王華俊.地電概論[M].北京：地震出版社，1984，3-95．

［2］鄭兆苾，汪雪泉.蒙城臺地電場資料分析[J].地震，2005，25(2)：41-47．

［3］林向東，徐平，等.地電場觀測中幾種常見干擾[J].華北地震科學，2007，25(1)，16-22．

［責任編輯：李書培］

倪曉寅（1981—），男，工程師，主要從事地磁地電數據分析研究工作。