Broadband sparse-array blind deconvolution using frequency-difference beamforming

Synthetic time reversal (STR) is a technique for blind deconvolution of receiving-array recordings of sound from an unknown source in an unknown multipath environment. It relies on generic features of multipath sound propagation. In prior studies, the pivotal ingredient for STR, an estimate of the source-signal's phase (as a function of frequency ω), was generated from conventional beamforming of the received-signal Fourier transforms, P(j)(ω), 1 ≤ j ≤ N, where N is the number of array elements. This paper describes how STR is implemented even when the receiving-array elements are many wavelengths apart and conventional beamforming is inadequate. Here, the source-signal's phase is estimated by beamforming P(j)(*)(ω(1))P(j)(ω(2)) at the difference frequency ω(2) - ω(1). This extension of STR is tested with broadband signal pulses (11-19 kHz) and a vertical 16-element receiving array having a 3.75-m-spacing between elements using simple propagation simulations and measured results from the FAF06 experiment involving 2.2 km of down slope propagation from 46 to 92 m water depth. The cross-correlation coefficient between the source-broadcast and STR-reconstructed-signal waveforms for the simulations and experiments are 98% and 91%-92%, respectively. In addition, frequency-difference beamforming can be used to determine signal-path-arrival angles that conventional beamforming cannot.