搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用随机多项式展开的海底声学参数反演方法

李风华 王翰卓

引用本文:
Citation:

利用随机多项式展开的海底声学参数反演方法

李风华, 王翰卓

Geo-acoustic inversion using polynomial chaos expansion

Li Feng-Hua, Wang Han-Zhuo
PDF
HTML
导出引用
  • 为了提高地声反演算法的计算效率, 探索克服地声反演结果多值性问题, 本文利用宽带、多收发位置的传播损失数据结合传播损失在地声参数先验搜索区间内的随机多项式展开系数矩阵, 反演得到海底纵波声速、吸收率和密度比重. 使用随机多项式展开近似传播损失时, 展开系数的自变量为声波频率、收发位置等参数, 随机多项式的自变量为表示声速、吸收率、比重在各自搜索区间内均匀分布的随机变量. 传播损失的展开系数通过嵌入随机多项式的声学宽角抛物方程结合盖辽金投影、最小角度回归算法计算求得. 在低频、一定声传播水平距离以内和地声参数搜索区间长度适中时, 使用随机多项式展开近似传播损失的相对误差在1%以下. 仿真发现, 在浅海环境中使用低频、一定声传播水平距离以内的传播损失数据, 在接收信号信噪比较高、声源和水听器相对位置误差较小时, 选择合适的随机多项式展开截断幂次可较准确地反演海底声速、吸收率和密度比重, 且计算效率比网格遍历搜索方法提高一个数量级以上.
    In order to improve the computational efficiency of algorithms while exploring the method to overcome the ambiguity problems in underwater geo-acoustic inversion, we use the data of transmission losses at the broadband sound frequencies and multiple propagating distances with the matrix of polynomial chaos expansion coefficients of transmission losses to invert the speed (c), attenuation (α) of compression sound wave and the density ratio of seabed to seawater (ρ) in their prior searching intervals. When approximating the transmission loss with the polynomial chaos expansion, the expansion coefficients are the functions of parameters including sound frequency, source and hydrophone’s position while the polynomial bases are functions of the above geo-acoustic parameters which are uniformly distributed in their respective intervals. The expansion coefficients are calculated by embedding the orthogonal polynomial bases into the acoustic wide-angle parabolic equation model. After that, the coefficients are deduced using the Galerkin projection and least angel regression. Under the situations of low sound frequency, short or medium sound propagation distance and short or medium length of intervals of geo-acoustic parameters, the polynomial chaos expansion can approximate the transmission losses accurately with the relatively error less than 1%. In the simulation case, with the high signal to noise ratio and the low errors of relative distances between source and receivers, the geo-acoustic parameters can be inverted accurately when the appropriate truncated powers are chosen. And the time cost is reduced by at least an order of magnitude compared with that of traversal grids searching procedure.
      通信作者: 李风华, lfh@mail.ioa.ac.cn
      Corresponding author: Li Feng-Hua, lfh@mail.ioa.ac.cn
    [1]

    尚尔昌 2019 应用声学 038 468Google Scholar

    Shang E C 2019 J. Appl. Acoust. 038 468Google Scholar

    [2]

    Gerstoft P 1994 J. Acoust. Soc. Am. 95 770Google Scholar

    [3]

    Gerstoft P 1998 J. Acoust. Soc. Am. 104 808

    [4]

    Dosso S E 2002 J. Acoust. Soc. Am. 111 129Google Scholar

    [5]

    Dosso S E, Nielsen P L 2002 J. Acoust. Soc. Am. 111 143Google Scholar

    [6]

    Sambridge M 1999 Geophys. J. Int. 138 479Google Scholar

    [7]

    Holland C W, Osler J 2000 J. Acoust. Soc. Am. 107 1263Google Scholar

    [8]

    李整林, 鄢锦, 李风华 2002 声学学报 06 487Google Scholar

    Li Z L, Yan J, Li F H 2002 Acta Acoustica 06 487Google Scholar

    [9]

    李梦竹, 李整林, 周纪浔, 张仁和 2019 物理学报 68 094301Google Scholar

    Li M Z, Li Z L, Zhou J X, Zhang R H 2019 Acta Phys. Sin. 68 094301Google Scholar

    [10]

    Li Z L, Zhang R H, Yan J 2004 IEEE. J. Oceanic. Eng. 29 973

    [11]

    Li Z L, Zhang R H 2004 Chin. Phys. Lett. 21 1100Google Scholar

    [12]

    Li Z L, Zhang R H 2007 Chin. Phys. Lett. 24 471Google Scholar

    [13]

    Li Z L, Li F H 2010 Chin. J. Oceanol. Limnol. 28 990Google Scholar

    [14]

    Xiu D B, Karniadakis G E 2001 SIAM. J. Sci. Comput. 24 619Google Scholar

    [15]

    Finette S 2005 J. Acoust. Soc. Am. 117 997Google Scholar

    [16]

    Finette S 2006 J. Acoust. Soc. Am. 120 2567Google Scholar

    [17]

    Khine Y Y, Creamer D B, Finette S 2011 J. Comput. Acoust. 18 397Google Scholar

    [18]

    程广利, 张明敏 2013 声学学报 38 294

    Cheng G L, Zhang M M 2013 Acta Acoustica 38 294

    [19]

    过武宏, 笪良龙, 赵建昕 2013 应用声学 32 464

    Guo W H, Da L L, Zhao J X 2013 Appl. Acoust. 32 464

    [20]

    笪良龙, 过武宏, 赵建昕 2015 声学学报 40 137

    Da L L, Guo W H, Zhao J X 2015 Acta Acoustica 40 137

    [21]

    笪良龙, 崔宝龙, 过武宏 2017 声学学报 42 25

    Da L L, Cui B L, Guo W L 2017 Acta Acoustica 42 25

    [22]

    LePage K D 2007 9th International Conference on Information Fusion Florence, Italy, July 10–13, 2007 pp1–5

    [23]

    James K R, Dowling D R, 2011 J. Acoust. Soc. Am. 129 589Google Scholar

    [24]

    张鹏, 吴立新 2020 应用声学 40 422

    Zhang P, Wu L X 2020 Appl. Acoust. 40 422

    [25]

    Finette S 2009 J. Acoust. Soc. Am. 126 2242Google Scholar

    [26]

    Gerdes F, Finette S 2012 J. Acoust. Soc. Am. 132 2251Google Scholar

    [27]

    Collins M D 1993 J. Acoust. Soc. Am. 93 1736Google Scholar

    [28]

    Creamer D B 2008 Wave Random Complex 18 197Google Scholar

    [29]

    Choi S K, Canfield R A, Grandhi R V 2004 AIAA J. 42 1191Google Scholar

    [30]

    Blatman G, Sudret B 2011 J. Comput. Phys. 230 2345Google Scholar

    [31]

    Creamer D B 2006 J. Acoust. Soc. Am. 119 1919Google Scholar

    [32]

    Jensen F B, Kuperman W A, Porter M B 2011 Computational Ocean Acoustics (New York: Springer) pp337–455

  • 图 1  不同截断幂次下, 在“频率f-距离r”平面上随机多项式展开平面波声压相对误差1%的等值线. 其中声速的范围是1645−1655 m/s, 吸收率的范围是0.55−0.65 dB/λ

    Fig. 1.  Isolines of 1% relative error about sound pressure for plane wave expanded by the polynomial chaos in frequency-range space. The intervals of sound speed and attenuation are 1645−1655 m/s and 0.55−0.65 dB/λ, respectively.

    图 2  Pekeris波导中100 m深度不同水平距离上的声传播损失, 圈点为传播损失的部分极大值点

    Fig. 2.  Acoustic transmission loss at different horizontal ranges with the depth of 100 m in Pekeris waveguide. Circled points are partial local maximum points.

    图 3  Pekeris波导中100 m深度不同水平距离上随机多项式展开近似声传播损失的验证集误差, 圈点为部分误差极大值点

    Fig. 3.  Validation set errors of acoustic transmission loss expanded by polynomial chaos at different horizontal ranges with the depth of 100 m in Pekeris waveguide. Circled points are partial local maximum points.

    图 4  不同随机多项式截断幂次N、不同深度方向差分网格数Z下, 遍历法与随机多项式展开法计算复杂度之比

    Fig. 4.  Calculation complexity ratio of the ergodic method and the polynomial chaos expansion method under different polynomial truncated power N and difference grid number in depth direction Z.

    图 5  仿真计算使用的海水声速剖面

    Fig. 5.  Sound speed profile used in the simulation case.

    图 6  不同随机多项式展开截断幂次下, (a)海底声速、(b)吸收率、(c)比重反演结果的误差

    Fig. 6.  Errors of geo-acoustic inversion results for (a) sound speed, (b) attenuation and (c) ratio of density under different truncated powers of polynomial chaos.

    图 7  随机多项式展开截断幂次$ N=1 $时的海底声速的反演结果

    Fig. 7.  Geo-acoustic inversion results of sound speed of seabed when the truncated power of polynomial chaos N =1.

    图 8  信噪比为15 dB, 收发水平距离误差为 ± 20 m时, (a)模拟声速、(b)吸收率、(c)比重的反演结果及其置信区间

    Fig. 8.  Geo-acoustic inversion results and its confidential intervals of (a) sound speed, (b) sound attenuation and (c) ratios of density between seabed and sea water when the signal to noise ratio is 15 dB and error of horizontal distance is ± 20 m.

    表 1  Pekeris波导的水文环境和声源参数

    Table 1.  Hydrological conditions of Pekeris waveguide and acoustic source parameters.

    海水深度/m 海水声速/(m·s–1) 声源频率/Hz 发射深度/m 接收深度/m
    100 1500 50 100 100
    下载: 导出CSV

    表 2  接收信号不同信噪比下反演误差均值

    Table 2.  Average inversion errors under different signal to noise ratios.

    信噪比/dB 0 5 10 15 20
    海底声速$ \overline {{\Delta}c}/({\rm{m}}\cdot {\rm{s}}^{-1}) $ 2.63 1.76 1.13 0.71 0.50
    海底吸收率$ \overline {{\Delta}\alpha }/({\rm{dB}}\cdot {\lambda }^{-1}) $ 0.09 0.04 0.02 0.02 0.01
    海底比重$ \overline {{\Delta}\left(\rho \right)} $ 0.10 0.05 0.03 0.02 0.01
    下载: 导出CSV

    表 3  收发水平距离不同误差下反演误差均值

    Table 3.  Average inversion errors under different horizontal distances.

    收发水平距离误差$ {\Delta}r/{\rm{m}} $ 50 30 20 10
    海底声速$ \overline {{\Delta}c}/({\rm{m}}\cdot {\rm{s}}^{-1}) $ 2.49 1.64 1.21 0.65
    海底吸收率$ \overline {{\Delta}\alpha }/({\rm{dB}}\cdot {\lambda }^{-1}) $ 0.04 0.03 0.02 0.01
    海底比重$ \overline {{\Delta}\left(\rho \right)} $ 0.06 0.04 0.03 0.01
    下载: 导出CSV

    表 4  使用遍历法与随机多项式展开法反演海底参数所需计算时间

    Table 4.  Calculation time of geo-acoustic inversion using method of traversing and polynomial chaos expansion.

    遍历法计算时间/min 随机多项式展开法反演计算时间/min
    截断幂次N
    1 2 3 4 5 6 7 8 9 10
    7047 2 9 10 77 199 574 2254 6130 18783 57547
    下载: 导出CSV
  • [1]

    尚尔昌 2019 应用声学 038 468Google Scholar

    Shang E C 2019 J. Appl. Acoust. 038 468Google Scholar

    [2]

    Gerstoft P 1994 J. Acoust. Soc. Am. 95 770Google Scholar

    [3]

    Gerstoft P 1998 J. Acoust. Soc. Am. 104 808

    [4]

    Dosso S E 2002 J. Acoust. Soc. Am. 111 129Google Scholar

    [5]

    Dosso S E, Nielsen P L 2002 J. Acoust. Soc. Am. 111 143Google Scholar

    [6]

    Sambridge M 1999 Geophys. J. Int. 138 479Google Scholar

    [7]

    Holland C W, Osler J 2000 J. Acoust. Soc. Am. 107 1263Google Scholar

    [8]

    李整林, 鄢锦, 李风华 2002 声学学报 06 487Google Scholar

    Li Z L, Yan J, Li F H 2002 Acta Acoustica 06 487Google Scholar

    [9]

    李梦竹, 李整林, 周纪浔, 张仁和 2019 物理学报 68 094301Google Scholar

    Li M Z, Li Z L, Zhou J X, Zhang R H 2019 Acta Phys. Sin. 68 094301Google Scholar

    [10]

    Li Z L, Zhang R H, Yan J 2004 IEEE. J. Oceanic. Eng. 29 973

    [11]

    Li Z L, Zhang R H 2004 Chin. Phys. Lett. 21 1100Google Scholar

    [12]

    Li Z L, Zhang R H 2007 Chin. Phys. Lett. 24 471Google Scholar

    [13]

    Li Z L, Li F H 2010 Chin. J. Oceanol. Limnol. 28 990Google Scholar

    [14]

    Xiu D B, Karniadakis G E 2001 SIAM. J. Sci. Comput. 24 619Google Scholar

    [15]

    Finette S 2005 J. Acoust. Soc. Am. 117 997Google Scholar

    [16]

    Finette S 2006 J. Acoust. Soc. Am. 120 2567Google Scholar

    [17]

    Khine Y Y, Creamer D B, Finette S 2011 J. Comput. Acoust. 18 397Google Scholar

    [18]

    程广利, 张明敏 2013 声学学报 38 294

    Cheng G L, Zhang M M 2013 Acta Acoustica 38 294

    [19]

    过武宏, 笪良龙, 赵建昕 2013 应用声学 32 464

    Guo W H, Da L L, Zhao J X 2013 Appl. Acoust. 32 464

    [20]

    笪良龙, 过武宏, 赵建昕 2015 声学学报 40 137

    Da L L, Guo W H, Zhao J X 2015 Acta Acoustica 40 137

    [21]

    笪良龙, 崔宝龙, 过武宏 2017 声学学报 42 25

    Da L L, Cui B L, Guo W L 2017 Acta Acoustica 42 25

    [22]

    LePage K D 2007 9th International Conference on Information Fusion Florence, Italy, July 10–13, 2007 pp1–5

    [23]

    James K R, Dowling D R, 2011 J. Acoust. Soc. Am. 129 589Google Scholar

    [24]

    张鹏, 吴立新 2020 应用声学 40 422

    Zhang P, Wu L X 2020 Appl. Acoust. 40 422

    [25]

    Finette S 2009 J. Acoust. Soc. Am. 126 2242Google Scholar

    [26]

    Gerdes F, Finette S 2012 J. Acoust. Soc. Am. 132 2251Google Scholar

    [27]

    Collins M D 1993 J. Acoust. Soc. Am. 93 1736Google Scholar

    [28]

    Creamer D B 2008 Wave Random Complex 18 197Google Scholar

    [29]

    Choi S K, Canfield R A, Grandhi R V 2004 AIAA J. 42 1191Google Scholar

    [30]

    Blatman G, Sudret B 2011 J. Comput. Phys. 230 2345Google Scholar

    [31]

    Creamer D B 2006 J. Acoust. Soc. Am. 119 1919Google Scholar

    [32]

    Jensen F B, Kuperman W A, Porter M B 2011 Computational Ocean Acoustics (New York: Springer) pp337–455

  • [1] 刘广凯, 全厚德, 康艳梅, 孙慧贤, 崔佩璋, 韩月明. 一种随机共振增强正弦信号的二次多项式接收方法. 物理学报, 2019, 68(21): 210501. doi: 10.7498/aps.68.20190952
    [2] 江鹏飞, 林建恒, 孙军平, 衣雪娟. 考虑噪声源深度分布的海洋环境噪声模型及地声参数反演. 物理学报, 2017, 66(1): 014306. doi: 10.7498/aps.66.014306
    [3] 路文龙, 谢军伟, 王和明, 盛川. 基于多项式调频Fourier变换的信号分量提取方法. 物理学报, 2016, 65(8): 080202. doi: 10.7498/aps.65.080202
    [4] 臧鸿雁, 柴宏玉. 一个二次多项式混沌系统的均匀化及其熵分析. 物理学报, 2016, 65(3): 030504. doi: 10.7498/aps.65.030504
    [5] 郭晓乐, 杨坤德, 马远良. 一种基于简正波模态频散的远距离宽带海底参数反演方法. 物理学报, 2015, 64(17): 174302. doi: 10.7498/aps.64.174302
    [6] 冯菊, 廖成, 张青洪, 盛楠, 周海京. 蒸发波导中的时间反演抛物方程定位法. 物理学报, 2014, 63(13): 134101. doi: 10.7498/aps.63.134101
    [7] 曾喆昭, 雷妮, 盛立锃. 不确定混沌系统的多项式函数模型补偿控制. 物理学报, 2013, 62(15): 150506. doi: 10.7498/aps.62.150506
    [8] 范洪义, 楼森岳, 潘孝胤, 笪诚. 涉及Hermite多项式的二项式定理和Laguerre多项式的负二项式定理. 物理学报, 2013, 62(24): 240301. doi: 10.7498/aps.62.240301
    [9] 屈科, 胡长青, 赵梅. 利用时域波形快速反演海底单参数的方法. 物理学报, 2013, 62(22): 224303. doi: 10.7498/aps.62.224303
    [10] 杨珺, 孙秋野, 杨东升. 基于多项式模型的混沌系统平方和算法脉冲控制. 物理学报, 2012, 61(20): 200511. doi: 10.7498/aps.61.200511
    [11] 刘伟伟, 任煜轩, 高红芳, 孙晴, 王自强, 李银妹. 泽尼克多项式校正全息阵列光镊像差的实验研究. 物理学报, 2012, 61(18): 188701. doi: 10.7498/aps.61.188701
    [12] 杨坤德, 马远良. 利用海底反射信号进行地声参数反演的方法. 物理学报, 2009, 58(3): 1798-1805. doi: 10.7498/aps.58.1798
    [13] 闫 华, 魏 平, 肖先赐. 基于Bernstein多项式的自适应混沌时间序列预测算法. 物理学报, 2007, 56(9): 5111-5118. doi: 10.7498/aps.56.5111
    [14] 周永道, 马 洪, 吕王勇, 王会琦. 基于多元局部多项式方法的混沌时间序列预测. 物理学报, 2007, 56(12): 6809-6814. doi: 10.7498/aps.56.6809
    [15] 徐 进, 王文祥, 岳玲娜, 宫玉彬, 魏彦玉. 脊加载椭圆波导的多项式表示分析方法. 物理学报, 2007, 56(11): 6393-6397. doi: 10.7498/aps.56.6393
    [16] 吴 楚. 多项式角动量代数的代数表示及实现. 物理学报, 2006, 55(6): 2676-2681. doi: 10.7498/aps.55.2676
    [17] 马少娟, 徐 伟, 李 伟. 基于Laguerre多项式逼近法的随机双势阱Duffing系统的分岔和混沌研究. 物理学报, 2006, 55(8): 4013-4019. doi: 10.7498/aps.55.4013
    [18] 马少娟, 徐 伟, 李 伟, 靳艳飞. 基于Chebyshev多项式逼近的随机 van der Pol系统的倍周期分岔分析. 物理学报, 2005, 54(8): 3508-3515. doi: 10.7498/aps.54.3508
    [19] 郭奇志, 谭维翰, 孟义朝. 推广的Tschebyshev多项式及其在解强耦合波导方程中的应用. 物理学报, 2003, 52(12): 3092-3097. doi: 10.7498/aps.52.3092
    [20] 叶开沅. 应力函数或扭曲函数爲三次多项式形式的扭转问题. 物理学报, 1953, 9(4): 255-274. doi: 10.7498/aps.9.255
计量
  • 文章访问数:  3091
  • PDF下载量:  63
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-19
  • 修回日期:  2021-03-23
  • 上网日期:  2021-06-07
  • 刊出日期:  2021-09-05

/

返回文章
返回