2D Inversion of magnetic and gravity data: a case study on Golgohar mine

Document Type : Research Paper


1 Faculty of Mining, Petroleum &Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran.

2 Institute of Geophysics University of Tehran, Tehran, Iran.


Because of the limitations of manipulating single geophysical data sets to interpret subsurface anomalies for many cases, it is required to combine geophysical data to decrease the ambiguity and non-uniqueness of the interpretation. Integration interpretation of two different geophysical data sets is one of the most common ways to integrate geophysical data and in this paper, we want to utilize the combination of gravity and magnetic data for the Golgohar mine in Iran. This mining case is located in the Sanandaj-Sirjan zone in the province of Kerman. Gravity and magnetic data are interpreted using a MATLAB code written based on the damped weighted minimum length solution for which the model weighting is the product of the multiplying of compactness and depth weighting constraints. At first, the inversion algorithm is applied to the synthetic case to investigate its reliability for practical application on real data. Reconstructed models from the noise-contaminated synthetic data are suggestive of the productivity of the inversion algorithm. Ultimately, the algorithm is applied for the interpretation of the real data and the inversion results of both data sets show a high correlation between the magnetite anomaly position horizontally and vertically. The results represent an anomaly with the depth ranging approximately from 25 to 130 m with the horizontal extension of about 120 m from 280 to 400 m relative to the start of the interesting profile.


Main Subjects

[1] Martinez, C., Li, Y., Krahenbuhl, R., & Braga, M. A. (2013). 3D inversion of airborne gravity gradiometry data in mineral exploration: A case study in the Quadrilátero Ferrífero, Brazil. Geophysics, 78(1), B1-B11
[2] Nettleton, L. L. (1976). Gravity and Magnetics in Oil Prospecting. McGraw-Hill Companies.
[3] Al-Farhan, M., Oskooi, B., Ardestani, V. E., Abedi, M., & Al-Khalidy, A. (2019). Magnetic and gravity signatures of the Kifl oil field in Iraq. Journal of Petroleum Science and Engineering, 183, 106397. 
[4] Paterson, N. R., & Reeves, C. V. (1985). Applications of gravity and magnetic surveys: The state-of-the-art in 1985. Geophysics, 50(12), 2558-2594.
[5] Hinze, W. J. (1990). The role of gravity and magnetic methods in engineering and environmental studies. In Geotechnical an Environmental Geophysics: Volume I: Review and Tutorial (pp. 75-126). Society of Exploration Geophysicists.
[6] Ward, S. H. (Ed.). (1990). Geotechnical and Environmental Geophysics: Volume I, Review and Tutorial: Volume I: Review and Tutorial. Society of Exploration Geophysicists.
[7] Sharma, P. V. (1987). Magnetic method applied to mineral exploration. Ore geology reviews, 2(4), 323-357.
[8] Cheng, Z. L., Zuo, Q. H., & Wang, J. F. (2013). Magnetic prospecting on magnet mining exploration area of Laos bokeo province. Progress in Geophysics, 3.
[9] Eppelbaum, L. V. (2011). Study of magnetic anomalies over archaeological targets in urban environments. Physics and Chemistry of the Earth, Parts A/B/C, 36(16), 1318-1330.
[10] Varfinezhad, R., Oskooi, B., & Fedi, M. (2020). Joint Inversion of DC Resistivity and Magnetic Data, Constrained by Cross Gradients, Compactness and Depth Weighting. Pure and Applied Geophysics, 1-19.
[11] Eventov, L. (1997). Applications of magnetic methods in oil and gas exploration. The Leading Edge, 16(5), 489-492.
[12] Zaher, M. A., Saibi, H., Mansour, K., Khalil, A., & Soliman, M. (2018). Geothermal exploration using airborne gravity and magnetic data at Siwa Oasis, Western Desert, Egypt. Renewable and Sustainable Energy Reviews, 82, 3824-3832.
[13] Oni, A. G., Eniola, P. J., Olorunfemi, M. O., Okunubi, M. O., & Osotuyi, G. A. (2020). The magnetic method as a tool in groundwater investigation in a basement complex terrain: Modomo Southwest Nigeria as a case study. Applied Water Science, 10(8), 1-18.
[14] Fregoso, E., & Gallardo, L. A. (2009). Cross gradients joint 3D inversion with applications to gravity and magnetic data. Geophysics, 74(4), 31–42.
[15] Joulidehsar, F., Moradzadeh, A., & Doulati Ardejani, F. (2018). An improved 3D joint inversion method of potential field data using cross-gradient constraint and LSQR method. Pure and Applied Geophysics, 175, 4389–4409
[16] Orlando, L. (2005). Joint interpretation of geophysical data for archaeology. A case study. Subsurface Sensing Technologies and Applications, 6(2), 235–250.
[17] Ogaya, X., Alcalde, J., Marza´n, I., Ledo, J., Queralt, P., Marcuello, A., et al. (2016). Joint interpretation of magnetotellurics, seismic, and well-log data in Hontomı´n (Spain). Solid Earth, 7, 1–15.
[18] Singh, A., Mishra, P. K., & Sharma, S. P. (2019). 2D cooperative inversion of direct current resistivity and gravity data: A case study of uranium-bearing target rock. Geophysical Prospecting, 67(3), 696-708
[19] Moradi, M., Oskooi, B., Pushkarev, P., Smirnov, M., & Oghaz, H. E. (2019). Cooperative inversion of magnetotelluric and seismic data on Shurab diapirs in Central Iran. Environmental Earth Sciences, 78(11), 1-14.
[20] Gallardo, L. A., & Meju, M. A. (2011). Structure coupled multiphysics imaging in geophysical sciences. Reviews of Geophysics, 49, RG1003.
[21] Varfinezhad, R., Hafizi, M. K., & Hashemi, H. (2015). Application of different inverse methods for combination of vS and vGPR data to estimate porosity and water saturation. Journal of the Earth and Space Physics, 41(4), 89-94.
[22] Wang, K. P., Tan, H. D., & Wang, T. (2017). 2D joint inversion of CSAMT and magnetic data based on cross-gradient theory. Applied Geophysics, 14(2), 279–290.
[23] Blakely, R. J. (1996). Potential Theory In Gravity And Magnetic Applications. Cambridge university press.
[24] Last, B. J., & Kubik, K. (1983). Compact gravity inversion. Geophysics, 48(6), 713-721.
[25] Guillen, A., & Menichetti, V. (1984). Gravity and magnetic inversion with minimization of a specific functional. Geophysics, 49(8), 1354-1360.
[26] Barbosa, V. C. F., & Silva, J. B. (1994). Generalized compact gravity inversion. Geophysics, 59(1), 57-68.
[27] Li, Y., & Oldenburg, D. W. (1998). 3-D inversion of gravity data. Geophysics, 63(1), 109-119.
[28] Ansari, A., Pasyar, A., Ghorbani, A. (2022). Solving 2-D gravity inversion problems using a PDE model in geophysics exploration. International Journal of Mining and Geo-Engineering, doi: 10.22059/ijmge.2021.317502.594888
[29] Boulanger, O., & Chouteau, M. (2001). Constraints in 3D gravity inversion. Geophysical prospecting, 49(2), 265-280.
[30] Cella, F., & Fedi, M. (2012). Inversion of potential field data using the structural index as weighting function rate decay. Geophysical Prospecting, 60(2), 313-336.
[31] Paoletti, V., Ialongo, S., Florio, G., Fedi, M., & Cella, F. (2013). Self-constrained inversion of potential fields. Geophysical Journal International, 195(2), 854-869.
[32] Tikhonov, A. N., & Arsenin, V. Y. (1977). Solutions of ill-posed problems. New York, 1-30.
[33] Menke, W. (2012). Geophysical data analysis: discrete inverse theory, MATLAB edn. New York: Academic.
[34] Li, Y., & Oldenburg, D. W. (1996). 3-D inversion of magnetic data. Geophysics, 61(2), 394-408.
[35] Fedi, M. (2007). DEXP: A fast method to determine the depth and the structural index of potential fields sources. Geophysics, 72(1), I1-I11. 
[36] R. Varfinezhad and V.E. Ardestani (2021). Kernel-based regularisation parameter and source-dependent depth weighting in gravity data inversion. Bulletin of Geophysics and Oceanography, pp. 249-266.
[37] Sahandi MR, Delavar ST, Sadeghi M, Jafari A, Moosavi A (2005).  Geological map of Iran (1:1,000,000) Geological survey of Iran.