地基ELF線天線在地-電離層殼體中產(chǎn)生的場

崔玉國 王元新

(中國電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點實驗室,青島 266107)

?

地基ELF線天線在地-電離層殼體中產(chǎn)生的場

崔玉國 王元新

(中國電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點實驗室,青島 266107)

為了研究地震產(chǎn)生的機理,需要分析和研究地基極低頻(Extreme-Low Frequency,ELF)電磁波的傳播特性．地基任意ELF線天線產(chǎn)生的電磁場可由ELF點源產(chǎn)生的電磁場進(jìn)行疊加．通過地基ELF點源電磁場的表達(dá)式導(dǎo)出了ELF線天線在地-電離層殼體中電磁場的表達(dá)式．利用加速收斂算法得到了單條和兩條ELF線天線產(chǎn)生電磁場的場強空間分布圖,用色標(biāo)表示場強值在空間分布的強弱．電場水平分量的場強值要比垂直分量的場強值小,磁場水平分量的場強值要比垂直分量的場強值大．而且兩條線天線產(chǎn)生的場強值要比單條線天線產(chǎn)生的場強值大.

點源;地-電離層殼體;空間分布圖;色標(biāo)

DOI 10.13443/j.cjors.2016051901

引 言

近年來的觀測研究表明,地震前幾天或幾小時將會出現(xiàn)大量的低頻電磁輻射異常,它提供了明顯的地震前兆信息[1]．地震電磁輻射是在地震的孕育過程中,由地震輻射源產(chǎn)生的電磁信號．它的出現(xiàn)為探索地震預(yù)測提供了新的方法和依據(jù)[2]，它已經(jīng)成為地震日常監(jiān)測預(yù)報特別是短臨預(yù)報最有效的手段之一．地震孕育區(qū)的范圍很大,產(chǎn)生的場極其復(fù)雜,其地下場源實際上是體電流或線電流．雖然地震前會出現(xiàn)全頻段的輻射,但只有低頻輻射到達(dá)地表,所以地震電磁輻射前兆主要能量集中在超低頻(Super-Low Frequency,SLF)和極低頻(Extreme-Low Frequency,ELF)頻段[3]．

為了深入研究地震電磁輻射前兆的規(guī)律,必須進(jìn)一步分析和研究SLF和ELF頻段電波的傳播規(guī)律．對于ELF頻段,電波波長已可與地球周長相比擬,此時,地面與電離層之間的傳播空間不能看作一個“波導(dǎo)”,而只能看作一個“殼體”,現(xiàn)有的SLF電波傳播理論不能用于ELF頻段的電波傳播描述．近幾年很多學(xué)者做了大量的理論研究工作,試圖了解ELF頻段的電波傳播特性,但所得電磁場表達(dá)式中含有難以解析計算的廣義索末菲積分[4-5]．

1999年Donald E.Barrick[6]提出了一種ELF地上垂直電偶極子在地-電離層殼體中的球諧級數(shù)算法,并就地面和電離層都是理想導(dǎo)體條件下,導(dǎo)出了電磁場的解析表達(dá)式,從而算出了電磁場沿地面的變化規(guī)律,但作者沒有進(jìn)一步計算出地面與電離層都是非理想導(dǎo)電時的結(jié)果．有關(guān)ELF頻段地面與電離層都是非理想導(dǎo)電時地上水平電偶極子和垂直電偶極子在地-電離層殼體中的計算結(jié)果本文作者已在前幾篇文章[7-8]中做了詳細(xì)介紹．本文是在Barrick方法和已有工作的基礎(chǔ)上對地基線天線進(jìn)行研究,其產(chǎn)生的場可由點源產(chǎn)生的場進(jìn)行疊加[9]．

1 地基ELF點源在地-電離層殼體中產(chǎn)生的場

1.1 理論計算物理模型

對于ELF頻段,電波波長已與地球周長可比擬,作為最簡單的模型,我們把地面與電離層都理想化為球面均勻銳邊界媒質(zhì),故本文的物理模型如圖1所示．

取球坐標(biāo)如圖1所示,地面與電離層都理想化為一個具有一定表面阻抗的反射壁,分別記為Δg和Δi,它們距球心的距離分別為ra(6 370 km)和rc=ra+h,h是電離層離地高度．激勵源理想化為一個水平電偶極子,其電流距為Idl,放置在(r=rb,θ=0°, φ=0°)處．

圖1 地-電離層殼體坐標(biāo)示意圖

1.2 點源在地層和大氣層中產(chǎn)生場的球諧級數(shù)表達(dá)式

波導(dǎo)中的電磁場可由一對位函數(shù)來表示,即

(1)

(2)

(3)

除了場源所在處外,位函數(shù)應(yīng)滿足

(4)

(5)

由水平電偶極子產(chǎn)生的初始場和次級擾動場可以得到位函數(shù)的表達(dá)式為:

(6)

(7)

將位函數(shù)的表達(dá)式(6)、(7)代入式(1)、 (2)和(3),得到波導(dǎo)中場分量的表達(dá)式為:

(8)

(9)

在r=ra處,電磁場應(yīng)滿足地面的阻抗邊界條件,即

(10)

在r=rc處, 電離層的阻抗邊界條件為

(11)

由上述阻抗邊界條件可以得到各系數(shù),其表達(dá)式可參考文獻(xiàn)[8]．

1.3 加速收斂算法

電場分量Eθ的加速收斂算法的公式為

(12)

2 地基ELF任意線天線在地-電離層殼體中產(chǎn)生的場

2.1 理論計算物理模型

如圖2所示,記天線坐標(biāo)為Ms(xs,ys,0),觀察點坐標(biāo)為Mr(xr,yr, 0),規(guī)定地面上正東方向為x方向,正北方向為y方向,天線上一小單元記為dls,它的方向為:ixdxs+iydys．

圖2 天線單元及坐標(biāo)轉(zhuǎn)換示意圖

觀察點Mr(xr,yr, 0)相對于天線單元中心的距離和方位(正東為x方向)應(yīng)為:

r=[(xr-xs)2+(yr-ys)2]1/2;

(13)

(14)

dls天線單元的走向為

(15)

觀察點相對于dls單元的方位為

φs=φ-ψs.

(16)

2.2 任意線天線在地-電離層殼體中產(chǎn)生場的表達(dá)式

因為dls的走向沿天線有變化,故iθ和iφ對于不同的天線單元方向不一致,在考慮整個天線上電流產(chǎn)生的場時,需轉(zhuǎn)換至直角坐標(biāo)系．

球坐標(biāo)轉(zhuǎn)換到直角坐標(biāo)需作如下變換:

(17)

由此可以得到:

(18)

eφ(s)sinφ)I(s)dls;

(19)

(20)

hz(s)=

3 計算實例

3.1 地基單線天線在地面上產(chǎn)生電磁場的場強空間分布

計算如圖3所示地上線天線在地面上產(chǎn)生電磁場的場強空間分布,取線天線長度100 km,電流為200 A,天線單元如圖2所示．起始點坐標(biāo)A(2,3),天線AB段長度為60 km,天線BC段長度為40 km,角度ψs1=60°,ψs2=30°,位于地下20 km處,頻率為f=10 Hz和f=1 Hz．低電離層離地高度h=70 km,地球半徑a=6 370 km,地面導(dǎo)電率σg=0.001 s/m,電離層導(dǎo)電率σi=0.000 01 s/m．

圖3 單條線天線示意圖

在圖4中,我們分別利用加速收斂算法計算了地基單條ELF線天線產(chǎn)生的電磁場四個分量的場強空間分布．在每個分布圖中,雖然色標(biāo)值是不同的,但都用來表示場強在空間的分布強弱，暗紅色表示場強值最強,深藍(lán)色表示的場強值最弱．從圖4可以看出：電場水平分量的場強值要比垂直分量的場強值?。淮艌鏊椒至康膱鰪娭狄却怪狈至康膱鰪娭荡螅?/p>

圖4 單條線天線產(chǎn)生的電磁場分量的場強空間分布圖

3.2 地基雙線天線在地面上產(chǎn)生電磁場的場強空間分布

計算如圖5所示地基兩條線天線在地面上產(chǎn)生電磁場的場強空間分布,取線天線長度100 km,電流為200 A,天線單元如圖2所示．起始點坐標(biāo)A(2,3),D(-3,4),天線AB段和DE段長度為60 km,天線BC段和EF段長度為40 km,角度ψs1=ψs3=60°,ψs2=ψs4=30°,位于地下20 km處,頻率為f=10 Hz和f=1 Hz．低電離層離地高度h=70 km,地球半徑a=6 370 km,地面導(dǎo)電率σg=0.001 s/m,電離層導(dǎo)電率σi=0.000 01 s/m．

圖5 兩條線天線示意圖

圖6 兩條線天線產(chǎn)生的電磁場分量的場強空間分布圖

在圖6中,我們分別利用加速收斂算法計算了地基兩條ELF線天線產(chǎn)生電磁場四個分量的場強空間分布．在每個分布圖中,用色標(biāo)值表示場強在空間的分布強弱．暗紅色表示場強值最強,深藍(lán)色表示的場強值最弱．從圖6可以看出：電場水平分量的場強值要比垂直分量的場強值??；磁場水平分量的場強值要比垂直分量的場強值大。對比圖4可以看出兩條線天線產(chǎn)生的場強值要比單條線天線產(chǎn)生的場強值大．

4 結(jié) 論

為了研究更接近于實際的地震輻射源在地面上激起的電磁場的空間分布規(guī)律,本文利用加速收斂算法分析和研究了地基ELF線天線在地面上產(chǎn)生電磁場的空間分布規(guī)律．由于天線單元沿天線有變化,故在考慮整個天線上電流產(chǎn)生的場時,需轉(zhuǎn)換至直角坐標(biāo)系．地基任意ELF線天線在地-電離層殼體中產(chǎn)生的電磁場可由地基ELF點源產(chǎn)生的場進(jìn)行疊加,求得ELF點源在地-電離層殼體中產(chǎn)生的電磁場的球諧級數(shù)表達(dá)式,從而可以求得地基任意線天線在半無限大地面空間中任意場點產(chǎn)生的電磁場的場強值．

雖然低頻電磁輻射異常是地震前兆之一,ELF輻射源與地震輻射源也有很大的相關(guān)性,但是通過本文的研究還無法對地震進(jìn)行預(yù)測預(yù)報,目前地震輻射源的機理有待深入研究,還有很多工作要做．本文可為后續(xù)地震輻射源機理和預(yù)報方法研究打下良好基礎(chǔ)．

[1]丁鑒海, 申旭輝, 潘威炎, 等.地震電磁前兆研究進(jìn)展[J].電波科學(xué)學(xué)報, 2006, 21(5)：791-801.

DING J H, SHEN X H, PAN W Y, et al.Seismo-electromagnetism precursor research progress[J].Chinese journal of radio science, 2006, 21(5)：791-801.(in Chinese)

[2]于海雁, 周洪娟, 喬曉林.汶川 MS 8.0 地震前 ELF 異常電磁輻射傳播分析[J].地震學(xué)報, 2010, 32(6)：641-648.

YU H Y, ZHOU H J, QIAO X L.Study on wave propagation of ELF emission anomaly before MS8.0 Wenchuan earthquake[J].Acta seismologica sinica, 2010, 32(6)：641-648.(in Chinese)

[3]丁鑒海, 盧振業(yè), 余素榮.地震地磁學(xué)概論[M].合肥：中國科學(xué)技術(shù)大學(xué)出版社, 2011.

[4]潘威炎.長波超長波極長波傳播[M].成都：電子科技大學(xué)出版社, 2004.

[5]熊皓,等.無線電波傳播[M].北京：電子工業(yè)出版社, 2000.

[6]BARRICK D E.Exact ULF/ELF dipole field strengths in the earth-ionosphere cavity over the Schumann resonance region：idealized boundaries[J].Radio science, 1999, 34(1)：209-227.

[7]WANG Y X, FAN W S, PAN W Y.Spherical harmonic series solution of fields excited by vertical electric dipole in earth-ionosphere cavity[J].Frontiers of electrical and electronic engineering in China, 2008, 3(1)：61-69.

[8]王元新, 彭茜, 潘威炎, 等.SLF/ELF水平電偶極子在地-電離層波導(dǎo)中的場[J].電波科學(xué)學(xué)報, 2007, 22(5)：728-734.

WANG Y X, PENG Q, PAN W Y, et al.The fields excited by SLF/ELF horizontal electric dipole in earth-ionosphere cavity[J].Chinese journal of radio science, 2007, 22(5)：728-734.(in Chinese)

[9]GRADSHTEYN I S, RYZHIK I M.Table of integrals, series, and products[M].NewYork：Academic Press, 1980.

[10]王元新, 金榮洪, 耿軍平.地下ELF線天線在地-電離層殼體中產(chǎn)生的場[J].微波學(xué)報, 2011, 27(1)：26-31.

WANG Y X, JIN R H, GENG J P.Fields excited by the ELF underground linear antenna in earth-ionosphere cavity[J].Journal of microwaves, 2011, 27(1)：26-31.(in Chinese)

Fields excited by groundsill ELF linear antenna in earth-ionosphere cavity

CUI Yuguo WANG Yuanxin

(NationalKeyLaboratoryofElectromagneticEnvironment,ChinaResearchInstituteofRadiowavePropagation,Qingdao266107,China)

To study the mechanism of the earthquake production, we need to study the propagation characteristics of the groundsill extreme-low frequency(ELF) electromagnetic wave.The electromagnetic fields excited by the groundsill ELF arbitrary linear antenna may be the sum of the electromagnetic fields excited by the ELF dot sources.The expressions of electromagnetic fields excited by the groundsill ELF linear antenna in earth-ionosphere cavity are derived from the expressions of the electromagnetic fields excited by the ELF dot sources.The field intensity spatial distribution maps of the electromagnetic fields excited by the ELF single linear antenna and two linear antennas are obtained by the speed numerical convergence algorithm.The strong and weak field intensity values in the spatial distribution are expressed by the color code.The field intensity values of the horizontal electric field component are less than those of the vertical component, whereas, the field intensity values of the horizontal magnetic field component are larger than those of the vertical component.Furthermore, the field intensity values excited by two linear antennas are larger than those excited by the single linear antenna.

dot source;earth-ionosphere cavity;spatial distribution figure;color code

崔玉國, 王元新.地基ELF線天線在地-電離層殼體中產(chǎn)生的場[J].電波科學(xué)學(xué)報,2016,31(5):851-857.

10.13443/j.cjors.2016051901

CUI Y G, WANG Y X.Fields excited by groundsill ELF linear antenna in earth-ionosphere cavity[J].Chinese journal of radio science,2016,31(5):851-857.(in Chinese).DOI：10.13443/j.cjors.2016051901

2016-05-19

中電集團重點實驗室專項基金(A171501023)

TN011

A

1005-0388(2016)05-0851-07

崔玉國 (1984-),男,山東人,碩士,中國電波傳播研究所工程師，研究方向為電磁波傳播、相控陣系統(tǒng)設(shè)計及大型天線陣外場測試技術(shù).

王元新 (1978-),男,山東人,博士,中國電波傳播研究所高級工程師,研究方向為電波傳播、電磁場理論及應(yīng)用．

聯(lián)系人：王元新 E-mail：wyx781215@163.com