Broad-band geoacoustic inversion in shallow water from waveguideimpulse response measurements on a single hydrophone: theory andexperimental results

The paper discusses an inversion method that allows the rapid determination of in situ geoacoustic properties of the ocean bottom without resorting to large acoustic receiving apertures, synthetic or real. The method is based on broad-band waterborne measurements and modeling of the waveguide impulse response between a controlled source and a single hydrophone. Results from Yellow Shark '94 experiments in Mediterranean shallow waters using single elements of a vertical array are reviewed, inversion of the bottom parameters is performed with an objective function that includes the processing gain of a model-based matched filter (MBMF) receiver relative to the conventional matched filter. The MBMF reference signals incorporate waveguide Green's functions for known geometry and water column acoustic model and hypothesized bottom geoacoustic models. The experimental inversion results demonstrated that, even for complex environmental conditions, a single transmission of a broad-band (200-800 Hz) coded signal received at a single depth and a few hundred forward modeling runs were sufficient to correctly resolve the bottom features. These included the sound speed profile, attenuation, density, and thickness of the top clay sediment layer, and sound speed and attenuation of the silty clay bottom. Exhaustive parameter search proved unequivocally the low-ambiguity and high-resolution properties of the MBMF-derived objective. The single-hydrophone results compare well with those obtained under identical conditions from matched-field processing of multitone pressure fields sampled on the vertical array. Both of these results agree with expectations from geophysical ground truth. The MBMF has been applied successfully to a field of advanced drifting acoustic buoys on the Western Sicilian shelf, demonstrating the general applicability of the inversion method presented     

A time series analysis of sound propagation in a strongly multipathshallow water environment with an adiabatic normal mode approach

Measured time series were generated by small omnidirectional explosive sources in a shallow water area. A bottom-mounted hydrophone recorded sound signals that propagated over a sloping bottom. The time series in the 250-500 Hz band were analyzed with a broad-band adiabatic normal mode approach. The measured waveforms contain numerous bottom interacting multipaths that are complicated by the subbottom structure that contains high-velocity layers near the water-sediment interface. Several of the details of the geoacoustic structure and the depth of the water column at the receiver are inferred from comparisons of the measured data to simulated time series. The sensitivity of broad-band matched-field ambiguity surfaces in the range-depth plane for a single receiver to selected waveguide parameters is examined. A consistent analysis is made where the simulated time series are compared to the measured time series along with the single-receiver matched-field localization solutions for ranges out to 5 km. In this range interval, it was found that the peak cross-correlation between the measured and simulated time series varied between 0.84 and 0.69. The difference between the GPS range and the range obtained from the matched-field solution varied from 0 to 63 m. The geoacoustic structure obtained in the analysis consists of an 8-m low-velocity sediment layer over an 8-m high-velocity layer followed by a higher velocity, infinite half-space     

Matched-field inversion for geoacoustic model parameters in shallowwater

Matched-field inversion is used to, estimate geoacoustic properties from data obtained in an experiment with a vertical line array (VLA). The experiment was carried out using broad-band sources (shots) in water depths of about 200 m on the continental shelf off Vancouver Island. The data were processed to obtain spectral components of the field for frequencies near the bubble frequency for the shot. The ocean bottom in this region consists of a layer of mainly sandy sediments (about 100 m thick) overlying older consolidated material. Consequently, the inversion was designed to estimate the parameters of a two-layer elastic sediment model. In the inversion, an adaptive global search algorithm was used to investigate the multidimensional space of geoacoustic models in order to determine the set of values corresponding to the best replica field. Convergence is driven by adaptively guiding the search to regions of the parameter space associated with above-average values of the matched field correlation between the measured and replica fields. The geoacoustic profile estimated by the inversion consisted of a 125-m layer with compressional speed ~1700 m/s and shear speed ~400 m/s, overlying a layer with compressional speed ~1900 m/s. This model is consistent with the results from conventional seismic experiments carried out in the same region     

Geoacoustic seafloor exploration with a towed array in a shallowwater area of the Strait of Sicily

Acoustic propagation in shallow water is greatly dependent on the geoacoustic properties of the seabottom. This paper exploits this dependence for estimating geoacoustic sediment properties from the bottom acoustic returns of known signals received on a hydrophone line array. There are two major issues in this approach: one is the feasibility of acoustic inversion with a limited aperture line array, the other is related to the knowledge of the geometry of the experimental configuration. To test the feasibility of this approach, a 40-hydrophone-4-m spaced towed array together with a low-frequency acoustic source, was operated at a shallow water site in the Strait of Sicily. In order to estimate the array deformation in real time, it has been equipped with a set of nonacoustic positioning sensors (compasses, tilt-meters, pressure gauges). The acoustic data were inverted using two complementary approaches: a genetic algorithm (GA) like approach and a radial basis functions (RBF) inversion scheme. More traditional methods, based on core sampling, seismic survey and geophone data, together with Hamilton's regression curves, have also been employed on the same tracks, in order to provide a ground truth reference environment. The results of the experiment, can be summarized as follows: 1) the towed array movement is not negligible for the application considered and the use of positioning sensors are essential for a proper acoustic inversion, 2) the inversion with GA and RBF are in good qualitative agreement with the ground truth model, and 3) the GA scheme tends to have better stability properties. On the other hand, repeated in version of successive field measurements requires much less computational effort with RBF     

浅海聚焦定位中简化地声模型的应用研究

为了提高定位算法的环境宽容性,聚焦法将环境参数纳入了寻优空间。聚焦法虽然降低了对环境测量的要求,但是反演参数的增加也增加了反演的复杂性。基于海底反射特性,用两个参数对海底进行建模。通过标准的反演测试问题对简化地声模型在浅海聚焦定位中的有效性进行了分析。结果表明:基于简化地声模型的聚焦定位是可行的。在获得正确定位结果的同时,随着地声参数个数的减少,匹配场处理的便捷性得到了提高。文中引入的简化地声模型是聚焦问题中参数最少的地声模型,它可以有效减少聚焦定位参数维数以提升反演的便捷性。同时,简化地声模型在参数敏感性和耦合性上有较好的表现,这些优点可以保证定位结果的稳健性。     

Combined Geoacoustic Inversion of Propagation and Reverberation Data

Sonar performance predictions in shallow water are strongly dependent on good knowledge of the geoacoustic and scattering properties of the seabed. One technique to extract information about the bottom is to use a towed source and a towed horizontal array. This towed system has been shown to be applicable for characterizing the bottom properties locally by inversion of the acoustic signals received directly on the towed array at short ranges. The same towed system has also been applied to extract bottom properties from long-range reverberation data providing effective bottom properties over a large area. However, independent geoacoustic inversion of the short-range propagation and long-range reverberation data can introduce low sensitivity and uncertainty in the extracted bottom properties. An attempt to resolve this low sensitivity and ambiguity is made by a simultaneous geoacoustic inversion of short-range propagation and long-range reverberation data with the intention of constraining the possible solutions of the bottom properties.      

Effective parameters for matched field geoacoustic inversion in range-dependent environments

Using the well-established technique of geoacoustic inversion, one can estimate a set of acoustic sea-bed parameters from sonar array data. Simultaneously, one can search for geometric parameters such as range, water depth, and hydrophone depth. When the technique is applied in a range-dependent environment, there is a potentially much larger set of parameters to match, unless one has perfect knowledge of the bathymetry. From the point of view of optimization, one needs to handle uncertainties in bathymetry without hugely increasing the amount of computation. A simple time-domain view (which is shown to be equivalent to the adiabatic approximation) suggests that it is sufficient to use a range-independent model with an empirical "effective" depth even when the bottom is not flat. In fact, there is a set of effective environments that will suffice; one can choose whichever is the most convenient. The success of this concept is demonstrated with some test cases from a recent Geoacoustic Inversion Techniques Workshop.     

A hybrid simplex genetic algorithm for estimating geoacousticparameters using matched-field inversion

Matched-fieId inversion (MFI) undertakes to estimate the geometric and geoacoustic parameters in an ocean acoustic scenario by matching acoustic field data recorded at hydrophone array with numerical calculations of the field. The model which provides the best fit to the data is the estimate of the actual experimental scenario. MFI provides a comparatively inexpensive method for estimating ocean bottom parameters over an extensive area. The basic components of the inversion process are a sound propagation model and matching (minimization) algorithm. Since a typical MFI problem requires a large number of computationally intensive sound propagation calculations, both of these components have to be efficient. In this study, a hybrid inversion algorithm which uses a parabolic equation propagation model and combines the downhill simplex algorithm with genetic algorithms is introduced. The algorithm is demonstrated on synthetic range-dependent shallow-water data generated using the parabolic equation propagation model. The performance for estimating the model parameters is compared for realistic signal-to-noise ratios in the synthetic data     

Geoacoustic model parameter estimation using a bottom-moored hydrophone array

This paper describes results from an experiment carried out to investigate geoacoustic inversion with a bottom-moored hydrophone array located in the shallow waters of the Timor Sea off the northern coast of Australia. The array consisted of two arms in a V shape, horizontally moored at a site that was essentially flat over a large area. Hydrophone positions were estimated using an array element localization (AEL) technique that established relative uncertainties of less than 1 m on the seafloor. The data used for geoacoustic inversion were from experiments with continuous wave (CW) tones in the 80- to 195-Hz band transmitted from a towed projector. A hybrid search algorithm determined the set of geoacoustic model parameters that maximized the Bartlett fit (averaged coherently spatially at each tone and incoherently over frequency) between the measured and modeled data at the array. Due to the long range experimental geometry, the inversion was sensitive to attenuation in the sediment. The inverted geoacoustic profile performed well in a simple test for localizing the sound source at other sites in the vicinity of the array. Range-depth localization performance for the horizontal array was comparable to that for an equivalent vertical array.     

消频散变换反演海底地声参数

提出1种将消频散变换应用到海底地声参数反演的方法。对单一水听器接收声压信号进行消频散处理后,根据群延时差建立代价函数,反演得到主要海底参数,最后根据贝叶斯统计理论给出了待反演地声参数的边缘后验概率密度。对单层波导进行仿真证明这种新方法的有效性。     

Geoacoustic inversion of tow-ship noise via near-field-matched-field processing

This paper discusses geoacoustic inversion from tow-ship noise data acquired via a horizontal towed array. Through simulations and experimental results, it is shown that even very quiet ships radiate sufficient noise power to enable self-noise inversion of basic geoacoustic parameters such as effective bottom velocity. The experimental results presented are particularly encouraging in view of the high level of interference shown to be tolerated from nearby shipping.     

Estimates of geoacoustic model parameters from inversions of horizontal and vertical line array data

This paper describes results from geoacoustic inversion of low-frequency acoustic data recorded at a receiving array divided into two sections, a sparse bottom laid horizontal array (HLA) and a vertical array (VLA) deployed in shallow water. The data are from an experiment conducted by the Norwegian Defence Research Establishment (FFI) in the Barents Sea, using broadband explosives (shot) sources. A two-layer range-independent geoacoustic model, consistent with seismic profiles from the area, described the environment. Inversion for geoacoustic model parameters was carried out using a fast implementation of the hybrid adaptive simplex simulated annealing (ASSA) inversion algorithm, with replica fields computed by the ORCA normal mode code. Low-frequency (40-128 Hz) data from six shot sources at ranges 3-9 km from the array were considered. Estimates of sediment and substrate p-wave velocities and sediment thickness were found to be consistent between independent inversions of data from the two sections of the array.     

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.     

Reverberation vertical coherence and sea-bottom geoacoustic inversion in shallow water

Optimal array-processing techniques in the ocean often require knowledge of the spatial coherence of the reverberation. A mathematical model is derived for the reverberation vertical coherence (RVC) in shallow water (SW). A method for analysis of RVC data is introduced. Measured reverberation cross-correlation coefficients as a function of time and frequency, obtained during the Asian Seas International Acoustic Experiment (ASIAEX) in the East China Sea, are reported. SW reverberation from a single shot provides a continuous spatial sampling of the surrounding sound field up to several tens of kilometers and holds valuable information on the geoacoustic properties of the sea floor over this distance. SW reverberation data can, therefore, be used as the basis for a quick and inexpensive method for geoacoustic inversion and has the obvious advantage that acquiring the data in situ requires only a single platform. This paper considers the use of the vertical coherence of the reverberation as the starting point for such an inversion. Sound speed and attenuation in the sea bottom at the ASIAEX site are obtained over a frequency range of 100-1500 Hz by finding values that provide the best match between the measured and predicted RVC.     

Geoacoustic characterization of a range-dependent ocean environment using towed array data

This paper examines geoacoustic inversion over a range-dependent multiple-layer seabed using a towed acoustic source and towed horizontal array. The approach is based on combining the results of a series of short-range, range-independent inversions to form a range-dependent representation of the environment. The data were collected in the Strait of Sicily during the MAPEX 2000 experiment. Issues such as the resolvability of multilayer structure and the sensitivity of various geoacoustic parameters are investigated by inversion of simulated data and by comparison of the MAPEX 2000 inversion results to a high-resolution seismic profile and to sediment core measurements. It appears that two, and in some cases possibly three, sediment layers can be resolved.     

Geoacoustic Information Content of Horizontal Line Array Data

This paper examines the effectiveness of horizontal line arrays (HLAs) for matched-field inversion (MFI) by quantifying geoacoustic information content for a variety of experiment and array factors, including array length and number of sensors, source range and bearing, source-frequency content, and signal-to-noise ratio (SNR). Emphasis is on bottom-moored arrays, while towed arrays are also considered, and a comparison with vertical line array (VLA) performance is made. The geoacoustic information content is quantified in terms of marginal posterior probability distributions (PPDs) for model parameters estimated using a fast Gibbs sampler approach to Bayesian inversion. This produces an absolute, quantitative estimate of the geoacoustic parameter uncertainties which can be directly compared for various experiment and array factors.     

An Equivalent Transform Method for Evaluating the Effect of Water-Column Mismatch on Geoacoustic Inversion

An "equivalent transform method" for quantitatively evaluating the effect of water-column mismatch on geoacoustic inversions is presented. This method is based on the concept of error transferral from one medium to another and is derived from linear perturbative inverse theory. To illustrate the method, geoacoustic inversions using adiabatic mode data, including wave numbers, group velocities, and travel times, are considered. In the test cases, both linear and nonlinear internal waves are considered as the causes of the water-column mismatch, and the inversion errors due to the water-column mismatch in each case are discussed. In the case of linear internal waves, range-averaged inversion errors are largely eliminated at the full cycle distance of the internal wave; however, in the case of nonlinear internal waves, the range-averaged inversion errors are accumulated but scaled down with increasing range. Furthermore, the inversion errors produced by water-column mismatch will be large if the modes used for inversion are sensitive to the mismatch; for instance, using higher order modes might lead to increased error even though it would normally improve inversion results. Using lower frequency signal, which penetrates deeper into the bottom, extends the effective depth of the inversion solutions, but this also allows the water-column mismatch to transfer error into the deeper bottom     

Geoacoustic inversion of range-dependent data with added Gaussian noise

A rotated coordinates inversion algorithm is used on subsets of the Inversion Techniques 2001 Geoacoustic Workshop data, to which white Gaussian noise is added. The resulting data sets are equivalent to noisy broad-band signals received on a horizontal line array (HLA) during a single integration time interval. The inversions are performed using a technique called systematic decoupling using rotated coordinates (SDRC), which expands the original idea of rotated coordinates by using multiple sets of rotated coordinates, each corresponding to a different set of bounds, to systematically decouple the unknowns in a series of efficient simulated annealing inversions. The cost function minimized in the inversion is based on the coherent broad-band correlation between data and model cross spectra, which increases the coherence gain of the signal relative to incoherent noise. Using the coherent broad-band cost function with sparse HLA-like data sets, the SDRC inversion method yields good estimates for the sensitive environmental parameters for signal-to-noise ratios as low as -15 dB.     

Application of an improved self-starter to geoacoustic inversion

The self-starter is improved using the operator of the split-step Pade solution. In addition to providing greater stability and being applicable closer to the source, the improved self-starter is an efficient forward model for geoacoustic inversion. It is necessary to solve only O(10) tridiagonal systems of equations to obtain the acoustic field on a vertical array located O(10) wavelengths from a source. This experimental configuration is effective for geoacoustic inverse problems involving unknown parameters deep in the ocean bottom. For problems involving depth-dependent acoustic parameters, the improved self-starter can be used to solve nonlinear inverse problems involving O(10) unknown sediment parameters in less than a minute on the current generation of workstations     

Estimating geoacoustic parameters using matched-field inversion methods

In this paper, we use matched-field inversion methods to estimate the geoacoustic parameters for three synthetic test cases from the Geoacoustic Inversion Techniques Workshop held in May 2001 in Gulfport, MS. The objective of this work is to use a sparse acoustic data set to obtain estimates of the parameters as well as an indication of their uncertainties. The unknown parameters include the geoacoustic properties of the sea bed (i.e., number of layers, layer thickness, density, compressional speed, and attenuation) and the bathymetry for simplified range-dependent acoustic environments. The acoustic data used to solve the problems are restricted to five frequencies for a single vertical line array of receivers located at one range from the source. Matched-field inversion using simplex simulated annealing optimization is initially used to find a maximum-likelihood (ML) estimate. However, the ML estimate provides no information on the uncertainties or covariance associated with the model parameters. To estimate uncertainties, a Bayesian formulation of matched-field inversion is used to generate posterior probability density distributions for the parameters. The mean, covariance, and marginal distributions are determined using a Gibbs importance sampler based on the cascaded Metropolis algorithm. In most cases, excellent results were obtained for relatively sensitive parameters such as wave speed, layer thickness, and water depth. The variance of the estimates increase for relatively insensitive parameters such as density and wave attenuation, especially when noise is added to the data.