搜索

文章查询

x

留言板

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

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

利用海洋环境噪声空间特性估计浅海海底分层结构及地声参数

李赫 郭新毅 马力

利用海洋环境噪声空间特性估计浅海海底分层结构及地声参数

李赫, 郭新毅, 马力
PDF
HTML
导出引用
导出核心图
  • 海洋环境噪声场中包含了海洋中的诸多信息, 海底地声参数是影响海洋环境噪声场空间分布的主要因素之一. 对于不同的海底分层结构, 海底反射损失会根据沉积层厚度和各层声速呈现出不同的临界角和干涉条纹结构. 本文利用Harrison能流理论, 从理想反射系数出发, 分别考虑了声速、密度、衰减系数、沉积层厚度等几种参数对无沉积层和单层沉积层中反射系数的影响, 并对单层沉积层海底的反射系数进行了化简, 结合互易原理解释了反射损失条纹结构的形成机理. 中国黄海某海区试验结果表明, 利用海洋环境噪声空间方向谱获得的海底反射损失, 可以提取海底反射临界角和干涉条纹信息, 由此可估计出海底分层结构、声速和沉积层厚度等海底参数信息.
      通信作者: 郭新毅, guoxinyi@mail.ioa.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 11874384)资助的课题
    [1]

    Wenz G M 1962 J. Acoust. Soc. Am. 34 1936

    [2]

    Kuperman W A, Ingenito F 1980 J. Acoust. Soc. Am. 67 1988

    [3]

    HarrisonC H 1996 J. Acoust. Soc. Am. 99 2055

    [4]

    林建恒 2002 博士学位论文 (北京: 中国科学院研究生院)

    Lin J H 2002 Ph. D. Dissertation (Beijing: Graduate University of Chinese Academy of Science) (in Chinese)

    [5]

    张乾初 2017 博士学位论文 (北京: 中国科学院大学)

    Zhang Q C 2017 Ph. D. Dissertation (Beijing: University of Chinese Academy of Science) (in Chinese)

    [6]

    郭新毅, 铁广鹏, 李凡, 马力 2014 物理 43 723

    Guo X Y, Tie G P, Li F, Ma L 2014 Physics 43 723

    [7]

    江鹏飞, 林建恒, 孙军平, 衣雪娟 2017 物理学报 66 014306

    Jiang P F, Lin J H, Sun J P, Yi X J 2017 Acta Phys. Sin. 66 014306

    [8]

    蒋光禹, 孙超, 谢磊, 刘雄厚 2019 物理学报 68 024302

    Jiang G Y, Sun C, Xie L, Liu X H 2019 Acta Phys. Sin. 68 024302

    [9]

    Zeng J, Chapman N R, Bonnel J 2013 J. Acoust. Soc. Am. 134 EL394

    [10]

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

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

    [11]

    周建波, 朴胜春, 刘亚琴, 张海刚, 屈科, 李晓曼, 张明辉 2019 声学学报 44 337

    Zhou J B, Piao S C, Liu Y Q, Zhang H G, Qu K, Li X M, Zhang M H 2019 Acta Acoust. 44 337

    [12]

    江鹏飞, 林建恒, 马力, 殷宝友, 蒋国建 2016 声学学报 41 59

    Jiang P F, Lin J H, Ma Li, Yin B Y, Jiang G J 2016 Acta Acoust. 41 59

    [13]

    骆文于 2002 博士学位论文 (北京: 中国科学院声学研究所)

    Luo W Y 2002 Ph. D. Dissertation (Beijing: Institute of Acoustics, Chinese Academy of Science) (in Chinese) (in Chinese)

    [14]

    布列霍夫斯基Л М 著 (山东省海洋学院海洋物理系, 中国科学院声学研究所水声研究室 译) 1983 海洋声学 (北京: 科学出版社) 第516−520页

    Бреховских Л М (translated by Department of Oceanophysics Shandong College of Oceanology, Laboratory of Underwater Acoustic Institute of Acoustics Chinese Academy of Science) 1983 Fundamentals of Ocean Acoustics (Beijing: Science Press) pp516−520 (in Chinese)

    [15]

    Harrison C H, Simons D G 2002 J. Acoust. Soc. Am. 112 1377

    [16]

    Harrison C H 2004 J. Acoust. Soc. Am. 115 1505

    [17]

    Siderius M, Harrison C H, Porter M B 2006 J. Acoust. Soc. Am. 120 1315

    [18]

    Muzi L, Siderius M, Quijano J E, Dosso S E 2015 J. Acoust. Soc. Am. 137 481

    [19]

    Muzi L, Siderius M, Nielsen P L 2016 J. Acoust. Soc. Am. 140 1513

    [20]

    Schmidt H 2004 OASES version 3.1 User Guide and Reference Manual (Cambridge, MA: Massachusetts Institute of Technology)

    [21]

    Jensen F B, KupermanW A, Porter M B, Schmidt H 2011 Computational Ocean Acoustics (2nd Ed.) (New York: Springer) pp38−43

    [22]

    Zeng J, Zhao Z D, Peng D Y, Ma L, Shang E C 2019 J. Acoust. Soc. Am. 145 EL243

  • 图 1  浅海中利用垂直阵接收噪声空间结构

    Fig. 1.  Spatial structure of receiving ocean ambient noise using vertical line array in shallow water

    图 2  (a) 真实反射损失BL; (b) 根据噪声垂直空间指向性利用OASN模块仿真得到的$\widehat {BL}$; (c) 1800 Hz下两种方法的比较, 实线为BL, 虚线为$\widehat {BL}$

    Fig. 2.  (a) True BL; (b) $\widehat {BL}$ computed by vertical directionality of ocean ambient noise using OASN; (c) two methods compare under 1800 Hz, the full line is BL, the imaginary line is $\widehat {BL}$

    图 3  1800 Hz下不同海底分层结构下的$\widehat {BL}$(实线为阵元遍布全海深, 点划线为42阵元, 间隔均为0.2 m) (a) 海底为无限大液体声学半空间; (b) 海底为单层沉积层

    Fig. 3.  The $\widehat {BL}$ of different structure of sub-bottom under 1800 Hz: (a) Infinite acoustic half space; (b) 1 layer of sediment. The full line corresponds to the condition that the elements set across the sea. The imaginary line corresponds to the condition that 42 elements set at the depth of 30 m

    图 4  海底分层模型 (a)液体无限大声学半空间海底; (b) 存在一层沉积层海底

    Fig. 4.  Model of sub-bottom stratification: (a) Infinite acoustic half space; (b) 1 layer of sediment

    图 5  Pekeris分支割线、极点和积分的复波数平面

    Fig. 5.  Pekeris branch cut map

    图 6  不同谱域内噪声场的指向性和$\widehat {BL}$ (a), (b)连续谱部分噪声场; (c), (d) 离散谱部分噪声场; (e), (f)全谱域噪声场

    Fig. 6.  Vertical directivity and $\widehat {BL}$ of ocean ambient noise in different spectral domain: (a), (b) Near field continuous spectrum; (c), (d) far field discrete spectrum; (e), (f) full spectrum domain

    图 7  不同海底参数下的BL (a)衰减系数$\alpha $; (b)声速c; (c)密度$\rho $

    Fig. 7.  The BL under different parameters of sub-bottom: (a) Attenuation coefficient $\alpha $; (b) sound speed c; (c) density $\rho $

    图 8  单层沉积层海底反射模型

    Fig. 8.  Reflection model of sub-bottom with 1 layer of sediment

    图 9  (a) (8)式计算的真实BL; (b) (10)式计算化简后的BL; (c) 1800 Hz下两条BL曲线

    Fig. 9.  (a) True BL computed by equation (8); (b) BL after simplified computed by equation (10); (c) the curves of BL of two methods under 1800 Hz

    图 10  (10)式计算化简后的BL (横轴利用(7)式转为${\rm{sin}}{\theta _2}$)

    Fig. 10.  Curve of BL after simplified computed by equation (10) with the horizontal axis converted to ${\rm{sin}}{\theta _2}$ by equation (7)

    图 11  (a) 42阵元设置在水下30 m处$\widehat {BL}$; (b) 1800 Hz下$\widehat {BL}$曲线

    Fig. 11.  (a) $\widehat {BL}$ with 42 elements at the depth of 30 m; (b) the curve of $\widehat {BL}$ under 1800 Hz.

    图 12  实验实测海水声速剖面

    Fig. 12.  Sound speed profile of the experiment

    图 13  第9组数据处理结果 (a) $\widehat {BL}$; (b) 1600 Hz下的$\widehat {BL}$曲线

    Fig. 13.  Results of the 9th set of data: (a) $\widehat {BL}$; (b) curve of $\widehat {BL}$ under 1600 Hz

    图 14  第11组数据处理结果 (a) $\widehat {BL}$; (b) 1600 Hz下$\widehat {BL}$曲线

    Fig. 14.  Results of the 11th set of data: (a) $\widehat {BL}$; (b) curve of $\widehat {BL}$ under 1600 Hz

    图 15  数据计算结果统计分布直方图 (a)沉积层声速; (b)基底声速; (c)沉积层厚度

    Fig. 15.  Statistical distribution histogram of 20 sets of data: (a) Sound speed of the sediment; (b) sound speed of the basement; (c) thickness of the sediment

    图 16  浅地层剖面仪探测试验海区海底结构(图(b)中用黑点标明海底分界面)

    Fig. 16.  Sub-bottom structure of the experimental area detected by sub-bottom profiler. The interface is marked by black spot in panel (b)

    表 1  海水及海底声学参数

    Table 1.  Parameters of ocean and sub-bottom

    介质
    类型
    声速
    c/m·s–1
    衰减系数
    α/dB·λ–1
    对水
    密度ρ
    深度
    H/m
    海水 1500 0 1 40
    沉积层 1600 0.2 1.5 3
    基底 2000 0.6 2
    下载: 导出CSV
  • [1]

    Wenz G M 1962 J. Acoust. Soc. Am. 34 1936

    [2]

    Kuperman W A, Ingenito F 1980 J. Acoust. Soc. Am. 67 1988

    [3]

    HarrisonC H 1996 J. Acoust. Soc. Am. 99 2055

    [4]

    林建恒 2002 博士学位论文 (北京: 中国科学院研究生院)

    Lin J H 2002 Ph. D. Dissertation (Beijing: Graduate University of Chinese Academy of Science) (in Chinese)

    [5]

    张乾初 2017 博士学位论文 (北京: 中国科学院大学)

    Zhang Q C 2017 Ph. D. Dissertation (Beijing: University of Chinese Academy of Science) (in Chinese)

    [6]

    郭新毅, 铁广鹏, 李凡, 马力 2014 物理 43 723

    Guo X Y, Tie G P, Li F, Ma L 2014 Physics 43 723

    [7]

    江鹏飞, 林建恒, 孙军平, 衣雪娟 2017 物理学报 66 014306

    Jiang P F, Lin J H, Sun J P, Yi X J 2017 Acta Phys. Sin. 66 014306

    [8]

    蒋光禹, 孙超, 谢磊, 刘雄厚 2019 物理学报 68 024302

    Jiang G Y, Sun C, Xie L, Liu X H 2019 Acta Phys. Sin. 68 024302

    [9]

    Zeng J, Chapman N R, Bonnel J 2013 J. Acoust. Soc. Am. 134 EL394

    [10]

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

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

    [11]

    周建波, 朴胜春, 刘亚琴, 张海刚, 屈科, 李晓曼, 张明辉 2019 声学学报 44 337

    Zhou J B, Piao S C, Liu Y Q, Zhang H G, Qu K, Li X M, Zhang M H 2019 Acta Acoust. 44 337

    [12]

    江鹏飞, 林建恒, 马力, 殷宝友, 蒋国建 2016 声学学报 41 59

    Jiang P F, Lin J H, Ma Li, Yin B Y, Jiang G J 2016 Acta Acoust. 41 59

    [13]

    骆文于 2002 博士学位论文 (北京: 中国科学院声学研究所)

    Luo W Y 2002 Ph. D. Dissertation (Beijing: Institute of Acoustics, Chinese Academy of Science) (in Chinese) (in Chinese)

    [14]

    布列霍夫斯基Л М 著 (山东省海洋学院海洋物理系, 中国科学院声学研究所水声研究室 译) 1983 海洋声学 (北京: 科学出版社) 第516−520页

    Бреховских Л М (translated by Department of Oceanophysics Shandong College of Oceanology, Laboratory of Underwater Acoustic Institute of Acoustics Chinese Academy of Science) 1983 Fundamentals of Ocean Acoustics (Beijing: Science Press) pp516−520 (in Chinese)

    [15]

    Harrison C H, Simons D G 2002 J. Acoust. Soc. Am. 112 1377

    [16]

    Harrison C H 2004 J. Acoust. Soc. Am. 115 1505

    [17]

    Siderius M, Harrison C H, Porter M B 2006 J. Acoust. Soc. Am. 120 1315

    [18]

    Muzi L, Siderius M, Quijano J E, Dosso S E 2015 J. Acoust. Soc. Am. 137 481

    [19]

    Muzi L, Siderius M, Nielsen P L 2016 J. Acoust. Soc. Am. 140 1513

    [20]

    Schmidt H 2004 OASES version 3.1 User Guide and Reference Manual (Cambridge, MA: Massachusetts Institute of Technology)

    [21]

    Jensen F B, KupermanW A, Porter M B, Schmidt H 2011 Computational Ocean Acoustics (2nd Ed.) (New York: Springer) pp38−43

    [22]

    Zeng J, Zhao Z D, Peng D Y, Ma L, Shang E C 2019 J. Acoust. Soc. Am. 145 EL243

  • [1] 黎雪刚, 杨坤德, 张同伟, 邱海宾. 基于拖曳倾斜线列阵的海底反射损失提取方法. 物理学报, 2009, 58(11): 7741-7749. doi: 10.7498/aps.58.7741
    [2] 屈科, 胡长青, 赵梅. 利用时域波形快速反演海底单参数的方法. 物理学报, 2013, 62(22): 224303. doi: 10.7498/aps.62.224303
    [3] 周天, 李海森, 朱建军, 魏玉阔. 利用多角度海底反向散射信号进行地声参数估计. 物理学报, 2014, 63(8): 084302. doi: 10.7498/aps.63.084302
    [4] 江鹏飞, 林建恒, 孙军平, 衣雪娟. 考虑噪声源深度分布的海洋环境噪声模型及地声参数反演. 物理学报, 2017, 66(1): 014306. doi: 10.7498/aps.66.014306
    [5] 杨坤德, 马远良. 利用海底反射信号进行地声参数反演的方法. 物理学报, 2009, 58(3): 1798-1805. doi: 10.7498/aps.58.1798
    [6] 张鹏, 李整林, 吴立新, 张仁和, 秦继兴. 深海海底反射会聚区声传播特性. 物理学报, 2019, 68(1): 014301. doi: 10.7498/aps.68.20181761
    [7] 李倩倩, 阳凡林, 张凯, 郑炳祥. 不确定海洋环境中基于贝叶斯理论的声源运动参数估计方法. 物理学报, 2016, 65(16): 164304. doi: 10.7498/aps.65.164304
    [8] 夏麾军, 马远良, 刘亚雄. 海洋环境噪声场对称性分析及噪声消除方法. 物理学报, 2016, 65(14): 144302. doi: 10.7498/aps.65.144302
    [9] 胡治国, 李整林, 张仁和, 任云, 秦继兴, 何利. 深海海底斜坡环境下的声传播. 物理学报, 2016, 65(1): 014303. doi: 10.7498/aps.65.014303
    [10] 李晟昊, 李整林, 李文, 秦继兴. 深海海底山环境下声传播水平折射效应研究. 物理学报, 2018, 67(22): 224302. doi: 10.7498/aps.67.20181480
    [11] 赵现斌, 严卫, 王迎强, 陆文, 马烁. 基于海面散射模型的全极化合成孔径雷达海洋环境探测关键技术参数设计仿真研究. 物理学报, 2014, 63(21): 218401. doi: 10.7498/aps.63.218401
    [12] 刘宗伟, 孙超, 向龙凤, 易锋. 不确定海洋环境中的模态子空间重构稳健定位方法. 物理学报, 2014, 63(3): 034304. doi: 10.7498/aps.63.034304
    [13] 陈桂波, 汪宏年, 姚敬金, 韩子夜. 各向异性海底地层海洋可控源电磁响应三维积分方程法数值模拟. 物理学报, 2009, 58(6): 3848-3857. doi: 10.7498/aps.58.3848
    [14] 侯倩男, 吴金荣. 浅海小掠射角的海底界面声反向散射模型的简化. 物理学报, 2019, 68(4): 044301. doi: 10.7498/aps.68.20181475
    [15] 李梦竹, 李整林, 周纪浔, 张仁和. 一种低声速沉积层海底参数声学反演方法. 物理学报, 2019, 68(9): 094301. doi: 10.7498/aps.68.20190183
    [16] 李扬, 郭树旭. 基于稀疏分解的大功率半导体激光器1/f噪声参数估计的新方法. 物理学报, 2012, 61(3): 034208. doi: 10.7498/aps.61.034208
    [17] 王龙昊, 秦继兴, 傅德龙, 李整林, 刘建军, 翁晋宝. 深海大接收深度海底混响研究. 物理学报, 2019, 68(13): 134303. doi: 10.7498/aps.68.20181883
    [18] 刘宗伟, 孙超, 杜金燕. 不确定海洋声场中的检测性能损失环境敏感度度量. 物理学报, 2013, 62(6): 064303. doi: 10.7498/aps.62.064303
    [19] 郭晓乐, 杨坤德, 马远良. 一种基于简正波模态频散的远距离宽带海底参数反演方法. 物理学报, 2015, 64(17): 174302. doi: 10.7498/aps.64.174302
    [20] 汤立国, 许肖梅, 刘胜兴. 海底爆破辐射声场的理论及数值研究. 物理学报, 2008, 57(7): 4251-4257. doi: 10.7498/aps.57.4251
  • 引用本文:
    Citation:
计量
  • 文章访问数:  228
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-28
  • 修回日期:  2019-08-13
  • 上网日期:  2019-11-26
  • 刊出日期:  2019-11-01

利用海洋环境噪声空间特性估计浅海海底分层结构及地声参数

  • 1. 中国科学院声学研究所, 中国科学院水声环境特性重点实验室, 北京 100190
  • 2. 中国科学院大学, 北京 100049
  • 通信作者: 郭新毅, guoxinyi@mail.ioa.ac.cn
    基金项目: 国家自然科学基金(批准号: 11874384)资助的课题

摘要: 海洋环境噪声场中包含了海洋中的诸多信息, 海底地声参数是影响海洋环境噪声场空间分布的主要因素之一. 对于不同的海底分层结构, 海底反射损失会根据沉积层厚度和各层声速呈现出不同的临界角和干涉条纹结构. 本文利用Harrison能流理论, 从理想反射系数出发, 分别考虑了声速、密度、衰减系数、沉积层厚度等几种参数对无沉积层和单层沉积层中反射系数的影响, 并对单层沉积层海底的反射系数进行了化简, 结合互易原理解释了反射损失条纹结构的形成机理. 中国黄海某海区试验结果表明, 利用海洋环境噪声空间方向谱获得的海底反射损失, 可以提取海底反射临界角和干涉条纹信息, 由此可估计出海底分层结构、声速和沉积层厚度等海底参数信息.

English Abstract

    • 海洋环境噪声是海洋中永恒存在的声场, 主要由风浪、降雨、舰船、海洋生物、工业等因素形成, 其中风成噪声在各个海域普遍存在, 其频段覆盖几百至几十千赫兹[1]. 一方面海洋环境噪声是影响声呐工作性能的主要因素之一, 只有对其进行充分的了解研究才能充分提高水下设备的工作性能, 所以国内外学者纷纷对海洋环境噪声进行建模研究[2-5]. 另一方面, 海洋环境噪声场中也包含了诸多水体、海底、海面等环境信息[6], 因此利用海洋环境噪声进行参数反演, 获取海水、海底等环境信息也成为时下的热点问题之一[7,8].

      海底地声参数反演主要分为主动反演和被动反演两部分. Zeng等[9]利用爆炸声源对中国黄海海域海底衰减系数进行反演; 李梦竹等[10]提出了一种适用于低声速沉积层的海底参数声学反演方法. 相比于主动反演, 利用海洋环境噪声进行被动反演无需主动发射声源信号, 在实验中仅需数据接收装置, 大幅度节省了实验消耗以及工作任务, 并且对海洋生物没有任何影响. 周建波等[11]利用矢量水听器通过海洋环境噪声提取了声场格林函数; 江鹏飞等[12]根据海洋环境噪声的空间指向性在不同掠射角范围内对海底参数的敏感度不同, 提出了分步反演的方法; 骆文于[13]利用中国东海环境噪声数据反演了液态半空间海底声速、密度、衰减系数; 早在1982年, 布列霍夫斯基[14]就曾利用海洋环境噪声场的各向异性, 来获取海底反射信息; Harrison在2002年进一步系统地研究了海洋环境噪声中的能流理论, 利用上行波和下行波的比值提取了海底反射损失, 并对海底进行了反演[15], 并利用希尔伯特变换补全反射系数的虚部, 通过计算垂直方向的冲击响应利用漂流浮标对海底进行成像[16]; Siderius等[17]利用海洋环境噪声提取了格林函数用来表征两点间的冲击响应, 对海底分层进行成像; Muzi等[18,19]在Harrison的基础上利用互谱密度矩阵的特性对阵元进行了理论扩展处理, 提高了利用海洋环境噪声提取的海底损失曲线的分辨率.

      本文在Harrison能流理论基础上, 从理想环境反射系数出发, 解释了不同海底参数对海底反射损失(bottom-loss, BL)曲线的影响, 并利用实测海洋环境噪声数据对海底分层结构、声速及沉积层厚度进行估计. 本文首先引入Harrison能流理论估计海底反射损失的方法, 分析了不同海底分层情况下的曲线特点; 然后利用简正波理论, 讨论了浅海波导中近场连续谱和远场离散谱在噪声估计海底反射损失中的贡献; 并从理想反射系数出发, 分别讨论了不同海底参数对无沉积层、单一海底沉积层两种情况的影响, 结合互易原理解释了反射损失条纹结构的产生机理, 并根据噪声获得的海底反射损失曲线给出海底声速、沉积层厚度的估计方法; 最后介绍了 2016年冬中国黄海某海区海上试验, 对实测海洋环境噪声数据进行处理, 估计其海底分层结构、海底声速和沉积层厚度.

    • 从浅海理想波导出发, 假设噪声源均匀分布在海面, 利用垂直阵接收环境噪声信号, 其空间结构如图1所示. 波导中的声波可分为两部分: 接收阵上方的来波称为下行波, 从噪声源直达接收阵; 下方的称为上行波, 经过海底一次反射后到达接收阵. 在此近场情况下不考虑多次反射, 后文将给出解释.

      图  1  浅海中利用垂直阵接收噪声空间结构

      Figure 1.  Spatial structure of receiving ocean ambient noise using vertical line array in shallow water

      声速剖面为等声速的理想情况下, 表面声源的出射角${\theta _{\rm{s}}}$、海底掠射角${\theta _{\rm{b}}}$与接收角${\theta _{\rm{r}}}$均相等, 在此均用掠射角${\theta _{\rm{b}}}$表示. 对于某一角频率$\omega $下, 海底反射损失BL有如下定义:

      $ BL\left( {{\theta _{\rm{b}}},\omega } \right) = - 10{\rm{lo}}{{\rm{g}}_{10}}\left( {{{\left| {R\left( {{\theta _{\rm{b}}},\omega } \right)} \right|}^2}} \right). $

      下文用$\widehat {BL}$表示利用噪声估计出的海底反射损失, BL表示理想反射系数得到的海底反射损失, 以便区分. 这里R为平面波的海底反射系数, 其幅度可以通过波导中的上行波与下行波能量比值求得. 对垂直阵接收的信号进行波束形成以获得不同俯仰角噪声能量, 计算公式如下:

      $ {\left| {\hat R\left( {{\theta _{\rm{b}}},\omega } \right)} \right|^2} = \frac{{B\left( { - {\theta _{\rm{b}}},\omega } \right)}}{{B\left( {{\theta _{\rm{b}}},\omega } \right)}}, $

      其中B为波束形成得到的不同角度和角频率下的接收能量, 即垂直方向能量谱,

      $ B\left( {{\theta _{\rm{b}}},\omega } \right) = E\left( {{{{w}}^{\rm{H}}}{{{C}}_\omega }{{w}}} \right) = {{{w}}^{\rm{H}}}E\left( {{{{p}}_\omega } \cdot {{p}}_\omega ^{\rm{H}}} \right){{w}}, $

      其中E代表期望; H为共轭转置; w是波束形成的加权向量; p为对应角频率$\omega $的接收数据, 在M阵元接收时, 为M维向量; ${{{C}}_\omega }$为噪声互谱密度矩阵.

      为说明此方法计算得到BL的有效性, 先考虑极端情况, 假设阵元间隔为0.2 m的垂直阵遍布在声速为1500 m/s、水深为40 m的全海深中, 忽略海底横波的影响, 海底参数如表1所列.

      介质
      类型
      声速
      c/m·s–1
      衰减系数
      α/dB·λ–1
      对水
      密度ρ
      深度
      H/m
      海水 1500 0 1 40
      沉积层 1600 0.2 1.5 3
      基底 2000 0.6 2

      表 1  海水及海底声学参数

      Table 1.  Parameters of ocean and sub-bottom

      互谱密度矩阵直接利用OASES软件中的OASN模块计算获得[20], w设置为常规波束形成加权向量. 图2(a)为根据理想情况下真实反射系数公式求得的BL[21], 后文将对该公式进行讨论; 图2(b)是根据噪声垂直空间指向性利用OASN模块仿真得到的$\widehat {BL}$; 图2(c)提取出了1800 Hz下两种方法的掠射角BL曲线. 图2(a)图2(b)中的条纹与图2(c)中的曲线峰值成对应关系, 阵元遍布全海深的理想状态下, 图2(c)中两条曲线的峰值一一对应, 仅在幅度上存在一定偏差, 因此利用Harrison能流理论估计海底反射损失是可行的.

      图  2  (a) 真实反射损失BL; (b) 根据噪声垂直空间指向性利用OASN模块仿真得到的$\widehat {BL}$; (c) 1800 Hz下两种方法的比较, 实线为BL, 虚线为$\widehat {BL}$

      Figure 2.  (a) True BL; (b) $\widehat {BL}$ computed by vertical directionality of ocean ambient noise using OASN; (c) two methods compare under 1800 Hz, the full line is BL, the imaginary line is $\widehat {BL}$

      而在实际海上试验中, 阵元个数是有限的, 很难达到上述遍布全海深的理想状态. 有限的阵元个数会带来不同掠射角方向分辨率的下降, 这样会导致图2(b)中条纹信息一定程度上的缺失, 同时也使图2(c)中峰值的模糊, 但仍然有很多特征可以利用. 图3给出了两种不同海底分层结构的$\widehat {BL}$, 保留表1中的参数, 图中实线为间隔0.2 m的水听器遍布全海深的结果, 虚线为间隔0.2 m的42元水听器阵, 第一个阵元设置在水深30 m处, 频率均为1800 Hz. 海底模型如图4所示, 其中图4(a)为无限大液态声学半空间海底, 图4(b)为带有一层沉积层的海底.

      图  3  1800 Hz下不同海底分层结构下的$\widehat {BL}$(实线为阵元遍布全海深, 点划线为42阵元, 间隔均为0.2 m) (a) 海底为无限大液体声学半空间; (b) 海底为单层沉积层

      Figure 3.  The $\widehat {BL}$ of different structure of sub-bottom under 1800 Hz: (a) Infinite acoustic half space; (b) 1 layer of sediment. The full line corresponds to the condition that the elements set across the sea. The imaginary line corresponds to the condition that 42 elements set at the depth of 30 m

      图  4  海底分层模型 (a)液体无限大声学半空间海底; (b) 存在一层沉积层海底

      Figure 4.  Model of sub-bottom stratification: (a) Infinite acoustic half space; (b) 1 layer of sediment

      当海底不存在沉积层时, $\widehat {BL}$曲线中只存在一次突变, 即存在一处临界角, 介质中仅存在一层分界面, 小掠射角时会发生全内反射, 而当掠射角大于临界角时, 发生折射, 能量泄漏传播到海底, 导致$\widehat {BL}$突增; 而存在一层沉积层时, 介质中有两个分界面, 在每层介质物理性质差异明显的情况下, 会产生两个临界角, 这是由于当掠射角略微大于第一个临界角时, 折射到沉积层中的声线会发生偏转, 通常情况下下层介质声速大于上层介质声速, 声线会向上偏转, 使得声波在沉积层中的掠射角变小, 而随着入射波的掠射角再次增大, 折射波的掠射角会达到第二个临界角, 因此图3(b)中会存在两次幅度上的突变, 而三层及多层海底结构复杂, 会引起第二个临界角不明显, 本文不予考虑. 可以看出, 利用噪声垂直方向谱获得的$\widehat {BL}$曲线可以用来估计海底的分层结构.

    • 在浅海波导中, 点声源在大掠射角范围内发出的声波很难进行远距离传播, 因为这一部分经过海底反射存在能量损失, 反射波携带着一定海底信息, 这一特性对于均匀分布在表面的噪声源同样适用. 在分层介质中, 点源声场的波函数的积分表达式为

      $ p\left( {r,z} \right) = \frac{1}{2}\int \nolimits_{ - \infty }^\infty {{G}}\left( {{z_1},{z_0},{k_r}} \right){\rm{H}}_0^{\left( 1 \right)}\left( {{k_r}r} \right){k_r}{\rm{d}}{k_r}, $

      式中${\rm{H}}_0^{\left( 1 \right)}$为0阶第一类Bessel函数. 该式可以利用积分的闭合围线留数定理求解, 回路积分等于回路所围的极点留数之和以及绕分支割线上下岸的积分. 本文选用Pekeris割线方法进行简单说明, 如图5所示. ${k_1}$, ${k_2}$两点将${k_r}$的正轴分为三个部分, 当$0 < {k_r} < {k_2}$时, ${k_r}$解为复数, 此谱域内声波辐射到海底, 存在声能泄漏, 称为连续谱. 当${k_2} < $ ${k_r} < {k_1} $时, ${k_r}$解为实数, 这部分谱含有离散极点, 相应的模式为无损模式, 称为离散谱.

      图  5  Pekeris分支割线、极点和积分的复波数平面

      Figure 5.  Pekeris branch cut map

      为了描述准确, 这一部分在K/I模型下进行讨论[2]. Kuperman给出海面噪声源均匀情况下噪声互谱密度的波数积分表达式

      $\begin{split} & {C_\omega }\left( {r,{z_1},{z_2}} \right) \\ =\, & \frac{{8{{\text{π}}^2}p{q^2}}}{{{k^{2p}}}}\int \nolimits_0^k {\left( {{k^2} - k_r^2} \right)^{p - 1}} \\ & \times {{g}}\left( {{k_r},{z_1},z'} \right){{{g}}^*}\left( {{k_r},{z_2},z'} \right){{\rm{J}}_0}\left( {{k_r}r} \right){k_r}{\rm{d}}{k_r},\end{split}$

      和简正波表达式

      $ \begin{split} & {C_\omega }\left( {r,{z_1},{z_2}} \right) \\ =\, & \frac{{{\text{π}}{q^2}}}{{2{\rho ^2}{k^2}}}\mathop \sum \nolimits_m \frac{{{{\left[ {{\varPsi _{\rm{m}}}\left( {z'} \right)} \right]}^2}{\varPsi _m}\left( {{z_1}} \right){\varPsi _m}\left( {{z_2}} \right){{\rm{J}}_0}\left( {{\kappa _m}r} \right)}}{{{\alpha _m}{\kappa _m}}},\end{split}$

      式中g为点源格林函数, ${\varPsi _m}$为模式函数. 由于接收器为垂直阵, 式中r = 0. 对于近场连续谱部分采用(5)式波数积分表达式进行计算, 而远场离散谱部分则采用(6)式简正波表达式, 将其分别代入(3)式中求得不同谱域内的$\widehat {BL}$. 图6给出了无限大液态声学半空间海底情况下, 不同谱域内噪声场的垂直指向性和对应的$\widehat {BL}$, 图6(a)图6 (b)是利用(5)式计算近场连续谱部分, 即$0 < {k_r} < {k_2}$, 这一部分的贡献主要来自大掠射角, 其对应的是大于临界角的部分, 而这部分经过海底后存在能量泄漏, 返回到接收器阵的声波可以表达出海底反射的信息; 图6(c)图6(d)是(6)式计算得到的远场离散谱部分, 此时${k_2} < {k_r} < {k_1}$, 其贡献主要来自于小掠射角, 对应小于临界角的部分, 这一角度内发生全内反射, 海底没有声能吸收, 对$\widehat {BL}$的大掠射角条纹部分没有贡献; 图6(e)图6 (f)为全谱域结果, $\widehat {BL}$与仅考虑近场条件下的$\widehat {BL}$基本一致. 在近场条件下大掠射角部分反射损失较高, 所以本文不考虑多次海底反射的情况.

      图  6  不同谱域内噪声场的指向性和$\widehat {BL}$ (a), (b)连续谱部分噪声场; (c), (d) 离散谱部分噪声场; (e), (f)全谱域噪声场

      Figure 6.  Vertical directivity and $\widehat {BL}$ of ocean ambient noise in different spectral domain: (a), (b) Near field continuous spectrum; (c), (d) far field discrete spectrum; (e), (f) full spectrum domain

    • 前文中介绍了利用噪声获得$\widehat {BL}$的有效性, 而有限阵元带来的影响尚不能明确, 因此从理想反射系数出发, 来解释这一过程. 首先考虑图4(a)中的无限大液态声学半空间海底模型.

      影响BL的海底地声参数有3个, 分别是海底衰减系数$\alpha $、声速c和密度$\rho $. 图7别为不同海底参数下的BL. 由图7(a)可见, $\alpha $仅在掠射角小于临界角${\theta _{\rm{c}}}$时, 对BL有影响, 在建立浅海远程传播模型中起到了至关重要的作用[21]. 但第3节中, 通过对近场连续谱和远场离散谱的分解已经证明: 噪声估计的$\widehat {BL}$主要来自于近场连续谱的贡献, 在近场条件下, 衰减系数$\alpha $又是一个不敏感参数, 因此可以忽略衰减系数$\alpha $带来的影响. 在分界面上掠射角${\theta _1}$, ${\theta _2}$利用Snell定律建立关系

      图  7  不同海底参数下的BL (a)衰减系数$\alpha $; (b)声速c; (c)密度$\rho $

      Figure 7.  The BL under different parameters of sub-bottom: (a) Attenuation coefficient $\alpha $; (b) sound speed c; (c) density $\rho $

      $ {k_1}{\rm{cos}}{\theta _1} = {k_2}{\rm{cos}}{\theta _2}, $

      透射波水平射出, 即${\theta _2} = 0^\circ $时存在临界角效应, 此时, ${\rm{cos}}{\theta _1} ={{{c_1}}}/{{{c_2}}}$, 在海水声速${c_1}$已知的情况下, 反射系数的临界角${\theta _1}$完全由${c_2}$决定. 除此之外, 海底声速对于反射系数的幅度也有一定影响. 图7(b)为不同海底声速下的BL. 海底密度$\rho $仅对BL的幅值起到影响, 与临界角不存在任何关系, 如图7(c).

    • 下面考虑存在一层沉积层时的反射情况, 反射系数的表达式为

      $ R = \frac{{{R_{12}} + {R_{23}}{\rm{exp}}\left( {2{\rm{j}}\varphi } \right)}}{{1 + {R_{12}}{R_{23}}{\rm{exp}}\left( {2{\rm{j}}\varphi } \right)}}, $

      其中, ${R_{12}},{R_{23}}$分别为两个分界面上的反射系数, $\varphi = {k_2}h{\rm{sin}}{\theta _2}$表示声波穿过沉积层的垂直相移. (8)式的分子和分母中均含有${\rm{exp}}\left( {2{\rm{j}}\varphi } \right)$一项, 难以讨论其各个参数对R的影响, 特别是图2反射系数中影响条纹的因素, 为了让各个参数在两式中的作用明显表达, 需要对其进行一定程度上的化简, 下面从单一沉积层的反射模型开始考虑, 如图8.

      图  8  单层沉积层海底反射模型

      Figure 8.  Reflection model of sub-bottom with 1 layer of sediment

      声波在沉积层中传播, 会形成多次返回海水中的反射波, 反射系数的原始表达为各个反射波的叠加

      $ \begin{split} R = \, &{R_{12}} + {T_{12}}{R_{23}}{T_{21}}\exp \left( {2{\rm{j}}\varphi } \right) \\ &+ {T_{12}}R_{23}^2{R_{21}}{T_{21}}\exp \left( {4{\rm j}\varphi } \right) + \cdots , \end{split}$

      式中, ${T_{12}}$表示海水进入沉积层的透射系数, ${T_{21}}$代表沉积层进入海水中的透射系数. 近场条件中, 海面风成噪声源源级较低, 可忽略海底多次反射后返回海水中反射波带来的影响, 仅考虑图8中前两条返回声线, 这种情况下反射系数化简为

      $ R = {R_{12}} + {T_{12}}{R_{23}}{T_{21}}\exp \left( {2{\rm{j}}\varphi } \right). $

      与(8)式相比, (10)式非常简洁, 便于讨论不同参数与R的对应关系. 图9分别为(8)和(10)两式计算的海底反射损失BL和1800 Hz下的BL曲线, 依然沿用表1中参数. 从图9可以看出, 两式计算的结果基本一致, 仅在幅值上存在细微差异. 从级数的角度解释, 大掠射角近场条件下, R为一纯小数, (9)式中第三项之后均为第二项的高阶小量, 因此仅在幅值上存在一定影响. 综上, 可以利用(10)式化简后的反射系数R代替(8)式, 对单层沉积层海底反射进行讨论.

      图  9  (a) (8)式计算的真实BL; (b) (10)式计算化简后的BL; (c) 1800 Hz下两条BL曲线

      Figure 9.  (a) True BL computed by equation (8); (b) BL after simplified computed by equation (10); (c) the curves of BL of two methods under 1800 Hz

      ${\theta _{{\rm{c}}1}}$${\theta _{{\rm{c}}2}}$表示两个临界角, 当掠射角$\theta < {\theta _{{\rm{c}}1}}$时, 没有透射波, (10)式仅存在${R_{12}}$一项, 此时R为模近似1的复数; 当${\theta _{{\rm{c}}1}} < \theta < {\theta _{{\rm{c}}2}}$时, (10)式中${R_{12}},\;{T_{12}}$, ${T_{21}}$均为实数, ${R_{23}}$为模近似1的复数, ${R_{12}}$为随$\theta $单调递增的负数, 其模单调递减, ${T_{12}} = $$1 + {R_{12}} $为正且单调递增, ${T_{21}}$同理, 因此对于此部分振荡完全由$\exp \left( {2{\rm{j}}\varphi } \right)$一项引起; 当$\theta > {\theta _{{\rm{c}}2}}$时, (10)式中除$\exp \left( {2{\rm{j}}\varphi } \right)$一项以外均为随$\theta $单调改变的实数, 其本身不存在极值点, 因此这一部分与${\theta _{{\rm{c}}1}} < \theta < {\theta _{{\rm{c}}2}}$时一样, BL曲线中的振荡均由$\exp \left( {2{\rm{j}}\varphi } \right)$一项引起. 在计算BL过程中, 取R的模, 所以$\exp \left( {2{\rm{j}}\varphi } \right)$项中仅有实部对BL曲线有贡献, 其实部通过Euler公式可知

      $ \cos \left( {2\varphi } \right) = \cos \left( {2{k_2}h{\rm{sin}}{\theta _2}} \right), $

      因此BL曲线中的振荡由声速和沉积层厚度h共同决定. 在海底声速确定的情况下, 对于某一掠射角${\theta _0}$, BL曲线对于频率f的振荡周期可由(11)式确定, 把频率f看作变量, (11)式表示为

      $ {\rm{cos}}\left( {2{\text{π}} \cdot \frac{{2h{\rm{sin}}{\theta _2}}}{{{c_2}}} \cdot f} \right), $

      其振荡周期为${{{c_2}}}/{{2h{\rm{sin}}{\theta _2}}}$, 在某一频率${f_{\rm{c}}}$范围内振荡次数为${{2h{\rm{sin}}{\theta _2}}}/{{{c_2}}} \cdot {f_{\rm{c}}}$, 振荡次数即BL的极值点数, 这与图9(a)图9(b)BL的条纹个数是对应的.

      根据互易原理, 可以看作不同距离的表面噪声源到达同一接收点的路径, 图8中前两条反射声线看作入射声线, 原来的入射声线看作汇聚后的出射声线. 两条入射声线的相位差等于声线在沉积层中走过路线的相位变化, 即垂直相位变化$2\varphi $. 在相位相差$2{\text{π}} $的整数倍时发生相长干涉,

      $ \begin{split} & 2\varphi = 2{k_2}h = 2{\text{π}}n,\\ & {f_0} = \frac{{{c_2}n}}{{2h{\rm{sin}}{\theta _2}}},~~n = 1,2,3, \cdots , \end{split} $

      (13)式${f_0}$即为(12)式中f的峰值点, 与其振荡周期相对应, 同样也解释了图9BL条纹的形成机理.

      而在某一频率确定的情况下, (12)式中${\rm{sin}}{\theta _2}$可用(7)式确定, 把${\rm{sin}}{\theta _2}$作为变量, (12)式的振荡周期为${{{c_2}}}/{{(2fh)}}$, 从图10可以看出BL在大掠射角范围内关于${\rm{sin}}{\theta _2}$均匀振荡.

      图  10  (10)式计算化简后的BL (横轴利用(7)式转为${\rm{sin}}{\theta _2}$)

      Figure 10.  Curve of BL after simplified computed by equation (10) with the horizontal axis converted to ${\rm{sin}}{\theta _2}$ by equation (7)

    • BL曲线的临界角效应和周期振荡这两种特征可以分别用来确定单沉积层海底的分层声速和沉积层深度. 在利用噪声获得$\widehat {BL}$时, 由于阵元的个数限制, 会引起$\widehat {BL}$条纹的模糊, 图11(a)仿真了42阵元, 第一阵元设置在30 m处的$\widehat {BL}$, 其他参数依旧沿用表1, 与图2阵元遍布全海深的$\widehat {BL}$相比, 条纹显得模糊很多. 图11(b)单独提取出了1800 Hz下的$\widehat {BL}$曲线, 该曲线中依然可以看到较为明显的临界角信息, 即幅值发生跃变的点, 通过两处临界角${\theta _{{\rm{c}}1}},{\theta _{{\rm{c}}2}}$可以利用上文所述的Snell定律, 在折射角为0时确定沉积层声速

      图  11  (a) 42阵元设置在水下30 m处$\widehat {BL}$; (b) 1800 Hz下$\widehat {BL}$曲线

      Figure 11.  (a) $\widehat {BL}$ with 42 elements at the depth of 30 m; (b) the curve of $\widehat {BL}$ under 1800 Hz.

      $ {c_2} = \frac{{{c_1}}}{{{\rm{cos}}{\theta _{{\rm{c}}1}}}} $

      和基底声速

      $ {c_3} = \frac{{{c_1}}}{{{\rm{cos}}{\theta _{{\rm{c}}2}}}}. $

      而在${\theta _{{\rm{c}}2}}$后, 存在一个明显“凸点”, 这个“凸点”不一定是这条曲线本身的极值点, 但实际上对应的是理想情况下反射损失BL的一个极值点, 由于有限阵元带来的模糊, 使得这个极值点变为了“凸点”. 假设“凸点”的位置对应为${\theta _{\rm{t}}}$, 那么利用${\theta _{\rm{t}}}$${\theta _{{\rm{c}}2}}$便可以得到(12)式在某一固定频率下关于$\sin {\theta _2}$的半周期

      $\begin{split} & \sin \left[ {{\rm{arcos}}\left( {\frac{{{c_2}}}{{{c_1}}}{\rm{cos}}{\theta _{\rm{t}}}} \right)} \right] - \sin \left[ {{\rm{arcos}}\left( {\frac{{{c_2}}}{{{c_1}}}{\rm{cos}}{\theta _{{\rm{c}}2}}} \right)} \right] \\ =\, & \frac{1}{2} \cdot \frac{{{c_2}}}{{2fh}},\\[-15pt]\end{split}$

      从而可以估计出沉积层厚度$\hat h$,

      $ \hat h \!=\! \frac{{{c_2}}}{{4f \!\left\{\!{{\rm{sin}}\left[\!{{\rm{arcos}}\left(\!{\dfrac{{{c_2}}}{{{c_1}}}{\rm{cos}}{\theta _{\rm{t}}}}\!\right)}\!\right]\! - \!{\rm{sin}}\!\left[\!{{\rm{arcos}}\left(\!{\dfrac{{{c_2}}}{{{c_1}}}{\rm{cos}}{\theta _{{\rm{c}}2}}}\!\right)}\!\right]}\!\right\}}}. $

      图11(b)可以看出, ${\theta _{\rm{t}}}$约为46°, ${\theta _{{\rm{c}}2}}$约为41°, 利用(17)式估计出沉积层厚度$\hat h$约为2.84 m, 与实际仿真设置的3 m相比偏小, 这是因为有限阵元带来的模糊, 即空间采样点不足导致$\widehat {BL}$曲线中的极值点之间的分辨率下降. 也可以利用理论阵元合成的方法对噪声的互谱密度矩阵进行扩展, 以提高$\widehat {BL}$曲线的分辨率[17,18].

    • 实验数据于2016年12月在中国黄海某海区采集, 当地实际水深35 m, 冬季海水近似看作等声速, 声速为1487 m/s, 实测声速剖面如图12所示. 实验采用阵元间隔为0.2 m的25元船挂阵采集噪声数据, 将实测数据分为20组, 每组时长5 min, 对每组数据分别进行处理. 图13图14分别给出了第9组和第11组数据处理结果, 1600 Hz下两图的$\widehat {BL}$曲线中明显存在两次突变, 因此存在两个临界角, 符合单层沉积层海底模型, 可以利用上文所述的方法对海底分层结构及地声参数进行估计. 图中已标注出${\theta _{{\rm{c}}1}},{\theta _{{\rm{c}}2}},{\theta _{\rm{t}}}$的位置, ${\theta _{{\rm{c}}1}}$, ${\theta _{{\rm{c}}2}}$的选取依据主要是根据$\widehat {BL}$曲线中幅值突变的程度, 选取其突变的初始角度, ${\theta _{\rm{t}}}$选取的是${\theta _{{\rm{c}}2}}$后的第一个“凸点”, 然后利用(14), (15), (17)三式对该海区海底声速和沉积层厚度进行估计. 图15为20组数据计算结果的统计分布直方图, 其中图15(a)为沉积层声速, 图15(b)为基底声速, 图15(c)为沉积层厚度. 沉积层声速的估计平均值为1544 m/s, 这与同海区的反演结果1531 m/s相近[22]; 基底声速的估计平均值为1935 m/s, 相比沉积层声速其波动范围略大, 主要集中在1900—1960 m/s之间, 目前暂无其他方法对这一结果进行对比; 沉积层厚度的估计平均值为2 m, 考虑到有限阵元带来的模糊影响, 实际沉积层厚度可能略高于估计值.

      图  12  实验实测海水声速剖面

      Figure 12.  Sound speed profile of the experiment

      图  13  第9组数据处理结果 (a) $\widehat {BL}$; (b) 1600 Hz下的$\widehat {BL}$曲线

      Figure 13.  Results of the 9th set of data: (a) $\widehat {BL}$; (b) curve of $\widehat {BL}$ under 1600 Hz

      图  14  第11组数据处理结果 (a) $\widehat {BL}$; (b) 1600 Hz下$\widehat {BL}$曲线

      Figure 14.  Results of the 11th set of data: (a) $\widehat {BL}$; (b) curve of $\widehat {BL}$ under 1600 Hz

      图  15  数据计算结果统计分布直方图 (a)沉积层声速; (b)基底声速; (c)沉积层厚度

      Figure 15.  Statistical distribution histogram of 20 sets of data: (a) Sound speed of the sediment; (b) sound speed of the basement; (c) thickness of the sediment

      为验证所估计沉积层厚度的可靠性, 试验船搭载了浅地层剖面仪(sub-bottom profiler), 该地层剖面仪发射调频信号, 利用声波在水中和水下沉积层内传播和反射的特性来探测海底地层结构. 图16给出了从试验点向东北方向扫描的海底地形, 35 m深处为海底表层, 在其下方存在沉积层, 图16(b)中用黑色标出, 厚度约为2.2 m, 噪声估计出的沉积层厚度与此结果基本符合.

      图  16  浅地层剖面仪探测试验海区海底结构(图(b)中用黑点标明海底分界面)

      Figure 16.  Sub-bottom structure of the experimental area detected by sub-bottom profiler. The interface is marked by black spot in panel (b)

    • 本文根据Harrison能流理论, 假设噪声源均匀分布在海面, 利用垂直阵接收的海洋环境噪声数据对海底反射损失进行估计, 讨论了无限大液体声学半空间海底和存在一层沉积层两种情况的$\widehat {BL}$曲线特点, 并将近场连续谱和远场离散谱两部分的贡献分开考虑, 从理想反射系数出发, 推导了各个海底参数对BL的影响.

      结果表明, 当海底不存在沉积层时, $\widehat {BL}$曲线中仅存在一处幅度突变点, 即一个临界角, 而存在一层沉积层时, 由于有两个分界面, 所以产生了两个临界角; 远场离散谱部分对$\widehat {BL}$几乎没有任何贡献, 主要由近场连续谱部分反映处海底信息; 不同参数对理想反射系数影响的研究表明, 沉积层声速已知的情况下, BL曲线中的振荡主要由沉积层厚度h决定, 利用互易原理两条入射声线会发生相长干涉, 也同样解释了BL条纹的形成机理. 利用$\widehat {BL}$曲线中的临界角信息可以分别对沉积层声速、基底声速进行估计, 而利用“凸点”${\theta _{\rm{t}}}$与第二个临界角${\theta _{{\rm{c}}2}}$可以确定$\widehat {BL}$曲线在某一频率下关于${\rm{sin}}{\theta _2}$的振荡的周期, 从而估计出沉积层的深度, 由于有限阵元个数带来的模糊影响, 这一估计出的结果可能会略小于实际值. 实测海洋环境噪声数据处理结果表明, 利用该方法估计得到的沉积层声速与海底沉积层厚度的结果可信度较高.

参考文献 (22)

目录

    /

    返回文章
    返回