专利名称：一种射频毁损仪的工作参数实时调控方法
技术领域：
本发明涉及一种射频毁损仪的工作参数实时调控的方法，尤其是利用近红外光谱 技术采集到的约化散射系数μ ‘ s对射频毁损仪的工作参数进行实时调控的方法。
背景技术：
射频(Radio frequency, RF)是一种频率为50kHz 2000kHz的电磁波，其电流高 度集中，并随高频产生离子振动磨擦，从而在组织内产生热量，使组织溶解、气化。毁损术对 于当今的帕金森病、各类肿瘤(原发性肝肿瘤、乳腺肿瘤、脑胶质瘤、前列腺瘤等)的毁损效 果较为满意。毁损程度与靶点定位精度、毁损灶大小、热凝温度、烧灼时间以及刺激参数的 选取等因素密切相关，目前国内临床医师针对毁损灶大小的毁损参数(毁损温度T和毁损 时间t)的选取都是通过先前的经验，术中的毁损效果也是看病人的反映来估计，缺乏实时 监控毁损程度的方法，这导致了手术治愈率不高且多复发的情况。临床上比较理想的射频毁损术中毁损程度的监控方法主要是术中MRI实时导航 和超声术中导航。术中MRI利用毁损靶点的温度场分布来进行毁损程度的监测，可以得出 毁损体积，评价毁损效果，得到国内外的普遍认可，但存在几个重要的问题1、术中MRI设 备价格昂贵，考虑到使用环境射频毁损设备也必须特别设计与制作，目前在我国绝大多数 医院的射频毁损手术没有这个条件。2、术中MRI实时监控的是温度，通过温度场分布来估 计毁损效果，不同温度点的毁损时间，其依据仍然缺少实时组织监控参数(如毁损靶点区 域的组织凝固程度等)，也就是缺少毁损程度的实时定量描述，或者说温度场分布不等于毁 损程度模型。近红外光谱(Near Infrared Spectroscopy，NIRS)技术的研究很早就已出现。 1977年，Jobsis在《Science》杂志上首次报道了血红蛋白和细胞色素在特定近红外区的吸 收特性，并发现氧合血红蛋白和还原血红蛋白分别在735nm和850nm处有两个吸收峰，其变 化可以反映血红蛋白的载氧情况。该报道引起了广泛重视。此后，很多研究小组从各自不 同的角度将这种方法推广开来，在不同的领域对其进行了深入的研究。目前，在生物医学领域NIRS技术被广泛应用于研究生物组织的光学性质、组织类型识别、诊断癌症、药物的传输过程监测以及组织结构和成分研究、组织内部散射粒子大 小、密度、近红外吸收药物的效应等。近10年来，随着光导纤维及传感技术的发展，近红外 光谱技术在医学领域得到了广泛的应用。由于生物体中不同的组织对近红外光具有不同的 吸收和散射特性，因此近红外光对不同组织和组织的变化具有较强的区分能力。根据这种 特性，可以利用近红外光谱法测量生物组织的某些光学参数(如约化散射系数)，从而与组 织的某些生理参数进行有意义的关联，建立光学参数和这些生理参数的关联数学模型。
发明内容
发明目的本发明的目的是针对目前射频损毁治疗仪实际工作中存在的盲目性高、精确性低的问题，提出一种利用近红外光谱技术，对射频毁损仪的工作参数进行实时调控的方法。技术方案本发明为实现上述发明目的，采用如下技术方案一种射频毁损仪的工作参数实时调控方法，包括以下步骤步骤一，采用射频毁损仪对目标体进行毁损时，分别设置η组对应的毁损温度T和 毁损时间t进行毁损，同时采用生物组织光学参数测量系统采集与每组毁损温度T和毁损 时间t相对应的实时约化散射系数μ ‘ s ；通过采集上述η组毁损温度Τ、毁损时间t以及 约化散射系数μ ‘ s的实时对应数据，建立约化散射系数μ ‘ s与毁损温度T和毁损时间 t的数学模型μ ‘ s = f1(T, t)(1)其中η为自然数；步骤二，在步骤一中，当射频毁损结束后，同时采集目标体毁损灶体积V与毁损温 度Τ、毁损时间t的实时对应数据，获得目标体毁损的体积V与毁损温度T、毁损时间t的数 学模型f2 V = f2(T, t)(2)步骤三，根据数学模型和数学模型f2，得出约化散射系数μ ‘ s与目标体毁损 的体积V的数学关系式(3)V = g(u ‘ s)(3)步骤四，在射频毁损仪正式采集数据前，通过核磁共振图像获得预期毁损目标体 的体积Vtl，根据数学模型f2确定毁损温度T和毁损时间t ；步骤五，通过数学关系式(3)获得与预期毁损目标体的体积Vtl相对应的预期约化 散射系数P ‘ sth ；步骤六，在射频毁损过程中利用生物组织光学参数测量系统，实时采集获得约化 散射系数μ ‘ s的值，当μ ‘ s的值达到预期值μ ‘ sth时，即控制射频毁损仪结束毁损。进一步的，本发明的射频毁损仪的工作参数实时调控方法的步骤一、步骤二中，T 的取值范围为50°C 100°C，t的取值范围为5s 120s。进一步的，本发明的射频毁损仪的工作参数实时调控方法的步骤四中，当预期毁 损目标体的体积\已知，而毁损温度T和毁损时间t均未知时，毁损温度T和毁损时间t的 确定具体步骤为步骤A、设定m组毁损时间t，根据数学模型f2，求得m组对应的毁损温度T的值， 得出m组关于预期毁损目标体的体积Vtl的毁损方案，T取值范围为50°C 100°C，t的取值 范围为5s 120s，m为自然数；步骤B、根据射频毁损仪可设置的温度、时间的最小精度，对步骤A的m组毁损方案 进行筛选，排除那些超出射频毁损仪设置精度范围的毁损方案，得出有效毁损方案；步骤C，对步骤B得到的有效毁损方案进行二次筛选，筛选出同一毁损时间下毁损 温度最低，或者同一毁损温度下毁损时间最短的方案，将其作为首选方案。进一步的，本发明的射频毁损仪的工作参数实时调控方法的步骤六中，当μ' s < μ ‘ sth时，根据μ ‘ s的曲线的态势，采取以下两种处理方法①、当μ ‘ s的曲线处于走平阶段时，即在设定的初始毁损温度T1下，毁损了、时长后，不管毁损时间再如何增大，毁损体积均不再增大，μ ' s不可能达到预期值μ' sth， 此时应重新设置毁损温度为T2，且T2 > T1,重新返回步骤四，采用射频毁损仪重新设置好毁 损温度T2和毁损时间t2后，继续进行毁损，同时监测y ‘ s的变化趋势；②、当μ' s的曲线处于上升阶段时，即在设定的初始毁损温度T1下，毁损了、 时长后，μ' s还没有达到极限值，同时毁损体积也还没有到达预期值；若延长毁损时间， μ ‘ s仍然将继续上升，直至达到预期值μ ‘ sth大小，对应的毁损体积也会继续增大，直至 无限接近目标体积Vtl，此时应保持毁损温度不变，即T2 = T1，延长毁损时间，即t2 >、，继续 进行毁损，直至μ' s= μ ‘ sth。 有益效果1、实时性强、准确性高。实验表明，毁损程度，即毁损体积的大小与μ ‘ s的变化 情况密切相关，因此μ ‘ s作为毁损程度的实时监控因子具有很好的效果。2、成本低。本发明提供的方法成本大大降低，生物组织光学参数测量系统的搭建 相对方便，实现起来也相对比较容易。3、可操作性强。本发明提出的方法，可以实现射频毁损仪的欠操作控制和过操作 控制，可操作性强，安全性高。
图1是本发明的流程图。图2是本发明的实施例中μ ‘ s变化情况。
具体实施例方式下面结合附图对本发明的技术方案进行详细说明附图1是本发明的详细流程图。采用射频毁损仪进行毁损操作前，通过其他医学 成像手段，如核磁共振图像获得目标毁损灶体积大小Vo，根据数学模型(2) :V = f2(T，t)设 定毁损温度和毁损时间分别为1\、、。毁损开始时，利用近红外光谱采集系统观测约化散射 系数μ' 3的变化。若其大小达到预期值μ' sth，则可认为达到毁损目的，可立即结束毁 损；若其大小未达到预期值μ ‘ sth，则进入方案修正环节，具体的情况分析和修正方案可 见具体实施方式
中步骤五。如图2所示，为约化散射系数μ ‘ s变化情况。附图2中三条曲线1、2、3是约化 散射系数μ' s的三种变化情况，圆点为μ ‘ sth的值。图中左侧为射频毁损电极在目标体 形成的毁损体剖面示意图，D为内置双光纤的射频毁损电极，A、B、C三个区域分别为术前测 量所得病灶区(体积为V、实际毁损形成的区域和毁损电极内置光纤的前视距离。以下对本发明进行详细说明。步骤一，定位。术前根据相关医学检测手段获得射频毁损灶的中心位置坐标、横 径、纵径大小等参数，经椭球体积计算公式计算出目标毁损灶的体积Vtl，并利用步进电机将 毁损电极尖端定位至中心坐标位置处。步骤二，确定毁损温度和毁损时间。根据毁损灶的体积V0和数学模型(2) :V = f2(T，t)，可求出毁损温度T和毁损时间t的各种毁损方案组合。从这些毁损方案组合中选 择最佳组合，现举一例进行说明。
如表1所示，当Vtl = 25. 5mm3时，表1给出了方案的筛选过程。初始方案是根据上 述方法求解得出对应于毁损时间t为5s、10s、15s、20s、……120s时各个毁损温度T的值， 现针对一具体的射频毁损实例加以说明，其他各类射频毁损术中毁损方案的选择均可参照 此技术路线。根据初始方案1，当给予的毁损时间是20s时，欲毁损Vtl = 25. 5mm3的体积， 则应设置毁损温度为99. 7°C;同理，根据初始方案3，当给予的毁损时间是30s时，欲毁损Vtl =25. 5mm3的体积，则应设置毁损温度为86. 2V。温度T变化范围为50°C 100°C，毁损时 间t变化范围为5s 120s，而毁损时间为5s、10s、15s时对应的毁损温度均超出最高温度 100°C，因此，初始方案中毁损时间为5s、10s、15s的方案不满足条件，并未给出。由于温度、时间值必须符合射频仪可设置的温度、时间的最小精度，因此，利用相 关算法对21种初始方法进行一次筛选，得出6种优选方法，这类方案属于有效方案，为方案 的选择提供了参考依据。为不造成过毁损，在该实例中首选方案通常采取取值相对保险的方案，因此可利 用算法对方案进行二次筛选。二次筛选是从优选方案中筛选出同一毁损时间下毁损温度最 低，或者同一毁损温度下毁损时间最短的方案，将其作为首选方案的一种，如当毁损温度均 为70°C时，相对保险的毁损方案是毁损时间最短的一组，即方案1、2、3、4为可操作的首选 方案。表1方案筛选过程
步骤三，计算。实时毁损流程见附图1。根据步骤二所得参数，计算目标毁损灶体 积Vtl (单位mm3)，由数学模型(2) :V = f2 (Τ, t)求得对应于V。的T和t的值，同时根据公式 (3)计算得出对应于该体积的约化散射系数的预期值μ ‘ sth。步骤四，监测。根据步骤三得到的T和t，制定首选毁损方案1\、、(毁损方案的合 理选择见步骤一)毁损病灶区。正式开始毁损时，记录约化散射系数μ ‘ s的变化情况。
步骤五，分析。实时监测的P ‘ s曲线见附图2，圆点为y ‘ sth的值，图中左侧为 射频毁损电极在目标体形成的毁损灶剖面示意图，D为内置双光纤的射频毁损电极，A、B、C 三个区域分别为术前测量所得病灶区(体积为、实际毁损形成的区域和毁损电极内置光 纤的前视距离。根据本发明的技术路线，比较实时监测获得的P ‘ 3与步骤三中计算得出的 U ‘ sth大小，有两种显著的结果其一，在首选方案，即毁损温度为，毁损时间为、的情况 下，ii' U ‘ sth，见附图2中的曲线3，此时的i! ‘ s达到了预期值，说明当前毁损形成 的毁损区域刚好无限接近或者超出了目标毁损灶大小，再继续毁损则很有可能烧伤组织中 的神经或其他重要部位，应立即结束毁损，这实现了对毁损的过操作控制，此种情况可以视 为一旦达成，则基本完成毁损；其二，毁损温度为1\，毁损时间达到、时，y ‘ s< u ‘ sth， 此时说明当前毁损形成的毁损区域还未达到目标毁损灶大小，仍需继续毁损，实现了对毁 损过程的欠操作控制。现针对欠操作控制的修正方案进行详细说明P ‘ s< U ‘ sth时，说明约化散射 系数s没有达到预期值，见附图2中的曲线1、2，其中曲线1为y ‘ 3处于走平阶段时 的情况，此时表明在设定的初始毁损温度下，毁损了、时长后，不管毁损时间再如何增 大，毁损灶体积均不再增大，即P ‘ s不可能达到预期值sth。针对性的修正方案为根 据优选方案重新设置毁损温度为T2，且T2 > 1\，T2所对应的毁损时间t2由模型t2 = f(T2, V0)得出。设置好参数值后，继续毁损，观测修正后P ‘ s的变化趋势。约化散射系数P ‘ 3没有达到预期值的另一种情况如附图2中的曲线2，此时 U ‘ s处于上升阶段。这说明在设定的初始毁损温度下，毁损了、时长后，s还没有 达到极限值。这种情况下的毁损体积还没有到达极限，若延长毁损时间，P ‘ s仍然有可能 上升，甚至达到预期值P ‘ sth大小，其对应的毁损体积也有可能继续增大，甚至无限接近 目标体积％。针对性的修正方案为保持毁损温度不变，即t2 = ，延长毁损时间，即t2 > 、，直至 ii ‘ s = u ‘ sth。
权利要求
一种射频毁损仪的工作参数实时调控方法，其特征在于包括以下步骤步骤一，采用射频毁损仪对目标体进行毁损时，分别设置n组对应的毁损温度T和毁损时间t进行毁损，同时采用生物组织光学参数测量系统采集与每组毁损温度T和毁损时间t相对应的实时约化散射系数μ′s；通过采集上述n组毁损温度T、毁损时间t以及约化散射系数μ′s的实时对应数据，建立约化散射系数μ′s与毁损温度T和毁损时间t的数学模型f1μ′s＝f1(T，t)(1)其中n为自然数；步骤二，在步骤一中，当射频毁损结束后，同时采集目标体毁损灶体积V与毁损温度T、毁损时间t的实时对应数据，获得目标体毁损的体积V与毁损温度T、毁损时间t的数学模型f2V＝f2(T，t)(2)步骤三，根据数学模型f1和数学模型f2，得出约化散射系数μ′s与目标体毁损的体积V的数学关系式(3)V＝g(μ′s)(3)步骤四，在射频毁损仪正式采集数据前，通过核磁共振图像获得预期毁损目标体的体积V0，根据数学模型f2确定毁损温度T和毁损时间t；步骤五，通过数学关系式(3)获得与预期毁损目标体的体积V0相对应的预期约化散射系数μ′sth；步骤六，在射频毁损过程中利用生物组织光学参数测量系统，实时采集获得约化散射系数μ′s的值，当μ′s的值达到预期值μ′sth时，即控制射频毁损仪结束毁损。
2.根据权利要求1所述的射频毁损仪的工作参数实时调控方法，其特征在于所述步 骤一、步骤二中，T取值范围为50°C 100°C，t的取值范围为5s 120s。
3.根据权利要求1所述的射频毁损仪的工作参数实时调控方法，其特征在于所述步 骤四中，当预期毁损目标体的体积\已知，而毁损温度T和毁损时间t均未知时，毁损温度 T和毁损时间t的确定具体步骤为步骤A、设定m组毁损时间t，根据数学模型f2，求得m组对应的毁损温度T的值，得出 m组关于预期毁损目标体的体积％的毁损方案，T取值范围为50°C 100°C，t的取值范围 为5s 120s，m为自然数；步骤B、根据射频毁损仪可设置的温度、时间的最小精度，对步骤A的m组毁损方案进行 筛选，排除那些超出射频毁损仪设置精度范围的毁损方案，得出有效毁损方案；步骤C，对步骤B得到的有效毁损方案进行二次筛选，筛选出同一毁损时间下毁损温度 最低，或者同一毁损温度下毁损时间最短的方案，将其作为首选方案。
4.根据权利要求1所述的射频毁损仪的工作参数实时调控方法，其特征在于在步骤 六中，当y ‘ s< u ‘ sth时，根据P ‘ s的曲线的态势，采取以下两种处理方法①、当P ‘ s的曲线处于走平阶段时，即在设定的初始毁损温度1\下，毁损了t1时长后，不管毁损时间再如何增大，毁损体积均不再增大，P ‘ s不可能达到预期值sth，此 时应重新设置毁损温度为T2，且T2 > ，重新返回步骤四，采用射频毁损仪重新设置好毁损 温度T2和毁损时间t2后，继续进行毁损，同时监测P ‘ s的变化趋势；②、当P ‘ s的曲线处于上升阶段时，即在设定的初始毁损温度1\下，毁损了、时长 后，s还没有达到极限值，同时毁损体积也还没有到达预期值；若延长毁损时间，P ‘ s 仍然将继续上升，直至达到预期值P ‘ sth大小，对应的毁损体积也会继续增大，直至无限 接近目标体积％，此时应保持毁损温度不变，即T2 = ，延长毁损时间，即t2 >、，继续进行 毁损，直至P ‘ s= u ‘ sth。
全文摘要
本发明公开了一种射频毁损仪的工作参数实时调控方法，属于近红外光谱数据处理技术领域。本发明根据射频毁损仪采集n组毁损温度T和毁损时间t对应值，同时通过光学参数测量系统实时采集得到与其同步的n组约化散射系数μ′s，得出关系式μ′s＝f1(T，t)，然后再通过同步采集到的毁损体积V的数据，得到毁损体积V与毁损温度T和毁损时间t的关系式V＝f2(T，t)；通过以上两个关系式得出约化散射系数μ′s与毁损体积V的关系式V＝g(μ′s)；本发明通过关系式V＝g(μ′s)，根据预期毁损体积V0设定相对应的预期约化散射系数μ′sth，当μ′s＝μ′sth时，射频毁损过程结束；当μ′s＜μ′sth时，根据μ′s的曲线态势调整射频毁损仪的参数，以解决以前射频毁损仪毁损盲目性高、精确性低的问题。
文档编号A61B18/18GK101836892SQ20101016422
公开日2010年9月22日 申请日期2010年5月6日 优先权日2010年5月6日
发明者李韪韬, 肖笛, 钱志余 申请人:南京航空航天大学