A fast algorithm for regularized focused 3D inversion of gravity data using randomized singular-value decomposition

Saeed Vatankhah, Rosemary Renaut, Vahid Ebrahimzadeh Ardestani

Research output: Contribution to journalArticle

Abstract

We develop a fast algorithm for solving the under-determined 3D linear gravity inverse problem based on randomized singular-value decomposition (RSVD). The algorithm combines an iteratively reweighted approach for L 1 -norm regularization with the RSVD methodology in which the large-scale linear system at each iteration is replaced with a much smaller linear system. Although the optimal choice for the low-rank approximation of the system matrix with m rows is q = m, acceptable results are achievable with q«m. In contrast to the use of the iterative LSQR algorithm for the solution of linear systems at each iteration, the singular values generated using RSVD yield a good approximation of the dominant singular values of the large-scale system matrix. Thus, the regularization parameter found for the small system at each iteration is dependent on the dominant singular values of the large-scale system matrix and appropriately regularizes the dominant singular space of the large-scale problem. The results achieved are comparable with those obtained using the LSQR algorithm for solving each linear system, but they are obtained at a reduced computational cost. The method has been tested on synthetic models along with real gravity data from the Morro do Engenho complex in central Brazil.

Original languageEnglish (US)
Pages (from-to)G25-G34
JournalGeophysics
Volume83
Issue number4
DOIs
StatePublished - Jul 1 2018

Fingerprint

Singular value decomposition
Linear systems
Gravitation
decomposition
gravity
matrix
Large scale systems
inverse problem
Inverse problems
methodology
inversion
cost
Costs

Keywords

  • 3d
  • Algorithm
  • Gravity
  • Inversion
  • Modeling

ASJC Scopus subject areas

  • Geochemistry and Petrology

Cite this

A fast algorithm for regularized focused 3D inversion of gravity data using randomized singular-value decomposition. / Vatankhah, Saeed; Renaut, Rosemary; Ardestani, Vahid Ebrahimzadeh.

In: Geophysics, Vol. 83, No. 4, 01.07.2018, p. G25-G34.

Research output: Contribution to journalArticle

@article{315206ff59d8452e87258ad62eea8655,
title = "A fast algorithm for regularized focused 3D inversion of gravity data using randomized singular-value decomposition",
abstract = "We develop a fast algorithm for solving the under-determined 3D linear gravity inverse problem based on randomized singular-value decomposition (RSVD). The algorithm combines an iteratively reweighted approach for L 1 -norm regularization with the RSVD methodology in which the large-scale linear system at each iteration is replaced with a much smaller linear system. Although the optimal choice for the low-rank approximation of the system matrix with m rows is q = m, acceptable results are achievable with q«m. In contrast to the use of the iterative LSQR algorithm for the solution of linear systems at each iteration, the singular values generated using RSVD yield a good approximation of the dominant singular values of the large-scale system matrix. Thus, the regularization parameter found for the small system at each iteration is dependent on the dominant singular values of the large-scale system matrix and appropriately regularizes the dominant singular space of the large-scale problem. The results achieved are comparable with those obtained using the LSQR algorithm for solving each linear system, but they are obtained at a reduced computational cost. The method has been tested on synthetic models along with real gravity data from the Morro do Engenho complex in central Brazil.",
keywords = "3d, Algorithm, Gravity, Inversion, Modeling",
author = "Saeed Vatankhah and Rosemary Renaut and Ardestani, {Vahid Ebrahimzadeh}",
year = "2018",
month = "7",
day = "1",
doi = "10.1190/geo2017-0386.1",
language = "English (US)",
volume = "83",
pages = "G25--G34",
journal = "Geophysics",
issn = "0016-8033",
publisher = "Society of Exploration Geophysicists",
number = "4",

}

TY - JOUR

T1 - A fast algorithm for regularized focused 3D inversion of gravity data using randomized singular-value decomposition

AU - Vatankhah, Saeed

AU - Renaut, Rosemary

AU - Ardestani, Vahid Ebrahimzadeh

PY - 2018/7/1

Y1 - 2018/7/1

N2 - We develop a fast algorithm for solving the under-determined 3D linear gravity inverse problem based on randomized singular-value decomposition (RSVD). The algorithm combines an iteratively reweighted approach for L 1 -norm regularization with the RSVD methodology in which the large-scale linear system at each iteration is replaced with a much smaller linear system. Although the optimal choice for the low-rank approximation of the system matrix with m rows is q = m, acceptable results are achievable with q«m. In contrast to the use of the iterative LSQR algorithm for the solution of linear systems at each iteration, the singular values generated using RSVD yield a good approximation of the dominant singular values of the large-scale system matrix. Thus, the regularization parameter found for the small system at each iteration is dependent on the dominant singular values of the large-scale system matrix and appropriately regularizes the dominant singular space of the large-scale problem. The results achieved are comparable with those obtained using the LSQR algorithm for solving each linear system, but they are obtained at a reduced computational cost. The method has been tested on synthetic models along with real gravity data from the Morro do Engenho complex in central Brazil.

AB - We develop a fast algorithm for solving the under-determined 3D linear gravity inverse problem based on randomized singular-value decomposition (RSVD). The algorithm combines an iteratively reweighted approach for L 1 -norm regularization with the RSVD methodology in which the large-scale linear system at each iteration is replaced with a much smaller linear system. Although the optimal choice for the low-rank approximation of the system matrix with m rows is q = m, acceptable results are achievable with q«m. In contrast to the use of the iterative LSQR algorithm for the solution of linear systems at each iteration, the singular values generated using RSVD yield a good approximation of the dominant singular values of the large-scale system matrix. Thus, the regularization parameter found for the small system at each iteration is dependent on the dominant singular values of the large-scale system matrix and appropriately regularizes the dominant singular space of the large-scale problem. The results achieved are comparable with those obtained using the LSQR algorithm for solving each linear system, but they are obtained at a reduced computational cost. The method has been tested on synthetic models along with real gravity data from the Morro do Engenho complex in central Brazil.

KW - 3d

KW - Algorithm

KW - Gravity

KW - Inversion

KW - Modeling

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

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

U2 - 10.1190/geo2017-0386.1

DO - 10.1190/geo2017-0386.1

M3 - Article

AN - SCOPUS:85045904448

VL - 83

SP - G25-G34

JO - Geophysics

JF - Geophysics

SN - 0016-8033

IS - 4

ER -