Spectral induced polarization forward modeling in the framework of finite difference scheme

Document Type : Research Paper


Institute of Geophysics, University of Tehran, Tehran, Iran



Dealing with numerous reviews and widespread inquiries, it has been concluded that much more information and interpretive parameters are accessible regarding the subsurface structures when using a particular frequency range in the spectral induced polarization (SIP) measurements. Therefore, the interpretation uncertainty would diminish which causes studies with more valid and authentic outcomes. This could be achieved by using a comprehensive and general model which is appropriate for representing electrical features variation in terms of frequency, known as the Cole-Cole model. By using the SIP method and applying a defined broad of frequencies, it would be conceivable to describe items such as medium properties, spectral behavior of the studied area, and the intensity of each single parameter. The widespread use of the SIP method requires accurate and fast modeling and inversion algorithms. An integral part of every geo-electrical data inversion is an accurate and efficient forward modeling resulting in numerical simulation of responses for a given physical property model. In other words, like every other geophysical method, a reliable spectral-induced polarization inversion is highly dependent on the accuracy of the forward problem. Forward modeling is accomplished over a 2D earth structure to generate complex resistivity data by simulating current flow into the earth's surface and solving the Poisson equation containing complex values. In this contribution, a finite difference algorithm is applied to solve the complex partial differential equations (PDEs) restricted by a mixed boundary condition. A spatial Fourier transform of the PDEs, with respect to a defined range of wavenumbers, is carried out along the strike direction to elucidate 3D source characteristics. Eventually, it is necessary to conduct an inverse Fourier transform to obtain potential solutions in the spatial domain. To verify the accuracy of the proposed numerical algorithm, some synthetic models are simulated and the forward responses, including resistance and phase values with respect to a specific frequency spectrum, are calculated. Furthermore, a comparison between our numerical results and those of Geotomo geo-electrical software is provided.


Main Subjects

[1] Pelton, W. H, Rijo, L, Swift, J. R, 1978, Inversion of two-dimensional resistivity and induced polarization data. Geophysics, 681-904.
[2] Olsson, P.-I., Dahlin, T., Fiandaca, G., Auken, E., 2015. Measuring time-domain spectral induced polarization in the on-time: decreasing acquisition time and increasing signalto- noise ratio. J. Appl. Geophys. 123, 316–321.
[3] Revil, A., Karaoulis, M., Johnson, T., and Kemna., A., 2012, Review:Some low-frequency electrical methods for subsurface characterization and monitoring in hydrogeology: Hydrogeology Journal, 20(4), 617–658.
[4] Dey, A., Morrison, H.F., 1979, Resistivity modeling for arbitrary shaped two-dimensional structures, Geophysical Prospecting 27, 1020–1036.
[5] 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
[6] Lowery, T, Allen, M, B, Shive, P, N, 1989. Singularity removal: a refinement of resistivity modelling techniques. Gepphysics, 680-789.
[7] Weller A., Seichter M. and Kampke A. 1996. Induced-polarization modelling using complex electrical conductivities. Geophysical Journal International 127, 387-398
[8] Routh, P.S., Oldenburg, D.W., Li, Y., 1998. Regularized inversion of spectral IP parameters from complex resistivity data. In: 68th SEG Meeting, New Orleans, USA, Expanded Abstracts, pp. 810–813.
[9] Loke, M.H., Chambers, J.E., 2006, Inversion of 2D spectral induced polarization imaging data, Geophysical Prospecting, 54, 287–301.
[10] Kemna, A., Binley, A., Ramirez, A. L., Daily, W. D., 2000, Complex resistivity tomography for environmental applications.
[11] Rucker, C, Gunther, T, Spitzer, K, 2006, Three dimensional inversion of DC resistivity data incorporating topography. Geophysics Journal 2006.
[12] Xu, K. J., 2007, Study on 2.5D complex resistivity electromagnetic forward and in version: PhD Thesis, Jilin University.
[13] Farias, V. J., Maranhao, C. H., Andrade, P, O, Induced polarization forward modelling using finite element method and the fractal model, Elsevier. 2007.
[14] Rucker, C, Gunther, T, Spitzer, K, 2011, Three dimensional modelling of DC resistivity data incorporating topography. Geophysics Journal 2011.
[15] Zhang, Z. Y., Tan, H. D., Wang, K. P., Lin, C. H., Zhang,B. and Xie, M. B., 2016, Two-dimensional inversion of spectral induced polarization data using MPI parallel.
[16] Günther, T., Martin, Tina, 2016, Spectral two-dimensional inversion of frequency-domain induced polarization data, Journal of Applied Geophysics 135, 436–448.
[17] Fiandaca, G., Ramm, J., 2013, Resolving spectral information from time domain induced polarization data through 2-D inversion, Geophys. J. Int 192, 631–646.
[18] Kemna, Andres., 2004, Tomographic inversion of complex resistivity, theory and application.
[19] Ghanati, R., Azadi, Y., Fakhimi, R., 2020. RESIP2DMODE: A MATLAB-Based 2D Resistivity and Induced Polarization Forward Modeling Software, Iranian Journal of Geophysics 13 (4), 60-78.