3. 摄动源与轨道演化
二体问题给出的开普勒轨道只是理想模型：轨道六根数恒定不变，卫星沿固定椭圆无休止地运行。真实航天器却始终处在多种摄动力的共同作用下，轨道根数随时间缓慢漂移。轨道参数估计的核心任务之一，就是把这些摄动源正确地写进动力学模型，让滤波器能够区分"观测噪声引起的抖动"和"真实摄动引起的长期演化"。本节按摄动力量级和建模方式，逐一分析四类主要摄动源对轨道参数的影响，并说明它们在估计器中的实现形式。
一、摄动力量级与轨道高度的关系
摄动加速度相对于地球中心引力的比值，决定了它对轨道演化的影响速度。以地球引力加速度为参考，各摄动源的典型量级如下表所示（低轨 500 km、中轨 20000 km、地球同步轨道 42164 km 三种高度对比）：
| 摄动源 | 低轨 500 km | 中轨 20000 km | GEO 42164 km | 主要影响 |
|---|---|---|---|---|
| 地球非球形 J2 | ~1×10⁻³ | ~1×10⁻³ | ~1×10⁻³ | Ω、ω 长期漂移 |
| 大气阻力 | ~1×10⁻⁶～10⁻⁴ | 可忽略 | 可忽略 | a 衰减、e 减小 |
| 太阳光压 | ~1×10⁻⁷ | ~1×10⁻⁷ | ~1×10⁻⁷ | e 矢量缓慢旋转 |
| 第三体引力（日月） | ~1×10⁻⁷ | ~1×10⁻⁶ | ~1×10⁻⁵ | e、i 长周期变化 |
可以看到，J2 项在所有高度都占主导，而大气阻力只对低轨有意义，第三体引力则随轨道半径增大而相对增强。估计器建模时不必把所有摄动都写成同等精度，而应按任务轨道高度和精度需求取舍。
二、地球非球形引力：J2 项的长期漂移
地球引力位展开式中，J2 项（扁率项）是最主要的非球形摄动。它对轨道根数的影响具有明确的解析形式，其中升交点赤经 Ω 和近地点幅角 ω 的长期变化率为：
dΩ/dt = -1.5 * J2 * n * (Re/p)^2 * cos(i)
dω/dt = 0.75 * J2 * n * (Re/p)^2 * (5cos²i - 1)
其中 n 为平均角速度，Re 为地球赤道半径，p = a(1-e²) 为半通径，i 为轨道倾角。这两个公式说明：Ω 的漂移方向由 cos(i) 的符号决定，顺行轨道（i < 90°）升交点西退，逆行轨道东进；当 i = 63.43° 时 5cos²i - 1 = 0，ω 的长期漂移为零，这就是著名的"临界倾角"，常用于设计冻结轨道。
在估计器中，J2 通常以解析加速度的形式直接加入动力学方程。设卫星位置矢量为 r，J2 摄动加速度为：
a_J2 = -1.5 * J2 * μ * Re² / r⁵ * [ (1 - 5(z/r)²)·x,
(1 - 5(z/r)²)·y,
(3 - 5(z/r)²)·z ]
这一项在 ECI 坐标系中计算，无需额外状态量，直接叠加到二体加速度上即可。对于精度要求更高的任务，可继续加入 J3、J4 甚至 8×8 阶次的球谐项，但计算量随之上升。
三、大气阻力：低轨衰减与弹道系数
大气阻力加速度为 a_drag = -0.5 · ρ · (Cd·A/m) · |v_rel| · v_rel，其中 ρ 为大气密度，Cd 为阻力系数，A 为迎风面积，m 为质量，v_rel 为卫星相对大气的速度。阻力方向始终与相对速度相反，因此它持续消耗轨道能量，导致半长轴 a 单调衰减、偏心率 e 逐渐减小，最终轨道圆化并再入。
建模的难点不在公式本身，而在大气密度 ρ 的获取。ρ 随高度、太阳活动 F10.7 指数、地磁指数 Ap 以及地方时变化，常用模型包括指数模型、MSISE90 和 NRLMSISE00。在估计器中，ρ 的误差往往比动力学积分误差更大，因此工程上常把弹道系数 B = Cd·A/m 作为可估参数，与轨道状态一起在线估计，以吸收密度模型偏差。
四、太阳光压与第三体引力：中高轨的主要摄动
太阳光压加速度为 a_SRP = -P·Cr·A/m·(AU/|r_s|)²·û，其中 P ≈ 4.56×10⁻⁶ N/m² 为 1 AU 处的太阳辐射压，Cr 为反射系数（1 表示全吸收，2 表示全反射），û 为太阳指向卫星的单位矢量。光压对轨道半长轴和倾角影响很小，但会持续改变偏心率矢量，使 e 缓慢旋转，对 GEO 卫星的经度保持和冻结轨道设计有实际影响。
第三体引力来自太阳和月球。在 ECI 坐标系中，其加速度为 a_3rd = μ_k·[(r_k - r)/|r_k - r|³ - r_k/|r_k|³]，其中 r_k 为第三体位置。日月摄动对低轨影响微弱，但对 GEO 和高轨卫星，月球摄动可引起倾角以约 0.85°/年的速率增长，是 GEO 轨道倾角保持必须补偿的主要来源。
五、估计器中的建模策略与参数选择
把上述摄动写入估计器时，需要做三件事：确定摄动项集合、确定各项的计算精度、确定哪些量作为可估参数。工程上常用的策略是：
- 低轨（200～800 km）：J2 + J3 + 大气阻力，弹道系数 B 作为可估参数，密度用 NRLMSISE00 或指数模型。J2 用解析式，阻力用相对速度计算。
- 中轨（10000～25000 km）：J2 + J4 + 日月第三体 + 太阳光压，Cr·A/m 可作为可估参数。大气阻力可忽略。
- GEO 及高轨：J2 + 日月第三体 + 太阳光压，第三体位置由 JPL DE 星历或简化解析星历提供，光压面积和反射系数按卫星构型建模。
一个常见的误区是把所有摄动项都写成"可估参数"。实际上，J2 的系数是已知的地球物理常数，不需要估计；大气密度模型的偏差可以通过 B 来吸收；而第三体位置由星历给出，也不应作为状态量。真正需要在线估计的，通常只有弹道系数 B、反射系数 Cr 以及少数与卫星姿态和构型相关的面积参数。