Deconvolution denoising in the f-x domain has some defects when facing situations like complicated geology structure, coherent noise of steep dip angles, and uneven spatial sampling. To solve these problems, a new filtering method is proposed, which uses the generalized S transform which has good time-frequency concentration criterion to transform seismic data from the time-space to time-frequency-space domain （t-f-x）. Then in the t-f-x domain apply Empirical Mode Decomposition （EMD） on each frequency slice and clear the Intrinsic Mode Functions （IMFs） that noise dominates to suppress coherent and random noise. The model study shows that the high frequency component in the first IMF represents mainly noise, so clearing the first IMF can suppress noise. The EMD filtering method in the t-f-x domain after generalized S transform is equivalent to self-adaptive f-k filtering that depends on position, frequency, and truncation characteristics of high wave numbers. This filtering method takes local data time-frequency characteristic into consideration and is easy to perform. Compared with AR predictive filtering, the component that this method filters is highly localized and contains relatively fewer low wave numbers and the filter result does not show over-smoothing effects. Real data processing proves that the EMD filtering method in the t-f-x domain after generalized S transform can effectively suppress random and coherent noise of steep dips.
We designed the window function of the optimal Gabor transform based on the time-frequency rotation property of the fractional Fourier transform. Thus, we obtained the adaptive optimal Gabor transform in the fractional domain and improved the time-frequency concentration of the Gabor transform. The algorithm first searches for the optimal rotation factor, then performs the p-th FrFT of the signal and, finally, performs time and frequency analysis of the FrFT result. Finally, the algorithm rotates the plane in the fractional domain back to the normal time-frequency plane. This promotes the application of FrFT in the field of high-resolution reservoir prediction. Additionally, we proposed an adaptive search method for the optimal rotation factor using the Parseval principle in the fractional domain, which simplifies the algorithm. We carried out spectrum decomposition of the seismic signal, which showed that the instantaneous frequency slices obtained by the proposed algorithm are superior to the ones obtained by the traditional Gabor transform. The adaptive time frequency analysis is of great significance to seismic signal processing.
