Dispersion analysis of an optimally accurate 3-D finite difference scheme for the elastic case

Abstract

The accuracy, efficiency and stability of numerical schemes for computing synthetic seismograms are important factors when the schemes are used for forward modeling and waveform inversion. We develop a method for evaluating the dispersion error directly for each volume element of the numerical grid using the criterion for optimal accuracy (eq. 2.20 of Geller & Takeuchi, 1995, GJI). This allows us to also conduct dispersion analyses of the OPT scheme for cylindrical coordinates. We conduct dispersion analysis and an evaluation of grid anisotropy for a stable and optimally accurate finite difference scheme (referred to below as "OPT") for a Cartesian grid (Mizutani et al., 2008). For 15 grid points per wavelength in a homogeneous medium, the maximum value of the dispersion error of the phase velocity for the OPT scheme implemented by a predictor-corrector (P-C) algorithm is about 0.02 %, which is 30 times less than that for a scheme which uses only the predictor step (e.g. conventional time marching). We find that the grid anisotropy for the P-C scheme is almost the same order as that of the predictor-only scheme, but the absolute value of the dispersion error is, as noted above, about 30 times smaller for the P-C scheme.