Explicit coupling of acoustic and elastic wave propagation in finite difference simulations

06/18/2020
by   Longfei Gao, et al.
0

We present a mechanism to explicitly couple the finite-difference discretizations of 2D acoustic and isotropic elastic wave systems that are separated by straight interfaces. Such coupled simulations allow the application of the elastic model to geological regions that are of special interest for seismic exploration studies (e.g., the areas surrounding salt bodies), while with the computationally more tractable acoustic model still being applied in the background regions. Specifically, the acoustic wave system is expressed in terms of velocity and pressure while the elastic wave system is expressed in terms of velocity and stress. Both systems are posed in first-order forms and discretized on staggered grids. Special variants of the standard finite-difference operators, namely, operators that possess the summation-by-parts property, are used for the approximation of spatial derivatives. Penalty terms, which are also referred to as the simultaneous approximation terms, are designed to weakly impose the elastic-acoustic interface conditions in the finite-difference discretizations and couple the elastic and acoustic wave simulations together. With the presented mechanism, we are able to perform the coupled elastic-acoustic wave simulations stably and accurately. Moreover, it is shown that the energy-conserving property in the continuous systems can be preserved in the discretization with carefully designed penalty terms.

READ FULL TEXT

page 17

page 19

research
12/27/2020

Energy-conserving 3D elastic wave simulation with finite difference discretization on staggered grids with nonconforming interfaces

In this work, we describe an approach to stably simulate the 3D isotropi...
research
09/01/2022

Strongly imposing the free surface boundary condition for wave equations with finite difference operators

Acoustic and elastic wave equations are routinely used in geophysical an...
research
06/21/2023

𝔼^𝐅𝐖𝐈: Multi-parameter Benchmark Datasets for Elastic Full Waveform Inversion of Geophysical Properties

Elastic geophysical properties (such as P- and S-wave velocities) are of...
research
07/16/2022

Accurate simulations of nonlinear dynamic shear ruptures on pre-existing faults in 3D elastic solids with dual-pairing SBP methods

In this paper we derive and analyse efficient and stable numerical metho...
research
04/22/2020

Scaling through abstractions – high-performance vectorial wave simulations for seismic inversion with Devito

[Devito] is an open-source Python project based on domain-specific langu...
research
07/20/2019

Non-stiff narrow-stencil finite difference approximations of the Laplacian on curvilinear multiblock grids

The Laplacian appears in several partial differential equations used to ...

Please sign up or login with your details

Forgot password? Click here to reset