采用能量守恒和高阶Padé近似的三维水声抛物方程模型

为了充分考虑海底地形随三维空间变化的海洋环境中水平方位角耦合效应对声传播的影响,建立了一种三维柱坐标系下流体高阶抛物方程算法。该算法采用泰勒近似将二维方根算子分裂成一维方根算子,并采用分裂步进的高阶Pade近似将一维方根算子写成微分算子有理分式连乘的形式,进而应用Galerkin离散化方法来处理微分算子,最终将微分方程写成矩阵方程的形式;采用能量守恒近似来处理海底边界,以考虑复杂海底对于声传播的影响;采用交替方向隐式格式,实现了三维声场的步进计算。楔形和海底山等典型海域声场仿真计算表明,相比于已有的声场计算模型,三维柱坐标系下高阶抛物方程模型可以更加精确地计算楔形海域和海底山区域的三维声场,实现水平方位全空间声场计算。     

抛物方程的多重非均匀网格模型及其应用

提出了抛物方程的多重非均匀网格模型,以准确求解三维空间存在多辐射源的电波传播问题。通过对不同辐射源建立不同的坐标系,并对其仿真空间采用不同的非均匀网格划分,构建了抛物方程的多重非均匀网格模型。在此基础上,实现了三维多辐射源问题的并行计算。实例仿真了空间存在四个辐射源的电波传播特性。结果表明,抛物方程的多重非均匀网格模型能够准确求解多源的空间电磁场分布特性,且在该算例中,并行技术使得抛物方程的计算速度提升了2.41倍,极大地提高了抛物方程对三维多源问题的求解效率。     

抛物方程的多重非均匀网格模型及其应用

提出了抛物方程的多重非均匀网格模型,以准确求解三维空间存在多辐射源的电波传播问题。通过对不同辐射源建立不同的坐标系,并对其仿真空间采用不同的非均匀网格划分,构建了抛物方程的多重非均匀网格模型。在此基础上,实现了三维多辐射源问题的并行计算。实例仿真了空间存在四个辐射源的电波传播特性。结果表明,抛物方程的多重非均匀网格模型能够准确求解多源的空间电磁场分布特性,且在该算例中,并行技术使得抛物方程的计算速度提升了2.41倍,极大地提高了抛物方程对三维多源问题的求解效率。     

复等效深度耦合简正波模型

为了考虑海底地形随距离变化的非水平分层介质中割线积分对声场的贡献,提出了复等效深度耦合简正波模型。该耦合简正波模型由介质运动方程和连续性方程推导得到了耦合微分方程组,此方程组满足海底地形随距离变化情况下的边界条件且仅包含一个耦合矩阵,并通过引入复等效深度理论处理连续谱和离散谱之间的相互耦合。仿真计算表明,复等效深度耦合简正波模型提高了波导简正波本征值位于割线枝点附近情况下声传播损失的计算精度,充分考虑了波导简正波、非波导简正波和割线积分对声场的贡献,可快速而准确地计算非水平分层介质中的声场。     

三维浅海下弹性结构声辐射预报的有限元-抛物方程法

利用多物理场耦合有限元法对结构和流体适应性强、抛物方程声场计算高效准确的特点,提出了三维浅海波导下弹性结构声振特性研究的有限元-抛物方程法.该方法采用多物理场耦合有限元理论建立浅海下结构近场声辐射模型,计算局域波导下结构声振信息,并提取深度方向上复声压值作为抛物方程初始值;然后采用隐式差分法求解抛物方程以步进计算结构辐射声场.重点介绍了该方法对浅海下结构声辐射计算的准确性、高效性以及快速收敛性后,对Pekeris波导中有限长弹性圆柱壳的声振特性进行了分析.研究得出,当圆柱壳靠近海面(海底)时,其耦合频率比自由场下的要高(低),当潜深达到一定范围时,与自由场耦合频率基本趋于一致;在低频远场,结构辐射场与同强度点源声场具有一定的等效性,且等效距离随着频率增加而增加;由于辐射声场受结构振动模态、几何尺寸和简正波模式影响,结构辐射场传播的衰减规律按近场声影响区、球面波衰减区、介于球面波和柱面波衰减区、柱面波衰减区四个扩展区依次进行.     

南海海域跨海沟环境的声场会聚特性

利用南海北部海域的声传播实验数据分析了跨海沟条件下的声传播特性。结果表明:由于负梯度声速剖面和海底地形的共同作用,导致声能量在开始随着海沟深度变化向更深层位置上弯曲传播,传播损失在海沟中心位置附近达到最大。在海深逐渐变浅的距离上,由于海底反射使得声能量逐渐会聚,声传播损失比单纯陆坡变到深海环境下要减小20 dB以上。当会聚效果小于扩散和海底反射等引起的损失时,声传播损失达到最小,之后随着距离增大而增大;在海沟最后一段,当海深在8 km范围内从850 m突然变浅到311 m,声场能量逐渐会聚效果再次凸显,使得声传播损失减小10 dB。利用射线理论和抛物方程近似数值分析,解释了海沟环境下的由于地形变化引起的声场会聚传播特性。     

陆架斜坡海域上坡波导环境中声能量急剧下降现象及其定量分析

深度较浅的声源其辐射声波在陆架斜坡海域上坡传播时,在斜坡顶端会出现声能量急剧下降现象.利用射线声学模型分析了造成这一现象的原因,并根据抛物方程声场模型计算的深海和浅海平均传播损失定义了"声能量急剧下降距离",定量分析了声源位置对该现象的影响.结果表明:声源深度对"声能量急剧下降距离"影响较大,而声源与斜坡底端水平距离对其影响较小;当声源深度变大时,部分掠射角较小的声线最终能够达到斜坡顶端,致使"声能量急剧下降距离"增大,继续增加声源深度,将导致上坡声能量急剧下降现象消失.利用抛物方程声场模型对陆架斜坡海域上坡声传播进行数值仿真,结合"声能量急剧下降距离"的定义,计算并比较了声源位置不同时该距离的变化,数值计算结果验证了理论分析.     

深海海底山环境下声传播水平折射效应研究

声波在深海海底山环境中传播时,海底山会对声传播产生重要影响.2016年在南海深海进行了一次海底山环境下的声传播实验,观测到了由海底山引起的三维声传播效应,本文利用BELLHOP射线理论解释了海底山环境下的三维声传播机理.结果表明：声波在传播过程中与海底山作用后破坏了深海会聚区结构,导致传播损失增大,在海底山后形成具有明显边界的声水平折射区,利用二维声传播模型无法解释实验现象,海底山后声水平折射区实验测量的声场结构与N×2D模型计算结果存在明显差异,实验的传播损失比N×2D模型计算结果大10 dB.通过三维射线模型分析N×2D模型计算结果与实验结果存在明显差异产生的原因,发现由于声波水平折射作用,部分声线无法到达接收器,使得三维声传播效应对海底山后一定角度范围内声场影响较为明显.因此,深海海底山会引起明显的三维水平折射效应,应在水下目标探测和定位等应用中给予重视.     

空气中声源激发浅海水下声场传播特性的实验研究

为了获得空气中远距离声源激发水下声场的精细结构,2013年3月,声场声信息国家重点实验室在南海海域进行了一次空气中声源激发水下声场的实验。采用汽笛作为空气声源,海底放置水听器作为接收,在实验过程中,发射船由距离水听器2.4 km处行驶至9.8 km。本文对该次实验数据进行分析,获得了收发距离远达9.8 km、频率分别为128 Hz和256 Hz的声传播损失曲线,该曲线随传播距离变化存在清晰的震荡结构.利用波数积分方法计算实验环境下的水下声场理论值,并对获得的声场传播特性进行了较好的物理解释。     

冬季冲绳海域地形及黑潮对声传播的影响

冲绳海域地形复杂且冬季存在较强的黑潮海洋锋,利用数值实验研究斜坡地形和海洋锋同时存在时由浅海至深海的声传播特性。海洋模式数值预报环境数据表明,分布于冲绳海槽斜坡上方的海洋锋导致该海域上层水体声速存在水平变化,纬度越高,水平变化越大。利用抛物方程声场模型计算声传播损失,通过简正模态分析存在表面声道的环境中声能量分布,利用声线轨迹图解释海底斜坡和海洋锋对声传播的影响。结果表明:声源频率低于表面声道截止频率时,声传播主要受海底地形影响;声源频率高于表面声道截止频率时,位于表面声道内的声源激发的声能量主要在表面声道传播,部分声能量从表面声道泄漏沿斜坡向深海传播,位于表面声道深度以下的声源激发的声能量主要沿斜坡向深海传播,斜坡地形导致表面声道下方至共轭深度这一深度范围呈现为声影区,海洋锋的存在可导致表面声道传播损失变化明显,影响程度与声源深度有关。     

Computer modeling of sound propagation in the ocean with fractal inhomogeneities

The wave-field computer code based on the wide-angle parabolic equation is modified and adapted to the problems of sound scattering in a medium with anisotropic inhomogeneities of fractal type. To verify the computer code, a model numerical experiment on determining the angular dependence of the scattered sound field is performed for different anisotropy coefficients of the sound speed inhomogeneities. The comparison of the computed data with the theoretical dependences shows their rather good agreement and indicates that the computer code can be applied to calculations of sound propagation in the ocean with fine-structure inhomogeneities possessing fractal properties.     

Three-dimensional coupled-mode model and characteristics of low-frequency sound propagation in ocean waveguide with seamount topography

Large-scale topography, such as a seamount, substantially impacts low-frequency sound propagation in an ocean waveguide, limiting the application of low-frequency acoustic detecting techniques. A three-dimensional (3D) coupled-mode model is developed to calculate the acoustic field in an ocean waveguide with seamount topography and analyze the 3D effect. In this model, a correction is introduced in the bottom boundary, theoretically making the acoustic field satisfy the energy conservation. Furthermore, a large azimuth angle calculation range is obtained by using the operator theory and higher-order Padé approximation. Additionally, the model has advantages related to the coupling mode and parabolic equation theory. The couplings corresponding to the effects of range-dependent environment are fully considered, and the numerical implementation is kept feasible. After verifying the accuracy and reliability of the model, low-frequency sound propagation characteristics in the seamount environment are analyzed. The results indicate lateral variability in bathymetry can lead to out-of-plane effects such as the horizontal refraction phenomenon, while the coupling effect tends to restore the abnormal sound field and produces acoustic field diffraction behind the seamount. This model effectively considers the effects of the horizontal refraction and coupling, which are proportional to the scale of the seamount.     

Incorporating source directionality into outdoor sound propagation calculations

Many outdoor sound sources, such as aircraft or ground vehicles, exhibit directional radiation patterns. However, long-range sound propagation algorithms are usually formulated for omnidirectional point sources. This paper describes two methods for incorporating directional sources into long-range sound propagation algorithms. The first is the equivalent source method (ESM), which determines a compact distribution of omnidirectional point sources reproducing a given directivity pattern in the far field. This method can be used with any propagation algorithm because it explicitly reconstructs a source function as a set of point sources with certain amplitudes and positions. The second is a directional starter method (DSM), which is developed specifically for the parabolic equation (PE) algorithms. This method derives narrow- or wide-angle directional starter fields, corresponding to a given source directivity pattern, without reconstructing the equivalent source distribution. Although the ESM can also be used for the PE, the DSM is simpler and can be more convenient, especially if the sound propagation is calculated only for one or a few azimuthal directions. While these two methods are found to produce generally distinct starter fields, they nonetheless yield identical directivity patterns.     

深远海线性内波条件下声场时间相关性的空间分布

深海声场特定的干涉结构导致其时间相关性的空间起伏,研究这种空间特性可以为水声信号的探测与处理提供重要参考。利用抛物方程声场仿真模型,联合Monte-Carlo数值方法计算分析了深远海线性内波条件下声场时间相关性的空间分布特性。与现有的研究相比,给出了时间相关性的距离和深度起伏特征。结果表明,当接收达到一定距离,声场时间相关性的空间分布具有与声场干涉条纹类似的结构,声场干涉越强,时间相关性越好。此外,声源频率和声速标准差的变化会引起时间相关性空间分布规律的改变,且会聚区传播模式下的改变强于深海声道传播模式。     

Stochastic differential equation analysis for sound scattering by random internal waves in the ocean

The scattering of a weakly divergent narrow sound beam by random inhomogeneities of a fluctuating ocean is considered in the coupled-mode approximation. The random index of sound refraction is described using the Garrett-Munk internal wave spectrum. The problem is solved using the stochastic differential equations for the first-and second-order statistical moments of the acoustic field. The equations are formulated according to the cumulant expansion method. The existence of weakly divergent narrow sound beams in long-range sound propagation was one of the last discoveries of L.M. Brekhovskikh, to which he attached much importance. The concentration of sound into narrow beams away from the axis of the underwater sound channel was first observed experimentally and then explained by Brekhovskikh and his former students Goncharov, Kurtepov, and Petukhov. In the present paper, the scattered field intensity of a sound beam is calculated for different frequencies and source depths. Analytical expressions are obtained for the coefficients of the differential equation. The intermode energy transfer that accompanies the long-range propagation of a weakly divergent sound beam is analyzed. A comparison with the conventionally used Monte Carlo simulation in the parabolic equation approximation is performed.     

三维绝热简正波-抛物方程理论及应用

复杂海域通常存在环境参数的水平变化,这会导致声波在传播过程中发生水平折射,呈现出三维效应.利用绝热简正波-抛物方程理论进行三维声场建模,在垂直方向上使用标准简正波模型KRAKEN求解本征值和本征函数,水平方向上使用宽角抛物方程模型RAM求解简正波幅度.该模型物理意义清晰,计算效率高,但由于忽略了各号简正波之间的耦合,只适用于环境参数水平变化缓慢的问题.使用该模型分析了内波环境和大陆架楔形波导中的声波水平折射现象,结果表明,声波的水平折射将水平平面分为不同区域,每个区域内的声场结构明显不同.此外,声强在水平平面内的分布与声源频率和简正波号数有关,这种依赖关系是导致声信号频谱变化、波形畸变以及声场时空扰动的主要原因.     

Sound propagation and scattering in media with random inhomogeneities of sound speed, density and medium velocity

This review presents both classical and new results of the theory of sound propagation in media with random inhomogeneities of sound speed, density and medium velocity (mainly in the atmosphere and ocean). An equation for a sound wave in a moving inhomogeneous medium is presented, which has a wider range of applicability than those used before. Starting from this equation, the statistical characteristics of the sound field in a moving random medium are calculated using Born-approximation, ray, Rytov and parabolic-equation methods, and the theory of multiple scattering. The results obtained show, in particular, that certain equations previously widely used in the theory of sound propagation in moving random media must now be revised. The theory presented can be used not only to calculate the statistical characteristics of sound waves in the turbulent atmosphere or ocean but also to solve inverse problems and develop new remote-sensing methods. A number of practical problems of sound propagation in moving random media are listed and the further development of this field of acoustics is considered.     

A High-Efficiency Spectral Method for Two-Dimensional Ocean Acoustic Propagation Calculations

The accuracy and efficiency of sound field calculations highly concern issues of hydroacoustics. Recently, one-dimensional spectral methods have shown high-precision characteristics when solving the sound field but can solve only simplified models of underwater acoustic propagation, thus their application range is small. Therefore, it is necessary to directly calculate the two-dimensional Helmholtz equation of ocean acoustic propagation. Here, we use the Chebyshev–Galerkin and Chebyshev collocation methods to solve the two-dimensional Helmholtz model equation. Then, the Chebyshev collocation method is used to model ocean acoustic propagation because, unlike the Galerkin method, the collocation method does not need stringent boundary conditions. Compared with the mature Kraken program, the Chebyshev collocation method exhibits a higher numerical accuracy. However, the shortcoming of the collocation method is that the computational efficiency cannot satisfy the requirements of real-time applications due to the large number of calculations. Then, we implemented the parallel code of the collocation method, which could effectively improve calculation effectiveness.     

Sound propagation and scattering in media with random inhomogeneities of sound speed,density and medium velocity

Abstract

This review presents both classical and new results of the theory of sound propagation in media with random inhomogeneities of sound speed, density and medium velocity (mainly in the atmosphere and ocean). An equation for a sound wave in a moving inhomogeneous medium is presented, which has a wider range of applicability than those used before. Starting from this equation, the statistical characteristics of the sound field in a moving random medium are calculated using Born-approximation, ray, Rytov and parabolic-equation methods, and the theory of multiple scattering. The results obtained show, in particular, that certain equations previously widely used in the theory of sound propagation in moving random media must now be revised. The theory presented can be used not only to calculate the statistical characteristics of sound waves in the turbulent atmosphere or ocean but also to solve inverse problems and develop new remote-sensing methods. A number of practical problems of sound propagation in moving random media are listed and the further development of this field of acoustics is considered.     

Development of a hybrid computer model for simulating the complicated virtual sound field in enclosures

A hybrid computer model, SOFIS, has been developed for the simulation of an enclosed virtual sound field in an arbitrary shaped enclosure. It can be used to calculate the impulse response and acoustical parameters of different positions in a virtual enclosure. This paper describes the way in which SOFIS models the sound source, the receiver and the sound propagation throughout the enclosed space with curved surfaces or barriers. A phase tracing method and the calculation of acoustic indexes are also discussed in this paper.