Sparse Gabor Transform and Its Application in Seismic Data Analysis

Considering the limitation of the Gabor transform due to the uncertainty principle, where time and frequency resolution cannot both be maximized simultaneously, we propose a post-processing strategy for the time-frequency spectrum to mitigate this limitation and improve time-frequency concentration. In the frequency band of seismic data, the time and frequency window sizes of the Gabor transform are fixed, implying that the Gabor transform’s time-frequency spectrum is formed by the 2D convolution of a high-resolution time-frequency spectrum with a Gaussian-shaped point spread function (PSF). Therefore, based on compressed sensing theory, we use sparse constraints to the time-frequency spectrum and solve the 2D deconvolution of the Gabor transform’s time-frequency spectrum using the alternating direction method of multipliers algorithm to eliminate the influence of the time-frequency window function. The PSF used for deconvolution is determined by the variances of the time and frequency windows, and by altering the size of the PSF, we can obtain frequency-sparse Gabor transform (FSGT) and time-sparse Gabor transform (TSGT). Simulation signals demonstrate the effectiveness of this post-processing strategy. For actual data, we prove that the sparse Gabor transform can enhance time-frequency concentration and improve the accuracy of seismic data analysis by integrating thin layer identification and frequency-dependent amplitude variation with offset attributes.