Ilk_inv: a Matlab based algorithm for rapid computation of pseudo-3D density contrast distribution by using Bouguer gravity data

Abstract

The new generation Matlab-based algorithm provides a rapid estimation of density contrast distribution. The 3D assumption, which is based on the 1D equation, is used. Therefore, the output is called pseudo-3D instead of 3D. The algorithm uses singular value decomposition and the median filter to produce pseudo-3D results. The success of the method is tested by theoretical and field studies. For synthetic studies, ingle-source models produce reasonable outputs, compared to the true density contrast value. However, the multiple source model shows slight deviations which are ±0.3 g/cm3, with respect to the true density contrast value. The acceptable results are observed for the Bouguer anomaly of the eastern Mediterranean region. The resolution matrix indicates that the inversion process is biased due to the generalised inverse. The algorithm provides a quite different qualitative interpretation perspective to the interpreter.

Keywords

• singular value decomposition
• pseudo- 3d density map
• gravity modelling.