Three-dimensional gravity modelling and focusing inversion using rectangular meshes

• Published in: Geophysical Prospecting Vol. 59; no. 5; pp. 966 - 979
• Main Author: Commer, Michael
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Publishing Ltd 01.09.2011
• Subjects: ALGORITHMS; Conjugate gradients; DECAY; Density; ENVIRONMENTAL SCIENCES; FLEXIBILITY; FOCUSING; GEOSCIENCES; Gradiometry; Gravitation; Gravity; IMPLEMENTATION; Imposition; Inversion; Inversions; KERNELS; Mathematical models; RESOLUTION; SIMULATION; TRANSFORMATIONS; VECTORS; WEIGHTING FUNCTIONS
• ISSN: 0016-8025; 1365-2478
• DOI: 10.1111/j.1365-2478.2011.00969.x

ABSTRACT Rectangular grid cells are commonly used for the geophysical modeling of gravity anomalies, owing to their flexibility in constructing complex models. The straightforward handling of cubic cells in gravity inversion algorithms allows for a flexible imposition of model regularization constraints, which are generally essential in the inversion of static potential field data. The first part of this paper provides a review of commonly used expressions for calculating the gravity of a right polygonal prism, both for gravity and gradiometry, where the formulas of Plouff and Forsberg are adapted. The formulas can be cast into general forms practical for implementation. In the second part, a weighting scheme for resolution enhancement at depth is presented. Modelling the earth using highly digitized meshes, depth weighting schemes are typically applied to the model objective functional, subject to minimizing the data misfit. The scheme proposed here involves a non‐linear conjugate gradient inversion scheme with a weighting function applied to the non‐linear conjugate gradient scheme's gradient vector of the objective functional. The low depth resolution due to the quick decay of the gravity kernel functions is counteracted by suppressing the search directions in the parameter space that would lead to near‐surface concentrations of gravity anomalies. Further, a density parameter transformation function enabling the imposition of lower and upper bounding constraints is employed. Using synthetic data from models of varying complexity and a field data set, it is demonstrated that, given an adequate depth weighting function, the gravity inversion in the transform space can recover geologically meaningful models requiring a minimum of prior information and user interaction.