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.     

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

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

Mid-frequency geoacoustic inversion using bottom loss data from the Shallow Water 2006 Experiment

Geoacoustic inversion work has typically been carried out at frequencies below 1 kHz, assuming flat, horizontally stratified bottom models. Despite the relevance to Navy sonar systems many of which operate at mid-frequencies (1-10 kHz), limited inversion work has been carried out in this frequency band. This paper is an effort to demonstrate the viability of geoacoustic inversion using bottom loss data between 2 and 5 kHz. The acoustic measurements were taken during the Shallow Water 2006 Experiment off the coast of New Jersey. A half-space bottom model, with three parameters density, compressional wave speed, and attenuation, was used for inversion by fitting the model to data in the least-square sense. Inverted sediment sound speed and attenuation were compared with direct measurements and with inversion results using different techniques carried out in SW06. Inverted results of the present work are consistent with other measurements, considering the known spatial variability in this area. The observations and modeling results demonstrate that forward scattering from topographical changes is important at mid-frequencies and should be taken into account in sound propagation predictions and geoacoustic inversion. To cope with fine-scale topographic variability, measurement technique such as averaging over tracks may be necessary.     

Geoacoustic Inversion Based on Dispersion Characteristic of Normal Modes in Shallow Water

A geoacoustic inversion method based on dispersion characteristic of normal modes is presented. An adaptive time-frequency analysis technique with a high resolution in both time and frequency domains is applied to derive the dispersion characteristic of normal modes from the broadband propagation signal. The bottom acoustic parameters are inverted by matching the calculated group delays of normal modes with the experimental data. Finally, some experimental results which could validate the inversion method are given.     

Data error covariance in matched-field geoacoustic inversion

Many approaches to geoacoustic inversion are based implicitly on the assumptions that data errors are Gaussian-distributed and spatially uncorrelated (i.e., have a diagonal covariance matrix). However, the latter assumption is often not valid due to theory errors, and can lead to reduced accuracy for geoacoustic parameter estimates and underestimation of parameter uncertainties. This paper examines the effects of data error (residual) covariance in matched-field geoacoustic inversion. An inversion approach is developed based on a nonparametric method of estimating the full covariance matrix (including off-diagonal terms) from the data residuals and explicitly including this covariance in the misfit function. Qualitative and quantitative statistical tests for Gaussianity and for correlations in complex residuals are considered to validate the inversion results. The approach is illustrated for Bayesian geoacoustic inversion of broadband, vertical-array acoustic data measured in the Mediterranean Sea.     

Optimized constraints for the linearized geoacoustic inverse problem

A geoacoustic inversion scheme to estimate the depth-dependent sound speed characteristics of the shallow-water waveguide is presented. The approach is based on the linearized perturbative technique developed by Rajan et al. [J. Acoust. Soc. Am. 82, 998-1017 (1987)]. This method is applied by assuming a background starting model for the environment that includes both the water column and the seabed. Typically, the water column properties are assumed to be known and held fixed in the inversion. Successful application of the perturbative inverse technique lies in handling issues of stability and uniqueness associated with solving a discrete ill-posed problem. Conventionally, such problems are regularized, a procedure which results in a smooth solution. Past applications of this inverse technique have been restricted to cases for which the water column sound speed profile was known and sound speed in the seabed could be approximated by a smooth profile. In this work, constraints that are better suited to specific aspects of the geoacoustic inverse problem are applied. These techniques expand on the original application of the perturbative inverse technique by including the water column sound speed profile in the solution and by allowing for discontinuities in the seabed sound speed profile.     

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.     

Bayesian geoacoustic inversion of ship noise on a horizontal array

This paper applies geoacoustic inversion to low-frequency narrow-band acoustic data from a quiet surface ship recorded on a bottom-moored horizontal line array in shallow water. A Bayesian matched-field inversion method is employed which quantifies geoacoustic uncertainties and allows for meaningful comparison of inversion results from different data sets. Geoacoustic inversion results for ship-noise data are compared with inversion results for multitone data from a towed controlled source collected in the same experiment, and with independent geophysical measurements. To increase the information content of low-level ship-noise data, the effect of including multiple, independent data segments in the inversion is investigated and shown to significantly reduce geoacoustic parameter uncertainties. Geoacoustic uncertainties are also shown to depend on ship range and orientation, with increased uncertainties for long ranges and for the ship stern oriented away from the array.     

Effects of incoherent and coherent source spectral information in geoacoustic inversion

This paper examines the effect on matched-field geoacoustic inversion of including source spectral information, as can be available in controlled-source acoustic surveys. Source information can consist of relative or absolute knowledge of the source amplitude and/or phase spectra, and can allow frequency-coherent processing of spatial acoustic-field data. A number of multi-frequency acoustic processors, appropriate for specific types of source information, are defined based on the likelihood function for complex acoustic-field data with Gaussian noise. The information content of the various processors is quantified in terms of marginal probability distributions and highest-probability density intervals for the unknown geoacoustic and geometric parameters, which define the accuracy expected in inversion. Marginal distributions are estimated using a fast Gibbs sampler approach to Bayesian inversion, which provides an efficient, unbiased sampling of the multi-dimensional posterior probability density. The analysis is illustrated for incoherent and coherent processors corresponding to several types of source knowledge ranging from complete information to no information, and the results are considered as a function of the spatial and frequency sampling of the acoustic fields.     

Statistical geoacoustic inversion from vertical correlation of shallow water reverberation

Statistical geoacoustic inversion results based on reverberation vertical correlation （RVC） data from 300 Hz to 800 Hz are presented. The data were obtained during Yellow Sea Experiment-2005 （YSE-05）. An uncertainty analysis is carried out. It is found that the inversion sea bottom sound velocities decrease when the frequency increases. So it is difficult to determine the sea bottom sound velocity. In order to solve this problem, a two-layer bottom model is assumed, and a multi-frequency inversion approach based on Genetic Algorithm is proposed. The approach is demonstrated using YSE-05 experiment data. Both RVC and normal modes depth functions （NMDFs） calculated using the inverted geoacoustic parameters are in good agreement with the experimental measurements.     

Subbottom profiling using a ship towed line array and geoacoustic inversion

Bottom profiling traditionally uses broadband signals received on a line array at long ranges to estimate the bottom layer structure and thickness. In this paper, a subbottom profiling method is developed and applied to a ship-towed line array using the same ship towed source to estimate the subbottom layer structure and thickness. A ship towed line-array system can be used to estimate bottom properties using geoacoustic inversion and can cover a wide area in a short time. It needs some prior information about the subbottom structure and layer thickness, without which the solution can be ambiguous and even erratic when resolving parameters over a wide area. It is shown that the required subbottom information can be obtained from the time-angle relation by beamforming the same acoustic signal data used for geoacoustic inversion. The time-angle analysis is used to expose the prevalent physics intrinsic to geoacoustic inversion. One finds that the tau-p relation of the bottom and the bottom reflection coefficients, sampled at discrete angles associated with bottom and multiple surface-bottom returns, are often adequate, for many practical applications, to uniquely determine the geoacoustic bottom at low (< or =1 kHz) frequencies.     

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.     

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.     

A phase regulated back wave propagation technique for geoacoustic inversion

An inversion method based on the concept of back wave propagation (BWP) is described in this paper for estimation of geoacoustic parameters from acoustic field data. A phase-regulation technique is introduced to increase the sensitivity of the method for geoacoustic model parameters having low sensitivity. The case of data consisting of signal plus additive noise is also addressed. It is shown theoretically that the sensitivity can be increased by a factor alpha using the phase regulation procedure, and that the spatial resolution of signal energy that is concentrated by BWP at the known source position is increased when a increases. This result suggests an effective criterion for use in the inversion, based on the spatial distribution of signal energy around the true source location. The basis for the criterion is the spatial variance of the back-propagated pressure field in a window around the known source location. A multistep search process is proposed to avoid using a complicated multidimensional search process. Inversion results from both simulations and experimental data are given. The real data were taken from the Pacific Shelf experiment carried out in shallow water off the West Coast of Vancouver Island in the Northeast Pacific Ocean.     

负跃层浅海声传播问题的解析解及其在水下脉冲声传播研究中的应用

提出一个基于波数积分理论的解析解,可用于计算包含海面下均匀层、负跃层、海底上方均匀层以及液态均匀海底的负跃层浅海波导中的声场。对声源位于水体不同层中的情况给出详细的推导,给出深度格林函数的解析解,并介绍了如何利用其计算点源和线源这两种情况下的声场.通过一个典型的负跃层浅海算例,验证了本文所提方法的精度,以及处理点源及线源问题的能力.此外,将本文提出的方法应用于水下脉冲声传播研究。针对某次实验的实际海洋环境,利用本文方法研究接收到的脉冲信号的波形特点,并与利用射线理论分析得到的相邻脉冲到达时间间隔等理论结果进行比较·结果表明,利用本方法得到的相邻脉冲到达时间数值结果与理论值高度一致。本方法非常适用于需要大量重复计算声场的情况,如水下脉冲声传播、海底声学参数反演等问题的计算机仿真研究。      

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.     

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.     

地声参数及传播损失不确定性估计与建模

地声参数的不确定性对水声传播具有重要的影响。通过贝叶斯理论建立水声环境不确定性推理模型,理论推导了地声参数的似然函数以及地声参数和传播损失的后验概率密度,并采用MCMC(Markov Chain Monte Carlo)进行了仿真计算,给出了地声参数的二维后验联合概率密度和一维边缘概率密度,在此基础上对传播损失的不确定性进行了估计,得到了传播损失80%的可信区间。仿真和实验结果表明,该方法适用于地声参数反演和不确定性估计,并能获取因地声参数不确定性导致的传播损失不确定性估计。     

Quantifying the uncertainty of geoacoustic parameter estimates for the New Jersey shelf by inverting air gun data

This paper describes geoacoustic inversion of low frequency air gun data acquired during an experiment on the New Jersey shelf. Hybrid optimization and Bayesian inversion techniques based on matched field processing were applied to multiple shots from three air gun data sets recorded by a vertical line array in a long-range shallow water geometry. For the Bayesian inversions, full data error covariance matrix was estimated from a set of consecutive shots that had high temporal coherence and small spatial variation in source position. The effect of different data error information on the geoacoustic parameter uncertainty estimates was investigated by using the full data error covariance matrix, a diagonalized version of the full error covariance, and a diagonal matrix with identical variances. The comparison demonstrated that inversion using the full data error information provided the most reliable parameter uncertainty estimates. The inversions were highly sensitive to the near sea floor geoacoustic parameters, including sediment attenuation, of a simple single-layer geoacoustic model. The estimated parameter values of the model were consistent with depth averaged values (over wavelength scales) of a high resolution geoacoustic model developed from extensive ground truth information. The interpretation of the frequency dependence of the estimated attenuation is also discussed.     

Tabu for matched-field source localization and geoacoustic inversion

Tabu is a global optimization technique that has been very successful in operations research. In this paper, a Tabu-based method is developed for source localization and geoacoustic inversion with underwater sound data; the method relies on memory to guide the multiparameter search. Tabu is evaluated through a comparison to simulating annealing. Both methods are tested by inverting synthetic data for various numbers of unknown parameters. Tabu is found to be superior to the simulated annealing variant implemented here in terms both of accuracy and efficiency. Inversion results from the SWellEX-96 data set are also presented.