A robust spatial filtering technique for multisource localization and geoacoustic inversion

Geoacoustic inversion and source localization using beamformed data from a ship of opportunity has been demonstrated with a bottom-mounted array. An alternative approach, which lies within a class referred to as spatial filtering, transforms element level data into beam data, applies a bearing filter, and transforms back to element level data prior to performing inversions. Automation of this filtering approach is facilitated for broadband applications by restricting the inverse transform to the degrees of freedom of the array, i.e., the effective number of elements, for frequencies near or below the design frequency. A procedure is described for nonuniformly spaced elements that guarantees filter stability well above the design frequency. Monitoring energy conservation with respect to filter output confirms filter stability. Filter performance with both uniformly spaced and nonuniformly spaced array elements is discussed. Vertical (range and depth) and horizontal (range and bearing) ambiguity surfaces are constructed to examine filter performance. Examples that demonstrate this filtering technique with both synthetic data and real data are presented along with comparisons to inversion results using beamformed data. Examinations of cost functions calculated within a simulated annealing algorithm reveal the efficacy of the approach.     

一种地声参数的联合反演方法

根据地声参数对不同声场物理量影响不同,提出了一种利用简正波频散特征结合声传播损失反演地声参数的联合反演方法。首先,考虑到简正波的频散特性(群速度)对海底的密度和声速较为敏感,而对海底吸收系数不敏感,利用自适应时频分析方法,获得不同频率不同号数简正波的到达时间差,以此作为代价函数,采用全局优化算法,反演出海底密度和海底声速的分层结构,并用概率统计的方法评价反演结果的有效性。反演出海底密度和声速后,利用实验船辐射噪声得到随距离连续的声传播损失来反演出海底吸收系数。最后,把反演的参数很好的用于声源匹配定位验证了反演结果的有效性。     

A hybrid scheme for geoacoustic inversion

A hybrid inversion scheme of geoacoustic parameters, which combines dispersion characteristics of normal modes with transmission loss, has been proposed based on the facts that different parameters have different effects on sound field. First, considering the dispersion characteristics of normal modes were sensitive to the density and the sound speed of seabottom, and not sensitive to the seabottom attenuation coefficient, the group time delays of different modes and different frequencies were acquired through adaptive optimal kernel time-frequency representation, then the group time delays were taken as the cost function, and using a global optimal algorithm, the seabottom density and the profile of seabottom sound speed were deduced. Successively, the validity of inversion results was evaluated by the a posteriori probability. Last, the attenuation coefficient of sediment is inverted using the transmission loss data recorded from the ship radiation noise. The validity of the total inversion results has been evaluated by making use of the matched field processing source localization.     

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.     

Sequential trans-dimensional Monte Carlo for range-dependent geoacoustic inversion

This paper develops a sequential trans-dimensional Monte Carlo algorithm for geoacoustic inversion in a strongly range-dependent environment. The algorithm applies advanced Markov chain Monte Carlo methods in combination with sequential techniques (particle filters) to carry out geoacoustic inversions for consecutive data sets acquired along a track. Changes in model parametrization along the track (e.g., number of sediment layers) are accounted for with trans-dimensional partition modeling, which intrinsically determines the amount of structure supported by the data information content. Challenging issues of rapid environmental change between consecutive data sets and high information content (peaked likelihood) are addressed by bridging distributions implemented using annealed importance sampling. This provides an efficient method to locate high-likelihood regions for new data which are distant and ∕ or disjoint from previous high-likelihood regions. The algorithm is applied to simulated reflection-coefficient data along a track, such as can be collected using a towed array close to the seabed. The simulated environment varies rapidly along the track, with changes in the number of layers, layer thicknesses, and geoacoustic parameters within layers. In addition, the seabed contains a geologic fault, where all layers are offset abruptly, and an erosional channel. Changes in noise level are also considered.     

Theoretical framework for geoacoustic inversion by adjoint method

Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new physical framework for geoacoustic retrievals. A parabolic approximation of wave equation with non-local boundary condition is used as the forward propagation model. The expressions of the corresponding tangent linear model and the adjoint operator are derived, respectively, by variational method. The analytical expressions for the gradient of the cost function with respect to the control variables can be formulated by the adjoint operator, which in turn can be used for optimization by the gradient-based method.     

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.     

利用海底反射信号进行地声参数反演的方法

针对现有反演方法的缺点,提出了一种基于海底反射信号的地声参数高分辨反演方法.它利用短距离声源在不同深度上发射宽带线性调频信号,采用垂直阵进行接收,首先通过匹配滤波方法提取多径到达信息,然后利用海底反射损失曲线,反演海底表层的声速和密度,最后利用浅底层反射信号估计沉积层参数.由于海水中直达波受到内波的强烈影响,选择海底表面反射作为参考,用以可靠地计算浅底层反射的相对到达时间和幅度,从而估计出沉积层的厚度、速度和衰减系数.通过海上实验,验证了利用浅底层反射信号反演参数的有效性. 关键词： 海底参数 反演 浅底层反射信号     

High-resolution geoacoustic inversion in shallow water: a joint time- and frequency-domain technique

High-resolution geoacoustic data are required for accurate predictions of acoustic propagation and scattering in shallow water. Since direct measurement of geoacoustic data is difficult, time-consuming, and expensive, inversion of acoustic data is a promising alternative. However, the main problem encountered in geoacoustic inversion is the problem of uniqueness, i.e., many diverse geoacoustic models can be made to fit the same data set. A key, and perhaps unique, aspect of this approach is the combination of data analysis in both the space-time and the space-frequency domains. This combination attempts to ameliorate the uniqueness problem by exploiting as much independent data as possible. In order to meet the stringent requirements of high spatial resolution and uniqueness, an entire method has been developed including a new measurement technique, processing/analysis technique, and inversion strategy. These techniques are described and then illustrated with a shallow-water data set. Sound-speed gradients in the upper few meters of the sub-bottom appear to be much higher (one order of magnitude) than generally assumed. And, although often ignored, a large density gradient was observed in the top layer that played an acoustically significant role.     

Sequential Bayesian geoacoustic inversion for mobile and compact source-receiver configuration

Geoacoustic characterization of wide areas through inversion requires easily deployable configurations including free-drifting platforms, underwater gliders and autonomous vehicles, typically performing repeated transmissions during their course. In this paper, the inverse problem is formulated as sequential Bayesian filtering to take advantage of repeated transmission measurements. Nonlinear Kalman filters implement a random-walk model for geometry and environment and an acoustic propagation code in the measurement model. Data from MREA/BP07 sea trials are tested consisting of multitone and frequency-modulated signals (bands: 0.25-0.8 and 0.8-1.6 kHz) received on a shallow vertical array of four hydrophones 5-m spaced drifting over 0.7-1.6 km range. Space- and time-coherent processing are applied to the respective signal types. Kalman filter outputs are compared to a sequence of global optimizations performed independently on each received signal. For both signal types, the sequential approach is more accurate but also more efficient. Due to frequency diversity, the processing of modulated signals produces a more stable tracking. Although an extended Kalman filter provides comparable estimates of the tracked parameters, the ensemble Kalman filter is necessary to properly assess uncertainty. In spite of mild range dependence and simplified bottom model, all tracked geoacoustic parameters are consistent with high-resolution seismic profiling, core logging P-wave velocity, and previous inversion results with fixed geometries.     

深海低声速沉积层简正波用于海底参数反演*

现有深海地声参数反演研究通常基于全波场理论,存在计算量大、多值性以及需要准确的水文环境信息等问题。针对这一问题,提出一种基于简正波频散特征的深海低声速沉积层海底参数反演方法。在南海北部大陆坡海域的一次实验中,坐底式水听器(深度约1740 m)接收的沉底弹信号中观察到一种低频成分先到而高频成分后至的到达结构,由被限制在低声速沉积层的简正波组成。通过匹配沉积层简正波60-220 Hz频段内的到达时间差对实验海域的低声速海底参数进行反演,得到沉积层厚度为16.4 m,沉积层声速为1450.5 m/s,与底质采样数据比较吻合,且具有较高的可信度,验证了所提方法的有效性。由于沉积层密度对简正波频散特征的敏感性较差,其反演结果可信度不是很高,需要进一步研究。     

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.     

一种基于简正波模态频散的远距离宽带海底参数反演方法

在浅海环境中, 海底环境参数对声传播有着重要的影响. 由于利用单个宽带声源进行海底参数反演时, 随着距离的增大, 误差变大, 本文提出利用warping变换对在浅海波导中传播的, 不同距离上的两个宽带爆炸声源进行简正波的有效分离, 实现了宽带爆炸声源的远距离海底参数反演. 采用全局寻优遗传算法对提取出的模态频散到达时间差与理论计算的模态频散到达时间差进行匹配处理, 并结合随距离连续变化的声传播损失, 实现了利用单水听器进行海底参数的反演. 实验结果表明: 运用反演出的海底参数提取模态频散时间差和实测数据提取出的模态频散时间差吻合得较好; 而通过传播损失反演得到的海底衰减系数与频率呈指数关系. 最后, 对反演结果进行了后验概率分析, 并将本组爆炸声源的反演结果用于另一组不同距离上爆炸声源时仍然有效, 来评价反演结果的有效性.     

Spectral technique for measuring phase retardation of a wave plate

We present a technique and algorithm for measuring the phase retardation of a wave plate based on spectral transmission curve. Through accurately extracting the intersection points' wavelengths from the spectral transmission curve, the effective phase retardation, absolute phase retardation, order, and physical thickness of the wave plate can be measured simultaneously in a wide spectral range. Experimental results show that the proposed technique has many advantages, such as higher data utilization, simpler extraction algorithm, and no strict requirement for the directions of transmission axes of the polarizer and analyzer, and the fast axis of the wave plate.     

基于反相位脉冲技术的生物组织谐波成像研究

理论及实验研究了反相位脉冲技术在生物组织二次谐波成像中的应用。结合有限振幅声波的非线性传播理论,从理论上证明了反相位脉冲技术可有效抑制基波信号,同时二次谐波信号增强两倍,轴向及径向声场的实验测量结果基波被抑制30~50dB,二次谐波提高6dB,与理论相符。建立了相应的成像系统,对若干生物离体组织进行了基于反相位脉冲相位技术的二次谐波成像,并与常规的基波及二次谐波图像对比,进一步证明了该技术能有效提高图像的对比度和清晰度。     

A Bayesian method of seabed structure and geoacoustic parameters inversion using ship radiated noise

Aiming at the problem of inversion shallow seabed layered structure and geoacoustic parameter inversion with ship radiated noise as the reference sound source,a shallow sea multi-layered seabed geoacoustic parameter inversion method based on Bayesian theory is studied.In the inversion,the line spectrum component of the ship radiated noise is the research object,and then the nonlinear Bayesian inversion method is used to invert the seabed structure,the wave speed in the layer,the wave speed atten...     

A statistical geoacoustic inversion scheme based on a modified radial basis functions neural network

This paper addresses the task of recovering the geoacoustic parameters of a shallow-water environment using measurements of the acoustic field due to a known source and a neural network based inversion process. First, a novel efficient "observable" of the acoustic signal is proposed, which represents the signal in accordance with the recoverable parameters. Motivated by recent studies in non-Gaussian statistical theory, the observable is defined as a set of estimated model parameters of the alpha-stable distributions, which fit the marginal statistics of the wavelet subband coefficients, obtained after the transformation of the original signal via a one-dimensional wavelet decomposition. Following the modeling process to extract the observables as features, a radial basis functions neural network is employed to approximate the vector function that takes as input the observables and gives as output the corresponding set of environmental parameters. The performance of the proposed approach in recovering the sound speed and density in the substrate of a typical shallow-water environment is evaluated using a database of synthetic acoustic signals, generated by means of a normal-mode acoustic propagation algorithm.     

Joint time/frequency-domain inversion of reflection data for seabed geoacoustic profiles and uncertainties

This paper develops a joint time/frequency-domain inversion for high-resolution single-bounce reflection data, with the potential to resolve fine-scale profiles of sediment velocity, density, and attenuation over small seafloor footprints (approximately 100 m). The approach utilizes sequential Bayesian inversion of time- and frequency-domain reflection data, employing ray-tracing inversion for reflection travel times and a layer-packet stripping method for spherical-wave reflection-coefficient inversion. Posterior credibility intervals from the travel-time inversion are passed on as prior information to the reflection-coefficient inversion. Within the reflection-coefficient inversion, parameter information is passed from one layer packet inversion to the next in terms of marginal probability distributions rotated into principal components, providing an efficient approach to (partially) account for multi-dimensional parameter correlations with one-dimensional, numerical distributions. Quantitative geoacoustic parameter uncertainties are provided by a nonlinear Gibbs sampling approach employing full data error covariance estimation (including nonstationary effects) and accounting for possible biases in travel-time picks. Posterior examination of data residuals shows the importance of including data covariance estimates in the inversion. The joint inversion is applied to data collected on the Malta Plateau during the SCARAB98 experiment.     

Selection of spiral waves in excitable media with a phase wave at the wave back

Universal relationships between the medium excitability and the angular velocity and the core radius of rigidly rotating spiral waves in excitable media are derived for situations where the wave front is a trigger wave and the wave back is a phase wave. Two universal limits restricting the region of existence of spiral waves in the parameter space are demonstrated. The predictions of the free-boundary approach are in good quantitative agreement with results from numerical reaction-diffusion simulations performed on the Kessler-Levine model.