基于误差发散规律的低轨卫星大气阻力系数计算方法-紫金山天文台.PDF

第57卷第4期 天 文 学 报 Vol.57 No.4 2016年7月 ACTA ASTRONOMICA SINICA Jul., 2016 doi: 10.15940/j.cnki.0001-5245.2016.04.006 基于误差发散规律的低轨卫星大气阻力 系数计算方法 1,2,3 1,2 4 1,2 汪宏波 赵长印 柳仲贵 张 伟 (1 中国科学院紫金山天文台南京210008) (2 中国科学院空间目标与碎片观测重点实验室南京210008) (3 宇航动力学国家重点实验室西安710043) (4 北京跟踪与通信技术研究所北京100094) 摘要 低轨卫星轨道预报精度受到大气模型和大气阻力系数精度的制约, 给一些高精度 的空间和航天任务带来不利影响. 提出了一种基于沿迹方向误差发散规律的大气阻力系 数计算新方法. 首先通过理论推导给出低轨卫星轨道预报中沿迹误差发散的分析表达式, 定量描述初值误差和模型误差对沿迹误差的综合影响; 提出利用定轨段的基本信息, 优选 预报段所采用的阻力系数, 抑制沿迹误差的发散速率, 从而降低沿迹方向预报误差的最大 值, 提高短期预报精度. 以400 km 附近的GRACE-A卫星的全弧段星载GPS高精度资料为 基础, 检验了方法的精度和成功率. 结果表明: 相对于传统的定轨预报方法, 新方法能提 高24 h短期预报精度约45%, 成功率约71%, 总体有效率约86%; 方法对低、中、高等3种 太阳辐射水平均有效, 对于中低等级的地磁扰动也有效, 具备较好的应用价值. 关键词 天体力学: 轨道计算和定轨, 方法：数值, 高层大气阻力 中图分类号: P135; 文献标识码: A 1 引言 一直以来, 低轨卫星大气阻力摄动加速度的计算存在较大误差, 是轨道预报误差的 最主要来源. 大气阻力难以精确计算的原因有两个方面: 一是大气密度模型普遍不准, 目前还没有哪个模型能够准确描述和预测高层大气在不同时空、不同环境下的大气密 度 ; 二是大气阻力系数存在较大偏差, 从物理的角度看, 它与大气分子的成分和浓度、 卫星表面材料及温度、撞击角度和速度等因素有关, 很难被精确测定. 计算大气阻力系数一般有数值模拟和轨道拟合两种方式, 其中数值模拟是基于大气 分子与卫星表面撞击作用的理论模型, 输入各种条件参数进行大量的蒙特卡罗模拟, 从 2015-11-14收到原稿, 2015-12-03收到修改稿 宇航动力学国家重点实验室开放基金资助 whb@pmo.ac.cn 448 天 文 学 报 57 卷 而得到阻力系数 . 这种方法比较接近物理本质, 结果也相对精确, 但缺点是在实际任务 中很多输入条件是不确定的, 而且计算耗时巨大, 适用性和时效性较差. 通过轨道拟合方法计算大气阻力系数, 是指在定轨阶段将大气阻力系数C 作为待 ˙ 估参数, 与卫星状态参数(⃗r, ⃗r) 共同求解, 进而应用在轨道预报中. 这种方法原理简单, 计算方便, 能够取得稳健的预报精度, 在实际任务中得到了普遍应用. 但是对于部分低 轨卫星的高精度空间任务, 例如: 空间交会、碰撞预警、精密跟踪等, 目前的预报精度 还不能满足要求. 这主要是因为定轨段轨道拟合中, 阻力系数与卫星轨道状态值、大气 模型形成了自洽系统, 阻力系数C 是定轨阶段的最优估计, 它吸收了观测资料、大