趙海生 徐朝輝 高敬帆2) 許正文 吳健 馮杰 徐彬薛昆 李輝 馬征征

1)(中國電波傳播研究所,電波環(huán)境特性及模化技術(shù)重點實驗室,青島 266107)

2)(西安電子科技大學(xué)物理與光電工程學(xué)院,西安 710071)

1 引 言

1959年先驅(qū)者II發(fā)射時,Booker[1]首次利用垂測儀觀測到了電離層洞的存在,此后又多次觀測到火箭發(fā)射時產(chǎn)生的電離層洞,尤其在1973年Skylab發(fā)射期間觀測到了空前巨大的電離層空洞[2?4].Mendillo等[5,6]分析了電離層洞形成和快速消失的原因,認為是電離層中的氧離子(O+)與尾焰中氫氣(H2)和水(H2O)分子進行快速的離子-原子交換反應(yīng),生成分子離子,隨后迅速與電離層中的電子復(fù)合,造成電離層等離子體快速損耗.Mendillo等[7]深入研究了中性氣體釋放物在電離層中的動力學(xué)過程和釋放物在電離層近真空環(huán)境中的膨脹、冷凝、碰撞加熱和自由擴散等物理過程,建立了化學(xué)物質(zhì)釋放研究的理論基礎(chǔ).Anderson和Bernhardt[8]在考慮了真實的電離層化學(xué)反應(yīng)速度后,利用模式研究了在赤道電離層釋放H2后的電子損耗量化結(jié)果,先后建立了H2點源釋放的受擾電離層演化模型.H2和H2O等釋放物先與O+離子反應(yīng)后產(chǎn)生相應(yīng)的正離子,正離子再與電子進行復(fù)合反應(yīng)耗空電子,與H2和H2O吸附電子原理相同的釋放物還有二氧化碳(CO2)、氧氣(O2)、二氧化氮(NO2)等.

六氟化硫(SF6)氣體吸附電子的方式不同于H2和H2O,其直接與電子反應(yīng)產(chǎn)生相應(yīng)的負離子,因此SF6耗空電子的速度更快、效率更高.Paul和Bernhardt等對SF6氣體釋放開展了廣泛而深入的研究,先后研究了釋放物云團在不均勻大氣中的運動過程[9,10];研究了電離層擾動的等離子體和流體不穩(wěn)定性[11];研究了SF6在電離層中的化學(xué)反應(yīng)和動力學(xué)過程,建立點源釋放的仿真模型[12];研究了SF6釋放產(chǎn)生氣輝輻射增強的原理,提出了氣輝輻射強度的計算方法[13].趙海生等[14]建立了基于SF6釋放的面向試驗場景的三維精細化效應(yīng)仿真模型.

與中性氣體吸附電子產(chǎn)生電離層洞的原理不同,堿金屬、堿土金屬等金屬蒸氣在電離層中通過光致電離產(chǎn)生高密度電子云團,被稱為電子密度增強類化學(xué)物質(zhì).Hunton[15]研究了Ba的光致電離機制,以及Ba原子在電離層中的動力學(xué)過程.Koons和Rocdcr[16]從試驗現(xiàn)象出發(fā),研究了Ba電離層釋放激發(fā)電離層電場和等離子體波的物理機制.Schunk和Szuszczewicz[17,18]首先采用一維物理模型系統(tǒng)研究了不同電離層背景條件下,釋放電子密度增強類釋放物形成的人工等離子體云的擴散和經(jīng)典掃雪機效應(yīng).此后,隨著理論研究的深入和計算技術(shù)的發(fā)展,開展了大量三維等離子體云物理建模研究[19?30],建立了更為精細的等離子體云模型,研究了釋放物在電離層中的動力學(xué)過程,計算了釋放物密度的時空分布.模型從連續(xù)性方程和動量方程出發(fā),研究了背景中性風場力、慣性力、地磁場以及電離層電場對等離子體云運動過程的影響,仿真了等離子體云的演化過程.

電離層化學(xué)物質(zhì)釋放技術(shù)從生產(chǎn)至今已經(jīng)開展了大量的理論和試驗研究,也取得了豐碩的成果.但是,對電離層化學(xué)物質(zhì)釋放的早期試驗效應(yīng),特別是中性氣體釋放的早期試驗效應(yīng)的研究很少.在電離層釋放電子吸附類中性氣體能夠引起電離層電子密度耗空,在釋放之后快速形成電離層洞;同時,由于釋放中性氣體的快速膨脹,擠壓背景等離子體,在電離層洞的外邊緣產(chǎn)生殼狀電子密度增強結(jié)構(gòu),電離層洞和電子密度增強結(jié)構(gòu)同時存在是釋放早期試驗效應(yīng)的顯著特征.本文研究了電離層中性氣體釋放的早期試驗效應(yīng),建立了早期試驗效應(yīng)物理模型,仿真了釋放早期電子密度的時空演化過程,同時采用射線追蹤方法研究了釋放后10 s和120 s不同頻率信號在擾動區(qū)的傳播效應(yīng).電離層化學(xué)物質(zhì)釋放對深入研究電離層動力學(xué)過程及耦合機制、研究電離層不穩(wěn)定激發(fā)機制和過程、構(gòu)建新的不穩(wěn)定性理論體系具有重要的意義.

2 釋放物在電離層中的膨脹過程

釋放物的擴散過程滿足連續(xù)性方程和動量方程.

連續(xù)方程為

式中,n是釋放氣體密度,v是釋放氣體速度,P和L分別是化學(xué)產(chǎn)生和損失項.連續(xù)方程來源于在假設(shè)粒子流連續(xù)的情況下的玻爾茲曼方程.

釋放物速度可通過動量方程得到:

式中,g=?g?az是重力加速度,?az是向上速度的單位矢量,T是氣體溫度,ub是背景大氣風速,υ是釋放物總和碰撞頻率.由于釋放物粒子自由運動周期遠大于碰撞周期,加速度可以忽略,即dv/dt=0.那么,速度的表達式為

式中,D是擴散系數(shù),H=kT/(mg)是釋放氣體云團標高.氣體速度由兩部分組成:密度和溫度梯度產(chǎn)生的擴散速度,背景大氣運動產(chǎn)生的飄移速度.

將(3)式代入(1)式得到一般擴散方程:

如果D和T是常量,忽略重力、大氣運動、產(chǎn)生率和損失率,(4)式可歸納為簡化的擴散方程.

在背景大氣平面分層的假設(shè)下,忽略背景大氣的運動、釋放物分子的熱擴散、釋放物質(zhì)的膨脹過程,釋放物擴散方程可以表示為

式中,?是釋放物反應(yīng)損失系數(shù).選取笛卡爾坐標系,以釋放點為原點,z軸為垂直于地面的方向.假設(shè)大氣參量服從指數(shù)分布,有D(z)=D0exp(z/Hb),其中D0為釋放點(z=0)的擴散系數(shù),Hb是背景大氣標高.(5)式的最后一項表示注入氣體與電離層的化學(xué)反應(yīng)損耗.通過譜方法求解微分方程(5),可得到[10]

表達式(6)可用于計算自由擴散階段任意時刻中性釋放物密度分布.

3 早期釋放效應(yīng)數(shù)值模擬理論

在釋放之初,釋放物快速膨脹(約3000 m/s),膨脹速度大于背景大氣聲速(約600 m/s),在釋放物云團的邊界處形成了激波,波陣面附近背景氣體因壓縮而密度升高,電子密度也隨之升高,形成電子密度聚集區(qū).同時釋放物與背景等離子體化學(xué)反應(yīng),形成電子密度空洞,隨著釋放物的膨脹,電離層洞隨之擴大.之后釋放物擴散速度逐漸降低到背景聲速之下,激波消失,聚集區(qū)的電子密度逐漸恢復(fù)到背景水平;而形成的電離層洞繼續(xù)擴大,洞的深度逐漸減小.以下從等離子體連續(xù)性方程和動量方程出發(fā),推導(dǎo)電子密度增強區(qū)和電離層洞的形成及演化理論.

釋放區(qū)等離子體連續(xù)性方程可以表示為

等離子體動量方程可以表示為

式中,i=1代表背景電離層正離子O+,i=2代表釋放物云的負離子e為電子;ni,ui,mi,ei,pi分別代表密度、飄移速度、質(zhì)量、電量和壓強;E,B,j分別代表電場、磁場和總電流密度;為基于點源釋放的密度變化,為基于點源釋放的動量變化.

對于背景電離層O+,假設(shè)在釋放之初,其密度不隨時間變化對于釋放物負離子其密度變化可以表示為其中,σ為產(chǎn)生率,ns為釋放物密度.

對于O+離子,其動量變化可以表示為

由于地磁場對帶電離子的控制作用,不同帶電離子垂直于磁場方向的速度差很小,因此相關(guān)的碰撞項可以忽略不計,但是沿磁場方向的碰撞項不能忽略.忽略帶電離子的慣性,假設(shè)地磁場恒定,那么可以得到不同離子的漂移速度,表達式如下:

對于離子O+,a1⊥,a1‖,s1⊥,s1||分別定義為:

在(20)—(29)式中引進了電勢Φ,其定義為Φ=??(kTe/e)ln(ne),那么電場可以表示為

式中,E0為背景電場.假設(shè)在整個計算過程中E0保持恒定,可以得到以下關(guān)于Φ的方程:

式中,σ⊥為Pederson電導(dǎo)率,σ‖為平行電導(dǎo)率,Q為電荷總量.

通過求解(31)式能夠得到釋放區(qū)的電勢分布,從而求得釋放區(qū)電場分布.通過(16)式求得任意離子的飄移速度分布,將飄移速度代入連續(xù)性方程(7),可以求得任意時刻電子密度分布.

4 數(shù)值模擬

根據(jù)第3節(jié)理論推導(dǎo)結(jié)果,本文通過虛構(gòu)的試驗場景(試驗參數(shù)如表1所列),數(shù)值模擬了釋放早期,電離層洞的形成以及釋放物快速膨脹擠壓背景等離子體,引起背景電子密度升高的演化過程.根據(jù)釋放參數(shù)和SF6在電離層中的物理化學(xué)過程,通過求解釋放物擴散過程和等離子體擴散方程,建立SF6釋放的效應(yīng)仿真模型,模型算法設(shè)計流程如圖1所示.

效應(yīng)仿真流程按照以下5個步驟進行:

1)根據(jù)試驗具體情況進行釋放參數(shù)設(shè)置,包括釋放參數(shù)和火箭飛行參數(shù),同時,根據(jù)釋放物膨脹和化學(xué)反應(yīng)過程,計算得到初始時刻空間電荷分布;

2)將空間電荷分布代入等離子體泊松方程,求解得到釋放區(qū)域空間電場分布;

3)依據(jù)計算得到的空間電場分布、國際參考磁場模型、國際參考電離層模型IRI2012、標準大氣模型,代入等離子體動量方程,求解等離子體飄移速度;

4)根據(jù)等離子體飄移速度代入等離子體連續(xù)性方程求解得到等離子體密度分布(包括正負離子和電子);

5)將計算得到的等離子體密度分布代入步驟3),重復(fù)執(zhí)行步驟2)—4),計算得到任意時刻電子密度分布.

圖1 電離層中性氣體釋放的早期試驗效應(yīng)仿真流程圖Fig.1.The simulation fl owchart of the early time effects produced by neutral gas ionospheric chemical release.

為了驗證效應(yīng)仿真模型,本文采用虛擬試驗場景開展了數(shù)值模擬研究,仿真參數(shù)設(shè)置如表1所列,試驗時間201509301800LT代表當?shù)貢r間2015年9月30日18:00.

圖2描述了在300 km高度釋放40 kg中性氣體SF6,釋放后0—120 s釋放物密度分布的演化過程.在釋放之初的2 min內(nèi)釋放物密度逐漸降低,從釋放后1 s的最大密度2.3×109cm?3快速下降到釋放后120 s的3.4×105cm?3,釋放物分布范圍快速增加,釋放后5 s釋放物分布范圍為直徑約10 km的球體,釋放后120 s釋放物分布直徑達80 km以上.釋放物密度降低主要取決于釋放物快速膨脹體積增加,另外,釋放物與背景電子化學(xué)反應(yīng)產(chǎn)生相應(yīng)的負離子,也造成釋放物密度下降.

圖3描述了在300 km高度釋放40 kg中性氣體SF6,釋放后0—120 s擾動區(qū)域電子密度分布的演化過程.

從仿真結(jié)果可以看出,在釋放之初由于釋放物膨脹速度遠大于背景聲速,在釋放物邊界處形成了殼狀電子密度聚集區(qū),隨后膨脹速度逐漸降低,邊界處電子密度逐漸恢復(fù)到背景水平.釋放后1 s釋放物邊界處的最大電子密度達2.46×106cm?3,約為背景電子密度的2.8倍,隨后電子密度逐漸降低;釋放后30 s,邊界處的最大電子密度約1.58×106cm?3,約為背景電子密度的1.7倍.釋放后120 s,邊界處的最大電子密度約1.12×106cm?3,約為背景電子密度的1.2倍.整個電子密度聚集過程持續(xù)時間大約120 s,之后聚集區(qū)電子密度恢復(fù)到背景水平.在釋放后2 min內(nèi),電離層洞的尺度逐漸增加,釋放后5 s釋放物分布范圍為直徑約10 km的球體,釋放后120 s釋放物分布直徑達70 km以上,同時電離層洞的耗空深度緩慢降低,釋放后1 s電離層洞100%耗空,到釋放后120 s電離層洞耗空深度約為95%.

表1 主要仿真參數(shù)Table 1.Parameters for the simulation.

圖2 釋放物密度分布隨時間的變化Fig.2.The distribution of release density varies with time.

圖3 釋放早期釋放區(qū)電子密度演化過程Fig.3.The electron density evolution process in release zone in the early time.

基于仿真得到的釋放初期擾動效應(yīng)的演化過程,采用射線追蹤方法研究了電波信號在擾動區(qū)的傳播過程.射線追蹤是研究電離層電波傳播的有效方法之一,尤其對高頻以上頻段電波信號,射線追蹤方法能夠準確描述電波的傳播路徑.射線追蹤可以分為數(shù)值法射線追蹤和解析法射線追蹤兩種,解析法射線追蹤不能應(yīng)用于非均勻的電離層,其應(yīng)用受到局限;數(shù)值法射線追蹤可以應(yīng)用于非均勻電離層,本文采用二維射線追蹤方法研究電波信號在擾動區(qū)的傳播效應(yīng).由于射線追蹤算法十分成熟,在此不再贅述,詳細算法請查閱文獻[33,34].

圖4以釋放后10 s擾動區(qū)的電子密度分布為背景,研究了12—15 MHz電波信號在擾動區(qū)的傳播過程.釋放后10 s電離層洞的直徑約30 km,電子密度增強區(qū)厚度約10 km,最大電子密度值約1.82×106cm?3.

從圖4可以看出,11 MHz無線電波全部被電子密度增強區(qū)反射,12—14 MHz無線電波部分穿越電子密度增強區(qū),部分被電子密度增強區(qū)反射,穿越電子密度增強區(qū)的電波信號產(chǎn)生了顯著的聚焦效應(yīng),聚焦效應(yīng)隨電波頻率的減小而增強.聚焦效應(yīng)是由于空洞區(qū)域介電系數(shù)大于背景電離層介電系數(shù)產(chǎn)生的,根據(jù)電離層色散介質(zhì)的特性,頻率越小的電波介電系數(shù)增加越大,聚焦效應(yīng)越明顯.

圖5以釋放后120 s擾動區(qū)的電子密度分布為背景,研究了9—12 MHz電波信號在電子密度擾動區(qū)的傳播過程.釋放后120 s電離層洞的直徑約70 km,電子密度增強區(qū)厚度約50 km,最大電子密度約1.12×106cm?3.

圖4 短波信號在強反射區(qū)中的傳播Fig.4.The transmission of short wave signals in strong re fl ection areas.

圖5 短波信號在強反射區(qū)中的傳播Fig.5.The transmission of short wave signals in strong re fl ection areas.

從圖5可以看出,9 MHz無線電波在電離層洞內(nèi)經(jīng)過多次繞射,部分信號穿越電離層洞,部分信號反射回地面,而傳播到電子密度增強區(qū)的無線電波,被電子密度增強區(qū)反射回地面.10—11 MHz無線電波在電離層洞內(nèi)產(chǎn)生了聚焦效應(yīng),而后穿越電離層洞,而傳播到電子密度增強區(qū)的無線電波,部分信號穿越增強區(qū),另一部分信號反射回地面.12 MHz無線電波在電離層洞內(nèi)聚焦,在電子密度增強區(qū)散焦后全部穿越擾動區(qū).

在一次利用探空火箭開展的空間物理試驗中,位于試驗區(qū)的電離層垂直探測儀在火箭發(fā)射后觀測到了電子密度突然增強現(xiàn)象[35].該火箭為固體火箭,火箭噴焰的主要產(chǎn)物為H2O和CO2等,主要為電子密度耗空類化學(xué)物質(zhì).火箭發(fā)射前電離層臨界頻率為7.3 MHz,發(fā)射后154 s在F層高度出現(xiàn)一條強反射描跡,最高反射頻率超過12 MHz(如圖6(b)[35]所示).本次試驗觀測到強反射描跡極有可能是釋放初期,釋放氣體快速膨脹擠壓背景等離子體,引起邊界處等離子體密度升高產(chǎn)生的.

圖6 人工強反射區(qū)仿真電離圖和實測電離圖對比Fig.6.The comparison between the simulated ionogram and the measured ionogram in the arti fi cial strong re fl ection zone.

根據(jù)無線電波在電離層擾動區(qū)的傳播規(guī)律和本文虛擬的試驗場景,采用射線追蹤方法反演得到了釋放后30 s的垂測電離圖(如圖6(a)所示).反演得到的電離圖與實測電離圖吻合較好,初步驗證了本文建立的釋放早期的試驗效應(yīng)物理模型的正確性.由于缺少空間主動試驗的實測數(shù)據(jù),本文無法開展數(shù)值模擬結(jié)果與實測結(jié)果的詳細對比研究,隨著我國空間主動試驗的開展,該模型將進一步完善和改進.

5 總結(jié)與討論

在電離層高度釋放中性氣體能夠引起電離層電子密度耗空,在釋放之后快速形成電離層洞;同時,由于釋放中性氣體的快速膨脹,擠壓背景等離子體,在電離層洞的外邊緣產(chǎn)生“殼狀”電子密度增強結(jié)構(gòu),電離層洞和電子密度增強結(jié)構(gòu)同時存在是釋放早期試驗效應(yīng)的顯著特征.本文研究了電離層中性氣體釋放的早期試驗效應(yīng),建立了早期試驗效應(yīng)物理模型,仿真了釋放早期電子密度的時空演化過程.釋放后1 s增強區(qū)最大電子密度達2.46×106cm?3,約為背景電子密度的2.8倍,隨后邊界處的電子密度逐漸降低;釋放后30 s,最大電子密度約1.58×106cm?3,約為背景電子密度的1.7倍;釋放后120 s,增強區(qū)的最大電子密度約1.12×106cm?3,約為背景電子密度的1.2倍.在釋放后2 min內(nèi),電離層洞的尺度逐漸增加,釋放后5 s釋放物分布范圍為直徑約10 km的球體,釋放后120 s釋放物分布直徑達70 km以上,同時電離層洞的耗空深度緩慢降低,釋放后1 s電離層洞100%耗空,到釋放后120 s電離層洞耗空深度約為95%.

在300 km高度釋放40 kg SF6,以釋放后10 s和120 s的效應(yīng)仿真結(jié)果為研究場景,采用射線追蹤方法研究了不同頻率電波信號在擾動區(qū)的傳播效應(yīng).釋放后10 s電子密度增強效應(yīng)十分顯著,增強區(qū)厚度約10 km,最高可反射14 MHz的無線電波信號.釋放后120 s電子密度增強現(xiàn)象減弱,增強區(qū)厚度增加,最高可反射11 MHz的無線電波,9—12 MHz無線電波在擾動區(qū)發(fā)生了復(fù)雜的繞射、聚焦和散焦等傳播效應(yīng).根據(jù)電離層垂測儀工作原理和射線追蹤理論,反演得到了釋放后30 s的電離層垂測電離圖,電離圖上出現(xiàn)了明顯了附加增強描跡,最高反射頻率達11.6 MHZ.反演結(jié)果與一次火箭噴焰的實際觀測數(shù)據(jù)吻合較好,初步驗證了本模型的正確性.

電離層中性氣體釋放早期試驗效應(yīng)物理模型,重點研究釋放后2 min內(nèi)電離層擾動效應(yīng)的演化過程,是電離層中性氣體釋放物理模型的重要改進,更為準確地描述了釋放后2 min內(nèi)擾動區(qū)電子密度的時空演化過程.

[1]Booker H G 1961J.Geophys.Res.66 1073

[2]Mendillo M,Hawkins G S,Klobuchar J A 1975Science187 343

[3]Mendillo M,Hawkins G S,Klobuchar J A 1975J.Geophys.Res.80 2217

[4]Zinn J,Sutherland C D,Stone S N,Duncan L M,Behnke R 1982J.Atmos.Terr.Phys.44 1143

[5]Mendillo M,Jeffrey M,Forbes J M 1978J.Geophys.Res.83 151

[6]Mendillo M,Forbes J M 1982J.Geophys.Res.87 8273

[7]Mendillo M,Smith S,Coster A,Erickson P,Baumgardner J,Martinis C 2008SPACE WEATHER6 S09001

[8]Anderson D N,Bernhardt P A 1978J.Geophys.Res.83 4777

[9]Paul A,Bernhardt P A 1979J.Geophys.Res.84 4341

[10]Paul A,Bernhardt P A 1979J.Geophys.Res.84 793

[11]Bernhardt P A 1982J.Geophys.Res.87 7539

[12]Bernhardt P A 1984J.Geophys.Res.89 3929

[13]Bernhardt P A,Weber E J,Moore J G,Baumgardner J,Mendillo M J 1986Geophys.Res.Lett.91 8937

[14]Zhao H S,Xu Z W,Wu Z S,Feng J,Wu J,Xu B,Xu T,Hu Y L 2016Acta Phys.Sin.65 209401(in Chinese)[趙海生,許正文,吳振森,馮杰,吳健,徐彬,徐彤,胡艷莉2016物理學(xué)報65 209401]

[15]Hunton D E 1993Geophys.Res.Lett.20 563

[16]Koons H C,Rocdcr J L 1995J.Geophys.Res.100 5801

[17]SchunkR W,Szuszczewicz E P 1988J.Geophys.Res.93 12901

[18]Schunk R W,Szuszczewicz E P 1991J.Geophys.Res.96 1337

[19]Drake J F,Mulbrandon M,Huba J D 1988Phys.Fluids31 3412

[20]Zalesak S T,Drake J F,Huba J D 1988Radio Sci.23 591

[21]Zalesak S T,Drake J F,Huba J D 1990Geophys.Res.Lett.17 1597

[22]Ma T Z,Shunk R W 1991J.Geophys.Res.96 5793

[23]Ma T Z,Shunk R W 1993J.Geophys.Res.98 323

[24]Gatsonis N A,Hastings D E 1991J.Geophys.Res.96 7623

[25]Shuman N S,Hunton D E,Viggiano A A 2015J.Chem.Rev.115 4542

[26]Schunk R W,Szuszczewicz E P 1991J.Geophys.Res.96 1337

[27]Doles J H,Zabusky N J,Perkins F W 1976J.Geophys.Res.81 5987

[28]Zhao H S,Feng J,Xu Z W,Wu Z S 2016J.Geophys.Res.121 10508

[29]Zabushky N J,Doles J H,Perkins F W 1973J.Geophys.Res.78 711

[30]Lloyd K H,Haerendel G 1973J.Geophys.Res.78 7389

[31]Bernhardt P A 1984J.Geophys.Res.89 3929

[32]Bernhardt P A 1987J.Geophys.Res.92 4617

[33]Haselgrove J 1955Proc.Phys.Soc.London23 355

[34]John M Kelso 1968Radio Sci.3 1

[35]Zhang R F 1986Proceedings of the International Symposium on Space PhysicsBeijing,China,November 10–14,1986 pp5–100