[论文解读] Analysis of time-profiles with in-beam PET monitoring in charged particle therapy

A. Kraan, S. Muraro+15 more|arXiv (Cornell University)|Sep 28, 2018

Radiation Therapy and Dosimetry参考文献 19被引用 3

一句话总结

本研究通过分析质子照射后短时间内正电子发射核素（11C、15O、10C）的时间轮廓，评估了质子治疗中在射野内PET监测的性能。结合FLUKA蒙特卡罗模拟与130 MeV质子照射的幻影实验数据，结果表明测量值与模拟值在空间分布和衰变速率上具有良好一致性，但实验数据中10C的含量显著高于蒙特卡罗预测，尤其在高密度聚乙烯幻影中更为明显，突显出当前核反应模型在短寿命核素模拟方面的潜在局限性。

ABSTRACT

Background: Treatment verification with PET imaging in charged particle therapy is conventionally done by comparing measurements of spatial distributions with Monte Carlo (MC) predictions. However, decay curves can provide additional independent information about the treatment and the irradiated tissue. Most studies performed so far focus on long time intervals. Here we investigate the reliability of MC predictions of space and time (decay rate) profiles shortly after irradiation, and we show how the decay rates can give an indication about the elements of which the phantom is made up. Methods and Materials: Various phantoms were irradiated in clinical and near-clinical conditions at the Cyclotron Centre of the Bronowice proton therapy centre. PET data were acquired with a planar 16x16 cm$^2$ PET system. MC simulations of particle interactions and photon propagation in the phantoms were performed using the FLUKA code. The analysis included a comparison between experimental data and MC simulations of space and time profiles, as well as a fitting procedure to obtain the various isotope contributions in the phantoms. Results and conclusions: There was a good agreement between data and MC predictions in 1-dimensional space and decay rate distributions. The fractions of $^{11}$C, $^{15}$O and $^{10}$C that were obtained by fitting the decay rates with multiple simple exponentials generally agreed well with the MC expectations. We found a small excess of $^{10}$C in data compared to what was predicted in MC, which was clear especially in the PE phantom.

研究动机与目标

• 评估FLUKA蒙特卡罗模拟在预测130 MeV质子照射后短时间内组织中正电子发射核素（11C、15O、10C）的空间与时间衰变轮廓方面的可靠性。
• 通过多指数拟合从实验衰变曲线中提取相对核素丰度，实现对组织成分的推断。
• 研究实验测量与模拟结果中10C贡献的差异，特别是在富含碳材料（如高密度聚乙烯）中的表现。
• 展示利用短时间间隔PET数据在空间与时间上绘制不同核素分布（15O、11C、10C）的可行性，为未来的生物与剂量学应用提供支持。

提出的方法

• 使用130 MeV质子笔形束（10^10个质子，5秒照射）照射三种均质幻影（PMMA、高密度聚乙烯、水）及一个非均质Zebra幻影（PMMA/聚乙烯层状结构）。
• 采用基于LYSO的平面16×16 cm²探测系统（DoPET）获取时间分辨PET数据，并应用一种非传统的直接重建方法（SRA）实现空间与时间事件的同时重建。
• 利用FLUKA模拟粒子输运与光子传播，预测空间与时间上的活度分布，包括各核素特异的衰变贡献。
• 对照射后8–300秒的实验衰变曲线进行三指数函数拟合，以提取15O（t₁/₂=2 min）、11C（t₁/₂=20 min）和10C（t₁/₂=19 s）的相对丰度。
• 对比不同幻影中测量与模拟的一维z方向轮廓及衰变速率，并通过2 mm z方向分层方式生成各核素贡献的空间分布图。
• 验证了被忽略的核素（如5B、14O、13N）在选定时间窗口内对总活度的贡献均低于2%。

实验结果

研究问题

• RQ1 FLUKA模拟在预测130 MeV质子照射后幻影中11C、15O和10C的空间与时间衰变轮廓方面准确性如何？
• RQ2 照射后前5分钟内测量的衰变速率在多大程度上反映了蒙特卡罗模拟所预期的核素丰度？
• RQ3 为何实验数据中10C的含量持续高于FLUKA预测，特别是在高密度聚乙烯中？
• RQ4 时间分辨PET数据能否用于在非均质幻影中绘制不同正电子发射核素（15O、11C、10C）的空间分布？
• RQ5 观测到的10C差异对在射野内PET中组织成分分析与生物建模有何影响？

主要发现

• 测量值与FLUKA预测的一维z方向活度轮廓之间表现出良好一致性，Zebra幻影的分层结构因PMMA与PE中15O生成差异而清晰可辨。
• 实验数据中10C的相对丰度始终高于FLUKA模拟结果，高密度聚乙烯幻影中分别为20.7%（实验）与14.3%（MC）。
• 在水和PMMA幻影中，10C丰度在实验数据中也有所升高（分别为1.6% vs. 2.4% MC 和 8.8% vs. 6.6% MC），但程度低于聚乙烯幻影。
• 15O的丰度被FLUKA准确预测（水中为89.4%实验值 vs. 91.5% MC），证实了对富氧组织相互作用的建模准确性。
• 核素的空间分布图显示明显差异：15O在PMMA区域丰富但在PE中缺失，而10C分布均匀但衰变更迅速，与其中短半衰期一致。
• 实验数据中10C的分布图与MC模拟存在差异，尤其在幻影深层更为显著，提示在致密材料中短寿命核素生成或输运的建模可能存在局限性。