共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper applies a Bayesian formulation to range-dependent geoacoustic inverse problems. Two inversion methods, a hybrid optimization algorithm and a Bayesian sampling algorithm, are applied to some of the 2001 Inversion Techniques Workshop benchmark data. The hybrid inversion combines the local (gradient-based) method of downhill simplex with the global search method of simulated annealing in an adaptive algorithm. The Bayesian inversion algorithm uses a Gibbs sampler to estimate properties of the posterior probability density, such as mean and maximum a posteriori parameter estimates, marginal probability distributions, highest-probability density intervals, and the model covariance matrix. The methods are applied to noise-free and noisy benchmark data from shallow ocean environments with range-dependent geophysical and geometric properties. An under-parameterized approach is applied to determine the optimal model parameterization consistent with the resolving power of the acoustic data. The Bayesian inversion method provides a complete solution including quantitative uncertainty estimates and correlations, while the hybrid inversion method provides parameter estimates in a fraction of the computation time. 相似文献
2.
Fialkowski L.T. Lingevitch J.F. Perkins J.S. Dacol D.K. Collins M.D. 《Oceanic Engineering, IEEE Journal of》2003,28(3):370-379
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. 相似文献
3.
4.
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. 相似文献
5.
Fallat M.R. Nielsen P.L. Dosso S.E. Siderius M. 《Oceanic Engineering, IEEE Journal of》2005,30(1):198-206
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. 相似文献
6.
Chapman N.R. Chin-Bing S. King D. Evans R.B. 《Oceanic Engineering, IEEE Journal of》2003,28(3):320-330
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. 相似文献
7.
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. 相似文献
8.
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 相似文献
9.
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 相似文献
10.
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. 相似文献
11.
Phenomenological and global optimization inversion 总被引：1，自引：0，他引：1
Gerstoft P. Hodgkiss W.S. Kuperman W.A. Heechun Song 《Oceanic Engineering, IEEE Journal of》2003,28(3):342-354
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. 相似文献
12.
采用结合匹配法和经验正交函数法的射线声层析反演方法,针对流场水平分布不均匀特点,使用距离等效分段方法对反演算法进行改进,利用三个断面声层析数据,对胶州湾口潮流场的垂直和水平分布进行反演。与传统方法比较,大大降低与实测流速垂直分布间的偏差,平均偏差小于0.02 m/s,流场垂直分布间的相关系数提高到0.85以上。分析了声传播断面上不同的流场结构,以及涨潮期间在胶州湾团岛附近出现的涡流现象,并计算得到湾口海水流量变化。结果表明,改进的射线声层析反演方法可有效地用于水平分布不均匀流场监测,仅采用少量声学观测站位,即可获得大范围的复杂流场三维信息,有利于近海海洋资源开发、海洋环境保护和船舶航行安全。 相似文献
13.
《Oceanic Engineering, IEEE Journal of》2009,34(4):516-525
14.
传统的四维变分数据同化方法在同化观测资料的同时可以对数值模式参数进行优化,然而传统的四维变分方法需要针对不同的数值模式编写特有的伴随模式,因此算法的可移植性差,同时计算时耗费大量资源。本文提出了一种新的基于解析四维集合变分的参数优化方法,该方法以迭代搜索得到的模式参数为基准展开扰动并构建样本集合,由此显式地计算协方差矩阵,并得到代价函数极小值的解析解,从而避免了伴随模式的使用。基于Lorenz-63模型对该方法进行单参数和多参数数值试验和优化效果检验,并在不同的同化时间窗口长度和观测采样间隔情况下,采用传统四维变分方法与之进行对比,结果显示,新方法表现出与传统四维变分相同的优化性能,都能有效收敛到真值,而新方法不需要计算伴随模式,可移植性好。本文还测试了不同的集合成员个数和模式参数真值的情况下新方法的同化效果,结果表明,新方法对集合样本个数及模型参数真值不敏感,采用较少的集合样本即可完成数据同化。 相似文献
15.
经验模态分解(empirical mode decomposition,简称EMD)是反演海洋内波参数的有效方法之一,但由于EMD存在模态混叠等问题,对海洋内波进行参数反演时会产生一定误差。相比较于EMD,变分模态分解(variational mode decomposition,简称VMD)能够有效地抑制模态混叠现象。为了更好地对海洋内波进行参数反演,提出了一种基于VMD对合成孔径雷达(synthetic aperture radar,简称SAR)遥感图像中的内波参数进行自动反演的方法。该方法先对SAR图像进行Canny处理,获取图像中的内波条纹信息,再根据内波传播方向自动选取灰度剖面;然后利用集合经验模态分解(ensemble empirical mode decomposition,简称EEMD)信号特征自适应分解模态函数的特点,再将分解得到的有效模态数作为VMD中参数K的参考值;最后利用VMD分解后的数据进行内波参数反演。试验结果表明:通过对Canny预处理后的条纹信息进行灰度剖面自动选取,解决了人为选取剖面所可能导致的误差;通过对剖面信号进行VMD处理不仅解决了EMD模态混叠的问题,成功地反演出内波的前导波振幅,而且所反演的结果与EEMD反演参数以及实测资料数据吻合得很好。 相似文献
16.
针对机载合成孔径雷达(SAR)对海探测特点,采用多入射角法从SAR数据本身得到与海浪参数反演区域时空匹配的同步海面风速和风向,并结合线性变换关系,计算得到海浪初猜谱对应的仿真SAR图像谱,将仿真SAR图像谱和观测SAR图像谱输入代价函数中进行迭代运算,通过非线性方程的解算得到最适海浪谱;采用交叉谱法去除海浪传播180°方向模糊,最终得到海浪参数。论文提出的基于同步风场的机载SAR海浪参数反演方法,充分利用了机载SAR海洋环境探测的优势,解决了传统SAR海浪参数反演中初猜谱构造依赖外部风场的问题,机载同步飞行试验的海浪参数反演结果与浮标观测值的有效波高、波向的均方根误差分别为0.23 m和13.23°,验证了该方法的有效性,可为机载SAR海浪参数反演业务化提供支持。 相似文献
17.
Parameter estimation is defined as the process to adjust or optimize the model parameter using observations. A long-term problem in ensemble-based parameter estimation methods is that the parameters are assumed to be constant during model integration. This assumption will cause underestimation of parameter ensemble spread,such that the parameter ensemble tends to collapse before an optimal solution is found. In this work, a two-stage inflation method is developed for parameter estimation, which ... 相似文献
18.
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 相似文献
19.
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. 相似文献
20.
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 相似文献