IGUG: A MATLAB package for 3D inversion of gravity data using graph theory

被引:10
|
作者
Vatankhah, Saeed [1 ]
Ardestani, Vahid Ebrahimzadeh [1 ]
Niri, Susan Soodmand [1 ]
Renaut, Rosemary Anne [2 ]
Kabirzadeh, Hojjat [3 ]
机构
[1] Univ Tehran, Inst Geophys, Tehran, Iran
[2] Arizona State Univ, Sch Math & Stat Sci, Tempe, AZ USA
[3] Univ Calgary, Dept Geomat Engn, Calgary, AB, Canada
关键词
Gravity; 3D inversion; Graph theory; Equidistance function; Mobrun;
D O I
10.1016/j.cageo.2019.03.008
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present an open source MATLAB package, IGUG, for 3D inversion of gravity data. The algorithm implemented in this package is based on methodology that was introduced by Bijani et al. (2015). A homogeneous subsurface body is modeled by an ensemble of simple point masses. The model parameters are the Cartesian coordinates of the point masses and their total mass. The set of point masses, assumed to each have the same mass, is associated to the vertices of a weighted complete graph in which the weights are computed by the Euclidean pairwise distances separating vertices. Kruskal's algorithm is used to solve the minimum spanning tree (MST) problem for the graph, yielding the reconstruction of the skeleton of the body described by the model parameters. The algorithm is stabilized using an equidistance function that restricts the spatial distribution of point masses and favors a homogeneous distribution for the subsurface structure. The non-linear global objective function for the model parameters comprises the data misfit term and the equidistance stabilization function. A regularization parameter lambda is introduced to balance the two terms of the objective function, and reasonable physically-relevant bound constraints are imposed on the model parameters. A genetic algorithm is used to minimize the bound constrained objective function for a fixed lambda, subject to the bound constraints. A new diagnostic approach is presented for determining a suitable choice for lambda, requiring a limited number of solutions for a small set of lambda. This contrasts the use of the L-curve which was suggested for estimating a suitable lambda in Bijani et al. (2015). Simulations for synthetic examples demonstrate the efficiency and effectiveness of the implementation of the algorithm. It is verified that the constraints on the model parameters are not restrictive, even with less realistic bounds acceptable approximations of the body are still obtained. Included in the package is the script GMD.m which is used for generating synthetic data and for putting measurement data in the format required for the inversion implemented within IGUG.m. The script Diagnostic_Results.m is included within IGUG.m for analyzing and visualizing the results, but can also be used as a standalone script given import of prior results. The software can be used to verify the simulations and the analysis of real data that is presented here. The real data set uses gravity data from the Mobrun ore body, north east of Noranda, Quebec, Canada.
引用
收藏
页码:19 / 29
页数:11
相关论文
共 50 条
  • [1] 3D inversion of borehole gravity data using cokriging
    GENG Meixia
    HUANG Danian
    XU Bowen
    Global Geology, 2014, 17 (04) : 225 - 230
  • [2] 3D Joint inversion of gravity and gravity tensor data
    ZHAO Simin
    GAO Xiuhe
    QIAO Zhongkun
    JIANG Dandan
    ZHOU Fei
    LIN Song
    Global Geology, 2018, 21 (01) : 55 - 61
  • [3] 3D Monte Carlo geometry inversion using gravity data
    Wei, Xiaolong
    Sun, Jiajia
    Sen, Mrinal
    GEOPHYSICS, 2024, 89 (03) : G29 - G44
  • [4] 3D stochastic inversion of gravity data using cokriging and cosimulation
    Shamsipour, Pejman
    Marcotte, Denis
    Chouteau, Michel
    Keating, Pierre
    GEOPHYSICS, 2010, 75 (01) : I1 - I10
  • [5] 3D inversion of airborne gravity-gradiometry data using cokriging
    Geng, Meixia
    Huang, Danian
    Yang, Qingjie
    Liu, Yinping
    GEOPHYSICS, 2014, 79 (04) : G37 - G47
  • [6] Fast 3D inversion of gravity data using Lanczos bidiagonalization method
    Reza Toushmalani
    Hakim Saibi
    Arabian Journal of Geosciences, 2015, 8 : 4969 - 4981
  • [7] Fast 3D inversion of gravity data using Lanczos bidiagonalization method
    Toushmalani, Reza
    Saibi, Hakim
    ARABIAN JOURNAL OF GEOSCIENCES, 2015, 8 (07) : 4969 - 4981
  • [8] 3D modelling of Trompsburg Complex (in South Africa) using 3D focusing inversion of gravity data
    Rezaie, Mohammad
    Moradzadeh, Ali
    Kalate, Ali Nejati
    Aghajani, Hamid
    Kahoo, Amin Roshandel
    Moazam, Sahar
    JOURNAL OF AFRICAN EARTH SCIENCES, 2017, 130 : 1 - 7
  • [9] 3D Gravity Inversion using Tikhonov Regularization
    Toushmalani, Reza
    Saibi, Hakim
    ACTA GEOPHYSICA, 2015, 63 (04) : 1044 - 1065
  • [10] 3D Gravity Inversion using Tikhonov Regularization
    Reza Toushmalani
    Hakim Saibi
    Acta Geophysica, 2015, 63 : 1044 - 1065