造影剂微泡声学特性测量方法、装置、设备及存储介质制造方法及图纸

本发明专利技术公开了造影剂微泡声学特性测量方法、装置、设备及存储介质，其中方法包括：标定阶段，获取发射探头发射的脉冲信号和接收探头接收的回波信号，并利用脉冲信号和回波信号计算得到接收探头的接收传递函数；测量阶段，控制发射探头在预设频率范围内扫频以发射超声激励信号，且控制接收探头接收仿体血管中微泡的散射信号；从散射信号中提取得到次谐波信号电压幅度；利用接收探头的接收传递函数将次谐波信号电压幅度转换为不同频率下的实际次谐波散射压；根据实际次谐波散射压确认次谐波最优驱动频率和次谐波散射压力大小。本发明专利技术通过获取接收探头的接收传递函数，并利用其获取实际的散射声压信号，从而提高微泡非线性散射特性的测量准确性。性的测量准确性。性的测量准确性。

【技术实现步骤摘要】
造影剂微泡声学特性测量方法、装置、设备及存储介质


[0001]本申请涉及超声造影剂特性测量
，特别是涉及一种造影剂微泡声学特性测量方法、装置、设备及存储介质。

技术介绍

[0002]近年来在临床医学超声诊断与生物组织成像的领域中，微泡型超声造影剂受到越来越多的关注，超声技术也被应用到多种治疗设备中，国内外大量的研究成果表明，利用声波激励微气泡时所产生的非线性振动和散射特性，可以提高超声治疗的效率、实施血管内溶栓治疗，而通过微纳气泡携带治疗药物或基因，使用超声作为介导手段，能够进行抗肿瘤药物靶向递送、基因定位转染或递送等方面的治疗，因此，将微气泡和超声结合起来进行一些重大疾病的治疗已成为国内外医学界所关注的热点之一。
[0003]目前，衡量超声造影剂微泡声学特性主要是通过衰减系数和背向散射系数两个重要指标进行评估，而由于探头的传递函数通常难以获得，因此一般通过将样本散射信号的功率谱除以参考信号的功率谱，以消除探头传递函数的影响，从而得到超声造影剂微泡的衰减系数和背向散射系数。但是，当声波在介质中传播时，气泡的受迫振动会引起声散射，体现强的非线性声特性，而衰减系数和背向散射系数仅能反映低声压下微泡的线性散射特性(基波共振频率和散射强度)，而不能反映微泡在高声压下的非线性散射特性，如二次谐波或者次谐波成分的共振频率、散射强度等。由于没有获得探头的传递函数，而无法得到实际的散射声压信号，这导致较强的频率成分可能并不在探头的中心频率，使得探头对该频率点的响应比较弱，而将其他频率成分误认为微泡的非线性共振频率，最终导致对超声造影剂的非线性声学特性评估不准确。

技术实现思路

[0004]有鉴于此，本申请提供一种造影剂微泡声学特性测量方法、装置、设备及存储介质，以解决现有超声造影剂声学特性评估不准确、不全面的问题。
[0005]为解决上述技术问题，本申请采用的一个技术方案是：提供一种造影剂微泡声学特性测量方法，包括：在标定阶段，获取发射探头发射的脉冲信号和接收探头接收的回波信号，并利用脉冲信号和回波信号计算得到接收探头的接收传递函数；在测量阶段，控制发射探头在预设频率范围内扫频以发射超声激励信号，且控制接收探头接收仿体血管中微泡的散射信号；从散射信号中提取得到次谐波信号电压幅度；利用接收探头的接收传递函数将次谐波信号电压幅度转换为不同频率下的实际次谐波散射压；根据实际次谐波散射压确认次谐波最优驱动频率和次谐波散射压力大小。
[0006]作为本申请的进一步改进，在标定阶段，获取发射探头发射的脉冲信号和接收探头接收的回波信号，并利用脉冲信号和回波信号计算得到接收探头的接收传递函数，包括：获取发射探头发射的多个脉冲信号，并对多个脉冲信号进行求和平均，得到平均脉冲信号；对平均脉冲信号进行快速傅里叶变化，得到电压信号幅度谱；将电压信号幅度谱转换为声
压信号幅度谱；根据声压信号幅度谱计算得到发射探头表面的第一声压谱；获取接收探头接收的多个回波信号，并对多个回波信号进行求和平均，得到平均回波信号；对平均回波信号进行快速傅里叶变化，得到平均回波信号的复数谱；根据发射探头表面的第一声压谱计算得到回波信号在接收探头表面产生的声压的第二声压谱；根据平均回波信号的复数谱和第二声压谱计算得到接收探头的接收传递函数。
[0007]作为本申请的进一步改进，第一声压谱的计算过程表示为：根据经典的声学理论，t时刻，声场中一点r处的瞬时声压为：t时刻，声场中一点r处的瞬时声压为：其中，ρ为介质的密度，为速度势，通过引入一个δ函数的卷积，则该阵元在点r处产生的速度势为其中，阵元的空间脉冲响应函数为：其中，S为阵元的面积，v0(t)为表面法向速度，r
′
为空间点r与阵元上的微元dS之间的距离，c介质中的声速，则其中，声场中一点的声压为：其中，声场中一点的声压为：为声场中一位置处质点的速度，则发射探头表面的声压为：聚焦发射探头焦点处的空间脉冲响应函数：其中，为聚焦探头的深度，R为发射探头的曲率半径，a表示发射探头的半径，表示傅里叶变换，对进行快速傅里叶变换得到：进行快速傅里叶变换得到：其中，f表示频率，j表示复数，ω＝2πf表示角频率表示，表示，表示；第二声压谱的计算过程表示为：其中，P0()为第一声压谱，D
ref
(2
ref
,)为从发射探头表面到参考平面再回到探头表面的声耦合函数，表示接收探头表面的平均声压，当接收探头为聚焦探头时，Z
ref
＝，则有：D
ref
(2,ω)＝
‑
{1
‑
exp(
‑
jG
p
)[J0(G
p
)+jJ1(
p
)]}，其中，J0、J1分别为0阶和1阶的贝塞尔函数，G
p
＝a2/2为声压在焦点上的增益，表示，exp表示指数函数，表示波数表示，当G
p
>时，有：
[0008]作为本申请的进一步改进，接收传递函数的计算公式表示为：其中，T
R
(f)为接收传递函数，V
out
(f)为平均回波信号的复数谱，P
R,0
(f)为第二声压谱。
[0009]作为本申请的进一步改进，在标定阶段，获取发射探头发射的脉冲信号和接收探头接收的回波信号，并利用脉冲信号和回波信号计算得到接收探头的接收传递函数，包括：获取发射探头发射的多个脉冲信号，并对多个脉冲信号进行求和平均，得到平均脉冲信号；对平均脉冲信号进行快速傅里叶变化，得到电压信号幅度谱；将电压信号幅度谱转换为声压信号幅度谱；根据声压信号幅度谱计算得到发射探头表面的第一声压谱；获取施加于发射探头表面的激励电压的复数谱；根据第一声压谱和激励电压幅度谱计算得到发射探头的发射传递函数；当发射探头和接收探头为超声换能器时，根据发射传递函数和预先获取的
互易常数计算得到接收探头的接收传递函数，互易常数为超声换能器的参数。
[0010]作为本申请的进一步改进，发射探头的发射传递函数的计算公式表示为：其中，T
T
(f)为发射传递函数，P
T,0
(f)为第一声压谱，V
in
(f)为发射探头表面的激励电压的复数谱；接收传递函数的计算公式表示为：其中，R
load
为连接到探头上的负载的电阻，A为探头的有效表面积，ρ为水的密度，c为水中的声速，const为互易常数。
[0011]作为本申请的进一步改进，在测量阶段，控制发射探头在预设频率范围内扫频以发射超声激励信号，且控制接收探头接收仿体血管中微泡的散射信号，包括：获取多组预设的声压和压力环境参数；分别在不同的声压和压力环境参数条件下，控制发射探头在预设频率范围内扫频以发射超声激励信号，且控制接收探头接收仿体血管中微泡的散射信号。
[0012]为解决上述技术问题，本申请采用的又一个技术方案是：提供一种造影剂微泡声学特性测量装置，包括：标定模块，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种造影剂微泡声学特性测量方法，其特征在于，所述方法包括：在标定阶段，获取发射探头发射的脉冲信号和接收探头接收的回波信号，并利用脉冲信号和回波信号计算得到接收探头的接收传递函数；在测量阶段，控制所述发射探头在预设频率范围内扫频以发射超声激励信号，且控制所述接收探头接收仿体血管中微泡的散射信号；从所述散射信号中提取得到次谐波信号电压幅度；利用所述接收探头的接收传递函数将所述次谐波信号电压幅度转换为不同频率下的实际次谐波散射压；根据所述实际次谐波散射压确认次谐波最优驱动频率和次谐波散射压力大小。2.根据权利要求1所述的造影剂微泡声学特性测量方法，其特征在于，所述在标定阶段，获取发射探头发射的脉冲信号和接收探头接收的回波信号，并利用脉冲信号和回波信号计算得到接收探头的接收传递函数，包括：获取所述发射探头发射的多个脉冲信号，并对所述多个脉冲信号进行求和平均，得到所述平均脉冲信号；对所述平均脉冲信号进行快速傅里叶变化，得到电压信号幅度谱；将所述电压信号幅度谱转换为声压信号幅度谱；根据所述声压信号幅度谱计算得到所述发射探头表面的第一声压谱；获取接收探头接收的多个回波信号，并对所述多个回波信号进行求和平均，得到平均回波信号；对所述平均回波信号进行快速傅里叶变化，得到平均回波信号的复数谱；根据所述发射探头表面的第一声压谱计算得到所述回波信号在所述接收探头表面产生的声压的第二声压谱；根据所述平均回波信号的复数谱和所述第二声压谱计算得到所述接收探头的接收传递函数。3.根据权利要求2所述的造影剂微泡声学特性测量方法，其特征在于，所述第一声压谱的计算过程表示为：根据经典的声学理论，t时刻，声场中一点r处的瞬时声压为：其中，ρ为介质的密度，为速度势，通过引入一个δ函数的卷积，则该阵元在点r处产生的速度势为度势为其中，阵元的空间脉冲响应函数为：其中，S为阵元的面积，v0(t)为表面法向速度，r
′
为空间点r与阵元上的微元dS之间的距离，c介质中的声速，则其中，声场中一点的声压为：其中，声场中一点的声压为：为声场中一位置处质点的速度，则所述发射探头表面的声压为：聚焦发射探头焦点处的空间脉冲响应函数：其中，为聚焦探头的深度，R为发射探头的曲率半径，a表示发射探头的半径，表示傅里叶变换，对为聚焦探头的深度，R为发射探头的曲率半径，a表示发射探头的半径，表示傅里叶变换，对进行快速傅里叶变换得到：
其中，f表示频率，j表示复数，ω＝2πf表示角频率；所述第二声压谱的计算过程表示为：所述第二声压谱的计算过程表示为：其中，P0(ω)为所述第一声压谱，D
ref
(2Z
ref
,ω)为从发射探头表面到参考平面再回到探头表面的声耦合函数，表示接收探头表面的平均声压，当接收探头为聚焦探头时，Z
ref
＝R，则有：D
ref
(2R,ω)＝
‑
{1
‑
exp(
‑
jG
p
)[J0(G
p
)+jJ1(G
p
)]}，其中，J0、J1分别为0阶和1阶的贝塞尔函数，G
p
＝ka2/2R为声压在焦点上的增益，表示，exp表示指数函数，表示波数，当G
p
...

【专利技术属性】
技术研发人员：郑海荣，李飞，黄来鑫，孟龙，
申请(专利权)人：中国科学院深圳先进技术研究院，
类型：发明
国别省市：