袁俊军1,2,3, 李凯1,3, 唐成盼1,3, 周善石1,3, 胡小工1,3, 曹建峰4
1. 中国科学院上海天文台, 上海 200030;
2. 中国科学院大学, 北京 100049;
3. 上海市空间导航与定位技术重点实验室, 上海 200030;
4. 北京航天飞行控制中心, 北京 100094
基金项目:国家自然科学基金(12103077;12173072;41804030;41874039)
摘要:LEO卫星精密轨道预报是LEO导航增强系统中重要的技术环节之一, 本文使用多种算法来实现不同任务需求下的轨道预报。对于在地面处理系统实现的LEO轨道预报, 算法1采用定轨预报同时处理的策略, 算法2将离散轨道点进行动力学拟合再进行积分外推。GRACE-C卫星预报5、10、15min的URE平均精度分别为5.25、5.67、6.25cm; HY2A卫星为7.83、8.69、9.66cm; SWARM-A卫星为8.88、9.22、9.63cm; SWARM-B卫星为8.49、8.98、9.63cm。对于计算条件受限的LEO星上轨道预报, 本文利用单个轨道点及简单动力学模型进行轨道积分外推的算法。该算法主要考虑地球中心引力及非球形引力摄动, 因此地球重力场阶次对轨道预报精度产生较大影响。平均高度为500km的LEO卫星选取60阶重力场, 高度为1000km的LEO卫星选取30阶重力场, 可实现预报10min轨道优于10cm的预报精度。
关键词:LEO导航增强 轨道预报 动力学拟合 轨道积分
YUAN Junjun, LI Kai, TANG Chengpan, et al. Accuracy analysis of LEO satellites orbit prediction for precise position service[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(5): 640-647. DOI: 10.11947/j.AGCS.2022.20210473
(3)
表 1 LEO星下轨道定轨策略Tab. 1 LEO Orbit determination strategy on ground
模型 | 参数 | 算法1 | 算法2 |
测量模型 | 观测量 | 伪距、相位 | 坐标 |
定轨弧长/h | 12 | 12 | |
预报时长/min | 5、10、15 | 5、10、15 | |
权重 | 相位: 伪距为100∶1 | 等权 | |
GPS产品 | CODE产品 | - | |
截止高度角/(°) | 5 | - | |
动力学模型 | 重力场 | EIGEN_GL04C, 120×120 | EIGEN_GL04C, 120×120 |
N体 | JPL DE405 | JPL DE405 | |
相对论 | IERS 2010 | IERS 2010 | |
固体潮 | IERS 2010 | IERS 2010 | |
海潮 | FES 2004 | FES 2004 | |
太阳光压 | 2 h估计一组 | 2 h估计一组 | |
大气阻力 | 2 h估计一组 | 2 h估计一组 | |
经验力 | 详细讨论 | 详细讨论 |
表选项
图 1 星下轨道预报流程Fig. 1 The process of LEO orbit prediction on ground |
图选项 |
图 2 不同经验力分段策略下LEO卫星预报5 min轨道预报精度Fig. 2 LEO orbit prediction accuracy for 5 min under different empirical force segmentation strategies |
图选项 |
表 2 不同经验力分段策略下LEO轨道预报5 min精度统计Tab. 2 Statistical accuracy of LEO orbit prediction for 5 min under different empirical force segmentation strategies
卫星 | 1个轨道周期 | 2个轨道周期 | 3个轨道周期 | |||||
算法1 | 算法2 | 算法1 | 算法2 | 算法1 | 算法2 | |||
GRCC | 7.07 | 7.07 | 5.25 | 5.24 | 5.88 | 5.87 | ||
HY2A | 8.76 | 8.75 | 8.01 | 8.02 | 7.84 | 7.81 | ||
SWRA | 12.08 | 12.08 | 8.87 | 8.88 | 9.00 | 9.01 | ||
SWRB | 12.48 | 12.48 | 8.49 | 8.48 | 8.62 | 8.62 |
表选项
表 3 LEO星下轨道预报精度统计Tab. 3 Statistical accuracy of LEO orbit prediction on ground
卫星 | 策略 | R方向 | T方向 | N方向 | URE | |||||||||||
5 min | 10 min | 15 min | 5 min | 10 min | 15 min | 5 min | 10 min | 15 min | 5 min | 10 min | 15 min | |||||
GRCC | 算法1 | 2.71 | 2.81 | 2.93 | 3.47 | 3.99 | 4.66 | 1.62 | 1.67 | 1.77 | 5.25 | 5.67 | 6.26 | |||
算法2 | 2.71 | 2.81 | 2.93 | 3.47 | 3.99 | 4.66 | 1.62 | 1.67 | 1.77 | 5.24 | 5.66 | 6.24 | ||||
HY2A | 算法1 | 4.45 | 4.85 | 5.24 | 5.17 | 6.17 | 7.29 | 2.02 | 1.93 | 1.86 | 7.84 | 8.70 | 9.68 | |||
算法2 | 4.43 | 4.82 | 5.20 | 5.17 | 6.16 | 7.28 | 1.99 | 1.90 | 1.84 | 7.81 | 8.67 | 9.64 | ||||
SWRA | 算法1 | 2.79 | 3.07 | 3.35 | 7.09 | 7.37 | 7.77 | 2.52 | 2.55 | 2.61 | 8.87 | 9.21 | 9.61 | |||
算法2 | 2.77 | 3.06 | 3.33 | 7.09 | 7.39 | 7.80 | 2.55 | 2.57 | 2.63 | 8.88 | 9.23 | 9.64 | ||||
SWRB | 算法1 | 3.65 | 4.13 | 4.55 | 5.55 | 6.00 | 6.72 | 3.34 | 3.33 | 3.33 | 8.49 | 8.98 | 9.63 | |||
算法2 | 3.65 | 4.13 | 4.55 | 5.55 | 6.00 | 6.72 | 3.34 | 3.33 | 3.33 | 8.48 | 8.97 | 9.62 |
表选项
表 4 LEO卫星摄动力量级分析Tab. 4 The analysis of LEO satellites perturbation force level
摄动力 | GRCC | HY2A | SWRA | SWRB |
非球形摄动 | 1.83E-03 | 1.55E-03 | 1.86E-03 | 1.83E-03 |
N体摄动 | 1.69E-04 | 2.83E-04 | 2.29E-04 | 2.58E-04 |
相对论 | 1.94E-09 | 1.81E-09 | 1.95E-09 | 1.93E-09 |
大气摄动 | 7.91E-10 | 2.14E-10 | 2.68E-09 | 8.99E-10 |
光压摄动 | 1.45E-09 | 9.44E-09 | 1.05E-09 | 6.70E-10 |
经验力 | 2.74E-09 | 3.27E-09 | 2.71E-09 | 4.37E-09 |
表选项
图 3 不同重力场阶数下LEO卫星星上轨道预报精度Fig. 3 The LEO orbit prediction accuracy on board with different Earth gravity field's orders |
图选项 |
表 5 LEO星上轨道预报精度统计Tab. 5 Statistical accuracy of LEO orbit prediction on board
卫星 | 阶次 | R方向 | T方向 | N方向 | URE | |||||||
5 min | 10 min | 5 min | 10 min | 5 min | 10 min | 5 min | 10 min | |||||
GRCC | 10阶 | 0.27 | 0.83 | 0.22 | 0.51 | 0.23 | 0.66 | 0.47 | 1.31 | |||
30阶 | 0.04 | 0.08 | 0.03 | 0.05 | 0.04 | 0.09 | 0.07 | 0.15 | ||||
60阶 | 0.02 | 0.04 | 0.02 | 0.03 | 0.01 | 0.04 | 0.03 | 0.07 | ||||
HY2A | 10阶 | 0.11 | 0.33 | 0.08 | 0.16 | 0.08 | 0.28 | 0.18 | 0.51 | |||
30阶 | 0.02 | 0.03 | 0.02 | 0.04 | 0.02 | 0.04 | 0.04 | 0.08 | ||||
60阶 | 0.01 | 0.03 | 0.02 | 0.04 | 0.02 | 0.04 | 0.04 | 0.07 | ||||
SWRA | 10阶 | 0.32 | 0.98 | 0.19 | 0.51 | 0.30 | 0.82 | 0.54 | 1.53 | |||
30阶 | 0.05 | 0.12 | 0.04 | 0.08 | 0.05 | 0.14 | 0.09 | 0.22 | ||||
60阶 | 0.01 | 0.03 | 0.04 | 0.05 | 0.02 | 0.04 | 0.05 | 0.08 | ||||
SWRB | 10阶 | 0.30 | 0.84 | 0.18 | 0.43 | 0.19 | 0.68 | 0.45 | 1.29 | |||
30阶 | 0.03 | 0.10 | 0.05 | 0.08 | 0.03 | 0.08 | 0.07 | 0.17 | ||||
60阶 | 0.02 | 0.04 | 0.04 | 0.06 | 0.02 | 0.03 | 0.05 | 0.09 |
往期推荐
资讯
○《国土空间规划城市时空大数据应用基本规定》(征求意见稿)公开征求意见
会议
《测绘学报》
○ 测绘学报 | 张烁:嫦娥五号探测器月面采样封装任务的定位精度
○《测绘学报》同济专刊 | 王权:陆地卫星遥感监测体系及应用前景
○《测绘学报》同济专刊 | 孙伟伟:空谱协同多尺度顶点成分分析的高光谱影像端元提取
○《测绘学报》同济专刊 | 乔刚:全球气候变化背景下极地冰盖关键参数遥感观测验证
《测绘通报》
《北京测绘》
《测绘科学技术学报》
《地球信息科学学报》
《测绘工程》
《中国空间科学技术》
《卫星应用》
《Journal of Geodesy and Geoinformation Science》
○《测绘学报(英文版)》专刊征稿 | 地图学与地球空间信息教育:理论与实践
○ 《测绘学报(英文版)》专刊征稿 | 用于三维地理信息的摄影测量和计算机视觉
《Satellite Navigation》
○ 杨飞博士:GNSS天顶对流层精化模型的构建与分析| SANA佳文速递
○ 2022征文| SatNav“普适定位、室内导航与基于位置服务”专题
○ 李昕博士:多频多星座PPP-RTK原理及性能| SANA佳文速递
《自然资源遥感》
○ 《自然资源遥感》征稿:“海岸带空间资源及生态健康遥感监测”专栏
《Journal of Geovisualization and Spatial Analysis》
○《Journal of Geovisualization and Spatial Analysis》入驻“智绘科服”融媒体平台!
○ 高被引论文推荐 | Journal of Geovisualization and Spatial Analysis
○ JGSA论文推荐 | 地理信息科学研究在过去几十年中都在关注什么?
《全球定位系统》
○《全球定位系统》入选《世界期刊影响力指数(WJCI)报告(2021科技版)》
《导航定位与授时》
Send to Author