-
超声是检测不透明液体中气泡的有效方法, 声散射模型是超声反演技术的核心. 经典气泡散射模型通常是基于球形假设及ka ≤ 1 (a为气泡半径, k为入射波的波数), 然而实际应用中这些假设并不总能得到满足. 本研究针对非球形气泡及ka偏离假设情况, 提出一种超声反演定量方案. 建立不受ka约束的球形气泡级数背散射模型, 将其与经典Medwin (
$ka\ll1$ )和Anderson (ka≈1)散射模型进行对比, 发现ka偏离引发的散射截面效应仅体现在散射高阶共振峰位置及大小差异上. 据此提出: 可通过曲线拟合法解决散射截面σbs/(πa2)与ka间的多值问题, 同时用当量半径a* 对非球形气泡尺寸进行量化. 具体首先利用非球形气泡背散射信号的频域信息测定其散射截面σbs, 再根据σbs与非球形当量半径a* 间拟合曲线进行反演, 同时利用回波时域信息勾画气泡形状轮廓. 反演结果通过高速摄影定量结果进行检验. 结果表明: 气泡沿之字形路径上升过程中产生非球形形变, 当9 ≤ kr0 ≤ 35时, 反演得到的非球形气泡当量半径与高速摄影定量半径r0的相对误差小于45%, 对应的最大绝对误差约为1 mm. 说明该方法在一定测量精度范围可用于非球形气泡的声反演.Ultrasonic detection is an effective method to quantify bubbles in opaque liquid, and acoustic scattering model is the key in ultrasonic inversion technique. Classical scattering models are usually based on the spherical assumption, and ka is much less than 1. However, these conditions are not always satisfied in practical applications. In this study, a quantitative strategy of ultrasonic inversion is proposed for non-spherical bubbles and ka deviation assumption. A series of solution models for a spherical gas bubble is established without considering the ka constraint, and it is compared with the classical Medwin ($ka\ll1 $ ) and Anderson (ka ≈ 1) models. The difference in scattering cross section σbs betweem them is only at the higher order formants of scattering, so the fitted line can be used to solve the multi-valued problem between σbs and ka. For a non-spherical bubble, σbs is determined by the frequency domain backscattering signal, the size is characterized by the equivalent radius a*, and the inversion is performed by fitted curve from series solution model. Ultrasonic quantitative results are examined by high-speed photography. Results show that during the bubbles rising along a zigzag path, they develop non-spherical bubbles, their scattering cross sections are measured by the frequency domain scattering signal obtained at a position of ultrasonic measurement, and the equivalent radius is inverted by the series solution fitting curve. The deviation of the result from the actual result r0 is about 1mm (relative error less than 45%) when 9≤kr0≤35. This method can be used for implementing the acoustic inversion of non-spherical bubbles in a certain range of measurement accuracy.-
Keywords:
- ultrasonic detection /
- scattering cross section /
- non-spherical bubbles /
- equivalent radius
[1] Leighton T G, Dogan H, Fox P, Mantouka A, Best A I, Robb G B R, White P R 2021 J. Acoust. Soc. Am. 150 2705
Google Scholar
[2] Judd A G, Hovland M, Dimitrov L I, García Gil S, Jukes V 2002 Geofluids 2 109
Google Scholar
[3] Kracht W, Moraga C 2016 Miner. Eng. 98 122
Google Scholar
[4] Liu J, Gao Q, Tang Z, Xie Y, Gui W, Ma T, Niyoyita J P 2020 IEEE Trans. Instrum. Meas. 69 9618
Google Scholar
[5] Guédra M, Inserra C, Gilles B 2017 Ultrason. Sonochem. 38 298
Google Scholar
[6] Buckey J C, Knaus D A, Alvarenga D L, Kenton M A, Magari P J 2005 Acta Astronaut. 56 1041
Google Scholar
[7] Wen W, Zong G H, Bi S S 2014 Rev. Sci. Instrum. 85 065106
Google Scholar
[8] Jarmo I, Tuomas E, Heikki M, Lasse L, Jari K, Heikki K 2014 19th Iberoamerican Congress Puerto Vallarta, Mexico Puerto Vallarta, Mexico, November 2−5, 2014 p38
[9] Ilonen J, Juránek R, Eerola T, Lensu L, Dubská M, Zemčík P, Kälviäinen H 2018 Pattern Recognit. Lett. 101 60
Google Scholar
[10] Ohta J, Doyama N, Wakabayashi D, Suzuki H 2013 T. Jpn. Soc. Mech. Eng. Part B 79 2397
Google Scholar
[11] Bradley P B, Seth J P 1992 Phys. Rev. Lett. 69 03839
Google Scholar
[12] Ren W, Jin N, Zhang J 2022 Ultrasonics 124 106740
Google Scholar
[13] Padilla A M, Loranger S, Kinnaman F S, Valentine D L, Weber T C 2019 J. Geophys. Res. Ocean. 124 2472
Google Scholar
[14] Padilla A M, Weber T C 2021 J. Acoust. Soc. Am. 149 2504
Google Scholar
[15] Clay C S, Medwin H 1977 Acoustical Oceanography: Principles and Applications (New York: Wiley) pp461–466
[16] Anderson V C 1950 J. Acoust. Soc. Am. 22 426
Google Scholar
[17] Zheng B L, Poojitha D Y M 2000 J. Hydraul. Eng. 126 852
Google Scholar
[18] Dong X, Su M, Cai X 2012 Particuology 10 117
Google Scholar
[19] Spiekhout S, Voorneveld J, van Elburg B, Renaud G, Segers T, Lajoinie G P R, Versluis M, Verweij M D, de Jong N, Bosch J G 2022 J. Acoust. Soc. Am. 151 3993
Google Scholar
[20] Sage K A, George J, Oberall H 1979 J. Acoust. Soc. Am. 65 1413
Google Scholar
[21] Sam A, Gomez C O, Finch J A 1996 Int. J. Miner. Process. 47 177
Google Scholar
[22] 郑晖, 林树青 2008 超声检测 (北京: 中国劳动社会保障社) 第75页
Zheng H, Lin S Q 2008 Ultrasonic Inspection (Beijing: China Labor and Social Security Publishing House) p75 (in Chinese)
[23] Weber T C, Ward L G 2015 J. Acoust. Soc. Am. 138 2169
Google Scholar
-
图 1 气泡测量装置(1-高速摄影机, 2-水槽, 3-流量控制阀, 4-进气口, 5-超声换能器, 6-声卡, 7-工控机)
Fig. 1. Bubble measuring device. 1-high-speed camera, 2-tank, 3-flow control valves, 4-air inlet, 5-ultrasonic transducer, 6-sound card, 7-controller.
图 3 气泡散射模型对比 (a) 0.001≤ka≤1; (b) 1<ka≤50
Fig. 3. Comparison of scattering models for bubbles: (a) 0.001≤ka≤1; (b) 1<ka≤50.
图 4 上升气泡的形态变化 (h = 120 mm) (a) r0 = 0.5 mm; (b) r0 = 1.3 mm; (c) r0 = 2.5 mm
Fig. 4. Shape change of rising bubbles (h = 120 mm): (a) r0 = 0.5 mm; (b) r0 = 1.3 mm; (c) r0 = 2.5 mm.
图 5 气泡“之”字形上升轨迹 (a) r0 = 4.2 mm, h = 120, 160 mm; (b) r0 = 4.5 mm, h = 120, 160 mm; (c) r0 = 4.2 mm, h = 160 mm
Fig. 5. Zigzag rising trajectory of the bubbles: (a) r0 = 4.2 mm, h = 120, 160 mm; (b) r0 = 4.5 mm, h = 120, 160 mm; (c) r0 = 4.2 mm, h = 160 mm.
图 6 自同一水深处上升气泡轨迹 (h = 120 mm) (a) r0 = 0.5 mm; (b) r0 = 1.3 mm; (c) r0 = 2.5 mm
Fig. 6. Rising bubble trajectory with same depth (h = 120 mm): (a) r0 = 0.5 mm; (b) r0 = 1.3 mm; (c) r0 = 2.5 mm.
图 7 气泡散射时域信号及采集方法 (a)信号采集示意图; (b)不同角度的气泡时域散射信号
Fig. 7. Time domain scattering signal of bubble and data acquisition method: (a) Schematic diagram of signal acquisition; (b) time domain scattering signals of bubbles at different angles.
图 8 频域散射信号(a)和散射截面校准因子(b)
Fig. 8. Frequency domain scattering signal (a) and calibration factor of scattering cross section (b).
图 11 超声与高速摄影定量结果对比
Fig. 11. Comparison of quantitative results of ultrasound and high-speed camera.
表 1 测量位置处各泡的等效半径
$ r^* $ Table 1. Equivalent radius r* of each bubble at the measurement position.
r0/mm 0.5 1.0 1.3 2.5 3.2 4.2 4.5 4.8 5.5 r*/mm 0.5 1.1 1.4 2.4 3.2 4.6 4.4 4.9 5.9 χ 1.5 2.7 3.2 2.7 1.9 1.5 3.6 1.6 2.2 
下载: 导出CSV
表 2 测量位置处气泡的形状及形变率
Table 2. Shape and deformation rate of bubbles at the measurement position.
r0/mm 0.5 1.3 2.5 水平(χ*) 
竖直(χ)
下载: 导出CSV
-
[1] Leighton T G, Dogan H, Fox P, Mantouka A, Best A I, Robb G B R, White P R 2021 J. Acoust. Soc. Am. 150 2705
Google Scholar
[2] Judd A G, Hovland M, Dimitrov L I, García Gil S, Jukes V 2002 Geofluids 2 109
Google Scholar
[3] Kracht W, Moraga C 2016 Miner. Eng. 98 122
Google Scholar
[4] Liu J, Gao Q, Tang Z, Xie Y, Gui W, Ma T, Niyoyita J P 2020 IEEE Trans. Instrum. Meas. 69 9618
Google Scholar
[5] Guédra M, Inserra C, Gilles B 2017 Ultrason. Sonochem. 38 298
Google Scholar
[6] Buckey J C, Knaus D A, Alvarenga D L, Kenton M A, Magari P J 2005 Acta Astronaut. 56 1041
Google Scholar
[7] Wen W, Zong G H, Bi S S 2014 Rev. Sci. Instrum. 85 065106
Google Scholar
[8] Jarmo I, Tuomas E, Heikki M, Lasse L, Jari K, Heikki K 2014 19th Iberoamerican Congress Puerto Vallarta, Mexico Puerto Vallarta, Mexico, November 2−5, 2014 p38
[9] Ilonen J, Juránek R, Eerola T, Lensu L, Dubská M, Zemčík P, Kälviäinen H 2018 Pattern Recognit. Lett. 101 60
Google Scholar
[10] Ohta J, Doyama N, Wakabayashi D, Suzuki H 2013 T. Jpn. Soc. Mech. Eng. Part B 79 2397
Google Scholar
[11] Bradley P B, Seth J P 1992 Phys. Rev. Lett. 69 03839
Google Scholar
[12] Ren W, Jin N, Zhang J 2022 Ultrasonics 124 106740
Google Scholar
[13] Padilla A M, Loranger S, Kinnaman F S, Valentine D L, Weber T C 2019 J. Geophys. Res. Ocean. 124 2472
Google Scholar
[14] Padilla A M, Weber T C 2021 J. Acoust. Soc. Am. 149 2504
Google Scholar
[15] Clay C S, Medwin H 1977 Acoustical Oceanography: Principles and Applications (New York: Wiley) pp461–466
[16] Anderson V C 1950 J. Acoust. Soc. Am. 22 426
Google Scholar
[17] Zheng B L, Poojitha D Y M 2000 J. Hydraul. Eng. 126 852
Google Scholar
[18] Dong X, Su M, Cai X 2012 Particuology 10 117
Google Scholar
[19] Spiekhout S, Voorneveld J, van Elburg B, Renaud G, Segers T, Lajoinie G P R, Versluis M, Verweij M D, de Jong N, Bosch J G 2022 J. Acoust. Soc. Am. 151 3993
Google Scholar
[20] Sage K A, George J, Oberall H 1979 J. Acoust. Soc. Am. 65 1413
Google Scholar
[21] Sam A, Gomez C O, Finch J A 1996 Int. J. Miner. Process. 47 177
Google Scholar
[22] 郑晖, 林树青 2008 超声检测 (北京: 中国劳动社会保障社) 第75页
Zheng H, Lin S Q 2008 Ultrasonic Inspection (Beijing: China Labor and Social Security Publishing House) p75 (in Chinese)
[23] Weber T C, Ward L G 2015 J. Acoust. Soc. Am. 138 2169
Google Scholar
下载:
计量
- 文章访问数: 3593
- PDF下载量: 87
- 被引次数: 0
