紅山臺(tái)ELF電磁場(chǎng)功率譜密度的時(shí)間分布特征

宋 昭，羅 娜，趙志遠(yuǎn)，高登平，王秀敏，解 真

(紅山基準(zhǔn)地震臺(tái)，河北 邢臺(tái) 054001)

紅山臺(tái)ELF電磁場(chǎng)功率譜密度的時(shí)間分布特征

宋 昭，羅 娜，趙志遠(yuǎn)，高登平，王秀敏，解 真

(紅山基準(zhǔn)地震臺(tái)，河北 邢臺(tái) 054001)

利用紅山地震臺(tái)ELF極低頻觀測(cè)數(shù)據(jù)處理得到不同觀測(cè)分量、不同頻率功率譜數(shù)據(jù)的時(shí)間序列，并通過對(duì)數(shù)據(jù)進(jìn)行作圖分析，得出紅山臺(tái)電磁場(chǎng)功率譜密度的時(shí)間分布具有如下特征： ①電磁場(chǎng)功率譜密度均具有夏強(qiáng)冬弱的季節(jié)分布，呈正弦波形態(tài),夏季變化幅度是冬季的1～10倍不等；②不同頻率上，各測(cè)項(xiàng)功率譜密度無多年趨勢(shì)性變化特征，年際變化幅度小于年內(nèi)1～3個(gè)數(shù)量級(jí)。

ELF；電磁場(chǎng)；功率譜密度；時(shí)間分布；紅山臺(tái)

0 引言

地震的孕育和發(fā)生是一個(gè)極其復(fù)雜的地球物理過程。由于其發(fā)生在地球內(nèi)部,更增加了地震預(yù)報(bào)研究的難度[1]。過去幾十年，通過地面和空間觀測(cè)，記錄到地震前大量的電磁異?，F(xiàn)象，證明地震電磁效應(yīng)是確實(shí)存在的。巖石實(shí)驗(yàn)和理論模型研究也為地震電磁現(xiàn)象的存在提供了一定的試驗(yàn)和解釋依據(jù)。電磁場(chǎng)觀測(cè)不同于僅測(cè)量電場(chǎng)或磁場(chǎng)的傳統(tǒng)電法或磁法，它同時(shí)處理電場(chǎng)和磁場(chǎng)，并利用電磁感應(yīng)原理研究地球電磁場(chǎng)和地質(zhì)問題[2]。

“十五”期間，全國(guó)范圍內(nèi)共安裝了12套ELF大地電磁測(cè)深儀，觀測(cè)電磁場(chǎng)時(shí)間序列、電磁場(chǎng)頻譜、視電阻率、阻抗相位和電阻率結(jié)構(gòu)隨時(shí)間的變化等。ELF大地電磁測(cè)深觀測(cè)產(chǎn)出了大量數(shù)據(jù)，在地震預(yù)測(cè)研究和相關(guān)領(lǐng)域得到了廣泛運(yùn)用[3]。隨著ELF儀器的廣泛布設(shè)，統(tǒng)計(jì)并分析全國(guó)范圍內(nèi)電磁場(chǎng)功率譜密度背景值的時(shí)空分布有著重要意義。它不僅能夠使高質(zhì)量的觀測(cè)數(shù)據(jù)得到充分利用，還將對(duì)儀器的運(yùn)行、維護(hù)提供有益幫助，為前兆事業(yè)的發(fā)展起到一定的推動(dòng)作用。

紅山臺(tái)ELF儀器自2009年7月份正式觀測(cè)以來，運(yùn)行情況較為穩(wěn)定，產(chǎn)出了高質(zhì)量的觀測(cè)數(shù)據(jù)。本文選取了該儀器2011—2015年期間的數(shù)據(jù)進(jìn)行分析處理，試圖找到紅山臺(tái)ELF電磁場(chǎng)功率譜密度的時(shí)間分布背景動(dòng)態(tài)。在提高基層臺(tái)站對(duì)本臺(tái)大量觀測(cè)數(shù)據(jù)的利用、更好地認(rèn)識(shí)觀測(cè)數(shù)據(jù)的同時(shí)，為獲取更大范圍內(nèi)的電磁場(chǎng)背景動(dòng)態(tài)分布乃至識(shí)別電磁前兆觀測(cè)異常奠定基礎(chǔ)。

1 觀測(cè)臺(tái)站和儀器簡(jiǎn)介

1.1 觀測(cè)臺(tái)站概況

紅山基準(zhǔn)臺(tái)(以下簡(jiǎn)稱紅山臺(tái))臺(tái)基穩(wěn)固、觀測(cè)環(huán)境優(yōu)良，遠(yuǎn)離城鎮(zhèn)和村莊等干擾源，加上臺(tái)站工作人員對(duì)觀測(cè)系統(tǒng)維護(hù)到位，使得紅山臺(tái)ELF電磁場(chǎng)記錄數(shù)據(jù)較為完整且數(shù)據(jù)質(zhì)量高，利于分析利用。

1.2 觀測(cè)儀器指標(biāo)

2007年5月由中國(guó)地震局地質(zhì)研究所在紅山臺(tái)敷設(shè)了ELF電磁觀測(cè)，使用俄羅斯生產(chǎn)的ACF-4M型ELF電磁觀測(cè)儀器，觀測(cè)4道水平天然場(chǎng)的電磁信號(hào)，即2道相互垂直電場(chǎng)(E1、E2)和2道相互垂直磁場(chǎng)(H1、H2)，觀測(cè)頻率范圍0.1～800 Hz。電場(chǎng)傳感器(Pb-PbCL2 不極化電極)和磁場(chǎng)傳感器(磁感應(yīng)線圈)采用與正南北、正東西呈45°交角的方式埋入地下，以消除交流供電干擾。傳感器連接線采用埋地(50 cm 深度)敷設(shè)，電極線和磁電纜線埋設(shè)相距2 m以上，不平行。儀器布設(shè)情況如圖1，儀器參數(shù)詳見表1。

圖1 紅山臺(tái)ELF布設(shè)圖

紅山臺(tái)ELF的觀測(cè)分為3個(gè)頻率段：第1頻段(D1)，頻率范圍是0.1～40 Hz，采樣率為160 Hz； 第2頻段(D2)，觀測(cè)頻率1～400 Hz，采樣率1 600 Hz；第3頻段(D3)，觀測(cè)頻率1～800 Hz，采樣率3 200 Hz。觀測(cè)儀器每天進(jìn)行2個(gè)時(shí)段的觀測(cè)：第1觀測(cè)時(shí)段為0：00—6：00，第2時(shí)段為6：30—12：30。每1時(shí)段的觀測(cè)頻段順序均為D3、D2、D1。其中，D3頻段觀測(cè)20分鐘，D2頻段觀測(cè)1小時(shí)，D1頻段觀測(cè)4小時(shí)。所有通道的前置放大器增益均為1。

紅山臺(tái)ELF觀測(cè)儀記錄的是天然電磁場(chǎng)。天然電磁場(chǎng)在極低頻頻段的場(chǎng)源主要為：雷電產(chǎn)生的電磁場(chǎng)、高空電離層擾動(dòng)產(chǎn)生的電磁場(chǎng)、地球上各種人文干擾以及來自震前的電磁輻射等[4]。

表1 紅山臺(tái)ELF儀器布設(shè)參數(shù)

2 數(shù)據(jù)計(jì)算及處理

2.1 從原始時(shí)域到頻率域

DFT即離散傅里葉變換算法，可以將信號(hào)進(jìn)行時(shí)域和頻率域的相互轉(zhuǎn)換。對(duì)于有限長(zhǎng)離散信號(hào)x(n),n=0,1,…,N-1,其DFT定義為：

(1)

FFT是DFT變換的快速算法，它利用了WN的對(duì)稱性和周期性，將N點(diǎn)DFT分解為2個(gè)N/2點(diǎn)的DFT，再分解為4個(gè)N/4點(diǎn)的DFT等一直分解下去，則2M點(diǎn)的DFT均可分解為2點(diǎn)的DFT以簡(jiǎn)化計(jì)算。

ELF的原始觀測(cè)資料是包含多種頻譜的時(shí)間序列數(shù)據(jù)。數(shù)據(jù)處理時(shí)將每天記錄到的原始時(shí)間序列數(shù)據(jù)通過FFT變換，分別得到4個(gè)觀測(cè)項(xiàng)的頻率序列。為更好地從頻率域中描述電磁場(chǎng)信號(hào)，一般通過電磁場(chǎng)的能量來描述電磁場(chǎng)的強(qiáng)度，即功率譜的方法。

2.2 抽取固定頻率組成時(shí)間序列

對(duì)所得功率譜序列按D1頻段0.5 Hz、1 Hz、4 Hz、8 Hz、18 Hz、22 Hz、32 Hz、39 Hz 8個(gè)頻率，D2頻段64 Hz、82 Hz、128 Hz、178 Hz、225 Hz、288 Hz、330 Hz、380 Hz 8個(gè)頻率，D3頻段420 Hz、470 Hz、520 Hz、580 Hz、625 Hz、680 Hz、725 Hz、780 Hz 8個(gè)頻率，共24個(gè)頻率點(diǎn)進(jìn)行提取，得到以上頻率電、磁場(chǎng)每個(gè)觀測(cè)分量的自功率譜數(shù)據(jù)的時(shí)間序列數(shù)據(jù)。

3 數(shù)據(jù)選取和作圖分析

使用紅山臺(tái)ELF儀器觀測(cè)時(shí)段中地磁場(chǎng)相對(duì)平靜的第1觀測(cè)時(shí)段(00：00—06：00)原始時(shí)間序列，運(yùn)用FFT方法對(duì)所選數(shù)據(jù)進(jìn)行計(jì)算最終得到了2011—2015年電磁場(chǎng)4個(gè)觀測(cè)項(xiàng)、24個(gè)固定頻率點(diǎn)的自功率譜密度的時(shí)間分布。由于本文只分析數(shù)據(jù)變化的背景動(dòng)態(tài)和一般趨勢(shì)，功率譜數(shù)據(jù)在分析繪圖前去掉了因觀測(cè)系統(tǒng)故障或人工測(cè)試等原因產(chǎn)生的錯(cuò)誤數(shù)據(jù)和個(gè)別突跳點(diǎn)，以減小其對(duì)曲線形態(tài)的壓制和對(duì)計(jì)算結(jié)果的影響。

3.1.1 同頻率功率譜密度分布

將紅山臺(tái)ELF儀器2015年330 Hz的電磁場(chǎng)H1、H2、E1、E2四個(gè)測(cè)項(xiàng)自功率譜密度繪圖(圖2)，并計(jì)算其各分量在冬、夏不同季節(jié)的平均值、變化幅度和變異系數(shù)(表2)。通過圖2可以直觀地看出，頻率為330 Hz 4個(gè)測(cè)項(xiàng)的電磁場(chǎng)功率譜密度的時(shí)間分布變化趨勢(shì)一致，均具有季節(jié)特征：即曲線夏季升高，冬季回落(文中季節(jié)按照冬季(11-12月,1-2月)，春、秋季(3-4月,9-10月)和夏季(5-8月)的方法劃分)，且以年為周期由夏季到冬季表現(xiàn)出逐漸減小的變化趨勢(shì)，全年呈正弦波形態(tài)。觀察表2發(fā)現(xiàn)，夏季各測(cè)項(xiàng)數(shù)據(jù)的觀測(cè)值、變化幅度和離散度均較大，而冬季較小。夏季變化幅度是冬季變化幅度的1～10倍不等。紅山臺(tái)其他年份和其他頻率數(shù)據(jù)均具有類似特征。

3.1.2 不同頻率功率譜密度分布

通過各分量之間功率譜密度分布比較，發(fā)現(xiàn)磁場(chǎng)H2分量表現(xiàn)出了更加明顯的季節(jié)特征。于是，從紅山臺(tái)2015年ELF電磁場(chǎng)數(shù)據(jù)中選取H2分量D1低頻段的1 Hz、8 Hz，D2中間頻段的128 Hz、288 Hz以及D3高頻段的420 Hz、680 Hz，共6個(gè)頻率點(diǎn)的功率譜時(shí)間序列進(jìn)行繪圖(圖3)。同樣，計(jì)算所選6個(gè)頻率的季節(jié)變化參數(shù)(表3)。通過圖3可以發(fā)現(xiàn)，對(duì)于H2分量來說，所選各頻率功率譜數(shù)據(jù)的總體變化趨勢(shì)一致，均表現(xiàn)出季節(jié)分布。由表3可知， H2高頻數(shù)據(jù)的功率譜平均值比低頻點(diǎn)小1～2個(gè)數(shù)量級(jí)，低頻點(diǎn)夏季變幅是冬季變幅的1～2倍而高頻率則能夠達(dá)到7～10倍。由圖3更加能夠直觀地看出，高頻率數(shù)據(jù)曲線的弧度大且功率譜強(qiáng)度小。分析圖表可以得出結(jié)論：雖然隨著頻率的增高，磁場(chǎng)H2分量的功率譜強(qiáng)度在逐漸減弱，但高頻數(shù)據(jù)卻表現(xiàn)出了更明顯的季節(jié)變化特征。

圖2 2015年紅山臺(tái)ELF電磁場(chǎng)(330 Hz)各分量功率譜曲線

表2 2015年紅山臺(tái)ELF電磁場(chǎng)(330 Hz)4分量功率譜季節(jié)變化參數(shù)

出現(xiàn)季節(jié)變化現(xiàn)象的原因，考慮是夏季地磁場(chǎng)擾動(dòng)較大，地磁場(chǎng)的擾動(dòng)歸因于日地能量的耦合，它反映的是地磁場(chǎng)所包含的能量的變化[5]。這種現(xiàn)象可能是太陽紫外輻射產(chǎn)生的電離層傳播效應(yīng)，也可能是受夏季雷電天氣多發(fā)、太陽活動(dòng)強(qiáng)烈、太陽風(fēng)和磁層及地磁脈動(dòng)等相互作用影響的結(jié)果[6]。而電磁場(chǎng)密不可分，兩者之間的關(guān)系遵循麥克斯韋方程，形態(tài)應(yīng)具有一定相關(guān)性。這也同引文中地電場(chǎng)觀測(cè)資料日變幅的季節(jié)變化特征[7]相吻合。

K指數(shù)作為“磁情指數(shù)”，以3小時(shí)為1時(shí)段, 用0，1，…，9的數(shù)字分10級(jí)來表示地球變化磁場(chǎng)受太陽微粒輻射而形成的擾動(dòng)變化的程度[8]。本文選取了紅山臺(tái)2015年全年的K指數(shù)，將每天0—6時(shí)的2個(gè)K指數(shù)數(shù)據(jù)累加，得到與所選ELF第一觀測(cè)時(shí)段數(shù)據(jù)相對(duì)應(yīng)的K指數(shù)繪圖(圖4a)。由圖4可見，紅山臺(tái)2015年每天0—6時(shí)K指數(shù)和具有同功率譜密度相似的季節(jié)分布。而且，相關(guān)研究發(fā)現(xiàn)ELF電磁場(chǎng)數(shù)據(jù)的變化滯后于地磁擾動(dòng)，表明ELF電磁場(chǎng)功率譜密度的季節(jié)分布受到了磁層活動(dòng)、地磁擾動(dòng)的影響[8]。另外，紅山臺(tái)2015年地磁總場(chǎng)強(qiáng)度F日均值的一階差分曲線也呈現(xiàn)相似的季節(jié)分布(圖4b)，可見，紅山臺(tái)電磁場(chǎng)功率譜密度的季節(jié)變化特征與磁場(chǎng)擾動(dòng)強(qiáng)弱程度的季節(jié)分布相關(guān)。

圖3 2015年紅山臺(tái)ELF觀測(cè)儀H2分量(1 Hz、8 Hz、128 Hz、288 Hz、420 Hz、680 Hz)功率譜曲線

表3 2015年紅山臺(tái)ELF電磁場(chǎng)H2分量功率譜季節(jié)變化參數(shù)

3.2 年際分布

分析紅山臺(tái)ELF電磁場(chǎng)功率譜密度的年際分布特征，即將時(shí)間單位以年計(jì)，計(jì)算年際變化。表4列出了H1分量所有24個(gè)頻點(diǎn)2011—2015年的年際和年內(nèi)變化幅度。由表4可見， 功率譜數(shù)據(jù)的年變幅比年際變化幅度大1～3個(gè)數(shù)量級(jí)。相較于年內(nèi)變化，功率譜數(shù)據(jù)的年際變化可忽略不計(jì)。紅山臺(tái)ELF功率譜其他分量觀測(cè)數(shù)據(jù)均有類似特征。

圖4 2015年紅山臺(tái)地磁K指數(shù)累加值與地磁總場(chǎng)F日均值一階差分曲線

表4 2011—2015年紅山臺(tái)ELF磁場(chǎng)H1分量24個(gè)頻點(diǎn)功率譜年際和年內(nèi)變化幅度對(duì)比

表4中H1分量功率譜年際、年內(nèi)變化趨勢(shì)整體隨頻率升高而減小，但在4～32 Hz頻率段似乎出現(xiàn)年內(nèi)變化幅度較小的異常值。考慮舒曼諧振波的第一階在7.8 Hz左右[9]，其前四階諧振波恰好分別在頻率為7.83 Hz、13.8 Hz、19.7 Hz、25.7 Hz附近，上述異?，F(xiàn)象可能與其相關(guān)。于是繪制紅山臺(tái)ELF數(shù)據(jù)(0～35 Hz)頻譜圖(圖5)。由圖5可見，紅山臺(tái)磁場(chǎng)功率譜數(shù)據(jù)能夠記錄到舒曼諧振的前4階諧振波，頻譜圖分別在頻率為7.83 Hz、13.8 Hz、19.7 Hz、25.7 Hz附近有功率譜增強(qiáng)的現(xiàn)象，且頻率越小，變化幅度越大。同時(shí)可以發(fā)現(xiàn)，電磁場(chǎng)4個(gè)分量在D1頻段的功率譜強(qiáng)度為H1、E1分量?jī)?yōu)于H2、E2，但磁場(chǎng)H2相較H1分量對(duì)諧振波現(xiàn)象有更好的響應(yīng)。這可能源于觀測(cè)系統(tǒng)差異或地下結(jié)構(gòu)的不均勻分布，可將此作為今后進(jìn)一步的研究方向。

圖5 紅山臺(tái)ELF電磁場(chǎng)數(shù)據(jù)(0～35 Hz)頻譜圖

上述分析表明：0～35 Hz頻率范圍內(nèi)，磁場(chǎng)H2能夠更好的響應(yīng)舒曼諧振波。因此，表5給出磁場(chǎng)H2分量位于D1頻率段，且接近諧振波的6個(gè)頻率點(diǎn)年際和年內(nèi)變化幅度結(jié)果。表5反應(yīng)出磁場(chǎng)H2分量位于舒曼諧振前4階諧振波附近的頻率點(diǎn)，其功率譜密度的年際分布變化幅度比年變幅小1個(gè)數(shù)量級(jí)，同樣沒有表現(xiàn)出多年變化趨勢(shì)。

表5 2011—2015年紅山臺(tái)ELF磁場(chǎng)H2分量6個(gè)頻點(diǎn)功率譜年際和年內(nèi)變化幅度對(duì)比

將磁場(chǎng)H2分量中1 Hz、8 Hz、22 Hz 3個(gè)頻率點(diǎn)2011—2015年功率譜密度的年變幅、年均值繪圖，將5年均值變化擬合直線，以突出多年變化長(zhǎng)期趨勢(shì)，并標(biāo)出擬合直線斜率以反映多年變幅(圖6)。圖6可直觀反映出各頻點(diǎn)數(shù)據(jù)年際分布均勻，多年變化幅度相比平均年變幅小1～2個(gè)量級(jí)。數(shù)據(jù)多年變化僅有微小的上升或下降浮動(dòng)，無明顯的趨勢(shì)性變化。

圖6 2011—2015年紅山臺(tái)磁場(chǎng)H2分量(1 Hz、8 Hz、22 Hz)功率譜年變幅和多年趨勢(shì)擬合

4 結(jié)論

本文通過分析紅山臺(tái)ELF極低頻電磁場(chǎng)功率譜密度的時(shí)間分布，得出如下結(jié)論：

1)各頻段的電磁場(chǎng)功率譜密度均具有正弦波形態(tài)的季節(jié)變化特征，表現(xiàn)形式為夏強(qiáng)冬弱。夏季變化幅度是冬季變化幅度的1～10倍不等。這種變化特征與磁場(chǎng)擾動(dòng)強(qiáng)弱的季節(jié)分布相關(guān)。

2)電磁場(chǎng)功率譜密度的年際分布均勻，其年變幅比年際變化幅度大1～3個(gè)數(shù)量級(jí)。相對(duì)于年內(nèi)的季節(jié)變化來說，年際變化幅度可忽略。電磁場(chǎng)數(shù)據(jù)不存在多年趨勢(shì)性變化。

3)紅山臺(tái)電磁場(chǎng)功率譜D1頻率段數(shù)據(jù)能夠記錄到舒曼諧振的前4階諧振波現(xiàn)象。其中，磁場(chǎng)H2分量表現(xiàn)出更好的季節(jié)分布特征和對(duì)舒曼諧振波更好的響應(yīng)。

[1] 郝建國(guó), 張?jiān)聘? 潘懷文, 等. 震前極低頻電磁異常及其頻譜特征[J]. 地震學(xué)報(bào), 1995, 17(1): 81-88.

[2] 趙國(guó)澤, 湯吉, 鄧前輝, 等. 人工源超低頻電磁波技術(shù)及在首都圈地區(qū)的測(cè)量研究[J]. 地學(xué)前緣, 2003, 10(S1): 248-257.

[3] 李細(xì)順, 高登平, 王靜, 等. 極低頻電磁觀測(cè)數(shù)據(jù)噪聲處理方法研究[J]. 山西地震, 2013(4): 18-22.

[4] 孫維懷, 湯吉, 張平, 等. 云南ELF電磁觀測(cè)資料的異常特征分析[J]. 地震地質(zhì), 2012, 34(3): 488-499.

[5] 齊瑋, 王秀芳, 王仁明, 等. 地磁場(chǎng)K變化的信息熵[J]. 地球物理學(xué)報(bào), 2011, 54(3): 780-786.

[6] 廖懷哲, Egeland A. 日側(cè)極隙緯度低頻波及其與地磁擾動(dòng)的關(guān)系[J]. 地球物理學(xué)報(bào), 1995, 38(3): 296-304.

[7] 郭建芳, 李非, 張秀霞, 等. 地電場(chǎng)日變幅與地電暴分析[J]. 地震地磁觀測(cè)與研究, 2010, 31(3): 18-23.

[8] 程安龍, 周錦屏, 黃蔚北, 等. 中國(guó)地磁臺(tái)網(wǎng)的K指數(shù)測(cè)量結(jié)果分析[J]. 地震地磁觀測(cè)與研究, 1998, 19(2): 48-55.

[9] 范曄, 湯吉, 趙國(guó)澤, 等. 電磁監(jiān)測(cè)臺(tái)站觀測(cè)的舒曼諧振背景變化[J]. 地球物理學(xué)報(bào), 2013, 56(7): 2369-2377.

Time-distribution Characteristic of Power-spectrum Density of ELF Electromagnetic Observations in Hongshan Station

SONG Zhao, LUO Na, ZHAO Zhi-yuan, GAO Deng-ping, WANG Xiu-min, XIE Zhen

(Hongshan Standard Seismic Station, Xingtai 054001, China)

In the paper, power-spectrum time series of the extremely low frequency electromagnetic field in different item recorded by Hongshan standard seismic station is calculated and studied using diagraph analysis. The result shows that, 1) all power-spectrum densities in each frequency and observation item have sine shaped seasonal distribution characteristics, mean that, datum are dispersed in summer and dense in winter. The variation in summer is 1 to 10 times as much as in winter.2) power-spectrum densities in each frequency and observation item have no trend change for years, and its interannual variability is 1～3 orders of magnitude less than its annual variability.

ELF; electromagnetic field; power-spectrum density; time-distribution; Hongshan station

宋昭，羅娜，趙志遠(yuǎn)，等.紅山臺(tái)ELF電磁場(chǎng)功率譜密度的時(shí)間分布特征[J]．華北地震科學(xué)，2017，35(1):67-72.

2016-10-13

河北省地震科技星火計(jì)劃項(xiàng)目：ELF電磁場(chǎng)功率譜數(shù)據(jù)的時(shí)空分布特征研究

宋 昭(1989—)，女，河北邢臺(tái)人，助理工程師，現(xiàn)主要從事前兆觀測(cè).E-mail：229243906@qq.com

P315.7

A

1003-1375(2017)01-0067-06

10.3969/j.issn.1003-1375.2017.01.011