Azimuthal amplitude difference inversion constrained by azimuth velocity anisotropy

Amplitude variation with incident angle and azimuth (AVAZ) inversion is usually used for estimating subsurface medium properties from azimuthal stacked seismic data for fracture prediction. And weak anisotropy parameters play an important role in the fracture prediction of shale reservoirs. However, the low contribution of the anisotropic parameters to the reflection coefficient and too many condition numbers of forward solver lead to insufficient inversion stability of the anisotropic parameters. Therefore, a new azimuthal-amplitude-difference inversion method with the azimuth-velocity-anisotropic constraints is proposed to improve the stability of inversion and obtain a reasonable estimation of anisotropic parameters. The azimuthal amplitude difference AVAZ forward solver in the anisotropic medium is initially established from a rewritten approximate reflectivity equation and variable substitution. Sensitivity analysis of estimated parameters of the rewritten approximate reflectivity indicates that our approach is sufficiently suitable for seismic inversion of anisotropic parameters. Furthermore, an azimuthal amplitude difference Bayesian AVAZ inversion approach with azimuth-velocity-anisotropic constraint is developed to improve the stability of the inversion. We add Gaussian noise to the synthetic seismic records and assume the Cauchy distributions as a prior constraint on the model parameters. The anisotropic information obtained from the azimuth-velocity-anisotropic inversion is utilized as a constraint for improving the stability of inversion. Synthetic data tests reveal that the anisotropic parameters can be estimated stably even with moderate noise. Field data tests illustrate the feasibility and reliability of the presented inversion method for estimating anisotropic parameters in a shale reservoir with vertical or near-vertical fractures.

Automated summary

AVAZ inversion is used to estimate subsurface medium properties from azimuthal stacked seismic data for fracture prediction. However, the low contribution of anisotropic parameters and the instability of inversion due to the condition numbers of the forward solver limit the accuracy of the results. Therefore, a new azimuthal-amplitude-difference inversion method with azimuth-velocity-anisotropic constraints is proposed to improve stability and obtain a reasonable estimation of anisotropic parameters.