GCH_gravinv: A MATLAB-based program for inverting gravity anomalies over sedimentary basins

This paper presents a Matlab-based program GCH_gravinv including an easy-to-use graphical user interface (GUI) for determining the depth to the basement of a sedimentary basin derived from inverting its gravity anomalies by an iterative procedure. The developed code uses an advanced iterative rapid...

Full description

Saved in:
Bibliographic Details
Published in:Computers & geosciences 2018-11, Vol.120, p.40-47
Main Authors: Pham, Luan Thanh, Oksum, Erdinc, Do, Thanh Duc
Format: Article
Language:English
Subjects:
Basement depth
Chintalapudi
Gravity inversion
Matlab
Sedimentary basin
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:This paper presents a Matlab-based program GCH_gravinv including an easy-to-use graphical user interface (GUI) for determining the depth to the basement of a sedimentary basin derived from inverting its gravity anomalies by an iterative procedure. The developed code uses an advanced iterative rapid algorithm based on the combination of the FFT-based algorithm of Granser and the space domain technique of Cordell and Henderson. As an advantage, this combination performs a fast computation with high precision and the inversion scheme does not require a mean depth or low-pass filtering. Given the required density contrast varying exponentially with depth and a pre-assigned criterion for the termination of the iterative procedure, the GUI allows the user an instant view of the convergence between the observed and inverted anomalies. The feasibility of the proposed algorithm is demonstrated on synthetic data from a 3D model where the obtained results coincide well with the actual depths. The inversion of a real gravity anomaly over Chintalapudi subbasin (India) is presented to compute the basement depths as a practical example. The obtained results from the real data application are in agreement with published information available for the study area. •We developed a MATLAB program to invert gravity anomalies over sedimentary basins.•The algorithm combines the advantages of both the FFT-based and space domain techniques.•It can be applied to the case where density contrast varies exponentially with depth.•The inversion scheme is independent of mean depth and low-pass filtering.•The code provided good result over 3D synthetic model as well as on real gravity data.
ISSN:0098-3004
DOI:10.1016/j.cageo.2018.07.009