TY - JOUR

T1 - Application of optimal control to inversion of self-potential data: theory and synthetic examples

AU - Malovichko, M. S.

AU - Tarasov, A. V.

AU - Yavich, N. B.

AU - Titov, K. V.

N1 - Publisher Copyright: IEEE

PY - 2021

Y1 - 2021

N2 - Last decades, there has been an increased interest in the use of the self-potential (SP) method in hydrogeophysics. In response to this strong interest, we develop a novel approach to the inversion of SP data. Mathematically, the SP inverse problem is the source identification problem for the Poisson equation. Our approach substantially differs from the standard regularization approach, which explicitly includes the forward-problem operator into the cost functional. We formulated the inverse problem is as an optimal control problem and then translate it into a variational system. The system is approximated in suitable finite-element spaces giving rise to an algebraic problem with the saddle-point structure. In contrast to the standard approach, which leads to a dense linear system, our method results in a system with a sparse block matrix. It can be efficiently solved by either direct sparse solvers or preconditioned iterative solvers. In this paper, we present the formulation of the problem and its finite-element approximation. We discuss the iterative solution and preconditioning strategies. Our software implementation is based on an industrial finite-element package. We also present a numerical experiment with node-based linear finite elements on tetrahedral grids. Our results suggest that the proposed approach may serve as a rapid and reliable tool for large-scale SP inverse problems. Moreover, the same technique can easily be extended to a wide range of geophysical linear inverse problems, such as inversions of magnetic and gravity data.

AB - Last decades, there has been an increased interest in the use of the self-potential (SP) method in hydrogeophysics. In response to this strong interest, we develop a novel approach to the inversion of SP data. Mathematically, the SP inverse problem is the source identification problem for the Poisson equation. Our approach substantially differs from the standard regularization approach, which explicitly includes the forward-problem operator into the cost functional. We formulated the inverse problem is as an optimal control problem and then translate it into a variational system. The system is approximated in suitable finite-element spaces giving rise to an algebraic problem with the saddle-point structure. In contrast to the standard approach, which leads to a dense linear system, our method results in a system with a sparse block matrix. It can be efficiently solved by either direct sparse solvers or preconditioned iterative solvers. In this paper, we present the formulation of the problem and its finite-element approximation. We discuss the iterative solution and preconditioning strategies. Our software implementation is based on an industrial finite-element package. We also present a numerical experiment with node-based linear finite elements on tetrahedral grids. Our results suggest that the proposed approach may serve as a rapid and reliable tool for large-scale SP inverse problems. Moreover, the same technique can easily be extended to a wide range of geophysical linear inverse problems, such as inversions of magnetic and gravity data.

KW - Geologic measurements

KW - Geophysical measurements

KW - Geophysics

KW - Inverse problems

KW - Jacobian matrices

KW - KKT system

KW - Sparse matrices

KW - Standards

KW - inverse problem

KW - optimal control

KW - self-potential

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

U2 - 10.1109/TGRS.2021.3121538

DO - 10.1109/TGRS.2021.3121538

M3 - Article

AN - SCOPUS:85118278520

JO - IEEE Transactions on Geoscience and Remote Sensing

JF - IEEE Transactions on Geoscience and Remote Sensing

SN - 0196-2892

ER -