spyro: a Firedrake-based wave propagation and full-waveform-inversion finite-element solver

In this article, we introduce spyro, a software stack to solve wave propagation in heterogeneous domains and perform full waveform inversion (FWI) employing the finite-element framework from Firedrake, a high-level Python package for the automated solution of partial differential equations using the...

Full description

Saved in:
Bibliographic Details
Published in:Geoscientific Model Development 2022-11, Vol.15 (23), p.8639-8667
Main Authors: Roberts, Keith J, Olender, Alexandre, Franceschini, Lucas, Kirby, Robert C, Gioria, Rafael S, Carmo, Bruno S
Format: Article
Language:English
Subjects:
Basis functions
Differential equations
Equations of state
Finite element method
Geophysics
Inverse problems
Lumping
Material properties
Mathematical analysis
Mathematical models
Methods
P waves
Partial differential equations
Propagation
Seismic response
Seismic velocities
Seismic waves
Software
Velocity
Wave propagation
Waveforms
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:In this article, we introduce spyro, a software stack to solve wave propagation in heterogeneous domains and perform full waveform inversion (FWI) employing the finite-element framework from Firedrake, a high-level Python package for the automated solution of partial differential equations using the finite-element method. The capability of the software is demonstrated by using a continuous Galerkin approach to perform FWI for seismic velocity model building, considering realistic geophysics examples. A time domain FWI approach that uses meshes composed of variably sized triangular elements to discretize the domain is detailed. To resolve both the forward and adjoint-state equations and to calculate a mesh-independent gradient associated with the FWI process, a fully explicit, variable higher-order (up to degree k=5 in 2D and k=3 in 3D) mass-lumping method is used. We show that, by adapting the triangular elements to the expected peak source frequency and properties of the wave field (e.g., local P-wave speed) and by leveraging higher-order basis functions, the number of degrees of freedom necessary to discretize the domain can be reduced. Results from wave simulations and FWIs in both 2D and 3D highlight our developments and demonstrate the benefits and challenges with using triangular meshes adapted to the material properties.
ISSN:1991-9603
1991-959X
1991-962X
1991-9603
1991-962X
DOI:10.5194/gmd-15-8639-2022