Iterative methods for 3D implicit finite-difference migration using the complex Padé approximation

Autor: Carlos A. N. da Costa, Francisco A. Silva Netoé, Jessé C. Costa, Itamara S. Campos, Jörg Schleicher, Amélia Novais
Rok vydání: 2013
Předmět:
Zdroj: Journal of Geophysics and Engineering. 10:045011
ISSN: 1742-2140
1742-2132
DOI: 10.1088/1742-2132/10/4/045011
Popis: Conventional implementations of 3D finite-difference (FD) migration use splitting techniques to accelerate performance and save computational cost. However, such techniques are plagued with numerical anisotropy that jeopardises the correct positioning of dipping reflectors in the directions not used for the operator splitting. We implement 3D downward continuation FD migration without splitting using a complex Pade approximation. In this way, the numerical anisotropy is eliminated at the expense of a computationally more intensive solution of a large-band linear system. We compare the performance of the iterative stabilized biconjugate gradient (BICGSTAB) and that of the multifrontal massively parallel direct solver (MUMPS). It turns out that the use of the complex Pade approximation not only stabilizes the solution, but also acts as an effective preconditioner for the BICGSTAB algorithm, reducing the number of iterations as compared to the implementation using the real Pade expansion. As a consequence, the iterative BICGSTAB method is more efficient than the direct MUMPS method when solving a single term in the Pade expansion. The results of both algorithms, here evaluated by computing the migration impulse response in the SEG/EAGE salt model, are of comparable quality.
Databáze: OpenAIRE