基于笔射束核阵列的三维适形放射治疗剂量计算方法研究

摘要

肿瘤放射治疗学与肿瘤外科学、肿瘤内科学一起构成恶性肿瘤治疗的三大主要手段,约有60％～70％的肿瘤病人在不同的治疗阶段需要接受放射治疗。放射治疗的目的是达到治疗的最优化,即靶区所受剂量尽可能大,而正常组织和危及器官尽可能少受或免受照射。计算机三维放射治疗计划系统是当今放射治疗的灵魂。本文针对放射治疗计划系统中的剂量计算模型进行了系统而深入的研究。
　　 伴随着放射肿瘤学的发展,关于精确放射治疗技术的研究从未停止过。1949年,瑞典科学家Leksell首先提出放射外科学的构想,利用立体定向技术,使用大剂量的射束一次性摧毁需治疗的病灶,此后这项治疗方法被命名为立体定向放射外科(SRS),并在1955年应用于临床。上世纪60年代初期,日本科学家用一种电机驱动光栅控制放射线束的形状,使得射束的截面正好吻合病灶的投影。70年代初期,美国科学家开始认识到调强放疗(IMRT),的重要性并对其进行了初步研究。80年代末90年代初,美、德等国家相继开发出了商用的X刀；瑞典开发了第三代γ刀系统；同时期带有机械手的“Cyber刀”出现。1995年Spirou等人提出了使用动态多叶光栅来实现调强放疗,并用其治疗前列腺癌。进入新世纪以来,随着放射物理学、放射生物学、临床肿瘤学和医学影像学等相关学科的发展,放射治疗技术领域发生了巨大变革。精确放射治疗技术就是以“精确定位、精确计划、精确治疗”为特征的新的放射治疗技术的统称。精确放射治疗技术使得在肿瘤治疗上实现高精度、高剂量、高疗效和低损伤(三高一低)的现代放疗模式成为可能。
　　 精确放射治疗技术较之常规放射治疗技术有很大的进步,主要表现在:
　　 (1)精确放射治疗技术一般采用CT、MRI、PET等先进的影像技术作为确定靶区与重要器官的依据,而常规放射治疗技术一般采用X光模拟定位图象；
　　 (2)精确放射治疗技术一般具有三维治疗计划系统(TPS),需对影像数据进行三维重建和三维显示,并在此基础上进行三维计划设计,显示三维剂量分布,而常规放射治疗技术一般只能进行二维治疗计划设计；
　　 (3)精确放射治疗技术一般采用面膜、真空负压垫、头环等固定装置结合三维坐标系统,可对病灶进行单次或多次重复的准确定位,而常规放射治疗技术一般采用体表勾画病灶等比较粗糙的定位方法；
　　 (4)精确放射治疗技术一般采用多个非共面照射野或多个非共面照射弧的照射方式,而常规放射治疗技术一般采用单野照射或两野、三野照射；
　　 (5)精确放射治疗技术可使高剂量区域与靶区在三维形状上更加一致,并保证靶区内及表面的剂量处处相等。也即实现剂量上的适形,常规照射治疗技术则无法达到。
　　 剂量计算的结果作为评价治疗计划的依据,同时也是治疗计划优化与逆向治疗计划设计的基础。在计算机放射治疗计划系统迅速发展的同时,其核心内容之一的剂量计算模型和方法也在不断地发展。这种发展主要表现为从基于修正的剂量计算模式向基于模型的剂量计算模式的转变以及二者的相互渗透。基于修正的剂量计算模式是以参考条件下标准野的实验测量深度剂量、离轴比、散射因子等为基础,加之必要修正后得到实际放射物理条件和放射治疗对象剂量分布的剂量计算方法。随着计算机存储容量和运算速度的不断提高,其它更高精度的剂量计算方法不断地被用于计算机放射治疗计划系统中。这类计算方法一般以建立新的剂量计算模型为着眼点,全局考虑放射线与人体的相互作用,从而使基于理论计算的剂量分布与基于实验测量的剂量分布更吻合。文中对这两大类方法进行了较为详细的介绍,特别是集中对几种基于模型的剂量计算方法进行了详细的分析和研究。其中包括:将模体中任意一点的剂量分成原射线剂量和散射线剂量两部分的剂量计算方法被称作原射线剂量和散射线剂量分离法。此类计算包括了Clarkson法、微分散射空气比法、Day氏法等。
　　 相比于上述几种方法,学术界普遍认为模特卡罗法为目前剂量计算的金标准。该方法进行剂量计算的原理为用统计学方法来模拟大量单能光子在输运过程中与物质的作用。光子进入某种介质(如人体组织)后,通过与介质中原子的相互作用,而传递电离辐射的部分或全部能量。该方法的先进之处主要表现在它能适用于任意几何形状的体模和各种次级粒子(如光子辐射产生的电子)上。
　　 另一个比较精确的剂量计算模型是卷积法。在这种模型中,剂量计算是通过将放射线的光通分布与一个点扩展函数或卷积核相卷积来实现的。卷积法适用于任意的射野光通分布,所以这种方法特别适宜于不规则野的计算。另外,因为卷积可通过快速傅立叶变换来完成,这使卷积运算速度得到了显著地提高。
　　 通过对上述几种基于模型的剂量计算方法进行了分析,选定采用卷积模型完成剂量计算,并且对该模型的原理进行了深入的研究,在该方法已有的理论基础上进行了改进,提高了模型的计算精确度。
　　 卷积模型作为当前应用于放射治疗计划系统中的主流计算模型之一。当暂不考虑卷积核随深度的变化情况时,简化为平面内的注量分布与笔射束核的卷积。通过计算可以得到某一深度处于射束方向垂直的平面内的剂量分布。对于卷积运算中所需要用到的笔射束核,根据Fermi理论,在小角度散射的情况下,可以表示为两个相互独立的正交方向上的一维笔射束核的乘积形式。
　　 一维笔射束核可以通过实际测量的宽平行束的剂量分布数据求解得到。再利用笔射束核重建各种形状射野的剂量分布,具有一定的计算精度的优势,减少了额外工作量。该方法的优点是:测量数据综合了照射过程的各种散射因素,模型不存在理想化和简化的问题。利用笔射束核重建测量野自身或受照条件与之相近的射野的剂量分布效果较好。但采用单一笔射束核重建任意射野剂量分布仍存在一定问题。
　　 进入模体内的光子与介质相互作用产生次级电子将能量传递给介质,在建成区附近达到电子平衡。卷积模型中假设笔射束核具有空间不变性的这一前提,对直线加速器而言在达到电子平衡的区域外难以得到保证。在直线加速器的剂量分布中,散射线的剂量贡献占有极其重要的一部分。均匀模体中散射线剂量受射野大小和深度,以及同一方向上准直器的近端和远端边界的影响,其变化较为复杂。在较大深度处,模体内散射占主要；较浅深度处,受准直器散射影响较多。由于各种干扰因素的存在,造成射线能谱发生改变,使实际笔射束核的剂量分布并非是一个空间不变量,受到笔射束核为空间不变量这一假设的限制,对其进行修正,会减弱卷积模型可以通过快速傅里叶变换完成的这一优点。利用单一的笔射束核对任意射野的剂量分布进行重建,因为对射线能谱和射线质的改变考虑不够全面,造成计算结果和实际测量值之间有时会存在较大误差。本文通过对多个由测量数据反卷积得到的笔射束核做进一步处理,根据射野大小和计算平面深度做成阵列。在计算过程中,按照目标射野的照射条件,拟合出近似的笔射束核,在保证计算速度的同时,最大限度的保证了模型的计算精确度。实验证实,基于笔射束核阵列的剂量计算能够取得比较满意的结果,是一种可行的解决方案。
　　 在论文的最后,我们对研究工作中的遗留问题进行了讨论,并对今后的工作进行了展望。