中华放射医学与防护杂志  2019, Vol. 39 Issue (3): 197-201   PDF    
呼吸运动状态对动态调强放疗剂量分布影响的研究
吴先想1 , 牛振洋2 , 费振乐1 , 刘苓苓3 , 崔相利3 , 魏敏1 , 王磊2     
1. 安徽医科大学生命科学学院生物医学工程系, 合肥 230032;
2. 中国人民解放军联勤保障部队第901医院放疗科, 合肥 230031;
3. 中国科学院合肥肿瘤医院, 合肥 230031
[摘要] 目的 探讨不同幅度、周期、方向的呼吸运动对动态调强放疗(IMRT)计划中靶区剂量分布的影响。方法 选取30例肺癌病例,按靶区体积大小分为A(72.0~200.2 cm3)、B(271.7~380.0 cm3)、C(498.9~684.9 cm3)3组,每组10例,平均体积分别为151.5、327.1和583.3 cm3。使用呼吸运动模拟平台带动含二维电离室矩阵的模体沿枪靶方向运动。分别转动准直器至0°和90°,在不同呼吸运动幅度(0、4、8、12和15 mm)与周期(3、4和5 s)下,采集模体等中心层面剂量。其中周期为4 s测量5次,以绝对剂量及γ通过率(3 mm/3%)为指标,分析采集剂量与治疗计划系统(TPS)输出的剂量分布差异。结果 在两个方向上,呼吸运动降低了靶区边缘内侧剂量,提高了靶区边缘外侧剂量。呼吸运动周期之间的γ通过率差异最大达3.54%(t=2.301,P < 0.05)。当呼吸运动幅度超过8 mm时,γ通过率 < 90%,且随幅度增大而减小。静态与呼吸运动之间γ通过率的差值和靶区体积呈负相关,A、B、C 3组的平均γ通过率依次增大。5次叠加剂量的γ通过率高于单次剂量平均γ通过率,且差异有统计学意义(t=-9.36~-5.95,P < 0.05)。结论 动态IMRT靶区剂量分布主要受呼吸运动幅度及自身体积影响,部分幅度下呼吸运动周期对剂量分布有影响。多次剂量实施后,可消除部分单次剂量实施误差。医师需要根据呼吸运动幅度对靶区进行合理外扩,同时优化呼吸运动方向上靶区边缘组织受量。对于靶区体积过小以及呼吸运动幅度过大的患者,应采取呼吸管理技术提高靶区剂量实施的精准性。
[关键词] 呼吸运动     靶区体积     γ通过率     剂量分布     动态调强放疗    
The effect of respiratory motion states on dynamic IMRT dose distribution
Wu Xianxiang1, Niu Zhenyang2, Fei Zhenle1, Liu Lingling3, Cui Xiangli3, Wei Min1, Wang Lei2     
1. Department of Biomedical Engineering, College of Life Sciences, Anhui Medical University, Hefei 230032, China;
2. Department of Radiation Oncology, No. 901 Hospital of PLA Joint Logistics Support, Hefei 230031, China;
3. Cancer Hospital, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
Fund programs: Social Benefit Technology Application Research and Project Connexion of Anhui Province(1704f0804051)
Corresponding author: Fei Zhenle, Email:feizhenle@21.cn.com
[Abstract] Objective To investigate the effect of respiratory movement of different amplitude, period and direction on the dose distribution of target area in dynamic intensity modulated radiation therapy. Methods A total of 30 cases of lung cancer were selected and divided into three groups according to the volume size of the target area, including groups A (72.0-200.2 cm3), B (271.7-380.0 cm3) and C (498.9-684.9 cm3). The average volume was 151.5, 327.1 and 583.3 cm3, respectively. Breathing motion simulation platform was used to drive the mode body with two-dimensional ionization chamber matrix along the Gun-Target direction, then turn the collimator to 0° and 90°, respectively. The doses were collected at the central level in different amplitudes of 0, 4, 8, 12 and 15 mm, periodic respiratory movement at the intervals of 3, 4 and 5 s and respiratory motion measurement with a cycle of 4 s 5 times. The difference of dose distribution between the collected dose and TPS output was analyzed by taking the absolute dose and γ-passing rate (3 mm/3%) as indicators. Results In the two-sided upward, respiratory movement reduced the dose at the medial edge of the target area and increased the dose at the lateral edge of the target area. The difference of γ-passing rate between respiration cycle was up to 3.54% (t=2.301, P < 0.05), and when the respiration movement was more than 8 mm, the γ-passing rate was less than 90% and decreased with the increase of amplitude. The difference of γ-passing rate between static and respiratory motion was negatively correlated with the volume of target area, and the average γ-passing rate of A, B and C three groups increased gradually. The γ-passing rate of 5 composited dose was higher than that of single dose, and the difference was statistically sigificant(t=-9.36--5.95, P < 0.05). Conclusions The dose distribution of dynamic IMRT target area is mainly influenced by respiration range and its own volume, and the respiration cycle has an effect on dose distribution under partial amplitude. After implementing the multiple doses, some single dose implementation errors can be eliminated. Physicians need to expand the target area reasonably according to the range of respiratory movement, and optimize the amount of marginal tissue in the target area in the direction of respiratory movement. For patients with small target volume and large respiratory movement, respiratory management technology should be adopted to improve the accuracy of target dose implementation.
[Key words] Respiratory motor     Target volume     γ-passing rate     Dose distribution     Dynamic intensitymodulated radiation therapy    

在胸腹部肿瘤中,肿瘤的位置随患者呼吸而变化。根据TG76号报告[1],隔膜附近肿瘤的位置受呼吸影响较大,肺下叶肿瘤、肝部肿瘤和胰腺肿瘤运动在头脚方向上幅度最大。Shimizu等[2]认为肺部肿瘤呼吸运动幅度在头脚、左右、前后方向分别可达15.9、10.0和14.6 mm; Ozhasoglu和Murphy[3]认为静息状态下呼吸频率在16~20次/min。动态调强放疗(IMRT)中,照射野分为有限个子野,通过控制单个子野的跳数,提高靶区处方剂量的适形度和均匀性[4]。在照射过程中,呼吸导致肿瘤位置改变,并且相邻子野切换时,照射野中叶片连续运动,每个子野照射范围不确度增大,呼吸运动和不断变化的照射野之间的相互作用将导致剂量分布发生改变。Duan等[5]研究认为,在单个照射野中,呼吸运动导致的点剂量误差最大可达47.8%,多个射野剂量融合后的总体剂量差异最大3.5%。本研究基于呼吸运动模体,模拟患者在实际照射过程中呼吸运动和照射野之间的相互作用,通过对比分析实际剂量分布和治疗计划系统(TPS)输出剂量分布差异,研究呼吸运动对剂量分布的影响。

资料与方法

1.病例资料:选取2016年11月至2018年6月期间在联勤保障部队第901医院放疗中心接受治疗的肺癌患者30例,年龄49~82岁(中位年龄62岁),中央型18例,边缘型12例。根据患者靶区体积大小分为A、B、C 3组,每组10例,靶区平均体积分别为151.5 cm3(72.0~200.2 cm3)、327.1 cm3(271.7~380.0 cm3)、583.3 cm3(498.9~684.9 cm3)。

2.计划设计:使用热塑膜和放疗定位板固定患者体位,采用荷兰飞利浦大孔径CT模拟定位机对患者进行CT扫描。将扫描的CT图像传至瑞典医科达Monaco 5.11系统进行三维重建,医师根据国际辐射单位与测量委员会(ICRU) 83号报告[6]勾画靶区临床靶区(CTV),外扩0.5~1.0 cm获得计划靶区(PTV)。设计5野动态IMRT计划,射线类型选择6 MV X射线。PTV处方剂量为60 Gy,共30次,单次剂量2 Gy,靶区及敏感组织所受剂量均满足临床要求。将设计好的IMRT计划移植到验证模体上,创建机架归零的质量保证(QA)计划,输出等中心层面的冠状位剂量分布。

3.数据采集:将QUASAR可程控模拟呼吸运动平台放置在治疗床板上,其运动方向沿枪靶方向,准直器角度为0°,将验证模体置于该平台上,初始运动位置保持在0位,调整治疗床,使德国PTW公司探测器阵列(Octavius 729)的有效测量点位于等中心处。按照文献[7]方法,测量前,探测器阵列需预热(>6 Gy)。使用平台编辑正弦波模拟呼吸运动,其中幅度分别为0(静态)、4、8、12和15 mm,周期分别为3、4和5 s。其中周期为4 s时,需测量5次,并将剂量做叠加。将采集的剂量分布与TPS输出剂量做γ分析,标准为3 mm/3%,通过率>90%为剂量验证通过[8-10]。将30例患者计划的静态γ通过率与不同周期、幅度状态下的呼吸运动γ通过率做差值,得到Wijk。将准直器角度更改为90°(模拟平行于叶片运动方向的呼吸运动),重复上述步骤。

$ {W_{ijk}} = {S_i} - {M_{ijk}} $ (1)

式中,Si为体积为i的患者计划静态γ通过率;Mijk为体积为i的患者在j幅度和周期k下的γ通过率。

4.统计学处理:数据用x±s表示。采用SPSS 19.0软件进行分析,在A,B,C 3组中,对同幅度不同周期的3组数据分别进行Shapiro-Wilk正态性检验数据服从正态分布。若服从正态分布则进行配对t检验,否则进行Wilcoxon符号秩和检验。P < 0.05为差异有统计学意义。将Wijk与靶区体积做Pearson相关性分析,r>0.6为强相关,r=0.4~0.6为中度相关,r=0~0.4为弱相关,r < 0为负相关。

结果

1.呼吸运动周期与γ通过率之间的关系:呼吸运动垂直于叶片运动方向时,A组中,呼吸运动幅度为4 mm时,周期3和4 s的γ通过率之间差异有统计学意义(Z=-2.06,P < 0.05);呼吸运动幅度为8 mm时,周期3和4 s的γ通过率之间差异有统计学意义(Z=-1.99,P < 0.05)。C组中,呼吸运动幅度为4 mm时,周期3和5 s的γ通过率之间差异为0.39%,有统计学意义(t=2.49,P < 0.05)。呼吸运动平行于叶片运动方向,A组中,呼吸运动幅度为12 mm时,周期5 s和周期3、4 s的γ通过率之间相对差异有统计学意义(t=2.506、2.301,P < 0.05)。结果见表 12

表 1 3组肺癌患者治疗计划在不同幅度、不同周期、垂直于叶片方向的呼吸运动条件下γ通过率比较(%,x±s) Table 1 Comparison of γ-passing rate in three groups of treatment plans under different respiratory motions with different amplitudes, different periods, and perpendicular to the direction of the leaves (%, x±s)

表 2 3组肺癌患者治疗计划在不同幅度、不同周期、平行于叶片方向的呼吸运动条件下γ通过率比较(%,x±s) Table 2 Comparison of γ-passing rate in three groups of treatment plans under different respiratory motions with different amplitudes, different periods, and parallel to the direction of the leaves (%, x±s)

2.呼吸运动幅度和γ通过率的关系:由表 1表 2可见,两种呼吸运动方向下,A、B、C 3组中,呼吸运动幅度为4、8 mm时,平均γ通过率缓慢降低,但仍>90%;幅度超过8 mm时,γ通过率将显著降低;幅度为12 mm时,平均γ通过率 < 85%;幅度达15 mm时,平均 < 70%。

3.靶区体积与γ通过率的关系:Wijk和两方向的呼吸运动幅度全呈负相关, r < 0;垂直于叶片运动方向,幅度为12 mm时,Wijk和靶区体积呈中度负相关,r=-0.44;幅度为15 mm时,Wijk和靶区体积呈强负相关,r=-0.66~-0.64。平行于叶片运动方向,幅度为4 mm时,Wijk和靶区体积呈中度负相关,r=-0.47~-0.44;其余多数呈弱负相关。由表 1表 2可见,随着A、B、C 3组靶区平均体积的增大,γ通过率逐渐升高。

4.剂量叠加与分次剂量平均γ通过率的关系:由表 3可见,在30例患者计划中,呼吸运动周期为4 s时,在两种方向下的呼吸运动中,5次剂量叠加的γ通过率将高于5次剂量的平均γ通过率,但差异无统计学意义(P>0.05)。

表 3 30例肺癌患者治疗计划在不同幅度、不同方向的呼吸运动上5次测量剂量的平均和叠加γ通过率比较(%) Table 3 Comparison of the average of the five measured doses and the composited γ-passing rate of 30 treatment plans in respiratory motions of different amplitudes and directions (%)

5.呼吸运动对γ分析的影响:对30例计划的剂量分布进行γ分析发现,两种呼吸运动方向下,γ分析未通过的点中,热点分布于沿呼吸运动方向的靶区边缘外侧,冷点分布于沿呼吸运动方向的靶区边缘内侧,剂量梯度较大处容易产生剂量实施误差,靶区内侧高剂量区域几乎不受影响,并且冷点与热点随着呼吸运动幅度增大而增加。

讨论

呼吸运动导致实际照射剂量与计划输出剂量之间产生差异,主要是由呼吸运动的模糊效应和肿瘤运动与叶片运动之间的相互影响效应所导致[11-14]。模糊效应在肿瘤的运动方向上对剂量分布产生影响,这种效应通常会给靶区边缘带来过量照射。在动态IMRT中,多叶光栅在对射束调制过程中,叶片在高剂量区运动缓慢,而在低剂量区快速运动,实现高剂量梯度的目的,然而该调制过程加剧了剂量分布的不确定性。在临床中,由于患者在实际治疗过程中,呼吸幅度和周期都是在变化之中,通常对靶区进行外扩或采用慢速CT、4DCT技术来确定肿瘤运动范围[1],只要整个疗程中不出现明显的脱靶,靶区处方剂量就不会受到影响,这样就避免或降低了模糊效应及相互影响效应对靶区剂量分布的影响。

呼吸运动幅度对靶区剂量影响较大,且主要影响呼吸运动方向上靶区边缘的剂量分布,这与费振乐[13]、Xie等[15]研究结果相似,但其研究只考虑单方向的呼吸运动,本次研究加入了靶区体积、呼吸运动周期、多分次剂量叠加以及平行于叶片运动方向呼吸运动对靶区剂量分布的影响。本次研究发现,当呼吸运动垂直于叶片运动时,呼吸周期对剂量分布的影响较小。呼吸运动平行于叶片运动方向时,相互影响效应增大,不同周期间γ通过率差异较高,同时呼吸周期与出束时间相比非常短[15],多数周期之间差异无统计学意义。受模糊效应的影响,随着呼吸幅度的增大,偏移出照射野的靶区体积占靶区总体积的百分比增大,靶区外进入照射野的体积增多,靶区边缘内侧冷点和外侧热点数量升高。Wijk和靶区体积大小呈负相关,即随着靶区体积的增大,剂量分布差异缩小。由于呼吸运动方向上靶区边缘长度存在差异,研究中两个方向上的γ通过率与靶区体积相关性大小并无一致性。当呼吸幅度超过8 mm时,两方向上的呼吸运动平均γ通过率将 < 90%,剂量验证将不能通过,患者的处方剂量受到严重影响。多次照射后,叠加剂量的γ通过率整体高于单次剂量的平均γ通过率。单次照射后,靶区内剂量分布有微小的差异,但多次照射后,该差异会被弥补,靶区剂量因呼吸运动而被均匀化,γ通过率将升高。

临床中,通常采用对靶区外扩的方法防止靶区脱靶,但更多的正常组织将进入照射野,危及器官剂量将会升高。本次研究将为肺部肿瘤治疗提供临床参考意见。靶区外扩时,医师不但要考虑靶区在不同方向上的呼吸运动幅度[16],还需要考虑靶区体积的大小。

总之,通过分析不同呼吸运动的状态对动态IMRT剂量分布的影响,为肺部肿瘤的放射治疗提供个性化指导,尽量降低呼吸运动带来的剂量实施误差,提高剂量实施的精准性。

利益冲突
作者贡献声明 吴先想负责采集数据,论文撰写,对结果进行统计与分析;牛振洋指导实验过程;费振乐设计实验思路,指导论文的撰写和修改;刘苓苓、崔相利、魏敏、王磊负责计划设计
参考文献
[1]
Keall PJ, Mageras GS, Balter JM, et al. The management of respiratory motion in radiation oncology report of AAPM Task Group 76[J]. Med Phys, 2006, 33(10): 3874-3900. DOI:10.1118/1.2349696
[2]
Shimizu S, Shirato H, Ogura S, et al. Detection of lung tumor movement in real-time tumor-tracking radiotherapy[J]. Int J Radiat Oncol Biol Phys, 2001, 51(2): 304-310. DOI:10.1016/S0360-3016(01)01641-8
[3]
Ozhasoglu C, Murphy MJ. Issues in respiratory motion compensation during external-beam radiotherapy[J]. Int J Radiat Oncol Biol Phys, 2002, 52(5): 1389-1399. DOI:10.1016/S0360-3016(01)02789-4
[4]
Ezzell GA, Galvin JM, Low D, et al. Guidance document on delivery, treatment planning, and clinical implementation of IMRT:report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee[J]. Med Phys, 2003, 30(8): 2089-2115. DOI:10.1118/1.1591194
[5]
Duan J, Shen S, Fiveash JB, et al. Dosimetric and radiobiological impact of dose fractionation on respiratory motion induced IMRT delivery errors:a volumetric dose measurement study[J]. Med Phys, 2006, 33(5): 1380-1387. DOI:10.1118/1.2192908
[6]
Hodapp N. The ICRU report 83:prescribing, recording and reporting photon-beam intensity-modulated radiation therapy (IMRT)[J]. Strahlenther Onkol, 2012, 188(1): 97-99. DOI:10.1007/s00066-011-0015-x
[7]
Li JG, Yan G, Liu C. Comparison of two commercial detector arrays for IMRT quality assurance[J]. J Appl Clin Med Phys, 2009, 10(2): 2942.
[8]
牛振洋, 费振乐, 段宗锦, 等. 探测器阵列在TPS模型调整中的应用[J]. 中华放射肿瘤学杂志, 2017, 26(7): 806-809.
Niu ZY, Fei ZL, Duan ZJ, et al. Application of detector array in treatment planning system modling adjustment[J]. Chin J Radiat Oncol, 2017, 26(7): 806-809. DOI:10.3760/cma.j.issn.1004-4221.2017.07.018
[9]
Low DA, Dempsey JF. Evaluation of the gamma dose distribution comparison method[J]. Med Phys, 2003, 30(9): 2455-2464. DOI:10.1118/1.1598711
[10]
Low DA, Moran JM, Dempsey JF, et al. Dosimetry tools and techniques for IMRT[J]. Med Phys, 2011, 38(3): 1313-1338. DOI:10.1118/1.3514120
[11]
Court L, Wagar M, Berbeco R, et al. Evaluation of the interplay effect when using RapidArc to treat targets moving in the craniocaudal or right-left direction[J]. Med Phys, 2010, 37(1): 4. DOI:10.1118/1.3263614
[12]
Chan C, Lang S, Rowbottom C, et al. Intensity-modulated radiotherapy for lung cancer:current status and future developments[J]. J Thorac Oncol, 2014, 9(11): 1598-1608. DOI:10.1097/JTO.0000000000000346
[13]
费振乐. 呼吸运动幅度对静态IMRT剂量分布的影响研究[J]. 中华放射肿瘤学杂志, 2016, 25(5): 508-512.
Fei ZL. Effect of amplitude of respiratory movements on dose distribution of static intensity-modulated radiotherapy[J]. Chin J Radiat Oncol, 2016, 25(5): 508-512. DOI:10.3760/cma.j.issn.1004-4221.2016.05.019
[14]
费振乐, 牛振洋, 李志杰, 等. 模拟呼吸运动对三维放疗剂量分布的影响[J]. 中华放射肿瘤学杂志, 2012, 21(6): 554-556.
Fei ZL, Niu ZY, Li ZJ, et al. Simulation of respiratory motion effect on dose distributions in three dimensional radiotherapy[J]. Chin J Radiat Oncol, 2012, 21(6): 554-556. DOI:10.3760/cma.j.issn.1004-4221.2012.06.021
[15]
Xie K, Sun H, Lin T, et al. IMRT dose verification considering passing rate and respiratory motion[J]. Oncol Lett, 2018, 16(1): 963-969. DOI:10.3892/ol.2018.8724
[16]
田玲玲, 毛荣虎, 李定杰, 等. 呼吸运动对肺癌立体定向放疗靶区边缘剂量的影响[J]. 中华放射肿瘤学杂志, 2017, 26(12): 1426-1429.
Tian LL, Mao RH, Li DJ, et al. Effect of respiratory motion on radiation dose to the margin of target volume in stereotactic body radiotherapy for lung cancer[J]. Chin J Radiat Oncol, 2017, 26(12): 1426-1429. DOI:10.3760/cma.j.issn.1004-4221.2017.12.015