留言板

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

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

定位归算在嫦娥四号中继星变轨阶段的应用

张宇 郭丽 王广利 李金岭 黄勇 李培佳 郑为民 王文彬 舒逢春 刘庆会

张宇, 郭丽, 王广利, 李金岭, 黄勇, 李培佳, 郑为民, 王文彬, 舒逢春, 刘庆会. 定位归算在嫦娥四号中继星变轨阶段的应用[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
引用本文: 张宇, 郭丽, 王广利, 李金岭, 黄勇, 李培佳, 郑为民, 王文彬, 舒逢春, 刘庆会. 定位归算在嫦娥四号中继星变轨阶段的应用[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
ZHANG Yu, GUO Li, WANG Guangli, LI Jinling, HUANG Yong, LI Peijia, ZHENG Weimin, WANG Wenbin, SHU Fengchun, LIU Qinghui. Application of Positioning Reduction in the Orbital Transfer Stage of the Chang'E-4 Relay Satellite[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
Citation: ZHANG Yu, GUO Li, WANG Guangli, LI Jinling, HUANG Yong, LI Peijia, ZHENG Weimin, WANG Wenbin, SHU Fengchun, LIU Qinghui. Application of Positioning Reduction in the Orbital Transfer Stage of the Chang'E-4 Relay Satellite[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365

定位归算在嫦娥四号中继星变轨阶段的应用

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

国家自然科学基金 11873076

国家自然科学基金 11473056

详细信息
    作者简介:

    张宇, 硕士生, 主要研究方向为VLBI定位。zhangyu@shao.ac.cn

    通讯作者: 郭丽, 博士, 副研究员。lguo@shao.ac.cn
  • 中图分类号: P228;P129

Application of Positioning Reduction in the Orbital Transfer Stage of the Chang'E-4 Relay Satellite

Funds: 

The National Natural Science Foundation of China 11873076

The National Natural Science Foundation of China 11473056

More Information
    Author Bio:

    ZHANG Yu, postgraduate, specializes in VLBI positioning. E-mail:zhangyu@shao.ac.cn

    Corresponding author: GUO Li, PhD, associate researcher. E-mail:lguo@shao.ac.cn
  • 摘要: 2018-06嫦娥四号中继星发射任务圆满完成,为2018-12的嫦娥四号月球背面降落任务打下良好的基础,其中甚长基线干涉测量技术(very long baseline interferometry,VLBI)测轨分系统在卫星的定位、变轨等方面起到了重要的作用。统计了嫦娥四号中继星任务实时阶段VLBI观测的时延、时延率和定位结果的精度,讨论了定位归算在嫦娥四号中继星轨控弧段实时轨迹监测中的应用。在中继星进入Halo轨道的过程中,能够实时监测到定位归算结果给出的轨道根数变化,以实时监测变轨情况,并通过事后与定轨结果比对,判断定位结果的正确性与准确性,从而为后续飞行任务提供了一种新的Halo轨道实时监测的方法。
  • 图  1  时延数据定轨后残差统计

    Figure  1.  Residual Statistics of Delay Data After Orbit Determination

    图  2  时延率数据定轨后残差统计

    Figure  2.  Residual Statistics of Delay Rate Data After Orbit Determination

    图  3  观测代码s8525a定轨与定位归算地心坐标结果差值

    Figure  3.  Difference Between Orbit Determination and Positioning Reduction of the Geocentric Coordinates Result of Observation Code s8525a

    图  4  观测代码s8525a定轨与定位归算速度结果差值

    Figure  4.  Difference Between Orbit Determination and Positioning Reduction of the Velocity Result of Observation Code s8525a

    图  5  观测代码s8521a轨道半长轴变化图像

    Figure  5.  Orbital Semi-major Axis Change Image of Observation Code s8521a

    图  6  观测代码s8521a轨道倾角变化图像

    Figure  6.  Orbital Inclination Change Image of Observation Code s8521a

    图  7  观测代码s8525a轨道偏心率变化图像

    Figure  7.  Orbital Eccentricity Change Image of Observation Code s8525a

    图  8  观测代码s8607a轨道偏心率变化图像

    Figure  8.  Orbital Eccentricity Change Image of Observation Code s8607a

    图  9  观测代码s8614a轨道半长轴变化图像

    Figure  9.  Orbital Semi-major Axis Change Image of Observation Code s8614a

    图  10  观测代码s8614a轨道偏心率变化图像

    Figure  10.  Orbital Eccentricity Change Image of Observation Code s8614a

    图  11  观测代码s8614a轨道倾角变化图像

    Figure  11.  Orbital Inclination Change Image of Observation Code s8614a

    图  12  观测代码s8614a实时定位归算与事后定轨结果的坐标差值

    Figure  12.  Coordinate Difference of Observation Code s8614a Between Real-Time Positioning Reduction and Post-Time Orbit Determination

    图  13  观测代码s8614a实时定位归算与事后定轨结果的速度差值

    Figure  13.  Velocity Difference of Observation Code s8614a Between Real-Time Positioning Reduction and Post-Time Orbit Determination

    表  1  CE-4中继星任务与CE-2、CE-3、CE-5T1奔月段定轨残差均方根比较

    Table  1.   Root Mean Square Comparison of Orbit Determination Residuals of CE-4 Relay Satellite Mission, CE-2, CE-3 and CE-5T1 Lunar Missions

    卫星名称 信号波段 时延/ns 时延率/
    (ps•s-1
    CE-2 实时S波段 4.93 0.64
    事后S波段 4.04 0.63
    事后带宽综合 1.28 0.45
    CE-3 实时X波段 0.78 0.93
    事后X波段 0.85 0.90
    CE-5T1 实时X波段 1.43 0.65
    实时S波段 3.44 1.13
    事后X波段 1.41 0.59
    事后S波段 2.24 1.07
    CE-4ZJ 实时S波段 0.91 0.67
    下载: 导出CSV

    表  2  定位计算的误差统计

    Table  2.   Error Statistics of Positioning Calculation

    观测代码 总点数 赤经误差/as 赤纬误差/as
    s8521a 2 827 0.007 0.009
    s8522a 3 170 0.004 0.007
    s8523a 2 662 0.006 0.007
    s8524a 1 417 0.012 0.014
    s8525a 3 495 0.014 0.013
    s8526a 3 112 0.004 0.007
    s8527a 2 652 0.003 0.003
    s8528a 2 640 0.013 0.013
    s8529a 2 683 0.005 0.007
    s8530a 2 359 0.015 0.024
    s8531a 2 337 0.007 0.017
    s8601a 2 412 0.019 0.017
    s8602a 1 842 0.011 0.017
    s8603a 1 994 0.006 0.010
    s8604a 1 958 0.002 0.005
    s8605a 1 811 0.002 0.005
    s8606a 1 804 0.006 0.007
    s8607a 2 222 0.007 0.008
    s8608a 1 811 0.071 0.012
    s8609a 1 767 0.004 0.008
    s8610a 1 811 0.004 0.007
    s8612a 925 0.082 0.041
    s8613a 1 535 0.052 0.050
    s8614a 2 171 0.007 0.014
    s8615a 878 0.006 0.010
    s8620a 913 0.008 0.014
    下载: 导出CSV

    表  3  CE-4任务实时阶段主要轨道机动事件

    Table  3.   Main Orbital Maneuver Events of CE-4 Mission Real-Time Phase

    观测代码 时间 事件
    s8521a 2018-05-21 14:54:20—14:54:32 中途修正
    s8525a 2018-05-25 13:32:10—13:46:57 近月制动
    s8527a 2018-05-27 13:46:10—16:46:21 中途修正
    s8607a 2018-06-08 00:00:00—00:00:15 捕获机动
    s8614a 2018-06-14 03:00:00—03:06:01 捕获机动
    下载: 导出CSV

    表  4  第1次中途修正变轨前后轨道根数记录

    Table  4.   Oribital Element Records of First Track Correction Before and After Orbital Transfer

    时间 半长轴/km 偏心率 轨道倾角/(°)
    14:54:18201 4660.96728.496
    14:54:20201 4680.96728.541
    14:54:22201 4780.96728.505
    14:54:24201 3730.96828.466
    14:54:26201 3430.96728.607
    14:54:28201 2540.96828.676
    14:54:30201 2410.96728.799
    14:54:32201 0830.96828.793
    14:54:34201 0290.96828.824
    下载: 导出CSV

    表  5  观测代码s8525a变轨前后偏心率记录

    Table  5.   Eccentricity Record of Observation Code s8525a Before and After Orbital Transfer

    时间 偏心率
    13:32:051.303
    13:33:051.288
    13:34:051.310
    13:35:051.224
    13:36:051.201
    13:37:051.176
    13:38:051.151
    13:39:051.130
    13:40:051.108
    13:41:051.087
    13:42:051.066
    13:43:051.046
    13:44:051.025
    13:45:051.037
    13:46:050.996
    13:47:050.971
    13:48:050.966
    13:49:050.966
    下载: 导出CSV

    表  6  观测代码s8614a变轨前后轨道根数记录

    Table  6.   Orbital Element Records of Observation Code s8614a Before and After Orbital Transfer

    时间 半长轴/ km 偏心率 轨道倾角/(°)
    03:00:1037 0850.91477.126
    03:00:4036 8860.92475.295
    03:01:0936 6380.93772.287
    03:01:3936 5330.94270.621
    03:02:0836 3920.95067.767
    03:02:3836 2330.95862.265
    03:03:0736 1460.96360.575
    03:03:3736 0480.96856.106
    03:04:0635 9490.97352.139
    03:04:3635 9000.97545.309
    03:05:0535 8510.97838.921
    03:05:3535 8130.98031.840
    03:06:0435 7980.98122.643
    03:06:3435 7840.98221.723
    下载: 导出CSV
  • [1] 郭丽.基于VLBI跟踪观测的空间飞行器瞬时状态参量归算[D].上海: 中国科学院上海天文台, 2007 http://d.wanfangdata.com.cn/Thesis/Y1620865

    Guo Li. Reduction of the Instantaneous State Vectors of Spacecraft Based on VLBI Tracking Data[D]. Shanghai: Shanghai Astronomical Observatory, Chinese Academy of Sciences, 2007 http://d.wanfangdata.com.cn/Thesis/Y1620865
    [2] 李金岭, 刘鹂, 郑为民, 等.定位归算在嫦娥二号任务实时阶段的应用[J].中国科学:物理学力学天文学, 2011, 41(7):889-895 http://www.cnki.com.cn/Article/CJFDTotal-JGXK201107012.htm

    Li Jinling, Liu Li, Zheng Weimin, et al. The Application of Positioning Reduction in the Real-Time Stage of the Chang'E-2 Project[J]. Scientia Sinica:Physical Mechanics & Astronomy, 2011, 41(7):889-895 http://www.cnki.com.cn/Article/CJFDTotal-JGXK201107012.htm
    [3] 平劲松, 王明远, 史弦, 等.嫦娥一号绕月探测器轨道投入过程实时监测判定的原理与技术实现[J].空间科学学报, 2011, 31(3):330-337 http://d.old.wanfangdata.com.cn/Periodical/kjkxxb201103008

    Ping Jinsong, Wang Mingyuan, Shi Xian, et al. Real Time Monitoring the Orbit Insertion of Chang'E-1 Lunar Mission[J]. China Journal of Space Science, 2011, 31(3):330-337 http://d.old.wanfangdata.com.cn/Periodical/kjkxxb201103008
    [4] 杜兰, 张中凯, 于亮, 等.圆型限制性三体问题下Halo-绕月轨道的星间测距定轨[J].测绘学报, 2013, 42(2):184-190 http://d.old.wanfangdata.com.cn/Periodical/chxb201302004

    Du Lan, Zhang Zhongkai, Yu Liang, et al. SST Orbit Determination of Halo-LMO Constellation in CRTBP[J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(2):184-190 http://d.old.wanfangdata.com.cn/Periodical/chxb201302004
    [5] 高珊, 周文艳, 梁伟光, 等.地月拉格朗日L2点中继星轨道分析与设计[J].深空探测学报, 2017, 4(2):122-129 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=sktcxb201702004

    Gao Shan, Zhou Wenyan, Liang Weiguang, et al. Trajectory Analysis and Design for Relay Satellite Using Lagrange L2 Point of Earth-Moon System[J]. Journal of Deep Space Exploration, 2017, 4(2):122-129 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=sktcxb201702004
    [6] 孙超, 唐玉华, 李翔宇, 等.地-月L2点中继星月球近旁转移轨道设计[J].深空探测学报, 2017, 4(3):264-269, 275 http://www.cqvip.com/QK/71905X/201703/7000242832.html

    Sun Chao, Tang Yuhua, Li Xiangyu, et al. Design of Earth-Moon L2 Halo Orbit Transfer Trajectory for Relay Satrllites Using Lunar Flybys[J]. Journal of Deep Space Exploration, 2017, 4(3):264-269, 275 http://www.cqvip.com/QK/71905X/201703/7000242832.html
    [7] Richardson D L. Analytical Construction of Periodic Orbits About the Collinear Points[J]. Celestial Mechanics, 1980, 22(3):241-253 doi:  10.1007/BF01229511
  • [1] 魏二虎, 任晓斌, 刘经南, 李连艳, 聂桂根, 李岩林.  利用VLBI观测量对月球天平动参数解算及着陆器定位和测速的改进 . 武汉大学学报 ● 信息科学版, 2022, 47(4): 483-491. doi: 10.13203/j.whugis20200311
    [2] 李金岭, 赵达新, 柳聪, 张会, 虞林峰, 赵融冰, 王锦清, 范庆元.  佘山13 m口径射电望远镜指向扫描数据解析 . 武汉大学学报 ● 信息科学版, 2020, 45(2): 159-166. doi: 10.13203/j.whugis20180415
    [3] 张志斌, 王广利, 宋世泽, 张阿丽, 王宏.  归心测量中轮轨式VLBI天线的水平钢轨形变建模 . 武汉大学学报 ● 信息科学版, 2019, 44(10): 1498-1504. doi: 10.13203/j.whugis20180050
    [4] 李斐, 邵先远, 曲春凯, 鄢建国, 郝卫峰, 叶茂, 杨轩, 金炜桐.  利用2006—2015年VLBI数据进行地球定向参数解算与分析 . 武汉大学学报 ● 信息科学版, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
    [5] 孙和平, 刘清超, 崔小明, 徐建桥.  利用重力和VLBI技术检测地球液态核的动力学效应 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 2058-2063. doi: 10.13203/j.whugis20180110
    [6] 李金岭, 黄飞, 孙中苗.  Ka波段天测与测地VLBI观测的利弊分析 . 武汉大学学报 ● 信息科学版, 2017, 42(2): 257-262. doi: 10.13203/j.whugis20140723
    [7] 魏二虎, 汤深权, 刘建栋, 董翠军.  基于联邦卡尔曼滤波的月球车联合定位方法 . 武汉大学学报 ● 信息科学版, 2017, 42(10): 1431-1437. doi: 10.13203/j.whugis20160137
    [8] 魏二虎, 刘文杰, WEI Jianan, 金双根, 刘经南.  VLBI和GPS观测联合解算地球自转参数和日长变化 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 66-71,92. doi: 10.13203/j.whugis20130435
    [9] 刘鹂, 李金岭.  中国VLBI网测站参数解析 . 武汉大学学报 ● 信息科学版, 2014, 39(3): 262-266. doi: 10.13203/j.whugis20120005
    [10] 李金岭, 张津维.  利用GPS测量监测VLBI天线参考点的仿真分析 . 武汉大学学报 ● 信息科学版, 2013, 38(12): 1387-1391.
    [11] 张志斌, 王广利, 刘祥, 唐正宏.  中国VLBI网观测地球定向参数能力分析 . 武汉大学学报 ● 信息科学版, 2013, 38(8): 911-915.
    [12] 李金岭, 郭丽, 孙中苗.  关于脉冲星参数的我国VLBI网测量 . 武汉大学学报 ● 信息科学版, 2013, 38(4): 386-389.
    [13] 魏二虎, 张奇.  利用1985~2010年VLBI数据进行日长变化分析 . 武汉大学学报 ● 信息科学版, 2012, 37(9): 1020-1023.
    [14] 李金岭, 魏二虎, 孙中苗, 钱志瀚.  关于我国天测与测地VLBI网络未来建设的讨论 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 670-673.
    [15] 魏二虎, 李雪川, 易慧, 刘经南.  利用2008年VLBI观测数据进行EOP参数的解算与分析 . 武汉大学学报 ● 信息科学版, 2010, 35(8): 988-990.
    [16] 魏二虎, 刘经南, 潘培鲸.  近三年VLBI观测的数据处理及分析 . 武汉大学学报 ● 信息科学版, 2008, 33(12): 1275-1278.
    [17] 朱新慧, 孙付平, 秦勇.  组合GPS和VLBI数据建立板块运动模型 . 武汉大学学报 ● 信息科学版, 2005, 30(7): 604-608.
    [18] 张捍卫, 郑勇, 杜兰.  河外射电源地面VLBI观测的相对论时延模型 . 武汉大学学报 ● 信息科学版, 2003, 28(5): 585-588.
    [19] 钟胜坚, 舒逢春, 李金岭, 吴徳, 何旋, 童锋贤.  利用中国VGOS站开展国际联测的EOP精度仿真分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210204
    [20] 郭丽, 李金岭, 王广利, 黄逸丹, 王文彬, 郑为民.  嫦娥五号探测器月面定位与精度分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): -. doi: 10.13203/j.whugis20210226
  • 加载中
图(13) / 表(6)
计量
  • 文章访问数:  677
  • HTML全文浏览量:  127
  • PDF下载量:  55
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-08-28
  • 刊出日期:  2020-02-05

定位归算在嫦娥四号中继星变轨阶段的应用

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

    国家自然科学基金 11873076

    国家自然科学基金 11473056

    作者简介:

    张宇, 硕士生, 主要研究方向为VLBI定位。zhangyu@shao.ac.cn

    通讯作者: 郭丽, 博士, 副研究员。lguo@shao.ac.cn
  • 中图分类号: P228;P129

摘要: 2018-06嫦娥四号中继星发射任务圆满完成,为2018-12的嫦娥四号月球背面降落任务打下良好的基础,其中甚长基线干涉测量技术(very long baseline interferometry,VLBI)测轨分系统在卫星的定位、变轨等方面起到了重要的作用。统计了嫦娥四号中继星任务实时阶段VLBI观测的时延、时延率和定位结果的精度,讨论了定位归算在嫦娥四号中继星轨控弧段实时轨迹监测中的应用。在中继星进入Halo轨道的过程中,能够实时监测到定位归算结果给出的轨道根数变化,以实时监测变轨情况,并通过事后与定轨结果比对,判断定位结果的正确性与准确性,从而为后续飞行任务提供了一种新的Halo轨道实时监测的方法。

English Abstract

张宇, 郭丽, 王广利, 李金岭, 黄勇, 李培佳, 郑为民, 王文彬, 舒逢春, 刘庆会. 定位归算在嫦娥四号中继星变轨阶段的应用[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
引用本文: 张宇, 郭丽, 王广利, 李金岭, 黄勇, 李培佳, 郑为民, 王文彬, 舒逢春, 刘庆会. 定位归算在嫦娥四号中继星变轨阶段的应用[J]. 武汉大学学报 ● 信息科学版, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
ZHANG Yu, GUO Li, WANG Guangli, LI Jinling, HUANG Yong, LI Peijia, ZHENG Weimin, WANG Wenbin, SHU Fengchun, LIU Qinghui. Application of Positioning Reduction in the Orbital Transfer Stage of the Chang'E-4 Relay Satellite[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
Citation: ZHANG Yu, GUO Li, WANG Guangli, LI Jinling, HUANG Yong, LI Peijia, ZHENG Weimin, WANG Wenbin, SHU Fengchun, LIU Qinghui. Application of Positioning Reduction in the Orbital Transfer Stage of the Chang'E-4 Relay Satellite[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
  • 探月工程嫦娥四号(Chang’E-4,CE-4)包括中继星和着巡组合体(着巡体)。着巡体在月球背面着陆,开展巡视探测任务。中继星在地月L2点绕飞,建立月地通信链路。CE-4中继星“鹊桥”于2018-05-21 05:28利用CZ-4C火箭在西昌发射,着巡体于2018-12-08利用CZ-3B火箭在西昌发射。CE-4的目标是实现人类首次月球背面软着陆并进行探测。由于月球自转与地球公转周期相同,月球背面始终不会与地球相对,因此需要发射一颗中继星承担月球背面与地球之间信号传输中转的作用。CE-4中继星经过发射段、地月转移段、月球-地月L2点转移段和地月L2点工作段4个飞行时段,最终围绕在地月L2点的Halo轨道上运行,执行地月之间信息中转传输的任务。在整个飞行过程中,中继星经过中途修正、近月制动、捕获机动等多次变轨后,于2018-06-14成功进入围绕地月L2点的Halo轨道。

    甚长基线干涉测量技术(very long baseline interferometry,VLBI)测轨分系统是探月工程测控系统的组成部分,需要开展中继星2万km远的地月转移段、月球-地月L2点转移段和地月L2点使命轨道段的测量工作。VLBI测轨分系统将综合利用VLBI测量结果、测距测速数据实现对中继星各飞行段的定轨定位和轨道预报。定位归算为VLBI测轨分系统的软件配置项之一[1]。本文统计分析了CE-4中继星多次变轨前后VLBI时延及时延率观测量和定位归算的结果,并将定位归算结果转化为轨道根数,通过轨道根数的变化情况实时判断中继星的变轨情况,另外分析了轨道根数较传统通过视向速度判断卫星变轨情况方法的优势之处。

    • CE-4中继星任务以S频段数据为主,并联合测距测速和VLBI两种测控手段对探测器进行轨道确定。CE-4中继星的VLBI观测和嫦娥三号(CE-3)一样,利用的是差分单程定位技术。在实时模式下,对卫星和其附近位置精确已知的射电源交替观测,修正相关系统误差,从而获得高精度的观测量。

      图 1图 2分别为时延和时延率数据定轨后残差统计结果。从图 1图 2可知,此次CE-4中继星的时延和时延率数据1 d弧长定轨后,残差平均值分别在0.91 ns和0.67 ps/s的水平。

      图  1  时延数据定轨后残差统计

      Figure 1.  Residual Statistics of Delay Data After Orbit Determination

      图  2  时延率数据定轨后残差统计

      Figure 2.  Residual Statistics of Delay Rate Data After Orbit Determination

      表 1为CE-4中继星任务与CE-2、CE-3、CE-5T1奔月段定轨残差均方根的比较。从表 1可以看到,CE-4中继星任务(CE-4ZJ)的VLBI时延测量数据定轨后残差相较CE-3的X频段数据精度略有下降,但精度高于嫦娥二号(CE-2)和嫦娥五号实验(CE-5T1)任务。

      表 1  CE-4中继星任务与CE-2、CE-3、CE-5T1奔月段定轨残差均方根比较

      Table 1.  Root Mean Square Comparison of Orbit Determination Residuals of CE-4 Relay Satellite Mission, CE-2, CE-3 and CE-5T1 Lunar Missions

      卫星名称 信号波段 时延/ns 时延率/
      (ps•s-1
      CE-2 实时S波段 4.93 0.64
      事后S波段 4.04 0.63
      事后带宽综合 1.28 0.45
      CE-3 实时X波段 0.78 0.93
      事后X波段 0.85 0.90
      CE-5T1 实时X波段 1.43 0.65
      实时S波段 3.44 1.13
      事后X波段 1.41 0.59
      事后S波段 2.24 1.07
      CE-4ZJ 实时S波段 0.91 0.67

      表 2为此次实时任务期间定位计算获得的赤经和赤纬误差。项目要求3测站测角精度为赤经0.08~0.37 as、赤纬0.12~0.40 as;3测站以上赤经0.05~0.09 as、赤纬0.07~0.12 as,从表 2可以看出,实时定位结果满足项目精度要求。

      表 2  定位计算的误差统计

      Table 2.  Error Statistics of Positioning Calculation

      观测代码 总点数 赤经误差/as 赤纬误差/as
      s8521a 2 827 0.007 0.009
      s8522a 3 170 0.004 0.007
      s8523a 2 662 0.006 0.007
      s8524a 1 417 0.012 0.014
      s8525a 3 495 0.014 0.013
      s8526a 3 112 0.004 0.007
      s8527a 2 652 0.003 0.003
      s8528a 2 640 0.013 0.013
      s8529a 2 683 0.005 0.007
      s8530a 2 359 0.015 0.024
      s8531a 2 337 0.007 0.017
      s8601a 2 412 0.019 0.017
      s8602a 1 842 0.011 0.017
      s8603a 1 994 0.006 0.010
      s8604a 1 958 0.002 0.005
      s8605a 1 811 0.002 0.005
      s8606a 1 804 0.006 0.007
      s8607a 2 222 0.007 0.008
      s8608a 1 811 0.071 0.012
      s8609a 1 767 0.004 0.008
      s8610a 1 811 0.004 0.007
      s8612a 925 0.082 0.041
      s8613a 1 535 0.052 0.050
      s8614a 2 171 0.007 0.014
      s8615a 878 0.006 0.010
      s8620a 913 0.008 0.014
    • CE-4任务实时阶段的主要轨道机动事件如表 3所示,定位归算用于监测轨道机动的实时变化。图 3图 4分别是观测代码为s8525a的中继星地心坐标和速度定轨结果与定位归算结果差值的变化。可以看出,在非轨道机动时间段,定位归算结果与定轨结果符合性较好,这说明定位归算结果的可靠性较高;当处于轨道机动时间段,定位归算结果与定轨结果有显著差异。这是因为在轨控阶段,卫星的受力情况复杂,难以进行准确模拟,变轨期间也不能以定轨长弧段方式给出实时轨道信息,因此轨控阶段定轨精度低于非轨控阶段定轨精度。定位归算是基于卫星与地面测站网的相对几何关系,与卫星具体的受力情况无关,所以定位归算精度不受是否轨控的影响[2]。在轨控阶段使用定位归算结果实时判断变轨情况是定位归算的一种特殊应用。

      表 3  CE-4任务实时阶段主要轨道机动事件

      Table 3.  Main Orbital Maneuver Events of CE-4 Mission Real-Time Phase

      观测代码 时间 事件
      s8521a 2018-05-21 14:54:20—14:54:32 中途修正
      s8525a 2018-05-25 13:32:10—13:46:57 近月制动
      s8527a 2018-05-27 13:46:10—16:46:21 中途修正
      s8607a 2018-06-08 00:00:00—00:00:15 捕获机动
      s8614a 2018-06-14 03:00:00—03:06:01 捕获机动

      图  3  观测代码s8525a定轨与定位归算地心坐标结果差值

      Figure 3.  Difference Between Orbit Determination and Positioning Reduction of the Geocentric Coordinates Result of Observation Code s8525a

      图  4  观测代码s8525a定轨与定位归算速度结果差值

      Figure 4.  Difference Between Orbit Determination and Positioning Reduction of the Velocity Result of Observation Code s8525a

      将定位归算结果的地心坐标和速度矢量转化为轨道根数,在非轨控阶段,轨道根数会保持稳定或缓慢变化;当进入变轨时间段,轨道根数会突然有明显的变化;变轨结束后,又渐渐趋于稳定,由此可以实时判断卫星变轨的时间及情况。图 5图 6分别是观测代码为s8521a的中继星轨道半长轴和轨道倾角在变轨前后的变化图像,图 7图 8分别是观测代码为s8525a、s8607a的中继星轨道偏心率在变轨前后的变化图像。变化图像表明,虽然由于数据的数量和精度的影响,结果存在一定波动,但可以清楚地看出,在非轨控阶段轨道根数较为平稳,当进入轨控阶段轨道会有明显变化,轨控结束后又趋于平稳。由此可以用定位归算结果转化的轨道根数实时判断CE-4中继星的变轨情况。

      图  5  观测代码s8521a轨道半长轴变化图像

      Figure 5.  Orbital Semi-major Axis Change Image of Observation Code s8521a

      图  6  观测代码s8521a轨道倾角变化图像

      Figure 6.  Orbital Inclination Change Image of Observation Code s8521a

      图  7  观测代码s8525a轨道偏心率变化图像

      Figure 7.  Orbital Eccentricity Change Image of Observation Code s8525a

      图  8  观测代码s8607a轨道偏心率变化图像

      Figure 8.  Orbital Eccentricity Change Image of Observation Code s8607a

      2018-05-21,中继星进行了第1次中途修正,表 4为中继星变轨时刻前后的轨道根数记录统计。从表 4可以看出,变轨时刻的轨道根数都有较明显的变化,说明确实按计划完成了中继星的轨道修正任务。2018-05-25中继星进行了近月制动,表 5为中继星近月制动时轨道偏心率的变化记录。从表 5可以看出,轨道偏心率从1.303逐渐减小到0.966并稳定,由偏心率大于1的双曲线轨道转变为偏心率小于1的椭圆轨道,说明中继星成功被月球引力场捕获并制动。

      表 4  第1次中途修正变轨前后轨道根数记录

      Table 4.  Oribital Element Records of First Track Correction Before and After Orbital Transfer

      时间 半长轴/km 偏心率 轨道倾角/(°)
      14:54:18201 4660.96728.496
      14:54:20201 4680.96728.541
      14:54:22201 4780.96728.505
      14:54:24201 3730.96828.466
      14:54:26201 3430.96728.607
      14:54:28201 2540.96828.676
      14:54:30201 2410.96728.799
      14:54:32201 0830.96828.793
      14:54:34201 0290.96828.824

      表 5  观测代码s8525a变轨前后偏心率记录

      Table 5.  Eccentricity Record of Observation Code s8525a Before and After Orbital Transfer

      时间 偏心率
      13:32:051.303
      13:33:051.288
      13:34:051.310
      13:35:051.224
      13:36:051.201
      13:37:051.176
      13:38:051.151
      13:39:051.130
      13:40:051.108
      13:41:051.087
      13:42:051.066
      13:43:051.046
      13:44:051.025
      13:45:051.037
      13:46:050.996
      13:47:050.971
      13:48:050.966
      13:49:050.966

      国内学者对探月工程的前几颗探月卫星的重要弧段进行变轨监测,通过观测卫星视向速度的变化计算卫星总速度变化量,不同的总速度变化量对应不同的变轨后轨道,最后与初始设计的变轨后轨道比较,进而确定卫星是否变轨成功。这种方法需要变轨前预先模拟不同卫星总速度变化量对应的变轨轨道,变轨时基于多普勒测速数据计算卫星总速度变化,并对结果进行插值[3],增加了额外的计算任务量,且当没有多普勒测距数据时也无法使用,无法判断实际变轨后的轨道与计划轨控后轨道的实际差异。使用定位归算结果转换的轨道根数判断变轨情况更加便捷,利用轨道根数的变化反映变轨情况也更加直观,不用增加额外计算,在应用实时VLBI(electronic VLBI,e-VLBI)后,VLBI的数据处理和传输速度更快,计算结果相对观测量滞后一般在20 s左右,能达到实时的要求。

    • 忽略太阳的引力,CE-4中继星处于地-月引力场中,对这类问题的轨道计算最常使用圆形限制性三体问题模型进行计算。圆形限制性三体问题模型用于质量可以忽略的小质点卫星在两个以圆轨道互相旋转环绕的两个大质量天体中受到引力作用下的运动问题[4]。在旋转的二体引力场中共存在5个平动点,也叫拉格朗日点。平动点是两个大质量天体引力场中的平衡点。在地月引力场中,卫星在平动点受到的引力平衡可以在基本不消耗推进剂的前提下保持与双星系统相同的公转角速度,且平动点可与地月两大天体组成稳定不变的空间构型,由于平动点特殊的动力学特性和相对地月系的固定位置,使得平动点在停泊中转、中继通信等深空探测应用中有很大的研究价值[5]

      平动点L2点位于地月连线的延长线上,距离月球大约6.5万km,在L2点附近轨道运行的中继星可以保持与地月不变的相对位置,并且可以持续与地球和月球背面不间断进行通信,还可以降低轨道转移消耗燃料,尽量延长中继星在轨工作寿命,因此是作为地月通信中继或月球背面探测的最佳位置[6]。Halo轨道是围绕L2点的周期性轨道,CE-4中继星最终将进入Halo轨道运行,以保证将要发射的于月球背面登陆的探测器与地球之间的通信联系。Halo轨道无法得到准确的表达式,只能通过摄动法获得三阶近似解析解[7]。本文将实时定位归算的结果转化为轨道根数,通过轨道根数的变化判断是否成功变轨,用变轨后的定轨结果与定位归算结果比较判断是否正确进入Halo轨道。

      2018-06-14 03:00,中继星进行了第2次捕获机动,之后将变轨进入Halo轨道开始周期运行。图 9图 10图 11分别为观测代码s8614a变轨前后轨道半长轴、偏心率、轨道倾角的变化,表 6为观测代码s8614a变轨前后轨道根数的记录。可以明显看出,03:00左右,中继星轨道根数出现了明显变化,几分钟后开始保持稳定,说明中继星在预定时间顺利完成了变轨任务。

      图  9  观测代码s8614a轨道半长轴变化图像

      Figure 9.  Orbital Semi-major Axis Change Image of Observation Code s8614a

      图  10  观测代码s8614a轨道偏心率变化图像

      Figure 10.  Orbital Eccentricity Change Image of Observation Code s8614a

      图  11  观测代码s8614a轨道倾角变化图像

      Figure 11.  Orbital Inclination Change Image of Observation Code s8614a

      表 6  观测代码s8614a变轨前后轨道根数记录

      Table 6.  Orbital Element Records of Observation Code s8614a Before and After Orbital Transfer

      时间 半长轴/ km 偏心率 轨道倾角/(°)
      03:00:1037 0850.91477.126
      03:00:4036 8860.92475.295
      03:01:0936 6380.93772.287
      03:01:3936 5330.94270.621
      03:02:0836 3920.95067.767
      03:02:3836 2330.95862.265
      03:03:0736 1460.96360.575
      03:03:3736 0480.96856.106
      03:04:0635 9490.97352.139
      03:04:3635 9000.97545.309
      03:05:0535 8510.97838.921
      03:05:3535 8130.98031.840
      03:06:0435 7980.98122.643
      03:06:3435 7840.98221.723

      图 12图 13分别是观测代码s8614a变轨前后定位归算与定轨结果的坐标、速度比较。可以看出,变轨前后定位归算的结果与定轨结果的符合度都很好,说明定位归算结果正确,使用定位归算结果能够判断中继星正确进入Halo轨道。

      图  12  观测代码s8614a实时定位归算与事后定轨结果的坐标差值

      Figure 12.  Coordinate Difference of Observation Code s8614a Between Real-Time Positioning Reduction and Post-Time Orbit Determination

      图  13  观测代码s8614a实时定位归算与事后定轨结果的速度差值

      Figure 13.  Velocity Difference of Observation Code s8614a Between Real-Time Positioning Reduction and Post-Time Orbit Determination

    • 本文统计的CE-4中继星任务实时阶段的观测资料和定位结果表明,在时延、时延率、定位精度上都满足任务要求。在中继星的地月转移段,月球-地月L2点转移段中,飞行器要经过多次中途轨道修正、近月制动、捕获机动等轨控弧段,通过定位归算得到的轨道根数变化可以实时高精度监测中继星的变轨情况,验证了定位归算在任务轨道机动时的特殊功用。

      通过定位归算,可以在第一时间准确判断出中继星的轨道修正和近月制动任务的正确实施,为后续的探月飞行任务提供了重要的帮助。通过分析CE-4中继星变轨进入Halo轨道前后的轨道根数变化以及将变轨后定位归算结果与定轨结果对比,可以实时准确地判断中继星变轨进入Halo轨道的时间及入轨情况,为后续飞行任务在Halo轨道处的监控提供了一种新的有效方法。

参考文献 (7)

目录

    /

    返回文章
    返回