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

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

Uncertainty estimation in simultaneous Bayesian tracking and environmental inversion

This paper develops a Bayesian approach for two related inverse problems: tracking an acoustic source when ocean environmental parameters are unknown, and determining environmental parameters using acoustic data from an unknown (moving) source. The formulation considers source and environmental parameters as unknown random variables constrained by noisy acoustic data and by prior information on parameter values (e.g., physical limits for environmental properties) and on inter-parameter relationships (limits on radial and vertical source speed). The goal is not simply to estimate parameter values, but to rigorously determine parameter uncertainty distributions, thereby quantifying the information content of the data/prior to resolve source and environmental parameters. Results are presented as marginal posterior probability densities (PPDs) for environmental parameters and joint marginal PPDs for source ranges and depths. Given the numerically intensive inversion, an efficient Markov-chain Monte Carlo importance-sampling approach is developed which combines Metropolis and heat-bath Gibbs' sampling, employs efficient proposal distributions based on a linearized PPD approximation, and considers nonunity sampling temperatures to ensure a complete parameter search. The approach is illustrated with two simulated examples representing tracking a quiet submerged source and geoacoustic inversion using noise from an unknown ship of opportunity. In both cases, source, seabed, and water-column parameters are unknown.     

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

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

声场重构的声辐射模态阶数分析

利用声辐射模态法对椭球状声源辐射声场进行重构,分析了模态阶数对相对误差的影响,以及模态阶数与频率的关系;同时深入探讨了声源形状变化对模态阶数选取的影响。研究发现利用较少的模态阶数即可得到很好的重构效果,且声源形状的变化对重构所需模态阶数的影响较小,因此可将该方法用于任意形状声源的重构,体现了此方法在声场重构中的优越性。     

水下声学滑翔机用于东印度洋声传播特性分析及声源估计

常规实验方法无法同步获取深海大尺度声学和水文数据,水下滑翔机可作为同步观测平台解决该问题.首先利用在东印度洋北部海域水下滑翔机同步获取的声传播和水文实验数据,分析了水下滑翔机的自噪声谱级和实验海区声传播特性,然后推算并修正了滑翔机水下运动轨迹,利用第一影区水下滑翔机接收声传播信号的脉冲多途到达时间差对声源进行测距与定深。潜标接收噪声与滑翔机自噪声谱级对比表明,水下滑翔机在海洋中无动力运动时的系统自噪声接近于潜标观测的海洋环境噪声。滑翔机实测的声传播损失与模型计算结果吻合较好,第一影区水下声源测距定深结果与实际位置较为一致,测距与定深的相对误差均小于5%。利用加载水听器的水下滑翔机可以实现水文环境数据与声学信号的同步观测,对深海声传播特性测量及定位算法研究具有重要意义。      

Model-oriented ocean tomography using higher frequency, bottom-mounted hydrophones

A tomographic scheme is presented that ingests ocean acoustic measurements into an ocean model using data from bottom-mounted hydrophones. The short distances between source-receiver pairs (1-10 km) means arrival times at frequencies of 8-11 kHz are readily detectable and often distinguishable. The influence of ocean surface motion causes considerable variability in acoustic travel times. Techniques are presented for measuring travel times and removing the variability due to surface waves. An assimilation technique is investigated that uses differences in measured and modeled acoustic travel times to impose corrections on the oceanographic model. Equations relating travel time differences to oceanographic variables are derived, and techniques are presented for estimating the acoustic and ocean model error covariance matrices. One test case using a single source-receiver pair shows that the tomographic information can have an impact on constraining the solution of the ocean circulation model but can also introduce biases in the predictions. A second test case utilizes knowledge of a bias in a model-predicted variable to limit grid cells that are impacted by the tomographic data. In this case, using the tomographic data results in significant improvements in the model predictions without introducing any biases.     

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.     

A generalized polynomial chaos based ensemble Kalman filter with high accuracy

As one of the most adopted sequential data assimilation methods in many areas, especially those involving complex nonlinear dynamics, the ensemble Kalman filter (EnKF) has been under extensive investigation regarding its properties and efficiency. Compared to other variants of the Kalman filter (KF), EnKF is straightforward to implement, as it employs random ensembles to represent solution states. This, however, introduces sampling errors that affect the accuracy of EnKF in a negative manner. Though sampling errors can be easily reduced by using a large number of samples, in practice this is undesirable as each ensemble member is a solution of the system of state equations and can be time consuming to compute for large-scale problems. In this paper we present an efficient EnKF implementation via generalized polynomial chaos (gPC) expansion. The key ingredients of the proposed approach involve (1) solving the system of stochastic state equations via the gPC methodology to gain efficiency; and (2) sampling the gPC approximation of the stochastic solution with an arbitrarily large number of samples, at virtually no additional computational cost, to drastically reduce the sampling errors. The resulting algorithm thus achieves a high accuracy at reduced computational cost, compared to the classical implementations of EnKF. Numerical examples are provided to verify the convergence property and accuracy improvement of the new algorithm. We also prove that for linear systems with Gaussian noise, the first-order gPC Kalman filter method is equivalent to the exact Kalman filter.     

Inverse radiation problem of temperature distribution in one-dimensional isotropically scattering participating slab with variable refractive index

In this paper, an inverse analysis is performed for estimation of source term distribution from the measured exit radiation intensities at the boundary surfaces in a one-dimensional absorbing, emitting and isotropically scattering medium between two parallel plates with variable refractive index. The variation of refractive index is assumed to be linear. The radiative transfer equation is solved by the constant quadrature discrete ordinate method. The inverse problem is formulated as an optimization problem for minimizing an objective function which is expressed as the sum of square deviations between measured and estimated exit radiation intensities at boundary surfaces. The conjugate gradient method is used to solve the inverse problem through an iterative procedure. The effects of various variables on source estimation are investigated such as type of source function, errors in the measured data and system parameters, gradient of refractive index across the medium, optical thickness, single scattering albedo and boundary emissivities. The results show that in the case of noisy input data, variation of system parameters may affect the inverse solution, especially at high error values in the measured data. The error in measured data plays more important role than the error in radiative system parameters except the refractive index distribution; however the accuracy of source estimation is very sensitive toward error in refractive index distribution. Therefore, refractive index distribution and measured exit intensities should be measured accurately with a limited error bound, in order to have an accurate estimation of source term in a graded index medium.     

Maximum entropy approach to statistical inference for an ocean acoustic waveguide

A conditional probability distribution suitable for estimating the statistical properties of ocean seabed parameter values inferred from acoustic measurements is derived from a maximum entropy principle. The specification of the expectation value for an error function constrains the maximization of an entropy functional. This constraint determines the sensitivity factor (β) to the error function of the resulting probability distribution, which is a canonical form that provides a conservative estimate of the uncertainty of the parameter values. From the conditional distribution, marginal distributions for individual parameters can be determined from integration over the other parameters. The approach is an alternative to obtaining the posterior probability distribution without an intermediary determination of the likelihood function followed by an application of Bayes' rule. In this paper the expectation value that specifies the constraint is determined from the values of the error function for the model solutions obtained from a sparse number of data samples. The method is applied to ocean acoustic measurements taken on the New Jersey continental shelf. The marginal probability distribution for the values of the sound speed ratio at the surface of the seabed and the source levels of a towed source are examined for different geoacoustic model representations.     

强非线性时间演化声速剖面的序贯反演

受海面波浪起伏、降雨和内波等海洋动力学过程的影响,浅水声速剖面的时间演化具有高度非线性,针对该问题提出使用改进的粒子滤波方法进行声速剖面序贯反演.该方法通过建立声速剖面的经验正交模型(EOF)以及描述声速剖面时间演化特征的状态空间模型,将声速剖面反演问题建模为状态跟踪问题,利用不敏粒子滤波(UPF:Uncented Particle Filter)算法进行声速剖面序贯反演。仿真试验通过实测声速剖面数据和先验地声参数信息产生接收声场数据,再利用模拟声场数据估计声速剖面的时间变化.结果表明,相比于集合卡尔曼滤波(EnKF:Ensemble Kalman Filter),在计算效率等同的情形下,该方法可以在状态参数的时间跳变点保持良好的跟踪性能,一定程度上克服了现有反演算法在跳变点发散的问题,可以有效提高声速剖面反演精度,尤其在声速剖面时变性较强时具有显著优势.      

基于分布源边界点法的局部近场声全息技术

为了克服基于分布源边界点法的近场声全息技术在小全息孔径条件下造成的重建误差问题,提出了基于分布源边界点法的局部近场声全息技术.该技术运用分布源边界点法,采用测得的较小全息面上的声压数据来外推较大全息面上的声压数据,然后用外推的数据进行全息重建.仿真和实验结果验证了采用该技术在小全息孔径条件下进行声源局部重建的有效性.     

声速剖面反演对基阵倾斜失配的敏感性研究

通过推导基阵倾斜失配给声传播时间计算带来的误差,分析了这种参数失配对基于声传播时间的声速剖面反演结果的敏感性。研究结果表明,小角度倾斜时,由基阵倾斜所带来的声速剖面反演误差与倾斜俯仰角近似成正比,与声传播时间近似成反比。当基阵倾斜方位与声源一基阵连线方向一致时带来的误差最大,而垂直时误差最小。理论计算与仿真、实验结果在变化趋势和数量级上具有较好的一致性,为反演误差预报提供了参考。      

Visualization of acoustic radiation from a vibrating bowling ball.

This paper presents visualization of acoustic radiation from a vibrating bowling ball using the Helmholtz equation least squares (HELS) method. In conducting the experiments, the ball is excited by a vibration shaker using stationary random signals. The radiated acoustic pressures are measured using two microphones and taken as input to the HELS formulations. The reconstructed acoustic pressures on the bowling ball surface are compared with those measured at the same locations. Also shown are comparisons of the reconstructed and measured acoustic pressure spectra at various locations on the bowling ball surface. Results demonstrate that the accuracy of reconstruction based on measurements over a conformal surface is much higher than that over a finite planar surface. This is because the latter often extends beyond the near-field region, making the accuracy of measurements inconsistent. Nevertheless, satisfactory reconstruction of acoustic pressure fields over the entire bowling ball surface can still be obtained based on the measurements taken over a finite planar surface on one side of the source. In a similar manner, the normal component of the surface velocity is reconstructed. Once these acoustic quantities are determined, the time-averaged acoustic intensity is calculated. Also presented are the formulations for estimating a priori the numbers of expansion functions and measurements required by the HELS method and the guidelines for determining the reconstruction error and optimum measurement locations, given the overall dimensions of the source and the highest frequency of interest in reconstruction.     

On the selection of loads in the multiload method for measuring the acoustic source parameters of duct systems

The in-duct source can be characterized by two acoustical parameters such as the source strength and the source impedance, which permit the prediction of radiated sound pressure or insertion loss of the whole duct system. One-port acoustic characteristics of an in-duct source can be measured by the multiload method using an overdetermined set of open pipes or side-branch pipes with different lengths as applied loads. The input data, viz. load pressure and load impedance, are usually contaminated by measurement error in the actual measurements, which result in errors in the calculated source parameters. In this paper, the effects of the errors in the input data on the results have been studied numerically, varying the number of loads and their impedances in order to determine what combination of the loads will yield the best result. It is noted that, frequently, only a set of open pipes is used when applying the multiload method to the internal combustion engine sources. A set of pipe lengths, which cause the calculated results to be least sensitive to the input data error, can be found when using open pipe loads. The present work is intended to produce guidelines for preparing an appropriate load set in order to obtain accurate source properties of fluid machines.     

用DIS数字化信息系统测量金属棒中声速的有效方法

本文设计了一种将驻波法与声波传感器相结合测定金属棒中声速的有效方法.将金属棒固定在中点,沿轴向摩擦激励纵振动并在棒中形成驻波而发声,用DIS数字化信息系统接收声波信号,得到波形图,由此测得棒中声波的基频和声速.本文分别测定了铝、钢、黄铜、紫铜4种不同金属材料棒中的声速.结果表明,同种材料、不同长度金属棒中声速的测量结果...     

On the causes of negative source impedance in the measurement of intake and exhaust noise sources

In order to evaluate the acoustic field within a duct system some knowledge is required about the acoustic source, which is characterized by frequency-dependent values of both strength and impedance. However, in using various experimental techniques for various types of source to obtain these values, negative real parts of the source impedance have been measured. Negative source resistance is physically implausible. In some instances, one measurement method gives positive values of source resistance where another gives negative values. This paper reviews the possible causes of the problem, with reference to experimental and theoretical results, in an attempt to clarify the issue. Various possible origins of a negative source resistance have been investigated. It is noted that, for the time being, the violation of the assumption of a time-invariant source and the defect in the inherent algorithm of the load method are the most probable origins of the negative source resistance.     

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.     

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

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

南海北部海域声速剖面水平变化环境声学层析

为获取水平变化的海水声速分布,在等效平均海水声速反演基础上,提出了一种声速剖面水平变化环境声学层析方法。首先利用经验正交函数结合匹配场处理方法反演不同距离下的等效平均声速剖面,进而通过逐步迭代重构出随水平距离变化的声速场分布,并通过数值仿真及海上实验两方面验证了声速剖面水平变化环境声学层析方法的可行性及其准确性。结果表明:声速剖面水平变化环境声学层析结果与实验测量声速剖面均方根误差小于0.8 m/s,利用声传播信号基本能较为准确地重构出随距离水平变化的海水声速分布,为进一步层析中尺度三维结构及声呐水声环境保障应用奠定了基础。