Improving the use of the randomized singular value decomposition for the inversion of gravity and magnetic data

06/26/2019
by   Saeed Vatankhah, et al.
0

The large-scale focusing inversion of gravity and magnetic potential field data using L_1-norm regularization is considered. The use of the randomized singular value decomposition methodology facilitates tackling the computational challenge that arises in the solution of these large-scale inverse problems. As such the powerful randomized singular value decomposition is used for the numerical solution of all linear systems required in the algorithm. A comprehensive comparison of the developed methodology for the inversion of magnetic and gravity data is presented. These results indicate that there is generally an important difference between the gravity and magnetic inversion problems. Specifically, the randomized singular value decomposition is dependent on the generation of a rank q approximation to the underlying model matrix, and the results demonstrate that q needs to be larger, for equivalent problem sizes, for the magnetic problem as compared to the gravity problem. Without a relatively large q the dominant singular values of the magnetic model matrix are not well-approximated. The comparison also shows how the use of the power iteration embedded within the randomized algorithm is used to improve the quality of the resulting dominant subspace approximation, especially in magnetic inversion, yielding acceptable approximations for smaller choices of q. The price to pay is the trade-off between approximation accuracy and computational cost. The algorithm is applied for the inversion of magnetic data obtained over a portion of the Wuskwatim Lake region in Manitoba, Canada

READ FULL TEXT

page 12

page 21

page 22

page 23

page 24

page 25

page 26

page 27

research
09/04/2019

Regularized Linear Inversion with Randomized Singular Value Decomposition

In this work, we develop efficient solvers for linear inverse problems b...
research
12/12/2019

Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory

A fast non-convex low-rank matrix decomposition method for potential fie...
research
08/19/2021

PyParSVD: A streaming, distributed and randomized singular-value-decomposition library

We introduce PyParSVD[https://github.com/Romit-Maulik/PyParSVD], a Pytho...
research
01/10/2020

Generalized L_p-norm joint inversion of gravity and magnetic data using cross-gradient constraint

A generalized unifying approach for L_p-norm joint inversion of gravity ...
research
08/27/2023

Learning end-to-end inversion of circular Radon transforms in the partial radial setup

We present a deep learning-based computational algorithm for inversion o...

Please sign up or login with your details

Forgot password? Click here to reset