[发明专利]一种月球探测器精密定轨方法

说明书

技术领域

本发明属于月球探测技术领域，涉及一种月球探测器精密定轨方法，具体地说，涉及一种利用甚长基线干涉测量(VLBI)观测、同时估计地球自转参数和力学模型参数的月球探测器精密定轨方法。

背景技术

月球探测器的精确位置信息对实现探月计划、科学目标和研究地月环境非常重要。VLBI技术被广泛应用于深空探测领域。国际上许多深空探测器的实时测轨过程中都成功使用了VLBI技术，我国探月工程中对嫦娥系列探测器的跟踪与测控也应用了VLBI技术，相关文献中均有描述。

VLBI技术的原理如图1所示。由于探测器距地面两测站的距离不同，因此探测器上发出的同一信号到达两测站的时间是不同的，测量这个时间差，建立探测器与两测站的几何测量关系，通过大量的观测值以及已知的测站位置就可以计算得到探测器的位置。

由于地面测站位于地球上，其坐标一般在协议地固参考系下描述，而探测器由于位于深空，其坐标根据其位置的不同又定义在不同的参考系统下，如环月探测器的坐标一般在J2000月心惯性系下描述。因此要通过VLBI技术建立地面测站与探测器的几何观测连接，就必须将它们定义在统一的参考系中。在深空领域中，常将两者统一在太阳质心参考框架(BCRS)下。这样就涉及到一个坐标系转换的问题。图2描述了由地固坐标系转换到太阳质心坐标系的过程。

在图2中，由地固坐标系转换到J2000地心惯性系可用下式表示：

其中，[GCRS]表示在J2000惯性坐标系下的坐标，[ITRS]为协议地固系下坐标，Q(t)为岁差章动矩阵，R(t)为地球旋转矩阵，W(t)为极移矩阵。其中，R(t)与W(t)由地球自转参数(ERP)决定，因此，ERP参数的精度高低将直接影响着转换过程的精度高低，从而间接影响了观测方程的建模精度，并最终影响到探测器的定轨精度。相关文献研究表明，不同ERP先验精度下“嫦娥一号”等深空探测器的定轨结果最大将会存在公里级的差异(魏二虎，史青，严韦，李雪川.ERP精度对“嫦娥一号”差分VLBI定位精度的影响.武汉大学学报·信息科学版，2011，36(11)：1324-1327)。

ERP参数目前由国际地球自转服务(IERS)发布，分为预报产品(Bulletin A)和事后精密产品(Bulletin B)两种，事后精密产品精度较高，可满足目前深空探测器的定轨要求，但由于是事后延迟发布，因此时效性低，在探测器的实时性定轨中不能使用。因此目前的深空探测器实时定轨任务使用的是精度较低的预报产品，影响了深空探测器实时定轨精度的提高。

目前，有学者提出了将ERP参数引入到观测方程中与探测器位置参数一同估计，在提高ERP参数精度的同时提高探测器位置的解算精度，有效地降低了ERP参数的预报精度对定位精度的影响。但这种方法假设了在一个较短的弧段内探测器轨道根数没有变化，限制了轨道确定的长度；而且基于几何定轨的原理，没有动力学信息，因此得到的轨道没有预报作用。

对月球探测器的每一种观测技术的观测模型中，涉及很多几何和物理模型参数的先验信息，如果直接利用而不改进其先验信息的精度，就会传播到所确定的月球探测器位置中。月球探测器在轨飞行时，其位置的先验信息一般分别通过经验的动力学模型积分获得。由于动力学经验模型参数不准确，影响了给出的数值积分位置的精度，从而影响了采用观测量确定探测器位置的精度。目前观测模型中不同参考系间的连接参数直接采用本观测系统以外观测系统所提供的先验信息，其精度水平也会影响月球探测器定位精度。如何改进这些几何和物理模型参数的先验精度，是目前月球探测器精密定位研究急需突破的关键问题和难点。对于月球探测器定轨，目前国际上使用的基本模型有几何观测模型和动力学定轨模型。这两类模型各有优劣。

发明内容

为了克服现有技术中存在的缺陷，本发明提供一种月球探测器精密定轨方法，在建立VLBI观测数据的观测模型基础上，提出了将观测模型与动力学模型进行有机结合，首次提出和建立一种同时估计参考系连接参数(ERP参数)、探测器初始状态参数、探测器飞行段动力学定轨的摄动模型参数的新模型，避免现有模型中直接引用这些几何和物理模型参数的先验精度对探测器定位精度影响的缺陷，为月球探测器精密定位建立基础。

其技术方案如下：

一种月球探测器精密定轨方法，包括以下步骤：

步骤1：将观测模型和动力学模型进行有机结合，建立新的结合模型分析月球探测器飞行段进行动力学定轨所涉及的摄动力因素，分析各摄动力对月球探测定轨的影响量级，研究摄动力模型。