Compressive simultaneous full-waveform simulation

Felix J. Herrmann, Yogi A. Erlangga, Tim T Y Lin

Research output: Contribution to journalArticle

84 Citations (Scopus)

Abstract

The fact that the computational complexity of wavefield simulation is proportional to the size of the discretized model and acquisition geometry and not to the complexity of the simulated wavefield is a major impediment within seismic imaging. By turning simulation into a compressive sensing problem, where simulated data are recovered from a relatively small number of independent simultaneous sources, we remove this impediment by showing that compressively sampling a simulation is equivalent to compressively sampling the sources, followed by solving a reduced system. As in compressive sensing, this reduces sampling rate and hence simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by recovering the full wavefield with a program that promotes sparsity. Depending on the wavefield's sparsity, this approach can lead to significant cost reductions, particularly when combined with the implicit preconditioned Helmholtz solver, which is known to converge even for decreasing mesh sizes and increasing angular frequencies. These properties make our scheme a viable alternative to explicit time-domain finite differences.

Original languageEnglish
JournalGeophysics
Volume74
Issue number4
DOIs
Publication statusPublished - 2009
Externally publishedYes

Fingerprint

waveforms
Sampling
sampling
simulation
Cost reduction
cost reduction
Computational complexity
mesh size
Imaging techniques
mesh
acquisition
Geometry
costs
harmonics
geometry
Costs
cost

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geophysics

Cite this

Herrmann, F. J., Erlangga, Y. A., & Lin, T. T. Y. (2009). Compressive simultaneous full-waveform simulation. Geophysics, 74(4). https://doi.org/10.1190/1.3115122

Compressive simultaneous full-waveform simulation. / Herrmann, Felix J.; Erlangga, Yogi A.; Lin, Tim T Y.

In: Geophysics, Vol. 74, No. 4, 2009.

Research output: Contribution to journalArticle

Herrmann, FJ, Erlangga, YA & Lin, TTY 2009, 'Compressive simultaneous full-waveform simulation', Geophysics, vol. 74, no. 4. https://doi.org/10.1190/1.3115122
Herrmann, Felix J. ; Erlangga, Yogi A. ; Lin, Tim T Y. / Compressive simultaneous full-waveform simulation. In: Geophysics. 2009 ; Vol. 74, No. 4.
@article{a364311e88064d98868985dc70379728,
title = "Compressive simultaneous full-waveform simulation",
abstract = "The fact that the computational complexity of wavefield simulation is proportional to the size of the discretized model and acquisition geometry and not to the complexity of the simulated wavefield is a major impediment within seismic imaging. By turning simulation into a compressive sensing problem, where simulated data are recovered from a relatively small number of independent simultaneous sources, we remove this impediment by showing that compressively sampling a simulation is equivalent to compressively sampling the sources, followed by solving a reduced system. As in compressive sensing, this reduces sampling rate and hence simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by recovering the full wavefield with a program that promotes sparsity. Depending on the wavefield's sparsity, this approach can lead to significant cost reductions, particularly when combined with the implicit preconditioned Helmholtz solver, which is known to converge even for decreasing mesh sizes and increasing angular frequencies. These properties make our scheme a viable alternative to explicit time-domain finite differences.",
author = "Herrmann, {Felix J.} and Erlangga, {Yogi A.} and Lin, {Tim T Y}",
year = "2009",
doi = "10.1190/1.3115122",
language = "English",
volume = "74",
journal = "Geophysics",
issn = "0016-8033",
publisher = "Society of Exploration Geophysicists",
number = "4",

}

TY - JOUR

T1 - Compressive simultaneous full-waveform simulation

AU - Herrmann, Felix J.

AU - Erlangga, Yogi A.

AU - Lin, Tim T Y

PY - 2009

Y1 - 2009

N2 - The fact that the computational complexity of wavefield simulation is proportional to the size of the discretized model and acquisition geometry and not to the complexity of the simulated wavefield is a major impediment within seismic imaging. By turning simulation into a compressive sensing problem, where simulated data are recovered from a relatively small number of independent simultaneous sources, we remove this impediment by showing that compressively sampling a simulation is equivalent to compressively sampling the sources, followed by solving a reduced system. As in compressive sensing, this reduces sampling rate and hence simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by recovering the full wavefield with a program that promotes sparsity. Depending on the wavefield's sparsity, this approach can lead to significant cost reductions, particularly when combined with the implicit preconditioned Helmholtz solver, which is known to converge even for decreasing mesh sizes and increasing angular frequencies. These properties make our scheme a viable alternative to explicit time-domain finite differences.

AB - The fact that the computational complexity of wavefield simulation is proportional to the size of the discretized model and acquisition geometry and not to the complexity of the simulated wavefield is a major impediment within seismic imaging. By turning simulation into a compressive sensing problem, where simulated data are recovered from a relatively small number of independent simultaneous sources, we remove this impediment by showing that compressively sampling a simulation is equivalent to compressively sampling the sources, followed by solving a reduced system. As in compressive sensing, this reduces sampling rate and hence simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by recovering the full wavefield with a program that promotes sparsity. Depending on the wavefield's sparsity, this approach can lead to significant cost reductions, particularly when combined with the implicit preconditioned Helmholtz solver, which is known to converge even for decreasing mesh sizes and increasing angular frequencies. These properties make our scheme a viable alternative to explicit time-domain finite differences.

UR - http://www.scopus.com/inward/record.url?scp=67650327365&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=67650327365&partnerID=8YFLogxK

U2 - 10.1190/1.3115122

DO - 10.1190/1.3115122

M3 - Article

AN - SCOPUS:67650327365

VL - 74

JO - Geophysics

JF - Geophysics

SN - 0016-8033

IS - 4

ER -