THREE-DIMENSIONAL FINITE ELEMENT MODELING OF ELECTRICAL WELL LOGGING
LI Chang-Wei1,2, XIONG Bin1, Lü Yu-Zeng1
1. College of Earth Sciences, Guilin University of Technology, Guilin 541004, China;
2. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
Abstract:Mesh techniques were studied for the reason of large computational domain and model complexity in three-dimensional finite element modeling of electrical well logging. The authors employed a radial grid with parallel tri-prism to improve mesh quality and reduce the number of elements; Accurate element integration of parallel tri-prism was implemented by an affine transformation in order to shorten elements integral computational time; In combination with the unstructured mesh techniques, modeling under complicated condition was implemented. The computational efficiency and accuracy of the authors' schemes were examined by numerical experiments of several synthetic models. The modeling results under different situations are discussed.
[1] Coggon J H. Electromagnetic and electrical modeling by the finite element method [J].Geophysics,1971,36(2):132-151.[2] Sasaki Y. 3D Resistivity inversion using finite-element method [J]. Geophysics,1994,59(11):1839-1848.[3] Zhou B, Greenhalgh S A. Finite element three-dimensional direct current resistivity modelling: accuracy and efficiency considerations[J].Geophys J Int. 2001,145:679-688.[4] Pain C C, Herwanger J V, Worthington M H,et al. Effective multidimensional resistivity inversion using finite-element technique[J].Geophys J Int. 2002,151(3):710-728.[5] Li Y, Spitzer K. Three-dimensional DC resistivity forward modeling using finite elements in comparison with finite-difference solutions[J]. Geophys J Int,2002,151(3):924-934.[6] 阮百尧,熊彬,徐世浙.三维地电断面电阻率测深有限元数值模拟 [J].地球科学,2001,26(1):73-77.[7] 阮百尧,熊彬.电导率连续变化的三维电阻率测深有限元模拟 [J].地球物理学报,2002,45(1):131-138.[8] Zienkiewicz O C, Taylor R L, Zhu J Z. The finite element method its basis & fundamentals (sixth ed.)[M].Butterworth-Heinemann, Oxford,2005.[9] 岳建华,刘志新.井-地三维电阻率成像技术[J].地物理学进展,2005,20(2):407-411.[10] 王志刚,何展翔,刘昱.井地直流电法三维数值模拟及异常规律研究[J].工程地球物理学报,2006,3(2):87-92.[11] 徐凯军,李桐林.垂直有限线源三维地电场有限差分正演研究 [J].吉林大学学报(地球科学版), 2006, 36(1): 137-147.[12] 柯敢攀, 黄清华. 井地电法的三维正反演研究 [J]. 北京大学学报(自然科学版),2009,45(2): 264-272.[13] Loke M H. Topographic modeling in resistivity imaging inversion.in Proceedings of the 62nd EAGE Conference & Technical Exhibition(Extended Abstract,D-2),2000.[14] Kerry K, Weiss C. Adaptive finite element modeling using unstructed grids, the 2D magnetotelluric example[J].Geophysics,2006,71:G291-294.[15] Rücker C, Günther T,Spitzer K. Three-dimensional modelling and inversion of dc resistivity data incorporating topography - I. Modelling [J].Geophys J Int, 2006,166:495-505.[16] Blome M, Maurer H R, Schmidt K. Advances in three-dimensional geoelectric forward solver techniques [J].Geophys J Int.,2009,176:740-752.[17] 徐世浙.地球物理中的有限单元法 [M].北京:科学出版社,1994:178-188.[18] 《数学手册》编写组.数学手册 [M].北京:人民教育出版社,1979:990-998.[19] 吴小平,汪彤彤.利用共轭梯度算法的电阻率三维有限元正演 [J]. 地球物理学报,2003,46(3):428-432.[20] Saad Y. Iterative Methods for Sparse Linear Systems[M]. PWS, Bostion,1996.[21] Barrett R, Berry M, Chan T, et al. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods [M]. SIAM, Philadelphia,1994.