「细胞剂量模型计算」– 以蒙地卡罗模拟方法研究细胞微剂量学.docVIP

「細胞劑量模型計算」– 以蒙地卡羅模擬方法研究細胞微劑量學 by 劉晉昇醫用物理師 台北榮總核醫部與清華大學醫學物理劑量實驗室合作開發之細胞劑量計算程式，是以蒙地卡羅數值模擬方法，可以計算單能量電子、阿伐粒子及各種核種之細胞劑量S值與微劑量學物理量，分別與美國核醫學會發表細胞S值表及美國橡樹嶺國家實驗室(Oak Ridge National Laboratory)所發表之蒙地卡羅套裝軟體 PENELOPE Code 模擬結果比較。在本報告中醫用物理師將介紹下列內容： 蒙地卡羅方法基本概念(Basic Concepts of Monte Carlo Simulation) 「蒙地卡羅方法」？它是一種數值方法，利用亂數取樣 (random sampling) 來模擬解決數學問題。原子核衰變、游離粒子在物質中行進產生反應，可以描述成是一個機率問題，也就是所謂Stochastic Processes。蒙地卡羅方法被廣泛應用在放射科學中，如放射治療、放射診斷與核子醫學領域。 首先利用所有可能發生結果，定義出一機率密度函數(Probability Density Function，PDF)。再將此機率密度函數累加成累積機率函數(Cumulative Probability Function，CPF)，並乘以一常數調整其值為1，此稱為歸一化(Normalization)。這也將正確反應出所有事件出現的總機率為1的機率特性，這也為亂數取樣與實際問題模擬建立起連結。也就是說我們將電腦所產生均勻分佈於 [0, 1] 之間的亂數，透過所欲模擬的過程所具有機率分佈函數，模擬出實際問題最可能結果。 (有關基本原理可參閱92年10月核醫會訊) 醫學體內劑量(Medical Internal Radiation Dosimetry, MIRD) 美國核醫學會發展的醫學體內輻射劑量(MIRD)計算程式，其基本假設為藥物平均分佈於器官或組織內, 自射源器官到標靶器宮的劑量是均勻，其所得劑量屬於巨觀的、平均的吸收劑量。為了計算身體一標靶區域的吸收劑量，必須精確地測定人體組織中與時間有關放射性強度變化，並代入MIRD 方程式計算可得。為了解整體研究計劃架構，在本節中將簡略介紹 MIRD 體內劑量方程式 ，其中為平均吸收劑量、為累積放射性強度， 而 S則是每單位累積強度平均吸收劑量。 吸收劑量定義為每單位質量吸收之能量，這也是用來表示放射性藥物在人體中所造成輻射吸收劑量的基礎。放射強度為放射性核種每單位時間衰變次數，單位以居禮(Ci)或貝克(Bq)表示，而累積放射性強度則定義為射源器官區域內放射性強度隨時間變化曲線下面積，即總衰變次數，以微居禮-小時為單位。 而標靶器官所吸收劑量也可用攝入總放射性強度A0表示，此時射源器官放射性強度滯留時間定義為 ，此式可視為總放射性強度 A0 在 ? 滯留時間內釋放出與累積放射性強度相同總衰變次數，因此。此方程式才是真正MIRD臨床所使用最基本的公式。這式子中影響吸收劑量因素有兩個資料類型： 一為隨時間變化因子(生理動力學的)，例如累積放射性強度與滯留時間 ?。另一為與時間無關因子(物理性的)，例如單位累積強度平均吸收劑量 S。 其中隨時間變化因子綜合了感興趣區域(ROI, regions of interest)放射性強度的攝入與滯留特性，它也包含核種的物理半衰期與生物半衰期。與時間無關物理性因子 S 則反應出核種釋放出粒子束種類與能量、射源與標靶區域尺寸大小和形狀、與其間之幾何距離和吸收組織成份。 細胞S值(Cellular S values) 美國核醫學會所發展細胞劑量模型，假設細胞模、細胞核由兩均質同心圓構成，其密度為 1 g/cm3，細胞模半徑 RC、細胞核半徑 RN，典型細胞直徑大小從6至20微米，而細胞核直徑大小從4至18微米。放射性活度假設以下列四種方式均勻分佈在細胞各區域：整個細胞(throughout the cell, C)、細胞質中(cytoplasm, Cy)、細胞模表面(cell surface, CS)、細胞核中(cell nucleus, N)。 首先，我們必須將巨觀MIRD公式改寫成微觀的、並適合細胞使用之S值公式。再利用蒙地卡羅數值模擬方法，假設放射性核種均勻分佈在四種區域情形，以國際輻射防護委員會(ICRP) No.49報告中所提供阿伐粒子與貝它粒子之能量射程關係式，編寫Fortran程式計算沉積於標靶細胞區域能量，從而得出細胞 S值，以便往後計算細胞劑量所用。 微劑量學物理量Microdosimetric Quantities 微劑量學一般常用物理量分析包括：比能(specific energy)、平均弦長(mean chord length)、線能(lineal energy)三項