基于反向运动学精密单点定位的GPS Block IIF卫星食季改进偏航姿态模型

《Satellite Navigation》：An improved yaw attitude model for GPS Block IIF satellites during eclipse seasons

　　

专业术语翻译（摘要部分） 本研究利用反向运动学精密单点定位（Reverse Kinematic Precise Point Positioning, RKPPP）偏航估计作为观测驱动参考，改进了全球定位系统（Global Positioning System, GPS）Block IIF卫星在食季的运控偏航姿态模型。研究人员利用一套长达7.5年的全球GPS跟踪数据集，重建了食季机动期间的偏航姿态时间序列，并将其与当前运控姿态模型进行了比较。比较结果揭示了四类反复出现的差异。其中三类发生在正午转向（noon turn）期间：（1）偶发的旋转方向不匹配；（2）食前偏航演化不一致，导致机动起始时刻出现旋转方向的突变；（3）明显的机动提前启动，表现为较大的偏航差异。第四类差异发生在午夜转向（midnight turn）期间，当β角接近零时，模型预测的旋转方向与RKPPP推断结果不一致。针对上述问题，研究人员为正午转向引入了卫星特定的β角偏差规则，并在食前采用经验修正的β角来调整偏航速率并延迟机动启动。对于小绝对值β角下的午夜转向，研究人员提出了一种基于机动起始时刻β角符号及β角变化率符号的判定规则。将该修正模型应用于网络轨道与钟差估计后，与姿态机动相关的相位残差特征显著减少，偏航转向期间的卫星钟差估计更加平滑且一致性更好。

论文解读文章

一、研究背景与问题

全球导航卫星系统（Global Navigation Satellite System, GNSS）卫星姿态的精确建模是精密定轨（Precise Orbit Determination, POD）、卫星钟差估计、观测值特定信号偏差（Observable-specific Signal Biases, OSBs）确定以及高精度精密单点定位（Precise Point Positioning, PPP）的关键前提。姿态误差会改变发射机有效相位中心几何关系，从而影响水平相位中心偏移的改正，同时还会影响载波相位缠绕效应以及太阳辐射压建模中的光照几何关系。在名义偏航控制（yaw-steering）模式下，GPS卫星保持太阳帆板指向太阳、导航天线指向地球。然而，当太阳相对于轨道面的高度角（即β角）趋近于零时，名义偏航控制律要求极高的偏航速率；在轨道正午和午夜附近，所需偏航速率可能超出星载硬件能力。此外，在地球食期间，太阳敏感器可能无法提供可靠姿态参考，实际姿态偏离名义偏航控制律，此时需要启动专门的正午转向、午夜转向和阴影控制模型。若在这些非名义工况下仍采用名义姿态模型，将产生系统误差，进而影响卫星钟差和相位OSB估计的精度，最终降低PPP解的精度。

GPS姿态建模经历了从GYM—94/95框架到Block IIR卫星动态偏航机动策略，再到Block IIF卫星固定偏航速率的演化过程。对于Block IIF卫星，现有运控食季姿态模型在正午转向期间采用固定偏航速率0.11°/s，并引入约0.7°的β角偏差来确定旋转方向；午夜转向则采用影锥进入与退出时刻之间的线性插值。然而，将该模型应用于精密定轨和钟差估计时，在部分正午和午夜转向机动期间仍观测到残余不一致性，表明现有模型仍有改进空间。因此，有必要利用观测驱动的方法估计Block IIF卫星的真实偏航姿态，并据此改进姿态机动模型。

二、研究内容与主要结论

本研究利用反向运动学精密单点定位（Reverse Kinematic Precise Point Positioning, RKPPP）技术，对12颗GPS Block IIF卫星在2017年11月1日至2025年4月20日（约7.5年）期间的食季偏航姿态进行了估计。研究识别出约3,935次正午转向事件和4,115次午夜转向事件，并将RKPPP估计的偏航姿态时间序列与现有运控模型进行了系统比较，发现四类反复出现的差异。

第一类差异发生在正午转向期间，共识别出88例模型旋转方向与RKPPP估计不一致的事件。研究发现，采用统一的?0.7°β角偏差并不能适用于所有Block IIF卫星，因此提出了卫星特定的β角偏差规则，各卫星的偏差值由长期RKPPP偏航姿态序列经验确定。应用该规则后，3,935次正午转向中仅剩6例旋转方向仍与RKPPP估计不一致，且这些例外均发生在β角非常接近偏差阈值的边界情形，绝对偏差仅约0.01°至0.05°。

第二类差异涉及食前偏航演化不一致。对于需要反向旋转的事件，旧模型在正午转向开始前保持名义偏航控制，仅在食影区间内施加反向旋转，导致机动起始时刻出现旋转方向的突变。研究人员将基于β角偏差的反转判定准则扩展至整个正午转向弧段，在正午转向进入前即触发反转，使模型偏航角呈现单一连续旋转，与RKPPP估计更加吻合。

第三类差异为机动提前启动。统计显示，3,935次正午转向事件中，177次事件的模型与RKPPP估计的最大偏航差超过40°，930次介于20°至40°之间。大偏差在负β角情形下更为频繁，且各卫星的最大偏航差均出现在接近其卫星特定β角偏差处。研究人员在食前偏航控制律中引入经验修正的β角（β_new），通过分段映射压缩β角绝对值，从而降低食前偏航速率并延迟机动启动。应用修正后，所有事件的偏航差均降至40°以下，但仍存在95次事件偏航差大于20°。

第四类差异发生在午夜转向期间，当|β|接近零时，33次事件的模型旋转方向与RKPPP估计不一致。研究人员提出基于机动起始时刻β角符号和β角变化率符号的判定规则来修正旋转方向。

将修正后的模型应用于网络轨道与钟差估计后，正午转向期间的POD后验相位残差显著降低，例如G066和G070卫星的残差95分位数从13.8 cm和9.6 cm分别降至3.6 cm和3.4 cm；卫星钟差序列在机动期间的跳变显著减小，变得更加平滑和一致。

三、主要技术方法

本研究采用的核心技术方法包括：其一，反向运动学精密单点定位（RKPPP）技术，该方法利用逐历元观测几何估计卫星天线相位中心（Antenna Phase Center, APC）的时变位置，结合已知的机体坐标系相位中心偏移（Phase Center Offset, PCO）和由POD获得的质心（Center of Mass, CoM）轨迹，逐历元反演偏航角。其二，迭代处理策略，即先以现有食季姿态模型初始化POD，再将RKPPP估计的偏航角反馈作为下一次迭代POD的规定姿态，重复直至偏航估计收敛；收敛性分析表明第四次迭代后两条独立初始化路径的RMS差异不超过3°。其三，卫星特定β角偏差的确定方法，基于RKPPP参考对各卫星全部正午转向事件的旋转方向分类，经验确定各卫星的反转触发阈值。其四，分段β角映射修正，通过经验修正的β角调整食前偏航速率。数据处理采用武汉大学研发的PANDA软件包，使用169个全球分布的IGS跟踪站观测数据。

四、研究结果

正午转向模型改进。 针对正午转向期间的三类差异分别提出修正方案。对于旋转方向不匹配问题，通过卫星特定β角偏差规则，将不匹配事件从88例降至6例，且残余例外均属边界情形。对于食前偏航演化不一致问题，将反转判定准则扩展至整个正午转向弧段，消除了机动起始时刻的旋转方向突变。对于机动提前启动问题，通过分段β角映射修正食前偏航速率，使所有事件的偏航差降至40°以下。在代表性算例中，修正后模型的偏航角与RKPPP估计的差异保持在约10°以内，POD相位残差从13.8 cm和9.6 cm分别降至3.6 cm和3.4 cm，钟差序列的跳变从7.4 cm和9.2 cm分别降至4.9 cm和2.0 cm。

午夜转向模型改进。 针对|β|接近零时33次事件的旋转方向不匹配问题，提出基于机动起始时刻β角符号和β角变化率符号的判定规则，修正了模型的旋转方向。

五、结论与讨论

本研究通过RKPPP技术对GPS Block IIF卫星食季偏航姿态进行了长达7.5年的观测驱动估计，系统揭示了现有运控姿态模型在正午转向和午夜转向期间的四类系统性缺陷，并提出了相应的修正方案。卫星特定β角偏差规则取代了统一偏差，分段β角映射有效延迟了机动启动，基于β角符号和β角变化率符号的判定规则修正了午夜转向的旋转方向。这些改进在精密定轨和钟差估计中显著降低了姿态相关误差，减少了相位残差特征，使卫星钟差估计在偏航转向期间更加平滑和一致。该研究对精密定位和授时应用以及未来的GNSS数据重处理具有积极意义。论文发表在《Satellite Navigation》。需要指出的是，所提出的β角偏差值和分段映射参数均属于经验性修正，并非卫星制造商姿态控制系统的物理常数，可能隐含吸收了Block IIF姿态控制实现中未公开的特征。