高强度聚焦超声经颅脑肿瘤治疗焦域的仿真研究*

为了研究脑组织和脑肿瘤组织在HIFU治疗时形成焦域的特性及坏死肿瘤组织、不同治疗剂量参数以及多次治疗时焦点间距和时间间隔对HIFU形成焦域的影响,本文以脑胶质瘤患者为例,利用患者头颅CT图像数据建立HIFU经颅治疗的仿真模型,基于Westervelt声波非线性传播方程和Pennes生物热传导方程进行HIFU经颅治疗的仿真研究。结果表明,脑组织和脑肿瘤组织内形成HIFU焦域的差异较小,坏死肿瘤组织对HIFU温度场分布有较小影响;辐照声强越大,焦点温升达到同一温度所用时间越短,焦域长短轴越短,颅骨处温升越低;当焦点间距在一定范围内时,第一次辐照形成的温度场分布对第二个焦点温升达到同一温度所需时间影响较大;两次辐照时间间隔对颅骨处温升和两次聚焦形成60?C以上温度分布影响较小。     

高强度聚焦超声间歇式治疗中焦域温度分布的仿真研究

探究不同治疗深度、组织类型和治疗模式对高强度聚焦超声焦域温度场的影响.采用有限元法建立高强度聚焦超声辐照组织的二维轴对称仿真模型,通过Westervelt方程和Pennes生物热传导方程计算高强度聚焦超声焦域的声场和温度场分布.仿真结果表明:随着治疗深度的增加,焦域内温度逐渐降低,有效温升面积减小;不同组织类型在相同条...     

基于换能器驱动信号特征的高强度聚焦超声焦域损伤实时监测

高强度聚焦超声(high intensity focused ultrasound,HIFU)焦域的实时监测是聚焦超声临床治疗面临的关键问题,目前临床常采用B超图像强回声的变化实现焦域组织损伤的监测,而B超图像出现的强回声大多与焦域处的空化及沸腾气泡相关,无法准确、实时地监测治疗状态.HIFU治疗中焦域组织会伴随温度升高、空化、沸腾和组织特性等变化,换能器表面的声学负载也在不断变化,针对该问题,本文构建了换能器电压电流实时检测平台,通过测量换能器电学参数来感知焦域组织的状态变化.以离体牛肝组织作为HIFU辐照对象,并将相位差变化的结果与离体牛肝组织损伤的结果进行了对照,实验结果表明,在HIFU辐照过程中,换能器电压与电流的相位会出现由相对平稳到大幅波动的过程,此时停止辐照可见焦域出现明显损伤,而此时B超图像灰度无明显变化;此外,当焦域出现空化时,其波动幅度与范围将较之更大.此方法可为HIFU焦域组织损伤监测提供一种新的研究方案和手段.     

基于分数导数研究高强度聚焦超声的非线性声场

在高强度聚焦超声(high intensity focused ultrasound, HIFU) 的研究中, 生物组织的衰减和色散性质会对声能量的空间分布产生影响. 本文提出应用分数导数修正非线性Khokhlov-Zabolotskaya-Kuznetsov (KZK)方程, 研究生物组织中非线性HIFU声场. 对三种生物仿体的衰减和声速色散的理论实验研究表明分数导数应用的可行性, 在此基础上通过数值仿真分析研究了衰减及声速随频率的变化对HIFU焦域分布的影响. 研究结果表明, 在计算强非线性聚焦超声时, 由于高次谐波的强色散作用, 引入分数导数来解决生物组织特殊的衰减以及色散问题可进一步提高HIFU治疗的安全性. 关键词： 分数导数 声衰减 色散 高强度聚焦超声     

HIFU治疗中换能器驱动电功率变化机制及规律*

高强度聚焦超声(HIFU)治疗中的驱动电功率对治疗效率起着非常关键的作用,驱动电功率控制的精准性势必会影响治疗的效率和安全性。前期研究表明：HIFU治疗过程中焦域瞬态物理特性的变化会导致换能器的负载阻抗发生变化,进而影响HIFU驱动电功率,但驱动电功率与焦域瞬态物理特性之间的影响关系及规律尚不明确。该文基于电压、电流传感器、空化检测探头和温度传感器等器件,构建了一种HIFU治疗中驱动电功率实时监测及焦域声空化、温度检测系统。基于该实验研究系统,以离体牛心组织作为HIFU辐照对象,分别研究了HIFU焦域温度变化、声空化及组织损伤与驱动电功率之间的变化关系及规律。研究结果表明：当焦域温度升高时,驱动电功率缓慢上升,驱动电功率与温度变化有良好的相关性;当空化产生时,驱动电功率出现明显的波动;当组织出现损伤时,驱动电功率呈陡然下降的变化。三种情景下,驱动电功率变化有明显区别,这有望为区分HIFU治疗过程中焦域处发生损伤和空化以及实时监测靶组织损伤程度提供一种新的解决方案。     

微泡对高强度聚焦超声焦域影响的研究

微泡对高强度聚焦超声(HIFU)治疗具有增效作用,而HIFU治疗中不同声学条件下微泡对HIFU治疗焦域的影响尚不清楚。本文基于声传播方程、Yang-Church气泡运动方程、生物热传导方程、时域有限差分法(FDTD)、龙格-库塔(RK)法数值仿真研究输入功率、激励频率和气泡初始半径对HIFU在含气泡体模中形成焦域的影响,并利用含Sono Vue造影剂的仿组织体模研究进行实验验证。结果表明,增大输入功率、气泡初始半径和升高激励频率均可增大焦域,随着输入功率的增大,焦域形状可能发生变化,而随着激励频率升高和气泡初始半径的增大,焦域会向远离换能器的方向移动。     

变厚度聚焦换能器对声焦域轴向长度影响的研究*

针对不同厚度的病变组织,改变声焦域轴向长度能提高高强度聚焦超声在临床治疗过程中的安全性和有效性。基于多频超声波叠加原理,该文提出了变厚度(多频)聚焦换能器,并设计了两种类型变厚度聚焦换能器。根据瑞利积分法推导了变厚度聚焦换能器声场,计算和分析了变厚度聚焦换能器的声焦域轴向长度,并与等厚度(单频)聚焦换能器声焦域轴向长度进行对比。结果显示,变厚度聚焦换能器中心到边缘的厚度变化趋势与声焦域轴向长度变化相关,中间薄两边厚换能器声焦域轴向长度缩短,中间厚两边薄换能器声焦域轴向长度变长,且实验验证了理论的正确性。研究结果可为变厚度聚焦换能器声场研究和高强度聚焦超声的临床治疗提供参考。     

格子玻尔兹曼方法对不同张角聚焦声束的建模

高强度聚焦超声(HIFU)是一种新型的无创治疗肿瘤新技术,其中换能器声场数值计算能够为HIFU治疗提供重要的依据。传统非线性KZK和SBE模型广泛应用于换能器声场数值计算,但依然存在某些不足。我们采用一种介观尺度的新型流体力学方法,即格子Boltzmann方法(LBM),基于2维9离散速度(D2Q9)格子构建了轴对称多弛豫参数LBM模型,并通过调节弛豫参数分析其对模型的影响;利用该模型对两个具有不同张角的球面聚焦换能器的声场进行数值模拟,并与KZK和SBE模型的计算结果进行比较。结果表明LBM模型能够很好地描述超声波的激发和传播机制,从流体力学的角度描述聚焦声场的分布,具有清晰的物理意义,且计算过程不受换能器张角的限制,在换能器声场的理论分析和模拟计算及其在HIFU治疗中的应用有着积极的意义。      

基于瑞利积分叠加法的聚焦换能器声场特性研究

基于经典的瑞利积分,提出考虑非线性传播、各次谐波衰减的声场瑞利积分线性叠加算法.分析其原理,并以凹球面聚焦换能器为例,采用该算法研究媒质衰减和非线性传播特性影响声焦域的规律.并与Khokhov-Zabolotskaya-Kuznesov(KZK)数值算法,以及实验结果进行比较,验证瑞利积分线性叠加算法描述高强度聚焦超声(HIFU)声场的有效性.     

时间反演聚焦经颅磁声电刺激仿真与实验研究*

无创脑神经调控技术是生物医学领域的研究热点,经颅磁声电刺激是利用静磁场和声场的耦合而产生的感应电场作用于神经组织,对大脑的目标位置进行刺激和调控的一项技术。颅骨的存在使超声在传播过程中发生相位畸变和幅值衰减,聚焦区域偏离,难以实现精准聚焦。该文基于时间反演法,模拟颅内点声源发射脉冲以及超声传播过程,计算各个阵元接收到的时间差,按照后到先发的原则发射脉冲进行聚焦刺激。与传统相控阵聚焦相比,焦点偏移现象基本得到解决,焦域横向、纵向分辨率均有所提高,提高了声束聚焦精度和感应电场峰值。通过搭建实验平台,将两种聚焦方法所测得的声场归一化处理,验证了时间反演法能补偿焦点偏移,并通过实验证实了超声换能器声场和产生感应电场分布存在较高的一致性。基于真实颅脑结构的虚拟点源时间反演聚焦可以实现无创、精准、灵活的经颅磁声电刺激,有助于推动精准神经调控技术的发展。     

Noninvasive treatment efficacy monitoring and dose control for high-intensity focused ultrasound therapy using relative electrical impedance variation

As an effective therapeutic modality, high-intensity focused ultrasound(HIFU) can destroy tumour tissues by thermocoagulation with less metastasis, but it is still limited by inaccurate non-invasive temperature monitoring and efficacy evaluation. A model of electrical impedance measurement during HIFU therapy was established using the temperatureimpedance relationship. Based on the simulations of acoustic pressure, temperature, and electrical conductivity, the impedance of the phantom was calculated and experimentally demonstrated for different values of acoustic power values and treatment time. We proved that the relative impedance variation(RIV) increases linearly with the increasing treatment time at a fixed acoustic power, and the relative impedance variation rate shows a linear relationship with the acoustic power.The RIV and treatment time required for HIFU treatment efficacy are inversely proportional to the acoustic power and the square of acoustic power, respectively. The favourable results suggest that RIV can be used as an efficient indicator for noninvasive temperature monitoring and efficacy evaluation and may provide new strategy for accurate dose control of HIFU therapy.     

Influence of dynamic tissue properties on temperature elevation and lesions during HIFU scanning therapy: Numerical simulation

When large tumors are treated,ablation of the entire volume of tumors requires multiple treatment spots formed by high intensity-focused ultrasound(HIFU)scanning therapy.The heating effect of HIFU on biological tissue is mainly reflected in temperature elevation and tissue lesions.Tissue property parameters vary with temperature and,in turn,the distribution of temperature as well as the heating effects change accordingly.In this study,an HIFU scanning therapy model considering dynamic tissue properties is provided.The acoustic fields and temperature fields are solved combining the Helmholtz wave equation with Pennes bio-heat transfer equation based on the finite element method(FEM)to investigate the effects of various tissue properties(i.e.,the attenuation coefficient,acoustic velocity,thermal conductivity,specific heat capacity,density,and blood perfusion rate)on heating performance.Comparisons of the temperature distribution and thermal lesions under static and dynamic properties are made based on the data of tissue property parameters varying with temperature.The results show that the dynamic changes of thermal conductivity,specific heat capacity,and acoustic velocity may account for the decrease of temperature elevation in HIFU treatment,while the dynamic changes of attenuation coefficient,density,and blood perfusion rate aggravate the increase of temperature on treatment spots.Compared with other properties,the dynamic change of attenuation coefficient has a greater impact on tissue temperature elevation.During HIFU scanning therapy,the temperature elevation and tissue lesions of the first treatment spot are smaller than those of the subsequent treatment spots,but the temperature on the last treatment spot drops faster during the cooling period.The ellipsoidal tissue lesion is not symmetrical;specifically,the part facing toward the previous treatment spot tends to be larger.Under the condition of the same doses,the temperature elevation and the size of tissue lesions under dynamic properties present significant growth in comparison to static properties.Besides,the tissue lesion begins to form earlier with a more unsymmetrical shape and is connected to the tissue lesion around the previous treatment spot.As a result,lesions around all the treatment spots are connected with each other to form a closed lesion region.The findings in this study reveal the influence of dynamic tissue properties on temperature elevation and lesions during HIFU scanning therapy,providing useful support for the optimization of treatment programs to guarantee higher efficacy and safety.     

斜位电流注入磁共振电阻抗成像（MREIT）与硬件设计

现有的磁共振电阻抗成像(MREIT)方法只适用于主磁场平行水平面的磁共振系统,本文中提出了一种新的MREIT方法即斜位电流注入磁共振电阻抗成像,同时适用于主磁场垂直于水平面的磁共振系统. 硬件系统以磁共振主控计算机为控制中心,将MRI与EIT有机结合于一体：(1)能够测量注入电流在成像体内部感应磁场的磁感应强度B(x,y,z);(2)向成像体注入电流;(3)满足时序匹配要求;(4)测量边界电压. 文章还展示了利用此装置所得的初步实验结果.     

Acoustic power measurement of high intensity focused ultrasound in medicine based on radiation force

How to measure the acoustic power of HIFU is one of the most important tasks in its medical application. In the paper a whole series of formula for calculating the radiation force related to the acoustic power radiated by a single element focusing transducer and by the focusing transducer array were given. Various system of radiation force balance (RFB) to measure the acoustic power of HIFU in medicine were designed and applied in China. In high power experiments, the dependence of radiation force acting the absorbing target on the target position at the beam axis of focusing transducer was fined. There is a peak value of "radiation force" acting the absorbing target in the focal region when the acoustic power through the focal plane exceeds some threshold. In order to avoid this big measurement error caused by the 'peak effect' in focal region, the distance between the absorbing target of RFB and the focusing transducer or transducer array was defined to be equal to or less than 0.7 times of the focal length in the National Standard of China for the measurements of acoustic power and field characteristics of HIFU. More than six different therapeutic equipments of HIFU have been examined by RFB for measuring the acoustic power since 1998. These results show that RFB with the absorbing target is valid in the acoustic power range up to 500W with good linearity for the drive voltage squared of focusing transducer or array. The uncertainty of measurement is within +/-15%.     

微泡对高强度聚焦超声声压场影响的仿真研究*

微泡对高强度聚焦超声(HIFU)治疗焦域具有增效作用,而HIFU治疗中不同声学条件下微泡对HIFU形成声压场的影响尚不清楚。本文基于气液混合声波传播方程、Keller气泡运动方程、时域有限差分(FDTD)法和龙格-库塔(RK)法数值仿真研究输入声压、激励频率、气泡初始空隙率和气泡初始半径对HIFU形成声压场的影响。研究结果表明,随着输入声压的增大,焦点处声压升高但焦点处最大声压与输入声压的比值减小,焦点位置几乎不变;随着激励频率和气泡初始半径的增大,焦点处声压升高且焦点位置向远离换能器方向移动;随着气泡初始空隙率的增大,焦点处声压降低且焦点位置向换能器方向移动。     

球形集声器在生物组织中形成的组织损伤

球形集声器可在亚波长焦域内形成高强度声压,在高强度聚焦超声治疗中具有潜在应用前景.本文结合非线性声传播理论及生物传热学理论,研究球形集声器在生物组织中形成的组织损伤.实验中采用430 kHz,内径为240 mm的球形集声器对肝组织作用,结果表明:集声器表面声压为53 kPa时作用2 s,可以形成小于波长尺度的组织损伤.理论计算结果与实验结果符合得较好,并且理论模型可优化球形集声器的开口孔径.研究结果表明,球形集声器可应用于肿瘤的精细超声治疗.     

Experimental measurement of the acousto-electric interaction signal in saline solution

Acoustic pressure alters local electrical conductivity in tissues and solutions. This work concerns the measurement of electrical conductivity in a liquid which is subjected to an acoustic pressure field created by a focused transducer. Measurements were made with four electrodes positioned in the ultrasonic focal zone, and the signal concerned is referred to as the acousto-electric interaction signal. A solution of sodium chloride in a measurement cell was subjected to ultrasound pressures of upto 1 MPa. It was shown that it is possible to quantitate the acousto-electric interaction signal once the ultrasonic vibration potential due to the Debye effect has been subtracted. The acousto-electric interaction signal was shown to be directly proportional to both the applied acoustic pressure and current. For the measurement cell used in this work, the interaction factor was found to be 5.3 microVmA(-1) MPa(-1).