A matrix-free reformulation of the multi-parameter descent and conjugate-gradient method for isotropic elastic iterative reverse-time migration

A matrix-free reformulation of the multi-parameter descent and conjugate-gradient method for isotropic elastic iterative reverse-time migration

Mulder, W.A. (TU Delft Applied Geophysics and Petrophysics; Shell Global Solutions International B.V.)

2021

Multi-parameter inversion of linear systems appears in many problems. The focus here is on isotropic elastic iterative reverse-time migration for three position-dependent subsurface model parameters, which amounts to data fitting of processed seismic data with synthetics from the Born approximation of the elastic wave equation. In that case, the matrix of the linear system is the hessian. As it is impractical to form, a matrix-free formulation is needed, which is readily derived for the gradient descent method. For single-parameter inversion, the conjugate-gradient (CG) method is generally more efficient than simple descent. However, the multiple-parameter CG method has a significantly higher cost than the descent method. Here, first a matrix-free data-domain reformulation is derived. Then, its performance is compared to the simple descent method to see of its faster convergence justifies the higher cost. A comparison on a marine 2-D toy problem with a salt body and sea-bottom receivers shows that the multiple-parameter descent method wins in terms of efficiency if the number of iterations is limited and that the single-parameter CG method is even faster.

http://resolver.tudelft.nl/uuid:24959b0b-6eb2-4082-8aae-c5047e881767

https://doi.org/10.3997/2214-4609.202112725

EAGE

2022-04-18

Proceedings of the 82nd EAGE Annual Conference & Exhibition: October 18-21, 2021

82nd EAGE Annual Conference & Exhibition
, 2021-10-18 → 2021-10-21, Amsterdam, Netherlands

Accepted Author Manuscript

Institutional Repository

conference paper

© 2021 W.A. Mulder