Varian Trilogy直线加速器6MV光子线治疗头的蒙特卡罗模拟

宋炳文，潘建南，吴思宇，朱凤盈，梁嘉乐

南方医科大学顺德医院放射治疗科，广东佛山528300

前言

作为肿瘤治疗的主要三大手段之一，放射治疗的目标是给予肿瘤足够精确的照射剂量，同时使正常组织和危及器官尽可能少地受到照射［1］。在现代放射治疗实践中，承担这一重要任务的主要为治疗计划系统（Treatment Planning System，TPS）［2］。其中剂量计算模块是TPS的核心，直接关系到剂量投送到患者身上的准确性，目前公认的最为准确的剂量计算方法为蒙特卡罗方法，蒙特卡罗剂量算法也是目前医学物理研究的热点之一［3］。蒙特卡罗剂量计算的关键是构建准确的辐射源模型，而目前传统光子线放射治疗的主要辐射源为医用电子直线加速器［4-8］。本研究利对Varian Trilogy直线加速器6 MV光子线治疗头进行蒙特卡罗模拟，为今后进一步的研究奠定基础。

1 材料与方法

1.1 MCNP5简介

MCNP5是由美国洛斯阿拉莫斯国家实验室开发的一款通用程序包，主要用于计算三维复杂几何结构中的中子、光子、电子或者耦合中子/光子/电子输运问题，其拥有丰富的粒子反应截面数据库以及多种减方差技巧［9］。此外，MCNP5 提供了多种几何描述方法，特别是用于描述规则几何体的“Macrobodies”非常方便地描述了加速器治疗头中的部件，此外，通过对面或者体进行交、并、补等逻辑运算，MCNP5可以方便地构建复杂几何结构，并提供多种计数方式，因而广泛应用于核科学领域。目前医学物理中常用的几种大型通用蒙特卡罗程序有EGS4、MCNP、GEANT4 等，相比于其他的蒙特卡罗程序，MCNP5在几何构造、减方差以及源设置上更加全面［4］，因此本研究采用MCNP5 程序来对加速器治疗头进行建模和剂量计算。

1.2 治疗头建模

作为放射治疗中外照射的主要放射源，医用电子直线加速器治疗头中射束产生的主要工作原理：由电场加速的电子经过偏转磁铁的偏转后，打到金属靶上，高速运动的电子与靶物质中的原子发生相互作用，电子在突然减速过程中发生轫致辐射产生X射线，产生的光子束经过初级准直器和均整块后，由次级准直器和多叶光栅形成所需的射野形状［1-2］。其中初级准直器位于靶的下方，中间一般为圆锥形镂空，轫致辐射产生的光子束经过初级准直器准直后形成锥形束，初级准直器的材料一般为高Z 材料，以阻挡和吸收光子。由于最初产生的射束分布式中央强度高，边缘强度较弱，因而在初级准直器的下方加上中央厚边缘薄的锥形状均整器，使最终光子束的强度分布是平坦并且对称的。均整器的材料一般也是钨、铅等高Z材料［4-8］。

随着计算机技术的提高以及蒙特卡罗程序的不断发展，医用电子直线加速器的模拟也经历了一个由简单到复杂的发展过程，从最初的点源模型到后来的射束部件模型，再到如今的详细加速器模型［4,10］。此前，已经有很多的研究人员使用各种不同的蒙特卡罗程序对不同厂商的直线加速器治疗头进行了模拟，不同的学者提出了不同的模拟方法，目前常用的模拟方法是先通过模拟电子打靶获得出射X射线的能谱，再将能谱用于下一步的剂量计算，甚至在国际原子能机构官网上可以下载使用蒙特卡罗程序建立的一些不同型号的加速器治疗头的相空间文件，这些相空间文件已经包含了次级准直器及其以上的所有部件，并且用户可以直接使用EGS 或者GEANT4 程序进行读取，从而减少模拟的复杂性［11-15］。考虑到本单位所使用的直线加速器型号以及后续能够更方便地进行复杂的剂量计算，参考Bednarz 等使用MCNPX 建立Varian Clinac 2100C型直线加速器详细的6 MV和18 MV光子线治疗头的思路，确定本研究的主要模拟步骤如图1所示［4-10］。其中Varian Trilogy直线加速器治疗头主要射束部件的几何参数和材料参数由厂商（Varian Medical Systems）提供，屏蔽部件的参数主要来自于Bryan等［4］和Kry等［5］描述，百分深度剂量曲线和横向剂量分布曲线来自于科室新机器验收时水箱的测量数据。

图1 治疗头蒙特卡罗模拟流程图Fig.1 Flow chart of Monte Carlo simulation of treatment head

医用直线加速器治疗头的蒙特卡罗模拟主要包括治疗头的几何建模以及初始电子束参数的调整。按照射束经过的依次顺序，医用直线加速器治疗头的主要部件包括靶、靶衬底、初级准直器、均整器、次级准直器、多叶光栅（Multi-Leave Collimators,MLC）以及屏蔽材料等超过20 个部件被建立在MCNP5 中，构建好的治疗头模型如图2所示，其中坐标原点取在靶下表面的中心，正Z 轴方向为靶指向患者的方向，正Y 方向为臂架指向患者的方向，正X 轴方向为当患者头部仰卧时患者左侧的方向。由于验收时测量的主要为方形野，为了简化模拟难度，在此次模拟中MLC 被等效成两个长方体钨金属块。对于沿弧形移动的两对次级准直器，在MCNP5 中可以使用TR 卡将次级准直器移动到相应射野大小时的位置。

1.3 MCNP5参数设置

对于MCNP5中物理参数的设置，程序可以通过调用“MODE P card”来使用详细的光子输运模型对光子与物质的相互作用进行模拟，该模型将相互作用过程处理为光电效应、电子对产生以及康普顿散射等［9］。电子输运的物理过程计算通过调用“MODE E card”实现。为了减少计算结果的统计涨落，节约模拟时间，使用“CUT card”将光子和电子的能量截断分别设置为“0.01 MeV”和“0.1 MeV”。此外，由于直线加速器治疗头模型涉及到的部件众多，特别是对次级准直器进行移动时，如果部件的几何边界发生重叠，会因为边界碰撞而导致MCNP5在运行期间大量粒子丢失，最终出现致命错误导致输出结果不能通过所有的统计检查而全部置为零，可使用“VOID card”使所有栅元的材料变成真空，防止几何内发生任何粒子碰撞［4-5］。

图2 加速器治疗头模型剖面图Fig.2 Profiles of accelerator treatment head model

1.4 初始电子束参数调整

准确模拟出治疗头的束流参数，需要知道该直线加速器初始电子束的参数。由于每台机器都不相同，并且同型号的直线加速器在不同的单位进行安装调试时，束流参数也不完全一样，因此还需根据实际使用的加速器束流测量参数对加速器治疗头的初始电子束的参数进行调整［10-11］,调整的部分主要为入射电子束的平均能量以及半高宽（Full Width at Half Maximum,FWHM），其中对比的测量数据来自于新机器验收时的水箱测量参数，为了便于与计算值对比，对测量数据进行了一定的插值处理。对于治疗头模型中百分深度剂量曲线和横向剂量分布的计算，采用了一个尺寸接近于实际水箱的（50×50×50）cm3的水模体，源皮距为100 cm，射野中心位于水模体的上表面中心。

2 结果与分析

计算百分深度剂量曲线时，由于要考虑光子在水模体中的剂量建成效应，采用了*F8能量沉积计数卡，计数的体素尺寸沿射野中心轴方向为（4×4×2）mm3。横向剂量分布的计算使用了1型“pedep”电子网格计数，不同深度处计数的网格大小同样为（4×4×2）mm3。百分深度剂量曲线的模拟计算中考虑了不同尺寸的射野，分别为（5×5）、（10×10）、（20×20）、（40×40）cm2。不同尺寸射野的横向剂量分布计算考虑了不同的深度，分别为1.5、5、10 cm。模拟的初始电子数为5×108个，所有的模拟均在Intel Core i7-6700K（主频4 GHz，四核心八线程）处理器上运行，采用多线程并行方式，每个程序运行时间约为60 h，所有模拟结果的统计误差均在5%以内，可以认为是可信的［9］，所有的计数结果均归一到最大吸收剂量处。

相关研究表面，直线加速器初始电子束基本为高斯分布，水模体中的剂量分布强烈依赖于入射到加速器靶上面的电子束的平均能量以及径向强度分布，深度剂量曲线强烈依赖于初始电子束的径向分布，并且对于最初的平均能量相当不敏感，横向剂量分布曲线取决于初始电子束的平均能量［2,10-11］。因此，本研究在模拟过程中通过几次迭代手动调整初始电子束的平均能量与径向分布，直到计算的百分深度剂量曲线和横向剂量分布与实际测量值的差异性小于2%［12-13］。优化前后的初始电子束参数如表1所示。

表1 优化前后初始电子束参数对比Tab.1 Comparison of initial electron beam parameters before and after optimization

不同射野大小下的百分深度剂量曲线结果如图3所示，计算值和测量值的最大差异性约为2%，其中（10×10）cm2射野下的最大剂量点深度约为1.5 cm，与实际测量值基本相符。

横向剂量分布结果如图4所示，在射野较大[（20×20）、（40×40）cm2]时，剂量曲线呈现明显的马鞍形趋势，射野范围内计算值和测量值的最大差异性约为2%，半影区和射野外的剂量分布差异性相对大一些。

3 讨论与结论

本研究基于MCNP5 程序对Varian Trilogy 直线加速器6 MV 光子线治疗头的几何结构和材料成分进行了建模，构建了一个详细的加速器治疗头模型，并对程序中的一些物理参数进行设置使模拟过程尽量接近粒子的实际输运过程。研究发现，水模体中的深度剂量曲线对于初始电子束的平均能量并不敏感，强烈依赖于径向分布，并且FWHM 越大，同一深度处的百分深度剂量增加，横向剂量分布曲线取决于初始电子束的平均能量，通过对入射电子束的平均能量以及FWHM 进行调试，使得束流参数与实际测量值达到比较好的符合度，研究结果达到了目前同类研究中的模拟精度［12-20］。因此，该模型可以用于更进一步的剂量计算以及机房和操作人员的辐射防护研究。

直线加速器治疗头是一个很复杂的系统，其中包含了很多不同的部件。值得注意的是本研究的模拟时间相对于其他研究略长［10］，这主要是因为模拟了电子的打靶过程，实际上电子通过轫致辐射产生光子的效率较低，导致本研究中需要对更多的初始电子进行抽样才能达到在统计误差内的结果，如何提高模拟效率也是后续研究中值得关注的问题。本研究中由于实验条件的限制，对于MLC 进行了等效处理，这种情况在研究方形野的情况下比较合适，后续的研究中将尝试对MLC进行精确的模拟以适应不规则射野的剂量计算条件。

图4 横向剂量曲线计算值与测量值的对比Fig.4 Comparison of calculated and measured values of lateral dose distribution