基于混合笔束模型的三维电子剂量计算研究-核技术及应用专业论文.docxVIP

基于混合笔束模型的三维电子剂量计算研究核技术及应用专业 基于混合笔束模型的三维电子剂量计算研究 核技术及应用专业 研究生勾成俊 指导教师罗正明 在人类藉化学、免疫、抗生素疗法三大防治手段，基本上控制了以往主要 威胁人类健康的传染病之后，诸如肿瘤、心脑血管等非传染性疾病逐渐成为人 类健康的主要敌人，其中尤以肿瘤的危害为甚。世界卫生组织的报告表明：二 十世纪后期近三十年以来癌症发病一直呈上升的趋势。而我国相关的统计表 明：近年来，肿瘤已经成为死亡的主要原因。因此，开展对肿瘤治疗手段的研 究成为一个追切的课题。目前比较成熟的肿瘤治疗手段包括手术治疗、放射治 疗和化疗。其中放射治疗占了相当大的比例。 放射治疗特别适用于局部的具有一定外形的肿瘤。实施放射治疗的先决条 件是确定肿瘤的位置和判断肿瘤的生长范围。这可以由计算机断层扫描(凹)、 核磁共振(NMR)或正电子发射断层扫描(PET)等技术手段来提供病人的三维 解剖信息。在此基础上．医生确定出肿瘤区(CTV)和敏感组织(OARs)。再和 医学物理学家一起定义照射靶区(PTV)。然后，利用放射治疗计划系统作病人 的治疗计划。 在放射治疗的临床应用上，剂量计算是一个基本的问题。无论从数学方法 还是从计算效率的角度，精确地计算电子束在人体中的剂量的空间分布都是一 件十分困难的工作。目前使用的计算电子柬剂量的算法包括Monte Carlo方法 和解析算法。前者包括MMC(Macro Monte Carlo)方法和VMC(Voxe Monte Carlo) 方法等。而解析算法是在Fermi—Eyges多次散射理论的基础上发展起来的．比 如PBRA(EIectrod]3enci卜beam redefinition algoEithm以及本文所使用 的混合笔束模型。 的混合笔束模型。 高能电子输运的Fermi—Eyges理论是Fermi在研究宇宙射线穿透薄层介质 时用 个直观的方法推导出来的。Fermj-Eyges理论利用如下的假设：(1)电 子在与介质的原子或原子的电子发生相互作用时，主要经历弹性散射，电子只 小范围改变运动方向而不损失能量。这种假设忽略了非弹性散射和轫致辐射的 影r晒。c2)}乜子在介质中运动时，其总数是不变的。实际上，d{于电子在介质 I}Jj垂动时，将弓}起介质电子拍电离或激发，从而入射电子将逐渐损失能量。当 入射电予能量减少到一定值时，电予实际上停留在介质中。因而随着深度的增 加，入刺电子的总数目要逐渐减少。同时，轫致辐射也将引起入射电子数目的 减少。t∞由于假定入射电子不损失能量。因而在确定的深度处，所有的入射 电子具有相同的能量。这种假设只能限于浅层情况。在较大深度处，由于电子 能齄损失歧离和径迹弯曲效应，电子能谱有很大的展宽，因而不能再被视为单 能． 出于这螳原因，Fermi—Eyges理论主要存在两个不足，即不能给出电子数 掰∞损失以及电子注量的累积效应。这两个问题的存在使Fermi—Eyges理论难 咀满足现代放射治疗技术对剂量汁算精度的要求，因而很少被直接川’F剂量计 尊。}logstrom，Shiu e￡a J．发展的笔柬模型得到了比较满意的结粜。笔柬模 型利粥了实验测定的百分深度剂量作为{赶子笔束的中心轴剂量，而其侧向分布 划鹿用Fermi--Eyges理论的解，目前在放射治疗中得到广泛的应用。但是，由 于该算法篱要很多测量数据，因丽不便于应用到非均匀介质。一 电予输运取群模型是罗正明研究员在1985年提出的用来求解电子在介质中 的输运行为的一个模型。该模型首先建立电子的微分注量分布函数在介质中的 ∞“”mm、力程。然后电子的微分注量分稚函数被分为向前电子和扩散电：F的 微分江盛分布，：二者通过扩散电子源互相联系。相应地，介质中的电子被分为 向前f乜予群和扩敫电子群，其输运方程破分解为向前电子的输运方程和扩散电 子的输运专程。向前电子群只由在小角方向运动的电子组成。小角是指电子的 运动方向与人射电子的方向之问只有一个小的角度。由于向前电子群的运动方 囊偈离入射力向不大，因而在求解向荫电子的微分注量分布的Bol妇，f}a11n力程 时，可以引入窄能谱近似和连续慢化近似。而扩敌电子群由在各个方向运动的 时，可以引入窄能谱近似和连续慢化近似。而扩敌电子群由在各个方向运动的 电子组成。扩散电子的角度分布可以认为基本上是满足各向同性的。这样，在 求解扩散电子的微分注量分布的BoltzmaHn方程州1，_可以利用较低阶的球谐矩 来近似扩散电子的微分注量分布函数。将扩散电子的微分注量分布函数离散化 后，可以利用Lax—Wendroff格式或特征线方法来获得扩散电子微分注量分布的 数值艇。在分别获摄向前电子和扩散电子的微分注量分布后，就可以进一步获 得我们所需要的电