目的 应用蒙特卡罗方法建立6711型125I放射性粒子模型,计算其在不同骨介质中的剂量分布,为进一步分析125I放射性粒子治疗头颈部肿瘤的剂量分布提供依据。方法 应用MCNP4C软件(截面库DLC-200)建立6711型125I放射性粒子模型,计算其剂量率常数及在水介质中的径向剂量函数值,并与美国医学物理学家协会(American Association of Physicists in Medicine,AAPM)2004年发布的TG43UI报告推荐值对比验证125I模型。利用该模型,采用ICRU44报告推荐的成人不同骨介质的元素组成及密度编写MCNP材料卡,计算125I模型在不同骨介质中的剂量分布。结果 应用MCNP4C软件建立的125I放射性粒子模型的几何形态及规格符合AAPM的TG43UI要求。模拟得到了125I模型剂量率常数Λ=0.977 78 cGy/(h·U)及在水介质中的径向剂量函数值,并得出125I模型在不同骨介质中的径向剂量分布。在同样的介质深度,125I粒子的剂量沉积能力自骨皮质、骨松质、软骨、水、黄骨髓、红骨髓依次递减。在介质深度分别为0.01、0.1、1 cm处,125I粒子在骨皮质中的剂量沉积分别是在水中的12.90倍、9.72倍、0.30倍。结论 本研究建立的125I放射性粒子模型可靠性好,可用于后续在骨介质中的模拟计算。125I放射性粒子在不同的骨介质中的剂量分布不等同于水,能量主要沉积在骨介质表面,因此在制定靶区毗邻骨组织的125I粒子植入计划时应考虑骨介质的影响。
Objective: To present a theoretical analysis of how the presence of bone in interstitial brachytherapy affects dose rate distributions with MCNP4C Monte Carlo code and to prepare for the next clinical study on the dose distribution of interstitial brachytherapy in head and neck neoplasm.Methods: Type 6711,125I brachytherapy source was simulated with MCNP4C Monte Carlo code whose cross section library was DLC-200. The dose distribution along the transverse axis in water and dose constant were compared with the American Association of Physicists in Medicine (AAPM) TG43UI update dosimetry formalism and current literature. The validated computer code was then applied to simple homogeneous bone tissue model to determine the affected different bone tissue had on dose distribution from125I interstitial implant.Results:125I brachytherapy source simulated with MCNP4C Monte Carlo code met the requirements of TG43UI report. Dose rate constant, 0.977 78 cGy/(h·U), was in agreement within 1.32% compared with the recommended value of TG43UI. There was a good agreement between TG43UI about the dosimetric parameters at distances of 1 to 10 cm along the transverse axis of the125I source established by MCNP4C and current published data. And the dose distribution of125I photon emitting source in different bone tissue was calculated. Dose-deposition capacity of photons was in decreasing order: cortical bone, spongy bone, cartilage, yellow bone marrow, red bone marrow in the same medium depth. Photons deposited significantly in traversal axis among the phantom material of cortical bone and sponge bone relevant to the dose to water. In the medium depth of 0.01 cm, 0.1 cm, and 1 cm, the dose in the cortical bone was 12.90 times, 9.72 times, and 0.30 times of water respectively.Conclusion: This study build a125I source model with MCNP4C Monte Carlo code, which is validated, and could be used in subsequent study. Dose distribution of photons in different bone medium is not the same as water, and its main energy deposits in bone medium surface, so we should consider the effect of bone medium when we design the target area adjacent to the bone tissue in125I sources implantation plan.
组织间植入放射性粒子近距离治疗肿瘤具有局部控制率高、微创、高度适形及毒副反应小等优势[1], 在前列腺、头颈部、肺、肝、胰腺等部位的实体肿瘤治疗中已取得良好的临床疗效[2, 3]。
目前, 在临床上进行的剂量计算中, 均是把人体组织模拟为水模体, 但与软组织相比较, 光子在骨介质中的剂量沉积与水介质有明显差异[4], 125I放射性粒子所释放的γ 射线即为光子。在前列腺、肺、肝等部位的肿瘤治疗中, 由于靶区均为软组织, 且距离骨组织较远, 所以骨组织对其剂量分布的影响较小, 而当125I放射性粒子治疗头颈部肿瘤时, 靶区内常存在或毗邻大量骨组织, 因此在制定治疗计划、设计处方剂量、疗效评估时, 明确125I放射性粒子在骨组织内剂量分布的影响对进一步提高疗效、同时降低毒副反应有着重要意义。
蒙特卡罗(Monte Carlo)模拟也称统计试验方法, 是一种基于概率论解决物理和数学问题的统计学方法, 可有效处理不同介质中的光子输运问题及不同组织中密度和原子序数的变化对于剂量沉积的影响, 其数值结果已经成为比较其他理论方法的非试验金标准[5, 6]。
本研究通过蒙特卡罗模拟方法建立125I粒子剂量场模型, 并利用该模型计算125I放射性粒子在人体骨介质(骨皮质、骨松质、黄骨髓、红骨髓以及软骨)中的剂量分布情况, 为临床中对邻近骨组织的靶区进行125I放射性粒子布源及剂量计算提供参考依据。
美国医学物理学家协会(American Association of Physicists in Medicine, AAPM)于1995年发布了指导近距离治疗剂量计算的TG-43报告, 随着近距离治疗的不断发展, 于2004年发布了其修订版TG43UI, 2007年再次对其进行了补充。
本研究的剂量计算所用公式[7]如下:
公式(1)和(2)所用的极坐标系统如图1。在极坐标系统下, Λ 为剂量率常数, 指的是参考点处的剂量率, 单位cGy/(h· U)。
模拟粒子源为125I放射性粒子(6711型, 宁波君安药业科技有限公司), 半衰期为59.40 d, 发射的光子能量分别为27.202、27.472、30.980、31.710、35.492 keV, 对应的每次衰变发出的光子数分别为0.406 0、0.757 0、0.202 0、0.043 9、0.066 8。采用蒙特卡罗MCNP4C程序, 在极坐标系统下, 建立125I放射性粒子蒙特卡罗模型。放射性粒子中介质的物理参数如表1所示。
依据公式(1), 在垂直于粒子源中心轴方向1 cm处取两个体积为0.04 cm3的栅元, 一个用于填充空气材料, 另一个填充水材料, 使用* F4卡和质量能量吸收系数计算沉积在两个栅元中的剂量, 两者的比值即为剂量率常数[8]。
将125I放射性粒子源置于半径为30 cm的液态水球中心, 在源中垂轴的正、负方向上各布置1组计算点, 每组包含16个计算点, 分别距源中心0.10、0.15、0.25、0.50、0.75、1.00、1.50、2.00、3.00、4.00、5.00、6.00、7.00、8.00、9.00和10.00 cm。计算上述这些位置上的吸收剂量, 将正、负方向上距离相同两点处的数值进行平均后作为该距离处的吸收剂量值。依据公式(2)求得各距离上径向剂量函数值, 单次计算时间为250~550 min, 模拟的光子数为108~109个, 将所得结果与TG43UI推荐数值及部分国内已发表的数据进行比较, 以验证本次模拟建模的可靠性。
建立单一125I粒子源对不同骨介质的照射模型。以半径为30 cm建立真空球形模型, 125I放射性粒子置于模型中心, γ 射线产生的光子在真空运行中不发生衰减和吸收, 直至撞击骨面。骨组织规格:5 cm× 5 cm× 4 cm (厚), 距粒子源1 cm。取较大边界条件, 保证光子在传输过程中基本不会逃逸出组织, 得到光子沿不同厚度骨介质的能量密度分布。光子入射骨组织后随机发生吸收和散射反应, 直到能量被吸收低于阈值或溢出组织。记录每个光子被吸收的能量和位置, 大量光子在组织中随机反应的统计结果就形成有界骨组织中的剂量分布。采用ICRU44报告推荐的成人骨皮质、骨松质、黄骨髓、红骨髓以及软骨的元素及密度组成[9](表2), 以同样的模型, 将骨组织的材料卡转为水介质的材料卡进行计算, 比较光子在骨组织和液态水中剂量分布的差异。模拟光子数为108~109个, 分析其在不同骨介质中的剂量分布。
MCNP4C所建立的6711型125I放射性粒子模型如图2所示, 最外层为钛壳, 钛壳与内部的银棒以空气相隔, 125I核素分布在银棒表面。
模拟计算得到125I模型在水介质中的剂量率常数Λ =0.977 78 cGy/(h· U)。125I在水介质中的径向剂量函数计算值与TG43UI报告理论值[7]以及孙亮等[10]的报道值的对比结果见表3。
125I放射性粒子在不同骨介质和液态水中不同距离的单光子剂量沉积见图3, 反映了125I放射性粒子剂量场下γ 射线在不同介质中不同深度的能量沉积情况。图中横坐标为光子进入介质的深度, 纵坐标为相对剂量分布, 所有介质内探测单元统计性误差均控制在5%以内。结果表明, 在不同介质中, 125I粒子辐射场的能量分布存在明显差异。同样介质深度下, 光子的沉积能力自骨皮质、骨松质、软骨、水、黄骨髓、红骨髓依次递减。125I剂量场下, 光子在进入介质表面极浅表的位置时主要发生光电效应, 光子转化为电子在介质中迅速沉积, 随着进入介质深度的增加, 其能量沉积呈指数衰减趋势, 主要发生康普顿散射(Compton scattering)。
骨介质与水介质在相同深度的剂量沉积差异较大, 在介质深度0.01 cm、0.1 cm、1 cm处, 125I粒子在骨皮质中的剂量沉积分别是在水中的12.90倍、9.72倍、0.30倍。由实验数据可推算出, 当单个活度为0.8 mCi的125I放射性粒子植入骨皮质或水中, 在距离125I放射性粒子为0.01 cm、0.1 cm、1 cm处时, 骨皮质所沉积的剂量为110.2 Gy、66.9 Gy、0.4 Gy, 水介质所沉积的剂量为16.4 Gy、13.2 Gy、2.8 Gy。
125I放射性粒子近距离治疗是治疗肿瘤的重要手段, 剂量的精确分布是其达到安全性和有效性的重要保障, 也是当前放射物理学面临的一个重要课题。临床上广泛应用的治疗计划系统(treatment planning system, TPS)将人体骨组织、软组织等效为液态水, 对不同组织间的介质差异未进行设计。在前列腺、肝、肺、胰腺等部位的近距离治疗中, 因靶区为软组织且毗邻骨组织较远, 其差异尚可接受, 但头颈部含有大量的骨组织, 相对肌肉、腺体等组织对电离辐射的吸收、散射的影响明显与水介质有区别[11, 12]。因此, 在对头颈部肿瘤靶区的术前设计中如何评估骨组织对剂量分布的影响, 进而优化术前设计, 是125I放射性粒子治疗头颈部肿瘤亟待解决的问题。
125I放射性粒子在骨介质中的剂量分布问题尚未见文献报道。由于人体解剖结构复杂, 个体差异大, 加之125I放射性粒子释放γ 射线能量较低, 仅为 27.202~35.492 keV, 采用热释光等剂量元件对125I粒子在骨介质中的能量沉积情况进行实验测定受到诸多限制, 故本实验选择蒙特卡罗模拟方法进行研究。
本研究中建立的125I放射性粒子蒙特卡罗模型, 其几何形态和规格符合TG43UI推荐值。模拟计算得出的剂量率常数与TG43UI推荐值[0.965 cGy/(h· U)]相差1.32%, 在水介质中的径向剂量函数数值非常接近TG43UI推荐值[7]及孙亮等[10]的研究结果(表3), 表明本研究所建立的125I放射性粒子蒙特卡罗模型可靠, 可用于后续研究125I粒子源在不同介质中的剂量分布。Belley等[4]的研究表明, 同等介质中, 光子源的平均能量决定了其在介质中的穿透和能量沉积能力, 不同组织介质中(骨、肌肉和水)光子的沉积能力存在明显差异, 且这种差异在光子能量为0.01~0.1 MeV的骨介质与水介质中明显[1], 并且可推测125I粒子在骨组织和水介质中剂量分布差异明显。本研究中125I粒子源在不同介质中剂量沉积的数据表明, 同样的介质深度, 125I放射性粒子所释放的γ 射线的沉积能力自骨皮质、骨松质、软骨、水、黄骨髓、红骨髓依次递减。125I放射性粒子能量沉积作用在骨介质表面处最为明显, 且骨皮质、骨松质的能量沉积迅速, 很快出现光子能量沉积的峰值, 随穿透深度增大, 与其他介质能量沉积的差异持续存在。临床治疗过程中, 随着植入粒子数量的增加, 这种差异的累积效应也将更明显。
本研究应用蒙特卡罗模拟建模的方式, 建立的125I放射性粒子模型较为精确地模拟了125I粒子源在水介质中的径向剂量分布, 模型可靠性好, 可用于后续其他组织的模拟计算; 采用不同的元素组成和密度模拟不同的介质, 得到实验难以测量的125I放射性粒子在骨介质中的剂量分布结果; 初步分析了125I放射性粒子在骨介质及水介质中剂量分布的差异, 可帮助临床医生对涉及骨组织靶区的术前设计进行优化; 同时, 本研究建立的模型及得到的相关剂量参数, 为后续研究125I放射性粒子治疗不同部位肿瘤的剂量分布提供了基础。
本研究在进行骨介质计算时, 其元素组成及密度采用ICRU44报告推荐值, 且为均质性模型, 而在临床应用中, 不同部位骨的厚度、组成不同, 多为多种骨介质混杂, 且元素分布、密度分布并不均匀。因此, 125I放射性粒子在人体的不同部位骨组织内的剂量分布还需要进一步研究。随着近距离治疗的不断发展, 对剂量计算的精度要求越来越高, 125I放射性粒子在骨介质中剂量分布的研究发展前景广阔。
(志谢:感谢北京大学口腔医院颌面外科的范祎老师对本研究的大力支持, 使研究得以顺利开展, 特此致谢!)
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|