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

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

浅海小掠射角的海底界面声反向散射模型的简化

在浅海,尤其是负梯度声速剖面和海面较为平静的浅海波导,海底界面反向散射是浅海混响的主要来源.经验散射模型只适用于分析浅海混响平均强度衰减特性,而基于物理机理建立的反向散射模型克服了这一缺陷,但同时也引入了其受地声模型约束的问题.本文结合了海底反射系数的三参数模型,对浅海远场海底反向散射模型进行了简化,以减少地声模型的输入参数.理论分析了海底反射系数的相移参数可以描述海底对声场的散射作用,无需任何海底地声参数的先验知识.通过对海底反向散射模型近似简化,结果表明在临界角附近和甚小掠射角范围内的海底粗糙界面反向散射模型的角度特性和强度特性受海底沉积层的影响不同:在临界角附近,海底反向散射的角度特性受海底反射系数的相移参数加权,而其散射系数则近似与相移参数无关;对于甚小掠射角,海底反向散射的角度特性近似与海底反射系数的相移参数无关,其散射系数则近似与相移参数的4次方成正比.     

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

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

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

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

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

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

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

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

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

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

浅海环境下数据同化声层析方法研究

讨论了适用于浅海的基于声速局部测量模型、声传播模型及声压场局部测量模型的声学数据同化方法,并给出了具体的执行算法。该算法根据最优化准则把局部声速和声压场的先验测量信息,以及声传播模型进行了有效地融合,提高了海洋环境参数估计的精度,为浅海声层析提供了新思路。同时,利用实验浅海声速测量数据,通过经验正交函数对实际海洋声速剖面进行了估计,并分析了各类噪声对水体声速场及海底声学参数估计精度的影响,验证了该执行算法的有效性。      

低频混响信号中的简正波衰减系数提取方法

提出一种从低频混响信号中提取简正波衰减系数的方法。利用简正波过滤技术对垂直阵混响信号进行分析,获得单阶简正波混响声场。假设海底反向散射矩阵可分离,从单阶简正波平均混响强度中提取出有效海底反向散射矩阵元素,最后利用不同距离上的有效海底反向散射矩阵元素计算出简正波的衰减系数。利用该方法从混响信号中提取出的简正波衰减系数预报的声传播损失和相同海域实测声传播损失一致。该简正波衰减系数提取方法有效避免了海底散射衰减和简正波传播衰减耦合的问题,同时对海底参数反演和水声环境的快速评估也具有重要的意义。      

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

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

Geoacoustic Inversion Based on a Vector Hydrophone Array

We propose a geoacoustic inversion scheme employing a vector hydrophone array based on the fact that vector hydrophone can provide more acoustic field information than traditional pressure hydrophones. Firstly, the transmission loss of particle velocities is discussed. Secondly, the sediment sound speed is acquired by a matchedfield processing （MFP） procedure, which is the optimization in combination of the pressure field and vertical particle velocity field. Finally, the bottom attenuation is estimated from the transmission loss difference between the vertical particle velocity and the pressure. The inversion method based on the vector hydrophone array mainly has two advantages： One is that the MFP method based on vector field can decrease the uncertain estimation of the sediment sound speed. The other is that the objective function based on the transmission loss difference has good sensitivity to the sediment attenuation and the inverted sediment attenuation is independent of source level. The validity of the inverted parameters is examined by comparison of the numerical results with the experimental data.     

Sequential Parameter Estimation Using Modal Dispersion Curves in Shallow Water

Existing sequential parameter estimation methods use the acoustic pressure of a line array as observations. The modal dispersion curves are employed to estimate the sound speed profile(SSP) and geoacoustic parameters based on the ensemble Kalman filter. The warping transform is implemented to the signals received by a single hydrophone to obtain the dispersion curves. The experimental data are collected at a range-independent shallow water site in the South China Sea. The results indicate that the SSPs are well estimated and the geoacoustic parameters are also well determined. Comparisons of the observed and estimated modal dispersion curves show good agreement.     

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.     

Geoacoustic inversion with ships as sources

Estimation of geoacoustic parameters using acoustic data from a surface ship was performed for a shallow water region in the Gulf of Mexico. The data were recorded from hydrophones in a bottom mounted, horizontal line array (HLA). The techniques developed to produce the geoacoustic inversion are described, and an efficient method for geoacoustic inversion with broadband beam cross-spectral data is demonstrated. The performance of cost functions that involve coherent or incoherent sums over frequency and one or multiple time segments is discussed. Successful inversions for the first sediment layer sound speed and thickness and some of the parameters for the deeper layers were obtained with the surface ship at nominal ranges of 20, 30, or 50 water depths. The data for these inversions were beam cross-spectra from four subapertures of the HLA spanning a little more than two water depths. The subaperture beams included ten frequencies equally spaced in the 120-200 Hz band. The values of the geoacoustic parameters from the inversions are validated by comparisons with geophysical observations and with the parameter values from previous inversions by other invesigators, and by comparing transmission loss (TL) measured in the experiment with modeled TL based on the inverted geoacoustic parameters.     

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.     

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.     

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.     

Model-based seafloor characterization employing multi-beam angular backscatter data--a comparative study with dual-frequency single beam

Sediment geoacoustic inversion results are estimated employing a multi-beam (MB) echo-sounding system operable at 95 kHz. To characterize the western continental shelf of India (off Goa) seafloor, MB backscatter signals were acquired along with grab sediment samples. The substrate type and roughness of the site were estimated using the composite roughness scattering model with the measured backscatter values. The seafloor parameters, namely mean grain size (M(φ)); roughness spectrum strength (w(2)) and exponent (γ(2)); and sediment volume parameter (σ(2)), for coarse and fine grain sediments are estimated by employing the MB system. These parameters have also been estimated at two other frequencies (33 and 210 kHz) and are compared to the ground truth data to provide sufficient support in validating the model results and increasing the understanding of the shelf seafloor processes. Distinct interclass separations between the sediment provinces are evident from the estimated mean grain size M(φ) and water-sediment interface roughness w(2). The seafloor parameters for coarse and fine grain sediments derived from the 95 kHz MB data are consistent with the sediment sample data as well as with the inversion results obtained using backscatter data at 33 and 210 kHz from the same locations.