基于beamnrc程序的均整器射野劑量分布研究

基于beamnrc程序的均整器射野劑量分布研究

現(xiàn)在，醫(yī)療直線模擬器通常使用4.20mav能量的x線作為輻射源，并使用x射線梁的劑量分布不均勻。為了獲得治療領域的35.40cm大小，滿足一定的平坦度和對稱性要求，可以進行臨床治療。由于x射線穿過不同的材料，因此具有不同的衰減率。因此，產(chǎn)生光電效應、康普頓效應或電子效應的相互作用的概率是不同的，因此需要調整平均裝置材料的選擇和幾何參數(shù)的配置是一個值得研究的問題。蒙特卡羅(MonteCarlo)方法非常適合于解決物質中粒子的輸運問題,可以用來模擬粒子與物質的相互作用,被稱作“理論上的實驗”.由于MonteCarlo方法具有解析方法無法比擬的優(yōu)點,因而在核科學、統(tǒng)計物理、分子動力學、醫(yī)學等自然科學領域得到了廣泛的應用.BEAMnrc是一種通用的MonteCarlo程序包,可用來模擬光子、電子等粒子輸運問題.BEAMnrc全面考慮了電子、光子在通過物質時發(fā)生的各種相互作用,通過對大量粒子歷史的追蹤記錄,得到治療野平面的相空間數(shù)據(jù).BEAMdp和DOSXYZnrc程序可以用來分析相空間數(shù)據(jù),得到粒子在治療野平面上的能量注量、能譜分布角分布以及模體中的離軸劑量分布、離軸比(off-axialratio,OAR)、百分深度劑量(percentagedepthdose,PDD)分布等信息.張潔熹等利用MonteCarlo程序BEAM對XHA9醫(yī)用加速器治療模式進行了模擬,計算了水模中的PDD分布曲線,分析了各部件產(chǎn)生的電子對治療電子束的貢獻,確定了限光筒產(chǎn)生的大量低能電子是導致出現(xiàn)標稱能量低、表面劑量偏高的主要原因.通過修改限光筒、散射箔等部件,分步優(yōu)化了治療頭結構.Sawkey等將測量所得的西門子加速器機頭參數(shù)作為模擬參數(shù)輸入,利用EGSnrc/BEAMnrc程序對6MV和18MV醫(yī)用加速器進行了模擬,并將模擬計算得到的X射線劑量分布曲線與實驗測量得到的劑量分布曲線進行了對比,二者的一致性非常好,說明MonteCarlo模擬計算結果非常準確,可將這種方法應用于臨床.本研究利用BEAMnrc程序構建了加速器治療頭模型,采用嘗試錯誤法設計了X射線均整器.首先對幾種單一材料均整器下的射野劑量分布進行了模擬計算,通過分析各種材料對X射線衰減程度的特征,選擇了兩種金屬作為均整器材料.反復調整幾何參數(shù)進行計算,最終使得X射線在通過所設計的均整器后射野劑量分布均勻,滿足了國家標準規(guī)定的均整度要求.此外,本研究還討論了X射線劑量利用率、粒子角分布及帶電粒子污染問題.1模型和方法1.1材料及幾何模型的建立本研究所建立的加速器治療頭模型如圖1所示.采用16MeV電子束打靶產(chǎn)生的X射線作為輻照源,電子束半徑為0.001cm,垂直入射打靶,靶幾何形狀為長方體,長和寬均為1.6cm,厚為1cm,分3層,各層厚度由上而下依次為0.16,0.64和0.20cm,材料依次為金、水和不銹鋼.模擬粒子數(shù)目為107個,取標稱源皮距(sourceskindistance,SSD)為100cm,射野為10cm×10cm光子束照射野.初級準直器、次級準直器等部件按廠家給定參數(shù)設置.受治療頭內部空間條件的限制,均整器的幾何外形給定為底面半徑為3.0cm,高為3.2cm的圓錐(見圖2),其他模擬參數(shù)按程序默認值設置,均整器的材料及其幾何參數(shù)需要多次模擬后確定.1.2均整器材料參數(shù)的確定圖3為模擬計算的流程圖.治療頭的主要參數(shù)設定如1.1節(jié)所述,均整器參數(shù)的確定采取嘗試錯誤法,即預先設定其材料和幾何參數(shù),在一次計算完成后對得到的數(shù)據(jù)進行分析,然后修改均整器的材料及幾何參數(shù)再次計算.重復上述步驟,直至得到滿意的結果為止.主要步驟如下:①選取鉛作為均整器材料,其他參數(shù)如1.1節(jié)所述,進行仿真模擬;②分別選取鋁、銅、鎢作為均整器材料,其余參數(shù)同步驟①,分3次進行計算;③分析以上兩步計算結果,設計復合材料均整器,幾何外形如圖4所示.采用兩錐或多錐嵌套結構,即在上述固定外形的圓錐形均整器內部依次插入一個或多個小圓錐,內外圓錐的材料參數(shù)根據(jù)各材料對X射線的衰減程度和均整程度的不同而定.通過反復調整內部圓錐的錐度及底面半徑,調節(jié)各材料的相對厚度,直到X射線劑量分布均勻且符合相關技術標準為止.2結果和分析2.1x-pb-2均整器設計及結果離軸劑量分布與離軸比是衡量射野平坦度和對稱性好壞的重要標準.由圖1可以看出,治療頭關于X,Y軸輪換對稱,故只計算了X軸方向的離軸比及離軸劑量分布.圖5為單一材料均整器對應的離軸比曲線.可以看出,低原子序數(shù)的銅、鋁均整器對X射線的均整效果不明顯;高原子序數(shù)的鉛、鎢對X射線有一定的均整度,其離軸比和劑量離軸分布曲線在距離射野中心軸附近有一定的平坦區(qū)域,對稱性較好,但邊緣處的相對劑量波動超過±5%,未能達到國際電子委員會(IEC)規(guī)定的±3%標準.圖6為各均整器對應的離軸劑量分布曲線.可以看出,相對鋁材料而言,重金屬鎢材料的均整器對應的X射線劑量利用率較低.考慮到低原子序數(shù)金屬對X射線衰減程度較小,銅材料的均整器也產(chǎn)生了一定的均整效果,這里利用鋁、銅兩種金屬制作了兩層相嵌套的復合均整器,外錐幾何尺寸如1.1節(jié)所述,材料為鋁,內部小圓錐材料為銅,采用了多種不同尺寸分別進行計算.圖6和圖7記錄了其中的兩組曲線(標記為Al-Cu-1和Al-Cu-2).從這兩組離軸劑量分布圖和離軸比可以看出,只采用輕金屬作為均整器材料很難將其離軸比和離軸劑量分布曲線中的凸起部分“消平”,要想實現(xiàn)X射線在治療射野內的劑量平坦度要求,必須在射野中心軸附近適當加入高原子序數(shù)材料,通過重金屬對X射線的衰減來降低射野中心軸附近的劑量值.采用重材料鉛錐嵌入鋁錐中心處(見圖4),外錐尺寸如1.1節(jié)所述,內部鉛錐底面半徑為1.6cm,高為2.6cm.由于鉛材料可使透射X射線能量有較大幅度的降低,而包圍鉛錐的鋁材料對X射線的衰減程度較小,因此,這種設計不僅可以“消去”各離軸分布曲線中心處的凸起部分,還不會過分降低射野邊緣區(qū)域的X射線劑量,進而在滿足射野劑量平坦度要求的同時,達到提高X射線劑量利用率的目的.由圖6和圖7中Al-Pb-1的離軸劑量分布和離軸比可以看出,射野中心軸附近有較寬的平坦部分,說明鋁-鉛復合材料均整器的均整效果已經(jīng)得到了顯著提高,但射野邊緣劑量值的陡降程度不夠理想,導致均整區(qū)域內最大吸收劑量點與最小吸收劑量點吸收劑量的比值超出了國標中規(guī)定的不大于106%的標準.經(jīng)過幾次調整鉛錐高度及底面半徑大小,當鉛錐底面半徑調整為1.2cm,高為2.2cm時,得到了標記為Al-Pb-2的離軸分布曲線和離軸比,如圖6和圖7所示.圖中可見,整個均整區(qū)域內,Al-Pb-2均整器下的X射線劑量分布平坦,且對稱性好,區(qū)域邊緣處劑量值急劇下降,均整區(qū)域內最大吸收劑量點與最小吸收劑量點吸收劑量的比值約為103%,完全滿足GB15213—94的要求.對比圖6中的各組離軸劑量分布曲線還可以看出,Al-Pb-2均整器下X射線劑量利用率較高.為驗證模體不同深度處的劑量分布是否理想,計算了模體表面下5和15cm深度處X軸方向的離軸比和離軸劑量分布,結果如圖8和圖9所示.該結果與10cm深度處的結果具有一致性.5cm深度處X射線劑量分布均勻,治療射野內劑量平坦度和對稱性非常好,邊緣區(qū)域劑量分布也滿足要求;15cm深度處的離軸比和劑量離軸分布曲線在均整區(qū)域內出現(xiàn)了小幅度波動,這主要是由于統(tǒng)計誤差造成的,更精確的計算需要繼續(xù)增大入射粒子數(shù).2.2u3000能量注量分布能量注量(energyfluence)是以進入輻射場內某點處單位截面積球體的粒子總動能描述輻射場性質的一個量.能量注量離軸分布可以從能量角度反映出治療射野劑量平坦度和對稱性問題.圖10為源皮距SSD為100cm的治療野平面上各均整器對應的能量注量離軸分布曲線.可以看出,X射線經(jīng)各均整器后得到的能量注量分布差別很大,其規(guī)律與離軸劑量分布類似:低原子序數(shù)的鋁、銅材料均整器下,能量注量在靠近射野中心軸附近較高,射野邊緣較低;高原子序數(shù)的鉛、鎢材料均整器下,能量注量在射野內有較寬的平坦寬度,但射野邊緣處劇降效果不明顯,而且整體上能量注量值相對較低,反映出X射線劑量利用率較低.復合材料均整器(Al-Pb-2)對應的能量注量分布在射野內非常平整,均整區(qū)域內射線能量分布均勻,且能量注量值較高,說明X射線劑量利用率較高,邊緣部分能量注量值急劇下降.2.3輻射防護.圖11為X射線經(jīng)過復合材料均整器(Al-Pb-2)均整后的粒子角分布曲線,插圖為放大后的電子和正電子角分布曲線.由圖可見,光子角分布集中在0°～4°之間,前向性較好,利于臨床治療使用.正負電子角分布范圍廣,角度較大,這種分布會使其偏離出治療射野外,不能進入模體,不會造成射野內帶電粒子污染.射野外的帶電子污染可通過輻射防護設計來解決.2.4電子和正電子能譜分布圖12為X射線經(jīng)過復合材料均整器(Al-Pb-2)后在治療野平面上的粒子能譜分布圖,插圖為放大后的電子和正電子能譜分布圖.由圖可見,光子能量在0～16MeV范圍內呈連續(xù)分布,峰值出現(xiàn)在1MeV附近,射野內光子平均能量約為3.81MeV,電子和正電子能譜分布在數(shù)值上與光子能譜分布相差近3個數(shù)量級,說明帶電粒子污染非常小.2.5x-ct/sdsd圖13為X射線經(jīng)復合材料均整器(Al-Pb-2)均整后得到的水模體內深度分布曲線,SSD為100cm.可以看出明顯的X射線劑量建成區(qū),最大劑量點約在模體下3.5cm處,表面劑量小于37%,滿足國標中規(guī)定的中心軸表面劑量不得超過最大劑量的60%的要求.由于最大劑量點不在表面,因此,在治療深度腫瘤時可以很好地保護皮膚.3模擬仿真模擬本研究利用BEAMnrc程序設計了復合材料均整器.由以上分析可以看出,重金屬材料對X射線的衰減程度較大,因此,在經(jīng)過了幾次對計算結果分析后,選擇了兩種金屬相結合作為均整材料.采用嘗試錯誤法多次調節(jié)材料各區(qū)域的相對厚度,即根據(jù)X射線束流分布特點,增大射野中心軸附近區(qū)域重金屬材料的厚度,減小邊緣區(qū)域重金屬材料的厚度,最終使得X射線分布均勻,在均整區(qū)域內得到了平坦的劑量