CN109085247A - 一种超声造影剂空化泡群共振状态测量方法 - Google Patents

Abstract

本发明涉及一种超声造影剂空化泡群共振状态测量方法。本发明是在水介质中加入超声造影剂后，首先用高频换能器激励泡群，泡群会因为受到声场的作用产生稳态空化，从而导致测量声波幅值发生变化（从示波器上可观察波形的振动情况），然后用低频换能器激励泡群达到共振状态，此时的声波幅值会发生明显变化（信号剧烈的抖动），所以可通过测量声波幅值变化得到泡群是否达到共振状态。本发明将有助于通过准确的监测超声造影剂泡群的共振状态，在焦域中产生超声治疗所需要的高声压，减少高强度聚焦声波在含微气泡液体中传播时的衰减和色散，这对超声医疗诊断、治疗和实验测量方面有着非常重要的意义。

Description

一种超声造影剂空化泡群共振状态测量方法

技术领域

本发明属于超声空化和声学技术领域的结合，具体涉及一种超声造影剂空化泡群共振状态测量方法。

背景技术

近年来在临床医学超声诊断与生物组织成像领域中，微泡型超声造影剂受到越来越多的关注，超声技术也被应用到多种治疗设备中，国内外大量的研究成果表明，利用声波激励微气泡时所产生的非线性振动和稳态空化，可以提高超声治疗的效率、实施血管内溶栓治疗，而通过微纳气泡携带治疗药物或基因，使用超声作为介导手段，能够进行抗肿瘤药物靶向递送、基因定位转染或递送等方面的治疗，因此，将微气泡和超声结合起来进行一些重大疾病的治疗已成为国内外医学界所关注的热点之一。

当液体介质中含有大量气泡时，入射声波与微气泡振动是一个相互影响的过程，微气泡的存在改变了介质的物理特性，增加了介质的散射、吸收和非线性系数，同时在焦域附近，泡群共振对焦域中声压分布将产生较大影响，会导致泡群中微气泡产生剧烈的非线性振动及空化。为了减少高强度聚焦声波在含气泡液体中传播时的衰减和色散，通过高频聚焦脉冲在焦域中形成稳态空化泡群，进而使用低频聚焦超声激发泡群共振，在焦点处可以形成粉碎结石所需要的高声压。为研究激发泡群共振时所需要的声场条件，分析在共振状态下泡群内部声压分布和微气泡振动规律，获得共振条件下泡群内部声压分布规律及与激励声压幅度、频率之间的关系，所获得的焦域声压分布会影响对治疗机理和治疗效果的判断。所以寻求一种空化泡群共振状态测量方法显得尤为重要。

发明内容

本发明针对现有技术的不足，提供了一种超声造影剂空化泡群共振状态测量方法。

本发明的技术解决方案为：

测量泡群共振状态时需要确定达到气泡空化的条件和共振频率的大小，而超声波能否产生空化现象和空化阈值有很大关系，对于分布函数为f(R₀)的气泡群而言，泡群的出现会对入射声波产生强烈的散射和阻挡作用，由于大量气泡和驱动声场的相互作用，泡群中气泡的振动明显区别于单个气泡的情况，气泡粒径变化也存在不同的状况，其声波幅值不能简单地视为是单个气泡的线性叠加，所有粒径大小气泡形成的声波幅值总和为：

其中，P₀为标准静压力，P_a为入射波声压幅值，圆频率ω＝2πf，f为入射波频率，α'为加入到水介质中超声造影剂的衰减系数，计算方法为：

其中，x表示发射信号和接收信号之间的距离；p和p_x分别表示水和含超声造影剂介质中接收信号幅值；α为已知超声使用频率和温度下纯净液体的衰减系数。R为气泡粒径大小，其随时间t和初始粒径R₀的变化关系式为：

T表示超声造影剂的平均消除半衰期，ψ(R₀,t)表示体积分数随气泡初始粒径R₀和时间t变化的函数，其表达式为：

f(R₀)表示气泡个数分布受粒径大小影响的函数，其表达式为：

k为波数，根据气泡群的散射特性，得到泡群的波数表达式为：

其中，阻尼系数δ＝2πR(R₀,t)/λ＝2πfR(R₀,t)/c，N为气泡的个数，在实验所选区域不太大的情况下，气泡的个数N和体积分数ψ的关系近似满足:

并且满足f₀为气泡的谐振频率，根据气泡动力学得到含气微泡的谐振频率是由气体、液体以及泡膜的性质共同决定的，可表示为：

其中，微泡内气体多方指数γ＝1.07，液体密度ρ＝998kg·m^-3，液体中静压P₀＝1.013×10⁵Pa，液体表面张力系数σ＝0.076N·m^-1，泡膜的弹力系数S_p＝4.2N·m^-1。

根据实验测量需求(高低频换能器一起工作)，得到气泡群达到共振状态时的声波幅值表达式为：

其中，P_a1为加在高频换能器上的入射声波幅值，ω₁＝2πf₁，f₁为入射波高频频率，P_a2为加在低频换能器上的入射声波幅值，ω₂＝2πf₂，f为入射波低频频率。通过高频聚焦脉冲信号在焦域中形成稳态空化泡群，然后使用低频聚焦超声激发泡群发生共振，在焦点处形成高声压，通过程控计算机中的LABVIEW程序可实时显示焦域的声波幅值变化，从而可判断泡群是否达到共振状态。

本发明中的实验装置包括信号发生器、功率放大器、聚焦换能器、高精度三维超声扫描控制机构、探针水听器、数字示波器、程控计算机等。其中，两个聚焦换能器固定在水槽相连的两侧，避免了同时移动换能器和水听器时，导致三者的声轴达不到共声轴的状态，水听器安装在高精度三维超声扫描控制机构上，由程控计算机控制水听器运动，采集不同位置的信号，水听器作为声波接收端与数字示波器相连，两个换能器分别连着两个功率放大器，再连着两个信号发生器，信号发生器的触发信号作为同步信号接入到数字示波器，数字示波器与程控计算机相连，对信号进行读取和存储。

进一步说，所述的两个信号发生器采用的是具备可以同时发送双通道信号功能的普源精电DG4062信号发生器，并可实现同相位操作的信号发生器。

进一步说，所述的换能器采用的一个是直径为900mm，谐振频率为990kHz，聚焦距离为105mm的球壳型换能器，另一个是直径为500mm，谐振频率为100kHz，聚焦距离为50mm的聚焦换能器，两个换能器固定在水槽相连的两端，二者的焦域相对着，避免了同时移动换能器和水听器时，导致三者达不到共声轴的状态。

进一步说，所述的水听器采用的是探针水听器，满足实验需求。

进一步说，所述的放大器是采用的一个是国产的JYH-200中频功率放大器，可放大频率低于400kHz的换能器；一个是进口的RPR-400功率放大器，最高功率输出达8KW，满足实验需求。

进一步说，所述的示波器采用的是普源精电的DS2102双通道示波器，它具备波形触发、内存读取、存储功能，实时采样率最高达2GSa/s，带宽达到200MHz，并具有串口通信，可直接与程控计算机相连，可同时接收两个换能器发射的信号。

进一步说，所述的高精度三维超声扫描控制机构，其***最小的有效控制步进精度为0.02mm，可实现水下三维运动，并可准确实现水下定位。

本发明的有益效果：综上所述，本发明的测量方法，可以针对不同种类的超声造影剂进行测量；水槽的尺寸设计为20cm×20cm×40cm，即保证了探针水听器可以在水槽内进行三维运动，以便找到泡群空化焦域，也保证了超声造影剂的浓度不至于太低，减小微泡向实验区域外扩散，使得微泡在实验测量区域内充分发挥作用，避免了超声造影剂的浪费；水槽上固定的两个高低频换能器即可以单独工作也可以同时工作，单独工作时可以达到稳态空化的状态，同时工作时可以更好的找到气泡群共振的状态，避免了单个换能器激励时的高功率，大大的减小了信号发生器的输入，使得很小的信号激励就能达到共振的状态，实验测量方法可以准确、实时的找到空化泡群共振的状态且操作简便、快捷。该测量方法的提出，将有助于通过使用适当的激励就能准确的监测超声造影剂泡群的共振状态，在焦域中产生超声治疗所需要的高声压，减少高强度聚焦声波在含微气泡液体中传播时的衰减和色散，这对超声医疗诊断、治疗和实验测量方面有着非常重要的意义。

附图说明

图1本发明测量装置示意图。

图2输出脉冲信号暂态稳态示意图。

图中：1、数字示波器，2、水槽，3、程控计算机，4、高精度三维超声扫描控制机构，5、探针水听器，6、聚焦换能器b，7、信号发生器a，8、功率放大器a，9、聚焦换能器a，10、信号发生器b，11、功率放大器b。

具体实施方式

下面结合附图对本发明作进一步的说明：

本发明的测量原理：在水介质中加入超声造影剂后，首先用高频换能器激励泡群，泡群会因为受到声场的作用产生稳态空化，从而导致测量声波幅值发生变化(从示波器上可观察波形的振动情况)，然后用低频换能器激励泡群达到共振状态，此时的声波幅值会发生明显变化(信号剧烈的抖动)，所以可通过测量声波幅值变化得到泡群是否达到共振状态。

如图1所示，首先在水槽2中加入经溶氧脱氧处理过的除气水，然后将配置好的超声造影剂混悬液(SonoVue)加入到水介质中。信号发生器7产生的猝发正弦脉冲信号作为同步信号接入数字示波器1的外部触发，用于捕获声压波形和计算聚焦换能器9和探针水听器5之间的距离，同时利用信号发生器1产生的30个周期、中心频率为990kHz的连续信号通过功率放大器8(频率范围0.5MHz～45MHz，增益50dB)激励聚焦换能器9(中心频率为990kHz，直径为900mm)发射声波，在水槽中形成辐射声场，探针水听器5作为接收器件将接收到声信号转化为电信号，在数字示波器1上进行显示。探针水听器5安装在高精度三维超声扫描控制机构4上，由程控计算机3控制探针水听器5运动，寻找声场空化焦域(尽量保证换能器和水听器共声轴，以增加***的接收灵敏度，提高对微弱信号的接收能力)，然后设置参数观察数字示波器1上的信号变化情况，在声波达到稳态空化后，再用信号发生器10通过功率放大器11激励聚焦换能器6发射信号，观察空化泡群的共振情况，然后将两个信号发生器的信号作同相位处理，在数字示波器上得到空化焦域的声波信号，并通过程控计算机3中的LABVIEW程序串口对数字示波器1上的信号进行读取和存储，得到空化焦域的声波数值，通过对信号的处理，得到不同情况下的空化泡群的共振状态特性。

(1)信号激励部分

由于实验过程中需要为示波器提供与发射信号同步的触发信号，用于捕获声压波形和计算发射波和接收波之间的时延，故选用的信号发生器为普源精电的DG4062双通道信号发生器，它具备双通道可以同时发送信号的功能，满足实验的要求。但是信号发生器的激励电压幅度有限，不足以激励换能器产生空化现象，因此需要对其进行放大。实验选用国产的JYH-200中频功率放大器，对低频换能器器进行激励；还选用RITEC的RPR-4000信号的功率放大器对高频换能器进行激励，它具有自发自收、门限激励、漏电保护等功能，最高功率输出达8KW，满足实验要求。同时为了达到空化气泡共振状态，使用两个型号一样的信号发生器(进行同相位操作)通过放大进行激励换能器。发射信号时为了保护换能器的同时避免产生混响场，激励信号一般采用猝发正弦脉冲信号，它既有脉冲信号性质又有稳定状态，可以在反射信号到来之前完成数据采集，避免形成混响场。

(2)数据采集部分

数据采集部分包括数字示波器和程控计算机两个部分。数字示波器选用的是普源精电的DS2102双通道示波器，它具备波形触发、内存读取、存储功能，实时采样率最高达2GSa/s，带宽达到200MHz，并具有串口通信，可通过软件编程来实现示波器的参数设置、状态查询和读取数据等，满足实验需求，数字示波器作为中转部件将探针水听器接收的信号传输给程控计算机，在计算机上进行显示和存储。

当正弦脉冲激励换能器时，换能器被驱动，开始阶段需要一段时间才能达到稳定的工作状态，当激励信号结束后，换能器还需要进行一段时间的阻尼振动。因此输出的脉冲信号的前部和尾部会有暂态现象出现，如图2所示。诸如有效声压、声工作频率等参数都是按换能器稳态工作状态时定义的，因此信号的暂态部分不可用且激励脉冲的宽度不能小于暂态过程，根据实际情况，实验给了30个周期信号，采集信号时读取稳态部分的周期信号(减小信号的不稳定性)，然后通过对波形的处理得到声波幅值，进而可分析出空化泡群共振状态特性。

(3)运动控制部分

实验中的运动控制部分主要是由计算机里面的LABVIEW程序控制着高精度三维扫描运动控制机构里的步进电机运动，步进电机作为运动的执行机构，接收计算机通过发送的脉冲信号，并将接收的电脉冲转换为自身的角位移，从而带动探针水听器在水槽中进行三维运动，完成各处信号的采集。

综上所述，本发明的测量方法，适合对不同种类的超声造影剂进行测量；水槽上的两个高低频换能器即可以单独工作也可以同时工作，从而更好的达到泡群稳态空化和空化泡群共振状态；水槽的尺寸设计即保证了探针水听器可以在水槽内进行三维运动，以便找到泡群空化焦域，也保证了超声造影剂的浓度不至于太低，减小微泡向实验区域外进行扩散，使得微泡在实验测量区域内充分发挥作用，避免了超声造影剂的浪费。该测量方法的提出，将有助于通过准确的监测超声造影剂泡群的共振状态，在焦域中产生超声治疗所需要的高声压，减少高强度聚焦声波在含微气泡液体中传播时的衰减和色散，这对超声医疗诊断、治疗和实验测量方面有着非常重要的意义。

Claims (1)

1.一种超声造影剂空化泡群共振状态测量方法，其特征在于：

测量泡群共振状态时需要确定达到气泡空化的条件和共振频率大小，对于分布函数为f(R₀)的气泡群而言，泡群的出现会对入射声波产生强烈的散射和阻挡作用，由于大量气泡和驱动声场的相互作用，导致泡群中气泡的振动明显区别于单个气泡的情况，气泡粒径变化也存在不同的状况，其声波幅值不能简单地视为是单个气泡的线性叠加，所有粒径大小气泡形成的声波幅值总和为：

其中，P₀为标准静压力，P_a为入射波声压幅值，圆频率ω＝2πf，f为入射波频率，α'为加入到水介质中超声造影剂的衰减系数：

其中，x表示发射信号和接收信号之间的距离；p和p_x分别表示水和含超声造影剂介质中接收信号幅值；α为已知超声使用频率和温度下纯净液体的衰减系数；R为气泡粒径大小，其随时间t和初始粒径R₀的变化关系式为：

T表示超声造影剂的平均消除半衰期，ψ(R₀,t)表示体积分数随气泡初始粒径R₀和时间t变化的函数，其表达式为：

f(R₀)表示气泡个数分布受粒径大小影响的分布函数，其表达式为：

k为波数，根据气泡群的散射特性，得到泡群的波数表达式为：

其中，阻尼系数δ＝2πR(R₀,t)/λ＝2πfR(R₀,t)/c，N为气泡的个数，气泡的个数N和体积分数ψ的关系近似满足：

并且满足f₀为气泡的谐振频率，根据气泡动力学得到含气微泡的谐振频率是由气体、液体以及泡膜的性质共同决定的，表示为：

其中，微泡内气体多方指数γ＝1.07，液体密度ρ＝998kg·m^-3，液体中静压P₀＝1.013×10⁵Pa，液体表面张力系数σ＝0.076N·m^-1，泡膜的弹力系数S_p＝4.2N·m^-1；

从而得到气泡群达到共振状态时的声波幅值表达式为：

其中，P_a1为加在高频换能器上的入射声波幅值，ω₁＝2πf₁，f₁为入射波高频频率，P_a2为加在低频换能器上的入射声波幅值，ω₂＝2πf₂，f为入射波低频频率；通过高频聚焦脉冲信号在焦域中形成稳态空化泡群，然后使用低频聚焦超声激发泡群发生共振，在焦点处形成高声压，通过程控计算机中的LABVIEW程序实时显示焦域的声波幅值变化，从而可判断泡群是否达到共振状态。