《一种医用电子直线加速器光子束能谱获取方法.pdf》由会员分享，可在线阅读，更多相关《一种医用电子直线加速器光子束能谱获取方法.pdf（10页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103336295 A (43)申请公布日 2013.10.02 CN 103336295 A *CN103336295A* (21)申请号 201310245158.0 (22)申请日 2013.06.19 G01T 1/36(2006.01) (71)申请人 南京航空航天大学 地址 210016 江苏省南京市白下区御道街 29 号 (72)发明人 周正东 陈元华 (74)专利代理机构 南京经纬专利商标代理有限 公司 32200 代理人 许方 (54) 发明名称 一种医用电子直线加速器光子束能谱获取方 法 (57) 摘要 本发明涉及一种医用电子直线加速器光子束 能谱。
2、获取方法， 通过测量待测加速器发出的光子 束在水模体中产生的百分深度剂量值向量， 基于 光子束能谱的单峰特性构建获取待测光子束能谱 的约束优化重建模型， 并通过对该模型的求解， 获 得该模型的最优解， 即获得待测光子束能谱 ； 本 发明设计的医用电子直线加速器光子束能谱获取 方法， 针对光子束能谱重建模型增加了能谱的约 束条件， 采用全局优化算法进行优化求解， 使得本 方法所获得的解与能谱的固有单峰特性相一致， 从而使得获得的优化解更加稳定， 大大减小了与 实际能谱的偏差， 具有高精度、 高效率的优点， 且 成本低廉、 简便快捷。 (51)Int.Cl. 权利要求书 1 页 说明书 5 页 附。
3、图 3 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书1页 说明书5页 附图3页 (10)申请公布号 CN 103336295 A CN 103336295 A *CN103336295A* 1/1 页 2 1. 一种医用电子直线加速器光子束能谱获取方法， 包括如下步骤 ： 步骤 1. 测量医用电子直线加速器所射出的待测光子束对应其射野中心轴在水模体内 不同深度的百分深度剂量值向量 ； 步骤 2. 将该待测光子束的标称能量 E 平均分成 m 个等间隔能量的能量区域， 取各能量 区域中标称能量的均值为该能量区域的标称能量， 记录为单能光子束 E1、 E2、 Em ； 。
4、步骤 3. 分别模拟计算各个单能光子束对应其射野中心轴在水模体内不同深度的百分 深度剂量值矩阵 ； 其特征在于， 还包括如下步骤 ： 步骤 4. 构建获取待测光子束能谱的约束优化模型如下 ： min F=|WA-D| 其中， Ei、 Ej、 Ep E1、 E2、 Em， 且 Ep为能谱峰值点的标称能量， D 为待测光子束对 应其射野中心轴在水模体内不同深度的百分深度剂量值向量， A 为各个单能光子束对应其 射野中心轴在水模体内不同深度的百分深度剂量值矩阵， W 为标称能量为 Ei的单能光子束 占标称能量为 E 的待测光子束的权重的向量 ； 求出该模型的最优解 W=(W1、 Wi、 Wm)， 即。
5、为待测光子束的能谱。 2. 根据权利要求 1 所述一种医用电子直线加速器光子束能谱获取方法， 其特征在于 ： 所述步骤 1 中， 所述待测光子束垂直射向水模体。 3. 根据权利要求 1 所述一种医用电子直线加速器光子束能谱获取方法， 其特征在于 ： 所述步骤 3 中， 使用蒙特卡罗方法分别模拟计算各个单能光子束对应其射野中心轴在水模 体内不同深度的百分深度剂量值矩阵。 4. 根据权利要求 1 所述一种医用电子直线加速器光子束能谱获取方法， 其特征在于 ： 所述步骤 4 中， 所述获取待测光子束能谱的模型的最优解采用遗传算法或模拟退火算法或 蚁群算法或上述方法的混合算法进行求解。 权 利 要 求。
6、 书 CN 103336295 A 2 1/5 页 3 一种医用电子直线加速器光子束能谱获取方法 技术领域 0001 本发明涉及一种医用电子直线加速器光子束能谱获取方法。 背景技术 0002 肿瘤放射治疗， 是利用电离辐射进行治疗， 主要是通过 射线和 射线治疗肿 瘤。 其中， 对于某些肿瘤只需要放疗， 即可治愈 ； 对于另外一些肿瘤， 无论是手术治疗还是化 疗， 为了降低风险， 在实施过程中往往都需要进行辅助放疗， 放射治疗几乎可用于治疗全身 各个部位的多种癌症， 是治疗癌症的重要手段之一， 随着放疗技术的发展， 其临床应用日益 广泛。 0003 但放射治疗也会带来副作用， 提高肿瘤的控制率。
7、并降低正常组织的并发症概率， 一直是放疗技术研究的焦点。国际辐射单位度量委员会 (International Commission on Radiation Units and Measurements,ICRU)24号与42号报告指出， 临床上对于肿瘤靶区的 剂量给予， 要求剂量不准确度必须小于 5% ； 而高精度的 collapsed cone 剂量计算方法需要 根据医用加速器治疗头所产生的射束能谱进行计算， 只有射束的能谱信息足够精确， 才能 保证剂量计算的精确性， 因此获得医用加速器 X 射线的能谱对于研究放疗患者体内的剂量 分布非常重要。 0004 但是， 加速器治疗头产生的光子束能。
8、量高， 其能谱很难直接测量 ； 根据容易获得的 剂量测量数据， 建立一定的数学模型快速重建出精确的能谱引起了诸多学者的关注。罗正 明等人提出了 “一种测量医用加速器能谱的方法” 的专利， 通过松弛迭代法求解线性方程组 得到加速器光子束的能谱。 但是， 其重建模型没有包含能谱的约束条件， 该方法所获得的解 往往不稳定， 与实际能谱偏差往往较大。 发明内容 0005 本发明所要解决的技术问题是提供一种具有高效率、 高精度的医用电子直线加速 器光子束能谱获取方法。 0006 本发明为了解决上述技术问题采用以下技术方案 ： 本发明设计了一种医用电子直 线加速器光子束能谱获取方法， 包括如下步骤 ： 0。
9、007 步骤 1. 测量医用电子直线加速器所射出的待测光子束对应其射野中心轴在水模 体内不同深度的百分深度剂量值向量 ； 0008 步骤 2. 将该待测光子束的标称能量 E 平均分成 m 个等间隔能量的能量区域， 取各 能量区域中标称能量的均值为该能量区域的标称能量， 记录为单能光子束 E1、 E2、 Em ； 0009 步骤 3. 分别模拟计算各个单能光子束对应其射野中心轴在水模体内不同深度的 百分深度剂量值矩阵 ； 0010 步骤 4. 构建获取待测光子束能谱的约束优化模型如下 ： 0011 min F=|WA-D| 说 明 书 CN 103336295 A 3 2/5 页 4 0012 。
10、0013 其中， Ei、 Ej、 Ep E1、 E2、 Em， 且 Ep为能谱峰值点的标称能量， D 为待测光子 束对应其射野中心轴在水模体内不同深度的百分深度剂量值向量， A 为各个单能光子束对 应其射野中心轴在水模体内不同深度的百分深度剂量值矩阵， W 为标称能量为 Ei的单能光 子束占标称能量为 E 的待测光子束的权重的向量 ； 求出该模型的最优解 W=(W1、 Wi、 Wm)， 即为待测光子束的能谱。 0014 作为本发明的一种优选技术方案 ： 所述步骤 1 中， 所述待测光子束垂直射向水模 体。 0015 作为本发明的一种优选技术方案 ： 所述步骤 3 中， 使用蒙特卡罗方法分别模拟。
11、计 算各个单能光子束对应其射野中心轴在水模体内不同深度的百分深度剂量值矩阵。 0016 作为本发明的一种优选技术方案 ： 所述步骤 4 中， 所述获取待测光子束能谱的模 型的最优解采用遗传算法或模拟退火算法或蚁群算法或上述方法的混合算法进行求解。 0017 本发明所述一种医用电子直线加速器光子束能谱获取方法采用以上技术方案与 现有技术相比， 具有以下技术效果 ： 0018 （1） 本发明设计的医用电子直线加速器光子束能谱获取方法具有高精度、 高效率 的优点， 且成本低廉、 简便快捷 ； 0019 （2） 本发明设计的医用电子直线加速器光子束能谱获取方法中， 针对获得的待测 光子束能谱的特点， 。
12、建立了求解能谱的约束优化模型， 采用遗传算法或模拟退火算法或蚁 群算法或上述方法的组合算法进行求解， 可获得全局最优解， 具有可靠性高的优点。 附图说明 0020 图 1 是本发明设计的医用电子直线加速器光子束能谱获取方法的流程图 ； 0021 图 2 是本发明实际应用过程中待测光子束在水中的百分百深度剂量值数据曲线 图 ； 0022 图 3 是本发明实际应用过程中采用蒙特卡罗方法模拟计算得到 3.125MeV 单能光 子束在水中的百分百深度剂量值数据曲线图 ； 0023 图 4 是本发明实际应用过程中获得 6MV 加速器的光子能谱。 具体实施方式 0024 下面结合说明书附图对本发明的具体实。
13、施方式作进一步详细的说明。 0025 如图 1 所示， 本发明设计了一种医用电子直线加速器光子束能谱获取方法， 包括 如下步骤 ： 0026 步骤 1. 测量医用电子直线加速器所射出的待测光子束对应其射野中心轴在水模 体内不同深度的百分深度剂量值向量 ； 0027 步骤 2. 将该待测光子束的标称能量 E 平均分成 m 个等间隔能量的能量区域， 取各 能量区域中标称能量的均值为该能量区域的标称能量， 记录为单能光子束 E1、 E2、 Em ； 0028 步骤 3. 分别模拟计算各个单能光子束对应其射野中心轴在水模体内不同深度的 百分深度剂量值矩阵 ； 说 明 书 CN 103336295 A 。
14、4 3/5 页 5 0029 步骤 4. 基于能谱的固有单峰特性， 构建获取待测光子束能谱的约束优化模型如 下 ： 0030 min F=|WA-D| 0031 0032 其中， Ei、 Ej、 Ep E1、 E2、 Em， 且 Ep为能谱峰值点的标称能量， D 为待测光子 束对应其射野中心轴在水模体内不同深度的百分深度剂量值向量， A 为各个单能光子束对 应其射野中心轴在水模体内不同深度的百分深度剂量值矩阵， W 为标称能量为 Ei的单能光 子束占标称能量为 E 的待测光子束的权重的向量 ； 求出该模型的最优解 W=(W1、 Wi、 Wm)， 即为待测光子束的能谱。 0033 本发明设计的医。
15、用电子直线加速器光子束能谱获取方法具有高精度、 高效率的优 点， 且成本低廉、 简便快捷。 0034 作为本发明的一种优选技术方案 ： 所述步骤 1 中， 所述待测光子束垂直射向水模 体。 0035 作为本发明的一种优选技术方案 ： 所述步骤 3 中， 使用蒙特卡罗方法分别模拟计 算各个单能光子束对应其射野中心轴在水模体内不同深度的百分深度剂量值矩阵。 0036 作为本发明的一种优选技术方案 ： 所述步骤 4 中， 所述获取待测光子束能谱的模 型的最优解采用遗传算法或模拟退火算法或蚁群算法或上述方法的混合算法进行求解。 0037 本发明设计的医用电子直线加速器光子束能谱获取方法中， 针对获得的。
16、待测光子 束能谱的约束优化模型， 采用遗传算法或模拟退火算法或蚁群算法或上述方法的组合算法 进行求解， 可获得全局最优解， 具有可靠性高的优点。 0038 本发明设计的医用电子直线加速器光子束能谱获取方法在实际应用过程当中， 首 先设定待测光子束的射野和标称能量 E， 其中， 针对医用电子直线加速器进行如下设置， 源 皮距为 100cm， 机架角和机头角均为 0， 调整准直器形成待测大小的射野 ； 将三维水箱置 于加速机头下方， 使机头中轴垂直三维水箱内水模体的表面并通过其几何中心， 待测光子 束由医用电子直线加速器射向水模体 ； 0039 紧接着， 通过设置在水中的剂量计， 分别测量该待测光。
17、子束对应其射野中心轴水 模体内不同深度 d1、 d2、 dn处的百分深度剂量值 P1、 P2、 Pn， 记录为该待测光子束百 分深度剂量值向量 D=（P1、 P2、 Pn） ； 0040 再将该待测光子束的标称能量 E 平均分成 m 个等间隔能量的能量区域， 取各能量 区域中标称能量的均值为该能量区域的标称能量， 记录为单能光子束 E1、 E2、 Em ； 0041 采用蒙特卡罗方法模拟计算标称能量为 E1、 E2、 Em的单能光子束对应其射野中 心轴不同深度 d1、 d2、 dn处的百分深度剂量值矩阵 0042 由于待测光子束由电子在靶中发生轫致辐射而产生， 其能谱连呈现连续的谱分 布， 将。
18、待测光子束能量区间分割的很小， 即 m 较大时， 可认为每个能量区间的标称能量近似 相等 ； 能量为Ei的单能光子束在加速器机头射出待测光子束中所占的权重Wi代表该能量区 说 明 书 CN 103336295 A 5 4/5 页 6 间的能谱， 其中 i 1、 2、 m。 0043 步骤 5. 构建获取待测光子束能谱的约束优化模型如下 ： 0044 min F=|WA-D| 0045 0046 其中， W 为标称能量为 Ei的单能光子束占标称能量为 E 的待测光子束的权重的 向量 ； 求出该模型的最优解 W=(W1、 Wi、 Wm)， 即为待测光子束的能谱， 其中， i 1、 2、 m， Ep。
19、为能谱峰值点的标称能量。其中， 所述获取待测光子束能谱的模型的最优解 采用遗传算法或模拟退火算法或蚁群算法或上述方法的组合算法进行求解。 0047 其中， 采用遗传算法获取待测光子束能谱的模型的最优解， 包括如下步骤 ： 0048 步骤随机初始化能谱个体数目为 K 的群体， 群体内每个能谱个体以实数编码形 式进行表示， 各个能谱个体满足待测光子束能谱的模型， 各个能谱个体分别为 ： 0049 0050 其中， g 1、 2、 G， G 为能谱优化求解问题的变量个数， k 1、 2、 K,ag,bg 为所述待测光子束能谱的模型的解的取值范围 ； 0051 步骤根据 F=|WA-D| 获得群体中各。
20、个能谱个体的适应度值， 以表示进化的好 坏程度， 为下一步骤的选择操作提供了数据和方向 ； 0052 步骤根据各个能谱个体的适应度值， 从父代中选择两个能谱个体， 适应度函数 值越大， 被选中的概率越大 ； 0053 步骤将选择出的两个能谱个体按杂交概率进行杂交， 产生两个新的能谱个体 ； 0054 步骤针对新的能谱个体按变异概率进行变异 ； 0055 步骤接受新的能谱个体并判断是否完成新群体的生成 ； 如果没有， 则返回到步 骤 ； 0056 步骤检验停止演化的条件， 若满足停止条件则停止演化， 停止演化时， 最优的能 谱个体解就是待测光子束的能谱， 其中， 每一次进化过程就产生新一代的群体。
21、， 群体内的每 个能谱个体表示一个可能的能谱解 ； 否则跳转至步骤重新进行演化。 0057 将获得的光子束能谱保存到数据文件， 当进行 collapsed cone 剂量计算时， 将该 光子束能谱数据导入放射治疗计划系统， 即可实现高精度剂量计算的优点。 0058 在实际应用过程当中， 医用电子直线加速器选用 6MV 的西门子加速器 （不仅仅限 于6MV,不仅仅限于西门子加速器） ， 设定好待测光子束的射野和标称能量后， 待测光子束由 6MV 的西门子加速器垂直射向水面， 待测光子束在水中的百分百深度剂量值数据如图 2 所 示 ； 之后， 分别采用蒙特卡罗方法模拟计算0MeV到6MeV的单能光。
22、子束在水中的百分百深度 剂量值数据， 其中， 3.125MeV 单能光子束在水中的百分百深度剂量值数据如图 3 所示 ； 构建 获取待测光子束能谱的模型 ； 采用遗传算法获取待测光子束能谱的模型的最优解， 即获得 待测光子束的能谱如图 4 所示 ； 最后将获得的待测光子束的能谱保存于数据文件中， 以便 该光子束能谱用于放射治疗计划系统中。 0059 本发明设计的医用电子直线加速器光子束能谱获取方法， 针对光子束能谱重建模 型增加了能谱的约束条件， 采用全局优化算法进行优化求解， 使得本方法所获得的解与能 说 明 书 CN 103336295 A 6 5/5 页 7 谱的固有单峰特性相一致， 从而使得获得的优化解更加稳定， 大大减小了与实际能谱的偏 差， 具有高精度、 高效率的优点， 且成本低廉、 简便快捷。 0060 上面结合附图对本发明的实施方式作了详细说明， 但是本发明并不限于上述实施 方式， 在本领域普通技术人员所具备的知识范围内， 还可以在不脱离本发明宗旨的前提下 做出各种变化。 说 明 书 CN 103336295 A 7 1/3 页 8 图 1 说 明 书 附 图 CN 103336295 A 8 2/3 页 9 图 2 说 明 书 附 图 CN 103336295 A 9 3/3 页 10 图 3 图 4 说 明 书 附 图 CN 103336295 A 10 。