地面甚低頻輻射滲透進(jìn)電離層的數(shù)值模擬分析

趙庶凡 廖力 張學(xué)民 申旭輝

(1.中國地震局地震預(yù)測研究所,北京 100036;2.中國地震局地球物理研究所,北京 100081;3.中國地震局地殼應(yīng)力研究所,北京 100085)

?

地面甚低頻輻射滲透進(jìn)電離層的數(shù)值模擬分析

趙庶凡1廖力2張學(xué)民1申旭輝3

(1.中國地震局地震預(yù)測研究所,北京 100036;2.中國地震局地球物理研究所,北京 100081;3.中國地震局地殼應(yīng)力研究所,北京 100085)

利用全波解模型進(jìn)行數(shù)值模擬實(shí)驗(yàn),研究了具有不同輻射源參數(shù)(輻射頻率和功率)的地面甚低頻(Very Low Frequency,VLF)輻射源在不同地磁場參數(shù)(地磁場強(qiáng)度和傾角)和電離層參數(shù)(電子密度和碰撞頻率)條件下激發(fā)的地球-電離層波導(dǎo)和電離層中的電磁場能量空間分布,并重點(diǎn)研究了電離層D/E區(qū)對電磁輻射能量的吸收.模擬結(jié)果發(fā)現(xiàn):VLF輻射在波導(dǎo)中衰減只受輻射源頻率的影響,不受輻射功率、地磁場參數(shù)和電離層參數(shù)變化影響,波導(dǎo)中的衰減隨頻率的增大而減小,而電離層D/E區(qū)吸收隨頻率增大而增大,兩者總衰減量隨頻率增大而增大.輻射源功率對電離層D/E區(qū)的吸收也無影響.地磁強(qiáng)度和地磁傾角越大,電離層D/E區(qū)吸收越小;電離層碰撞頻率和電子密度越大,電離層中能量衰減越大.

地基VLF人工源;全波解;電離層吸收

DOI 10.13443/j.cjors.2016021701

引 言

全球分布的地面甚低頻(Very Low Frequency, VLF)人工源產(chǎn)生的電磁輻射能滲透進(jìn)入電離層乃至磁層,對磁層輻射帶高能粒子產(chǎn)生影響導(dǎo)致其沉降[1].衛(wèi)星觀測的VLF人工源信噪比變化源于地震導(dǎo)致電離層參數(shù)的異常增大或減小[2].因此構(gòu)建地面VLF輻射滲透進(jìn)電離層的傳播模型,研究不同地面VLF輻射源在電離層中激發(fā)的電磁能量分布特征,以及電離層參數(shù)等對能量分布的影響,有助于進(jìn)一步研究粒子沉降和地震電離層電磁異常的物理機(jī)制.此外,還將有助于中國電磁試驗(yàn)衛(wèi)星電磁場觀測數(shù)據(jù)的可靠性驗(yàn)證和數(shù)據(jù)質(zhì)量評價(jià).

當(dāng)VLF波穿透電離層傳播時(shí),由于電子與中性大氣成份的頻繁碰撞會導(dǎo)致VLF波能量產(chǎn)生損耗,計(jì)算這種能量損耗多基于Helliwell[3]的吸收曲線(Helliwell absorption curve)進(jìn)行估算,該模型使用電離層電子密度模型的總體平均值給出2 kHz和20 kHz波的能量吸收估計(jì)值,而實(shí)際電離層D/E區(qū)電子密度隨緯度和季節(jié)變化,因此使用Helliwell吸收曲線時(shí)需要研究的問題是吸收值隨電離層D/E區(qū)電子密度變化將如何變化.文獻(xiàn)[4]利用火箭探測的電離層D區(qū)電子密度數(shù)據(jù),采用全波解方法(Full-wave method)計(jì)算,證實(shí)了Helliwell估計(jì)的吸收值偏大.

在VLF頻段,電離層D/E區(qū)介質(zhì)參數(shù)變化在垂直向上一個波長范圍內(nèi)變化劇烈,在水平方向上變化很平緩,所以將它理想化為水平分層模型,采用全波解的方法進(jìn)行求解.全波方法從文獻(xiàn)[5-6]的理論發(fā)展而來,至今已有大量的水平分層介質(zhì)場的計(jì)算方法[7-10].文獻(xiàn)[11]對這些方法進(jìn)行了總結(jié),指出這些方法大多在保持?jǐn)?shù)值計(jì)算穩(wěn)定性方面存在問題,當(dāng)耗散波的解(垂直波數(shù)虛部大的解)淹沒了我們感興趣的波模時(shí),需要采用特殊的技巧來維持計(jì)算的穩(wěn)定性[5,9].文獻(xiàn)[10]提出采用格蘭-施密特正交化方法克服傳播矩陣迭代過程中產(chǎn)生的數(shù)值溢出,可以對電離層中的總波場進(jìn)行自然分解得到左旋不可傳播模(L模式)和右旋可傳播模(R模式)的特征波波模.Helliwell經(jīng)典吸收曲線估算的就是電離層對R模式的吸收.文獻(xiàn)[12]的全波解方法源于文獻(xiàn)[7]波導(dǎo)模式計(jì)算方法的啟發(fā),基于文獻(xiàn)[9]提出的反射系數(shù)和波模振幅遞歸計(jì)算方法提出的,具有很好的內(nèi)在穩(wěn)定性可以克服“數(shù)值淹沒”問題,用于計(jì)算任意結(jié)構(gòu)的源在地球-電離層波導(dǎo)中任意水平距離和方向上的場,以及計(jì)算電離層中以哨聲波模向上傳播的電磁場強(qiáng)度.

作者利用文獻(xiàn)[10]的全波解方法計(jì)算了從電離層底邊界入射的單位平面波在電離層中的衰減,研究了不同電離層中激發(fā)電磁響應(yīng)的時(shí)變特征,模擬結(jié)果與DEMETER衛(wèi)星數(shù)據(jù)分析結(jié)果較為一致[13].此外,作者在模型中加入離子將文獻(xiàn)[13]的全波方法推廣到超低頻(Ultra Low Frequency, ULF)頻段,對震前電離層參數(shù)異常引起的衛(wèi)星高度的電場變化特征進(jìn)行了模擬研究[14].本文將利用文獻(xiàn)[12]的全波解方法構(gòu)建地面VLF輻射的空間傳播計(jì)算模型,計(jì)算球面波激發(fā)的電磁場沿磁子午面的二維分布,進(jìn)而模擬分析具有不同輻射源參數(shù)(輻射頻率和功率)的地面VLF輻射源在不同地磁場參數(shù)(地磁場強(qiáng)度和傾角)和電離層參數(shù)(電子密度和碰撞頻率)下激發(fā)電磁場分布的差異并總結(jié)其規(guī)律.

1 計(jì)算方法

高電離層電子密度變化比較緩慢,但是低電離層,特別是在VLF波段,電子密度在垂直向上一個波長范圍內(nèi)變化劇烈,在水平方向上變化很平緩,所以將它理想化為水平分層模型.這里將電離層看成是隨高度變化的不均勻的水平分層各向異性冷等離子體,對低電離層進(jìn)行分層,每層的邊界為zi(i=0,1,…,N),z0為地表,z1為電離層底邊界,zN為電離層頂邊界.每一層中認(rèn)為電離層參數(shù)如電離層介電張量沒有變化[15].此外,由于VLF頻率遠(yuǎn)高于離子回旋頻率,因此通常情況不需考慮離子的效應(yīng),本文只考慮電子的影響.本文將電離層磁導(dǎo)率設(shè)為1,場的時(shí)諧因子為ejωt.

低電離層每一層i的介電系數(shù)張量形式上相同,隨高度變化:

εi=ε0εri=ε0(I+Mi).

(1)

式中：εri為相對介電張量;矩陣Mi是第i層(zi

(2)

圖1 計(jì)算使用的坐標(biāo)系

任意指向的源都滿足

(3)

式中,li、mi、ni分別是源的指向在x、y、z方向上的方向余弦.對于地面VLF人工源,可以將其理想化為位于地表的垂直電偶極子源,即li=0, mi=0, ni=1,所以麥克斯韋方程組變?yōu)?/p>

(4)

通過傅里葉變換式對麥克斯韋方程式(4)展開并消去z方向上的分量得到考慮源項(xiàng)的矩陣形式方程組:

(5)

本文采用全波解方法[12]求解方程組(5).在本人其他文章中已有詳細(xì)闡述及驗(yàn)證[16],這里不做具體介紹.

2 模擬結(jié)果與討論

在計(jì)算時(shí)分別將65 km作為電離層底邊界,65～250 km為水平分層各向異性電離層,分層厚度為1 km.利用IRI模型[17]計(jì)算電子數(shù)密度在65～250 km的剖面.碰撞頻率采用Cummer[18]提出的模型.利用國際地磁場(International Geomagnetic Reference Field, IGRF)模型[19]計(jì)算地磁場強(qiáng)度和地磁場傾角.當(dāng)觀察點(diǎn)離場源水平距離不太遠(yuǎn),本文模型較為合理;若水平距離很大則必須考慮地球曲率,同時(shí)參照文獻(xiàn)[20]的模擬實(shí)驗(yàn),故本文模型水平方向上的計(jì)算范圍在±1 500 km范圍以內(nèi).如前文所述坐標(biāo)系選取,模擬結(jié)果中水平軸即x軸指向地磁北,縱軸即z軸指向垂直向上方向.

2.1 不同輻射頻率對電磁響應(yīng)分布的影響

VLF人工源或者閃電激發(fā)的電磁波具有不同頻率,因此需要考慮不同輻射頻率的電離層電磁響應(yīng).本小節(jié)模擬輻射功率為1 000 kW的地面垂直電偶極子,輻射頻率分別為20,10,5 kHz激發(fā)的電磁波能量的空間分布.圖2左列是三個不同頻率的輻射源激發(fā)的坡印廷能流在磁子午面內(nèi)的二維分布,右列是提取的地球-電離層波導(dǎo)中(地表0 km和地球-電離層波導(dǎo)上邊界65 km)以及電離層中(120 km和250 km)的能流水平分布.

由圖2左列可見輻射源的頻率對地球-電離層波導(dǎo)中以及電離層中的能流分布的形態(tài)和大小有較大影響.隨著輻射頻率降低,波導(dǎo)和電離層中的能量減小.滲透進(jìn)電離層的能流聚集成束向特定方向傳播,能量束在輻射源處的磁力線指向一側(cè),但是隨著頻率減小,原本集中在輻射源一側(cè)的能量束變成兩束.

從右列圖綠線和黑線所示波導(dǎo)的上下邊界處的能量分布可見波導(dǎo)中的場有清晰的波模干涉現(xiàn)象,且波導(dǎo)上下邊界的能量變化趨勢剛好相反.從電離層中120 km與250 km處的坡印廷變化曲線可見,兩個高度曲線形態(tài)及大小基本一致,只是水平位置向左錯動,射線仰角為向地理北57°,與地磁場傾角-55°基本一致,說明250 km處的能流是由120 km處的場沿地磁場方向向上映射的,也就是說VLF電磁波在穿透120 km以下的低電離層后,將以哨聲模形式基本無損耗地沿磁力線向高電離層傳播.隨著頻率降低,由第三行圖可見5 kHz時(shí)能流的變化形態(tài)從輻射源附近的單峰結(jié)構(gòu)變成兩個場強(qiáng)峰值,對照左列即電離層中的能量束從單束變成雙束.

圖2 不同頻率輻射源的坡印廷能流隨高度和水平傳播距離的分布(從上到下依次為20, 10, 5 kHz)

為了更清晰地了解輻射源頻率對激發(fā)的能量空間分布的影響,計(jì)算了30,25,20,15,10,5 kHz頻率的能量空間分布,并從計(jì)算結(jié)果中提取不同高度的坡印廷能流的極值,給出其隨高度的變化,如圖3所示.由圖可見電磁輻射從地表向電離層高度傳播的過程中能量不斷衰減,主要能量損耗集中在地球-電離層波導(dǎo)中(0 ～65 km)以及120 km以下的低電離層D/E區(qū),當(dāng)電磁波穿透電離層D/E區(qū)后,傳播過程中能量基本不再衰減.

圖3 不同頻率輻射源的坡印廷能流極值隨高度的變化

2.2 不同輻射功率對電磁響應(yīng)分布的影響

VLF人工源或者閃電激發(fā)的電磁輻射具有不同的能量,本節(jié)計(jì)算不同輻射功率下的電離層電磁響應(yīng).由于電流矩是計(jì)算中使用的激勵源項(xiàng),因此需要確定電流矩與輻射功率的對應(yīng)關(guān)系.對于垂直電偶極子天線,輻射功率通過以下公式算出對應(yīng)的電流矩.真空中,電流矩Il和輻射功率P之間的關(guān)系式如下:

(6)

當(dāng)?shù)厍虮砻婵醋魇抢硐雽?dǎo)體時(shí),必須考慮鏡像電流并且輻射只在地表以上的半空間,因此可以指定電流矩為

(7)

本節(jié)模擬輻射頻率為20 kHz,輻射功率為1 000,500,100 kW的坡印廷能流空間分布,如圖4所示.可見輻射源的功率不影響波導(dǎo)中以及電離層中的能流空間分布形態(tài),只影響能量的大小;功率越低,波導(dǎo)和電離層中的能量顯著減小.

圖4 不同輻射功率下偶極子源的坡印廷能流空間分布(從上到下依次為1 000，500，100 kW)

圖5是不同輻射功率下,不同高度能流密度的最大值隨高度的變化曲線,可見輻射功率越低,在地球-電離層波導(dǎo)和電離層中激發(fā)的能流越小.波導(dǎo)中能量的衰減和電離層D/E區(qū)吸收,與輻射源功率無關(guān).

圖5 不同輻射功率下每一高度坡印廷能流極值隨高度的變化曲線

2.3 不同地磁場強(qiáng)度和地磁傾角對電磁響應(yīng)的影響

VLF人工源和閃電是全球分布的,不同位置的地磁強(qiáng)度和地磁傾角也不同,本節(jié)將模擬VLF人工源在不同地磁強(qiáng)度和地磁傾角下滲透到電離層的電磁場分布.首先將地磁場傾角固定為-55°,計(jì)算地磁場強(qiáng)度分別為5.3×10-3,5.3×10-4,5.3×10-5,5.3×10-6T時(shí)坡印廷能流的空間分布,如圖6和7所示.可見地磁場強(qiáng)度大小對地球-電離層波導(dǎo)中能量的大小和分布沒有影響,主要影響電離層D/E區(qū)的吸收即滲透進(jìn)入電離層的能量.結(jié)果表明電離層中能量滲透的集中區(qū)域在輻射源處的磁力線指向一側(cè),隨著地磁場強(qiáng)度的減小,電離層中的能量越小,即電磁輻射越難進(jìn)入電離層.

由圖7可見,當(dāng)?shù)卮艌龊苄?.3×10-6T時(shí),電離層D/E區(qū)對電磁輻射的吸收作用非常強(qiáng)烈.地磁場強(qiáng)度增加時(shí),電離層D/E區(qū)吸收變小,但是當(dāng)?shù)卮艌鋈≈翟酱蟪^5.3×10-4T后,滲透進(jìn)電離層中的能量大小基本不再隨地磁場的增大而繼續(xù)增大.這可能是由于當(dāng)?shù)卮艌鰪?qiáng)度為5.3×10-6T時(shí)電子回旋頻率為1.5×105s-1,遠(yuǎn)小于電離層的碰撞頻率(量級107s-1),此時(shí)電離層碰撞起主導(dǎo)作用導(dǎo)致能量被吸收無法滲透進(jìn)入電離層.當(dāng)?shù)卮艌鰪?qiáng)度為5.3×10-4T即回旋頻率為1.5×107s-1時(shí),與電離層D/E區(qū)碰撞頻率(量級107s-1)可相比擬,碰撞不起主導(dǎo)作用,隨著地磁場強(qiáng)度繼續(xù)增加D/E區(qū)的吸收保持不變.

圖6 不同地磁場強(qiáng)度下偶極子源的坡印廷能流空間分布(從上到下依次為5.3×10-3，5.3×10-4，5.3×10-5，5.3×10-6 T)

圖7 不同地磁場強(qiáng)度下每一層坡印廷能流最大值隨高度的變化曲線

圖8 不同地磁場傾角下偶極子源的坡印廷能流空間分布(從上到下依次為88°，-9°，-30°，-77°)

圖9 不同地磁場傾角下每一層坡印廷能流最大值隨高度的變化曲線

接著將地磁場強(qiáng)度設(shè)定為5.3×10-5T,計(jì)算地磁場傾角分別為88°,-9°,-30°,-77°時(shí)地表垂直電偶極子激發(fā)的電磁場空間分布,如圖8和圖9所示.由圖8可見在計(jì)算范圍內(nèi)地球-電離層波導(dǎo)中的能流分布的形態(tài)和大小不受地磁場傾角的影響.但是,地磁場傾角對滲透進(jìn)電離層中的能量分布形態(tài)起著重要的作用,模擬結(jié)果得到不同地磁場傾角下能量束的角度分別是89°,0°,-34°,-78°,可見地磁場方向基本決定了電離層中電磁能量束的方向,即能流傾向于沿著磁力線導(dǎo)管傳播,尤其在緯度較高的地區(qū).在我們的模擬試驗(yàn)中,在磁赤道附近,能量基本無法進(jìn)入電離層全部被限制在地球-電離層波導(dǎo)內(nèi),還發(fā)現(xiàn)在極區(qū)附近電離層中的能量束由輻射源一側(cè)的一束變成以輻射源為中心對稱的兩束.

此外由圖9可見,地磁場傾角對波導(dǎo)中高度向衰減基本沒有影響,主要影響電離層D/E區(qū)的吸收;地磁場傾角越小,電離層D/E區(qū)對電磁輻射的吸收作用越強(qiáng)烈.

2.4 不同電子密度剖面和碰撞頻率剖面對電磁響應(yīng)的影響

眾所周知,電子密度的分布特征具有區(qū)域性和時(shí)間性[16,25],本節(jié)模擬不同電子密度剖面對電磁響應(yīng)空間分布的影響.模擬參數(shù)設(shè)定地面的垂直電偶極子源的輻射功率為1 000 kW,輻射頻率為20 kHz,仍選取地磁場強(qiáng)度為5.3×10-5T,地磁場傾角為-55°.將之前模擬時(shí)使用的電子密度稱為“正?！?研究電子密度“正?！?在“正?！被A(chǔ)上增大一個量級,和在“正?！被A(chǔ)上減少一個量級3種剖面下,不同高度的電磁能量最大值變化,模擬時(shí)碰撞頻率保持不變,結(jié)果如圖10所示.類似地,在模擬碰撞頻率的影響時(shí),電子密度保持不變,模擬碰撞頻率“正?！?在“正?！被A(chǔ)上增大一個量級,和在“正?！被A(chǔ)上減少一個量級3種情況,得到同樣的結(jié)果.電子密度和碰撞頻率對波導(dǎo)中能量大小和分布沒有影響,主要影響電離層D/E區(qū)的吸收.電子密度和碰撞頻率越小,電離層D/E區(qū)對電磁輻射的吸收作用越小.

圖10 不同電子密度剖面下每一高度坡印廷能流極值隨高度的變化曲線

2.5 在地球-電離層波導(dǎo)和低電離層中的衰減

由前述幾節(jié)可知,不同特性的VLF輻射源在不同地磁場參數(shù)、電離層參數(shù)下,在地球-電離層波導(dǎo)中高度向的衰減,以及在低電離層(D/E區(qū))的吸收不盡相同,因此最終滲透進(jìn)電離層的能量也不同.我們將地-電離層波導(dǎo)上下邊界能流密度最大值之差定義為波導(dǎo)中高度向的衰減,將波導(dǎo)上邊界和120 km高度的能流密度最大值之差定義為電離層D/E區(qū)的吸收,得到的兩種衰減隨輻射頻率、輻射功率、地磁場大小和傾角、電離層電子密度和碰撞頻率的變化如圖11所示.

圖11 地球-電離層和電離層D/E區(qū)的衰減隨輻射頻率和功率、地磁場大小和傾角、電離層電子密度和碰撞頻率的變化

輻射功率、地磁場參數(shù)和電離層參數(shù)變化對波導(dǎo)中高度向電磁波能量的衰減基本無影響,波導(dǎo)中衰減主要受輻射源頻率的影響,頻率越大,波導(dǎo)中能量高度向的衰減越小,而電離層D/E區(qū)的衰減越大.輻射源頻率小于15 kHz左右時(shí),D/E區(qū)的吸收較小,大于15 kHz左右時(shí),D/E區(qū)吸收比波導(dǎo)中的衰減大.輻射源功率不但對波導(dǎo)中的衰減沒有影響,對D/E區(qū)吸收也無影響,D/E區(qū)的吸收比波導(dǎo)中衰減大約8 dB.當(dāng)?shù)卮艌鰪?qiáng)度很小的時(shí)候,D/E區(qū)的吸收極大,電磁波無法穿透低電離層,隨著地磁強(qiáng)度增大,D/E區(qū)中的吸收逐漸減小,但當(dāng)?shù)卮艔?qiáng)度增大到10-5T的量級之后,D/E區(qū)中的吸收就不再隨地磁場強(qiáng)度增加而減小.D/E區(qū)吸收在磁赤道處非常大,隨地磁傾角增大而減小.電離層D/E區(qū)吸收隨電子密度、碰撞頻率的增大而增大.

為了研究波導(dǎo)中高度向的衰減和D/E區(qū)吸收引起的總衰減量,給出輻射源在地表激發(fā)的能流以及經(jīng)過波導(dǎo)衰減和D/E區(qū)吸收后進(jìn)入電離層中的能流大小,如圖12所示.由圖可見輻射頻率越高,在地表激發(fā)的能流越大,滲透到電離層的能流也越大,但是由兩者的差值可見,總的能流損耗也略有增大.

圖12 不同輻射頻率在地表和電離層中產(chǎn)生的最大能流

3 結(jié) 論

本文建立了地面VLF電磁輻射滲透電離層傳播的二維全波解模型,可計(jì)算地面各種類型的VLF輻射源在空間激發(fā)的電磁能量分布,用于進(jìn)一步研究VLF人工源對輻射帶高能粒子沉降,探索地震電離層電磁異常的物理機(jī)制.利用該模型研究了不同輻射源參數(shù)(輻射頻率和功率)的地面VLF輻射源,在不同地磁場參數(shù)和電離層參數(shù)條件下激發(fā)的電磁能量分布特征,并重點(diǎn)關(guān)注了電離層D/E區(qū)吸收,得出以下結(jié)論:

1) 輻射源的頻率在控制地球-電離層波導(dǎo)中以及電離層中的能流分布的形態(tài)和大小方面占有重要作用.頻率越低,波導(dǎo)中能量高度向的衰減越大,而電離層D/E區(qū)對其吸收作用卻越小.輻射功率與能量分布呈線性關(guān)系.

2) 地磁場強(qiáng)度和地磁傾角,電離層電子密度和碰撞頻率對地球-電離層波導(dǎo)中的能量的大小和分布形態(tài)都沒有影響,主要影響電離層D/E區(qū)的吸收,說明在研究VLF人工源在波導(dǎo)中的近距離傳播時(shí),忽略地磁場以及電離層參數(shù)等波導(dǎo)上邊界參數(shù)的影響是合理的,而研究滲透進(jìn)電離層的傳播必須考慮地磁場的影響.

3) 地磁場傾角對滲透進(jìn)電離層中的能量分布形態(tài)以及電離層D/E區(qū)的吸收起著重要的作用.地磁傾角越小,電離層中的能量束相對輻射源的偏角越小,電離層對電磁波的吸收增大,因此緯度越高，電磁輻射越容易穿透進(jìn)入電離層,靠近磁赤道時(shí),能量基本無法進(jìn)入電離層.

4) 電離層電子密度和碰撞頻率越小,電離層D/E區(qū)對電磁輻射的吸收作用越小,進(jìn)入電離層的能量越大,所以夜間電磁輻射更容易穿透進(jìn)入電離層.

本文的二維模型是在磁子午面內(nèi),主要研究電磁輻射在經(jīng)度向(南北)上的分布特征,為更清晰了解電磁輻射在整個三維空間的分布特征，還需要建立三維全波模型.

[1]UMAN M A,HORNSTEIN J.The lightning discharge[J].Physics today,1987, 42 (9)：443-445.

[2]CHRISSAN D A, FRASER-SMITH A C.Seasonal variations of globally measured ELF/VLF radio noise[J].Radio science, 1996, 31(5)：1141-1152.

[3]郄秀書, 周筠珺, 袁鐵.衛(wèi)星觀測到的全球閃電活動及其地域差異[J].地球物理學(xué)報(bào),2003, 46(6):743-750.

QIE X S, ZHOU Y J, YUAN T.Global lighting activities and their regional differences observed from the satellite[J].Chinese journal of geophysics, 2003, 46(6)：743-750.(in Chinese)

[4]GRAF K L, INAN U S, PIDDYACHIY D, et al.DEMETER observations of transmitter-induced precipitation of inner radiation belt electrons[J].Journal of geophysical research atmospheres,2009, 114(A7):155-164.

[5]MOLCHANOV O A, ROZHNOI A, SOLOVIEVA M, et al.Global diagnostics of the ionospheric perturbations related to the seismic activity using the VLF radio signals collected on the DEMETER satellite[J].Natural hazards, 2006, 6：745-753.

[6]HELLIWELL R A.Whistlers and related ionospheric phenomena[M].California：Stanford University Press, 1965.

[7]TAO X,BORTNIK J,FRIEDRICH M.Variance of transionospheric VLF wave power absorption[J].Journal of geophysical research：space physics (1978—2012), 2010, 115:A07303.

[8]HINES C O,PITTEWAY M L V,HINES C O.The reflection and ducting of atmospheric acoustic-gravity waves[J].Canadian journal of physics, 1965, 43 (12)：2222.

[9]ARANTES D S, SCARABUCCI R R.Full-wave analysis and coupling effects in a crossover region[J].Radio science, 1975, 10 (8/9)：801-81.

[10]WAIT J R.Ground-wave attenuation function for a Gaussian modulated carrier signal[J].Electronics letters, 1970, 6 (4)：106-107.

[11]YAGITANI S,NAGANO I,MIYAMURA K, et al.Full wave calculation of ELF/VLF propagation from a dipole source located in the lower ionosphere[J].Radio science, 1994, 29(1)：39-54.

[12]NYGREN T.A method of full wave analysis with improved stability[J].Planetary &space science, 1982, 30：427-430.

[13]NAGANO I,MAMBO M,HUTATSUISHI G.Numerical calculation of electromagnetic waves in an anisotropic multilayered medium[J].Radio science, 1975, 10 (6)：611-617.

[14]BUDDEN K G.The propagation of radio waves-the theory of radio waves of low power in the ionosphere and magnetosphere[M].Cambridge University Press, 1988.

[15]LEHTINEN N G, INAN U S.Radiation of ELF/VLF waves by harmonically varying currents into a stratified ionosphere with application to radiation by a modulated electrojet[J].Journal of geophysical research, 2008, 113(A6)：A06301.

[16]趙庶凡, 張學(xué)民, 趙正予, 等.NWC通信臺在電離層中激發(fā)電磁響應(yīng)的時(shí)變特征[J].地球物理學(xué)報(bào), 2015,58(7)：2263-2273.

ZHAO S F, ZHANG X M, ZHAO Z Y, et al.Temporal variations of electromagnetic responses in the ionosphere excited by the NWC communication station [J].Chinese journal of geophysics, 2015, 58(7)：2263-2273.(in Chinese)

[17]ZHAO S F, ZHANG X M, ZHAO Z Y, et al.The numerical simulation on ionospheric perturbations in electric field before large earthquakes[J].Annales geophysicae, 2014, 32(12):1487-1493.

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

PAN W Y.LF VLF ELF wave propagation [M].Chengdu：University of Electronic Science and Technology Press, 2004.(in Chinese)

[19]趙庶凡.地面VLF輻射滲透進(jìn)電離層和磁層傳播研究及其應(yīng)用[D].武漢：武漢大學(xué), 2015.

ZHAO S F.Study of the transionospheric propagation of terrestrial VLF radiations and up to the magnetosphere[D].Wuhan：Wuhan University, 2015.(in Chinese)

[20]BILITZA D.International reference ionosphere 2000[J].Radio science, 2001, 36(2):317.DOI:10.1029/2000RS002432.

[21]CUMMER S A.Modeling electromagnetic propagation in the Earth-ionosphere waveguide[J].IEEE transactions on antennas and propagation,2000, 48(9)：1420-1429.DOI：10.1109/8.898776.

[22]BARTON C E.International Geomagnetic Reference Field：the seventh generation[J].Journal of geomagnetism and geoelectricity,1997, 49：123-148.

[23]LEHTINEN N G, INAN U S.Full-wave modeling of transionospheric propagation of VLF waves[J].Geophysical research letters, 2009, 36 (L03104).

[24]OZAKI M, YAGITANI S, NAGANO I, et al.Ionospheric penetration characteristics of ELF waves radiated from a current source in the lithosphere related to seismic activity[J].Radio science, 2009, 44(1)：292-310.

[25]COHEN M B, INAN U.Terrestrial VLF transmitter injection into the magnetosphere[J].Journal of geophysical research：space physics (1978-2012), 2012, 117:A08310.

趙庶凡 (1985-),女,安徽人,中國地震局地震預(yù)測研究所助理研究員,博士,研究方向?yàn)殡姴▊鞑?

Full wave calculation of ground-based VLF radiation penetrating into the ionosphere

ZHAO Shufan1LIAO Li2ZHANG Xuemin1SHEN Xuhui3

(1.InstituteofEarthquakeScience,ChinaEarthquakeAdministration,Beijing100036,China;2.InstituteofGeophysics,ChinaEarthquakeAdministration,Beijing100081,China;3.InstituteofCrustalDynamics,ChinaEarthquakeAdministration,Beijing100085,China)

The spatial distribution of electromagnetic field excited by ground-based VLF transmitters with different radiation frequency and power under different geomagnetic parameters and ionospheric parameters are stimulated using the full-wave model constructed in this paper.Focusing on the absorption in the D/E region of the ionosphere, we can conclude that the attenuation of VLF radiation in the waveguide is only affected by the wave frequency in the near field which decreases with the increase of the wave frequency, while the D/E region absorption and total attenuation increase with wave frequency.The variation of the radiation power has no effect on the attenuation in the waveguide and the D/E region absorption.The D/E region absorption decreases with the increase of the geomagnetic field and the dip angle of the geomagnetic field.The D/E region absorption increases with ionospheric electron density and collision frequency.

ground-based VLF transmitter;full-wave solution;ionosphere absorption

趙庶凡, 廖力, 張學(xué)民,等.地面甚低頻輻射滲透進(jìn)電離層的數(shù)值模擬分析[J].電波科學(xué)學(xué)報(bào),2016,31(5):825-833.

10.13443/j.cjors.2016021701

ZHAO S F, LIAO L, ZHANG X M，et al.Full wave calculation of ground-based VLF radiation penetrating into the ionosphere[J].Chinese journal of radio science,2016,31(5):825-833.(in Chinese).DOI：10.13443/j.cjors.2016021701

2016-02-17

國家國際科技合作對俄科技合作專項(xiàng)(2014DFR21280)； 中國地震局地震預(yù)測研究所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2015IES0101)； 亞太空間合作項(xiàng)目以及ISSI-BJ項(xiàng)目聯(lián)合資助

P352

A

1005-0388(2016)05-0825-09

聯(lián)系人：趙庶凡 E-mail：zsf2008bj@126.com