A Study of Spatial Processing Gain in Underwater Acoustic Communications

Spatial processing, including beamforming and diversity combining, is widely used in communications to mitigate intersymbol interference (ISI) and signal fading caused by multipath propagation. Beamforming suppresses ISI (and noise) by eliminating multipath (and noise) arrivals outside the signal beam. Beamforming requires the signals to be highly coherent between the receivers. Diversity combining combats ISI as well as signal fading by taking advantage of the independent information in the signal. Classical (spatial) diversity requires that signals are independently fading, hence are (spatially) uncorrelated with each other. In the real world, the received signals are neither totally coherent nor totally uncorrelated. The available diversity is complex and not well understood. In this paper, we study the spatial processing gain (SPG) as a function of the number of receivers used, receiver separation, and array aperture based on experimental data, using beamforming and multichannel combining algorithms. We find that the output symbol signal-to-noise ratio (SNR) for a multichannel equalizer is predominantly determined by the array aperture divided by the signal coherence length, with a negligible dependence on the number of receivers used. For a given number of receivers, an optimal output symbol SNR (OSNR) is achieved by spacing the receivers equal to or greater than the signal coherence length. We model the SPG in decibels as the sum of the noise suppression gain (NSG, equivalent to signal-to-noise enhancement) and the ISI suppression gain (ISG, equivalent to signal-to-ISI enhancement) both expressed in decibels; the latter exploits the spatial diversity and forms the basis for the diversity gain. Data are interpreted using the modeled result as a guide. We discuss a beam-domain processor for sonar arrays, which yields an improved performance at low-input SNR compared to the element-domain processor because of the SNR enhancement from beamforming many sensors.     

A new approach to the design of broad-band element space antenna array processors

In this paper, a new approach to the design of broad-band time-domain element space antenna array processors is presented. The basic approach is based on the idea of minimizing the mean-square deviation between the desired look-direction response and the response of the processor over a frequency band of interest. With this approach, three types of presteering can be handled: no presteering, coarse presteering, and exact presteering. The elimination of presteering lime delays is important in a digital implementation of antenna array processors. The relationship that the new processor has to other broad-band processors is explored and the significance of the parameters associated with the new processor is established. Furthermore, the approach presented enables various types of errors and mismatches between signal model and actual scenario to be incorporated in the problem formulation.     

Optimization of conventional beamformer shading weights for conformal velocity sonar

A numerical optimization technique that uses sonar array noise measurements is used to determine conventional shading weights that maximize the broadband deflection coefficient at the output of the optimal square-law detector, across a frequency band of interest. This process maintains the structure of the conventional processor while providing performance improvement typical of adaptive techniques. The performance of the optimized time domain delay-and-sum beamformer is compared with that of the traditional beamformer that uses conventionally chosen shading weights. Application of this method to conformal velocity sonar array data is shown to provide large improvements in performance over heuristic designs.     

Differences between passive-phase conjugation and decision-feedback equalizer for underwater acoustic communications

Passive-phase conjugation (PPC) uses passive time reversal to remove intersymbol interferences (ISIs) for acoustic communications in a multipath environment. It is based on the theory of signal propagation in a waveguide, which says that the Green's function (or the impulse-response function) convolved with its time-reversed conjugate, summed over a (large-aperture) vertical array of receivers (denoted as the Q function) is approximately a delta function in space and time. A decision feedback equalizer (DFE) uses a nonlinear filter to remove ISI based on the minimum mean-square errors (mmse) between the estimated symbols and the true (or decision) symbols. These two approaches are motivated by different principles. In this paper, we analyze both using a common framework. We note the commonality and differences, and pros and cons, between the two methods and compare their performance in realistic ocean environments, using simulated and at-sea data. The performance measures are mean-square error (mse), output signal-to-noise ratio (SNR), and bit-error rate (BER) as a function of the number of receivers. For a small number of receivers, the DFE outperforms PPC in all measures. The reason for poor PPC performance is that, for a small number of receivers, the Q function has nonnegligible sidelobes, resulting in nonzero ISI. As the number of receivers increases, the BER for both processors approaches zero, but at a different rate. The modeled performance differences (in mse and SNR) between PPC and DFE are in general agreement with the measured values from at-sea data, providing a basis for performance prediction.     

Passive systems theory with narrow-band and linear constraints: Part I--Spatial diversity

This paper studies passive problems where the receiver extracts from the source radiated signature information concerning the parameters defining the relative source/receiver geometry. A model encompassing the fundamental global and local characteristics for passive positioning and navigation is presented. It considers narrow-band signals, imposes linear constraints on the geometry, and exhibits explicitly the symmetry between the space and time aspects. The analysis concentrates on questions of global geometry identifiability, emphasizing the passive global range acquisition. The maximum-likelihood processor is analyzed by studying the ambiguity structure associated with inhomogeneous passive narrowband tracking. Bounds on the global and local mean-square error performance are studied and tested via Monte Carlo simulations. By considering two limiting geometries, a distant and a close observer, simple approximate expressions for the mean-square errors are presented and compared to the exact bounds. Herein the study is restricted to stationary geometries where the source is located by an extended array (spatial diversity). Subsequent papers generalize the study to moving sources (temporal diversity) and to coupled geometries.     

Passive systems theory with narrow-band and linear constraints: Part III -- Spatial/Temporal diversity

This is the last of a series of three papers studying the theory of passive systems. The model assumes that i) the narrow-band signals are transmitted through a Rayleigh channel, ii) the observing array is geometrically linear, and iii) the source motion is deterministic. Ranging techniques based on synchronized measurements of the travel time delay are precluded by the incoherent phase model considered. The paper explores alternative methods that process the phase modulations induced on the signal by the extended geometry and relative dynamics. The present work applies maximum likelihood theory to design the receiver, being concerned with the global identifiability of all parameters defining the relative source/receiver geometry and dynamics. The emphasis is placed on the passive range global acquisition. In contradistinction with the previous papers, where the time stationarity (Part I) or the space homogeneity (Part II) lead to a one-dimensional processor, here the receiver involves processing over both domains. The paper considers the issues of space/time factorability and coupling arising in nonhomogeneous passive tracking. The cross coupling, resulting in more complex filters, improves the receiver acquisition capability. Resorting to Taylor's series type studies, the paper quantifies these improvements, as well as the receiver's mean square error performance, in terms of intuitively satisfying analytical expressions.     

The TAP III beamforming system

The Three-Array Processor (TAP III) beamforming system incorporating both wide-band time-domain beamforming and narrow-band frequency-domain beamforming is described. This paper briefly develops the beamforming theory and shows how the fast Fourier transform (FFT) is applied to accomplish frequency-domain beamforming. The frequency-domain beamformer operates in the frequency domain to form beams and power spectrum data over narrow frequency bands of interest. A real-time digital filtering technique is used to extract the narrow bands of interest from the broad-band input signal. The frequency-domain beamformer accomplishes real-time digital filtering and beamforming by using a high-speed array processor to do the complex calculations and data handling required by the algorithm. The time-domain beamformer operates in parallel with the frequency-domain beamformer to form up to 16 broad-band beams in the time domain. A programmable all-pass digital filter is used to create the fine time delays required by the time-domain beamformer.     

Performance investigation of a two-raft-type wave energy converter with hydraulic power take-off unit

Raft-type wave energy converter (WEC) is a multi-mode wave energy conversion device, using the relative pitch motion to drive its hydraulic power take-off (PTO) units for capturing energy from the ocean waves. The hydraulic PTO unit as its energy conversion module plays a significant role in storing large qualities of energy and making the output power smooth. However, most of the previous investigations on the raft-type WECs treat the hydraulic PTO unit as a linear PTO unit and do not consider the dynamics of the hydraulic circuit and components in their investigations. This paper is related to a two-raft-type WEC consisting of two hinged rafts and a hydraulic PTO unit. The aim of this paper is to make an understanding of the dynamics of the hydraulic PTO unit and how these affect the performance of the two-raft-type WEC. Therefore, a combined hydrodynamic and hydraulic PTO unit model is proposed to investigate and optimize the performance of the two-raft-type WEC; and based on the simulation of the combined model, the relationships between the optimal power capture ability, the optimal magnitude of the hydraulic PTO force and the wave states are numerically revealed. Results show that an approximately square wave type hydraulic PTO force is produced by the hydraulic PTO unit, which causes the performance of the two-raft-type WEC not to be sinusoidal and the energy capturing manner different from that of the device using a linear PTO unit; moreover, there is an optimal magnitude of the hydraulic PTO force for obtaining an optimal power capture ability, which can be achieved by adjusting the parameters of the hydraulic PTO unit; in regular waves, the optimal power capture ability as well as the optimal magnitude of the hydraulic PTO force normalized by the wave height presents little relationship with the wave height, mainly depends on the wave period; in irregular waves, the trends of the optimal power capture ability and the normalized optimal magnitude of the hydraulic PTO force against the peak wave periods at different significant wave heights are generally identical and show a good correlation. All means that the hydraulic PTO unit of the two-raft-type WEC can be tuned to the wave states, and these would provide a valuable guidance for the optimal design of its hydraulic PTO unit.     

The AWSUM filter: a 20-dB gain fluctuation-based processor

The advanced WISPR summation (AWSUM) filter, a natural extension of the WISPR filter for higher filter order numbers, is presented and its performance is compared to the performance of the WISPR filter and the conventional summation processor. It is shown that the AWSUM filter achieves substantial gains in various measures of processor performance above those of the other two processors in spatial and spectral resolution, minimum detectable level (MDL), clutter reduction, and signal-to-noise ratio (SNR) gain. The important processing parameters are shown to be the percentage of overlap of the voltage time series and the number of FFT's averaged. SNR gains in excess of 20 dB were shown to be achievable for low-fluctuation amplitude tonals using measured data     

Geoacoustic Information Content of Horizontal Line Array Data

This paper examines the effectiveness of horizontal line arrays (HLAs) for matched-field inversion (MFI) by quantifying geoacoustic information content for a variety of experiment and array factors, including array length and number of sensors, source range and bearing, source-frequency content, and signal-to-noise ratio (SNR). Emphasis is on bottom-moored arrays, while towed arrays are also considered, and a comparison with vertical line array (VLA) performance is made. The geoacoustic information content is quantified in terms of marginal posterior probability distributions (PPDs) for model parameters estimated using a fast Gibbs sampler approach to Bayesian inversion. This produces an absolute, quantitative estimate of the geoacoustic parameter uncertainties which can be directly compared for various experiment and array factors.     

Approximations to directivity for linear, planar, and volumetricapertures and arrays

Even after decades of sonar design, approximations to the directivity factor (DF) or index of an array, are often used inappropriately. Many of the approximations commonly used provide accurate directivity approximations for only the simplest of array geometries. As the array's size, shape, weighting, and complexity increase, there is a renewed need for better directivity approximations. Directivity is defined as the ratio of the output signal-to-noise (SNR) of an array to the input SNR at an omnidirectional element in a spherically isotropic noise field. Calculation of directivity is obtained by integrating the magnitude-squared response of the array over all angles of incidence. In spherical coordinates, these arrival angles are denoted by an azimuthal angle &thetas; and a polar angle φ. Hence, calculation of the directivity requires a two-fold integration over the angular space defined by the azimuthal and polar angles. For complex, large-size arrays consisting of thousands of array elements, directivity calculations using numerical integration procedures can be time consuming, even on state-of-the-art computing systems. This report provides a number of accurate formulas for estimating the directivity of linear, planar, and volumetric apertures and arrays, which are allowed to have arbitrary shading coefficients, steering angles, and directional array element responses     

一种适用于小信噪比条件的运动声源目标定位方法

本文把长时间积分与匹配场结合起来进行处理,对来自匹配场处理器的能量进行长时间积分,为解决弱信号和衰落信号下的目标检测和定位问题提供了一条新的途径。仿真结果表明：在小信噪比（－３５ｄＢ）情况下仍可对运动声源目标精确定位。     

输出控制方法对Wells透平性能输出特性影响的研究

文中利用实验方法系统研究了负载控制技术对OW C波能转换装置的影响。研究结果表明,当透平完成自起动过程后,在装置输出最佳工况点处,输入波况对平均最佳进气迎角的影响很小,其值基本处于8~10°之间。因此,为实现变工况下波能装置的最佳性能输出,只需保证透平的平均进气迎角始终处于上述范围。反之,利用文中给出的平均最佳进气迎角也可确定透平的最佳转速,进而确定最佳工况点。     

Passive sonar detection and localization by matched velocityfiltering

This paper presents a new bearings-only method of detecting and tracking low signal-to-noise ratio (SNR) wideband targets on a constant course and velocity trajectory. A track-before-detect strategy based on matched velocity filtering is adopted using spatial images constructed from a sequence of power bearing map (PBM) estimates accumulated during a track. To lower the threshold SNR for detection, a discrete bank of matched velocity filters integrates the PBM images over a range of hypothesized trajectories, such an approach eliminates the need to estimate the number of targets since signal detection is determined by comparing the output of each matched filter (MF) to a decision threshold. The distribution of the MF output is derived based on a single point target in diffuse noise assumption. Receiver operating characteristic curves show a definite detection gain under low SNR conditions for matched velocity filtering (track-before-detect) over detection from a single PBM     

Signal excess in K-distributed reverberation

Active sonar systems have recently been developed using larger arrays and broad-band sources to counter the detrimental effects of reverberation in shallow-water operational areas. Increasing array size and transmit waveform bandwidth improve the signal-to-noise ratio-and-reverberation power ratio (SNR) after matched filtering and beamforming by reducing the size of the range-bearing resolution cell and, thus, decreasing reverberation power levels. This can also have the adverse effect of increasing the tails of the probability density function (pdf) of the reverberation envelope, resulting in an increase in the probability of a false alarm. Using a recently developed model relating the number of scatterers in a resolution cell to a K-distributed reverberation envelope, the effect of increasing bandwidth (i.e., reducing the resolution cell size) on detection performance is examined for additive nonfluctuating and fluctuating target models. The probability of detection for the two target models is seen to be well approximated by that for a shifted gamma variate with matching moments. The approximations are then used to obtain the SNR required to meet a probability of detection and false-alarm performance specification (i.e., the detection threshold). The required SNR is then used to determine that, as long as the target and scatterers are not over-resolved, decreasing the size of the resolution cell always results in an improvement in performance. Thus, the increase in SNR obtained by increasing bandwidth outweighs the accompanying increase in false alarms resulting from heavier reverberation distribution tails for K-distributed reverberation. The amount of improvement is then quantified by the signal excess, which is seen to be as low as one decibel per doubling of bandwidth when the reverberation is severely non-Rayleigh, as opposed to the expected 3-dB gain when the reverberation is Rayleigh distributed.     

Using matched-field processing to estimate shallow-water bottomproperties from shot data taken in the Mediterranean Sea

It is extremely difficult to determine shallow ocean bottom properties (such as sediment layer thicknesses, densities, and sound speeds). However, when acoustic propagation is affected by such environmental parameters, it becomes possible to use acoustic energy as a probe to estimate them. Matched-field processing (MFP) which relies on both field amplitude and phase can be used as a basis for the inversion of experimental data to estimate bottom properties. Recent inversion efforts applied to a data set collected in October 1993 in the Mediterranean Sea north of Elba produce major improvements in MFP power, i.e., in matching the measured field by means of a model using environmental parameters as inputs, even using the high-resolution minimum variance (MV) processor that is notoriously sensitive and usually results in very low values. The inversion method applied to this data set estimates water depth, sediment thickness, density, and a linear sound-speed profile for the first layer, density and a linear sound-speed profile for a second layer, constant sound speed for the underlying half space, array depth, and source range and depth. When the inversion technique allows for the array deformations in range as additional parameters (to be estimated within fractions of a wavelength, e.g., 0.1 m), the MFP MV peak value for the Med data at 100 Hz can increase from 0.48 (using improved estimates of environmental parameters and assuming a vertical line array) to 0.68 (using improved estimates of environmental parameters PLUS improved phone coordinates). The ideal maximum value would be 1.00 (which is achieved for the less sensitive Linear processor). However, many questions remain concerning the reliability of these inversion results and of inversion methods in general     

Improved time-delay estimates of underwater acoustic signals usingbeamforming and prefiltering techniques

Passive sonar systems that localize broadband sources of acoustic energy estimate the difference in arrival times (or time delays) of an acoustic wavefront at spatially separated hydrophones, The output amplitudes from a given pair of hydrophones are cross-correlated, and an estimate of the time delay is given by the time lag that maximizes the cross correlation function. Often the time-delay estimates are corrupted by the presence of noise. By replacing each of the omnidirectional hydrophones with an array of hydrophones, and then cross-correlating the beamformed outputs of the arrays, the author shows that the effect of noise on the time-delay estimation process is reduced greatly. Both conventional and adaptive beamforming methods are implemented in the frequency domain and the advantages of array beamforming (prior to cross-correlation) are highlighted using both simulated and real noise-field data. Further improvement in the performance of the broadband cross-correlation processor occurs when various prefiltering algorithms are invoked     

Tidal Turbine Array Optimization Based on the Discrete Particle Swarm Algorithm

In consideration of the resource wasted by unreasonable layout scheme of tidal current turbines, which would influence the ratio of cost and power output, particle swarm optimization algorithm is introduced and improved in the paper. In order to solve the problem of optimal array of tidal turbines, the discrete particle swarm optimization (DPSO) algorithm has been performed by re-defining the updating strategies of particles’ velocity and position. This paper analyzes the optimization problem of micrositing of tidal current turbines by adjusting each turbine’s position, where the maximum value of total electric power is obtained at the maximum speed in the flood tide and ebb tide. Firstly, the best installed turbine number is generated by maximizing the output energy in the given tidal farm by the Farm/Flux and empirical method. Secondly, considering the wake effect, the reasonable distance between turbines, and the tidal velocities influencing factors in the tidal farm, Jensen wake model and elliptic distribution model are selected for the turbines’ total generating capacity calculation at the maximum speed in the flood tide and ebb tide. Finally, the total generating capacity, regarded as objective function, is calculated in the final simulation, thus the DPSO could guide the individuals to the feasible area and optimal position. The results have been concluded that the optimization algorithm, which increased 6.19% more recourse output than experience method, can be thought as a good tool for engineering design of tidal energy demonstration.     

Matched field inversion for geoacoustic model parameters usingadaptive simulated annealing

A method is described for the estimation of geoacoustic model parameters by the inversion of acoustic field data using a nonlinear optimization procedure based on simulated annealing. The cost function used by the algorithm is the Bartlett matched-field processor (MFP), which related the measured acoustic field with replica fields calculated by the SAFARI fast field program. Model parameters are perturbed randomly, and the algorithm searches the multidimensional parameter space of geoacoustic models to determine the parameter set that optimizes the output of the MFP. Convergence is driven by adaptively guiding the search to regions of the parameter space associated with above-average values of the MFP. The performance of the algorithm is demonstrated for a vertical line array in a shallow water enviornment where the bottom consists of homogeneous elastic solid layers. Simulated data are used to determine the limits on estimation performance due to error in experimental geometry and to noise contamination. The results indicate that reasonable estimates are obtained for moderate conditions of noise and uncertainty in experimental geometry     

Digital signal processing for precision wide-swath bathymetry

A mathematical model is formulated which accurately represents the envelope function of bottom return signals received from a number of spatial directions comprising a wide swath. The bottom return signals are processed utilizing a digital nonrecursive matched filter whose coefficients are tapered using a Tukey window. High-speed convolution employing the fast Fourier transform is examined for implementation of the digital matched filtering operation. Computer simulation of the signal processing system indicates that, even in the presence of considerable background and fluctuation noises, the processor provides an output signal having a well-defined peak. The error in time of arrival is found to be less than 3 ms, corresponding to an error in depth of less than 0.1 percent, for an average signal-to-noise ratio of 15 dB and a vertical ocean depth of 12 000 ft (3.7 km). These performance figures apply to the most difficult case of mapping at angles ofpm 45degoff vertical.