空间碎片环境演化模型 — 地球卫星种群不稳定性

论文信息： J.-C. Liou and N. L. Johnson, "Earth Satellite Population Instability: Underscoring the Need for Debris Mitigation," NASA Johnson Space Center, 2008.

研究机构： NASA 轨道碎片项目办公室 (NASA Orbital Debris Program Office, Johnson Space Center)

1. 引言：Kessler 综合征与碎片环境演化

空间碎片问题并非一夜之间出现。自 1957 年 Sputnik 1 进入轨道以来，人类的空间活动已在近地轨道（LEO）积累了超过 9,000 个可编目的在轨物体，总质量超过 500 万公斤。然而，最令人忧虑的并非当前已被跟踪的碎片数量，而是碎片数量的未来演化趋势。

1978 年，NASA 科学家 Donald J. Kessler 与 Burton G. Cour-Palais 在一篇开创性论文中提出了一个令人不安的预测：当近地轨道某些高度区间的物体空间密度超过某一临界空间密度（critical spatial density）时，碰撞产生的新碎片数量将超过因大气阻力衰减而离轨的碎片数量。此时，即使人类停止一切发射活动，碎片数量也会因碰撞级联（collisional cascading）而持续增长——这一过程后来被称为 Kessler 综合征（Kessler Syndrome）。

Liou 与 Johnson 于 2008 年发表的这项研究，利用 NASA 轨道碎片项目办公室开发的高保真三维演化模型 LEGEND（LEO-to-GEO Environment Debris model），首次系统性地论证了一个严峻的事实：即使从现在起完全停止所有发射，仅凭现有在轨物体之间的随机碰撞，2050 年后 LEO 碎片数量仍将持续增长。

2. LEGEND 模型：碎片环境的三维演化模拟

2.1 模型概述

LEGEND 是一个高保真三维物理模型，能够模拟从 LEO 到 GEO（地球静止轨道）的碎片环境历史演变与未来预测。与早期的简化分析模型不同，LEGEND 具备以下核心能力：

历史环境重建：基于美国空间监视网络（SSN）的编目数据，校准模型以匹配 1957 年至模拟起始时刻的已知碎片环境；
未来爆炸模拟：基于历史爆炸统计规律，采用蒙特卡洛方法模拟在轨物体的未来爆炸事件；
碰撞概率与碎片生成：计算轨道物体之间的碰撞概率，并利用 NASA 标准解体模型（Standard Breakup Model）生成碰撞产生的碎片云；
轨道衰减：考虑大气阻力、太阳辐射压等摄动力，模拟碎片的自然衰减过程。

2.2 蒙特卡洛方法

LEGEND 的未来预测采用蒙特卡洛方法。在每个时间步内，程序执行以下步骤：

对每一个完整的在轨物体（intact object），根据其类型（火箭箭体、航天器）和历史爆炸率，估计其在该时间步内发生爆炸的概率 $P_{\text{exp}}$；
生成一个 $[0, 1]$ 区间均匀分布的随机数 $r$，若 $r < P_{\text{exp}}$，则触发爆炸事件；
对每一对可能发生碰撞的目标-炮弹组合，计算其碰撞概率 $P_{\text{col}}$，以类似方式判定是否发生碰撞；
若发生碰撞，判断碰撞类型（灾难性碰撞或非灾难性碰撞），并生成相应碎片。

关键判据：当碰撞能量与目标质量的比值超过 $40 \text{ J/g}$ 时，定义为灾难性碰撞（catastrophic collision），目标完全解体。这一阈值源自地面超高速撞击实验。

由于蒙特卡洛方法的随机性，LEGEND 对每个场景执行 50 次独立模拟（每次模拟覆盖 200 年），最终结果取统计平均值，标准误差约 5%。

3. 核心发现：种群不稳定性

3.1 关键假设

为揭示碰撞级联的 "pure" 效应，研究设置了保守假设：

零未来发射：2004 年 12 月后，不再发射任何火箭箭体或航天器；
无主动离轨：现有航天器不执行为期 25 年的离轨处置机动；
自然衰减：仅依赖大气阻力等自然力清除碎片。

在这样的极端假设下，如果碎片数量仍然增长，则说明 Kessler 综合征已经或即将不可逆转。

3.2 全 LEO 种群演化

Fig 1

图 1：碎片种群的增长预测。 LEO 区域（200–2000 km）中大于 10 cm 物体的有效数量。有效数量定义为物体在每个轨道周期内停留在 200–2000 km 高度区间的时间比例。"完整物体"指尚未发生解体的火箭箭体和航天器。（来源：LEGEND 模拟，Liou & Johnson, 2008）

图 1 揭示了碎片演化的三个关键阶段：

1957–2000 年（历史阶段）：总碎片数量快速增长，主要由爆炸碎片（绿色曲线）和完整物体（蓝色曲线）贡献。爆炸碎片的峰值出现在 1990 年前后，反映了冷战时期大量火箭箭体在轨爆炸的历史事实；
2000–2050 年（过渡阶段）：爆炸碎片开始缓慢减少（火箭箭体消耗殆尽），完整物体数量因自然衰减而下降。然而，碰撞碎片（红色曲线）从接近零开始加速增长；
2050–2210 年（碰撞主导阶段）：碰撞碎片在 2060 年左右超过完整物体数量，在 2080 年左右超过爆炸碎片数量，成为 LEO 碎片的主导来源。到 2210 年，仅碰撞碎片一项就接近 7,000 个（大于 10 cm），推动总碎片数突破 11,000。

关键结论：即使在零发射的最保守假设下，碰撞级联驱动的碎片增长依然是不可避免的。

3.3 空间密度分布的演变

Fig 2

图 2：碎片环境的未来空间密度分布。 三个不同年份（2004、2104、2204）中，直径大于 10 cm 物体的空间密度随高度分布。900–1000 km 高度区间的急剧增长源自该区域预测的高碰撞活动率。

图 2 给出了更精细的空间分布视角。三条曲线分别对应 2004 年（基线）、2104 年和 2204 年的空间密度分布：

基线（2004 年，蓝色）：空间密度在 800 km 和 1400 km 附近各有一个峰值，这是历史上最常使用的轨道高度；
100 年后（2104 年，红色）：900–1000 km 高度区间出现显著的空间密度尖峰（接近 $3.8 \times 10^{-8} \text{ obj/km}^3$），远超基线水平；
200 年后（2204 年，绿色）：同一区间的空间密度继续攀升至近 $5.0 \times 10^{-8} \text{ obj/km}^3$。

900–1000 km 高度区间之所以成为碎片增长的"热点"，原因有三：
1. 该区域大气极其稀薄，碎片轨道寿命长达数百年至数千年；
2. 历史上大量航天器和火箭箭体部署于此；
3. 高空间密度导致碰撞概率显著升高，触发级联效应。

3.4 "红色区域"：900–1000 km

Fig 3

图 3："红色区域"。 900–1000 km 高度区间内大于 10 cm 物体的有效数量。该区域的种群增长是驱动整个 LEO 碎片环境恶化的主要原因。

图 3 聚焦于 900–1000 km 这个最危险的高度区间。与全 LEO 趋势相比，该区域表现出更早、更剧烈的增长：

总碎片数量从 2010 年开始即呈现加速增长态势（全 LEO 趋势直到 2050 年后才明显）；
碰撞碎片的贡献从一开始就占据主导地位；
到 2210 年，该区域大于 10 cm 的碎片数量超过 3,000 个。

Liou 与 Johnson 将这一高度区间称为"红色区域"（the red zone），因为它是整个 LEO 碎片环境中最不稳定、最危险的部分。

4. 数学框架：临界空间密度

4.1 Kessler 的原始判据

Kessler 在 1991 年的后续研究《Collisional cascading: The limits of population growth in low Earth orbit》中将临界空间密度概念形式化。考虑高度 $h$ 处的一个体积元 $V$，其中包含 $N$ 个碎片物体。碎片的平均碰撞截面为 $A_c$（数量级为 $10^{-6} \text{ km}^2$），平均相对速度为 $v_{\text{rel}}$。

碎片数量的演化方程可写为：

$$\frac{dN}{dt} = \dot{N}_{\text{launch}} + \dot{N}_{\text{explosion}} + \dot{N}_{\text{collision}} - \dot{N}_{\text{decay}}$$

其中碰撞产生率与空间密度的平方成正比：

$$\dot{N}_{\text{collision}} \propto \rho^2 \cdot A_c \cdot v_{\text{rel}}$$

这里 $\rho = N/V$ 为空间密度。关键洞察在于：衰减率 $\dot{N}_{\text{decay}}$ 在大气稀薄的高轨道区域很小，而碰撞产生率随 $\rho^2$ 增长。一旦 $\rho$ 超过临界值 $\rho_c$，即使 $\dot{N}_{\text{launch}} = 0$ 且 $\dot{N}_{\text{explosion}} = 0$，也必然有 $dN/dt > 0$。

4.2 LEGEND 中的碰撞概率计算

在 LEGEND 模型中，两物体在单位时间内的碰撞概率采用"气体动力学"类比：

$$P_{ij} = s_i s_j \cdot V_{\text{imp}} \cdot \frac{dt}{V_{\text{bin}}}$$

其中：
- $s_i, s_j$ 为物体 $i$ 和 $j$ 的碰撞截面半径；
- $V_{\text{imp}}$ 为两物体的平均撞击速度；
- $V_{\text{bin}}$ 为空间网格体积。

LEGEND 将 LEO 区域划分为数百个高度-纬度-经度网格单元，在每个单元内独立计算碰撞概率，从而实现了空间分辨的碎片演化模拟。

5. 对碎片减缓策略的启示

Liou 与 Johnson 的研究得出了几个对空间碎片减缓政策具有深远影响的结论：

5.1 仅靠减缓措施不够

长期以来，国际航天界主要关注的是减缓（mitigation）策略，即限制新碎片的产生。这些措施包括：
- 25 年离轨规则（LEO 航天器任务结束后 25 年内再入大气层）；
- 钝化措施（任务结束后排空剩余推进剂和电池能量，防止在轨爆炸）；
- 减少任务相关碎片（如分离螺栓、镜头盖等）的释放。

然而，LEGEND 的模拟结果表明，即使上述减缓措施被完全执行，也仅能延缓但不能阻止碎片种群的持续增长。原因在于，现有轨道上已经积累了足够的"种子"碎片，碰撞级联已经开始。

5.2 主动清除的必要性

论文明确指出，主动碎片清除（Active Debris Removal, ADR）应成为未来碎片治理的核心方向。具体建议包括：

优先清除大质量物体：大型火箭箭体和废弃航天器（质量超过 1,000 kg）是未来碰撞碎片的主要来源；
聚焦"红色区域"：900–1000 km 高度区间应作为 ADR 的优先目标区域；
尽快行动：每延迟一年，碰撞级联效应就会进一步增强，清除成本也随之增加。

6. 总结与展望

Liou 与 Johnson (2008) 的这项研究是空间碎片环境演化领域的分水岭。通过 LEGEND 模型的高保真蒙特卡洛模拟，研究给出了一个无法回避的结论：即使立即停止所有发射活动，Kessler 综合征的触发已经不可避免。

该研究的主要贡献可归纳为：

定量论证了碎片种群的不稳定性：零发射假设下，碰撞碎片将在 2050 年后驱动总碎片数持续增长；
识别了"红色区域"：900–1000 km 高度区间是碰撞级联最先触发、最严重的区域；
推动了政策转变：从仅关注"减缓"到呼吁"主动清除"。

时至今日（2026 年），随着 Starlink、OneWeb 等巨型星座的快速部署，LEO 的空间密度进一步增加。Liou 与 Johnson 的警告非但没有过时，反而比以往任何时候都更加紧迫。对于中国航天工程专业的研究生而言，深入理解碎片环境演化机制，不仅是学术研究的需要，更是未来轨道设计、碰撞规避和空间交通管理实践中的必备知识。

参考文献：

Kessler, D. J., & Cour-Palais, B. G. (1978). Collision frequency of artificial satellites: The creation of a debris belt. Journal of Geophysical Research: Space Physics, 83(A6), 2637–2646.
Kessler, D. J. (1991). Collisional cascading: The limits of population growth in low earth orbit. Advances in Space Research, 11(12), 63–66.
Liou, J.-C., & Johnson, N. L. (2006). Risks in space from orbiting debris. Science, 311(5759), 340–341.
Liou, J.-C., & Johnson, N. L. (2008). Earth satellite population instability: Underscoring the need for debris mitigation. NASA Johnson Space Center.
Liou, J.-C., Hall, D. T., Krisko, P. H., & Opiela, J. N. (2004). LEGEND — A three-dimensional LEO-to-GEO debris evolutionary model. Advances in Space Research, 34(5), 981–986.

【经典论文解读】空间碎片环境演化模型 — 地球卫星种群不稳定性

目录