Amplitude variation with incidence and azimuth stepwise inversion with coherence-attribute constraints for anisotropic parameters

With the development of 5D (3D + offset + azimuth) seismic technology, the stable acquisition of anisotropy information from wide-azimuth seismic data has become a key scientific problem in the seismic inversion of fractured reservoirs. The amplitude variation with incidence and azimuth (AVAZ) inversion method using wide-azimuth seismic data is an effective way to predict the anisotropic information of the subsurface medium. However, the conventional AVAZ inversion method suffers from too many parameters to be estimated, large variation in contribution, and inversion instability. Therefore, an AVAZ inversion method with coherence-attribute constraints is developed to solve the problem of unstable inversion of anisotropic parameters. First, we use seismic coherence attributes to build a fracture-probability-distribution model containing anisotropic information of the subsurface medium, which can be used to simulate large-scale subsurface fractures and faults. Then, it is added to the objective function as an anisotropic information constraint to improve the reliability and stability of the anisotropic inversion. Furthermore, an AVAZ inversion method in a Bayesian framework is implemented by using wide-azimuth seismic data. Gaussian distribution and a smoothing background model are added to the objective function to improve the reasonableness and stability of the inversion. In addition, we develop a stepwise optimization inversion method for isotropic and anisotropic parameters, prioritizing the inversion of parameters that contribute significantly to the reflection coefficient, and then using the results of the previous inversion as the initial values for the next inversion step to achieve multiparameter inversion. This method can reduce the number of the estimated parameters and thus improve the stability of the inversion of the anisotropic parameters. Field data examples indicate that this method produces suitable inversion results even at moderate levels of noise. Therefore, we can conclude that the proposed method has good applicability and stability in predicting the anisotropy parameters of fractured shale reservoirs.