Single-receiver geoacoustic inversion using modal reversal

This paper introduces a single-receiver geoacoustic-inversion method based on dispersion analysis and adapted to low-frequency impulsive sources in shallow-water environments. In this context, most existing methods take advantage of the modal dispersion curves in the time-frequency domain. Inversion is usually performed by matching estimated dispersion curves with simulated replicas. The method proposed here is different. It considers the received modes in the frequency domain. The modes are transformed using an operator called modal reversal, which is parameterized using environmental parameters. When modal reversal is applied using parameters that match the real environment, dispersion is compensated for in all of the modes. In this case, the reversed modes are in phase and add up constructively, which is not the case when modal reversal is ill-parameterized. To use this phenomenon, a criterion that adds up the reversed modes has been defined. The geoacoustic inversion is finally performed by maximizing this criterion. The proposed method is benchmarked against simulated data, and it is applied to experimental data recorded during the Shallow Water 2006 experiment.     

Particle filtering for dispersion curve tracking in ocean acoustics

A particle filtering method is developed for dispersion curve extraction from spectrograms of broadband acoustic signals propagating in underwater media. The goal is to obtain accurate representation of modal dispersion which can be employed for source localization and geoacoustic inversion. Results are presented from the application of the method to synthetic data, demonstrating the potential of the approach for accurate estimation of waveguide dispersion characteristics. The method outperforms simple time-frequency analysis providing estimates that are very close to numerically calculated dispersion curves. The method also provides uncertainty information on modal arrival time estimates, typically unavailable when traditional methods are used.     

Geoacoustic inversion in a dispersive waveguide using warping operators

This paper presents a single receiver geoacoustic inversion method adapted for low-frequency impulsive sources. It is applied to light bulb data collected during the Shallow Water 2006 experiment. The inversion is carried out by extracting dispersion curves from the received signal, and comparing them to simulated replicas. To achieve dispersion curve estimation in the time-frequency domain, modal separability is improved using a signal processing method called warping. The inversion scheme allows for a reliable estimation of the New Jersey Shelf sediment properties (compressional sound speed and density). It also provides an accurate estimation of the source/receiver range.     

地声参数及传播损失不确定性估计与建模

地声参数的不确定性对水声传播具有重要的影响。通过贝叶斯理论建立水声环境不确定性推理模型,理论推导了地声参数的似然函数以及地声参数和传播损失的后验概率密度,并采用MCMC(Markov Chain Monte Carlo)进行了仿真计算,给出了地声参数的二维后验联合概率密度和一维边缘概率密度,在此基础上对传播损失的不确定性进行了估计,得到了传播损失80%的可信区间。仿真和实验结果表明,该方法适用于地声参数反演和不确定性估计,并能获取因地声参数不确定性导致的传播损失不确定性估计。     

一种基于简正波模态消频散变换的声源距离深度估计方法

针对浅海环境中传播的低频宽带水声脉冲信号,基于简正波水平波数差和波导不变量之间的关系,本文提出了一种利用距离-频散参数二维平面聚焦测距与匹配模态能量定深的目标声源定位方法.首先,通过将由频散参数和波导不变量表示的前几阶模态相速度与由环境模型计算的相速度进行对比分析,从而估计出前几阶模态的频散参数和环境的波导不变量.其次,利用估计出的频散参数值和波导不变量对接收信号进行消频散变换处理,只有当接收信号的距离参数等于目标声源距离时,各号简正波的幅度均达到最大值,在距离-频散参数二维平面上,出现声压聚焦的现象,利用此现象可以估计目标声源的距离.不仅如此,消频散变换后的接收信号,前几阶模态在时域上明显地分离开来,可以准确地估计出前几阶模态的能量,采用多模态能量匹配的方式,可以估计出目标声源的深度.最后,通过对仿真和冬季获得的气枪信号数据处理结果验证了本文方法的有效性.     

Geoacoustic inversion based on reflection model of effective density fluid approximation

In order to obtain the physical and geoacoustic properties of marine sediments,an inverse method using reflection loss of different grazing angles is presented.The reflection loss is calculated according to the reflection model of effective density fluid approximation.A two-step hybrid optimization algorithm combining differential evolution and particle swarm optimization along with Bayesian inversion is employed in estimation of porosity,mean grain size,mass density and bulk modulus of grains.Based on the above physical parameters,geoacoustic parameters,including sound speed and attenuation,are further calculated.According to the numerical simulations,we can draw a conclusion that all the parameters can be well estimated with the exception of bulk modulus of grains.In particular,this indirect inverse method for bottom geoacoustic parameters performs high accuracy and strong robustness.The relative errors are 0.092%and 17%,respectively.Finally,measured reflection loss data of sandy sediments at the bottom of a water tank is analyzed,and the estimation value,uncertainty and correlation of each parameter are presented.The availability of this inverse method is verified through comparison between inverse results and part of measured parameters.     

利用多角度海底反向散射信号进行地声参数估计

针对现有海底地声参数估计方法的不足,提出了利用相控参量阵浅地层剖面仪接收的多角度海底反向散射信号进行地声参数估计的方法,首先利用正下方和斜入射方向上沉积层上、下表面的差频反向散射信号进行沉积层厚度和声速估计,然后利用正下方沉积层上、下表面两个不同频率的差频信号的反向散射信号估计沉积层衰减系数,最后利用正下方沉积层上表面原频反向散射信号估计沉积层阻抗,计算沉积层密度从而解决和声速的耦合性,通过水池试验验证了该方法的有效性。     

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.     

基于等效密度流体近似反射模型反演海底参数

为了获取海底沉积物的物理和地声特性,根据等效密度流体近似反射模型得到的不同掠射角下的海底反射损失,利用差分进化算法和粒子群算法相结合的两级混合优化算法以及Bayesian反演方法对海底沉积物的孔隙度、平均颗粒粒度、颗粒质量密度以及颗粒体积弹性模量进行反演,再根据这4个物理参数的估计值进一步计算出海底地声参数,包括声速和衰减。通过反演结果与仿真真值的比较可以看出,除颗粒体积弹性模量外,得到的估计结果是令人满意的,特别是这种间接反演海底地声参数的方法对地声参数的估计具有较高的精确性和较强的稳健性,相对误差分别为0.092%和17%。最后,对实验室水池池底沙质沉积物的反射损失数据进行处理,给出了沙质沉积物各参数的估计值、不确定性和相关性,并通过反演结果与部分实测参数值的比较验证了反演方法的可行性。      

Passive geoacoustic inversion with a single hydrophone using broadband ship noise

An inversion scheme is proposed, relying upon the inversion of the noise of a moving ship measured on a single distant hydrophone. The spectrogram of the measurements exhibits striations which depend on waveguide parameters. The periodic behavior of striations versus range are used to estimate the differences of radial wavenumber between couples of propagative modes at a given frequency. These wavenumber differences are stacked for several frequencies to form the relative dispersion curves. Such relative dispersion curves can be synthesized using a propagation model feeded with a bottom geoacoustic model. Inversion is performed by looking for the bottom properties that optimize the fit between measured and predicted relative dispersion curves. The inversion scheme is tested on simulated data. The conclusions are twofold: (1) a minimum 6 dB signal to noise ratio is required to obtained an unbiased estimate of compressional sound speed in the bottom with a 3 m s(-1) standard deviation; however, even with low signal to noise ratio, the estimation error remains bounded and (2) in the case of a multi-layer bottom, the scheme produces a single depth-average compressional sound speed. The inversion scheme is applied on experimental data. The results are fully consistent with a core sample measured around the receiving hydrophone.     

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

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

基于模式能量比的海底声衰减系数反演

简正波的分离性能对简正波频散曲线、模式能量比等物理量的提取至关重要,这些物理量可以用来反演海底地声参数。本文分析了短时傅里叶变换(STFT)和warping变换两种方法在简正波分离上的性能的差异,并利用该两种方法提取模式能量比来反演海底声衰减系数。仿真实验和2018年冬季北海实验数据处理表明,warping变换在简正波分离上具有较大的优势,尤其是在近距离接收上。     

Geoacoustic Inversion Based on Dispersion Characteristic of Normal Modes in Shallow Water

A geoacoustic inversion method based on dispersion characteristic of normal modes is presented. An adaptive time-frequency analysis technique with a high resolution in both time and frequency domains is applied to derive the dispersion characteristic of normal modes from the broadband propagation signal. The bottom acoustic parameters are inverted by matching the calculated group delays of normal modes with the experimental data. Finally, some experimental results which could validate the inversion method are given.     

Parametric Analysis of Complex Dispersion Curves for Flexural Lamb Waves in Layered Plates in the Low-Frequency Range

A combined asymptotical and iteration method is used to study dispersion curves for the case of dynamic bending of isotropically layered plates. Based on the explicit limit formulation of dispersion equation, asymptotics of roots are derived in closed form for large values of root moduli. The influence of elastic and geometric parameters of layers are analyzed. The existence of critical values of geometric parameters that correspond to change of the type of asymptotics is demonstrated. The errors of asymptotics are estimated, and an iterative method is proposed for calculating the exact values of roots in statics. A low-frequency long-wave asymptotics of complex dispersion curves is derived; its accuracy is the higher the lower the frequency and the greater the number of the curve are. It is also proved that each complex curve has a long flat segment, the length of which increases simultaneously with the number of curve. The dispersion curves themselves are also calculated by another specific iterative procedure. The fundamental bending mode is analyzed together with its purely imaginary sister. The existence of the additional purely imaginary curve at low frequency is proved. Examples of calculating the static roots and the dispersion curves for subcritical and supercritical values of geometrical parameters are presented, and the efficiency of the algorithm is estimated.     

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.     

Inversion of guided-wave dispersion data with application to borehole acoustics

The problem of inferring unknown geometry and material parameters of a waveguide model from noisy samples of the associated modal dispersion curves is considered. In a significant reduction of the complexity of a common inversion methodology, the inner of two nested iterations is eliminated: The approach described does not employ explicit fitting of the data to computed dispersion curves. Instead, the unknown parameters are adjusted to minimize a cost function derived directly from the determinant of the boundary condition system matrix. This results in an efficient inversion scheme that, in the case of noise-free data, yields exact results. Multimode data can be simultaneously processed without extra complications. Furthermore, the inversion scheme can accommodate an arbitrary number of unknown parameters, provided that the data have sufficient sensitivity to these parameters. As an important application, we consider the sonic guidance condition for a fluid-filled borehole in an elastic, homogeneous, and isotropic rock formation for numerical forward and inverse dispersion analysis. We investigate numerically the parametric inversion with errors in the model parameters and the influence of bandwidth and noise, and examine the cases of multifrequency and multimode data, using simulated flexural and Stoneley dispersion data.     

椭圆芯手征光纤的基模特性研究

对椭圆芯手征光纤进行了严格的求解，得到了模式场的解析解；推了了模式特征方程；并对基模的偶模和奇模特征方程进行数值求解，研究了手征参数对不同椭圆比下基模的色散曲线，双折射，模间色散等特性的影响。     

利用拖船自噪声进行浅海环境参数贝叶斯反演

研究了以拖船自噪声为参考声源的浅海环境参数反演问题,并针对反演结果不确定性快速量化评估问题,提出了一种基于自适应重要性抽样的贝叶斯反演新方法。反演利用了拖船自噪声低频线谱成分,并采用混合高斯推荐函数自适应推荐声场模型样本,使得样本集中于参数高概率密度区域,实现后验概率密度快速收敛计算。仿真试验结果表明:拖船自噪声反演能够准确估计水深、沉积层及阵列参数等。所提自适应重要性抽样贝叶斯反演方法的计算效率优于快速吉布斯抽样方法。利用试验数据处理验证,反演得到试验海域声学环境参数,计算传播损失与各阵元接收线谱强度变化吻合,说明反演最优环境模型能准确表征声场传播特征。      

Computing Eigen value equation and modal dispersion characteristics of an elliptical Bragg waveguide

Using the method of separation of variables in the elliptical coordinate system, a recursive formula for the electromagnetic fields in a new type of Bragg waveguide having elliptical core cross section with multilayered cladding is derived. The eigen equation is written in the form of Mathieu and the modified Mathieu functions and a dispersion relation is obtained for various modes supported by the proposed Bragg waveguide. The cutoff frequencies for several lower order even-odd modes have been calculated and their propagation characteristics are plotted. The results show that the dispersion curves are discontinuous and modes can exist only in particular wavelength bands. The effects of elliptical eccentricity on the mode cutoff values and mode transmission are addressed. Finally, the modal birefringence in the said waveguide is also estimated.     

基于等效密度流体模型的地声参数分析

目前的地声参数反演多采用液态海底模型,但是实际海底为多孔弹性海底。该文在等效密度流体模型基础上,通过计算液态海底和多孔弹性海底的反射系数及传播损失,给出了等效密度流体模型和液态海底模型的等效性分析。数值仿真结果表明在低频情况下,多孔弹性海底给出的快纵波声速与等效密度流体模型给出的声速以及等效密度流体的实部与真实的海底密度基本一致。将等效密度流体模型近似看作液态海底模型进行反射系数和传播损失计算,在小掠射角和远距离时,计算结果表明与多孔弹性海底计算结果具有较好一致性,从而以此为依据确定海底为液态,进行地声参数反演。