吳正新，孫慧斌,，何承發(fā)，童永彭，馬玉剛，陸景彬，荀明珠，呼延奇，蔡震波，朱志鵬（. 吉林大學 物理學院，長春 00；. 深圳大學 核技術研究所，深圳 58060；. 中國科學院 新疆理化技術研究所，烏魯木齊 800；. 北京空間飛行器總體設計部，北京 0009）

?

載人航天器艙內輻射環(huán)境及劑量分析

吳正新1，孫慧斌1,2，何承發(fā)3，童永彭2，馬玉剛1，
陸景彬1，荀明珠3，呼延奇4，蔡震波4，朱志鵬2
（1. 吉林大學 物理學院，長春 130012；2. 深圳大學 核技術研究所，深圳 518060；
3. 中國科學院 新疆理化技術研究所，烏魯木齊 830011；4. 北京空間飛行器總體設計部，北京 100094）

摘要：構建載人航天器模型，利用Geant4軟件計算了銀河宇宙射線中的質子穿過載人航天器后在水模體中的總吸收劑量以及次級粒子吸收劑量，并且統(tǒng)計了不同種類的次級粒子在水模體中的總數量分布。計算結果表明，次級粒子引起的輻射劑量所占比例為50%，且二次電子總數量達到了107量級。

關鍵詞：吸收劑量；銀河宇宙射線；載人航天器；屏蔽材料；蒙特卡羅方法；數值分析

http://www.bisee.ac.cnE-mail: htqhjgc@126.comTel:(010)68116407, 68116408, 68116544

0　引言

隨著航天技術的發(fā)展，未來的載人深空探索活動將集中于月球以及火星甚至更遠的空間區(qū)域，航天器將會遭遇與近地軌道不同的輻射環(huán)境。由于失去了地磁場的保護，航天器直接暴露在銀河宇宙射線與太陽粒子輻射環(huán)境中。太陽粒子事件具有偶然性，與太陽活動有關，相對于宇宙射線來說更容易防護[1]；而銀河宇宙射線中，質子占約87%[2]，這些質子會危及航天器壽命以及航天員的生命健康，因而研究高能質子的輻射屏蔽效應非常有必要。瑞典D. Mancusi以及意大利A. Bertucci等人[3-4]對能量為1GeV的質子做了大量實驗，通過研究質子在航天器材料中的劑量深度分布來考察質子的輻射屏蔽效應。本文主要研究高能質子穿過航天器外殼屏蔽材料后在水模體中的劑量深度分布：水模體的吸收劑量代表了穿過航天器艙壁材料之后初始射線以及次級射線在艙內的總能量沉積，而航天器外殼屏蔽材料中的劑量深度分布代表的是射線在外殼及屏蔽材料中的能量沉積，并不能反映艙內航天員的輻射損傷。因此，本文通過三維的蒙特卡羅程序Geant4建立簡易載人航天器模型時，在屏蔽材料之后設置了一個水模體來模擬人體組織，以獲得水模體中的劑量深度分布曲線，然后采用積分的方法可以獲得高能質子在水模體中沉積的總能量以及次級粒子引起的能量沉積。

1　Geant4程序與載人航天器模型

1.1Geant4程序介紹

Geant4[5-8]是基于C++面向對象程序設計的一款開源軟件，它是由歐洲核子中心與日本高能所主導，大約100多位科學家參與完成的一款蒙特卡羅程序包，其主要特點是采用面向對象以及類的概念構建數據模型，通過編寫C++代碼調用合適的物理包，最后統(tǒng)計所需要的數據。Geant4包括許多內置的和用戶編制的精細物理模型，常用的粒子與物質相互作用的物理過程有電磁相互作用和強相互作用：電磁相互作用包括多次散射、電離、光電效應、軔致輻射等過程；強相互作用包括彈性散射、非彈性散射和原子核嬗變。Geant4模擬的粒子種類繁多，幾乎包括了自然界中的所有粒子?？傊?，相對于其他蒙特卡羅軟件，Geant4具有無可比擬的優(yōu)勢，因而在高能物理、重離子治療以及空間輻射效應等研究中都有著廣泛的應用。

1.2載人航天器模型

航天員的人體輻射劑量效應分析是一個十分復雜的問題，將真正的人體作為對象研究輻射粒子的能量沉積，目前并不是一個可行的辦法。為簡化起見，按照國際上通用的方法，用水模體代替人體復雜的生物組織。水是人體的重要組成部分，構成人體中大約70%的重量和血液中約83%的重量，研究表明，水在質量密度、電子密度和原子序數上都與肌肉組織非常接近，因而將水作為等效人體組織是可行的。

模擬計算時選取的模型為簡易圓柱體殼形航天器[9]（見圖1），是典型空間航天器通用模型，其內半徑為50cm，高為50cm。其中綠色代表多層航天器外殼材料，藍色代表水模體（為直徑30cm、高30cm的圓柱體）。表1是簡易航天器外殼材料的成分[9]，圖1還給出銀河宇宙射線質子穿過航天器外殼材料的粒子徑跡?？梢钥闯?，當高能量的質子與飛行器外殼相互作用時，會激發(fā)出大量的次級射線，其中綠線、紅線、藍線分別表示帶正電、負電以及不帶電的粒子。

表1　簡易載人航天器外殼材料成分Table1　The composition of spacecraft’s external structural material

圖1　簡易載人航天器模型及穿越粒子徑跡Fig. 1　Model of manned spacecraft and the pass-through charged particles

銀河宇宙射線質子通量占了整個銀河宇宙射線的87%，因而選取質子作為典型宇宙射線，其能譜由美國NRL主導編制的 CREME96模型獲得，計算了太陽活動周期為1年的銀河宇宙射線質子微分能譜（將1977年太陽活動極小年時作為參考時期），如圖2所示[2]。模擬計算時，調用Geant4中的通用源GPS功能。

圖2　銀河宇宙質子能譜Fig. 2　GCR proton energy spectrum

2　模擬結果與討論

2.1模擬結果

本文在模擬計算中，調用了Geant4中電磁相互作用物理包以及強相互作用物理包：選取的電磁作用為Opition3物理包，該模型對光子、電子、離子等與靶原子作用都有較好的輸運；強相互作用調用Geant4中自帶的QGSP_BIC_HP[10]物理包，該物理包囊括了原子核內級聯、核粒子蒸發(fā)、復合核、衰變等基本的核反應過程，能跟蹤次級質子、中子以及其他帶電粒子。因而當高能質子與飛行器外殼材料發(fā)生相互作用時，該物理包能較好地處理質子、中子等等次級粒子的輸運，從而保證計算結果的正確性。圖3為計算出的銀河宇宙射線質子穿過載人航天器艙后在水模體中的總劑量深度分布及次級粒子劑量深度分布；圖4為水模體總的吸收劑量與次級粒子吸收劑量直方分布；圖5是銀河宇宙射線質子穿過航天器在水模體中產生的不同次級粒子數量。

圖3　GCR質子穿過航天器在水模體中的總劑量與次級粒子劑量深度分布Fig. 3　Dose-depth distribution of total and secondary particles in the water phantom when GCR protons pass through the spacecraft

圖 4　水模體吸收總劑量與次級粒子劑量直方分布Fig. 4　Histogram of total dose and secondary particle dose distributions in the water phantom

圖5　GCR質子穿過航天器在水模體中產生的不同次級粒子數量Fig. 5　The different secondary particles generated in the water phantom when GCR protons pass through the spacecraft

從圖3中可以看出，總的吸收劑量隨著水模體深度的增加變化平緩，次級粒子引起的劑量分布與總劑量分布曲線類似。從圖4中可得出，次級粒子引起的吸收劑量所占比例達到50%，因此這部分吸收劑量不可忽視。從圖5中可以看出，銀河宇宙射線質子穿過航天器后，與航天器外殼材料相互作用，產生次級粒子數量最多可達到107量級。

2.2討論

高能質子與航天器外殼材料碰撞涉及的物理過程復雜，主要包括3種相互作用：電磁相互作用，放射性衰變作用，強相互作用。

1）電磁相互作用

主要是指初始粒子與靶核作用，當粒子靠近靶核時，速度降低，傳遞能量，進而損失能量。

2）放射性衰變作用

入射粒子通常為銀河宇宙射線或者太陽質子，當粒子能量大于1GeV時，可能產生次級中子與次級介子；不帶電的介子會衰變成光子，光子又能產生電子；而帶電的介子會衰變?yōu)棣套樱套永^而又衰變?yōu)檎撾娮印?/p>

3）強相互作用

空間環(huán)境中含有的射線種類復雜，特別是當能量較高的宇宙射線穿過靶核時，會發(fā)生核反應，產生各種不同的次級粒子。例如當高能量的質子或者中子與靶核反應，會產生靶核碎片；靶核碎片與飛船屏蔽材料發(fā)生非彈性散射，繼而產生大量的中子；中子不斷地累加，進入飛船內會嚴重威脅航天員的健康?？傊?，進入航天器艙內的次級輻射環(huán)境極其復雜。銀河宇宙射線質子入射航天器后，與航天器外殼材料相互作用發(fā)生核反應，因其能量較高，產生大量的次級粒子，與總吸收劑量相比，次級粒子引起的吸收劑量也較大。另外，因為次級質子與中子隨著屏蔽材料的厚度增加而增加，所以對高通量的質子較難屏蔽。

3　結束語

本文利用Geant4構建了簡易三維圓柱體空殼型航天器，模擬計算了高能量的銀河宇宙射線質子與航天器外殼材料相互作用，通過調用Geant4中電磁相互作用與強相互作用物理包，得到了高能量的宇宙質子穿過載人航天器后在水模體中的總劑量以及次級粒子引起的劑量深度分布。計算結果表明，次級粒子引起的劑量占總吸收劑量比例為50%，在水模體的電子總數量達到了107量級。

本文建立的載人航天器模型可為我國將來的深空輻射探測提供借鑒，也可為空間輻射效應分析及輻射防護提供一定的數據參考。

致謝

本工作得到了美國NASA馬歇爾空間飛行中心林子威教授以及復旦大學核科學與技術系鄭川老師的幫助與支持，在此一并表示衷心感謝。

參考文獻（References）

[1]Spillantini P, Casolino M, Durante M, et al. Shielding from cosmic radiation for interplanetary missions: active and passive methods[J]. Radiation Measurements, 2007, 42(1): 14-23

[2]Bernabeu J, Casanova I. Geant4-based radiation hazard assessment for human exploration missions[J]. Advances in Space Research, 2007, 40(9): 1368-1380

[3]Mancusi D, Bertucci A, Gialanella G, et al. Comparison of aluminum and lucite for shielding against 1GeV protons[J]. Advances in Space Research, 2007, 40(4): 581-585

[4]Bertucci A, Durante M, Gialanella G, et al. Shielding of relativistic protons[J]. Radiation and Environmental Biophysics, 2007, 46(2): 107-111

[5]Asai M, Dotti A, Verderi M, et al. Recent developments in Geant4[J]. Annals of Nuclear Energy, 2015, 82: 19-28

[6]Bao G. The simulation of radiation effects to astronauts due to solar energetic particles in deep space[J]. Acta Astronautica, 2012, 71: 32-37

[7] Bao G. The simulation of GCRs at near earth interplanetary locations for radioprotection[J]. Science China: Physics, Mechanics and Astronomy, 2010, 53(9): 1739-1745

[8]Ersmark T, Carlson P, Daly E, et al. Geant4 Monte Carlo simulations of the galactic cosmic ray radiation environment on-board the International Space Station/ Columbus[J]. IEEE Transactions on Nuclear Science, 2007, 54(5): 1854-1862

[9]Martinez L M, Kingston J. Space radiation analysis: radiation effects and particle interaction outside the Earth's magnetosphere using GRAS and GEANT4[J]. Acta Astronautica, 2012, 72: 156-164

[10]Lin Z W, Adams J H, Barghouty A F, et al. Comparisons of several transport models in their predictions in typical space radiation environments[J]. Advances in Space Research, 2012, 49(4): 797-806

（編輯：閆德葵）

Analysis of radiation environment and its dose within sealed cabin of manned spacecraft

Wu Zhengxin1, Sun Huibin1,2, He Chengfa3, Tong Yongpeng2, Ma Yugang1, Lu Jingbin1, Xun Mingzhu3, Hu Yanqi4, Cai Zhenbo4, Zhu Zhipeng2
(1. College of Physics, Jilin University, Changchun 130012, China;
2. Institute of Nuclear Technology, Shenzhen University, Shenzhen 518060, China;
3. Xinjiang Technical Institute of Physics and Chemistry, China Academy of Sciences, Urumqi 830011, China; 4. Beijing Institute of Spacecraft System Engineering, Beijing 100094, China)

Abstract:By using the Geant4 toolkit, we have modeled the interaction process of the most abundant galactic cosmic ray particle fluxes(protons)injected into the water phantom in a manned spacecraft model. The total absorbed dose due to the protons and the absorbed dose due to the secondary particles generated by the radiation are estimated, as well as the total number of different kinds of secondary particles. Final results indicate that the proportion of the absorbed dose due to the secondary particles generated accounts for 50% of the total dose, and the total amount of the secondary electrons reaches the order of 107.

Key words:absorbed dose; GCR; manned spacecrafts; shield material; Monte Carlo method; numerical analysis

作者簡介：吳正新（1987—），男，博士研究生，研究方向為空間輻射效應的蒙特卡羅模擬；E-mail:wuzx13@mails.jlu.edu.cn。通信作者：孫慧斌（1957—），男，教授，博士生導師，研究方向為核物理及核技術應用；E-mail: hbsun@szu.edu.cn。

基金項目：國家重大科學儀器設備開發(fā)專項基金（批準號：2012YQ240121）；國家自然科學基金（編號：11575118；11375117）

收稿日期：2015-10-05；修回日期：2016-03-12

DOI:10.3969/j.issn.1673-1379.2016.02.007

中圖分類號：TL99

文獻標志碼：A

文章編號：1673-1379(2016)02-0154-04