《Aerospace》：An SSA-to-ADR Uncertainty Propagation Pipeline for Tether-Net Debris Capture Mission Planning
编辑推荐：
绳网捕获与空间态势感知（Space Situational Awareness, SSA）一直作为两个独立的学科领域发展，而捕获动力学研究几乎普遍预设一个确定性的目标状态。本文提出一种端到端的SSA到主动碎片清除（Active Debris Removal,
绳网捕获与空间态势感知（Space Situational Awareness, SSA）一直作为两个独立的学科领域发展，而捕获动力学研究几乎普遍预设一个确定性的目标状态。本文提出一种端到端的SSA到主动碎片清除（Active Debris Removal, ADR）不确定性传播管道，将绳网捕获视为一个以观测数据质量为条件的概率决策问题。该管道耦合了三个阶段：从双行轨道根数（TLE）和光度测光光变曲线数据表征碎片状态协方差；通过状态转移矩阵和阻力驱动的过程噪声将该协方差传播至弹射时刻；以及在弹射参数空间上构建捕获概率曲面。将该框架应用于Telkom-3（NORAD 38744）在11个历元上的数据，结果表明SSA数据质量通过预测状态的有效期（validity horizon）而非末端导航精度来制约捕获性能。两个分支相差两到三个数量级：平动有效期为8至9小时（8.0–9.3 h），而转动有效期在三个分析历元上从115 s降至23 s，与此同时自旋周期从20.3 s缩短至11.8 s。在30 s观测到弹射延迟下，捕获概率相应地从0.90降至0.36。因此，光度测光支持必须在弹射前数分钟而非数小时获得，且弹射时机需精确到预测有利转动相位约1秒以内。
论文解读
一、研究背景与问题
低地球轨道（LEO）中失效卫星和碎片碎片的累积已达到临界阈值，凯斯勒综合征（Kessler syndrome）对运行中的空间基础设施构成系统性风险。主动碎片清除（Active Debris Removal, ADR）被视为必要补充手段，其中专用航天器对高优先级目标进行物理捕获并使其离轨。执行ADR任务需要两种在实践中独立发展的能力：空间态势感知（Space Situational Awareness, SSA）——通过地基光学和雷达观测表征目标的轨道与转动状态；以及捕获动力学工程——设计负责物理回收的硬件与控制系统。两者的相互依赖是操作性的而非偶然性的——捕获尝试的质量受制于任务前表征的保真度——然而尚无已发表的框架将两者定量耦合为单一的任务规划流程。
在研发中的捕获技术中，绳网系统因其三个特性而受到持续关注：第一，网可被动耗散目标的刚体转动角动量；第二，网顺应目标形状而非要求指定接口，因此捕获可容忍显著的几何不确定性；第三，绳连接是柔顺的，可容纳捕获后运输中两质心之间较大且不完全已知的偏移。这些特性对无法提供遥测、无法机动协助捕获或稳定姿态的失效卫星尤为重要。
然而，现有研究存在一个普遍局限：目标位置、速度和姿态均被预设为确定性且精确已知的输入。这在操作上不可行，因为失效LEO卫星携带的TLE衍生位置不确定性达公里级，且随预测间隔的平方增长；其姿态状态必须从有限节奏和精度的地基光度或雷达观测中推断。已知翻滚速率允许确定性响应，而未知的翻滚速率则使捕获成功成为概率而非结果。所缺失的是从SSA网络实际掌握的目标信息——以元素集年龄、光学过境节奏和光度测量噪声导出的状态协方差表示——到应据此选择的弹射参数的正式路径。
二、研究内容与主要结论
本文提出一种端到端的SSA到ADR不确定性传播管道，将绳网捕获视为以可用观测质量为条件的概率决策问题。管道耦合三个阶段：第一阶段从TLE导出的轨道不确定性和光度光变曲线自旋状态估计表征完整的碎片状态协方差（位置、速度、自旋速率和主轴方向）；第二阶段通过状态转移矩阵和阻力驱动的过程噪声将其传播至网弹射历元；第三阶段将传播后的不确定性转化为捕获概率——即网周长在接触时包围目标投影轮廓至少指定比例的采样碎片状态比例——在弹射参数空间上评估以生成捕获概率曲面（Capture Probability Surface）。
框架以Telkom-3（NORAD 38744）为方法学测试案例加以开发和验证。该卫星于2012年8月因上面级故障被困于低椭圆轨道，2021年2月再入。其适用性依赖于两个罕见同时可用的数据集：跨越2016年12月至2020年7月的11个历元TLE历史记录，以及来自瑞士齐默瓦尔德天文台（AIUB）的同期光度自旋状态测量。
主要结果表明，状态估计的平动与转动两个分支并非同等约束。每个分支携带一个有效期（validity horizon），即预测状态保持足够精确以据此规划的区间，两者相差两到三个数量级：平动有效期为数小时，而转动有效期为数十秒且随目标自旋加速而缩短。转动分支因此构成约束，随之而来的观测要求是光度延迟而非跟踪精度。
三、关键技术方法
研究人员构建了三阶段耦合管道。第一阶段利用TLE元素集历史记录通过局部线性拟合获取平均运动导数，结合23.0%的阻力系数变异系数建立沿轨不确定性模型，径向和法向分量按固定比例设定；转动状态采用齐默瓦尔德天文台多历元光变曲线分析的自旋周期和主轴方向估计。第二阶段通过沿参考轨迹积分变分方程获得状态转移矩阵，采用二体和J2带谐加速度的引力梯度雅可比，过程噪声由阻力弥散驱动。第三阶段采用蒙特卡洛积分，以圆透镜面积公式计算包围分数，在弹射参数空间上以625个网格点、每点5×10?次实现评估捕获概率。轨道确定更新采用平方根信息滤波器的批处理最小二乘估计器，模拟齐默瓦尔德和蒂马乌两个站点的角度观测。
四、研究结果
4.1 转动状态及其不确定性
光度测定的自旋周期从2016年12月的28.212 s单调递减至2020年7月的11.830 s，而1σ估计不确定度从0.20 s增长至0.50 s。两趋势叠加：自旋周期缩短1.7倍而σP上升1.7倍，使σP/P增加4.9倍，固定质量的光度解在观测活动末期失去可用相位信息的速度比初期快近5倍。
4.2 位置不确定性增长
沿轨不确定性在一天元素集年龄下为E1、E2、E3分别6.62、8.91和8.73 km。排序并非随历元单调，因为位置不确定性由阻力的不可预测部分驱动，而该部分追踪太阳活动而非仅追踪轨道高度。这些值与3 m网半宽之间超过三个数量级的差距定义了两尺度问题。
4.3 跨观测场景的轨道确定
形式协方差在任一跟踪可用时从公里级坍缩至米级——沿轨标准差改善因子达700至3700。增加蒂马乌站使几何观测计数仅增加22%至41%，却将后验沿轨不确定性改善2.75至4.06倍，表明第二站的回报由其与第一站的几何分离而非观测时间决定。
4.4 协方差传播作为一致性检验
跟踪场景在一天规划地平线上以0.11%以内的精度复现解析表达式，验证了传播实现保持相位模式的性质。未跟踪场景S1达到34.8至42.4 km，为跟踪案例的4.8至5.3倍，主要由分配的历元协方差的径向分量驱动。由于阻力不确定性过程噪声在单日内主导，S3相对S2的历元时间优势在弹射历元到达前已被抹除。
4.5 SSA有效期：转动分支构成约束
转动有效期比平动有效期在每个历元短两到三个数量级，且随自旋加速单调缩短：E1为114.9 s，E2为53.7 s，E3为23.3 s。平动有效期在8至9小时几乎平坦。对Telkom-3而言，SSA数据有效性的约束条件在整个过程中是转动性的，且在三个分析历元间收紧4.9倍。即使在极端容差下，两分支仍保持1.5至2个数量级的差距。操作后果直接：支持捕获尝试的光度解必须在弹射前数分钟而非数小时获得，且弹射需精确到预测有利相位约1秒以内。
4.6 捕获概率曲面
在45 s观测到弹射延迟下，峰值捕获概率从E1的0.821降至E2的0.514和E3的0.251，由接触相位弥散增长驱动。每个历元从约0.95的共同相位已知值衰减至约0.20的共同均匀相位下限，历元维度仅控制衰减速率：边界由轮廓几何和末端导航精度设定，两者均无法由任何观测活动改变，而给定延迟下的位置完全由σP/P设定。操作包络在每个阈值下随历元单调收缩；在0.80阈值下仅第一历元非零，至第三历元完全消失。
4.7 光度表征何时重要
相位知识收益图表明，对低纵横比目标，相位知识收益有限，因为均匀相位下捕获概率已高；对高纵横比目标，收益显著。收益随网半宽增大而增加，因为更大的网对相位窗口的敏感性降低。对于Telkom-3的部署构型（纵横比9.16），相位知识在3 m网下收益约0.4，在6 m网下收益约0.7。
五、讨论与结论
讨论部分指出，捕获概率对包围阈值和导航残差的敏感性在均匀相位极限下均较弱。参数化曲面的操作包络面积对展开规律的具体形状不敏感，因为所有捕获概率均在最优锥半角处评估，该值对展开规律不变。研究人员还讨论了缩减阶展开模型的验证边界：已验证配置涵盖1.0和3.2 m网边长对应0.15和10–50 kg目标，6 m网和1903 kg目标超出该包络，一致性检查针对函数形式而非绝对量级。
研究结论为：SSA数据质量通过预测状态的有效期而非末端导航精度制约捕获性能。平动与转动分支的有效期相差两到三个数量级：平动为8至9小时，转动从115 s降至23 s。捕获概率在30 s延迟下从0.90降至0.36。光度测光支持必须在弹射前数分钟而非数小时获得，弹射时机需精确到预测有利转动相位约1秒以内。通过将光度SSA输出作为弹射参数选择的正式概率输入，该框架将ADR任务规划从确定性力学模拟转向以可量化观测质量为条件的概率任务分析。
