RESIP2DMODE: A MATLAB-Based 2D Resistivity and Induced Polarization Forward Modeling Software

نوع مقاله: مقاله تحقیقی‌ (پژوهشی‌)

نویسندگان

1 Assistant Professor, Institute of Geophysics, University of Tehran, Tehran, Iran

2 M.Sc. Student, Institute of Geophysics, University of Tehran, Tehran, Iran

چکیده

Forward modeling is an integral part of every geophysical modeling resulting in the numerical simulation of responses for a given physical property model. This Forward procedure is helpful in geophysics both as a means to interpret data in a research setting and as a means to enhance physical understanding in an educational setting. Calculation of resistivity and induced polarization forward responses is carried out using simulation of the current flow into the earth’s surface through solving the Poisson’s equation. In this contribution, a finite-difference algorithm is applied to discretize the simulated models restricted by a mixed boundary condition. To account for the 3D source characteristic, a spatial Fourier transform of the partial differential equations with respect to a range of wave numbers is performed along the strike direction. Then, an inverse Fourier transformation is conducted to obtain the potential solutions in the spatial domain. The present package provides a user-friendly interface designed to understand and handle for various conventional electrical configurations in the frame of the MATLAB programming language. To verify the program, initial responses of some simple models are compared with those of analytic solutions, which proved satisfactory in terms of accuracy. For further evaluation, the code is also examined on some complicated models.
 

کلیدواژه‌ها


عنوان مقاله [English]

RESIP2DMODE: A MATLAB-Based 2D Resistivity and Induced Polarization Forward Modeling Software

نویسندگان [English]

  • Reza Ghanati 1
  • Yosra Azadi 2
  • Razieh Fakhimi 2
1 Assistant Professor, Institute of Geophysics, University of Tehran, Tehran, Iran
2 M.Sc. Student, Institute of Geophysics, University of Tehran, Tehran, Iran
چکیده [English]

Forward modeling is an integral part of every geophysical modeling resulting in the numerical simulation of responses for a given physical property model. This Forward procedure is helpful in geophysics both as a means to interpret data in a research setting and as a means to enhance physical understanding in an educational setting. Calculation of resistivity and induced polarization forward responses is carried out using simulation of the current flow into the earth’s surface through solving the Poisson’s equation. In this contribution, a finite-difference algorithm is applied to discretize the simulated models restricted by a mixed boundary condition. To account for the 3D source characteristic, a spatial Fourier transform of the partial differential equations with respect to a range of wave numbers is performed along the strike direction. Then, an inverse Fourier transformation is conducted to obtain the potential solutions in the spatial domain. The present package provides a user-friendly interface designed to understand and handle for various conventional electrical configurations in the frame of the MATLAB programming language. To verify the program, initial responses of some simple models are compared with those of analytic solutions, which proved satisfactory in terms of accuracy. For further evaluation, the code is also examined on some complicated models.
 

کلیدواژه‌ها [English]

  • RESIP2DMODE
  • Resistivity
  • Chargeability
  • Forward modeling
  • MATLAB
  • finite-difference
Akca, I., 2016, LRIS2D: a MATLAB package for the 2D inversion of DC resistivity/IP data. Acta Geophysica, 64(2), 443–462.

Befus, K. M., 2017, pyres: A Python Wrapper for Electrical Resistivity Modeling with R2, Journal of Geophysics and Engineering, 15, 338–346.

Dahlin, T., Rosqvist, H., Leroux, V., 2010. Resistivity-IP for landfill applications. First Break, 28(8), 101–105.

Demkowicz, L., 2007, Computing with Hp-Adaptive Finite Elements. Chapman and Hall/CRC.

Dey, A. and Morrison, H. F., 1979a, Resistivity modeling for arbitrarily shaped three-dimensional structures. Geophysics, 44(4), 753–780.

Gazoty, A., Fiandaca, G., Pedersen, J., Auken, E., Christiansen, A.V., 2012, Mapping of land-fills using time-domain spectral induced polarization data: the Eskelund case study. Near Surf. Geophys., 10, 575–586.

Günther, T., 2004, Inversion Methods and Resolution Analysis for the 2D/3D Reconstruction of Resistivity Structures from DC Measurements. PhD thesis, Freiberg University of Mining and Technology.

Günther, T., Rücker, C., and Spitzer, K., 2006, 3-d modeling and inversion of DC resistivity data incorporating topography - Part II: Inversion. Geophys. J. Int., 166(2):, 506-517.

Hughes, T. J. R., 1987, The Finite Element Method. Linear Static and Dynamic Finite Element Analysis. Prentice Hall, Englewood Cliffs, NJ.

Loke, M. H, 2014, RES2DMOD ver. 3.01: Rapid 2D resistivity forward modeling using the finite-difference and finite-element methods, Geotomo Software. https://www.geotomosoft.com.

Lowry, T., Allen, M. B., and Shive, P. N., 1989, Singularity removal: A refinement of resistivity modeling techniques. Geophysics, 54(6), 766–774.

McGillivray, P. R., 1992, Forward modeling and inversion of DC resistivity and MMR data. Ph.D. thesis: University of British Columbia.

Mendez-Delgado, S., Gomez-Trevino, E., and Perez-Flores, M. A., 1999, Forward modelling of direct current and low-frequency elec-tromagnetic fields using integral equations: Geophysical Journal International, 137, 336–352.

Pidlisecky, A., Knight, R., 2008, FW2_5D: A MATLAB 2.5-D electrical resistivity modeling code. Computers & Geosciences, 34, 1645–1654.

Qiang, J. K., Han, X., and Dai, S.K., 2013, 3-D DC resistivity inversion with topography
based on regularized conjugate gradient method. International Journal of Geophysics, Article ID 931876.

Rücker, C., Günther, T., and Spitzer, K., 2006, Three-dimensional modelling and inversion of dc resistivity data incorporating Topography-I. Modelling. Geophysical Journal International, 166, 495–505.

Rücker, C., Günther, T., Wagner, F.M., 2017, pyGIMLi: an open-source library for mod-elling and inversion in geophysics. Comput. Geosci., 109, 106–123.

Spitzer, K., and Wurmstich, B., 1995, Speed and accuracy in 3d resistivity modeling. In 3D EM Symposium, Ridgefield, Connecticut.

Telford, W., Geldart, L., Sheriff, R., Keys, D., 1991, Applied Geophysics. Cambridge University Press.

Wu, X. P., 2003, A 3-D finite-element algorithm for DC resistivity modeling using the shifted incomplete Cholesky conjugate gradient method. Geophysical Journal International, 154, 947–956.

Yuan, Y., J. Qiang, J. Tang, Z., Ren, and Xiao, X., 2016, 2.5D direct-current resistivity forward modelling and inversion by finite-element–infinite-element coupled method. Geophysical Prospecting, 64, 767–779.

Zhang, J., Mackie, R. L., and Madden, T. R., 1995, 3-d resistivity forward modeling and inversion using conjugate gradients. Geophysics, 60(5), 1313–1325.