An iterative multilevel method for computing wavefields in frequency-domain seismic inversion

Yogi A. Erlangga, Felix J. Herrmann

Research output: Contribution to journalArticle

14 Citations (Scopus)

Abstract

We describe an iterative multilevel method for solving linear systems representing forward modeling and back propagation of wavefields in frequency-domain seismic inversions. The workhorse of the method is the so-called multilevel Krylov method, applied to a multigrid-preconditioned linear system, and is called multilevel Krylov-multigrid (MKMG) method. Numerical experiments are presented for 2D Marmousi synthetic model for a range of frequencies. The convergence of the method is fast, and depends only mildly on frequency. The method can be considered as the first viable alternative to LU factorization, which is practically prohibitive for 3D seismic inversions.

Original languageEnglish
Pages (from-to)1956-1960
Number of pages5
JournalSEG Technical Program Expanded Abstracts
Volume27
Issue number1
DOIs
Publication statusPublished - Jan 2008
Externally publishedYes

Fingerprint

Iterative methods
Linear systems
inversions
linear systems
Factorization
Backpropagation
multigrid methods
factorization
propagation
back propagation
Experiments
forward modeling
inversion
method
experiment

Keywords

  • Algorithm
  • Finite difference
  • Inversion
  • Modeling
  • Wave equation

ASJC Scopus subject areas

  • Geophysics
  • Geotechnical Engineering and Engineering Geology

Cite this

An iterative multilevel method for computing wavefields in frequency-domain seismic inversion. / Erlangga, Yogi A.; Herrmann, Felix J.

In: SEG Technical Program Expanded Abstracts, Vol. 27, No. 1, 01.2008, p. 1956-1960.

Research output: Contribution to journalArticle

@article{083bacbeefbb43afbc577d8a08b7da97,
title = "An iterative multilevel method for computing wavefields in frequency-domain seismic inversion",
abstract = "We describe an iterative multilevel method for solving linear systems representing forward modeling and back propagation of wavefields in frequency-domain seismic inversions. The workhorse of the method is the so-called multilevel Krylov method, applied to a multigrid-preconditioned linear system, and is called multilevel Krylov-multigrid (MKMG) method. Numerical experiments are presented for 2D Marmousi synthetic model for a range of frequencies. The convergence of the method is fast, and depends only mildly on frequency. The method can be considered as the first viable alternative to LU factorization, which is practically prohibitive for 3D seismic inversions.",
keywords = "Algorithm, Finite difference, Inversion, Modeling, Wave equation",
author = "Erlangga, {Yogi A.} and Herrmann, {Felix J.}",
year = "2008",
month = "1",
doi = "10.1190/1.3059279",
language = "English",
volume = "27",
pages = "1956--1960",
journal = "SEG Technical Program Expanded Abstracts",
issn = "1052-3812",
publisher = "Society of Exploration Geophysicists",
number = "1",

}

TY - JOUR

T1 - An iterative multilevel method for computing wavefields in frequency-domain seismic inversion

AU - Erlangga, Yogi A.

AU - Herrmann, Felix J.

PY - 2008/1

Y1 - 2008/1

N2 - We describe an iterative multilevel method for solving linear systems representing forward modeling and back propagation of wavefields in frequency-domain seismic inversions. The workhorse of the method is the so-called multilevel Krylov method, applied to a multigrid-preconditioned linear system, and is called multilevel Krylov-multigrid (MKMG) method. Numerical experiments are presented for 2D Marmousi synthetic model for a range of frequencies. The convergence of the method is fast, and depends only mildly on frequency. The method can be considered as the first viable alternative to LU factorization, which is practically prohibitive for 3D seismic inversions.

AB - We describe an iterative multilevel method for solving linear systems representing forward modeling and back propagation of wavefields in frequency-domain seismic inversions. The workhorse of the method is the so-called multilevel Krylov method, applied to a multigrid-preconditioned linear system, and is called multilevel Krylov-multigrid (MKMG) method. Numerical experiments are presented for 2D Marmousi synthetic model for a range of frequencies. The convergence of the method is fast, and depends only mildly on frequency. The method can be considered as the first viable alternative to LU factorization, which is practically prohibitive for 3D seismic inversions.

KW - Algorithm

KW - Finite difference

KW - Inversion

KW - Modeling

KW - Wave equation

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

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

U2 - 10.1190/1.3059279

DO - 10.1190/1.3059279

M3 - Article

VL - 27

SP - 1956

EP - 1960

JO - SEG Technical Program Expanded Abstracts

JF - SEG Technical Program Expanded Abstracts

SN - 1052-3812

IS - 1

ER -