Numerical Modeling of Wave Propagation in Poroelastic Media with a Staggered-Grid Finite-Difference Method

2013 ◽  
Vol 275-277 ◽  
pp. 612-617
Author(s):  
Wen Sheng Zhang ◽  
Li Tong
Keyword(s):  
Numerical Modeling ◽  
Wave Propagation ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Staggered Grid ◽  
Poroelastic Medium ◽  
Numerical Computations ◽  
Poroelastic Media ◽  
Biot’S Equations

In this paper, wave propagation in poroelastic medium is simulated with a staggered-grid finite-difference method. The formulation is discretized based on the second-order Biot’s equations rather than the corresponding velocity-stress form. In order to eliminate boundary reflections, the PML method is applied. Numerical computations are implemented and the results show the correctness and effectiveness of the schemes presented in this paper.

A High-Order Finite-Difference Method on Staggered Curvilinear Grids for Seismic Wave Propagation Applications with Topography

2021 ◽  
Author(s):  
Ossian O’Reilly ◽  
Te-Yang Yeh ◽  
Kim B. Olsen ◽  
Zhifeng Hu ◽  
Alex Breuer ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Vertical Axis ◽  
Cartesian Grid ◽  
Staggered Grid ◽  
Graphic Processing Units ◽  
Curvilinear Grids ◽  
Curvilinear Grid

ABSTRACT We developed a 3D elastic wave propagation solver that supports topography using staggered curvilinear grids. Our method achieves comparable accuracy to the classical fourth-order staggered grid velocity–stress finite-difference method on a Cartesian grid. We show that the method is provably stable using summation-by-parts operators and weakly imposed boundary conditions via penalty terms. The maximum stable timestep obeys a relationship that depends on the topography-induced grid stretching along the vertical axis. The solutions from the approach are in excellent agreement with verified results for a Gaussian-shaped hill and for a complex topographic model. Compared with a Cartesian grid, the curvilinear grid adds negligible memory requirements, but requires longer simulation times due to smaller timesteps for complex topography. The code shows 94% weak scaling efficiency up to 1014 graphic processing units.

scholarly journals Performance of 3D Wave Field Modeling Using the Staggered Grid Finite Difference Method with General-Purpose Processors

Energies ◽  
2020 ◽  
Vol 13 (17) ◽  
pp. 4573
Author(s):  
Anna Franczyk ◽  
Damian Gwiżdż ◽  
Andrzej Leśniak
Keyword(s):  
Numerical Modeling ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Wave Field ◽  
Difference Method ◽  
General Purpose ◽  
Staggered Grid ◽  
Computational Time ◽  
Computing Environment ◽  
General Purpose Processors

This paper aims to provide a quantitative understanding of the performance of numerical modeling of a wave field equation using general-purpose processors. In particular, this article presents the most important aspects related to the memory workloads and execution time of the numerical modeling of both acoustic and fully elastic waves in isotropic and anisotropic mediums. The results presented in this article were calculated for the staggered grid finite difference method. Our results show that the more realistic the seismic wave simulations that are performed, the more the demand for memory and the computational capacity of the computing environment increases. The results presented in this article allow the estimation of the memory requirements and computational time of wavefield modeling for the considered model (acoustic, elastic or anisotropic) so that their feasibility can be assessed in a given computing environment and within an acceptable time. Understanding the numerical modeling performance is especially important when graphical processing units (GPU) are utilized to satisfy the intensive calculations of three-dimensional seismic forward modeling.

Sign in / Sign up

Export Citation Format

Share Document