技术摘要：
本发明提供了一种基于铯喷泉钟和氢原子钟组的地方原子时生成方法，用以解决现有技术中由于对氢原子钟频率漂移计算误差而导致的地方原子时标长期稳定度及准确度不高的问题。所述地方原子时生成方法采集各台原子钟与主钟的钟差数据和主钟相对于铯喷泉钟的频率差值；再读  全部
背景技术：
时间是宇宙的基本结构之一，日常生活中的时间是一个狭义的范畴，在地球上被 划分为24个不同的时区。不同的国家根据自己所在的时区，对日常时间进行了不同的规定。 所确定的时间标准，因国家的不同而不同，称为地方原子时，被广泛应用于日常生活、天文 导航、大地测量和航天器跟踪等。 古代有晨钟暮鼓在日常生活中报时，随着科学技术的进步，现在各国在确定时间 标准时，广泛采用铯原子钟和/或氢原子钟来建立原子时。但是，铯原子钟和/或氢原子钟等 实际的物理时钟存在缺陷，如每台单独的钟可能突然失效。因此通常采用一台新钟代替失 效的钟，或利用一个钟组，通过计算给出一个时间尺度。在计算这种组合时间时，必须减小 由于新钟加入或停用某个钟而造成对时间尺度的影响。钟组内钟越多，加入或剔除一个钟 的影响越小。 利用原子钟组进行平均原子时间标准(以下简称平均原子时标)的计算，可以获得 更好的长期稳定度和准确度，目前可以采用的计算方法包括加权平均、卡尔曼滤波、小波分 析等，其中，经典加权平均方法(ALGOS)应用最为广泛，国际原子时计算即采用ALGOS算法， 该算法的主要特点是权重选取主要依据原子钟的历史特性，充分发挥较好原子钟的性能， 抑制较差原子钟的性能，使计算得到的时间尺度在指定时间间隔的稳定度理论上优于钟组 内任何一台原子钟，并且钟组中加入或剔除一个钟时，仍能保持较好的性能。 ALGOS方法的基本原理是，利用守时原子钟之间的钟差，通过数据预处理和原子钟 性能分析对各守时原子钟进行赋权，根据加权平均的方法给出综合原子时(TA)的纸面值。 ALGOS方法仅考虑了原子钟的钟速变化影响，适应于铯钟的原子时计算，而氢钟相对于理想 时标存在着明显的频率漂移。频率漂移是指原子钟的频率随时间的变化率。因此，由于计算 得到的时间尺度与各台原子钟抵偿结果存在相关性，导致原子时计算结果存在频率漂移， 计算所得的时间尺度的长期稳定度将受氢原子钟的频率漂移影响而明显变差。 现有技术中，ALGOS方法中对氢原子钟漂移率的修正方案为：以各台原子钟加权平 均计算出的时间尺度为基准，计算氢原子钟的漂移率，并在加权平均方法中予以抵偿。该方 法可准实时的计算氢原子钟的频率漂移。但是，2012年国际计量局BIPM的一项研究表明，由 于加权平均计算出的时间尺度与各台氢原子钟存在极大的相关性，因此计算得到的时间尺 度自身带有频率漂移，以此时间尺度为基准计算得到的氢原子钟漂移率，并非其实际漂移 率，进而在时间尺度计算中未能完全抵偿频率漂移，导致时间尺度的频率漂移恶化。时间尺 度的频率漂移恶化，进一步导致氢原子钟频率漂移计算的误差增大。因此，原子时计算结果 长期存在频率漂移，导致原子时的频率与真实频率之间的偏差不断扩大，依据此方法所测 定的地方原子时间标准长期稳定度和准确度不高。 4 CN 111610710 A 说　明　书 2/6 页
技术实现要素：
本发明实施例为了解决现有技术中由于对氢原子钟频率漂移计算误差而导致的 地方原子时标长期稳定度及准确度不高的问题，提出了一种基于铯喷泉钟和氢原子钟的地 方原子时生成方法，采用非相关抵偿，计算结果中不存在频率漂移，结合氢原子钟短期稳定 度与铯原子喷泉钟长期稳定度的双重优点，使生成的时间尺度同时具备短期稳定度与长期 稳定度。 为了实现上述目的，本发明实施例采用如下技术方案。 一种基于铯喷泉钟和氢原子钟组的地方原子时生成方法，所述地方原子时生成方 法利用铯喷泉钟频率基准进行非相关频率漂移计算，基于铯喷泉钟的时间尺度算法，产生 时间尺度及地方原子时。 优选地，所述地方原子时生成方法包括如下步骤： 步骤S1，将原子钟组的秒信号接入多通道时间间隔计数器，由数据采集与处理工 控机采集数据并存入时频数据库，得到各台原子钟与主钟的钟差数据； 步骤S2，将主钟的5MHz信号接入铯喷泉频率基准，得到主钟相对于铯喷泉钟的频 率差值，将所述频率差值存入时频数据库； 步骤S3，采用基于铯喷泉钟的时间尺度算法，读取时频数据库中的钟差数据与铯 喷泉钟的频率差值，计算各台原子钟相对于铯喷泉钟的频率差值，利用最小二乘法计算各 台原子钟的非相关频率漂移； 步骤S4，在时间尺度计算中除去所述非相关频率漂移，得到基于铯喷泉钟的时间 尺度及地方原子时。 优选地，所述步骤S3的计算过程，具体包括如下步骤： 步骤S301，铯喷泉钟对主钟的准确度评定； 步骤S302，铯喷泉钟对其他各台原子钟的准确度评定，同时评定多台原子钟； 步骤S303，根据铯喷泉钟对主钟及其他各台原子钟的准确度评定得到准连续运行 铯喷泉钟与氢原子钟组的频率差值，并进一步计算各台原子钟相对于铯喷泉钟的频率漂 移。 优选地，所述步骤S301中铯喷泉钟对主钟的准确度评定，进一步包括：读取时频数 据库中的所述步骤S2存入的频率差值，对频率差值进行数据处理，消除铯喷泉钟系统频率 差值，得到铯喷泉钟对主钟的准确度评定； 所述步骤S302中铯喷泉钟对其他各台原子钟的准确度评定，进一步包括：每日提 取主钟与其他各台原子钟的钟差数据并进行处理，得到当日主钟与各台原子钟的频差，再 根据所述铯喷泉钟对主钟的准确度评定，得到铯喷泉钟对其他各台原子钟的准确度评定。 优选地，所述步骤S301和步骤S302中的评定，采用最小二乘法拟合对频率差值和 钟差数据进行处理。 优选地，所述系统频率差值，包括二级塞曼频移、自旋交换/冷碰撞频移、微波功率 频移、Majorana跃迁频移、黑体辐射频移、引力频移。 优选地，所述铯喷泉钟对主钟的准确度评定公式(1)为： 所述当日主钟与各台原子钟的频差公式(2)为： 5 CN 111610710 A 说　明　书 3/6 页 fi＝fMC-fClock-i   (2) 所述铯喷泉钟对其他各台原子钟的准确度评定公式(3)为： 式(1)、(2)和(3)中， 为铯喷泉钟与主钟的频率差值，fNIM5-M为铯喷泉钟频率值， fMC为主钟频率值，fi为当日主钟与第i台原子钟的频率差值，fClock-i为第i台原子钟频率值， 为铯喷泉钟与第i台原子钟的频率差值。 优选地，所述步骤S4中的时间尺度TA，进一步通过公式(4)进行计算： 式(4)中，Hm为主钟的钟面读数，Pi为第i个原子钟的频率漂移权重，ξi为第i台钟相 对于主钟的实测钟差， 为第i台钟相对于MAT的起点时间修正值，帽子“～”表示采用上个 计算周期计算结果，为预报值； 为第i台钟相对于MAT的钟速率，也是上个计算周期计算结 果；t0为起点时刻，t为当前时刻，n为原子钟组中原子钟的个数。 优选地，其中，权重Pi的计算，利用各台原子钟相对于TA的钟差改正所计算得到的 稳定度求得，计算公式为： 其中 为第i台原子钟的天稳定度，根据阿伦方差进行计算。 本发明具有如下有益效果： 本发明实施例基于铯喷泉钟和氢原子钟的地方原子时生成方法，基于激光冷却铯 原子喷泉钟，计算得到氢原子钟组漂移率，并进行非相关抵偿方法，使计算得到的时间尺度 没有因时间尺度与各台原子钟的相关性而产生的频率漂移，结合了氢原子钟短期稳定度与 铯原子喷泉钟长期稳定度的双重优点，使生成的时间尺度同时具备短期稳定度与长期稳定 度。 附图说明 为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用 的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本 领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的 附图。 图1为本发明实施例基于铯喷泉钟和氢原子钟的地方原子时生成方法原理示意 图； 图2为本发明实施例基于铯喷泉钟和氢原子钟的地方原子时生成方法流程示意 图。 6 CN 111610710 A 说　明　书 4/6 页