Lawrence Berkeley National Laboratory
Efficient discontinuous finite difference meshes for 3-D Laplace-Fourier domain seismic wavefield modelling in acoustic media with embedded boundaries
- Author(s): Alsalem, HJ
- Petrov, P
- Newman, G
- Um, E
- Rector, J
- et al.
Published Web Locationhttps://doi.org/10.1093/gji/ggz361
© 2019 The Author(s). Published by Oxford University Press on behalf of The Royal Astronomical Society. Simulation of acoustic wave propagation in the Laplace?Fourier (LF) domain, with a spatially uniform mesh, can be computationally demanding especially in areas with large velocity contrasts. To improve efficiency and convergence, we use 3-D second- and fourth-order velocitypressure finite difference (FD) discontinuous meshes (DM). Our DM algorithm can use any spatial discretization ratio between meshes. We evaluate direct and iterative parallel solvers for computational speed, memory requirements and convergence. Benchmarks in realistic 3-D models and topographies show more efficient and stable results for DM with direct solvers than uniform mesh results with iterative solvers.