2.5D IP finite element numerical simulation in the time domain based on quadratic interpolation
Zhang Zhiyong1,2, Zhou Feng2, Li Zelin2,3
1. School of Geosciences and Info-Physics, Central South University, Changsha, Hunan 410083, China;
2. School of Nuclear Engineering and Geophysics, East China Institute of Technology, Nanchang, Jiangxi 330013, China;
3. Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan, Hubei 430074, China
Abstract:This paper presents an induced polarization numerical simulation in the time domain. We adopt the 2.5D finite element method based on the triangular element discrete and quadratic interpolated function. A cole-cole model is adopted to descript the complex resistivity of the discrete element. Then a boundary value problem of 2.5D complex resistivity is deduced. At last, a filtering algorithm based on Guptasarma is used to transform IP from the frequency-domain to the time domain. We design a fast algorithm suitable for 2.5D IP forward modeling in the time domain and a direct method for linear equations on symbolic analysis is used in this algorithm. The algorithm firstly processes approximately the boundary condition. So the obtained stiffness matrixes have only some relation with frequency and wave number, not with source point positions. Which the workload of matrix decomposition is reduced in the linear equations solution. Secondly, this algorithm achieves fast matrix LDLT decomposition with matrix rearrangement and element fill-in based on graph theory. The calculation scheme is optimized according to stiffness matrix and LDLT characteristics, so only one symbolic analysis is needed for the same domain discrete and only one LDLT decomposition is needed for all source points with the same frequency and wave number. Finally apparent resistivity, apparent polarization, and apparent dispersion rate of the synthetic model are analyzed.
张志勇, 周峰, 李泽林. 二次插值的时间域激电2.5维有限元数值模拟[J]. 石油地球物理勘探, 2015, 50(5): 999-1006.
Zhang Zhiyong, Zhou Feng, Li Zelin. 2.5D IP finite element numerical simulation in the time domain based on quadratic interpolation. OGP, 2015, 50(5): 999-1006.
Zhang Saizhen, Li Yingxian, Zhou Jiping et al. Induced polarization (IP) method in oil exploration—The cause of IP anomaly and it's relation to the oil reservoir. Chinese Journal of Geophysics,1986,29(6):597-614.
Liu Song. Effect of IP method for direct oil prospection and approaches for improving its application.Geological Science and Technology Information,1992,11(4):81-86.
[5]
王晓红.激电法在西北某油田上的应用效果.物探与化探,1996,20(1):64-70.
Wang Xiaohong. The effects of induced polarization method in a certain oil field of northwest China. Geophysical and Geochemical Exploration,1996,20(1):64-71.
[6]
郑秀芬,李金铭.油气电法勘探综述.地球物理学进展,1997,12(2):89-96.
Zheng Xiufen, Li Jinming. A review on the electrical exploration for oil-gas accumulation. Progress in Geophsics,1997,12(2):89-96.
Xu Chuanjian, Xu Zisheng, Yang Zhicheng et al. The effect of complex resistivity(CR) method to the exploration of oil and gas reservoirs.OGP,2004,39(S):31-35.
Li Zhengchen, Yan Weilin, Wang Wei. To identify oil and water layers in Daqing oilfield through complex resistivity well logging. World Well Logging Technology, 2007,22(1):23-26.
Zhang Hui, Wu Weijie, Liu Junli et al. The appliance prospect of complex electric resistance rate logging in land deposit oilfield. Fault-block Oil & Gas Field,2001,9(6):32-35.
Xiao Zhanshan,Xu Shizhe,Luo Yanzhong et al. Study on numeric simulation of complex resistivity logging.OGP,2007,42(3):343-347.
[14]
Snyder D D. A method for modeling the resistivity and IP response of two-dimensional bodies.Geophysics,1976,41(5):997-1015.
[15]
Oldenburg D W, Li Yaoguo. Inversion of induced polarization data. Geophysics,1994,59(9):1327-1341.
[16]
Weller A, Seichter M, Kampke A. Induced-polarization modeling using complex electrical conductivities.Geophys J Int,1996,127(2):387-398.
[17]
Routh P S, Oldenburg D W, Li Yaoguo. Regularized inversion of spectral IP parameters from complex resistivity data. SEG Technical Program Expanded Abstracts,1998,17:810-813.
[18]
Loke M H, Chambers J E, Ogilvy R D. Inversion of 2D spectral induced polarization imaging data. Geophysical Prospecting,2006,54(5):287-301.
[19]
Farias V J da C, Maranh C H de M, Rocha B R P da et al. Induced polarization forward modelling using finite element method and the fractal model. Applied Mathematical Modelling,2010,34(7):1849-1860.
Guptasarma D. Computation of the time-domain re-sponse of a polarizable ground.Geophysics,1982,47(11): 1574-1576.
[23]
Yuval, Oldenburg D W. Computation of Cole-Cole parameters from IP data. Geophysics,1997,62(2):436-448.
[24]
Hönig M, Tezkan B. 1D and 2D Cole-Cole-inversion of time-domain induced-polarization data. Geophysical Prospecting, 2007, 55(3):117-133.
[25]
Fiandaca G, Auken E, Christiansen A V et al.Time-domain-induced polarization: Full-decay forward modeling and 1D laterally constrained inversion of Cole-Cole parameters. Geophysics,2011,77(3):213-225.
Liu Bin, Li Shucai, Li Shuchen et al. 3D electrical resistivity inversion with least-squares method based on inequality constraint and its computation efficiency optimization. Chinese Journal of Geophysics,2012,55(1): 260-268.
Zhang Zhiyong, Liu Qingcheng. Study on 2D MT numerical simulation using FEM based on bi-tree Grid. OGP,2013,48(3):482-487.
[28]
徐世浙.地球物理中的有限单元法. 北京:科学出版社, 1994: 159-170.
[29]
Erdogan E, Demirci I, Candansayar M E. Incorporating topography into 2D resistivity modeling using finite-element and finite-difference approaches. Geophysics,2008,73(3): F135-F142.
[30]
刘崧.谱激电法.湖北武汉:中国地质大学出版社,1998.
[31]
Ellis R G, Oldenberg D W. The pole-pole 3D DC-resistivity inverse problem: a conjugate gradient approach. Geophys J Int,1994,119:187-194.
Zhang Jifeng, Tang Jingtian, Wang Ye et al. Magnetotellurics fast forward based on preconditioning conjugate gradient. Journal of Central South University (Science and Technology), 2010,41(5):1877-1882.