A 3D Helmholtz iterative solver with a semi-analytical preconditioner for acoustic wavefield modeling in seismic exploration problems
Numerical methods and programming, Tome 15 (2014) no. 3, pp. 514-529
Cet article a éte moissonné depuis la source Math-Net.Ru
An approach to the iterative solution of the 3D acoustic wave equation in a frequency domain is proposed, substantiated, and verified numerically. Our method is based on Krylov-type linear solvers, similarly to several other iterative solver approaches. The distinctive feature of our method is the use of a right preconditioner obtained as the solution of the complex dumped Helmholtz equation in a 1D medium, where velocities vary only with depth. The actual Helmholtz operator is represented as a perturbation of the preconditioner. As a result, a matrix-by-vector multiplication of the preconditioned system can be efficiently evaluated via 2D FFT in x and y directions followed by the solution of a number of ordinary differential equations in z directions. While solving ODE's. it is possible to treat the 1D velocity function as a piecewise constant one and to search for the exact solution as a superposition of upgoing and downgoing waves. This approach allows one not to use explicit finite-difference approximations of derivatives at all. The method has excellent dispersion properties in both lateral and vertical directions.
Keywords:
Helmholtz equation, iterative methods, preconditioners, acoustic waves, seismic exploration.
@article{VMP_2014_15_3_a11,
author = {D. A. Neklyudov and I. Yu. Silvestrov and V. A. Tcheverda},
title = {A {3D} {Helmholtz} iterative solver with a semi-analytical preconditioner for acoustic wavefield modeling in seismic exploration problems},
journal = {Numerical methods and programming},
pages = {514--529},
year = {2014},
volume = {15},
number = {3},
language = {ru},
url = {http://geodesic.mathdoc.fr/item/VMP_2014_15_3_a11/}
}
TY - JOUR AU - D. A. Neklyudov AU - I. Yu. Silvestrov AU - V. A. Tcheverda TI - A 3D Helmholtz iterative solver with a semi-analytical preconditioner for acoustic wavefield modeling in seismic exploration problems JO - Numerical methods and programming PY - 2014 SP - 514 EP - 529 VL - 15 IS - 3 UR - http://geodesic.mathdoc.fr/item/VMP_2014_15_3_a11/ LA - ru ID - VMP_2014_15_3_a11 ER -
%0 Journal Article %A D. A. Neklyudov %A I. Yu. Silvestrov %A V. A. Tcheverda %T A 3D Helmholtz iterative solver with a semi-analytical preconditioner for acoustic wavefield modeling in seismic exploration problems %J Numerical methods and programming %D 2014 %P 514-529 %V 15 %N 3 %U http://geodesic.mathdoc.fr/item/VMP_2014_15_3_a11/ %G ru %F VMP_2014_15_3_a11
D. A. Neklyudov; I. Yu. Silvestrov; V. A. Tcheverda. A 3D Helmholtz iterative solver with a semi-analytical preconditioner for acoustic wavefield modeling in seismic exploration problems. Numerical methods and programming, Tome 15 (2014) no. 3, pp. 514-529. http://geodesic.mathdoc.fr/item/VMP_2014_15_3_a11/