^125Ⅰ微型源在水模体中能量分布的蒙特卡罗模拟.pdfVIP

第 36卷 第2期 西 南 师 范 大 学 学 报 (自然科学版) 2011年 4月 Vo1．36 No．2 JournalofSouthwestChinaNormalUniversity(NaturalScienceEdition) Apr． 2011 文章编号 ：1000—5471(2011)02—0052—06 125I微型源在水模体中能量分布的蒙特卡罗模拟 刘 漪 ， 王春燕。 1．北京工商大学机械工程学院，北京 100048；2．首都医科大学 生物医学工程学院，北京 100069 摘要：介绍了 射线辐射场能量分布 的蒙特卡罗模拟计算方法 ；以立方体水为介质模型，利用 EGS4软件模拟 了7 光子在三维水模体中的输运过程 ，并对微型源”I在水模体 中的辐射场能量分布进行了模拟计算 ；研究了 I源 的 排列方式、几何特性以及介质 的均匀性对辐射场分布 的影响． 关 键 词 ：蒙特卡罗方法 ；辐射场；模拟计算 中图分类号 ：R814．2 文献标志码 ：A 恶性肿瘤是危害人们生命的重大疾病之一，目前治疗恶性肿瘤 的最有效手段之一是放射治疗．放射治 疗的 目的是消灭大量的癌细胞，然而，临床医学表明，在放射过程中，伴随着癌细胞的消失 ，肿瘤周围的正 常组织不可避免地会受到一定程度的辐射剂量，从而增加了正常细胞癌变的机会．因此，靶 区(肿瘤)所受 剂量及其分布是治疗过程中必须考虑的问题 ，剂量太小会影响治疗效果 ，剂量过大又会对人体造成伤 害 ．如何提高治疗效果 ，最大限度地消除癌细胞而保持肿瘤周 围正常组织避免辐射 的伤害，已成为当前 国内外医学物理学研究的热点问题． 近几年来 ，以三维适形和调强适形放射为代表的治疗技术在临床医学中得到了广泛的应用 ，这些新技 术的特点是 ：通过对肿瘤病变部位的精确模拟定位、对进人人体辐射剂量分布的精确计算 以及正确选用放 射源 ，来进行治疗方案的设计与实施，以提高靶区(肿瘤)的照射剂量 ，减少肿瘤周围正常组织所受到的放 射损伤，达到最佳的治疗效果_2．在治疗方案的实施过程 中，关键的问题是如何精确计算辐射剂量 ，常规 的计算方法是在假设介质半无限大的条件下求解粒子在人体组织中输运方程的近似解 ，对于非均匀介质的 问题，通常采用半经验方法处理，由此导致的计算误差很大 ．随着计算机的诞生，一种新的计算方 法——蒙特卡罗 (MonteCarlo)模拟计算方法逐渐发展起来 ，蒙特卡罗计算方法是一种统计试验方法 ，它直 接从实际物理模型出发 ，采用统计抽样技巧模拟粒子在人体组织中的输运过程 ，通过沉积能量的方法 ，计 算出辐射场能量 (剂量)的分布．这种方法具有计算精确度高，和问题 的几何结构、维数无关等特点 ，得到 了国际社会普遍的认可 ． 1 蒙特卡罗模拟计算方法 蒙特卡罗方法是一种随机模拟方法 ，又称计算机模拟方法．随着计算技术 的迅速发展 ，蒙特卡罗方法 计算技术不断完善 ，应用领域不断扩大．与一般数值计算方法的最大区别是，蒙特卡罗计算方法能够逼真 收稿 日期 ：2010—04—01 基金项 目：“材料物理学术创新团队”北京市属高等学校人才强教计划资助项 目(PHR201007122)}北京市 自然科学基金资助项 目 (2102014)． 作者简介 ：刘 漪 (1958一)，女 ，上海人 ，副教授 ，主要从事大学物理本科教学及应用物理方 面的科研工作． 第2期 刘 漪 ，等 ： I微 型源在水模体中能量分布的蒙特卡罗模拟 53 地描述随机事物(如 7光子的输运过程)的特点及模拟真实的物理实验过程 ，用它来部分取代物理实验，不 仅可 以间接地得到 Boltzmann输运方程的解 ，甚至可 以得到物理实验难 以得到的结果 ，解决一般数值计算 方法难 以解决的问题 ．因此蒙特卡罗计算方法被广泛地应用于核物理、高能物理、生物、医学物理学等 领域 ，用 以解决模拟粒子的各类运输 问题． 1．1 蒙特卡罗统计样本的建立及统计方法 粒子在介质中的运输问题，涉及粒子与介质之间的相互作用过程 ，虽然它的反应机制复杂 ，但却遵循 着一定的统计规律．蒙特卡罗方法解决此类 问题 的基本思想 ：首先建立统计模型 、确定统计模型的参数 ， 其次将实际问题的模拟