CN103776524A - 一种测量强聚焦超声非线性声场分布的系统及其方法 - Google Patents

一种测量强聚焦超声非线性声场分布的系统及其方法 Download PDF

Info

Publication number
CN103776524A
CN103776524A CN201410059590.5A CN201410059590A CN103776524A CN 103776524 A CN103776524 A CN 103776524A CN 201410059590 A CN201410059590 A CN 201410059590A CN 103776524 A CN103776524 A CN 103776524A
Authority
CN
China
Prior art keywords
partiald
sigma
theta
sound
shell transducer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410059590.5A
Other languages
English (en)
Other versions
CN103776524B (zh
Inventor
范庭波
陈涛
章东
胡济民
张崴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University filed Critical Nanjing University
Priority to CN201410059590.5A priority Critical patent/CN103776524B/zh
Publication of CN103776524A publication Critical patent/CN103776524A/zh
Application granted granted Critical
Publication of CN103776524B publication Critical patent/CN103776524B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种测量强聚焦超声非线性声场分布的系统及其方法,属于超声声场测量技术领域。本发明步骤为:一、获得球壳换能器轴向及其焦平面径向的声压分布曲线;二、根据瑞利积分计算球壳换能器轴向及焦平面径向的声压分布曲线,将计算结果与测量结果比较,获得有效参数;三、将有效参数代入非线性声传播模型,计算R21,再改变P0,计算不同P0对应的R21,得R21与P0的关系曲线;四、对球壳换能器施加激励,测量焦点波形,计算得R21′,插值得到R21′对应的P0′;五、将P0′带入非线性声传播模型,得到非线性声场分布。本发明有效解决了水听器在测量强非线性声场带宽不足的问题,使用相对值R21作为判断依据,消除了水听器灵敏度不高带来的测量隐患。

Description

一种测量强聚焦超声非线性声场分布的系统及其方法
技术领域
本发明涉及超声声场测量技术领域,更具体地说,涉及一种测量强聚焦超声非线性声场分布的系统及其方法。
背景技术
近十多年来,高强度聚焦超声(HIFU)作为一种新兴的非侵入式治疗肿瘤方法得到人们极大的关注。HIFU治疗肿瘤的主要原理是将体外超声能量聚焦于体内肿瘤组织内,靶区组织吸收声能量,温度在短时间内上升至65℃以上,发生不可逆凝固性坏死。为确保治疗的安全性和有效性,治疗前必须准确描述HIFU声场分布。由于HIFU治疗时声能量很高,焦点声压甚至能达到100MPa量级;强非线性导致声波波形严重畸变,甚至形成冲击波,频谱展宽。强声压和宽频带是HIFU声场测量中的两个挑战,要求声场测量系统既要能承受强声压又要有很宽的频带,而目前一直没有理想的声场测量系统能够满足上述两点。
传统的测量声场方法(例如使用PVDF水听器测量)无法满足HIFU声场的极端条件,目前比较常用的方法是使用光纤水听器,利用光纤端口的声-光效应(如压力变化导致光反射系数变化)测量HIFU声场。虽然光纤水听器灵敏度不如传统PVDF水听器高且噪声较大,但是其频带宽且鲁棒性好。然而,在实际测量中发现,当冲击波产生后,光纤水听器带宽仍然不能满足测量要求,测得的峰值正压常常小于实际值。此外,研究者们还提出了一些非侵入式的通过监测HIFU声场中质点位移的测量声场技术,如通过磁共振MRI(magneticresonance imaging)、粒子速度成像PIV(particle image velocimetry)等,但这些方法还很不成熟,且造价昂贵。
另一方面,HIFU非线性声传播模型目前已经较为成熟。Westervelt方程描述了精确到二阶的非线性声场。对于平面和弱聚焦换能器(半张角小于16°),在近轴近似下,Westervelt方程可以化简为著名的Khokhlov–Zabolotskaya–Kuznetsov(KZK)方程;对于强聚焦换能器,在椭球坐标系下,Westervelt方程可以化简为Spheroidal beam equation(SBE)模型。在已知边界条件,即换能器表面声压分布后,非线性声场模型可以准确的预测HIFU声场分布。
目前,有研究者提出了结合实验测量与数值模型来确定高强度聚焦超声非线性声场的方法。Canney等人于2008年提出了该方法:首先根据线性声场测量和数值计算结果得到换能器有效参数;其次根据测量的焦点声压和换能器聚焦增益得到换能器表面声压P0;假设激励换能器的电信号V0与换能器表面声压P0成线性关系,那么在任意电信号V激励下,换取表面声压p=P0(V/V0);最后将表面声压p作为边界条件,可以计算得到在该激励下的声场分布。然而,此方法的缺点在于当激励电信号很强时,电信号与换能器表面声压的线性关系无法得到保证。
发明内容
1.发明要解决的技术问题
本发明的目的在于克服现有高强度聚焦超声非线性声场的测量方法,在强非线性条件下,测量系统带宽不够,导致声场测量准确度不高的问题,提供了一种测量强聚焦超声非线性声场分布的系统及其方法,本发明提供的技术方案测量准确度高、容易实施,测量过程中对水听器的灵敏度要求较低,避免了因水听器灵敏度低导致检测结果不准确的隐患。
2.技术方案
为达到上述目的,本发明提供的技术方案为:
本发明的一种测量强聚焦超声非线性声场分布的方法,其测量步骤为:
步骤一、在线性声场条件下,测量球壳换能器轴向以及该球壳换能器焦平面径向的声压,获得球壳换能器轴向及其焦平面径向的声压分布曲线;
步骤二、根据瑞利积分计算球壳换能器轴向及该球壳换能器焦平面径向的声压分布曲线,将计算结果与步骤一获得的测量结果相比较,调节球壳换能器的张口半径a和几何焦距d,使得计算获得的球壳换能器轴向及其焦平面径向声压分布曲线与步骤一的测量结果相吻合,获得球壳换能器有效参数a和d;
步骤三、将步骤二获得的有效参数代入非线性声传播模型,计算焦点的二次谐波与基波的比值R21,再改变球壳换能器表面声压幅值P0,计算不同表面声压幅值P0对应的R21,得R21与P0的关系曲线;
步骤四、对球壳换能器施加单频电信号激励,测量焦点波形,对该焦点波形进行傅里叶变换,计算得R21′,根据步骤三得到的R21与P0的关系曲线进行插值,得到R21′对应的球壳换能器表面声压P0′;
步骤五、将步骤四获得的表面声压P0′带入步骤三所述的非线性声传播模型,计算得到步骤四施加的电信号激励对应的非线性声场分布。
更进一步地,步骤一测量声压的条件为:对球壳换能器施加低能量激励,使得球壳换能器焦点处基波幅度与二次谐波幅度的差值不低于40dB,忽略声传播过程中的非线性效应,将声场近似为线性声场。
更进一步地,步骤二涉及的瑞利积分公式为:
p ( M ) = - jk 2 π ∫ ∫ Source p 0 ( N ) exp ( - jkr ) r dS
式中,p(M)为声场空间任意一点M的声压,p0(N)为球壳换能器表面任意一点N处声压;k=ω0/c0为波数,其中,ω0=2πf0为声波的角频率,f0为声波的频率,c0为声速;r为M、N两点间距离,Source表示球壳换能器表面区域。
更进一步地,步骤二要求计算获得的球壳换能器轴向及其焦平面径向声压分布曲线与步骤一测量结果的吻合程度为:测量结果与理论计算结果的主瓣偏差低于5%。
更进一步地,步骤三在椭球坐标系下计算声场分布,椭球坐标系中一点(σ,η,φ)与直角坐标系(x,y,z)的变换关系为:
Figure BDA0000468404390000032
z=bση
式中,-∞<σ<+∞,0≤η≤1,
Figure BDA0000468404390000033
2b为椭球焦距;为简化声传播方程,作变换η=cosθ,0≤θ≤π/2;
所述的非线性声传播模型为SBE模型,即:
&PartialD; 2 P &PartialD; &tau; s &PartialD; &sigma; + 1 2 sin 2 &theta; &sigma; ( 1 + &sigma; 2 ) &PartialD; 2 P &PartialD; &tau; s &PartialD; &theta; + &epsiv; &sigma; 2 + sin 2 &theta; &sigma; ( 1 + &sigma; 2 ) ( &PartialD; 2 P &PartialD; &theta; 2 + cot &theta; &PartialD; P &PartialD; &theta; ) + E &sigma; &PartialD; P &PartialD; &tau; s = - &sigma; 2 + sin 2 &theta; &sigma; ( &alpha;b &PartialD; 3 P &PartialD; &tau; s 3 + b 2 l D &PartialD; 2 P 2 &PartialD; &tau; s 2 ) E , ( &sigma; < &sigma; 0 < 0 )
&PartialD; 2 P &PartialD; &tau; p &PartialD; &sigma; - &sigma; sin 2 &theta; 1 + &sigma; 2 &PartialD; 2 P &PartialD; &tau; p &PartialD; &theta; - &epsiv; ( 2 - cos &theta; ) 1 + &sigma; 2 ( &PartialD; 2 P &PartialD; &theta; 2 + cot &theta; &PartialD; P &PartialD; &theta; ) = ( &alpha;b &PartialD; 3 P &PartialD; &tau; p 3 + b 2 l D &PartialD; 2 P 2 &PartialD; &tau; p 2 ) E , ( &sigma; &GreaterEqual; &sigma; 0 )
其中,σ<σ0<0部分为声场靠近声源表面的球面波近似部分,σ≥σ0部分为声场远离声源表面的平面波近似部分,σ0为球面波区域与平面波区域分界面;式中,P=p/P0为归一化声压,p为点
Figure BDA0000468404390000036
处声压,P0为球壳换能器表面声压;τs=ω0[t+b(σ2+sin2θ)1/2/c0]为球面波弛豫时间,τp=ω0[t-bσcosθ/c0]为平面波弛豫时间,t为球壳换能器激发声波在传播介质中的传播时间;E=(σ2+cos2θ)/(1+σ2);α为声传播介质声衰减系数;lD=ρ0c0 3/(βP0ω0)为冲击波形成距离,ρ0为声传播介质密度,β为声传播介质非线性系数,ρ0,α,β均为常数。
更进一步地,步骤三所用的SBE模型的边界条件为:
P ( &sigma; = - &sigma; max , &theta; &tau; s ) = sin &tau; s &theta; &le; &theta; 0 0 &theta; > &theta; 0
式中,σmax=d/b, &theta; 0 = tan - 1 ( tan &alpha; 0 1 + 1 / &sigma; max 2 ) , &alpha; 0 = sin - 1 ( a d ) .
本发明的一种测量强聚焦超声非线性声场分布的系统,包括示波器、计算机、信号发生器、功率放大器、三维运动平台、水听器、球壳换能器和水槽,所述的示波器、计算机、信号发生器、功率放大器和球壳换能器依次相连,功率放大器将信号发生器的输出信号放大,驱动球壳换能器激发声场;所述的球壳换能器固定于水槽内部,水槽内盛装有水;所述的三维运动平台控制水听器移动,测量球壳换能器轴向及该球壳换能器焦平面径向声压分布;所述的示波器采集水听器的输出信号;所述的计算机用于调控整个测量系统,计算机完成的功能包括:控制信号发生器输出,控制三维运动平台运动以及存储数据。
3.有益效果
采用本发明提供的技术方案,与已有的公知技术相比,具有如下显著效果:
(1)本发明的一种测量强聚焦超声非线性声场分布的方法,通过引入焦点处二次谐波与基波比值R21反映换能器的边界条件,再根据非线性声传播模型,得到对应激励下的声场分布;由于仅需要测量二次谐波和基波,有效的解决了强非线性条件下,测量系统带宽不够的问题;此外由于R21为相对值,测量过程中对水听器的灵敏度要求较低,避免了因水听器灵敏度低导致检测结果不准确的隐患;
(2)本发明的一种测量强聚焦超声非线性声场分布的系统,结构简单、便于控制,测量准确度高、容易实施。
附图说明
图1是本发明中强聚焦超声非线性声场分布的测量流程图;
图2是本发明的一种测量强聚焦超声非线性声场分布的系统结构示意图;
图3是本发明进行实例测量使用的强聚焦球壳换能器的示意图;
图4是本发明通过测量得到的球壳换能器线性声场轴向分布和焦平面径向分布以及数值仿真结果图;
图5是本发明根据SBE模型计算得到的R21与球壳换能器表面声压P0的关系曲线图;
图6是本发明球壳换能器表面声压与信号发生器输出信号的关系曲线图;
图7中的(a)~(f)分别是本发明在三种激励下测量得到的以及使用本发明的测量方法得到的焦点波形、相应频谱对比图;
图8中的(a)~(c)分别是本发明在三种激励下球壳换能器焦平面上径向声压分布图。
示意图中的标号说明:
1、示波器;2、计算机;3、信号发生器;4、功率放大器;5、三维运动平台;6、水听器;7、球壳换能器;8、水槽。
具体实施方式
为进一步了解本发明的内容,结合附图和实施例对本发明作详细描述。
实施例1
结合图2,本实施例测量强聚焦超声非线性声场分布的系统,包括示波器(Agilent54830B)1、计算机2、信号发生器(Agilent33250A)3、功率放大器(A150,ENI,Rochester)4、三维运动平台(Newport ESP7000)5、水听器(HNA-0400,ONDA或者HFO-660,ONDA)6、球壳换能器7和水槽8,所述的示波器1、计算机2、信号发生器3、功率放大器4和球壳换能器7依次相连,功率放大器4将信号发生器3的输出信号放大,驱动球壳换能器7激发声场;所述的球壳换能器7固定于水槽8内部,水槽8内盛装有水;所述的三维运动平台5控制水听器6在三维空间(x,y,z)内移动,测量球壳换能器7轴向及该球壳换能器7焦平面径向声压分布。所述的示波器1采集水听器6的输出信号。所述的计算机2用于调控整个测量系统,完成功能包括:控制信号发生器3输出,控制三维运动平台5运动以及存储数据。
本实施例的一种测量强聚焦超声非线性声场分布的方法(测量流程参见图1),其测量步骤为:
步骤一、给球壳换能器7施加低能量的激励,使得焦点处二次谐波幅度比基波幅度至少低40dB,此时高次谐波幅度远小于基波幅度,忽略声传播中的非线性效应,将声场看作线性声场。将球壳换能器7固定在水槽8中,控制三维运动平台5带动水听器6逐点测量球壳换能器7轴向以及该球壳换能器7焦平面(焦平面即指垂直于球壳换能器7对称轴且经过焦点的平面)径向的声压,获得球壳换能器7轴向及其焦平面径向的声压分布曲线。
步骤二、根据瑞利积分计算球壳换能器7轴向及该球壳换能器7焦平面径向的声压分布曲线,涉及的瑞利积分公式为:
p ( M ) = - jk 2 &pi; &Integral; &Integral; Source p 0 ( N ) exp ( - jkr ) r dS
式中,j为单位虚数,p(M)为声场空间任意一点M的声压,p0(N)为球壳换能器7表面任意一点N处声压;k=ω0/c0为波数,其中,ω0=2πf0为声波的角频率,f0为声波的频率,c0为声速;r为M、N两点间距离,dS为球壳换能器7表面面积元,Source表示球壳换能器7表面区域,该表面区域级指球壳换能器7用于电-声转化的的压电陶瓷部分。
将计算结果与步骤一获得的测量结果相比较,一般情况下,由于球壳换能器7封装等原因,球壳换能器7表面振动并不是理想的同相位同振幅振动,这导致根据球壳换能器7标称参数计算得到的声场与实际测量声场有差异。微调球壳换能器7的张口半径a和几何焦距d(参看图3),使得计算获得的球壳换能器7轴向及其焦平面径向声压分布曲线与步骤一的测量结果相吻合,吻合程度要求主瓣的偏差低于5%,此时的球壳换能器7参数a和d为球壳换能器7的有效参数。
结合图3,本实施例中球壳换能器7的标称几何尺寸张口半径a为5.0cm,几何焦距d为10.0cm。球壳换能器7中心频率为1.1MHz。利用本实施例球壳换能器7有效参数的确定方法,信号发生器3输出10mv信号,经功率放大器4放大后驱动球壳换能器7,此时由于驱动信号很小,球壳换能器7激发的声场满足线性声场条件,测量球壳换能器7轴向和焦平面径向声压分布。根据球壳换能器7标称参数计算轴向及焦平面径向声场分布,与测量结果比较,发现焦点位置与测量结果不符合,调整球壳换能器7几何焦距为10.3cm,再次计算轴向和焦平面径向声压分布,此时计算结果与实验结果吻合很好。因此球壳换能器7有效参数为:张口半径5.0cm,几何焦距10.3cm。如图4所示,实际测量的球壳换能器7轴向声压和焦平面径向声压分布与计算结果吻合很好。
步骤三、将步骤二获得的有效参数代入非线性声传播模型,计算焦点的二次谐波与基波的比值R21,再改变球壳换能器7表面声压幅值P0,计算不同表面声压幅值P0对应的R21,得R21与P0的关系曲线。具体为:对于本实施例使用的大张角强聚焦球壳换能器7,在椭球坐标系下计算声场分布。椭球坐标系为以正交的椭球和双曲面为坐标轴的坐标系,椭球坐标系中一点(σ,η,φ)与直角坐标系(x,y,z)的变换关系为:
Figure BDA0000468404390000061
z=bση
式中,-∞<σ<+∞,0≤η≤1,
Figure BDA0000468404390000062
σ代表椭球面坐标,2b为椭球焦距,η代表双曲面坐标,
Figure BDA0000468404390000063
代表原点与点(x,y,z)的连线在x-y平面的投影与x轴的夹角;为简化声传播方程,作变换η=cosθ,0≤θ≤π/2。大张角强聚焦球壳换能器7的声场将分为两个部分:靠近球壳换能器7表面的球面波近似部分和远离球壳换能器7表面的平面波近似部分。椭球坐标系下非线性声传播方程,即Spheroidal beam equation(SBE)模型表示为:
&PartialD; 2 P &PartialD; &tau; s &PartialD; &sigma; + 1 2 sin 2 &theta; &sigma; ( 1 + &sigma; 2 ) &PartialD; 2 P &PartialD; &tau; s &PartialD; &theta; + &epsiv; &sigma; 2 + sin 2 &theta; &sigma; ( 1 + &sigma; 2 ) ( &PartialD; 2 P &PartialD; &theta; 2 + cot &theta; &PartialD; P &PartialD; &theta; ) + E &sigma; &PartialD; P &PartialD; &tau; s = - &sigma; 2 + sin 2 &theta; &sigma; ( &alpha;b &PartialD; 3 P &PartialD; &tau; s 3 + b 2 l D &PartialD; 2 P 2 &PartialD; &tau; s 2 ) E , ( &sigma; < &sigma; 0 < 0 )
&PartialD; 2 P &PartialD; &tau; p &PartialD; &sigma; - &sigma; sin 2 &theta; 1 + &sigma; 2 &PartialD; 2 P &PartialD; &tau; p &PartialD; &theta; - &epsiv; ( 2 - cos &theta; ) 1 + &sigma; 2 ( &PartialD; 2 P &PartialD; &theta; 2 + cot &theta; &PartialD; P &PartialD; &theta; ) = ( &alpha;b &PartialD; 3 P &PartialD; &tau; p 3 + b 2 l D &PartialD; 2 P 2 &PartialD; &tau; p 2 ) E , ( &sigma; &GreaterEqual; &sigma; 0 )
上式中,声场分为靠近球壳换能器7表面(即σ<σ0<0)的球面波近似部分,和远离球壳换能器7表面(即σ≥σ0)的平面波近似部分,σ0为球面波区域与平面波区域分界面。对于轴对称的球壳换能器7,声场分布与无关;P=p/P0为归一化声压,p为点处声压,P0为球壳换能器7表面声压;τs=ω0[t+b(σ2+sin2θ)1/2/c0]为球面波弛豫时间,τp=ω0[t-bσcosθ/c0]为平面波弛豫时间,ω0为声波的角频率,t为球壳换能器激发声波在传播介质中的传播时间,c0为声传播速度;ε=1/(2kb),k=ω0/c0为波数;E=(σ2+cos2θ)/(1+σ2);α为声传播介质的声衰减系数;lD=ρ0c0 3/(βP0ω0)为冲击波形成距离,ρ0为声传播介质的密度,β为声传播介质的非线性系数。ρ0,c0,α,β为声传播介质的常数,对于本实施例使用的纯水,ρ0=1000kg/m3,c0=1495m/s,α=0.025Np/m/MHz2,β=3.5。
SBE模型的边界条件为:
P ( &sigma; = - &sigma; max , &theta; &tau; s ) = sin &tau; s &theta; &le; &theta; 0 0 &theta; > &theta; 0 ;
式中,σmax=d/b,d为球壳换能器7几何焦距; &theta; 0 = tan - 1 ( tan &alpha; 0 1 + 1 / &sigma; max 2 ) , &alpha; 0 = sin - 1 ( a d ) , a为球壳换能器7张口半径。
将步骤二得到的球壳换能器7有效参数代入边界条件,数值求解SBE模型。本实施例通过有限差分算法计算SBE模型,具体过程为:第一步、通过傅里叶级数展开,将时域声压P变换到频域其中Cn为n次谐波幅度;第二步、将傅里叶级数展开式代入SBE模型,可以得到关于Cn的抛物线型偏微分方程;第三步、对于各次谐波,利用隐式后向差分算法,将偏微分方程离散化,再进行求解得到Cn。由于声场的该计算方法已经较为成熟,其具体计算过程就不再赘述,可参见相关文献。计算焦点处二次谐波与基波的比值R21
R21=C2/C1
改变球壳换能器7表面声压幅值P0,计算不同球壳换能器7表面声压幅值P0对应的R21,得到R21与P0的曲线。高强度聚焦超声的传播非线性与球壳换能器7表面声压幅度、传播距离以及传播介质的声学参数有关,当确定了传播介质,在焦点的波形非线性畸变程度(用焦点处的二次谐波C2与基波C1衡量)仅仅与球壳换能器7表面声压幅度有关,且表面声压幅度越强,焦点波形的非线性畸变越明显,即R21越大,因此R21随P0增加而单调增加(参见图5)。
步骤四、对球壳换能器7施加任意幅度的单频电信号激励,激发声场,利用光纤水听器测量焦点波形p(t),对该焦点波形进行傅里叶变换:
S ( &omega; ) = &Integral; - &infin; &infin; p ( t ) e - j&omega;t dt
计算得R21′:
R21′=S(2ω0)/S(ω0)
利用焦点二次谐波与基波比值R21与球壳换能器7表面声压P0的单调关系,测定R21即可获得P0;根据步骤三得到的R21与P0的关系曲线进行插值,得到R21′对应的球壳换能器7表面声压P0′。
步骤五、将步骤四获得的表面声压P0′带入步骤三所述的SBE模型,计算得到步骤四施加的电信号激励对应的非线性声场分布。如图6所示,信号发生器3输出信号幅值与球壳换能器7表面声压幅值的关系曲线。依据图6,在任意激励下的声场分布均可由SBE模型计算得到。
如图7所示,在表面声压分别为81kPa、259kPa、350kPa时,测量和计算得到的焦点波形,以及各自相应频谱。在表面声压为81kPa和259kPa时,实验测量的波形和理论计算的结果吻合很好;在表面声压为350kPa时,冲击波产生,实验测量结果和计算结果存在较大差异。对焦点波形进行傅里叶变换,比较各次谐波。示波器采样频率为100MHz,分析50MHz内的各次谐波幅值。由图7可见,在表面声压350kPa时,50MHz内的各次谐波幅值吻合的很好,这说明波形的差异主要由于测量系统带宽不够。另外,图7中在表面声压81kPa时,频谱在高于10MHz后,存在明显差异;主要原因是此时高次谐波太小,淹没在系统噪声中。
如图8所示,在表面声压分别为81kPa,259kPa,350kPa时,焦平面径向声压分布。计算结果和实验结果吻合较好。
综上所述,本实施例提供的强聚焦超声非线性声场分布测量方法,通过引入焦点二次谐波与基波比值R21反映球壳换能器7边界条件,再根据非线性声传播模型,得到对应激励下的声场分布。由于仅需要测量二次谐波和基波,有效的解决了强非线性下,测量系统带宽不够的问题;此外由于R21为相对值,测量中不需要使用水听器6的灵敏度,因此避免了水听器6灵敏度不高的隐患;本实施例的测量方法操作简便,实用性强,使用的测量系统结构简单、便于控制,测量准确度高、容易实施。

Claims (7)

1.一种测量强聚焦超声非线性声场分布的方法,其测量步骤为:
步骤一、在线性声场条件下,测量球壳换能器轴向以及该球壳换能器焦平面径向的声压,获得球壳换能器轴向及其焦平面径向的声压分布曲线;
步骤二、根据瑞利积分计算球壳换能器轴向及该球壳换能器焦平面径向的声压分布曲线,将计算结果与步骤一获得的测量结果相比较,调节球壳换能器的张口半径a和几何焦距d,使得计算获得的球壳换能器轴向及其焦平面径向声压分布曲线与步骤一的测量结果相吻合,获得球壳换能器有效参数a和d;
步骤三、将步骤二获得的有效参数代入非线性声传播模型,计算焦点的二次谐波与基波的比值R21,再改变球壳换能器表面声压幅值P0,计算不同表面声压幅值P0对应的R21,得R21与P0的关系曲线;
步骤四、对球壳换能器施加单频电信号激励,测量焦点波形,对该焦点波形进行傅里叶变换,计算得R21′,根据步骤三得到的R21与P0的关系曲线进行插值,得到R21′对应的球壳换能器表面声压P0′;
步骤五、将步骤四获得的表面声压P0′带入步骤三所述的非线性声传播模型,计算得到步骤四施加的电信号激励对应的非线性声场分布。
2.根据权利要求1所述的一种测量强聚焦超声非线性声场分布的方法,其特征在于:步骤一测量声压的条件为:对球壳换能器施加低能量激励,使得球壳换能器焦点处基波幅度与二次谐波幅度的差值不低于40dB,忽略声传播过程中的非线性效应,将声场近似为线性声场。
3.根据权利要求2所述的一种测量强聚焦超声非线性声场分布的方法,其特征在于:步骤二涉及的瑞利积分公式为:
p ( M ) = - jk 2 &pi; &Integral; &Integral; Source p 0 ( N ) exp ( - jkr ) r dS
式中,p(M)为声场空间任意一点M的声压,p0(N)为球壳换能器表面任意一点N处声压;k=ω0/c0为波数,其中,ω0=2πf0为声波的角频率,f0为声波的频率,c0为声速;r为M、N两点间距离,Source表示球壳换能器表面区域。
4.根据权利要求3所述的一种测量强聚焦超声非线性声场分布的方法,其特征在于:步骤二要求计算获得的球壳换能器轴向及其焦平面径向声压分布曲线与步骤一测量结果的吻合程度为:测量结果与理论计算结果的主瓣偏差低于5%。
5.根据权利要求2或3所述的一种测量强聚焦超声非线性声场分布的方法,其特征在于:步骤三在椭球坐标系下计算声场分布,椭球坐标系中一点
Figure FDA0000468404380000012
与直角坐标系(x,y,z)的变换关系为:
Figure FDA0000468404380000021
z=bση
式中,-∞<σ<+∞,0≤η≤1,2b为椭球焦距;为简化声传播方程,作变换η=cosθ,0≤θ≤π/2;
所述的非线性声传播模型为SBE模型,即:
&PartialD; 2 P &PartialD; &tau; s &PartialD; &sigma; + 1 2 sin 2 &theta; &sigma; ( 1 + &sigma; 2 ) &PartialD; 2 P &PartialD; &tau; s &PartialD; &theta; + &epsiv; &sigma; 2 + sin 2 &theta; &sigma; ( 1 + &sigma; 2 ) ( &PartialD; 2 P &PartialD; &theta; 2 + cot &theta; &PartialD; P &PartialD; &theta; ) + E &sigma; &PartialD; P &PartialD; &tau; s = - &sigma; 2 + sin 2 &theta; &sigma; ( &alpha;b &PartialD; 3 P &PartialD; &tau; s 3 + b 2 l D &PartialD; 2 P 2 &PartialD; &tau; s 2 ) E , ( &sigma; < &sigma; 0 < 0 )
&PartialD; 2 P &PartialD; &tau; p &PartialD; &sigma; - &sigma; sin 2 &theta; 1 + &sigma; 2 &PartialD; 2 P &PartialD; &tau; p &PartialD; &theta; - &epsiv; ( 2 - cos &theta; ) 1 + &sigma; 2 ( &PartialD; 2 P &PartialD; &theta; 2 + cot &theta; &PartialD; P &PartialD; &theta; ) = ( &alpha;b &PartialD; 3 P &PartialD; &tau; p 3 + b 2 l D &PartialD; 2 P 2 &PartialD; &tau; p 2 ) E , ( &sigma; &GreaterEqual; &sigma; 0 )
其中,σ<σ0<0部分为声场靠近声源表面的球面波近似部分,σ≥σ0部分为声场远离声源表面的平面波近似部分,σ0为球面波区域与平面波区域分界面;式中,P=p/P0为归一化声压,p为点
Figure FDA0000468404380000027
处声压,P0为球壳换能器表面声压;τs=ω0[t+b(σ2+sin2θ)1/2/c0]为球面波弛豫时间,τp=ω0[t-bσcosθ/c0]为平面波弛豫时间,t为球壳换能器激发声波在传播介质中的传播时间;E=(σ2+cos2θ)/(1+σ2);α为声传播介质声衰减系数;lD=ρ0c0 3/(βP0ω0)为冲击波形成距离,ρ0为声传播介质密度,β为声传播介质非线性系数,ρ0,α,β均为常数。
6.根据权利要求5所述的一种测量强聚焦超声非线性声场分布的方法,其特征在于:步骤三所用的SBE模型的边界条件为:
P ( &sigma; = - &sigma; max , &theta; &tau; s ) = sin &tau; s &theta; &le; &theta; 0 0 &theta; > &theta; 0
式中,σmax=d/b, &theta; 0 = tan - 1 ( tan &alpha; 0 1 + 1 / &sigma; max 2 ) , &alpha; 0 = sin - 1 ( a d ) .
7.一种测量强聚焦超声非线性声场分布的系统,包括示波器(1)、计算机(2)、信号发生器(3)、功率放大器(4)、水听器(6)、球壳换能器(7)和水槽(8),其特征在于:还包括三维运动平台(5),所述的示波器(1)、计算机(2)、信号发生器(3)、功率放大器(4)和球壳换能器(7)依次相连,功率放大器(4)将信号发生器(3)的输出信号放大,驱动球壳换能器(7)激发声场;所述的球壳换能器(7)固定于水槽(8)内部,水槽(8)内盛装有水;所述的三维运动平台(5)控制水听器(6)移动,测量球壳换能器(7)轴向及该球壳换能器(7)焦平面径向声压分布;所述的示波器(1)采集水听器(6)的输出信号;所述的计算机(2)用于调控整个测量系统,计算机(2)完成的功能包括:控制信号发生器(3)输出,控制三维运动平台(5)运动以及存储数据。
CN201410059590.5A 2014-02-21 2014-02-21 一种测量强聚焦超声非线性声场分布的方法 Active CN103776524B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410059590.5A CN103776524B (zh) 2014-02-21 2014-02-21 一种测量强聚焦超声非线性声场分布的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410059590.5A CN103776524B (zh) 2014-02-21 2014-02-21 一种测量强聚焦超声非线性声场分布的方法

Publications (2)

Publication Number Publication Date
CN103776524A true CN103776524A (zh) 2014-05-07
CN103776524B CN103776524B (zh) 2016-02-10

Family

ID=50569043

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410059590.5A Active CN103776524B (zh) 2014-02-21 2014-02-21 一种测量强聚焦超声非线性声场分布的方法

Country Status (1)

Country Link
CN (1) CN103776524B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106813774A (zh) * 2017-01-10 2017-06-09 中国计量大学 一种高强度聚焦超声的声强和声功率测量方法
CN108175442A (zh) * 2018-01-15 2018-06-19 中国计量大学 一种超声造影剂非线性参数测量方法
CN108760899A (zh) * 2018-08-03 2018-11-06 广东工业大学 一种超声换能器模型制作方法和装置
CN109891230A (zh) * 2016-11-07 2019-06-14 球波株式会社 用于测量气体浓度的系统、方法和计算机程序产品
CN109931995A (zh) * 2019-02-01 2019-06-25 天津大学 一种获取气体超声流量计内部声场的方法
CN110243454A (zh) * 2019-07-08 2019-09-17 西安交通大学 粘弹媒质中基于双频叠加超声脉冲的微泡增效动力学实验系统及方法
CN110440906A (zh) * 2018-05-04 2019-11-12 重庆海扶医疗科技股份有限公司 超声换能器的声场声强分布检测方法及其装置
CN110954207A (zh) * 2019-11-15 2020-04-03 重庆医科大学 一种聚焦超声焦点声波结构的检测装置及检测方法
CN111342908A (zh) * 2020-02-18 2020-06-26 中国人民解放军空军工程大学 聚束声场处理方法、装置和电子设备
CN112033518A (zh) * 2020-08-18 2020-12-04 中国船舶重工集团公司第七一五研究所 一种基于正弦函数拟合的水听器灵敏度标定方法
CN112067112A (zh) * 2020-09-21 2020-12-11 中国石油大学(华东) 一种水下超声空化场特征自动测量与可视化系统、方法
CN112967709A (zh) * 2021-02-20 2021-06-15 河北工业大学 一种基于无损脑神经刺激聚焦声场的相控方法
CN113179476A (zh) * 2021-04-12 2021-07-27 科大讯飞(苏州)科技有限公司 配置参数的获取方法、配置方法、电子设备及存储装置
CN113884173A (zh) * 2021-09-27 2022-01-04 浙江省计量科学研究院 一种高强度聚焦超声换能器焦点峰值声压联合定征方法
CN117872273A (zh) * 2024-03-11 2024-04-12 厦门市盛迅信息技术股份有限公司 基于人工智能的多环境声场声线识别方法和系统
CN118225225A (zh) * 2024-05-24 2024-06-21 山东科技大学 一种基于声学超材料的水下低频声长期监测系统及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4843597A (en) * 1986-01-13 1989-06-27 Gjessing Dag K T System for marine seismic exploration
CN102141427A (zh) * 2010-12-02 2011-08-03 中国船舶重工集团公司第七一五研究所 一种利用激光测振仪检测液体媒介中声场参数的方法
CN102662007A (zh) * 2012-05-23 2012-09-12 北京理工大学 相控阵超声换能器声场扫查方法
CN102818850A (zh) * 2012-08-31 2012-12-12 中国船舶重工集团公司第七一五研究所 一种大面积水声材料透声性能角谱的测量方法
CN102944298A (zh) * 2012-11-06 2013-02-27 重庆大学 端部敏感型光纤光栅高强度聚焦超声传感器及系统
CN103234619A (zh) * 2013-04-25 2013-08-07 重庆大学 光纤法布里-珀罗超声水听器及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4843597A (en) * 1986-01-13 1989-06-27 Gjessing Dag K T System for marine seismic exploration
CN102141427A (zh) * 2010-12-02 2011-08-03 中国船舶重工集团公司第七一五研究所 一种利用激光测振仪检测液体媒介中声场参数的方法
CN102662007A (zh) * 2012-05-23 2012-09-12 北京理工大学 相控阵超声换能器声场扫查方法
CN102818850A (zh) * 2012-08-31 2012-12-12 中国船舶重工集团公司第七一五研究所 一种大面积水声材料透声性能角谱的测量方法
CN102944298A (zh) * 2012-11-06 2013-02-27 重庆大学 端部敏感型光纤光栅高强度聚焦超声传感器及系统
CN103234619A (zh) * 2013-04-25 2013-08-07 重庆大学 光纤法布里-珀罗超声水听器及系统

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109891230B (zh) * 2016-11-07 2021-05-14 球波株式会社 用于测量气体浓度的系统、方法和计算机程序产品
CN109891230A (zh) * 2016-11-07 2019-06-14 球波株式会社 用于测量气体浓度的系统、方法和计算机程序产品
CN106813774A (zh) * 2017-01-10 2017-06-09 中国计量大学 一种高强度聚焦超声的声强和声功率测量方法
CN108175442A (zh) * 2018-01-15 2018-06-19 中国计量大学 一种超声造影剂非线性参数测量方法
CN110440906A (zh) * 2018-05-04 2019-11-12 重庆海扶医疗科技股份有限公司 超声换能器的声场声强分布检测方法及其装置
CN108760899A (zh) * 2018-08-03 2018-11-06 广东工业大学 一种超声换能器模型制作方法和装置
CN108760899B (zh) * 2018-08-03 2021-01-29 广东工业大学 一种超声换能器模型制作方法和装置
CN109931995A (zh) * 2019-02-01 2019-06-25 天津大学 一种获取气体超声流量计内部声场的方法
CN109931995B (zh) * 2019-02-01 2020-07-14 天津大学 一种获取气体超声流量计内部声场的方法
CN110243454A (zh) * 2019-07-08 2019-09-17 西安交通大学 粘弹媒质中基于双频叠加超声脉冲的微泡增效动力学实验系统及方法
CN110954207A (zh) * 2019-11-15 2020-04-03 重庆医科大学 一种聚焦超声焦点声波结构的检测装置及检测方法
CN111342908A (zh) * 2020-02-18 2020-06-26 中国人民解放军空军工程大学 聚束声场处理方法、装置和电子设备
CN111342908B (zh) * 2020-02-18 2021-09-10 中国人民解放军空军工程大学 聚束声场处理方法、装置和电子设备
CN112033518A (zh) * 2020-08-18 2020-12-04 中国船舶重工集团公司第七一五研究所 一种基于正弦函数拟合的水听器灵敏度标定方法
CN112067112A (zh) * 2020-09-21 2020-12-11 中国石油大学(华东) 一种水下超声空化场特征自动测量与可视化系统、方法
CN112967709A (zh) * 2021-02-20 2021-06-15 河北工业大学 一种基于无损脑神经刺激聚焦声场的相控方法
CN113179476A (zh) * 2021-04-12 2021-07-27 科大讯飞(苏州)科技有限公司 配置参数的获取方法、配置方法、电子设备及存储装置
CN113884173A (zh) * 2021-09-27 2022-01-04 浙江省计量科学研究院 一种高强度聚焦超声换能器焦点峰值声压联合定征方法
CN117872273A (zh) * 2024-03-11 2024-04-12 厦门市盛迅信息技术股份有限公司 基于人工智能的多环境声场声线识别方法和系统
CN117872273B (zh) * 2024-03-11 2024-05-31 厦门市盛迅信息技术股份有限公司 基于人工智能的多环境声场声线识别方法和系统
CN118225225A (zh) * 2024-05-24 2024-06-21 山东科技大学 一种基于声学超材料的水下低频声长期监测系统及方法

Also Published As

Publication number Publication date
CN103776524B (zh) 2016-02-10

Similar Documents

Publication Publication Date Title
CN103776524B (zh) 一种测量强聚焦超声非线性声场分布的方法
Bessonova et al. Focusing of high power ultrasound beams and limiting values of shock wave parameters
CN109085247B (zh) 一种超声造影剂空化泡群共振状态测量方法
CN109745077B (zh) 基于聚焦超声声振信号的弹性特性检测方法
Jin et al. An online impedance analysis and matching system for ultrasonic transducers
Baker Nonlinear effects in ultrasound propagation
CN110243950B (zh) 基于幅值调制的声辐射力脉冲激励方法
CN103143126A (zh) 生物组织非线性hifu声场确定的方法
CN103575377B (zh) 测量参量声场中差频波空间分布特性的方法
Filipczyński et al. Nonlinear and linear propagation of diagnostic ultrasound pulses
CN117874703A (zh) 基于超声激发组织位移场的焦域声吸收能量反演方法及系统
CN106768264B (zh) 基于焦域轴向最大振速的聚焦超声功率测量系统及方法
Chen et al. Acoustic characterization of high intensity focused ultrasound fields generated from a transmitter with a large aperture
Chen et al. Modeling of shock wave generated from a strong focused ultrasound transducer
Li et al. Determining the responsivity of air-coupled piezoelectric transducers using a comparative method: Theory and experiments
Fan et al. Acoustic characterization of high intensity focused ultrasound field generated from a transmitter with large aperture
Tao et al. Accurate acoustic power measurement for low-intensity focused ultrasound using focal axial vibration velocity
CN109044485B (zh) 一种高声压高声强超声波测量系统
Shou et al. Calibration of a focusing transducer and miniature hydrophone as well as acoustic power measurement based on free-field reciprocity in a spherically focused wave field
Li et al. Nonlinear acoustic-power measurement based on fundamental focal axial vibration velocity for high-intensity focused ultrasound
Karzova et al. Nonlinear effects in ultrasound fields of diagnostic-type transducers used for kidney stone propulsion: Characterization in water
Wójcik et al. Fast prediction of pulsed nonlinear acoustic fields from clinically relevant sources using time-averaged wave envelope approach: Comparison of numerical simulations and experimental results
Ming-He et al. Nonlinear effect on focusing gain of a focusing transducer with a wide aperture angle
Fan et al. A modeling approach to predict acoustic nonlinear field generated by a transmitter with an aluminum lens
Liu et al. The transmission of finite amplitude sound beam in multi-layered biological media

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant