水平阵信号压缩感知用于简正波分离

针对水平阵信号简正波分离过程中常规波束形成分辨率低以及warping模态滤波不适用于复杂声信号的问题,提出水平阵信号压缩感知用于简正波分离的方法。利用压缩感知在方位估计中的高分辨特性,通过估计水平阵接收信号在频率方位角上的二维分布,分离得到各阶简正波的方位谱,并逆Fourier变换得到时域波形。仿真孔径1 km、阵元间隔10 m水平阵接收20~200 Hz伪随机声信号和脉冲声信号,所提方法分离出的各阶简正波与理论波形的相关系数在0.97~1.0。对2011年北黄海声学实验中的海底28元水平阵接收的气枪信号,在合成至1 km孔径后使用压缩感知方法分离简正波,其与warping模态滤波分离得的前5阶简正波相关系数在0.82~0.93。仿真与实验都说明了水平阵信号压缩感知简正波分离方法的有效性。      

warping变换提取单模态反演海底衰减系数

为了获得浅海海底地声模型参数,利用warping变换方法分离出单模态简正波.对于接收深度固定、定深爆炸声源情况,以简正波理论为基础定义了距离归一化的简正波传播损失,并且其随传播的距离呈线性关系,故可通过此变化规律得到声压值实部的衰减因子,进而可求得海底地声模型参数:海底衰减系数.为验证此方法的有效性,仿真了warping变换提取单模态简正波的过程,同时将warping变换提取的单模态简正波与数值计算的结果进行比较验证;并针对某次黄海试验数据进行了处理,得到在150—550 Hz频带范围内海底衰减随频率的变化规律为α=0.581f_k~(1.86)(dB/m).通过与其他学者在相同海域试验结果的对比验证,变化规律基本相同.此外不同模态间反演相同频点的衰减系数接近也较好地支撑了结果.     

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

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

利用气枪声源数据的地声参数反演*

本文通过对东海采集的气枪声源数据进行处理,基于地声参数对不同声场观测量的敏感性差异,利用分步反演的策略进行分析。首先选择距离接收阵七个不同的距离点数据,通过warping变换准确提取第三阶至第八阶简正波的频散曲线;利用海底声速对频散曲线敏感的特性来反演声速;由Hamilton经验公式求得海底密度;通过传播损失拟合获得海底衰减。不同距离点数据反演的海底声速与密度一致性较高。实验提取的频散曲线和和反演参数仿真结果、实验获取的传播损失与反演参数获得传播损失均拟合较好。     

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

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

一种地声参数的联合反演方法

根据地声参数对不同声场物理量影响不同,提出了一种利用简正波频散特征结合声传播损失反演地声参数的联合反演方法。首先,考虑到简正波的频散特性(群速度)对海底的密度和声速较为敏感,而对海底吸收系数不敏感,利用自适应时频分析方法,获得不同频率不同号数简正波的到达时间差,以此作为代价函数,采用全局优化算法,反演出海底密度和海底声速的分层结构,并用概率统计的方法评价反演结果的有效性。反演出海底密度和声速后,利用实验船辐射噪声得到随距离连续的声传播损失来反演出海底吸收系数。最后,把反演的参数很好的用于声源匹配定位验证了反演结果的有效性。      

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

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

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

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

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

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

2005黄海实验混响垂直相关统计反演海底参数

处理并分析了2005黄海混响实验数据.在单层海底假设下利用300～800 Hz频段混响垂直相关数据反演了海底参数,并分析了参数反演结果的不确定性.结果发现:单频独立反演的声速值随频率增加而减小.类似的结果以前也有报道,导致难以有效地确定海底的平均声速.为了解决这一问题,本文假设海底具有两层分层结构,并提出了一种基于遗传算法的多频联合海底参数反演方法.同一海域声传播数据的简正波分离结果以及混响垂直相干实验值与基于反演参数计算的理论值吻合良好,验证了分层海底假设及参数反演结果的有效性.     

A passive source ranging method based on the frequency warping transform of the vertical intensity flux in shallow water

The phase of cross-correlation function of two different normal modes contains source range information, which can be extracted by warping transform due to the dispersive characteristics of the shallow water waveguide. The autocorrelation function of the received pressure or particle velocity contains both modal autocorrelation component(MAC) and modal cross-correlation component(MCC), with the former part usually treated as interference for source ranging. Because the real part of the vertical intensity flux(RPVIF) only contains MCC, a passive impulsive source ranging method based on the frequency warping transform of RPVIF with a single vector receiver in shallow water is presented. Using a waveguide-invariant-based frequency warping operator, the cross-correlation components of two different modes in the vertical intensity flux are warped into separable impulsive sequences, the time delays of which are subsequently used for source ranging. The advantages of source ranging based on warping the vertical intensity flux compared with warping the pressure autocorrelation function are pointed out, and the experiment results are also presented.     

Modal analysis of broadband acoustic receptions at 3515-km range in the North Pacific using short-time Fourier techniques

In 1995-1996 the Acoustic Thermometry of Ocean Climate (ATOC) experiment provided an opportunity to study long-range broadband transmissions over a series of months using mode-resolving vertical arrays. A 75-Hz source off the California coast transmitted broadband pulses to receiving arrays in the North Pacific, located at ranges of 3515 and 5171 km. This paper develops a short-time Fourier transform (STFT) processor for estimating the signals propagating in the lowest modes of the ocean waveguide and applies it to analyze data from the ATOC experiment. The STFT provides a convenient framework for examining processing issues associated with broadband signals. In particular, this paper discusses the required frequency resolution for mode estimation, analyzes the broadband performance of two standard modal beamforming algorithms, and explores the time/frequency tradeoffs inherent in broadband mode processing. Short-time Fourier analysis of the ATOC receptions at 3515 km reveals a complicated arrival structure in modes 1-10. This structure is characterized by frequency-selective fading and a high degree of temporal variability. At this range the first ten modes have equal average powers, and the magnitude-squared coherence between the modes is effectively zero. The coherence times of the peaks in the STFT mode estimates are on the order of 5.5 min. An analysis of mean arrival times yields modal dispersion curves and indicates that there are statistically significant shifts in travel time over 5 months of ATOC transmissions.     

Ultrasonic flaw detection using threshold modified S-transform

Interference noising originating from the ultrasonic testing defect signal seriously influences the accuracy of the signal extraction and defect location. Time–frequency analysis methods are mainly used to improve the defects detection resolution. In fact, the S-transform, a hybrid of the Short time Fourier transform (STFT) and wavelet transform (WT), has a time frequency resolution which is far from ideal. In this paper, a new modified S-transform based on thresholding technique, which offers a better time frequency resolution compared to the original S-transform is proposed. The improvement is achieved by the introduction of a new scaling rule for the Gaussian window used in S-transform. Simulation results are presented and show correct time frequency information of multiple Gaussian echoes under low signal-to-noise ratio (SNR) environment. In addition, experimental results demonstrate better and reliable detection of close echoes drowned in the noise.     

Passive source range estimation with a single receiver in shallow water

An approach for long-range passive impulsive source ranging with a single receiver in shallow water is proposed,which utilizes the frequency spectrum of the warped signal autocorrelation function via warping transform.For an ideal waveguide,there are invariable frequency features both in the frequency spectrum of the warped signal corresponding to modal cut-off frequencies and the warped signal autocorrelation function due to modal interference.These intrinsic frequency features can be used to passive source ranging.So,the approximate relationship between the frequency of warped signal at an unknown source range and the intrinsic frequency extracted by the time warping transform is derived.These rules can be generalized to an actual shallow water waveguide.Employing an acoustic model to offer the invariable frequency spectrum features,the impulsive signal data collected by a single hydrophone in the North Yellow Sea in December 2011 are analyzed to verify the proposed source ranging approach.The estimated ranges are in good agreement with the ranges measured by GPS,and the mean relative error of range estimation is less than 10%.     

Bayesian geoacoustic inversion using wind-driven ambient noise

This paper applies Bayesian inversion to bottom-loss data derived from wind-driven ambient noise measurements from a vertical line array to quantify the information content constraining seabed geoacoustic parameters. The inversion utilizes a previously proposed ray-based representation of the ambient noise field as a forward model for fast computations of bottom loss data for a layered seabed. This model considers the effect of the array's finite aperture in the estimation of bottom loss and is extended to include the wind speed as the driving mechanism for the ambient noise field. The strength of this field relative to other unwanted noise mechanisms defines a signal-to-noise ratio, which is included in the inversion as a frequency-dependent parameter. The wind speed is found to have a strong impact on the resolution of seabed geoacoustic parameters as quantified by marginal probability distributions from Bayesian inversion of simulated data. The inversion method is also applied to experimental data collected at a moored vertical array during the MAPEX 2000 experiment, and the results are compared to those from previous active-source inversions and to core measurements at a nearby site.     

VORTEX-INDUCED VIBRATIONS OF TWO SIDE-BY-SIDE EULER-BERNOULLI BEAMS

Vortex-induced vibration of two side-by-side elastic beams in a cross flow is numerically studied. The two beams are identical and fixed at both ends. In the numerical approach, the Euler-Bernoulli beam theory is used to model the beam vibration, and the laminar Navier-Stokes equations are solved to give the flow field. The flow equations are resolved using a finite element method and the flow-induced forces are calculated at every time step in order to correctly reflect the fluid-beam interaction. The beam response is calculated using the modal analysis method. Free vibrations of the two beams with three pitch ratios, T/D=1·13, 1·7 and 3·0, where T is the gap between the centers of the two beams and D is the beam diameter, are simulated at Re=800. Results obtained are compared with experimental measurements and other numerical results obtained assuming a two-degree-of-freedom (2-d.o.f.) model. The agreement is good in general. Correlation analysis is carried out, showing that the phase relation is different for differentT /D. The short-time Fourier transform (STFT) method is used to carry out the spectral analysis, along with the conventional auto-regressive moving averaging (ARMA) method for comparison. The STFT analysis shows that the time evolution of fluid force and beam vibration for T/D=1·13 and 3·0 are stationary. For theseT /D ratios, the STFT results are consistent with the ARMA results, but give a clearer picture of the higher order harmonics. For T/D=1·7, the time evolution is non-stationary. The STFT analysis shows that there are three types of frequency spectrum for the fluid force, with one, two, and three dominant frequencies respectively. The spectra intermittently change in a random way during the evolution. The ARMA results, though consistent with previous experiments, can only reveal a particular feature of the three different types of spectrum. This suggests that the STFT method is more appropriate to analyze the spectra of non-stationary time series in the study of flow-induced vibrations.