Efficient depth extrapolation of waves in elastic isotropic media |

We start with the wave equation governing the displacement in an arbitrary heterogenous isotropic elastic medium in the *Navier form* (see (Segall, 2010)):

where denote the components of a displacement field, is the shear modulus, is Poisson's ratio for the medium, and is the density. In this paper we consider a heterogenous elastic medium under the assumption of

where is the Lamé coefficient (see (Mavko et al., 2009),(Segall, 2010)), we get

where are horizontal wave numbers and is the frequency. The left-hand side of system 3 is the result of an ordinary differential operator applied to a vector-function and parameterized by horizontal wave numbers. In the present form equations 3 cannot be used for computationally efficient explicit depth extrapolation in a heterogeneous medium; however, these equations can be used for modeling displacements by solving a series of boundary-value problems (see (Maharramov, 2012)). In (Maharramov, 2012) it was suggested that equations 3 might be factorized in such a way as to allow solving them by alternating

where

(5) |

and are matrices with components that are complex-valued functions of the horizontal wave numbers, is the identity matrix. Performing the multiplication in equation 4 and using equation 3, we obtain:

where

Combining equations 6 and 7, we get the following equation for the operators and :

where

Equations 4, 8 in combination with equations 7 and 9 suggest the following procedure for extrapolating solutions to system 1 in depth:

- Solve the system of matrix equations 8 for , for each pair of horizontal wave numbers and two reference values of each elastic parameter and ;
- Evaluate
- Solve

by downward continuing in depth, using the formula

- Perform each step of the depth extrapolation for four combinations of the reference elastic parameters, then apply the inverse Fourier transform to the four fields and interpolate at each spatial point of the depth slice using true as e.g. in the PSPI method (see (Biondi, 2005)).
- After reaching the desired maximum depth, find the solution
by
*upward*extrapolation:

- Repeat the above steps for each frequency component

is not in the interior of the left half-plane, and the spectrum of

is not in the interior of the right half-plane. While the above algorithm tries to mimic two-way wave propagation, it is effectively just an approximation to the propagation process as it ignores the interaction between the up and down-going wave at intermediate depth steps. A less accurate alternative would be to downward-continue the wave field using equation 11 in a way similar to the one-way depth extrapolation using the scalar square-root equation (see (Claerbout, 1985),(Biondi, 2005)). The latter approach would be unable to image any dips beyond 90 , however, it would reduce the cost of extrapolation by a further factor of 2. Note the cost of solving equation 10 in depth is roughly three times that of solving the scalar square-root equation.

The above analysis may be extended to the case of an arbitrary anisotropic elastic medium. The fact that the components of the pseudo-differential operator matrices are not given in an analytical form, but are only computed numerically, does not limit their applicability.

Factorization of system 3 in the elastostatic case was one of the approaches mentioned by the author in Maharramov (2012). However, the one-way extrapolation technique is mostly useful for elastodynamic problems as the passband of the factorized depth extrapolators (e.g., as in equation 11) narrows down to zero with the temporal frequency passing to the zero static limit.

maximageigenval
The phase of a phase-shift operator corresponding to the maximum imaginary part of the eigenvalues of operator 15. Multicomponent ``phase-shift'' is defined by three such scalar phase-shift operators and a
matrix
of equation 16.
Figure 1. | |
---|---|

Note that equation 1 uses elastic parameterization that degenerates into a singularity if the shear modulus is equal to zero. This is not causing any problems with purely acoustic wave extrapolation as the singularity is effectively removed from equations 3 by the substitution in equation 2.

Efficient depth extrapolation of waves in elastic isotropic media |

2012-10-29