4. 轨道漂移的物理机制
轨道保持的核心任务,是让航天器在任务寿命内始终满足轨道约束。要做到这一点,必须先回答一个基础问题:轨道根数为什么会随时间变化?答案来自摄动。真实航天器并不处在理想二体引力场中,地球非球形引力、大气阻力、太阳光压、第三体引力等都会使轨道根数发生持续或周期性的偏移。把这些偏移按时间尺度分解,就得到长期项、长周期项和短周期项三类分量。理解这三类项的物理来源与作用规律,是设计保持策略、选择控制时机和估算推进剂消耗的前提。
一、三类项的物理含义与时间尺度
设某轨道根数 σ(t) 在摄动下的变化可写成:
σ(t) = σ₀ + σ̇·t + A_lp·sin(ω_lp·t + φ_lp) + A_sp·sin(ω_sp·t + φ_sp)
其中 σ₀ 为初始值,σ̇·t 是长期项,随时间线性累积;第二项是长周期项,其周期与近地点幅角 ω 或升交点赤经 Ω 的慢变化相关,通常为几天到几个月;第三项是短周期项,周期与轨道周期同量级或为其整数分之一,通常为几十分钟到几小时。
三类项的关键区别不在幅度大小,而在是否累积。短周期项和长周期项在一个周期内正负抵消,不产生净漂移;长期项则每圈都朝同一方向叠加,即使每圈增量极小,经过成百上千圈后也会显著偏离。因此轨道保持真正要对抗的是长期项,而周期项主要影响瞬时精度和测轨采样。
二、地球非球形引力:J2 项的主导作用
对近地轨道,地球扁率引起的 J2 摄动是最主要的漂移来源。J2 对轨道根数的影响具有明确的选择性:
- 升交点赤经 Ω:产生长期漂移,速率与轨道半长轴 a、倾角 i、偏心率 e 相关。对近圆轨道,Ω̇ ≈ -1.5·J2·(R_e/p)²·n·cos i,其中 p = a(1-e²),n 为平均角速度。倾角小于 90° 时 Ω 持续西退,这正是太阳同步轨道需要特定倾角的原因。
- 近地点幅角 ω:同样产生长期漂移,速率约为 Ω̇·(2 - 2.5·sin²i)/cos i 量级,决定近地点在轨道面内的旋转方向。
- 半长轴 a、偏心率 e、倾角 i:J2 不产生长期项,只引起短周期和长周期振荡。这意味着纯 J2 不会让轨道整体升高或降低,但会让 a、e、i 在标称值附近周期性摆动。
以高度约 700 km、倾角 98° 的太阳同步轨道为例,J2 引起的 Ω̇ 约为 +0.9856°/天,恰好补偿地球绕日公转的平均角速度,使升交点地方时保持不变。若倾角偏离设计值 0.1°,Ω̇ 的变化量约为 0.0015°/天,一年累积约 0.55°,对地方时要求严格的遥感任务已不可忽略。
三、大气阻力与太阳光压
大气阻力对低轨航天器的作用表现为持续的切向减速,直接导致半长轴 a 长期衰减,进而使轨道周期缩短、平均高度下降。其衰减速率与面质比、大气密度和速度相关,在太阳活动高年可增大数倍。阻力同时使偏心率缓慢减小,对近圆轨道影响有限,但对大椭圆轨道会明显改变近地点高度。
太阳光压对半长轴的影响取决于面质比和光照条件。对高面质比航天器,光压可产生不可忽略的长期项,尤其在轨道面与太阳方向夹角变化时,会调制偏心率矢量,形成以年为周期的长周期项。对低轨小卫星,光压通常小于大气阻力;对地球同步轨道和深空探测器,光压则成为主要摄动之一。
四、第三体引力与共振效应
太阳和月球的第三体引力对高轨目标影响显著。对地球同步轨道,日月引力引起轨道倾角的长期漂移,典型速率约 0.75°/年,若不控制,十几年后倾角将超出 5° 的保持窗口。同时,地球非球形引力的田谐项在同步轨道高度与轨道周期形成共振,可激发偏心率的长周期增长,周期约数十年。
对特定轨道,如 12 小时周期的半同步轨道,J2 与日月引力可能形成组合共振,使某些根数出现更复杂的长期与长周期耦合。这类效应在任务设计阶段就应通过数值积分评估,而不是等到在轨后再处理。
五、对轨道根数影响的综合判断
把上述机制汇总,可以得到一张按根数分类的影响关系图。下图以近地轨道为主,标出各摄动源对六个轨道根数的主要作用类型:L 表示长期项,LP 表示长周期项,SP 表示短周期项,— 表示影响可忽略。
从图中可以读出几个对工程实践有直接指导意义的结论。第一,近地轨道保持的首要矛盾是 Ω 和 ω 的长期漂移,以及大气阻力引起的 a 衰减。第二,e 和 i 在近地轨道上以周期项为主,控制需求相对较低;但在同步轨道上,i 的长期漂移成为主要控制对象。第三,同一摄动源对不同根数的作用类型不同,不能笼统地说“摄动使轨道漂移”,而应逐根数分析。
六、从机制到保持策略
明确三类项之后,轨道保持策略就有了清晰的着力点。对长期项,需要在漂移累积到容差边界前施加反向修正,修正量按 σ̇·Δt 估算,控制时机通常选在漂移方向与推力方向匹配的弧段。对长周期项,可在设计阶段通过选择轨道参数使周期振幅落在容差内,或在振幅峰值附近预留控制余量。对短周期项,一般不做专门控制,而是在测轨定轨和门限判断中通过滤波或取平均来抑制其影响。
一个实用的判断流程是:先由任务约束确定各根数的容差带,再用摄动模型计算各根数的长期漂移率,得到“不控制可维持时间”;若该时间小于任务周期,则必须安排保持机动。以 700 km 太阳同步轨道为例,若 Ω 容差为 ±0.1°,J2 引起的 Ω̇ 偏差约 0.0015°/天,则约 67 天就会触及边界,需要每两个月左右进行一次轨道面修正。这个量级可以直接用于推进剂预算和测控计划编排。
需要强调的是,长期项与周期项的划分依赖于所取的时间尺度和参考模型。同一物理效应在不同轨道高度、不同倾角下可能表现出不同的主导类型。因此在实际任务中,应以数值积分或解析摄动模型给出的根数时间序列为依据,而不是仅凭经验判断。只有把摄动机制、根数响应和任务容差三者对齐,才能把轨道保持从“被动纠偏”变成“可预测、可规划”的常规操作。