一种低声速沉积层海底参数声学反演方法

软泥底环境下沉积层参数的声学反演是国际水声领域的一个研究热点.浅海中,当高声速基底和海水之间存在一层低声速(小于海水声速)的沉积层时,小掠射角情况下不同频率声传播损失会出现周期性增大现象.基于此现象,提出一种适用于低声速沉积层的海底参数声学反演方法.首先,推导给出小掠射角情况下传播损失周期增大的频率间隔与沉积层声速、厚度及近海底海水声速之间的解析表达式;其次,利用一次黄海实验中软泥底环境下的宽带声传播信号,提取了小掠射角下传播损失增大的频率周期;再次,把该解析表达式作为约束条件,结合Hamilton密度与声速的经验公式,采用匹配场处理反演给出沉积层的声速、密度、厚度及基底的声速、密度;然后,利用声传播损失数据反演得到泥底环境下不同频率的声衰减系数,通过拟合发现泥底声衰减系数随频率近似呈线性关系;最后,给出了双层海底模型和半无限大海底模型等效性的讨论.反演结果为低声速沉积层海底声传播规律研究与应用提供了海底声学参数.     

采用小生境遗传算法反演浅海声速剖面研究*

针对浅海声速剖面反演问题,采用小生境遗传算法,结合声线搜索的最快本征声线匹配反演,实现浅海负跃层条件下的声速剖面估计。利用经验正交函数对声速剖面的多参数不确定性降维,依据声场计算模型获取的最快特征声线传播时延与观测声传播时长进行匹配,采用小生境遗传优化算法,获取最优经验正交函数估计,实现声速剖面反演。按上述方法反演处理浅海声传播实验数据,结果表明,该方法针能够有效反演浅海声速剖面,并且显著优于传统遗传算法反演结果。     

南海北部负跃层环境下海底参数声学反演

利用2015年南海北部声学实验的实验数据,对海底声学参数进行反演。在综合分析声速起伏及海底模型对海底反演影响的基础上,通过选择等效的海水声速剖面和海底模型,改进了多参量联合地声反演方法,使得其不仅能解决反演中的多值性问题,还能适用于负跃层起伏环境下的声学反演。即当用匹配场反演海底声速和用传播损失反演海底衰减时,如果温跃层内有内波等随机起伏,可使用传播路径上平均声速剖面和单层等效海底进行声场计算。反演得到的海底声速和密度结果与海底采样测量符合较好,拟合给出了海底衰减系数随频率的非线性经验关系式。反演结果可为南海北部声传播规律研究与应用提供海底参数。      

浅海环境噪声垂直指向性探测海底参数

浅海环境噪声垂直指向性和海底地声参数密切相关,可以利用它来反演海底参数。本文在水平分层介质噪声场模型下,讨论了海底参数对噪声垂直指向性的敏感度,并假设半无限大液态海底,采用效率较高的自适应单纯形模拟退火算法全局寻优,根据中国某海区实测噪声数据,估计了海底声速、密度及衰减,分析了反演结果的不确定性。鉴于声速和密度在不同频率的反演结果不够一致,文中提出多频海底参数联合反演方法,提高了反演海底声学参数的精度。     

沉积层声学特性对浅海低频声场分布的影响

已有对浅海低频水声场的讨论多是关注声能量在水体中的分布特性,对水体下沉积层、基底中低频声传播的同步研究相对较少。本文基于波动方程,在柱坐标系下推导了一种浅海水体/海底统一波导下低频声能流的计算模型,在此基础上结合具体仿真算例与波动理论阐述了不同沉积层声学参数对声场能量分布的影响规律及机理。仿真结果表明,在沉积层纵波声速>水中声速>沉积层横波声速的前提下,沉积层中密度与纵波声速数值越大,声能量越趋于保留在水体中而不向海底泄漏,横波声速的影响正好相反;沉积层厚度增加到一定量后,基底对流体层中声传播的影响可近似忽略不计。      

一种基于内潮动力特征的浅海声速剖面构建新方法

为了降低反演参数空间的维数,常利用正交经验函数(EOF)来构建声速剖面.然而,EOF方法的样本依赖性使之难以用于缺乏现场实测数据的海域.本文提出一种全新的利用历史数据而不依靠现场实时数据即可获得的声速剖面展开基函数.基于水质子流体静力方程和物态方程,推导了在缺乏实时测量的情况下从历史数据获得水动力模式基函数(HMB)的办法.利用WOA13季节平均温盐数据获得代表内潮动力特征的HMB进行分析.较之EOF, HMB及其对应的投影系数与海洋动力特征直接相关并具有明确的物理含义.基于东中国海实验获得的CTD (conductance-temperature-depth)及宽带爆炸声源声信号数据,利用声速剖面重构以及匹配场声层析对HMB进行了分析,并与EOF进行对比研究.结果表明:HMB可以以较好的精度构建浅海声速剖面.在对现场实时测量依赖更小的情况下,基于HMB方法的声场预报及声层析结果与EOF方法一样好.HMB的获取更简单且直接关联海水的物理特性,该方法可在实时测量样本不足的海域有效替代EOF进行海洋动力现象的声学监测.     

浅海复杂环境下等效声速剖面的构建方法*

针对浅海复杂环境声速剖面水平变化情况下的声传播损失预报及目标定位问题,提出了一种基于遗传算法的等效声速剖面重构算法。首先,将声速剖面进行时间和空间上的分解,从而将声速剖面抽象为对声速剖面前三阶正交函数系数的反演;其次,利用遗传算法,以先验声速剖面集为基础,进行参数反演。仿真结果表明,在浅海复杂条件下,传播损失预报受声速剖面及海洋参数的影响,不能直接运用接收或发射位置处的声速剖面进行传播损失预报,否则会对预报结果造成误差。通过构建具有声传播累积效应的等效声速剖面可以提高匹配场定位精度,完成目标定位,且在构建等效声速剖面时,接收位置处即本地声速剖面所占权重较大。     

自回归状态空间模型下时变声速剖面跟踪方法

讨论了一种适用于浅海的时变声速剖面跟踪方法。该方法将时变水体声速剖面的反演问题建模为由描述声速剖面时变特性的状态方程与包含声压场局部测量信息的测量方程组成的状态-空间模型,提出利用自回归分析拟合方法将声速场扰动建模为高阶自回归演化模型,并通过集合卡尔曼滤波序贯地估计时间演化的海洋声速场。利用2001年亚洲海实验环境与声速测量数据,仿真分析了基于高阶自回归演化模型的时变声速剖面集合卡尔曼滤波估计方法。结果表明,相比于利用传统随机游走状态演化模型的估计方法,该改进方法可有效降低声速的跟踪误差,并且在较低信噪比条件下仍具有较好的跟踪性能。      

三维声线追踪的正三棱锥前向伸展算法

提出一种在三维空间中求取声路径的高速声线追踪算法。该算法在声线前进方向的局部微小空间作一小正三棱锥,以使其内部的声速场可近似为线性;由于在线性声速场中声线轨迹总是在一个平面内,所以可使用二维空间中的三角形前向伸展算法确定声路径。完成了该算法的理论推导,并通过数值模型对其进行了检验。实验数据显示该算法有较高的运算速度和精确度,比传统三维算法更加简洁可靠,可望应用于三维超声层析像重建的散射修正。     

舰船辐射噪声贝叶斯方法反演浅海底层结构及地声参数EI北大核心CSCD

针对以舰船辐射噪声为参考声源的浅海海底分层结构及地声参数反演问题,研究了一种基于贝叶斯理论的浅海多层海底地声参数反演方法。反演中以舰船辐射噪声的线谱成分为研究对象,进而采用非线性贝叶斯反演方法反演浅海底层结构、层中声速、声速衰减和密度,并对反演结果的不确定性进行分析。反演结果的最大后验概率估计值和边缘概率分布分别通过拨正模拟退火算法和Metropolis-Hastings采样法在各参数先验区间内计算获得,并根据贝叶斯信息准则确定最佳海底分层结构。海上实验表明:根据该方法反演获得海底分层结构及地声参数,计算得到的声压场与实测舰船辐射噪声传播损失误差不超过10%,反演结果能够准确表征实验海区海底特征。反演结果不确定性分析表明:海底纵波声速、横波声速以及密度的不确定性更小,对声压场变化更加敏感,反演结果更有效、准确。     

Tracking of time-evolving sound speed profiles with an auto-regressive state-space model

An approach for time-evolving sound speed profiles tracking in shallow water is discussed.The inversion of time-evolving sound speed profiles is modeled as a state-space estimation problem,which includes a state equation for predicting the time-evolving sound speed profile and a measurement equation for incorporating local acoustic measurements.In the paper,auto-regression(AR) method is introduced to obtain a high-order AR evolution model of the sound speed field time variations,and the ensemble Kalman filter is utilized to track the sound speed field.To validate the approach,the accuracy in sound speed estimation is analyzed via a numerical implementation using the ASIAEX experimental environment and the sound velocity measurement data.Compared with traditional approaches based on the state evolution represented as a random walk,simulation results show the proposed AR method can effectively reduce the tracking errors of sound speed,and still keep good tracking performance at low signal-to-noise ratios.     

The impact of water column variability on horizontal wave number estimation and mode based geoacoustic inversion results

The influence of water column variability on low-frequency, shallow water geoacoustic inversion results is considered. The data are estimates of modal eigenvalues obtained from measurements of a point source acoustic field using a horizontal aperture array in the water column. The inversion algorithm is based on perturbations to a background waveguide model with seabed properties consistent with the measured eigenvalues. Water column properties in the background model are assumed to be known, as would be obtained from conductivity, temperature, and depth measurements. The scope of this work in addressing the impact of water column variability on inversion is twofold. Range-dependent propagation effects as they pertain to eigenvalue estimation are first considered. It is shown that mode coupling is important even for weak internal waves and can enhance modal eigenvalue estimates. Second, the effect of the choice of background sound speed profile in the water column is considered for its impact on the estimated bottom acoustic properties. It is shown that a range-averaged sound velocity profile yields the best geoacoustic parameter estimates.     

Study of a novel range-dependent propagation effect with application to the axial injection of signals from the Kaneohe source

A novel range-dependent propagation effect occurs when a source is placed on the seafloor in shallow water with a downward refracting sound speed profile, and sound waves propagate down a slope into deep water. Under these conditions, small grazing-angle sound waves slide along the bottom downward and outward from the source until they reach the depth of the sound channel axis in deep water, where they are detached from the sloping bottom and continue to propagate outward near the sound channel axis. This "mudslide" effect is one of a few robust and predictable acoustic propagation effects that occur in range-dependent ocean environments. As a consequence of this effect, a bottom mounted source in shallow water can inject a significant amount of acoustic energy into the axis of the deep ocean sound channel that can then propagate to very long ranges. Numerical simulations with a full-wave range-dependent acoustic model show that the Kaneohe experiment had the appropriate source, bathymetry, and sound speed profiles that allows this effect to operate efficiently. This supports the interpretation that some of the near-axial acoustic signals, received near the coast of California from the bottom mounted source located in shallow water in Kaneohe Bay, Oahu, Hawaii, were injected into the sound channel of the deep Pacific Ocean by this mechanism. Numerical simulations suggest that the mudslide effect is robust.     

The effects of ignored seabed variability in geoacoustic inversion

Geoacoustic inversion using a matched-field inversion algorithm is a well-established technique for estimating the geoacoustic parameters of the seabed. This paper demonstrates how parameter estimation can be affected by unknown or wishfully ignored random range dependence of the true environment when the inversion model is--for practicality--assumed to be range independent. Simulations with controlled statistics were carried out using a simple shallow water model: an isospeed water column over a homogeneous elastic halfspace. The inversion parameters included water depth, compressional speed in the seabed, seabed density, and compressional wave attenuation. On average the environment is range independent: some parameters are constant while other parameters are random with range-independent means and variances. A Parabolic Equation underwater acoustic propagation model is used to calculate the simulated data fields for the range-dependent environment as well as to calculate the model fields for the range-independent inversion model. The Adaptive Simplex Simulated Annealing inversion algorithm is used to estimate the best-fit solution. It is found that ignoring the variability of even a single geoacoustic parameter leads to significant and correlated uncertainty (bias and variance) in the estimation of all inverted parameters. Results are presented for range variation of compressional sound speed and water depth.     

浅海波导中目标散射的简正波-Kirchhoff近似混合方法

提出了计算浅海波导中复杂目标散射的数值方法:简正波-Kirchhoff近似混合方法。通过把目标散射的Kirchhoff近似方法和简正波声传播模型相结合,可对大尺寸复杂目标在浅海波导中的散射声场进行计算。以浅海波导中刚性球体散射的解析解为标准解验证了本方法,说明简正波-Kirchhoff近似混合方法是有一定计算精度的工程预报方法。数值计算Pekeris浅海波导中球体目标散射声场在深度-频率平面和距离-频率平面上的二维干涉结构及其与自由空间中的差异。进而通过FFT获得目标时域回波随深度的分布图,分析浅海波导中目标姿态、声速剖面对Benchmark潜艇目标散射的影响。      

Sequential Parameter Estimation Using Modal Dispersion Curves in Shallow Water

Existing sequential parameter estimation methods use the acoustic pressure of a line array as observations. The modal dispersion curves are employed to estimate the sound speed profile(SSP) and geoacoustic parameters based on the ensemble Kalman filter. The warping transform is implemented to the signals received by a single hydrophone to obtain the dispersion curves. The experimental data are collected at a range-independent shallow water site in the South China Sea. The results indicate that the SSPs are well estimated and the geoacoustic parameters are also well determined. Comparisons of the observed and estimated modal dispersion curves show good agreement.     

浅海时变声速环境下的自适应匹配场定位算法实现*

针对浅海环境下声速剖面失配引起的匹配场处理器失配问题,提出了一种自适应匹配场定位算法在声速剖面时变环境下的实现方式。将先验声速剖面集简化为经验正交函数表示,结合蒙特卡洛方法与环境扰动约束算法对当下时刻的目标声源进行匹配场定位。本文以某次试验获取的连续20小时的声速剖面数据为研究对象,通过仿真试验对该算法进行验证,结果表明：在先验声速剖面集的半小时之后,利用自适应算法的距离和深度定位成功率较常规匹配场算法有较大提升,其中,深度正确定位概率相对较低。