# AVO equations

Series | Investigations in Geophysics |
---|---|

Author | Öz Yilmaz |

DOI | http://dx.doi.org/10.1190/1.9781560801580 |

ISBN | ISBN 978-1-56080-094-1 |

Store | SEG Online Store |

Consider the two elastic half-space layers in Figure 11.2-3e. The Zoeppritz equations (**14**) can be solved for the reflected and refracted *P*- and *S*-wave amplitudes, *A*_{1}, *B*_{1}, *A*_{2}, and *B*_{2}. However, our interest in exploration seismology is largely the angle-dependency of the *P*-to-*P* reflections given by the coefficient *A*_{1}. Specifically, we wish to infer or possibly estimate elastic parameters of reservoir rocks from reflection amplitudes and relate these parameters to reservoir fluids.

The exact expression for *A*_{1} derived from the solution of the Zoeppritz equations (**13**) is complicated and not intuitive in terms of its practical use for inferring petrophysical properties of reservoir rocks. The first approximation to the Zoeppritz equation for *P*-to-*P* reflection amplitude is given by Bortfeld ^{[1]} as

**(**)

**(**)

**(**)

The arrangement of the two terms on the right-hand side of equation (**14**) is based on separating the acoustic (the first term) and the elastic (the second term) effects on reflection amplitudes. As such, equation (**14**) does not explicitly indicate angle- or offset-dependence of reflection amplitudes; therefore, its practical implementation for AVO analysis has not been considered.

Instead, we shall use the approximation provided by Aki and Richards ^{[2]} as the starting point for deriving a series of practical AVO equations. Now that we only need to deal with the *P*-to-*P* reflection amplitude *A*_{1}, we shall switch to the conventional notation by replacing *A*_{1} with *R*(*θ*) as the angle-dependent reflection amplitude for AVO analysis.

By assuming that changes in elastic properties of rocks across the layer boundary are small and propagation angles are within the subscritical range, the exact expression for *R*(*θ*) given by the Zoeppritz equation can be approximated by ^{[2]}

**(**)

where *α* = (*α*_{1} + *α*_{2})/2, average *P*-wave velocity and Δ*α* = (*α*_{2} − *α*_{1}), *β* = (*β*_{1} + *β*_{2})/2, average *S*-wave velocity and Δ*β* = *β*_{2} − *β*_{1}, *ρ* = (*ρ*_{1} + *ρ*_{2})/2, average density and Δ*ρ* = *ρ*_{2} − *ρ*_{1}, and *θ* = (*φ*_{1} + *φ*_{2})/2, average of the incidence and transmission angles for the *P*-wave (Figure 11.0-2e).

Figure 11.2-8 shows the angle-dependent reflection amplitude associated with an interface with contrast in *P*- and *S*-wave velocities and densities based on the exact Zoeppritz equation, the Bortfeld approximation described by equation (**14**), and the Aki-Richards approximation described by equation (**15**) ^{[3]}. Note that these approximations closely follow the exact solution within the range of angles of incidence that are achievable by the recording of seismic data used in exploration. For very shallow reflectors, however, the approximate solutions would deviate from the exact solution significantly at very wide angles of incidence.

Note that the Aki-Richards approximation to the angle-dependent reflection amplitude *R*(*θ*) given by equation (**15**) has three parts in terms of Δ*α*/*α* which describes the fractional change in *P*-wave velocity across the layer boundary and hence may be referred to as the *P*-wave reflectivity, Δ*β*/*β* which describes the fractional change in the *S*-wave velocity across the layer boundary and hence may be referred to as the *S*-wave reflectivity, and Δ*ρ*/*ρ* which describes the fractional change in density across the layer boundary.

In practice, we do not observe the separate effects of *P*-wave reflectivity Δ*α*/*α, S*-wave reflectivity Δ*β*/*β* and fractional change in density Δ*ρ*/*ρ* on the reflection amplitudes *R*(*θ*). Instead, we observe changes in reflection amplitudes as a function of angle of incidence. In fact, it is the elastic parameters such as the *P*-wave reflectivity Δ*α*/*α, S*-wave reflectivity Δ*β*/*β* and fractional change in density Δ*ρ*/*ρ* that we wish to estimate from the *observed* angle-dependent reflection amplitudes. To use the Aki-Richards equation (**15**) in the inversion of reflection amplitudes for these elastic parameters, we first need to recast it in successive ranges of angle of incidence. This change of philosophy in arranging the terms in the Aki-Richards equation (**15**) was first introduced by Shuey ^{[4]} and led to practical developments in AVO analysis. The new arrangement in terms of successive ranges of angle of incidence is given by

**(**)

**Figure 11.2-3**Reflection and refraction of an incident*P*-wave at a layer boundary. Medium parameters:*ρ*is density,*α*is*P*-wave velocity,*β*is*S*-wave velocity. (a) Reflected*P*-wave; (b) reflected*S*-wave; (c) refracted*P*-wave; (d) refracted*S*-wave; (e) raypaths associated with the incident*P*wave, and reflected and refracted*P*- and*S*-waves. The radius of the circular wavefront associated with Huygens’ secondary source at*A*on the layer boundary is*CB*for the reflected wave, (*β*_{1}/*α*_{1}) for the reflected*S*-wave, (*α*_{2}/*α*_{1})*CB*for the refracted*P*-wave, and (*β*_{2}/*α*_{1}) for the refracted*S*-wave. The relationship between the angles in (e) is given by Snell’s law (equation**10**).

Another practical matter of concern is that the Aki-Richards equation (**15**) or any of its modifications that we shall derive in this section describe the *modeled* reflection amplitudes as a function *of angle of incidence*. However, the *observed* reflection amplitudes are available from CMP data as a *function of offset*. A need then arises either to transform the model equation for the reflection amplitudes from angle to offset coordinates ^{[6]} or to actually transform the CMP data from offset to angle coordinates. While the first approach is theoretically appealing, the practical schemes are based on the latter approach. We have already discussed such a transformation in the radon transform — the Radon transform using the linear moveout equation or its robust variation in the form of slant stacking. Figure 11.2-9 shows Zoeppritz amplitude curves as a function of angle of incidence and offset ^{[6]}.

Based on the theoretical conjecture made earlier by Koefoed ^{[7]} that the elastic property that is most directly related to angular dependence of reflection coefficient *R*(*θ*) is Poisson’s ratio *σ*, Shuey ^{[4]} introduced a variable transformation from *S*-wave velocity *β* to *σ*. The relationship between the two variables is given by equation (**L-49**) which we rewrite below as

**(**)

to perform the necessary differentiation

**(**)

The compressional-wave velocity *α* is given by equation (**3b**) and the shear-wave velocity *β* is given by equation (**L-47**) which is rewritten below as

**(**)

**(**)

where *μ* is Lamé’s constant.

We also define the *P*-wave reflection amplitude *R _{P}* at normal incidence as

**(**)

Substitute equations (**17a**), (**17b**), and (**18**) into the Aki-Richards equation (**16**) and perform some algebraic simplification to obtain

**(**)

Define a new term *H*

**(**)

and by way of equation (**18**) note that

**(**)

Next combine equations (**20**) and (**21a**) to obtain

**(**)

Finally, substitute equations (**21a**) and (**21b**) into the second term on the right-hand side of equation (**19**) to obtain

**(**)

where

**(**)

Equation (**22**) is known as Shuey’s three-term AVO equation. The first term *R _{P}* is the reflection amplitude at normal incidence. At intermediate angles (0 <

*θ*< 30 degrees), the third term may be dropped, thus leading to a two-term approximation

**(**)

where

**(**)

Equation (**24**) is known as Shuey’s two-term AVO equation. In practice, amplitudes picked along a moveout-corrected event on a CMP gather plotted against sin^{2} *θ* can be fitted to a straight line. The slope of the line gives the *AVO gradient attribute* and the ordinate at zero angle gives the *AVO intercept attribute*. The AVO gradient given by equation (**25**) is directly related to *change in Poisson’s ratio* Δ*σ*, which in turn, is directly related to *fluid saturation in reservoir rocks*. The AVO intercept attribute represents the reflectivity *R _{P}* at normal incidence. Therefore, the AVO intercept attribute, in lieu of conventional stack, can be used as input to derive the

*acoustic impedance attribute*(acoustic impedance estimation), which is indirectly related to

*porosity in reservoir rocks*.

Shown in Figure 11.2-10a is a portion of a section derived from 2-D prestack time migration. The objective is to identify fluid-saturated reservoir zones at the apex and the flanks of the structural closure. This image section is derived from the stacking of common-reflection-point (CRP) gathers associated with the prestack time-migrated data. The CRP gather in Figure 11.2-10b shows three events with amplitude variations with offset which are plotted in Figure 11.2-10c. By using Shuey’s equation (**24**), the AVO gradient and intercept sections are computed from the CRP gathers as shown in Figures 11.2-11 and 11.2-12, respectively. Note that the gradient section exhibits a group of AVO anomalies in the vicinity of the structural apex, possibly indicating fluid-saturated reservoir rocks.

At large angles of incidence beyond 30 degrees, the third term in equation (**22**) gradually becomes dominant. Note that this term is related directly to fractional change in *P*-wave velocity, Δ*α*/*α*. So, not only the reflection traveltimes at far offsets (normal moveout) corresponding to large angles of incidence, but also the reflection amplitudes at large angles of incidence make the biggest contribution to the resolution needed to estimate the changes in *P*-wave velocities.

The two-term equation (**24**) can be specialized for a specific value of Poisson’s ratio, *σ* = 1/3 and *H*_{0} = −1 so that equation (**25**) takes the form

**(**)

which can be solved for the change in Poisson’s ratio across a layer boundary

**(**)

This is the AVO attribute equation for estimating changes in Poisson’s ratio ^{[8]}. Actually, as described by equation (**27**), this attribute is the scaled sum of the AVO intercept *R _{P}* and AVO gradient

*G*attributes.

By recasting the first-order approximation to the Zoeppritz equation, Wiggins ^{[9]} and Spratt ^{[10]} derived a practical expression for *S*-wave reflectivity. First, drop the term with tan^{2} *θ* in equation (**16**) and rearrange the remaining terms to obtain

**(**)

such that, much like the definition for the *P*-wave reflectivity *R _{P}* given by equation (

**18**), an expression for the

*S*-wave reflectivity

*R*

_{S}

**(**)

can be explicitly inserted back into equation (**28**) to get

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle R(\theta)=R_P+\left(R_P-8\frac{\beta^2}{\alpha^2}R_S\right)\sin^2\theta + \left(2\frac{\beta^2}{\alpha^2}-\frac{1}{2}\right)\frac{\Delta\rho}{\rho}\sin^2\theta.}****(**)

Set *β*/*α* = 0.5 to make the last term on the right-hand side vanish and obtain ^{[10]}

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle R(\theta)=R_P+(R_P-2\,R_S)\sin^2\theta.}****(**)

This equation is of the form given by equation (**24**) where

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle G=R_P-2\,R_S,}****(**)

which can be rewritten explicitly in terms of the shear-wave reflectivity *R _{S}*

**(**)

This is the AVO attribute equation for estimating the shear-wave reflectivity. Given the AVO intercept *R _{P}* and AVO gradient

*G*attributes, simply take half of the difference between the two attributes to derive the shear-wave reflectivity

*R*as described by equation (

_{S}**33**).

Figure 11.2-13 shows the reflection amplitudes as a function of angle as predicted by equation (**31**) for three combinations of *R _{P}* and

*R*. Note that equation (

_{S}**31**) is a good approximation to the

*P*-to-

*P*reflection amplitudes as predicted by the exact Zoeppritz equation up to nearly 30 degrees of angles of incidence.

**Figure 11.2-11**The AVO gradient section defined by*G*in equation (**24**); (b) close-up of (a) in the vicinity of the structural apex.**Figure 11.2-12**(a) The AVO intercept section defined by*R*in equation (_{P}**24**); (b) close-up of (a) in the vicinity of the structural apex.**Figure 11.2-13***P*-to-*P*reflection amplitude as a function of angle of incidence computed by using the exact Zoeppritz equation (the three curves labeled as*A*) and the approximate form given by equation (**31**) (the three curves labeled as*B*)^{[10]}.

Return to the Aki-Richards equation (**15**) and consider the case of *N*-fold CMP data represented in the domain of angle of incidence. Note that the reflection amplitude *R*(*θ*) is a linear combination of three elastic parameters — *P*-wave reflectivity Δ*α*/*α, S*-wave reflectivity Δ*β*/*β* and fractional change in density Δ*ρ*/*ρ*.

Smith and Gidlow ^{[3]} argue in favor of solving for only two of the three parameters by making use of the empirical relation between density *ρ* and *P*-wave velocity *α* ^{[5]}:

**(**)

where *k* is a scalar. This relation holds for most water-saturated rocks. Differentiate to get

**(**)

Now substitute equation (**34b**) into the original Aki-Richards equation (**15**) and combine the terms with Δ*α*/*α*

**(**)

Simplify the algebra to obtain the desired two-parameter model equation for prestack amplitude inversion (Section L.6)

**(**)

Redefine the coefficients *a _{i}* and

*b*and write the discrete form of equation (

_{i}**36**)

**(**)

where

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle a_i=\frac{5}{8}-\frac{1}{2}\frac{\beta^2}{\alpha^2}\sin^2\theta_i+\frac{1}{2}\tan^2\theta_i}****(**)

and

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle b_i=-4\frac{\beta^2}{\alpha^2}\sin^2\theta_i.}****(**)

Consider the case of *N*-fold CMP data represented in the domain of angle of incidence. We have, for each CMP location and for a specific reflection event associated with a layer boundary, *N* equations of the form as equation (**37**) and two unknowns — Δ*α*/*α* and Δ*β*/*β*. We have more equations than unknowns; hence, we encounter yet another example of a generalized linear inversion (GLI) problem (Section J.1). The objective is to determine the two parameters such that the difference between the *modeled* reflection amplitudes *R _{i}* represented by equation (

**37**) and the

*observed*reflection amplitudes

*X*is minimum in the least-squares sense

_{i}^{[3]}.

**Figure 11.2-14**(a)*P*-wave and (b)*S*-wave reflectivity sections derived from the Smith-Gidlow inversion of prestack amplitudes (equation**39**).

The GLI solution in matrix form is given by (Section L.6)

**(**)

This matrix equation is solved for the two parameters — *P*-wave reflectivity Δ*α*/*α* and *S*-wave reflectivity Δ*β*/*β*. Note from the definitions of the coefficients *a _{i}* and

*b*given by equations (

_{i}**38a**,

**38b**) that you have to choose a value for the ratio

*β*/

*α*to compute the two reflectivities.

Shown in Figure 11.2-14 are the *P*- and *S*-wave reflectivity sections that were derived from prestack amplitude inversion applied to the CRP gathers associated with the data shown in Figure 11.2-10a. Note the AVO anomalies along the flanks of the structural closure. It is apparent that some of the AVO anomalies in the *P*-wave reflectivity section stand out more distinctively as compared with those in the *S*-wave reflectivity section. This may be related to the fact that changing from brine to gas causes a small change in *S*-wave reflectivity, while it causes a significant change in *P*-wave reflectivity. Such inference is supported by the laboratory measurements of *P*- and *S*-wave reflectivities at a gas-brine interface in sandstones ^{[10]}. Figure 11.2-15 shows *P*-wave reflectivities for a group of reservoir sandstones in which the pore fluid changes from brine to gas. For rocks with low induration (weak rocks), there is a large change in *P*-wave reflectivity, while for rocks with high induration the change becomes less significant (Figure 11.2-15a). The change in *S*-wave reflectivity is negligible for rocks with low and high induration (Figure 11.2-15b). A demonstrative field data example of *P*- and *S*-wave reflectivity contrasts that arise from a change in lithology and fluid saturation is shown in Figure 11.2-16.

Refer to equation (**17b**) and note that the difference between the *P*-wave and *S*-wave reflectivities is related to change in Poisson’s ratio Δ*σ* — a direct hydrocarbon indicator. In fact, Smith and Gidlow ^{[3]} have coined the term *pseudo-Poisson reflectivity* to describe the difference between the *P*-wave and *S*-wave reflectivities

**(**)

Note from equation (**17b**) that the pseudo-Poisson reflectivity is not exactly the same as what may be referred to as the proper Poisson reflectivity Δ*σ*/*σ*. If we define the ratio

**(**)

by differentiation, we can derive the difference relation given by equation (**40**) and thus show that

**(**)

Castagna ^{[11]} defined a straight line in the plane of *S*-wave velocity versus *P*-wave velocity as shown in Figure 11.2-17. This is called the *mudrock line* and is represented by the equation

**(**)

where the scalar coefficients *c*_{0} and *c*_{1} are empirically determined for various types of rocks. Suggested values for these scalars are *c*_{0} = 1360 and *c*_{1} = 1.16 for water-saturated clastics ^{[11]}. Gas-bearing sandstones lie above the mudrock line, and carbonates lie below the mudrock line as shown on a crossplot of *P*- and *S*-wave velocities sketched in Figure 11.2-17. The crossplot shown in Figure 11.2-18 is based on log measurements at a gas-producing well ^{[12]}. Note the separation of gas-sandstone cluster from the water-sandstone and shale clusters in the manner as sketched in Figure 11.2-17.

A way to quantify the prospectivity of the reservoir rock of interest is by defining a *fluid factor attribute* that indicates the position of the rock property with respect to the mudrock line. First, apply differentiation to both sides of equation (**42**) and note that

**(**)

Then, define the fluid factor Δ*F*

**(**)

If Δ*F* is close to zero, it means that you have water-saturated rock. If you have a gas-saturated sandstone, Δ*F* will be negative at the top and positive at the base of the reservoir unit ^{[3]}.

Figure 11.2-19 shows the pseudo-Poisson and fluid-factor sections that were derived from prestack amplitude inversion applied to the CRP gathers associated with the data shown in Figure 11.2-10a. Note the distinctive AVO anomalies along the flanks of the structural closure.

The two parameters Δ*α*/*α* and Δ*β*/*β*, estimated by using the least-squares solution given by equation (**39**), represent fractional changes in *P*- and *S*-wave velocities. As such, they are related to *P*- and *S*-wave reflectivities, Δ*I _{P}*/

*I*and Δ

_{P}*I*/

_{S}*I*, respectively, where

_{S}*I*and

_{P}*I*are the

_{S}*P*- and

*S*-wave impedances given by

**(**)

and

**(**)

From Section L.6, the *P*- and *S*-wave reflectivities are given by

**(**)

and

**(**)

Assuming that the density obeys Gardner’s relation given by equation (**34a**), the *P*-wave reflectivity Δ*I _{P}*/

*I*is simply the fractional change of the

_{P}*P*-wave velocity Δ

*α*/

*α*scaled by a constant, whereas the

*S*-wave reflectivity Δ

*I*/

_{S}*I*is a linear combination of the fractional changes of the

_{S}*P*- and

*S*-wave velocities, Δ

*α*/

*α*and Δ

*β*/

*β*, respectively.

A direct estimation of the *P*- and *S*-wave reflectivities given by equations (**L-46a**,**L-46b**) can be made by using an alternative formulation of prestack amplitude inversion which is based on transforming the Aki-Richards equation (**15**) to the new variables Δ*I _{P}*/

*I*and Δ

_{P}*I*/

_{S}*I*

_{S}^{[13]}. Solve equation (

**46a**) for Δ

*α*/

*α*and equation (

**46b**) for Δ

*β*/

*β*, and substitute into the Aki-Richards equation (

**15**). Following the algebraic simplification, we obtain (Section L.6)

**(**)

Goodway ^{[13]} have implemented a specific form of equation (**47**) to derive the AVO attributes Δ*I _{P}*/

*I*and Δ

_{P}*I*/

_{S}*I*. For a specific value of

_{S}*α*/

*β*= 2 and small angles of incidence for which tan

*θ*≈ sin

*θ*, the third term in equation (

**47**) vanishes. Compare equations (

**18**) and (

**46a**), and equations (

**29**) and (

**46b**), and note that equation (

**47**) with the remaining terms takes the form

**(**)

The resulting equation then is solved for the *P*- and *S*-wave reflectivities, Δ*I _{P}*/

*I*= 2

_{P}*R*and Δ

_{P}*I*/

_{S}*I*= 2

_{S}*R*, respectively. Again, consider the case of

_{S}*N*-fold CMP data represented in the domain of angle of incidence. In discrete form, equation (

**48**) can be rewritten as

**(**)

where *i* is the trace index and the coefficients *a _{i}* and

*b*are given by

_{i}

**(**)

and

**(**)

The reflection amplitude *R _{i}* in equation (

**49**) is a linear combination of the two parameters,

*R*and

_{P}*R*. As for the Smith-Gidlow equation (

_{S}**37**), the two parameters,

*R*and

_{P}*R*, can be estimated for a specific event from CMP data using the least-squares solution given by

_{S}

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{pmatrix} \sum\nolimits_i^N\,a^2_i&\sum\nolimits_i^N\,a_ib_i\\ \sum\nolimits_i^N\,a_ib_i&\sum\nolimits_i^N\,b^2_i \end{pmatrix} \begin{pmatrix} R_P\\ R_S \end{pmatrix}= \begin{pmatrix} \sum\nolimits_i^N\,a_iX_i\\ \sum\nolimits_i^N\,b_iX_i \end{pmatrix}.}****(**)

Following the estimation of the *P*-wave reflectivity *R _{P}* and the

*S*-wave reflectivity

*R*using the least-squares solution given by equation (

_{S}**51**), the

*P*-wave impedance

*I*and the

_{P}*S*-wave impedance

*I*can be computed by integration.

_{S}By using the impedance attributes *I _{P}* and

*I*, Goodway

_{S}^{[13]}compute two additional AVO attributes in terms of Lamé’s constants scaled by density —

*λρ*and

*μρ*. Substitute equation (

**3b**) into (

**45a**) to get the relation

**(**)

and substitute equation (**17c**) into equation (**45b**) to get the relation

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mu\rho=I^2_S.}****(**)

Note from equations (**52a**,**52b**) that

**(**)

Figure 11.2-20 shows the *μρ* and *λρ* AVO attribute sections which were derived from prestack amplitude inversion applied to the CRP gathers associated with the data shown in Figure 11.2-10a. As in the case of the pseudo-Poisson and fluid-factor AVO attribute sections shown in Figure 11.2-19, note the distinctive AVO anomalies along the flanks of the structural closure. Goodway ^{[13]} convincingly demonstrates that separation of gas-bearing sands from tight sands and carbonates is much better with the crossplot of the Lamé attributes (*μρ, λρ*) in contrast with the crossplot of the *P*- and *S*-wave reflectivities or impedances (Figure 11.2-21).

Figure 11.2-22 outlines a summary of the AVO equations that are described in this section based on the various approximations to the Aki-Richards equation (**15**). Start with Shuey’s arrangement of the terms in the Aki-Richards equation, which itself is an approximation to the exact expression for the *P*-to-*P* reflection amplitudes given by the Zoeppritz equations. Then, apply a transformation from *S*-wave velocity to Poisson’s ratio, and drop the third term to get a simple expression for the *P*-to-*P* reflection amplitude that is a linear function of sin^{2} *θ*. This linear approximation yields the AVO intercept and gradient attributes.

Alternatively, you may drop the third term in the original Shuey arrangement given by equation (**16**) and apply Wiggins’ rearrangement of the terms. Then, assume the specific case of the *P*-wave velocity to be twice the *S*-wave velocity to obtain the equation that yields, once again, the AVO intercept and gradient attributes. The Wiggins AVO intercept and gradient attributes, unlike the Shuey counterparts, directly lead to deriving the *P*-wave and *S*-wave reflectivity attributes, albeit only under the assumption that the *P*-wave velocity is twice the *S*-wave velocity.

Returning to the three-term Aki-Richards equation (**15**), you may reduce it to the two-term Smith-Gidlow equation (**36**) by using Gardner’s relation of density and *P*-wave velocity. Then, for a specified ratio of *P*- to *S*-wave velocity, perform inversion of prestack amplitudes associated with target events on CMP or CRP gathers to obtain the *P*- and *S*-wave reflectivity attributes. In addition, using these attributes to derive two more AVO attributes — psuedo-Poisson and fluid factor.

Finally, you may recast the Aki-Richards equation (**15**) in terms of *P*- and *S*-wave reflectivities (equations **46a**,**46b**) and assume that *α*/*β* = 2 to obtain the two-term Goodway et al. equation (**48**). Again, perform inversion of prestack amplitudes associated with target events on CMP or CRP gathers to obtain the *P*- and *S*-wave impedance attributes. Subsequently, by using the relations between impedances and Lamé’s constants, you can compute two additional AVO attributes — *λρ* and *μρ*.

The AVO equations derived in this section are all expressed in terms of angle of incidence (equations **24**, **31**, and **36**). In practice, the mapping of amplitudes associated with a reflection event on a CMP gather from offset to angle of incidence needs to be performed. Assume that the CMP gather is associated with a horizontally layered earth model so that the event moveout on the CMP gather is given by the hyperbolic equation

**(**)

where *t* is the two-way traveltime from the source to the flat reflecting interface back to the receiver, *t*_{0} is the two-way zero-offset time, *x* is the offset, and *v _{rms}* is the rms velocity down to the reflector. The ray parameter

*p*is given by the stepout

*dt*/

*dx*measured along the hyperbolic moveout trajectory (the slant-stack transform). Apply differentiation to equation (

**53**) to get

**(**)

Since the ray parameter *p* is also expressed as the ratio of sin *θ*/*v _{int}*, where

*v*is the interval velocity above the reflecting interface, it follows that

_{int}

**(**)

By using equation (**55**), reflection amplitudes on a CMP gather can be transformed from offset to angle domain, and subsequently used in the AVO equations (**24**), (**31**), (**36**), and (**48**).

## References

- ↑ Bortfeld (1961), Bortfeld, R., 1961, Approximation to the reflection and transmission coefficients of plane longitudinal and transverse waves: Geophys. Prosp., 9, 485–502.
- ↑
^{2.0}^{2.1}Aki and Richards (1980), Aki, K. I. and Richards, P. G., 1980, Quantitative seismology: W. H. Freeman and Co. - ↑
^{3.0}^{3.1}^{3.2}^{3.3}^{3.4}^{3.5}Smith and Gidlow, 1987, Smith, G. C. and Gidlow, P. M., 1987, Weighted stacking for rock property estimation and detection of gas: Geophys. Prosp., 35, 993–1014. - ↑
^{4.0}^{4.1}Shuey (1985), Shuey, R. T., 1985, A simplification of the Zoeppritz equations: Geophysics, 50, 609–614. - ↑
^{5.0}^{5.1}Gardner et al. (1974), Gardner, G. H. F., Gardner, L. W., and Gregory, A. R., 1974, Formation velocity and density — The diagnostic basis for stratigraphic traps: Geophysics, 39, 770–780. - ↑
^{6.0}^{6.1}^{6.2}Demirbag and Coruh, 1988, Demirbag, E. and Coruh, C., 1988, Inversion of Zoeppritz equations and their approximations: 58th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1199–1203. - ↑ Koefoed (1955), Koefoed, O., 1955, On the effect of Poission’s ratio of rock strata on the reflection coefficients of plane waves: Geophys. Prosp., 3, 381–387.
- ↑ Hilterman, 1983, Hilterman, F. J., 1983, Seismic lithology: presented as a continuing education course at the 53rd Ann. Internat. Mtg., Soc. Expl. Geophys.
- ↑ Wiggins et al. (1984), Wiggins, R., Kenny, G. S., and McClure, C. D., 1984, A method for determining and displaying the shear-wave reflectivities of a geologic formation: European patent application (Mobil Oil Corporation).
- ↑
^{10.0}^{10.1}^{10.2}^{10.3}^{10.4}^{10.5}Spratt et al. (1984), Spratt, R. S., Goins, N. R., and Fitch, T. J., 1984, Pseudo-Shear — The analysis of AVO:*in*Offset-dependent reflectivity — theory and practice: Soc. Expl. Geophys. - ↑
^{11.0}^{11.1}^{11.2}Castagna et al., 1985, Castagna, J. P., Bazle, M. L., and Eastwood, R. L., 1985, Relationship between compressional-wave and shear-wave velocities in elastic silicate rocks: Geophysics, 50, 571–581. - ↑
^{12.0}^{12.1}^{12.2}Fatti et al., 1994, Fatti, J. L., Smith, G. C., Vail, P. J., Strauss, P. J., and Levitt, P. R., 1994. Detection of gas in sandstone reservoirs using AVO analysis: A 3-D seismic case history using the Geostack technique: Geophysics, 59, 1362–1376. - ↑
^{13.0}^{13.1}^{13.2}^{13.3}Goodway et al., 1998, Goodway, B., Chen, T., and Downton, J., 1998, AVO and prestack inversion: Presented Ann. Mtg. Can. Soc. of Expl. Geophys.

## See also

- Analysis of amplitude variation with offset
- Reflection and refraction
- Reflector curvature
- Processing sequence for AVO analysis
- Derivation of AVO attributes by prestack amplitude inversion
- Interpretation of AVO attributes
- 3-D AVO analysis