2020年6月21日日環食低頻時碼信號監測

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

2020年6月21日日環食低頻時碼信號監測

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

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

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

日環食；低頻時碼；電離層

0 引言

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

低頻時碼的傳播主要有地波和天波兩種方式。地波沿地球表面以繞射的方式到達接收點，天波經電離層反射傳播至接收點[3]。電離層在垂直方向呈現出分層結構，一般分為D層、E層和F層，其中，F層又分為F1層和F2層。夜間D層、E層和F1層基本消失，只剩下F2層仍可用于無線電傳播[4]。這種電離層的晝夜的變化，使低頻時碼信號在夜間強于白晝[5]，呈現出“晝弱夜強”的特點。太陽活動會影響電離層，進而對低頻時碼信號傳播產生影響。

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

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

1 2020年6月21日日環食及觀測

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

表1 日環食時刻及食甚時食分

注：本文出現時間皆為北京時間（UTC+8h）。

本次監測分別在西安臨潼及商洛洛南設置觀測點，對商丘BPC臺站發播的68.5 kHz低頻時碼信號進行監測，監測系統如圖2所示。

圖2 觀測系統

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

2 數據處理

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

圖3 采集數據后處理過程

采集數據的整體有效性受外部不穩定因素影響，如遠程連接觀測設備導致采集中斷等。以保證數據可靠性為前提，在對現有數據進行進一步處理前，應先篩出有效數據。有效數據應符合以下幾個條件：

① 連續每秒記錄數據，數據丟失量不超過整體數據的3%；

② 連續時間內，相同數據值不超過5條；

③ 在以天為單位的連續數據中，無數據中斷現象。

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

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

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

表2 觀測點日環食時刻表

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

3 觀測結果分析

通過24 h內信號的變化曲線可以看出，白天時段（08:00—18:00）信號波動幅度較為平穩，夜間時段（20:00—24:00，00:00—5:00）及晝夜過渡段（18:00—20:00）信號波動幅度較大。這是因為信號經電離層反射，白天時段內電離層相對穩定，夜間時段及晝夜過渡時段內的電離層處于變化狀態。其中， 05:00—08:00為BPC發播臺停機檢修時間。兩地監測68.5 kHz信號皆呈現“晝弱夜強”的整體趨勢，符合低頻時碼傳播隨電離層不同發生的晝夜變化，且在日環食期間信號發生異常時刻相近、時長相似。

3.1 臨潼觀測站

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

圖6 臨潼觀測日環食期間68.5 kHz信號變化

3.2 洛南觀測站

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

圖7 洛南觀測日環食期間68.5 kHz信號變化

4 結語

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

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

[2] 顧卓璟. 低頻時碼接收系統的研究設計[D]. 南京: 南京大學, 2012.

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

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

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

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

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

[8] 李文浩. 低頻授時設備優化及管理研究—基于商丘低頻臺[D]. 鄭州: 鄭州大學, 2020.

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

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

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

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

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

[14] 潘煉德, 苗永瑞. 日環食對低頻及甚低頻電波傳播的影響[J]. 電波科學學報, 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日日環食低頻時碼信號監測[J]. 時間頻率學報, 2021, 44(3): 215-221.

2021-04-17；

2021-05-19

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