-
在载人航天活动中,航天员面临空间辐射环境对其身体造成的损伤风险。随着我国载人航天活动从短期驻留向中长期驻留以及从低地球轨道(LEO)向月球基地和深空发展探索,评估航天员的辐射风险并提出有效的防护方案尤为重要[1-2]。参考并深入研究现有的航天员辐射剂量相关数据和结论,可为制定切实可行的航天员辐射防护方案提供重要支撑。
执行LEO任务的航天员所接收的辐射剂量主要源于地球辐射带(ERB)和银河宇宙线(GCR)粒子[3-5]。此外,太阳质子事件和太阳宇宙线等辐射源也会对航天员产生剂量贡献,但相对于总剂量而言其贡献较小[6]。关于ERB和GCR对航天员辐射剂量的贡献,国内外开展了大量研究工作。2001年国际空间站(ISS)搭载的Increment-2人体躯干模型[7]实验测量结果表明,ISS舱内航天员的有效剂量约为0.41 mSv·d-1,其中GCR的有效剂量贡献率约为58.1%,这个数值与Cucinotta等[8]通过HZETRN/QMSFRG空间辐射传输模型计算获得的GCR总剂量贡献率相比仅有3%左右的偏差。2014年El-Jaby等 [9]通过PHITS蒙特卡罗程序计算出在51.6°,450 km轨道处,ISS舱内航天员约75%~85%的总有效剂量当量来自GCR,且GCR中70%的剂量当量贡献来自于Z>2的重离子。由于轨道高度、GCR模型和屏蔽厚度等因素,El-Jaby等的计算结果高于Cucinotta等的计算结果。美国NASA的“火星科学实验室”计算出在10 g·cm-2铝屏蔽下,ISS舱内航天员6~12个月的有效剂量为50~100 mSv,其中1/3来自ERB质子,2/3来自GCR [10]。许峰等[11]利用MRI中国男性体素模型和Geant4蒙特卡罗粒子输运程序并结合我国“神舟五号”飞船的轨道参数,计算得到在1 g·cm-2 等效铝屏蔽下,GCR对航天员有效剂量的贡献率仅为11.1%(只考虑了GCR中质子的剂量贡献)。张斌全等[12]利用Fortran 90辐射剂量计算程序并结合CNMAN人体数字模型,计算出在5 g·cm-2等效铝屏蔽下GCR在航天员单个器官内的吸收剂量率为40 μGy·d-1左右,当量剂量率为200 μSv·d-1,其计算结果比Puchalska[13]等的计算结果偏小20%~50%。
由上可知,国内外对于GCR和ERB粒子对航天员吸收的剂量分析结果存在差异。计算GCR和ERB粒子对航天员吸收剂量的贡献需要综合考虑辐射环境、轨道参数以及人体模型等因素。本文以中国空间站(CSS)和ISS任务为研究对象,利用课题组建立的中国男性体素模型[14],创建Geant4蒙特卡罗程序,对不同轨道下航天员的吸收剂量、当量剂量和有效剂量进行对比分析,旨在评估航天员潜在的辐射风险并提供防护依据。
-
目前从事中长期载人航天活动的CSS和ISS均运行于LEO,空间站航天员平均在轨工作90天,最长工作时间达到438天[15],在轨开展空间科学实验、测试、交会对接等技术试验任务。CSS和ISS的轨道参数如表1所示[16],由表可见空间站处于ERB内带下边缘,受地磁场的保护,可屏蔽大部分宇宙线粒子。
表 1 CSS与ISS运行轨道参数
Table 1. Operational orbital parameters of CSS and ISS
空间站 轨道
倾角/(°)轨道
高度/km远地点/
km近地点/
kmCSS 42.4 393 450 340 ISS 51.6 408 410 400 -
ERB分为内带和外带,空间站处于内带,内带由电子和质子组成。通常将5 g·cm-2 或10 g·cm-2的等效铝屏蔽作为模拟空间站的典型屏蔽厚度[17],根据课题组之前的研究结果[18],经过5 g·cm-2 的等效铝屏蔽后,ERB电子在航天员体内的剂量贡献基本可忽略不计,即ERB剂量贡献主要来源于质子。
采用AP8MIN和AP8MAX模型分别计算太阳活动极小年和极大年ERB质子的微分通量,获得CSS和ISS在屏蔽前/后的微分能谱如图1所示。可以看出,ISS舱内的质子通量在50~60 MeV能量区间内最大,并且太阳极小年的质子通量大约为太阳极大年的2倍。ISS舱内的质子通量计算结果与Dobynde等[4] 的计算结果相比,规律一致。另外, ISS舱外低能质子的通量明显高于CSS的;高能质子的通量两者差别不大,CSS舱外略高。这是由于ISS轨道倾角大,飞行纬度高于CSS,所以其轨道环境中低能质子比较多;ERB质子通量与太阳活动周期同步变化,在太阳极小年呈现高通量,在太阳极大年呈现低通量。
-
采用CREME96模型计算空间站舱外太阳活动极小年和极大年的GCR环境,轨道参数与计算ERB质子模型时所采用的参数相同,同时考虑地磁场的屏蔽作用。图2是计算得到的不同太阳活动期CSS和ISS舱外H、He、C、Mg和Fe等部分GCR重离子的微分能谱。
图 2 不同太阳活动期空间站舱外部分GCR粒子的微分能谱
Figure 2. Differential energy spectra of GCR particles outside the space station during different solar activity periods
从图2可以看出,H和He的微分通量峰值比C、Mg、Fe要高出2~3个数量级,且ISS舱外所有GCR粒子在低能段(<250 MeV)的微分通量比CSS高出数倍。此外,太阳活动极大年时GCR高能段(>20 GeV)粒子通量与极小年相比明显降低。这是由于太阳活动极大年时太阳耀斑活动增强,耀斑释放大量能量所产生的辐射带电粒子会与GCR相互作用,降低高能段GCR粒子的通量。并且太阳活动极大年时太阳风的增强也会导致GCR粒子的散射增强,从而降低了高能段GCR粒子的通量。
-
本文采用课题组搭建的中国成年男性体素模型,创建基于蒙特卡罗算法的Geant4程序,开展仿真计算。体素模型基于一名中国男性全身临床计算机断层扫描(CT)数据建立[14],其特征为:身高176 cm、体重65 kg,全身共659层CT影像,每层的分辨率为512×512像素,单个体素大小为0.97 mm×0.97 mm×3 mm,累积体素为1.72×108 个,模型共包含17种人体组织器官以及各器官的密度和元素占比[14]。与国内早年建立的MRI[19]、CVP[20]和VCH[21]等中国成年男性体素模型相比,该男性体素模型分辨率更高,见表2。
-
吸收剂量、当量剂量和有效剂量是用于衡量航天员吸收辐射的不同剂量指标:吸收剂量为单位质量物质中受到辐射后的能量沉积值;当量剂量为不同组织或器官的吸收剂量乘以相应辐射种类的辐射权重因子,反映了不同类型辐射对人体组织的相对生物效应;有效剂量为不同组织或器官的当量剂量乘以相应的组织或器官权重因子,反映了不同组织或器官对辐射源的基准风险量[22]。
模拟中将航天器简化为一个直径4 m、高度6 m的封闭圆柱体,等效铝屏蔽厚度分别为2、5 和10 g·cm-2 ,航天器内部填充空气,外部为真空,入射源设置为半径5 m的球形面源,粒子从4π立体角各向同性入射到靶物质。人体模型位于航天器中心位置,如图3所示,入射粒子微分能谱分别参见图1和图2。
考虑次级粒子的影响,程序记录初级和次级粒子在航天员每个体素中的能量沉积,每个组织器官的吸收剂量由能量沉积后经过质量加权算法得到,计算当量剂量和有效剂量的辐射权重因子和组织权重因子采用ICRP 103号报告中的WR 和WT 值[23]。
-
重点考虑男性体素模型中眼睛、皮肤、红骨髓、甲状腺等17个人体组织器官,分别将1.1和1.2节中计算获得的ERB质子和GCR各粒子的微分能谱数据以离散点的形式由Geant4中的GPS粒子源类对ERB质子和GCR进行抽样计算,发射源粒子包括ERB质子和原子序数从1~92的GCR粒子,不同辐射粒子在模拟时均单独发射106 个,计算完成后再根据不同辐射粒子的真实通量进行换算得到不同辐射源在人体组织或器官中的辐射剂量。
-
在太阳活动极小年和极大年条件下,经5 g·cm-2等效铝屏蔽后,ERB质子、GCR质子、α粒子和其他重离子在航天员组织器官内的吸收剂量分别如图4 (a)~(d)所示。
图 4 不同太阳活动时期航天员器官或组织吸收剂量
Figure 4. Absorbed doses by astronaut organs or tissues during different solar activity periods
从图4中可以看出:太阳活动极小年时ERB质子和GCR粒子对航天员辐射剂量大;太阳活动极大年时ERB质子和GCR粒子对航天员辐射剂量小。这是由于太阳活动极小年时的质子通量比极大年高约1个数量级所致。从图4 (a)可以看出,在ERB质子辐照条件下,ISS舱内的大部分航天员组织器官所吸收的辐射剂量明显高于CSS的,只有个别器官如软组织和膀胱是在CSS内辐射剂量较高。这是由于从空间站所接收的ERB粒子能谱来看,ISS在低能区的粒子通量高于CSS的,而在高能区CSS比ISS的高。而高能粒子能量沉积在深层器官,因此个别深层器官吸收的辐射剂量在CSS舱内要高于ISS舱内。此外,在ERB质子辐照条件下,航天员表层器官所吸收的辐射剂量高,深层器官所吸收的辐射剂量低,即皮肤所吸收的辐射剂量最高,其次是眼睛、肝、结肠等。在CSS舱内,ERB质子对皮肤的辐射剂量率约为0.18 mGy·d-1,而ISS舱内约为0.45 mGy·d-1,后者约为前者的2.5倍。从图4(b)~(d)可以看出,在GCR辐照条件下,ISS舱内的航天员组织器官所吸收的辐射剂量明显高于CSS的,与ERB质子相比,在航天员最表层器官皮肤内产生的辐射剂量反而没有眼睛、肝、结肠等接近表层的器官剂量高。这与GCR的能谱分布有关:能量低的GCR粒子通量并不是最高,通量最高的GCR粒子能量普遍在几百到上千MeV/n,这些高能粒子一般可以穿透皮肤,沉积到航天员较深层器官中。
表3为不同太阳活动期ERB质子和GCR粒子在航天员体内的吸收剂量分布。可以看出,ISS舱内的航天员吸收剂量高于CSS的,不同粒子在航天员体内的剂量贡献排序一致:ERB质子>GCR质子>GCR α粒子>GCR其余重离子。此外,航天员组织器官约80%的吸收剂量来自于ERB质子,而GCR粒子在航天员体内辐射剂量中占比最高的是质子,约占14%;其次为α粒子,约占5%;其他重离子占1%左右。
表 3 不同太阳活动时期ERB和GCR粒子在航天员体内的吸收剂量
Table 3. Absorbed doses of ERB and GCR particles by an astronaut during different solar activity periods
单位:mGy·d-1 粒子 CSS ISS 太阳活动
极大年太阳活动
极小年太阳活动
极大年太阳活动
极小年ERB 质子 1.260 0.570 2.390 0.910 GCR 质子 0.210 0.110 0.370 0.160 α粒子 0.075 0.041 0.160 0.048 其他重离子 0.012 0.007 0.047 0.010 合计 1.557 0.728 2.967 1.128 基于已发表的ISS人体模型实验剂量探测数据,本文根据太阳活动极小年计算的ISS航天员部分器官的吸收剂量率与MATROSHKA人体模型实验[24] (太阳活动极小年附近)获得的器官剂量探测结果的平均值进行比较,如图5所示。
图 5 本文计算的ISS航天员部分器官的吸收剂量与模型实验测量值的比较
Figure 5. Comparison of the calculated absorbed doses in this article and the experimental measurement for selected organs of an ISS astronaut model
可看出,计算结果与人体模型实验测量结果相比一致性良好,相对误差均在30%之内,尤其是肾的吸收剂量与实验结果比较吻合,其他组织器官有所差别。这主要是由于本文模拟所使用的辐射环境模型是AP8和CREME96,是基于20世纪的卫星观测数据推导得到的经验公式模型,与空间站外真实辐射环境相比有所偏差[6]。同时,不同人体模型具有不同的器官位置、体积和形状,这种差异会影响辐射粒子与器官组织相互作用后能量的沉积和散射,从而造成辐射环境对特定个体器官的辐射剂量差异[2]。
-
ISS和CSS舱内航天员组织器官的当量剂量如图6所示。在当量剂量的计算中引入了辐射权重因子WR,不同粒子的辐射权重因子如表4所示,从表中可以看出GCR α粒子和重离子在进入航天员组织器官后产生的当量剂量是吸收剂量的20倍,而ERB质子进入航天员组织器官后的当量剂量是吸收剂量的2倍,两者相差10倍。
图 6 不同太阳活动期航天员器官或组织的当量剂量
Figure 6. Equivalent doses for organs or tissues of an astronaut during different solar activity periods
表 4 不同辐射粒子源的辐射权重因子
Table 4. Radiation weight factors for different particle radiation sources
辐射类型 辐射权重因子 >2 MeV质子(不包括反冲质子) 2 α粒子,裂变碎片,重核 20 从图6可以看出,α粒子在航天员各个器官的当量剂量比GCR质子高了3~4倍,其他重离子在航天员各个器官的当量剂量几乎与GCR质子持平。
表5是不同太阳活动期ERB质子和GCR粒子对航天员的当量剂量贡献分布,在太阳活动极小年期间,CSS舱内航天员组织器官大约55%的当量剂量来源于ERB质子;ISS舱内航天员组织器官大约49%的当量剂量来源于ERB质子。太阳活动极大年期间,CSS舱内航天员组织器官大约53%的当量剂量来源于ERB质子;ISS舱内航天员组织器官大约55%的当量剂量来源于ERB质子。综上,相较在航天员吸收剂量中的占比,GCR当量剂量贡献均在50%左右,占比显著上升。
表 5 不同太阳活动期ERB和GCR粒子在航天员体内的当量剂量贡献分布
Table 5. Distributions of equivalent dose contribution of ERB and GCR particles in an astronaut during different solar activity periods
单位:mSv·d-1 粒子 CSS ISS 太阳活动
极大年太阳活动
极小年太阳活动
极大年太阳活动
极小年ERB 质子 2.52 1.27 4.78 1.82 GCR 质子 0.42 0.22 0.74 0.32 α粒子 1.50 0.78 3.20 0.96 其他重离子 0.24 0.14 0.20 9.66 合计 4.68 2.41 9.66 3.30 -
从辐射防护角度来看,空间辐射风险评估的基准量是有效剂量,通过将当量剂量HT与组织权重因子WT相乘,可得到ERB质子和GCR粒子对航天员总有效剂量贡献的分布,结果如图7所示。通过模拟计算获得(在5 g·cm-2 等效铝屏蔽下):太阳活动极小年时,CSS舱内航天员的总有效剂量为385 μSv·d-1,GCR对航天员的总有效剂量贡献占比约为46%;ISS舱内航天员的总有效剂量为787 μSv·d-1,GCR对舱内航天员吸收的总有效剂量贡献占比较CSS略大,约为51%。太阳活动极大年时,CSS舱内航天员吸收的总有效剂量为192 μSv·d-1,GCR对航天员的总有效剂量贡献占比约为51%;ISS舱内航天员总有效剂量为266 μSv·d-1,GCR对航天员的总有效剂量贡献占比约为45%。
图 7 不同太阳活动期ERB和GCR粒子在航天员体内的日均有效剂量
Figure 7. Daily effective doses of ERB and GCR particles in an astronaut during different solar activity periods
本文除计算5 g·cm-2 典型铝屏蔽厚度下的航天员剂量外,还计算了2 g·cm-2 和10 g·cm-2 铝屏蔽厚度下航天员的有效剂量,如图8所示。通过不同屏蔽厚度下的航天员年有效剂量与NASA LEO男性航天员的年有效剂量限值[25]对比来看,当铝屏蔽厚度为2 g·cm-2 时,ISS男性航天员在太阳活动极小年吸收的年有效剂量超过辐射风险限值;当铝屏蔽厚度≥5 g·cm-2 时,CSS和ISS男性航天员在太阳活动极小年和极大年的年有效剂量均远低于NASA男性航天员辐射风险限值,表明将5 g·cm-2的铝屏蔽应用于LEO空间站辐射环境中,可兼顾减少材料成本与良好的防护作用。
-
通过计算发现:对于不同太阳活动期,航天员吸收的约80%和约20%的辐射剂量分别来自于ERB质子和GCR粒子的贡献,并且在ERB质子辐照下,航天员外层器官的辐射剂量高于深层器官,皮肤的最高;而在GCR粒子辐照下,皮肤表层的辐射剂量没有接近表层器官的剂量高;本文的计算结果与MATROSHKA人体模型实验测量结果相比,器官吸收剂量的相对误差均在30%之内。此外,航天员吸收的当量剂量和有效剂量50%左右来自于ERB质子,其他来自于GCR粒子的贡献。通过本文计算的不同屏蔽厚度下的航天员年有效剂量与NASA LEO男性航天员的年有效剂量限值对比来看,CSS和ISS男性航天员在不同太阳活动期的年有效剂量在5 g·cm-2的铝屏蔽下均低于辐射风险限值,表明5 g·cm-2的铝屏蔽能够在空间站辐射环境中起到良好的防护作用。
本研究仍存在局限性和不足:所借助的AP8和CREME96模型实现的ERB和GCR粒子的能谱输入是静态能谱模型,所用数据是20世纪的卫星观测数据,LEO空间中的部分区域没有被覆盖,与真实辐射环境间存在一定的偏差;另外本文简化了空间站核心舱的结构,而舱体结构的差异可能导致带电粒子在屏蔽材料中的作用过程不同,这会给航天员辐射剂量的评估带来不确定性。因此,下一步拟使用更接近空间站舱外真实情况的粒子能谱和空间站核心舱实际使用的材料和结构进行模拟,以更有效地评估航天员的真实辐射暴露风险。
Analysis of the dose effect of space radiation to astronauts on the space station
-
摘要: 为有效实施对航天员的空间辐射防护,采用计算机断层扫描数据(CT)建立了精细化男性体素模型;基于Geant4建立蒙特卡罗程序计算空间站轨道辐射经过舱壁(5 g·cm-2等效铝屏蔽)后在体模中的辐射剂量;分析了不同辐射粒子在体素模型组织或器官中的吸收剂量、当量剂量和有效剂量。计算结果表明:航天员体内吸收剂量大约80%来源于地球辐射带(ERB)质子;大约14%来源于银河宇宙线(GCR)质子;α粒子的剂量贡献占比约为5%;其余重离子的剂量贡献占比在1%左右。另外,航天员吸收的当量剂量和有效剂量50%左右来自于ERB质子,另50%左右来自于GCR粒子。计算结果将有助于评估航天员在空间站舱内的潜在辐射风险并提供辐射防护参考。Abstract: In order to effectively protect astronauts from space radiation, a finely detailed male voxel phantom was created using computed tomography (CT) data. A Monte Carlo program based on Geant4 was established to calculate orbital radiation doses in the phantom after the radiation passing through the module wall of the space station (equivalent to 5 g·cm-2 of aluminum shielding). The absorbed dose, equivalent dose, and effective dose of various radiation particles in the voxel phantom's tissues or organs were analyzed. The results show that approximately 80% of the absorbed dose in astronauts is from Earth radiation belt (ERB) protons; approximately 14% is from galactic cosmic ray (GCR) protons; and alpha particles accounted for about 5% of the dose contribution; the contribution of other heavy ions to the absorbed dose was around 1%. In addition, approximately 50% of equivalent dose and effective dose received by astronauts came from ERB protons, while the other 50% came from GCR. The calculation results will aid in evaluating the potential radiation risks received by astronauts aboard the space station and provide a reference for radiation protection.
-
Key words:
- space station /
- astronauts /
- radiation protection /
- Earth radiation belt /
- galactic cosmic ray /
- radiation dose /
- Geant4 software
-
表 1 CSS与ISS运行轨道参数
Table 1. Operational orbital parameters of CSS and ISS
空间站 轨道
倾角/(°)轨道
高度/km远地点/
km近地点/
kmCSS 42.4 393 450 340 ISS 51.6 408 410 400 表 2 中国成年男性不同体素模型分辨率对比
Table 2. Comparison of resolution of different voxel phantoms for a Chinese adult male
表 3 不同太阳活动时期ERB和GCR粒子在航天员体内的吸收剂量
Table 3. Absorbed doses of ERB and GCR particles by an astronaut during different solar activity periods
单位:mGy·d-1 粒子 CSS ISS 太阳活动
极大年太阳活动
极小年太阳活动
极大年太阳活动
极小年ERB 质子 1.260 0.570 2.390 0.910 GCR 质子 0.210 0.110 0.370 0.160 α粒子 0.075 0.041 0.160 0.048 其他重离子 0.012 0.007 0.047 0.010 合计 1.557 0.728 2.967 1.128 表 4 不同辐射粒子源的辐射权重因子
Table 4. Radiation weight factors for different particle radiation sources
辐射类型 辐射权重因子 >2 MeV质子(不包括反冲质子) 2 α粒子,裂变碎片,重核 20 表 5 不同太阳活动期ERB和GCR粒子在航天员体内的当量剂量贡献分布
Table 5. Distributions of equivalent dose contribution of ERB and GCR particles in an astronaut during different solar activity periods
单位:mSv·d-1 粒子 CSS ISS 太阳活动
极大年太阳活动
极小年太阳活动
极大年太阳活动
极小年ERB 质子 2.52 1.27 4.78 1.82 GCR 质子 0.42 0.22 0.74 0.32 α粒子 1.50 0.78 3.20 0.96 其他重离子 0.24 0.14 0.20 9.66 合计 4.68 2.41 9.66 3.30 -
参考文献
[1] CHANCELLOR J C, BLUE R S, CENGEL K A, et al. Limitations in predicting the space radiation health risk for exploration astronauts[J]. Npj Microgravity, 2018, 4(1): 1-11 doi: 10.1038/s41526-017-0036-6 [2] NORBURY J W, SLABA T C, AGHARA S, et al. Advances in space radiation physics and transport at NASA[J]. Life Sciences in Space Research, 2019, 22: 98-124 doi: 10.1016/j.lssr.2019.07.003 [3] CUCINOTTA F A, WILSON J W, WILLIAMS J R, et al. Analysis of Mir-18 results for physical and biological dosimetry: radiation shielding effectiveness in LEO[J]. Radiation Measurements, 2000, 32(3): 181-191 doi: 10.1016/S1350-4487(99)00273-5 [4] DOBYNDE M I, SHPRITS Y Y. Radiation environment created with GCRs inside a spacecraft[J]. Life Sciences in Space Research, 2020, 24: 116-121 doi: 10.1016/j.lssr.2019.09.001 [5] 赵磊, 尚钰轩, 袁爽, 等. 载人深空探索中空间辐射防护技术的研究进展[J]. 科学通报, 2019, 64(20): 2087-2103 ZHAO L, SHANG Y X, YUAN S, et al. Research progress in space radiation protection technology in manned deep space exploration[J]. Science Bulletin, 2019, 64(20): 2087-2103 [6] CUCINOTTA F A, SCHIMMERLING W, WILSON J W, et al. Uncertainties in estimates of the risks of late effects from space radiation[J]. Advances in Space Research, 2004, 34(6): 1383-1389 doi: 10.1016/j.asr.2003.10.053 [7] YASUDA H. Effective dose measured with a life size human phantom in a low Earth orbit mission[J]. Journal of Radiation Research, 2009, 50(2): 89-96 doi: 10.1269/jrr.08105 [8] CUCINOTTA F A, KIM M H Y, WILLINGHAM V, et al. Physical and biological organ dosimetry analysis for International Space Station astronauts[J]. Radiation Research, 2008, 170(1): 127-138 doi: 10.1667/RR1330.1 [9] EL-JABY S, TOMI L, SIHVER L, et al. Method for the prediction of the effective dose equivalent to the crew of the International Space Station[J]. Advances in Space Research, 2014, 53(5): 810-817 doi: 10.1016/j.asr.2013.12.022 [10] NELSON G A. Space radiation and human exposures, a primer[J]. Radiation Research, 2016, 185(4): 349-358 doi: 10.1667/RR14311.1 [11] XU F, ZENG Z, JIA X H. Estimation of doses to Taikonauts in SZ-5/6/7 missions using a male voxel model[C]//Tenth International Conference on Photonics and Imaging in Biology and Medicine (PIBM 2011). SPIE, 2012, 8329: 17-23 [12] 张斌全, 余庆龙, 梁金宝, 等. 航天员受银河宇宙线辐射的剂量计算[J]. 北京航空航天大学学报, 2015, 41(11): 2044-2051 ZHANG B Q, YU Q L, LIANG J B, et al. Dose calculation of astronauts exposed to galactic cosmic ray radiation[J]. Journal of Beijing University of Aeronautics and Astronautics, 2015, 41(11): 2044-2051 [13] PUCHALSKA M, SIHVER L, SATO T, et al. Simulations of MATROSHKA experiment outside the ISS using PHITS[J]. Advances in Space Research, 2012, 50(4): 489-495 doi: 10.1016/j.asr.2012.04.027 [14] 张昭, 方美华, 杨航. 不同等效屏蔽厚度下体模器官空间辐射剂量计算的研究[J]. 载人航天, 2022, 28(3): 285-290 doi: 10.3969/j.issn.1674-5825.2022.03.001 ZAHNG Z, FANG M H, YANG H. Study on the calculation of radiation dose in phantom organs with different equivalent shielding thickness[J]. Manned Spaceflight, 2022, 28(3): 285-290 doi: 10.3969/j.issn.1674-5825.2022.03.001 [15] 刘建忠, 王勇, 姚小丽, 等. 空间辐射剂量测量简介[J]. 中国辐射卫生, 2010, 19(4): 458-461 doi: 10.13491/j.cnki.issn.1004-714x.2010.04.017 LIU J Z, WANG Y, YAO X L, et al. Introduction to space radiation dose measurement[J]. China Radiation Health, 2010, 19(4): 458-461 doi: 10.13491/j.cnki.issn.1004-714x.2010.04.017 [16] 翟睿琼, 姜海富, 田东波, 等. 空间站原子氧环境仿真研究[J]. 装备环境工程, 2014, 11(3): 35-39 ZHAI R Q, JIANG H F, TIAN D B, et al. Simulation study on atomic oxygen environment of space station[J]. Equipment Environmental Engineering, 2014, 11(3): 35-39 [17] 沈自才, 夏彦, 杨艳斌, 等. 航天器空间辐射防护材料与防护结构[J]. 宇航材料工艺, 2020, 50(2): 1-7 SHEN Z C, XIA Y, YANG Y B, et al. Space radiation protection materials and structures for spacecraft[J]. Aerospace Material Technology, 2020, 50(2): 1-7 [18] GUO Y P, YAN F Z, FANG M H, et al. Monte Carlo calculation of Earth’s trapped radiation exposure to Chinese female astronauts onboard China Space Station[J]. Chinese Physics B, 2022, 13(1): 25-31 [19] JIA X H, XU F, HUANG Z X, et al. Estimates of space radiation exposure to astronauts using male voxel model based on MRI[J]. Space Medicine & Medical Engineering, 2008: 299-303 [20] LI J L, QIU R, ZHANG Z, et al. Organ dose conversion coefficients for external photon irradiation using the Chinese voxel phantom (CVP)[J]. Radiation Protection Dosimetry, 2009, 135(1): 33-42 doi: 10.1093/rpd/ncp087 [21] ZHANG G Z, LUO Q M, ZENG S Q, et al. The development and application of the visible Chinese human model for Monte Carlo dose calculations[J]. Health Physics, 2008, 94(2): 118-125 doi: 10.1097/01.HP.0000285256.48498.b4 [22] 曾志, 李君利, 贾向红, 等. 空间辐射剂量及屏蔽效应研究[J]. 清华大学学报: 自然科学版, 2008, 48(3): 391-394 ZENG Z, LI J L, JIA X H, et al. Study on space radiation dose and shielding effect[J]. Journal of Tsinghua University (Natural Science Edition), 2008, 48(3): 391-394 [23] 邬仁耀, 耿长冉, 田锋, 等. 基于蒙特卡罗方法的航空机组人员辐射剂量评估[J]. 辐射防护, 2022, 42(5): 442-449 doi: 10.3969/j.issn.1000-8187.2022.5.fsfh202205008 WU R Y, GENG C R, TIAN F, et al. Radiation dose assessment for aviation crew based on Monte Carlo method[J]. Radiation Protection, 2022, 42(5): 442-449 doi: 10.3969/j.issn.1000-8187.2022.5.fsfh202205008 [24] PUCHALSKA M, BILSKI P, BERGER T, et al. NUNDO: a numerical model of a human torso phantom and its application to effective dose equivalent calculations for astronauts at the ISS[J]. Radiation and Environmental Biophysics, 2014, 53(4): 719-727 doi: 10.1007/s00411-014-0560-7 [25] CUCINOTTA F A, HU S, SCHWADRON N A, et al. Space radiation risk limits and Earth-Moon-Mars environmental models[J]. Space Weather, 2010, 8(12): 77-85 -