Seismic signal denoising using thresholded variational mode decomposition

Noise reduction is important for signal analysis. In this paper, we propose a hybrid denoising method based on thresholding and data-driven signal decomposition. The principle of this method is to reconstruct the signal with previously thresholded intrinsic mode functions (IMFs). Empirical mode decomposition (EMD) based methods decompose a signal into a sum of oscillatory components, while variational mode decomposition (VMD) generates an ensemble of modes with their respective centre frequencies, which enables VMD to further decrease redundant modes and keep less residual noise in the modes. To illustrate its superiority, we compare VMD with EMD as well as its derivations, such as ensemble EMD (EEMD), complete EEMD (CEEMD), improved CEEMD (ICEEMD) using synthetic signals and field seismic traces. Compared with EMD and its derivations, VMD has a solid mathematical foundation and is less sensitive to noise, while both make it more suitable for non-stationary seismic signal decomposition. The determination of mode number is key for successful denoising. We develop an empirical equation, which is based on the detrended fluctuation analysis (DFA), to adaptively determine the number of IMFs for signal reconstruction. Then, a scaling exponent obtained by DFA is used as a threshold to distinguish random noise and signal between IMFs and the reconstruction residual. The proposed thresholded VMD denoising method shows excellent performance on both synthetic and field data applications.