[论文解读] Modelling beam transport and biological effectiveness to develop treatment planning for ion beam radiotherapy

L. Grzanka|arXiv (Cornell University)|Oct 1, 2014

Radiation Therapy and Dosimetry参考文献 1被引用 5

一句话总结

本博士论文通过整合粒子输运模拟与细胞存活的生物物理模型，开发了一套全面的碳离子放射治疗计划计算框架。该研究引入了用于轨道结构建模的 libamtrack 库，实现了相对生物效应（RBE）的 Katz 模型，并展示了适用于临床应用的优化深度-剂量分布与细胞存活曲线。

ABSTRACT

Radiation therapy with carbon ions is a novel technique of cancer radiotherapy, applicable in particular to treating radioresistant tumours at difficult localisations. Therapy planning, where the medical physicist, following the medical prescription, finds the optimum distribution of cancer cells to be inactivated by their irradiation over the tumour volume, is a basic procedure of cancer radiotherapy. The main difficulty encountered in therapy planning for ion radiotherapy is to correctly account for the enhanced radiobiological effectiveness of ions in the Spread Out Bragg Peak (SOBP) region over the tumour volume. In this case, unlike in conventional radiotherapy with photon beams, achieving a uniform dose distribution over the tumour volume does not imply achieving uniform cancer cell inactivation. In this thesis, an algorithm of the basic element (kernel) of a treatment planning system (TPS) for carbon ion therapy is developed. The algorithm consists of a radiobiological part which suitably corrects for the enhanced biological effect of ion irradiation of cancer cells, and of a physical beam transport model. In the radiobiological component, Katz's track structure model of cellular survival is applied, after validating its physical assumptions and improving some aspects of this model. The Katz model offers fast and accurate predictions of cell survival in mixed fields of the primary carbon ions and of their secondary fragments. The physical beam model was based on available tabularized data, prepared earlier by Monte Carlo simulations. Both components of the developed TPS kernel are combined within an optimization tool, allowing the entrance energy-fluence spectra of the carbon ion beam to be selected in order to achieve a pre-assumed uniform (flat) depth-survival profile over the SOBP region, assuring uniform cancer cell inactivation over the tumour depth.

研究动机与目标

• 开发一种基于物理的模型，用于描述粒子束在组织样介质中的输运行为。
• 整合径向能量沉积分布与线性能量转移（LET）计算，以实现精确的微剂量学建模。
• 实现并验证适用于混合场辐射暴露的 Katz 细胞存活模型。
• 优化碳离子束在治疗计划中的深度-剂量分布与细胞存活曲线。
• 提供一个计算基础设施（libamtrack），以实现高精度、高效率且可扩展的重离子治疗计划。

提出的方法

• 利用蒙特卡罗模拟生成碳离子束的能量通量谱。
• 开发并实现了 libamtrack 软件库，用于轨道结构计算，包括电子能量-射程关系与径向剂量分布（D(r)）公式。
• 应用 Katz 模型，采用积分与缩放两种版本，基于靶标尺寸与辐射品质计算细胞失活概率。
• 采用插值算法处理能量谱数据文件（SPC 格式），以实现高效的剂量与 RBE 计算。
• 通过自定义算法优化深度-剂量分布，实现均匀、临床适用的剂量分布。
• 将模型与实验数据进行对比，并与成熟的模型（如局部效应模型 LEM）进行基准测试。

实验结果

研究问题

• RQ1 在水与组织等效介质中，碳离子的径向剂量分布与 LET 值能否被准确建模？
• RQ2 Katz 模型的实现能在多大程度上准确预测混合场照射场景下的细胞存活？
• RQ3 能否通过束流输运模拟与剂量优化算法生成优化的深度-剂量分布？
• RQ4 Katz 模型的积分版本与缩放版本在预测生物效应方面有何差异？
• RQ5 libamtrack 库能否作为离子束放射治疗治疗计划的稳健且可扩展的基础？

主要发现

• libamtrack 库成功实现了电子能量-射程关系与径向剂量分布公式，其结果与实验数据高度一致。
• libamtrack 中的径向剂量分布（D(r)）公式在 1–100 nm 径向范围内与实验测量结果保持一致。
• Katz 模型的缩放版本在预测混合场照射下细胞存活方面优于积分版本，尤其在低剂量区域表现更优。
• 通过优化算法实现了目标体积内平坦的剂量分布，证明了其临床可行性。
• 基于模型的 RBE 计算结果强烈依赖于 LET 与细胞类型，不同细胞系中碳离子的 RBE 值在 2.0 至 4.0 之间。
• 该框架通过结合蒙特卡罗束流输运与生物物理建模，实现了高效、高精度的治疗计划，其结果经 published 实验与模拟数据验证。