Fourth-Order Staggered-Grid Finite-Difference Seismic Wavefield Estimation Using a Discontinuous Mesh Interface (WEDMI)

Shiying Nie, Yongfei Wang, Kim B. Olsen, & Steven M. Day

Accepted August 30, 2017, SCEC Contribution #7911

In a realistic geological structure with a large contrast in seismic wave speed between shallow and deep regions, simulation of seismic wave propagation using a spatially uniform grid can be computationally very demanding, due to overdiscretization of the high-speed material. Thus, numerical methods that allow for coarser discretization of the faster regions have the potential to be much more efficient. Discontinuous mesh (DM) methods, operating by exchanging wakefield information between media partitions discretized with two different grid\ spacings, provide a convenient way to improve such efficiency issues. Unfortunately, discontinuous staggered-grid finite-difference (FD) methods typically suffer from inherent stability problems, in particular in strongly heterogeneous media, arising from numerical noise generated at the overlap of the two regions with different grid spacing.We developed a 3D fourth-order velocity-stress staggered-grid FD DM method (AWP-DM) to model seismic wakefield estimation using a discontinuous mesh interface (WEDMI) between fine and coarse meshes. Benchmarks in models with realistic 3D velocity variations and finite-fault sources across the grid interface show stable results for a number of timesteps, exceeding the need dictated by current high-frequency ground-motion simulations. In the case of a factor-of-three ratio between the coarse and fine grid sizes, this method is capable of producing a level of accuracy comparable to that from the uniform fine-grid scheme, using at least 7–8 grid points per minimum S wavelength inside the mesh overlap zone.

Key Words
Ground Motion, finite differences, discontinuous mesh

Citation
Nie, S., Wang, Y., Olsen, K. B., & Day, S. M. (2017). Fourth-Order Staggered-Grid Finite-Difference Seismic Wavefield Estimation Using a Discontinuous Mesh Interface (WEDMI). Bulletin of the Seismological Society of America, (accepted).


Related Projects & Working Groups
SEISM2, Ground Motion Prediction (GMP)