等声速环境中目标前向声散射简正波耦合的垂直阵空域响应特征

水下目标的前向声散射会引起声波传播过程中的简正波耦合效应,使得接收声场结构发生变化,研究目标前向散射引起的垂直阵列上空域响应变化特征,可实现对直达波强干扰背景下的前向散射检测。通过将垂直阵波束形成技术分别用于信道中目标散射场理论模型计算数据和湖上实验验证数据,分析了等声速环境中目标前向声散射简正波耦合的垂直阵空域响应特征。结果表明,目标靠近接收端时前向散射引起的声波垂直达到结构与无目标时相比差异显著,高阶简正波向低阶简正波转化导致信号到达时延宽度展宽,采用指向水平方向的窄波束可显著提取目标前向散射引起的接收声波变化特征。      

浅海波导中目标散射的简正波-Kirchhoff近似混合方法

提出了计算浅海波导中复杂目标散射的数值方法:简正波-Kirchhoff近似混合方法。通过把目标散射的Kirchhoff近似方法和简正波声传播模型相结合,可对大尺寸复杂目标在浅海波导中的散射声场进行计算。以浅海波导中刚性球体散射的解析解为标准解验证了本方法,说明简正波-Kirchhoff近似混合方法是有一定计算精度的工程预报方法。数值计算Pekeris浅海波导中球体目标散射声场在深度-频率平面和距离-频率平面上的二维干涉结构及其与自由空间中的差异。进而通过FFT获得目标时域回波随深度的分布图,分析浅海波导中目标姿态、声速剖面对Benchmark潜艇目标散射的影响。      

水平变化波导中的简正波耦合与能量转移

针对海底地形水平变化对声场能量传播和声场干涉结构的影响,对简正波之间的耦合和能量转移进行了研究.建立了一种二维大步长格式的耦合简正波模型和三维楔形波导耦合简正波模型,以便快速有效地分析简正波之间的耦合和能量转移.基于耦合简正波模型,阐述了前向声场能量在水平变化波导中传播时的转移过程.并根据射线简正波理论,解释了海底地形变化对声场能量分布的影响机理.水平变化波导中声场的仿真计算表明,当本征值虚部发生剧烈变化时声场存在着较强的简正波耦合和能量转移,且海底地形变化将导致声场能量的水平传播方向偏转至海水深度增加的方向.在声场能量转移和传播方向变化中,声场的能量趋于保留在波导中而不向海底泄漏.同时,声场能量分布受到类似于压缩或稀疏的作用,从而形成椭圆状的干涉结构.     

用于浅海有源声呐目标深度估计的匹配相位处理EI北大核心CSCD

针对浅海波导中有源声呐目标深度估计问题,提出了通过单个声源发射低频宽带信号,垂直双接收水听器接收目标回波,利用宽带目标回波比值的相位特征进行匹配相位处理的方法,不仅消除了目标散射特性对有源声呐目标深度估计的影响,而且仅需计算单程的信道传递函数,运算量小,有利于实时处理的实现需求。首先通过理论推导,定义了一个隐含目标深度、且与目标散射特性无关的声场特征——单程传播向量宽带干涉结构,在Pekeris波导条件下仿真分析了简正波阶数与单程传播向量宽带干涉结构的幅度和相位的关系,发现选取较多的简正波阶数贡献的单程传播向量宽带干涉结构进行匹配相位处理可以提高目标深度估计的性能。进一步仿真分析表明,在发射信号选取一定的时间长度和带宽的条件下,利用全部阶简正波且信噪比大于-10 dB时,方法的深度估计误差在5 m之内。最后分析了信号时间长度和处理带宽对有源声呐目标深度估计性能的影响,以及海底声速、海深和声速剖面失配时方法的鲁棒性。     

水下任意刚性散射体对Bessel波的散射特性分析

本文利用Bessel波的谐波展开式, 采用T矩阵方法的推导思路, 建立了水下任意刚性散射体在Bessel波照射下的声散射场计算公式. 以水下刚性椭球体和两端附连半球的刚性圆柱体为例, 计算了在不同波锥角β 下的反向散射形态函数, 同时, 依据镜反射波和绕行波的干涉物理模型, 给出了预报Bessel波照射下的反向散射形态函数峰峰间隔值的计算模型. 仿真结果表明本文提出的Bessel波照射下反向散射形态函数峰峰间隔值预报方法是准确有效的, 同时也说明, 本文建立的基于T矩阵法计算水下任意刚性散射体在Bessel波束下的声散射场方法是有效的, 这拓展了T矩阵法的应用领域.     

浅海内波影响下的波导不变量变化特性分析

针对浅海内波引起波导不变量变化的问题,利用声场波导不变量的概率分布并结合声场简正波的理论,研究了内波活动下波导不变量的时变性,给出了波导不变量变化的机理和规律.具体结论是,在负跃层波导中,声场的波导不变量的最大概率取值具有明显的频变特性.内波环境下,当声传播方向与内波波阵面平行时,接收声场简正波的幅度变化不大,但是简正波的相慢度差和群慢度差的变化却能引起波导不变量最大概率取值的变化;而当声传播方向与内波波阵面垂直时,内波引起的简正波耦合同样会导致波导不变量的最大概率取值的明显变化.     

Pekeris波导中简正波的复声强及其应用

在Pekeris波导模型下,关注了简正波的矢量场,讨论了简正波水平复声强和垂直复声强的表述,并分析了其特征.单阶简正波在水平方向是行波,相应的水平复声强仅为有功的;在垂直方向为驻波,相应的垂直复声强仅为无功的.而多阶简正波相互干涉,因此总声场的复声强既有有功分量,也有无功分量,其中只有有功分量参与声能的输运,但无功分量是反映声场信息的重要组成部分.通过对垂直(交互)复声强无功分量和水平交互复声强有功分量的数值分析,对于甚低频率的点源声场,发现当声源深度变化时,上述声场分量的正负号呈有规变化,当接收传感器置 关键词： 目标深度分类 复声强 矢量场 Pekeris波导     

一种水平变化可穿透波导中声传播问题的耦合简正波方法

通过利用标准简正波程序KRAKEN计算本地简正波解及耦合矩阵, 进一步发展了求解水平变化波导中声场的全局矩阵耦合简正波方法(Luo et al., "A numerically stable coupled-mode formulation for acoustic propagation in range-dependent waveguides," Sci. China-Phys. Mech. Astron. 55, 572 (2012)), 使得该方法可以处理具有可穿透海底及随深度变化声速剖面等实际问题, 并提供声场的完全双向解. 本文还给出了双层波导中耦合矩阵的解析表达式, 并利用其验证了本方法中耦合矩阵数值算法的精度. 最后, 利用改善后的全局矩阵耦合简正波模型(DGMCM)计算了美国声学学会(ASA)提出的可穿透楔形波导标准问题, 将所得数值解与参考解比较, 结果表明DGMCM方法可以精确处理水平变化波导中声传播实际问题. 关键词： 耦合简正波理论 全局矩阵方法 可穿透楔形波导     

SBS相位共轭补偿现象的实验研究与图示化分析

通过实验观察了受激布里渊散射相位共轭波前补偿现象,分析了单横模激光、波前弱畸变及强畸变激光受激布里渊散射现象的差异.由于工作介质中强烈的热畸变所引起的波前破碎效应,高功率固体激光器的输出光束一般都不是理想的高斯型单横模,而是具有弱畸变或强畸变波前结构的光束,其包含多个子光束.只有总能量足够大的激光束被聚焦到布里渊介质中时,其包含的所有子光束才能各自按单横模受激布里渊散射规律产生后向散射光,再合并而实现完全的相位共轭波前反转,否则就只能实现部分相位共轭波前反转,即部分相位共轭补偿.     

浅海中利用单水听器的声源被动测距

针对浅海声波导中远距离脉冲声源被动测距问题,提出了一种利用单水听器接收信号自相关函数进行warping变换的声源被动测距方法。理想水下声波导中,接收信号warping变换输出的傅里叶变换频谱中具有不变性频率特征,即与声源距离无关的各简正波截止频率;信号自相关函数中不同简正波相干成分也存在不变性频率特征;推导了未知声源距离时特征频率提取值与不变性频率特征之间的近似关系式。这些规律可推广到实际浅海声波导,并用于声源被动测距。利用声场计算模型来提供具有不变性频率特征的频谱,对2011年12月北黄海海域水声实验中单水听器接收的脉冲声数据进行了处理,验证了方法的有效性,测距结果和实际距离符合良好,平均测距误差在10%以内。      

Efficiency of Spatial Matched Filtering of Shallow-Water Sound Waveguide Modes

Peculiarities of the mode shadow formation in a randomly inhomogeneous oceanic waveguide are studied using vertical radiating and receiving arrays. It is shown that the maximum mode shadow contrast can be reached for optimal parameters of the radiating and receiving arrays if characteristics of the waveguide at the points of location of the source and receiver are known. A decrease in the mode shadow contrast due to scattering of few-mode signals by random inhomogeneities in oceanic waveguides is analyzed.     

Numerical simulation of the EM scattering and Doppler spectrum of a low flying target above dynamic oceanic surface by using the FEM-DDM method

The Doppler (DP) spectrum of a moving target such as a ship, an iceberg or an air-plane above dynamic oceanic surfaces[1—3] is one of the most important subjects for ra-dar oceanic surveillance, target tracking and oceanic remote sensing[4—7]. With the ad-vancement of oceanic remote sensing and radar surveillance, the experimental observa-tion and theoretical modeling of oceanic clutter have been extensively studied, e.g. by using oceanic field measurement, wave tank experiment and some stud…     

Spectral and modal formulations for the Doppler-shifted field scattered by an object moving in a stratified medium

Spectral and normal mode formulations for the three-dimensional field scattered by an object moving in a stratified medium are derived using full-field wave theory. The derivations are based on Green's theorem for the time-domain scalar wave equation and account for Doppler effects induced by target motion as well as source and receiver motion. The formulations are valid when multiple scattering between the object and waveguide boundaries can be neglected, and the scattered field can be expressed as a linear function of the object's plane wave scattering function. The advantage of the spectral formulation is that it incorporates the entire wave number spectrum, including evanescent waves, and therefore can potentially be used at much closer ranges to the target than the modal formulation. The normal mode formulation is more computationally efficient but is limited to longer ranges. For a monochromatic source that excites N incident modes in the waveguide, there will be roughly N2 distinct harmonic components in the scattered field. The Doppler shifts in the scattered field are highly dependent upon the waveguide environment, target shape, and measurement geometry. The Doppler effects are illustrated through a number of canonical examples.     

Acoustical scattering by arrays of cylinders in waveguides

Multiple scattering of acoustic waves in a planar horizontal waveguide by finite-length cylinders is considered. Cylinder height equals the waveguide depth, and both are vertically constrained by the pressure-release boundaries. An analytically exact solution is obtained via normal mode expansion method in conjunction with the concept of the T matrix. The problem is decomposed into an infinite number of two-dimensional multiple scattering problems, modulated by waveguide mode shapes. Examples are presented for an isovelocity waveguide in which the medium is uniform and the waveguide depth is constant. It is found that, in numerical computations, including one or two evanescent modes captures the essence of the evanescent modes. Multiple scattering in the waveguide is compared with the corresponding two-dimensional case. It is concluded that, in low frequencies, the wave patterns in the two cases are very similar, with a shift in the frequency. The similarity diminishes when there are multiple propagating modes. Despite the mode mixing, some key features in the scattering as observed in the two-dimensional problem remain observable in the waveguide case.     

运动目标辐射声场干涉结构映射域特征研究

浅海低频声场呈稳定而显著的干涉特征,并蕴含声源状态和波导特性等信息.本文研究运动目标辐射声场干涉结构的简化映射以及映射域能量分布对目标运动状态的指示特征.理论分析了声压场和矢量声场空(时)频干涉谱图的二维傅里叶变换映射特征,推导了匀速运动目标时频干涉谱图映射域能量脊斜率与波导不变量以及距变率、航向角的关系式,证明了映射域脊斜率符号、脊斜率绝对值变化等对目标来袭或远离以及目标威胁程度的指示,并进行了数值仿真和海上试验研究.实测结果与理论、仿真分析有较好的一致性.研究结果表明:二维傅里叶变换可将声压场和矢量场时(空)频谱图干涉结构简化,匀速运动目标辐射声场干涉结构映射域的能谱脊斜率、距变率、航向角与波导不变量有解析关系式,声强谱、动能密度谱、声强流谱等声场干涉结构经映射后更为一致,映射域脊能简明的指示目标运动状态和威胁程度.     

A waveguide invariant adaptive matched filter for active sonar target depth classification

This paper addresses depth discrimination of a water column target from bottom clutter discretes in wideband active sonar. To facilitate classification, the waveguide invariant property is used to derive multiple snapshots by uniformly sub-sampling the short-time Fourier transform (STFT) coefficients of a single ping of wideband active sonar data. The sub-sampled target snapshots are used to define a waveguide invariant spectral density matrix (WI-SDM), which allows the application of adaptive matched-filtering based approaches for target depth classification. Depth classification is achieved using a waveguide invariant minimum variance filter (WI-MVF) which matches the observed WI-SDM to depth-dependent signal replica vectors generated from a normal mode model. Robustness to environmental mismatch is achieved by adding environmental perturbation constraints (EPC) derived from signal covariance matrices averaged over the uncertain channel parameters. Simulation and real data results from the SCARAB98 and CLUTTER09 experiments in the Mediterranean Sea are presented to illustrate the approach. Receiver operating characteristics (ROC) for robust waveguide invariant depth classification approaches are presented which illustrate performance under uncertain environmental conditions.     

Dynamics of an optical signal in a multimode waveguide structure

A mathematical model of linear dynamic distortions of average signal power in a planar optical waveguide stemming from the intermodal coupling and dispersion of waveguide modes caused by macroscopic bends and irregularities of the optical waveguide surface has been developed on the basis of pulse intermodal scattering matrix. It is shown that the elements of this matrix are determined at a waveguide segment whose length greatly exceeds the correlation scale of the field of irregularities but is much smaller than the normalization length of waveguide modes and are described by rectangular functions with a duration proportional to the relative velocity of the corresponding modes. The dependences of the pulse intermodal scattering matrix on the intensities of absorption, Rayleigh scattering, and optical signal scattering from waveguide irregularities have been established. It is shown that, at waveguide lengths exceeding the normalization length, the optical signal shape becomes almost stable and can be expressed in terms of the pulse characteristic of an integrating circuit. It has been established that irregularities of the film in a bent optical waveguide lead not only to additional attenuation but also to enhancement of multimode distortion. A matrix method is proposed for taking into account these distortions, which is based on the effect of equalization (scrambling) of the signal energy over the spectrum of the bent waveguide modes. The model developed was used to process the experimental data on the dynamics of a signal in a weakly guiding multimode polymer optical fiber in order to estimate the parameters of attenuation and irregularities of the waveguide surface.     

Time-frequency analysis and modeling of the backscatter of categorized dolphin echolocation clicks for target discrimination

A set of dolphin echolocation clicks collected from an Atlantic bottlenose dolphin in Kaneohe Bay, Hawai'i from a previous experiment is examined in terms of their time and frequency characteristics. The center frequency and rms bandwidth are calculated for the clicks and these are clustered into four classes by using a model based on the Bayesian information criterion. The echo signatures are attained from a solid, elastic homogeneous sphere for each class of clicks from an acoustic scattering model. The results from the scattering model are compared to experimental values. The joint time-frequency content of the resulting echo signals is obtained by the reduced interference distribution (RID). The RIDs are plotted and examined for each signal class for four spherical targets of different material compositions. RID correlation values are obtained for a standard target versus comparison targets by using a time-frequency correlator. The results suggest that dolphins may discriminate by auditory inspection of the time-frequency information returned by the targets. The modification of the outgoing clicks and examination of time-frequency target information may be fundamental to a dolphin's ability to identify and discriminate targets.