删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

浅海起伏海面下气泡层对声传播的影响

本站小编 Free考研考试/2021-12-29

摘要:分析了起伏海面下风浪引起的气泡层对海面反射损失和对声传播的影响. 一方面, 气泡层会改变原来水中的声速剖面; 另一方面, 气泡层会对声波产生散射和吸收作用. 考虑以上两方面的因素, 分析了不同风速下气泡层对海面反射损失和声传播损失的影响, 仿真发现, 在风速大于10 m/s时, 对于2 kHz以上频率时气泡层对小掠射角下海面反射损失的影响不可忽视. 在给定的水声环境中, 当声源深度和接收深度都为7 m时, 风速为16 m/s的风浪下生成的气泡层, 在10 km处对3 kHz的声传播损失的影响达到8.1 dB. 当声源深度和接收深度都为18 m时, 风速为16 m/s的风浪下生成的气泡层, 在10 km处对3 kHz的声传播损失的影响达到4 dB.
关键词: 起伏海面/
气泡层/
声传播

English Abstract


--> --> -->
理想状态下的平整海面对声波的反射系数为–1, 这种良好的反射作用不会引起海面反射损失; 而海面风浪作用引起的起伏海面对声波同时具有反射和散射的作用[1-4], 会引起海面反射损失[5-7], 同时风浪搅拌引起的海面附近的气泡层既能改变原来的水中声速剖面, 还会对声波具有散射和吸收作用[8-10], 因此起伏海面下的气泡层也是影响声传播的重要因素之一.
国内外****先后研究了强风浪时起伏海面引起的信号相关性的改变[11]、信号到达时间的起伏[12,13]等现象. 2016年, Jones等[14]以起伏海面下的声传播模型——Ramsurf模型为基准比较了不同风速下的海面反射损失模型; 2018年, Zou和Badiey[15]指出了起伏海面对经过一次海面反射的宽带信号声压幅度的影响; 2019年, Liu和Li[16]研究了不同风速的起伏海面下浅海声场的波动特性, 但都没有同时考虑强风浪引起的气泡层对声传播的影响.
强风浪时需要同时考虑起伏海面和近海面气泡层对声传播的影响. 2005年, Ainslie[17]基于HN气泡群模型考虑了风浪下气泡层对海面反射损失的影响, 包含了风力引起的气泡层对入射声波的折射效应以及吸收和散射效应, 但他只给出了风速为10 m/s时的单位体积内的气泡总体积计算式, 在具体计算时对于粗糙海面的反射系数计算和气泡引起的声场衰减系数计算也相对简化.
本文将起伏海面下气泡层对声传播的影响研究工作更推进一步, 在具体计算时, 采用适用性更广的小斜率散射理论来计算粗糙海面的反射系数, 并且考虑了气泡层引起的声波衰减系数对频率的依赖关系; 另外, 在分析气泡层对海面反射损失的影响基础上, 基于起伏海面下的声传播模型—Ramsurf模型, 考虑气泡层引起的声波折射效应和吸收散射效应, 将水中的声速剖面进行修正, 将气泡层引起的声衰减系数作为水体衰减系数输入到Ramsurf模型中, 分析了起伏海面下气泡层对声传播损失的影响. 本文全面考虑了风力引起的起伏海面和气泡层两种因素, 直接着眼于声传播损失的计算, 可以直接应用于声纳方程中.
2
2.1.海面声反射系数求解
-->粗糙海面散射(见图1)的理论方面, 目前较为热点的是小斜率近似散射理论[18,19]. 根据小斜率近似理论, 可将海面粗糙度表示为幂率粗糙度谱: 入射到粗糙海面$z = f(x)$上的平面波为$\exp ({\rm{i}}{{{k}}_{\rm{i}}} \cdot {{r}})$, 其中 ${{r}} = (x, z)$, ${{{k}}_{\rm{i}}} = ({k_{{\rm{i}}x}}, {k_{{\rm{i}}z}}) = ({k_{{\rm{i}}x}}, - {\kappa _{{\rm{i}}z}})$, 且有${\kappa _{{\rm{i}}z}} = $${\left[ {{k^2} - k_{{\rm{i}}x}^2} \right]^{{1 / 2}}} > 0$, $k = {\omega / c}$是辐射波数.
图 1 一维粗糙海面散射示意图
Figure1. Scattering geometry of one-dimensional rough sea surface.

该表面下方的散射场${\psi _{\rm{s}}}$可以用T (传递)矩阵表示为
${\psi _{\rm{s}}}({{r}},{k_{{\rm{i}}x}}) = \int {{\rm{d}}{k_{{\rm{s}}x}}} \exp \left[ {{\rm{i}}{{{k}}_{\rm{s}}} \cdot {{r}}} \right]{{T}}\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right).$
对于一维表面$f(x)$, 仅考虑水平和竖直方向上的位移, Voronovich提出将传递矩阵T表示为
$\begin{split} & {{T}}({k_{{\rm{s}}x}},{k_{{\rm{i}}x}}) \\=\, & A({k_{{\rm{s}}x}},{k_{{\rm{i}}x}})\int {\rm{d}}x \exp {\left. {\left[ {{\rm{i}}\left( {{{{k}}_{\rm{i}}} - {{{k}}_{\rm{s}}}} \right) \cdot {{r}}} \right]} \right|_{z = f(x)}} \\ &\times \varPhi \left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}},x} \right), \end{split}$
其中函数$\varPhi \left( {{k_{{\rm{s}}x}}, {k_{{\rm{i}}x}}, x} \right)$与起伏表面$f(x)$有关.
传递矩阵T可以用小斜率序列表示:
${{T}} = {{{T}}_0} + {{{T}}_1} + {{{T}}_2} + \cdot \cdot \cdot ,$
其中${{{T}}_0}$可认为是一阶斜率函数, ${{{T}}_1}$是二阶斜率函数, ${{{T}}_2}$是三阶斜率函数, 依此类推. 等式右边的前两项可以写为
${{{T}}_0}\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right) \!=\! - \frac{{2{\kappa _{{\rm{i}}z}}}}{{{v_z}}}\frac{1}{{2{\text{π}}}}\!{\int \!{{\rm{d}}x \! \cdot \exp \left. {\left[ {{\rm{i}}{{v}} \cdot {{r}}} \right]} \right|} _{z = f(x)}},$
$\begin{split}{{{T}}_1}\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right) =\, & - \frac{{{\rm{i}}{\kappa _{{\rm{i}}z}}}}{{{v_z}}}\frac{1}{{2{\text{π}}}}{\int {{\rm{d}}x \cdot \exp \left. {\left[ {{\rm{i}}{{v}} \cdot {{r}}} \right]} \right|} _{z = f(x)}}\\ & \times \int {{\rm{d}}{K_1}} \cdot \exp \left[ {{\rm{i}}{K_1}x} \right]F({K_1}) \\ & \times \left[ {k{\beta _{1 + {\rm{i}}}} + k{\beta _{{\rm{s}} - 1}} - {v_z}} \right],\\[-10pt]\end{split}$
式中, ${{v}} \!=\! {{{k}}_{\rm{i}}}\! - \!{{{k}}_{\rm{s}}} = ({v_x}, - {v_z})$, ${v_x} \!= \!{k_{{\rm{i}}x}} - {k_{{\rm{s}}x}}$, ${v_z} = $$ {\kappa _{{\rm{i}}z}} +{k_{{\rm{s}}z}}$; $k{\beta _{1 + {\rm{i}}}} = {\big[{{k^2} - {{\left( {{K_1} \!+\! {k_{{\rm{i}}x}}} \right)}^2}} \big]^{{1 / 2}}}$, 并且有${\rm{Im}} \left[ {{\beta _{1 + {\rm{i}}}}} \right] > 0$; $k{\beta _{{\rm{s}} - 1}} \!=\! {\big[ {{k^2} \!- \!{{\left( {{k_{{\rm{s}}x}}\! -\! {K_1}} \right)}^2}} \big]^{{1 / 2}}}$, 并且有${\rm{Im}} \left[ {{\beta _{{\rm{s}} - 1}}} \right]\! >\! 0$; $F(K)$$f(x)$的傅里叶变换. 则对于一维表面双基地散射截面可以写成如下的形式:
$\begin{split} \sigma \delta \left( {{k_{{\rm{i}}x}} - k_{{\rm{i}}x}'} \right) =\, &\frac{{k_{{\rm{s}}z}^2}}{k}\left[\left\langle {T\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right)T{{\left( {{k_{{\rm{s}}x}},k_{{\rm{i}}x}'} \right)}^ * }} \right\rangle\right. \\ & \left. - {{\left\langle {T\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right)T\left( {{k_{sx}},k_{{\rm{i}}x}'} \right)} \right\rangle }^ * } \right], \end{split}$
散射截面表示成
$\sigma = {\sigma _{00}} + {\sigma _{01}} + {\sigma _{11}} + \cdot \cdot \cdot ,$
其中${\sigma _{00}}$为二阶散射截面, ${\sigma _{01}}$, ${\sigma _{11}}$则分别对应第三和第四阶, 即
${\sigma ^{(2)}} = {\sigma _{00}},$
${\sigma ^{(3)}} = {\sigma ^{(2)}} + {\sigma _{01}},$
${\sigma ^{(4)}} = {\sigma ^{(3)}} + {\sigma _{11}}.$
此方法与扰动理论式是一样的, 并且包含了变量$F(K)$, $f(x)$存在于指数函数${\left. {\exp \left[ {{\rm{i}}{{v}} \cdot {{r}}} \right]} \right|_{z = f(x)}}$中, 对传递矩阵进行平均得到
$\begin{split}& \left\langle {{T_0}\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right)} \right\rangle \\ =\, & - \frac{{2{\kappa _{{\rm{i}}z}}}}{{{v_z}}}\frac{1}{{2{\text{π}}}}\int {\exp [{\rm{i}}{v_x}x]} \left\langle {\exp [ - {\rm{i}}{v_z}f(x)]} \right\rangle {\rm{d}}x \\= \, &- \exp [{{ - {\chi ^2}} / 2}]\delta \left( {{k_{{\rm{s}}x}} - {k_{{\rm{i}}x}}} \right),\end{split} $
$\begin{split} & \left\langle {{T_1}\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right)} \right\rangle \\ =\, & - \frac{{{\rm i}{\kappa _{{\rm{i}}z}}}}{{{v_z}}}\frac{1}{{{{\left( {2{\text{π}}} \right)}^2}}}\int {\exp [{\rm{i}}{v_x}x]} {\rm{d}}x \\ &\times \int {\exp [{\rm{i}}{K_1}x]} \left( {k{\beta _{1 + {\rm{i}}}} + k{\beta _{{\rm{s}} - 1}} - {v_z}} \right){\rm{d}}{K_1}\\&\times \int {\exp [ - {\rm{i}}{K_1}{x'}]} \left\langle {\exp [ - {\rm{i}}{v_z}f(x)]} \right\rangle {\rm{d}}{x'}\\ =\, & {\kappa _{{\rm{i}}z}}\exp [{{ - {\chi ^2}} / 2}]{\rm{\delta }}\left( {{k_{{\rm{s}}x}} - {k_{{\rm{i}}x}}} \right)\\ & \times \int {{\rm{d}}{K_1}} \;W\left( {{K_1}} \right)g\left( {{K_1}} \right),\end{split}$
其中$\chi = {v_z}h$, $g({K_1}) = k{\beta _{{\rm{i + }}1}} + k{\beta _{{\rm{s}} - 1}} - {v_z}$, $W({K_1})$为表面粗糙度幅度谱, 它是粗糙面相关函数的傅里叶变换. $h$为粗糙海面的波高均方根.
对于$\left\langle {{T}} \right\rangle $的小斜率近似, 使用前两项已相对精确:
$\left\langle {{T}} \right\rangle = \left\langle {{{{T}}_0}} \right\rangle + \left\langle {{{{T}}_1}} \right\rangle ,$
$\left\langle {{{{T}}_0}} \right\rangle = \exp \left[ { - \frac{{v_z^2{h^2}}}{2}} \right]{\rm{\delta }}\left( {{k_{sx}} - {k_{{\rm{i}}x}}} \right),$
$\begin{split} \left\langle {{{{T}}_1}} \right\rangle =\, & {\kappa _{{\rm{i}}z}}\exp \left[ { - \frac{{v_z^2{h^2}}}{2}} \right]{\rm{\delta }}\left( {{k_{{\rm{s}}x}} - {k_{{\rm{i}}x}}} \right)\\ & \times \int {{\rm{d}}{K_1}} W\left( {{K_1}} \right)g\left( {{K_1}} \right).\end{split}$
计算反射系数的方法如下:
$\left\langle {{{T}}\left( {{k_{{\rm{s}}x}},{k_{{\rm{i}}x}}} \right)} \right\rangle = R\left( {{k_{{\rm{i}}x}}} \right){\rm{\delta }}\left( {{k_{{\rm{i}}x}} - {k_{{\rm{s}}x}}} \right).$
因此可以得到小斜率近似下的反射系数为
$\begin{split} R\left( {{k_{{\rm{i}}x}}} \right) \approx\, & \exp \left[ { - \frac{{v_z^2{h^2}}}{2}} \right]\\ &\times\left[ { - 1 + {\kappa _{{\rm{i}}z}}\int {{\rm{d}}{K_1}} W\left( {{K_1}} \right)g\left( {{K_1}} \right)} \right], \end{split}$
其中${k_{{\rm{i}}x}} \equiv {k_{{\rm{s}}x}} = k\sin {\theta _{\rm{i}}}$, ${\kappa _{{\rm{i}}z}} \equiv \left| {{k_{{\rm{s}}z}}} \right| = k\cos {\theta _{\rm{i}}}$, ${v_z} = $$2 k\cos {\theta _{\rm{i}}} $. 其中$W({K_1})$作为表面粗糙度幅度谱, 与海面波浪谱有关. 目前已有多种形式的海面波浪谱, 本文将采用国际上最为广泛应用的PM谱[20,21].
对于PM谱, 一维粗糙度幅度谱$W({K_1})$[6]
$ W({K_1}) = \left\{ \begin{aligned}& \left[\frac{a}{{8{{\left| {{K_1}} \right|}^3}}} \right] \cdot \frac{2}{{\text{π}}}\left[\arcsin \left( {\frac{{\left| {{K_1}} \right|}}{{{K_{\rm{L}}}}}} \right)\right.\\ & \left. - \frac{{\left| {{K_1}} \right|\sqrt {K_{\rm{L}}^2 - K_1^2} }}{{K_{\rm{L}}^2}}\right], ~ \left| {{K_1}} \right| < {K_{\rm{L}}}, \\ &\left[ {\frac{a}{{8{{\left| {{K_1}} \right|}^3}}}} \right],\qquad\quad\quad\quad\left| {{K_1}} \right| > {K_{\rm{L}}},\end{aligned} \right. $
其中$a = 8.1 \times {10^{ - 3}}$, $g$是重力加速度, $U$为海面上方19.5 m高处的风速, ${K_{\rm{L}}}$是与风速有关的量: ${K_{\rm{L}}} = {{\left( {\sqrt \beta g} \right)} / {{U^2}}}$, $\beta = 0.74$.
2
2.2.气泡层结构的声反射系数求解
-->3
2.2.1.气泡群建模
-->Ainslie[17]对Hall-Novarini (HN)气泡群模型[22]进行了修正, 假设气泡半径大小在10—1000 μm之间, 则气泡群的分布函数为
$n(a) = (1.6 \times {10^{10}})G(a,z){\left( {\frac{{{v_{10}}}}{{13}}} \right)^3} \times \exp \left[ { - \frac{z}{{L({v_{10}})}}} \right],$
其中$a$为气泡半径, 单位为μm; $z$为深度, 单位为m, ${v_{10}}$为海面上方10 m处的风速, 单位为${\rm{m/s}}$, $L({v_{10}})$是与风速有关的量, 单位为${\rm{m}}$; $G(a, z)$是与气泡半径和水深有关的量:
$L({v_{10}}) = \left\{ \begin{aligned}& 0.4, \quad\quad\quad\quad\quad\quad\quad\quad\quad {v_{10}} \leqslant 7.5,\\& 0.4 + 0.115({v_{10}} - 7.5),\quad {v_{10}} > 7.5,\end{aligned} \right.$
$G(a,z) = \left\{ \begin{aligned}& {\left[ {{a_{{\rm{ref}}}}/a} \right]^4},~\qquad {a_{\min }} \leqslant a \leqslant {a_{{\rm{ref}}}}\left( z \right),\\& {\left[ {{a_{{\rm{ref}}}}/a} \right]^{x(z)}},~\quad{a_{{\rm{ref}}}}(z) < a \leqslant {a_{\max }}.\end{aligned} \right.$
$G(a, z)$表达式中的两个变量${a_{{\rm{ref}}}}\left( z \right)$$x\left( z \right)$与水深$z$的关系分别为:
${a_{{\rm{ref}}}}\left( z \right) = 54.4 + 1.984 \times {10^{ - 6}}z,$
$x\left( z \right) = 4.37 + {\left( {\frac{z}{{2.55}}} \right)^2}.$
1)气泡层对声波的折射效应
Ainslie[17]给出了修正后的含气泡层的水中声速:
$\frac{{c_{\rm{w}}^2}}{{{c_{\rm{m}}}{{\left( z \right)}^2}}} = 1 + \frac{{{\rho _{\rm{w}}}c_{\rm{w}}^2}}{{{\kappa _0}\left( z \right)P(z)}}U(z),$
其中${c_{\rm{w}}}$是不含气泡的水中声速, ${c_{\rm{m}}}$是修正后的含气泡水中声速. 其中$P(z)$是水中的绝对流体静力压强, $U(z)$是单位体积内气泡的总体积, Ainslie仅给出了风速为10 m/s时的$U(z)$表达式, 本文推导了不同风速下的$U(z)$表达式:
$\begin{split} U(z) =\, & \int_{{a_{\min }}}^{{a_{\max }}} {V(a)n(a){\rm{d}}a} = \frac{4}{3}{\text{π}}\int_{{a_{\min }}}^{{a_{\max }}} {{a^3}} n(a){\rm{d}}a \\ = \, &1.6 \times {10^{10}} \times {\left( {\frac{{{v_{10}}}}{{13}}} \right)^3} \times \exp \left( { - \frac{z}{L}} \right) \\ &\times \frac{4}{3}{\text{π}}\left\{ {\int_{{a_{\min }}}^{{a_{{\rm{ref}}}}} {\frac{{a_{{\rm{ref}}}^4\left( z \right)}}{a}{\rm{d}}a + \int_{{a_{{\rm{ref}}}}}^{{a_{\max }}} {\frac{{a_{{\rm{ref}}}^{x(z)}\left( z \right)}}{{{a^{x(z) - 3}}}}{\rm{d}}a} } } \right\} \\=\, &1.6 \times {10^{10}} \times {\left( {\frac{{{v_{10}}}}{{13}}} \right)^3} \times \exp \left( { - \frac{z}{L}} \right) \\ &\times \frac{4}{3}{\text{π}}a_{{\rm{ref}}}^4\left( z \right)\left\{\ln \frac{{{a_{{\rm{ref}}}}\left( z \right)}}{{{a_{\min }}}} \right. \\ &+\left. \frac{{1 - {{\left[ {{{{a_{{\rm{ref}}}}(z)} / {{a_{\max }}}}} \right]}^{x(z) - 4}}}}{{x(z) - 4}} \right\}.\end{split} $
在频率为2 kHz、水中声速${c_{\rm{w}}}$为1490 m/s时(温度取10 ℃, 盐度取35 psu), 海表面风速分别为7, 10, 13 和16 m/s时的水中声速剖面的变化如图2所示, 海表面0—10 m以内的水中声速有不同程度的改变, 在风速为13 m/s时, 海面表层水中声速的改变达到32 m/s左右.
图 2 气泡层引起的水中声速剖面的改变
Figure2. variation of the sound speed profile due to bubbles layer.

2)气泡层对声波的吸收散射效应
水中气泡对声能量具有散射和吸收的作用, 引起声能量的衰减. 衰减系数的计算方法为[8]
${\alpha _{\rm{b}}}\left( z \right) = 4.34\int_{{a_{\min }}}^{{a_{\max }}} {\frac{{4{\text{π}}{a^2}n\left( a \right)\left[ {{\delta / {(ka)}}} \right]}}{{{{\left( {{Y^2} - 1} \right)}^2} + {\delta ^2}}}} {\rm{d}}a, $
其中$Y = {{{f_{\rm{r}}}} / f}$, $f$是声频率, ${f_{\rm{r}}}$是半径为$a$的气泡在水深为$z$处的共振频率, 另外, $\delta $是阻尼常数. 图3是频率为2 kHz时不同海表风速下风浪引起的气泡层的水中衰减系数计算结果.
图 3 气泡层引起的声能量的衰减系数
Figure3. Extinction coefficient due to bubbles layer.

3
2.2.2.气泡层多层结构的声反射系数求解
-->研究气泡层对海面反射损失的影响时, 由于气泡层对声波的折射效应和消减效应随水深而变化, 因此考虑将气泡层划分为在垂直方向上的多层结构, 这里需要用到一种求解多层结构的反射系数的递归方法[23].
图4所示, 将气泡混合层划分为垂直方向上的多层结构, 包括最底层的layer 1和最顶层的layer m, 以及中间的m – 2个隔层, 最上面3层的反射系数为
图 4 气泡层的多层结构示意图
Figure4. Layered structure geometry of the bubbles layer.

${R_{(m - 2)m}} = \frac{{{R_{(m - 2)(m - 1)}} + {R_{(m - 1)m}}\exp (2{\rm{i}}{\phi _{m - 1}})}}{{1 + {R_{(m - 2)(m - 1)}}{R_{(m - 1)m}}\exp (2{\rm{i}}{\phi _{m - 1}})}},$
计算出${R_{(m - 2)m}}$, 下移一层至m – 3层, 计算下一个反射系数, 即
${R_{(m - 3)m}} = \frac{{{R_{(m - 3)(m - 2)}} + {R_{(m - 2)m}}\exp (2{\rm{i}}{\phi _{m - 2}})}}{{1 + {R_{(m - 3)(m - 2)}}{R_{(m - 2)m}}\exp (2{\rm{i}}{\phi _{m - 2}})}}.$
由此递归算法, 最后可求得总的反射系数为
${V_{\rm{c}}} = {R_{1m}} = \frac{{{R_{12}} + {R_{2m}}\exp (2{\rm{i}}{\phi _2})}}{{1 + {R_{12}}{R_{2m}}\exp (2{\rm{i}}{\phi _2})}}, $
其中${\phi _j} \equiv {k_j}{h_j}\sin {\theta _j}$是声场穿过厚度为${h_m}$的薄层的路径产生的与角度有关的垂直相移. 其中${k_j}$是第j层的波数, ${h_j}$是第j层的厚度, ${\theta _j}$是第j层的入射波的掠射角.
考虑小掠射角入射的声波, 水中声速为1490 m/s等声速剖面, 分别在风速为7, 10和13 m/s下计算海面反射损失(分别对应图5, 图6图7), 以$ - \ln \left| {{V_{\rm{c}}}} \right|$计算海面反射损失(单位为奈培).
图 5 风速为7 m/s时起伏海面下气泡层对海面反射损失的影响 (a) 无气泡层; (b) 考虑气泡层对声波的消减效应; (c) 考虑气泡层对声波的折射效应; (d) 同时考虑气泡层对声波的折射效应和消减效应
Figure5. Effects of the bubbles layer underneath the rough sea surface on reflection loss in nepers with a wind speed of 7 m/s: (a) Rough sea surface; (b) rough sea surface + scattering and absorption effect of the bubbles layer; (c) rough sea surface + refraction effect of the bubbles layer; (d) rough sea surface + scattering and absorption effect of the bubbles layer + refraction effect of the bubbles layer.

图 6 风速为10 m/s时起伏海面以及气泡层对海面反射损失的影响 (a) 无气泡层; (b) 考虑气泡层对声波的消减效应; (c) 考虑气泡层对声波的折射效应; (d) 同时考虑气泡层对声波的折射效应和消减效应
Figure6. Effects of the bubbles layer underneath the rough sea surface on reflection loss in nepers with a wind speed of 10 m/s: (a) Rough sea surface; (b) rough sea surface + scattering and absorption effect of the bubbles layer; (c) rough sea surface + refraction effect of the bubbles layer; (d) rough sea surface + scattering and absorption effect of the bubbles layer + refraction effect of the bubbles layer.

图 7 风速为13 m/s时起伏海面以及气泡层对海面反射损失的影响 (a) 无气泡层; (b) 考虑气泡层对声波的消减效应; (c) 考虑气泡层对声波的折射效应; (d) 同时考虑气泡层对声波的折射效应和消减效应
Figure7. Effects of the bubbles layer underneath the rough sea surface on reflection loss in nepers with a wind speed of 13 m/s: (a) Rough sea surface; (b) rough sea surface + scattering and absorption effect of the bubbles layer; (c) rough sea surface + refraction effect of the bubbles layer; (d) rough sea surface + scattering and absorption effect of the bubbles layer + refraction effect of the bubbles layer.

图5(a)是不考虑气泡层影响时起伏海面的海面反射损失计算, 图 5(b)图5(a)的基础上加入了气泡层对声波的消减效应, 不同风速下气泡层引起的声波衰减系数如图3所示; 图5(c)图5(a)的基础上加入了气泡层对声波的折射效应, 不同风速下气泡层引起的声速剖面的变化如图2所示, 将海表以下0—10 m的水体分为100层, 利用上述递归方法逐层计算, 最顶层的反射系数利用(17)式求得; 图5(d)是在图5(a)的基础上同时考虑了气泡层对声波的折射效应和消减效应. 后续图6图7的子图考虑的因素排序与图5一致.
观察图5, 在风速为7 m/s时, 起伏海面的粗糙度较小, 此时由于风浪搅拌引起的气泡层对声速剖面结构的改变以及对声波的散射作用都较小, 因此, 海面反射损失相对较小; 图5(c)图5(a)相比, 气泡层的折射效应对于3 kHz以上的海面反射损失影响较为明显. 观察图6, 图6(c)图6(a)相比, 在频率大于2 kHz时, 起伏海面下气泡层因素很大程度地改变了海面反射损失; 图6(c)图6(d)相差较小, 此时, 相对于气泡层对声波的散射和吸收效应, 气泡层对声波的折射效应起主要作用. 可见, 风速大于10 m/s时, 在频率大于2 kHz时, 海面粗糙度和气泡层对小掠射角下声传播的影响不可忽视.
观察图7, 相对于图5图6, 由于风速的增加, 图7中起伏海面下的海面反射损失以及气泡层对海面反射损失的影响都较大, 并且在风速为13 m/s时, 气泡混合层对1 kHz以上小掠射角下的海面反射损失的影响已经较为明显.
本文提出的起伏海面下含有气泡层时的声传播计算方法的总体思路是: 以起伏海面下的Ramsurf声传播模型为模型基础, 利用Monte-Carlo方法[24]产生一维PM谱起伏海面作为海面边界的输入参数; 在考虑气泡层的影响时, 将不含气泡的水中声速剖面用(20)式进行修正得到含有气泡层的水中声速剖面(如图2所示), 作为模型中的声速剖面参数; 同时利用(21)式计算气泡层引起的声波衰减系数(如图3所示)作为模型的水中衰减参数.
水声环境为80 m水深, 声源深度SD为7 m, 接收深度RD为7 m; 不含气泡的水中声速为1490 m/s等声速, 考虑气泡层影响时用(20)式进行修正. 分别在风速为7, 10, 13和16 m/s下计算起伏海面下不考虑气泡层和起伏海面下考虑气泡层时的传播损失曲线. 用宽带非相干方法以1/3倍频程的带宽计算中心频率为3 kHz的传播损失曲线(频率间隔取10 Hz), 每个频率采用40次Monte-Carlo方法的平均, 以下图8(a)图8(d)分别是风速为7, 10, 13和16 m/s时的声传播损失曲线, 其中红色曲线是平整海面下计算所得, 绿色曲线是起伏海面下不考虑气泡混合层因素计算所得, 蓝色曲线是起伏海面下考虑气泡层因素计算所得.
图 8 不同风速时平整海面、起伏海面不考虑气泡层因素和起伏海面考虑气泡层因素时的传播损失曲线比较(SD = 7 m, RD = 7 m) (a) v = 7 m/s; (b) v = 10 m/s; (c) v = 13 m/s; (d) v = 16 m/s
Figure8. Comparison of the transmission loss curves with different wind speeds (SD = 7 m, RD = 7 m): (a) v = 7 m/s; (b) v = 10 m/s; (c) v = 13 m/s; (d) v = 16 m/s.

观察图8(a), 风速为7 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差5 dB, 而绿色曲线和蓝色曲线几乎重合, 这是由于气泡层引起的衰减系数以及声速剖面的结构上的改变都较小, 因此, 气泡层对声传播损失的影响较小. 观察图8(b), 风速为10 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差4.5 dB, 绿色曲线和蓝色曲线相差0.4 dB. 观察图8(c), 风速为13 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差6.7 dB, 绿色曲线和蓝色曲线相差2.6 dB. 观察图8(d), 风速为16 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差11.2 dB, 绿色曲线和蓝色曲线相差8.1 dB. 可见, 风速越大, 起伏海面和气泡层对声传播的影响越大.
保持其他条件不变, 将声源深度SD和接收深度RD均设为18 m, 利用同样的方法分别在风速为7, 10, 13和16 m/s下计算起伏海面下不考虑气泡层和起伏海面下考虑气泡层时的传播损失曲线, 以下图9(a)图9(d)分别是风速为7, 10, 13和16 m/s时的声传播损失曲线.
图 9 不同风速时平整海面、起伏海面不考虑气泡层因素和起伏海面考虑气泡层因素时的传播损失曲线比较(SD = 18 m, RD = 18 m) (a) v = 7 m/s; (b) v = 10 m/s; (c) v = 13 m/s; (d) v = 16 m/s
Figure9. Comparison of the transmission loss curves with different wind speeds (SD = 18 m, RD = 18 m.): (a) v = 7 m/s; (b) v = 10 m/s; (c) v = 13 m/s; (d) v = 16 m/s.

观察图9(a), 风速为7 m/s时, 与图8(a)类似, 绿色曲线和蓝色曲线几乎重合, 风速较小时, 气泡层对声传播损失的影响较小. 观察图9(b), 风速为10 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差8.6 dB, 绿色曲线和蓝色曲线相差0.8 dB. 观察图9(c), 风速为13 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差9.9 dB, 绿色曲线和蓝色曲线相差2.5 dB. 观察图9(d), 风速为16 m/s时, 在距离为10 km处, 红色曲线和绿色曲线相差10 dB, 绿色曲线和蓝色曲线相差4 dB. 根据HN气泡群模型风浪引起的气泡层只改变了10 m水深范围以内的声速和声衰减系数, 但由于波导是一个整体, 风浪引起的气泡层仍然会影响发射深度和接收深度大于10 m时的声传播.
大风浪下在近起伏海面处会产生气泡层, 气泡层不但会改变原来的水中声速剖面, 还会因为散射吸收等引起声波的衰减, 本文提出了起伏海面下含有气泡层时的声传播计算方法, 并且从海面反射损失和声传播损失两个方面分析了起伏海面下气泡层对声传播的影响, 仿真实验表明, 在风速为10 m/s以上时, 气泡层对于频率大于2 kHz时小掠射角下的海面反射损失的影响较为显著; 在风速为13 m/s以上时, 气泡层对于频率大于1 kHz时小掠射角下的海面反射损失的影响较为显著. 在给定的水声环境中, 当声源深度和接收深度都为7 m时, 风速为13 m/s的风浪下生成的气泡层, 在10 km处对3 kHz的声传播损失的影响达到2.6 dB; 风速为16 m/s的风浪下生成的气泡层, 在10 km处对3 kHz的声传播损失的影响达到8.1 dB. 当声源深度和接收深度都为18 m时, 风速为13 m/s的风浪下生成的气泡层, 在10 km处对3 kHz的声传播损失的影响达到2.5 dB; 风速为16 m/s的风浪下生成的气泡层, 在10 km处对3 kHz的声传播损失的影响达到4 dB.
本文提出的起伏海面下含有气泡层时的声传播计算方法, 可为高海况下大风浪天气时的声场预报等提供模型基础和算法支撑, 不同风速下的仿真数据结果可为目标定位、水声环境监测等提供数据参考; 另外, 本文全面直接着眼于起伏海面下风浪引起的气泡层对声传播损失的影响, 可以直接应用于声纳方程中. Ramsurf模型在计算时只考虑前向散射. 对于强起伏界面下后向散射为主的情况, 是下一步的工作方向.
相关话题/传播 计算 结构 气泡 海面

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 缺陷组合嵌入VO<sub>2</sub>薄膜结构的可调太赫兹吸收器
    摘要:提出一种多缺陷组合嵌入VO2薄膜结构的可调太赫兹吸收器,它由上表面金属图案层、基体和底层金属板三层结构组成,在上表面和基体之间嵌入二氧化钒介质.计算结果表明在f=4.08THz和f=4.33THz两频点吸收率分别为99.8%和99.9%.通过改变外界环境温度可控制二氧化钒相变,从而使两个频点吸 ...
    本站小编 Free考研考试 2021-12-29
  • 天线方向系数的一类计算逼近方法
    摘要:天线的方向系数是天线的核心性能指标之一,准确计算方向系数是高性能天线应用的核心要求.本文基于平面近场测试理论、实测数据和快速傅里叶变换算法,系统阐述基于近场测试来数值计算天线方向系数的原理,并进行深入的误差分析.本文选择一种方向图函数和方向系数已知的被测天线,来检验所讨论的误差评估方案.评估分 ...
    本站小编 Free考研考试 2021-12-29
  • 有限步传播范围期望指标判别节点传播影响力
    摘要:在线社交网络逐渐成为人们不可或缺的重要工具,识别网络中具有高影响力的节点作为初始传播源,在社会感知与谣言控制等方面具有重要意义.本文基于独立级联模型,给出了一个描述有限步传播范围期望的指标—传播度,并设计了一种高效的递推算法.该指标在局部拓扑结构信息的基础上融合了传播概率对影响力进行刻画,能够 ...
    本站小编 Free考研考试 2021-12-29
  • 二维材料/铁电异质结构的研究进展
    摘要:二维材料是一类具有原子层厚度的层状材料,拥有独特的电学、磁学、光学和力学性能.以石墨烯和过渡金属硫族化合物为代表的二维材料展现出迁移率高、能带可调、可见光透过率高等特点,是近年来微纳科学领域的前沿热点.将二维材料与各种功能材料,如SiO2绝缘体、半导体、金属、有机化合物等结合,可以深化和拓宽二 ...
    本站小编 Free考研考试 2021-12-29
  • 基于密度泛函理论的外电场下C<sub>5</sub>F<sub>10</sub>O的结构及其激发特性
    摘要:采用密度泛函(DFT)B3LYP/6-311g(d)对C5F10O分子进行几何结构优化,研究外加电场(0-0.03a.u.,1a.u.=5.142×1011V/m)对分子的几何结构、能量、前线轨道能级、红外光谱的影响.在相同基组下,采用TD-DFT方法计算和分析C5F10O的轨道成分和激发特性 ...
    本站小编 Free考研考试 2021-12-29
  • 具有内参考热补偿功能的三层膜结构微球腔折射率传感器
    摘要:光学微腔在高灵敏度传感中有着重要的应用前景,而在传感中热漂移是制约其走向实用的重要因素.本文提出了一种镀有三层膜结构的微球腔,可以在实现高灵敏度折射率传感的同时,具备内参考热补偿功能.该结构由内到外分别涂覆折射率为高、低、高的薄膜,内外两高折射率层可以分别支持各自的回音壁模式,称之为内层模式和 ...
    本站小编 Free考研考试 2021-12-29
  • He离子辐照对石墨烯微观结构及电学性能的影响
    摘要:本文采用5.4keV不同剂量的He离子辐照单层石墨烯,利用X射线光电子能谱(XPS)、拉曼光谱(Raman)和半导体参数分析仪表征辐照前后石墨烯微观结构和电学性能变化.研究结果表明:随着辐照剂量增大,单层石墨烯的缺陷密度逐渐增加,当辐照剂量增至1.6×1013He+/cm2,石墨烯开始由纳米晶 ...
    本站小编 Free考研考试 2021-12-29
  • 含石墨烯分界面有耗分层介质的传播矩阵
    摘要:给出一种适用于含导电界面的有耗分层介质的传播矩阵方法.利用相位匹配原理给出斜入射时有耗介质波矢量的实部和虚部,二者方向不同使得在介质中传播非均匀平面波.根据边界条件,推导了跨越石墨烯界面的传播矩阵,以及“无限薄”石墨烯层的反透射系数解析式.最终将传播矩阵方法推广应用于含石墨烯界面的有耗分层介质 ...
    本站小编 Free考研考试 2021-12-29
  • TiAl电子态结构的<i>ab initio</i>计算
    摘要:应用完全活动基自洽场方法,结合N电子价态微扰近似(NEVPT2),对TiAl金属二聚体的基态和若干最低电子激发态的势能曲线进行了计算.完全活动空间由Al的3个价电子(3s23p1)轨道和Ti的4个价电子(3d24s2)轨道构成,计算基组选用Karlsruhegroup的价分裂全电子基组def2 ...
    本站小编 Free考研考试 2021-12-29
  • 空位及氮掺杂二维ZnO单层材料性质:第一性原理计算与分子轨道分析
    摘要:采用基于密度泛函理论的第一性原理计算方法,系统地研究了带缺陷的二维类石墨烯结构的ZnO(graphenelike-ZnO,g-ZnO)的几何结构、电子结构、磁性性质和吸收光谱性质.研究的缺陷类型包括锌原子空位(VZn_g-ZnO)、氧原子空位(VO_g-ZnO)、氮原子取代氧原子(NO_g-Z ...
    本站小编 Free考研考试 2021-12-29