Range-dependent geoacoustic inversion: results from the Inversion Techniques Workshop

Model-based geoacoustic inversion in range-dependent underwater environments is a challenging task constrained by data quality (synthetic or measured) and propagation-model efficiency and accuracy. The Inversion Techniques Workshop (ITW), held in Gulfport, MS, May 15-18, 2001, was organized for the acoustics community to present state-of-the-art numerical geoacoustic inversion capabilities in range-dependent shallow-water environments. The organizers defined five range-dependent test cases (three synthetic and two experimental cases). Two of the synthetic cases were adopted for geoacoustic inversion in this paper. The first test case (TC1) is a monotonic down-slope bathymetry problem and the adiabatic normal-mode model PROSIM was applied for the inversion in this case. The second test case (TC3) is a flat-bottom case with an intrusion. The forward model used in this case was RAMGEO. The global optimization package SAGA was used for geoacoustic inversion of the synthetically generated reference solutions for TC1 and TC3. In general, the geoacoustic inversion results are in good agreement with the true solutions provided by the organizers. The results obtained demonstrate the feasibility of performing geoacoustic inversion in synthetic range-dependent shallow-water environments. However, results show that the propagation model choice in the inversion is strongly dependent on the specific range-dependent environment.     

Matched field inversion for geoacoustic model parameters usingadaptive simulated annealing

A method is described for the estimation of geoacoustic model parameters by the inversion of acoustic field data using a nonlinear optimization procedure based on simulated annealing. The cost function used by the algorithm is the Bartlett matched-field processor (MFP), which related the measured acoustic field with replica fields calculated by the SAFARI fast field program. Model parameters are perturbed randomly, and the algorithm searches the multidimensional parameter space of geoacoustic models to determine the parameter set that optimizes the output of the MFP. Convergence is driven by adaptively guiding the search to regions of the parameter space associated with above-average values of the MFP. The performance of the algorithm is demonstrated for a vertical line array in a shallow water enviornment where the bottom consists of homogeneous elastic solid layers. Simulated data are used to determine the limits on estimation performance due to error in experimental geometry and to noise contamination. The results indicate that reasonable estimates are obtained for moderate conditions of noise and uncertainty in experimental geometry     

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.      

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.     

Using the adaptive simulated annealing algorithm to estimate ocean-bottom geoacoustic properties from measured and synthetic transmission loss data

This paper presents the results obtained using the adaptive simulated annealing (ASA) algorithm to invert the test cases from the Geoacoustic Inversion Techniques Workshop held in May 2001. The ASA algorithm was chosen for use in our inversion software for its speed and robustness when searching the geoacoustic parameter solution space to minimize the difference between the observed and the modeled transmission loss (TL). Earlier work has shown that the ASA algorithm is approximately 15 times faster than a modified Boltzmann annealing algorithm, used in prior versions of our TL inversion software, with comparable fits to the measured data. Results are shown for the synthetic test cases, 0 through 3, and for the measured data cases, 4 and 5. The inversion results from the synthetic test cases showed that subtle differences between range-dependent acoustic model version 1.5, used to generate the test cases, and parabolic equation (PE) 5.0, used as the propagation loss model for the inversion, were significant enough to result in the inversion algorithm finding a geoacoustic environment that produced a better match to the synthetic data than the true environment. The measured data cases resulted in better fits using ASTRAL automated signal excess prediction system TL 5.0 than using the more sophisticated PE 5.0 as a result of the inherent range averaging present in the ASTRAL 5.0 predictions.     

Geoacoustic Uncertainties From Viscoelastic Inversion of Seabed Reflection Data

This paper applies nonlinear Bayesian inversion to seabed reflection data to estimate viscoelastic parameters of the upper sediments. The inversion provides maximum a posteriori probability (MAP) parameter estimates with uncertainties quantified in terms of marginal probability distributions, variances, and credibility intervals; interparameter relationships are quantified by correlations and joint marginal distributions. The inversion is applied to high-resolution reflectivity data from two sites in the Strait of Sicily. One site is characterized by low-speed sediments, resulting in data with a well-defined angle of intromission; the second is characterized by high-speed sediments, resulting in a critical angle. Data uncertainties are quantified using several approaches, including maximum-likelihood (ML) estimation, treating uncertainties as nuisance parameters in the inversion, and analysis of experimental errors. Statistical tests are applied to the data residuals to validate the assumed uncertainty distributions. Excellent results (i.e., small uncertainties) are obtained for sediment compressional-wave speed, compressional attenuation, and density; shear parameters are less well determined although low shear-wave speeds are indicated. The Bayesian analysis provides a quantitative comparison of inversion results for the two sites in terms of the resolution of specific geoacoustic parameters, and indicates that the geoacoustic information content is significantly higher for intromission data     

Geoacoustic inversion of broad-band acoustic data in shallow water on a single hydrophone

The inversion of broad-band low-frequency acoustic signals received on sparse arrays can lead to robust and efficient estimations of sea-bed properties. This paper describes a shallow-water geoacoustic inversion scheme based on the use of a model-based matched-impulse response on a single hydrophone. Results from the INTIMATE'96 experiment on the Portuguese shelf break are reviewed. In order to minimize the effects of strong time variability due to internal tides, only the time-stable waterborne bottom-surface reflected arrivals are exploited. A quasi-linear inversion algorithm is first applied to refine the geometry of the experiment. Then, inversion of bottom parameters is performed with an objective function that only makes use of the bottom-surface reflected arrivals' amplitudes. The experimental results show that broad-band transmissions (300-800 Hz) received on a single hydrophone, combined with the use of a simple eigenray code, are sufficient to correctly resolve geometrical parameters and bottom features. The analysis of the reflection coefficients both on simulated and real data helps to understand the validity of the inverted parameters and to derive the basis of an equivalent medium concept for geoacoustic inversion based on a "through-the-sensor" approach.     

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

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

Performance comparison between vertical and horizontal arrays for geoacoustic inversion

Most of the research on model-based geoacoustic inversion techniques has concentrated on data collected using moored vertical receiver arrays. However, there are many advantages to considering geoacoustic inversion using a towed horizontal array. Towed arrays are easily deployed from a moving platform; this mobility makes them well suited for surveying large areas for sea-bed properties. Further, if a model-based geoacoustic inversion scheme uses both a towed source and array, the separation between the two can be kept short, which reduces the requirement for range-dependent modeling. Range-independent modeling is used for inverting all the horizontal array data considered in this paper. Using the Inversion Techniques Workshop Benchmark Test Cases, the performance of a horizontal (simulated towed) and vertical arrays are compared and found to be very similar. However, it will be shown that, for Benchmark Test Case 3, where the bathymetry is flat and a hidden bottom intrusion exists, a towed horizontal array is ideal for determining the range-dependent sea-bed properties. The practical advantages of using a towed array are clear and the purpose of this paper is to show that the performance is similar (and in some cases better) than using moored vertical arrays.     

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.     

An investigation of algorithm-induced variability in geoacoustic inversion

This paper investigates the inherent variability in the results of matched-field geoacoustic inversion algorithms. This algorithm-induced variability must be considered when interpreting inversion results in terms of environmental changes as a function of time or space. Fast simulated annealing (FSA), genetic algorithms (GA), and a hybrid algorithm (adaptive simplex simulated annealing; ASSA) are compared by performing multiple inversions of benchmark synthetic data (noise free and noisy) and acoustic data measured over both low- and high-speed sea-bed sediments in the MAPEX 2000 experiment. ASSA produced the lowest variability in inversion results for all cases, followed by GA and FSA. For the high-speed MAPEX 2000 case, the variability is essentially negligible, while for the low-speed case the variability is significant as compared with environmental variations reported in the literature.     

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.     

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.     

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     

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     

Phenomenological and global optimization inversion

This paper discusses geoacoustic inversion results based on benchmark range-dependent data using SAGA, a global inversion package, and using phenomenological inversions. In phenomenological inversions, physical and signal-processing approaches are used to enhance the data to extract specific features. The global optimization approach is carried out on complex-valued vertical array data, transmission loss data, and reverberation data. The importance of checking the solution is emphasized by inspecting the match with the data and the error estimates and by checking the solution using data that has not been used in constructing the solution. The results show that we are able to estimate the geoacoustic parameters and that these parameters could be used to predict the field for different frequencies and/or source-receiver geometry than used in the inversion.     

Active rapid geoacoustic characterization using a seismic survey source

A survey of received acoustic energy levels from a seismic profiler were performed in Long Beach Harbor, CA, for compliance with the Marine Mammal Protection Act (MMPA). In addition to direct acoustic measurements, a rapid geoacoustic inversion algorithm was applied to the data to estimate the sediment properties acoustically. This inversion algorithm has matching criteria based on time spread, range-frequency interference patterns, and the range dependence of transmission loss. Self-consistency was checked by comparing acoustic measurements with predictions based on the inversion. With an estimated geoacoustic profile, predictions of received levels as a function of position in the range-dependent environment of Long Beach Harbor were then performed.     

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     

A New Multistep Algorithm for Geoacoustic Inversion via Subspace Extraction

This communication presents a new multistep matched-field algorithm for geoacoustic inversion by subspace extraction with a threshold. In this algorithm, according to the varying sensitivities of geoacoustic parameters, parameters are separated into several subsets (or subspaces). Then, inversions are carried out in each sensitive subspace using an optimization algorithm, and for each inversion, a sub-subspace is extracted where values of objective functions are lower than a given threshold. Finally, in all the extracted sub-subspaces combined with the subspace of insensitive parameters, an inversion is performed for all parameters to find the optimal solution. After the extracting process, the search space is greatly reduced, and generally, the true parameter values will not be excluded from the sub-subspace if a reasonable threshold is designed. Thus, higher efficiency and accuracy can be obtained when compared with other algorithms. Simulation is carried out on synthetic data and results indicate that the new algorithm's performance is significantly superior to those of other algorithms.      

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