Feature Extraction for Acoustic Scattering from a Buried Target

Elastic acoustic scattering is important for buried target detection and identification. For elastic spherical objects, studies have shown that a series of narrowband energetic arrivals follow the first specular one. However, in practice, the elastic echo is rather weak because of the acoustic absorption, propagation loss, and reverberation, which makes it difficult to extract elastic scattering features, especially for buried targets. To remove the interference and enhance the elastic scattering, the de-chirping method was adopted here to address the target scattering echo when a linear frequency modulation (LFM) signal is transmitted. The parameters of the incident signal were known. With the de-chirping operation, a target echo was transformed into a cluster of narrowband signals, and the elastic components could be extracted with a band-pass filter and then recovered by remodulation. The simulation results indicate the feasibility of the elastic scattering extraction and recovery. The experimental result demonstrates that the interference was removed and the elastic scattering was visibly enhanced after de-chirping, which facilitates the subsequent resonance feature extraction for target classification and recognition.


Introduction
The detection and recognition of underwater objects buried in the seafloor continue to be a challenge (Chen and Zhou 2012). For a buried target, the common methods can be broadly divided into two categories: imaging detection and non-imaging detection. With imaging sonar (such as synthetic aperture sonar) and image processing, a buried target can be detected and recognized. Studies have shown that the targets and clutter can be sorted according to their acoustic imaging characteristics (Waters et al. 2012b). With localized image symmetry and target strength correlation in 2D feature space, unexploded ordnance (UXO) targets can be separated from at least some non-UXO targets (Bucaro et al. 2008;Bucaro et al. 2012). Compared to imaging detection, non-imaging detection is easier for real-time process, and its requirements on the detection platform are less stringent. In addition, the detection range can be enlarged by transmitting a signal with a lower frequency and smaller grazing angle. Fishell (Fischell and Schmidt 2015) classified seabed spherical and cylindrical targets according to their rigid scattering features. Non-imaging detection implements buried target detection based on its elastic scattering. According to the resonance scattering theory, the scattering from an elastic target contains two terms: background (rigid scattering) and resonance spectrum (elastic scattering) (Junger 1952;Flax et al. 1978). The resonance spectrum reflects the material and internal construction of the target (Gaunaurd and Werby 1991;Morse et al. 1998;Tesei et al. 2002), which are crucial for buried target detection and classification. The elastic contribution to the backscattering from a target remains significant with respect to the specular echo, even when the object is deeply buried (Tesei et al. 2008). The background term has to be removed to obtain the resonance spectrum; Article Highlights • For buried target detection and classification, the weak elastic scattering characteristics need to be improved.
• The de-chirping method is adopted here to address the buried target scattering echo.
• The simulation and experimental results are presented and compared with the MIIR result.
* Xiukun Li lixiukun@hrbeu.edu.cn 1 however, this is rather difficult in practice because the transmission loss is considerable during wave insonification and propagation in seafloor sediment (Wan et al. 2006). In addition, with an active sonar, the presence of reverberation also impedes target detection. There are two popular methods for elastic scattering extraction, namely time inverse technology and the method of isolation and identification of resonance (MIIR). The time inverse technology for buried target detection was introduced by Waters (Waters et al. 2009). The procedure consists of exciting the target object with a broadband pulse, sampling the return with a finite time window, reversing the signal in time, and using this reversed signal as the source waveform for the next interrogation. It has been indicated that the spectrum of the returns rapidly converges to the dominant mode in the backscattering response, and the signal-to-noise of the target echo is enhanced for a target buried in different depths (Waters et al. 2012a;Waters and Barbone 2012;Yu et al. 2014). The MIIR, which was introduced by G. Maze (Maze et al. 1988;Maze et al. 1989), extracts elastic scattering in the time domain by transmitting a narrow pulse. Experimental results have shown that the resonance spectrum is visible with the MIIR, and it is effective for buried target detection (Maze 1991;Décultot et al. 2010).
It is well known that resonance characteristics are significant for target classification; however, the elastic scattering needs to be enhanced as it is usually disturbed by the propagation loss and the interference of reverberation. In this study, to remove the reverberation and specular echo and highlight the elastic resonant characteristics, de-chirping method (Caputi 1971) was adopted when a linear frequency modulation (LFM) signal was transmitted. An echo could be transformed into a series of signals with different frequencies, and the elastic ones could be filtered out with a band-pass filter and then recovered by remodulation. The experimental result shows that the resonance was visibly enhanced after processing, which is beneficial to the next resonance feature extraction for target classification and recognition.

Simulation and Analysis of Buried Target Scattering
The target in this study was a stainless steel sphere with a radius of 0.079 cm, which was buried in the sand-water mixture at a depth of 2 cm. Table 1 lists the selected physical properties for the numerical simulations. A COMSOL model (Hu et al. 2016) was built to obtain the theoretical scattering of the target. Considering the limitation of hardware, the calculation frequency range was set from 60 to 120 kHz with an interval of 100 Hz. The frequency response is shown in Fig. 1. The blue line is the normal series solution of the free-field scattering with a frequency interval of 5 Hz, and the red one is the COMSOL result of the buried target when the attenuation coefficient of sand was 0.5 dB/λ.
Compared to the normal series result of the target in free field, there was no marked change in the frequency response of that in the buried condition, except a slight intensity decrease due to the absorption of sand. A Hanning window was used as a band-pass filter, and the waveform of the inverse Fourier transformation (Anderson et al. 2012) is shown in Fig. 2, where the left one is the waveform of the free-field scattering, and the right one is that of the buried target. For both free-field and buried target scattering, series of narrowband energetic arrivals follow the first specular one, i.e., elastic echo, while their amplitudes rapidly drop over time. Then, the MIIR was employed, and the specular and the elastic scattering were separated in the time domain, the spectra of which are shown in Fig. 3. In the figure, the upper panel gives the result of free-field scattering, where the left is the spectrum of the specular echo, and the right is that of the elastic ones; similarly, the bottom panel is the result of the buried situation. For free-field scattering, the specular echo spectrum was consistent with the incident signal, while there were more details of elastic scattering after the MIIR. The same results were obtained for the buried situation.
The above analysis was performed without considering the incident signal, which is important for signal processing. In underwater studies, LFM signals are very popular because of their pulse compression. Suppose the frequency modulation range of the transmitted LFM signal was from 60 to 120 kHz, and the pulse width was 100 sampling points; the corresponding waveform and elastic spectrum are shown in Fig. 4. The Stainless steel 5940 3100 7900 Fig. 1 Simulation of the scattering from a sphere in free field and in a buried condition aliasing of the rigid scattering and the elastic ones is expected to be noticed when we increase the pulse width of the incident signal. Figure 5 presents the scattering waveforms when the pulse widths were 1000 and 2000 sampling points. With a larger pulse width, the specular echo and elastic scattering were mixed in the time domain. In this situation, it was not easy to extract the exact elastic echo in the time domain with the MIIR; thus, the de-chirping method was employed.

De-chirping Processing
In engineering, a target can be considered as a linear timeinvariant system and the target echo as the system response to the incident signal. In the far field, according to the highlight model (Tang 1994), a target echo can be written as the superposition of a series of subwaves, i.e., highlights. The frequency modulation properties of specular echo and elastic scattering are basically the same , although they have different frequency responses. To simplify the analysis, every highlight is considered as a time-delay copy of the incident signal, without considering the amplitude response and phase jump. An incident LFM signal can be written as where k = (f 1 − f 0 )/T is the slope of frequency modulation, f 0 ∼ f 1 is the frequency range, and T is the pulse width. A single scattering component can be described as where τ i is the corresponding time delay. An LFM signal with a fixed time and the same frequency and frequency modulation are taken to execute a conjugate convolution with the signal to be processed. Select a reference point between the target and receiver, which is R ref away from the receiver, then the time delay is τ ref = 2R ref /c, where c is the sound speed. The reference signal is  (2) and we will get which can be seen as a single-frequency signal, and the frequency is proportional to the time delay . If τ ref = 0, then Eq. (4) can be rewritten as.
This way, scattering with different time delays from a buried target can be seen as a series of single-frequency signals with different frequencies. Therefore, through a band-pass filter, the elastic components of the target scattering can be extracted, and the rigid scattering and interference can be removed. The procedure consists of the following steps: 1) De-chirping the observed signal: With the known parameters of the incident signal, the scattering signal is dechirped as mentioned above. 2) Band-pass filtering the de-chirped signal: The elastic components are filtered out by several band-pass filters, and the center frequencies of the filters are determined based on the peak frequencies of the de-chirped signal. 3) Recovering the filtered component: The elastic scattering components are recovered by remodulating the filtered signal with the transmitted signal.
To verify the proposed method, we consider an LFM signal whose frequency range was 60 to 120 kHz and pulse width was 100 sampling points, and the result of de-chirping is shown in Fig. 6. Figures 6 (a)

Experimental Signal Processing
To further verify the effectiveness of the de-chirping method, an experimental signal was processed. The experiment was conducted in a pool at Shanghai Jiaotong University, and the experimental setup and the parameters are shown in Fig. 7 (Hu et al. 2016). A stainless steel sphere was buried in sand at a depth of 2 cm. The transmitter and receiver were joined into a monostatic device, which could move horizontally through the guide rail above the target. Before burying the sphere in the sand-water mixture, the backscattering in the free field was collected. The signal mentioned in Section 3 was transmitted, and the scattering echo is shown in Fig. 8. Figure 8(a) is the waveform of the backscattering echo and the envelope after matched filtering. As the object was suspended in the pool with a rope, the second arrival was the scattering of the rope. Figure 8(b) is the spectrum, (c) is the elastic scattering obtained by the MIIR, including the rope scattering, and (d) is the pure elastic scattering obtained by the de-chirping method. The resonance peaks in the spectrum of the entire echo are very close to the theoretical result shown in Fig. 1.