Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Implementation of the Kirchhoff integral for elastic waves in staggered-grid modeling schemes

Journal Article · · Geophysics; (United States)
DOI:https://doi.org/10.1190/1.1443576· OSTI ID:6684026
 [1]
  1. IKU Petroleum Research, Trondheim (Norway)

Implementation of boundary conditions in finite-difference schemes is not straightforward for the elastic wave equation if a staggered grid formulation is used. Reverse time migration of VSP data requires a proper description of the recording surface so as not to excite false P- and S-waves. Such contributions may cause artifacts in the imaging procedure. The boundary conditions for the elastic stress tensor can be implemented numerically in a staggered coarse grid modeling scheme by using band-limited spatial delta-functions and band-limited first-order derivatives of these spatial delta-functions. A representation theorem for elastic waves is derived to test the implementation of the spatial part of the boundary condition. The implementation is tested in a 2-D numerical experiment for a closed, but curved, boundary S enclosing a volume V. The test condition is that within the volume V, the difference between the forward modeled field and the retropropagated field should be equal to zero. Both P- and S-waves are properly recovered in a 2-D reverse time modeling example. The numerical artifacts related to the proposed spatial approximation of the boundary condition are found to be negligible.

OSTI ID:
6684026
Journal Information:
Geophysics; (United States), Journal Name: Geophysics; (United States) Vol. 59:12; ISSN GPYSA7; ISSN 0016-8033
Country of Publication:
United States
Language:
English