3. 轨道摄动源识别

二体问题给出的开普勒轨道只是理想模型：航天器只受一个质量均匀球体的中心引力。真实航天器同时受到多种附加作用力，这些力相对中心引力的偏差统称为轨道摄动。轨道预报的第一步不是急着积分，而是识别“有哪些摄动源、各自来自什么物理机制、在给定轨道上量级多大”，因为不同摄动项对轨道根数的影响方向、时间尺度和空间分布完全不同。本节按物理来源逐项分析，并给出量级比较方法。

一、摄动加速度的统一表达与量级基准

设航天器位置矢量为 r，中心引力加速度为 a₀ = −μ_Er/r³，其中 μ_E ≈ 3.986×10¹⁴ m³/s²。所有摄动项写成附加加速度 a_p，则运动方程为 r̈ = a₀ + Σa_p。工程上常用相对量级 a_p/a₀ 来比较，因为 a₀ 随高度迅速衰减，同一摄动源在不同轨道上的相对重要性差别很大。例如 500 km 高度 a₀ ≈ 8.4 m/s²，而 36000 km 高度 a₀ ≈ 0.22 m/s²，相差近 40 倍。

二、地球非球形引力：最主要的保守摄动

地球并非理想球体，赤道隆起使引力场呈扁球状。用球谐展开时，J₂ 项（约 1.0826×10⁻³）远大于其他带谐项。J₂ 摄动加速度量级约为 a₀·(3/2)J₂(R_E/r)²，在 500 km 高度约 1.1×10⁻² m/s²，相对量级约 1.3×10⁻³。它的物理效应是使升交点赤经 Ω 和近地点幅角 ω 产生长期漂移，对近地轨道尤为显著。J₂ 是保守力，不改变轨道能量，但会持续改变轨道面朝向和拱线方向。更高阶项 J₃、J₄ 及田谐项量级依次递减，通常在精密预报中保留到 8×8 或更高阶。

三、大气阻力：低轨道的耗散性摄动

大气阻力加速度为 a_D = −(1/2)ρv_rel²(C_DA/m)v̂_rel，其中 ρ 为当地大气密度，C_D 为阻力系数（典型 2.2），A/m 为面质比。它沿相对速度反方向，属于耗散力，持续降低轨道半长轴和偏心率。在 400 km 高度，ρ 约 10⁻¹¹～10⁻¹² kg/m³，若 A/m = 0.01 m²/kg，a_D 约 10⁻⁶～10⁻⁷ m/s²，相对量级约 10⁻⁷～10⁻⁸。虽然量级远小于 J₂，但它是长期累积效应，数天到数周即可造成公里级位置偏差。大气密度受太阳活动、地磁指数和地方时影响，可变化一个数量级以上，因此阻力是低轨预报中不确定性最大的摄动源之一。

四、第三体引力：远距离的潮汐式作用

第三体引力指月球和太阳对航天器的直接引力与对地球的间接引力之差，即潮汐加速度。对地月系统，太阳质量虽大但距离远，月球虽小但距离近，两者量级接近。第三体摄动加速度近似为 2μ₃r/d³，其中 d 为第三体到地心距离。在 GEO 轨道上，太阳和月球引起的加速度约 10⁻⁶～10⁻⁷ m/s²，相对量级约 10⁻⁶。对高轨和深空轨道，第三体引力是主导摄动，会引起轨道偏心率、倾角和升交点赤经的周期性或长期变化。对近地轨道，其量级通常小于 J₂ 和大气阻力。

五、太阳光压与潮汐：小量级但不可忽略

太阳光压加速度为 a_SRP = −(P_SR·A/m·C_R)·ŝ，其中 P_SR ≈ 4.56×10⁻⁶ N/m² 为 1 AU 处太阳光压，C_R 为反射系数（1～2）。对 A/m = 0.01 m²/kg 的航天器，a_SRP 约 5×10⁻⁸ m/s²，相对量级约 10⁻⁸。它主要影响高面质比航天器（如太阳帆、大型太阳翼卫星）和 GEO 以上轨道。地球固体潮、海潮和大气潮则通过改变地球引力场分布产生附加摄动，量级通常比 J₂ 小两到三个数量级，在厘米级精密定轨中才需精细建模。

六、量级比较与工程判断

把上述摄动项放在同一张“量级谱”上，可以得到清晰的工程判断：对 500 km 近地轨道，J₂ 相对量级约 10⁻³，大气阻力约 10⁻⁷～10⁻⁸，第三体约 10⁻⁸，光压约 10⁻⁸，潮汐约 10⁻⁹。因此近地轨道预报中，J₂ 决定轨道面的长期演化，大气阻力决定半长轴的衰减速率，其余项在短期预报中可作为小量修正。对 GEO 轨道，J₂ 相对量级降至约 10⁻⁵，第三体引力和太阳光压升至 10⁻⁶～10⁻⁷，成为主导摄动，大气阻力几乎可忽略。对深空轨道，第三体引力和光压进一步主导，需采用行星星历和光压模型联合积分。

七、摄动源识别检查表

实际做轨道预报时，可按以下检查表逐项确认是否纳入模型，避免漏项或过度建模：

识别摄动源之后，下一步才是选择积分器、步长和力模型阶数。量级比较的意义在于：它决定了哪些项必须精细建模、哪些项可以简化，从而在预报精度和计算效率之间取得平衡。对同一航天器，不同任务阶段（如近地点、远地点、地影区）摄动主导项也可能切换，需要在预报过程中动态判断。