Total variation regularization of the 3-D gravity inverse problem using a randomized generalized singular value decomposition

Saeed Vatankhah, Rosemary Renaut, Vahid E. Ardestani

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

We present a fast algorithm for the total variation regularization of the 3-D gravity inverse problem. Through imposition of the total variation regularization, subsurface structures presenting with sharp discontinuities are preserved better than when using a conventional minimumstructure inversion. The associated problem formulation for the regularization is nonlinear but can be solved using an iteratively reweighted least-squares algorithm. For small-scale problems the regularized least-squares problem at each iteration can be solved using the generalized singular value decomposition. This is not feasible for large-scale, or even moderate-scale, problems. Instead we introduce the use of a randomized generalized singular value decomposition in order to reduce the dimensions of the problem and provide an effective and efficient solution technique. For further efficiency an alternating direction algorithm is used to implement the total variation weighting operator within the iteratively reweighted least-squares algorithm. Presented results for synthetic examples demonstrate that the novel randomized decomposition provides good accuracy for reduced computational and memory demands as compared to use of classical approaches.

Original languageEnglish (US)
Pages (from-to)695-705
Number of pages11
JournalGeophysical Journal International
Volume213
Issue number1
DOIs
StatePublished - Apr 1 2018

Fingerprint

Singular value decomposition
inverse problem
Inverse problems
Gravitation
decomposition
gravity
gravitation
iteration
Mathematical operators
discontinuity
inversions
formulations
Data storage equipment
operators

Keywords

  • Asia
  • Gravity anomalies and Earth structure
  • Inverse theory
  • Numerical approximations and analysis

ASJC Scopus subject areas

  • Geophysics
  • Geochemistry and Petrology

Cite this

Total variation regularization of the 3-D gravity inverse problem using a randomized generalized singular value decomposition. / Vatankhah, Saeed; Renaut, Rosemary; Ardestani, Vahid E.

In: Geophysical Journal International, Vol. 213, No. 1, 01.04.2018, p. 695-705.

Research output: Contribution to journalArticle

@article{4ec1ddb678f94ea281f53ce2e86830af,
title = "Total variation regularization of the 3-D gravity inverse problem using a randomized generalized singular value decomposition",
abstract = "We present a fast algorithm for the total variation regularization of the 3-D gravity inverse problem. Through imposition of the total variation regularization, subsurface structures presenting with sharp discontinuities are preserved better than when using a conventional minimumstructure inversion. The associated problem formulation for the regularization is nonlinear but can be solved using an iteratively reweighted least-squares algorithm. For small-scale problems the regularized least-squares problem at each iteration can be solved using the generalized singular value decomposition. This is not feasible for large-scale, or even moderate-scale, problems. Instead we introduce the use of a randomized generalized singular value decomposition in order to reduce the dimensions of the problem and provide an effective and efficient solution technique. For further efficiency an alternating direction algorithm is used to implement the total variation weighting operator within the iteratively reweighted least-squares algorithm. Presented results for synthetic examples demonstrate that the novel randomized decomposition provides good accuracy for reduced computational and memory demands as compared to use of classical approaches.",
keywords = "Asia, Gravity anomalies and Earth structure, Inverse theory, Numerical approximations and analysis",
author = "Saeed Vatankhah and Rosemary Renaut and Ardestani, {Vahid E.}",
year = "2018",
month = "4",
day = "1",
doi = "10.1093/gji/ggy014",
language = "English (US)",
volume = "213",
pages = "695--705",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "1",

}

TY - JOUR

T1 - Total variation regularization of the 3-D gravity inverse problem using a randomized generalized singular value decomposition

AU - Vatankhah, Saeed

AU - Renaut, Rosemary

AU - Ardestani, Vahid E.

PY - 2018/4/1

Y1 - 2018/4/1

N2 - We present a fast algorithm for the total variation regularization of the 3-D gravity inverse problem. Through imposition of the total variation regularization, subsurface structures presenting with sharp discontinuities are preserved better than when using a conventional minimumstructure inversion. The associated problem formulation for the regularization is nonlinear but can be solved using an iteratively reweighted least-squares algorithm. For small-scale problems the regularized least-squares problem at each iteration can be solved using the generalized singular value decomposition. This is not feasible for large-scale, or even moderate-scale, problems. Instead we introduce the use of a randomized generalized singular value decomposition in order to reduce the dimensions of the problem and provide an effective and efficient solution technique. For further efficiency an alternating direction algorithm is used to implement the total variation weighting operator within the iteratively reweighted least-squares algorithm. Presented results for synthetic examples demonstrate that the novel randomized decomposition provides good accuracy for reduced computational and memory demands as compared to use of classical approaches.

AB - We present a fast algorithm for the total variation regularization of the 3-D gravity inverse problem. Through imposition of the total variation regularization, subsurface structures presenting with sharp discontinuities are preserved better than when using a conventional minimumstructure inversion. The associated problem formulation for the regularization is nonlinear but can be solved using an iteratively reweighted least-squares algorithm. For small-scale problems the regularized least-squares problem at each iteration can be solved using the generalized singular value decomposition. This is not feasible for large-scale, or even moderate-scale, problems. Instead we introduce the use of a randomized generalized singular value decomposition in order to reduce the dimensions of the problem and provide an effective and efficient solution technique. For further efficiency an alternating direction algorithm is used to implement the total variation weighting operator within the iteratively reweighted least-squares algorithm. Presented results for synthetic examples demonstrate that the novel randomized decomposition provides good accuracy for reduced computational and memory demands as compared to use of classical approaches.

KW - Asia

KW - Gravity anomalies and Earth structure

KW - Inverse theory

KW - Numerical approximations and analysis

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

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

U2 - 10.1093/gji/ggy014

DO - 10.1093/gji/ggy014

M3 - Article

AN - SCOPUS:85052851334

VL - 213

SP - 695

EP - 705

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 1

ER -