In prestack imaging, the first step is the extrapolation of the source and received wave fields into subsurface. The source and received wave fields are extrapolated recursively for each depth level as follows:

(46) |

(47) |

(48) |

(49) |

Therefore, one of the most important properties of the extrapolation operator is the unitary that makes equation () valid. Even though many algorithms based on the wave equation are accepted as unitary operators, each algorithm has a different property in terms of closeness to unitary. Among the many extrapolation algorithms which allows lateral velocity changes like finite-difference Claerbout (1985), PSPI Gazdag and Sguazzero (1984) and split-step Fourier Stoffa and Fokkema (1990). I chose the split-step Fourier method Stoffa and Fokkema (1990) for this thesis because of its pseudounitary property (Appendix A).

To show the closeness to unitary of the split-step Fourier extrapolation, I tested a plane wave extrapolation followed by its adjoint operator in the Marmousi velocity model. The input is a plane wave which is horizontally aligned spikes at time equals zero (Figure (a)). Figure (b) shows synthesis operator obtained by extrapolating a plane wave from depth level 800 m to surface. Figure (c) shows source wave field obtained by applying its adjoint operator to the synthesis operator which is extrapolating the synthesis operator at surface to depth level 800 m. Figure is amplitude of the source wave field (Figure (c)) along the plane wave.

Figure 3

Figure 4

2/5/2001