The migration algorithm used for this project uses a variation of the split-step Fourier domain migration method, developed by Stoffa Stoffa et al. (1990). Simple Fourier domain migration based on the single square root operator (SSR) is also known as phase-shift migration Gazdag (1978). In phase-shift migration, a wavefield at the surface is downward-continued, assuming the wavefield is generated by the subsurface reflectors exploding Claerbout and Black (1997). The single square root operator marches the wavefield down, by shifting in the domain:

(2) |

The imaging principal sums over all , giving the wavefield at t=0, which is the migrated image Yilmaz (1987).
The split-step methods developed by Stoffa use an approximation of the SSR, where v is replaced by an
average velocity of a given depth step, *v*_{ref}. A new term is added to correct for the difference of the
average velocity and the actual medium velocity *v*_{m}. The equation looks like this variation of
equation (2) (*k*_{z} is shown):

(3) |

In equation (3), the first term is calculated in the domain, the second in the
domain. So at every depth step, a 2-D FFT transforms the wavefield into space. The
wavefield is downward-propagated using the phase-shift operator and v=*v*_{ref}. A 1-D IFFT transforms into
space, where the correction term is applied. Once again, summing over all frequency gives the
image. This method is exact if *v*_{ref}=*v*_{m}, or if layers are flat. If lateral velocity variations are
significant or units are dipping steeply, a better approximation of the SSR equation is needed. To handle
these conditions, we use the double square root equation (DSR). The DSR equation is a higher order approximation
of the SSR equation. The DSR equation as displayed in IEI Claerbout (1984) is

(4) |

where,

(5) |

(6) |

The DSR equation inherently includes both migration and NMO. Therefore, with input CMPs and a velocity function, a good image is created. The actual migration program used in this project is one developed by Biondo Biondi that uses a variation of the split-step method with the DSR equation Biondi (1998). The DSR equation can better handle lateral velocity variations as well as steeper dips (as compared to the SSR equation). The use of the split-step method with the DSR makes the algorithm more effective. In this way, multiple reference velocities can be used, and sharp lateral velocity variations, such as salt in contact with sediment, can be handled. For this project, three reference velocities were used. Biondi's program also has an input salt velocity, so that the reference velocities in the sediment will not be anomalously high if the salt is at the same depth step. So for each depth step, three wavefields were generated, then an interpolation in the space domain was used as needed.

Figure 5

Figure 6

4/27/2000