首页 > 论文范文 > 自然科学论文 > 天文学论文 > 天文观测论文 > 关于GPS实时钟差确定算法的研究

关于GPS实时钟差确定算法的研究

2020-07-13 175 上传者：管理员

摘要：针对目前常用的GPS实时钟差确定算法主要基于实时观测数据流实现,使产品可靠性和完整性受限于网络质量,及广播星历和IGU-P精度较低的问题,提出并实现基于小时观测文件拼接的实时钟差确定算法。该算法基于小时观测文件数据,通过将钟差估计与自适应超短期钟差预报相结合以确定实时钟差。30d实时在线结果显示,基于小时观测文件确定的实时钟差精度约为0.25ns,优于广播星历和IGU-P,与IGSRTS提供的实时钟差精度相当,且能保证产品的可靠性和完整性。

关键词：
可靠性
天文观测
实时
小时观测文件
钟差确定
钟差预报
加入收藏

随着GNSS系统的快速发展,基于广域(全球)跟踪站网络的精密实时应用已拓展到导航和地学领域的多个方面,如实时精密单点定位、地震海啸等灾害的实时监测、空间气象实时监测、对地观测、卫星实时精密定轨等。导航卫星的实时精密轨道与钟差确定技术作为GNSS精密实时服务的核心技术,已成为热门的前沿性研究课题之一。

目前,IGS提供的产品可供实时应用的包括广播星历、超快产品的预报部分(IGU-P)和实时服务。其中,广播星历卫星钟差精度(STD)约为2.5ns,IGU-P卫星钟差精度约为1.5ns,均不能满足高精度实时应用的需求。从2013-04-01起,IGS实时服务(RTS)正式发布运行,旨在通过提供广播星历的改正参数以满足高精度实时应用的需求。然而,IGS未公布RTS实时钟差产品参考时间的确切含义,同时IGSRTS产品可能受网络影响而中断。近年来,国际上CODE、JPL、Thales、BKG、GFZ[1,2,3]等研究机构均围绕实时产品展开深入研究。与此同时,我国相关学者[4,5,6,7,8,9]在实时钟差确定方面也取得一定进展。目前常用的实时钟差确定方法为非差法、历元间差分法和混合法,这些方法均基于实时观测数据流而实现,采用滤波方法进行参数估计。上述3种方法均依赖于网络,同时需要较长的收敛时间,这会在一定程度上影响产品的完整性和可靠性。

本文提出基于小时观测文件拼接的实时钟差确定方法,该方法基于文件形式的小时观测数据(小时文件或15min文件),而不是实时观测数据流,可削弱对网络的依赖性;同时,生成的实时钟差以文件形式播发,无需专用软件接收。该方法采用最小二乘法进行参数估计,不需要参数收敛时间。因此,该方法较传统的实时钟差确定方法,可提高实时产品的可靠性和完整性。

1、基于小时观测文件拼接的实时钟差算法确定

基于小时观测文件拼接的实时钟差确定方法,采用钟差估计和预报的方式确定实时钟差。该算法包含2个关键技术,即基于弧段观测数据的钟差估计和自适应超短期钟差预报。

根据更新频率不同,IGS提供的观测数据可分为天文件、小时文件、15min文件和实时数据流。其中,天文件更新延迟较长,实时数据流强依赖于网络且需要专用软件(BNC)接收,而小时文件和15min文件的快速、可靠更新使基于小时观测文件拼接的实时钟差确定具备可能性。

基于小时观测文件拼接的实时钟差确定可分为5个步骤,具体流程如下:

1)小时观测文件下载,可以按需选择小时文件或15min观测文件;

2)观测文件拼接,形成一定长度的观测弧段数据;

3)基于弧段观测数据,按照更新周期窗口滑动,估计钟差;

4)按需确定训练序列长度和预报弧段长度;

5)超短期钟差预报,形成实时钟差。

本文选择小时观测文件应用于实时钟差确定,如更新周期小于1h,则选择15min观测文件。

1.1基于弧段观测数据的钟差估计

与传统精密卫星钟差估计方法一致[7],基于弧段观测数据的钟差估计一般采用非差消电离层组合观测值,其相位和伪距观测值的误差方程可表示为:

式中,k为测站号,j为卫星号,i为相应的观测历元,c为真空中光速;Δtk(i)为接收机钟差,Δtj(i)为卫星钟差;δρjk,trop为对流层延迟影响,εjk,Φ(i)、εjk,P(i)为多路径、观测噪声等未模型化的误差影响;Pjk(i)、Φjk(i)为已消除电离层影响的相应卫星、测站和历元的组合观测值,vjk,Φ、vjk,P为其观测误差;λ为相应的波长,ρjk为信号发射时刻的卫星位置到信号接收机位置之间的几何距离。

基于弧段观测数据的钟差估计采用BerneseGNSSsoftwareV5.2软件[10],估计流程如图1所示。在该方法中,输入文件主要包括RINEX格式的弧段观测文件、初始站坐标(CRD)、地球自转参数(ERP)和轨道等文件,输出CLK格式的钟差文件。估计策略如表1所示。

图1基于弧段观测数据的钟差估计流程

表1基于弧段观测数据的钟差估计策略

1.2自适应超短期钟差预报

本文采用自适应超短期钟差预报方法,其流程如图2所示。

图2自适应超短期钟差预报流程

1.2.1数据预处理

钟差预报要求具有可靠、稳定的钟差建模数据。如果钟差数据存在粗差、跳变等异常现象,这会对钟差建模产生影响,导致钟差预报精度大大降低。数据预处理主要包括相位数据与频率数据的转换、粗差探测及修复和钟跳探测及修复等。

首先将钟差相位数据转换为频率数据,利用中位数法进行粗差探测,然后对剔除的粗差进行置零处理[11]。但如果数据距离零值较远,置零处理会产生新的粗差。因此,对于粗差剔除,本文采用内插的方法进行补齐。在钟差模型拟合预报过程中,算法需要可以自适应地处理钟跳,避免其污染拟合模型,具体算法可参考文献[12]。

1.2.2超短期钟差预报

对于短期钟差预报而言,最常用的为一次多项式模型和二次多项式模型,方法选择取决于实际卫星的频漂特性是否明显[13,14]。然而,具体到GPS系统,星载原子钟类型多样,将所有卫星钟归为固定不变的模型并不合理,为此,黄观文[15]提出采用拟合数据段的残差标准差进行自适应钟差模型确定。具体判断准则如下:

式中,RMSLine为一次多项式钟差拟合模型的残差标准差,RMSQuadratic为二次多项式钟差拟合模型的残差标准差。

2、基于小时观测文件拼接的实时钟差策略确定

基于小时观测文件拼接的实时钟差确定可采用钟差估计和预报的方式。首先基于钟差估计部分展开策略优选,然后确定基于小时观测文件拼接的实时钟差策略。

2.1影响因素分析

从式(1)、(2)可以看出,钟差估计的影响因素有轨道、ERP和观测数据,其中观测数据包括测站数目和观测弧段长度。基于这些影响因素,设置3种实验方案。方案1:设置不同长度的观测弧段估计钟差,并与IGS精密产品进行比较;方案2:基于方案1,只改变测站数量,估计钟差并与IGS精密产品进行比较;方案3:基于方案1,只改变轨道产品和ERP,估计钟差并与IGS精密产品进行比较。本文对钟差估计结果的精度评价采用二次差法。

2.1.1观测弧段长度影响

方案1采用全球均匀分布的120个IGS测站2016年第178～207天共30d的观测数据,并将STD、ERP及CRD固定为CODE相关产品,设置观测弧段为24h、12h、6h、3h、1h和0.5h,估计GPS卫星钟差。为便于和IGS钟差结果进行比较,估计卫星钟差时采用5min的时间间隔,且在钟差比较时选择GPS10为参考星。GPS04卫星状态为“不可用”,因此无观测数据及钟差结果。

图3、4为实验结果,第202～203天存在卫星机动,故不作统计。

图3方案1实验结果

图4方案1基于卫星的实验结果统计

实验表明,基于该算法的钟差估计,观测弧段设置越短,钟差估计结果精度越差,原因可能为在钟差计算的最后一步参数估计时,将模糊度参数估计为浮点解。弧段越短,数据量越少,模糊度估计越不精确,导致钟差结果精度越差。

从图3可以看出,绝大多数卫星的钟差估计结果遵循观测弧段越短、精度越差的规律。少数卫星钟差精度较差,这可能与卫星自身特性有关,如发射年限、卫星型号及卫星钟类型等因素。弧段长度与钟差精度关系为:当弧段分别设置为24h、12h、6h、3h、1h、0.5h时,大部分卫星钟差精度分别小于0.05ns、0.15ns、0.20ns、0.25ns、0.30ns、0.35ns。

2.1.2测站数目影响

为了比较测站数目对不同观测弧段钟差估计的影响,方案2与方案1数据时间段、输入设置、采样频率、观测弧段设置相同,但测站数目减少为60,生成的钟差数据如图5所示。从图5可以看出,当测站数目减少为方案1的一半时,仍然满足观测弧段越短钟差估计结果精度越差的规律。

图5方案2实验结果

比较方案1与方案2可知,当测站数目减少约一半时,与IGS相比,钟差精度明显降低,当弧段为24h、12h、6h、3h、1h、0.5h时,精度分别降低100%、73%、47%、48%、52%、58%。

2.1.3轨道和ERP影响

为了比较轨道对不同观测弧段钟差估计的影响,方案3将ERP和STD固定为比CODE产品精度略差的IGU-P相关产品,观测数据、测站数目、采样频率和弧段设置均保持与方案1一致。2016年第178～207天的统计结果如图6所示,其中第189～191天由于IGU-P轨道数据精度较差或缺失,故不作统计。

图6方案3实验结果

从图6可以看出,输入不同的轨道产品时,仍然满足观测弧段越短钟差估计结果精度越差的规律。对比图3和图6可知,由于轨道结果作为输入值参与钟差估计,所以当输入的轨道产品精度降低时,钟差精度也随之降低,但随着观测弧段减小,影响逐渐减弱。当弧段为24h、12h、6h、3h、2h、1h、0.5h时,精度分别降低120%、27%、5%、4%、4%、7%。

综上所述,观测弧段长度、测站数目、轨道及ERP精度都是影响钟差估计的重要因素。在基于弧段观测数据的钟差估计算法中,观测弧段越短钟差估计精度越差的原因可能为在利用非差法计算钟差时需要估计模糊度参数,而当观测弧段较短时,数据量较少,模糊度估计不准确,导致钟差结果精度较低。另外,在基于弧段观测数据的钟差估计算法中,轨道、ERP作为已知量输入,这也会对钟差造成一定影响。

当观测弧段长度一定时,提高钟差结果精度的关键应该为增加测站数量,其次再考虑轨道、ERP等因素,但钟差的估计效率会随着测站数量的增加而降低。因此在实际应用中,应综合考虑钟差估计效率与钟差精度,以求两者平衡。

2.2基于小时文件拼接的实时钟差策略确定

观测弧段长度、测站数目、轨道及ERP精度是影响钟差估计的重要因素,并且在超短期钟差预报时,训练序列长度、预报弧段长度将直接影响基于小时文件拼接确定的实时钟差的精度。同时,实时钟差确定效率也需要综合考虑。

综上所述,图7为基于小时文件拼接的实时钟差确定算法策略,观测弧段长度设置为24h,计算弧段长度为1h,训练序列长度为3h,预报弧段长度为2h,更新周期为1h。小时观测文件下载时间节点设置为每个整点的20min,实时钟差确定时间节点设置为每个整点的30min,实时钟差确定完成时间约为每个整点的45min。

图7基于小时文件拼接的实时钟差确定算法策略示意图

3、结果与分析

基于小时观测文件拼接的实时钟差确定自2017-04-26(DOY116)开始在线测试,测试基于IGS发布的小时观测文件和IGU轨道产品开展,测站数量约为100。将2017-04-26(DOY116)～2017-05-28(DOY148)的实时钟差结果与IGS最终产品进行比对分析,DOY121、DOY127和DOY144由于IGU轨道精度较差,故不作统计。统计每天及每个SESSION的实时钟差结果(图8、9)。30d的在线测试统计结果显示,基于小时观测文件拼接确定的实时钟差精度(STD)约为0.25ns,与IGSRTS精度相当,优于目前成熟的实时钟差产品广播星历。

图8实时钟差精度统计结果(按天统计)

图9实时钟差精度统计结果(按SESSION统计)

目前,GPS卫星搭载5种不同类型的卫星钟,按照不同的卫星钟类型,各选择1颗卫星(表2),按天统计实时钟差精度(图10)。图11为32颗卫星30d的统计结果,从图中可以看出,实时钟差精度从高到低分别为GPS01、GPS02、GPS07、GPS08。由于搭载短期稳定性能较差的Cs原子钟,GPS08与其他几颗卫星的实时钟差结果相差较大,出现明显的分层现象;而其他卫星均搭载Rb原子钟,实时钟差相差较小。

表2选取卫星列表

图10不同类型卫星钟实时钟差统计结果(按天统计)

图11每颗卫星30d结果

结果表明,基于小时观测文件拼接确定实时钟差,一方面可保证实时钟差的精度,得到精度约为0.25ns的GPS实时钟差,与IGSRTS精度相当,远优于广播星历和IGU-P;另一方面,该方法不再依赖实时数据流形式的观测数据,可摆脱网络质量对传统实时钟差确定算法的制约。基于小时观测文件拼接的实时钟差确定方法中包含钟差预报技术,而钟差预报必须考虑到原子钟本身具备的物理特性,且不同类型的原子钟的物理特性差异较大[10]。因此,基于小时观测文件拼接确定的实时钟差精度与卫星钟类型密切相关。

4、结语

本文针对GNSS系统实时性应用的发展需求,以及目前实时钟差产品及算法受网络质量制约的现象,提出基于小时观测文件拼接的实时钟差确定算法。该方法基于由小时观测文件拼接形成的弧段观测数据,结合钟差估计与超短期钟差预报确定实时钟差。基于该算法可得到STD约为0.25ns的GPS实时钟差产品,与IGSRTS精度相当。

相较于目前常用的GPS实时钟差确定算法,基于小时观测文件拼接的实时钟差算法具有以下特点:1)该算法不受网络质量的制约,并且可采用最小二乘法进行参数估计,无需收敛时间,可改善产品的可靠性和完整性;2)该算法生成的钟差产品以文件形式播发,无需专用软件进行接收和修正,便于推广应用;3)相较于目前成熟的GPS实时产品广播星历和IGU-P,该算法可大幅提高GPS实时钟差产品的精度。

本文只对GPS实时钟差展开研究,鉴于目前GNSS,尤其是我国自主研发的BDS-3系统的快速发展,下一步将开展多GNSS系统实时钟差的深入研究。

参考文献：

[4]丁文武.实时精密单点定位系统的若干关键技术研究[D].武汉:中国科学院测量与地球物理研究所,2012.

[5]楼益栋.导航卫星实时精密轨道与钟差确定[D].武汉:武汉大学,2008.

[6]楼益栋,施闯,周小青,等.GPS精密卫星钟差估计与分析[J].武汉大学学报:信息科学版,2009,34(1):88-91.

[7]蔡华,赵齐乐,楼益栋,等.精密卫星钟差确定系统的实现与精度分析[J].武汉大学学报:信息科学版,2009,34(11):1293-1296.

[8]李星星,徐运,王磊.非差导航卫星实时/事后精密钟差估计[J].武汉大学学报:信息科学版,2010,35(6):661-664.

[12]张成军,贾学东.接收机钟跳对GPS定位的影响及探测方法[J].测绘通报,2009(12):7-9.

[13]于合理,郝金明,刘伟平,等.卫星钟差超短期预报模型分析[J].大地测量与地球动力学,2014,34(1):161-164.

[15]黄观文.GNSS星载原子钟质量评价及精密钟差算法研究[D].西安:长安大学,2012.

武美芳,苏行,张喆.基于小时观测文件拼接的GPS实时钟差确定算法[J].大地测量与地球动力学,2020,40(07):730-735.

基金:中国科学院“西部青年学者”项目(XAB2018B19).

网友评论

我要评论

期刊名称：天文学报

期刊人气：1799

期刊详情

主管单位：中国科学院

主办单位：中国天文学会

出版地方：江苏

专业分类：科学

国际刊号：0001-5245

国内刊号：32-1113/P

邮发代号：2-818

创刊时间：1953年

发行周期：双月刊

期刊开本：16开

见刊时间：一年半以上

论文导航

91学术

关于GPS实时钟差确定算法的研究

相关论文

典型模拟火星土壤探究进程

关键词

月球钍元素发电技术实践探究

关键词

机器学习基础上日冕仪图像分类方法探究

关键词

冕流电流片的白光和紫外观测探究分析

关键词

地球表面非构造季节性运动的物理因素探究

关键词

适用于SKA低频孔径阵列的接收模块设计研究

关键词

关于空间天气预警中太阳风分类研究的前沿动态

关键词

有噪中度量子时代及黑洞之光的开始

关键词