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

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

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

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

一种海洋环境噪声分步反演地声参数方法

分析不同海底参数对环境噪声垂向空间特性的敏感度,根据海底密度对环境噪声垂直指向性的小掠射角部分不敏感,而对等效海底损失的大掠射角部分相对敏感的特点,提出了一种海洋环境噪声分步反演地声参数方法:先用环境噪声垂直指向性小掠射角部分反演海底声速、衰减;之后利用大掠射角部分来反演海底密度。仿真算例和海上实验数据处理结果验证了该方法的有效性。结果表明:分步反演加强了匹配物理量对海底参数敏感部分的关注,在反演结果精度(海底密度)以及反演效率上都有所改善,具有一定可行性。      

Mid-frequency geoacoustic inversion using bottom loss data from the Shallow Water 2006 Experiment

Geoacoustic inversion work has typically been carried out at frequencies below 1 kHz, assuming flat, horizontally stratified bottom models. Despite the relevance to Navy sonar systems many of which operate at mid-frequencies (1-10 kHz), limited inversion work has been carried out in this frequency band. This paper is an effort to demonstrate the viability of geoacoustic inversion using bottom loss data between 2 and 5 kHz. The acoustic measurements were taken during the Shallow Water 2006 Experiment off the coast of New Jersey. A half-space bottom model, with three parameters density, compressional wave speed, and attenuation, was used for inversion by fitting the model to data in the least-square sense. Inverted sediment sound speed and attenuation were compared with direct measurements and with inversion results using different techniques carried out in SW06. Inverted results of the present work are consistent with other measurements, considering the known spatial variability in this area. The observations and modeling results demonstrate that forward scattering from topographical changes is important at mid-frequencies and should be taken into account in sound propagation predictions and geoacoustic inversion. To cope with fine-scale topographic variability, measurement technique such as averaging over tracks may be necessary.     

Complex reflection phase gradient as an inversion parameter for the prediction of shallow water propagation and the characterization of sea-bottoms

In this paper a quantity is proposed, referred to as the complex reflection phase gradient, whose use in a matched field inversion procedure allows for the rapid extraction of first order geo-acoustic information about the sea-bottom. It is based on the observation that at low grazing angles the reflection phase and bottom loss for a wide range of sea-bottom types commonly exhibits an approximate linear relationship to the vertical component of the acoustic wave number at the seabed. The real part of this quantity specifies the rate at which the reflection phase varies with vertical acoustic wave number while the imaginary part quantifies the rate of change of bottom loss. Despite being defined with just two real parameters it is shown that it provides an accurate prediction of the sound field for a wide range of bottom types. In addition, its measurement permits an estimate to be made for the input impedance to the seabed in the zero grazing angle limit and, in the case of a homogeneous elastic half-space of known density, the compressional and shear wave speed. The main advantage of the two-parameter seabottom representation is that each parameter is readily inverted from comparatively few acoustic pressure measurements. The usefulness of the technique is illustrated by the results from computer simulated acoustic pressure measurements made at just eleven sensors in a simple shallow water channel, and results from a 10 cm deep laboratory channel at frequencies between 10 kHz and 75 kHz.     

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

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

浅海海底反射系数幅值参数的反演

理论分析了一种通过混响强度衰减特性获取海底反射系数的幅值参数的方法.将海底反射系数的幅值参数和相位参数引入到全波动混响模型中,为海底反射系数的反演提供理论基础。理论分析和数值仿真表明,在小掠射角条件下,利用混响强度衰减特性反演海底反射系数幅值参数的可行性和准确性。该反演方法只需要输入4个变量:本地混响强度的衰减特性,反射系数的相位参数,海深以及海深处的声速,同时要求混响数据具有一定的混响噪声比(大于6 dB)才能够使反演结果准确可信。根据本地静态混响实验数据成功反演得到海底反射系数的幅值参数.      

Explicit results for wave scattering and transmission through a rough fluid-fluid interface

This paper describes an analysis of reflection and transmission of acoustic waves from an imperfectly reflecting, rough fluid-fluid interface within the Kirchhoff approximation. It presents the results of calculations of the coherent and diffuse field. This work is motivated by the fact that few explicit results of the characteristics of the scattered and transmitted wave field exist in the literature for this problem. For the problem of interest, the surface reflection coefficient depends at each point upon the local angle between the incident wave and the rough surface. For surfaces with statistically independent local surface position and gradient, coherent field calculations show that the correction to constant reflection coefficient analyses is simply a multiplicative factor that depends upon the statistical characteristics of the surface gradient, sound speed and density ratio across the surface. This multiplicative factor is interpreted as an average reflection or transmission coefficient, <R> and <T>, respectively. The principle differences between these results and constant reflection coefficient analyses occur when waves impinge upon regions with higher sound speeds, where total internal reflection may occur. While the wave characteristics of smooth or constant reflection coefficient surfaces change abruptly in the vicinity of the angle of total internal reflection, the average reflection coefficient exhibits a much smoother dependence upon angle of incidence or sound speed ratio for rough surfaces. It is also shown that the direction of maximum diffuse scattering moves relative to its value were the reflection coefficient constant.     

Determination of the inhomogeneous structure of a medium from its plane wave reflection response,part II: A numerical approximation

This paper is concerned with a standard one dimensional inverse scattering problem: given the reflection response of an unknown inhomogeneous medium for plane waves under normal or oblique incidence, determine its sound speed and density structures. The problem is solved by means of a simple numerical technique which involves only fast Fourier transform operations and numerical integration of ordinary differential equations. Three cases are specifically considered: (a) sound speed is unknown, density is known; (b) sound speed is known, density is unknown; (c) sound speed and density are to be determined simultaneously. Numerical simulations performed on reflection coefficients computed in Part I for a limited band of frequencies lead to accurate reconstructions of the original structures of various media.     

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

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

Hybrid geoacoustic inversion of broadband Mediterranean Sea data

This paper describes an acoustic experiment (PROSIM'97) carried out to investigate inversion for seabed properties at a site off the west coast of Italy where previous acoustic and geophysical studies have been performed. Acoustic fields were measured at a vertical hydrophone array due to a swept-frequency source towed over weakly range-dependent bathymetry. Based on the known geology, the seabed is modeled as a sediment layer overlying a semi-infinite basement with unknown model parameters consisting of the sediment thickness, sediment and basement sound speeds, source range and depth, water depth at the source and array, and array tilt. A hybrid inversion algorithm is applied to determine the model values that minimize the mismatch with the measured acoustic fields. Multiple data sets are analyzed to examine the consistency of the inversion results. It is found that the low sound speed of the sediment layer, together with a large uncertainty in bathymetry, leads to strong correlations between the water depths and sediment thickness. This precludes reliable estimation of these parameters individually; however, the total depth to the basement can be estimated reliably. In addition, the basement speed and geometric parameters are estimated consistently, and all parameters compare favorably with the geophysical ground-truth information and with previous inversion results.     

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

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

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

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

浅海中沉积层参数对多途到达时差反演声速剖面的影响

本文以简正波模型和声线模型为基础，通过计算机仿真，研究了浅海多途到达时延差结构与沉积层密度、声速、传播损失系数等参数的关系，针对100m左右浅海和600-1300Hz的信号频率，对两种模型仿真结果进行了比较，结果表明，声线模型仿真得到的多途到达时延差结构与简正波模型仿真的结果基本一致，沉积层参数失配对基于声线模型的多途到达时差反演SSP的影响为：反演声速剖面的均方根误差的变化小于1.45m/s，平均误差的变化小于2.60m/s，在所使用的频段内，浅海沉积层参数对多途到达时差反演浅海声速剖面的影响不大，工程应用中可以忽略不计。     

气泡体积分数对沙质沉积物低频声学特性的影响

由于有机物质分解等原因,实际的海底沉积物中存在气泡,气泡的存在会显著影响沉积物低频段的声学特性,因此研究气泡对沉积物低频段声速的影响机理具有重要意义.考虑到外场环境的不可控性,在室内水池中搭建了大尺度含气非饱和沙质沉积物声学特性获取平台,在有界空间中应用多水听器反演方法首次获取了含气非饱和沙质沉积物300—3000 Hz频段内的声速数据(79—142 m/s),并同时利用双水听器法获取了同一频段的数据(112—121 m/s).在声波频率远低于沉积物中最大气泡的共振频率时,根据等效介质理论,将孔隙水和气泡等效为一种均匀流体,改进了水饱和等效密度流体近似模型.模型揭示了气泡对沉积物低频段声学特性的影响规律,理论上解释了沉积物中声速的降低.通过分析模型预报声速对模型参数的敏感性,根据测量得到的声速分布反演得到了沉积物不同区域的气泡体积分数,气泡体积分数从1.07%变化到2.81%.改进的模型为沉积物中气泡体积分数估计提供了一种新方法.     

Acoustic inversions for measuring boundary layer suspended sediment processes

Although sound has been applied to the study of sediment transport processes for a number of years, it is acknowledged that there are still problems in using the backscattered signal to measure suspended sediment parameters. In particular, when the attenuation due to the suspension becomes significant, the uncertainty associated with the variability in the scattering characteristics of the sediments in suspension can lead to inversion errors which accumulate as the sound propagates through the suspension. To study this attenuation propagation problem, numerical simulations and laboratory experiments have been used to assess the impact unpredictability in the scattering properties of the suspension has on the acoustically derived suspended sediments parameters. The results clearly show the commonly applied iterative implicit inversion can lead to calculated sediment parameters, which become increasingly erroneous with range, as the sound propagates through the suspension. To address this problem an alternative approach to the iterative implicit formulation is investigated using a recently described dual frequency inversion. This approach is not subject to the accumulation of errors and has an explicit solution. Here the dual frequency inversion is assessed and calculated suspended sediment parameters are compared with those obtained from the iterative implicit inversion.     

利用海洋环境噪声估计浅海海底表层声速

为提高反演效率,提出一种快速估计浅海海底表层声速的方法。根据噪声能流理论,垂直阵接收的海洋环境噪声数据能够用于无源提取海底反射损失,反射损失曲线中具有明显的临界角效应,从而接估计海底表层声速。以射线模型为基础,推导了噪声提取的反射损失与理论值之间的差异,并讨论实际阵列波束形成在不同角度和频率下的性能。考虑到非等声速环境下声线会发生弯曲,需对角度进行修正以提高方法的广泛适用性。不同频率的临界角与有效深度之间存在对应关系,黄海某海区数据处理结果表明,在临界角不变的有效深度内,海底表层可以视为等声速层,该海区海底表层0.5 m内声速估计结果为1547 m/s,与有源反演结果相近。      

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.     

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

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