一种放射性支架空间剂量分布的计算方法.pdf

1、10申请公布号CN102339360A43申请公布日20120201CN102339360ACN102339360A21申请号201110291941122申请日20110930G06F19/0020060171申请人成都君晟科技有限公司地址610000四川省成都市高新区紫荆北路19号72发明人张会雄时曼娜陈俊如帅先高珊54发明名称一种放射性支架空间剂量分布的计算方法57摘要放射性支架通过放射性核素发出的射线抑制血管平滑肌细胞VSMC的增生,从而抑制内膜的增厚,是预防PTCA经皮穿刺冠状动脉腔内成形术后血管再狭窄的有效工具。放射剂量模拟非常重要，因为这是对放射治疗疗效的预测，同时也是保护病人的

2、前提。本发明提供一种放射性支架空间剂量分布的计算方法，该方法中将支架的丝用一系列线性方程表示，简化了计算过程。本发明的有益效果是，计算过程相对简单，网格形状只需改变方程即可适应不同支架，剂量的模拟计算简单方便，不容易出现计算错误。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书4页附图3页CN102339373A1/2页21一种放射性支架空间剂量分布计算方法，其特征是将支架网线简化为一系列线性方程来处理，主要包括以下步骤，（1）钢丝网的坐标表示把支架展开，侧面是一个942MM20MM的矩形（长20MM，宽942MM），如图1所示,钢丝网的坐标位置可以在直角坐标

3、系（XZ平面）中表示，原点是（0，0）。2用简单的数学几何考虑，如果把长分为20等分，间隔1毫米，钢丝网之间相交的角度是45，所有的钢丝可以用58条方程表示,用N来标记这些等式，可以分成3组（EQ1）一旦和N被随机产生了，一个肯定在钢丝网上的放射源点R也就随之产生了（2）判断BETA射线穿过支架时是否穿过0根，1根或两根钢丝从支架背面射出的电子在在组织中释放能量之前会穿过支架的正面,算出BETA射线的轨迹和支架表面的交点位置，并计算交点和每根钢丝的距离，如果有一个或两个距离小于00015厘米，（不锈钢钢丝的直径是0003厘米），就认为电子将穿过一根或两根钢丝,根据图2图3计算交点BEQ2L4L

4、2L3,IE,B和所有不锈钢钢丝的距离EQ3（3）计算等效距离要正确计算支架材料的自吸收，很需要确定从各个方向穿过钢丝的电权利要求书CN102339360ACN102339373A2/2页3子所走的等效路程，一个方法是用蒙特卡罗，一系列的点，它们的圆心角是,，在一根钢丝的圆周侧面上被随机选取，在它们和测量剂量的场点（测试点）之间画一条线（如图4与图5所示），将所有直线穿过钢丝时所走的路程求和，再求平均，于是各个方向的BETA射线穿过钢丝时走过的等效路程就被计算出来了，（EQ4）D是测量剂量的场点和钢丝中心的距离，钢丝的半径是R，等效路程可以计算出来了；（4）计算支架近场的剂量分布放射源是分散成

5、一系列分布在支架表面的具有单位积累活度的源点，再把各自的贡献求和，吸收剂量就可以计算出来，以上步骤确定了放射点源的具体位置；源点和场点（计算吸收剂量）之间的距离；射线是否穿过钢丝；射线从各个方向穿过钢丝所走过的等效距离这些参量以后，就可以分开两种情况1）从支架正面放射源射出的射线，不经过钢丝，用点核函数公式；2）从支架背面放射源射出的射线，经过钢丝，用点核函数公式，计算两种情况下分别在场点产生的吸收剂量，最后把所有的源点产生的吸收剂量求和，那么某一确定场点的吸收剂量就能计算出来。权利要求书CN102339360ACN102339373A1/4页4一种放射性支架空间剂量分布的计算方法技术领域00

6、01一种放射性支架空间剂量分布的计算方法，属于放射医学领域，主要应用于计算放射性支架的空间剂量分布。背景技术0002放射性支架通过放射性核素发出的射线抑制血管平滑肌细胞VSMC的增生,从而抑制内膜的增厚,是预防PTCA经皮穿刺冠状动脉腔内成形术后血管再狭窄的有效工具。放射剂量模拟非常重要，因为这是对放射治疗疗效的预测，同时也是保护病人的前提。0003在WVPRESTWICH等人的早期工作中（PRESTWICH,WV,JNUNES,ANDCSKWOK,BETADOSEPOINTKERNELSFORRADIONUCIDESOFPOTENTIALUSEINRADIOIMMUNOTHERAPYTHEJ

7、OURNALOFNUCLEARMEDICINE,198930P10361046），给出了发射BETA射线核素的标度剂量点核函数（DPK），式中的参数是以LEVENBERGMARQUARDT算法为基础的非线性最小二乘法拟合出来的。在WVPRESTWICH等人之后的工作中（WVPRESTWISH,TJKENNETT,ANDFWKUS,THEDOSEDISTRIBUTIONPRODUCEDBYAP32COATEDSTENTMEDICALPHYSICS,1995223P313320），用这个函数计算了P32电镀表面的均匀圆柱支架剂量的空间分布。后来，CJANICKI等人（JANICKI,C,ETAL,

8、RADIATIONDOSEFROMAPHOSPHOROUS32IMPREGNATEDWIREMESHVASCULARSTENTMEDICALPHYSICS,1997243P437445），对WVPRESTWICH的剂量点核函数作了修正，加入了照射时间的影响因素，模拟了PALMAZSCHATZ支架（侧面为钢丝网）近场的剂量分布，并用放射线剂量胶片测量验证。然而，支架材料对放射能量的自吸收是产生放射剂量在血管里的分布的更加准确预测的一个不可忽视的因素。当BETA射线从支架背面的放射源射出，在射线到达增生组织并释放能量之前，要穿过支架正面的钢丝网，钢丝网要吸收一部分的放射能量，这时，自吸收就产生了。

9、REYNAERT等人（REYNAERT,N,ETAL,MONTECARLOCALCULATIONSOFDOSEDISTRIBUTIONSAROUND32PAND198AUSTENTSFORINTRAVASCULARBRACHYTHERAPYMEDICALPHYSICS,1999268P14841491）用EGS4MONTECARLO模拟方法计算了磷32PALMAZSCHATZ支架的三维剂量分布。发现在距支架中心轴2MM（即距支架壁125MM）处，由于支架本身对射线的吸收，使支架的剂量下降了23。两年后，REYNAERT等人（REYNAERT,NANDUOHFELI,SELFABSORPTION

10、CORRECTIONFORP32,AU198ANDRE188STENTSDOSEPOINTKERNELCALCULATIONSVERSUSMONTECARLOMEDICALPHYSICS,2001289P18831897）又将考虑了自吸收的磷32ACSMULTILINK和PALMAZSCHATZ支架放射剂量EGS4，MCNP4BMONTECARLO模拟和DPK计算和放射线剂量胶片测量，三种方法结果作对比。0004在以上方法中，均匀圆柱支架是理想情况，实际应用中支架都是网格结构的。PALMAZSCHATZ支架和ACSMULTILINK支架的网格结构都很复杂，给吸收剂量的计算带来很大的难度，因为复

11、杂的结构不容易被理解和准确描述，所以加大了计算误差的可能。说明书CN102339360ACN102339373A2/4页5发明内容0005本发明提供一种放射性支架空间剂量分布计算方法，该方法中将支架的丝用一系列线性方程表示，简化了计算过程。0006本发明的有益效果是，计算过程相对简单，网格形状只需改变方程即可适应不同支架，剂量的模拟计算简单方便，不容易出现计算错误。0007附图说明图1是所选支架的侧面展开图，本图中我选择模拟的支架规格是20MM长，直径是3MM，因为这是其中一种流行的在动脉中使用广泛的支架。把支架展开，侧面是一个942MM20MM的矩形（长20MM，宽942MM）。图2是支架前

12、视图，图3是支架俯视图。图4是一个理想的放射源点与一根不锈钢钢丝，图5是由此放射源放出的射线穿过此不锈钢钢丝的不同路径。图6是取样测量点的在考虑与不考虑支架自吸收时的剂量对比，图7是支架的轮廓图，图8考虑和忽略自吸收时感兴趣点剂量的差异。0008此发明的实施方案包括以下步骤，（1）钢丝网的坐标表示本实施方案中，我选择模拟的支架规格是20MM长，直径是3MM，因为这是其中一种流行的在动脉中使用广泛的支架。把支架展开，侧面是一个942MM20MM的矩形（长20MM，宽942MM），如图1所示。钢丝网的坐标位置可以在直角坐标系（XZ平面）中表示，原点是（0，0）。用简单的数学几何考虑，如果把长分为2

13、0等分，间隔1毫米，钢丝网之间相交的角度是45，所有的钢丝可以用58条方程表示。用N来标记这些等式，可以分成3组（EQ1）一旦和N被随机产生了，一个肯定在钢丝网上的放射源点R也就随之产生了。0009此例中N是在158取随机，是在02PI之间随机选取。支架网线之间相交角度作为一个可变量，以适应不同的支架。如果相交角度是45度之外的其他角度，以上式子将比此实施方案中的EQ1复杂。0010（2）判断BETA射线穿过支架时是否穿过0根，1根或两根钢丝从支架背面射出的电子在在组织中释放能量之前会穿过支架的正面。算出BETA射线的轨迹和支架表面的交点位置，并计算交点和每根钢丝的距离，如果有一个或两个距离小

14、于00015厘米，（不锈钢钢丝的直径是0003厘米），就认为电子将穿过一根或两根钢丝。0011根据图2图3计算交点BEQ2说明书CN102339360ACN102339373A3/4页6L4L2L3,IE,B和所有不锈钢钢丝的距离EQ3（3）计算等效距离要正确计算支架材料的自吸收，很需要确定从各个方向穿过钢丝的电子所走的等效路程。一个方法是用蒙特卡罗。一系列的点，它们的圆心角是，在一根钢丝的圆周侧面上被随机选取，在它们和测量剂量的场点（测试点）之间画一条线（如图4与图5所示），将所有直线穿过钢丝时所走的路程求和，再求平均，于是各个方向的BETA射线穿过钢丝时走过的等效路程就被计算出来了。（EQ

15、4）D是测量剂量的场点和钢丝中心的距离，钢丝的半径是R，等效路程可以计算出来了。0012（4）计算支架近场的剂量分布放射源是分散成一系列分布在支架表面的具有单位积累活度的源点，再把各自的贡献求和，吸收剂量就可以计算出来。以上步骤确定了放射点源的具体位置；源点和场点（计算吸收剂量）之间的距离；射线是否穿过钢丝；射线从各个方向穿过钢丝所走过的等效距离这些参量以后，就可以分开两种情况1）从支架正面放射源射出的射线，不经过钢丝，用点核函数公式；2）从支架背面放射源射出的射线，经过钢丝，用点核函数公式，计算两种情况下分别在场点产生的吸收剂量，最后把所有的源点产生的吸收剂量求和，那么某一确定场点的吸收剂量

16、就能计算出来。0013模拟结果如图6所示。0014支架初始活度（放射性同位素）是37KBQ（）；20毫米长，直径3毫米的圆柱说明书CN102339360ACN102339373A4/4页7体支架；放射时间是143D（磷32的半衰期），感兴趣点离支架中轴16MM处，（相当于离支架表面01MM处）。轴向。0015绿色点表示模拟时，不考虑自吸收；而红色点表示模拟时考虑自吸收。0016支架边缘除外，不考虑自吸收时剂量峰值是20GY，谷值是10GY。0017考虑自吸收时剂量峰值是14GY，谷值是88GY支架初始活度（放射性同位素）是37KBQ（）；20毫米长，直径3毫米的圆柱体支架；放射时间是143D（磷32的半衰期），感兴趣点离支架中轴16MM处，（相当于离支架表面01MM处）。轴向。0018模拟结果显示，自吸收约占初始放射能量的215，或者说自吸收占据了215的剂量。因为沿支架轴向取一系列的感兴趣点测量，自吸收都是215，这个常数。A不考虑自吸收时的剂量；B考虑自吸收时的剂量；C自吸收；公式C（AB）/A。说明书CN102339360ACN102339373A1/3页8图1图2图3图4说明书附图CN102339360ACN102339373A2/3页9图5图6说明书附图CN102339360ACN102339373A3/3页10图7图8说明书附图CN102339360A