Nonlinear seismic amplitude versus offset inversion using the exact Zoeppritz equation

TL;DR

This paper introduces an efficient seismic amplitude versus offset inversion method using an explicit analytical gradient of the exact Zoeppritz equation, improving accuracy and computational speed for elastic property estimation.

Contribution

It develops an adjoint-state-based analytical gradient formulation for the exact Zoeppritz equation, enhancing seismic inversion accuracy and efficiency over previous numerical approaches.

Findings

Effective in estimating elastic properties from synthetic and real data.

Provides stable and reliable property models with varying noise levels.

Improves computational efficiency by avoiding numerical gradient approximations.

Abstract

The amplitude-variation-with-offset inversion techniques are formulated to estimate elastic properties by fitting modeled seismic responses to observed data. Solving inverse seismic problems requires minimizing a target objective function for which gradient-based methods are frequently adopted. However, the efficiency and accuracy of these methods depend significantly on the approach used to compute the gradient of the target function. This work presents an explicit analytical gradient formulation of the exact Zoeppritz equation, discretized for multilayer media and derived using the adjoint-state method. The resulting expressions provide the gradient of a convolution-based objective function with respect to P-wave velocity, S-wave velocity, and density. The adjoint state-based solution improves computational efficiency by avoiding numerical approximations while maintaining high accuracy in calculating the gradient for seismic inversion. Additionally, using the exact Zoeppritz equation helps overcome the limitations associated with weak elastic property contrasts across subsurface layers. The least squares target function is minimized using a nonlinear limited-memory quasi-Newton algorithm. We demonstrate the effectiveness of the analytical gradient solution of the exact Zoeppritz equations in seismic inversion problems involving P-wave and S-wave velocity and density models. The inversion methodology is validated using 1D well logs-based and 2D synthetic seismic data with varying noise levels. Then it is applied to a 2D field data set from the Troll oil and gas field in the Norwegian North Sea. The results demonstrate that the proposed inversion framework provides stable and reliable estimates of elastic property models.