- Journal: Geophysics
- Source code: pinga-lab/paper-planting-densities
- Supplement: 10.6084/m9.figshare.91574
- PDF: paper-planting-anomalous-densities-2012.pdf
- doi: 10.1190/geo2011-0388.1

Uieda, L., and V. C. F. Barbosa (2012), Robust 3D gravity gradient inversion by planting anomalous densities, Geophysics, 77(4), G55-G66, doi:10.1190/geo2011-0388.1

- How two gravity-gradient inversion methods can be used to reveal different geologic features of ore deposit — A case study from the Quadrilátero Ferrífero (Brazil)
- Imaging iron ore from the Quadrilátero Ferrífero (Brazil) using geophysical inversion and drill hole data

The inversion method proposed in this paper is implemented in the open-source
Fatiando a Terra Python library
as the `fatiando.gravmag.harvester`

module.
The module was introduced in
version 0.1
of the library.

The following is an animation of the growth algorithm during the inversion of synthetic data. The video is available at figshare: 10.6084/m9.figshare.91469

We have developed a new gravity gradient inversion method for estimating a 3D density-contrast distribution defined on a grid of rectangular prisms. Our method consists of an iterative algorithm that does not require the solution of an equation system. Instead, the solution grows systematically around user-specified prismatic elements, called “seeds,” with given density contrasts. Each seed can be assigned a different density-contrast value, allowing the interpretation of multiple sources with different density contrasts and that produce interfering signals. In real world scenarios, some sources might not be targeted for the interpretation. Thus, we developed a robust procedure that neither requires the isolation of the signal of the targeted sources prior to the inversion nor requires substantial prior information about the nontargeted sources. In our iterative algorithm, the estimated sources grow by the accretion of prisms in the periphery of the current estimate. In addition, only the columns of the sensitivity matrix corresponding to the prisms in the periphery of the current estimate are needed for the computations. Therefore, the individual columns of the sensitivity matrix can be calculated on demand and deleted after an accretion takes place, greatly reducing the demand for computer memory and processing time. Tests on synthetic data show the ability of our method to correctly recover the geometry of the targeted sources, even when interfering signals produced by nontargeted sources are present. Inverting the data from an airborne gravity gradiometry survey flown over the iron ore province of Quadrilátero Ferrífero, southeastern Brazil, we estimated a compact iron ore body that is in agreement with geologic information and previous interpretations.

@article{uieda2012, title = {Robust 3D gravity gradient inversion by planting anomalous densities}, volume = {77}, issn = {00168033}, doi = {10.1190/geo2011-0388.1}, number = {4}, journal = {Geophysics}, author = {Uieda, Leonardo and Barbosa, Valéria C. F.}, year = {2012}, pages = {G55--G66}, }