一种基于简正波模态消频散变换的声源距离深度估计方法

针对浅海环境中传播的低频宽带水声脉冲信号,基于简正波水平波数差和波导不变量之间的关系,本文提出了一种利用距离-频散参数二维平面聚焦测距与匹配模态能量定深的目标声源定位方法.首先,通过将由频散参数和波导不变量表示的前几阶模态相速度与由环境模型计算的相速度进行对比分析,从而估计出前几阶模态的频散参数和环境的波导不变量.其次,利用估计出的频散参数值和波导不变量对接收信号进行消频散变换处理,只有当接收信号的距离参数等于目标声源距离时,各号简正波的幅度均达到最大值,在距离-频散参数二维平面上,出现声压聚焦的现象,利用此现象可以估计目标声源的距离.不仅如此,消频散变换后的接收信号,前几阶模态在时域上明显地分离开来,可以准确地估计出前几阶模态的能量,采用多模态能量匹配的方式,可以估计出目标声源的深度.最后,通过对仿真和冬季获得的气枪信号数据处理结果验证了本文方法的有效性.     

一种基于单水听器的浅海水下声源被动测距方法

针对浅海波导中宽带脉冲声源的被动测距问题,本文在群延迟理论的基础上,与warping变换处理相结合,提出了一种适用于浅海波导中宽带声源的单水听器被动测距方法.利用warping变换可以实现对脉冲声源接收信号各阶简正波的分离提取,对分离后的简正波进行时频分析处理可以得到各阶简正波到达时刻和频率之间的关系,即各阶简正波的频散曲线,从而得到任意两阶简正波到达接收水听器的时延差.海底相移参数P是描述海底地声参数的一个重要参量,包含了海底地声参数信息,在海底环境参数未知而P已知的情况下,利用P和简正波水平波数之间的关系可以求得任意两阶简正波的?S_(g,mn)(群慢差).根据群延迟理论,利用得到的任意两阶简正波的时延和?S_(g,mn)可实现利用单水听器对水下声源进行被动测距.本文提出的测距方法测量简单、计算方便,具有较强的实用意义.数值仿真和海上实验数据处理结果的测距误差都在10%以内,证明了该方法的有效性.     

利用简正波频散特征映射方法进行海洋环境参数反演

浅海波导中传播的脉冲信号可以通过频散特征映射方法(warping变换)进行简正波的有效分离,这使得频散曲线的获取更加容易和准确.利用简正波频散特征映射方法分析2010年南海声层析实验中单个阵元接收的声学信号,有效分离出各号简正波,并进一步通过Wigner-Ville分布处理提取出波导中各号简正波的频散曲线.将提取的频散曲线与理论计算到达时差进行匹配,实现了利用单水听器信号的海洋环境参数反演.实验结果表明:将简正波分离与Wigner-Ville分布相结合可有效消除Wigner-Ville分布在估计多分量信号时的交叉干扰,提取的频散特征可很好地用于同时反演海水和海底声速剖面。为评价反演结果的有效性,对反演结果进行了后验概率分析,并将一小时内的海水声速剖面反演结果与温度链实测值作对比,二者符合较好.      

一种基于简正波模态频散的远距离宽带海底参数反演方法

在浅海环境中, 海底环境参数对声传播有着重要的影响. 由于利用单个宽带声源进行海底参数反演时, 随着距离的增大, 误差变大, 本文提出利用warping变换对在浅海波导中传播的, 不同距离上的两个宽带爆炸声源进行简正波的有效分离, 实现了宽带爆炸声源的远距离海底参数反演. 采用全局寻优遗传算法对提取出的模态频散到达时间差与理论计算的模态频散到达时间差进行匹配处理, 并结合随距离连续变化的声传播损失, 实现了利用单水听器进行海底参数的反演. 实验结果表明: 运用反演出的海底参数提取模态频散时间差和实测数据提取出的模态频散时间差吻合得较好; 而通过传播损失反演得到的海底衰减系数与频率呈指数关系. 最后, 对反演结果进行了后验概率分析, 并将本组爆炸声源的反演结果用于另一组不同距离上爆炸声源时仍然有效, 来评价反演结果的有效性.     

利用简正模态相位关系的浅海声源深度分辨方法

针对浅海波导声源深度分辨问题,本文提出了一种利用简正波的相位关系进行声源深度分辨的方法。该方法通过对垂直阵接收到的信号进行简正波分解,得到简正波系数,在已知声源距离的情况下通过对简正波系数项的处理得到各阶简正波的相位关系,根据相位关系估计声源的深度区间。该方法在阵元数较少时依然可以分辨声源深度,而且对宽带信号、窄带信号、连续信号和脉冲信号均适用。理论分析和数值仿真的结果说明:只有近水面处声源激发的简正波满足相位一致的特性,利用此特性可分辨水面水下目标。仿真和海试结果证明:在海水声速随深度不变或缓变的水文环境下,该方法可以准确的分辨声源深度区间,并且在深度采样不充分的情况下仍然适用。      

联合简正波水波和底波频散特性的贝叶斯地声参数反演

大多数基于浅海简正波模态频散数据的地声参数反演方法无法对深层底质声学参数进行可靠估计,究其原因是仅利用了简正波水波频散特征,忽略了与深层底质声学参数密切相关的底波频散特征,因此,本文在分析了包含水波和底波的浅海宽带数据的基础上,结合底波频散特征对深层底质声学参数变化更加敏感的物理特性,实现了基于完整简正波频散特性的贝叶斯地声参数反演,并针对简正波宽带声场模型计算复杂度较高的现实问题,利用变分贝叶斯蒙特卡罗方法的推断优势,完成了未知参数的可靠估计和快速后验分析.仿真和海上实验结果表明:联合简正波水波和底波频散特征数据的贝叶斯地声参数反演,不仅可以有效估计深层底质声学参数,而且降低了其他相关环境参数的估计不确定性.     

消频散阵不变量距离估计方法

针对传统阵不变量方法在低信噪比和远距离测距性能差的问题,提出了消频散阵不变量测距方法。该方法首先对垂直阵接收信号进行消频散波束形成,消除频散影响,压缩脉冲宽度,提高简正波到达时间分辨率;然后,对得到的角度-时间谱采用反卷积处理提高简正波俯仰角度分辨率。通过消频散和反卷积处理,减小了简正波相对到达时间和俯仰角在角度-时间谱图中的聚焦范围,有利于准确的提取阵不变量和估计目标声源的距离。计算机仿真证明了该方法在低信噪比条件下的鲁棒性;2009年崂山湾实验数据处理结果表明:相比于传统阵不变量测距方法,该方法距离估计误差减小2.5%。      

一种基于模态匹配的浅海波导中宽带脉冲声源的被动测距方法

针对浅海波导中宽带脉冲声源的被动测距问题,本文在模态匹配和匹配场处理定位方法的基础上,提出了一种适用于具有液态半无限空间海底的浅海波导中声源的单水听器被动测距方法.利用warping变换可以对脉冲声源接收信号的各阶简正波实现有效分离,由此得到各阶简正波的频域信号.海底相移参数是描述海底地声参数的一个重要参量,包含了海底地声参数信息,而各阶简正波的水平波数可以通过含有海底相移参数的表达式来表达.此外,由于声速剖面对简正波的各阶水平波数具有相近的影响,因此通过对任意两阶简正波进行联合处理,可以近似消除声速剖面对简正波水平波数差的影响.任意两阶简正波的水平波数差只近似用于海底相移参数、海深以及波导中平均声速三个参数有关,可以简单、快速地计算相应拷贝场,然后通过建立代价函数并对简正波模态进行匹配,可以实现对水下脉冲声源的被动测距.与传统的模态匹配定位方法相比,本文提出的方法既不需要使用水听器阵,又可以简单、快速地计算出拷贝场.数值仿真和海上实验数据处理结果的测距误差都在10%以内,证明了该方法的有效性.     

利用自相关函数warping变换的浅海声源深度判别

针对浅海波导声源深度判别问题,提出了一种利用warping变换提取接收信号简正波相关项的深度判别方法。对接收信号自相关函数做warping变换,分离其简正波相关项,利用能量占主导的简正波相关项特征频率之间的比例关系,确定其号数,从而分辨水面声源和水下声源。公式推导和数值仿真结果说明:不同深度声源激发的能量占主导的简正波相关项特征频率不同,可以用于判别声源深度。海试实验结果证明:在海水声速随深度不变或缓变的水文环境,低频条件下,该方法可以分离简正波相关项并确定其号数,在未知声源距离的情况下有效进行声源深度判别,且无须声源相对接收器运动。      

浅海波导中简正波干涉特征频率增强

针对浅海波导中水平线列阵接收的低频宽带低信噪比信号，提出了一种利用多拍信号相干累加来提高干涉简正波特征频谱信噪比的方法。这种方法主要针对未知距离和声源形式的运动声源，对不同时刻（或接收距离）处阵列输出信号进行WARPING变换，基于干涉简正波特征频率不变性原理，通过距离和径向速度比值的搜索来使得不同时刻或距离处信号自相关函数WARPING变换后频谱具有近似相同的特征频率，进而通过相干累加来增强信号干涉简正波特征频率。仿真和海上实验数据分析均验证了方法的性能。     

Nonlinear synthesis of infrasound propagation through an inhomogeneous, absorbing atmosphere

An accurate and efficient method to predict infrasound amplitudes from large explosions in the atmosphere is required for diverse source types, including bolides, volcanic eruptions, and nuclear and chemical explosions. A finite-difference, time-domain approach is developed to solve a set of nonlinear fluid dynamic equations for total pressure, temperature, and density fields rather than acoustic perturbations. Three key features for the purpose of synthesizing nonlinear infrasound propagation in realistic media are that it includes gravitational terms, it allows for acoustic absorption, including molecular vibration losses at frequencies well below the molecular vibration frequencies, and the environmental models are constrained to have axial symmetry, allowing a three-dimensional simulation to be reduced to two dimensions. Numerical experiments are performed to assess the algorithm's accuracy and the effect of source amplitudes and atmospheric variability on infrasound waveforms and shock formation. Results show that infrasound waveforms steepen and their associated spectra are shifted to higher frequencies for nonlinear sources, leading to enhanced infrasound attenuation. Results also indicate that nonlinear infrasound amplitudes depend strongly on atmospheric temperature and pressure variations. The solution for total field variables and insertion of gravitational terms also allows for the computation of other disturbances generated by explosions, including gravity waves.     

Identification of a dynamical model for an acoustic enclosure using the wavelet transform

This paper presents results of research work in the identification of a dynamical model for an acoustic enclosure, a duct with rectangular cross-section, closed ends, and side-mounted speaker enclosures. An acoustic enclosure is excited randomly and random decrement functions are built to convert the random responses to free acoustic responses. It is shown that the estimation of resonance frequencies is possible using the wavelet transform of the system’s free response. Using a particular form of the son wavelet function, results are improved compared to those obtained with the traditionally Morlet wavelet function. An optimal value of a parameter of the son wavelet function is obtained by minimization of the wavelet entropy. The accuracy of this new technique is confirmed by applying it to a numerical example, and to an acoustic enclosure. The advantage of using the wavelet transform method over the Fourier-based modal analysis that would normally be used for the enclosure modes problem is established.     

A k-space Green's function solution for acoustic initial value problems in homogeneous media with power law absorption

An efficient Green's function solution for acoustic initial value problems in homogeneous media with power law absorption is derived. The solution is based on the homogeneous wave equation for lossless media with two additional terms. These terms are dependent on the fractional Laplacian and separately account for power law absorption and dispersion. Given initial conditions for the pressure and its temporal derivative, the solution allows the pressure field for any time t>0 to be calculated in a single step using the Fourier transform and an exact k-space time propagator. For regularly spaced Cartesian grids, the former can be computed efficiently using the fast Fourier transform. Because no time stepping is required, the solution facilitates the efficient computation of the pressure field in one, two, or three dimensions without stability constraints. Several computational aspects of the solution are discussed, including the effect of using a truncated Fourier series to represent discrete initial conditions, the use of smoothing, and the properties of the encapsulated absorption and dispersion.     

Acoustic radiation of cylindrical elastic shells subjected to a point source: investigation in terms of helical acoustic rays.

The paper deals with the acoustic radiation of a cylindrical elastic shell with no internal loading surrounded by a fluid medium when its external surface is subjected to a point source. The problem is addressed via the use of the spatial Fourier transform. An expression is obtained for the radiated pressure that is evaluated for the far field using both the stationary phase method and the fast Fourier transform (FFT). The acoustic field calculated from the FFT is much more complicated than that obtained by using only the stationary phase method. In agreement with the geometrical theory of diffraction (GTD), alternative interpretations of the radiated field in terms of helical acoustic rays allows one to understand the reason for this result. The outstanding phenomenon underlined by the use of the FFT is the emergence of an infinite number of spatial dispersion curves associated with each leaky wave propagating in shells when excited by a point source.     

Reverberation generated by sequential underwater explosions

In order to obtain long-time reverberation, an innovative experiment is operated, using a sequence of explosions with a interval time. The relationship between reverberation generated by single explosion and a sequence of explosions is discussed. The data are obtained from the experiment, then analyzed with method of wavelet transform, conclusion is obtained that power generated by sequence of explosions with different interval time spreads uniformly at low frequency range, unequally at high frequency range.     

The wavelet response as a multiscale NDT method

We analyze interfaces by using reflected waves in the framework of the wavelet transform. First, we introduce the wavelet transform as an efficient method to detect and characterize a discontinuity in the acoustical impedance profile of a material. Synthetic examples are shown for both an isolated reflector and multiscale clusters of nearby defects. In the second part of the paper we present the wavelet response method as a natural extension of the wavelet transform when the velocity profile to be analyzed can only be remotely probed by propagating wavelets through the medium (instead of being directly convolved as in the wavelet transform). The wavelet response is constituted by the reflections of the incident wavelets on the discontinuities and we show that both transforms are equivalent when multiple scattering is neglected. We end this paper by experimentally applying the wavelet response in an acoustic tank to characterize planar reflectors with finite thicknesses.     

Multi-Symplectic Wavelet Collocation Method for Maxwell's Equations

In this paper, we develop a multi-symplectic wavelet collocation method for three-dimensional (3-D) Maxwell's equations. For the multi-symplectic formulation of the equations, wavelet collocation method based on autocorrelation functions is applied for spatial discretization and appropriate symplectic scheme is employed for time integration. Theoretical analysis shows that the proposed method is multi-symplectic, unconditionally stable and energy-preserving under periodic boundary conditions. The numerical dispersion relation is investigated. Combined with splitting scheme, an explicit splitting symplectic wavelet collocation method is also constructed. Numerical experiments illustrate that the proposed methods are efficient, have high spatial accuracy and can preserve energy conservation laws exactly.     

Bending wavelet for flexural impulse response

The work addresses the definition of a wavelet that is adapted to analyze a flexural impulse response of a beam or plate that can be modeled with the Euler-Bernoulli bending theory. The wavelet gives the opportunity to directly analyze the dispersion characteristics of a pulse. The aim is to localize a source or to measure material parameters. An overview of the mathematical properties of the wavelet is presented. An algorithm for the optimal extraction of the dispersion characteristics with the use of genetic algorithms is outlined. The application of the wavelet is shown in an example and experiment.     

Design aspects of acoustic sensor networks for environmental noise monitoring

An increase in public awareness of noise pollution and the impact of noise on human health has led to the need for enhanced insight in complex noise situations. This insight is commonly obtained either by brief measurements or by evaluation of a simplified acoustic model. Both of these approaches however have limitations in complex noise situations. Noise monitoring can be an appropriate and cost efficient measure to obtain more insight, for it allows to measure at many locations for long periods of time. Monitoring can for example be used to improve the accuracy of models or to assess and respond to changes in the acoustic situation. Several monitoring approaches, for various applications, have been developed or are under development. In the first part of this paper an overview is given of current developments in acoustic monitoring networks and their key aspects. Sensor networks for environmental noise monitoring are here divided into four different categories, distinguished by five aspects: hardware costs, scalability, flexibility, reliability and accuracy. These five aspects determine the range of applications for which a network is suited. In the second part of this paper a monitoring network developed in-house is used to further illustrate the relevance of these aspects. This network was designed to facilitate research into the field of acoustic monitoring networks and is used to experiment with and learn from a broad field of applications.     

Schlieren成像和声场可视化

Schlieren技术是利用声场引起透明媒质光折射率的变化而实现声场可视化的光学成像技术。它具有对声场无干扰、快速、瞬态成像的特点。本文利用二维光学Fourier变换分析了Schlieren技术的成像原理,在采用连续激光和高速ICCD的Schlieren成像系统中,实验研究了平面波声场和线聚焦声场中换能器光学校准方法和声压的定量检测技术。发展声场瞬态和动态成像技术,观测了声波的聚焦过程和固-液界面的声场分布和变化。这些结果表明Schlieren技术是一种有效的声场可视化和定量检测的光学成像技术。