蒙特卡罗方法在肿瘤放射物理学中的应用.pdf

蒙特卡罗方法在肿瘤放射物理学中的应用 蒙特卡罗方法 又称随机抽样技术或统计试验方法， 是基于计算机模拟物理的 思想， 抓住物理过程的数量和几何特征， 利用数值方法加以仿真， 能够比较逼真地描述事物 的特点及物理过程。 该方法是求解辐射输运问题和粒子能量沉积的一种相当成熟、 实用和有 效的方法，所以在肿瘤放射物理学尤其是远距离放射治疗中的应用已越来越广泛。 粒子与物质相互作用时服从统计学规律， 发生作用的位置、 作用的形式 (如对光子而言， 有光电效应、康普顿效应、电子对效应 )，发生作用后粒子可能被吸收或散射，散射粒子的 运动方向和能量、 两次作用位置间的距离等参数均是随机变量。 蒙特卡罗方法可以模拟粒子 与物质相互作用的全过程，通过模拟 10 万甚至 100 万个粒子的输运过程，就可以比较精确 地计算出粒子束与物质相互作用的宏观特征， 如注量分布、 吸收剂量分布。 用蒙特卡罗方法 解粒子输运问题一般包括 3 个过程： (1)源分布抽样过程，产生粒子的初始状态； (2)空间、 能量和运动方向的随机游动过程，产生粒子的运动状态序列； (3)记录贡献与分析结果，记 录每个粒子对所求量的贡献并分析所求量的误差 [1] 。 蒙特卡罗方法的优点是可以处理粒子输运的各种复杂情况，尤其是一些难以进行实验 测量的情况。下面就蒙特卡罗方法在肿瘤放射物理学中的应用综述如下。 1 外照射源模拟 以加速器产生的 X 射线的输运过程为例，蒙特卡罗技术是用随机抽样技术去模拟两个 过程： (1)由加速器 X 射线靶产生的 X 射线和一级准直器、均整器产生的散射 X 射线组成的 初始射线和散射线的能谱及离轴分布； (2)初始射线及散射线光子在介质中的输运过程。 目前已有蒙特卡罗程序，如 EGS4和 EGS-BEAM等，不仅能模拟均匀介质中上述两种 过程，而且能够精确模拟不均匀介质中，光子和经一级和二级碰撞产生的次级电子的径迹， 和计算它们的能量沉积， 对一个治疗射野， 蒙特卡罗模拟跟踪上亿个光子事件的能量沉积过 程。 在最常用的 BEAM 程序中，将粒子的输运过程分成几个区间，通过对上百万个粒子的 初始状态、 碰撞位置、 碰撞性质以及碰撞后粒子的能量与方向的确定， 得到每个区间或截面 的宏观特征， 并记录结果。 例如单纯的光子输运过程， 输运之前确定其初始状态 (空间位置、 能量的大小、运动方向 )，并根据光子与物质作用的宏观截面的大小，确定下次作用事件的 位置和作用性质等， 由光子与物质作用的形式， 在进入下一个区间时确定是否继续跟踪此光 子。 BEAM 程序提供加速器机头的组成模块，以文件形式输出，并提供一些降低方差技术， 使用者可根据加速器的具体情况而灵活应用 [2] 。 2 外照射时模体内辐射场模拟 模拟粒子束进入模体能量沉积的过程，包括粒子通过加速器与模体表面之间的空气 间隙的作用过程。 最常用的软件是 SPRXYZ和 DOSXYZ，SPRXYZ是在 SPRRZ和 DOSXYZ的基 础上产生的计算法则。 Rogers 等利用 SPRXYZ对管电压标称值为 10MV 和 20MV 的 X 射线 的水、空气阻止本领比 (sw ，，air) 作了计算，结果与 Andreo 和 Nahum 的计算结果一致。 结果显示水、 空气阻止本领比 (sw ，，air) 随深