磁暴下的LEO卫星精密轨道确定.pdfVIP

磁暴下的LEO卫星精密轨道确定 陈建荣王家松 (西安卫星测控中心·西安·710043) 摘要对低轨卫星(LEO)，大气阻尼摄动是主要的定轨误差源。尤其在发生 磁暴时，求解一个大气阻尼因子的定轨方法已不能充分吸收大气密度计算不准所 造成的定轨误差，因而在标校统一S波段(USB)的测量系统差和随机差时往往 计算失真．本文提出了一种求解折线型Cd因子的新方法，克服了动力学模型不 准所带来的定轨．误差，通过和独立的GPS数据比较，定孰精度有明显提高，同时 给出的测量系统差和随机差更加真实可信． 关键词轨道确定；磁暴；Cd因子；大气密度 l 引 言 随着卫星应用技术的发展，对卫星轨道确定精度的要求也越来越高。例如我国第二 代“海洋卫星”要求径向定轨精度达到分米级。定轨精度与测量数据的数量和质量、测量 数据的时间和地理分布，以及卫星动力学模型精度密切相关。当测量精度较高、覆盖较好， 而卫星动力学模型精度较低时，我们可以在定轨中增加估计大气阻尼因子、太阳光压系数， 或引入周期性的经验力，甚至采用简化动力学方法来补偿力模型的不足(文献【l，2】)。同 样如果测量精度较低，而卫星动力学模型精度较高，我们可以在定轨中增加估计测量数据 的系统差。实际定轨中必须综合考虑测量精度和动力学模型精度，通过不断调试，制定最 佳的定轨策略。 通常的卫星标校方法，认为卫星在已知摄动力模型作用下，按照理论力学轨道飞行， 轨道确定过程中计算观测量和实测值之差(即残差)假定主要是由于测量误差所导致。在 定轨过程中通常求解卫星的初始位置和速度、一个大气阻尼因子Cd(补偿大气阻力摄动计 算误差)和每圈每站每种测量类型一个常数系统差。这种方法对中高轨卫星效果很好，但 对低轨卫星，由于动力学模型误差导致的轨道计算误差已经上升为主要矛盾(假定USB的 动比较剧烈(即磁暴)的时候，这种方法应用效果不佳。 由于对上大气层的变化和组成缺乏足够的了解，所以目前所有的大气密度模型都不是 很精确。对LEO卫星，如果定轨弧段选取太长，则在积分运动方程时，大气阻尼误差的累 积效应将很明显；如果定轨弧段太短，则参加定轨的测量数据太少，也就忽略了大量数据 的统计特性，所以我们通常选取的定轨弧长一般为O．8天～3天。以200公里高的YG卫星 454 A采用1．5天数据定轨，算例B采片j0．8天数据定轨。我们发现两组算例在发生磁暴的时 候，标校的M1站和w站测量系统差存在着明显差异，而这种差异在未发生磁暴的时候并 不突出。 表1 YG卫星任务中，两组算例给出的部分站测量系统差和随机差(9051]91圈) 圈次 站名 算例 测距／m 测速／(,：nVO 随机差 系统差 随机差 系统差 91 Ml A 37．3 122．6 4．6 77．3 B 1．8 ．56．9 6．1 ．12．8 90 w A 15．2 ．7．0 4．0 40．6 B 1．3 17．4 6．1 ．24．1 表l表明磁暴的发生给轨道确定带来了困难。如何在发生磁暴的情况下，仍然保证低 轨卫星(LEO)定轨精度的连续性和测量标校值的可信度，是我们需要解决的难题。本文 提出了一种求解折线型Cd因子的新方法，解决了这个难题。这种新方法不仅在发生磁暴 情况下可以大幅度提高LEO卫星的定轨精度，而且在非磁暴情况下，也有很好的应用效果。 2求解折线型Cd因子的方法 通常的定轨方法中，全弧段仅解一个Cd，这只能吸收大气密度计算误差的平均部分， 而大气密度快速波动部分带来的误差则很难吸收。新方法中采用解多个大气阻尼因子来克 服大气阻力计算误差。比如对2天定轨弧段，可以求解5个大气阻尼因子，定轨弧段起始 插值获得。图l示意了Cd随时间的变化： 图1 计算多个cd因子的示意图 其他任意时刻对应的Cd，用公式表达如下 啪)2糟啪小告‰1