首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
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.  相似文献   

2.
This paper presents an adaptive hybrid algorithm to invert ocean acoustic field measurements for seabed geoacoustic parameters. The inversion combines a global search (simulated annealing) and a local method (downhill simplex), employing an adaptive approach to control the trade off between random variation and gradient-based information in the inversion. The result is an efficient and effective algorithm that successfully navigates challenging parameter spaces including large numbers of local minima, strongly correlated parameters, and a wide range of parameter sensitivities. The algorithm is applied to a set of benchmark test cases, which includes inversion of simulated measurements with and without noise, and cases where the model parameterization is known and where the parameterization most be determined as part of the inversion. For accurate data, the adaptive inversion often produces a model with a Bartlett mismatch lower than the numerical error of the propagation model used to compute the replica fields. For noisy synthetic data, the inversion produces a model with a mismatch that is lower than that for the true parameters. Comparison with previous inversions indicates that the adaptive hybrid method provides the best results to date for the benchmark cases  相似文献   

3.
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  相似文献   

4.
A method for estimating properties of the ocean bottom such as bathymetry and geoacoustic parameters such as sound speed, density and attenuation, using matched-field inversion is considered. The inversion can be formulated as an optimization problem by assuming a discrete model of unknown parameters and a bounded search space for each parameter. The optimization then involves finding the set of parameter values which minimizes the mismatch between the measured acoustic field and modeled replica fields. Since the number of possible models can be extremely large, the method of simulated annealing, which provides an efficient optimization that avoids becoming trapped in suboptimal solutions, has been used. The matching fields are computed using a normal mode model. In inversions for range-dependent parameters, the adiabatic approximation is employed. This allows mode values to be precomputed for a grid of parameter values and stored in look-up tables for fast reference, which greatly improves computational efficiency. Synthetic inversion examples are presented for realistic range-independent and range-dependent environments  相似文献   

5.
A perturbative inversion method for estimating sediment compressional-wave-speed profiles from modal travel-time data is extended to include range-dependent environments. The procedure entails dividing a region into range-independent sections and obtaining estimates of the sediment properties for each region. Inversion results obtained using synthetic data show that range-dependent properties can be obtained if an experiment is designed to include multiple source/receiver combinations. This approach is applied to field data collected during the 2006 Shallow Water Experiment (SW06). The sediment compressional-wave-speed profiles resulting from analysis of the field data are evaluated by comparing acoustic fields predicted based on the inversion to acoustic fields measured during a different experiment conducted in the same region. The model is also compared to seismic reflection survey data collected during SW06. Resolution and variance estimated for the inversion results are also presented.   相似文献   

6.
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  相似文献   

7.
The problem of rapid environmental assessment in a range-dependent environment is addressed. For rapid assessment, the exact geoacoustic parameters are not required, nor is it a requirement that the exact structure of the acoustic field (location of peaks and s) be matched by an acoustic prediction model. The parameters that are relevant are the overall transmission loss (incoherent TL), the time spread (/spl tau/), and the slopes of the range/frequency interference patterns (/spl beta/, the waveguide invariant). The rapid geoacoustic characterization algorithm uses a homogeneous single-sediment layer overlying a hard acoustic basement model to optimally match the predicted acoustic observables with those estimated from data. The approach is presented here and is applied to the range-dependent benchmark cases TC1 and TC2 from the Inversion Techniques Workshop held in Gulfport, MS, in May 2001. The technique successfully reproduces the acoustic observables and estimates the sediment sound-speed, density, and attenuation profiles, as well as the sediment thickness.  相似文献   

8.
An efficient method for geoacoustic inversions in a range-dependent ocean waveguide is implemented and tested with synthetic data. This method combines a simulated annealing search with an optimal coordinate rotation that increases the efficiency of navigating parameter landscapes for which parameter coupling is important. The coordinate rotation associated with the parameter couplings also provides information about which parameters are resolvable for a particular inversion frequency and array geometry. Using this information, results from several single-frequency inversions can be combined to obtain an estimate for the sediment parameters.  相似文献   

9.
Phenomenological and global optimization inversion   总被引:1,自引:0,他引:1  
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.  相似文献   

10.
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.  相似文献   

11.
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.  相似文献   

12.
Inversion methods have been developed over the past decade to extract information about unknown ocean-bottom environments from acoustic field data. This paper summarizes results from the Office of Naval Research/Space and Naval Warfare Systems Command (SPAWAR) Geoacoustic Inversion Techniques Workshop, which was designed to benchmark present-day inversion methods. The format of the workshop was a blind test to estimate unknown geoacoustic profiles by inversion of synthetic acoustic field data. The fields were calculated using a high-angle parabolic approximation and verified using coupled normal modes for three range-dependent shallow-water test cases: a monotonic slope; a shelf break; and a fault intrusion in the sediment. Geoacoustic profiles were generated to simulate sand, silt, and mud sediments in these environments. Several different approaches for inverting the acoustic field data were presented at the workshop: model-based matched-field methods; perturbation methods; methods using transmission loss data; and methods using horizontal array information. An effective inversion must provide both an estimate of the bottom parameters and a measure of the uncertainty of the estimated values. New methods were presented at the workshop to formalize the measure of uncertainty in the inversion. Comparisons between the different inversions are discussed in terms of a metric-based transmission loss calculated using the inverted profiles. The results demonstrate the effectiveness of present-day inversion techniques and indicate the limits of their capabilities for range-dependent waveguides.  相似文献   

13.
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.  相似文献   

14.
It is extremely difficult to determine shallow ocean bottom properties (such as sediment layer thicknesses, densities, and sound speeds). However, when acoustic propagation is affected by such environmental parameters, it becomes possible to use acoustic energy as a probe to estimate them. Matched-field processing (MFP) which relies on both field amplitude and phase can be used as a basis for the inversion of experimental data to estimate bottom properties. Recent inversion efforts applied to a data set collected in October 1993 in the Mediterranean Sea north of Elba produce major improvements in MFP power, i.e., in matching the measured field by means of a model using environmental parameters as inputs, even using the high-resolution minimum variance (MV) processor that is notoriously sensitive and usually results in very low values. The inversion method applied to this data set estimates water depth, sediment thickness, density, and a linear sound-speed profile for the first layer, density and a linear sound-speed profile for a second layer, constant sound speed for the underlying half space, array depth, and source range and depth. When the inversion technique allows for the array deformations in range as additional parameters (to be estimated within fractions of a wavelength, e.g., 0.1 m), the MFP MV peak value for the Med data at 100 Hz can increase from 0.48 (using improved estimates of environmental parameters and assuming a vertical line array) to 0.68 (using improved estimates of environmental parameters PLUS improved phone coordinates). The ideal maximum value would be 1.00 (which is achieved for the less sensitive Linear processor). However, many questions remain concerning the reliability of these inversion results and of inversion methods in general  相似文献   

15.
A common problem in sonar system prediction is that the ocean environment is not well known. Utilizing probabilistic based results from geoacoustic inversions we characterize parameters relevant to sonar performance. This paper describes the estimation of transmission loss and its statistical properties based on posterior parameter probabilities obtained from inversion of ocean acoustic array data. This problem is solved by first finding an ensemble of relevant environmental model parameters and the associated posterior probability using a likelihood based inversion of the acoustic array data. In a second step, each realization of these model parameters is weighted with their posterior probability to map into the transmission loss domain. This approach is illustrated using vertical-array data from a recent benchmark data set and from data acquired during the Asian Seas International Acoustics Experiment (ASIAEX) 2001 in the East China Sea. The environmental parameters are first estimated using a probabilistic-based geoacoustic inversion technique. Based on the posterior probability that each of these environmental models fits the ocean acoustic array data, each model is mapped into transmission loss. This enables us to compute a full probability distribution for the transmission loss at selected frequencies, ranges, and depths, which potentially could be used for sonar performance prediction.  相似文献   

16.
This paper presents the results of a perturbative inverse approach applied to the range-dependent acoustic data provided as part of the Geoacoustic Inversion Techniques Workshop. The method is based on the Hankel transform relationship between complex-pressure field data measured on a horizontal array and the depth-dependent Green's function for a horizontally stratified medium. The input data to the inversion algorithm are discrete values of horizontal wavenumbers estimated for locally range-independent subapertures of the acoustic data. Inversion results are presented with emphasis placed on inverting for compressional wave speed as a function of both range and depth in the bottom.  相似文献   

17.
This paper proposes a method, based on the Biot model, for estimating the physical and acoustic properties of surficial ocean sediments from normal incidence reflection data acquired by a chirp sonar. The inversion method estimates sediment porosity from reflection coefficient measurements and, using the estimated porosity and the measured change in fast wave attenuation with frequency, estimates the permeability of the top sediment layer. The spectral ratio of echoes from the interface at the base of the upper sediment layer and from the sediment-water interface provides a measure of the change in attenuation with frequency. Given the porosity and permeability estimates, the Kozeny-Carman equation provides the mean grain size and the inversion method yields the acoustic properties of top sediment layer. The inversion technique is tested using chirp sonar data collected at the 1999 Sediment Acoustics Experiment (SAX-99) site. Remote estimates of porosity, grain size, and permeability agree with direct measurements of those properties.  相似文献   

18.
19.
不确定海洋环境中基于贝叶斯理论的多声源定位算法   总被引:2,自引:0,他引:2  
环境参数失配导致定位性能大幅度下降是匹配场定位所面临的难题之一。应用贝叶斯理论对环境聚焦,是当前解决该难题的研究热点。环境聚焦方法的实质是将未知环境参数和声源位置联合优化估计,当出现多个目标时,估计的参数会随着声源个数成倍增加,因此不得不利用有限的观测信息来实现众多参数的估计。本文采用最大似然比方法,获得信号源谱和误差项的最大似然估计,实现这些敏感性较弱参数的间接反演,有效降低了反演参数维数和定位算法复杂度。针对遗传算法的早熟和稳定性差的问题,改进了似然函数的经验表达式。将多维后验概率密度在参数起伏变化范围内积分,得到反演参数的一维边缘概率分布,求解最优值的同时进行反演结果的不确定性分析。本文仿真了位于相同距离、不同深度的两个声源,使用仿真实验验证了提出算法的有效性。  相似文献   

20.
In October 1997, the EnVerse 97 shallow-water acoustic experiments were jointly conducted by SACLANT Centre, TNO-FEL, and DERA off the coast of Sicily, Italy. The primary goal of the experiments was to determine the sea-bed properties through inversion of acoustic data. Using a towed source, the inversion method is tested at different source/receiver separations in an area with a range-dependent bottom. The sources transmitted over a broadband of frequencies (90-600 Hz) and the signals were measured on a vertical array of hydrophones. The acoustic data were continuously collected as the range between the source and receiving array varied from 0.5 to 6 km. An extensive seismic survey was conducted along the track providing supporting information about the layered structure of the bottom as well as layer compressional sound speeds. The oceanic conditions were assessed using current meters, satellite remote sensing, wave height measurements, and casts for determining conductivity and temperature as a function of water depth. Geoacoustic inversion results taken at different source/receiver ranges show sea-bed properties consistent with the range-dependent features observed in the seismic survey data. These results indicate that shallow-water bottom properties may be estimated over large areas using a towed source fixed receiver configuration  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号