Fully scalable solver for frequency-domain visco-elastic wave equations in 3D heterogeneous media: A controllability approach

We develop a controllability strategy for the computation of frequency-domain solutions of the 3D visco-elastic wave equation, in the perspective of seismic imaging applications. We generalize the controllability results for such equations beyond the sound-soft scattering (obstacle) problem. We deta...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Journal of computational physics Jg. 468; S. 111514
Hauptverfasser: Tang, Jet Hoe, Brossier, Romain, Métivier, Ludovic
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Cambridge Elsevier Inc 01.11.2022
Elsevier Science Ltd
Elsevier
Schlagworte:
Algorithms
Computational physics
Computer Science
Conjugate gradient
Controllability
Discretization
Domain decomposition methods
Elastic wave equations
Elastic waves
Frequency domain analysis
Full waveform inversion
Geophysics
Iterative methods
Modeling and Simulation
Parallelization
Physics
Runge-Kutta method
Spectral element
Viscoelasticity
Wave equations
Waveforms
Controllability
Conjugate gradient
Full waveform inversion
Parallelization
Elastic wave equations
Spectral element
ISSN:0021-9991, 1090-2716
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We develop a controllability strategy for the computation of frequency-domain solutions of the 3D visco-elastic wave equation, in the perspective of seismic imaging applications. We generalize the controllability results for such equations beyond the sound-soft scattering (obstacle) problem. We detail the conjugate gradient implementation and show how an inner elliptic problem needs to be solved to compute the Riesz representative of the gradient at each iteration. We select a spectral-element spatial discretization and a fourth-order Runge-Kutta time discretization. We implement the controllability method in the framework of the SEM46 full waveform modeling and inversion software, to inherit for its excellent scalability which relies on an efficient domain decomposition algorithm. We perform a series of numerical experiments to validate the approach and illustrate its scalability up to more than fifteen hundred cores. In this case, with an elapsed time of less than 50 minutes, we solve a problem on a cubic domain containing up to 160 wavelengths in each direction, involving more than 1.7 billion unknowns. •Efficient numerical solution of frequency-domain visco-elastic wave equation.•Controllability methods (CM) extended to general visco-elastic problem beyond sound-soft scattering problems.•CM combined with spectral element method, explicit time scheme and parallelization through domain decomposition.•CM offer robust, non-intrusive approach for time-harmonic wave equations.•CM achieve strong scalability on massively parallel architectures.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2022.111514