Seismic pre-stack AVA inversion scheme based on lithology constraints

Abstract

Seismic pre-stack AVA inversion using the Zoeppritz equation and its approximations as a forward engine yields P- and S-wave velocities and density. Due to the presence of seismic noise and other factors, the solution to seismic inversion is generally ill-posed and it is necessary to add constraints to regularize the algorithm.Moreover, since pre-stack inversion is a nonlinear problem, linearized optimization algorithms may fall into false local minima. The simulated annealing (SA) algorithm, on the other hand, is capable of finding the global optimal solution regardless of the initial model. However, when applied to multi-parameter pre-stack inversion, standard SA suffers from instability. Thus, a nonlinear pre-stack inversion method is proposed based on lithology constraints. Specifically, correlations among the elastic parameters are introduced to establish constraints based on a Bayesian framework, with special intention of mitigating the ill-posedness of the inversion problem as well as addressing the lithological characteristics of the formations. In particular, to improve the stability, a multivariate Gaussian distribution of elastic parameters is incorporated into the model updating the SA algorithm.We apply the algorithm to synthetic and field seismic data, indicating that the proposed method has a good resolution and stability performance.