王广善,刘桂梅,李晓波,李祥勇,张旭
首都医科大学附属北京胸科医院放射治疗科,北京101149
国内有机构对不同加速器的束流特性进行对比研究,发现束流特性之间的相关性,为以后束流匹配奠定了基础[1-4]。国外很多放疗部门购买同一厂家相同系列的医用直线加速器,物理师利用模体对束流特性进行测量,对相同射线能量进行束流匹配,方便患者治疗[5-9]。近几年国内有放疗部门对束流匹配也做了相关研究[10-11],尝试进行束流匹配,但数量较少。我科拥有两台瓦里安直线加速器(简称治疗机),一台为Clinac iX(简称iX),一台为高端机Trilogy(简称TR),两台同为较成熟的C 系列型号,按照购置合同要求及束流匹配规程[12-13],对两台治疗机的6 MV X射线进行束流匹配,以方便一台治疗机出现故障时能将6 MV 计划切换至另一台治疗机执行,而不耽误患者治疗。本研究就两台治疗机6 MV 治疗计划进行分析,验证其数据传递的一致性及计划执行的精准性。基础数据的测量、比较结果及多叶准直器(Multi-leaf Collimator,MLC)DLG(Dosimetric Leaf Gap)和透射因子已在临床使用之前全部完成,此处因篇幅原因暂未列出,对于实际病例及治疗技术的验证更具临床意义。国外及国内束流匹配加速器均取得了较好的临床应用结果[8-11,14-16]。
iX 治疗机安装于2012年,其配有6、15 MV X 射线和4、6、9、12、15、18 MeV 电子线,各能量剂量率均为100~600 MU/min,三级准直系统,包括X、Y 钨门准直及60对千禧(Millennium)多叶准直器,其中中间40 对叶片在等中心宽度为0.5 cm、两边各10 对叶片在等中心的宽度为1.0 cm,Extract IGRT(Image-Guided Radiotherapy)碳纤维床板,可做三维适形、固定野调强、电子治疗。TR 治疗机安装于2018年,配有6、10 MV X 射线和6、9、12、15、18、22 MeV 电子线,剂量率、准直系统、治疗床与iX 配置相同。两台直线加速器使用同一套治疗计划系统(Treatment Planning System,TPS)Eclipse 13.6。
IBA DOSE1 剂量仪、CC13 电离室(灵敏体积0.125 cm3)、固体水(30×30)cm2、IBA MatriXX 二维电离室矩阵(简称Mt,空气电离室,灵敏体积0.08 cm3,面积24.4 cm×24.4 cm,1 020 个电离室,相邻间距0.762 cm)及配套分析软件OmniPro I'MRT、加速器自带的PortalVision 平板(简称PV,aSi1000 非晶硅,面积40 cm×30 cm,分辨率1 024×768)及Portal Dosimetry图像分析软件、统计分析软件SPSS 17.0。
病例选取,第一组:随机选取适合三维适形(3-Dimensional Conformal Radiotherapy,3D-CRT)治疗的病例,其中头颈部病人3例、胸部10例、腹盆部2例,共15例,按照临床要求对每例患者先设计一个基于iX的计划iX1[17],然后将该计划的治疗机改为TR,生成另一个计划TR1。第二组:随机选取适合固定野调强(Intensity Modulated Radiation Therapy,IMRT)治疗的病例,其中头颈部病人2例、胸部9例、腹盆部4例,共15例,设计一个基于iX的调强计划iX2,然后将该计划的治疗机改为TR,生成另一个调强计划TR2。
对每一个计划设计点剂量验证计划和面剂量验证计划,对每一个验证计划在各自的治疗机下执行,分别验证点剂量和面剂量。利用DOSE1 剂量仪、CC13 电离室和固体水采集点剂量;由MatriXX 平面矩阵和PV 平板分别独立进行面剂量验证。然后两台治疗机交叉执行计划,再次进行点剂量和面剂量验证。平面剂量验证采用γ 通过率分析法,选择3 mm/2%比较标准。
对所得数据采用SPSS17.0 软件进行分析,数据用均数±标准差表示,两组结果采用配对样本t检验,P<0.05为差异有统计学意义。
Eclipse计划系统自带计划评估及计划比较功能。图1为Eclipse计划比较界面,从3个不同层面可以看出二者剂量分布几乎相同。将每个病例iX1和TR1、iX2和TR2做对比,按照本研究的方案设计读取计划中的指标值,在剂量分布方面比较两计划之间的差异,比较的计划参数包括靶区PTV的D98%、D2%、Dmax和危及器官的相关指标,计算结果列于表1、表2、表3。两个计划相同指标的评价方法按公式(TR1-iX1)/iX1×100%或(TR2-iX2)/iX2×100%得出。
在PTV剂量方面,D98%、D2%、Dmax剂量差异最大平均值标准差分别为-0.52%±0.30%、0.53%±0.45%、-0.55%±0.17%,均小于1.00%,表现出非常好的一致性。危及器官的受量差异最大平均值出现在腹盆部IMRT 计划的膀胱D50%,其值为0.94%±0.84%,小于1.00%。两个束流匹配计划的靶区和危及器官剂量学差异能够达到临床可接受的范围,与部位和技术类型没有相关性。
图1 同一病例不同治疗机计划DVH及层面剂量比较Fig.1 Comparison of DVH and dose levels in different treatment machine plans of the same case
表1 头颈部iX和TR治疗机计划PTV和危及器官剂量差异比较(%,± s)Tab.1 Comparison of dosimetric difference between Linac iX and TR on PTV and OAR in head&neck (%,Mean±SD)
表1 头颈部iX和TR治疗机计划PTV和危及器官剂量差异比较(%,± s)Tab.1 Comparison of dosimetric difference between Linac iX and TR on PTV and OAR in head&neck (%,Mean±SD)
计划类型PTV D98%PTV D2%PTV Dmax 左晶体Dmax 右晶体Dmax 脑干Dmax 脊髓Dmax MU 3D-CRT IMRT-0.52±0.30-0.40±0.63-0.38±0.12-0.35±0.82-0.55±0.17-0.08±0.25-0.42±2.14-0.22±2.48-0.66±2.97-0.56±2.36-0.77±0.77-0.62±0.85-0.77±0.77-0.62±0.85-0.31±0.28-0.28±0.37
表2 胸部iX和TR治疗机计划PTV和危及器官剂量差异比较(%,± s)Tab.2 Comparison of dosimetric difference between Linac iX and TR on PTV and OAR in thorax (%,Mean±SD)
表2 胸部iX和TR治疗机计划PTV和危及器官剂量差异比较(%,± s)Tab.2 Comparison of dosimetric difference between Linac iX and TR on PTV and OAR in thorax (%,Mean±SD)
计划类型PTV D98%PTV D2%PTV Dmax 左肺V5Gy 右肺V5Gy 双肺V20Gy 双肺V5Gy 脊髓Dmax 心脏V30Gy MU 3D-CRT IMRT-0.36±0.23 0.02±0.50-0.08±0.34 0.05±0.61-0.11±0.38-0.05±0.77 0.54±0.41 0.84±0.65 0.20±0.38 0.79±0.55-0.55±0.76 0.54±0.62 0.92±0.74 0.93±0.44 0.93±0.24 0.71±0.04 0.41±0.39 0.29±0.69-0.37±0.43-0.20±0.40
表3 腹盆部iX和TR治疗机计划PTV和危及器官剂量差异比较(%,± s)Tab.3 Comparison of dosimetric difference between Linac iX and TR on PTV and OAR in abdomen&pelvis (%,Mean±SD)
表3 腹盆部iX和TR治疗机计划PTV和危及器官剂量差异比较(%,± s)Tab.3 Comparison of dosimetric difference between Linac iX and TR on PTV and OAR in abdomen&pelvis (%,Mean±SD)
计划类型PTV D98%PTV D2%PTV Dmax 膀胱D50%小肠D2cc 直肠D50%脊髓Dmax 左股骨头V40Gy 右股骨头V40Gy MU 3D-CRT IMRT-0.21±0.63 0.41±0.40-0.34±0.35 0.53±0.45-0.38±0.46 0.40±0.75 0.11±0.30 0.94±0.84-0.16±0.21 0.34±0.38 0.01±0.05 0.45±0.83-0.01±0.03 0.00±0.00-0.86±0.26 0.15±0.31 0.08±0.24 0.08±0.18 0.93±0.16 0.12±0.45
将模体中测量点剂量与TPS 同一位置的点剂量进行比较,以TPS 为基准计算百分比偏差。8 种组合方式点剂量验证结果列于图2。数据显示,iX 治疗机执行iX1、iX2、TR1、TR2剂量偏差百分比平均值及标准差分别为1.35%±0.61%、1.23%±0.53%、1.21%±0.40%、1.66%±0.44%,最大偏差2.36%。TR 治疗机执行TR1、TR2、iX1、iX2 计划通过率平均值及标准差分别为 -1.06%±0.75%、-1.26%±0.60%、-1.36%±0.67%、-0.87%±0.67%,最大偏差-2.80%。对以上8 种治疗机与计划的组合分别编号A1、A5、A4、A8、A3、A7、A2、A6,见图2。A1 和A2、A3 和A4、A5 和A6、A7 和A8分别进行两样本配对t检验,P值分别0.006、0.010、0.026、0.039,差异有统计学意义。
相对剂量采用MatriXX平面矩阵和PV平板两种验证工具独立验证,它们在使用之前均经过校正和刻度,保证其准确性。γ 分析采用3 mm/2%标准,MatriXX 矩阵对iX 和TR 治疗机验证的结果列于图3,iX治疗机执行iX1、iX2、TR1、TR2计划通过率平均值及标准差分别为97.93%±0.96%、97.29%±1.07%、99.02%±0.99%、98.40%±0.76%,最小值为95.40%。TR 治疗机执行TR1、TR2、iX1、iX2 计划通过率平均值及标准差为97.81%±0.86%、97.47%±0.80%、98.33%±0.92%、98.25%±0.81%,最小值为95.90%。对以上8 种治疗机与计划的组合分别编号B1、B5、B4、B8、B3、B7、B2、B6,见图3。B1 和B2、B3 和B4、B5 和B6、B7 和B8 分别做两样本配对t检验,P值分别为0.012、0.009、0.006、0.009,差异有统计学意义。
图2 iX与TR治疗机执行本机及交叉执行计划绝对剂量验证结果Fig.2 Absolute dose verification results of local and cross execution plans of Linac iX and TR
图3 MatriXX对iX与TR治疗机执行本机计划及交叉执行计划相对剂量验证结果Fig.3 MatriXX′s local plan and cross plan relative dose verification results for Linac iX and TR
PV 对iX 和TR 治疗机的IMRT 计划再次做独立验证,结果列于图4,iX 治疗机执行iX2、TR2 计划通过率平均值及标准差分别为96.92%±0.60%、96.97%±0.49%,最小值为96.20%。TR 治疗机执行iX2、TR2计划通过率平均值及标准差为97.40%±0.63%、97.28%±0.80%,最小值为95.90%。对以上4 种治疗机与计划的组合分别编号C1、C4、C2、C3 见图4。C1和C2、C3和C4分别做两样本配对t检验,P值分别为0.009、0.008。
图4 PV对iX与TR治疗机执行本机计划及交叉执行计划相对剂量验证结果Fig.4 PV′s local plan and cross plan relative dose verification results for Linac iX and TR
本科室在安装Trilogy的同时并与Clinac iX 进行束流匹配调试。束流匹配的最大优点在于,其中一台治疗机突然故障停机时,患者可以随时切换到另一台治疗机继续治疗,而不需重新设计计划,保证患者治疗疗程完整。另外,同一系列的加速器在束流匹配后,若其中一台治疗机故障,另一台可以作为参考,查看参数设置,互换零件,极大地提高了维修效率,为临床治疗争取时间。
有研究表明,剂量传递过程中5%的差异就会导致肿瘤控制率10%~20%的变化和正常组织20%~30%出现放射副作用的可能性[18]。有研究者提出必须满足剂量处方点的误差≤3.5%,整个计划靶区其他点的剂量误差≤5%,才能满足患者靶区剂量的准确[19]。从计划质量的比较来看,剂量分布剂量曲线重叠。因为在TR 调试过程中,束流参数要一一匹配精准,比如射线质D20/D10相差0.20%、射线对称性小于1.00%等。另一方面,同一套计划系统Eclipse13.6也保证了算法与建模的一致性,保证基础数据准确无误,治疗计划便符合预期结果。Seonghee 等[14]同样对束流匹配的60 例不同部位的治疗计划进行比较,得到的结果是PTV 和危及器官剂量差异都在1.00%以内。Krishnappan 等[15]对6 个未做束流匹配但同一生产厂商的加速器进行计划比较发现,选取的病例使用3D-CRT、IMRT 和容积旋转弧形治疗的计划最大剂量偏差为2.60%,平均偏差小于1.90%,并且剂量偏差与使用技术类型没有相关性。
图2数据表明,两台治疗机在执行所有计划方面均能准确地传递MU 和调制参数,iX 和TR 的点剂量偏差均值都小于3.00%,差异有统计学意义(P<0.05),两台治疗机对计划的执行略有不同,但差异在可接受范围内。iX 实际测量剂量比TPS 计算偏高,TR 实际测量剂量比TPS计算偏低。原因可能来自两方面,一是治疗机本身方形射野输出剂量准确,执行计划时剂量率、MLC 等调制参数共同配合下剂量传递引入误差;二是受温度、湿度、气压等因素影响的验证工具的精度及摆位误差。Amin等[20]在研究用于测量束流匹配加速器输出剂量的电离室系统稳定性和一致性中发现,不确定的来源主要来自于测量摆位、刻度测量和实际测量之间的温度差异,这些因素综合产生了0.6%的测量不确定度。Ashokkumar 等[9]对3 台束流匹配的加速器进行了剂量学研究,验证结果在±3%的容差值以内。
在相对剂量验证方面,本研究采用比较严格的3 mm/2%的γ 分析标准,图3和图4的数据表明,MatriXX和PV验证通过率结果平均值分别在97%和96%以上,iX最小值95.40%,TR最小值95.90%,这和国外一些研究得到了一致的结果[8-9,14],结果经过t检验,差异有统计学意义(P<0.05),说明交换执行计划确实存在统计学意义上的差异,造成这种差异的原因有很多,比如两台治疗机的MLC运动精度、剂量率控制等。Oliver 等[21]报道,由于1 mm 的随机MLC 到位误差,剂量传递的最大变化约为1.21%。最终测量结果是很多因素综合产生的,但差异甚微,可以接受。两种独立的验证方式均有力地证明了束流匹配的可靠性与可行性。分析中发现PV 通过率偏低,并且比MatriXX 通过率低,这与MLC 叶片间隙漏射线比叶片透射线多关系密切,对于所有射野机头角度都为0°的计划,每个射野MLC叶片间隙漏射线叠加,PV平板分辨率很高,能响灵敏,因此通过率低的区域出现在MLC叶片间隙处。如果每个计划中的每个射野机头角度不同,则可减弱MLC 叶片间隙漏射线对通过率的影响。Seonghee 等[14]利用EPID 对束流匹配的60 例不同部位的治疗计划进行验证,2%/2 mm和2%/3 mm标准下的γ通过率分别高于93%和97%。全面的束流特性、计划差异和预治疗验证结果均在可接受的临床范围内,和TPS 计算展示了很好的一致性。
两台经过束流匹配的治疗机之间交换执行头颈、胸部和腹盆计划的剂量学分析在临床可接受的范围内。这些结果为在繁忙的临床环境中通过与束流匹配的治疗机交换治疗患者提供了重要的质控保证,而无需重新设计计划。虽然两台匹配束流的加速器交叉执行计划可行可靠,但每台治疗机有独立的束流数据,日常治疗时计划必须在对应的治疗机上执行,以免不可预料的错误发生。只有在治疗机故障等万不得已的时刻,才采取交叉治疗的方案,治疗前必须做锥形束计算机断层扫描,并与计划CT 进行比较,以匹配最佳治疗区域。另外,交叉实施之前,必须确保计划匹配的验证[16]。