This study proposed the variational problem of 2.5D finite element forward modeling of the complex resistivity method and detailed the process of solving stiff matrix of finite equations. The Fractal model was introduced as a research model for studying the equivalent induced polarization anomalies of spectra. Furthermore, the complex conductivity and complex potential of a grid unit were linearly interpolated. Then, to obtain anomalous complex potential, finite element linear equations were solved using the biconjugate gradient stabilized method with incomplete LU decomposition. The results of three typical geoelectric models validated the correctness and accuracy of the algorithm proposed in this study. Furthermore, this study analyzed the abnormal response characteristics of 2.5D complex resistivity under different frequencies.
龙秀洁, 陈汉波, 莫亚军, 区小毅, 卢胜辉. 基于Fractal模型的复电阻率法2.5D有限元数值模拟[J]. 物探与化探, 2022, 46(4): 887-896.
Long Xiu-Jie, Chen Han-Bo, Mo Ya-Jun, Ou Xiao-Yi, Lu Sheng-Hui. Fractal model-based 2.5 D finite element modeling of complex resistivity method. Geophysical and Geochemical Exploration, 2022, 46(4): 887-896.
Yang Z W, Zhen W, Li X B, et al. The development and prospect of the spectral induced polarization method[J]. Geophysical and Geochemical Exploration, 2015, 39(1):22-28.
[2]
罗延钟, 张桂青. 频率域激电法原理[M]. 北京: 地质出版社, 1988.
[2]
Luo Y Z, Zhang G Q. Principle of frequency domain IP method[M]. Beijing: Geological Publishing House, 1988.
Zhang G Q, Cui X W, Luo Y Z. The Determination of intrinsic parameters by inversing IP apparent spectrum[J]. Geology and Exploration, 1987, 23(4) : 48-54.
Liu S, Guan S Y, Gao P F. Joint sip inversion for estimation of intrinsiccole-cole parameters of apolarizable ellipsoid[J]. Chinese Journal of Geophysice, 1994, 37(S1) : 542-551.
[6]
Hohmann G W. Three-dimensional induced polarization and electromagnetic modelling[J]. Geophysics, 1975, 40 (2) :309-324.
doi: 10.1190/1.1440527
[7]
Weller A, Seichter M, Kampke A. Induced-polarization modelling using complex electrical conductivities[J]. Geophys. J. Int., 1996, 127(1):387-398.
doi: 10.1111/j.1365-246X.1996.tb04728.x
[8]
Oldenburg D W, Li Y. Inversion of induced polarization data[J]. Geophysics, 1994, 59 (9) :1327-1341.
doi: 10.1190/1.1443692
[9]
Rocha B R P, Habashy T M. Fractal geometry, porosity and complex resistivity. I: From rough pore interfaces to hand specimens, in: Developments in Petrophysics[M], London: Geological Society Publishing House, London, 1995.
Liu J X, Jiang P F, Tong X Z, et al. Application of BICGSTAB algorithm with incomplete LU decomposition preconditioning to two-dimensional magnetotelluric forward modeling[J]. Journal of Central South University: Science and Technology, 2009, 40(2):484 -491.