近距离粒子植入放射治疗，是将封装有放射源的粒子，通过装置植入到肿瘤内部，对病灶进行持续的照射，达到消融肿瘤的手术。与外照射相比较，此种手术方法能够有效的控制肿瘤，并减少对周围正常组织器官的照射。

临床上，剂量计划的实施是借助徒手穿刺、共面模板和3D非共面模板完成［1-2］。徒手穿刺定位不精确，使得肿瘤内的125I粒子分布情况很难完全符合术前剂量规划的粒子分布，从而产生剂量的冷区和热区。共面模板借助一个带有多排等间距平行孔的模板辅助穿刺针导向定位，并依照术前计划结果布置放射性粒子。这种手术方法简单快速，能够应对大部分体部肿瘤的粒子植入手术，并取得了很好的治疗效果［3-4］。但模板孔的方向和孔的间距固定，面对多病灶、大肿瘤和骨骼或重要器官遮挡情况，其剂量规划方法很难使得剂量分布满足临床要求。有研究表明，使用多个模板能增加粒子植入的灵活性，提高放射治疗效果［5］，是3D非共面模板开发的基础。非共面模板能够有效避开骨骼、重要器官的遮挡，使得剂量有效覆盖整个靶区。研究表明，使用3D非共面模板的手术剂量精准，可达术前剂量设计要求［6］。

传统的粒子植入手术剂量规划方法是基于平面的二维规划［7］。术前剂量规划需要物理医师操作TPS在合适位置放置粒子，通过观察等剂量曲线或曲面对重建靶区的剂量覆盖情况，合理调整粒子的位置和数量，最终满足DVH曲线评估。此过程，即使有丰富经验的医师指导，剂量结果也不一定最优。

本研究提出了一种新颖的基于非共面空间针的自动布源的面向体部肿瘤的剂量规划系统，用于帮助物理医师准确、快速、有效地完成术前剂量规划。

SAP10-2视野检测程序发现视野损害前青光眼黄斑区视功能损伤的能力…………………… 汪星朦 孙兴怀戴毅 等（6）748

1 系统设计理论基础

1.1 空间针的设计

空间针具有位置灵活特性，能从任意角度和位置进行穿刺，使粒子能够放置在靶区的任意位置，为剂量有效的覆盖靶区提供了保证。同时也能躲避靶区周围的重要器官，减低对患者的创伤。空间针的布置方法如图1所示。首先，将针尖放置到靶区，得到在横断面P1内并且与Y轴平行针初始位置；然后，在P1面内旋转得θ1到针的中间位置；最后，将中间针在P2平面内旋转θ2得到针最终位置。

  

图1 空间针在软件中的布置

式中：Dj表示PTV内感兴趣点剂量值，Dp代表处方剂量，ω1是PTV的权重因子；Dk表示心脏内感兴趣点剂量值，Hp是心脏的安全剂量值，ω2是心脏的安全系数的权重因子；D1表示脊髓内感兴趣点剂量值，Sp是脊髓的安全剂量值，ω3是脊髓的安全系数的权重因子。

 

  

图2 空间针针尾的坐标变换

1.2 剂量计算模型

剂量计算模型选取6711型125I放射性粒子，粒子活度选取 0.5 mCi（1 mCi=3.7×107Bq）。基于 TG-43 U1报告［8］的二维线源计算公式，将其扩展到三维剂量计算得到公式（2）。

 

为有效遏制围标，要对参加投标的单位数量进行限制。大型水利工程项目由于资金数额巨大、施工技术复杂，投标单位数量一般不多，而越小的项目，投标单位越多，竞争越激烈，增加了发生多种违规的可能。因此对于中小型水利工程项目，要根据投资金额、招标费用合理确定投标单位数量，减少企图围标者的围标机会。

 

1.3 模拟退火模型

选取骨继发性恶性肿瘤病例为实验材料，此病例已完成临床手术，植入了63颗粒子和消耗22根针。实验仿真选取脊髓为危机器官，使用活度为0.5 mCi的6711型125I粒子，设置处方剂量为120 Gy。DICOM（Digital Imaging and Communications in Medicine）图像总体素为 499×499×170，算法提取的感兴趣区域为 99×99×90。

式中：κB为玻尔兹曼常数；n是温度下降次数；Tn控制退火温度下降的温度参数；α控制退火温度的下降速度，速度下降越慢，模拟退火算法结果越精确，运算时间越长，反之亦然。

 

假设针的总长度为L，针尖为C点，针尾的坐标为T点。针尾的坐标变换可以表示成如图2所示。图中的坐标系定义为 o-xyz、o'uvw、o'-u'v'w'和 o'-u''v''w''。针的位置调整可以等效成坐标系变换。穿刺针首先在 o-xyz坐标系内移动到（x，y，z），后绕 w轴旋转 θ1，最后绕 u'轴旋转 θ2。最终坐标如公式（1）所示。

由图6①剂量规划靶区和重要器官的三维重建可知肿瘤已经侵蚀到了脊髓，所以空间针位置的布置应尽量避免影响脊髓；图6②描述了剂量规划所使用的空间针的分布、三维粒子的位置、等剂量曲面对靶区的覆盖情况以及生成的三维虚拟模板；图6③展示了靶区和脊髓的DVH曲线；图6④是将虚拟模板3D打印出的非共面模板。软件优化算法运行15 s，共使用22根针，消耗63颗放射性粒子。如表1所示，剂量评估参数V100＝90.24，V150＝48.09%以及V200＝14.63%均满足临床剂量标准。由于危机器官被肿瘤侵蚀，最大照射剂量达到96 Gy，超过了安全剂量指标。但其它剂量指标D90＝1.26 Gy，V100＝2.48%和V200＝2.06%，满足剂量要求。

模拟退火算法过程如下：根据Metropolis原则，设当前i状态的状态方程为F（M）i，之后对当前状态进行扰动，然后得到了一个新状态j以及它的状态方程为F（M）j。此时对两个状态方程进行比较，如果 F（M）j＜F（M）i，粒子分布结果将会从状态 i跳转至状态 j；如果 F（M）j≥F（M）i，此时粒子分布状态将会以一定的概率 Pi进行跳转［9］。

 

看来，对方相貌虽凶，到底还是有些忌惮师父的。青辰暗暗松了口气，正要从树后走出，却忽然发现，在对方转身迈出几步之后，那垂落在腰间的左手，掌心似乎有什么东西动了一下。

1.4 放射剂量评估方法

采用剂量体积直方图（dose-volume histogram，DVH）进行治疗计划评估，它描述的是相应体积的靶区与所接受的照射剂量的对应关系。参考美国物理医师协会报告和国内相关文章，设置评估参数如下：90%体积的靶区至少接受100%处方剂量照射量，即：D90＞100%；100%处方剂量至少照射的靶区体积的90%，即：V100＞90%；150%处方剂量照射的靶区体积不超过50%，即：V150＜50%；200%处方剂量照射的靶区体积不超过 25%，即：V200＜25%［10］。当满足以上评估参数时，剂量规划结果满足临床治疗标准。

1.5 DVH评估方法

多个剂量结果的优劣评估，主要依据剂量评估参数V100，V150和V200作为基准。如图3所示，S1和S2面积越大靶区接受的剂量越充分，规划结果越好；S3和S4面积越小，高剂量区域越小，规划果越好。由图所示推出DVH评价方法如下式：

 

式中：ε1，ε2，ε3和 ε4为权重因子，Vp表示处方剂量百分比是j时所对应的体积。

  

图3 DVH曲线评估方法

2 自动剂量规划的主要步骤

综上，提出了自动三维剂量优化算法，主要步骤如图4。

为加强施工现场管理，江苏省出台了《江苏省重点水利工程建设管理规定》，按照规定要求，重点水利工程建设必须实行项目法人责任制、招标投标制、建设监理制、合同管理制和竣工验收制。在建设运行管理过程中，实现集中管理和分级管理相结合的体制。在现场施工管理工作中，全面落实现场管理的各项制度，抓好现场运行管理，保证工程建设有序地开展。

3 实验平台与系统的搭建

综上设计的算法流程以及临床医师的经验，运用 VTK（visualization toolkit）、ITK（Insight Segmentation and Registration Toolkit）和 Visual Studio 2010 等工具，设计了面向体部肿瘤的近距离粒子植入计划软件。图5。

忽略粒子之间的剂量影响，多颗粒子在剂量空间剂量场的影响呈线性叠加关系。所以N颗粒子在空间 P（x，y，z）点的总剂量叠加值为：

4 计算机仿真与实验验证

4.1 术前剂量计划

图像选用层厚是5 mm的CT。手动分割靶区并完成重建之后，提取包裹靶区体积的最小立方体，并外扩50 mm作为剂量优化的计算范围；设置单颗粒子计算范围为50 mm×50 mm×50 mm；除此之外，穿刺针的布置间距为5～10 mm，穿刺针上可植入粒子位置点间距是5 mm。

1.3.1 目标函数 放射治疗剂量规划的目标是使处方剂量覆盖整个靶区，而且邻近器官所接受的放射剂量最小。然而，靶区剂量目标和危机器官剂量目标相互制约，即靶区剂量越高，周围器官接受的放射剂量越大，反之亦然。为了得到最优的剂量分布，目标函数综合考虑了靶区和危及器官的剂量覆盖情况，并根据器官的重要程度为每种组织器官设置相应的权重因子，定义如下：

  

图4 自动剂量规划算法流程图

 

（1）对肿瘤进行手动分割，重建出计划靶区。查询诺模图得出满足处方剂量需要的最少粒子数目。布置空间针，使之覆盖靶区所有位置，并计算出每根针上的粒子剂量叠加运算模型。结合空间针和PTV边界，计算出所有空间针上植入粒子的候选位置。（2）在初始粒子数目、粒子候选位置和粒子剂量叠加模型都准备完毕后，进入改进模拟退火算法，表示当前分布方案。①初始随机产生粒子布置方案。②添加扰动产生新的粒子布置方案。对方案进行剂量计算，将高剂量区域的一颗粒子，转移到低剂量区。③计算新旧两种方案的状态方程，若差值＜0或＞0且概率Pi＞0，1之间随机产生的数，则此时从状态i跳转到状态j，将设为当前粒子分布方案。 ④如果温度未达到最低温度则进行步骤②，否则返回当前状态最优解。然后进行剂量评估，如果评估通过则将此解作为最终结果。否则判断粒子数目是否达到上限，若果没有进行步骤①；反之，将所有粒子数目对应的最优结果，进行DVH评估，将评估式Q最小的粒子布置方案作为最终结果输出

1.3.2 模拟退火算法 剂量优化模拟退火算法将目标函数作为控制退火温度缓慢下降的参数，结合以一定的概率跳出的中间解，能够有效防止算法陷入局部最优。

还有安南（越南）自古以来就使用中国钱币。最早开始铸钱是在丁朝大瞿越大平年间(970)，仿照中国宋朝方孔圆钱的形制铸成的，钱文多有宋体之意。从黎仁宗延宁起，到昭宗光绍年间（1454—1521），则是安南铸钱的成熟时期，无论铸钱技术，还是工艺水平，都达到了较高水平。光绍以后，安南古钱的制作便逐步走向下坡。安南曾经大量仿制过中国古钱，与中国铜钱同名的有15种。

  

图5 体部肿瘤近距离粒子植入剂量规划系统（BTPS）

 

表1 靶区和脊髓的剂量参数

  

?

4.2 算法精度验证

随机选取10例体部肿瘤病例，物理医师分别使用Oncentra Seeds和本文的BTPS进行剂量规划。运用控制变量法，确保靶区大小、粒子型号、处方剂量（120 Gy）和粒子位置等变量一致，计算BTPS靶区的剂量参数与Oncentra Seeds的相对剂量差异，计算结果如表2。

从BTPS与Oncentra Seeds的对比结果可以看出，所有剂量参数的偏差小于3.6%，关键剂量评估参数D90、V100、V150和V200最大偏差分别为2.4%、1.2%、2.3%和3.0%。所以BTPS具备临床应用的准确性和安全性。

4.3 实验讨论

仿真实验从耗材上与实际手术方案一致，重要剂量参数满足临床要求，计算速度达到实时规划要求，验证了基于空间针的三维自动布源方法的临床适应能力。BTPS与商用软件的剂量对比验证了剂量算法的相对准确性和安全性。为进一步验证剂量算法的绝对准确性和安全性，后续需要对BTPS进行放射粒子剂量实验。

二、影响一国需求结构的最主要因素是平均收入水平。高收入国家对技术水平高、加工程度深、价值较大的高档商品的需求较大，而低收入国家则以低档商品的消费为主，以满足基本生活需求。所以，收入水平可以作为衡量两国需求结构或偏好相似程度的指标。例如高尔夫球在欧美是普及运动，但在发展中国家却不是代表性需求。

  

图6 剂量规划情况示意

 

①靶区和脊髓的三维重建；②三维剂量规划结果可视化展示；③DVH曲线和主要剂量参数；④3D打印非共面模板

 

表2 两种软件剂量算法的实验比较

  

?

在基于平面的二维术前剂量规划过程中，熟练的物理医师大约需要消耗20 min，应用本文提出的方法大约需要消耗6 min，明显提高了术前剂量规划速度。穿刺针的布置现依靠物理医师手动安放，这个过程比较耗时。为了减少布针的时间，可以通过保持穿刺针间距固定，且尽量使同一平面内的针保持平行，后续工作需要在穿刺针的布置方面进行深入研究和实验。

当课题通过教育科研部门鉴定后，并不意味着课题研究的结束，而应采取灵活多样的推广形式，及时进行成果的推广和应用，使更多的学校、教师了解成果，并汲取优秀成果的精华。把科研成果运用到教育实践中去，让其发挥出应有的作用。

基于空间针剂量规划每根针的方向都可能不同，理论上需要多个平面模板来实施。过多平面模板增加了手术操作复杂性，所以本软件通过提取裁剪患者表面皮肤，安放穿刺针导向柱，制作三维虚拟模板，并使用3D打印机实现三维模板制造。怎样使三维模板精准复位，确保剂量完全依照术前规划实施，仍需要进一步研究。

本文针对临床实际，研发了一种基于空间针的面向体部肿瘤近距离粒子植入剂量规划系统，分析了三维自动布源的理论基础，介绍了三维自动布源的算法和工作流程，通过计算机仿真实验重点验证了算法的速度和临床适应能力，经过软件剂量算法对比实验验证算法的准确性和安全性。实验结果表明该体部肿瘤近距离粒子植入剂量规划系统达到了临床要求，实现了精确、快速和有效的剂量规划。

［参 考 文 献］

［1］ 柴树德.我国放射性粒子植入治疗支气管肺癌现状与技术创新［J］.山东大学学报·医学版，2017，55：4-7.

［2］ 屈 佳.超声引导经皮徒手穿刺125I放射性粒子植入治疗中央型肝癌的临床研究［D］.吉林大学，2016.

［3］ Li J，Xie Q，Wang W，et al.CT-guided implantation of125I seeds（permanent brachytherapy） for metastatic tumors of the hepatic portal system： effectiveness and safety in 13 patients［J］.Brachytherapy， 2016， 15： 224-230.

［4］ Liu B， Zhou T， Geng J， etal.Percutaneouscomputed tomography-guided iodine-125 seeds implantation for unresectable pancreatic cancer［J］.Indian J Cancer， 2015， 52（Suppl 2）： 52-69.

［5］ Chaswal V， Thomadsen BR， Henderson DL.Template independent prostate brachytherapy simulation using adjoint sensitivity based treatment planning optimization technique ［J］.Brachytherapy， 2010， 9（Suppl 1）： S41.

［6］ 郭福新，姜玉良，吉 喆，等.3D打印非共面模板辅助CT引导125I粒子植入治疗锁骨上复发转移癌的剂量学研究［J］.北京大学学报·医学版，2017，49：506-511.

［7］ 唐富龙，任菊娜，吴 娟，等.125I放射性粒子单平面布源剂量学研究［J］.介入放射学杂志， 2014， 23： 619-622.

［8］ Rivard MJ，Coursey BM，Dewerd LA，et al.Update of AAPM Task Group No.43 Report：a revised AAPM protocol for brachytherapy dose calculations［J］.Med Phys， 2004， 31： 633-674.

［9］ 蔡东阳，卓子寒，王 婕，等.基于模拟退火算法的磁感应治疗热籽分布［J］.清华大学学报·自然科学版，2014，54：153-158.

［10］Kirisits C，Siebert FA，Baltas D，et al.Accuracy of volume and DVH parameters determined with different brachytherapy treatment planning systems［J］.Radiother Oncol， 2007， 84： 290.