基于电阻抗层析成像的高强度聚焦超声温度监测技术-郭各朴.pdf

《基于电阻抗层析成像的高强度聚焦超声温度监测技术-郭各朴.pdf》由会员分享，可在线阅读，更多相关《基于电阻抗层析成像的高强度聚焦超声温度监测技术-郭各朴.pdf（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、基于电阻抗层析成像的高强度聚焦超声温度监测技术郭各朴宿慧丹丁鹤平马青玉Noninvasive temperature monitoring for high intensity focused ultrasound therapy based on electricalimpedance tomographyGuo Ge-Pu Su Hui-Dan Ding He-Ping Ma Qing-Yu引用信息Citation: Acta Physica Sinica , 66, 164301 (2017) DOI: 10.7498/aps.66.164301在线阅读View online: http

2、:/dx.doi.org/10.7498/aps.66.164301当期内容View table of contents: http:/ you may be interested in表面粗糙度对固体内部超声背散射的影响Effects of surface roughness on diffuse ultrasonic backscatter in the solids物理学报.2016, 65(21): 214301 http:/dx.doi.org/10.7498/aps.65.214301球形集声器在生物组织中形成的组织损伤Tissue lesion induced by a sphe

3、rical cavity transducer物理学报.2014, 63(4): 044301 http:/dx.doi.org/10.7498/aps.63.044301基于分数导数研究高强度聚焦超声的非线性声场Study of nonlinear acoustic field of high intensity focused ultrasound by the fractional wave物理学报.2013, 62(5): 054301 http:/dx.doi.org/10.7498/aps.62.054301万方数据物理学报Acta Phys. Sin. Vol. 66, No.

4、16 (2017) 164301基于电阻抗层析成像的高强度聚焦超声温度监测技术 郭各朴宿慧丹丁鹤平马青玉y(南京师范大学物理科学与技术学院,南京210023)(2017年4月8日收到; 2017年6月7日收到修改稿)作为一种对正常组织无损伤且不易引起癌细胞转移的非入侵肿瘤治疗手段,高强度聚焦超声(HIFU)治疗过程中焦域的温度监测是实现剂量精准控制的关键.本文基于生物组织的温度-电阻抗的关系,将电阻抗层析成像(EIT)和HIFU治疗相结合,提出了一种利用组织焦平面的表面电压实现电阻抗重构的检测技术.建立了HIFU治疗和EIT综合系统模型,在考虑组织的声吸收条件下,对三维Helmholtz方程在

5、柱坐标下的声场计算进行了二维简化,并引入Pennes生物热传导方程来计算HIFU焦域的声压和温升分布特性;引入生物组织的温度-电阻抗关系,基于麦克斯韦电磁场理论,建立了具有温度分布HIFU焦域的电流和电压计算模型,利用恒流注入的边界条件实现电场计算,获得焦平面的表面电压分布.在数值计算中,利用实验聚焦换能器参数,模拟了在固定声功率下组织焦域的声场和温度场分布,以及中心和偏心聚焦条件下不同治疗时刻的电导率分布;然后通过对称电极的循环电流注入,计算了组织模型焦平面内的电流密度和电势分布,获得了焦平面圆周分布的表面电极电压;进一步采用修正的牛顿-拉夫逊算法,利用32 32的表面电极电压实现了焦平面内

6、电导率分布的重建.结果表明,基于温度-电阻抗关系的EIT电导率重建技术不但能准确定位HIFU焦域中心,还能恢复HIFU治疗中焦域的温度分布,证明了EIT用于HIFU治疗中温度监测的可行性,为其疗效评估和剂量控制提供了一种无创电阻抗测量和成像新方法.关键词:高强度聚焦超声,电阻抗层析成像,温度监测,表面电极电压PACS: 43.80.Ev, 43.35.Yb DOI: 10.7498/aps.66.1643011引言高强度聚焦超声1 4(high intensity focusedultrasound, HIFU)是一种具有广阔应用前景的无创肿瘤治疗技术,它利用超声波在组织中的穿透性和易聚焦性,

7、将换能器发射的超声波汇聚到肿瘤靶区,利用组织的声热效应产生65 C以上的高温,实现肿瘤组织在短时间内凝固性坏死从而达到治疗肿瘤的目的.在HIFU治疗过程中,既要杀灭肿瘤细胞,又不损伤周围的正常组织,准确的温度控制是关键,其中实时温度和疗效监测对HIFU的临床应用具有重要的意义.为了避免在HIFU治疗中插入测温探针,减少探针对HIFU声场的影响并降低癌细胞转移概率5,国内外学者提出了多种无创测温技术.微波测温技术利用组织温度和热辐射的关系,通过体外测量体内的热辐射来推测体内温度,但渗透深度有限,测量精度较差.磁共振成像测温技术6通过温度相关的扩散系数、质子共振频率或弛豫时间的测量实现组织温度图像

8、的重建,具有无创伤、无电离辐射、高温度分辨率的优点,但其时间分辨率不高.通过将超声测温技术7 11和HIFU系统进行融合,利用声速和回波时移以及非线性等参数实现温度国家自然科学基金(批准号: 11474166, 11604156)、江苏省自然科学基金(批准号: BK20161013)、国家博士后基金(批准号:2016M591874)和江苏高校优势学科资助的课题.通信作者. E-mail: 2017中国物理学会Chinese Physical Society http:/164301-1万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 164301监

9、测,但它们的温度变化较小,测量精度较低.近年来B超12被用来进行HIFU定位引导,监测治疗前后组织的供血变化,但是还不能实现高精度的温度监控和实时疗效评价.研究表明,生物组织也是一种导电体,在低频信号( I2.在一维声传播条件下,声场中单位时间内单位体积的媒质所吸收的声能量可以由声强表示:Q = I2dS I1dSdzdS = dIdz: (6)在三维空间的声传播中,热源则可以由声强的空间梯度28计算,Q = I: (7)在平面波近似下,声波沿z方向传播,若入射超声强度I0,传播距离z处的声强可以表示为I = I0 exp( 2 az),其中 a为吸收系数,代入(7)式得到28:Q(z) =

10、dIdz = 2 aI0 exp( 2 az) = 2 aI: (8)在生物组织中,入射声能量的损失一般由声衰减系数 表征,其为吸收系数 a和散射系数 s的总和.实际计算中很难区分组织对声能量的吸收和散射分量,一般直接将衰减系数等效为吸收系数28,即a.在忽略介质中的热耗散和热传导,传播距离z处的超声热源表示为28Q = 2 aI 2 I = 2 tctp2: (9)164301-3万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 164301在声场计算中, HIFU焦点处的声强可以由声压各次谐波幅值表示28:I = p2tct = 21n=1jCn

11、j2/ tct; (10)其中, t和ct分别为组织密度和声速, Cn是n次谐波的声压幅值.在本研究中所用超声功率不大,忽略非线性影响, (10)式可以简化为用基波声压计算.在不考虑血液流动影响的前提下,将组织所吸收的热量应用到Pennes生物热传导方程28;29中,tCtTt = Kt2T + Q; (11)其中Kt是组织热导率, T为组织温度, T0为初始温度.由于HIFU的热效应,在焦点处形成中心温度最高、周围温度较低、具有明显的温度梯度分布的椭球状(mm)的温升焦域.为了定量分析组织电阻抗随温度的变化,将TIVF14 16应用到电导率中,得到组织的温度-电导率分段函数为(T) =8:0

12、:73; (T = 293K);0:73 1:0204T T0; (293K 343K) TIVF = 2%: (12)可见,随着HIFU治疗时间的增长,生物组织焦域的温度逐渐升高,电导率随温度的升高而增大,在70 C组织凝固时产生快速的电导率提升,在焦域内形成电导率的梯度分布,这为EIT在HIFU中的温度监测提供了物理基础.基于EIT的HIFU焦域温度监测技术的正问题研究包括声场、温度和电导率分布,焦平面内的电势分布以及边界电极的电压分布.图1(b)显示了HIFU治疗中组织模型(半径为R)焦平面的电导率分布,在边界圆周上均匀设置N个电极,用对称电极1/17电流注入,计算焦平面内的电流密度和电

13、势分布,得到32个边界电极的电压值;然后利用对称电极的循环电流注入,重复以上过程,得到边界电极电压Vij(i;j = 1;2; ;N),其中i和j表示激励电极和测量电极的序号.EIT正问题求解是在已知HIFU焦域的内部电阻抗分布的前提下,且具有特殊边界条件的电场计算30;31,组织模型可以等效为导体,其电场分布满足麦克斯韦方程组32,在10100 kHz的低频电流激励下,忽略介电常数的影响,得到场域的表达式33:8: : (T) = 0;1 : = 0;2 : (T)n = Jn;(13)其中, (T)为场域内电导率分布, T为节点的温度;1为第一类边值条件,表示已知位函数在场域边界上各点的值

14、, 为场域内电势分布函数, 0为给定的边界电位,初始时认为给定点边界电位为0; 2为第二类边值问题,表示已知位函数在场域边界上各点的法向导数值; Jn是给定边界注入的电流密度.对(13)式采用变分法求解33,建立拉普拉斯方程的泛函I() = 12( (T)d;根据格林定理,考虑模型的边界条件,并使整个场域的总能量最小,得到8: : minI()=min12 (T)()2d+Jnd ;1 : = 0:(14)在HIFU治疗中,焦域的电导率分布随着温升而改变,结合激励电极的边界条件,利用有限元算法可以计算出组织模型内各剖分单元和节点的电压,进一步获得边界电极测量电压Vij.基于EIT的HIFU焦域

15、温度监测的逆问题是通过已经获得的模型边界电极电压和电流激励模式,利用Vij重建出模型内的电阻抗分布.修正的牛顿-拉夫逊(MNR)算法18;34;35是通过不断迭代来改变电阻率分布,进而使目标函数(重建的边界电极电压和测量电压之间误差范数的平方)最小.如组织的电导率为 (T),则相应的电阻率分布为= 1/ (T) =(1 2 m 1 m)T.在HIFU治疗的正问题中,通过对称电极的循环电流注入得164301-4万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 164301到32 32边界电极测量电压Vij;在MNR电阻抗重建的逆问题中,假设相应对称电极

16、电流注入时计算得到边界电极电压为Uij( ),则目标函数为f( ) = 12Ni=1U( ) V 2= 12NiNj(Uij( ) Vij)2; (15)其中N = 32.通过不断迭代使目标函数最小,得到焦平面的稳定电阻率分布.为使目标函数最小,令f( ) = U( )(U( ) V ) = 0,其中U( )称为Jacobian矩阵36.对f( )泰勒级数展开,并保留线性项得到f f( k)+f( k) k,其中 k =k+1 k, f( k) = J( k)TJ( k)为Hessian矩阵, k = J( k)TJ( k) 1J( k)TU( k) V ,得到MNR的迭代公式为k+1 = k

17、 + k; (16)其中k是迭代次数.为了避免对病态J( k)TJ( k)求逆,在重建中引入Tikhonov正则化33修正,并将其补偿项表示为 的平方函数形式,得到重构电阻率分布的迭代公式为k+1 = k + k= k fJ( k)TJ( k) + Rg 1J( k)TU( k) V : (17)将MNR重建方法应用到HIFU焦域的电阻率重建中,先基于EIT正问题得到的边界电极测量电压Vij,假设模型焦平面内电阻率分布均匀,利用有限元模型计算得到边界电极电压Uij( ),建立目标函数f( ),通过求解不同电阻率分布下的雅克比矩阵和海森矩阵,获得新的迭代电阻率分布,直到目标函数小于预设值,此时的

18、电阻率分布 k即为重建出的组织模型焦平面的电阻率分布,进一步利用(T) = 1/ k可以计算出电导率分布.3数值计算如图1所示,由于换能器和声传播的对称性,在有限元37计算中, HIFU声场和温度场以及组织电导率分布采用二维轴对称的柱坐标模型,其中轴向z是声传播方向, r是半径方向,仿真区域为超声换能器、水域环境以及组织模型.通过调整换能器表面振速控制输出声功率,计算组织焦域的声场、温度场以及电导率分布.为了在保证计算精度的前提下提高计算速度,水域环境区域剖分网格尺寸为 Water/4,组织区域和焦域的剖分网格尺寸分别 Tissue/4和 Tissue/8.设置聚焦超声换能器的直径和焦距均为1

19、0 cm,中心频率1.13 MHz.直径和高度分别为32 mm和35 mm的圆柱形组织模型中心放在超声换能器的焦域处.在模型焦平面的表面均匀设置32个电极测量表面电极的电压.在HIFU治疗中,随着超声的作用,模型内焦域的温度升高,其电导率随之提高,形成电导率的梯度分布.模型采用仿组织透明凝胶38,其物理参数和人体组织较为接近.计算中水和凝胶组织以及人体组织的相关参数列于表1.结合HIFU系统的实验参数,将换能器表面振幅设定为10.2 nm,通过有限元计算得到焦域处的声压以及声强分布,定义焦平面内声压衰减6 dB的面积为有效截面面积,通过AcsI(x;y)dA计算得到声功率为15.68 W17,

20、并以此作为声源参数进行相关计算.表1温度为293 K (20 C)时的仿真参数Table 1. Parameters used in simulation at the tem-perature of 293 K (20 C).参数水凝胶人体组织密度/kg m 3 1000 1044 10001100速度/m s 1 1500 1568 14501640衰减系数/Np m 1 MHz 1 0.025 8.55 4.0317.27比热容/J kg 1 K 1 N/A 3710 36003890导热系数/W m 1 K 1 N/A 0.59 0.450.56电导率/S m 5.5 10 6 0.7

21、3 0.8相对介电常数81 40 40通过有限元仿真得到固定声功率15.68 W时HIFU焦域的轴向剖面和径向焦平面的声压分布,结果如图2(a)和图2(b)所示,可见HIFU焦域呈现椭球状,焦点处的声压最大,随着轴向和径向范围的扩大,声压逐渐降低.结合(11)式,计算不同治疗时间(t)组织焦域的温升,得到如图3所示的焦平面温度分布.在t = 1 s时,焦域半径小于0.5 mm,焦域中心的温度较低,未达到70 C.随着治疗时间的延长,焦域的能量逐渐积累,焦点及周围组织的温度不断升高,同时由于组织的热扩散,164301-5万方数据物理学报Acta Phys. Sin. Vol. 66, No. 1

22、6 (2017) 164301径向距离/mm112 14121086420-2-4-6-8-10-12-14110108106104102轴向距离/mm轴向距离/mm10098969492908886-15 -10 -5 0 5 10 15径向距离/mm-15 -10 -5 0 5 10 156X106 X10654321654321图2 (网刊彩色)固定声功率为15.68 W时, HIFU的(a)轴向剖面和(b)焦平面的声压分布Fig. 2. (color online) Pressure distributions of (a) axial prole and (b) focal plan

23、e for HIFU at a xedacoustic power of 15.68 W.径向焦平面温度分布/K2 s1 s 3 s径向焦平面电导率分布/S.m-11s 2s 3s径向焦平面电导率分布/S.m-11 s 2 s 3 s-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-404812

24、16-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-40481216-16-15 -10 -5 0 5 10 15-12-8-404812163503403303203103003.02.52.01.51.03.02.52.01.51.0(a)(b)(c)图3 (网刊彩色)不同治疗时刻HIFU焦平面的(a)温度和(b)电导率分布及(c)偏心聚焦的电导率分布Fig. 3. (color online) Focal distributions

25、 of (a) temperature and (b) conductivity for centric HIFU therapy,as well as (c) the conductivity distributions for eccentric HIFU therapy at dierent treatment times.164301-6万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 164301焦域面积不断扩大,在焦平面上形成圆形的焦斑.在t = 2 s时,焦域中心0.2 mm半径范围内的温度超过70 C,达到了治疗效果.在t = 3 s

26、时,焦域半径超过1.3 mm,其中超过在半径0.5 mm的范围内温度超过70 C.进一步将组织的温度-电阻抗关系应用到图3(a)中,得到如图3(b)所示的不同治疗时间的焦平面电导率分布.可见随着的延长,HIFU焦域的温度逐渐升高,电导率相应提高,形成焦域中心导电能力强,周围导电性能弱的梯度分布.为了证明研究的普适性,将焦点沿着径向移动8 mm形成HIFU偏心聚焦,得到了如图3(c)所示的不同治疗时间HIFU焦平面的电导率分布.由于组织媒质相同,模型内不同位置的声学特性和电阻抗特性相同,因此除了焦域位置的移动,焦平面上焦点处的电导率分布和图3(b)基本一致.在HIFU治疗的同时,通过对称电极的电

27、流注入,在模型内部形成相应的电场分布,其变化能反映焦域的温度变化.图4(a)和图4(b)分别显示了HIFU中心和偏心聚焦时,在电极1/17电流注入,t = 1, 2, 3 s时焦平面内的电流密度和电势分布.如图4(a)的虚线环所示,随着治疗时间的延长,焦域的温度和电导率升高,原来较为均匀分布的电流向焦域的中心汇聚,形成向上弯曲的电势等位线.对于如图4(b)所示的偏心聚焦,在焦域内也形成汇聚的电流和弯曲的电势等位线,其弯曲程度由焦域的电导率分布决定.为了进行基于EIT的电阻抗分布重建,采用如图1所示的三角形网格剖分39,其中剖分节点数为1225,剖分单元数为2320,并将有限元计算得到的焦平面电

28、导率分布和和电势分布导入到所建立的剖分网格中,将组织模型的焦平面圆周32等分,1s 2s 3s2.52.01.5 03.03.5 4(a)2.52.03.04.03.51.51s 2s 3s0(b)0 0.75 1.50 2.25 3.00 3.75 4.500 0.75 1.50 2.25 3.00 3.75 4.500 1 2 3 4 50 1 2 3 4 5图4 (网刊彩色)对称电极1/17电流注入, HIFU (a)中心和(b)偏心聚焦时焦平面的电流密度和电势分布Fig. 4. (color online) Distributions of current density and el

29、ectrical potential of the focal plane with (a)centric and (b) eccentric focusing for the current injection from the electrodes 1/17.164301-7万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 1643012 4 6 8 10 12 14 16 18 20 22 24 26 28 30 3200.51.01.52.02.53.07 81.741.761.781.801.821.841.86电压/V电压/V电压/V电压

30、/V电极编号电极编号电极编号电极编号中心1 s中心2 s中心3 s偏心1 s偏心2 s偏心3 s(a)2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 3200.51.01.52.02.53.017 181.581.601.621.641.661.681.70中心1 s中心2 s中心3 s偏心1 s偏心2 s偏心3 s(b)图5 (网刊彩色)对称电极(a) 1/17和(b) 9/25电流注入时,模型焦平面表面电极的电压分布Fig. 5. (color online) Electrical voltage distributions of the surface

31、electrodes in the focal plane with thecurrent injections from the symmetric electrodes of 1/17 and 9/25.获得循环对称电极电流注入时的边界电极测量电压Vij.图5(a)和图5(b)分别显示了对称电极1/17和9/25电流注入时32个边界电极上的测量电压分布.可见随着治疗时间的延长,焦域温度升高,电阻率降低,相同电流注入时边界电极电压降低.在改变电流注入电极后,边界电极的电压分布发生相应的变化.然而,由于相对于模型尺寸,焦域面积很小,电阻抗差异较小,因此电极电压分布随着治疗时间和温度以及位置的变

32、化不明显.图5(a)和图5(b)的内插局部放大图分别显示了两种情况下测量得到电极7/8和17/18的电压分布,可以看出相应电极的电压产生了幅度较小的变化(mV量级).在EIT的电阻抗重建中,根据HIFU治疗中组织焦域的电导率改变,设定初始均匀电导率为0.73 S/m,利用Matlab编程计算焦平面上的电极电位分布Uij( ),建立目标函数,并通过求解每次迭代电阻率分布下的雅克比矩阵和海森矩阵获得新的迭代电阻率分布,当迭代次数达到30次时,获得了较稳定的重建结果.图6(a)和图6(b)分别显示2s 3s1s3s1s 2s2.62.42.22.01.81.61.41.21.00.80.920.90

33、0.880.860.840.820.800.780.760.74(a)(b)图6 (网刊彩色) HIFU中心聚焦时, (a)模型焦平面的网格化电导率分布和(b)重建的电导率分布Fig. 6. (color online) (a) Grid distributions and (b) reconstructed distributions of conductivity in the focalplane for centric focusing.164301-8万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 164301了HIFU中心聚焦时模型焦平

34、面的网格化电导率分布和基于边界电极电压的重建电导率分布,为了进行重建效果比较, HIFU偏心聚焦时的相应结果在图7(a)和图7(b)中给出.可见重建结果能够较为准确地反映出治疗区域的位置和形状以及电导率的变化趋势,能够实现准确的焦点定位,重建焦域的尺寸和电导率随治疗时间的增长而变大,和模型计算结果的变化趋势一致,能够反映出不同治疗时间HIFU焦域的电导率差异.但是由于中心聚焦的焦域具有完全对称性,其边界电极电压的偏差较小,因此重建焦域的精度不如偏心聚焦效果好.由于在电导率的网格化处理中使用了三角形内部的数据平均,所得到的最大电导率小于有限元的仿真结果.同时由于重建算法的病态性,所得到的电导率小

35、于实际值.为了比较重建效果,将重建结果用对应时刻模型电导率的最大值进行归一化处理,得到如图8所示的不同治疗时间模型焦平面的网格化电导率径向分布和重建电导率的径向分布.虽然重建结果比实际电导率范围宽,但是仍然可以精确定位焦点位置,反映焦平面内焦域温度和电导率的分布;随着治疗时间的增加,焦域的电导率变化增大,重建的电导率分布和模型电导率分布趋势一致,能够较好地反映组织电导率随温度的变化特性.3s2s1s3s1s 2s2.62.42.22.01.81.61.41.21.00.81.41.21.00.81.31.10.9(a)(b)图7 (网刊彩色) HIFU偏心聚焦时, (a)模型焦平面的网格化电导

36、率分布和(b)重建的电导率分布Fig. 7. (color online) (a) Grid distributions and (b) reconstructed distributions of conductivity in the focalplane for eccentric focusing.-15 -10 -5 0 5 10 1500.20.40.60.81.0模型 1 s模型 2 s模型 3 s重建 1 s重建 2 s重建 3 s半径/mm归一化的电导率00.20.40.60.81.0归一化的电导率(a)-16 -12 -8 -4 0 4 8 12 16(b)模型 1 s模型

37、 2 s模型 3 s重建 1 s重建 2 s重建 3 s半径 /mm图8 (网刊彩色) HIFU (a)中心和(b)偏心聚焦时,不同治疗时间组织模型焦平面的网格化电导率归一化径向分布和重建电导率的归一化径向分布Fig. 8. (color online) Normalized radial grid distributions and reconstructed distributions of conductivity inthe focal plane for (a) centric and (b) eccentric focusing.164301-9万方数据物理学报Acta Phys.

38、 Sin. Vol. 66, No. 16 (2017) 1643014讨论本研究中由于HIFU的声功率不太高,超声非线性所产生温升对线性声场所产生的误差较小,为了简化温度场的计算,选用黏滞媒质中的线性方程求解,取得较好的计算效果.但实际中, HIFU声场一般为非线性,要获得更为准确的声场和温度场,甚至考虑声空化的影响,在进一步的研究中采用Westervelt40, KZK41;42和SBE43方程求解,会得到更为精确HIFU焦域的温度分布.理论和模拟结果证明,虽然测量的边界电极的电压变化较小,但利用EIT可以精准定位HIFU焦域的位置和焦域内电阻抗的变化趋势,重建图像中电导率的差异可以反映焦

39、域的温度分布特性,在HIFU的疗效监测中有着广泛的应用前景.另外,当前EIT逆问题求解算法主要包括动态重构算法(例如等位线反投影法)和静态重构算法(如牛顿-拉夫逊迭代算法)等.鉴于动态重构算法计算速度快和效率高的优点,本研究首先采用等位线投影法对32电极的测量数据进行了仿真计算,然而由于HIFU焦域的面积小,电导率变化范围大,而边界电极的电压变化小,因此重建算法对电阻抗的变化不敏感,重建图像的精度差.虽然修正的牛顿-拉夫逊算法是一种静态的重建算法,但具有良好的收敛性,重建结果较动态算法更好,但由于该算法重建问题的病态性会随剖分网格的增加而增加,电阻抗重建的网格剖分不能过密,限制了重建精度的提高

40、,因此需要进一步研究EIT重建算法,在较少的网格剖分和表面电极条件下,获得更高精度的重建图像.本研究使用循环对称电极电流注入,其电流密度和电势分布对位置较为敏感,当HIFU焦域在模型中心时,不同位置激励下的表面电极电压的差异较小,因此电阻抗重建的精度和分辨率较低.为了获得较大的表面电极电压差异,需要改进电极注入方式,实现了多种激励模式下的电阻抗分布重建.当前EIT理论较为成熟,算法简便,设备廉价,对生物组织伤害小,重建速度较快,本研究把EIT和HIFU相结合,实现电学系统和声学系统完全分离,对HIFU治疗系统干扰小,可以同步实现HIFU治疗和电阻抗测量,利用组织温度-电导率的关系,可以实现HI

41、FU治疗过程中焦域温度的反演.进一步研究中,将模型对称电极的RIV测量和边界电极电压监测相结合,可以有效提高焦点定位和焦域温度监测的准确性,为HIFU疗效评估和剂量控制提供实时电阻抗测量和成像新方法.5结论本文针对HIFU治疗中的无创温度监测难题,基于组织的温度-电阻抗关系,提出了一种基于EIT的HIFU焦域的温度和疗效监测技术.本研究将HIFU治疗的声学系统和电阻抗测量的电学系统有机结合,利用组织模型表面的电阻抗测量来实现HIFU焦域的电阻抗分布的重建,实现温度和疗效估计.建立了HIFU治疗和EIT测量系统模型,模拟了固定声功率(15.68 W)时组织焦域的声场和温度场分布,以及中心和偏心聚

42、焦的电导率分布.通过对称电极的旋转电流注入,计算了圆柱组织模型焦平面内的电流密度和电势分布,获得了焦平面的表面电极电压分布.在电阻抗重建中,采用修正的牛顿-拉夫逊算法,通过迭代计算实现了焦平面的电导率分布重建.结果表明,重建的电导率分布能准确反映HIFU焦域的位置和范围,体现随治疗时间增长的电导率变化,证明EIT在温度监测中应用的可行性.本研究为HIFU治疗中焦点的精确定位,焦域的实时温度监测和疗效评估以及剂量控制提供了一种无创电阻抗测量和成像新方法.参考文献1 Hutchinson L 2011 Nat. Rev. Clin. Oncol. 8 3852 Kennedy J E 2005 N

43、at. Rev. Canc. 5 3213 Qian S Y, Wang H Z 2001 Acta Phys. Sin. 50 501 (inChinese) 钱盛友,王鸿樟2001物理学报50 5014 Gavrilov L R 2013 J. Acoust. Soc. Amer. 133 43485 Jiang L X, Hu B 2006 Tech. Acoust. 25 43 (in Chinese)姜立新,胡兵2006声学技术25 436 Shen J, Shen J L, Zou J Z 2007 J. Ultrasound Clin. Med.9 486 (in Chinese

44、) 沈洁,申俊玲,邹建中2007临床超声医学杂志9 4867 Ye G, Smith P P, Noble J A 2010 Ultrasound Med. Biol.36 2348 Daniels M J, Varghese T, Madsen E L, Zagzebski J A2007 Phys. Med. Biol. 52 48279 Fan T B, Zhang D, Zhang Z, Ma Y, Gong X F 2008Chin. Phys. B 17 337210 Anand A, Kaczkowski P J 2004 J. Acoust. Soc. Amer.115 249

45、011 Ma Y, Zhang D, Gong X F, Liu X Z, Ma Q Y, Qiu Y Y2007 Chin. Phys. 16 2745164301-10万方数据物理学报Acta Phys. Sin. Vol. 66, No. 16 (2017) 16430112 Fan L Z, Luo F, Yu D Y, Liu X T, Zhang J, Xie M X2005 Clin. J. Med. Instru. 29 115 (in Chinese) 范良志,罗飞,喻道远,刘夏天,张静,谢明星2005中国医疗器械杂志29 11513 Gabriel C, Peyman A,

46、 Grant E H 2009 Phys. Med. Biol.54 486314 ZurbuchenU, Holmer C, Lehmann K S, Stein T, RogganA, Seifarth C, Buhr H J, Ritz J P 2010 Int. J. Hyper-thermia 26 2615 Griths H, Ahmed A 1987 Clin. Phys. Physiol. Meas.8 14716 Cai H, You F S, Shi X T, Fu F, Liu R G, Tang C, DongX Z 2010 Chin. Med. Equip. J.

47、31 8 (in Chinese) 蔡华,尤富生,史学涛,付峰,刘锐岗,汤池,董秀珍2010医疗卫生装备31 817 Su H D, Guo G P, Ma Q Y, Tu J, Zhang D 2017 Chin.Phys. B 26 05430218 Li K Q 2015 M. S. Dissertation (Nanjing: Nanjing Uni-versity of Posts and Telecommunications) (in Chinese)李凯强2015硕士学位论文(南京:南京邮电大学)19 Xu G X 2004 Ph. D. Dissertation (Chongq

48、ing:Chongqing University) (in Chinese) 徐管鑫2004博士学位论文(重庆:重庆大学)20 Zhang L 2013 Ph. D. Dissertation (Nanjing: NanjingUniversity of Science and Technology) (in Chinese) 张丽2013博士学位论文(南京:南京理工大学)21 Curra F P, Mourad P D, Khokhlova V A, Crum L A2000 IEEE Trans. Ultrason. Ferroelect. Freq. Control47 107722 Soneson J E, Myers M R 2010 IEEE Trans. Ultrason.Ferroelect. Freq. Control 57 245023 Myers M R, Soneson J E 2009 J. Acoust. Soc. Amer.126 42524 Soneson J E, Mye