河北紅山基準地震臺ELF電磁觀測資料映震能力分析

賈　華，王利兵，高登平，陳凱男，暢國平

(1.河北省地震局邯鄲中心地震臺，河北　邯鄲　056000；2.河北省地震局紅山基準地震臺，河北　邢臺　054000)

?

河北紅山基準地震臺ELF電磁觀測資料映震能力分析

賈華1，王利兵2，高登平2，陳凱男2，暢國平2

(1.河北省地震局邯鄲中心地震臺，河北邯鄲056000；2.河北省地震局紅山基準地震臺，河北邢臺054000)

摘要：為探究ELF電磁觀測中出現的電磁擾動與地震間存在的關系，文章選取河北省地震局紅山基準地震臺2011年至2012年ELF電磁觀測資料進行分析。發現電磁觀測中存在有與地震相關的異常電磁擾動，并將異常信號幅度與地震震級和震中距的關系作進一步的分析。

關鍵詞：功率譜；電磁擾動；ELF觀測

0引言

在中國，從20世紀90年代末開始，開展了基于大功率人工源超低頻/極低頻技術的方法研究，也開展了該方法在地震監測領域的應用研究。與天然場源的電磁波信號相比，人工可控源的大功率超低頻/極低頻電磁波技術具有信號強度大、受電離層的干擾小、傳播穩定可靠、覆蓋范圍廣、可連續定點觀測等優點，并可在數據記錄及后續處理過程中，對信號采用相關的分析技術，大大提高了測量數據的精度，即使在具有較強干擾背景的地區，仍能得到很高信噪比的觀測資料。計算獲得的地下介質電性結構具有很高的精度和穩定性，特別適合于地震預測等需要精確觀測數據的研究[1]。

“十五”期間，全國安裝12套超低頻電磁觀測儀器(ELF)，總體運行良好。2007年5月，由中國地震局地質研究所在紅山基準地震臺(以下簡稱紅山臺)安裝架設ELF電磁觀測儀器(俄羅斯生產的ACF—4M型)，每天連續觀測11小時，所選頻率范圍是：D1:0.1～80 Hz，D2:0.1～800 Hz，D3:0.1～1 600 Hz。

1ELF觀測數據處理

1.1ELF觀測設置

紅山臺ELF儀器觀測的是四分量的電磁場信號：二分量水平電場(正東西Ex和正南北Ey)，二分量的磁場(正東西Hx和正南北Hy)，電場傳感器(電極)和磁場傳感器埋入地下，連接線采用埋地(50 cm深度)敷設，前置放大器增益為1(所有通道均相同)，采用與正南北、正東西45°交角布設儀器，進行四道電磁場監測，觀測相互垂直的水平磁場和電場，頻率范圍為0.1～1 600 Hz，布設參數如表1所示。

表1　紅山臺ELF儀器布設參數

ELF—4M儀器具有計算機網絡設置的功能，使用專用軟件設置儀器的各項參數和自動記錄的開機時間及記錄長度。記錄數據傳輸到計算機后，使用軟件進行處理，對時間序列進行FFT(快速傅式變換)和儀器標定后，得到電場、磁場的自功率譜、互功率譜及視電阻率和阻抗相位曲線。

1.2ELF功率頻譜

ELF儀器的測量方法實際上是大地電磁測深法。利用電磁場電磁感應的趨膚效應，即高頻電磁場穿透淺、低頻電磁場穿透深的特點，通過觀測不同頻段電磁場的頻率，來達到測深目的，將地面上測得的電場、磁場分量，計算出視電阻率和阻抗相位的變化[2]。

對于一個確定性信號或隨機信號，都可以采用譜分析的方法得到信號的譜特征。由于相關函數和功率譜是一對傅里葉變換(FFT)，因而可以先在時間域上求自相關或互相關函數，然后進行傅里葉變換，得到其頻率域的自功率譜或互功率譜，從譜中可以得到帶內的關于信號的各種特征。

SE(f)=∫RE(τ)e-j2πfτdτ，

(1)

式中：RE(τ)是時間域的信號；SE(f)是頻率域功率譜。

對于每天觀測的電磁場信號，均采用以上方法進行計算，得到ELF觀測數據隨頻率的功率譜變化。

在大地電磁測深法中，視電阻率作為一個反映地下介質電性變化的地電參數，建立了表面阻抗與地下介質電阻率之間的關系。一般來說，視電阻率不是某層介質的真正電阻率，而是層狀介質電性參數分布的綜合反映，與地電斷面參數及觀測電磁場信號的周期有關。一定頻段的大地電磁場有一定的穿透深度和影響范圍，而視電阻率表示這一范圍內地電斷面的平均效應。顯然，高頻電磁場分量分布范圍小，平均效應涉及的范圍也小，而低頻電磁場分量的分布范圍大，平均效應涉及的范圍也大，即穿透深度也大[3]。

阻抗相位代表了與之相關的電場水平分量和與之相垂直的磁場水平分量之間的相位差，是頻率的函數，同樣也是反映介質電性特征的一個重要物理量。

2ELF觀測資料映震能力分析

2.1ELF異常對應的地震

在紅山臺2011年至2012年ELF觀測期間，ELF異常共對應中強地震異常7次(見表2)，還存在有地震而無異常對應的情況，可能是由于震級較小或者震中距較大的原因。

表2　2011年至2012年紅山臺ELF觀測對應的中強地震

2.2電磁場功率頻譜的變化分析

將2011年1月1日至2012年4月30日紅山地震臺ELF觀測各頻段的電場(E)和磁場(H)進行數據分析，選用每天連續記錄4 h的D1頻段數據，采用快速FFT，將連續的時間序列數據變換為頻譜序列，經傅里葉變換處理后，可得到電磁場各分量的自、互功率譜密度序列。

使用相應的頻率抽取軟件，對觀測資料處理后的頻譜進行頻率抽取，生成相應頻率的數據文件，繪制出隨時間變化的電磁場變化曲線S(E)和S(H)(見圖1)。通過曲線可知觀測臺站在各種頻率上的電磁場背景變化，從而基于背景的變化和地磁活動，判斷其觀測的電磁場是否存在異常變化[4]。

圖1　紅山臺ELF電磁場D1頻段1 Hz的功率譜曲線Fig.1　Power spectrum in D1 band of 1Hz of ELF electromagnetic field at Hongshan Station

在沒有太陽活動和磁暴的情況下，天然電磁場的變化在一定范圍內如果存在異常，則可能與人為干擾或地球構造運動事件相關。對于測點附近的人為干擾，可通過其他觀測手段或實地調查來判定[5]。

由于電場信號具有更強的局部效應，受臺站地下結構的影響較大，同時，臺站周圍的人為干擾對大地電場也有影響，使得ELF磁場信號往往較電場信號具有更好的信噪比。所以，對于某些地震，磁場信號可以記錄到而電場信號卻不能，因而對于電場功率譜信號也存在一些有震無異常的情況(見第44頁圖2)。

以上兩圖中，選用D1頻段的1 Hz和32 Hz的數據作為分析對象，從每天經過處理的功率譜密度中抽取1 Hz和32 Hz的數據文件，得到其在時間序列上的功率譜變化。用1 Hz和32 Hz頻率的電場、磁場自功率譜繪圖，得到該頻率隨時間變化的曲線，發現在多次地震臨震前存在明顯的電磁能量增強的特點。

2011年2月15日臺灣5.1級地震前，1 Hz南北磁場和東西磁場都比正常值高出兩個數量級，32 Hz磁場和電場功率譜的變化打破正常形態，長時間、大幅度下降。2012年3月4日新疆6.0級地震前，1 Hz和32 Hz功率譜都出現較大異常，磁場與電場均比正常值高出3～6個數量級。分析以上地震前的地磁Kp指數表明，在觀測到異常的時段并沒有出現較強的地磁活動，屬于磁靜日。這說明紅山臺ELF觀測到的電磁場異常變化不是來自外空間電磁活動，可能是由于地球內部變化引起的。

圖2　紅山臺ELF電磁場D1頻段32 Hz功率譜曲線Fig.2　Power spectrum in D1 band of 32Hz of ELF electromagnetic field at Hongshan Station

對觀測到的7次地震事件與ELF 資料的對應變化分析得出，ELF觀測的電磁場計算電場的自功率譜S(E)、磁場功率譜S(H)對于震中距在2 000 km以上，震級Ms>5.5的新疆地震，震前全部參量有異常，異常出現的時間為震前48 h。S(E)異常幅度達到正常時段平均值的5～6個數量級，S(H)異常幅度達到3～4個數量級。對于震中距在0～400 km，震級Ms>4.0的山西、河南地震，震前全部參量有異常，異常出現的時間為震前24 h，S(E)異常幅度達到0.5～0.7個數量級，S(H)異常幅度達到0.5～0.8個數量級(見表3)。

ELF功率譜各參量的變化幅度不僅與信號的頻率、布極方位有關，還與震級大小、震中距等有關。新疆6.0級地震及余震前后，各參量出現明顯跳動。其中，高頻(32 Hz)比低頻(1 Hz)變化幅度強，說明接收的電磁場信號強弱具有頻率可選性。山西4.2級與河南4.3級地震前，1 Hz、32 Hz磁場信號H2(EW向)比H1(SN向)變化強，電場信號E2(EW向)比E1(SN向)變化強，說明電磁場強弱具有方向性。

表3　紅山臺ELF資料功率譜異常幅度與對應地震

對以上地震的視電阻率變化特征進行分析，由于存在50 Hz工頻干擾，使整個曲線全部壓縮，不易看出震前、震時、震后的視電阻率變化，需要經過數字信號處理，才可得到明顯的變化結果，不在此詳述。

3結語

對ELF觀測資料的分析表明，震前的電磁效應明顯。電磁異常不僅在強震震源區可以記錄到，在遠場的記錄中也有反映。紅山臺記錄的是天然電磁場，所以會存在一些人為的電磁干擾或雷電產生的電磁場，因而電、磁場自功率譜曲線中出現的許多突跳并沒有對應地震發生，屬于正常現象。當在不同頻率的電磁場隨時間變化曲線出現異常后，確認資料處理是否準確無誤，周圍環境是否存在人為的電磁干擾，再與同時段臺站的其他觀測手段、其他觀測臺站的記錄進行比較，確認異常類型，從而達到利用ELF電磁觀測資料進行地震短臨預報的目的。

ELF電磁觀測數據有效的處理方法和地震短臨信息提取技術，對研究地震前兆異常與中強地震的關系，對未來地震監測、地震短臨預測方法的拓展具有積極意義。

參考文獻：

[1]王蘭煒，趙家騮，張世中.LF/ELF電磁接收機研究及觀測實驗[J].地震地質，2010，32(3)：482-491.

[2]高曙德，杜學彬，郝臻.低頻電磁波在汶川強余震時的同震效應[J].地震地磁觀測與研究，2009，30(6)：1-7.

[3]趙國澤，王立鳳，湯吉.地震監測人工源極低頻電磁技術(CSELF)新試驗[J].地球物理學報，2010，53(3)：479-486.

[4]趙國澤，湯吉，鄧前輝.人工源超低頻電磁波技術及在首都圈地區的測量研究[J].地學前緣，2003，10(8)：248-257.

[5]吳小平，柳建新，雒志峰.一種多尺度電磁測深方法[J].地球物理學進展，2010，25(5)：1 670-1 675.

Earthquake Reflecting Ability of ELF Electromagnetic Data at Hongshan Reference Seismological Station

JIA Hua1, WANG Li-bing2, GAO Deng-ping2, CHEN Kai-nan2, CHANG Guo-ping2

(1.Handan Central Seismic Station of Earthquake Administration of Hebei Province, Handan, Hebei 056000, China;2.Hongshan Reference Seismic Station of Earthquake Administration of Hebei Province, Xingtai, Hebei 054000, China)

Abstract:ELF electromagnetic data at Hongshan Reference Seismological Station from 2011 to 2012 are chose to research the relationship between electromagnetic disturbance in ELF electromagnetic observation and earthquake. It is found that there are abnormal electromagnetic disturbance associated with the earthquake in electromagnetic observation. And the relationship between abnormal signal amplitude and the earthquake magnitude and epicentral distance is analyzed.

Key words:Power spectrum; Electromagnetic disturbance; ELF observation

中圖分類號：P315.72

文獻標志碼：A

作者簡介：第一賈華(1983—)，女，河北省邯鄲人。2009年畢業于河北工業大學，碩士研究生，工程師。

基金項目：河北省地震局地震科技星火項目(DZ20140707009)。

收稿日期：2015-03-10

文章編號：1000-6265(2016)01-0042-04