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

by   Longfei Gao, et al.

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.



There are no comments yet.


page 17

page 19


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...

Upwind Summation By Parts Finite Difference Methods for Large Scale Elastic Wave Simulations In Complex Geometries

High-order accurate summation-by-parts (SBP) finite difference (FD) meth...

Elastic wave propagation in curvilinear coordinates with mesh refinement interfaces by a fourth order finite difference method

We develop a fourth order accurate finite difference method for the thre...

3D Acoustic-Elastic Coupling with Gravity: The Dynamics of the 2018 Palu, Sulawesi Earthquake and Tsunami

We present a highly scalable 3D fully-coupled Earth ocean model of e...

A comparative study of two-dimensional vocal tract acoustic modeling based on Finite-Difference Time-Domain methods

The two-dimensional (2D) numerical approaches for vocal tract (VT) model...

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...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.