Solving 2-D gravity inversion problems using a PDE model in geophysics exploration

Document Type : Research Paper


Faculty of Mining and Metallurgical Engineering, Yazd University, Yazd, Iran


Inverse modeling is one of the useful solutions to create a logical model with relationships between observed and measured values. In geophysical and subsurface investigations such as cavities or mineral explorations, solving inverse problems using problem physics in a partial differential equation (PDE) system is very important. In this research, COMSOL multiphysics’ optimization interface, combined with a PDE or physics interface, was used to solve inverse-modeling problems. Also, a framework is presented to solve undetermined inverse problems using COMSOL multiphysics’ optimization. COMSOL multiphysics does not include a gravity calculation module. However, since Poisson’s equation governs gravity and electrostatics, a gravity model can be created in the electrostatics module by changing the electrical permittivity value. We present a general adjoint state formulation that may be used in this framework and allows for faster calculation of sensitivity matrices in a variety of commonly encountered underdetermined problems. First of all, 2D inversion of gravity data has been run and validated in COMSOL multiphysics software using one synthetic model and synthetic data in a forward modeling process. Afterward, using real gravity data surveyed along a cross-section of the sinkholes in the NW of Abarkuh, the lateral structure and subsurface cavities were estimated. The inverted gravitational acceleration values, then cross-correlated with observed gravity data and available surface pieces of evidence such as sinkholes and circular structures. The results indicated that our COMSOL-based routines for the solution of PDE-based inverse problems using adjoint states, while high in computational speed, can be used in modeling a wide range of physical systems governed by the partial differential equation laws and also can accurately discriminate between low-density contrast regions and background.


[1]  Zhang, J., Wang, C.-Y., Shi, Y., Cai, Y., Chi, W.-C., Dreger, D., Cheng, W. B., and Yuan, Y.H., 2004. Three-dimensional crustal structure in central Taiwan from gravity inversion with a parallel genetic algorithm. Geophysics, 69, 917–924. 2004.
[2]  Farquharson, C.G., Mosher, C.R.W., 2009. Three-dimensional modeling of gravity data using finite differences. Journal of Applied Geophysics, 68, 417–422.
[3]  COMSOL Multiphysics, User’s Guide, 2017. Version 5.2, COMSOL AB, Stockholm, Sweden.
[4]  Butler, S.L., Sinha, G., 2012. Forward modeling of applied geophysics methods using Comsol and comparison with analytical and laboratory analog models. Computers & Geosciences, 42, 168–176.
[5]  Cardiff, M., Kitanidis, P.K., 2008. Efficient solution of nonlinear, underdetermined inverse problems with a generalized PDE model. Computers & Geosciences, 34, 1480–1491.
[6]  Kalavagunta, A., Weller, R.A., 2005. Accurate geometry factor estimation for the four-point probe method using COMSOL multiphysics. In: Proceedings of the Comsol Users Conference, Boston.
[7]  Park, J., Bjornara, T.I., Farrelly, B.A., 2010. Absorbing boundary domain for CSEM 3D modeling. In: Excerpt from the Proceedings of the COMSOL Conference 2010, Paris.
[8]  Volkmann, J., Mohnke, O., Klitzsh, N., Blaschek, R., 2008. Microscale modeling of the frequency-dependent resistivity of porous media. In: Proceedings of the COMSOL Conference 2008, Hanover.
[9]  Stoll, J.B., 2005. FE-modelling of electrical borehole tool responses. In: Proceedings of the COMSOL Multiphysics User’s Conference 2005, Frankfurt.
[10] Braun, M., Rommel, I., Yaramanci, U., 2005. Modeling of magnetic resonance sounding using finite elements (FEMLAB) for 2D resistivity extension. In: Proceedings of the COMSOL Multiphysics User’s Conference 2005, Frankfurt.
[11] Telford, W.M., Geldart, L.P., Sheriff, RE., 1990. Applied Geophysics, 2nd ed. Cambridge University Press, Cambridge.
[12] Jackson, J.D., 1998. Classical Electrodynamics, Wiley, New York.
[13] Kitanidis, P.K., 1995. Quasi-linear Geostatistical Theory for Inversing. Water Resources Research, vol. 31, no. 10, pp., 2411–2419.
[14] Madsen, K, Nielsen, H.B., and Tingleff, , 2004. Methods for Non-Linear Least Squares Problems, 2nd ed. Informatics and Mathematical Modelling (IMM), Technical University of Denmark (DTU).