CN102564571B - 基于红外测温的高强度聚焦超声声场测量方法 - Google Patents
基于红外测温的高强度聚焦超声声场测量方法 Download PDFInfo
- Publication number
- CN102564571B CN102564571B CN2012100303838A CN201210030383A CN102564571B CN 102564571 B CN102564571 B CN 102564571B CN 2012100303838 A CN2012100303838 A CN 2012100303838A CN 201210030383 A CN201210030383 A CN 201210030383A CN 102564571 B CN102564571 B CN 102564571B
- Authority
- CN
- China
- Prior art keywords
- sound
- ultrasonic
- absorption thing
- temperature rise
- ultrasonic absorption
- 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.)
- Active
Links
Images
Landscapes
- Radiation Pyrometers (AREA)
Abstract
本发明涉及一种基于红外测温的高强度聚焦超声声场测量方法,包括:第一步,水听器测量聚焦超声声源的声强值;第二步,放置超声吸收物,水听器再次测量声强值,将第二步中获得的声强值与第一步获得的声强值进行比较,得到的衰减物理量值;第三步,通过红外测温仪测量超声吸收物在聚焦超声声源设定时间照射后在表面形成的温升,通过计算求取表面温升的最大温升率,得到超声吸收物的声热参数综合特征值;第四步,通过使用红外测温仪测量在不同位置下超声吸收物表面的温升求取对应的最大温升率,将此最大温升率与衰减物理量值和声热参数综合特征值进行运算得到相对应不同位置的声强分布。与现有技术相比,本发明具有实现对HIFU声场三维快速测量等优点。
Description
技术领域
本发明涉及一种高强度聚焦超声声场测量方法,尤其是涉及一种基于红外测温的高强度聚焦超声声场测量方法。
背景技术
高强度聚焦超声三维声场分布是衡量HIFU治疗探头设计与制造工艺、治疗效果好坏的关键因素,也是衡量设备安全高效的重要指标,此外,HIFU三维声场分布也对治疗计划的制定有着重要的参考意义,所以高强度聚焦超声声场的测量是HIFU研究领域中的一个重要环节。
自上世纪90年代以来,国际电工学会(IEC)在超声声场参数测量方面制定了大量的标准,但是由于聚焦超声声场具有如下特点:
焦点体积通常很小(呈米粒状,横截面直径约1~3mm,声束轴长度约8~15mm);
焦域的最大声强≥10000Wcm-2;
声压超过10Mpa;
非线性谐波、声饱和。
鉴于高强度聚焦超声声场的复杂性,在2005年之前,HIFU声场检测标准的制定都在讨论中。2005年9月,我国国家标准化管理委员会批准发布了《高强度聚焦超声(HIFU)声功率和声场特性的测量》国家标准(GB/T 19890),并从2006年4月1日起实施,成为世界上首个HIFU测量标准。此外,2005年12月7日,国家武汉医用超声设备检测中心起草的《高强度聚焦超声(HIFU)治疗系统》行业标准(YY0592-2005)由国家药监局批准发布,自2006年12月1日起实施。
目前已有的关于聚焦声场特性和超声功率测量标准的基本思想和做法是:在小功率条件下(全功率输出的10%~20%),用水听器三维扫描法找到声压焦点位置,对聚焦声场进行扫描测量,并计算出焦点的几何参数;用水听器或辐射力法测量超声功率,再按照相关公式求出最大声强。
这些标准和方法的主要存在以下问题:
1、获取聚焦声场三维特征值的过程十分耗时;
2、焦域的声强能量过大,容易导致水听器的损坏;
3、换能器辐射面各点振动相位若不全一致,则焦点处声强不一定正比于声压平方;
4、HIFU设备在全功率输出时,超声在水中或组织中传播时会产生非线性和衍射效应,而非线性累积在水中和组织中是不同的,所以不能用水中测量得到的声场分布直接反映组织中的情况。
5、由于有限振幅波在传播过程中会产生波形畸变,滋生谐波、声衰减、声饱和等非线性现象,通过外插办法来推测HIFU声场的方法得到的结果并不理想。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于红外测温的高强度聚焦超声声场测量方法。
本发明的目的可以通过以下技术方案来实现:
一种基于红外测温的高强度聚焦超声声场测量方法,其特征在于,通过使用红外测温法测量置于超声声场内的超声吸收物表面温度分布,实现测量超声声场分布,具体包括以下步骤:
第一步,通过校准过的水听器测量聚焦超声声源在焦点位置的声强值;
第二步,在水听器和超声声源之间放置超声吸收物,水听器再次测量聚焦超声声源在经过超声吸收物后的声强值,将第二步中获得的声强值与第一步获得的声强值进行比较,得到超声吸收物对超声声源的衰减物理量值;
第三步,将超声吸收物的上表面平移至焦平面位置,通过红外测温仪测量超声吸收物在聚焦超声声源设定时间照射后在表面形成的温升,通过计算求取表面温升的最大温升率,并将这个温升率与第二步中获得的声强值进行运算得到超声吸收物的声热参数综合特征值;
第四步,将超声吸收物在超声声场范围内沿着声束轴的方向平行移动,通过使用红外测温仪测量在不同位置下超声吸收物表面的温升求取对应的最大温升率,将此最大温升率与第二步得到的衰减物理量值和第三步得到的声热参数综合特征值进行运算得到相对应不同位置的声强分布。
所述的超声吸收物与水听器之间距离为1~5mm。
所述的超声吸收物对超声声源的衰减物理量值计算如下:FAmplify=1/FAttenuation=I0/I1,其中FAttenuation为超声吸收物对超声声源的衰减物理量值,FAmplify为FAttenuation的倒数,I0和I1分别为第一步和第二步测量到的焦点位置声强值。
所述的最大温升率的计算如下:TCR0max=max(ΔT0/Δt),TCR1max=max(ΔT1/Δt),其中,TCR0max是低驱动电压下得到的超声吸收物表面的最大温升率,用于求取超声吸收物声热参数综合特征值;TCR1max是高驱动电压下得到的超声吸收物表面的最大温升率,用于求取超声声场中对应水平面的超声声强值,Δt为间隔时间;ΔT0为在低驱动电压下超声吸收物表面在间隔时间内的温度变化;ΔT1为高驱动电压下超声吸收物表面在间隔时间内的温度变化。
所述的声热参数综合特征值的计算如下:MTransform=I1/TCR0max,其中MTransform为表征超声吸收物的声热参数综合特征值,I1为测量第二步测量到的焦点位置声强值,TCR0max是低驱动电压下得到的超声吸收物表面的最大温升率。
所述的将此最大温升率与第二步得到的衰减物理量值和第三步得到的声热参数综合特征值进行运算具体为:IEstimate=FAmplify×MTransform×TCR1max,其中IEstimate为超声声场声强值;FAmplify为表征超声吸收物对超声声源的衰减物理量;MTransform为表征超声吸收物的声热参数综合特征值;TCR1max是高驱动电压下得到的超声吸收物表面的最大温升率。
与现有技术相比,本发明具有根据超声声场在被测物中的声场三维分布和利用能量热传导方程求解出的三维温度场分布,通过使用红外测温技术快速测量被测物表面温度,以实现对HIFU声场三维快速可视化测量。
附图说明
图1为本发明的流程示意图;
图2为本发明实现的硬件结构示意图;
图3为本发明构建的软件系统功能模块图。
其中图1、图2,1为水听器,2为红外测温仪,3为超声吸收物,4为超声探头。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图1、图2所示,一种基于红外测温的高强度聚焦超声声场测量方法,通过使用红外测温法测量置于超声声场内的超声吸收物表面温度分布,实现测量超声声场分布,具体包括以下步骤:
第一步,通过校准过的水听器1测量聚焦超声声源在焦点位置的声强值;
第二步,在水听器1和超声声源之间放置超声吸收物3,所述的超声吸收物3与水听器1之间距离为1~5mm,水听器1再次测量聚焦超声声源在经过超声吸收物3后的声强值,将第二步中获得的声强值与第一步获得的声强值进行比较,得到超声吸收物3对超声声源的衰减物理量值;
第三步,将超声吸收物3的上表面平移至焦平面位置,通过红外测温仪2测量超声吸收物3在聚焦超声声源设定时间照射后在表面形成的温升,通过计算求取表面温升的最大温升率,并将这个温升率与第二步中获得的声强值进行运算得到超声吸收物3的声热参数综合特征值;
第四步,将超声吸收物3在超声声场范围内沿着声束轴的方向平行移动,通过使用红外测温仪2测量在不同位置下超声吸收物3表面的温升求取对应的最大温升率,将此最大温升率与第二步得到的衰减物理量值和第三步得到的声热参数综合特征值进行运算得到相对应不同位置的声强分布。
Claims (2)
1.一种基于红外测温的高强度聚焦超声声场测量方法,其特征在于,通过使用红外测温法测量置于超声声场内的超声吸收物表面温度分布,实现测量超声声场分布,具体包括以下步骤:
第一步,通过校准过的水听器测量聚焦超声声源在焦点位置的声强值;
第二步,在水听器和超声声源之间放置超声吸收物,水听器再次测量聚焦超声声源在经过超声吸收物后的声强值,将第二步中获得的声强值与第一步获得的声强值进行比较,得到超声吸收物对超声声源的衰减物理量值;
第三步,将超声吸收物的上表面平移至焦平面位置,通过红外测温仪测量超声吸收物在聚焦超声声源设定时间照射后在表面形成的温升,通过计算求取表面温升的最大温升率,并将这个温升率与第二步中获得的声强值进行运算得到超声吸收物的声热参数综合特征值;
第四步,将超声吸收物在超声声场范围内沿着声束轴的方向平行移动,通过使用红外测温仪测量在不同位置下超声吸收物表面的温升求取对应的最大温升率,将此最大温升率与第二步得到的衰减物理量值和第三步得到的声热参数综合特征值进行运算得到相对应不同位置的声强值;
所述的超声吸收物对超声声源的衰减物理量值计算如下:FAmplify=1/FAttenuation=I0/I1,其中FAttenuation为超声吸收物对超声声源的衰减物理量值,FAmplify为FAttenuation的倒数,I0和I1分别为第一步和第二步测量到的焦点位置声强值;
所述的最大温升率的计算如下:TCR0max=max(ΔT0/Δt),TCR1max=max(ΔT1/Δt),其中,TCR0max是低驱动电压下得到的超声吸收物表面的最大温升率,用于求取超声吸收物声热参数综合特征值;TCR1max是高驱动电压下得到的超声吸收物表面的最大温升率,用于求取超声声场中对应水平面的超声声强值;Δt为间隔时间;ΔT0为在低驱动电压下超声吸收物表面在间隔时间内的温度变化;ΔT1为高驱动电压下超声吸收物表面在间隔时间内的温度变化;
所述的声热参数综合特征值的计算如下:MTransform=I1/TCR0max,其中MTransform为表征超声吸收物的声热参数综合特征值,I1为测量第二步测量到的焦点位置声强值,TCR0max是低驱动电压下得到的超声吸收物表面的最大温升率;
所述的将此最大温升率与第二步得到的衰减物理量值和第三步得到的声热参数综合特征值进行运算具体为:IEstimate=FAmplify×MTransform×TCR1max,其中IEstimate为超声声场声强值;FAmplify为表征超声吸收物对超声声源的衰减物理量;MTransform为表征超声吸收物的声热参数综合特征值;TCR1max是高驱动电压下得到的超声吸收物表面的最大温升率。
2.根据权利要求1所述的一种基于红外测温的高强度聚焦超声声场测量方法,其特征在于,所述的超声吸收物与水听器之间距离为1~5mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100303838A CN102564571B (zh) | 2012-02-10 | 2012-02-10 | 基于红外测温的高强度聚焦超声声场测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100303838A CN102564571B (zh) | 2012-02-10 | 2012-02-10 | 基于红外测温的高强度聚焦超声声场测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102564571A CN102564571A (zh) | 2012-07-11 |
CN102564571B true CN102564571B (zh) | 2013-11-13 |
Family
ID=46410600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012100303838A Active CN102564571B (zh) | 2012-02-10 | 2012-02-10 | 基于红外测温的高强度聚焦超声声场测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102564571B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106388861B (zh) * | 2014-10-21 | 2019-05-03 | 孙艳华 | 自动化程度高的便携移动式医用超声检测装置的使用方法 |
CN107874779B (zh) * | 2017-10-31 | 2020-12-18 | 江西中医药大学 | 一种高强度聚焦超声测量装置及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7611278B2 (en) * | 2003-12-02 | 2009-11-03 | White Box, Inc. | Infrared thermometers |
CN100365432C (zh) * | 2005-10-11 | 2008-01-30 | 汪柏年 | 一种基于超声波或声波连续声场鉴相原理的位移/距离测量方法及其装置 |
CN101382451B (zh) * | 2008-11-05 | 2010-09-22 | 北京理工大学 | 超声液浸换能器声场精密测量方法 |
-
2012
- 2012-02-10 CN CN2012100303838A patent/CN102564571B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102564571A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144356B (zh) | 非制冷红外焦平面探测器阵列热响应时间测试系统及方法 | |
Fu et al. | Impact source identification in a carbon fiber reinforced polymer plate by using embedded fiber optic acoustic emission sensors | |
CN102857850B (zh) | 一种高频超声发射器及阵列声学参数近场校准方法 | |
Wang et al. | Modelling nonlinear ultrasound propagation in absorbing media using the k-Wave toolbox: experimental validation | |
CN105806466B (zh) | 一种超声波功率测量系统 | |
CN106813774A (zh) | 一种高强度聚焦超声的声强和声功率测量方法 | |
CN101825497A (zh) | 一种基于热声效应的温度实时测量系统与方法 | |
CN106422090B (zh) | 一种医用超声设备声场快速校准方法 | |
CN104501939A (zh) | 一种利用单水听器测量非消声水池混响时间的反演方法 | |
CN108318123A (zh) | 一种水下声基阵互辐射阻抗的近场声层析测试方法 | |
CN102564571B (zh) | 基于红外测温的高强度聚焦超声声场测量方法 | |
CN107132279A (zh) | 一种基于阵列超声表面波的构件损伤概率成像定位方法 | |
US20130265856A1 (en) | Thermoacoustic sensor | |
CN107874779B (zh) | 一种高强度聚焦超声测量装置及方法 | |
CN105806471A (zh) | 一种超声空化噪声信号测量装置与方法 | |
Civale et al. | Focused ultrasound transducer spatial peak intensity estimation: a comparison of methods | |
JP2015230171A (ja) | 避雷素子及びそれを覆う碍管からなる避雷器における避雷素子の温度測定方法 | |
CN105136069A (zh) | 一种激光超声波技术测量圆柱体内锥形锥尖角度的方法 | |
CN103575381B (zh) | 基于动态光弹性法的超声换能器声场的测量方法 | |
CN105300887A (zh) | 一种激光超声无损检测材料织构的方法 | |
Uchida et al. | Development of high ultrasonic power measurement technique by calorimetric method using water as heating elements for high ultrasonic power standard | |
CN109000780B (zh) | 一种升降温式声速测量仪 | |
Koukoulas et al. | Uncertainty contributions in the optical measurement of free-field propagating sound waves in air and water | |
CN103207321B (zh) | 基于热声效应的脉冲微波辐射场分布的测量系统与方法 | |
CN103983344A (zh) | 一种超声光栅位相振幅的定量测量方法 |
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 |