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

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

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

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

浅海地声参数的变维贝叶斯反演

为解决海底沉积层分层结构未知时的反演问题,提出了一种变维粒子滤波方法,利用声场的互谱密度,估计沉积层分层结构以及地声参数。仿真结果表明：变维粒子滤波在沉积层分层结构未知时,能有效反演沉积层层数以及地声参数,粒子的并行计算能使其相较于可逆跳蒙特卡洛马尔可夫链(rjMCMC)更加高效。利用变维粒子滤波,对南海垂直线阵列接收到的线性调频信号进行处理,反演结果与rjMCMC反演得到的沉积层层数和地声参数结果相近,说明了此方法能有效估计沉积层层数的同时反演浅海地声参数,得到可靠的参数后验概率密度。     

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

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

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

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

深海近距离声场频率-距离干涉结构反演海底声学参数

地声反演一般存在多解性的问题,不易有效分离沉积层声速与厚度,本文提出了一种当声源与接收水听器分别位于海面与海底附近时匹配声场干涉周期反演深海沉积层声速与厚度的方法。该声场周期结构由直达波与海底沉积层反射回波干涉形成,其干涉周期受沉积层声速和厚度影响,通过匹配实测与理论计算的声场干涉周期可有效反演沉积层声速和厚度。开展了深海实验验证,利用水听器接收到的近距离船舶噪声的声场干涉条纹,反演得到的实验海域海底沉积层声速与该海域海底表层采样平均值差小于20 m/s。数值仿真与实验结果表明该方法可以有效反演沉积层的声速与厚度。      

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

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

利用时域波形快速反演海底单参数的方法

海底单参数反演由于反演维数少, 较之多维反演具有实施简单快捷的优点. 首先定义小掠射角下反射损失随掠射角变化率F_dB为单个参数描述海底性质, 构建单参数反演模型; 其次, 从简正波与射线理论出发, 对时域脉冲波形与海底反射性质之间的关系进行了推导, 并讨论衰减机制; 最后, 基于该关系设计出利用最小二乘法对时域波形进行反演的方法. 对北黄海实验数据进行了实际反演, 通过海底采样验证反演结果对海底性质估计的正确性, 同时基于反演结果的传播损失预报值与实测值符合良好. 该方法实施仅需要单个水听器, 待测物理量少, 所获结果较为可靠, 且能满足大多数应用的需要. 关键词： 单参数 地声反演 浅海 海底反射损失     

利用拖船辐射噪声的深海声学参数近场匹配场反演方法

针对深海声学参数难以通过远距离合作声源反演获取的问题,提出了利用拖船低频噪声近场匹配场反演方法。首先,利用聚焦波束形成计算拖曳阵接收拖船噪声的方向性,获得传播路径特征;然后,构建多参数反演模型,由波数积分声传播模型计算拷贝场,采用遗传算法对多频匹配场目标函数进行反演。同时,采用蒙特卡罗方法分析参数后验概率密度。仿真与试验结果表明:深海环境中拖曳阵接收拖船噪声主要来自海底反射路径,利用该特性反演得到海水深度、噪声源距离、阵列深度、沉积层厚度等参数,多频联合反演可以提高沉积层厚度等参数反演准确性。宽带匹配场处理表明,利用反演最优参数模型能准确给出拖船噪声源的空间位置。      

利用传播损失反演海底单参数

根据地声反演复杂性随着待反演参数的个数减少不断降低的原理,提出一种利用传播损失反演海底单个参数的方法。通过对海底声阻抗的推导,利用声速、密度和衰减系数拟合出海底反射损失对掠射角的斜率F。基于简正波理论推导了用F描述声场的公式,并据此设计出对传播损失数据进行最小二乘法处理的反演方法。得益于将待反演参数减少至一个,该方法只需单个水听器,避免了复杂测量及多维寻优。对东海实验数据进行了实际反演,并介绍了利用反演结果F进行传播损失预报和海底性质估计的步骤。所获结果与多参数混合反演方法及实测真值进行比较,其一致性验证了方法的有效性。      

Bayesian geoacoustic inversion in a dynamic shallow water environment

This paper presents results for matched field Bayesian geoacoustic inversion of multitonal continuous wave data collected on the New Jersey continental shelf. To account for effects of significant spatial and temporal variation of the water column sound speed, the sound speed profile was represented by empirical orthogonal functions. Data error information for the inversion was estimated from multiple time windows of the data. Inversion results for the sediment sound speeds at three ranges are in excellent agreement with the ground truth.     

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

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

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.     

Quantifying the uncertainty of geoacoustic parameter estimates for the New Jersey shelf by inverting air gun data

This paper describes geoacoustic inversion of low frequency air gun data acquired during an experiment on the New Jersey shelf. Hybrid optimization and Bayesian inversion techniques based on matched field processing were applied to multiple shots from three air gun data sets recorded by a vertical line array in a long-range shallow water geometry. For the Bayesian inversions, full data error covariance matrix was estimated from a set of consecutive shots that had high temporal coherence and small spatial variation in source position. The effect of different data error information on the geoacoustic parameter uncertainty estimates was investigated by using the full data error covariance matrix, a diagonalized version of the full error covariance, and a diagonal matrix with identical variances. The comparison demonstrated that inversion using the full data error information provided the most reliable parameter uncertainty estimates. The inversions were highly sensitive to the near sea floor geoacoustic parameters, including sediment attenuation, of a simple single-layer geoacoustic model. The estimated parameter values of the model were consistent with depth averaged values (over wavelength scales) of a high resolution geoacoustic model developed from extensive ground truth information. The interpretation of the frequency dependence of the estimated attenuation is also discussed.     

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.     

Subbottom profiling using a ship towed line array and geoacoustic inversion

Bottom profiling traditionally uses broadband signals received on a line array at long ranges to estimate the bottom layer structure and thickness. In this paper, a subbottom profiling method is developed and applied to a ship-towed line array using the same ship towed source to estimate the subbottom layer structure and thickness. A ship towed line-array system can be used to estimate bottom properties using geoacoustic inversion and can cover a wide area in a short time. It needs some prior information about the subbottom structure and layer thickness, without which the solution can be ambiguous and even erratic when resolving parameters over a wide area. It is shown that the required subbottom information can be obtained from the time-angle relation by beamforming the same acoustic signal data used for geoacoustic inversion. The time-angle analysis is used to expose the prevalent physics intrinsic to geoacoustic inversion. One finds that the tau-p relation of the bottom and the bottom reflection coefficients, sampled at discrete angles associated with bottom and multiple surface-bottom returns, are often adequate, for many practical applications, to uniquely determine the geoacoustic bottom at low (< or =1 kHz) frequencies.