Finite element 2D forward modeling of DC method with line source

Authors

XIE Haijun, College of Geology and Environment, Xian University of Science and Technology, Xian 710054, China
LI Zhiqiang, College of Geology and Environment, Xian University of Science and Technology, Xian 710054, China
LI Sheng, College of Geology and Environment, Xian University of Science and Technology, Xian 710054, China

Abstract

In direct current prospecting, the forward calculation of 3D point source is complex and the calculation quantity is large, and using the Fourier transform can converted into 2D point source, but it is more complicated than that of 2D line source. This paper is based on Matlab programming, using the first boundary condition, and gives numerical solution of two-dimensional finite element of DC method with line source, through comparing with theoretical value of homogeneous half-space, the correctness of calculation results was verified. Then the high and low resistivity geoelectric models were built, the response of single side power supply and double side power supply were compared and analyzed respectively, and the corresponding apparent resistivity profile was obtained. The results show that unilateral power supply can accurately reflect the vertical position of the abnormal body, while bilateral power supply can more accurately reflect the horizontal position and the shape of the abnormal body, providing theoretical basis for actual production.

Keywords

DC method with line source, the first boundary condition, 2D finite element, forward modeling, unilateral(bilateral) power supply

DOI

10.3969/j.issn.1001-1986.2019.01.030

Recommended Citation

XIE Haijun, LI Zhiqiang, LI Sheng, et al. (2019) "Finite element 2D forward modeling of DC method with line source," Coal Geology & Exploration: Vol. 47: Iss. 1, Article 31.
DOI: 10.3969/j.issn.1001-1986.2019.01.030
Available at: https://cge.researchcommons.org/journal/vol47/iss1/31

Reference

[1] 傅良魁. 电法勘探教程[M]. 北京:地质出版社,1983.

[2] 陈明生,严又生. 二维水平电偶极变频测深阻抗视电阻率的有限元正演计算[J]. 地球物理学报, 1987, 30(2):201-208. CHEN Mingsheng,YAN Yousheng. Finite element forward calculation of apparent resistivity of two-dimensional horizontal electric dipole frequency conversion sounding impedance[J]. Chinese Journal of Geophysics,1987,30(2):201-208.

[3] 汤井田,肖晓,杜华坤,等. ANSYS在直流电法正演中的应用[J]. 地球物理学进展,2006,21(3):987-992. TANG Jingtian,XIAO Xiao,DU Huakun,et al. Application of ANSYS in direct current method forward modeling[J]. Progress in Geophysics,2006,21(3):987-992.

[4] 王小龙,冯宏,田华光,等. 基于COMSOL MULTIPHYSICS的直流电法正演模拟[J]. 煤田地质与勘探,2011,39(5):76-80. WANG Xiaolong,FENG Hong,TIAN Huaguang,et al. Direct current forward modeling based on Comsol Multiphysics[J]. Coal Geology & Exploration,2011,39(5):76-80.

[5] 何继善,周正秀. 电阻率法地形改正中的点电源场与线电源场问题[J]. 地质与勘探,1975,11(8):45-51. HE Jishan,ZHOU Zhengxiu. Problems of point power field and line power field in resistivity topography correction[J]. Geology and Exploration,1975,11(8):45-51.

[6] 刘树才,刘志新,姜志海,等. 矿井直流电法三维正演计算的若干问题[J]. 物探与化探,2004,28(2):170-172. LIU Shucai,LIU Zhixin,JIANG Zhihai,et al. Some problems of three dimensional forward calculation of mine direct current method[J]. Geophysical and Chemical Exploration,2004,28(2):170-172.

[7] 彭艳华,柳建新,刘海飞,等. 直流电阻率2.5维有限元数值模拟[J]. 地球物理学进展,2014,29(4):1926-1931. PENG Yanhua,LIU Jianxin,LIU Haifei,et al. Numerical simulation of DC resistivity 2.5 dimensional finite element[J]. Progress in Geophysics,2014,29(4):1926-1931.

[8] 蒋鑫. 高密度电阻率法的正演数值模拟研究[D]. 长春:吉林大学,2013.

[9] 金聪,刘江平. 二维高密度电阻率法数值模拟与应用[J]. 地质与勘探,2014,50(5):984-990. JIN Cong,LIU Jiangping. Two-dimensional numerical simulation and application of the high-density resistivity method[J]. Geology & Exploration,2014,50(5):984-990.

[10] 严波,韩波. 二维直流电阻率倾斜各向异性自适应有限元正演[J]. 石油物探,2017,56(5):766-773. YAN Bo,HAN Bo. Two-dimensional DC resistivity tilt anisotropic adaptive finite element method[J]. Oil Geophysical Exploration,2017,56(5):766-773.

[11] 张宏萄. 电阻率有限元正演方法研究与应用[D]. 北京:中国地质大学(北京),2011.

[12] 张继锋,智庆全,李貅,等. 起伏地形电偶源2.5维有限元数值模拟[J]. 中国有色金属学报,2013,23(9):2540-2550. ZHANG Jifeng,ZHI Qingquan,LI Xiu,et al. 2.5D finite element numerical simulation for electric dipole source on ridge terrain[J]. The Chinese Journal of Nonferrous Metals,2013,23(9):2540-2550.

[13] 鲁杏,张胜业,崔先文. 基于异常场的2.5维电阻率有限元正演模拟[J]. 地球物理学进展,2014,29(6):2718-2722. LU Xing,ZHANG Shengye,CUI Xianwen. Finite element forward modeling of 2.5 dimensional resistivity based on anomalous field[J]. Progress in Geophysics,2014,29(6):2718-2722.

[14] 金建铭,王建国,葛德彪,等. 电磁场有限元方法[M]. 西安:西安电子科技大学出版社,1998.

[15] 周熙襄. 直流电法正演计算中的有限单元法(第三讲)[J]. 物化探电子计算技术,1982(4):65-76. ZHOU Xirang. The finite element method in the calculation of forward direct current method[J]. Computing Techniques for Geophysical and Geochemical Exploration,1982(4):65-76.

[16] 张钱江,戴世坤,陈龙伟,等. 多源条件下直流电阻率法有限元三维数值模拟中一种近似边界条件[J]. 地球物理学报, 2016,59(9):3448-3458. ZHANG Qianjiang,DAI Shikun,CHEN Longwei,et al. DC resistivity method under the condition of multi-source an approximate boundary conditions in finite element three-dimensional numerical simulation[J]. Chinese Journal of Geophysics, 2016,59(9):3448-3458.