两阶段几何-动力学方法用于实时空间编目碎片通量估计

《Acta Astronautica》：Two-stage geometric-dynamic approach for real-time space catalogued debris flux estimation

　　

空间碎片通量是量化碰撞风险的关键指标，然而现有方法面临根本性权衡：统计环境模型可提供长期全球平均通量场，但缺乏实时响应能力；全编目确定性轨道传播虽物理精度高，但计算成本过高。本研究提出一种面向编目碎片群的目标中心通量估计的两阶段几何-动力学框架，采用双行根数（TLE）数据。第一阶段，在目标周围构建地心空心球壳，并利用快速径向区间相交检测保守筛选出不相交轨道，借助几何稀疏性缩减候选集。第二阶段，仅对保留的候选碎片选择性执行简化通用摄动-4（SGP4）模型传播，评估时间分辨的控制球穿越事件，并计算年度化、表面归一化的通量密度。该方法的计算复杂度表现为一次线性筛选加上稀疏加权传播工作量，在保持全传播物理决策逻辑的同时显著加速运算。实验表明，该方法在检测的穿越事件及衍生的通量密度方面与全传播结果高度一致，并在所测试的计算条件下支持快速局部目标级风险评估。

研究背景与意义

近二十年来，近地空间中的驻留空间物体（RSO）数量呈指数增长，导致轨道拥堵与碰撞风险急剧上升。据欧洲航天局（ESA）2025年《空间环境报告》统计，目前地球轨道上已有超过44,300个编目碎片，速度超过7.5 km/s，毫米级颗粒即可对航天器造成严重损伤，并可能触发凯斯勒综合征（Kessler Syndrome）。空间碎片通量作为衡量碎片风险的核心指标，定义为单位时间穿过单位面积的碎片穿越次数，广泛应用于防护设计、撞击风险评估及模型校准。现有经验模型如NASA的ORDEM 3.1与ESA的MASTER-8，可提供基于观测与物理源模型的长期统计通量分布，但其依赖预计算数据库，缺乏空间态势感知（SSA）与空间交通管理（STM）所需的实时性。另一方面，基于简化通用摄动-4（SGP4）模型的全编目轨道传播可提供高物理保真度，但随着碎片数量增加，轨道-区域相交检测的计算复杂度达到O(n²)，导致运算成本过高，难以在日常操作中实施。因此，如何在保持物理精度的同时显著提升计算效率，成为当前空间碎片环境建模的重大挑战。

主要技术与方法

研究人员针对编目碎片群体提出两阶段框架。第一阶段引入基于地心空心球壳的几何预筛选机制，通过保守的径向区间相交检测剔除无法与目标高度层相交的轨道，从而利用空间稀疏性大幅缩减候选碎片数量。该步骤不同于传统成对近距离接近筛查，而是面向目标中心的通量估计任务。第二阶段对筛选后的碎片采用SGP4摄动模型进行选择性轨道传播，考虑地球扁率项J₂与大气阻力效应，计算目标控制球内的时间分辨通量。技术路线将几何筛选与时间动态解耦，实现亚公里级空间精度并保持TLE数据的可追踪性。实验基于公开TLE数据，以全SGP4传播（FSP）为基准，在不同仿真时长、时间分辨率与筛选阈值下进行性能评估。

研究结果

概念框架

研究人员将通量定义为编目碎片穿越目标控制球的时间积分穿越率，而非全尺寸环境通量。该框架兼顾计算效率与物理保真度，适用于局部目标中心的风险估计。

实验

在多种仿真条件下，两阶段框架均显著减少需要传播的碎片数量，同时保持与全传播方法在检测到的穿越事件及通量密度上的一致性。计算耗时随筛选阈值的优化进一步降低，验证了方法的加速能力与可扩展性。

结论

该方法在相同的TLE/SGP4假设下保留了全传播决策结果，同时通过几何筛选与选择性传播显著降低运算成本。它并非替代MASTER类统计模型，而是作为其补充工具，用于实时、确定性的目标中心通量分析。

讨论与结论翻译

本研究提出的两阶段几何-动力学框架，可在避免全编目传播高昂成本的同时，实现时间分辨的确定性局部通量评估。研究结果表明，这种重构方法在相同输入条件下能够保持全传播的决策结果，并显著提高计算效率。其意义在于为目标级空间态势感知提供了一种兼具物理精度与实时性的工具，可与统计环境模型联合使用，以提升轨道安全管理与可持续运行能力。