南海北部负跃层环境下海底参数声学反演

利用2015年南海北部声学实验的实验数据,对海底声学参数进行反演。在综合分析声速起伏及海底模型对海底反演影响的基础上,通过选择等效的海水声速剖面和海底模型,改进了多参量联合地声反演方法,使得其不仅能解决反演中的多值性问题,还能适用于负跃层起伏环境下的声学反演。即当用匹配场反演海底声速和用传播损失反演海底衰减时,如果温跃层内有内波等随机起伏,可使用传播路径上平均声速剖面和单层等效海底进行声场计算。反演得到的海底声速和密度结果与海底采样测量符合较好,拟合给出了海底衰减系数随频率的非线性经验关系式。反演结果可为南海北部声传播规律研究与应用提供海底参数。      

一种低声速沉积层海底参数声学反演方法

软泥底环境下沉积层参数的声学反演是国际水声领域的一个研究热点.浅海中,当高声速基底和海水之间存在一层低声速(小于海水声速)的沉积层时,小掠射角情况下不同频率声传播损失会出现周期性增大现象.基于此现象,提出一种适用于低声速沉积层的海底参数声学反演方法.首先,推导给出小掠射角情况下传播损失周期增大的频率间隔与沉积层声速、厚度及近海底海水声速之间的解析表达式;其次,利用一次黄海实验中软泥底环境下的宽带声传播信号,提取了小掠射角下传播损失增大的频率周期;再次,把该解析表达式作为约束条件,结合Hamilton密度与声速的经验公式,采用匹配场处理反演给出沉积层的声速、密度、厚度及基底的声速、密度;然后,利用声传播损失数据反演得到泥底环境下不同频率的声衰减系数,通过拟合发现泥底声衰减系数随频率近似呈线性关系;最后,给出了双层海底模型和半无限大海底模型等效性的讨论.反演结果为低声速沉积层海底声传播规律研究与应用提供了海底声学参数.     

基于大掠射角海底反射特性的深海地声参数反演

海底声学参数的获取对于海洋声学研究有着重要意义.通过推导分层吸收介质下的海底反射系数,理论分析了大掠射角条件下吸收系数对海底反射系数的影响.海底反射系数随频率振荡过程中,将其等于海水-沉积层界面反射系数模时所对应的频点定义为1/4振荡周期频率.在该频率下,沉积层吸收系数与基底地声参数的耦合程度小于其他频点.本文基于大掠射角下的海底反射特性,提出一种深海地声参数分步反演方法.首先,利用相关法提取得到海底反射系数的干涉周期,利用干涉周期反演了沉积层声速和厚度.声速的反演结果结合Hamilton经验公式反演密度.第二步,通过结合基底声速的穷举边界,给出沉积层吸收系数的假设值,利用1/4振荡周期频率下的海底反射系数对基底声速进行一维反演.最后利用半波层频率下的海底反射系数对沉积层吸收系数进行一维反演.大掠射角海底反射特性结合分步反演,实现了基底声速和沉积层吸收系数一定程度的解耦合.实验结果表明,在大掠射角测量条件下,该方法反演的地声参数可有效应用在一定范围内的传播损失预报.     

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

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

舰船辐射噪声贝叶斯方法反演浅海底层结构及地声参数EI北大核心CSCD

针对以舰船辐射噪声为参考声源的浅海海底分层结构及地声参数反演问题,研究了一种基于贝叶斯理论的浅海多层海底地声参数反演方法。反演中以舰船辐射噪声的线谱成分为研究对象,进而采用非线性贝叶斯反演方法反演浅海底层结构、层中声速、声速衰减和密度,并对反演结果的不确定性进行分析。反演结果的最大后验概率估计值和边缘概率分布分别通过拨正模拟退火算法和Metropolis-Hastings采样法在各参数先验区间内计算获得,并根据贝叶斯信息准则确定最佳海底分层结构。海上实验表明:根据该方法反演获得海底分层结构及地声参数,计算得到的声压场与实测舰船辐射噪声传播损失误差不超过10%,反演结果能够准确表征实验海区海底特征。反演结果不确定性分析表明:海底纵波声速、横波声速以及密度的不确定性更小,对声压场变化更加敏感,反演结果更有效、准确。     

利用气枪声源数据的地声参数反演*

本文通过对东海采集的气枪声源数据进行处理,基于地声参数对不同声场观测量的敏感性差异,利用分步反演的策略进行分析。首先选择距离接收阵七个不同的距离点数据,通过warping变换准确提取第三阶至第八阶简正波的频散曲线;利用海底声速对频散曲线敏感的特性来反演声速;由Hamilton经验公式求得海底密度;通过传播损失拟合获得海底衰减。不同距离点数据反演的海底声速与密度一致性较高。实验提取的频散曲线和和反演参数仿真结果、实验获取的传播损失与反演参数获得传播损失均拟合较好。     

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

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

利用随机多项式展开的海底声学参数反演方法

为了提高地声反演算法的计算效率,探索克服地声反演结果多值性问题,本文利用宽带、多收发位置的传播损失数据结合传播损失在地声参数先验搜索区间内的随机多项式展开系数矩阵,反演得到海底纵波声速、吸收率和密度比重.使用随机多项式展开近似传播损失时,展开系数的自变量为声波频率、收发位置等参数,随机多项式的自变量为表示声速、吸收率、比重在各自搜索区间内均匀分布的随机变量.传播损失的展开系数通过嵌入随机多项式的声学宽角抛物方程结合盖辽金投影、最小角度回归算法计算求得.在低频、一定声传播水平距离以内和地声参数搜索区间长度适中时,使用随机多项式展开近似传播损失的相对误差在1%以下.仿真发现,在浅海环境中使用低频、一定声传播水平距离以内的传播损失数据,在接收信号信噪比较高、声源和水听器相对位置误差较小时,选择合适的随机多项式展开截断幂次可较准确地反演海底声速、吸收率和密度比重,且计算效率比网格遍历搜索方法提高一个数量级以上.     

浅海环境噪声垂直指向性探测海底参数

浅海环境噪声垂直指向性和海底地声参数密切相关,可以利用它来反演海底参数。本文在水平分层介质噪声场模型下,讨论了海底参数对噪声垂直指向性的敏感度,并假设半无限大液态海底,采用效率较高的自适应单纯形模拟退火算法全局寻优,根据中国某海区实测噪声数据,估计了海底声速、密度及衰减,分析了反演结果的不确定性。鉴于声速和密度在不同频率的反演结果不够一致,文中提出多频海底参数联合反演方法,提高了反演海底声学参数的精度。     

由简正波群延时及幅度反演海底参数

由简正波过滤技术得到1号与2号简正波群延时及幅度。通过匹配简正波群延时差来反演海底的声速、密度等参数,衰减系数由简正波的幅度反演得到。遗传算法用来在参数域内搜索待反演的参数。在 200～500 Hz频率范围内海底衰减系数随频率的变化关系为α=0.29f1.9dB/m（f的单位为 kHz）。最后给出由反演出的参数计算得到的传播损失与实验的声传播损失数据的比较,结果符合的较好。     

Geoacoustic inversion for bottom parameters in a thermocline environment in the northern area of the South China Sea

Bottom acoustic parameters have important influence on the application of underwater acoustic propagation and source location.The acoustic parameters of the seabed in the northern of the South China Sea(SCS) were inversed using the experiment data from an acoustic experiment in 2015.Based on the comprehensive analysis of the influence of the sound speed fluctuation and the geoacoustic model on seabed inversion,the multi-parameter hybrid acoustic inversion scheme is improved by selecting the equivalent mean sound speed profile(SSP) and half-infinite liquid bottom model to save the inversion dimensions in the matched field processing(MFP) inversion.The inverted bottom sound speed and density are in good agreement with the core sampling measurements.The nonlinear empirical relationship of the attenuation coefficient with frequency is given out.The inversion results are meaningful to the sound propagation research and application in the northern area of the SCS.     

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 effects of ignored seabed variability in geoacoustic inversion

Geoacoustic inversion using a matched-field inversion algorithm is a well-established technique for estimating the geoacoustic parameters of the seabed. This paper demonstrates how parameter estimation can be affected by unknown or wishfully ignored random range dependence of the true environment when the inversion model is--for practicality--assumed to be range independent. Simulations with controlled statistics were carried out using a simple shallow water model: an isospeed water column over a homogeneous elastic halfspace. The inversion parameters included water depth, compressional speed in the seabed, seabed density, and compressional wave attenuation. On average the environment is range independent: some parameters are constant while other parameters are random with range-independent means and variances. A Parabolic Equation underwater acoustic propagation model is used to calculate the simulated data fields for the range-dependent environment as well as to calculate the model fields for the range-independent inversion model. The Adaptive Simplex Simulated Annealing inversion algorithm is used to estimate the best-fit solution. It is found that ignoring the variability of even a single geoacoustic parameter leads to significant and correlated uncertainty (bias and variance) in the estimation of all inverted parameters. Results are presented for range variation of compressional sound speed and water depth.     

Full wave-field reflection coefficient inversion

This paper develops a Bayesian inversion for recovering multilayer geoacoustic (velocity, density, attenuation) profiles from a full wave-field (spherical-wave) seabed reflection response. The reflection data originate from acoustic time series windowed for a single bottom interaction, which are processed to yield reflection coefficient data as a function of frequency and angle. Replica data for inversion are computed using a wave number-integration model to calculate the full complex acoustic pressure field, which is processed to produce a commensurate seabed response function. To address the high computational cost of calculating short range acoustic fields, the inversion algorithms are parallelized and frequency averaging is replaced by range averaging in the forward model. The posterior probability density is interpreted in terms of optimal parameter estimates, marginal distributions, and credibility intervals. Inversion results for the full wave-field seabed response are compared to those obtained using plane-wave reflection coefficients. A realistic synthetic study indicates that the plane-wave assumption can fail, producing erroneous results with misleading uncertainty bounds, whereas excellent results are obtained with the full-wave reflection inversion.     

射线理论下的浅海环境噪声场空间相干性分析及海底声学参数反演

基于射线理论,以简化海底模型为基础,获得了浅海噪声场空间相干性在频域上的表示方法,并用于海底声学参数反演。将风关海洋噪声看作分布在无限大平面上的点源,通过计算噪声场垂直方向上的能量分布,进一步获得空间两点噪声场频域相干性表达式。利用数值仿真验证了模型对于不同海底类型的敏感性。并结合噪声场空间两点的频域相干性海上实验数据,反演得到该处海域的底质类型,进一步获得海底声速、密度和声衰减系数。结果表明,反演得到的海底参数与实际情况符合较好,文中提出采用的这种计算方法能够有效的反映海底反射对于噪声场空间相干性的影响,并且可在一定程度上提高反演海底声学参数效率。      

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.     

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 geoacoustic inversion using wind-driven ambient noise

This paper applies Bayesian inversion to bottom-loss data derived from wind-driven ambient noise measurements from a vertical line array to quantify the information content constraining seabed geoacoustic parameters. The inversion utilizes a previously proposed ray-based representation of the ambient noise field as a forward model for fast computations of bottom loss data for a layered seabed. This model considers the effect of the array's finite aperture in the estimation of bottom loss and is extended to include the wind speed as the driving mechanism for the ambient noise field. The strength of this field relative to other unwanted noise mechanisms defines a signal-to-noise ratio, which is included in the inversion as a frequency-dependent parameter. The wind speed is found to have a strong impact on the resolution of seabed geoacoustic parameters as quantified by marginal probability distributions from Bayesian inversion of simulated data. The inversion method is also applied to experimental data collected at a moored vertical array during the MAPEX 2000 experiment, and the results are compared to those from previous active-source inversions and to core measurements at a nearby site.