A Rotated Staggered Grid (RSG) finite-difference elastic wave solver based on Devito

Listed in Datasets

By Oumeng Zhang1, Douglas R Schmitt1

Purdue University

An Optimized Rotated Staggered Grid (RSG) finite-difference elastic wave solver for modeling 2D/3D seismic wave propagation in complex anisotropic media, based on the open-source Python finite-difference package, Devito.

Version 1.0 - published on 14 Nov 2023 doi:10.4231/R8N7-HX28 - cite this Archived on 14 Dec 2023

Licensed under GNU General Public License 3.0

3 citation(s)

226 total view(s), 34 download(s)

saenger2004.png

Description

Built upon the open-source Python finite-difference framework, Devito, this Rotated Staggered Grid elastic wave solver is designed to model seismic wave propagation in intricate anisotropic media.

Key Components:

  • model.py: A refined file derived from Devito's original model.py. It introduces a new object, AnisoSeismicModel, tailored to support the input of 21 distinct elastic parameters for the grid.

  • wavesolver.py: Solver object provides the operators for forward modeling computations. It encapsulates both temporal and spatial discretizations pertinent to a specified problem setup.

  • operators.py: Construct Method for the forward modeling operator handling source excitation and monitoring component recorded by receivers.

  • stiffnessoperator.py: Serving as an operator for the stiffness tensor, this object accommodates isotropy, Vertical Transverse Isotropy (VTI), and 3D rotations. Users can either opt for manual input of the 21 elastic parameters or utilize the automated setup. When paired with a mask, it aids in crafting the AnisoSeismicModel.

Getting Started:

For a practical introduction and hands-on experience, refer to the tutorial Jupyter notebook Tutorial.ipynb. This guide walks users through the application of the package, showcasing wave propagation modeling within a two-layered structure (comprising an upper water layer and a lower TTI media) in 2D.

Open Research Archive (OpenResearchArchive.zip) contains the result and Python code to conduct the numerical tests in terms of stability, numerical dispersion and performance benchmark described in the manuscript: Zhang, O. and Schmitt, D. R. (2023) An optimized 2D/3D Finite-difference Seismic Wave Propagator Using Rotated Staggered Grid for Complex Elastic Anisotropic Structures. ESS Open Archive. November, 2023. File Readme.md details the content in each subdirectory.

Cite this work

Researchers should cite this work as follows:

Tags

  1. anisotropy
  2. EAPS (Dept)
  3. Modeling
  4. modeling and simulation
  5. seismic behavior

Notes

v1.0 - Compatible with up-to-data Devito version >4.6.2  <=4.8.3 (November 2023)

The Purdue University Research Repository (PURR) is a university core research facility provided by the Purdue University Libraries and the Office of the Executive Vice President for Research and Partnerships, with support from additional campus partners.