BNCT治疗脑胶质瘤剂量特性蒙特卡罗模拟研究

doi:10.27045/j.1000-8187.202603008

摘要/Abstract

摘要： 本文以西安脉冲反应堆中子源为基础,采用蒙特卡罗方法,对脑胶质瘤BNCT治疗中的束流整形组件(beam shaping assembly,BSA)设计及剂量特性开展了研究。在满足IAEA推荐束流指标的约束条件下,对慢化体、反射体、准直系统、热中子吸收体及γ吸收体等关键结构参数进行了优化,形成了以AlF₃慢化体为核心的超热中子束流方案。进一步构建Snyder椭球脑模型及含肿瘤组织模型,计算并分析了不同硼浓度条件下硼俘获剂量、热中子剂量、快中子剂量及γ剂量的空间分布特征。结果表明,采用约60 cm的AlF₃慢化体、约15 cm的Pb反射体、50 cm准直器、Gd热中子吸收体(0.3 cm/0.1 cm)及5 cm Bi γ吸收体的BSA组合,可获得满足治疗要求的超热中子通量;硼俘获剂量在肿瘤区形成显著峰值,并随硼浓度升高而逐渐增强,而非选择性剂量分量基本不受影响。通过总生物剂量评估,确定了不同硼浓度条件下的剂量率及相应治疗时间。上述结果为反应堆型BNCT束流设计及脑胶质瘤治疗的剂量优化提供了一定的数据支持。

关键词: 硼中子俘获治疗, 脑胶质瘤, 吸收剂量率, Snyder模型, 蒙特卡罗模拟

Abstract: In this study, based on the neutron source of the Xi′an Pulsed Reactor, Monte Carlo simulations were performed to systematically investigate the beam shaping assembly (BSA) design and dosimetric characteristics for BNCT treatment of brain glioma. Under the constraints of the beam quality criteria recommended by the IAEA, key components of the BSA—including the moderator, reflector, collimation system, thermal neutron absorber, and gamma filter—were optimized, resulting in an epithermal neutron beam configuration centered on an AlF₃ moderator. Furthermore, a Snyder ellipsoidal head phantom and a corresponding tumor-bearing model were constructed to calculate and analyze the spatial distributions of boron capture dose, thermal neutron dose, fast neutron dose, and gamma dose under different boron concentrations. The results indicate that a BSA configuration consisting of an AlF₃ moderator (~60 cm), a Pb reflector (~15 cm), a 50 cm collimator, Gd thermal neutron absorber (0.3 cm/0.1 cm), and a 5 cm Bi gamma filter can provide an epithermal neutron flux that satisfies therapeutic requirements. A pronounced peak of boron capture dose is formed within the tumor region and increases significantly with increasing boron concentration, while the non-selective dose components remain essentially unchanged. Based on the evaluation of the total biologically weighted dose, the corresponding dose rates and treatment times for different boron concentrations were determined. These results provide a reliable theoretical basis for reactor-based BNCT beam design and dosimetric optimization in the treatment of brain glioma.

Key words: boron neutron capture therapy, glioma, absorbed dose rate, snyder model, Monte Carlo simulation

中图分类号:

TL99

李波宁, 王晓娟, 方青龙, 原林, 王博宇, 李茁, 刘洋. BNCT治疗脑胶质瘤剂量特性蒙特卡罗模拟研究[J]. 辐射防护, 2026, 46(3): 250-263.

LI Boning, WANG Xiaojuan, FANG Qinlong, YUAN Lin, WANG Boyu, LI Zhuo, LIU Yang. Monte Carlo simulation study on dosimetry of BNCT treatment for glioma[J]. RADIATION PROTECTION, 2026, 46(3): 250-263.

参考文献

[1] 国家卫生健康委员会医政医管局,中国抗癌协会脑胶质瘤专业委员会,中国医师协会脑胶质瘤专业委员会.脑胶质瘤诊疗指南(2022版)[J].中华神经外科志,2022,38(8):757-777.
Medical Administration Bureau of National Health Commission, Professional Committee of Glioma of China Anti-Cancer Association, Professional Committee of Glioma of Chinese Medical Doctor Association. Guidelines for the Diagnosis and Treatment of Gliomas (2022 Edition)[J]. Chinese Journal of Neurosurgery, 2022, 38(8): 757-777.
[2] 邱有恒, 邓力, 应阳君. BNCT蒙特卡罗剂量计算的混合网格算法研究[J]. 高能物理与核物理, 2004, 28(6): 669-674.
QIU Youheng, DENG Li, YING Yangjun. Mix mesh sizes arithmetic research of BNCT dose rate calculations by Monte Carlo method[J]. High Energy Physics and Nuclear Physics, 2004, 28(6): 669-674.
[3] 江海燕. 硼中子俘获治疗中细胞损伤的Monte Carlo模拟与分析[D]. 合肥: 合肥工业大学, 2004.
JIANG Haiyan. Monte Carlo simulation and analysis of radiation damage to cells in boron neutron capture therapy[D]. Hefei: Hefei University of Technology, 2004.
[4] 孙爱扣, 陈珍平, 匡蓝珺, 等. 硼中子俘获治疗剂量蒙特卡罗计算程序Magic开发与验证[J]. 现代应用物理, 2023, 14(4): 39-46, 52.
SUN Aikou, CHEN Zhenping, KUANG Lanjun, et al. Development and verification of a Monte Carlo calculation program magic for dose of boron neutron capture therapy[J]. Modern Applied Physics, 2023, 14(4): 39-46, 52.
[5] 江新标, 朱广宁, 王道华, 等. 硼中子俘获治疗脑胶质瘤超热中子束的品质参数研究[J]. 中华神经外科疾病研究杂志, 2005, 4(3): 252-256.
JIANG Xinbiao, ZHU Guangning, WANG Daohua, et al. Study on quality parameters of epithermal neutron for boron neutron capture therapy of human gliomas[J]. Chinese Journal of Neurosurgical Disease Research, 2005, 4(3): 252-256.
[6] 江新标, 赵英, 曹锐锋, 等. 铀氢锆堆中子源辐照含硼C6胶质瘤细胞的生物学效应研究[J]. 核技术, 2005, 28(6): 458-462.
JIANG Xinbiao, ZHAO Ying, CAO Ruifeng, et al. Biological effects of boron-absorbed C6 glioma cells irradiated by UHZr reactor neutron source[J]. Nuclear Techniques, 2005, 28(6): 458-462.
[7] 刘书焕, 江新标, 仲云红, 等. 反应堆快中子实验装置辐射场参数测量[J]. 核技术, 2006, 29(9): 660-664.
LIU Shuhuan, JIANG Xinbiao, ZHONG Yunhong, et al. Measurement the parameters of the irradiation field for the fast neutron device in reactor[J]. Nuclear Techniques, 2006, 29(9): 660-664.
[8] Ganjeh Z A, Eslami-Kalantari M. Design and optimization of two-sided beam based on ⁷Li(p,n)⁷Be source using in BNCT for brain and liver tumors[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2019, 916: 290-295.
[9] Rasouli F S, Masoudi S F. Design and optimization of a beam shaping assembly for BNCT based on D-T neutron generator and dose evaluation using a simulated head phantom[J]. Applied Radiation and Isotopes, 2012, 70(12): 2755-2762.
[10] Fantidis J G, Antoniadis A. Optimization study for BNCT facility based on a DT neutron generator[J]. International journal of radiation research, 2014, 13(1): 13-24.
[11] Lee J C, CHEN Yiwei, Chuang K S, et al. The dosimetric impact of shifts in patient positioning during boron neutron capture therapy for brain tumors[J]. BioMed Research International, 2018, 2018: 5826174.
[12] Monshizadeh M, Kasesaz Y, Khalafi H, et al. MCNP design of thermal and epithermal neutron beam for BNCT at the Isfahan MNSR[J]. Progress in Nuclear Energy, 2015, 83: 427-432.
[13] Hu N, Tanaka H, Kakino R, et al. Evaluation of a treatment planning system developed for clinical boron neutron capture therapy and validation against an Independent Monte Carlo dose calculation system[J]. Radiation Oncology, 2021, 16(1): 243.
[14] Report I .Photon, electron, proton and neutron interaction data for body tissues[M].International Commission on Radiation Units and Measurements,1992.
[15] Sakamoto S, Kiger W S, Harling O K. Sensitivity studies of beam directionality, beam size, and neutron spectrum for a fission converter-based epithermal neutron beam for boron neutron capture therapy[J]. Medical Physics, 1999, 26(9): 1979-1988.
[16] Koivunoro H, Bleuel D L, Nastasi U, et al. BNCT dose distribution in liver with epithermal D-D and D-T fusion-based neutron beams[J]. Applied Radiation and Isotopes, 2004, 61(5): 853-859.
[17] 朱养妮, 江新标, 张强, 等. 硼中子俘获治疗人头开颅模型内宏观吸收剂量分布[J]. 原子能科学技术, 2009, 43(S2): 404-406.
ZHU Yangni, JIANG Xinbiao, ZHANG Qiang, et al. Analysis of absorbed dose distribution in equivalent model of head in case of operation for boron neutron capture therapy[J]. Atomic Energy Science and Technology, 2009, 43(S2): 404-406.
[18] Kawabata S, Miyatake S I, Kuroiwa T, et al. Boron neutron capture therapy for newly diagnosed glioblastoma[J]. Journal of radiation research, 2009, 50(1): 51-60.
[19] 鱼红亮, 郑传城, 孙亮. MCNP计算含肿瘤Snyder修正头部模型的硼中子俘获治疗剂量[J]. 原子能科学技术, 2010, 44(1): 89-94.
YU Hongliang, ZHENG Chuancheng, SUN Liang. Boron neutron capture therapy dose calculation for tumor modified Snyder head phantom using MCNP[J]. Atomic Energy Science and Technology, 2010, 44(1): 89-94.
[20] 刘栋, 孙光爱, 彭述明. 中子散射技术在化学研究中的应用[J]. 中国科学: 化学, 2023, 53(11): 2237-2249.
LIU Dong, SUN Guangai, PENG Shuming. Applications of neutron scattering techniques in investigations of chemistry[J]. Scientia Sinica Chimica, 2023, 53(11): 2237-2249.
[21] Rasouli F S, Masoudi S F. Simulation of the BNCT of brain tumors using MCNP code: beam designing and dose evaluation[J]. 2012, 9(3): 183-192.
[22] 邢星河, 赵亚丽, 鲁维, 等. 基于ICRU模型的辐射剂量计算[J]. 国际航空航天科学, 2020, 8(4): 125-133.
XING Xinghe, ZHAO Yali, LU Wei, et al. Radiation dose calculation based on ICRU model[J]. Journal of Aerospace Science and Technology, 2020, 8(4): 125-133.
[23] 沈江燕, 闫聪冲. 空间辐射粒子致航天员辐射剂量的蒙特卡罗模拟计算[J]. 辐射防护, 2023, 43(S1): 8-13.
SHEN Jiangyan, YAN Congchong. Monte Carlo simulations of radiation dose of astronauts caused by space radiation particles[J]. Radiation Protection, 2023, 43(S1): 8-13.
[24] Goorley J T, Kiger W S, Zamenhof R G. Reference dosimetry calculations for neutron capture therapy with comparison of analytical and voxel models[J]. Medical Physics, 2002, 29(2): 145-156.
[25] 邱有恒, 邓力, 应阳君, 等. 硼中子俘获治疗的蒙特卡罗方法模拟(英文)[J]. 高能物理与核物理, 2003, 27(10): 936-942.
QIU Youheng, DENG Li, YING Yangjun, et al. Monte Carlo simulation of boron neutron capture therapy[J]. High Energy Physics and Nuclear Physics, 2003, 27(10): 936-942.
[26] Ishiyama S. Deterministic parsing model of the Compound Biological Effectiveness ( CBE ) factor for intracellular ¹⁰boron distribution in boron neutron capture therapy[J]. Journal of Cancer Therapy, 2014, 5(14): 1388-1398.