蒙特卡罗方法在放射外科治疗计划中的应用.docxVIP

蒙特卡罗方法在放射外科治疗计划中的应用 人体组织的剂量分布计算方法基本上可以分为两类：半经验分析法和蒙特卡罗法。 半经验解析方法包括基于离轴比 (OAR, Off Axis Ratio) 经验公式的方法, 以及基于笔束核和点核的卷积/叠加 (Convolution/Superposition) 方法等。 式中参数都是针对某一种能量下某一类型射线的, 都是能量的函数, s是射野面积, d为深度, r为离轴半径;Mu是治疗机剂量监测器的计数 (Monitor Unit) ;S 对于复杂射野剂量分布的计算, 尤其是适形、调强放疗技术的应用, 一般采用笔束核 (Pencil Beam Kernel) 或点核 (Point Kernel) 来进行组织剂量分布的计算。在电子处于平衡状态的条件下, 组织内任意点r的剂量D (r) 可以表示为能量注量分布和能量沉积核的卷积。实际上, 在组织密度变化较大的区域、并不满足电子平衡条件, 这时可以使用比释总能 (Terma) 和点分布函数 (Point Spread Function) 概念进行剂量计算, 点r处的剂量表示为: 其中, h为点分布函数, T为微分比释总能 (Terma) 。 半解析剂量计算算法的精确性和适用范围是有限的, 基于测量数据的面核模型 (即离轴比公式) 只适用于开放野治疗情况, 而不适用于调强、适形放疗, 而且必须修正组织不均匀性对剂量分布的影响。笔束核可以由测量的面核通过微分运算获得, 或由蒙特卡罗粒子输运方法求解。它可以用于调强、适形放疗剂量分布的计算, 速度较快且适合射野强度分布的优化运算 (用在逆向计划中) , 但其适用的前提是要求满足电子平衡条件, 且对于直接使用测量的数据进行剂量分布计算时, 因为实际很难得到微束 (细的窄束, 亦称为笔形束) 的准确放射测定数据, 因而进一步引入了剂量计算的不确定性 而蒙特卡罗方法, 因为其处理复杂问题 (复杂几何、复杂的放射源布置等) 的能力而成为一种不可替代的方法。蒙特卡罗方法可以精确地对放射治疗过程所涉及到的物理过程进行建模, 而使用较少的近似。蒙特卡罗方法的最大不足之处在于计算强度较大, 耗时较长。然而, 随着计算机处理速度的大幅提升和成本的持续降低以及在蒙特卡罗方法中引入新的方差减小技巧, 一些蒙特卡罗程序系统已经开始在临床中应用开来。例如, PEREGRINE程序就是在劳伦斯利弗莫尔国家实验室专为放射治疗而开发的多粒子蒙特卡罗程序, 可以作为治疗计划系统的剂量计算模块使用。Hartmann等 本文首先对Varian Clinac 2100C 6MV X射线治疗头进行了模拟, 得到治疗平面上粒子的相空间信息, 并把计算结果与文献计算结果进行了比较。然后利用我们自己开发的基于蒙特卡罗方法的精确放疗计划系统 (ARTS, the Advanced/Accurate Radio-Therapy System) 对一个头部实例进行了分析研究, 并与Varian公司基于常规剂量计算方法的放疗计划软件CADPlan6.2计算得到的结果做了比较分析, 两种结果基本一致。另外, 本文还重点介绍了由患者的医学CT图像建立蒙特卡罗输运几何的方法和过程。 1 在平台模拟中的应用 OMEGA项目是加拿大国家研究委员会 (NRCC) 和位于Madison的威斯康星大学 (Univercity of Wisconsin) 的合作项目, 其目的是开发基于EGSnrc程序系统的软件包使用户能快速地建立模拟直线加速器及患者模体的应用程序。其所开发的BEAM程序 蒙特卡罗粒子输运方法需要输入射野内粒子能谱分布、角分布和空间分布等有关治疗射线束特性的详细信息。按照获得这些信息的不同方式可将目前这方面的研究分为两大类:一类是用蒙特卡罗方法模拟大量粒子穿过治疗机头的输运过程, 获得粒子在模体表面的相空间数据, 直接利用相空间数据作为蒙特卡罗剂量算法的输入;另一类是通过对模拟得到的相空间数据进行分析, 建立治疗机源模型, 并以源模型为基础用蒙特卡罗方法计算模体或病人几何的剂量分布。本文使用BEAM程序对Varian Clinac 2100C (图1, 钨靶, 均整器为不锈钢) 加速器治疗头进行模拟的主要步骤如下:1.建立模拟加速器治疗头的用户程序和所需的光子–电子与介质相互作用的截面数据;2.完成用户的输入, 主要包括加速器治疗头各部件模块的几何定义, 入射粒子束定义以及程序运行控制参数的选择;3.进行模拟计算和分析计算的结果, 并把模拟计算结果作为体模吸收剂量计算的源项输入。图2是模拟Varian Clinac 2100C 6 MV X射线 (开放野, 源皮距为100 cm, 射野大小为5 cm×5 cm) 得到粒子能量注量的离轴分布。能量注量在射野内分布平