Analysis of the stability and dispersion for a Riemannian acoustic wave equation



Título de la revista

ISSN de la revista

Título del volumen




The construction of images of the Earth's interior using methods as reverse time migration (RTM) or full wave inversion (FWI) strongly depends on the numerical solution of the wave equation. A mathematical expression of the numerical stability and dispersion for a particular wave equation used must be known in order to avoid unbounded numbers of amplitudes. In case of the acoustic wave equation, the Courant–Friedrich–Lewy (CFL) condition is a necessary but is not a sufficient condition for convergence. Thus, we need to search other types of expression for stability condition. In seismic wave problems, the generalized Riemannian wave equation is used to model their propagation in domains with curved meshes which is suitable for zones with rugged topography. However, only a heuristic version of stability condition was reported in the literature for this equation. We derived an expression for stability condition and numerical dispersion analysis for the Riemannian acoustic wave equation in a two-dimensional medium and analyzed its implications in terms of computational cost. © 2018 Elsevier Inc.


Palabras clave

Acoustic wave propagation, Acoustic waves, Acoustics, Dispersion (waves), Finite difference method, Numerical methods, Seismic prospecting, Wave equations, Computational costs, Finite difference scheme, Mathematical expressions, Neumann, Numerical dispersions, Numerical solution, Reverse time migrations, Stability condition, Stability criteria