IMF-Slices for GPR Data Processing Using Variational Mode Decomposition Method

Abstract

Using traditional time-frequency analysis methods, it is possible to delineate the time-frequency structures of ground-penetrating radar (GPR) data. A series of applications based on time-frequency analysis were proposed for the GPR data processing and imaging. With respect to signal processing, GPR data are typically non-stationary, which limits the applications of these methods moving forward. Empirical mode decomposition (EMD) provides alternative solutions with a fresh perspective. With EMD, GPR data are decomposed into a set of sub-components, i.e., the intrinsic mode functions (IMFs). However, the mode-mixing effect may also bring some negatives. To utilize the IMFs' benefits, and avoid the negatives of the EMD, we introduce a new decomposition scheme termed variational mode decomposition (VMD) for GPR data processing for imaging. Based on the decomposition results of the VMD, we propose a new method which we refer as "the IMF-slice". In the proposed method, the IMFs are generated by the VMD trace by trace, and then each IMF is sorted and recorded into different profiles (i.e., the IMF-slices) according to its center frequency. Using IMF-slices, the GPR data can be divided into several IMF-slices, each of which delineates a main vibration mode, and some subsurface layers and geophysical events can be identified more clearly. The effectiveness of the proposed method is tested using synthetic benchmark signals, laboratory data and the field dataset.

Data Processing and Interpretation of Antarctic Ice-Penetrating Radar Based on Variational Mode Decomposition

“…This study uses the VMD method [22][23][24][25]. It is a signal decomposition estimation method based on classical Wiener filtering [26,27], Hilbert transform [28] and mixing variational problem.…”

The attenuated Ricker wavelet basis for seismic trace decomposition and attenuation analysis

“…; Zhang et al . ). This paper highlights the seismic decomposition based on a new type of wavelet basis and the extended applications.…”

“…As the first issue discussed above, inspired by the nonstationary convolution in time-frequency domain (Margrave et al 2002), we propose a new type of wavelet basis in this paper based on the classic Ricker wavelet, where the seismic attenuation is considered and referred to as the quality factor Q. For the second issue, we have combined the multi-channel decomposition technique and the idea of OMP to form an improved algorithm, which we referred to as the multi-channel OMP in our recent research work, and the examples demonstrated sparser decomposition results Zhang et al 2018). This paper highlights the seismic decomposition based on a new type of wavelet basis and the extended applications.…”

Dip Filter and Random Noise Suppression for GPR B-Scan Data Based on a Hybrid Method in f - x Domain

“…As part of the long-term program focused on the algorithms and applications of the GPR data time-frequency decomposition, this paper can be seen as a follow-up study of our previous work [18]. The research targets of this paper are to: (1) review the method of the f -x EMD and test its effectiveness on GPR data; (2) extend the original f -x EMD to the dip filter in f -x domain by reconstructing the frequency slices from the linear combination of different IMFs; (3) form the hybrid scheme by combining noise suppression based on SSA with sliding windows with the dip filter in f -x domain.…”