📚 课程目录第2章空间辐射环境

2. 空间辐射环境

单粒子效应评估的第一步，不是选器件，而是把器件将要面对的空间辐射环境量化清楚。空间辐射环境由三类主要来源构成：银河宇宙线（GCR）、太阳粒子事件（SPE）和地球俘获带（辐射带）。三者的粒子种类、能谱形状、时间行为与空间分布差异极大，因此不能用一个统一的“总剂量”或“总通量”来概括，而必须分别建立能谱模型，再折算到具体轨道、具体屏蔽条件下的在轨通量。

三类辐射源的物理特征

银河宇宙线来自太阳系外，由高能带电粒子组成，其中约85%是质子、约14%是α粒子，其余为重离子（常以HZE表示，即高能重带电粒子）。GCR的特点是能量极高（可到10^20 eV量级）、通量低但穿透力强，且受太阳活动调制：太阳活动高年，日球层磁场增强，低能段GCR被削弱，通量下降；太阳活动低年则相反。对单粒子效应而言，GCR中的重离子是导致高可靠器件发生SEU、SEL的关键因素，尤其是铁、氪、氙等高LET粒子。

太阳粒子事件由太阳耀斑和日冕物质抛射驱动，以质子为主，也包含α粒子和少量重离子。SPE是偶发事件，持续时间从数小时到数天，通量可比GCR背景高几个数量级。其能谱通常用幂律或指数形式描述，低能段粒子多、高能段迅速衰减。SPE对单粒子效应的威胁在于：在事件期间，质子通量骤增，高能质子可通过核反应在器件敏感区产生次级重离子，从而触发SEU；同时SPE也会显著抬高总剂量。

地球俘获带由地磁场捕获的质子和电子组成，分为内带和外带。内带以高能质子为主，峰值能量可达数百MeV，主要分布在赤道附近低高度区域；外带以电子为主。对单粒子效应评估而言，俘获带质子是低地球轨道（LEO）任务中必须考虑的因素，尤其是南大西洋异常区（SAA），该区域磁场较弱，俘获带质子下沉到较低高度，使经过该区域的航天器质子通量显著升高。

能谱模型与在轨通量计算方法

把上述三类源转化为可用于评估的输入，需要经过“能谱模型 → 轨道调制 → 屏蔽传输 → 在轨通量”四个环节。能谱模型给出的是自由空间或特定参考条件下的微分通量，常用单位是 particles/(cm²·s·MeV/n)，其中 MeV/n 表示每核子能量，便于比较不同离子。

GCR常用CREME96或ISO-15390模型。CREME96提供太阳活动极大、极小和平均三种条件下的GCR能谱，并给出从氢到铀的多种离子通量。使用时先选定太阳活动相位，再按任务轨道计算地磁截止刚度Rc。Rc描述地磁场对带电粒子的偏转能力：垂直入射时，刚度低于Rc的粒子被阻挡。对于倾角较低的LEO轨道，Rc在赤道附近最大（约15 GV），在两极附近趋近于0，因此极区GCR通量更高。工程上常用Störmer公式或查表法得到沿轨道的Rc分布，再对GCR能谱做截止处理。

SPE模型常用King模型、August 1972模型或ESP模型。King模型给出不同置信水平（如95%、99%）下的质子积分通量，适合任务周期较长的保守评估；August 1972是历史上最大的SPE之一，常作为最坏情况参考。SPE能谱通常用幂律拟合：Φ(>E) = Φ0 · E^(-γ)，其中γ为谱指数，典型值在2~4之间。使用时需注意SPE是偶发事件，应结合任务时长和太阳活动周期确定事件发生概率，而不是简单取最大值。

俘获带质子常用AP-8或AP-9模型。AP-8按太阳活动极大/极小给出质子通量分布，AP-9则提供更精细的能量和空间分辨率。对于LEO轨道，需要沿轨道积分，得到经过SAA区域时的质子通量增强。一个典型教学示例：假设轨道高度550 km、倾角51.6°，任务周期1年，太阳活动中等。用AP-8模型沿轨道采样，得到能量大于10 MeV的质子平均通量约为 4×10^3 particles/(cm²·s)，而在SAA区域峰值可达 1×10^5 particles/(cm²·s) 量级。这些数值属于示例性假设，实际任务应使用最新模型和实测数据校核。

屏蔽传输环节决定有多少粒子能到达器件敏感区。常用方法是射线追踪法：从器件位置向各个方向发射射线，计算每条射线穿过航天器结构、设备壳体和PCB的等效铝厚度，再对能谱做截止。对于质子，还需考虑核反应产生的次级粒子。工程上可用SHIELDOSE-2或类似工具计算剂量-深度曲线，再结合能谱得到器件位置的粒子通量。

计算要点：在轨通量不是单一数值，而是能量、方向、时间的函数。评估单粒子效应时，通常需要的是LET谱而非质子谱，因为SEU由器件敏感区内的电离能量沉积决定。质子通过核反应产生次级重离子，其LET谱需用蒙特卡罗工具（如Geant4、MCNP）或半经验模型计算。

从通量到单粒子效应率

得到器件位置的在轨粒子通量后，下一步是结合器件截面数据计算单粒子效应率。对于重离子直接电离，常用方法是将LET谱与Weibull或幂律截面对接，积分得到错误率。对于质子，则需先通过核反应模型把质子通量转换为等效LET谱，再与截面对接。典型计算流程为：

输入：轨道参数、屏蔽厚度、太阳活动相位
步骤1：用AP-8/AP-9、CREME96、SPE模型得到自由空间能谱
步骤2：按轨道Rc和SAA穿越时间做地磁截止与时间加权
步骤3：用射线追踪法计算等效铝屏蔽厚度，做能谱截止
步骤4：对质子做核反应次级粒子计算，得到器件位置LET谱
步骤5：与器件截面数据卷积，得到SEU率（错误/天/器件）
输出：在轨SEU率、SEL率、总剂量贡献

需要强调的是，三类源对单粒子效应的贡献不同：GCR重离子直接电离是SEU的主要来源，SPE质子通过次级粒子在事件期间可能造成突发错误率上升，俘获带质子则在LEO轨道长期累积贡献。评估时应分别计算，再按任务剖面加权求和，而不是只取其中一项。对于高可靠任务，还需考虑最坏情况（太阳活动极大、SPE叠加、SAA穿越）下的错误率上限，作为设计裕度的依据。

常见误区：把GCR通量直接当作器件位置通量，忽略地磁截止和屏蔽；或把SPE峰值通量直接用于长期平均，导致过度保守。正确做法是先明确任务轨道、屏蔽配置和太阳活动相位，再选择对应模型和置信水平。

综上，空间辐射环境的量化是单粒子效应评估的输入基础。GCR、SPE和俘获带三类源在能谱、时间行为和空间分布上各有特点，必须分别建模、分别传输、分别加权。只有把在轨通量算准，后续的截面测试、错误率计算和加固设计才有可靠依据。