3D Inversion of Magnetic Data through Wavelet based Regularization Method

Document Type : Research Paper


1 Department of Mining Engineering, College of Engineering, University of Tehran, Iran

2 Institute of Geophysics, University of Tehran, Iran


This study deals with the 3D recovering of magnetic susceptibility model by incorporating the sparsity-based constraints in the inversion algorithm. For this purpose, the area under prospect was divided into a large number of rectangular prisms in a mesh with unknown susceptibilities. Tikhonov cost functions with two sparsity functions were used to recover the smooth parts as well as the sharp boundaries of model parameters. A pre-selected basis namely wavelet can recover the region of smooth behaviour of susceptibility distribution while Haar or finite-difference (FD) domains yield a solution with rough boundaries. Therefore, a regularizer function which can benefit from the advantages of both wavelets and Haar/FD operators in representation of the 3D magnetic susceptibility distributionwas chosen as a candidate for modeling magnetic anomalies. The optimum wavelet and parameter β which controls the weight of the two sparsifying operators were also considered. The algorithm assumed that there was no remanent magnetization and observed that magnetometry data represent only induced magnetization effect. The proposed approach is applied to a noise-corrupted synthetic data in order to demonstrate its suitability for 3D inversion of magnetic data. On obtaining satisfactory results, a case study pertaining to the ground based measurement of magnetic anomaly over a porphyry-Cu deposit located in Kerman providence of Iran. Now Chun deposit was presented to be 3D inverted. The low susceptibility in the constructed model coincides with the known location of copper ore mineralization.


Abedi, M., Gholami, A., Norouzi, G.H. and Fathianpour, N. (2013). “Fast inversion of magnetic data using Lanczos bidiagonalization method.” Journal of Applied Geophysics, Vol. 90, PP. 126-137.
[2] Abedi, M., Gholami, A. and Norouzi, G.H. (2014). “3D inversion of magnetic data seeking sharp boundaries: a case study for a porphyry copper deposit from Now Chun in central Iran.” Near Surface Geophysics, Vol. 12, PP. 657-666.
[3] Gholami, A. and Siahkoohi, H.R. (2010). “Regularization of linear and non-linear geophysical ill-posed problems with joint sparsity constraints.” Geophysical Journal International, Vol. 180, No. 2, PP. 871-882.
[4] Tikhonov, A.N. and Arsenin, V.Y. (1977). “Solutions of Ill-Posed Problems.” Winston, Washington, D.C.
[5] Li. Y. and Oldenburg, D.W. (1996). “3-D inversion of magnetic data.” Geophysics, Vol. 61, PP. 394–408.
[6] Pilkington, M. (2009). “3D magnetic data-space inversion with sparseness constraints.” Geophysics, Vol. 74, PP. L7-L15.
[7] Li. Y. and Oldenburg, D.W. (1998). “3-D inversion of gravity data.” Geophysics, Vol. 63, PP. 109-119.
[8] Boulanger, O. and Chouteau, M. (2001). “Constraints in 3D gravity inversion.” Geophysical Prospecting, Vol. 49, No. 2, PP. 265– 280.
[9] Portniaguine, O. and Zhdanov, M.S. (2002). “3-D magnetic inversion with data compression and image focusing.” Geophysics, Vol. 67, PP. 1532–1541.
[10] Chasseriau, P. and Chouteau, M. (2003). “3D gravity inversion using a model of parameter covariance.” Journal of Applied Geophysics, Vol. 52, No. 1, PP. 59-74.
[11] Caratori Tontini, F., Cocchi, L. and Carmisciano, C. (2006). “Depth-to-the-bottom optimization for magnetic data inversion: Magnetic structure of the Latium volcanic region, Italy.” Journal of Geophysical Research, Vol. 111, PP. 1-17.
[12] Pignatelli, A., Nicolosi, I. and Chiappini, M. (2006). “An alternative 3D source inversion method for magnetic anomalies with depth resolution.” Annals of Geophysics, Vol. 49, PP. 1021-1027.
[13] Malehmir, A., Thunehed, H. and Tryggvason, A. (2009). “A Case History: the Paleoproterozoic Kristineberg mining area, northern Sweden: Results from integrated 3D geophysical and geologic modeling, and implications for targeting ore deposits.” Geophysics, Vol. 74, PP. B9-B22.
[14] Namaki, L., Gholami, A. and Hafizi, M.K. (2011). “Edge-preserved 2-D inversion of magnetic data: an application to the Makran arc-trench complex.” Geophysical Journal International, Vol. 184, No. 3, PP. 1058-1068.
[15] Shamsipour, P., Chouteau, M. and Marcotte, D. (2011). “3D stochastic inversion of magnetic data.” Journal of Applied Geophysics, Vol. 73, No. 4, PP. 336-347.
Abedi et al. / Int. J. Min. & Geo-Eng., Vol.49, No.1, June 2015
[16] Cella, F. and Fedi, M. (2012). “Inversion of potential field data using the structural index as weighting function rate decay.” Geophysical Prospecting, Vol. 60, No. 2, PP. 313-336.
[17] Li, Y. and Oldenburg, D.W. (2003). “Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method.” Geophysical Journal International, Vol. 152, No. 2, PP. 251-265.
[18] Charbonnier, P., Blanc-Feraud, L., Aubert, G. and Barlaud, M. (1997). “Deterministic edge-reserving regularization in computed imaging.” IEEE Transactions on Image Processing, Vol. 6, PP. 298–310.
[19] Geman, S. and McClure, D.E. (1985). “Bayesian image analysis: an application to single photon emission tomography, in Proc. Statistical Computation Section.” PP. 12–18, Amer. Statistical Assoc., Washington, DC.
[20] Sacchi, M.D. and Ulrych, T.J. (1995). “High resolution velocity gathers and offset space reconstruction.” Geophysics, Vol. 60, PP. 1169-1177.
[21] Bertete-Aguirre, H., Cherkaev, E. and Oristaglio, M. (2002). “Non-smooth gravity problem with total variation penalization functional.” Geophysical Journal International, Vol. 149, No. 2, PP. 499–507.
[22] Farquharson, C.G. (2008). “Constructing piecewise-constant models in multidimensional minimum-structure inversions.” Geophysics, Vol. 73, No. 1, PP. K1-K9. [23] Gholami, A. and Siahkoohi, H.R. (2009). “Simultaneous constraining of model and data smoothness for regularization of geophysical inverse problems.” Geophysical Journal International, Vol. 176, No. 1, PP. 151–163.
[24] Abedi, M., Gholami, A. and Norouzi, G.H. (2013). “A stable downward continuation of airborne magnetic data: A case study for mineral prospectivity mapping in Central Iran.” Computers & Geosciences, Vol. 52, PP. 269–280.
[25] Abedi, M., Mosazadeh, K., Dehghani, H. and MadanchiZare, A. (2014). “Enhancing magnetic signals in unexploded ordnances (UXO) detection based on edge-preserved stable downward continuation method.” Journal of Mining & Environment, Vol. 5, No. 1, PP. 13-24.
[26] Routh, P.S., Qu, L., Sen, M.K. and Anno, P.D. (2007). “Inversion for nonsmooth models with physical bounds.” SEG, Expanded Abstracts, Vol. 26, No. 1, PP. 1795–1799.
[27] Donoho, D.L. (1995). “Nonlinear solution of linear inverse problems by wavelet-vaguelette decomposition.” Applied and Computational Harmonic Analysis, Vol. 2, No. 2, PP. 101–126.
[28] Donoho, D.L. and Johnsstone, I.M. (1998). “Minmax estimation via wavelet shrinkage.” Annals of Statistics, Vol. 26, PP. 879-921.
[29] Daubechies, I., Defriese, M. and De Mol, C. (2004). “An iterative thresholding algorithm for linear inverse problems with a sparsity constraint.” Communications on Pure and Applied Mathematics, Vol. LVII, PP. 1413–1457.
[30] Loris, I., Nolet, G., Daubechies, I. and Dahlen, F.A. (2007). “Tomographic inversion using l1-norm regularization of wavelet coefficients.” Geophysical Journal International, Vol. 170, PP. 359–370. [31] Lelièvre, P.G. and Oldenburg, D.W. (2009). “A 3D total magnetization inversion applicable when significant, complicated remanence is present.” Geophysics, Vol. 74, No. 3, PP. L21-L30.
[32] Li, Y., Shearer, S.E., Haney, M.M. and Dannemiller, N. (2010). “Comperhensive approaches to 3D inversion of magnetic data affected by remanent magnetization.” Geophysics, Vol. 75, PP. L1-L11.
[33] Bhattacharyya, B.K. (1964). “Magnetic anomalies due to prism-shaped bodies with arbitrary polarization.” Geophysics, Vol. 29, PP. 517–53.
[34] Rao, D.B. and Babu, N.R. (1991). “A rapid method for three-dimensional modeling of magnetic anomalies.” Geophysics, Vol. 56, PP. 1729-1737.
[35] Daubechies, I. (1988). “Orthonormal bases of compactly supported wavelets.” Communications on Pure and Applied Mathematics, Vol. 41, PP. 909-996.
[36] Aster, R.C., Borchers, B. and Thurber, C. (2003). “Parameter Estimation and Inverse Problems.” Academic Press, New York, NY.
[37] Bregman, L. (1967). “The relaxation method of finding the common points of convex sets and its application to the solution of problems in convex optimization.” USSR Computational Mathematics and Mathematical Physics, Vol. 7, PP. 200–217.
[38] Donoho, D.L. (1995). “De-noising by soft thresholding.” IEEE Transactions on Information Theory, Vol. 41, PP. 613–627.
Abedi et al. / Int. J. Min. & Geo-Eng., Vol.49, No.1, June 2015
[39] Daubechies, I. (1992). “Ten Lectures on Wavelets”. SIAM, Philadelphia.
[40] Mallat, S.G. (1998). “A Wavelet Tour of Signal Processing”. Academic Press, San Diego.
[41] Li, Y. and Oldenburg, D.W. (2000). “Joint inversion of surface and three-component borehole magnetic data.” Geophysics, Vol. 65, PP. 540-552.
[42] Elyasi, G.R. (2009). “Mineral Potential Mapping in Detailed Stage Using GIS in One of Exploration Prospects of Kerman Province.” Master of Science Thesis, University of Tehran (published in Persian).
[43] Clark, D.A. (1999). “Magnetic petrology of igneous intrusions-Implications for exploration and magnetic interpretation.” Exploration Geophysics, Vol. 20, PP. 5–26.
[44] John, D.A., Ayuso, R.A., Barton, M.D., Blakely, R.J., Bodnar, R.J., Dilles, J.H., Gray, Floyd, Graybeal, F.T., Mars, J.C., McPhee, D.K., Seal, R.R., Taylor, R.D. and Vikre, P.G. (2010). “Porphyry copper deposit model, chap. B of Mineral deposit models for resource assessment.” U.S. Geological Survey Scientific Investigations Report 2010–5070–B, PP. 169.
[45] Thoman, M.W., Zonge, K.L. and Liu, D. (2000). “Geophysical case history of North Silver Bell, Pima County, Arizona—A supergene-enriched porphyry copper deposit, in Ellis RB, Irvine R, Fritz F, eds., Northwest Mining Association 1998 Practical Geophysics Short Course Selected Papers on CD-ROM: Spokane, Washington. Northwest Mining Association, paper 4, PP. 42.
[46] Oldenburg, D.W., Li, Y. and Ellis, R.G. (1997). “Inversion of geophysical data over a copper gold porphyry deposit: A case history for Mt. Milligan.” Geophysics, Vol. 62, PP. 1419-1431.
[47] Telford, W.M., Geldart, L.P. and Sheriff, R.E. (2004). “Applied Geophysics”. Second Edition, Cambridge University Press.
[48] Hezarkhani, A. (2009). “Hydrothermal fluid geochemistry at the Chah-Firuzeh porphyry copper deposit, Iran, evidence from fluid inclusions.” Journal of Geochemical Exploration, Vol. 101, No. 3, 254–264