Seismic wavefield modelling in two‐phase media including undulating topography with the modified Biot/squirt model by a curvilinear‐grid finite difference method

In this paper, we simulated two-dimension numerical on the strong ground motion through the hybrid scheme based on the pseudo-spectral method (PSM) and finite difference method (FDM). We based on the same focal depth, and 2 different thick deposition layers are used as models to analyze the relationship between site situation and the peak displacement of strong ground motion. The results show that the hybrid PSM/FDM method for seismic wavefield simulation combines with advantages of the pseudospectral method and the finite difference method and makes up for the disadvantage of the pseudospectral method and the finite difference method, so this method can process well the calculation of the discontinuous medium surface, then the calculation accuracy is similar to the pseudospectral method. Through the wavefield simulation it is known that the range of the seismic wavefield the peak ground displacement (PGD) of the thicker deposition is larger and the influence of the secondary surface wave at the basin edge is more obvious. The thicker deposition amplitude of strong ground motion in the basin is larger and the duration is longer, and the reflected wave of which is more obvious and stronger. However, the difference of the site condition has little influence to strong ground motion in the horizontal direction.

Download Full-text

A cell-centered finite difference method for a degenerate elliptic equation arising from two-phase mixtures

Computational Geosciences ◽

10.1007/s10596-017-9649-9 ◽

2017 ◽

Vol 21 (4) ◽

pp. 701-712 ◽

Cited By ~ 3

Author(s):

Todd Arbogast ◽

Abraham L. Taicher

Keyword(s):

Elliptic Equation ◽

Finite Difference ◽

Finite Difference Method ◽

Difference Method ◽

Degenerate Elliptic Equation ◽

Two Phase ◽

Degenerate Elliptic ◽

A Cell

Download Full-text

2-D poroelastic wave modelling with a topographic free surface by the curvilinear grid finite-difference method

Geophysical Journal International ◽

10.1093/gji/ggz263 ◽

2019 ◽

Vol 218 (3) ◽

pp. 1961-1982 ◽

Cited By ~ 1

Author(s):

Yao-Chong Sun ◽

Hengxin Ren ◽

Xu-Zhen Zheng ◽

Na Li ◽

Wei Zhang ◽

...

Keyword(s):

Free Surface ◽

Finite Difference ◽

Finite Difference Method ◽

Difference Method ◽

Wave Modelling ◽

Curvilinear Grid

Download Full-text

VARIABLE GRID FINITE DIFFERENCE METHOD FOR TWO-DIMENSIONAL TWO-PHASE IMMISCIBLE FLOW

Acta Mathematica Scientia ◽

10.1016/s0252-9602(17)30591-x ◽

1998 ◽

Vol 18 (4) ◽

pp. 379-386

Author(s):

Wentao Sun

Keyword(s):

Finite Difference ◽

Finite Difference Method ◽

Difference Method ◽

Two Dimensional ◽

Immiscible Flow ◽

Two Phase

Download Full-text

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

Bulletin of the Seismological Society of America ◽

10.1785/0120210096 ◽

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.

Download Full-text