浅海声速剖面与移动声源的跟踪定位

在水平非均匀分布的浅海环境中,针对移动声源跟踪时,声速剖面的变化会对声场产生影响,提出了一种利用集合卡尔曼滤波算法的声速剖面跟踪反演和移动声源跟踪定位的方法。首先,将声速剖面进行距离和深度的参数化表示,从而将对声速剖面的跟踪转化为对声速剖面前3阶经验正交函数系数的跟踪;其次,通过将声源状态信息和声速剖面信息表示为状态变量,而将垂直线列阵接收到的声场信息作为测量值建立状态-测量模型,然后利用集合卡尔曼滤波方法对模型状态变量进行跟踪。仿真结果得出:声速剖面跟踪反演的均方根误差和移动声源跟踪定位的绝对误差都非常小,对声源的跟踪定位精度很高。并且通过增加集合样本数、增加接收信号信噪比以及增加接收阵元数目都可以提高跟踪定位结果精度。最后,利用东海实验数据对本方法进行了验证。      

自回归状态空间模型下时变声速剖面跟踪方法

讨论了一种适用于浅海的时变声速剖面跟踪方法。该方法将时变水体声速剖面的反演问题建模为由描述声速剖面时变特性的状态方程与包含声压场局部测量信息的测量方程组成的状态-空间模型,提出利用自回归分析拟合方法将声速场扰动建模为高阶自回归演化模型,并通过集合卡尔曼滤波序贯地估计时间演化的海洋声速场。利用2001年亚洲海实验环境与声速测量数据,仿真分析了基于高阶自回归演化模型的时变声速剖面集合卡尔曼滤波估计方法。结果表明,相比于利用传统随机游走状态演化模型的估计方法,该改进方法可有效降低声速的跟踪误差,并且在较低信噪比条件下仍具有较好的跟踪性能。      

采用小生境遗传算法反演浅海声速剖面研究*

针对浅海声速剖面反演问题,采用小生境遗传算法,结合声线搜索的最快本征声线匹配反演,实现浅海负跃层条件下的声速剖面估计。利用经验正交函数对声速剖面的多参数不确定性降维,依据声场计算模型获取的最快特征声线传播时延与观测声传播时长进行匹配,采用小生境遗传优化算法,获取最优经验正交函数估计,实现声速剖面反演。按上述方法反演处理浅海声传播实验数据,结果表明,该方法针能够有效反演浅海声速剖面,并且显著优于传统遗传算法反演结果。     

非等声速信道下的循环神经网络机动目标跟踪模型

针对水下非等声速信道声线弯曲导致传统滤波跟踪轨迹偏差的问题,提出一种循环神经网络的目标跟踪模型。该模型在缺乏声速剖面信息的情况下,通过数据驱动迭代训练,学习输入观测值与输出状态值之间的映射关系,实现目标位置和瞬态特征变化的精确获知。蒙特卡洛仿真实验结果表明,本文模型在非等声速信道下复杂机动场景中相较传统单模型滤波算法以及交互式多模型算法,水平距离跟踪精度分别提升4.06%,1.57%,深度估计精度分别提升0.87%,0.85%。本文模型相较于传统滤波方法具有更高的跟踪精度,并且能够在失配声速分布信道下进行迁移学习,提升模型在失配声场环境下的泛化性。      

深海上层声速剖面水平非均匀分布对声场会聚区特性预报的影响

针对深海上层声速剖面水平方向弱变化对声场会聚区特征的影响,通过对0°出射声线传播轨迹的分析,提出了只考虑声源及会聚区位置附近声速剖面水平变化实现深海声场会聚区特征准确预报的方法,该方法有助于降低深海声传播调查时声速剖面水平非均匀分布的测量工作量。通过模拟小振幅内波引起深海上层声速剖面水平非均匀变化,对该方法进行了仿真验证。南海声传播试验数据分析表明,当不考虑上层声速剖面的水平非均匀分布时,模型预报的声场会聚区位置及会聚区能量分布与实测结果有一定差异,该差异具有距离累积效应,到第二会聚区时位置累计误差可达3 km,当考虑会聚区宽度内声速剖面水平变化时,模型预报的声场会聚区位置与实测结果非常吻合,第二会聚区位置预报误差小于1 km。     

浅海复杂环境下等效声速剖面的构建方法*

针对浅海复杂环境声速剖面水平变化情况下的声传播损失预报及目标定位问题,提出了一种基于遗传算法的等效声速剖面重构算法。首先,将声速剖面进行时间和空间上的分解,从而将声速剖面抽象为对声速剖面前三阶正交函数系数的反演;其次,利用遗传算法,以先验声速剖面集为基础,进行参数反演。仿真结果表明,在浅海复杂条件下,传播损失预报受声速剖面及海洋参数的影响,不能直接运用接收或发射位置处的声速剖面进行传播损失预报,否则会对预报结果造成误差。通过构建具有声传播累积效应的等效声速剖面可以提高匹配场定位精度,完成目标定位,且在构建等效声速剖面时,接收位置处即本地声速剖面所占权重较大。     

一种声速剖面展开的正交基函数获取方法

为了评估经验正交函数(Empirical Orthogonal Functions)在声速剖面的声学反演应用中的可靠性和适用性,降低声速剖面基函数的获取对样本数据的依赖性,提出了一种声速剖面展开的正交基函数的获取方法.该方法通过分析引起声速剖面变化的海洋动力学过程,利用动力学方程求解得到的水动力简正波(Hydrodynamic Normal Modes)来给出声速剖面展开的基函数.HNMs方法获得的声速剖面正交基函数与EOFs方法的结果非常吻合,且对样本数据的依赖性小。HNMs方法物理意义清晰,能够给出决定声速剖面正交基函数的物理参数,使评估经验正交函数的准确性和适用性成为可能。      

基于强跟踪平方根容积卡尔曼滤波的纯方位目标运动分析方法

针对纯方位目标跟踪系统中模型状态简化、系统噪声统计特性未知、目标初始距离信息不准确导致的滤波收敛时间长和滤波精度不高的问题,以自主水下机器人(Autonomous Underwater Vehicle, AUV)跟踪水下动态目标为例,提出了一种基于强跟踪平方根容积卡尔曼滤波器(Strong Tracking Square Root Cubature Kalman Filter, STFSRCKF)的纯方位目标运动分析算法。该算法在滤波过程中,利用平方根容积卡尔曼滤波器(Square Root Cubature Kalman Filter, SRCKF)完成预测更新,对于SRCKF中的每个容积点采用强跟踪滤波器(Strong Tracking Filter, STF)进行更新,设计滤波增益以抑制噪声对系统状态估计的影响,有效提高了滤波的数值稳定性,减小了状态估计误差。通过仿真分析,比较了扩展卡尔曼滤波器(Extended Kalman Filter, EKF)、无迹卡尔曼滤波器(Unscented Kalman Filter, UKF)、平方根容积卡尔曼滤波器(Square-Root Cubature Kalman Filter, SRCKF)、STFSRCKF的算法性能,实验表明所提算法具有跟踪速度快,精度高等优点。     

混合退火粒子滤波器

针对非线性、非高斯系统状态的在线估计问题，提出一种新的基于序贯重要性抽样的粒子滤波算法. 在滤波算法中，用状态参数分解和退火系数来产生重要性概率密度函数，此概率密度函数综合考虑了转移先验、似然、噪声的统计特性以及最新的观察数据，因此更接近于系统状态的后验概率. 理论分析与仿真实验表明该粒子滤波器的性能明显优于标准的粒子滤波器和扩展卡尔曼滤波器. 关键词： 非线性 非高斯 粒子滤波 序贯重要性抽样     

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

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

集合卡尔曼滤波在时变声速剖面追踪中的性能分析

对集合卡尔曼滤波在时变海洋环境下的声速剖面追踪性能进行了分析。将南海实验背景下普林斯顿海洋模型预报的声速剖面正交分解为3阶系数组成的状态-空间形式,其状态转移方程建模为3阶自回归过程;基于卡尔曼反馈理论,利用适合于水平非均匀模型RAM仿真的观测声压场对系统状态进行校正,实现声速剖面的动态追踪。在水平均匀、水平非均匀和海底参数失配环境下的仿真结果均能较好地实现对声速剖面的追踪,验证了算法的可行性。同时对不同信噪比、粒子数、阵元数和海底参数失配等情况下的分析表明,观测信息量的增加可以有效抑制观测误差和模型误差的影响,相关结论得到了实验数据的验证。      

Tracking of time-evolving sound speed profiles with an auto-regressive state-space model

An approach for time-evolving sound speed profiles tracking in shallow water is discussed.The inversion of time-evolving sound speed profiles is modeled as a state-space estimation problem,which includes a state equation for predicting the time-evolving sound speed profile and a measurement equation for incorporating local acoustic measurements.In the paper,auto-regression(AR) method is introduced to obtain a high-order AR evolution model of the sound speed field time variations,and the ensemble Kalman filter is utilized to track the sound speed field.To validate the approach,the accuracy in sound speed estimation is analyzed via a numerical implementation using the ASIAEX experimental environment and the sound velocity measurement data.Compared with traditional approaches based on the state evolution represented as a random walk,simulation results show the proposed AR method can effectively reduce the tracking errors of sound speed,and still keep good tracking performance at low signal-to-noise ratios.     

浅海时变声速环境下的自适应匹配场定位算法实现*

针对浅海环境下声速剖面失配引起的匹配场处理器失配问题,提出了一种自适应匹配场定位算法在声速剖面时变环境下的实现方式。将先验声速剖面集简化为经验正交函数表示,结合蒙特卡洛方法与环境扰动约束算法对当下时刻的目标声源进行匹配场定位。本文以某次试验获取的连续20小时的声速剖面数据为研究对象,通过仿真试验对该算法进行验证,结果表明：在先验声速剖面集的半小时之后,利用自适应算法的距离和深度定位成功率较常规匹配场算法有较大提升,其中,深度正确定位概率相对较低。     

Effect of shallow water internal waves on ocean acoustic striation patterns

Contour plots of underwater acoustic intensity, mapped in range and frequency, often exhibit striations. It has been claimed that a scalar parameter 'beta', defined in terms of the slope of the striations, is invariant to the details of the acoustic waveguide. In shallow water, the canonical value is β=1. In the present paper, the waveguide invariant is modelled as a distribution rather than a scalar. The effects of shallow water internal waves on the distribution are studied by numerical simulation. Realizations of time-evolving shallow water internal wave fields are synthesized and acoustic propagation simulated using the parabolic equation method. The waveguide invariant distribution is tracked as the internal wave field evolves in time. Both random background internal waves and more event-like solitary internal waves are considered.     

Effect of shallow water internal waves on ocean acoustic striation patterns

Abstract

Contour plots of underwater acoustic intensity, mapped in range and frequency, often exhibit striations. It has been claimed that a scalar parameter ‘beta’, defined in terms of the slope of the striations, is invariant to the details of the acoustic waveguide. In shallow water, the canonical value is β=1. In the present paper, the waveguide invariant is modelled as a distribution rather than a scalar. The effects of shallow water internal waves on the distribution are studied by numerical simulation. Realizations of time-evolving shallow water internal wave fields are synthesized and acoustic propagation simulated using the parabolic equation method. The waveguide invariant distribution is tracked as the internal wave field evolves in time. Both random background internal waves and more event-like solitary internal waves are considered.     

Tracking-positioning of sound speed profiles and moving acoustic source in shallow water

An ensemble Kalman filter(EnKF) approach is proposed to perform sequential tracking of water column sound speed profile(SSP) using a moving acoustic source. First,the SSPs are discretized in depth and range, and are expressed by the empirical orthogonal functions(EOFs). Second, the acoustic source state information and the first three orders of EOF coefficients are expressed as the state variable, and the acoustic field information received by the vertical line array are the measured values. Successively, the state variables and measured values are used to establish the state-measure model. Last, the EnKF is utilized to track the state variables. The simulation results show that the root mean square error of SSP and the absolute error of source are all small, and thus the acoustic source tracking-positioning has high accuracy. Moreover, increasing the number of sample collection, the signal-to-noise ratio and the number of receiving elements can improve the tracking-positioning results. The method is verified using the experimental data of the East China Sea.     

A real-time plane-wave decomposition algorithm for characterizing perforated liners damping at multiple mode frequencies

Perforated liners with a narrow frequency range are widely used as acoustic dampers to stabilize combustion systems. When the frequency of unstable modes present in the combustion system is within the effective frequency range, the liners can efficiently dissipate acoustic waves. The fraction of the incident waves being absorbed (known as power absorption coefficient) is generally used to characterize the liners damping. To estimate it, plane waves either side of the liners need to be decomposed and characterized. For this, a real-time algorithm is developed. Emphasis is being placed on its ability to online decompose plane waves at multiple mode frequencies. The performance of the algorithm is evaluated first in a numerical model with two unstable modes. It is then experimentally implemented in an acoustically driven pipe system with a lined section attached. The acoustic damping of perforated liners is continuously characterized in real-time. Comparison is then made between the results from the algorithm and those from the short-time fast Fourier transform (FFT)-based techniques, which are typically used in industry. It was found that the real-time algorithm allows faster tracking of the liners damping, even when the forcing frequency was suddenly changed.     

Acoustic propagation through anisotropic internal wave fields: transmission loss, cross-range coherence, and horizontal refraction

Results of a computer simulation study are presented for acoustic propagation in a shallow water, anisotropic ocean environment. The water column is characterized by random volume fluctuations in the sound speed field that are induced by internal gravity waves, and this variability is superimposed on a dominant summer thermocline. Both the internal wave field and resulting sound speed perturbations are represented in three-dimensional (3D) space and evolve in time. The isopycnal displacements consist of two components: a spatially diffuse, horizontally isotropic component and a spatially localized contribution from an undular bore (i.e., a solitary wave packet or solibore) that exhibits horizontal (azimuthal) anisotropy. An acoustic field is propagated through this waveguide using a 3D parabolic equation code based on differential operators representing wide-angle coverage in elevation and narrow-angle coverage in azimuth. Transmission loss is evaluated both for fixed time snapshots of the environment and as a function of time over an ordered set of snapshots which represent the time-evolving sound speed distribution. Horizontal acoustic coherence, also known as transverse or cross-range coherence, is estimated for horizontally separated points in the direction normal to the source-receiver orientation. Both transmission loss and spatial coherence are computed at acoustic frequencies 200 and 400 Hz for ranges extending to 10 km, a cross-range of 1 km, and a water depth of 68 m. Azimuthal filtering of the propagated field occurs for this environment, with the strongest variations appearing when propagation is parallel to the solitary wave depressions of the thermocline. A large anisotropic degradation in horizontal coherence occurs under the same conditions. Horizontal refraction of the acoustic wave front is responsible for the degradation, as demonstrated by an energy gradient analysis of in-plane and out-of-plane energy transfer. The solitary wave packet is interpreted as a nonstationary oceanographic waveguide within the water column, preferentially funneling acoustic energy between the thermocline depressions.