中华放射医学与防护杂志  2015, Vol. 35 Issue (9): 665-669   PDF    
基于多叶准直器日志文件的鼻咽癌调强放疗的累积剂量验证
曹瑛1, 杨振2 , 邱小平1, 杨晓喻2, 唐慧敏1, 吕知平2, 张子健2, 雷明军2, 刘归2, 胡永梅2    
1. 421000 衡阳, 南华大学核科学技术学院;
2. 中南大学湘雅医院肿瘤科
[摘要]    目的 利用多叶准直器(MLC)日志文件验证鼻咽癌患者动态调强放疗的累积剂量,分析MLC叶片位置偏差与剂量学关系。方法 跟踪记录5例鼻咽癌患者的28个分次动态调强放疗,共计1 400个日志文件。通过Argus软件解析和提取MLC叶片的实际位置信息,并与MLC叶片的计划位置信息比较,计算每个MLC叶片的位置偏差和所有患者的整体平均叶片位置偏差。将MLC文件导入治疗计划系统(TPS)进行剂量重建计算,分别比较和分析每个病例重建计划和原计划的28次累积剂量差异及单次相对剂量差异。结果 单个叶片位置偏差<0.50 mm约占90%,0.50~1.00 mm之间的占4%~11%,1.00~1.50 mm约占8‰,整体平均叶片位置偏差为0.16 mm。与原始计划累积剂量相比,重建计划的叶片位置偏差导致靶区累积相对剂量差异 < ±0.1%,危及器官的累积相对剂量差异 < ±0.6%,差异均无统计学意义(P>0.05);而与原始计划的单次剂量相比,重建计划的靶区和脑干、脊髓和右眼晶状体单次剂量差异有统计学意义(z=-2.02~4.61和4.46、-4.51、2.07, P<0.05)。结论 基于动态调强日志文件的分次累积剂量验证,一定程度上减小了单次MLC叶片位置偏差对剂量的影响,验证结果可更加精确地反映患者的实际受照剂量。
[关键词]     鼻咽癌    调强放疗日志文件    位置偏差    剂量学    
The IMRT accumulative dosimetric verification of NPC using MLC log files
Cao Ying1, Yang Zhen2 , Qiu Xiaoping1, Yang Xiaoyu2, Tang Huimin1, Lyu Zhiping2, Zhang Zijian2, Lei Mingjun2, Liu Gui2, Hu Yongmei2    
1. School of Nuclear Science and Technology, University of South China, Hengyang 421000, China;
2. Department of Oncology Xiangya Hospital, Central South University, Changsha 410008, China
[Abstract]    Objective To verify the accumulative dose of intensity-modulated radiotherapy (IMRT) of nasopharyngeal carcinoma (NPC) using linac logfiles, and analyze the relationship between the leaf position errors of multileaf collimator (MLC) and dose distribution changes. Methods In total, 1 400 dynalog files of 28 fractions of 5 NPC patients were tracked and recorded. MLC leaf position information was compared with the planned MLC leaf position information to calculate the MLC leaf position errors of each leaf and the overall mean leaf position error of all patients, which were analyzed and extracted through Argus program. The MLC files were imported into treatment planning system (TPS) to recalculate the dose. Twenty-eight accumulative dose differences between original plans and reconstruction plans, and single relative dose differences were analyzed. Results About 90% of the leaf position errors were less than 0.50 mm, 4%-11% of the leaf position errors were between 0.50-1.00 mm, while there were about 8‰ leaf position errors within 1.00-1.50 mm and the average leaf position error was 0.16 mm. Compared with accumulated dose of original plans and that of reconstruction plans, the average relative dose discrepancies in targets and organs at risk (OARs) induced by leaf position errors were less than±0.1% and±0.6%, and the differences were not statistically significant (P>0.05). Compared with single dose of original plans, single dose differences in targets, brainstem, spinal cord and right lens of reconstruction plans were statistically significant (z=-2.02-4.61, 4.46, -4.51, 2.07, P<0.05). Conclusions The accumulative dose verification based on multi-fraction dynamic intensity-modulated MLC log files, to some extent, lessen the dose influence of the single-fraction leaf position errors and the verification results could manifest the practical irradiation dose more accurately.
[Key words]     Nasopharyngeal carcinoma    IMRT dynalog file    Position error    Dosimetry    

剂量验证是保证鼻咽癌调强放疗安全可靠的一个重要环节。利用多叶准直器(MLC)产生的日志文件(dynalog file)提取实际治疗时的叶片序列,检测静态、动态调强以及容积调强技术的机械稳定性及剂量验证,方便快捷[1, 2, 3, 4, 5, 6]。叶片位置偏差对剂量的影响呈线性相关[7]。随机误差对剂量的影响很小,系统误差对剂量的影响明显[8]。这些研究结果是基于单次(实际治疗的1个分次)剂量重建得出,而实际的放疗包括多个分次的累积。在整个疗程中,实际的随机误差不可避免地累积或抵消,系统误差可能累积,也可能由于日常的质量保证消除或减小。因此,本研究利用MLC日志文件对鼻咽癌动态调强放疗的全程累积剂量进行验证,明确全疗程中叶片位置偏差的大小,并分析其对靶区和危及器官剂量的影响。

资料与方法

1. 临床资料:随机选取5例鼻咽癌患者的调强放射治疗(IMRT)计划。靶区和危及器官由临床医生勾画。肿瘤靶区(GTV)包括原发灶(GTVnx)和转移淋巴结(GTVnd),分别外放5和0~3 mm作为计划原发灶靶区(PGTVnx)和计划淋巴结靶区(PGTVnd);临床靶区(CTV)定义了高危区(CTV1)和预防照射区(CTV2),分别外扩3 mm得到计划靶区(PTV1和PTV2)。靶区处方剂量按照同步推量方式给予,PGTVnx:7 095~7 392 cGy,33次;PGTVnd:6 996~7 095 cGy,33次;PTV1:6 105~6 270 cGy,33次;PTV2:5 040~5 152 cGy,28次。治疗分两个阶段,第1阶段28次,第2阶段舍弃PTV2再照射5次,由于第2阶段强度调制用的叶片比第1阶段少,为了更严谨地研究叶片位置偏差,只跟踪第1阶段28次日志文件。

2. 叶片位置偏差描述:利用美国瓦里安公司生产的Varian Trilogy直线加速器,配备千分(Millennium)120叶MLC,共A、B两组(平行组),每组各有60个叶片。调强放疗照射期间,系统间隔50 ms记录1次MLC叶片位置信息,并写入日志文件。该文件是按时序生成的一长串ASCII字符,可读性差,借助Varian Argus软件对其进行解析和数据提取,便于逐个叶片与计划 MLC文件比较。为了表征全疗程MLC的整体位置偏差,定义了整体平均叶片位置偏差 ± s,A、B组单个叶片平均位置偏差 ± sA ± sB,和每对叶片的平均叶片间隙差异△Gap,相关计算公式如下。

式中,x1xp分别表示日志文件记录的实际和计划叶片的绝对位置,cm;NmNcNfNp分别表示运动叶片数、控制点数、射野数、病例数;XLAXLBXPAXPB分别表示日志文件记录的A、B组实际和计划叶片的绝对位置,cm;△Gap>0、=0和<0分别表示实际叶片间隙缩小、不变和扩大。

3. 剂量验证:剂量验证的评价指标有95%靶区受到的最低剂量D95%、靶区平均剂量Dmean,脊髓、脑干、眼晶状体的最大剂量Dmax,以及靶区和危及器官剂量体积直方图(DVH)差异。Siebers等[9]提出当剂量差异达到3%~5%时,可能导致靶区欠量或并发症发生概率增加,因此,定义相对剂量差异来评估叶片位置偏差对剂量的影响程度:相对剂量偏差(%)=(DP-DL)/DL×100%。式中,DPDL分别为计划和实际剂量,cGy。

对单个病例,将每个射野的日志文件通过Argus解析后导出为MLC格式文件,用C+[KG-*4]+编程进行转换,生成Eclipse兼容格式的MLC文件。在Eclipse中复制原计划,并导入MLC文件,进行剂量重建计算,得到重建计划。需要注意的是,跳数和原计划保持不变(相同的有效位)、算法和计算网格(2.5 mm×2.5 mm)相同。每个病例各得到28个重建计划,并分别进行计划叠加,得到基于日志文件的28次累积剂量,和原治疗计划28次累积剂量做比较和分析。随机选取1例患者的4次日志文件,比较4个单次(第1次~第4次)和累积计划中靶区与危及器官相对剂量差异的大小。

4. 统计学处理:采用SPSS 17.0软件分别对重建计划与原计划的累积剂量和单次剂量进行Wilcoxon秩和检验。P < 0.05为差异有统计学意义。

结 果

1. MLC叶片位置偏差:5例患者各机架角的平均叶片位置偏差、单个病例的平均叶片位置偏差和整体平均叶片位置偏差列于表1。5例患者各机架角的平均叶片位置偏差在0.14~0.20 mm,无明显的角度依赖性;每例患者的平均叶片位置偏差在0.15~0.17 mm之间;整体平均叶片位置偏差为0.16 mm。各机架角MLC叶片位置偏差、两组各叶片的平均位置偏差和平均叶片间隙差异结果分别示于图12。5例患者各28次治疗,逐个叶片统计得出叶片位置偏差<0.50 mm约占90%,0.50~1.00 mm之间的占4%~11%,1.00~1.50 mm约占8‰,未发现>1.50 mm的情况。A、B组的平均位置偏差分布趋势一致,但A组叶片位置偏差基本是B组叶片位置偏差的两倍;ΔGap>0,表明整体上实际的叶片间隙有所减小;另外,由于MLC两端叶片未参与剂量调制,位置偏差明显减小。

表 1 5例患者在不同机架角度的叶片平均位置偏差(mm)

图1 各机架角MLC叶片位置偏差

图2 两组各叶片的平均位置偏差和平均叶片间隙差异

2. 累积剂量学差异:靶区与危及器官的累积相对剂量差异结果列于表2。5例患者靶区D95%Dmean的累积剂量学差异均在-0.05%~0.09%之间;危及器官累积相对剂量差异在-0.37%~0.51%之间,其中以病例2结果最为显著,脑干差异为0.51%,脊髓差异为-0.37%,左、右眼晶状体差异分别为0.04%、0.41%,剂量分布差异均无统计学意义(P>0.05)。危及器官的累积相对剂量差异比靶区稍大,但总体上差异都在±1%以内,剂量分布差异无统计学意义(P>0.05)。

表 2 5例患者靶区与危及器官的累积相对剂量差异(%)

3. 单次与累积相对剂量差异的比较:靶区与危及器官的单次和累积相对剂量差异结果列于表3。4个单次的靶区D95%Dmean相对剂量差异在-0.06%~0.10%之间,而累积剂量学差异均在-0.01%~0.06%之间,且剂量差异有统计学意义(z=-2.02~4.61,P < 0.05);4个单次的危及器官Dmax剂量差异在-0.25%~0.25%之间,而累积剂量学差异在-0.02%~0.11%之间,除左眼晶状体外,其他危及器官的Dmax差异有统计学意义(z=4.46、-4.51、2.07,P < 0.05)。单次计划相对剂量差异有正有负,累积计划在一定程度上减小了单次MLC叶片位置偏差对剂量的影响。

表 3 靶区与危及器官的单次和累积相对剂量差异(%)

4. DVH结果:某例患者两种计划DVH差异如图3所示。DVH图一致性很好,靶区较危及器官DVH差异稍大,但总体上每个结构的DVH改变<0.8%,原计划与实际计划的剂量分布基本吻合,说明叶片位置偏差对剂量的影响甚小。

图3 两种计划的DVH图差异(重建计划与原计划之差) 注:PGTVnx.计划原发灶靶区;PGTVnd.计划淋巴结靶区;PTV1.计划靶区1;PTV2.计划靶区2
讨 论

Luo等[7]的研究显示,0.20 mm的平均位置偏差将导致1%剂量改变。Mu等[10]分别模拟了随机误差和系统误差对简单和复杂鼻咽癌计划的影响,表明随机误差对复杂计划(子野数>50)靶区的影响在1%左右,对危及器官的影响稍大但平均不超过2%,而系统误差的影响非常明显。本研究采用了动态调强放疗技术,子野(或控制点)数目均远>50,同时,通过实际治疗的日志文件得到的0.16 mm整体平均叶片位置偏差,来源于随机误差和系统误差的叠加,而由此导致的靶区累积剂量改变<0.1%,危及器官累积剂量改变也<1%,低于单次计划的相对剂量差异,也远低于其他研究,甚至比单纯的随机误差引起的剂量改变还小[7, 10, 11]。因此,与单次治疗相比,较多分次的治疗确实可以在一定程度上减小叶片位置偏差的影响。本研究未发现叶片位置偏差的角度依赖性,说明重力作用影响很小,可以忽略。结果显示,危及器官累积相对剂量变化比靶区大,可能原因是,危及器官一般位于剂量跌落区,小野照射或窄条野照射的次数较多,叶片走位复杂,导致出现走位误差的次数增加,更多的叶片偏差对剂量分布的影响也会更大[12];在相同的网格分辨率条件下,[JP+1]危及器官体积较小,微小的叶片位置偏差可能带来较明显的剂量差异。同时发现,A、B两组的实际位置相对计划位置有一定的延迟,且A组更明显,说明这是一种系统误差,可能由于控制系统的延迟或每天开机时的MLC刻度产生。这种误差造成了实际叶片间隙的减小,但剂量差异可忽略。而在用Argus 软件解析日志文件时却发现每次治疗中总会出现某个别射野的实际控制点数比计划控制点数多1个,且最后一个控制点的MU都显示为0,说明在执行该控制点之前,已经停止了出束,结合这一发现,可以推断,本研究系统误差来源于控制系统的延迟的可能性非常高。一般来说,系统的延迟是控制系统本身固有的,很难消除,但从本研究结果来看,其影响是没有临床意义的。同时,为了避免严重的MLC叶片位置偏差,本研究中,加速器叶片允许偏差设置为2 mm,统计的1 400个日志文件中8‰的叶片位置偏差超过Klein等[13]建议的1 mm,但这对剂量分布造成的影响却可忽略不计。加速器叶片允许偏差不宜太小,过小容易造成连锁停止出束,反而容易引起剂量的改变,建议叶片偏差设定在1.50~2.00 mm之间。

值得注意的是,调强放疗实际照射的每个子野的MU也有可能变化。Luo等[7]的研究发现动态调强治疗中MU会有微小的改变,但其对剂量的影响极小,可以忽略。在本研究中并未发现MU的改变,因此,并没有涉及MU改变对剂量的影响。此外,动态调强放疗中,叶片运动速度是剂量率的函数,叶片的延迟会导致剂量率的改变,但在日志文件和记录验证系统中,没有记录每个子野(控制点)的实际剂量率。而本研究与之前的相关研究在剂量重建时,假设实际剂量率和计划剂量率保持不变,这可能给剂量重建的准确性带来一定的影响,影响的程度尚无法估计。但可通过观察叶片位置偏差的大小来定性评估剂量率改变的程度,本研究结果表明,叶片位置偏差的范围非常小,说明剂量率的改变也非常小。因此,如果要把基于日志文件的动态调强验证作为一种日常的验证手段,就要更加注意MLC本身的日常质控和维护保养,尽可能减小叶片位置偏差。基于此,建议该方法不应成为一个放疗中心调强验证的唯一手段,而应结合其他测量方法,保证调强放疗的安全性。

综上所述,基于动态调强日志文件的多分次累积剂量验证在一定程度上减小了单次MLC叶片位置偏差对剂量的影响,验证结果更加精确地反映了患者的实际受照剂量。

参考文献
[1] Rangaraj D, Zhu M, Yang D, et al. Catching errors with patient-specific pretreatment machine log file analysis[J]. Pract Radiat Oncol, 2013, 3(2): 80-90.
[2] 白瞳,朱健,尹勇. 静态调强放疗中多叶光栅到位精度数字化分析[J]. 中华放射医学与防护杂志, 2014, 34(11): 837-840.
[3] 余晖,张书旭. 日志文件在静态调强计划MLC到位精度检测中的应用[J]. 中国医疗设备, 2013, 28(1): 24-26.
[4] 张国前,张书旭,韩鹏慧. 动态调强放射治疗滑窗技术的剂量和机械稳定性评价[J]. 中国医疗设备,2011,26(10):21-24.
[5] 马阳光,张可,胡志辉,等. 利用加速器日志文件验证容积旋转调强放疗计划的可行性[J]. 中华放射医学与防护杂志, 2012, 32(3): 285-288.
[6] 石锦平,滕建建,谢秋英,等. 直线加速器日志文件对容积调强三维剂量验证研究[J]. 中华肿瘤防治杂志, 2014, 21(16): 1270-1273.
[7] Luo W, Li J, Price RA Jr, et al. Monte Carlo based IMRT dose verification using MLC log files and R/V output[J]. Med Phys, 2006, 33(7): 2557-2564.
[8] Bai S, Li G, Wang M, et al. Effect of MLC leaf position, collimator rotation angle, and gantry rotation angle errors on intensity-modulated radiotherapy plans for nasopharyngeal carcinoma[J]. Med Dosim, 2012, 38(2): 143-147.
[9] Siebers JV, Keall PJ, Wu Q, et al. Effect of patient setup errors on simulataneously integrated boost head and neck IMRT treament plans[J]. Int J Radiat Oncol Biol Phys, 2005, 63(2): 422-433.
[10] Mu G, Ludlum E, Xia P. Impact of MLC leaf position errors on simple and complex IMRT plans for head and neck cancer[J]. Phys Med Biol, 2008, 53(1): 77-88.
[11] Rangel A, Dunscombe P. Tolerance on MLC leaf position accuracy for IMRT delivery with a dynamic MLC[J]. Med Phys, 2009, 36(7): 3304-3309.
[12] 陈诚,周晓艺,韩光,等.静态调强计划中多叶准直器叶片到位误差的检测及其对剂量分布的影响[J]. 中华放射医学与防护杂志, 2015, 35(3): 210-213.
[13] Klein EE, Hanley J, Bayouth J, et al. Task Group 142 report: quality assurance of medical accelerators[J]. Med Phys, 2009, 39(6): 4197-4212.