搜索

x

留言板

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

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

浅海小掠射角的海底界面声反向散射模型的简化

侯倩男 吴金荣

侯倩男, 吴金荣. 浅海小掠射角的海底界面声反向散射模型的简化. 物理学报, 2019, 68(4): 044301. doi: 10.7498/aps.68.20181475
引用本文: 侯倩男, 吴金荣. 浅海小掠射角的海底界面声反向散射模型的简化. 物理学报, 2019, 68(4): 044301. doi: 10.7498/aps.68.20181475
Hou Qian-Nan, Wu Jin-Rong. Simplification of roughness bottom backscattering model at small grazing angle in shallow-water. Acta Phys. Sin., 2019, 68(4): 044301. doi: 10.7498/aps.68.20181475
Citation: Hou Qian-Nan, Wu Jin-Rong. Simplification of roughness bottom backscattering model at small grazing angle in shallow-water. Acta Phys. Sin., 2019, 68(4): 044301. doi: 10.7498/aps.68.20181475

浅海小掠射角的海底界面声反向散射模型的简化

侯倩男, 吴金荣

Simplification of roughness bottom backscattering model at small grazing angle in shallow-water

Hou Qian-Nan, Wu Jin-Rong
PDF
HTML
导出引用
  • 在浅海, 尤其是负梯度声速剖面和海面较为平静的浅海波导, 海底界面反向散射是浅海混响的主要来源. 经验散射模型只适用于分析浅海混响平均强度衰减特性, 而基于物理机理建立的反向散射模型克服了这一缺陷, 但同时也引入了其受地声模型约束的问题. 本文结合了海底反射系数的三参数模型, 对浅海远场海底反向散射模型进行了简化, 以减少地声模型的输入参数. 理论分析了海底反射系数的相移参数可以描述海底对声场的散射作用, 无需任何海底地声参数的先验知识. 通过对海底反向散射模型近似简化, 结果表明在临界角附近和甚小掠射角范围内的海底粗糙界面反向散射模型的角度特性和强度特性受海底沉积层的影响不同: 在临界角附近, 海底反向散射的角度特性受海底反射系数的相移参数加权,而其散射系数则近似与相移参数无关; 对于甚小掠射角, 海底反向散射的角度特性近似与海底反射系数的相移参数无关, 其散射系数则近似与相移参数的4次方成正比.
    Bottom backscattering due to roughness seafloor is the main source of shallow water reverberation, especially in the waveguide with downward reflection profile or a calm sea-surface. Empirical backscattering models with a simple form has an important limitation to analyzing other characteristics of reverberation except for the intensity characteristics, which originates from optics and describes the relationship between the bottom backscattering strength and scattering grazing angle of plane-wave in half-infinite space. In the shallow water, such a plane-wave backscattering model cannot be used due to frequency dispersion. The model of bottom backscattering based on physical scattering principle is made to relieve such a limitation, but thereby bringing about another restraint by a geoacoustics model. The bottom backscattering model, which is formulated during modeling the full-wave reverberation theory at small grazing angle in range-independent shallow water waveguide, is simplified by combining with bottom reflection coefficient model which is independent of the geoacoustics model. The bottom reflection coefficient model as referred to the proposed phase parameter P in this paper is equivalent to velocity and density of sediment to describe sound field interacted with sea-bottom. Therefore simplification of bottom backscattering model can be handled by the phase parameter without any knowledge of bottom geoacoustic parameters. The angular dependency and intensity dependency of bottom backscattering due to roughness seafloor at small grazing angle are studied more in depth through such a simplified model. Marking 2/P as the cut-off point, the grazing angle is divided into two stages. Near the critical angle, as grazing angle is greater than 2/P and less than critical grazing angle, the angular dependency of bottom backscattering due to roughness seafloor is weighted by phase parameter of bottom reflection coefficient, while the intensity dependency is independent of phase parameter. At each small grazing angle, as grazing angle is less than 2/P, the angular dependency of bottom backscattering due to roughness seafloor is proportional to incident and scattering grazing angle squared and irrespective of phase parameter of bottom reflection coefficient which is like the empirical bottom backscattering model, while the intensity dependency is proportional to the fourth power of phase parameter. So the bottom has different influences on the angular dependency and intensity dependency of bottom backscattering in different stages of grazing angle.
      通信作者: 吴金荣, wujinrong@mail.ioa.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 11374323, 11774375)资助的课题.
      Corresponding author: Wu Jin-Rong, wujinrong@mail.ioa.ac.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11374323, 11774375).

    海底作为浅海混响的主要散射源, 浅海海底界面混响是很多学者探究的课题. 浅海混响根据其形成过程可以分为传播和散射两个过程. 入射声传播和散射声传播过程同属信道声传播问题, 研究较为完善. 而散射过程则比较复杂, 散射强度的测量也较难实现, 尤其是小掠射角的反向散射的测量, 所以针对散射的研究相对比较受限. 在混响研究伊始, 各国学者均采用经验的散射模型描述海底对声场的散射作用. 由于其形式较为简单, 运算速度较快, 在主动声纳预报系统中有很大的优势[1-4], 至今仍被延续使用. 但是该类型的散射模型也存在着很大的局限性, 很难分析海底散射机理. 最初的经验散射模型是借鉴光学理论中提出的半无限自由空间中的Lambert散射定律给出的海底散射模型, 是海底散射强度与平面波掠射角之间的关系. 而在波导环境中, 由于频散效应的存在, 平面波散射理论不再适用于声场处理, 遂采用简正波理论, 避免了频散产生的多途现象. 然而简正模态的掠射角不同于平面波的掠射角, 因此对波导环境中的声散射问题不能够直接采用经验的Lambert散射模型. 近年来, 越来越多的学者开始从海底散射的物理机理着手建立散射模型. 国际上关于物理散射模型的建模方法包括有限元方法[5]、Kirchhoff近似方法[6]、微扰近似方法[7-10]. 国内关于物理散射模型的研究最具代表性的是中国科学院声学研究所的高天赋和尚尔昌. Shang等[11]在Gao[12]和Tang[13]的研究基础上提出了全波动混响理论, 并对该模型进行了一系列的发展. 2001年, Gao等[14]提出了依据混响数据反演海底反向散射矩阵的方法. Wu等[15,16]将海底反射模型引入到浅海混响模型中简化了混响衰减特性与海底沉积层参数之间的关系, 为海底地声参数的反演提供了一种新方法.

    本文在尚尔昌和吴金荣的研究基础上, 对反向散射模型进行进一步分析. 通过引入海底反射系数的相移参数, 描述海底对声场的散射作用, 分析海底反向散射模型的角度特性和强度特性. 在实际应用允许的误差范围内, 进一步分析海底粗糙界面反向散射模型的角度特性和强度特性, 并与经验散射模型进行对比, 说明两者之间的异同.

    在浅海波导环境中, 海底粗糙界面是海底混响的主要散射源. 在图1所示的浅海波导环境中, 平坦海底的平均海深为${H_0}$, 任意水平位置$r$的海底粗糙界面相对于平均海深的起伏高度为$\eta (r)$, 其量值远小于${H_0}$, 均值为0, 方差为$\sigma _\eta ^{\rm{2}}$. 海底粗糙界面对声场的散射强度远小于入射声场的强度, 即满足弱散射的条件. 海底为半无限均匀介质, 其声速、密度和声吸收系数分别为$ {c_b},{\rho _b}\,{\text{和}}\,\alpha _b^{(w)} $(声吸收系数的单位是${{{\rm{dB}}}/{\rm{\lambda }}}$). 海面为平坦的自由边界, 不考虑其对声场的散射作用.

    图 1  浅海波导环境
    Fig. 1.  Shallow water waveguide environment.

    全波动混响理论的基础思想是将海底粗糙界面对声场的散射作用看作是“二次声源”向外辐射声能量, 从而将海底粗糙界面的散射作为“声源”向波导环境中辐射声能量[11]. 该“声源”的声源强度与入射声强度成正比关系. 在水平均匀的波导环境中, 单位简谐点源入射声场满足波动方程为

    ${\nabla ^2}{u_{{\rm{w,b}}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}} + k_{{\rm{w,b}}}^{\rm{2}}{u_{{\rm{w,b}}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}} = 4{\rm{{\text π} \delta (}}{{{R}}_1} - {{{R}}_0}{\rm{)}}, $

    (1)

    海底粗糙界面的“二次声源”满足

    $ \begin{split} & {\left. {\frac{{\partial {u_{\rm{w}}}({ R_1},{ R_0})}}{{\partial z}} - \frac{{\partial {u_{\rm{b}}}({ R_1},{ R_0})}}{{\partial z}}} \right|_{z = {H_0}}}\\ = & V({R_1}){G^{\rm{i}}}({ R_1},{ R_0}), \end{split} $

    (2)

    $ \begin{split} & {\left. {{\rho _{\rm{w}}}{u_{\rm{w}}}({ R_1},{ R_0}) - {\rho _{\rm{b}}}{u_{\rm{b}}}({ R_1},{ R_0})} \right|_{z = {H_0}}}\\ = & p({ R_1}){G^{\rm{i}}}({ R_1},{ R_0}), \end{split} $

    (3)

    ${G^{\rm{i}}}({{{R}}_1},{{{R}}_0}) = \sqrt {\frac{{2{\rm{\pi }}i}}{{{k_{\rm{w}}}r}}} \sum\limits_{m = 1}^M {{\phi _m}({z_0}){\phi _m}{\rm{(}}{H_0}{\rm{)}}{{\rm{e}}^{{\rm{i}}{k_m}r - {\beta _m}r}}} , $

    (4)

    $V({{{R}}_1}) = \left(k_{\rm{w}}^{\rm{2}} - \frac{{k_{\rm{b}}^{\rm{2}}}}{\kappa }\right)\eta {\rm{(}}{r_1}{\rm{)}} + \left(1 - \frac{1}{\kappa }\right){\nabla _ \bot } \cdot {\rm{(}}\eta ({r_1}{\rm{)}}{\nabla _ \bot }{\rm{)}}, $

    (5)

    $p{\rm{(}}{{{R}}_1}{\rm{)}} = {\rm{(}}{\rho _{\rm{w}}} - {\rho _{\rm{b}}}{\rm{)}}\eta {\rm{(}}{r_1}{\rm{)}}\frac{\partial }{{\partial z}}, $

    (6)

    ${k_{\rm{w}}} = \frac{{2{{\text π}}f}}{{{c_{\rm{H}}}}},\;\;\;\;\;\;{k_{\rm{b}}} = \frac{{2{{\text π}}f}}{{{c_{\rm{b}}}}},\;\;\;\;\;\;\kappa = \frac{{{\rho _{\rm{b}}}}}{{{\rho _{\rm{w}}}}}, $

    (7)

    其中, ${{{R}}_0} = {\rm{(}}{r_{\rm{s}}},{z_0}{\rm{)}}$, ${{{R}}_1} = ({r_1},{H_0})$分别是声源和散射点的空间位置, 两者之间的水平距离差为$r = |{r_{\rm{s}}} - {r_1}|$; $ u_{\rm{w}}, u_{\rm{b}}$分别是水介质和沉积层介质中“二次源”的辐射声场(下脚标w和b分别代表水介质和沉积层介质); ${k_{{\rm{w,b}}}}$是水介质和沉积层介质中的波数, 与声源频率$f$和介质声速有关; ${c_{\rm{H}}} = $$ {c_{\rm{w}}}{\rm{(}}{H_0}{\rm{)}}$是海深处的声速, 沉积层介质和水介质的密度比用$\kappa $表示; ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$是波导环境中初级声场的格林函数. 根据简正波理论的思想, 在浅海波导环境中, ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$可以表示为简正模态${\phi _m}({H_0})$(实际上是${\phi _m}(z)$, $z$是接收深度, 在这里的接收深度为海深${H_0}$)叠加的形式, 如(4)式所示, km${\beta _m}$分别是简正模态本征值的实部和虚部. $V({{{R}}_1})$$p({{{R}}_1})$则分别是质点振速算子和压力算子, 作用到初级声场的格林函数${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$, 得到海底粗糙界面的“二次速度源”的声源强度$V{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$和“二次压力源”的声源强度$p{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$. 在瑞利参数较小的情况下, Bass给出了(5)式和(6)式的近似表达式. ${\nabla _ \bot }$为水平梯度算子,

    ${\nabla _ \bot } = \frac{\partial }{{\partial r}}.$

    (8)

    根据格林定理, 水平均匀波导环境中, “二次速度源”在接收点${{R}} = {\rm{(}}{r_{\rm{r}}},z{\rm{)}}$处的辐射声场可以通过格林函数$G_{{\rm{w,b}}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$给出形式解. 辐射声场的格林函数满足

    $\begin{aligned} & {\nabla ^2}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} + k_{\rm{w}}^{\rm{2}}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = 0\;\;\;z < {H_0}\\ & {\nabla ^2}G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} + k_{\rm{b}}^2G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = 0\;\;\;z > {H_0} \end{aligned},$

    (9)

    同时满足边界条件

    $\begin{gathered} {\left. {\frac{{\partial G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}}}{{\partial z}} - \frac{{\partial G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}}}{{\partial z}}} \right|_{z = {H_0}}} = {\rm{\delta (}}{{R}} - {{{R}}_1}{\rm{)}} \\ {\left. {\rho G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} - {\rho _{\rm{b}}}G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}} \right|_{z = {H_0}}} = 0 \\ \end{gathered} .$

    (10)

    在接收点 “二次速度源”的辐射声场满足形式解,

    $u_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = \int {G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}V{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)d}}{{{R}}_1}} .$

    (11)

    “二次压力源”的辐射声场同样有(11)式的形式解,

    $u_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = \int {G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}p{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)d}}{{{R}}_1}} ,$

    (12)

    其中, $G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$是“二次压力源”的辐射声场的格林函数. 根据边界处法向振速与压力之间的关系, 不难得到$G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$$G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$满足

    $G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = - \frac{1}{{{\rho _{\rm{w}}}}}\frac{\partial }{{\partial z}}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}.$

    (13)

    叠加“二次速度源”和“二次压力源”的辐射声场, 得到海底粗糙界面向波导环境中辐射的总声场

    ${u_{\rm{w}}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = u_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} + u_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}}.$

    (14)

    将(11)—(13)式代入(14)式, 得到海底粗糙界面散射声场的一般形式解

    $ \begin{split} & {u_{\rm{w}}}( R,{ R_1},{ R_0})\\ = & \int [G_{\rm{w}}^{\rm{V}}( R,{ R_1})V({ R_1})\\ & - \frac{1}{{{\rho _{\rm{w}}}}}\frac{\partial }{{\partial z}}G_{\rm{w}}^{\rm{V}}( R,{ R_1})p({ R_1})]{G^{\rm{i}}}({ R_1},{ R_0}){\rm{d}}{ R_1}. \end{split} $

    (15)

    实际上, 入射声场的格林函数、散射场的格林函数以及初级声场的格林函数有一致的形式解, 忽略他们之间的差异,

    ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{R}}{\rm{)}} = G_{\rm{w}}^{\rm{V}}{\rm{(}}{{{R}}_1},{{R}}{\rm{)}} = G{\rm{(}}{{{R}}_1},{{R}}{)},$

    (16)

    所以对(15)式进行简化得到

    $ \begin{split} & {u_{\rm{w}}}( R,{ R_1},{ R_0})\\ = & \int_S {[G( R,{ R_1})V({ R_1})} \\ & - \frac{1}{{{\rho _{\rm{w}}}}}\frac{{\partial G( R,{ R_1})}}{{\partial z}}p({ R_1})]G({ R_1},{ R_0}){\rm{d}}{ R_1}. \end{split} $

    (17)

    对于本地混响, 声源和接收均在同一水平位置. 为了方便计算, 通常令其位于过源点的垂直轴线上, 即${r_{\rm{s}}} = {r_{\rm{r}}} = 0$. 将(4)—(6)式代入(17)式, 并通过分部积分法得到简谐点源的海底粗糙界面的反向散射声场, 将其描述为简正模态的叠加形式[11],

    $ \begin{split} & {u_{\rm{w}}}(r,z)\\ = & \displaystyle\frac{{2\pi i}}{{{k_{\rm{w}}}r}}\sum\limits_m^M {\sum\limits_n^M {{\phi _m}({z_0}){\phi _m}({H_0})} } \cdot C_{mn}^\eta K_{mn}^\eta ({k_m},{k_n})\\ & \times {\phi _n}({H_0}){\phi _n}(z){{\rm{e}}^{ - ({\beta _m} + {\beta _n})r}} \end{split} $

    (18)

    其中

    $C_{mn}^\eta = k_{\rm{w}}^{\rm{2}} - \frac{{k_{\rm{b}}^2}}{\kappa } + {\rm{(}}1 - \frac{1}{\kappa }{\rm{)}}{k_m}{k_n} + \frac{{1 - \kappa }}{{{\kappa ^2}}}{\gamma _m}{\gamma _n}, $

    (19)

    $K_{mn}^\eta {\rm{(}}{k_m},k{}_n{\rm{)}} = \int {\eta {\rm{(}}r'{\rm{)}}{{\rm{e}}^{{\rm{i(}}{k_m} + {k_n}{\rm{)}}r'}}{\rm{d}}r'} .$

    (20)

    实际上(18)式是平坦波导环境中海底粗糙界面散射场的稳定解形式. 对于声源脉冲信号s(t), 其功率谱$F\left( \omega \right)$满足Fourie变换

    $F{\rm{(}}\omega {\rm{)}} = \int_{ - \infty }^{ + \infty } {s{\rm{(}}t{\rm{)}}{{\rm{e}}^{ - {\rm{i}}\omega t}}{\rm{d}}t},\tag{21a}$

    (21a)

    $s{\rm{(}}t{\rm{)}} = \frac{1}{{2{\rm{\pi }}}}\int_{ - \infty }^{ + \infty } {F{\rm{(}}\omega {\rm{)}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}\omega }.\tag{21b}$

    (21b)

    根据Fourie变换的性质, 得到入射信号为s(t)的海底散射声场为

    $ u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{(}}r,z{\rm{)}} = \frac{1}{{2{\rm{\pi }}}}\int_{ - \infty }^{ + \infty } {F{\rm{(}}\omega {\rm{)}}{u_{\rm{w}}}{\rm{(}}r,z{\rm{)}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}\omega } .$

    (22)

    将声源信号的功率谱函数$F\left( \omega \right)$和(18)式代入(22)式, 同时取水平波束在中心频率处的两阶泰勒级数展开, 进而得到脉宽为$\tau $声源信号s(t)的反向散射声场,

    $ \begin{split} & u_{\rm{w}}^{s(t)}(r,z)\\ = & \frac{{2\pi i}}{{{k_{\rm{w}}}r}}s(t - {t_0})\sum\limits_m^M {\sum\limits_n^M {{\phi _m}({z_0}){\phi _m}({H_0}) } } \\ & \times C_{mn}^\eta K_{mn}^\eta \cdot {\phi _n}({H_0}){\phi _n}(z){{\rm{e}}^{ - ({\beta _m} + {\beta _n})r}} \end{split} $

    (23)

    其中, t0是信号的传播时间, 与散射环的水平距离满足${t_0} \approx {{2r}/{{c_{\rm{H}}}}}$.

    混响是能够同一时刻被接收、来自各方向散射回波的叠加. 对于海底粗糙界面混响而言, 同一时刻接收到的散射回波来自于宽度为$\Delta r = {{{c_{\rm{H}}}\tau }/2}$的散射圆环, 其内径为${r_{{\rm{in}}}} = r - \displaystyle\frac{{{c_{\rm{H}}}\tau }}{4}$, 外径为${r_{{\rm{ex}}}} = r + $$ \displaystyle\frac{{{c_{\rm{H}}}\tau }}{4}$. 对(23)式计算强度, 并在散射环面积上进行积分得到海底粗糙界面混响平均强度

    $ {I_\eta }{\rm{(}}r,z{\rm{)}} = \int_A {\langle u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^{*}\rangle } {\rm{d}}A, $

    (24)

    其中, $ {\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^* $$u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}$的复共轭.

    假设散射环面积内的海底界面起伏高度满足各向同性的原则, 对(24)式进一步简化为

    $ {I_\eta }{\rm{(}}r,z{\rm{)}} = 2{\rm{\pi }}r\int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\langle u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^{*}\rangle {\rm{d}}r} . $

    (25)

    将(23)式代入(25)式, 并假设脉宽远小于混响时间t, 将声源强度近似为

    ${E_0} = {s^2}\left( {t - {t_0}} \right) \cdot \tau .$

    (26)

    在这里, 只考虑其非相干特性, 得到短脉冲的非相干混响平均强度

    $ \begin{split} {I_\eta }(r,z) = & {E_0} \cdot {\left(\frac{{2{\text{π}}}}{{{k_{\rm{w}}}r}}\right)^2}\\ & \times {\text{π}} rc\sum\limits_m^M {\sum\limits_n^M {\phi _m^2({z_0})\phi _m^2({H_0}) \cdot S_{mn}^\eta } } \\ & \times \phi _n^2({H_0})\phi _n^2(z){{\rm{e}}^{ - 2({\beta _m} + {\beta _n})r}}, \end{split} $

    (27)

    $S_{mn}^\eta = {\left[ {C_{mn}^\eta } \right]^2}{S^\eta }\left( {2{k_0}} \right),$

    (28)

    $ \begin{split} & {S^\eta }(2{k_0})\\ = & \int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\langle \eta (r')\eta (r'')\rangle {{\rm{e}}^{{\rm{i}}({k_m} + {k_n})r' - {\rm{i}}({k_m} + {k_n})r''}}} } \\ & \times {\rm{d}}r''{\rm{d}}r'\; = \;\sigma _\eta ^{\rm{2}}{P^\eta }(2{k_0}), \end{split} $

    (29)

    其中, ${S^\eta }(2{k_0})$是海底界面的粗糙度谱, 与海底沉积层的地声参数无关, 只是声源频率的函数; $S_{mn}^\eta $是海底粗糙界面引起的入射模态和散射模态之间的耦合系数, 反映海底粗糙界面的反向散射能力, 由海底界面粗糙度谱和地声参数决定. 定义

    $\Theta _{mn}^\eta = \phi _m^2({H_0})S_{mn}^\eta \phi _m^2({H_0})$

    (30)

    是区别于经验散射模型的物理散射模型. 它结合了波导环境特性和简正波理论的思想, 建立了受格林函数约束的散射模型, 明确了海底地声参数以及海底粗糙界面与海底反向散射的定量关系. 因此, 地声参数的准确度直接影响到海底反向散射函数的精确度. 通常情况下, 很难直接获取大面积的地声数据, 通过反演得到的地声参数又受地声模型的影响较大, 而且反演的未知参数较多, 导致结果存在很大的不确定性. 从反射系数的角度思考, 海底对声场的影响主要体现在海底反射系数. 海底反射系数的三参数模型与海底地声模型无关, 且参数较少, 利用海底反射系数代替地声参数作为输入参数, 能很大程度上简化海底粗糙界面的反向散射模型.

    尚尔昌在1979年提出了海底反射系数的三参数模型[17], 随后给出了海底小掠射角的反射系数幅值和相移与掠射角之间的量化关系, 该参数与地声模型无关. 海底反射系数的幅值用Q参数表示

    $\ln \left| {V\left( \theta \right)} \right| = - Q\theta ,\;\;\theta \to 0,$

    (31)

    海底反射系数的相移用P参数表示

    $\arg \left[ {V\left( \theta \right)} \right] = -{\text{π}} + P/\theta ,$

    (32)

    $\theta $是海底界面处入射声波的掠射角; Q参数反映海底声吸收效果, 体现声能量的衰减; P参数则描述各阶简正模态在海深处的能量, 即控制各阶简正模态的水平波数. 海底反射系数的相移参数P与地声参数之间的转换关系满足[18,19]

    $ {P^{{\rm{Low}}}} = \frac{{2\kappa }}{{\sqrt {1 - {\upsilon ^{ - 2}}} }},\;\;{\rm{ }}\theta < \theta ^{\rm{*}}, $

    (33)

    $ {{{P^{{\rm{Crit}}}} = {\text{π}}}/{\theta ^{*}}},\;\;{\rm{ }}\theta \to \theta ^{\rm{*}}, $

    (34)

    其中, $\upsilon = {{{c_{\rm{b}}}}/{{c_{\rm{H}}}}}$表示沉积层介质和水介质的声速比; $\kappa $是密度比, 如(7)式所示; $ \theta ^{*}$是临界掠射角.

    (33)式和(34)式分别给出了不同海底掠射角时海底反射系数的相移参数的不同形式. ${P^{{\rm{Low}}}}$是甚小掠射角处的海底反射系数相移参数, 是海底沉积层和水介质的声速比和密度比的函数; ${P^{Crit}}$是临界角附近海底反射系数的相移参数, 是临界角的函数.

    海底反射系数的幅值参数Q与地声参数之间的转换关系满足[18,19]

    $Q = 0.036 \cdot \frac{{{\upsilon ^{ - 1}}\kappa }}{{{{\left[ {1 - {\upsilon ^{ - 2}}} \right]}^{3/2}}}} \cdot \alpha _b^{\left( w \right)},$

    (35)

    其中, $\alpha _b^{\left( w \right)}$是海底沉积层的声吸收系数, 单位是$dB/\lambda $.

    根据Snell折射定律, $\upsilon $与临界掠射角之间满足

    ${\rm{cos(}}\theta ^{*}{\rm{)}} = {{{c_{\rm{H}}}}/{{c_{\rm{b}}}}} = {\upsilon ^{ - 1}}.$

    (36)

    将(36)式代入(34)式, 临界角附近的海底反射系数的相移参数${P^{Crit}}$与声速之间满足

    $\upsilon = {\rm{co}}{{\rm{s}}^{ - 1}}{\rm{(}}\frac{{\rm{\pi }}}{{{P^{{\rm{Crit}}}}}}{\rm{)}}.$

    (37)

    联合(37)式和(33)式得到甚小掠射角的海底反射系数的相移参数${P^{Low}}$与密度之间满足关系式

    $\kappa = \frac{{{P^{{\rm{Low}}}}}}{2}\sin \left(\frac{{\text{π}}}{{{P^{{\rm{Crit}}}}}}\right).$

    (38)

    (37)式和(38)式分别是海底反射系数的参数与等效的均匀半无限海底介质的地声参数之间的转换关系.

    常用的经验散射模型如(39)式所示:

    ${\Theta _E} = {\mu _E}{\sin ^l}\theta {\sin ^k}\varphi ,$

    (39)

    其中, ${\mu _{\rm{E}}}$是经验散射系数, 反映反向散射强度; $\theta $是入射掠射角; $\varphi $是散射掠射角. 当$l = k = 1$, 是经验的Lambert散射模型. 该类型的散射模型是可分离的经验海底反向散射模型, 不可分离的经验散射模型(参见附录)也被用来描述海底反向散射.

    (30)式的海底反向散射模型是根据波导环境建立的散射模型, 明确了入射简正模态、散射简正模态以及海底地声参数对反向散射的强度特性和角度特性的影响. 浅海波导环境中简正模态是受格林函数严格约束的稳定解, 与海底沉积层的介质参数和水文环境有关. 在掠射角小于临界角的条件下, 忽略各阶简正模态的水平波数(掠射角)之间的差异, 即${k_m} \approx {k_n} \approx {k_{\rm{w}}}$, 则(19)式可以近似描述为与模态无关的量.

    $ \begin{split} C_{mn}^\eta = & \left[1 \!- \!\frac{{k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^2\kappa }} \!+\! \left(1 - \frac{1}{\kappa }\right)\frac{{{k_m}{k_n}}}{{k_{\rm{w}}^2}} + \frac{{1 - \kappa }}{{{\kappa ^2}}}\frac{{{\gamma _m}{\gamma _n}}}{{k_{\rm{w}}^2}}\right]k_{\rm{w}}^2 \\ & \approx \;\left[1 - \frac{{k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^{\rm{2}}\kappa }} + \left(1 - \frac{1}{\kappa }\right) + \frac{{1 - \kappa }}{{{\kappa ^2}}}\frac{{k_{\rm{w}}^2 - k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^{\rm{2}}}}\right]k_{\rm{w}}^{\rm{2}} \\ & = \;\left[2 - 2{\kappa ^{ - 1}} + {\kappa ^{ - 2}} - {\left(\kappa \upsilon \right)^{ - 2}}\right]k_{\rm{w}}^2. \end{split} $

    (40)

    将(37)式和(38)式代入(40)式, 用海底反射系数的相移参数代替地声参数, 得到

    $C_{mn}^\eta = k_{\rm{w}}^2 \cdot \varsigma {\rm{(}}P{\rm{)}},$

    (41)

    其中

    $\varsigma {\rm{(}}P{\rm{)}} \approx 2 - \frac{4}{{P{\rm{sin(}}{{\text{π}}/P}{\rm{)}}}} + \frac{4}{{{P^2}}}.$

    (42)

    在这里忽略${P^{Low}}$${P^{Crit}}$之间的差异, 文献[18]中也说明了这种近似的合理性.

    同时, 海底反射系数的相移参数$P$决定了各阶简正模态在海底粗糙界面平均深度处的能量, 即

    ${\phi _m}({H_0}) = \sqrt {\frac{2}{{{H_0}}}} \sin \left(\frac{{P{\theta _m}}}{2}\right).$

    (43)

    将(41)式和(43)式代入(30)式得到远距离(小掠射角)条件下的海底反向散射核函数,

    ${\Theta _P} = \frac{4}{{H_0^2}}S(2{k_0}){[k_w^2 \cdot \varsigma (P)]^2} \cdot {\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right),$

    (44)

    其中$\theta = {\theta _m}$为入射掠射角, $\varphi = {\theta _n}$为散射掠射角.

    采用海底反射系数的相移参数$P$描述海底反射系数, 明确了海底对声散射的物理机理, 不同于经验反向散射模型.

    在远距离条件下, (44)式表明, 海底反向散射的角度特性受海底反射系数的相移参数$P$的影响, 不同于经验散射模型. 在临界角附近,

    $\mathop {\lim }\limits_{\theta \to \theta *} \frac{{P\theta }}{2} = \frac{{\text{π}}}{2}.$

    (45)

    $\displaystyle\frac{{P\theta }}{2}$为宗量, 不能对$\sin \left(\displaystyle\frac{{P\theta }}{2}\right)$进行小角度近似, 所以在临界角附近, (44)式所描述的反向散射模型随掠射角并不满足线性变化的关系.

    ${{P\theta }/2} = 1$为分界点, 当${{P\theta }/2} < 1$时, 为甚小掠射角范围, 小角度近似引起的海底反向散射强度小于3 dB, 在实际应用中可忽略不计; 当${{P\theta }/2} > 1$时, 为临界角附近, 小角度近似使海底反向散射强度出现大于3 dB的误差, 实际应用中不可忽略. 所以, 在随后的分析中, 对反向散射模型的强度特性和角度特性进行分段处理.

    ${{P\theta }/2} > 1$, 即${2/P} < \theta < \theta *$时, (44)式不能进行小角度近似, 其随角度的变化关系受P参数加权.

    $\begin{split}\Theta _P^{Crit} = & \mu _P^{Crit}{\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right) \\ & \propto {\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right),\end{split}$

    (46)

    其中

    $\mu _P^{Crit} = \frac{4}{{H_0^2}}S(2{k_0}){[k_0^2 \cdot \varsigma (P)]^2}.$

    (47)

    ${{P\theta }/2} < 1$, 即$\theta < {2/P}$时, 对(44)式的角度项进行小角度近似,

    $\sin \left(\frac{{P\theta }}{2}\right) \approx \frac{{P\theta }}{2}.$

    (48)

    将(48)式代入(44)式得到

    $\Theta _P^{Low} = \mu _P^{Low} \cdot {\theta ^2}{\varphi ^2} \propto {\theta ^2}{\varphi ^2},$

    (49)

    其中

    $\mu _P^{Low} = \frac{{{P^4}}}{{4H_0^2}}S(2{k_0}){[k_0^2 \cdot \varsigma (P)]^2}.$

    (50)

    对(39)式的经验散射模型进行小角度近似得到

    ${\Theta _E} \approx {\mu _E}{\theta ^l}{\varphi ^k} \propto {\theta ^l}{\varphi ^k}.$

    (51)

    $l = k = 2$时, 比较(51)式和(49)式具有相同的角度关系. 所以当掠射角$\theta < {2/P}$时, 基于物理散射机理的反向散射模型与$l = k = 2$的可分离经验散射模型具有一致的角度特性.

    通过4.1节的分析, 掠射角不同时, 海底反向散射的角度特性不同, 同时也反映了海底反向散射系数的差异.

    (47)式指出, 在临界角附近, 海底反向散射系数由海底反射系数的相移参数P、海底界面粗糙度谱, 海深以及声源信号的频率决定. 忽略临界角附近以及甚小掠射角的P参数之间的差异[18], 则(42)式中等式右边第三项相比前两项为小量, 可以忽略不计. 第二项中采用小角度近似,

    ${\rm{sin}}\frac{{\text{π}}}{P} \approx \frac{{\text{π}}}{P},$

    (52)

    使得(42)式中的第二项近似与P参数无关, 所以$ \varsigma (P) \approx 0.73 $. 图2也表明了$\varsigma (P)$P参数的变化不是很明显. 所以, 在临界角附近的海底反向散射系数近似与海底沉积层介质无关, 只是海底界面粗糙度谱、海深以及声源频率的函数.

    $ \mu _P^{Crit} \approx \frac{{2.11}}{{H_0^2}}S(2{k_0})k_0^4. $

    (53)
    图 2  $\zeta $P参数的变化
    Fig. 2.  $\zeta $ varied with P parameters near critical angle.

    对于甚小掠射角的海底反向散射系数满足(50)式. 结合$\varsigma (P)$近似为常数的结果, (50)式可以近似为与${P^4}$成正比, 即

    $ \mu _P^{Low} \approx {0.13}\frac{{{P^4}}}{{H_0^2}}S(2{k_0})k_0^4 \propto {P^4}.$

    (54)

    图3以第一类海底为例说明了这种近似结果的合理性. 图中给出了甚小掠射角的海底反向散射系数随$P$参数的变化曲线(仿真参数: 海深为50 m, 声源频率为600 Hz, 海底粗糙界面的标准差为$0.1\;{\rm{m}}$, 相关尺度为${\rm{10}}\;{\rm{m}}$, 计算得到Goff-Jordan谱[20]为-32.8 dB)以及对其进行$P$参数的四次方拟合结果, 两者的变化趋势基本一致. 所以, 在甚小掠射角的条件下, 海底反向散射系数与海底反射系数的相移参数$P$密切相关, 即与海底沉积层的声速比和密度比有关, 与声吸收系数无关.

    图 3  海底甚小掠射角反向散射系数对P参数的依赖
    Fig. 3.  Relationship between P and bottom backscattering coefficient at very low grazing angle.

    另外, 从混响平均强度的角度分析, 同样说明了这种近似的合理性. 海底界面粗糙度谱函数采用图3的仿真参数, 海底是第一类均匀半无限介质, 声速为$18{\rm{36}}\;{\rm{m/s}}$, 密度为${\rm{2}}{\rm{.03}}\;{\rm{g/c}}{{\rm{m}}^{\rm{3}}}$, 声吸收系数为$0.88\;{\rm{dB/}}\lambda $, 海深为50 m, 水介质为1500 m/s的等声速剖面的均匀水体. 声源频率为600 Hz, 发射深度为10 m, 接收深度为30 m. 计算水平散射距离从5 km到50 km的混响平均强度衰减曲线. 图4中实线是根据(30)式仿真的非近似混响平均强度衰减曲线; “”是根据(46)式和(53)式仿真的临界角的近似结果, 与非近似的混响平均强度衰减曲线在13 km以前完全重合; 随着水平距离的增加, 两者相差逐渐增大; “$ \bullet $”是根据(49)式和(54)式仿真的甚小掠射角的近似结果, 在20 km以前, 与非近似结果之间的差异随水平距离的增大而减小, 在20 km以后两者基本吻合. 在散射距离较近时, 掠射角相对较大, 临界角附近近似的海底反向散射模型与真实结果更接近; 在散射距离较远时, 掠射角相对较小, 甚小掠射角的近似结果与真实结果更接近.

    图 4  海底反向散射模型的比较
    Fig. 4.  Compare with bottom backscattering model with different approximate.

    在浅海混响平均强度模型中, 经验散射模型在混响特性分析中存在明显的局限性, 而现有的物理散射模型受地声模型的影响较大. 本文以全波动混响理论的物理散射模型为基础, 结合海底反射系数的三参数模型, 将海底反射系数的相移参数等效代替地声参数, 描述海底对声场的散射作用, 简化了海底反向散射模型. 通过理论分析, 明确了海底掠射角以${2/P}$为分界点, 海底反向散射的角度特性和强度特性对海底反射系数的相移参数存在不同的依赖关系. 在掠射角满足$\theta < {2/P}$时, 海底反向散射的角度特性近似描述为与可分离经验散射模型的角度特性一致, 与$P$参数无关; 而其散射系数则近似描述为与P4线性增强. 当海底掠射角满足${2/P} < \theta < \theta ^{*}$时, 海底反向散射的角度特性是受$P$参数加权的, 即受海底地声参数的影响; 而其散射系数则近似为与$P$参数无关. 所以不同掠射角范围, 海底对反向散射声场的强度特性和角度特性的贡献不同. 掠射角较大时海底对反向散射声场的影响主要体现在其角度特性; 掠射角非常小时, 海底的影响主要体现在其强度特性.

    附录A 经验反向散射模型

    经验反向散射模型有多种形式, 主要体现在其角度特性之间的差异. 最常用的海底散射模型是可分离的散射模型,

    $\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\theta {\rm{si}}{{\rm{n}}^k}\varphi, \tag{A1}$

    (A1)

    其中, ${\mu _{\rm{E}}}$是海底反向散射系数; $l,k$取不同的值对应不同的散射模型. 常见的Lambert散射模型则是(A1)式中$l = k = 1$时的散射模型. 该模型最初是Mackenzie[1]在处理深海海底混响模型时由光学的散射原理引入, 给出经验的海底反向散射模型, 随后根据实验数据与理论结果对比得到${\mu _{\rm{E}}}$的分贝值大约在-27 dB左右, 也就是说${\mu _{\rm{E}}} \approx {10^{ - 2.7}}$. 当$l = 1,k = 0$表示散射源强正比于$\mu {\rm{sin}}\theta $的均匀散射; 当$l = 0,k = 0$表示与角度无关的均匀散射; 同时也存在$l = k = 2$的反向散射模型.

    在海底界面大尺度不均匀波导环境中, 小掠射角和基尔霍夫近似的条件下的海底散射模型可以表示为

    $\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\left(\frac{{\theta + \varphi }}{2}\right) \;\;l = 0,\frac{1}{2},1,\frac{3}{2},....\tag{A2}$

    (A2)

    在海底界面小尺度不均匀波导环境中, 小掠射角和基尔霍夫近似的条件下的海底散射模型可以表示为

    $\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\left[{\rm cos}^{- 1}\left(\frac{{{\rm{cos}}\theta + {\rm{cos}}\varphi }}{2}\right)\right]\;\;l = 0,\frac{1}{2},1,\frac{3}{2},....\tag{A3}$

    (A3)

    另外一种不可分离的散射模型, 同样被用来描述海底反向散射过程

    $\Theta = {\mu _{\rm{E}}}\frac{{{\rm{sin}}\theta {\rm{sin}}\varphi }}{{{\rm{sin}}\theta + {\rm{sin}}\varphi }}.\tag{A4}$

    (A4)

    通常无论是可分离的经验反向散射强度还是不可分离的经验反向散射强度, 它们的共同点在于海底反向散射强度特性与角度特性相互分离, 相互独立.

    [1]

    张仁和, 李文华, 裘辛方, 金国亮 1995 声学学报 20 417

    Zhang R H, Li W H, Qiu X F, Jin G L 1995 Acta Acoust. 20 417

    [2]

    刘建军, 李风华, 张仁和 2006 声学学报 31 173Google Scholar

    Liu J J, Li F H, Zhang R H 2006 Acta Acoust. 31 173Google Scholar

    [3]

    姚万军, 蔡志明, 卫红凯 2009 声学学报 34 223Google Scholar

    Yao W J, Cai Z M, Wei H K 2009 Acta Acoust. 34 223Google Scholar

    [4]

    Zhou J X, Zhang X Z 2012 J. Acoust. Soc. Am. 131 2611Google Scholar

    [5]

    Isakson M J, Chotiros N P 2011 J. Acoust. Soc. Am. 129 1237

    [6]

    彭朝晖, 周纪浔, 张仁和 2004 中国科学G辑 物理学 力学 天文学 34 378

    Peng Z H, Zhou J X, Zhang R H 2004 Sci. China (Series G) 34 378

    [7]

    Ivakin A N 1998 J. Acoust. Soc. Am. 103 827Google Scholar

    [8]

    Ivakin A N 2016 J. Acoust. Soc. Am. 140 657Google Scholar

    [9]

    Moe J E, Jackson D R 1994 J. Acoust. Soc. Am. 96 1748Google Scholar

    [10]

    Tang D J, Jackson D R 2017 J. Acoust. Soc. Am. 142 2968Google Scholar

    [11]

    Shang E C, Gao T F, Wu J R 2008 IEEE J. Ocean Eng. 33 451Google Scholar

    [12]

    高天赋 1989 声学学报 14 126

    Gao T F 1989 Acta Acoust. 14 126

    [13]

    Tang D J 1997 Internal Conference on Shallow-Water Acoustics Beijing, China, April 21-25, 1997 p323

    [14]

    Gao T F, Shang E C, Tang D J 2001 Theoretical and Computational Acoustics Beijing, China, May 21-25, 2001 p67

    [15]

    Wu J R, Shang E C, Gao T F 2010 AIP Conf. Proc. 1272 314

    [16]

    Wu J R, Shang E C, Gao T F 2010 J. Comput. Acous. 18 209Google Scholar

    [17]

    尚尔昌 1979 海洋学报 1 58

    Shang E C 1979 Acta Ocean. Sin. 1 58

    [18]

    Zhao Z D, Ma L, Shang E C 2014 J. Comput. Acoust. 22 1440005Google Scholar

    [19]

    赵振东 2014 博士学位论文 (北京: 中国科学院大学)

    Zhao Z D 2007 Ph. D. Dissertation (Beijing: The University of Chinese Academy of Sciences) (in Chinese)

    [20]

    Goff J A, Jordan T H 1988 J. Geophys. Res. 93 13589Google Scholar

    期刊类型引用(3)

    1. 汪磊,黄益旺,郭霖,任超. 浅海粗糙海底声散射建模及声场特性. 物理学报. 2024(03): 91-98 . 百度学术
    2. 郑毅,秦志亮,于盛齐,马本俊,赵吉祥,刘梦婷. 海底声散射特性进展及展望. 哈尔滨工程大学学报. 2024(04): 748-757 . 百度学术
    3. 马国栋. 基于改进侧扫声呐法的水下抛石精准识别研究. 人民长江. 2022(04): 210-214 . 百度学术

    其他类型引用(0)

  • 图 1  浅海波导环境

    Fig. 1.  Shallow water waveguide environment.

    图 2  $\zeta $P参数的变化

    Fig. 2.  $\zeta $ varied with P parameters near critical angle.

    图 3  海底甚小掠射角反向散射系数对P参数的依赖

    Fig. 3.  Relationship between P and bottom backscattering coefficient at very low grazing angle.

    图 4  海底反向散射模型的比较

    Fig. 4.  Compare with bottom backscattering model with different approximate.

  • [1]

    张仁和, 李文华, 裘辛方, 金国亮 1995 声学学报 20 417

    Zhang R H, Li W H, Qiu X F, Jin G L 1995 Acta Acoust. 20 417

    [2]

    刘建军, 李风华, 张仁和 2006 声学学报 31 173Google Scholar

    Liu J J, Li F H, Zhang R H 2006 Acta Acoust. 31 173Google Scholar

    [3]

    姚万军, 蔡志明, 卫红凯 2009 声学学报 34 223Google Scholar

    Yao W J, Cai Z M, Wei H K 2009 Acta Acoust. 34 223Google Scholar

    [4]

    Zhou J X, Zhang X Z 2012 J. Acoust. Soc. Am. 131 2611Google Scholar

    [5]

    Isakson M J, Chotiros N P 2011 J. Acoust. Soc. Am. 129 1237

    [6]

    彭朝晖, 周纪浔, 张仁和 2004 中国科学G辑 物理学 力学 天文学 34 378

    Peng Z H, Zhou J X, Zhang R H 2004 Sci. China (Series G) 34 378

    [7]

    Ivakin A N 1998 J. Acoust. Soc. Am. 103 827Google Scholar

    [8]

    Ivakin A N 2016 J. Acoust. Soc. Am. 140 657Google Scholar

    [9]

    Moe J E, Jackson D R 1994 J. Acoust. Soc. Am. 96 1748Google Scholar

    [10]

    Tang D J, Jackson D R 2017 J. Acoust. Soc. Am. 142 2968Google Scholar

    [11]

    Shang E C, Gao T F, Wu J R 2008 IEEE J. Ocean Eng. 33 451Google Scholar

    [12]

    高天赋 1989 声学学报 14 126

    Gao T F 1989 Acta Acoust. 14 126

    [13]

    Tang D J 1997 Internal Conference on Shallow-Water Acoustics Beijing, China, April 21-25, 1997 p323

    [14]

    Gao T F, Shang E C, Tang D J 2001 Theoretical and Computational Acoustics Beijing, China, May 21-25, 2001 p67

    [15]

    Wu J R, Shang E C, Gao T F 2010 AIP Conf. Proc. 1272 314

    [16]

    Wu J R, Shang E C, Gao T F 2010 J. Comput. Acous. 18 209Google Scholar

    [17]

    尚尔昌 1979 海洋学报 1 58

    Shang E C 1979 Acta Ocean. Sin. 1 58

    [18]

    Zhao Z D, Ma L, Shang E C 2014 J. Comput. Acoust. 22 1440005Google Scholar

    [19]

    赵振东 2014 博士学位论文 (北京: 中国科学院大学)

    Zhao Z D 2007 Ph. D. Dissertation (Beijing: The University of Chinese Academy of Sciences) (in Chinese)

    [20]

    Goff J A, Jordan T H 1988 J. Geophys. Res. 93 13589Google Scholar

  • [1] 汪磊, 黄益旺, 郭霖, 任超. 浅海粗糙海底声散射建模及声场特性. 物理学报, 2024, 73(3): 034301. doi: 10.7498/aps.73.20231472
    [2] 黎章龙, 胡长青, 赵梅, 秦继兴, 李整林, 杨雪峰. 基于大掠射角海底反射特性的深海地声参数反演. 物理学报, 2022, 71(11): 114302. doi: 10.7498/aps.71.20211915
    [3] 李郝, 杨鑫, 张正平. THz波在不同角度磁化的非均匀磁化等离子体中的传输特性分析. 物理学报, 2021, 70(7): 075202. doi: 10.7498/aps.70.20201450
    [4] 聂丹丹, 冯晋霞, 戚蒙, 李渊骥, 张宽收. 基于光学参量振荡器的可调谐红外激光的强度噪声特性. 物理学报, 2020, 69(9): 094205. doi: 10.7498/aps.69.20191952
    [5] 张肖肖, 吴振森, 苏翔. 阻抗劈绕射对破碎波后向散射特性的影响. 物理学报, 2016, 65(21): 214101. doi: 10.7498/aps.65.214101
    [6] 相坤生, 程天海, 顾行发, 郭红, 陈好, 王颖, 魏曦, 包方闻. 基于多角度偏振载荷数据的中国典型地物偏振特性研究. 物理学报, 2015, 64(22): 227801. doi: 10.7498/aps.64.227801
    [7] 宋洪胜, 庄桥, 刘桂媛, 秦希峰, 程传福. 菲涅耳深区散斑强度统计特性及演化. 物理学报, 2014, 63(9): 094201. doi: 10.7498/aps.63.094201
    [8] 周天, 李海森, 朱建军, 魏玉阔. 利用多角度海底反向散射信号进行地声参数估计. 物理学报, 2014, 63(8): 084302. doi: 10.7498/aps.63.084302
    [9] 谢东海, 顾行发, 程天海, 余涛, 李正强, 陈兴峰, 陈好, 郭婧. 基于多角度偏振相机的城市典型地物双向反射特性研究. 物理学报, 2012, 61(7): 077801. doi: 10.7498/aps.61.077801
    [10] 姚呈康, 曾晓东, 曹长庆. 机械抖动棱镜式激光陀螺出射光强度特性. 物理学报, 2012, 61(9): 094216. doi: 10.7498/aps.61.094216
    [11] 鄂鹏, 韩轲, 武志文, 于达仁. 磁场强度对霍尔推力器放电特性影响的实验研究. 物理学报, 2009, 58(4): 2535-2542. doi: 10.7498/aps.58.2535
    [12] 王可嘉, 张清泉, 吕健滔, 杜泽明, 刘劲松. 二维无序介质中横磁模的谱线宽度随抽运强度的变化特性. 物理学报, 2008, 57(5): 2941-2945. doi: 10.7498/aps.57.2941
    [13] 余 辉, 江晓清, 杨建义, 戚 伟, 王明华. 窄光束在掠反射条件下的特性. 物理学报, 2008, 57(7): 4208-4213. doi: 10.7498/aps.57.4208
    [14] 程天海, 顾行发, 陈良富, 余 涛, 田国良. 卷云多角度偏振特性研究. 物理学报, 2008, 57(8): 5323-5332. doi: 10.7498/aps.57.5323
    [15] 刘 军, 陈晓伟, 刘建胜, 冷雨欣, 朱 毅, 戴 君, 李儒新, 徐至展. 负啁啾高强度飞秒脉冲在正常色散材料中传输特性研究. 物理学报, 2006, 55(4): 1821-1826. doi: 10.7498/aps.55.1821
    [16] 厉以宇, 顾培夫, 李明宇, 刘 旭, 杨 彗. 波状膜层全角度偏振分束特性的分析. 物理学报, 2006, 55(2): 910-913. doi: 10.7498/aps.55.910
    [17] 尚也淳, 刘忠立, 王姝睿. SiC Schottky结反向特性的研究. 物理学报, 2003, 52(1): 211-216. doi: 10.7498/aps.52.211
    [18] 宋远红, 王友年, 宫 野. 氢离子在固体表面掠角散射与能量损失的数值模拟. 物理学报, 1999, 48(7): 1275-1281. doi: 10.7498/aps.48.1275
    [19] 黄燕霞, 郭光灿. 依赖强度耦合Jaynes-Cummings模型中原子和场的Disentangled-states及其演化特性. 物理学报, 1999, 48(1): 49-57. doi: 10.7498/aps.48.49
    [20] 郑海菜, 苟增光, 崔明, 李秀英, 莫育俊. 固-液系统中表面增强喇曼散射强度ISERS对浓度c的依赖关系及水杨酸分子的吸附等温特性研究. 物理学报, 1992, 41(6): 1027-1031. doi: 10.7498/aps.41.1027
  • 期刊类型引用(3)

    1. 汪磊,黄益旺,郭霖,任超. 浅海粗糙海底声散射建模及声场特性. 物理学报. 2024(03): 91-98 . 百度学术
    2. 郑毅,秦志亮,于盛齐,马本俊,赵吉祥,刘梦婷. 海底声散射特性进展及展望. 哈尔滨工程大学学报. 2024(04): 748-757 . 百度学术
    3. 马国栋. 基于改进侧扫声呐法的水下抛石精准识别研究. 人民长江. 2022(04): 210-214 . 百度学术

    其他类型引用(0)

计量
  • 文章访问数:  7729
  • PDF下载量:  106
  • 被引次数: 3
出版历程
  • 收稿日期:  2018-08-02
  • 修回日期:  2018-12-25
  • 上网日期:  2019-02-01
  • 刊出日期:  2019-02-20

/

返回文章
返回