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

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Journal of computational physics Ročník 468; s. 111514
Hlavní autori: Tang, Jet Hoe, Brossier, Romain, Métivier, Ludovic
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Cambridge Elsevier Inc 01.11.2022
Elsevier Science Ltd
Elsevier
Predmet:
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
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí: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.
Bibliografia: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