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

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

现在，医疗直线模拟器通常使用4.20mav能量的x线作为辐射源，并使用x射线梁的剂量分布不均匀。为了获得治疗领域的35.40cm大小，满足一定的平坦度和对称性要求，可以进行临床治疗。由于x射线穿过不同的材料，因此具有不同的衰减率。因此，产生光电效应、康普顿效应或电子效应的相互作用的概率是不同的，因此需要调整平均装置材料的选择和几何参数的配置是一个值得研究的问题。蒙特卡罗(MonteCarlo)方法非常适合于解决物质中粒子的输运问题,可以用来模拟粒子与物质的相互作用,被称作“理论上的实验”.由于MonteCarlo方法具有解析方法无法比拟的优点,因而在核科学、统计物理、分子动力学、医学等自然科学领域得到了广泛的应用.BEAMnrc是一种通用的MonteCarlo程序包,可用来模拟光子、电子等粒子输运问题.BEAMnrc全面考虑了电子、光子在通过物质时发生的各种相互作用,通过对大量粒子历史的追踪记录,得到治疗野平面的相空间数据.BEAMdp和DOSXYZnrc程序可以用来分析相空间数据,得到粒子在治疗野平面上的能量注量、能谱分布角分布以及模体中的离轴剂量分布、离轴比(off-axialratio,OAR)、百分深度剂量(percentagedepthdose,PDD)分布等信息.张洁熹等利用MonteCarlo程序BEAM对XHA9医用加速器治疗模式进行了模拟,计算了水模中的PDD分布曲线,分析了各部件产生的电子对治疗电子束的贡献,确定了限光筒产生的大量低能电子是导致出现标称能量低、表面剂量偏高的主要原因.通过修改限光筒、散射箔等部件,分步优化了治疗头结构.Sawkey等将测量所得的西门子加速器机头参数作为模拟参数输入,利用EGSnrc/BEAMnrc程序对6MV和18MV医用加速器进行了模拟,并将模拟计算得到的X射线剂量分布曲线与实验测量得到的剂量分布曲线进行了对比,二者的一致性非常好,说明MonteCarlo模拟计算结果非常准确,可将这种方法应用于临床.本研究利用BEAMnrc程序构建了加速器治疗头模型,采用尝试错误法设计了X射线均整器.首先对几种单一材料均整器下的射野剂量分布进行了模拟计算,通过分析各种材料对X射线衰减程度的特征,选择了两种金属作为均整器材料.反复调整几何参数进行计算,最终使得X射线在通过所设计的均整器后射野剂量分布均匀,满足了国家标准规定的均整度要求.此外,本研究还讨论了X射线剂量利用率、粒子角分布及带电粒子污染问题.1模型和方法1.1材料及几何模型的建立本研究所建立的加速器治疗头模型如图1所示.采用16MeV电子束打靶产生的X射线作为辐照源,电子束半径为0.001cm,垂直入射打靶,靶几何形状为长方体,长和宽均为1.6cm,厚为1cm,分3层,各层厚度由上而下依次为0.16,0.64和0.20cm,材料依次为金、水和不锈钢.模拟粒子数目为107个,取标称源皮距(sourceskindistance,SSD)为100cm,射野为10cm×10cm光子束照射野.初级准直器、次级准直器等部件按厂家给定参数设置.受治疗头内部空间条件的限制,均整器的几何外形给定为底面半径为3.0cm,高为3.2cm的圆锥(见图2),其他模拟参数按程序默认值设置,均整器的材料及其几何参数需要多次模拟后确定.1.2均整器材料参数的确定图3为模拟计算的流程图.治疗头的主要参数设定如1.1节所述,均整器参数的确定采取尝试错误法,即预先设定其材料和几何参数,在一次计算完成后对得到的数据进行分析,然后修改均整器的材料及几何参数再次计算.重复上述步骤,直至得到满意的结果为止.主要步骤如下:①选取铅作为均整器材料,其他参数如1.1节所述,进行仿真模拟;②分别选取铝、铜、钨作为均整器材料,其余参数同步骤①,分3次进行计算;③分析以上两步计算结果,设计复合材料均整器,几何外形如图4所示.采用两锥或多锥嵌套结构,即在上述固定外形的圆锥形均整器内部依次插入一个或多个小圆锥,内外圆锥的材料参数根据各材料对X射线的衰减程度和均整程度的不同而定.通过反复调整内部圆锥的锥度及底面半径,调节各材料的相对厚度,直到X射线剂量分布均匀且符合相关技术标准为止.2结果和分析2.1x-pb-2均整器设计及结果离轴剂量分布与离轴比是衡量射野平坦度和对称性好坏的重要标准.由图1可以看出,治疗头关于X,Y轴轮换对称,故只计算了X轴方向的离轴比及离轴剂量分布.图5为单一材料均整器对应的离轴比曲线.可以看出,低原子序数的铜、铝均整器对X射线的均整效果不明显;高原子序数的铅、钨对X射线有一定的均整度,其离轴比和剂量离轴分布曲线在距离射野中心轴附近有一定的平坦区域,对称性较好,但边缘处的相对剂量波动超过±5%,未能达到国际电子委员会(IEC)规定的±3%标准.图6为各均整器对应的离轴剂量分布曲线.可以看出,相对铝材料而言,重金属钨材料的均整器对应的X射线剂量利用率较低.考虑到低原子序数金属对X射线衰减程度较小,铜材料的均整器也产生了一定的均整效果,这里利用铝、铜两种金属制作了两层相嵌套的复合均整器,外锥几何尺寸如1.1节所述,材料为铝,内部小圆锥材料为铜,采用了多种不同尺寸分别进行计算.图6和图7记录了其中的两组曲线(标记为Al-Cu-1和Al-Cu-2).从这两组离轴剂量分布图和离轴比可以看出,只采用轻金属作为均整器材料很难将其离轴比和离轴剂量分布曲线中的凸起部分“消平”,要想实现X射线在治疗射野内的剂量平坦度要求,必须在射野中心轴附近适当加入高原子序数材料,通过重金属对X射线的衰减来降低射野中心轴附近的剂量值.采用重材料铅锥嵌入铝锥中心处(见图4),外锥尺寸如1.1节所述,内部铅锥底面半径为1.6cm,高为2.6cm.由于铅材料可使透射X射线能量有较大幅度的降低,而包围铅锥的铝材料对X射线的衰减程度较小,因此,这种设计不仅可以“消去”各离轴分布曲线中心处的凸起部分,还不会过分降低射野边缘区域的X射线剂量,进而在满足射野剂量平坦度要求的同时,达到提高X射线剂量利用率的目的.由图6和图7中Al-Pb-1的离轴剂量分布和离轴比可以看出,射野中心轴附近有较宽的平坦部分,说明铝-铅复合材料均整器的均整效果已经得到了显著提高,但射野边缘剂量值的陡降程度不够理想,导致均整区域内最大吸收剂量点与最小吸收剂量点吸收剂量的比值超出了国标中规定的不大于106%的标准.经过几次调整铅锥高度及底面半径大小,当铅锥底面半径调整为1.2cm,高为2.2cm时,得到了标记为Al-Pb-2的离轴分布曲线和离轴比,如图6和图7所示.图中可见,整个均整区域内,Al-Pb-2均整器下的X射线剂量分布平坦,且对称性好,区域边缘处剂量值急剧下降,均整区域内最大吸收剂量点与最小吸收剂量点吸收剂量的比值约为103%,完全满足GB15213—94的要求.对比图6中的各组离轴剂量分布曲线还可以看出,Al-Pb-2均整器下X射线剂量利用率较高.为验证模体不同深度处的剂量分布是否理想,计算了模体表面下5和15cm深度处X轴方向的离轴比和离轴剂量分布,结果如图8和图9所示.该结果与10cm深度处的结果具有一致性.5cm深度处X射线剂量分布均匀,治疗射野内剂量平坦度和对称性非常好,边缘区域剂量分布也满足要求;15cm深度处的离轴比和剂量离轴分布曲线在均整区域内出现了小幅度波动,这主要是由于统计误差造成的,更精确的计算需要继续增大入射粒子数.2.2u3000能量注量分布能量注量(energyfluence)是以进入辐射场内某点处单位截面积球体的粒子总动能描述辐射场性质的一个量.能量注量离轴分布可以从能量角度反映出治疗射野剂量平坦度和对称性问题.图10为源皮距SSD为100cm的治疗野平面上各均整器对应的能量注量离轴分布曲线.可以看出,X射线经各均整器后得到的能量注量分布差别很大,其规律与离轴剂量分布类似:低原子序数的铝、铜材料均整器下,能量注量在靠近射野中心轴附近较高,射野边缘较低;高原子序数的铅、钨材料均整器下,能量注量在射野内有较宽的平坦宽度,但射野边缘处剧降效果不明显,而且整体上能量注量值相对较低,反映出X射线剂量利用率较低.复合材料均整器(Al-Pb-2)对应的能量注量分布在射野内非常平整,均整区域内射线能量分布均匀,且能量注量值较高,说明X射线剂量利用率较高,边缘部分能量注量值急剧下降.2.3辐射防护.图11为X射线经过复合材料均整器(Al-Pb-2)均整后的粒子角分布曲线,插图为放大后的电子和正电子角分布曲线.由图可见,光子角分布集中在0°～4°之间,前向性较好,利于临床治疗使用.正负电子角分布范围广,角度较大,这种分布会使其偏离出治疗射野外,不能进入模体,不会造成射野内带电粒子污染.射野外的带电子污染可通过辐射防护设计来解决.2.4电子和正电子能谱分布图12为X射线经过复合材料均整器(Al-Pb-2)后在治疗野平面上的粒子能谱分布图,插图为放大后的电子和正电子能谱分布图.由图可见,光子能量在0～16MeV范围内呈连续分布,峰值出现在1MeV附近,射野内光子平均能量约为3.81MeV,电子和正电子能谱分布在数值上与光子能谱分布相差近3个数量级,说明带电粒子污染非常小.2.5x-ct/sdsd图13为X射线经复合材料均整器(Al-Pb-2)均整后得到的水模体内深度分布曲线,SSD为100cm.可以看出明显的X射线剂量建成区,最大剂量点约在模体下3.5cm处,表面剂量小于37%,满足国标中规定的中心轴表面剂量不得超过最大剂量的60%的要求.由于最大剂量点不在表面,因此,在治疗深度肿瘤时可以很好地保护皮肤.3模拟仿真模拟本研究利用BEAMnrc程序设计了复合材料均整器.由以上分析可以看出,重金属材料对X射线的衰减程度较大,因此,在经过了几次对计算结果分析后,选择了两种金属相结合作为均整材料.采用尝试错误法多次调节材料各区域的相对厚度,即根据X射线束流分布特点,增大射野中心轴附近区域重金属材料的厚度,减小边缘区域重金属材料的厚度,最终使得X射线分布均匀,在均整区域内得到了平坦的剂量