A Practical Understanding of Pre- and Poststack Migrations
This volume, SEG Course Notes Series No. 13, is designed to give the practicing geophysicist an understanding of the principles of poststack migration, presented with intuitive reasoning rather than laborious math. Modeling is introduced as a natural process that starts with a geologic model and then builds seismic data. Migration is then described as the reverse process that uses seismic data to find the geologic model. Many other topics are covered relating to the quality of the migrated section, such as aliasing, rugged topography, or use of the correct velocity. Significant new material has been added in this revised edition of the original 1997 book, especially algorithms based on the phase-shift method, such as PSPI and the omegaX method.
Solving the wave equation
-
Published:January 01, 2007
Abstract
There are a number of ways to derive the wave equation, and many forms in which it is expressed. The approach taken in these notes will be based on the works by Claerbout [23], [294], Yilmaz [83], Stolt [21], and Brysk [100]. The main objective is to define the paths to the solutions used in seismic processing for modelling and migration. The main emphasis will be on the downward continuation method of seismic migration. Since this is a collection from the above authors, the notations used will reflect that used by the specific author to enable comparison to with their results. Consequently some parameters may change, specifically when moving from the solutions of the wave equations to the finite difference solutions.
It is assumed that the density is constant, and that the velocity varies in a vertical (z) and horizontal (x) manner. The 2-D seismic model is assumed to be 3-D with axis x, z, and time (t), with the pressure amplitude defined at any point in the volume as P(x, z, t). The zero offset or stacked section is defined by the surface P(x, z=0, t), and the geological cross-section or desired depth migration as P(x, z, t=0). The interval velocity v(x, z) is assumed to be isotropic (independent of direction). Note that some solutions use modified forms of the velocities such as RMS velocities (i.e. Kirchhoff time) and other use velocities defined in time as v(x, t).