Combination of the discontinuous Galerkin method with finite differences for simulation of seismic wave propagation
- Institute of Petroleum Geology and Geophysics SB RAS, Novosibirsk (Russian Federation)
- University of Stavanger (Norway)
We present an algorithm for the numerical simulation of seismic wave propagation in models with a complex near surface part and free surface topography. The approach is based on the combination of finite differences with the discontinuous Galerkin method. The discontinuous Galerkin method can be used on polyhedral meshes; thus, it is easy to handle the complex surfaces in the models. However, this approach is computationally intense in comparison with finite differences. Finite differences are computationally efficient, but in general, they require rectangular grids, leading to the stair-step approximation of the interfaces, which causes strong diffraction of the wavefield. In this research we present a hybrid algorithm where the discontinuous Galerkin method is used in a relatively small upper part of the model and finite differences are applied to the main part of the model.
- OSTI ID:
- 22570240
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 311; ISSN JCTPAH; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Unstructured discontinuous Galerkin for seismic inversion.