Statistical Synchrosqueezing Transform and Its Application to Seismic Thin Interbed Analysis

Synchrosqueezing transform (SST) benefits from an instantaneous frequency (IF) estimator in the time-frequency domain, providing an energy-concentrated time-frequency representation to describe the time-varying frequency of seismic signals. To enhance the concentration performance of SST, this paper theoretically proposes a statistical SST (SSST) by constructing a spectrum-weighted IF estimator in the short-time Fourier transform (STFT) domain. In this SSST, a linear chirp signal is introduced to better capture its chirp rate, and then the quadratic STFT spectra are weighted by a window-related weighting function. On this basis, two equations related to the IF and chirp rate are constructed to derive the spectrum-weighted IF estimator. Finally, the STFT coefficients are squeezed to the estimated IF trajectories by a frequency fixed-point iterative algorithm, thereby providing a more concentrated time-frequency representation for multi-component signals than existing advanced methods, while enabling to retrieve each component. Two synthetic examples and one field seismic data on thin interbeds are utilized to demonstrate the effectiveness of the proposed SSST and show its ability to highlight the time-varying frequency features of seismic signals, which is a promising seismic data analysis tool, such as characterizing thin interbed thickness variations.