A High-Resolution Spectral Analysis Method Based on Fast Iterative Least Squares Constraints

The prediction of reservoir and caprock thickness is important in geological evaluations for site selection for aquifer underground gas storage. Therefore, high-resolution seismic identification of reservoirs and caprocks is crucial. High-resolution time–frequency decomposition is one of the key methods for identifying sedimentary layers. Based on this, we propose a least squares constrained spectral analysis method using a greedy fast shrinkage algorithm. This method replaces the traditional Tikhonov regularization objective function with an L1-norm regularized objective function and employs a greedy fast shrinkage algorithm. By utilizing shorter window lengths to segment the data into more precise series, the method significantly improves the computational efficiency of spectral analysis while also enhancing its accuracy to a certain extent. Numerical models demonstrate that compared to the time–frequency spectra obtained using traditional methods such as wavelet transform, short-time Fourier transform, and generalized S-transform, the proposed method can achieve high-resolution extraction of the dominant frequencies of seismic waves, with superior noise resistance. Furthermore, its application in a research area in southern China shows that the method can effectively predict thicker sedimentary layers in low-frequency ranges and accurately identify thinner sedimentary layers in high-frequency ranges.