BNCT
硼中子俘獲治療(Boron Neutron Capture Therapy,簡稱BNCT)。硼中子俘獲治療通過在腫瘤細胞內(nèi)的原子核反應(yīng)來摧毀癌細胞。它的原理是這樣的:先給病人注射一種含硼的特殊化合物,這種化合物與癌細胞有很強的親和力,進入人體后,迅速聚集于癌細胞內(nèi),而其他組織內(nèi)分布很少。這種含硼化合物對人體無毒無害,對癌癥也無治療作用。這時,用一種中子射線進行照射,這種射線對人體的損傷不大,但中子與進入癌細胞里的硼能發(fā)生很強的核反應(yīng),釋放出一種殺傷力極強的射線,這種射線的射程很短,只有一個癌細胞的長度。所以只殺死癌細胞,不損傷周圍組織。這種有選擇的只殺死形狀復(fù)雜的癌細胞而不損傷正常組織的技術(shù),稱為硼中子俘獲治療技術(shù)。

基于Geant4的硼中子俘獲治療的能量沉積分布研究
周小平,張 巍
(中國原子能科學(xué)研究院,北京 102413)
●●●摘 要
硼中子俘獲治療(BNCT)是一種治療癌癥的有效手段。采用國際通用的Geant4開發(fā)程序包,以Snyder模型為對象,構(gòu)建中子輸運計算模型,評估在給定中子束流能譜和入射方向條件下,腦組織吸收不同10B質(zhì)量比的含硼藥物對癌變組織及正常組織的能量沉積的來源和分布情況。研究表明,10B引起的能量沉積呈先升高后降低的變化趨勢,在入射深度為3.0 cm處達到峰值;在入射深度為2~5 cm區(qū)間內(nèi),能量沉積主要來自10B;10B引起的能量沉積與10B質(zhì)量比呈線性關(guān)系;光子、質(zhì)子、電子及其他粒子引起的能量沉積占比很小,且不隨10B質(zhì)量比的變化而變化。研究結(jié)果對BNCT的方案設(shè)計、效果評估具有參考意義。
●●●引 言
硼中子俘獲治療(boron neutron capture therapy,BNCT)主要通過將具有選擇性的含硼藥物注射入人體血液,然后富集于腫瘤細胞中,受中子照射后,經(jīng)過10B(n,α)7Li反應(yīng)放出的高能量、短射程的α粒子和7Li粒子,來殺死腫瘤細胞,而不傷害正常組織[1-3]。
本文以修正的Snyder模型為對象,采用中心點算法[4],網(wǎng)格尺寸為2.0 mm×2.0 mm×2.0 mm,僅考慮5種材料(皮膚、顱骨、腦組織、癌變組織和空氣),假設(shè)網(wǎng)格中心點處的材料成分和密度為整個網(wǎng)格的材料成分和密度。采用Geant4軟件工具包模擬加速器束流中子在人體組織內(nèi)的輸運行為,統(tǒng)計各種次級粒子引起的能量沉積,分析束流中子對癌變組織及正常組織的輻射損傷。
Geant4是歐洲核子研究中心(European Organization for Nuclear Research,CERN)開發(fā)的面向?qū)ο缶幊痰能浖ぞ甙诿商乜宸椒M各種粒子在物質(zhì)中的運動行為,廣泛用于高能物理、核能實驗、醫(yī)學(xué)、加速器和空間物理學(xué)研究等領(lǐng)域,已被全球許多研究項目所使用。
●●●能量沉積
當(dāng)量劑量是組織或器官接受的平均吸收劑量乘以輻射權(quán)重因子后得到的乘積。其中,吸收劑量是指單位質(zhì)量物質(zhì)受輻射后吸收輻射的能量沉積,其國際制單位為Gy(1 Gy=1 J/kg)。
在BNCT治療中,能量沉積的來源包括:α粒子、7Li粒子、光子、電子、質(zhì)子、其他重核等。主要的能量沉積來自于10B(n,α)7Li反應(yīng)放出的高能量、短射程的α粒子和7Li粒子,其反應(yīng)關(guān)系式為:

其中:
不同次級粒子引起的能量沉積對吸收劑量的輻射權(quán)重因子不同[5-6],如表1所示。在Geant4模擬計算中,中子不直接產(chǎn)生能量沉積,而是通過其他次級粒子作用引起能量沉積。
表1 次級粒子的輻射權(quán)重因子
Table 1 Radiation weight factor of secondary particle

含硼藥物(10B)引起的能量沉積由α粒子(僅指由中子與10B反應(yīng)所產(chǎn)生的α粒子)和7Li粒子引起的能量沉積之和。
●●●模型介紹
1、Snyder模型
目前,理論計算普遍采用的頭部模型是修正的Snyder橢球模型[3-4]。它是由3個橢球組成,將整個頭部組織由外至內(nèi)分為頭皮、顱骨和腦組織三個部分,如圖1所示。其中,腦組織在+Z方向有1 cm的偏心。各分界面滿足下列3個橢球方程(1)、(2)、(3)。
(1)頭皮與空氣:

(2)顱骨與頭皮:

(3)腦組織與顱骨:

橢球模型被空氣包圍,整個空間被分為4個部分,相應(yīng)的生物組織分別為:空氣、皮膚、顱骨、腦組織,其密度分別為1.293×10-3、1.040、1.610、1.090 g/cm3。4種組織的元素成分由ICRU46給出,如表2所示。
圖1 Snyder模型剖面圖
Fig.1 the profile map of Snyder model

表2 人體組織元素成分
Table 2 The element composition of human tissue 單位:%

2、計算方案
假定入射中子能譜為加速器束流質(zhì)子轟擊鈹靶后經(jīng)慢化體慢化整形后超熱中子譜[10],如圖2所示。
圖2 束流中子能譜
Fig.2 Spectrum of neutron beam

假設(shè)中子從某個固定點、向固定方向發(fā)射,出射方向為-Z方向,初始出射位置為(0 cm,0 cm,10 cm)。
僅考慮腦組織吸收不同10B質(zhì)量比的含硼藥物,其他組織成分不變。這里的10B質(zhì)量比,是指10B原子占組織成分的質(zhì)量分?jǐn)?shù)。含硼藥物中10B質(zhì)量比分別為0、10、20、30、40 mg/kg。
同時考慮在(0,0,5 cm)位置有一個腫瘤組織,其10B質(zhì)量比為40 mg/kg,正常組織的10B質(zhì)量比為10 mg/kg,如圖3所示。其分界面滿足橢球方程(4)。

圖3 含癌變組織的Snyder模型(YZ)
Fig.3 Snyder model with canceration issue(YZ view)

●●●結(jié)果分析
Snyder模型統(tǒng)計了腦組織吸收不同含硼質(zhì)量比情況下的次級粒子能量沉積,得到了各次級粒子(α粒子,7Li粒子,電子,光子,質(zhì)子及其他粒子)沿入射中子方向(Z-方向)的單位中子引起的能量沉積分布情況[6-7]。初始計算粒子數(shù)目為2.0×109個,各組織區(qū)域的能量沉積分布平均方差不超過2%。
當(dāng)10B質(zhì)量比為0時,其次級粒子(電子,光子,質(zhì)子及其他粒子)引起的沉積能量和總的能量沉積的分布情況(其中光子引起的沉積能量對應(yīng)于次縱坐標(biāo)),如圖4所示。結(jié)果表明,光子引起的能量沉積低,可以忽略不計;質(zhì)子、其他重核引起的能量沉積次之,占比很低;當(dāng)腦組織含10B質(zhì)量比為0 mg/kg時,能量沉積主要由電子引起。
圖4 次級粒子能量沉積分布
Fig.4 Energy deposition distribution of secondary particles

電子、光子、質(zhì)子和其他粒子引起的能量沉積,不隨含硼藥物質(zhì)量比增加而增加,如圖5、圖6、圖7、圖8所示。
圖5 電子引起的能量沉積分布
Fig.5 Distribution of energy deposition caused by electron

沿入射中子方向(Z-)的總能量沉積和10B引起的能量沉積分布情況如圖9所示。結(jié)果表明,10B引起的能量沉積隨硼質(zhì)量比的增加而增加,且占總能量沉積的比例很大;能量沉積呈先上升后下降的變化趨勢,且在3.0 cm處,能量沉積達到峰值;當(dāng)入射深度在2.0~5.0 cm時,10B引起的能量沉積占比超過80%,峰值占比達到96.59%,如表3所示。
圖6 質(zhì)子引起的能量沉積分布
Fig.6 Distribution of energy deposition caused by protons

圖7 光子引起的能量沉積分布
Fig.7 Distribution of energy deposition caused by photons

圖8 其他粒子引起的能量沉積分布
Fig.8 Distribution of energy deposition caused by other particles

單位10B質(zhì)量比引起的能量沉積分布情況,如圖10。結(jié)果表明,單位10B質(zhì)量比引起的能量沉積隨含硼藥物質(zhì)量比略微下降,10B引起的能量沉積與含硼藥物質(zhì)量比呈線性正相關(guān),說明提高含硼藥物質(zhì)量比能有效增加癌變組織的輻射損傷效果。
圖9 不同含10B質(zhì)量比藥物對總的和10B引起的能量沉積的影響
Fig.9 Effects of different concentrations of boron-containing drugs on energy deposition caused by total and B-10

圖10 單位10B質(zhì)量比的能量沉積分布
Fig.10 Energy deposition distribution of B-10 per unit concentration

●●●結(jié) 論
本文采用Geant4開發(fā)程序包,網(wǎng)格化Snyder模型,構(gòu)建中子輸運計算模型,評估加速器束流中子照射不同質(zhì)量比含硼藥物的腦組織條件下的能量沉積的來源和分布情況。
表3 10B能量沉積在總的能量沉積中的占比
Table 3 The proportion of boron-10 energy deposition in total energy deposition

研究表明,癌變組織吸收一定質(zhì)量比的含硼藥物,在一定照射深度范圍,能有效提高癌變組織的能量沉積,從而達到治療效果;10B引起的能量沉積與10B質(zhì)量比呈線性關(guān)系,應(yīng)在后續(xù)BNCT研究中提高癌變組織的含硼藥物質(zhì)量比,同時抑制正常組織的含硼藥物質(zhì)量比。
研究結(jié)果對BNCT的方案設(shè)計、效果評估具有參考意義,同時表明,Geant4軟件開發(fā)包可以用于BNCT的方案設(shè)計和效果評估。
成都瑞庚醫(yī)療科技咨詢有限公司
公司網(wǎng)址丨www.fortisrecovery.com
公司郵箱丨rg@ruigeng.net
洽詢電話丨028-61985320
微信客服丨18982287599
微信公眾號丨瑞庚臺醫(yī)

赴臺醫(yī)療
疫苗接種
赴臺培訓(xùn)
赴臺考察
赴臺交流
赴臺健檢
臺灣養(yǎng)老
兩岸論壇
兩岸合作