基于多算法电子剂量场反演方法.pdfVIP

基于多算法的电子剂量场反演方法 1,2 1,2 3 1,2 李贵 郑华庆 吴爱东 吴宜灿 FDS 团队 1 中国科学院等离子体物理研究所 安徽合肥 230031 2 中国科学技术大学核科学技术学院 安徽合肥 230027 3 安徽省立医院放射治科 安徽合肥 230001 摘要 目的：研究精确放射治疗中的电子剂量场反演方法，并研究采用不同的反演算法来实现电子剂量场 反演模型的效果。方法：建立了非线性规划电子剂量场反演模型，提出利用体外透射剂量重建全空间电子 剂量场的方法；并通过包括Levenberg-Marquardt 、Quasi-Newton、Gradient 、Conjugate Gradient 、Newton 、 Principal-Axis 与 NMinimize 等 7 种算法实现。结果：本文设计了 7 个例题进行测试；测试结果表明除了 NMinimize 算法不能实现该模型(只有一算例，误差大)外，其他算法都能实现，反演误差可以小于0.009 （即 0.9% ）; Quasi-Newton 算法表现出很好的特性，精度快速，可以作为电子剂量场反演的首选算法。进一步 测试表明，利用接近电子射程的透射剂量会带来较大的反演误差 （尾差）。结论：对于尾差现象，原因可 能是能谱平均所带来的，因此改进的反演模型需要考虑改进对平均能量的近似。尽管如此，由于反演时间 较短，精度高，可以实现实时剂量场反演；并且每一种算法实现非线性规划反演模型都有特点，这为在实 际实现中提供多方法的选择可能性，可满足临床上对实时剂量验证的需要。 关键词 剂量反演；剂量验证；能谱；放射治疗 中图分类号 R730.55,R811 1 引言 在精确放射治疗中，精确剂量计算如蒙特卡罗模拟等方法可以达到很高的正向计算精 [1-3] [4, 5] 度 ，但是在实际的治疗操作中，由于不可控的误差存在 ，导致照射区域的剂量分布 偏移了计划剂量分布，因此需要实施实时的逆向剂量计算即剂量场反演，以便保证患者体 内肿瘤区域所得剂量符合计划的要求。对于精确放疗的剂量验证，目前的方法基本考虑光 子照射的剂量验证，但是电子照射的剂量验证也是必须的，否则不能保证“精确放疗”的要 求。另外，电子剂量场反演还可以应用于工业探测 （如损伤，材料均匀情况的了解），利 用高能电子射线照射工件，通过测量透射剂量反推工件内部剂量场分布，从而获得工件的 损伤与材料分布等情况。 电子照射的剂量验证的难点有两个：其一电子射程短，超过射程以外基本测不到电子 剂量分布；其二就是光子污染问题，与光子照射不同，光子剂量场可以很好控制电子污染， 而电子剂量场中的光子不能很好地消除。对于第一个问题，目前还没有很好的解决方案， 因为超出电子射程以外电子剂量场的就基本消失，体外测不到电子剂量场的一部分也就难 以反推其全部分布。不过由于电子在高能区如 20MeV( 电子放射治疗的能量在 50MeV 以 [6] 内 ) 以外就能穿透 10cm 以上水模，可以穿透人体一些部位 （如乳腺癌的治疗），这种能 量的电子照射可以实施剂量验证，因此本文把电子剂量场限定在高能区。对于第二个问题， 本文将把光子场也考虑在内，建立一个混合场模型。 本文的研究目的是发展一种新型的电子剂量场的反演方法，在混合笔形束模型的基础 上，建立非线性规划反演模型并利用多种非线性反演方法实现快速精确的电子束剂量场实 [7-10] 时反演，最后通过与蒙特卡罗模拟对比 进行了验证。 [基金项目] 国家“973”计划项目 （2006CB708307 ），安徽省自然科学基金项目(070413081) 。 1 [作者简介] 李贵，博士生，现从事精确放射治疗关键问题研究， Email:gli@ipp.ac.cn Tel:0551-5369742 2 基于混合笔形束的非线性规划反演模型 2.1 剂量场误差分析与反演模