蒙特卡罗方法在放疗计划中的应用.pdfVIP

第 29 卷 第 11 期 核 技 术 Vol. 29, No.1 2006 年 1 月 NUCLEAR TECHNIQUES January 2006 蒙特卡罗方法在放疗计划中的应用 陈朝斌 黄群英 吴宜灿 林 辉 （中国科学院等离子体物理研究所 合肥 230031 ） 摘要：本文就蒙特卡罗方法在放射外科治疗计划中的应用做了初步的研究和分析。首先利用 Beam 程序模拟了 Varian 公司 Clinac 2100C 6MV X 射线加速器治疗头，得到源皮距为 100 cm、射野大小为 5 cm×5 cm 的开放野 治疗平面上的粒子相– 空间信息，并与文献进行了比较。然后把粒子相–空间信息作为输入，利用 Dosxyz 程序 计算了体模中吸收剂量的空间分布，并与一基于解析剂量计算方法的商用治疗计划软件的结果进行了比较， 两种结果基本一致。另外，本文还重点介绍了由患者的医学 CT 图像建立蒙特卡罗输运几何的方法和过程。 关键词：放射治疗计划，蒙特卡罗方法，CT 图像，剂量计算，精确放疗计划系统（ARTS ） 中图分类号：R730.55 文献标识码：A 文章编号：03-252J 人体组织剂量分布计算方法基本上可以分为两 目前规则靶区常规放疗仍占大多数，多叶准直器 大类，即半经验解析方法和蒙特卡罗方法。 （MLC ）在接近规则野情况下仍可使用半经验解析 半经验解析方法包括基于离轴比（OAR,Off 剂量计算方法。这些方法使用的绝大部分数据都是 Axis Ratio ）经验公式的方法，以及基于笔束核和点 针对具体治疗机的实测数据，靶区内任一点 核的卷积/叠加 （Convolution/Superposition ）方法等。 r (x , y ,z ) 处的吸收剂量D (x, y , z ) 可表达为： D (x , y , z ) Mu=⋅S (s ) ⋅TMR(d ,s ) ⋅OAR(s ,r ) ⋅F 2 ⋅F ⋅F ⋅F ⋅C ⋅C （1） t w SAD s p c 式中参数都是针对某一种能量下某一类型射线的， 其中， h 为点分布函数， T 为微分比释总能 都是能量的函数，s 是射野面积，d 为深度，r 为离 （Terma ）。 轴半径；Mu 是治疗机剂量监测器的计数（Monitor 半解析剂量计算算法的精确性和适用范围是有 Unit ）；St (s)是总的散射修正因子；TMR 为组织最 限的，基于测量数据的面核模型（即离轴比公式） 大比；OAR 为离轴比；F2 为反平方修正因子；Fw 只适用于开放野治疗情况，而不适用于调强、适形 为楔形因子；FSAD 为源轴距（Source Axis Distance） 放疗，而且必须修正组织不均匀性对剂量分布的影 因子；F 为防护托盘的透射因子；C 为非均匀组织 响。笔束核可以由测量的面核通过微分运算获得， s p 密度修正因子；Cc 为人体曲面修正因子。 或由蒙特卡罗粒子输运方法求解。它可以用于调强、 对于复杂射野剂量分布的计算，尤其是适形、 适形放疗剂量分布的计算，速度较快且适合射野强 调强放疗技术的应用，一般采用笔束核（Pencil 度分布的优化运算（用在逆向计划中），但其适用 Beam Kernel ）或点核（Point Kernel ）来进行组织剂 的前提是要求满足电子平衡条件，且对于直接使用 量分布的计算。在电子处于平衡状态的条件下，组