中国辐射卫生  2021, Vol. 30 Issue (1): 94-98  DOI: 10.13491/j.issn.1004-714X.2021.01.021

引用本文 

孔栋, 魏贤顶, 惠琳, 孔旭东, 赵于天, 杨波. Bolus下空腔对各向异性分析算法计算精度影响[J]. 中国辐射卫生, 2021, 30(1): 94-98. DOI: 10.13491/j.issn.1004-714X.2021.01.021.
KONG Dong, WEI Xianding, HUI Lin, KONG Xudong, ZHAO Yutian, YANG Bo. Effect of cavity under Bolus to shallow dose calculation precision of anisotropic analytical algorithm[J]. Chinese Journal of Radiological Health, 2021, 30(1): 94-98. DOI: 10.13491/j.issn.1004-714X.2021.01.021.

通讯作者

杨波,E-mail:wuxiyangbo@163.com

文章历史

收稿日期:2020-08-27
Bolus下空腔对各向异性分析算法计算精度影响
孔栋 , 魏贤顶 , 惠琳 , 孔旭东 , 赵于天 , 杨波     
江南大学附属医院 肿瘤放疗科,江苏 无锡 214062
摘要目的 基于蒙特卡罗算法研究Bolus下不同大小空腔对各向异性分析算法(AAA)计算浅层组织剂量精度的影响。方法 在Eclipse中构建上表面位于加速器源轴距、中心与射野中心轴重合的30 cm × 30 cm × 30 cm水模体,其上表面构建不含或含有不同大小空腔的1 cm厚Bolus,采用AAA分别计算不同空腔下水模体中心轴和不同深度x轴上剂量分布;用Geant4构建加速器模型和相同水模体及补偿膜,验证加速器模型正确性,并计算水模体相同位置剂量分布;以Geant4结果为基准,对比不同空腔下AAA计算精度变化。结果 对于1 cm厚Bolus,空腔面积为2 cm × 2 cm时,其厚度 < 5 mm,对AAA精度影响在2%左右,空腔厚度进一步增加,AAA会过高评估浅表剂量,随深度增加对AAA影响迅速减小;空腔面积增加,对AAA计算中心区域影响逐渐减小直至基本无影响,影响较大区域逐渐外移,空腔面积 ≥ 射野时,对AAA影响不再变化,影响较大区域移至射野边缘附近。结论 计划制定时对于相较治疗子野面积较小的空腔,若其厚度 ≥ 5 mm应适当提高空腔下浅表剂量;若空腔面积较治疗子野过大或空腔厚度过大,则需重新定位。
关键词补偿膜    空腔    各向异性分析算法    Geant4    剂量    
Effect of cavity under Bolus to shallow dose calculation precision of anisotropic analytical algorithm
KONG Dong , WEI Xianding , HUI Lin , KONG Xudong , ZHAO Yutian , YANG Bo     
Department of Radiation Oncology, Affiliated Hospital of Jiangnan University, Wuxi 214062 China
Abstract: Objective To research the effect of cavity under Bolus to anisotropic analytical algorithm (AAA) on calculation precision of dose in shallow tissue based on Monte Carlo method; Methods A 30 cm × 30 cm × 30 cm water phantom with the upper surface was constructed which was located at the source-axis distance (SAD) of the medical linear accelerator and the center as well as coincided with the central axis of the radiation field in Eclipse treatment panning system. Above the water phantom, a water film of 1 cm thick with or without different cavities was constructed or. AAA was used to calculate the dose distribution on the central axis and the x-axis of different depth of the water model with different cavities respectively. The accelerator model, the same water phantom and the water film were constructed and the dose distributions of the same positions were calculated with Geant4. Based on the Geant4 calculation result, the calculation precision of AAA with different cavity were compared; Results For cavities with area of 2 cm × 2 cm, if the thickness is smaller than 0.5 cm, the AAA calculation error is about 2%. with the cavity thickness increase, the AAA would overestimate the dose in the shallow area under the cavity. With the cavity area increase, the area where AAA overestimate the shallow dose gradually moved out until near the edge of the radiation field, and the calculation error on the central area gradually reduced until there is basically no error. Conclusions The shallow dose would be increased according to the cavity size when planning with AAA; If there are cavities with large volume, it is better to reposition.
Key words: Bolus    Cavity    Anisotropic Analytical Algorithm (AAA)    Geant4    Dose    

高能X射线存在剂量建成区,治疗浅表肿瘤时需在皮肤表面加Bolus,由于人体体表凹凸不平,通常会在Bolus下形成一些空腔,并反映在定位CT上。Bolus下空腔对放疗中浅表剂量和百分深度剂量的影响多有研究[1-4],而这些空腔对治疗计划系统(TPS)影响的研究较少。各向异性分析算法(AAA)是Eclipse TPS中的一种剂量算法,组织的非均质性对该算法计算精度影响较大,尤其是当空腔存在时[5-7],虽然Eclipse之后又推出了新的算法[8-10],但一些医疗机构目前依然在用AAA。蒙特卡罗(MC)算法是剂量计算的金标准,在放疗及辐射防护中多有应用[6, 7, 11, 12]。本研究基于Geant4 MC软件包[13]计算结果对比了不同大小空腔对AAA计算浅表剂量的影响,为计划制定和评估提供参考,对Bolus应用有一定指导意义。

1 材料与方法 1.1 材料

Visual Studio 2013,Geant4蒙特卡罗软件包(v10_02_p02),瓦里安2300CD医用直线加速器,IBA蓝水箱剂量测量系统,Eclipse治疗计划系统(v8.6),OriginPro 8。

1.2 治疗计划系统(TPS)模型构建和数据获取

在Eclipse中构建层厚0.1 cm、大小30 cm × 30 cm × 30 cm、上表面位于加速器源轴距且与射线入射方向垂直的水模体,CT值指定为0。根据研究目的,在水模体上表面分别设置3组补偿膜,依次为对照组、厚度组和面积组。对照组(图1A)为紧贴水模体上表面、大小30 cm × 30 cm、厚 1 cm的补偿膜;厚度组(图1B)在对照组基础上对补偿膜进行变形,在补偿膜下形成面积2.0 cm × 2.0 cm、厚度依次为0.2 cm、0.3 cm、0.5 cm、0.7 cm和0.9 cm的空腔;面积组与厚度组类似,在补偿膜下形成厚度为0.9 cm,面积依次为2 cm × 2 cm、4 cm × 4 cm、8 cm × 8 cm、10 cm × 10 cm、12 cm × 12 cm和15 cm × 15 cm的空腔。空腔部分CT值指定为−1000,其他部分CT值指定为0。BODY和计算区域包含水模体、补偿膜和空腔,计算网格为0.1 cm × 0.1 cm × 0.1 cm,入射野面积为10 cm × 10 cm,跳数100 MU,采用AAA计算模体中剂量分布。以0.2 cm × 0.2 cm分辨率导出通过射野中心轴的横断面剂量分布。

图 1 空腔模型 Figure 1 Cavity models 注:A对照组模型;B厚度组和面积组模型。
1.3 加速器蒙特卡罗模型构建和相空间文件获取

本研究所用服务器CPU为Intel Xeon E5-2630,16核16线程,操作系统为64位win7,C++编译器为Visual Studio 2013,Geant4版本号10_02_p02。基于Varian公司提供模型,以Geant4构建能量6 MV、射野面积10 cm × 10 cm的2300CD医用直线加速器机头模型,算法物理模型为G4EmPenelopePhysics,截止射程1 mm,计算水中剂量分布,通过与IBA蓝水箱测量数据对比验证模型的准确性(图2)。在Geant4中分3组模拟,每组初始粒子1 × 108,每次产生30个相同粒子,记录到达Jaw上粒子信息作为相空间文件1,粒子数约1.36 × 108,该过程历时约208 h;以相空间文件1作为粒子源,依次读取其中粒子信息,每次产生10个相同粒子,记录到达Jaw下粒子信息作为相空间文件2,粒子数约5.4 × 107,该过程历时约15 h。

图 2 Geant4计算与测量数据对比 Figure 2 Comparison of Calculated and measured data 注:A.百分深度剂量分布;B.10 cm深度离轴比。
1.4 Geant4中不同空腔模型剂量计算

构建与TPS中相同的3组模型,并在水模体中心轴及0.1 cm、0.3 cm、0.5 cm、0.9 cm、1.5 cm深度处X方向上(−8 cm,8 cm)范围设置0.2 cm × 0.2 cm × 0.2 cm的体素以记录剂量沉积。每组模型均分3组进行计算,以相空间文件2作为粒子源,依次读取粒子信息,每次产生50个相同粒子,分别获取体素中剂量沉积,该过程历时约164 h。每组射野范围内剂量相对均方差小于2%。

1.5 数据处理及绘图

Geant4计算数据和TPS导出数据均保存为ASCII文件,将其拷贝到Excel中进行统计分析。为进一步减小模型带来的误差,本研究采用相对深度剂量比(relative depth dose ratio,RDDR)(式(1))和相对离轴剂量比(relative off-axis dose ratio,ROADR)(式(2))描述不同空腔对中心轴和不同深度处侧向profile带来的影响,其越接近100%说明影响越小,大于100%说明AAA计算值高于Geant4,小于100%说明AAA计算值低于Geant4。图片均由OriginPro 8软件绘制。

$RDDR = {{\left( {{{{D_{C,d,AAA}}} / {{D_{d,AAA}}}}} \right)} / {\left( {{{{D_{C,d,MC}}} / {{D_{d,MC}}}}} \right)}} \times 100\% $ (1)

式中, ${D_{C,d,AAA}}$ ${D_{d,AAA}}$ 分别为含空腔和不含空腔时模体深度d处AAA计算值; ${D_{C,d,MC}}$ ${D_{d,MC}}$ 分别为含空腔和不含空腔时模体深度d处Geant4计算值。

$ROADR = {{\left( {{{{D_{C,p,AAA}}} / {{D_{p,AAA}}}}} \right)}/ {\left( {{{{D_{C,p,MC}}} / {{D_{p,MC}}}}} \right)}} \times 100\% $ (2)

式中, ${D_{C,p,AAA}}$ ${D_{p,AAA}}$ 分别为含空腔和不含空腔时模体p处AAA计算值; ${D_{C,p,MC}}$ ${D_{p,MC}}$ 分别为含空腔和不含空腔时模体p处Geant4计算值。

2 结 果 2.1 对照组

无空腔时,Geant4与AAA计算百分深度剂量(PDD)差异在2%以内(图3A),0.1 cm深度处离轴比差异在射野范围(−5 cm,+5 cm)也均在2%以内(图3B)。因此无空腔存在时两种算法计算结果基本吻合。考虑到无Bolus时,Geant4计算结果与测量数据吻合很好(图2),因此该模型可用于评估Bolus下空腔对AAA计算精度的影响。

图 3 对照组Geant4与AAA计算结果对比 Figure 3 Calculation result of Geant4 and AAA in control group 注:A.百分深度剂量;B.0.1 cm深度离轴比。
2.2 空腔厚度对AAA计算精度影响

图4图5可知,对于面积为2 cm × 2 cm的空腔,厚度 ≤ 0.5 cm时,对AAA计算精度影响在2%左右,不超过3%;随空腔厚度进一步增加,相对深度剂量比和空腔下0.1 cm深度相对离轴剂量比迅速增加,但这种影响主要表现在浅表层,在0.5 cm深度处即使对于0.9 cm厚的空腔,影响也已降到2%左右。

图 4 不同厚度空腔下相对深度剂量比 Figure 4 Distributions of relative depth dose ratio with different thickness cavities

图 5 不同厚度空腔下相对离轴剂量比分布 Figure 5 Distributions of relative off-axis dose ratio with different thickness cavities 注:A. 0.1 cm深度处;B. 0.5 cm深度处。
2.3 空腔面积对AAA计算精度影响

图6图7可知,对于厚度为0.9 cm的空腔,随面积增加,在中心区域对AAA计算精度影响减小,影响较大区域逐渐外移,直至空腔面积 ≥ 射野时,影响较大区域位于射野内侧约1 cm范围内不再移动。影响仅限于浅表层,在0.5 cm深度处影响最大只有2%左右。

图 6 不同面积空腔下相对深度剂量比 Figure 6 Distributions of relative depth dose ratio with different area cavities
3 讨 论

放疗中Bolus的应用对提高患者表面剂量,保证浅表肿瘤剂量适形性和均匀性极其重要,胡作怀等[14]研究证明定位时在体模内贴皮肤加Bolus相较计划中添加虚拟Bolus更为合理,但定位时添加Bolus的好坏,即Bolus是否紧贴皮肤表面将对后续计划制定及最终治疗效果产生一定影响。

本文基于Geant4计算数据系统研究了Bolus下空腔厚度和面积对AAA计算精度的影响。可以看到,空腔厚度主要影响计算误差大小,尤其是空腔厚度 > 0.5 cm时,浅表剂量计算误差将迅速增加,AAA会明显高估浅表剂量,这与陈华等 [7]研究结果类似,但Bolus下空腔影响更大,这可能是因为本研究中空腔处尚未建立带电粒子平衡,此时又需进行二次建成,而AAA低估了电子不平衡带来的影响;空腔面积主要影响误差的分布,尤其是空腔面积 ≥ 射野面积时,误差主要发生在射野边缘附近浅表层,这可能是因为AAA未充分考虑电子发散的影响。本研究结果可为临床工作提供依据。虽然临床所遇空腔厚度和面积大多都较小,但射野的倾斜入射会变相增加空腔厚度,同时调强放疗中会涉及一系列面积较小的子野,相对这些子野,某些空腔面积还是偏大的。升级TPS算法或许能解决这个问题,但在此之前物理师在使用相关TPS制定放疗计划时应根据空腔大小适当增加空腔下浅表肿瘤剂量,医生在评估计划时也应考虑到这点,尤其是空腔较大时,对肿瘤区域剂量均匀性的过高要求极有可能意味着肿瘤浅表部分剂量的不足。当空腔特别大时,建议重新定位,但每次定位都会给患者带来额外剂量,应尽量避免这种情况的发生。

图 7 不同面积空腔下相对离轴剂量比分布 Figure 7 Distributions of relative off-axis dose ratio with different area cavitie 注:A. 0.1 cm深度处;B. 0.5 cm深度处。

本研究是针对1 cm厚Bolus和目前放疗中应用最多的6 MV X射线,对于其他厚度的Bolus及射线能量,空腔对AAA计算精度的影响可能会有一些变化,同时放疗TPS中有很多算法,空腔对这些算法的影响也各不相同,这有待进一步研究。

参考文献
[1]
Sharma SC, Johnson MW. Surface dose perturbation due to air gap between patient and bolus for electron beams[J]. Med Phys, 1993, 20(2 Pt 1): 377-378. DOI:10.1118/1.597079
[2]
Kong M, Holloway L. An investigation of central axis depth dose distribution perturbation due to an air gap between patient and bolus for electron beams[J]. Australas Phys Eng Sci Med, 2007, 30(2): 111-119. DOI:10.1007/BF03178415
[3]
李承军, 胡健, 张爱华, 等. 6 MV X线时组织等效补偿膜与人体空气间隙对表面剂量的影响[J]. 医疗卫生装备, 2011, 32(3): 83-84.
Li CJ, Hu J, Zhang AH, et al. Surface dose perturbation due to air gap between patient and bolus for 6MV X-ray[J]. Chin Med Equip J, 2011, 32(3): 83-84. DOI:10.3969/j.issn.1003-8868.2011.03.033
[4]
孔栋, 惠琳, 魏贤顶, 等. 补偿膜下空腔间隙对浅表组织剂量影响[J]. 中华放射肿瘤学杂志, 2018, 27(11): 1009-1013.
Kong D, Hui L, Wei XD, et al. Effect of cavity under bolus on the dose of superficial tissue[J]. Chin J Radiat Oncol, 2018, 27(11): 1009-1013. DOI:10.3760/cma.j.issn.1004-4221.2018.11.012
[5]
Codel G, Serin E, Pacaci P, et al. SU-F-SPS-01: accuracy of the small field dosimetry using acuros XB and AAA dose calculation algorithms of eclipse treatment planning system within and beyond heterogeneous media for trubeam 2.0 unit[J]. Medical Physics, 2016, 43(6Part4): 3350. DOI:10.1118/1.4955676
[6]
张彦秋, 邱小平, 杨振. 用蒙特卡罗模拟评估小野条件下肺介质中AAA算法的精度[J]. 中国医学物理学杂志, 2012, 29(1): 3099-3103.
Zhang YQ, Qiu XP, Yang Z. Evaluation the calculation accuracy of AAA algorithm for the situation with small fields in lung by Monte Carlo[J]. Chin J Med Phys, 2012, 29(1): 3099-3103. DOI:10.3969/j.issn.1005-202X.2012.01.005
[7]
陈华, 徐义果, 庄志邈, 等. 放疗计划系统中空腔边缘剂量计算准确性研究[J]. 中华放射肿瘤学杂志, 2017, 26(1): 69-73.
Chen H, Xu YG, Zhuang ZM, et al. Dose accuracy research on air cavity interface in treatment planning system[J]. Chin J Radiat Oncol, 2017, 26(1): 69-73. DOI:10.3760/cma.j.issn.1004-4221.2017.01.016
[8]
Rana S, Rogers K, Pokharel S, et al. Acuros XB algorithm vs. anisotropic analytical algorithm: a dosimetric study using heterogeneous phantom and computed tomography (CT) data sets of esophageal cancer patients[J]. J Cancer Ther, 2013, 4(1): 138-144. DOI:10.4236/jct.2013.41019
[9]
刘致滨, 石锦平, 张利文, 等. Acuros XB与AAA算法在肺癌调强放疗计划设计中的比较[J]. 中国医学物理学杂志, 2017, 34(10): 984-987.
Liu ZB, Shi JP, Zhang LW, et al. Comparison of Acuros XB algorithm and anisotropic analytical algorithm in intensity-modulated radiotherapy plan for lung cancer[J]. Chin J Med Phys, 2017, 34(10): 984-987. DOI:10.3969/j.issn.1005-202X.2017.10.003
[10]
苗利, 邱正帅, 程孝国. 基于Eclipse研究不同人体替代组织的深度剂量特性[J]. 中国辐射卫生, 2018, 27(6): 621-623.
Miao L, Qiu ZS, Cheng XG. The depth dose characters of different tissue substitutes based on Eclipse[J]. Chin J Radiol Heal, 2018, 27(6): 621-623. DOI:10.13491/j.issn.1004-714x.2018.06.029
[11]
Wang L, Ding GX. The accuracy of the out-of-field dose calculations using a model based algorithm in a commercial treatment planning system[J]. Phys Med Biol, 2014, 59(13): N113-N128. DOI:10.1088/0031-9155/59/13/n113
[12]
李东, 刘平, 季芳, 等. 基于蒙卡GEANT4研究某放疗场所的辐射水平分布规律[J]. 中国辐射卫生, 2018, 27(4): 413-416.
Li D, Liu P, Ji F, et al. Research on radiation distribution rule of a radiation therapy place based on GEANT4 Monte Carlo progress[J]. Chin J Radiol Heal, 2018, 27(4): 413-416. DOI:10.13491/j.issn.1004-714x.2018.04.034
[13]
[14]
胡作怀, 付建东, 陈芳, 等. 组织补偿膜添加方式引起的放疗剂量差异研究[J]. 中华放射肿瘤学杂志, 2016, 25(4): 388-390.
Hu ZH, Fu JD, Chen F, et al. Difference in radiotherapy dose caused by different ways of adding bolus[J]. Chin J Radiat Oncol, 2016, 25(4): 388-390. DOI:10.3760/cma.j.issn.1004-4221.2016.04.018