Aki–Richards equation

From SubSurfWiki
Jump to: navigation, search

The Aki–Richards equation is an important linear approximation of the Zoeppritz equations. It is valid for reflection angles up to about 40°. An even simpler approximation is the Shuey approximation.

The reflection coefficient for plane elastic waves varies with incidence angle, as described by the Zoeppritz equations. These equations are very complicated to solve, so Aki & Richards (1980[1]) gave linear approximations to them. Various authors have since derived slightly different versions of these approximations. The one you use might depend on the required accuracy, or what the software you're using happens to have implemented.

Avseth's formulation

Adapted from Avseth et al (2006)[2]:

 R(\theta) = W - X \sin^2 \theta + Y \frac{1}{\cos^2 \theta_\mathrm{avg}} - Z \sin^2 \theta

where

 W = \frac{1}{2} \frac{\Delta \rho}{\rho}
 X = 2 \frac{V^2_\mathrm{S}}{V^2_\mathrm{P1}} \frac{\Delta \rho}{\rho}
 Y = \frac{1}{2} \frac{\Delta V^2_\mathrm{P}}{V^2_\mathrm{P}}
 Z = 4 \frac{V^2_\mathrm{S}}{V^2_\mathrm{P1}} \frac{\Delta V^2_\mathrm{S}}{V^2_\mathrm{S}}

where the 'delta' expressions are, for example:

 \frac{\Delta \rho}{\rho} = \frac{\rho_2 - \rho_1}{(\rho_1 + \rho_2)/2}

where ρ1 is the density of the upper layer, and ρ2 is the same property for the lower layer, and VP1 is the P-wave velocity of the upper layer.

Notice that the θ in the third term is the mean of the incident and transmission angles. Sometimes this is approximated by the incident angle.

Alternate formulation

For small contrasts and small angles, the approximations can be given as:[citation needed]

 R_\mathrm{PP}(\theta) = \frac{1}{2} \left( \frac{\Delta V_\mathrm{P}}{V_\mathrm{P}} + \frac{\Delta \rho }{\rho} \right) + \frac{1}{2} \left[ \frac{\Delta V_\mathrm{P}}{V_\mathrm{P}} - 4 \gamma^2 \left( \frac{\Delta \rho}{\rho} + 2 \frac{\Delta V_\mathrm{S}}{V_\mathrm{S}} \right) \right] \theta^2
 R_\mathrm{PS} (\theta) = -\frac{1}{2} \left( \frac{\Delta \rho}{\rho} \right) + 2 \gamma \left[ \frac{\Delta \rho}{\rho} + 2 \frac{\Delta V_\mathrm{S}}{V_\mathrm{S}} \right] \theta

where θ is the PP incident angle in units of radians, and

V_\mathrm{P} = \frac{1}{2}[ V_\mathrm{P1} + V_\mathrm{P2}]
V_\mathrm{S} = \frac{1}{2}[ V_\mathrm{S1} + V_\mathrm{S2}]
\rho = \frac{1}{2}[ \rho_1 + \rho_2]
\gamma = \frac{V_\mathrm{S}}{ V_\mathrm{P}}. Values for V_\mathrm{S} and V_\mathrm{P} are generally taken as the average value over the region of interest.
\Delta V_\mathrm{P} = V_\mathrm{P2} - V_\mathrm{P1}\
\Delta V_\mathrm{S} = V_\mathrm{S2} - V_\mathrm{S1}\
\Delta \rho = \rho_2 - \rho_1\

Implementation in Python

Some of these approximations have been implemented in Agile's bruges Python library.

Avseth's formulation is called bruges.reflection.akirichards.

See also

References

  1. Aki, K, and PG Richards (1980). Quantitative Seismology: Theory and Methods. WH Freeman and Co.
  2. Avseth, P, T Mukerji, G Mavko (2006), Quantitative Seismic Interpretation, Cambridge University Press.