海基单站短弧光学跟踪GEO精密定轨

2022-06-04宋叶志谭龙玉张少伟曾春平

兵器装备工程学报 2022年5期

邵瑞，宋叶志，陈赟，谭龙玉，张少伟，叶钊，曾春平

(1.中国科学院上海天文台，上海 200030； 2.中国科学院大学，北京 100049；3.上海航天控制技术研究所，上海 201109； 4.航天东方红卫星有限公司，北京 100094)

1 引言

随着我国航天事业的飞速发展，越来越多的卫星进入太空，在导航、定位、通信、遥感等领域发挥着重要的作用。近年来各大国对空间资源的竞争愈发激烈，使得轨道资源愈发紧张，尤其是拥有重要战略地位的地球同步轨道资源。因此，为更加合理地规划轨道资源、避免卫星碰撞、保障太空安全，各国都对空间目标的精密定轨提出了更高的要求，除了要对已知的合作目标进行精密跟踪，还需要对非合作的空间目标进行监测。

传统的地基监测系统在时间、空间的监测上均存在监视盲区，无法满足我国当前对空间态势感知能力的高要求。自20世纪末，美国就致力于天基监测系统的建设，试图将地基系统与天基观测相结合，实现对天区的全面覆盖。相比美国，中国对天基系统的建设起步较晚，且在地基台站的建设上同时受到政治、地理等因素的限制，在空间覆盖上还需要其他的手段来进行补充。海基观测可以在陆地以外的海域进行，包括我国的领海以及公海，所以海基观测在一定程度上可以弥补地基观测的不足，我国天基观测系统建成后，海基观测依然可以作为一个重要组成部分，为航天测控提供有力的支持。

目前地基测控研究已经有了诸多成果，于涌等基于上海天文台佘山站1.56 m望远镜，利用单站CCD漂移扫描技术对北斗同步卫星进行定轨，单圈精度为数千米量级，多圈精度可达50 m。宋叶志等利用多台站，采用双频双程测距模式，对风云四号进行定轨，非变轨期间精度优于20 m。YongHuang等利用VLBI测轨数据以及C波段转发式测距数据对GEO进行定轨分析，VLBI时延测量精度为3.6 ns，定轨精度约为20 m。美国从1996年发射第一个携带天基光学传感器的低轨卫星，论证了单个平台对高轨带监视的可行性开始，在此后的20多年的时间里，不断发射卫星，构建SSBS天基太空监视系统。国内的天基监视系统尚未建设完成，周海银等设计SBSS系统覆盖区域的简易算法，对天基监视系统进行仿真，得到低轨目标精度10 m，高轨目标精度500 m左右的结论。周庆勇等基于天基光学监视平台的空间目标可视性以及轨道改进算法，对多平台联合定轨进行仿真。海基测控在国内起步较晚，康德勇等在2005年提出，海上移动站精密测轨的方法，并对低轨卫星进行了定轨。在2010年，分析了船位误差对船载外测数据的影响，随着距离的增加对定轨结果的影响逐渐减弱。朱伟康等建立测速数据的船姿船速修正模型等一系列优化途径，使得初轨计算半长轴外符合精度提高3倍。陈红英等使用最小二乘法解算双站目标坐标的数据处理方法，提高了弹道状态估计精度。倪兴等提出在航天测控数据处理中使用简化的大气折射公式，得到的轨道精度满足要求，同时满足实时数据处理数据的计算速度。

目前航天测量船的测控手段多为雷达测控，且观测对象多为低轨目标，应用于高轨光学监测的相对较少。本研究使用海基单站测量设备采集到的光学测角数据，利用数值方法，在观测弧段较短且不连续的情况下，对GEO卫星进行轨道改进。

2 海基天文光学测量方法

海基天文光学测量方法是利用CCD成像中的空间目标和背景恒星的相对位置，通过计算其理想坐标，得到其赤经赤纬。

以原点建立直角坐标系-,是CCD底片上的一点，为赤纬圈投影，与垂直，为赤经圈投影，正方向取赤经赤纬增加的方向，该坐标系为理想坐标系。望远镜视场平面与天球的交点为(,)。空间目标在该坐标系中的投影坐标(,)与其赤道坐标(,)的转换公式为：

(1)

理想坐标系和量度坐标系由于设备参数的原因在中心原点、坐标轴比例以及坐标轴平行等方面存在差异，可以通过多项式建立二者之间的关联。

底片模型参数是通过理想坐标和量度坐标之间建立方程，利用最小二乘原理求得。利用底片模型参数，可以得到观测目标在理想坐标系下的坐标，进而得到其观测坐标，即赤经、赤纬

(2)

3 海基测定轨原理

3.1 海基测定轨涉及的坐标系统

用于定轨的观测资料是在惯导地平坐标系-(舰船地理坐标系)下描述的，其坐标原点是舰船摇摆运动中心在舰船甲板上的投影，通过惯导三轴(北向、东向、方位)定义基准面。指向正东，指向正北，指向由,,按右手螺旋规则确定。此外还有惯导甲板坐标系-(舰船甲板坐标系)，其中心与与重合，指向与艏艉线一致，指向垂直，平行于甲板面，以右舷为正，指向由和按照右手螺旋规则确定；原始的观测资料是在测量坐标系-(设备直角坐标系)下描述，该坐标系的中心位于设备中心，三轴指向与-一致。

将原始观测资料从-下转换到-需要借助船载的形变测量设备和惯导测量设备。假设某个时刻，形变测量设备的观测值为(,,)，该值是设备中心到()之间几何矢量的形变；惯导设备的观测量为(,,)，该值分别是舰船的艏摇角、纵摇角和横摇角。

图1 船载坐标系示意图

由-转换到-，考虑船体形变，进行如下变换：

(3)

再由-转换到-，考虑船体姿态，进行如下变换：

(4)

其中旋转矩阵具体表达如下：

3.2 光学测角精密定轨方法

空间目标在围绕地球进行运行的过程中会受到许多力的作用。这些力大致可以分为两类，一类是保守力，诸如中心天体的引力，其他大型天体的引力等。另一类是耗散力，则系统能量损耗，对于卫星而言，主要是大气阻力、太阳光压、地球红外辐射等。

在惯性系中，空间目标的运动方程可以表示为：

(5)

观测量为赤经赤纬，其观测方程为：

(6)

Δ=-

(7)

式中：(,,),(,,)分别为目标卫星和观测平台在天球参考系下的坐标；,为大气折射改正；,为周年光行差改正；,为周日光行差改正。

设在时刻有一组观测量：

(8)

(9)

对式(9)进行变形可得：

(10)

(11)

(12)

式中，为待估参数的个数，=+6。

若令

(13)

可得到观测方程，即条件方程：

(14)

利用协方差分析理论可以研究不同误差源对估值精度的影响，将观测方程表达如下：

=++

(15)

式中：为求解参数；为考察参数；和分别待估参数和考察参数的偏导矩阵。定义灵敏矩阵为：

(16)

(17)

式中：为计算协方差；为权阵。另外定义摄动矩阵=π12，其意义是考察参数有1的先验误差时，该误差对待估参数的影响。摄动矩阵是评价解算的外符合精度的重要统计量。其中：

(18)

利用灵敏度矩阵和摄动矩阵，考察协方差阵可以表述为：

=+π=+

(19)

因为为正定矩阵，式(19)意味着考察参数的误差会降低待估参数的精度，并且协方差分析可定量的计算这种影响的大小。

4 数据处理与分析

4.1 仿真分析

对海基单站光学跟踪GEO卫星定轨进行仿真模拟。初始历元为2021年6月5日14时，台站(船载设备)的初始位置以及目标GEO卫星的初始星历如表1所示(船载设备的位置以经纬度、高程表示，GEO卫星的初轨以无奇点轨道根数表示)。表2所示为台站的初始位置。

表1 GEO目标的初始星历

无奇点轨道根数与开普勒轨道根数的关系式如下：

(20)

式中：为轨道半长轴；为轨道倾角；为轨道偏心率；为近地点幅角；为升交点赤经；为平近地点角。

表2 台站的初始位置

仿真的观测噪声为3 ″，该噪声为随机误差，服从高斯分布，采样率为3 s，模拟的观测弧段为2021年6月5日14∶00—15∶00以及16∶00—17∶30，观测不连续。

4.1.1 台站位置噪声为50 m

将台站位置在地固系下3个方向上的分量依次加上50 m噪声，使用时间序列化的台站位置和模拟的观测数据进行定轨，定轨得到的残差如图2所示。

图2 位置噪声为50 m的定轨残差示意图

改进后的轨道和仿真参考轨道在观测弧段内的重叠精度如图3所示。

图3 改进轨道与仿真轨道曲线

4.1.2 台站位置噪声为10 m

将台站位置在地固系下3个方向上的分量依次加上10 m噪声，同样使用时间序列化的台站位置和模拟的观测数据进行定轨，定轨得到的残差如图4所示。

图4 位置噪声为10 m的定轨残差示意图

改进后的轨道和参考轨道在观测弧段内的重叠精度如图5所示。

对比台站位置噪声振幅50 m和10 m的定轨结果可以看出：不同的噪声振幅对于定轨残差的影响较小，这是由较小的台站位置偏移和台站到目标的长距离所共同决定；不同的噪声振幅对于轨道重叠的精度影响有所差异，但是基本处于相同的量级，相比较于其他观测手段在相同条件下，精度也基本符合。

图5 改进轨道与仿真轨道曲线

4.2 实测数据处理

本文利用某次海上测控任务采集的数据进行分析，观测目标主要针对静止轨道卫星。选取观测弧段为2021年6月5日14∶00—15∶00以及16∶00—17∶30内的四颗GEO卫星进行数据解算，根据GEO卫星的NORAD 编号查询其对应的两行根数，转换成位置速度作为轨道改进的初始值，4颗GEO卫星在历元为2021年6月5日14时的初轨如表3所示(初轨使用无奇点根数表示)。