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

蒙特卡罗方法在肿瘤放射物理学中的应用 班级：XXXXX 学号：XXXXXXXXXX 姓名：XXXX 【摘要】蒙特卡罗方法。蒙特卡罗方法（MonteCarlomethod），也称统计模拟方法，是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数（或更常见的伪随机数）来解决很多计算问题的方法。蒙特卡罗方法又称随机抽样技术或统计试验方法，是基于计算机模拟物理的思想，抓住物理过程的数量和几何特征，利用数值方法加以仿真，能够比较逼真地描述事物的特点及物理过程。该方法是求解辐射输运问题和粒子能量沉积的一种相当成熟、实用和有效的方法，所以在肿瘤放射物理学尤其是远距离放射治疗中的应用已越来越广泛。 【关键词】蒙特卡罗方法；辐射输运；能量沉积【正文】粒子与物质相互作用时服从统计学规律，发生作用的位置、作用的形式(如对光子而言，有光电效应、康普顿效应、电子对效应)，发生作用后粒子可能被吸收或散射，散射粒子的运动方向和能量、2次作用位置间的距离等参数均是随机变量。蒙特卡罗方法可以模拟粒子与物质相互作用的全过程，通过模拟10万甚至100万个粒子的输运过程，就可以比较精确地计算出粒子束与物质相互作用的宏观特征，如注量分布、吸收剂量分布。用蒙特卡罗方法解粒子输运问题一般包括3个过程：源分布抽样过程，产生粒子的初始状态；空间、能量和运动方向的随机游动过程，产生粒子的运动状态序列；记录贡献与分析结果，记录每个粒子对所求量的贡献并分析所求量的误差。蒙特卡罗方法的优点是可以处理粒子输运的各种复杂情况，尤其是一些难以进行实验测量的情况，下面就蒙特卡罗方法在肿瘤放射物理学中的应用综述如下。 1外照射源模拟以加速器产生的X射线的输运过程为例，蒙特卡罗技术是用随机抽样技术去模拟两个过程：由加速器X射线靶产生的X射线和一级准直器、均整器产生的散射X射线组成的初始射线和散射线的能谱及离轴分布；初始射线及散射线光子在介质中的输运过程。目前已有蒙特卡罗程序，如EGS4和EGSBEAM等，不仅能模拟均匀介质中上述两种过程，而且能够精确模拟不均匀介质中，光子和经一级和二级碰撞产生的次级电子的径迹，和计算它们的能量沉积，对一个治疗射野，蒙特卡罗模拟跟踪上亿个光子事件的能量沉积过程。在最常用的BEAM程序中，将粒子的输运过程分成几个区间，通过对上百万个粒子的初始状态、碰撞位置、碰撞性质以及碰撞后粒子的能量与方向的确定，得到每个区间或截面的宏观特征，并记录结果。例如单纯的光子输运过程，输运之前确定其初始状态（空间位置、能量的大小、运动方向），并根据光子与物质作用的宏观截面的大小，确定下次作用事件的位置和作用性质等，由光子与物质作用的形式，在进入下一个区间时确定是否继续跟踪此光子。BEAM程序提供加速器机头的组成模块，以文件形式输出，并提供一些降低方差技术，使用者可根据加速器的具体情况而灵活应用。2外照射时模体内辐射场模拟模拟粒子束进入模体能量沉积的过程，包括粒子通过加速器与模体表面之间的空气间隙的作用过程。最常用的软件是SPRXYZ和DOSXYZ，SPRXYZ是在SPRRZ和DOSXYZ的基础上产生的计算法则。Rogers等利用SPRXYZ对管电压标称值为10 mV和20 mV的X射线的水、空气阻止本领比（sw：air）作了计算，结果与Andreo和Nahum的计算结果一致。结果显示水、空气阻止本领比（sw：air）随深度的变化可达到1.25%，变化趋势是在模体表面附近略大，在最大电离深度处大幅度减小，然后随着深度的增加均匀减小。利用DOSXYZ除了可以准确给出模体内离轴比和百分深度剂量以及射野的大小外，还可以求解除剂量以外的一些物理量，比如能谱和初始射线、散射线的平均能量。Rogers等用蒙特卡罗技术对管电压标称值为10 mV和20 mV能级的X射线的中心轴百分深度剂量作了模拟计算，计算结果与测量结果相差不到最大深度剂量的1.5%，并且对TRS277中的有效测量点用蒙特卡罗方法作了验证。结果表明有效测量点的修正比水、空气阻止本领的变化对中心轴百分深度剂量的影响更大，而且有效测量点下移0.6 r(电离室空腔内径)，与0.75 r相比，前者的测量结果与蒙特卡罗计算结果更为接近，而后者的差异可达到1.5。 3外照射治疗计划应用由于蒙特卡罗方法运算速度慢，还不能直接用于治疗计划的制定，而是对计划中所采用的计算规则和技术作理论验证，以及为新的治疗计划提供配置数据。放射治疗的基本目标是努力提高放射治疗的增益比，即最大限度地将放射线的剂量集中到病变（靶区）内，杀灭肿瘤细胞，而减少周围正常组织和器官的损伤。因此理想的放射治疗技术应按照肿瘤形状给靶区很高的致死剂量，而靶区周围的正常组织不受到照射。适形调强放射治疗是一种提高治疗增益比的较为有效的物理措施，但由于缺乏治疗验证措施和计划评估手段，限制了适形调强放射治疗的广泛应用。建立在CT影像基础上的蒙特卡罗三维计算，可给出各种适形野以及由多野结合、楔形板、组织补偿技术等造成的不规则野的剂量分布，并提供理论验证或各种配置数据。4腔内放疗源周围辐射场模拟主要应用在近距离放射治疗中，如研究125I在人体内的剂量分布，利用蒙特卡罗方法计算放射性液体球囊在血管组织内的剂量分布等。随着计算机技术的迅猛发展，蒙特卡罗方法在肿瘤放射物理学中的应用将会越加广泛。当然，蒙特卡罗方法在解决粒子输运问题中依然存在收敛速度慢以及概率性质的误差等问题。除此之外，经验证明，只有当系统的大小与粒子的平均自由程可以相比较时，即系统大小为十多个平均自由程时，结果较为满意。而对于大系统，算出结果偏低，但是数值方法适应于大系统的计算，得到结果较好，现在已经有人将两者结合起来使用，取得了一定效果。放射治疗在肿瘤治疗中起到了非常重要的作用,为了保证放射治疗的质量和效果,必须验证放射治疗辐射场的剂量分布.胶片与其他探测器相比,具有使用方便,空间分辨率高,能记录完整的二维剂量分布.模体中而不用考虑侧向电子不平衡问题等优点.因此胶片是验证先进放射治疗技术(立体定向放射外科手术(SRS)、动态楔型板和调强适形放射治疗(IMRT)等)的剂量分布不可缺少的工具，是研究使用胶片测量放射治疗辐射场的剂量分布.以电子直线加速器、X刀、钴60远距离治疗机和伽玛刀放射治疗系统为研究对象,使用胶片(XV、ED（略）和扫描水箱(RFA-300和RTS-200)测量以上放疗设备的放疗辐射场剂量分布,分析比较剂量分布的两个重要参数(半高宽和半影)来评价（略）的一致性.结果分析表明,测量6MV电子直线加速器在标准条件下的辐射场剂量分布参数半高宽和半影,XV、EDR、EBT与RFA-300扫描（略）果的最大绝对偏差分别是0.5、0.3、0.2mm;测量钴60治疗机标准条件下辐射场剂量分布参数半高宽和半影,XV、EDR、EBT与RTS-200扫描水箱的测量结果的最大绝对偏差分别是2.6mm、2.8mm、3.0mm。【参考文献】 1 胡逸民.肿瘤放射物理学M.北京：原子能出版社，1999： 495496. 2 徐志勇，张良安，苑淑渝.32P液体球囊在血管组织中的剂量分布 J.中华