摘要目的提出一种非匀速优化旋转治疗方法，从剂量学角度研究了其可行性。方法包括优化治疗剂量分布计算和实验两部分。在剂量计算中，设计了一个用于合理调节射束旋转速率的优化因子，并用半适形野取代单一野，以利于改善剂量分布。将优化方法加进基本数学模型中，以大量实验数据为基础，设计完成了剂量计算法和相应C程序。结果引入优化方法后，剂量分布得到显著改善。结论本文提出的优化方法能够在一定程度上解决中心横截面上的适形问题；这一初步结论可为下一步研究——采用变频技术控制机架旋转速率，将该方法运用到直线加速器和⁶⁰Co治疗机上提供了必要的剂量学参考。

Theoretical and Experimental Studies of Non-uniform Speed Optimized Rotational Radiotherapy

Ma Shaogang,Bai Sen,Zhang Yiyun,el al.

(The Affiliated Hospital of Medical College of Qingdao University,Qingdao 266003)

Abstract ObjectiveA method of non-uniform speed optimized rotation therapy was proposed,and the feasibility was studied from dosimetrical angle.Methods The study content is divided into two parts:the dose distribution calculation of the optimized therapy and the experiment. The optimization factor designed for rotation speed and quasi-conformal field was used to improve the dose distribution.The two optimizations are added into the basical math-model.An algorithm and the responding C program are designed on the basis of great amount of experimental data.Results The results of the dose calculation pre and post being optimized show that the dose distributions are remarkably improved by the optimization.ConclusionFrom the theoretical calculation and experimental results,by using non-uniform speed rotational irradiation we may partially solve the problem of conformation on the horizontal axis (mid-plane).By using the variable frequency technique to control the rotation speed our technique will be easily applied to accelerator or⁶⁰Co therapy machine.

Key WordsOptimized rotation therapyDosimetryIsodose curveDose distributionFilm dosimeter

本文提出了一种非匀速优化旋转治疗方法，以典型的偏中心肿瘤——肾癌为例，通过理论计算和实验证对该方法的剂量学可行性进行了研究。研究中采用了自行研制的逐行逐列扫描式等剂量线的跟踪和绘制程序。^［1］

1材料和方法

1.1优化旋转治疗剂量计算

本文针对典型的偏中心肿瘤——肾癌的治疗问题，设计了一个用于调节射速旋转速率的优化因子，并采用了半适形野优化方法，来共同改善治疗剂量分布。我们将优化方法融合到基本数学模型当中，以大量实验数据为基础，设计完成了剂量计算算法和相应C程序。

需要说明的是：为了研究方便，暂且采用组织等效的均匀体模，并对体模条件作了适当简化：把穿过人体横断面所形成的体模轮廓和靶区形状近似看成椭圆。该研究只涉及过肿瘤中心的人体横截面(同时又是照射野中心平面)上的剂量计算。另外，为了将计算时间减少到适当程度，又不影响计算结果，我们采用N_beam个固定射束来模拟计算连续旋转射束在体模内形成的剂量分布。

1.1.1优化方法

非匀速优化对于旋转照射，要想改变授予体模内任一点P点的剂量值有两种方法：改变剂量率和改变旋转速度。前者是文献［2，3］的优化方法，后者就是本文的优化方法。本文设计了一个优化因子F_opt，用于合理地调节旋转速度，改善剂量分布。F_opt应表示成以下三个变量的函数：①靶区形状：由旋转中心到靶区连缘的距离r₁描述；②体模轮廓：由旋转中心到体表的距离r₂描述；③组织密度：由有效路径d_eff描述(此项暂不考虑)。

F_opi=F_opt(r₁,r₂,d_eff)

(1)

●半适形野优化在机架旋转过程中的任一瞬间，可令等中心处的射野大小与该处的肿瘤的投影宽宽W_conf一致。肿瘤投影宽度是指，在通过等中心点垂直于射野中心轴的平面上把放射源假设成光源时肿瘤的投影宽度。

1.1.2数学模型

固定射束在体膜内任意点的剂量计算问题是旋转照射剂量计算的基础，下面介绍其数学模型。

本文采用基于实验测量的剂量计算数学模型，它需要大量的测量数据，需对几种射野在一些特定深度处，沿射野中心轴和与它垂直的截面方向作相对剂量测量，然后将测好的数据存入相应表格，任一点的剂量计算可通过查表和在表格中插值完成。

为了描述方便，选用与射束相关的坐标系，原点取射野中心轴在体模的入射点，X轴和Y轴分别平行与方野的两边，Z轴沿着射野中心轴指向体模。

本文所用数学模型的一个重要特点就是，通过分离变量简化计算过程。变量的分离过程是：先定义一个刻度剂量点，然后将剂量分布表达成刻度点剂量的百分量。刻度剂量点取在射野中心轴上的最大剂量建成点。于是，体膜内任意一点P(x,y,d)处的吸收剂量可表达为

D(x,y,z)=D_m·P(x,y,z)

(2)

其中，百分量

P(x,y,z)=PDD(Z,FSZo)·OAR((x,z，FSZo)·OAR(y,z,FSZo)

(3)

式中FSZo是等中心平面上的射野大小。本文只考虑X轴上的情况(y=0)，OAR(y,z,FSZo)=1。接着，考虑到非标称源皮距修正和人体曲面校正，上式还需乘上F因子和曲面校正因子BCC，

D(x,z)=Dm·PDD(d,PSZo)·OAR(x,z,FSZo)·F·BCC

(4)

引入优化算法：在上式中乘上优化因子F_opt，并令FSZo=W_conf,得到完整公式

D(x,z)=Dm·PDD(z,W_conf)·OAR(x,z,W_conf)·F·BCC·F_opt(r₁,r₂,d_eff)

(5)

1.1.3剂量计算算法

剂量计算算法选用了三个坐标系：直角坐标系XOY、直角坐标系X’SY’、极坐标系SX’。在进行剂量计算之前，预先将实验数据表格用自编的数据表存取程序存放好，以备计算程序中查表使用。剂量网格密度取格距0.5cm可满足研究需要^［4］。

图1剂量计算算法流程图

1.2优化旋转治疗模拟实验

1)准备材料：①慢感光胶片剂量计(Kodak X-Omat V)；②混有超细氧化物EVA体模：密度均匀(CT值=-20)，呈椭圆柱状，长轴30cm，短轴20cm；体模分两块，每块厚8cm，以便将胶片夹在其间。

2)做出胶片剂量计的剂量—黑度刻度曲线。

3)照射实验：在暗室中将胶片剂量计夹在两块体模中间，按照体模轮廓将胶片割齐，用闭光胶带沿两块体模间隙密封一圈。将体模放到治疗床上，调整体模位置，保证胶片平面与射束中心轴旋转平面重合，把旋转中心放在体模肿瘤中心位置。本实验用36个固定照射野模拟优化旋转治疗，等中心点的总吸收剂量取700cGy。每个机架角下的剂量权重(对应F_opt)，及野大小(W_conf)都是按照有关算法算出的。

4)胶片的黑度处理：本文在胶片的黑度处理中，使用了计算机多媒体摄像机和图象处理软件photostyler，大大提高了几何分辨率和黑度分辨率。用2)步得出的胶片刻度曲线将实验胶片的黑度分布换算成剂量分布，得图5。

2结果

图2是用现行旋转照射获得的计算结果，即只使用一种射野尺寸FSZo=5cm的匀速旋转照射。使用速率优化因子和半适形野优化后，可得图3所示的计算结果。优化旋转治疗模拟实验结果见图4。

图2现行旋转治疗剂量分布

图3非匀速优化旋转治疗剂量分布

图4非匀速优化旋转治疗模拟实