3. 空间辐射环境建模与仿真工具

要完成典型轨道的辐射环境仿真，核心不是“打开某个软件点几下”，而是先明确要算的物理量，再选择能给出该物理量的模型，最后把轨道参数、屏蔽厚度和任务时间正确映射到模型输入上。本节围绕四类常用工具展开：AP8/AP9 与 AE8/AE9 负责俘获带粒子通量，CREME96 负责银河宇宙线与太阳粒子事件的重离子环境及单粒子效应估算，SPENVIS 则作为集成平台把轨道传播、各模型调用和屏蔽计算串起来。

一、先分清四类工具的分工

AP8/AE8 是静态、长期平均的俘获带模型，输出的是给定轨道上质子（AP）和电子（AE）的积分或微分通量，适合做总剂量和位移损伤的初步评估。AP9/AE9 是新一代概率模型，除均值外还能给出不同置信水平（如 50%、95%）的通量，适合需要评估环境不确定度的任务。CREME96 处理的是高能重离子，输出 LET 谱和积分通量，用于单粒子效应和总剂量中的重离子贡献。SPENVIS 本身不产生新物理模型，而是把上述模型与轨道传播、屏蔽输运（如 SHIELDOSE-2、MULASSIS）组合成一条可复现的计算链。

二、AP8/AE8 与 AP9/AE9 的实际使用

以一条 800 km 太阳同步轨道为例，倾角约 98.7°，任务周期 3 年。若用 AP8/AE8，需要先确定太阳活动条件：AP8 提供太阳活动最大（MAX）和最小（MIN）两套静态通量，工程上常取两者包络或按任务期太阳活动相位加权。在 SPENVIS 中，先输入轨道六根数，选择“Trapped Radiation”模块，指定 AP8 的 MAX 或 MIN，再设定能量网格（例如质子 0.1–400 MeV，电子 0.04–7 MeV），平台会沿轨道积分给出积分通量谱。

AP9/AE9 的输入多了一个关键参数——置信水平。以 AE9 为例，选择“Mean”“95th percentile”或“99th percentile”，模型会基于轨道采样和太阳活动概率分布给出对应分位通量。对 800 km 轨道，电子通量对太阳活动相位较敏感，若任务要求高可靠，建议同时输出均值和 95% 分位，用两者差值评估环境不确定度。需要注意，AP9/AE9 的轨道采样点数、能量分档和屏蔽厚度都会影响结果，配置时应记录这些参数，便于复现。

三、CREME96 的重离子环境计算

CREME96 的典型用法是：选择银河宇宙线（GCR）或太阳粒子事件（SPE）模式，指定太阳活动调制参数（如 GCR 的太阳调制势 φ，常用 400–1200 MV），设定屏蔽厚度（铝等效，例如 100 mils ≈ 2.54 mm），然后输出 LET 谱和积分通量。以 100 mils 铝屏蔽为例，GCR 的 LET 谱在 1–10 MeV·cm²/mg 区间仍有可观通量，而 SPE 在太阳活动高年可能高出几个数量级。工程上常把 CREME96 输出的 LET 谱与器件截面曲线卷积，估算单粒子翻转率。

四、SPENVIS 集成仿真的操作要点

SPENVIS 的价值在于把上述模型串成一条链。实际操作顺序通常是：新建项目 → 输入轨道 → 选择环境模型（AP8/AE8 或 AP9/AE9、CREME96）→ 设置屏蔽厚度和材料 → 运行计算 → 导出通量谱和剂量结果。以 800 km 轨道、3 年任务、2.54 mm 铝屏蔽为例，可先运行 AP8/AE8 得到质子与电子积分通量，再用 SHIELDOSE-2 计算屏蔽后总剂量；同时调用 CREME96 得到重离子 LET 谱，用于单粒子效应评估。导出结果时，建议保留能量网格、屏蔽厚度、太阳活动条件和模型版本，这些是后续复现和对比的基础。

五、结果判读与常见问题

拿到通量谱后，先看量级是否合理：800 km 太阳同步轨道上，>10 MeV 质子积分通量在 3 年任务中通常在 10⁹–10¹⁰ cm⁻² 量级（示例范围，实际随太阳活动和屏蔽变化）。若结果明显偏离，优先检查轨道参数单位（km 还是 m）、能量单位（MeV 还是 keV）、屏蔽厚度单位（mils 还是 mm）以及太阳活动条件是否与任务期匹配。另一个常见问题是把 AP8/AE8 的 MAX 结果直接当作任务期平均值，导致剂量高估；正确做法是按任务期太阳活动相位加权，或直接用 AP9/AE9 的概率模型。

最后，所有模型结果都应视为工程估算而非实测。AP8/AE8 基于历史卫星数据，AP9/AE9 虽更新了数据但仍存在模型不确定度；CREME96 的 SPE 统计也依赖太阳活动假设。因此，仿真结果应结合在轨监测数据或地面试验进行校核，并在设计裕度中体现模型不确定度。