Gravity inversion method using L0-norm constraint with auto-adaptive regularization and combined stopping criteria

We present a gravity inversion method that can produce compact and sharp images to assist the modeling of non-smooth geologic features. The proposed iterative inversion approach makes use of L0-norm-stabilizing functional, hard and physical parameter inequality constraints and a depth-weighting func...

Full description

Saved in:
Bibliographic Details
Published in:Solid earth (Göttingen) 2023-02, Vol.14 (2), p.101-117
Main Authors: Mesay Geletu Gebre, Lewi, Elias
Format: Article
Language:English
Subjects:
Algorithms
Earth science
Geometry
Gravity
Gravity data
Inverse problems
Iterative methods
Methods
Noise levels
Parameters
Physical properties
Regularization
Stabilizing
Weighting functions
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We present a gravity inversion method that can produce compact and sharp images to assist the modeling of non-smooth geologic features. The proposed iterative inversion approach makes use of L0-norm-stabilizing functional, hard and physical parameter inequality constraints and a depth-weighting function. The method incorporates an auto-adaptive regularization technique, which automatically determines a suitable regularization parameter and error-weighting function that helps to improve both the stability and convergence of the method. The auto-adaptive regularization and error-weighting matrix are not dependent on the known noise level. Because of that, the method yields reasonable results even if the noise level of the data is not known properly. The utilization of an effectively combined stopping rule to terminate the inversion process is another improvement that is introduced in this work. The capacity and the efficiency of the new inversion method were tested by inverting randomly chosen synthetic and measured data. The synthetic test models consist of multiple causative blocky bodies, with different geometries and density distributions that are vertically and horizontally distributed adjacent to each other. Inversion results of the synthetic data show that the developed method can recover models that adequately match the real geometry, location and densities of the synthetic causative bodies. Furthermore, the testing of the improved approach using published real gravity data confirmed the potential and practicality of the method in producing compact and sharp inverse images of the subsurface.
ISSN:1869-9510
1869-9529
DOI:10.5194/se-14-101-2023