Next: PROBLEMS WITH THE INTEGRATION Up: Stationary Phase Approach To Previous: INTRODUCTION

# BLEISTEIN'S DERIVATION

Bleistein (1987) proposed an inversion operator to process the seismic data:

 (6.1)

where is a pair of parameters that determines the source and receiver coordinates:

The domain of integration is the set of values which are required to cover the source-receiver array. is regarded as a smoothly tapered version of the source wavelet. The functions and and are the solutions to the ray-theoretic eikonal equation and transport equation with the source at and observation point at , and is the seismic data. The term is the determinant, i.e.,

 (6.2)

The weight function is chosen to give the true ampltude image:

 (6.3)

The integration in equation 6.1 is carried out over the entire data space, that is, the set of that is required to cover the entire source-receiver array. This inversion operator is a Kirchhoff-type integration. Bleistein applied this operator to the upward reflected field which is given by the Kirchhoff approximation :

 (6.4)

where is the angular-dependent reflection coefficient and is the upward normal unit vector. Equation 6.4 was plugged into equation 6.1 to yield

 (6.5)

where

 (6.6)

Then Bleistein evaluated equation 6.5 at the stationary phase point which yields

 (6.7)

The condition for the stationary phase point is

 = = 0 , (6.8)

and

 = = 0, (6.9)

where is a pair of parameters which determines the reflector coordinates

 (6.10)

The integration in equation 6.5 reduces to equation 6.7 which is evaluated at the stationary point. The evaluation at the stationary point gives the following result:

 (6.11)

where is a band-limited delta function with its support on the reflector. This equation shows that the output image of inversion operator 6.1 is a singular function scaled by the geometric-optics angular-dependent reflection coefficients. The term reaches its peak amplitude at the reflector, and the peak amplitudes are proportional to the reflection coefficients.

Note that in fact it is the stationary phase approximation (i.e. equation 6.7) to the diffraction stack formula given in equation 6.1 that leads to the true amplitude image given in equation 6.11. In other words, equation 6.7 is more close to equation 6.11 than equation 6.1. Once we know the stationary point, equation 6.7 directly yields a result described by equation 6.11.

Next: PROBLEMS WITH THE INTEGRATION Up: Stationary Phase Approach To Previous: INTRODUCTION
Caldera OpenLinux User
1999-02-09