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

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

利用海底反射信号进行地声参数反演的方法

针对现有反演方法的缺点,提出了一种基于海底反射信号的地声参数高分辨反演方法.它利用短距离声源在不同深度上发射宽带线性调频信号,采用垂直阵进行接收,首先通过匹配滤波方法提取多径到达信息,然后利用海底反射损失曲线,反演海底表层的声速和密度,最后利用浅底层反射信号估计沉积层参数.由于海水中直达波受到内波的强烈影响,选择海底表面反射作为参考,用以可靠地计算浅底层反射的相对到达时间和幅度,从而估计出沉积层的厚度、速度和衰减系数.通过海上实验,验证了利用浅底层反射信号反演参数的有效性. 关键词： 海底参数 反演 浅底层反射信号     

利用拖船辐射噪声的深海声学参数近场匹配场反演方法

针对深海声学参数难以通过远距离合作声源反演获取的问题,提出了利用拖船低频噪声近场匹配场反演方法。首先,利用聚焦波束形成计算拖曳阵接收拖船噪声的方向性,获得传播路径特征;然后,构建多参数反演模型,由波数积分声传播模型计算拷贝场,采用遗传算法对多频匹配场目标函数进行反演。同时,采用蒙特卡罗方法分析参数后验概率密度。仿真与试验结果表明:深海环境中拖曳阵接收拖船噪声主要来自海底反射路径,利用该特性反演得到海水深度、噪声源距离、阵列深度、沉积层厚度等参数,多频联合反演可以提高沉积层厚度等参数反演准确性。宽带匹配场处理表明,利用反演最优参数模型能准确给出拖船噪声源的空间位置。      

利用海洋环境噪声估计浅海海底表层声速

为提高反演效率,提出一种快速估计浅海海底表层声速的方法。根据噪声能流理论,垂直阵接收的海洋环境噪声数据能够用于无源提取海底反射损失,反射损失曲线中具有明显的临界角效应,从而接估计海底表层声速。以射线模型为基础,推导了噪声提取的反射损失与理论值之间的差异,并讨论实际阵列波束形成在不同角度和频率下的性能。考虑到非等声速环境下声线会发生弯曲,需对角度进行修正以提高方法的广泛适用性。不同频率的临界角与有效深度之间存在对应关系,黄海某海区数据处理结果表明,在临界角不变的有效深度内,海底表层可以视为等声速层,该海区海底表层0.5 m内声速估计结果为1547 m/s,与有源反演结果相近。      

Geoacoustic inversion for bottom parameters in a thermocline environment in the northern area of the South China Sea

Bottom acoustic parameters have important influence on the application of underwater acoustic propagation and source location.The acoustic parameters of the seabed in the northern of the South China Sea(SCS) were inversed using the experiment data from an acoustic experiment in 2015.Based on the comprehensive analysis of the influence of the sound speed fluctuation and the geoacoustic model on seabed inversion,the multi-parameter hybrid acoustic inversion scheme is improved by selecting the equivalent mean sound speed profile(SSP) and half-infinite liquid bottom model to save the inversion dimensions in the matched field processing(MFP) inversion.The inverted bottom sound speed and density are in good agreement with the core sampling measurements.The nonlinear empirical relationship of the attenuation coefficient with frequency is given out.The inversion results are meaningful to the sound propagation research and application in the northern area of the SCS.     

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.     

Broadband sound propagation in shallow water and geoacoustic inversion

Part of an experiment to test a measurement package in a shallow water region in the Gulf of Mexico was designed to gather broadband acoustic data suitable for inversion to estimate seabed geoacoustic parameters. Continuous wave tow acoustic signals at multiple frequencies and broadband impulsive source signals were recorded on a horizontal line array in a high-noise environment. Simulated annealing with a normal mode forward propagation model is utilized to invert for a geoacoustic representation of the seabed. Several inversions are made from different data samples of two light bulb implosions, the measured sound speed profiles at the HLA and at the positions of the light bulb deployments, and for two different cost functions. The different cost functions, measured sound speed profiles, and measured time series result in different inverted geoacoustic profiles from which transmission loss is generated for comparison with measurements. On the basis of physical consistency and from the comparison of the transmission loss and time series, a best estimate geoacoustic profile is selected and compared to those obtained from previously reported inversions. Uncertainties in the sound speed profile are shown to affect the uncertainties of the estimated seabed parameters.     

Geoacoustic inversion with ships as sources

Estimation of geoacoustic parameters using acoustic data from a surface ship was performed for a shallow water region in the Gulf of Mexico. The data were recorded from hydrophones in a bottom mounted, horizontal line array (HLA). The techniques developed to produce the geoacoustic inversion are described, and an efficient method for geoacoustic inversion with broadband beam cross-spectral data is demonstrated. The performance of cost functions that involve coherent or incoherent sums over frequency and one or multiple time segments is discussed. Successful inversions for the first sediment layer sound speed and thickness and some of the parameters for the deeper layers were obtained with the surface ship at nominal ranges of 20, 30, or 50 water depths. The data for these inversions were beam cross-spectra from four subapertures of the HLA spanning a little more than two water depths. The subaperture beams included ten frequencies equally spaced in the 120-200 Hz band. The values of the geoacoustic parameters from the inversions are validated by comparisons with geophysical observations and with the parameter values from previous inversions by other invesigators, and by comparing transmission loss (TL) measured in the experiment with modeled TL based on the inverted geoacoustic parameters.     

Geoacoustic Inversion Based on a Vector Hydrophone Array

We propose a geoacoustic inversion scheme employing a vector hydrophone array based on the fact that vector hydrophone can provide more acoustic field information than traditional pressure hydrophones. Firstly, the transmission loss of particle velocities is discussed. Secondly, the sediment sound speed is acquired by a matchedfield processing （MFP） procedure, which is the optimization in combination of the pressure field and vertical particle velocity field. Finally, the bottom attenuation is estimated from the transmission loss difference between the vertical particle velocity and the pressure. The inversion method based on the vector hydrophone array mainly has two advantages： One is that the MFP method based on vector field can decrease the uncertain estimation of the sediment sound speed. The other is that the objective function based on the transmission loss difference has good sensitivity to the sediment attenuation and the inverted sediment attenuation is independent of source level. The validity of the inverted parameters is examined by comparison of the numerical results with the experimental data.     

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.     

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

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

深海近距离声场频率-距离干涉结构反演海底声学参数

地声反演一般存在多解性的问题,不易有效分离沉积层声速与厚度,本文提出了一种当声源与接收水听器分别位于海面与海底附近时匹配声场干涉周期反演深海沉积层声速与厚度的方法。该声场周期结构由直达波与海底沉积层反射回波干涉形成,其干涉周期受沉积层声速和厚度影响,通过匹配实测与理论计算的声场干涉周期可有效反演沉积层声速和厚度。开展了深海实验验证,利用水听器接收到的近距离船舶噪声的声场干涉条纹,反演得到的实验海域海底沉积层声速与该海域海底表层采样平均值差小于20 m/s。数值仿真与实验结果表明该方法可以有效反演沉积层的声速与厚度。      

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

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

Range-dependent seabed characterization by inversion of acoustic data from a towed receiver array

The MAPEX2000 experiments were conducted in the Mediterranean Sea in March, 2000 to determine seabed properties using a towed acoustic source and receiver array. Towed systems are advantageous because they are easy to deploy from a ship and the moving platform offers the possibility for estimating spatially variable (range-dependent) seabed properties. In this paper, seabed parameters are determined using a matched-field geoacoustic inversion approach with measured, towed array data. Previous research has successfully applied matched-field geoacoustic inversion techniques to measured acoustic data. However, in nearly all cases the inverted data were collected on moored, vertical receiver arrays. Results here show that seabed parameters can also be extracted by inverting acoustic measurements from a towed array of receivers, and these agree with those inverted using data received simultaneously on a vertical array. These findings imply that a practical technique could be developed to map range-dependent seabed parameters over large areas using a towed acoustic system. An example of such a range-dependent inversion is given using measurements from the MAPEX2000 experiments.     

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

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

利用气枪声源数据的地声参数反演*

本文通过对东海采集的气枪声源数据进行处理,基于地声参数对不同声场观测量的敏感性差异,利用分步反演的策略进行分析。首先选择距离接收阵七个不同的距离点数据,通过warping变换准确提取第三阶至第八阶简正波的频散曲线;利用海底声速对频散曲线敏感的特性来反演声速;由Hamilton经验公式求得海底密度;通过传播损失拟合获得海底衰减。不同距离点数据反演的海底声速与密度一致性较高。实验提取的频散曲线和和反演参数仿真结果、实验获取的传播损失与反演参数获得传播损失均拟合较好。     

不确定海洋声场中的检测性能损失环境敏感度度量

现有的检测算法在实际的不确定海洋环境中会出现失配情况, 进而导致检测性能下降, 但是定量分析这种检测性能损失的工作却很少见, 因此本文定义了检测性能损失敏感度函数, 并给出了基于蒙特卡罗方法的计算方法.检测性能损失敏感度是一个表征海洋环境不确定度的物理量, 它反映了海洋参数变动和检测性能损失之间的关系.利用地中海某处海洋环境进行仿真, 针对各个海洋环境参数, 计算了全海域的检测性能敏感度. 结果表明: 1)检测性能损失呈现了很强的空间性, 在20 km处, 水面声速有4 m/s的变化下, 检测损失最小为1%, 而最大达到60%.声信道中的检测性能较为稳定, 其在远距离上更是如此; 2)各个海洋环境参数对检测性能损失有不同的影响, 水体声速剖面和第一层海底介质的声速和厚度是对检测性能影响最大的量; 3)环境参数敏感度呈现很强的频率特性, 海底底质参数包括底质厚度、密度和吸收系数等随着频率的升高对检测性能的损失影响变小. 关键词： 检测性能损失 环境敏感度 不确定海洋环境     

Effects of sediment acoustic characteristics on low-frequency acoustic field distribution in shallow water

Most of the published literatures on low-frequency underwater sound propagation are focused on the sound propagation features in the water column,while studies on sound propagation features in the sediment layer or the semi-infinite basement are rare.In this paper,based on the wave equation,a computational model for sound energy flux in the sediment layer and the basement as well as in the water column is proposed under a cylindrical coordinate system.On this basis,the effects of various sediment acoustic parameters on the sound energy distribution and the corresponding mechanisms are elaborated through numerical examples and acoustic theory.Simulation results reveal that,in a situation where sediment P-wave speed>water sound speed> sediment S-wave speed,the greater the values of density and P-wave speed in sediment,the more likely it is that the sound energy remains in the water column without leaking to the sea floor.Conversely,the influence of the variation of S-wave speed is reversed.Basement influence on the sound propagation in the fluid layer is approximately negligible if the sediment layer is sufficiently thick.     

基于拖曳倾斜线列阵的海底反射损失提取方法

在浅海环境中,海底声学参数对水下声场的精确预报十分重要.现有的海底声学参数反演方法大多数是采用固定垂直阵进行的,其缺点是不能实施大面积、高分辨的走航式反演.在已提出的垂直阵海底反射损失反演方法的基础上,研究了基于三种线列阵阵形的海底反射损失提取方法及其特点,提出了基于拖曳倾斜线列阵与三个声源组合的走航式海底参数快速获取方法,并对该方法进行了误差分析.研究结果表明：覆盖同样的掠射角范围,垂直线列阵需与多个距离的声源组合,拖曳水平线列阵只需一个声源组合但需要采用较大的物理孔径,拖曳倾斜线列阵综合了垂直阵和水平 关键词： 水下声场预报 海底声学参数 海底反射损失 拖曳倾斜线列阵     

利用拖船自噪声进行浅海环境参数贝叶斯反演

研究了以拖船自噪声为参考声源的浅海环境参数反演问题,并针对反演结果不确定性快速量化评估问题,提出了一种基于自适应重要性抽样的贝叶斯反演新方法。反演利用了拖船自噪声低频线谱成分,并采用混合高斯推荐函数自适应推荐声场模型样本,使得样本集中于参数高概率密度区域,实现后验概率密度快速收敛计算。仿真试验结果表明:拖船自噪声反演能够准确估计水深、沉积层及阵列参数等。所提自适应重要性抽样贝叶斯反演方法的计算效率优于快速吉布斯抽样方法。利用试验数据处理验证,反演得到试验海域声学环境参数,计算传播损失与各阵元接收线谱强度变化吻合,说明反演最优环境模型能准确表征声场传播特征。