戴順，羅近濤，劉海文，李實鋒，李祎豐，李浩，劉長虹，劉偉，劉永朋

2020年6月21日日環(huán)食低頻時碼信號監(jiān)測

戴順1,2，羅近濤2，劉海文3，李實鋒2，李祎豐2，李浩2，劉長虹2，劉偉2，劉永朋2

（1. 華東交通大學(xué)，南昌 330013；2. 中國科學(xué)院 國家授時中心，西安 710600；3. 西安交通大學(xué)，西安 710049）

2020年6月21日日環(huán)食前后，在西安臨潼及商洛洛南搭建觀測系統(tǒng)，對商丘B(yǎng)PC發(fā)播臺68.5kHz信號進行了觀測。將采集數(shù)據(jù)可視化處理，對比分析多日信號波動情況，結(jié)果表明日食期間低頻時碼信號幅度發(fā)生波動，并在食甚后達到極值。根據(jù)電離層變化分析信號特點，從日出—日落效應(yīng)層面闡述波動原因，證實低頻時碼信號傳播受日食影響，為研究日食期間低頻時碼信號變化特點提供參照。

日環(huán)食；低頻時碼；電離層

0 引言

隨著現(xiàn)代通信、導(dǎo)航、航天等領(lǐng)域的發(fā)展，高精度授時技術(shù)的重要性日益凸顯。低頻時碼授時是在長波授時基礎(chǔ)上發(fā)展起來的一種授時方法，是國際電信聯(lián)盟（ITU）建議和鼓勵采用的授時方式[1]。該方法采用時間編碼方式，通過無線電發(fā)射裝置，將高精度原子鐘產(chǎn)生和保持的時間信息以無線電波為載體進行發(fā)播，因其覆蓋范圍廣、傳播損耗小、接收簡單，得到廣泛應(yīng)用[2]。

低頻時碼的傳播主要有地波和天波兩種方式。地波沿地球表面以繞射的方式到達接收點，天波經(jīng)電離層反射傳播至接收點[3]。電離層在垂直方向呈現(xiàn)出分層結(jié)構(gòu)，一般分為D層、E層和F層，其中，F(xiàn)層又分為F1層和F2層。夜間D層、E層和F1層基本消失，只剩下F2層仍可用于無線電傳播[4]。這種電離層的晝夜的變化，使低頻時碼信號在夜間強于白晝[5]，呈現(xiàn)出“晝?nèi)跻箯姟钡奶攸c。太陽活動會影響電離層，進而對低頻時碼信號傳播產(chǎn)生影響。

中國科學(xué)院國家授時中心于2007年在河南省商丘市虞城縣新建設(shè)了一座68.5 kHz低頻時碼發(fā)播臺，呼號為BPC[6]。其發(fā)射功率為100 kW，信號覆蓋半徑約為2 000 km[7]，日發(fā)播時長21 h[8]，北京時間早上5點至8點停機檢修維護。

2020年6月21日發(fā)生了一次日環(huán)食[9]，我國大部分地區(qū)可見偏食。日食發(fā)生時，月亮處在太陽和地球中間，遮擋住了部分射向地球的太陽光，在日食區(qū)會出現(xiàn)短暫的晝夜更替現(xiàn)象。

1 2020年6月21日日環(huán)食及觀測

2020年6月21日日環(huán)食中心食區(qū)橫跨非洲、亞洲、太平洋，如圖1（a）所示①http://www.heavens-above.com/SolarEclipselnteractive.aspx?jdmax=2459021.77865012，圖中圓圈標(biāo)記處為本次日環(huán)食最大食分地。在我國大陸，環(huán)食帶由西藏中部進入，從福建南部離開，主要覆蓋貴州東北部、湖南中部、江西南部等地，其他地區(qū)可見不同程度的偏食，如圖1（b）所示。我國部分城市的日食時刻信息如表1所示[10]。

表1 日環(huán)食時刻及食甚時食分

注：本文出現(xiàn)時間皆為北京時間（UTC+8h）。

本次監(jiān)測分別在西安臨潼及商洛洛南設(shè)置觀測點，對商丘B(yǎng)PC臺站發(fā)播的68.5 kHz低頻時碼信號進行監(jiān)測，監(jiān)測系統(tǒng)如圖2所示。

圖2 觀測系統(tǒng)

高增益環(huán)形天線接收到的低頻信號，通過饋線傳輸至射頻處理器件，進行下變頻濾波等處理，處理所得的信號通過SDRuno頻譜存儲軟件進行采集存儲，本文所用數(shù)據(jù)采集頻次為1s一次。數(shù)據(jù)采集系統(tǒng)通過網(wǎng)絡(luò)時間服務(wù)器進行時間同步，保障兩地觀測時間的一致性。

2 數(shù)據(jù)處理

臨潼觀測點于6月8日至7月5日期間對68.5 kHz信號進行監(jiān)測，洛南觀測點從6月12日至6月21日期間對68.5 kHz信號進行監(jiān)測，在7月6日至7月19日期間進行后續(xù)監(jiān)測。對獲得數(shù)據(jù)的處理過程如圖3所示。

圖3 采集數(shù)據(jù)后處理過程

采集數(shù)據(jù)的整體有效性受外部不穩(wěn)定因素影響，如遠程連接觀測設(shè)備導(dǎo)致采集中斷等。以保證數(shù)據(jù)可靠性為前提，在對現(xiàn)有數(shù)據(jù)進行進一步處理前，應(yīng)先篩出有效數(shù)據(jù)。有效數(shù)據(jù)應(yīng)符合以下幾個條件：

① 連續(xù)每秒記錄數(shù)據(jù)，數(shù)據(jù)丟失量不超過整體數(shù)據(jù)的3%；

② 連續(xù)時間內(nèi)，相同數(shù)據(jù)值不超過5條；

③ 在以天為單位的連續(xù)數(shù)據(jù)中，無數(shù)據(jù)中斷現(xiàn)象。

在采集數(shù)據(jù)時，系統(tǒng)采集響應(yīng)占用一定量的時間，造成約42 s采集41s數(shù)據(jù)的情況；且低頻時碼信號會受到干擾和噪聲的影響，出現(xiàn)信號突變現(xiàn)象。對數(shù)據(jù)進行平滑處理可提高信息量，減少誤差。數(shù)據(jù)平滑可采用多種方法，常用的方法包括自適應(yīng)平滑、簡易平滑、移動平均平滑[11]等。其中，移動平均平滑在修勻數(shù)據(jù)點的同時，能夠較為準(zhǔn)確地反映出變化趨勢。本次數(shù)據(jù)平滑處理采用移動平均平滑法，平滑參數(shù)根據(jù)經(jīng)驗值選為300。

采集數(shù)據(jù)具有晝夜強度差異，為對比1 d中不同時段的信號強弱，將數(shù)據(jù)點以二維圖的形式進行可視化處理；采集數(shù)據(jù)在不同日期同一時刻的強度具有相似性，為對比不同日期同一時間段的信號強度，根據(jù)日期將數(shù)據(jù)分成00:00—24:00的子數(shù)據(jù)，并用自主編寫的python程序?qū)⒍鄠€子數(shù)據(jù)集成到同一個二維圖中，以不同線條形式進行區(qū)分，以便事后分析。

根據(jù)觀測站的坐標(biāo)位置及自主編寫的python程序，計算得出觀測站日環(huán)食時刻信息如表2所示，獲得時刻值與日食計算器②http://interesting-sky.china-vo.org/astronomical-database/rishijisuanqi.html相近。將采集的數(shù)據(jù)進行處理，信號變化曲線如圖4和圖5所示。圖中實線為日環(huán)食當(dāng)天數(shù)據(jù)曲線，標(biāo)志初虧、復(fù)圓時刻為當(dāng)?shù)厝窄h(huán)食的變化時刻。

表2 觀測點日環(huán)食時刻表

圖4 臨潼觀測68.5 kHz信號變化

3 觀測結(jié)果分析

通過24 h內(nèi)信號的變化曲線可以看出，白天時段（08:00—18:00）信號波動幅度較為平穩(wěn)，夜間時段（20:00—24:00，00:00—5:00）及晝夜過渡段（18:00—20:00）信號波動幅度較大。這是因為信號經(jīng)電離層反射，白天時段內(nèi)電離層相對穩(wěn)定，夜間時段及晝夜過渡時段內(nèi)的電離層處于變化狀態(tài)。其中， 05:00—08:00為BPC發(fā)播臺停機檢修時間。兩地監(jiān)測68.5 kHz信號皆呈現(xiàn)“晝?nèi)跻箯姟钡恼w趨勢，符合低頻時碼傳播隨電離層不同發(fā)生的晝夜變化，且在日環(huán)食期間信號發(fā)生異常時刻相近、時長相似。

3.1 臨潼觀測站

在臨潼觀測日環(huán)食發(fā)生時段內(nèi)，68.5 kHz信號存在劇烈變化，該觀測站日食前后約一個月的其他有效數(shù)據(jù)中，同一時段內(nèi)信號平穩(wěn)無異常，如圖4及圖6所示。臨潼觀測點接收的68.5 kHz信號在日食期間的異常，表現(xiàn)為在15:18附近變?nèi)?，?6:23附近恢復(fù)；在15:59和16:05附近達到極值，極值時刻稍晚于當(dāng)?shù)厥成鯐r刻。在日環(huán)食過程中，太陽輻射各分量均有明顯的變化[12]，其有效輻射主要由地面低頻時碼輻射和大氣逆輻射所決定，且日食過程中入射到地面的太陽輻射量有一個作用過程。因此低頻時碼輻射在日食過程中的減弱較小，幅度的最大變化量與食甚時的食分值基本上成正相關(guān)，且最大減弱時刻并非食甚時刻，而是稍有滯后[12-14]。

圖6 臨潼觀測日環(huán)食期間68.5 kHz信號變化

3.2 洛南觀測站

在洛南觀測日環(huán)食發(fā)生時段內(nèi)，68.5 kHz信號同樣存在劇烈變化，該觀測站日食后約兩周的部分其他有效數(shù)據(jù)中，同一時段內(nèi)信號有極微弱的變化趨勢，但幅度均遠不及日食當(dāng)天，如圖5及圖7所示。洛南觀測點接收的68.5 kHz信號在日食期間的異常，表現(xiàn)為在15:32附近變強，在16:41附近恢復(fù)；在16:08附近達到極值，極值時刻稍晚于當(dāng)?shù)厥成鯐r刻。洛南觀測站地處秦嶺山脈中，海拔約為1 km。在不吸收電磁能量的真空中，電波的傳輸損耗只與頻率和距離有關(guān)，但實際上電波傳播的實際環(huán)境總是涉及到復(fù)雜的傳輸媒介，傳輸媒介的不均勻、多徑效應(yīng)、地貌特征等都會對信號產(chǎn)生影響。

圖7 洛南觀測日環(huán)食期間68.5 kHz信號變化

4 結(jié)語

2020年6月21日日環(huán)食前后，在西安臨潼和商洛洛南兩地對商丘B(yǎng)PC發(fā)播臺68.5 kHz信號進行觀測。對觀測數(shù)據(jù)分析表明，日食這一天象對低頻時碼有明顯的影響，作用機制為影響電離層進而影響低頻時碼信號傳輸。電離層隨經(jīng)緯度的不同呈現(xiàn)出復(fù)雜的空間變化，并且伴隨晝夜、季節(jié)、太陽黑子活動等現(xiàn)象發(fā)生不同的變化。因此，要進一步了解低頻時碼在空間的傳播特征及太陽活動對低頻時碼授時的影響，需要開展長期監(jiān)測及研究。

[1] 李國棟. 基于FPGA的低頻時碼數(shù)字接收機設(shè)計與實現(xiàn)[D]. 西安: 中國科學(xué)院研究生院(國家授時中心), 2007.

[2] 顧卓璟. 低頻時碼接收系統(tǒng)的研究設(shè)計[D]. 南京: 南京大學(xué), 2012.

[3] 靳建新. 長波天波授時可行性分析[J]. 全球定位系統(tǒng), 2012, 37(3): 5-8+15.

[4] 翟厚濤. VF/ELF電磁波在電離層中的傳播機理研究[D]. 杭州: 浙江大學(xué), 2007.

[5] 姜濤. 地—電離層波導(dǎo)中甚低頻波晝夜過渡期特性研究[J]. 通信電源技術(shù), 2014, 31(4): 43-45.

[6] 劉小花, 許林生, 華宇. 基于FPGA的便攜式BPC定時接收機設(shè)計[J]. 時間頻率學(xué)報, 2012, 35(2): 88-95.

[7] 白燕, 吳貴臣, 許林生. 低頻時碼接收機高精度定時的研究[J]. 電子測量與儀器學(xué)報, 2006, 20(3): 102-105.

[8] 李文浩. 低頻授時設(shè)備優(yōu)化及管理研究—基于商丘低頻臺[D]. 鄭州: 鄭州大學(xué), 2020.

[9] 梁志峰, 秦放, 崔方.“6·21”日環(huán)食對光伏發(fā)電及電網(wǎng)運行影響分析[J]. 電力系統(tǒng)自動化, 2021, 45(7): 1-7.

[10] 何玉囡, 王桂芬.2020—2049年日食的見食情況[J]. 紫金山天文臺臺刊, 1993, 12(2): 123-139.

[11] 胡昌林, 王蕾. 一種基礎(chǔ)模糊理論的自適應(yīng)數(shù)據(jù)平滑方法[J]. 現(xiàn)代雷達, 2016, 38(7): 49-51+86.

[12] 季國良. 1987年日環(huán)食的氣象效應(yīng)研究取得成果[J]. 中國科學(xué)基金, 1995(4): 55-56.

[13] “1987年9月23日日環(huán)食期間太陽輻射和大氣重力波的分析研究”課題組. 太陽輻射和地面氣象要素的日環(huán)食效應(yīng)[J]. 地球物理學(xué)報, 1990(4): 399-407.

[14] 潘煉德, 苗永瑞. 日環(huán)食對低頻及甚低頻電波傳播的影響[J]. 電波科學(xué)學(xué)報, 1990, 5(1): 38-43.

Observations of low frequency time-code signal during the annular eclipse of June 21, 2020

DAI Shun1,2, LUOJin-tao2, LIU Hai-wen3, LI Shi-feng2, LI Yi-feng2,LI Hao2, LIU Chang-Hong2, LIU Wei2, LIU Yong-peng2

(1. East China Jiaotong University, Nanchang 330013, China；2. National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China;3. Xi’an Jiaotong University, Xi’an 710049, China)

Around the annular eclipse on June 21, 2020, observing systems were installed in Lintong District of Xi’an City and Luonan County of Shangluo City to monitor the 68.5 kHz signal from the Shangqiu BPC station. The recorded data was processed and visualized. Analysises of the signal’s fluctuations were carried out, and the results showed that the amplitude of the low frequency time-code signal fluctuated during the solar eclipse and reached an extreme value after the maximum of the eclipse. Analysises of the signal were implemented based on the change of the ionosphere. The signal’s fluctuations can be caused by the sunrise-sunset effects, and it is suggested that the propagation of the low frequency time-code signal is affected by the solar eclipse. This work provides a reference for studying the low frequency time-code signal during solar eclipses.

annular eclipse; low frequency time-code; ionospheric

10.13875/j.issn.1674-0637.2021-03-0215-07

戴順, 羅近濤, 劉海文, 等.2020年6月21日日環(huán)食低頻時碼信號監(jiān)測[J]. 時間頻率學(xué)報, 2021, 44(3): 215-221.

2021-04-17；

2021-05-19

國家自然科學(xué)基金資助項目（11973046；91736207；42030105；11873050）