利用拖船自噪声进行浅海环境参数贝叶斯反演

研究了以拖船自噪声为参考声源的浅海环境参数反演问题,并针对反演结果不确定性快速量化评估问题,提出了一种基于自适应重要性抽样的贝叶斯反演新方法。反演利用了拖船自噪声低频线谱成分,并采用混合高斯推荐函数自适应推荐声场模型样本,使得样本集中于参数高概率密度区域,实现后验概率密度快速收敛计算。仿真试验结果表明:拖船自噪声反演能够准确估计水深、沉积层及阵列参数等。所提自适应重要性抽样贝叶斯反演方法的计算效率优于快速吉布斯抽样方法。利用试验数据处理验证,反演得到试验海域声学环境参数,计算传播损失与各阵元接收线谱强度变化吻合,说明反演最优环境模型能准确表征声场传播特征。      

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

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

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

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

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

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

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

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

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

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

利用浅海声场干涉条纹反演海底参数方法研究

浅海频散特性导致声场在距离-频率平面形成干涉条纹,该干涉条纹蕴含着声源以及海底信息。本文通过降噪、滤波等方法,提取船舶辐射噪声干涉条纹中距离域干涉准周期、波导不变量以及去扩展衰减系数(不考虑扩展损失后接收到噪声级随距离的衰减程度)三种特征量的观测值,基于无限均匀沉积层海底模型,采用简正波计算程序理论计算上述三种特征量作为拷贝,提出一种以模拟退火法作为搜索寻优策略的新的匹配反演海底参数的方法。海试数据处理结果表明,该方法反演估计所得海底密度、声速、衰减系数三参数用于浅海负跃层条件下的声场预报效果比较良好。     

匹配场处理舰船辐射噪声级估计方法

针对水声信道对舰船辐射噪声声传播的影响,进而导致声源级测量结果不准确的问题,提出了基于匹配场处理的舰船辐射噪声级估计方法。在海洋环境噪声为空间均匀高斯白噪声的假设下,当海洋环境参数已知、信噪比满足一定要求时,匹配场处理能有效地给出被测噪声源的位置信息及该位置处的能量响应。从能量估计角度出发,推导了声源位置处匹配场输出响应的能量修正因子计算公式,从理论上证明了匹配场处理在被测声源位置处输出响应与能量修正因子的乘积为真实声源级的最小方差无偏(MVU)估计。该方法首先选择合适的声场计算模型计算拷贝场向量,对接收到的辐射噪声信号进行匹配场处理,得出接收信号级和被测声源位置;其次利用该位置所对应的拷贝场向量替换能量修正因子公式中的真实信道传输函数以计算能量修正因子的估计值;最后由接收信号级与能量修正因子估计值相乘得出舰船辐射噪声声源级的MVU估计。针对典型的浅海水声信道,进行了计算机仿真试验,结果表明:该方法能有效地进行舰船辐射噪声测量,当信噪比满足一定要求时,测量得到的声源级与实际声源级相比,误差小于1 dB。      

基于空时级联滤波的拖船自噪声抵消方法

拖船自噪声属于近场多途强干扰信号,具有显著角度扩展现象,是影响拖曳阵探测性能的重要因素之一。在逆波束形成自噪声抵消方法基础上,结合空域和时域滤波,提出了空时级联滤波的拖船自噪声抵消方法。首先,由方位估计器估计拖船自噪声方位,经过波束形成空域滤波估计得到自噪声信号;然后,利用自噪声空域滤波输出设计最优维纳滤波器,将自噪声抵消在阵元级。基于简正波声场模型的仿真和实际海试数据的处理表明,该方法能够有效抵消拖船自噪声,性能优于逆波束形成方法。     

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

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

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.     

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.     

Bottom parameter behavior in shallow water

This paper discusses the matched field processing behavior of key geoacoustic parameters as a function of frequency (25-100 Hz), range (250 m-2 km), and sediment layer thickness (12, 22, 40 m with scenarios corresponding to possible Shallow Water'06 tests) in a systematic manner for simulated data. It considers ideal, as well as approximate knowledge of other conditions, such as source depth, phone depth, and an ocean sound-speed profile. Some of these approximations cause more trouble than others, but they are all still manageable at the low frequencies and short ranges of interest here. A single layer sediment is always assumed here, as is a vertical array and range independence. However, the approach is easily adaptable to other situations. Finally, we introduce a new way of processing broadband information and find that averaging over all available frequencies and ranges may not be optimal for maximizing sensitivity to, e.g., finding, geoacoustic parameters.     

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

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

Statistical geoacoustic inversion from vertical correlation of shallow water reverberation

Statistical geoacoustic inversion results based on reverberation vertical correlation （RVC） data from 300 Hz to 800 Hz are presented. The data were obtained during Yellow Sea Experiment-2005 （YSE-05）. An uncertainty analysis is carried out. It is found that the inversion sea bottom sound velocities decrease when the frequency increases. So it is difficult to determine the sea bottom sound velocity. In order to solve this problem, a two-layer bottom model is assumed, and a multi-frequency inversion approach based on Genetic Algorithm is proposed. The approach is demonstrated using YSE-05 experiment data. Both RVC and normal modes depth functions （NMDFs） calculated using the inverted geoacoustic parameters are in good agreement with the experimental measurements.     

Range-dependent seabed characterization by inversion of acoustic data from a towed receiver array

The MAPEX2000 experiments were conducted in the Mediterranean Sea in March, 2000 to determine seabed properties using a towed acoustic source and receiver array. Towed systems are advantageous because they are easy to deploy from a ship and the moving platform offers the possibility for estimating spatially variable (range-dependent) seabed properties. In this paper, seabed parameters are determined using a matched-field geoacoustic inversion approach with measured, towed array data. Previous research has successfully applied matched-field geoacoustic inversion techniques to measured acoustic data. However, in nearly all cases the inverted data were collected on moored, vertical receiver arrays. Results here show that seabed parameters can also be extracted by inverting acoustic measurements from a towed array of receivers, and these agree with those inverted using data received simultaneously on a vertical array. These findings imply that a practical technique could be developed to map range-dependent seabed parameters over large areas using a towed acoustic system. An example of such a range-dependent inversion is given using measurements from the MAPEX2000 experiments.