張世田，陳林如，王元新

?

2009年7月22日日食期間VLF信號(hào)相位和場(chǎng)強(qiáng)的觀測(cè)和分析

張世田，陳林如，王元新

（中國電波傳播研究所青島分所，青島 266071）

2009年7月22日日食期間，分別在重慶、桂林和青島利用阿爾法場(chǎng)強(qiáng)與相位監(jiān)測(cè)儀觀測(cè)了來自俄羅斯阿爾法導(dǎo)航系統(tǒng)位于Novosibirsk、Krasnodar和Khabarovsk 的3個(gè)發(fā)射臺(tái)3個(gè)頻率的甚低頻信號(hào)，并對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了分析。通過建立日食期間的電離層模型，結(jié)合實(shí)際傳播路徑地面電磁參數(shù)，推測(cè)全食區(qū)低電離層的等效高度；給出經(jīng)不同路徑傳播的VLF（甚低頻）信號(hào)相位突變與日食發(fā)生時(shí)刻的關(guān)系及原因；分析了日食期間經(jīng)不同路徑傳播的VLF信號(hào)的場(chǎng)強(qiáng)起伏變化，發(fā)現(xiàn)場(chǎng)強(qiáng)變化滯后于相位變化這一現(xiàn)象。

甚低頻(VLF)異常傳播；日食；傳播相位；阿爾法導(dǎo)航系統(tǒng)

日食期間是進(jìn)行某些科學(xué)觀測(cè)的絕好時(shí)機(jī)?；谌帐称陂g的觀測(cè)，可以進(jìn)行日地關(guān)系方面的課題研究，包括太陽對(duì)電離層的控制作用的研究。2009年7月22日，在我國西起西藏、東至浙江的廣泛區(qū)域出現(xiàn)了日全食。此次日全食在我國覆蓋范圍廣、可觀測(cè)時(shí)間長，除了全食帶，我國很多地區(qū)處于偏食區(qū)。國內(nèi)很多單位在日食期間組織安排了相關(guān)觀測(cè)。

為了深入研究日食對(duì)電離層的影響，中國電波傳播研究所對(duì)此次日食進(jìn)行了較大規(guī)模的電離層加密觀測(cè)。作者課題組承擔(dān)了利用阿爾法場(chǎng)強(qiáng)與相位監(jiān)測(cè)儀分別在重慶（106.26oN，29.68oE），桂林（109.94oN，25.39oE）及青島（120.31oN，36.13oE）監(jiān)測(cè)來自俄羅斯阿爾法導(dǎo)航系統(tǒng)位于新西伯利亞發(fā)射臺(tái)（Novosibirsk，俗稱主臺(tái)，84.43oE，55.75oN），克拉斯諾達(dá)爾發(fā)射臺(tái)（Krasnodar，俗稱西副臺(tái)，38.16oE，45.41oN）和哈巴羅夫斯克發(fā)射臺(tái)（Khabarovsk，俗稱東副臺(tái)，136.61oE，50.07oN）這3個(gè)臺(tái)的3個(gè)頻率的VLF（甚低頻）導(dǎo)航信號(hào)在日食期間相位和場(chǎng)強(qiáng)變化的任務(wù)，獲得了完整的實(shí)測(cè)數(shù)據(jù)。據(jù)此數(shù)據(jù)分析了日食期間我國上空低電離層的等效高度的變化；研究了日食發(fā)生與VLF信號(hào)相位變化的關(guān)系及其原因；還分析了日食期間VLF信號(hào)場(chǎng)強(qiáng)起伏量及場(chǎng)強(qiáng)變化與日食的關(guān)系。

1 信號(hào)傳播路徑及導(dǎo)航信號(hào)發(fā)射順序

2009年7月22日我國見食區(qū)、發(fā)射臺(tái)、觀測(cè)站分布和信號(hào)傳播路徑等情況如圖1所示。

圖1 日食帶、發(fā)射臺(tái)、觀測(cè)站分布和信號(hào)傳播路徑

日食剛發(fā)生時(shí)，太陽光照強(qiáng)度減弱，當(dāng)?shù)厣峡盏碗婋x層電子密度降低，低電離層等效高度升高。當(dāng)日全食發(fā)生時(shí)，光照被完全遮蔽，低電離層等效高度上升到最大高度，生光后低電離層等效高度連續(xù)下降，直至回到正常高度為止。日食期間信號(hào)傳播路徑和低電離層等效高度示意于圖2。

為了全面地研究日食期間傳播的甚低頻（VLF）信號(hào)相位和場(chǎng)強(qiáng)的變化，選取重慶、桂林和青島3地作為觀測(cè)點(diǎn)，每個(gè)觀測(cè)點(diǎn)可以接收3條電波傳播路徑上各3個(gè)頻率的甚低頻信號(hào)。發(fā)射臺(tái)到重慶觀測(cè)點(diǎn)的電波傳播路徑穿過偏食區(qū)止于全食帶；發(fā)射臺(tái)到桂林觀測(cè)點(diǎn)的電波傳播路徑穿過偏食區(qū)及全食帶后止于偏食區(qū)；發(fā)射臺(tái)到青島觀測(cè)點(diǎn)的電波傳播路徑僅位于偏食區(qū)內(nèi)。

實(shí)測(cè)所用的甚低頻信號(hào)場(chǎng)強(qiáng)和相位接收設(shè)備為“CJA-1阿爾法場(chǎng)強(qiáng)與相位監(jiān)測(cè)儀”，該設(shè)備可同時(shí)接收3個(gè)發(fā)射臺(tái)發(fā)射的3個(gè)頻率的電波信號(hào)，具有較高的相位測(cè)量和場(chǎng)強(qiáng)測(cè)量的精度。

圖2 日食期間信號(hào)傳播路徑和低電離層等效高度示意圖

3個(gè)發(fā)射臺(tái)（主臺(tái)、東副臺(tái)、西副臺(tái)）的信號(hào)發(fā)射順序（發(fā)射周期3.6 s）如圖3所示。

圖3 信號(hào)發(fā)射順序

注：1=11.905 kHz，2=12.649 kHz，3=14.881 kHz

在每個(gè)觀測(cè)點(diǎn)利用阿爾法場(chǎng)強(qiáng)與相位監(jiān)測(cè)儀可以同時(shí)觀測(cè)東副臺(tái)、主臺(tái)和西副臺(tái)發(fā)射的3個(gè)頻率信號(hào)的相位和場(chǎng)強(qiáng)。

2 VLF電波傳播原理

當(dāng)不考慮沿路徑上可能發(fā)生的“模轉(zhuǎn)換”影響時(shí)[1-2]，VLF電波傳播路徑上的電場(chǎng)強(qiáng)度分量[3]可以表示為

式（2）中，

由于低電離層等效高度的變化，場(chǎng)強(qiáng)將產(chǎn)生起伏，場(chǎng)強(qiáng)起伏標(biāo)準(zhǔn)偏差估計(jì)[4]公式

式（4）中

3 實(shí)驗(yàn)數(shù)據(jù)分析

3.1 主臺(tái)到觀測(cè)點(diǎn)的相位變化

各觀測(cè)點(diǎn)及主臺(tái)的見食時(shí)刻如表1所示，表中時(shí)間為北京時(shí)間。主臺(tái)至各觀測(cè)點(diǎn)的傳播路徑的初虧時(shí)刻就是各觀測(cè)點(diǎn)初虧時(shí)刻，日食復(fù)圓時(shí)，日食陰影已完全掃過整條傳播路徑。圖4，圖5和圖6分別為在重慶、桂林和青島3地接收到的主臺(tái)信號(hào)的相位變化曲線，其相位初虧時(shí)刻、相位復(fù)原時(shí)刻、最大相位突變時(shí)刻和相位最大突變量總結(jié)于表2中。

表1 觀測(cè)點(diǎn)及主臺(tái)見食時(shí)刻 （hh:mm）

圖5 主臺(tái)—桂林傳播信號(hào)的相位（信號(hào)頻率12.649 kHz）

圖6 主臺(tái)—青島傳播信號(hào)的相位（信號(hào)頻率12.649 kHz）

表2 日食期間主臺(tái)信號(hào)相位變化情況

對(duì)比表1和表2，可以看出（以觀測(cè)點(diǎn)重慶為例）：相位初虧時(shí)刻08:30遲于光學(xué)初虧時(shí)刻08:07，但早于光學(xué)食甚時(shí)刻09:15；相位復(fù)原時(shí)刻10:20早于光學(xué)復(fù)圓時(shí)刻10:32；最大相位突變時(shí)刻09:22出現(xiàn)在光學(xué)食甚（或生光）之后。在桂林和青島的觀測(cè)也得出了類似的結(jié)果。

相位初虧遲于觀測(cè)點(diǎn)光學(xué)初虧，這主要是：當(dāng)各觀測(cè)點(diǎn)發(fā)生光學(xué)初虧時(shí)，電離層變化比較緩慢，對(duì)相位的影響比較小，隨著時(shí)間的推移，電離層變化加劇，等效高度明顯升高，此時(shí)，傳播相位發(fā)生明顯變化，出現(xiàn)相位初虧。

相位復(fù)原早于觀測(cè)點(diǎn)光學(xué)復(fù)圓時(shí)刻，這主要是各觀測(cè)點(diǎn)光學(xué)復(fù)圓時(shí)，日食陰影已完全掃過整條傳播路徑。在接近各觀測(cè)點(diǎn)光學(xué)復(fù)圓時(shí)刻時(shí)，電離層對(duì)傳播相位的影響漸漸恢復(fù)正常，傳播相位基本恢復(fù)正常，相位復(fù)原。

3.2 東副臺(tái)到觀測(cè)點(diǎn)的相位變化

各觀測(cè)點(diǎn)及東副臺(tái)的見食時(shí)刻如表3所示。就東副臺(tái)而言，各觀測(cè)點(diǎn)初虧時(shí)刻即為東副臺(tái)至各觀測(cè)點(diǎn)的傳播路徑的初虧時(shí)刻，但觀測(cè)點(diǎn)日食復(fù)圓時(shí)，日食陰影并沒有完全掃過整條傳播路徑，東副臺(tái)仍位于日食陰影中。圖7和圖8分別是在重慶和桂林接收到的東副臺(tái)信號(hào)的相位變化曲線，其相位初虧時(shí)刻、相位復(fù)原時(shí)刻、最大相位突變時(shí)刻和相位最大突變量總結(jié)于表4。

表3 觀測(cè)點(diǎn)及東副臺(tái)見食表 （hh:mm）

圖8 東副臺(tái)—桂林傳播信號(hào)的相位（信號(hào)頻率11.905 kHz）

表4 日食期間東副臺(tái)相位變化情況

對(duì)比表3、表4可以看出（以重慶為例，其他觀測(cè)點(diǎn)得到類似的觀測(cè)結(jié)果）：相位初虧時(shí)刻08:40遲于光學(xué)初虧時(shí)刻08:07，相位復(fù)原時(shí)刻10:37早于日食陰影掃過東副臺(tái)時(shí)刻10:53，最大相位突變時(shí)刻09:37出現(xiàn)在光學(xué)食甚（或生光）時(shí)刻09:15之后，這與主臺(tái)到各觀測(cè)點(diǎn)相位突變現(xiàn)象一致，其原因亦相同。

3.3 西副臺(tái)到觀測(cè)點(diǎn)的相位變化

西副臺(tái)到各觀測(cè)點(diǎn)的傳播路徑較長（大于6 060 km），且大部分的路徑都在日食區(qū)內(nèi)。各觀測(cè)點(diǎn)初虧時(shí)刻基本上亦是各傳播路徑的初虧時(shí)刻，日食復(fù)圓時(shí)，日食陰影已完全掃過傳播路徑。圖9，圖10和圖11分別為在重慶、桂林及青島接收到的西副臺(tái)信號(hào)的相位變化曲線，其相位初虧時(shí)刻，相位復(fù)原時(shí)刻，最大相位突變時(shí)刻和相位最大突變量總結(jié)于表5。

圖9 西副臺(tái)—重慶傳播信號(hào)的相位（信號(hào)頻率12.649 kHz）

圖10 西副臺(tái)—桂林傳播信號(hào)的相位（信號(hào)頻率12.649 kHz）

圖11 西副臺(tái)—青島傳播信號(hào)的相位（信號(hào)頻率12.649 kHz）

表5 日食期間西副臺(tái)信號(hào)相位變化情況

由圖9至圖11可以看出，日食期間相位發(fā)生變化，但由于傳播路徑較長，各階波模中主要是基模起作用，故相位只是緩變。

對(duì)比表1、表5可以看出（以重慶為例，其他觀測(cè)點(diǎn)得到類似結(jié)果）：相位初虧時(shí)刻06:54遠(yuǎn)早于光學(xué)初虧時(shí)刻08:07，提前2 h左右；相位復(fù)原時(shí)刻12:00遠(yuǎn)遲于光學(xué)復(fù)圓時(shí)刻10:32，滯后1.5 h左右，這一現(xiàn)象的原因有待進(jìn)一步分析研究。

3.4 主臺(tái)到觀測(cè)點(diǎn)場(chǎng)強(qiáng)的變化

日食期間場(chǎng)強(qiáng)的起伏變化及場(chǎng)強(qiáng)變化時(shí)刻與日食時(shí)刻的關(guān)系在理論上很難給出，以前的日食觀測(cè)也沒有進(jìn)行這方面的研究[5-6]。為了研究日食期間場(chǎng)強(qiáng)的變化，對(duì)其進(jìn)行了觀測(cè)，日食期間在桂林和青島接收到的主臺(tái)信號(hào)的場(chǎng)強(qiáng)變化情況分別如圖12和圖13所示。

圖12 主臺(tái)—桂林傳播信號(hào)的場(chǎng)強(qiáng)（信號(hào)頻率14.881 kHz）

圖13 主臺(tái)—青島傳播信號(hào)的場(chǎng)強(qiáng)（信號(hào)頻率12.649 kHz）

由圖12可以看出：在初虧時(shí)刻08:10之后，場(chǎng)強(qiáng)幅度逐漸增大；在復(fù)圓時(shí)刻10:36前后，場(chǎng)強(qiáng)幅度出現(xiàn)“快速日出”現(xiàn)象（場(chǎng)強(qiáng)迅速減小，繼而增大）；場(chǎng)強(qiáng)在16:39迅速減小，出現(xiàn)“快速日落”現(xiàn)象（場(chǎng)強(qiáng)急劇減小到正常值）。

由圖13可以看出：在初虧時(shí)刻08:22場(chǎng)強(qiáng)幅度沒有明顯變化，至10:51場(chǎng)強(qiáng)幅度有所增加；在復(fù)圓時(shí)刻10:57前后，場(chǎng)強(qiáng)幅度出現(xiàn)“快速日出”現(xiàn)象；場(chǎng)強(qiáng)在16:51迅速減小，出現(xiàn)“快速日落”現(xiàn)象。

日食期間各觀測(cè)點(diǎn)接收到的東副臺(tái)和西副臺(tái)的信號(hào)場(chǎng)強(qiáng)變化不明顯，基本上無變化。

3.5 日食對(duì)傳播相位的影響分析

當(dāng)傳播路徑比較短時(shí)，由于接收點(diǎn)與發(fā)射臺(tái)的距離較近，除了基模以外還存在其他波模，因而日食期間會(huì)有明顯的相位突變發(fā)生，但是相位突變與日食發(fā)生時(shí)刻并不一致。

當(dāng)傳播路徑較長時(shí)，由于各階波模中主要是基模起作用，所以相位變化較緩不發(fā)生突變。相位變化的時(shí)刻與日食發(fā)生的時(shí)刻相差比較大，這主要是因?yàn)檎麠l傳播路徑低電離層等效高度發(fā)生變化需要較長時(shí)間所致。

為切實(shí)合理推測(cè)日食期間低電離層等效高度，本文采取了以下幾項(xiàng)措施：

1）結(jié)合中國電波傳播研究所地面導(dǎo)電率數(shù)據(jù)庫[7]計(jì)算相速、激勵(lì)因子相位等傳播參數(shù)；

2）考慮激勵(lì)因子相位的影響；

表6 日食期間低電離層等效高度和相位變化量的推算結(jié)果

3.6 日食對(duì)場(chǎng)強(qiáng)的影響分析

由圖12和圖13可以看出：日食期間場(chǎng)強(qiáng)變化不是很明顯，但是恰在日食復(fù)圓時(shí)刻會(huì)出現(xiàn)“快速日出”現(xiàn)象，繼而場(chǎng)強(qiáng)迅速增加，此時(shí)場(chǎng)強(qiáng)起伏標(biāo)準(zhǔn)偏差可以達(dá)到30~40 dB，而正常情況下，由（4）式推算可知場(chǎng)強(qiáng)起伏標(biāo)準(zhǔn)偏差為十幾個(gè)dB，日食期間場(chǎng)強(qiáng)起伏標(biāo)準(zhǔn)偏差可達(dá)30~40 dB的現(xiàn)象將持續(xù)6h左右，此后恢復(fù)正常。通過對(duì)比圖6和圖13也可發(fā)現(xiàn)：電離層復(fù)原對(duì)傳播相位的影響與對(duì)場(chǎng)強(qiáng)的影響并不一致，對(duì)場(chǎng)強(qiáng)的影響往往要晚一些、持續(xù)的時(shí)間更長一些，造成這一現(xiàn)象的原因不能簡單地從低電離層等效高度變化方面解釋，有待于進(jìn)一步分析研究。

4 結(jié)語

由以上實(shí)驗(yàn)數(shù)據(jù)和分析得出以下結(jié)論：

1）日食對(duì)VLF傳播相位的影響程度與傳播路徑有關(guān)：當(dāng)傳播路徑比較短時(shí)，相位將發(fā)生突變；當(dāng)傳播路徑較長時(shí)，相位緩慢變化不發(fā)生突變。相位的變化與日食帶和傳播路徑交點(diǎn)處的初虧、食甚、復(fù)圓并不同步，路徑越長現(xiàn)象越明顯。這次實(shí)驗(yàn)中，日食期間相位變化量在1.52~7.12μs之間。

2）日食期間低電離層等效高度會(huì)明顯升高，但是并不會(huì)升高到夜晚低電離層的等效高度，全食帶低電離層等效高度為78 km左右，較白天低電離層等效高度（70 km左右）升高8 km左右。

3）當(dāng)傳播路徑較短時(shí)，日食對(duì)甚低頻信號(hào)場(chǎng)強(qiáng)影響明顯，有30~40dB的起伏，但場(chǎng)強(qiáng)明顯增大的時(shí)刻遲于日食復(fù)圓時(shí)刻，其變化與相位變化亦不同步，造成這一現(xiàn)象的原因有待進(jìn)一步研究。

[1] WATT A D. VLF Radio Engineering[M]. New York: Pergamon Press in Oxford, 1967.

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

[3] 田育庶, 潘威炎, 余遠(yuǎn)超. 甚低頻(VLF)無線電系統(tǒng)電波場(chǎng)強(qiáng)與相位計(jì)算方法[S].

[4] 張世田, 陳恩平, 王元新, 等. 甚低頻在多模干涉區(qū)的場(chǎng)強(qiáng)起伏標(biāo)準(zhǔn)偏差估計(jì)[J] . 電波科學(xué)學(xué)報(bào), 2011, 26(2): 286-290.

[5] 趙協(xié)中.日食期間的甚低頻傳播效應(yīng)[J]. 電波科學(xué)學(xué)報(bào), 1990, 1(5): 23-31.

[6] 劉萬通. VLF相位突然異常與太陽x射線事件相關(guān)性的分析[J]. 空間科學(xué)學(xué)報(bào), 1987, 3(7): 185-189.

[7] ITU-R.ITU-R Rec. 832. World atlas of ground conductivities[S].USA: National Telecommunications and International Administration, 1982.

Observation and analysis of VLF phase/amplitudeduring a solar eclipse of July 22, 2009

ZHANG Shi-tian, CHEN Lin-ru, WANG Yuan-xin

(Qingdao Branch, Chinese Research Institute of Radiowave Propagation, Qingdao 266071, China)

The phase and field strength of VLF signals at three frequencies, which were transmitted from Novosibirsk, Krasnodar and Khabarovs stations of the Russian alpha hyperbolic navigational system, were measured with the alpha phase/field-strength monitor in Chongqing, Guilin and Qingdao respectively during a solar eclipse of July 22, 2009, and the measurements were analyzed. The ionospheric equivalent height in the total eclipse regions was derived by establishing the ionospheric model and considering the electromagnetic parameters for the propagation path. The relation between the sudden change in phase of VLF signal and the time when the solar eclipse took place were obtained for different propagation paths, and the cause of the relation were analyzed. Besides, the amplitude changes occurred during the solar eclipse for different propagation paths were analyzed and it was found that the phase change occurred earlier than the amplitude change.

VLF(very low frequency) abnormal propagation; eclipse; phase of propagation; Alpha navigation system

TN011

A

1674-0637(2011)01-0066-11

2010-11-13

國家重點(diǎn)實(shí)驗(yàn)室基金資助項(xiàng)目（9140C080405100C08）

張世田，男，碩士，高級(jí)工程師，主要從事電波傳播、電磁散射與逆散射等研究。