Interactive methods for 3-D finite-difference migration and modeling
, by Steve Cole
Splitting wave extrapolation operators introduces some inaccuracy in 3-D finite difference
modeling and migration. While this splitting can be avoided, traditional methods for doing
so are prohibitively expensive. Iterative methods such as Gauss-Seidel and successive
over-relaxation (SOR) are less expensive than direct solution methods and just as accurate.
An additional benefit is that iterative methods can be interrupted when a desired level of
accuracy is reached, making the methods even cheaper. A serious drawback is that the
iterative methods I use here have stability problems for the 45 degree approximation, which
is precisely where they are needed.
STANFORD