Anisotropic total variation pre-stack multitrace inversion based on Lp norm constraint

Due to the instability of extracting elastic parameters from conventional elastic impedance variation with incident angle (EVA) pre-stack inversion and the fact that single-trace inversion does not consider the correlation between seismic traces and the lateral continuity of inversion results. An anisotropic total variation multitrace inversion method based on the $L_P$ norm constraint is proposed on the basis of EVA inversion. In this paper, the $L_P$ norm is used as the regularization constraint, the longitudinal and transverse difference operators are introduced, and the alternating direction multiplier algorithm and the 2D Sylvester equation solver algorithm are used as inversion algorithms to achieve a multitrace inversion of the elastic parameters. The feasibility of the method proposed in this paper is verified by making full use of the advantages of both methods and conducting trial calculations on some blocks of the Marmousi2 model. The method was further validated by applying it to actual data from marine shales of the Longmaxi-Wufeng Formation in the southern Sichuan Basin. The inversion results were validated using a combination of brittleness indices and microseismic monitoring to confirm the feasibility of the method in this paper.