Compressive simultaneous full-waveform simulation

TitleCompressive simultaneous full-waveform simulation
Publication TypeConference
Year of Publication2009
AuthorsTim T.Y. Lin, Yogi A. Erlangga, Felix J. Herrmann
Conference NameSEG Technical Program Expanded Abstracts
KeywordsPresentation, SEG

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 is 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 allows for a reduction in sampling rate and hence in simulation costs. We demonstrate this principle for the time-harmonic Helmholtz solver. The solution is computed by inverting the reduced system, followed by a recovery of the full wavefield with a sparsity promoting program. Depending on the wavefield’s sparsity, this approach can lead to a significant cost reduction, in particular 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-difference.

Citation Keylin2009SEGcsf