中国辐射卫生  2020, Vol. 29 Issue (4): 362-365  DOI: 10.13491/j.issn.1004-714X.2020.04.010

引用本文 

李明生, 姜庆寰, 程金生. 利用Geant4计算生物样品照射时的剂量差异[J]. 中国辐射卫生, 2020, 29(4): 362-365. DOI: 10.13491/j.issn.1004-714X.2020.04.010.
LI Mingsheng, JIANG Qinghuan, CHENG Jinsheng. Calculation of the dose deviation of biological samples irradiation with Geant4[J]. Chinese Journal of Radiological Health, 2020, 29(4): 362-365. DOI: 10.13491/j.issn.1004-714X.2020.04.010.

基金项目

国家传染病防治重大专项(2018ZX10734401-014)

通讯作者

程金生,E-mail:chengjs3393@163.com

文章历史

收稿日期:2020-02-26
利用Geant4计算生物样品照射时的剂量差异
李明生 , 姜庆寰 , 程金生     
中国疾病预防控制中心辐射防护与核安全医学所 辐射防护与核应急中国疾控中心重点实验室,北京 100088
摘要目的 评估生物样品照射时,容器中间位置处剂量与容器内溶液总体剂量的差异。方法 利用Geant4程序模拟4种照射模式下,容器内的中心处剂量和容器内溶液总剂量,并计算两者的比值。结果 容器中心处剂量与溶液总体剂量的比值与射线能量、容器大小以及照射方式紧密相关。结论 生物样品照射实验时,应该根据现有条件选择合适的容器、射线以及照射模式以便中心剂量与总体剂量近似相等。
关键词生物样品照射    辐射剂量    蒙特卡罗模拟    
Calculation of the dose deviation of biological samples irradiation with Geant4
LI Mingsheng , JIANG Qinghuan , CHENG Jinsheng     
National institute of radiation protection China CDC, Beijing 100088 China
Abstract: Objective To irradiate biological samples and evaluate the difference between the dose in the middle of the container and the total dose of the solution in the container. Methods Simulate the central dose in the container and the total dose of the solution in the container under four irradiation modes with Geant4, and then calculate the ratio of the above two dose. Results The ratio is closely related to the ray energy, the size of the container, and the modes of irradiation. Conclusion when irradiating biological samples, appropriate containers, rays, and irradiating patterns should be selected according to the existing conditions so that the central dose is approximately equal to the overall dose.
Key words: Biological Sample Irradiation    Radiation Dose    Monte Carlo Simulation    
1 材料和方法 1.1 生物容器

生物样品照射中,常常是物理技术人员了解所照射位置处的空气吸收剂量或照射量,而对于各种规格容器中的样品总剂量无从判断;而对于生物技术人员,对于剂量的了解则更少,往往只能提出一个模糊的概念,比如需要多少照射剂量。这样的结果导致照射时,辐照技术人员只能给出容器中间位置处的空气吸收剂量,把该值作为整个容器内的吸收剂量值。但是由于生物容器的多样性,以及材质的影响,中间位置处的剂量与总体剂量会有不同,再结合射线能量以及照射方式不同导致的这个差异可能会更大。鉴于此采用了Geant4程序[1-2]模拟了4中常用规格生物容器,在不同射线和不同照射方式下容器内的剂量分布,并计算了中间位置处剂量与总体剂量的比值。

模拟实验采用了4中规格的生物容器,体积分别为1.5 mL、4 mL、15 mL和50 mL,其中1.5 mL容器的材质为聚苯乙烯;4 mL容器壁的材质为聚苯乙烯,帽为丁基橡胶;15 mL容器壁材质为聚苯乙烯,帽材质为聚乙烯;50 mL容器壁材质为PET塑料,帽材质为聚乙烯。

1.2 照射方式

采用了4种照射方式,分别为竖直照射和水平照射,其中每种照射方式对应着两种溶液状态,满管溶液和半管溶液。如图1所示。

图 1 照射方式示意图
1.3 射线能量

采用电子线和γ、X射线来照射容器,选取的能量分别为电子线[3-4]:8 MeV、10 MeV、12 MeV、15 MeV、18 MeV和22 MeV;γ和X射线[5-6]分别选取能量为:0.662 MeV、1.25 MeV、2 MeV、3 MeV和5 MeV,其中0.662 MeV为137Cs产生的γ射线,1.25 MeV为60Co产生的γ射线,2 MeV、3 MeV和5 MeV为加速器产生的X射线。

1.4 Geant4建立模型

本文采用Geant4来模拟建立容器模型,溶液为水,如图2所示,剂量利用Geant4的scoring mesh方法[7-10]来获得容液内的剂量分布曲线,如图3所示,1.5 mL容器竖直选取了22层,水平选取10层,每层1 mm;4 mL容器竖直选取了70层,水平选取10层,每层1 mm;15 mL容器竖直选取了100层,水平选取20层,每层1 mm;50 mL容器竖直选取了100层,水平选取20层,每层1 mm;mesh中心与容器中点相一致。总剂量是计算容器内溶液整体受到的剂量值[11]。最后取剂量曲线中间位置处的剂量与溶液总剂量相比(以下简称比值),这个比值反应了在生物样品照射中预期剂量与实际剂量的差异。如果中间位置处剂量与溶液总剂量相等,则比值为1,如果大于总剂量则比值大于1,如果小于总剂量则比值小于1。剂量的差异通过比值相对于1的大小来反映。

图 2 模型与实物对比图

图 3 水平和竖直模式下scoring mesh示意图
2 结果

利用不同射线对四种容器分别采用四种照射方式进行模拟计算,其中剂量值是与入射粒子数相关的相对值。1.5 mL容器比值与射线的关系详见图4;4 mL容器内比值与射线的关系详见图5;15 mL容器比值与射线关系详见图6;50 mL容器比值与射线关系详见图7表1表2表3表4分别与图4图5图6图7相对应,表中列出了在四种照射模式下相应容器的最大比值和最小比值以及平均比值和标准差。平均比值是在固定容器与照射方式的情况下,比值随着射线能量变化的情况,标准差反映了比值的波动大小。

表 1 1.5 mL容器下的比值情况

表 2 4 mL容器下的比值情况

表 3 15 mL容器下的比值情况

图 4 1.5 mL容器比值与射线以及照射方式的关系图

表 4 50 mL容器下的比值情况

图 5 4 mL容器比值与射线以及照射方式的关系图

图 6 15 mL容器比值与射线以及照射方式的关系图

图 7 50 mL容器比值与射线以及照射方式的关系图
3 讨论

上文总结了比值随着容器大小、照射方式与射线能量的变化趋势。从容器大小来说,为了分析容器对比值的影响,对每种容器4中照射模式下的平均比值进行了再平均,得到1.5 mL对应的比值为1.04 ± 0.18;4 mL对应的比值为0.98 ± 0.16;15 mL对应的比值为1.04 ± 0.34;50 mL对应的比值为1.18 ± 0.32。据此可以说,采用4 mL容器进行生物样品照射,比值最接近1,且标准差最小,容器太大或者太小,比值都不理想。

从照射方式来说,对不同照射方式下的平均值进行了再平均,其中竖直满管照射模式下比值为1.03 ± 0.28;竖直半管模式下比值为1.29 ± 0.28;水平满管照射模式下为1.07 ± 0.01;水平半管照射模式下为0.83 ± 0.02。因此可以说水平满管照射模式下比值最接近1,且标准差小,是生物样品照射的理想模式。

射线选择应该与照射模式一起考虑,射线分为电子线和X、γ射线,对电子线4种照射模式下的比值进行平均,其中竖直满管照射模式下比值为0.92 ± 0.41;竖直半管模式下比值为1.71 ± 0.39;水平满管照射模式下为1.04 ± 0.01;水平半管照射模式下为0.98 ± 0.02。对X、γ射线4种照射模式下的比值进行平均,其中竖直满管照射模式下比值为1.17 ± 0.20;竖直半管模式下比值为0.79 ± 0.20;水平满管照射模式下为1.11 ± 0.02;水平半管照射模式下为0.64 ± 0.03。因此在选择电子线照射时,应该采用水平半管模式;选择X、γ射线时,应该采用水平满管模式。

在进行生物样品照射实验时,应该综合考虑容器大小、射线类型和照射模式,选择合适的照射方式以便使得比值最接近1。

参考文献
[1]
J Allison, K Amako, J Apostolakis, et al. Geant4 developments and applications[J]. IEEE Trans Nucl Sci, 2006, 53(1): 270-278. DOI:10.1109/TNS.2006.869826
[2]
S. Agostinelli, J. Allison, K. Amaker, et al. Geant4: a Simulation Toolkit[J]. Nucl. Instrum. Meth. A, 2003, 506(3): 250-303. DOI:10.1016/S0168-9002(03)01368-8
[3]
洪天祺, 戴瑜. 0.5~2.5 MeV辐射加工电子加速器X射线发射率的MCNP5计算[J]. 中国辐射卫生, 2015, 24(3): 293-295. DOI:10.13491/j.cnki.issn.1004-714x.2015.03.045
[4]
邓磊, 周宁, 邱诗宝, 等. 江西省医用电子加速器性能检测与分析[J]. 中国辐射卫生, 2017, 26(5): 552-553. DOI:10.3969/j.issn.1004-714X.2017.05.015
[5]
徐洋洋, 邓磊, 杨亚新, 等. 放射治疗外照射的蒙特卡罗计算[J]. 中国辐射卫生, 2016, 25(1): 15-17.
[6]
李建国, 尹晶晶, 袁慧, 等. 不同剂量γ射线照射人肝细胞株对BMPR2基因表达的影响研究[J]. 中国辐射卫生, 2013, 22(1): 18-21. DOI:10.13491/j.cnki.issn.1004-714x.2013.01.058
[7]
Choi WQ, Cho KY, Yeo IS, et al. Performance profiling for brachytherapy applications[J]. Computer Physics Communication, 2018, 226: 180-186. DOI:10.1016/j.cpc.2017.12.022
[8]
Han M C, Ku Y, Lee H S, et al. New calculation method for 3D dose distribution in tetrahedral-mesh phantoms in Geant4[J]. Phys Med, 2019, 66: 97-103. DOI:10.1016/j.ejmp.2019.09.239
[9]
Jiang H, Paganetti H. Adaptation of GEANT4 to Monte Carlo dose calculations based on CT data[J]. Med Phys, 2004, 31(10): 2811-2818. DOI:10.1118/1.1796952
[10]
Landry G, Reniers B, Pignol J P, et al. The difference of scoring dose to water or tissues in Monte Carlo dose calculations for low energy brachytherapy photon sources[J]. Med Phys, 2011, 38(3): 1526-1533. DOI:10.1118/1.3549760
[11]
孙亮, 李士骏. 电离辐射剂量学基础[M]. 北京: 中国原子能出版社, 2014: 24-35.