Frequency-dependent AVO inversion applied to physically based models for seismic attenuation

Seismic inversion of amplitude versus offset (AVO) data in viscoelastic media can potentially provide high-resolution subsurface models of seismic velocities and attenuation from offset/angle seismic gathers. P- and S-wave quality factors (Q), whose inverse represent a measure of attenuation, depend on reservoir rock and pore fluid properties, in particular, saturation, permeability, porosity, fluid viscosity and lithology; however, these quality factors are rarely taken into account in seismic AVO inversion. For this reason, in this work, we aim to integrate quality factors derived from physically based models in AVO inversion by proposing a gradient descent optimization-based inversion technique to predict the unknown model properties (P- and S-wave velocities, the related quality factors and density). The proposed inversion minimizes the non-linear least-squares misfit with the observed data. The optimal solution is iteratively obtained by optimizing the data misfit using a second-order limited-memory quasi-Newton technique. The forward model is performed in the frequency-frequency-angle domain based on a convolution of broad-band signals and a linearized viscoelastic frequency-dependent AVO (FAVO) equation. The optimization includes the adjoint-state-based gradients with the Lagrangian formulation to improve the efficiency of the non-linear seismic FAVO inversion process. The inversion is tested on synthetic seismic data, in 1-D and 2-D, with and without noise. The sensitivity for seismic quality factors is evaluated using various rock physics models for seismic attenuation and quality factors. The results demonstrate that the proposed inversion method reliably retrieves the unknown elastic and an-elastic properties with good convergence and accuracy. The stability of the inverse solution especially seismic quality factors estimation relies on the noise level of the seismic data. We further investigate the uncertainty of the solution as a function of the variability of the initial models.