留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法

张万威 王甫红 龚学文 郭磊

张万威, 王甫红, 龚学文, 郭磊. 一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
引用本文: 张万威, 王甫红, 龚学文, 郭磊. 一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
ZHANG Wanwei, WANG Fuhong, GONG Xuewen, GUO Lei. A Centimeter-Level Real-Time Orbit Determination Method Using Space-Borne GPS Measurements Considering IGS-RTS Data Receiving Interruption[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
Citation: ZHANG Wanwei, WANG Fuhong, GONG Xuewen, GUO Lei. A Centimeter-Level Real-Time Orbit Determination Method Using Space-Borne GPS Measurements Considering IGS-RTS Data Receiving Interruption[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432

一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法

doi: 10.13203/j.whugis20200432
基金项目: 

国家自然科学基金 91638203

国家自然科学基金 62073044

武汉大学地球空间环境与大地测量教育部重点实验室开放基金 19-02-03

详细信息
    作者简介:

    张万威, 博士生, 主要从事卫星导航定位与卫星定轨理论与方法研究。wwzhang@sgg.whu.edu.cn

    通讯作者: 王甫红, 博士,教授。fhwang@sgg.whu.edu.cn
  • 中图分类号: P228

A Centimeter-Level Real-Time Orbit Determination Method Using Space-Borne GPS Measurements Considering IGS-RTS Data Receiving Interruption

Funds: 

The National Natural Science Foundation of China 91638203

The National Natural Science Foundation of China 62073044

the Open Research Fund of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University 19-02-03

More Information
    Author Bio:

    ZHANG Wanwei, PhD candidate, specializes in the theory and method of GNSS navigation and satellite orbit determination. E-mail: wwzhang@sgg.whu.edu.cn

    Corresponding author: WANG Fuhong, PhD, professor, E-mail: fhwang@sgg.whu.edu.cn
  • 摘要: 当前,在低轨导航增强、对地观测和科学应用等领域,低轨卫星对轨道参数的精度和实时性提出了更高的要求。利用国际GNSS服务组织实时服务(international global navigation satellite system service real-time service,IGS-RTS)播发的GPS卫星轨道与钟差改正数,针对极区实时改正数接收中断情况下实现厘米级星载GPS实时定轨的关键问题开展研究。首先分析了实时改正数及其短时外推引起的星历综合误差的变化特性;然后以此为依据,在星载GPS实时精密定轨数学模型中构建分段随机游走的伪模糊度参数随机模型,以减小星历综合误差对实时精密定轨的影响,从而实现厘米级精度的实时定轨。采用自主研制的实时精密定轨软件SATODS,使用GPS广播星历与CLK93实时产品,对为期一周的重力场恢复和气候实验(gravity recovery and climate experiment,GRACE)C卫星的GPS双频实测数据模拟在轨实时精密定轨处理。实验结果表明,在考虑极区改正数接收中断的情况下,所提实时精密定轨方法可以达到7.04 cm的位置精度以及0.20 mm/s的速度精度,所提方法具有可行性和有效性。
  • 图  1  CLK93产品的轨道与钟差精度

    Figure  1.  Orbit and Clock Accuracies of CLK93 Products

    图  2  各GPS卫星的轨道与钟差精度

    Figure  2.  Orbit and Clock Accuracies for Each GPS Satellite

    图  3  ELOS相关参数的变化曲线

    Figure  3.  Variation Curves of Parameters Related to ELOS

    图  4  5种方案的每天位置和速度精度统计

    Figure  4.  Daily Position and Velocity Accuracy Statistics of Five Schemes

    图  5  4种方案的三维位置误差比较

    Figure  5.  3D Position Error Comparison of Schemes 2, 3, 4 and 5

    图  6  方案3和方案4定轨位置误差对比

    Figure  6.  Comparison of Orbit Errors for Schemes 3 and 4

    表  1  实时精密定轨策略及设置

    Table  1.   Stratrgies of Real-Time Precise Orbit Determination

    模型参数 相关设置
    GPS观测数据 GPS双频伪距、载波相位无电离层组合
    EOP参数 IERS Bulletin A预报
    地球重力场模型 EGM2008 70×70
    N体引力 仅考虑日月,近似公式计算日月位置
    固体潮 简易固体潮模型
    大气阻力 改进后的Harris-Priester模型[14],固定表面积,估计大气阻力系数
    太阳光压 简单的球状模型,固定表面积,估计太阳光压系数
    经验加速度 RAC 3个方向的动力学模型补偿一阶高斯-马尔可夫随机模型
    参数估计方法 扩展卡尔曼滤波
    测量更新间隔/s 30
    下载: 导出CSV

    表  2  GRACE C卫星的实时精密定轨结果统计

    Table  2.   Accuracy Statistics of Real-Time Orbit Determination for GRACE C

    方案 位置精度/cm 速度精度/(mm·s-1
    R方向 A方向 C方向 3D R方向 A方向 C方向 3D
    1 10.3 26.5 12.3 30.9 0.298 0.133 0.161 0.363
    2 4.74 5.32 4.00 8.18 0.144 0.097 0.099 0.201
    3 3.74 4.85 3.45 7.04 0.139 0.096 0.103 0.197
    4 3.12 3.68 2.72 5.55 0.134 0.094 0.099 0.191
    5 2.85 2.85 2.54 4.78 0.131 0.093 0.099 0.188
    下载: 导出CSV
  • [1] Li Xingxing, Ma Fujian, Li Xin, et al. LEO Constellation-Augmented Multi-GNSS for Rapid PPP Convergence[J]. Journal of Geodesy, 2019, 93 (5): 749-764 doi:  10.1007/s00190-018-1195-2
    [2] Li Bofeng, Ge Haibo, Ge Maorong, et al. LEO Enhanced Global Navigation Satellite System (LeGNSS)for Real-Time Precise Positioning Services[J]. Advances in Space Research, 2019, 63 (1): 73-93 doi:  10.1016/j.asr.2018.08.017
    [3] Montenbruck O, Hauschild A, Andres Y, et al. Real-Time Orbit Determination for GNSS Radio Occultation Processing[J]. GPS Solutions, 2013, 17(2): 199-209 doi:  10.1007/s10291-012-0271-y
    [4] 王甫红, 龚学文, 刘万科. 顾及广播星历误差分离吸收的分米级星载GPS实时定轨新方法[J]. 武汉大学学报·信息科学版, 2015, 40(9): 1 230-1 236 doi:  10.13203/j%20.whu%20g%20is20130682

    Wang Fuhong, Gong Xuewen, Liu Wanke. A Novel Decimeter-Level Real-Time Orbit Determination Algorithm Using Space-Borne GPS Measurements with Separation and Absorption of Broadcast Ephemeris Error[J]. Geomatics and Information Science of Wuhan University, 2015, 40(9): 1 230-1 236 doi:  10.13203/j%20.whu%20g%20is20130682
    [5] Wang Fuhong, Gong Xuewen, Sang Jizhang, et al. A Novel Method for Precise Onboard Real-Time Orbit Determination with a Standalone GPS Receiver [J]. Sensors, 2015, 15(12): 30 403-30 418 doi:  10.3390/s151229805
    [6] 凌三力, 王甫红, 张万威. 风云三号C卫星星载GPS/BDS实时定轨分析[J]. 测绘地理信息, 2021, 46(3): 17-20 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202103005.htm

    Ling Sanli, Wang Fuhong, Zhang Wanwei. Analysis on the Real-Time Orbit Determination for FY3C Satellite Using Space-Borne GPS/BDS[J]. Journal of Geomatics, 2021, 46(3): 17-20 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202103005.htm
    [7] Wermuth M, Hauschild A, Montenbruck O, et al. TerraSAR-X Precise Orbit Determination with Real-Time GPS Ephemerides[J]. Advances in Space Research, 2012, 50(5): 549-559 doi:  10.1016/j.asr.2012.03.014
    [8] Sun Xiucong, Han Chao, Chen Pei. Precise RealTime Navigation of LEO Satellites Using a SingleFrequency GPS Receiver and Ultra-Rapid Ephemerides[J]. Aerospace Science and Technology, 2017, 67: 228-236 doi:  10.1016/j.ast.2017.04.006
    [9] Giordano P, Zoccarato P, Otten M, et al. P2OD: Real-Time Precise Onboard Orbit Determination for LEO Satellites[C]// The 30th International Technical Meeting of The Satellite Division of the Institute of Navigation, Portland, Oregon, America, 2017
    [10] 师一帅. 低轨卫星实时运动学精密定轨方法研究[D]. 郑州: 信息工程大学, 2018

    Shi Yishuai. Research on Real-Time Kinematic Precise Orbit Determination of Low-Earth-Orbit Satellites[D]. Zhengzhou: Informatiom Engineering University, 2018
    [11] Hauschild A, Tegedor J, Montenbruck O, et al. Precise Onboard Orbit Determination for LEO Satellites with Real-Time Orbit and Clock Corrections[C]// The 29th International Technical Meeting of the Satellite Division of the Institute of Navigation, Portland, Oregon, America, 2016
    [12] 王甫红, 凌三力, 龚学文, 等. 风云三号C卫星星载GPS/ BDS分米级实时定轨模型研究[J]. 武汉大学学报· 信息科学版, 2020, 45(1): 1-6 doi:  10.13203/j.whugis20180385

    Wang Fuhong, Ling Sanli, Gong Xuewen, et al. Decimeter-Level Orbit Determination for FY3C Satellite Based on Space-Borne GPS/BDS Measurements[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 1-6 doi:  10.13203/j.whugis20180385
    [13] 楼益栋, 施闯, 周小青, 等. GPS精密卫星钟差估计与分析[J]. 武汉大学学报·信息科学版, 2009, 34(1): 88-91 doi:  10.13203/j.whugis20180385

    Lou Yidong, Shi Chuang, Zhou Xiaoqing, et al. Realization and Analysis of GPS Precise Clock Products[J]. Geomatics and Information Science of Wuhan University, 2009, 34(1): 88-91 doi:  10.13203/j.whugis20180385
    [14] Montenbruck O, Gill E, Lutze F H. Satellite Orbits: Models, Methods, and Applications[J]. Applied Mechanics Reviews, 2002, 55(2): B27 doi:  10.1115/1.1451162
  • [1] 杨宇飞, 杨元喜, 徐君毅, 许扬胤, 赵昂.  低轨卫星对导航卫星星座轨道测定的增强作用 . 武汉大学学报 ● 信息科学版, 2020, 45(1): 46-52. doi: 10.13203/j.whugis20180215
    [2] 王甫红, 凌三力, 龚学文, 郭磊.  风云三号C卫星星载GPS/BDS分米级实时定轨模型研究 . 武汉大学学报 ● 信息科学版, 2020, 45(1): 1-6. doi: 10.13203/j.whugis20180385
    [3] 柳聪亮, 白伟华, 夏俊明, 孙越强, 孟祥广, 杜起飞.  低轨卫星星载GNSS反射事件模拟分析 . 武汉大学学报 ● 信息科学版, 2018, 43(6): 826-831, 839. doi: 10.13203/j.whugis20160161
    [4] 宫晓春, 王宇谱, 王宁, 翟树峰, 吕志平.  GPS卫星钟差改正数实时预报算法 . 武汉大学学报 ● 信息科学版, 2018, 43(6): 867-873. doi: 10.13203/j.whugis20160154
    [5] 张兵兵, 王正涛, 冯建迪, 邱耀东.  伪随机脉冲先验值对低轨卫星简化动力学定轨精度的影响 . 武汉大学学报 ● 信息科学版, 2018, 43(8): 1222-1227, 1241. doi: 10.13203/j.whugis20160168
    [6] 田英国, 郝金明, 陈明剑, 于合理, 衡培深.  GPS卫星钟差及观测数据采样间隔对LEO卫星定轨精度影响 . 武汉大学学报 ● 信息科学版, 2017, 42(12): 1792-1796. doi: 10.13203/j.whugis20150591
    [7] 王甫红, 龚学文, 刘万科.  顾及广播星历误差分离吸收的分米级星载gps实时定轨新方法 . 武汉大学学报 ● 信息科学版, 2015, 40(9): 1230-1236. doi: 10.13203/j .whu g is20130682
    [8] 马洋, 欧吉坤, 袁运斌, 霍星亮, 丁文武.  导航卫星天线相位中心变化估计及对LEO精密定轨影响 . 武汉大学学报 ● 信息科学版, 2015, 40(7): 894-900. doi: 10.13203/j.whugis20130626
    [9] 王甫红, 徐其超, 龚学文, 张蔚.  GAAF在星载GPS实时定轨中的应用研究 . 武汉大学学报 ● 信息科学版, 2014, 39(1): 47-51.
    [10] 张小红, 李盼, 左翔.  固定模糊度的精密单点定位几何定轨方法及结果分析 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1009-1013.
    [11] 刘俊宏, 谷德峰, 赖育网, 易东云.  DORIS实时定轨中的频偏估计方法研究 . 武汉大学学报 ● 信息科学版, 2012, 37(8): 956-959.
    [12] 王甫红.  高精度星载GPS实时定轨卡尔曼滤波模型 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 653-656.
    [13] 周晓青, 胡志刚, 张新远.  低轨卫星星载GNSS精密定轨的精度检核方法 . 武汉大学学报 ● 信息科学版, 2010, 35(11): 1342-1345.
    [14] 王正涛, 靳祥升, 党亚民, 姜卫平.  低轨卫星精密定轨的初轨向量与力模型参数数值积分误差分析 . 武汉大学学报 ● 信息科学版, 2009, 34(6): 728-731.
    [15] 施闯, 李敏, 楼益栋, 邹蓉.  利用区域基准站进行导航卫星近实时精密定轨研究 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 697-700.
    [16] 楼益栋, 施闯, 葛茂荣, 赵齐乐.  GPS卫星实时精密定轨及初步结果分析 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 815-817.
    [17] 韩保民, 朱秀英, 柳林涛, 曲国庆.  伪随机脉冲估计及其在简化动力学定轨中的应用 . 武汉大学学报 ● 信息科学版, 2007, 32(5): 466-469.
    [18] 韩保民.  精密卫星钟差加密方法及其对星载GPS低轨卫星定轨精度影响 . 武汉大学学报 ● 信息科学版, 2006, 31(12): 1075-1078.
    [19] 王甫红, 刘基余.  星载GPS载波相位测量的周跳探测方法研究 . 武汉大学学报 ● 信息科学版, 2004, 29(9): 772-774.
    [20] 林高宇, 王磊, 何飞扬, 宋晓迪, 郭际明.  Swarm低轨卫星星座的GPS接收机差分码偏差估计 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200479
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  254
  • HTML全文浏览量:  107
  • PDF下载量:  62
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-10
  • 刊出日期:  2021-11-05

一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法

doi: 10.13203/j.whugis20200432
    基金项目:

    国家自然科学基金 91638203

    国家自然科学基金 62073044

    武汉大学地球空间环境与大地测量教育部重点实验室开放基金 19-02-03

    作者简介:

    张万威, 博士生, 主要从事卫星导航定位与卫星定轨理论与方法研究。wwzhang@sgg.whu.edu.cn

    通讯作者: 王甫红, 博士,教授。fhwang@sgg.whu.edu.cn
  • 中图分类号: P228

摘要: 当前,在低轨导航增强、对地观测和科学应用等领域,低轨卫星对轨道参数的精度和实时性提出了更高的要求。利用国际GNSS服务组织实时服务(international global navigation satellite system service real-time service,IGS-RTS)播发的GPS卫星轨道与钟差改正数,针对极区实时改正数接收中断情况下实现厘米级星载GPS实时定轨的关键问题开展研究。首先分析了实时改正数及其短时外推引起的星历综合误差的变化特性;然后以此为依据,在星载GPS实时精密定轨数学模型中构建分段随机游走的伪模糊度参数随机模型,以减小星历综合误差对实时精密定轨的影响,从而实现厘米级精度的实时定轨。采用自主研制的实时精密定轨软件SATODS,使用GPS广播星历与CLK93实时产品,对为期一周的重力场恢复和气候实验(gravity recovery and climate experiment,GRACE)C卫星的GPS双频实测数据模拟在轨实时精密定轨处理。实验结果表明,在考虑极区改正数接收中断的情况下,所提实时精密定轨方法可以达到7.04 cm的位置精度以及0.20 mm/s的速度精度,所提方法具有可行性和有效性。

English Abstract

张万威, 王甫红, 龚学文, 郭磊. 一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
引用本文: 张万威, 王甫红, 龚学文, 郭磊. 一种顾及IGS-RTS数据接收中断的厘米级星载GPS实时定轨方法[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
ZHANG Wanwei, WANG Fuhong, GONG Xuewen, GUO Lei. A Centimeter-Level Real-Time Orbit Determination Method Using Space-Borne GPS Measurements Considering IGS-RTS Data Receiving Interruption[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
Citation: ZHANG Wanwei, WANG Fuhong, GONG Xuewen, GUO Lei. A Centimeter-Level Real-Time Orbit Determination Method Using Space-Borne GPS Measurements Considering IGS-RTS Data Receiving Interruption[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1620-1626. doi: 10.13203/j.whugis20200432
  • 近年来,随着空间任务日益增多,一些低轨卫星对轨道参数的精度和实时性提出了很高的要求。例如在低轨导航增强系统中,为了地面用户能实现分米-厘米级的实时精密定位,低轨卫星自身的实时轨道精度需要达到厘米级[1-2]。在全球导航卫星系统(global navigation satellite system,GNSS)无线电掩星测量中,实时厘米级定轨是无线电掩星测量数据实时处理的前提 [3]

    目前受GNSS广播星历的轨道误差与卫星钟误差等因素的制约,实时自主定轨只能达到0.4~1.0 m的位置精度[4-6],使用国际GNSS服务组织(international GNSS service,IGS)发布的超快速预报星历进行实时定轨也只能达到0.2~0.5 m的定轨精度[6-8],离实时厘米级的精度需求仍存在一定的差距。

    自2013-04-01起,IGS通过基于互联网的RTCM网络传输协议(networked transport of radio technical commission for maritime services via internet protocol,NTRIP)向全球播发IGS实时服务(IGS real-time service,IGS-RTS)数据流,其轨道精度约为5 cm,时钟精度约为0.3 ns,时间延迟小于30 s。日本准天顶卫星系统(quasi-zenith satellite system,QZSS)、中国北斗三号导航卫星系统(the third-generation BeiDou navigation satellite system,BDS-3)通过地球同步轨道(geostationary earth orbit,GEO)卫星播发实时改正数,可应用于区域的实时精密定位。一些商业机构如NavCom、Fuguro等也通过GEO卫星向全球授权用户提供实时改正数信息。这些区域或全球实时改正数产品的播发,使星载GNSS实时精密定轨成为可能。国内外学者基于全弧段的实时改正数产品,采用实时运动学定轨模式,对哨兵3号A星(Sentinel-3A)、欧洲空间局的蜂群(SWARM)等卫星的GPS实测数据开展实时精密定轨实验,轨道精度为10~20 cm[9-10]。由于GEO卫星存在覆盖盲区,以及低轨卫星接收天线仰角和增益等因素的限制,轨道倾角接近于90°的低轨极轨卫星在极区弧段将无法接收到GEO卫星播发的实时改正数信息。文献[11]基于Fuguro公司的实时轨道钟差产品,考虑GEO卫星的覆盖盲区引起的改正数接收中断问题,采用动力学实时精密定轨方法,SWARM卫星定轨精度可达8.5 cm,其解决极区改正数缺失的方法是根据历史数据在线预报,实时定轨精度严重依赖于在线预报的精度。

    本文针对极区实时改正数接收中断情况下实现厘米级星载GPS实时定轨的关键问题开展研究。首先分析了实时改正数及其短时外推引起的星历综合误差的变化特性;然后以此为依据,在星载GPS实时精密定轨数学模型中,构建分段随机游走的伪模糊度参数随机模型,减小星历综合误差对实时精密定轨的影响,从而实现厘米级精度的实时定轨;最后使用重力恢复与气候实验(gravity recovery and climate experiment,GRACE)C卫星GPS双频实测数据和CLK93实时产品验证了本文提出的实时精密定轨方法的有效性,并进行精度分析。

    • 对于低轨卫星星载GPS双频接收机,采用GPS双频无电离层(ionospheric free,IF)组合伪距和载波相位观测值消除电离层延迟误差,且不考虑对流层延迟等误差,其观测方程可简化为:

      PIF=ρ+c·δtr-c·δts+MP+εPLIF=ρ+c·δtr-c·δts+BIF+ML+εL ]]>

      式中,PIFLIF分别为GPS双频伪距和载波相位无电离层组合观测值,单位:m;ρ为星载接收机与GPS卫星之间的几何距离;δtrδts分别为星载接收机和GPS卫星的钟差参数;c为真空中的光速;BIF为组合载波相位的模糊度;MPML分别为组合伪距和载波相位的多路径误差;εPεL分别为组合伪距和载波相位的观测噪声。

    • 使用与GPS广播星历参数中星历/时钟数据期号相符且最近的实时改正数,用户可得到GPS卫星的精密轨道和钟差。假设当前历元信号发射t时刻由GPS广播星历计算的轨道和钟差分别为rbrdcsdtbrdcs,则GPS卫星精密轨道rs和精密钟差dts的计算如下:

      rs=rbrdcs-δr ]]>
      δr=ereaecδo ]]>
      δo=[δorδoaδoc]T+[dδordδoadδoc]T(t-t0) ]]>
      ea=v/vec=(r×v)/r×ver=ea×ec ]]>
      dts=dtbrdcs+δc/c ]]>
      δc=c0+c1(t-t0)+c2(t-t0)2 ]]>

      式中,rv分别为由GPS广播星历计算的在当前历元信号发射t时刻WGS84坐标系下GPS卫星位置、速度;ereaec分别为在GPS卫星轨道径向(radial,R)、切向(along,A)、法向(cross,C)上的单位向量;δo为GPS卫星在轨道坐标系下RAC 3个方向上的改正量;δc为GPS卫星钟差改正量;c为真空中的光速;(δor,δoa,δoc)(dδor,dδoa,dδoc)分别为t0参考时刻的GPS卫星在轨道坐标系下RAC 3个方向的轨道改正值及其变化率;(c0,c1,c2)为GPS卫星钟差改正的系数。这些轨道与钟差改正数信息可直接从实时产品中得到,在短时接收中断情况下,用最近时刻的改正数来外推当前时刻的改正数。

    • 由于实时改正数产品是由IGS分析中心等机构根据精密定轨得到的轨道与钟差参数经短时间预报产生。当前,由于轨道动力学模型精度较高,短时轨道预报精度可达到厘米级。但是GPS星座使用不同类型的原子钟,卫星钟差的预报精度存在一定差异。另外,在极区以及接收中断等情况下,直接使用最近时刻的改正数来外推当前时刻的改正数,易导致改正数的误差特性发生较大的变化。由实时改正数及其短时外推改正数计算的GPS卫星轨道与钟差,与真实的卫星轨道和钟差仍存在一定差异,这种差异在GPS信号传播路径(location of signal,LOS)上引起的星历综合误差ELOS可表示为[512]

      ELOS=ρ*-cδts*-ρ-cδts=dρ-c·dδts ]]>

      式中,ρ*cδts*分别为真实的低轨卫星与GPS卫星之间的几何距离和GPS钟差;ρcδts分别为由GPS广播星历和实时改正数计算的几何距离和GPS钟差;c为真空中的光速;dρc·dδts分别为由星历误差引起的在GPS信号传播路径上的几何距离误差与GPS钟误差。

      由于ELOS参数和式(1)中模糊度参数BIF的系数一致且均未知,无法准确分离,同时,为了减少后续实时定轨滤波模型中待估参数的数量,提高实时定轨计算效率,本文将这两个参数合并为一个参数,称为伪模糊度参数[512],即A=BIF+ELOS。因此,可将观测方程改写为:

      PIF=ρ+c·δtr-c·δts+MP+εPLIF=ρ+c·δtr-c·δts+A+ML+εL ]]>
    • 低轨卫星围绕地球运转的动力学方程可表示为:

      atotal=ag+ang+aw ]]>

      式中,atotal为低轨卫星所受的总加速度;ag为由保守力(包括地球中心引力和非球形引力、N体引力、地球固体潮汐和海洋潮汐摄动力等)引起的加速度;ang为由非保守力(包括大气阻力、太阳光压力等)引起的加速度;aw为人为引入的经验加速度,用于补偿无法模型化或错误模型的微小摄动力的影响,采用一阶高斯-马尔可夫随机模型对RAC 3个方向进行动力学模型补偿[5]

    • 以GPS双频无电离层组合伪距和载波相位作为主要观测值,结合低轨卫星简化动力学模型,用扩展卡尔曼滤波估计卫星的位置、速度等状态参数,实时定轨滤波模型中的状态方程和观测方程如下:

      Xk=Φk,k-1Xk-1+Wk-1Zk=HkXk-1+Vk ]]>

      式中,XkZk分别为扩展卡尔曼滤波状态量和观测量;Φk,k-1Hk分别为状态转移矩阵和观测矩阵;Wk-1Vk为系统噪声和观测噪声。其中,待估滤波状态量选取为:

      Xk=Y6×1T2×1P5×1An×1T ]]>

      式中,Y6×1=rvk时刻低轨卫星轨道参数(包括三维位置r和速度v);T2×1=δtrdδtr为GPS接收机钟差δtr和钟漂dδtr参数;P5×1=cdcraw为动力学模型相关参数(包括大气阻力系数cd、太阳光压系数cr、RAC 3个方向的经验加速度aw);An×1=A1A2An为所观测到的n颗GPS卫星待估的伪模糊度参数。

      在参数估计时,模糊度BIF为常量,伪模糊度参数A的随机模型与ELOS参数紧密相关,需要根据ELOS的变化特性来合理建模。通过伪模糊度参数的实时准确估计,减小ELOS对实时定轨的影响,可提高实时定轨精度。

    • 本文以IGS最终精密星历和精密钟差(30 s)产品为参考,分析IGS实时改正数产品的精度及其变化特性。利用BNC(BKG Ntrip Client)软件接收了CLK93(采样率为5 s)挂载点在2019-04-08—2019-04-14,即年积日(day of year,DOY)第98~104天,的IGS-RTS数据流。使用式(2)~式(7),以30 s为间隔计算CLK93实时精密轨道与钟差结果,并与IGS精密星历和精密钟差产品计算结果作比较,在计算过程中,要将IGS精密轨道从卫星质心改正到天线相位中心。

      将两者的轨道差值转换到R、A、C方向,并统计各方向的均方根误差(root mean square error,RMSE)。为消除卫星钟差的系统偏差影响,本文选取G01卫星作为基准星,采用求二次差的方法评估实时钟差的精度[13],统计钟差二次差的均值(MEAN)、RMSE、标准差(STD)。

      图 1为DOY第98~104天CLK93实时产品中所有卫星轨道和钟差每天的精度统计。由图 1可以看出,CLK93实时产品在RAC方向的轨道精度优于5 cm,钟差精度为0.2~0.5 ns。

      图  1  CLK93产品的轨道与钟差精度

      Figure 1.  Orbit and Clock Accuracies of CLK93 Products

      图 2给出了DOY第101天每颗GPS卫星轨道和钟差的统计结果,其中G04没有轨道和钟差改正信息。由图 2可以看出,不同卫星之间的轨道误差和卫星钟同步误差都存在一定的不一致性,其中,G08卫星轨道误差超过了10 cm,G03、G18、G21等卫星钟同步误差都超过了1 ns。实时产品误差将直接影响低轨卫星实时精密定轨的轨道精度,在实时定轨算法中,需要加以建模来减小其影响。

      图  2  各GPS卫星的轨道与钟差精度

      Figure 2.  Orbit and Clock Accuracies for Each GPS Satellite

    • 考虑到GEO覆盖盲区,本文假定低轨卫星星下点纬度大于65°的区域为极区范围,对于轨道高度约500 km的极轨卫星,通过极区所需时间大约为13 min。在极区改正数接收中断时,直接使用最后可用时刻的改正数来外推得到。

      以IGS精密星历和精密钟差(30 s)产品的计算结果代替式(4)的真实轨道与钟差,使用GPS广播星历和CLK93实时产品,计算GRACE C卫星观测的所有GPS卫星的dρc·dδtsELOSdELOS图 3给出了DOY第101天连续4次过极区期间的与ELOS相关参数的变化情况,每次选取1颗GPS卫星的连续跟踪弧段为例,其中各曲线黑色部分对应为极区弧段。

      图  3  与ELOS相关参数的变化曲线

      Figure 3.  Variation Curves of Parameters Related to ELOS

      图 3可知,在非极区弧段内,dρc·dδts的曲线变化平缓且光滑,但是由变化率dELOS曲线可知,ELOS仍受到小量级的随机误差的影响。在极区弧段内,因受到外推期间广播星历误差的影响,ELOS的波动明显大于非极区,其变化特性与广播星历相近。进入极区时,ELOS变化不明显,但在出极区时刻,ELOS出现显著的台阶变化,这是由于出极区后重新接收到了新的改正数,ELOS迅速减小。

      针对ELOS在极区和非极区的不同变化特性,本文采用分段随机游走过程对伪模糊度A进行建模:

      Ak+1=Ak+ωk ]]>

      当低轨卫星进入极区弧段时,A不用重新初始化,由于直接使用最后可用时刻的改正数来外推当前时刻的改正数,相当于在广播星历上进行了常值改正,其变化特性与广播星历相近,所以过程噪声要设置为广播星历对应的数值(广播星历ELOS变化特性与随机模型设置详见文献[4-512]);在极区弧段进入非极区弧段时,因为ELOS出现显著的台阶变化,所以A要重新初始化,其过程噪声设置与实时改正数的精度有关。本文将A参数的过程噪声方差设置为(1 mm/s·dt2,其中,dt为滤波中前后历元的时间间隔。

    • 为了验证本文提出的实时精密定轨方法的有效性,根据使用GPS星历的不同,本文设计了5种实时精密定轨实验方案,方案1:仅采用GPS广播星历;方案2:采用GPS广播星历和CLK93实时产品,不考虑极区改正数接收中断情况,即全弧段使用CLK93实时产品,伪模糊度参数采用随机常量模型,过程噪声方差设为0,认为改正数没有误差;方案3:采用GPS广播星历和CLK93实时产品,考虑极区改正数接收中断情况,伪模糊度的随机模型采用本文提出的分段随机游走过程;方案4:使用的数据与方案2相同,只是伪模糊度的随机模型选取不同,这里采用随机游走过程;方案5:采用IGS精密轨道和精密钟差(30 s)产品。5种方案定轨策略及设置均相同,如表 1所示。

      表 1  实时精密定轨策略及设置

      Table 1.  Stratrgies of Real-Time Precise Orbit Determination

      模型参数 相关设置
      GPS观测数据 GPS双频伪距、载波相位无电离层组合
      EOP参数 IERS Bulletin A预报
      地球重力场模型 EGM2008 70×70
      N体引力 仅考虑日月,近似公式计算日月位置
      固体潮 简易固体潮模型
      大气阻力 改进后的Harris-Priester模型[14],固定表面积,估计大气阻力系数
      太阳光压 简单的球状模型,固定表面积,估计太阳光压系数
      经验加速度 RAC 3个方向的动力学模型补偿一阶高斯-马尔可夫随机模型
      参数估计方法 扩展卡尔曼滤波
      测量更新间隔/s 30

      根据上述5种实验方案设置,利用自主研制的实时精密定轨软件SATODS对2019年DOY第98~104天共7 d的GRACE C卫星GPS双频实测数据模拟连续实时精密定轨数据处理,计算的轨道结果与JPL精密参考轨道求差,并统计位置与速度在R、A、C方向和三维(3D)的精度。

      图 4给出了5种实验方案每天的3D位置和速度精度统计。表 2给出了5种实验方案连续7天实时定轨的位置与速度精度。除了方案1仅采用广播星历外,方案2~4都使用CLK93实时产品后,实时定轨精度从30.9 cm提高到优于10 cm,实现了厘米级实时定轨。其中,方案4全弧段使用CLK93实时产品且考虑伪模糊度随机模型与方案5采用IGS精密产品的实时定轨精度相当。

      图  4  5种方案的每天位置和速度精度统计

      Figure 4.  Daily Position and Velocity Accuracy Statistics of Five Schemes

      表 2  GRACE C卫星的实时精密定轨结果统计

      Table 2.  Accuracy Statistics of Real-Time Orbit Determination for GRACE C

      方案 位置精度/cm 速度精度/(mm·s-1
      R方向 A方向 C方向 3D R方向 A方向 C方向 3D
      1 10.3 26.5 12.3 30.9 0.298 0.133 0.161 0.363
      2 4.74 5.32 4.00 8.18 0.144 0.097 0.099 0.201
      3 3.74 4.85 3.45 7.04 0.139 0.096 0.103 0.197
      4 3.12 3.68 2.72 5.55 0.134 0.094 0.099 0.191
      5 2.85 2.85 2.54 4.78 0.131 0.093 0.099 0.188

      图 5给出了DOY第101天方案2、3、5分别与方案4的3D位置误差对比。由图 5可以看出,全弧段都使用CLK93实时产品时,方案2的大部分弧段轨道误差要大于方案4,说明不同GPS卫星不同时段的实时改正数精度并不均匀,对实时定轨精度的影响较大。方案4采用随机游走过程的伪模糊度模型,大幅减弱了实时改正数精度不均匀的影响,定轨精度接近于方案5。与方案4相比,方案3考虑了极区实时改正数接收中断情况,采用分段随机游走过程来描述伪模糊度参数,可以消除大部分极区弧段的改正数外推误差的影响,定轨精度可以达到7.04 cm,但部分弧段较差,甚至达到了20 cm。

      图  5  4种方案的三维位置误差比较

      Figure 5.  3D Position Error Comparison of Schemes 2, 3, 4 and 5

      图 6给出了方案3与方案4部分精度较差弧段(对应图 5中虚线框内所示范围)在RAC方向和3D位置误差对比。由图 6可以看出,方案3在R、C方向受极区改正数外推的影响较小,而A方向受影响较大,位置误差曲线的波峰出现在LEO卫星出极区弧段,可能与部分极区弧段改正数的外推值不能准确消除广播星历的误差有关,极区弧段伪模糊度的随机模型切换到广播星历的随机

      图  6  方案3和方案4定轨位置误差对比

      Figure 6.  Comparison of Orbit Errors for Schemes 3 and 4

      模型参数设置,如果该弧段广播星历的误差不显著,定轨误差会得到较好保持,如果该弧段广播星历的误差较大,定轨误差将会出现逐渐增加的现象,定轨精度将会退化到方案1的分米级。但是在出极区后,由于重新初始化了伪模糊度参数,随着新改正数的使用,伪模糊度参数将逐渐重新收敛,非极区弧段的定轨误差又达到了厘米级。

    • 本文利用IGS播发的GPS卫星轨道与钟差改正数实时数据流,提出用分段随机游走过程来构建伪模糊度参数随机模型,减小星历综合误差对实时精密定轨精度的影响,从而实现厘米级星载GPS

      实时定轨。基于CLK93实时产品,设计了多种实时定轨实验方案,对GRACE C卫星GPS双频实测数据模拟在轨实时精密定轨处理,得出以下结论:(1)不考虑实时改正数中断的情况下,基于伪模糊度随机模型的实时精密定轨算法,定轨精度优于6 cm(3D RMSE),接近于IGS精密轨道与钟差产品的实时定轨精度;(2)针对极区弧段实时改正数中断情况下,基于本文提出的伪模糊度分段随机游走模型的实时精密定轨算法,定轨精度约为7 cm(3D RMSE),大幅减小实时改正数外推误差对定轨精度的影响。

      综上所述,本文提出基于IGS-RTS的厘米级星载GPS实时定轨的数学模型及其处理策略具有可行性和有效性,较好地解决了极区实时改正数中断时外推误差对实时定轨精度的影响,对于部分极区弧段实时定轨精度较差的影响因素,还需要进一步的研究。

参考文献 (14)

目录

    /

    返回文章
    返回