关于X和γ射线束楔形板剂量的修正算法.docxVIP

关于X和γ射线束楔形板剂量的修正算法 关于X和γ射线束楔形板剂量的修正算法 2021年10月 第40卷第5期四川大学学报(自然科学版) Journal of Sichuan University (Natural Science Edition ) Oct. 2021Vol. 40　No. 5文章编号:049026756(2021) 0520901204 关于X 和γ射线束楔形板剂量的修正算法 杨代伦, 罗正明 (四川大学原子核科学技术研究所?辐射物理及技术教育部重点实验室, 成都610064) 摘要:利用加速器治疗肿瘤时, 通常要在射线束的路径上加入楔形板以改变平野的剂量分布. 作者研究了用衰减系数法计算加楔形板后的剂量修正算法. 即通过测量楔形板的材料衰减系 数, 再根据楔形板剖面的尺寸大小, 建立插值算法以计算射线穿透楔形板的路径长度, 从而准 确地计算射线楔衰减因子修正平野的剂量. 关键词:治疗计划系统; 剂量算法; 楔形板剂量修正 中图分类号:R815, TL72　　　文献标识码:A 为了适应临床治疗的需要, 在用加速器治疗肿瘤时[1], 通常要在射线束的路径上加楔形板过滤器以改变平野(即不加楔形板) 的剂量分布. 楔形板是用高密度材料作成的, 它能使射野的输出剂量形成最佳分布, 以达到更好的治疗效果. 传统用的楔形板的楔形角有15°,30°,45°和60°4种. 随着楔形板的广泛运用, 希望楔形角是可变的, 因此实际用的楔形角为至15°至60°间, 而且楔剖面也有变化. 为了计算加楔形板的剂量分布, 一般在加速器出束下, 要实际测量每个楔形板在不同射野的百分深度剂量(PDD ) 和不同深度的离轴剂量(OAR ) . 在放射治疗计划系统[2]中, 通常用测量加和不加楔形板的剂量来计算楔形衰减因子, 再计算出加入楔形板后的剂量分布. 测量楔形板的剂量一般在水模体下进行, 不仅需要测量的数据很多, 而且数据处理、储存等工作量较大. 我们研究了用衰减系数法计算加楔形板后的剂量修正算法. 该方法的特点是:不需要作楔形板的大量PDD 和OAR 数据测量, 仅通过测量在不同能量和射野楔形板的材料衰减系数, 再根据楔形板剖面的尺寸大小, 计算出射线穿透楔形板的路径长度, 并准确地计算射线穿透楔形板的楔形因子, 最后用楔形因子修正平野的剂量. 我们所提出的计算射线穿透楔形板路径的插值算法, 既可适用于标准和非标准形状的楔形板, 也可使计算更简单. 其计算结果与美国AAPM (American Association of Physicists in Medicine ) 55#报告的测量相比较, 符合较好. 1　楔形板剂量修正算法 对于X 射线和γ射线束, 在加楔形板后人体内某一点P (x , y , z ) 处的剂量D w (x , y , z ) 可表示为 (1) D w (x , y , z ) =D (X , Y , Z ) ?f w 式中D (X , Y , Z ) 为平野(即不加楔形板) 的剂量, 可用矩阵模型、Clarkson 或C/S [2～4]方法计算, 作人体表面和组织非均匀修正[5]后求得, f w 为楔形修正因子, 采用衰减系数法有 F w =e -μw d (2) 式中μw 为楔形板材料的衰减系数, d 为在P (x , y , z ) 点的射线穿过楔形板的路径长度(楔形板厚度) . 楔形因子F w 定义为在射野中心轴某一深度d 0点处, 加和不加楔形板的剂量比, 即 F w =D d 0/D d 0w (3) 其中D w d 0, D d 0分别为加和不加楔形板时的剂量, 可通过测量得到. 如果要测量出楔形板材料的衰减系数μw , 则可利用以下关系式求得, 即 收稿日期:2021203224; 修回日期:2021204221 902四川大学学报(自然科学版) μw =-ln F w /d 第40卷(4) 2　楔形修正因子的计算 图1是楔形板示意图, L 为板长, W 为板宽, H 为板厚, σ为楔角(楔角为楔形板的几何角) . 楔形角的定义为:在射束中心轴上的皮表下10cm 深度处等剂量曲线与束轴垂线的夹角. 标准楔形板的楔形角为15°, 30°,45°和60°. 实际使用的楔形板可能有两种:一为标准楔形板, 其楔剖面的形状为三角形; 另一为非标准楔形板, 其楔剖面形状为多边形, 如图2. 传统使用的楔形板为标准楔形板. 随着楔形板用途的拓展, 为使楔形板同时具有挡块的功能,