孫愛(ài)扣,陳珍平,3?,匡藍(lán)珺,高克坤,劉程偉,于 濤

(1. 南華大學(xué) 核科學(xué)技術(shù)學(xué)院;2. 南華大學(xué) 先進(jìn)核能技術(shù)設(shè)計(jì)與安全教育部重點(diǎn)實(shí)驗(yàn)室：湖南衡陽(yáng)421001;3. 蘇州大學(xué) 放射醫(yī)學(xué)與輻射防護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 江蘇蘇州 215123)

惡性腫瘤現(xiàn)已成為導(dǎo)致人類死亡的頭號(hào)殺手,且發(fā)病率正在逐年增加。2022年,國(guó)家癌癥中心發(fā)布我國(guó)癌癥患者為406.4萬(wàn)人[1]。目前,癌癥的治療方式主要有手術(shù)治療、化學(xué)治療及放射治療。國(guó)內(nèi)將近50%～70%腫瘤患者在治療過(guò)程中接受放射治療。其中,硼中子俘獲治療(boron neutron capture therapy,BNCT)是一種嶄新的二元靶向放射治療手段,主要原理是通過(guò)將具有選擇性的含硼藥物注射入人體內(nèi),讓10B富集于腫瘤區(qū),用熱中子束照射10B富集區(qū)域,由于10B具有異常大的熱中子吸收截面,大量的10B會(huì)發(fā)生10B(n,α)7Li反應(yīng),釋放出α,7Li 2個(gè)具有高傳能線密度(linear energy transfer,LET)的粒子,這兩種粒子在約12～13 μm的范圍內(nèi)沉積能量,使沉積能量區(qū)域的細(xì)胞DNA產(chǎn)生難以修復(fù)的損傷,最終殺死腫瘤細(xì)胞,治療效果優(yōu)于常規(guī)治療手段[2]。

BNCT的治療包括計(jì)算和分析患者體內(nèi)的照射劑量分布,用以確定中子射束照射時(shí)間及方向角度的最優(yōu)化配置,并同時(shí)遵守危及正常組織和器官的劑量限制。BNCT中組織劑量基本來(lái)自于：(1)硼劑量。由10B(n,α)7Li反應(yīng)所獲得,且由于硼對(duì)熱中子的吸收截面非常大(3 840 b),釋放能量多,因而硼劑量對(duì)總劑量的貢獻(xiàn)很大,為主要?jiǎng)┝?(2)熱中子劑量。熱中子與人體中的14N原子發(fā)生核反應(yīng)生成反沖核14C和質(zhì)子,在人體細(xì)胞中沉積能量所產(chǎn)生的的劑量;(3)超熱中子及快中子劑量。統(tǒng)計(jì)的是除熱中子外所有中子由彈性散射而沉積的能量;(4)光子劑量。主要有伴隨入射中子束產(chǎn)生的γ射線、硼中子俘獲反應(yīng)中產(chǎn)生的射線及熱中子與人體內(nèi)的H原子發(fā)生俘獲反應(yīng)生成的γ射線3種來(lái)源。上述4種劑量單元的計(jì)算非常復(fù)雜,目前尚無(wú)精確劑量計(jì)算的相關(guān)經(jīng)驗(yàn)公式,現(xiàn)階段BNCT的治療計(jì)劃系統(tǒng)中主要采用蒙特卡羅模擬計(jì)算患者或體模的個(gè)體化模型劑量[3]。目前,國(guó)際上針對(duì)BNCT的治療計(jì)劃系統(tǒng)(treatment planning system,TPS)日漸成熟,開發(fā)了多種治療軟件,主要有：美國(guó)的NCTPlan和SERA、日本的JCDS、北京應(yīng)用物理與計(jì)算數(shù)學(xué)研究所研制的MCDB及中硼醫(yī)療公司研發(fā)的NeuMANTA等[4]。

本文基于成熟的蒙特卡羅方法研制了自主化BNCT蒙特卡羅劑量計(jì)算程序Magic。針對(duì)國(guó)際基準(zhǔn)題Snyder修正頭部模型,在此基礎(chǔ)上添加一個(gè)假設(shè)的腫瘤模型,建立1 mm精細(xì)體素模型來(lái)逼近解析模型,并在程序中采用網(wǎng)格中心點(diǎn)方法填充體素模型材料,同時(shí)對(duì)中子源、Kerma因子及S(α,β)熱散射等參數(shù)設(shè)置,使用MCNP和Magic程序進(jìn)行模型深度-劑量率曲線模擬對(duì)比,對(duì)結(jié)果進(jìn)行分析研究,初步驗(yàn)證了Magic程序正確性,并以Magic為基礎(chǔ),研究10B藥代動(dòng)力學(xué)模型下的硼濃度對(duì)BNCT劑量學(xué)影響規(guī)律,可為BNCT 劑量計(jì)算基準(zhǔn)化和臨床實(shí)踐應(yīng)用提供參考。

1 蒙特卡羅劑量計(jì)算程序Magic開發(fā)

目前,國(guó)際上BNCT治療計(jì)劃系統(tǒng)的劑量計(jì)算引擎主要使用MCNP作為核心計(jì)算傳輸工具,然而MCNP的禁止轉(zhuǎn)讓限制了其在中國(guó)的使用。劑量計(jì)算是整個(gè)BNCT治療計(jì)劃系統(tǒng)的核心部分,MCNP的限制將影響國(guó)內(nèi)自主化BNCT治療計(jì)劃系統(tǒng)的發(fā)展,同時(shí)MCNP程序體系龐大,主要用于核能、核工程及核技術(shù)等各領(lǐng)域的相關(guān)理論計(jì)算研究,并不適用于特殊研究需求的發(fā)展。鑒于此,本文基于蒙特卡羅粒子輸運(yùn)方法自主開發(fā)了BNCT治療計(jì)劃系統(tǒng)專用計(jì)算引擎—蒙特卡羅劑量計(jì)算程序Magic。Magic主要采用C++開發(fā),借助于C++豐富的標(biāo)準(zhǔn)庫(kù)和第三方庫(kù),覆蓋了多個(gè)可擴(kuò)展功能,便于未來(lái)的BNCT治療計(jì)劃系統(tǒng)的發(fā)展和維護(hù)。

Magic程序由7塊功能模塊組合構(gòu)成,包括幾何模塊(Geometry Module)、源項(xiàng)模塊(Source Module)、數(shù)據(jù)庫(kù)模塊(Database Module)、粒子輸運(yùn)模塊(Particle Transport Module)、計(jì)數(shù)模塊(Tallies Module)、輸出模塊(Output Module)和輔助模塊(Auxiliary Module),圖 1為Magic開發(fā)框架。幾何模塊主要采用面方程描述方法構(gòu)造復(fù)雜的劑量計(jì)算模型,基于組合幾何(CSG)方法進(jìn)行幾何描述,同時(shí)支持基于Lattice重復(fù)結(jié)構(gòu)的超精細(xì)毫米量級(jí)體素(voxel)建模,并對(duì)相應(yīng)組織區(qū)域從數(shù)據(jù)庫(kù)模塊中填充組織材料的核素成分及比例;源項(xiàng)模塊旨在設(shè)置粒子的初始條件和源的屬性,如源的空間分布、能量分布、角度分布及粒子類型等信息,并對(duì)相應(yīng)源參數(shù)支持多種概率分布抽樣(離散、連續(xù)、直方圖及混合),滿足模擬BNCT復(fù)雜中子源的要求;數(shù)據(jù)庫(kù)模塊儲(chǔ)存了基于ENDF/B-VIII.0評(píng)價(jià)數(shù)據(jù)庫(kù),利用具有層次結(jié)構(gòu),可處理大數(shù)據(jù)集等優(yōu)點(diǎn)的HDF5數(shù)據(jù)格式來(lái)存儲(chǔ)粒子輸運(yùn)中物理相關(guān)的各反應(yīng)截面數(shù)據(jù)[5],同時(shí)根據(jù)ICRU46文和ICRU63文內(nèi)容[6-7],創(chuàng)建了包括106種人體組織成分材料庫(kù)和相應(yīng)的Kerma因子庫(kù),并提供了支持修改和編輯相關(guān)數(shù)據(jù)的功能;粒子輸運(yùn)模塊主要處理程序模擬粒子從產(chǎn)生到消失的過(guò)程,中子與物質(zhì)相互作用主要涉及到彈性散射、非彈性散射及吸收,其中,吸收包含BNCT所關(guān)注劑量相關(guān)反應(yīng)(n,α),(n,p),(n,γ)等,在此過(guò)程中由中子產(chǎn)生的次級(jí)光子將被儲(chǔ)存在粒子庫(kù)中,為后續(xù)的光子輸運(yùn)做好準(zhǔn)備,光子與物質(zhì)相互作用主要考慮光電效應(yīng)、電子對(duì)效應(yīng)、相干散射及康普頓散射;粒子輸運(yùn)完成后,調(diào)用計(jì)數(shù)模塊,根據(jù)用戶的需求,設(shè)置不同的統(tǒng)計(jì)參數(shù),包括感興趣區(qū)域、能量區(qū)間和粒子類型等,以獲取相應(yīng)的注量,結(jié)合數(shù)據(jù)庫(kù)模塊中的Kerma因子庫(kù),得到劑量;最后通過(guò)輸出模塊對(duì)BNCT治療中關(guān)注的各種劑量,包括硼劑量、熱中子劑量、超熱中子劑量、快中子劑量、光子劑量及總相對(duì)生物效應(yīng)劑量進(jìn)行統(tǒng)計(jì),并提供可視化處理;輔助模塊為上述模塊提供相關(guān)函數(shù)庫(kù)及異常捕獲等。

Magic程序以粒子輸運(yùn)模塊為核心,與上游的幾何模塊、源項(xiàng)模塊及截面模塊相互連接,同時(shí)支持下游的計(jì)數(shù)模塊和輸出模塊。各模塊之間相互獨(dú)立,且模塊間通過(guò)接口進(jìn)行消息傳遞實(shí)現(xiàn)各模塊功能,利于程序維護(hù)和擴(kuò)展。

圖1 Magic開發(fā)框架

2 模型及方法描述

2.1 Snyder修正頭部模型

目前,國(guó)際上對(duì)于BNCT劑量理論計(jì)算普遍采用的基準(zhǔn)問(wèn)題為Snyder修正頭部模型,由3個(gè)橢球面來(lái)界定不同組織的邊界,將頭部由里到外分為腦組織、顱骨及頭皮組織3部分。整個(gè)頭部又被空氣包圍,通常稱為解析模型[3]。各分界面滿足的3個(gè)橢圓面解析方程可表示為

腦組織與顱骨邊界：

(1)

顱骨與頭皮組織邊界：

(2)

頭皮組織與空氣邊界：

(3)

但在實(shí)際情況中,由于患者腫瘤形狀、位置及大小各不相同,無(wú)法使用解析方程去描述模型,因此在國(guó)際臨床中,普遍使用連續(xù)均勻的立方體網(wǎng)格模型來(lái)逼近模型結(jié)構(gòu),通常稱為體素模型。目前Snyder修正體素模型網(wǎng)格尺寸通常取4,8,16 mm,當(dāng)構(gòu)建的體素網(wǎng)格越小,越能更好地趨近解析模型。故本文利用Magic構(gòu)建1 mm×1 mm×1 mm精細(xì)Snyder體素模型,共分224(x方向)×224(y方向)×224(z方向)=11 239 424個(gè)體素網(wǎng)格,驗(yàn)證正確性的同時(shí),也對(duì)Magic構(gòu)建千萬(wàn)網(wǎng)格數(shù)的承載力提出考驗(yàn)。模擬計(jì)算的同時(shí),在腦組織中加入半徑為1.5 cm的球面作為腫瘤邊界,模擬實(shí)際臨床中腦部膠質(zhì)腫瘤的情形。本文所模擬的含腫瘤Snyder修正模型如圖2所示。

含腫瘤Snyder修正模型中腦組織、顱骨及頭皮組織的元素成分和密度由Magic數(shù)據(jù)庫(kù)獲得,腫瘤詳細(xì)的組織成分參考文獻(xiàn)[8]。目前臨床試驗(yàn)要求,腫瘤細(xì)胞的10B濃度必須大于正常組織細(xì)胞的2.5倍,故本文設(shè)置頭皮組織及腦內(nèi)正常組織硼濃度為10-5,腫瘤組織內(nèi)的硼濃度為3×10-5,顱骨內(nèi)不含硼[9-10]。對(duì)解析模型中材料,可精確填充描述,但對(duì)網(wǎng)格模型,會(huì)因網(wǎng)格在兩種物質(zhì)交界面處難以處理,傳統(tǒng)處理方式通過(guò)得到每個(gè)立方體網(wǎng)格內(nèi)部不同材料精確的體積比填充,當(dāng)網(wǎng)格數(shù)量較多時(shí),材料填充極其耗時(shí)且會(huì)帶來(lái)一定偏差。為降低模擬問(wèn)題的復(fù)雜度,在Magic程序模型構(gòu)造中使用網(wǎng)格中心點(diǎn)的方法,根據(jù)每個(gè)網(wǎng)格的中心點(diǎn)位置確定材料密度,能大大降低材料數(shù)目,可由原來(lái)的286種材料降低至4種材料[11]。

(a)Tumor modified Snyder analytical model

2.2 相對(duì)生物學(xué)效應(yīng)劑量

由于不同種類的電離輻射會(huì)導(dǎo)致生物體產(chǎn)生不同的生物效應(yīng),BNCT涉及多個(gè)不同輻射劑量組分的混合場(chǎng)照射,因此,對(duì)BNCT治療須衡量每種劑量組分的相對(duì)生物學(xué)效應(yīng)劑量,評(píng)估其對(duì)正常組織細(xì)胞和腫瘤細(xì)胞的殺傷效應(yīng)。用實(shí)驗(yàn)測(cè)量的相對(duì)生物學(xué)效應(yīng)(relative biological effectiveness,RBE)來(lái)表征熱中子劑量、超熱中子及快中子劑量與光子劑量的生物學(xué)效應(yīng)高低。對(duì)于硼劑量的評(píng)估,須考慮α粒子和7Li離子的RBE及10B在組織和細(xì)胞中的微觀幾何位置對(duì)生物學(xué)效應(yīng)的綜合影響,被稱為復(fù)合生物學(xué)效應(yīng)(compound biological effectiveness,CBE)[12]。將上述4種劑量成分分別乘上各自對(duì)應(yīng)的加權(quán)因子得到總相對(duì)生物學(xué)效應(yīng)劑量H,表示為

H=WBDB+WγDγ+WnDn+WPDP

(4)

其中：WB為硼劑量的CBE因子;Wγ,Wn及WP分別為光子、快中子及氮俘獲產(chǎn)生熱中子劑量的RBE因子;DB,Dγ,Dn,DP分別為硼劑量、光子劑量、快中子劑量及氮俘獲產(chǎn)生的熱中子劑量。為使本文計(jì)算結(jié)果與國(guó)際上BNCT劑量結(jié)果具有可比性,選取腫瘤及健康組織不同劑量成分的RBE和CBE因子[13-16],如表1所列。

表1 腫瘤及健康組織不同劑量成分的RBE和CBE因子[13-16]

2.3 中子源設(shè)置

合適的中子源是BNCT治療腫瘤成功的關(guān)鍵要素之一,理想情況下,需選擇不同中子源來(lái)治療不同深度處的腫瘤。BNCT劑量計(jì)算中,要對(duì)入射源項(xiàng)粒子束的類型和形狀、粒子源的種類、粒子束能量及方向等參數(shù)進(jìn)行詳細(xì)描述。本文采取國(guó)際上針對(duì)Snyder修正模型所提出的寬譜超熱混合中子束,其中能量小于0.5 eV的熱中子占10%,0.5 eV～10 keV的超熱中子占89%,10 keV～2 MeV的快中子占1%,總體中子源能譜服從1/E(E為中子能量)分布,均勻分布在5 cm的圓盤面上,源強(qiáng)度為1010cm-2·s-1。

2.4 其他相關(guān)參數(shù)的設(shè)置

由于BNCT劑量計(jì)算時(shí)涉及中子-光子耦合輸運(yùn)過(guò)程,劑量可通過(guò)注量率與Kerma因子得到,表示為

(5)

其中：φ(j,E)為第j個(gè)體素網(wǎng)格,能量為E的歸一化注量率;V(j)為第j個(gè)體素網(wǎng)格的體積;K(j,E)為j個(gè)體素網(wǎng)格,能量為E的中子/光子劑量轉(zhuǎn)換因子,即Kerma因子。本文采用ICRU63給出的相對(duì)于ICRU46成年人腦的Kerma因子,但該Kerma因子所對(duì)應(yīng)的最低能量為0.025 3 eV,而對(duì)于BNCT的劑量來(lái)說(shuō),低于0.025 3 eV能量的中子對(duì)熱中子劑量的貢獻(xiàn)非常重要,故本文采取雙對(duì)數(shù)插值,將低于0.025 3 eV對(duì)應(yīng)的Kerma值外推到0.000 1 eV對(duì)應(yīng)的數(shù)據(jù),本文中光子的Kerma因子基于Seltzer[17]所計(jì)算的質(zhì)能吸收系數(shù)得到。

低能中子對(duì)BNCT的劑量貢獻(xiàn)很大,但當(dāng)中子能量降低到幾個(gè)電子伏時(shí),散射靶核的熱運(yùn)動(dòng)會(huì)對(duì)碰撞產(chǎn)生強(qiáng)烈的影響,從而對(duì)出射中子的能量及出射角度都會(huì)造成影響,最終會(huì)導(dǎo)致結(jié)果產(chǎn)生偏差。MCNP對(duì)于這種情況提供了S(α,β)熱散射模型,直接調(diào)用相應(yīng)核素的S(α,β)熱散射模型數(shù)據(jù)即可[18],本文Magic采取相同的S(α,β)熱散射模型處理Snyder材料中H的熱中子散射。

3 數(shù)值驗(yàn)證

3.1 MCNP與Magic計(jì)算結(jié)果比較

由于MCNP程序通常可作為各類蒙特卡羅計(jì)算的“標(biāo)準(zhǔn)”,常常用于標(biāo)定其他蒙特卡羅程序的準(zhǔn)確性測(cè)試,尤其是MCNP已成熟運(yùn)用在BNCT的治療計(jì)劃中。本文利用MCNP程序計(jì)算模擬Snyder解析模型的各劑量率作為標(biāo)準(zhǔn),使用Magic程序采用上文參數(shù)設(shè)置模擬同一解析模型,比較二者計(jì)算結(jié)果的差異,并在此基礎(chǔ)上繼續(xù)利用Magic程序?qū)? mm千萬(wàn)量級(jí)體素網(wǎng)格的Snyder模型進(jìn)行模擬,將模擬結(jié)果與MCNP解析模型結(jié)果進(jìn)行比較。本文重點(diǎn)在于探究Magic程序的正確性,對(duì)模型沿z軸統(tǒng)計(jì)相關(guān)劑量,其中每個(gè)統(tǒng)計(jì)點(diǎn)的計(jì)數(shù)網(wǎng)格大小設(shè)置為1.6 cm×1.6 cm×0.4 cm。在實(shí)際BNCT治療中有部分光子由中子源伴隨產(chǎn)生的,這部分光子是無(wú)法避免的,本文假設(shè)中子源無(wú)伴隨光子污染,只討論中子與體內(nèi)的核素發(fā)生反應(yīng)所產(chǎn)生的次級(jí)光子。

為減小蒙特卡羅程序的深穿透所帶來(lái)的結(jié)果影響,兩個(gè)蒙特卡羅程序輸運(yùn)模擬的粒子數(shù)均為1×108,確保樣本數(shù)量充足。圖3-圖7分別為硼劑量率、熱中子劑量率、超熱中子及快中子劑量率、次級(jí)光子劑量率和總的相對(duì)生物效應(yīng)劑量率及相對(duì)偏差隨深度的變化關(guān)系。相對(duì)偏差是以MCNP解析模型結(jié)果為標(biāo)準(zhǔn)值,Magic的解析/體素模型的結(jié)果與之比較得到的偏差百分比。

(a)Dose rate vs. depth

(b)Relative deviation vs. depth

(a)Dose rate vs. depth

(b)Relative deviation vs. depth

(a)Dose rate vs. depth

(b)Relative deviation vs. depth

(a)Dose rate vs. depth

(b)Relative deviation vs. depth

(a)Dose rate vs. depth

3.2 結(jié)果分析

由圖3(a)可見(jiàn),深度為0.5～1.5 cm,14～17 cm處為顱骨,由于顱骨不含硼,故沒(méi)有硼劑量率,這兩處的相對(duì)偏差為0;深度為4.5～7 cm時(shí),分布曲線上的陡然變化是由于該處為腫瘤區(qū)域,硼含量較高,導(dǎo)致硼劑量率顯著增加,有助于在臨床中更精確地確定腫瘤的位置。由圖3(b)可見(jiàn),相對(duì)偏差保持在BNCT臨床允許的5%以內(nèi)[19-20]。由圖3(a)、圖4(a)和圖6(a)可見(jiàn),分布曲線趨勢(shì)一致,這是由于次級(jí)光子來(lái)自于中子反應(yīng),中子劑量越多的地方產(chǎn)生的光子也越多,在頭部的淺表組織產(chǎn)生的損傷最大;由圖3(b)、圖4(b)和圖6(b)可見(jiàn),相對(duì)偏差也均保持在5%以內(nèi)。由圖5(a)可見(jiàn),因中子源發(fā)射大部分為超熱中子,一開始就發(fā)生中子慢化到熱中子,隨著深度的加深,劑量率逐漸降低,深度小于12 cm時(shí),相對(duì)偏差都基本保持在5%以內(nèi);深度大于12 cm時(shí),為深穿透問(wèn)題,蒙特卡羅程序計(jì)算時(shí)到達(dá)深處的粒子數(shù)較少,造成統(tǒng)計(jì)結(jié)果相對(duì)偏差較大。最后,將圖3-圖6的劑量率按照相應(yīng)的權(quán)重因子整合,得到總相對(duì)生物效應(yīng)劑量率,如圖7(a)所示。由圖7(b)可見(jiàn),相對(duì)偏差均保持在5%以內(nèi)?？傮w上,由圖3-圖7可見(jiàn),Magic的解析模型與體素模型的模擬計(jì)算結(jié)果和MCNP解析模型的結(jié)果符合較好,表明Magic的計(jì)算結(jié)果是正確的。

4 基于Magic的BNCT動(dòng)態(tài)硼濃度分析

中子注量率及硼濃度共同決定了BNCT沉積在組織中的劑量,因此在計(jì)算劑量時(shí),這二者對(duì)劑量評(píng)估的準(zhǔn)確性有很大的影響。目前在進(jìn)行理論BNCT的治療計(jì)劃評(píng)估時(shí),對(duì)腫瘤和正常組織的硼濃度設(shè)置都是時(shí)間固定不變的,但實(shí)際給藥后組織和血液中的硼濃度會(huì)因新陳代謝隨時(shí)間而變化,因此為臨床劑量評(píng)估的準(zhǔn)確性,有必要輸運(yùn)計(jì)算時(shí)建立10B的藥代動(dòng)力學(xué)模型。本文基于已驗(yàn)證正確性的蒙特卡羅劑量計(jì)算程序Magic定量研究動(dòng)態(tài)硼濃度下BNCT劑量學(xué)規(guī)律,為臨床應(yīng)用的安全性提供理論基礎(chǔ)。

目前,在BNCT含硼藥物BSH和BPA上,所使用的血硼濃度藥代動(dòng)力學(xué)連續(xù)模型為兩房室及三房室模型,而通常組織硼濃度是通過(guò)血硼濃度和靜態(tài)組織與血液硼濃度之比的乘積作為實(shí)際組織硼濃度的代替物來(lái)建模的[21]。本文模擬含腫瘤(膠質(zhì)腫瘤)Snyder模型,選取BPA-F介導(dǎo)的多形性膠質(zhì)母細(xì)胞瘤BNCT情況進(jìn)行分析研究,圖8為硼藥(BPA-F)靜脈注射后顱內(nèi)血液中10B的開放式二房室模型。

其中：中央室(central compartment)代表血液和注射分布良好的組織;周邊室(peripheral compartment)代表注射較差的組織,如大腦;I為藥物BPA-F注射的速率,mg·kg-1·min-1;C1,C2分別為各室的10B濃度;V1,V2分別為各室的分布容積;K12,K21,K10為速率常數(shù),min-1。

結(jié)合二室模型,并以哈佛-麻省理工學(xué)院I期臨床試驗(yàn)患者的藥代動(dòng)力學(xué)研究數(shù)據(jù)[22]做相應(yīng)推導(dǎo),得到中央室血硼濃度C1隨時(shí)間的變化關(guān)系,如圖9所示。根據(jù)圖9,在BPA-F介導(dǎo)的多形性膠質(zhì)母細(xì)胞瘤BNCT中,腫瘤與血液的硼濃度之比通常假定為3.5∶1,腦部與血液的硼濃度之比假定為1∶1,頭皮與血液的硼濃度之比假定為1.5∶1,從而利用血硼濃度變化曲線來(lái)預(yù)測(cè)其他組織中的硼濃度變化[22]。本文選取上述曲線硼濃度上升階段5個(gè)時(shí)間點(diǎn),下降階段10個(gè)時(shí)間點(diǎn),共15個(gè)時(shí)間點(diǎn),結(jié)合各組織倍數(shù)關(guān)系,利用Magic計(jì)算10B藥代動(dòng)力學(xué)模型下腫瘤的總相對(duì)生物效應(yīng)劑量率隨時(shí)間的變化關(guān)系,如圖10所示。

圖9 中央室血硼濃度隨時(shí)間的變化關(guān)系

圖10 10B藥代動(dòng)力學(xué)模型下腫瘤的總相對(duì)生物效應(yīng)劑量率隨時(shí)間的變化關(guān)系

圖10中陰影部分的時(shí)間段,可視為BNCT治療時(shí)持續(xù)照射的時(shí)間,紅色陰影部分為建立了10B的藥代動(dòng)力學(xué)模型,考慮動(dòng)態(tài)硼濃度變化后實(shí)際所受的劑量;藍(lán)色陰影處為目前理論BNCT穩(wěn)態(tài)下硼濃度所受的劑量。由圖10可直觀比較出二者的差異,對(duì)陰影面積(劑量值)進(jìn)行偏差分析,二者陰影區(qū)域相對(duì)偏差為11.78%。因此,在BNCT臨床治療過(guò)程中,若不考慮患者體內(nèi)硼濃度的動(dòng)態(tài)變化過(guò)程,則可能導(dǎo)致實(shí)際受照劑量與處方劑量存在較大偏差,影響治療效果。

4 結(jié)論

本文基于自主開發(fā)的BNCT蒙卡劑量計(jì)算程序Magic對(duì)國(guó)際Snyder修正模型進(jìn)行劑量計(jì)算分析,通過(guò)MCNP得到深度-劑量率結(jié)果為標(biāo)準(zhǔn)進(jìn)行對(duì)比,得出Magic程序計(jì)算的結(jié)果與參考文獻(xiàn)結(jié)果相符,硼劑量率、熱中子劑量率、超熱中子及快中子劑量率、次級(jí)光子劑量率和總的相對(duì)生物效應(yīng)劑量率的相對(duì)偏差均小于5%,滿足臨床治療對(duì)計(jì)算精度的要求,驗(yàn)證了Magic程序正確性,同時(shí)證明Magic程序具有構(gòu)建毫米量級(jí)千萬(wàn)網(wǎng)格數(shù)模型的承載力?；贛agic建立10B的開放式二房室藥代動(dòng)力學(xué)模型,對(duì)動(dòng)態(tài)硼濃度下BNCT劑量學(xué)進(jìn)行定量研究,根據(jù)血硼濃度變化曲線及各組織倍數(shù)關(guān)系,計(jì)算得出不同時(shí)間點(diǎn)下各硼濃度的腫瘤隨時(shí)間變化的總相對(duì)生物效應(yīng)劑量率變化擬合曲線。選取BNCT持續(xù)照射時(shí)間下,考慮動(dòng)態(tài)硼濃度變化后,實(shí)際的所受的真劑量值及穩(wěn)態(tài)下的硼濃度所受的假劑量值,計(jì)算得出二者的相對(duì)偏差為11.78%,可為臨床更加準(zhǔn)確的劑量計(jì)算及下一步完善含時(shí)蒙特卡羅劑量計(jì)算程序Magic提供參考。