The 3-D Global Prestack Seismic Inversion in the Time–Frequency Mixed Domain

Abstract

Traditional seismic inversion has the problems of low lateral resolution and poor anti-noise performance. To address these challenges, the advantages of inversion methods in time and frequency domains are combined; the Lp norm that can preserve more sparsity information is introduced into seismic inversion. Meanwhile, multi-trace simultaneous inversion has become the mainstream of inversion methods, which can effectively improve the quality of pre-stack seismic inversion. However, traditional multi-trace simultaneous inversion methods have high computational costs and are difficult to obtain the 3D result. Ensuring spatial continuity of inversion results is a challenging task. To address this issue, a 3D global pre-stack inversion method in the time-frequency mixed domain is proposed. The proposed method adds a frequency domain fidelity constraint the objective function of traditional time domain inversion to improve the resolution of the inversion results. To improve the computational efficiency of multi-trace simultaneous inversion, the anisotropic total variation regularization method under Lp norm constraint is employed. In addition, the 2D Sylvester equation is extended to 3D space to enable 3D global inversion. The method not only improves computational efficiency, but also ensures spatial continuity of the results. The stability and feasibility of the proposed method have been demonstrated through its application on both the 3D overthrust model and field pre-stack seismic data.