1. 引言
航天员在深空探索中面临着空间辐射的威胁,在其执行任务前应预先针对复杂的空间辐射环境进行风险评估,才能针对性地提出防护措施,降低航天员飞行过程中所受的空间辐射风险。按照来源不同,空间辐射可分为银河宇宙射线(GCR)和太阳质子事件(SPE),以及在近地轨道中的地磁捕获辐射(TPR)。GCR辐射粒子能量很高,但其在空间环境中相对稳定且注量率较低;SPE是随机事件,虽然一年仅出现几次,但SPE注量率极高,若无有效措施防护,航天员健康将面临严重威胁。无论是近地还是深空飞行,质子均是航天员受到的主要辐射粒子,在空间中占比80%以上,其中,GCR质子的能量相对分散,能量以10~1000 MeV为主,而SPE质子能量主要集中在10~100 MeV [1]。
空间质子能量分布宽,不同能量的质子在人体组织器官中会产生不同的吸收剂量。此外,不同组织器官处于人体不同深度位置且辐射敏感性不同,因此空间辐射对各组织器官的危害也不同。获得不同能量的质子在组织器官中的吸收剂量,研究不同能量的质子对人体组织的作用,可为辐射防护策略的制定提供输入,对于飞行前辐射风险评估十分必要。
由于直接获取人体器官剂量数据困难,航天国家使用国际辐射单位与测量委员会(ICRU)设定的组织等效球/板 [2]、人体模型 [3] [4] 等方法进行器官辐射效应研究。人体模型质量大、探测器取放复杂,在航天应用中有一定局限,而ICRU球/板与探测器整合构建的监测体系 [5] [6] 与人体模型相比载荷小、操作简便,是空间监测设备研发和校准的基础工具,从任务实际出发更适于空间辐射监测评估。然而实际制作ICRU球/板工艺、辐射源条件等因素会略有差异,不同监测体系的实测结果可能存在偏差,因此,基于ICRU球/板进行辐射剂量仿真计算,获得特定照射条件下的吸收剂量,是建立ICRU球/板与探测器构成的监测体系的前提,基于ICRU球/板进行辐射剂量仿真计算可为剂量监测有效性验证及校准、深部组织剂量预估提供参考和方法体系。本工作基于Geant4程序建立了ICRU模型,通过仿真计算分析了空间主要能量质子及地面常用辐射源γ射线与ICRU球和平板相互作用,模拟了不同深度等效人体组织内探测到的吸收剂量,为ICRU球/板与探测器整合构建的监测体系的有效性验证和空间辐射风险评估提供了参考依据。
2. 研究方法
2.1. ICRU模型
国际辐射单位与测量委员会(ICRU)设定的球模型和平板模型如图1所示。ICRU球模型直径为30 cm,平板模型尺寸为30 cm × 30 cm × 15 cm,二者材料均由71.39%的O,14.89%的C,10.00%的H,3.47%的N,0.15%的Na和0.10%的Cl组成,密度均为1 g/cm3。基于ICRU模型,可计算个人剂量dp(H)用于辐射监测。通常取深度为10 mm的剂量,作为躯干所受有效剂量的近似值;取深度为3 mm的剂量,作为晶状体所受当量剂量的近似值;取深度为0.07 mm的剂量,作为皮肤所受当量剂量的近似值。
![](//html.hanspub.org/file/9-2980145x10_hanspub.png)
Figure 1. Sphere (a) and slab (b) model recommended by ICRU
图1. ICRU设定的球模型(a)和平板模型(b)
2.2. 辐射源
1) 质子
入射质子与靶材料相互作用包括弹性和非弹性物理过程。与阻止介质的核外电子的非弹性碰撞,使核外电子激发或者电离,是质子能量损失的最主要形式。此外,与原子核发生弹性碰撞,原子核因获得动能而发生反冲产生原子位移缺陷,也是质子的能量损失方式之一,而质子与原子核的非弹性碰撞、与核外电子的弹性碰撞的能量损失较小可忽略不计 [7]。
质子穿入物质与材料中的电子和原子相互碰撞而损失能量,通常用单位距离上的平均能量损失(−dE/dx)对该过程进行描述。质子的能量损失主要为电离能量损失,遵循Bethe-Block公式 [8]:
Bethe-Block公式中m0是电子静止质量,ν是入射粒子的速度,z是带电粒子电荷,Z是阻止介质原子序数,
,c是光速,I是介质原子的平均电离电势。
质子在某物质中沿入射方向穿过的最大距离为质子的射程。图2是利用SRIM软件,计算的不同能量质子在人体组织等效材料内的射程。质子与物质相互作用时在射程末端形成布拉格(Bragg)峰,随着质子穿入介质的深度增加,不断产生电离,质子能量减少,运动速度降低;而慢速的质子又会引起更多的电离即Bragg峰,Bragg峰后粒子能量基本全部释放,电离能力突变降为0。
2) γ射线
γ射线本身不带电,因而用相互作用截面、能量转移和能量吸收的几率来描述其在物质中的行为。γ射线与物质相互作用最主要的方式为光电效应、康普顿效应、电子对效应等,γ射线能量主要转换为次级电子能量,通过计算次级电子能量沉积可得到γ射线能量沉积。
![](//html.hanspub.org/file/9-2980145x13_hanspub.png)
Figure 2. The range of protons in the tissue-equivalent materials
图2. 质子在人体组织等效材料内的射程
2.3. 方法程序
蒙特卡罗方法是使用随机数或伪随机数来解决大量随机样本中某一事件发生概率的方法 [9]。在辐射评估方面,参与核反应的辐射粒子在物质中输运的过程具有一定的统计规律性,适用蒙特卡罗方法分析。运用蒙特卡罗方法,通过在计算机上随机的模拟粒子的运动过程,并追踪这些粒子的运动信息,从而得出这些粒子运动的统计规律,分析粒子在物质中的输运问题。该方法的核心是随机数产生器目前蒙特卡洛方法中随机数并行化的技术已经较成熟,常见的蒙特卡洛程序包有MCNP、FLUKA、GEANT等。
![](//html.hanspub.org/file/9-2980145x14_hanspub.png)
Figure 3. Main module structure diagram of GEANT4
图3. GEANT4的主要模块结构图
GEANT4是欧洲核子中心(CERN)开发的一种蒙特卡洛程序,在核物理实验模拟、空间辐射、加速器模拟、放射医学等领域均有广泛应用。该程序程序代码开放,采用面向对象的C++语言编写,具备完整的粒子与物质相互作用过程的材料库,包括几何跟踪,探测器响应,可视化和用户界面等模块的软件包,是粒子物理与核物理实验模拟最好的工具之一。本工作选用GEANT4程序基于ICRU模型对吸收剂量进行模拟计算。GEANT4主要模块图如图3所示。其中Global模块包含了程序所需要的常数、单位、数值运算以及随机数产生等,Particle和Material模块描述粒子和材料的物理特性,Tracking和Track模块与粒子的运输过程有关,Process模块描述粒子与物质相互作用,可设置粒子的能量、动量以及种类、对参与过程的物质等,Run是对模拟事件的控制指令,Visualization模块可以将探测器、粒子径迹等转化为图形。在模拟时,首先要明确粒子会与物质会发生哪些重要的物理过程。在各种粒子与物质的相互作用的过程中,Geant4中主要有7种物理过程包括电与磁相互作用、强子之间相互作用过程、原子核的衰变过程、光轻子与强子之间的相互作用过程、光学的光子的过程、设置参数过程以及粒子输运过程。
2.4. 输入条件
本计算使用Geant4.9.6版本,根据ICRU平板模型和球模型的几何构造和元素组成建立了物理模型,采用质子和γ射线两种辐射源进行仿真计算,选用的物理模型有强相互作用G4QGSM模型,光电效应、康普顿效应、电子效应、EM输运等物理模型,粒子的截止阈值为1 mm [10]。
1) 基于空间质子能谱及地面质子加速器试验能段范围,质子能量选取10 MeV、20 MeV、40 MeV、50 MeV、60 MeV、80 MeV、100 MeV、120 MeV、150 MeV及200 MeV;基于空间舱内屏蔽后的γ能谱分布、地面钴60的γ平均能量1250 keV,γ射线能量选取1000 keV、1173 keV、1320 keV、1500 keV、2000 keV、2500 keV、3000 keV及3500 keV;
2) 质子、γ射线均采用入射角为cos分布、度数0~90度之间照射;
3) 质子注量选取SPE典型恶劣环境注量107 cm−2且为归一化常数,γ射线注量与质子一致;
4) ICRU模型深度选取0.07 mm、3 mm、10 mm、30 mm及100 mm,探测在该深度下的吸收剂量。
3. 结果与讨论
3.1. ICRU球模型中辐射剂量分布
表1列出了球模型中注量107 cm−2的归一化吸收剂量与质子能量、入射深度的对应关系。通过表1中深度为0.07、3、10 mm归一化吸收剂量与实际粒子注量的乘积,可近似得到特定能量质子照射下,皮肤、晶状体所受当量剂量以及躯干所受有效剂量,为辐射评估提供参考依据。
![](Images/Table_Tmp.jpg)
Table 1. The normalized absorbed dose (Gy·cm2) in sphere model of different detecting depths caused by protons of different energies
表1. 不同探测深度下不同能量质子造成的球模型中归一化吸收剂量(Gy·cm2)
图4为球模型中归一化吸收剂量随质子能量(a)和材料厚度(b)的变化。由图4(a)可见,随着质子能量增加,吸收剂量先增大后然后逐渐下降,在某个能量下达到峰值,且随着探测深度增加,峰值对应的质子入射能量逐渐增大。这与图2中质子在组织等效材料的射程有关,质子与等效人体组织相互作用不断产生电离而损失能量,能量沉积在靶物质中。当探测深度固定,质子能量小射程小时,探测深度大于质子射程,因此该能量的质子对吸收剂量贡献小;随着质子能量增加,探测深度与质子射程接近,吸收剂量达到峰值,而后随着质子能量继续增加,射程逐渐超过探测厚度,吸收剂量逐渐减小,且随着探测厚度增加,达到吸收剂量峰值需要的质子入射能量增大。结合表1数据,由图4(b)可见,质子能量在0~50 MeV时,随着探测深度增加吸收剂量减小,若等效为人体剂量,则人体内部器官相对皮肤受到辐射的剂量较小。当质子能量大于60 MeV,30 mm探测深度的吸收剂量超过0.07,3,10 mm探测深度的值,即深部器官相对皮肤受到的辐射剂量更大。
![](//html.hanspub.org/file/9-2980145x15_hanspub.png)
Figure 4. The proton energy (a) and material thickness (b) dependence of the normalized absorbed dose (Gy·cm2) in sphere model
图4. 球模型中归一化吸收剂量随质子能量(a)和材料厚度(b)的变化
表2列出了不同探测深度下不同能量γ射线造成球模型中的注量107 cm−2的归一化吸收剂量。γ射线为地面常用辐射源。通过表2中归一化吸收剂量与实际粒子注量的乘积得到特定注量下的吸收剂量,可作为ICRU球的地面试验验证的参考。
![](Images/Table_Tmp.jpg)
Table 2. The normalized absorbed dose (Gy·cm2) in sphere model of different detecting depths caused by γ ray of different energies
表2. 不同探测深度下不同能量γ射线造成球模型中的归一化吸收剂量(Gy·cm2)
图5为球模型中注量107 cm−2的归一化吸收剂量随γ射线能量(a)和探测厚度(b)的变化,与图4对比可见,γ射线与质子的作用效果不同。图5(a)中,当探测深度小于3 mm且质子能量较大时,吸收剂量随能量上升而下降,这是由于γ射线在物质中产生的次级电子还未达到平衡;当深度大于3 mm时,次级电子逐渐平衡,因此随入射γ射线能量增大,组织等效材料的吸收剂量也变大。图5(b)中曲线前半段吸收剂量随厚度上升,这是正面碰撞产生的次级电子的积累作用造成,后半段随着深度增加,由于电离损耗逐渐减小,吸收剂量下降。
![](//html.hanspub.org/file/9-2980145x16_hanspub.png)
Figure 5. The γ-ray energy (a) and material thickness (b) dependence of the normalized absorbed dose (Gy·cm2) in sphere model
图5. 球模型中归一化吸收剂量随γ射线能量(a)和探测厚度(b)的变化
3.2. ICRU平板模型中辐射剂量分布
平板模型地面模拟照射实验可行性高,是作为地面材料试验验证的有效参考。因此表3中列出了不同探测深度下不同能量质子造成的平板模型中注量107 cm−2的归一化吸收剂量,可作为ICRU平板实物地面照射的对照,为不同深度辐射剂量评估提供参考依据。与表1数据对比可见,相同辐照条件下,球模型吸收剂量大于平板模型,这是因为球模型是4π照射,次级粒子散射到同一探测深度的概率较大。通过表3中的数据作图,得到图6中平板模型中注量107 cm−2的归一化吸收剂量随质子能量(a)和材料厚度(b)的变化,通过对比可见,平板模型中吸收剂量的变化规律与图4球模型中的规律一致,均先变大,达到一定阈值后缓慢衰减。
![](Images/Table_Tmp.jpg)
Table 3. The normalized absorbed dose (Gy·cm2) in slab model of different detecting depths caused by protons of different energies
表3. 不同探测深度下不同能量质子造成的平板模型中的归一化吸收剂量(Gy·cm2)
![](//html.hanspub.org/file/9-2980145x17_hanspub.png)
Figure 6. The proton energy (a) and material thickness (b) dependence of the normalized absorbed dose (Gy·cm2) in slab model
图6. 平板模型中归一化吸收剂量随质子能量(a)和材料厚度(b)的变化
表4列出了平板模型中注量107 cm−2的归一化吸收剂量与γ射线能量、入射深度的对应关系。γ射线为地面常用辐射源,通过计算归一化吸收剂量与实际粒子注量的乘积得到特定注量下的吸收剂量,可作为ICRU平板实物地面试验验证的参考。与表2相比,相同辐照条件下,平板模型相同深度的吸收剂量与球模型基本一致,差异在5%以内,因此,平板模型的吸收剂量随γ射线能量和探测深度的变化规律与图4球模型中的一致。
![](Images/Table_Tmp.jpg)
Table 4. The normalized absorbed dose (Gy·cm2) in slab model of different detecting depths caused by γ-rays of different energies
表4. 不同探测深度下不同能量γ射线造成的平板模型归一化吸收剂量(Gy·cm2)
4. 结论
本工作建立了ICRU球模型和平板模型,基于Geant4软件获得了ICRU模型不同深度的吸收剂量与质子、γ射线能量的对应关系,建立了空间辐射剂量监测校准和剂量预估仿真计算体系。由仿真结果可见,平板模型和球模型中组织等效材料的吸收剂量随着质子能量和探测深度变化的规律一致,随着单能质子能量的增加,吸收剂量均先变大,达到一定阈值后缓慢衰减;由于4π方向入射,相同条件下球模型的吸收剂量比平板模型略高。对地面常用参考辐射源γ射线进行了仿真计算,结果表明不同能量γ射线、不同探测深度下,γ射线在平板模型和球模型中的作用规律和吸收的剂量基本一致。
致谢
本文是在航天员科研训练中心各位领导和专家支持下完成,向各位表示衷心的感谢!
基金项目
本项目由航天员科研训练中心计量课题、实验技术研究青年基金项目2017SY54C0401支持。
NOTES
*通讯作者。