Next: Discrete matrix formulation of Up: Wave-equation datuming operators Previous: Wave-equation datuming operators

# Kirchhoff datuming

Based on Green's second identity, the integral expression for a pressure wavefield within a source free volume bounded by surface S is

G is a Green's function satisfying the scalar wave equation, and is the inward pointing normal to the bounding surface (Figure ). For seismic data, the integral can be split into two parts so that S = S1 + S2 and the wavefield can be considered as arising from causal sources below S1, which represents the data collection surface. Since seismic data are recorded for a finite time, the surface S2 can be chosen in such a way that the contribution from this part of the integral is negligible. The choice of Green's function that vanishes on surface S1 reduces the integral to

When the undulations of the surface are small over a wavelength, the Green's function can be approximated by

where and . The subscript s indicates that the coordinate is on surface S1. is angular frequency and v is propagation velocity.

 greensthm Figure 1 Geometry for derivation of the Kirchhoff integral (after Berkhout, 1982).

The surface integral then becomes

The integral can be split into near-field and far-field contributions by writing
 (1)
where represents the traveltime from the surface S1 to the observation point .This equation can be recast in the time domain by using the relationships

to obtain
 (2)
Equations () and () are the 3-D frequency and time domain equations for the upward continuation of a scalar wavefield. For computational efficiency it is common to make the far-field approximation , and to drop the second term in the equations.

If the wavefield is invariant in one space coordinate, the frequency domain equation can be rewritten as

Integrating out the y dependence by the method of stationary phase yields
 (3)
Making the far-field approximation yields
 (4)
The equivalent time domain expression is
 (5)
For the 2-D versions of the integrals, .The term is a convolutional operator corresponding to frequency domain multiplication by . Rather than making the far-field approximation, Berryhill (1979) obtains a 2-D version in the time domain by integrating out the y dependence of equation (). His expression is similar to equation (), except that his convolutional operator depends on r. Berryhill's time domain formulation is accurate since it includes the near-field term, but in practice, the added computation of the delay-time variant convolution is often unnecessary.

Equations () through () are the Kirchhoff integral equations for upward continuation. Wave-equation datuming is performed by implementing these equations. Downward continuation is performed by means of the adjoint processes. The concept of adjoint operators will be more fully described in a subsequent section, but for now, it is only important to note that the adjoint forms of the integral equations will have the opposite sign in or in the time delay factor . The adjoint forms of the time domain integrals, along with the imaging condition, result in Schneider's (1978) Kirchhoff migration operators. Wiggins (1984) uses the far-field approximation in his implementation of Kirchhoff datuming and migration.

Next: Discrete matrix formulation of Up: Wave-equation datuming operators Previous: Wave-equation datuming operators
Stanford Exploration Project
2/12/2001