時間頻率專業數據集

時間頻率專業數據集
作者：章宇　王燕平　袁海波　董紹武
2020年6月23日

本作品收錄於《中國科學數據》

章宇, 王燕平, 袁海波, 董紹武.時間頻率專業數據集[J/OL]. 中國科學數據, 2020, 5(2). (2020-06-22). DOI: 10.11922/csdata.2019.0077.zh.

摘要&關鍵詞

摘要：標準時間頻率系統是國家重要的戰略基礎設施，在國民經濟和國家安全方面起着基礎保障作用。時間頻率數據資源是國家重要信息資源之一，推動着時間頻率中守時、授時和時間頻率應用的發展。在大數據融合的背景下，時間頻率專業數據集在中國科學院科學大數據工程項目支持下應運而生。本數據集通過依託中國科學院國家授時中心（NTSC）定期向國際計量局（BIPM）提交我國國家標準時間系統測量比對數據、國際計量局國際原子時（TAI）歸算數據、全球各守時實驗室保持的本地時間UTC (k)與協調世界時(UTC)的偏差數據，及各實驗室原子鐘在國際原子時歸算中所占的權重公報（w公報）、速率公報（r公報）等數據，通過數據字典匹配，生成TAI(k)、UTC (k)、參與TAI歸算的原子鐘權重、速率等數據集。此外，數據集擴展包括《時間頻率公報》電子版以及收集到的國際地球自轉服務機構（IERS）公布的地球自轉A、B、C、D公報相關內容。本數據集為2018–2019年數據，採用時頻科學數據規程控制數據質量。未來，隨着數據庫資源的不斷完善，時間頻率專業數據集可以為時頻科學研究及國家標準時間頻率用戶提供更加豐富準確的基礎科學數據支持。

關鍵詞：時頻體系；原子鐘權重；原子鐘速率；NTSC時頻公報；地球自轉公報；數據質量

Abstract & Keywords

Abstract: The standard time-frequency system is an important strategic infrastructure of China, and it plays a fundamental role in safeguarding the national economy and security. Time-frequency data is one of the important national information resources. The dataset of time and frequency system includes TAI (International Atomic Time), UTC (NTSC) (Coordinate Universal Time maintained by National Time Service Center of Chinese Academy of Sciences), atomic clock weights, atomic clock rates, the electronic version of the Time and Frequency Bulletin released by NTSC and the Earth Rotation Bulletins A, B, C, and D published by the International Earth Rotation Service (IERS), etc. The current dataset contains data from 2018 to 2019. The data quality is controlled through time-frequency scientific data processing procedures. In the future, with continuous improvement of database resources, the time-frequency system dataset can provide more abundant and accurate basic scientific data for time-frequency scientific research and national standard time-frequency users.

Keywords: time-frequency system; atomic clock weight; atomic clock rate; time-frequency dataset; earth rotation bulletins; data quality

數據庫（集）基本信息簡介

數據庫（集）名稱	時間頻率專業數據集
數據作者	章宇、王燕平、袁海波、董紹武
數據通信作者	袁海波（yuanhb@ntsc.cn）
數據時間範圍	2018–2019年
地理區域	全國
數據格式	*.xls*.txt
數據服務系統網址	http://www.sciencedb.cn/dataSet/handle/917
基金項目	中國科學院信息化專項科學大數據工程（XXH135）
數據庫（集）組成	數據集由7部分數據組成，主要包含原子鐘權重、原子鐘速率、NTSC時頻公報、地球自轉公報A報、地球自轉公報B報、地球自轉公報C報、地球自轉公報D報。

Dataset Profile

Title	A dataset of time and frequency system
Data corresponding author	Haibo Yuan (yuanhb@ntsc.cn)
Data authors	Yu Zhang, Yanping Wang, Haibo Yuan,Shaowu Dong
Time range	2018–2019
Geographical scope	Nationwide
Data format	*.xls*.txt
Data service system	<http://www.sciencedb.cn/dataSet/handle/917>
Source of funding	CAS informatization project during the Thirteenth Five-Year Plan (XXH135)
Dataset composition	This dataset consists of seven pieces of data, including atomic clock weights, atomic clock rates, NTSC time-frequency bulletin, and Earth Rotation Bulletins (A, B, C, and D).

引 言

時間是目前基本物理量中測量精度最高的物理量，在科學技術研究、發展與應用中具有基礎性作用，同時標準時間頻率信號具有實時服務的重要職能，因此時間頻率的作用不僅在基礎領域研究、工程技術領域的探索方面，對於國民經濟、社會民生等方方面面都有重大影響。中國科學院國家授時中心（NTSC，下簡稱國家授時中心）是我國的標準時間、標準頻率的產生、發播單位，在我國時頻科學領域有着舉足輕重的作用，為國民經濟建設做出了重大貢獻。國家授時中心產生並保持有我國時間頻率大量科學數據，匯集建設時頻科學數據集並進行發布共享，不僅對於基礎的時間頻率相關領域研究有重大意義，同時對於時頻設備的製造、導航技術的研究也能夠提供豐富的時頻基礎數據支持。在實現學科數據大融合的背景下，時間頻率特色數據集的建設對於促進科技創新和經濟社會發展有重大作用，填補了時間頻率學科科學數據集建設的空白。

國家授時中心有多台高性能原子鐘，是參與全球國際原子時歸算的重要實驗室^[1]。國家授時中心每周向國際計量局（BIPM，Bureau International Des Poids etMesures）發送每天UTC 0h的原子鐘比對數據，BIPM通過遠程時間比對鏈路獲得包括UTC (NTSC)在內的全球各守時實驗室原子鐘與國際比對中心站UTC (PTB)的時差，運用加權平均算法得到自由原子時（EAL），隨後經過聯合噴泉鍾校準得到國際原子時（TAI），進而得到UTC與各守時實驗室UTC (k)的偏差、各守時實驗室原子鐘一個月的平均速率值、原子鐘在國際原子時計算中的權重值，以T公報、w公報及r公報形式，每月月初通過BIPM網站發布，本月發布數據為上月計算結果。本論文依據BIPM公布的2018–2019年原子鐘的權重、速率報文，通過數據字典匹配關係得到UTC (NTSC)的原子鐘權重、速率信息，收集國家授時中心每月出版的時間頻率公報以及國際地球自轉服務機構（IERS）公布的地球自轉公報相關內容，對數據進行校驗、規範化處理，形成數據集。其中，數據字典即數據字段類型等信息，具體來講其結構包含BIPM交互數據TAI_TA編號、BIPM交互數據UTC_UTC編號、BIPM交互數據文件編號、科研機構（實驗室）編號、國家授時中心原子鐘自編號、原子鐘國際編號、原子鐘權重編號、原子鐘速率編號、原子鐘權重_儒略日、原子鐘速率_儒略日、原子鐘權重_權重、原子鐘速率_速率、原子鐘權重/速率是否標記為刪除等信息。

本文詳述了時間頻率特色數據集數據產生的過程、數據的結構、數據的預處理、同時驗證了數據集的有效性，對於數據集的價值也進行了一定的闡述。此外，對於時頻數據集的未來發展，如在多模導航系統兼容及互操作等方面應用進行了介紹。在生成數據集，管理數據集的過程中，同時形成數據管理規範，進而控制時頻科學數據質量，為更好地實現時頻科學數據共享提供強有力支撐，也為數據用戶提供完整的時頻科學數據服務提供保障。

1 數據採集和處理方法

UTC (NTSC)守時系統本地比對系統及遠程時間比對鏈路比對數據產生的原理如圖1所示。UTC (NTSC)每隔一個小時進行一次主鍾與各守時原子鐘的比對，並將所得的鐘差比對結果採集存儲於數據庫中，文件命名為mc及jm，更新頻率為每小時一次，即為源數據。每周向BIPM發送一次前一周測量的UTC早上八點比對數據，BIPM獲取到全球各實驗室數據，經過歸算得到權重、速率報文。我們在BIPM官網下載數據報文，並對數據報文信息提取、整合，隨後得到我們的原子鐘權重、速率數據集。

圖片

圖1 NTSC守時比對系統原理圖

此外，國家授時中心通過統計我國長波授時台（BPL）時間信號、短波授時台（BPM）的UTC、世界時UT1時號等相對於UTC (NTSC)的發播時間^[2]、採集BIPM公布的Circular-T公報協調世界時UTC與國家授時中心保持的協調時世界時UTC (NTSC)、綜合原子時UTC (JATC)偏差等形成時頻公報，通過獲取到時頻公報電子版進行整理形成時頻公報數據集。

通過在IERS官網下載地球自轉參數報文，完成地球自轉公報數據採集。

1.1 數據來源

1.1.1 原子鐘速率數據的產生過程

BIPM接收到國家授時中心在內的全球參與TAI歸算的各守時實驗室的鐘差數據、以GNSS PPP（全球衛星導航系統精密單點定位時間傳遞）和TWSTFT（衛星雙向時間頻率傳遞）等遠距離時間頻率比對技術得到的UTC (PTB)與UTC (NTSC)的相位偏差^[3]，以此二者為數據基礎利用ALOGOS算法^[4]可以計算得到EAL，對EAL經過校準得到TAI。由國際原子時TAI與世界時UT1^[5]協調獲得協調世界時UTC。

隨後，經過換算可以得到UTC與UTC (NTSC)的相位偏差，結合UTC (NTSC)和原子鐘clock(i)鐘差數據，可以求得UTC相對NTSC的各台原子鐘clock(i)相位差，對此數據進行最小二乘擬合可以得到NTSC各台原子鐘相對於UTC的速率。

1.1.2 原子鐘權重數據產生原理

首先，利用之前一個月迭代計算出的權重作為本次計算的權重，然後計算出鍾hi 相對於自由原子時EAL的相位差。

隨後，計算每個月的實際頻率值yi,I 與預報值Si,I 的絕對偏差：

\({\epsilon }_{i,I}=\left|{y}_{i,I}-{S}_{i,I}\right| \) （1）

其中i表示第i台鐘，I表示計算的時間間隔。隨後求出\({\epsilon }_{i,j}^{2}\)。

該算法認定新的測量數據更具有可靠性，所以利用式（2）求出方差：

圖片 （2）

其中，j表示計算間隔，M代表月份（範圍是5–12月，5個月是一台鐘參與UTC計算的最小周期，而一年12個月則是標準的觀測周期^[6]），N代表原子鐘個數。

則可以得到原子鐘的權重如式（3）：

圖片 （3）

由於該算法強調可預測性，當一台鐘可預測性較好的時候會占據較大的權重，一旦該鐘出現故障，則系統的可靠性就會大大降低，因此引入最大權限值的規則，即當原子鐘hi 的權重超過最大權限值的wmax =4/N的時候，值為w。

此外，當一台鐘的性能變差時，如檢驗到實際的鐘速與預測鐘速超過5 ns/day，即在計算間隔內預報鐘差與實際鐘差的計算結果超過150 ns，則該鐘權重為0。採用這種方法可剔除參與UTC計算1%的鐘數^[6]。

以上權重、速率值是每台原子鐘每月對應一個值。可通過BIPM官網下載相關報文，隨後人工上傳數據庫^[7]，數據庫經過數據字典匹配自動採集NTSC對應約化儒略日、對應鐘號的權重、速率值，至此，權重、速率數據集生成，可為原子鐘使用單位和個人，標準時間使用單位和個人共享查詢，隨後導出xlsx文件。

1.1.3 時頻公報數據的產生

國家授時中心多年來積累了大量的時間頻率科學數據，一方面，國家授時中心具有長期的國際合作經驗，在數據自動化傳輸方面，技術已經成熟可靠，能夠保障現有數據的完整性、準確性，另一方面國家授時中心作為我國標準時間產生和發播單位，產生大量高精度時間基準數據。國家授時中心依據嚴格的時頻科學數據質量控制規範以時間順序梳理匯總時頻數據，並以文本文件的形式保存，整合後以時間頻率公報形式，定期公開出版發行。本文時間頻率公報數據集是通過對Excel格式的《時間頻率公報》進行採集形成的，以便於用戶的下載以及使用。

1.1.4 地球自轉公報數據的產生

IERS通過統籌全世界各地的觀測台站數據，經過IERS中心進行匯總處理並進行分析，得出聯測地球自傳變化參數並發布，主要應用於地球變化的監測、地球定向參數標準及規範的提供等。本文地球自轉公報數據集主要收集地球自轉A、B、C、D公報，因其涉及國際閏秒的調整預報、國際協調世界時與國際原子時的偏差預報，世界時與協調世界時的偏差預報等時間頻率相關參數預報。通過對公報採集並整合，形成地球自轉公報數據集。

1.2 數據預處理

鐘差比對結果出現異常值，及由於斷電或者是採集程序的故障引起數據缺失，會使得一個正常的鐘獲得不合理的評估，使得較好的鐘反而獲得較小的權重。因此，在報送BIPM之前需要對數據進行預處理。

對於異常值剔除採用萊特準則，原理如下：

首先對實測鐘差數據進行擬合為Si ，則實測值與擬合數據的殘差為：

圖片 （4）

以殘差為基礎值，隨後計算方差：

圖片 （5）

其中，圖片為均值，假設在t1 時刻，滿足式(6)條件：

圖片 （6）

則視hi (t1 )為檢測出的異常值，隨後以擬合值替代異常值：

圖片 （7）

對於原子鐘比對結果的缺少值可利用線性插值或最小二乘擬合^[8]，其實質是求解二次多項式係數，二次項係數對應的是頻漂，一次項係數對應的是擬合鐘速。以此方法補齊缺失數據，使得數據完整性得到保障，進而提高計算結果的準確性。

1.3 數據規範化處理

對於採集到的數據進行人工的核驗工作，數據量小的時候採用遍歷的方法。如果對應原子鐘標號有數據缺失，再次進行導入報文操作。對於導入時系統出現錯誤提示，則需檢測報文是否完整、格式是否正確。數量大的時候，可以通過採取抽檢的方式或是依據文件約化儒略日日期來確認文件採集更新是否正常。需要明確的數據集更新頻率標準是BIPM的權重w公報每月更新一次，BIPM的速率r公報每月更新一次，NTSC時頻公報每月更新一次，地球自轉公報A報每周更新一次，地球自轉公報B報每月更新一次，地球自轉公報C報按需更新，地球自轉公報D報按需更新。

2 數據樣本描述

由上介紹，BIPM公布的r公報、w公報的數據來源是包括NTSC在內的全球各國和重要地區守時實驗室的原子鐘比對結果，以及通過遠距離比對技術溯源到國際時間比對匯聚中心（PTB）的比對結果，兩者結合獲得的。

BIPM原子鐘速率計算周期為一個月，每月初公布上月相關約化儒略日對應的速率計算結果，單位為ns/day。其中，不同數字編號對應了不同的鐘的類型，其中35-表示銫原子鐘、40-表示氫原子鐘，這兩種守時原子鐘占據了NTSC鍾類型的絕大多數，此外還有其他的鐘類型。報文形式如BIPM r19.08.txt，即圖2所示^[9]。

圖片

圖2 BIPM公布r報（部分）

本文原子鐘速率數據集由獲得的BIPM r公報經過數據字典匹配得到NTSC相應鐘號，對應約化儒略日的原子鐘速率信息生成，具體如圖3所示。

圖片

圖3 UTC (NTSC)原子鐘速率集（部分）

依據速率數據集，可分析原子鐘每月的速率變化情況，即頻漂情況。如果速率變化較大則說明原子鐘的性能下降，即原子鐘運行不穩定；反之，則說明原子鐘運行良好。

原子鐘權重的大小是國際原子時歸算時的參數和衡量原子鐘長期性能水平的標誌^[10]。根據原子鐘的性能確定其權重，可以充分發揮性能優良的原子鐘優勢，以及在原子時計算中的重要地位。BIPM原子鐘權重計算周期為一個月，每月初公布上月相關約化儒略日對應權重的計算結果。報文形式如BIPM w19.08.txt即圖4所示^[9]。

圖片

圖4 BIPM公布w報（部分）

本文原子鐘權重數據集由獲得的BIPM w公報經過數據字典匹配得到NTSC相應鐘號，對應約化儒略日的原子鐘權重信息生成，具體如圖5所示。

圖片

圖5 UTC (NTSC)原子鐘權重集（部分）

隨機選取NTSC的6台高性能原子鐘，其中氫鍾3台，銫鐘3台，分析2019年權重分布情況，相關結果如圖6所示。由圖可得銫鐘所占的權重普遍較小，氫鍾所占權重普遍較大。這是由於2014年後BIPM改進ALGOS算法更強調原子鐘的可預測性，如果原子鐘頻率變化幅度穩定，可預測性較強，則原子鐘的權重占比較高。氫鐘的短期穩定度較好，雖然有頻漂但漂移穩定，可預測性較強。銫鐘長期穩定性好，短期穩定性較差、頻漂不穩定難以預測，故氫鍾所占的比重整體高於銫鐘所占比。

圖片

圖6 氫、銫原子鐘的權重分布

時頻公報數據集數據結構除以上介紹外，還包括美國導航衛星全球定位系統（GPS）的時間信號等相對國家授時中心協調世界時系統UTC (NTSC)的時差數據，地方原子時TA(NTSC)與UTC (NTSC)的相位偏差，IERS提供的地球自轉參數的最終確定結果，同時還有相對IERS的UT1系統，利用我國短波授時系統給出了世界時時號的改正數，計算了以JYD為原點的地極坐標等^[2]。以此公報為基礎，整理並發布。

IERS的A公報主要給出地球方向參數、世界時UT1與協調世界時UTC差值及其每日間隔的誤差、對未來1年的預測值，地球自轉B公報主要提供了有關地球方向的最新信息，包括世界時間，地球極坐標和天極偏移，地球自轉C公報主要是提供UTC的閏秒信息，以及協調世界時UTC與國際原子時TAI的相對差值、地球自轉D公報則是給出UT1與UTC的差值，並在一定精度範圍內發送^[11]。

3 數據質量控制和評估

為保證數據質量，時間頻率數據集有效實現共享，結合《時頻科學數據質量控制規程》、《時頻科學數據庫數據上傳管理規定》（內部規程、規定）對數據集進行三階段控制，具體如下：

第一階段：數據格式規範劃分，分析原始時頻數據產生的規律、數據特點、文件類型，建立各種原始數據整理的規範格式，即形成時頻數據字典。

第二階段：數據導入原則，將時頻原始數據轉化成規範格式，提交到時頻數據庫中。一方面按照數據集數據增新的要求，按質按時提交新增數據到數據庫中；另一方面按倒時順序整理正式試運行前的數據，不斷追加到數據庫中，逐步提前數據集的起始數據的起點，最大限度地保護原始數據。

第三階段：數據集校驗，分為導入數據前人工校驗、數據庫校驗（報文格式、鍾號、MJD等校驗）、上傳數據後導出數據集抽檢。

4 數據使用方法和建議

數據集的使用方法和建議，從三個方面舉例說明：一是可開展對於守時實驗室守時能力的評估，守時實驗室參與計算TAI的原子鐘數量越多、原子鐘權重取值越大，說明守時實驗室的影響力越大，同時依據BIPM Circular-T公報中UTC與UTC (k)偏差數據，可評估守時實驗室守時能力的強弱，如依據國家授時中心時頻公報數據集可以得出2018–2019年UTC (NTSC)相對於UTC的相位偏差保持在±5 ns之內，進而可分析UTC (NTSC)的短期性能和長期性能，具體如圖7所示。結合BIPM Circular-T公報數據集可以比較得出，國家授時中心的守時水平已居於世界前列。二是可應用於對於原子鐘性能的評估，由BIPM守時原子鐘權重數據集可以直觀得出每台原子鐘的性能差異，也可以以BIPM守時原子鐘速率、權重數據集為基礎計算原子鐘的ALLAN方差[] ，評估原子鐘的長期穩定度和短期穩定度。此外還可開展異常值研究，如頻率、相位跳變和頻漂等。三是對於時間尺度^[12]的研究，利用時間頻率數據集資源，獲得原子鐘權重、速率數據和鐘差數據，開展原子時尺度研究，可驗證ALOGOS、指數濾波AT1或全新設計的時間尺度算法，為提高原子時的保持水平，保障其穩定性、獨立性，更好地駕馭控制本地標準時間的輸出，提高守時能力發揮重要作用。^[13][14]

圖片

圖7 UTC (NTSC)相對於UTC的相位偏差

5 數據價值

時頻專業數據集不僅包含國內時頻相關數據，而且包含國際重要時頻科學數據，數據來源較為廣泛，包括BIPM、IERS、NTSC等。同時在數據的採集加工過程中、數據質量的把控中形成了數據管理質量體系規範。

本數據集的發布與共享，可進一步促進我國時頻科學領域的研究，提高時頻科學數據應用水平。未來，隨着時間頻率數據集的不斷完善、共享，其意義將愈加深遠，如研究範圍將從時間頻率信號產生、時間尺度建立、時間頻率傳遞、時間頻率測量、時間頻率應用等各方面擴展延伸，可以開展時頻系統誤差分析研究，從而進一步促進我國守時、授時發播水平的提高，並為天基與地基多模導航技術研究提供重要的數據支撐。

參考文獻

1. 趙書紅. UTC (NTSC)控制方法研究[D]. 北京: 中國科學院研究生院（國家授時中心）, 2014.
2. ^ ^{2.0 2.1} 中國科學院國家授時中心. 時間頻率公報[M]. 西安: 中國科學院國家授時中心出版社, 1979.
3. 董紹武. 守時中的若干重要技術問題研究[D]. 北京: 中國科學院研究生院（國家授時中心）, 2007.
4. HANADO Y, IMAE M, AIDA M, et al. Algorithm of ensemble atomic time[J]. NICT,2003, 50:155-167.
5. 胡永輝,漆貫榮.時間測量原理[M].香港:香港亞太科學出版社,2000.
6. ^ ^{6.0 6.1} PANFILO G, ARIAS F. The Coordinated Universal Time (UTC)[J]. Metrologia,2019,56(4).
7. 王燕平,董紹武,袁海波. 時間頻率科學數據庫的建設和應用[C].全國時間頻率學術會議,四川成都,2009.
8. 費業泰. 誤差理論與數據處理[M]. 北京:機械工業出版社, 2015.
9. ^ ^{9.0 9.1} BIPM.Circular T[EB/OL].(2019–08) [2019–09].https://www.bipm.org/en/about-us/.
10. 胡錦倫,李樹洲,李大志.國際原子時進展中的原子鐘[J].宇航計測技術, 2002,22(05): 6-12.
11. IERS. A、B、C、D報[EB/OL]. (2018-2019) [2019–09].https://www.iers.org/IERS/EN/Home/home_node.html.
12. ARIAS E F, PANFILO G, PETIT G.Timescales at the BIPM[J].Metrologia,2011, 48:S145–S153.
13. 漆貫榮. 時間科學基礎[M]. 北京:高等教育出版社, 2006.
14. 袁海波. 原子鐘噪聲數字化模擬與鐘速動態預測模型研究[D]. 北京: 中國科學院研究生院（國家授時中心）, 2009.

數據引用格式

章宇, 王燕平, 袁海波, 董紹武.時間頻率專業數據集[DB/OL]. Science Data Bank, 2019. (2019-11-28). DOI: 10.11922/sciencedb.917.