Affiliation:
1. Shell International E&P, Rijswijk, the Netherlands. .
Abstract
A preconditioned iterative solver for the 3D frequency-domain wave equation applied to seismic problems is evaluated. The preconditioner corresponds to an approximate inverse of a heavily damped wave equation deduced from the (undamped) wave equation. The approximate inverse is computed with one multigrid cycle. Numerical results show that the method is robust and that the number of iterations increases roughly linearly with frequency when the grid spacing is adapted to keep a constant number of discretization points per wavelength. To evaluate the relevance of this iterative solver, the number of floating-point operations required for two imaging problems are roughly evaluated. This rough estimate indicates that the time-domain migration approach is more than one order of magnitude faster. The full-wave-form tomography, based on a least-squares formulation and a scale separation approach, has the same complexity in both domains.
Publisher
Society of Exploration Geophysicists
Subject
Geochemistry and Petrology,Geophysics
Reference30 articles.
1. Multifrontal parallel distributed symmetric and unsymmetric solvers
2. Aminzadeh, F., J. Brac, and T. Kunz, 1997, 3-D salt and overthrust models: SEG/EAGE 3-D Modeling Series no.1, SEG.
3. Robust Ordering of Sparse Matrices using Multisection
4. A perfectly matched layer for the absorption of electromagnetic waves
5. Briggs, W. L. , V. E. Henson, and S. F. McCormick, 2000, A multigrid tutorial, 2nd ed.: SIAM.
Cited by
108 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献