[發(fā)明專利]一種解析蒙特卡羅劑量計(jì)算方法

無效

• 申請?zhí)枺?/td>: 200910116115.6
• 申請日: 2009-01-22
• 公開（公告）號: CN101477203A
• 公開（公告）日: 2009-07-08
• 發(fā)明（設(shè)計(jì)）人: 鄭華慶;宋鋼;李貴;孟耀;蘭海洋;吳宜燦
• 申請（專利權(quán)）人: 中國科學(xué)技術(shù)大學(xué)
• 主分類號: G01T1/10
• 分類號: G01T1/10
• 代理公司: 安徽合肥華信知識產(chǎn)權(quán)代理有限公司
• 代理人: 余成俊

1.一種解析蒙特卡羅劑量計(jì)算方法，其特征在于包括以下步驟：

(1)獲得計(jì)算參數(shù)，包括以下內(nèi)容：

蒙特卡羅程序獲得的能量沉積核h：h是一系列單能或一定能譜的光子入射到某一均勻 介質(zhì)中，對應(yīng)于球坐標(biāo)(θ，R)的輻射能量值矩陣；

一系列單能或一定能譜的光子在某一介質(zhì)中的質(zhì)量衰減系數(shù)μ/ρ；

一系列單能或一定能譜的電子在某一介質(zhì)中的射程；

計(jì)算模型柵元劃分信息；

計(jì)算模型柵元物理密度信息；

計(jì)算模型柵元電子密度信息；

計(jì)算模型中計(jì)算區(qū)域邊界標(biāo)示信息；

計(jì)算模型中感興趣區(qū)域標(biāo)示文件；

獲得放射源的能譜信息；

射線在介質(zhì)表面形成的形狀，由用戶給出；

(2)利用下列MDSA模型進(jìn)行計(jì)算：

使用MDSA模型計(jì)算公式，如下：

D(r)=∫E∫∫∫VTE(r′)ρ(r′)ρ0c2h(E,c(r,r′)(r-r′))d3r′dE---(1)]]>

其中，T_E(r′)為微分比釋總能；ρ(r′)是作用柵元的物理密度；h為預(yù)先計(jì)算好的物理密 度為ρ₀的均勻介質(zhì)中的能量沉積核；c_(r，r′)代表了三維非均勻修正系數(shù)，是以某個作用柵元為 圓心，在垂直于作用柵元到沉積柵元連線方向，以帶電粒子在介質(zhì)中以射程為半徑的圓面為 底，以作用柵元到沉積柵元連線為高，所得的圓柱體內(nèi)的柵元電子密度的平均值：

c=c(r,r′)=∫01ρrel[r′-l(r′-r)]dl---(2)]]>

T_E(r′)使用子束入射到介質(zhì)表面的能量注量計(jì)算所有深度的微分比釋總能，即根據(jù)各個 小子束在介質(zhì)表面的劃分，隨著深度的增加，每個小野的尺寸按比例擴(kuò)大，保持小野數(shù)目不 變，相對于傳統(tǒng)DSA方法，減小各深度計(jì)算的采樣點(diǎn)數(shù)目，達(dá)到減小計(jì)算時間的目的；

在使用MDSA模型計(jì)算時，根據(jù)不同的計(jì)算區(qū)域，分別采用不同的計(jì)算點(diǎn)密度：感興 趣區(qū)域的計(jì)算點(diǎn)密度高，其他區(qū)域計(jì)算點(diǎn)密度低。