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

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

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

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

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.     

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

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

海底单参数频域反演与沉积物分类

为了快速获取海底声学特性,从降低反演参数维数入手,提出一种基于海底单参数(小掠射角下海底反射损失随掠射角的变化率)的频域反演方法,并用于沉积物分类。在海底单参数模型基础上,对海底反射信号进行频域相关处理,利用几何射线理论,推导获得海底单参数与归一化频域自相关系数下降至0.5时频移量之间的关系表达式。通过研究海底单参数与沉积物类型之间的等效关系,对海底沉积物进行分类。对波兰波罗的海实验数据进行实际反演,得到1000~2500 Hz频段内的海底单参数值,进而判定实验海域海底沉积物类型为极细砂。结果表明,根据单参数反演值预报的声传播损失与实验测量值符合较好,利用频域反演结果获得的海底沉积物类型与取样结果一致,验证了实验海域频域反演方法的有效性及分类结果的正确性。      

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.     

warping变换提取单模态反演海底衰减系数

为了获得浅海海底地声模型参数,利用warping变换方法分离出单模态简正波.对于接收深度固定、定深爆炸声源情况,以简正波理论为基础定义了距离归一化的简正波传播损失,并且其随传播的距离呈线性关系,故可通过此变化规律得到声压值实部的衰减因子,进而可求得海底地声模型参数:海底衰减系数.为验证此方法的有效性,仿真了warping变换提取单模态简正波的过程,同时将warping变换提取的单模态简正波与数值计算的结果进行比较验证;并针对某次黄海试验数据进行了处理,得到在150—550 Hz频带范围内海底衰减随频率的变化规律为α=0.581f_k~(1.86)(dB/m).通过与其他学者在相同海域试验结果的对比验证,变化规律基本相同.此外不同模态间反演相同频点的衰减系数接近也较好地支撑了结果.     

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

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

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.     

Propagation in a waveguide with range-dependent seabed properties

The ocean environment contains features affecting acoustic propagation that vary on a wide range of time and space scales. A significant body of work over recent decades has aimed at understanding the effects of water column spatial and temporal variability on acoustic propagation. Much less is understood about the impact of spatial variability of seabed properties on propagation, which is the focus of this study. Here, a simple, intuitive expression for propagation with range-dependent boundary properties and uniform water depth is derived. It is shown that incoherent range-dependent propagation depends upon the geometric mean of the seabed plane-wave reflection coefficient and the arithmetic mean of the cycle distance. Thus, only the spatial probability distributions (pdfs) of the sediment properties are required. Also, it is shown that the propagation over a range-dependent seabed tends to be controlled by the lossiest, not the hardest, sediments. Thus, range-dependence generally leads to higher propagation loss than would be expected, due for example to lossy sediment patches and/or nulls in the reflection coefficient. In a few instances, propagation over a range-dependent seabed can be calculated using range-independent sediment properties. The theory may be useful for other (non-oceanic) waveguides.     

Geoacoustic inversion for fine-grained sediments

Fine-grained water saturated sediments like silty clays have the curious property that the speed of sound through its bulk medium is lower than that of the interstitial pore fluid. When a fine-grained sediment is at the water-sediment interface, classical theory predicts that there is an angle at which the reflection coefficient is zero, and there is total transmission of sound into the seabed. This angle is called the angle of intromission and has been directly observed at the seafloor only rarely. Data from a new measurement technique show this phenomenon with remarkable clarity. The presence of the angle of intromission creates an opportunity for a direct (i.e., without search) inversion for the sediment sound speed and density. Though acoustic techniques generally do not estimate sediment density very precisely, this technique is quite sensitive to density. Geoacoustic inversion results from reflection measurements in the Straits of Sicily compare very favorably with independent "ground truth" data indicating that the method is robust.     

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

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

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.     

基于拖曳倾斜线列阵的海底反射损失提取方法

在浅海环境中,海底声学参数对水下声场的精确预报十分重要.现有的海底声学参数反演方法大多数是采用固定垂直阵进行的,其缺点是不能实施大面积、高分辨的走航式反演.在已提出的垂直阵海底反射损失反演方法的基础上,研究了基于三种线列阵阵形的海底反射损失提取方法及其特点,提出了基于拖曳倾斜线列阵与三个声源组合的走航式海底参数快速获取方法,并对该方法进行了误差分析.研究结果表明：覆盖同样的掠射角范围,垂直线列阵需与多个距离的声源组合,拖曳水平线列阵只需一个声源组合但需要采用较大的物理孔径,拖曳倾斜线列阵综合了垂直阵和水平 关键词： 水下声场预报 海底声学参数 海底反射损失 拖曳倾斜线列阵     

考虑噪声源深度分布的海洋环境噪声模型及地声参数反演

考虑到海洋环境噪声源深度分布不集中,建立了噪声源随深度分布的海洋环境噪声模型,分析了源深度对噪声场垂向特征的影响并从简正波角度予以解释,发现海底声阻抗和声源深度都显著影响由海洋环境噪声获得的等效海底反射损失大掠射角部分,进而将该模型用于地声参数反演.两段实测噪声数据200—525 Hz频段的反演结果表明:基于海洋环境噪声的地声参数反演最优值与声传播的反演结果相近;源平均深度最优值随频率增加有变小的趋势,说明随频率增加环境噪声主要贡献源逐渐由航船转为风浪;当海况大于3级时,400 Hz以上频段噪声源深度平均值很小,与Monahan气泡理论的描述一致.     

Minimization of finite beam effects in the determination of reflection and transmission coefficients of an elastic layer

A special data acquisition technique was applied to determine the acoustic plane-wave reflection and transmission properties of a plane-parallel aluminum plate. In this technique, the reflected and transmitted wavefield along a plane or line normal to the reflected or transmitted wave vector is recorded at equidistant receiver positions. The obtained traces are subsequently added up in the temporal domain to satisfy plane-wave conditions, thus effectively removing the effect of the limited beam of commonly used transducers. The agreement between plane-wave theory and experiment was found to be excellent, both in the temporal and in the frequency domain.     

Seabed geoacoustic characterization with a vector sensor array

This paper proposes a vector sensor measurement model and the related Bartlett estimator based on particle velocity measurements for generic parameter estimation, illustrating the advantages of the Vector Sensor Array (VSA). A reliable estimate of the seabed properties such as sediment compressional speed, density and compressional attenuation based on matched-field inversion (MFI) techniques can be achieved using a small aperture VSA. It is shown that VSAs improve the resolution of seabed parameter estimation when compared with pressure sensor arrays with the same number of sensors. The data considered herein was acquired by a four-element VSA in the 8-14 kHz band, during the Makai Experiment in 2005. The results obtained with the MFI technique are compared with those obtained with a method proposed by C. Harrison, which determines the bottom reflection loss as the ratio between the upward and downward beam responses. The results show a good agreement and are in line with the historical information for the area. The particle velocity information provided by the VSA increases significantly the resolution of seabed parameter estimation and in some cases reliable results are obtained using only the vertical component of the particle velocity.     

Near-specular acoustic scattering from a buried submarine mud volcano

Submarine mud volcanoes are objects that form on the seafloor due to the emission of gas and fluidized sediment from the Earth's interior. They vary widely in size, can be exposed or buried, and are of interest to the underwater acoustics community as potential sources of active sonar clutter. Coincident seismic reflection data and low frequency bistatic scattering data were gathered from one such buried mud volcano located in the Straits of Sicily. The bistatic data were generated using a pulsed piston source and a 64-element horizontal array, both towed over the top of the volcano. The purpose of this work was to appropriately model low frequency scattering from the volcano using the bistatic returns, seismic bathymetry, and knowledge of the general geoacoustic properties of the area's seabed to guide understanding and model development. Ray theory, with some approximations, was used to model acoustic propagation through overlying layers. Due to the volcano's size, scattering was modeled using geometric acoustics and a simple representation of volcano shape. Modeled bistatic data compared relatively well with experimental data, although some features remain unexplained. Results of an inversion for the volcano's reflection coefficient indicate that it may be acoustically softer than expected.     

An iterative implementation of rotated coordinates for inverse problems

A generalized inversion method is presented that uses a rotated coordinates technique [Collins and Fishman, J. Acoust. Soc. Am. 98, 1637-1644 (1995)] in simulated annealing to invert for both the location of an acoustic source and parameters that describe the ocean seabed. The rotated coordinates technique not only aids in the inversion process but also indicates the coupling of the source and environmental parameters and the relative sensitivities of the cost function to changes in the various parameters. The information obtained from the rotated coordinates provides insights into how the inversion problem can be effectively decoupled. An iterative process consisting of multiple simulated annealing runs that each use a different set of rotated coordinates is demonstrated. This multistep algorithm is called systematic decoupling using rotated coordinates and is especially helpful when inverting for a large number of unknown parameters. The cost function minimized in the inversion algorithm is model-data cross-hydrophone spectra summed coherently over frequency and receiver pairs. The results of applying this inversion method to simulated data are presented in this paper.     

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

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