On Godunov-type finite volume methods for seismic wave propagation

Juan B. Camargo,Pedro S. Peixoto,Felipe A. G. Silva

2024-07-26

Abstract:The computational complexity of simulating seismic waves demands continual exploration of more efficient numerical methods. While Finite Volume methods are widely acclaimed for tackling general nonlinear hyperbolic (wave) problems, their application in realistic seismic wave simulation remains uncommon, with rare investigations in the literature. Furthermore, seismic wavefields are influenced by sharp subsurface interfaces frequently encountered in realistic models, which could, in principle, be adequately solved with Finite Volume methods. In this study, we delved into two Finite Volume (FV) methods to assess their efficacy and competitiveness in seismic wave simulations, compared to traditional Finite Difference schemes. We investigated Gudunov-type FV methods: an upwind method called wave propagation algorithm (WPA), and a Central-Upwind type method (CUp). Our numerical analysis uncovered that these finite volume methods could provide less dispersion (albeit increased dissipation) compared to finite differences for seismic problems characterized by velocity profiles with abrupt transitions in the velocity. However, when applied to more realistic seismic models, finite volume methods yielded unfavorable outcomes compared to finite difference methods, the latter offering lower computational costs and higher accuracy. This highlights that despite the potential advantages of finite volume methods, such as their conservative nature and aptitude for accurately capturing shock waves in specific contexts, our results indicate that they are only advantageous for seismic simulations when unrealistic abrupt transitions are present in the velocity models.

Numerical Analysis,Geophysics

What problem does this paper attempt to address?

### Problems the Paper Aims to Solve

This paper aims to explore the effectiveness of Godunov-type finite volume methods in simulating seismic wave propagation and compare them with traditional finite difference methods. Specifically:

1. **Exploration of Numerical Methods**:
- The researchers evaluated two Godunov-type finite volume methods: the Wave Propagation Algorithm (WPA) and the Central-Upwind (CUp) method, to determine their effectiveness and competitiveness in seismic wave simulation.

2. **Background of the Problem**:
- Seismic wave simulation is a complex process that typically requires solving the elastic wave equation. Although finite difference methods are widely used due to their low computational cost and simplicity, they may fail when dealing with complex earth structures with large velocity contrasts, strong heterogeneity, or topographic attenuation.

3. **Research Motivation**:
- Finite volume methods have advantages in handling discontinuous media and can better capture phenomena such as shock waves. However, their application in seismic wave simulation is relatively rare. Therefore, this paper attempts to fill this research gap by evaluating the performance of finite volume methods in seismic wave simulation.

4. **Main Findings**:
- The research results indicate that for idealized models with sudden velocity changes, finite volume methods exhibit lower dispersion but higher dissipation effects. However, in more realistic seismic models, finite volume methods are less effective than finite difference methods, which offer lower computational costs and higher accuracy.

Through the above research, the paper aims to provide new insights into the field of seismic wave propagation simulation, particularly regarding the performance differences of various numerical methods in handling complex media with discontinuities.