2.5D DC resistivity forward modeling based on a local weak form element-free method
MA Changying1,2,3, YAN Lingling1,3, YAO Zhen'an1,3, LIU Jianxin2, ZHAO Wenxue1,3, ZHOU Cong1,3
1. Engineering Research Center for Seismic Disaster Prevention and Engineering Geological Disaster Detection of Jiangxi Province, East China University of Technology, Nanchang, Jiangxi 330013, China; 2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring(Central South University), Ministry of Education, Changsha, Hunan 410083, China; 3. School of Geophysics and Measurement-Control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China
Abstract:The traditional numerical simulation method relies on elements based on nodal connection information, so its application is limited. Although the element-free Galerkin method(EFGM) based on the global weak form does not require the division of elements based on nodal connection information, it still relies on the background quadrature elements divided in the global domain. In this paper, the partition of unity quadrature method is used to transform the global local integration into nodal local domain integration based on EFGM. Therefore, the background quadrature element divided in the global domain is not needed, and the dependence on the element is further reduced. As a result, 2.5-dimensional direct current(DC) resistivity forward modeling based on a local weak form elementfree method is realized with higher flexibility and accuracy. The proposed method, the EFGM, and the finite element method(FEM) are respectively used to simulate different geoelectric models, and the results are compared, which prove the correctness and effectiveness of the proposed method for the forward modeling of DC resistivity. In addition, the flexibility and adaptability of the proposed method are better than those of EFGM and FEM, with higher accuracy and low computational efficiency. In order to improve the accuracy of this method in simulating the undulating terrain, node encryption is carried out in the undulating terrain, and Gaussian points outside the terrain are not calculated. The results are basically consistent with those obtained by the EFGM and FEM.
麻昌英, 闫玲玲, 姚振岸, 柳建新, 赵文学, 周聪. 基于无单元局部弱式法的2.5维直流电阻率正演[J]. 石油地球物理勘探, 2023, 58(4): 1002-1016.
MA Changying, YAN Lingling, YAO Zhen'an, LIU Jianxin, ZHAO Wenxue, ZHOU Cong. 2.5D DC resistivity forward modeling based on a local weak form element-free method. Oil Geophysical Prospecting, 2023, 58(4): 1002-1016.
强建科, 罗延钟. 三维地形直流电阻率有限元法模拟[J]. 地球物理学报, 2007, 50(5):1606-1613.QIANG Jianke, LUO Yanzhong. The resistivity FEM numerical modeling on 3-D undulating topography[J]. Chinese Journal of Geophysics, 2007, 50(5):1606-1613.
[2]
TANG W, LI Y, SWIDINSKY A, et al. Three-dimensional controlled-source electromagnetic modelling with a well casing as a grounded source:a hybrid method of moments and finite element scheme[J]. Geophysical Prospecting, 2015, 63(6):1491-1507.
[3]
TANG W, LI Y, LIU J, et al. Three-dimensional controlled-source electromagnetic forward modeling by edge-based finite element with a divergence correction[J]. Geophysics, 2021, 86(6):E367-E382.
[4]
TANG W, CHEN W, DENG J, et al. Forward calculation of 3D controlled-source electromagnetic responses based on joint application of secondary field and coupled potential formulations[J]. Geophysics, 2022, 87(4):E253-E265.
[5]
张钱江, 戴世坤, 陈龙伟, 等. 多源条件下直流电阻率法有限元三维数值模拟中一种近似边界条件[J]. 地球物理学报, 2016, 59(9):3448-3458.ZHANG Qianjiang, DAI Shikun, CHEN Longwei, et al. An approximate boundary condition for FEM-based 3-D numerical simulation with multi-source direct current resistivity method[J]. Chinese Journal of Geophysics, 2016, 59(9):3448-3458.
[6]
XIONG B, LUO T, CHEN L. Direct solutions of 3-D magnetotelluric fields using edge-based finite element[J]. Journal of Applied Geophysics, 2018, 159:204-208.
[7]
XIE J, CAI H, HU X, et al. 3-D magnetotelluric inversion and application using the edge-based finite element with hexahedral mesh[J]. IEEE Transactions on Geoscience and Remote Sensing, 2021, 60:1-11.
[8]
戴世坤, 冉应强, 张莹, 等. 空间-波数域二维强磁场及其梯度张量数值模拟研究[J]. 石油地球物理勘探, 2023, 58(1):228-237.DAI Shikun, RAN Yingqiang, ZHANG Ying, et al. Numerical simulation of two-dimensional strong magnetic field and its gradient tensor in space-wavenumber domain[J]. Oil Geophysical Prospecting, 2023, 58(1):228-237.
[9]
李静和, 何展翔, 穆桐. 接触式激发极化法油气储层压裂监测有限元模拟[J]. 石油地球物理勘探, 2022, 57(3):719-727.LI Jinghe, HE Zhanxiang, MU Tong. Finite element simulation for fracturing monitoring of oil and gas reservoirs by using contact induced polarization method[J]. Oil Geophysical Prospecting, 2022, 57(3):719-727.
[10]
GUO R, WANG Y, EGBERT G D, et al. An efficient multigrid solver based on a four-color cell-block GaussSeidel smoother for 3D magnetotelluric forward modeling[J]. Geophysics, 2022, 87(3):E121-E133.
[11]
REN Z, TANG J. 3D direct current resistivity modeling with unstructured mesh by adaptive finite-element method[J]. Geophysics, 2010, 75(1):H7-H17.
[12]
韩骑, 胡祥云, 程正璞, 等. 自适应非结构有限元MT二维起伏地形正反演研究[J]. 地球物理学报, 2015, 58(12):4675-4684.HAN Qi, HU Xiangyun, CHENG Zhengpu, et al. A study of two-dimensional MT inversion with steep topography using the adaptive unstructured finite element method[J]. Chinese Journal of Geophysics, 2015, 58(12):4675-4684.
[13]
蔡红柱, 熊彬, ZHDANOV Michael. 电导率各向异性的海洋电磁三维有限单元法正演[J]. 地球物理学报, 2015, 58(8):2839-2850.CAI Hongzhu, XIONG Bin, ZHDANOV M. Three-dimensional marine controlled-source electromagnetic modelling in anisotropic medium using finite element method[J]. Chinese Journal of Geophysics, 2015, 58(8):2839-2850.
[14]
吴小平, 刘洋, 王威. 基于非结构网格的电阻率三维带地形反演[J]. 地球物理学报, 2015, 58(8):2706-2717.WU Xiaoping, LIU Yang, WANG Wei. 3D resistivity inversion incorporating topography based on unstructured meshes[J]. Chinese Journal of Geophysics, 2015, 58(8):2706-2717.
[15]
YAN B, LI Y, LIU Y. Adaptive finite element modeling of direct current resistivity in 2-D generally anisotropic structures[J]. Journal of Applied Geophysics, 2016, 130:169-176.
[16]
CHOU T K, CHOUTEAU M, DUBÉ J S. Intelligent meshing technique for 2D resistivity inverse problems[J]. Geophysics, 2016, 81(4):IM45-IM56.
[17]
李勇, 吴小平, 林品荣, 等. 电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟[J]. 地球物理学报, 2017, 60(5):1955-1978.LI Yong, WU Xiaoping, LIN Pinrong, et al. Three-dimensional modeling of marine controlled-source electromagnetism using the vector finite element method for arbitrary anisotropic media[J]. Chinese Journal of Geophysics, 2017, 60(5):1955-1978.
[18]
CAI H, HU X, XIONG B, et al. Finite-element time-domain modeling of electromagnetic data in general dispersive medium using adaptive Pad series[J]. Computers & Geosciences, 2017, 109(C):194-205.
[19]
曹晓月, 殷长春, 张博, 等. 面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟[J]. 地球物理学报, 2018, 61(6):2618-2628.CAO Xiaoyue, YIN Changchun, ZHANG Bo, et al. A goal-oriented adaptive finite-element method for 3D MT anisotropic modeling with topography[J]. Chinese Journal of Geophysics, 2018, 61(6):2618-2628.
[20]
CAI H, LONG Z, LIN W, et al. 3D multinary inversion of controlled-source electromagnetic data based on the finite-element method with unstructured mesh[J]. Geophysics, 2021, 86(1):E77-E92.
[21]
周峰, 张志勇, 陈煌, 等. 基于非结构网格的三种CSEM有限元三维正演系统分析[J]. 石油地球物理勘探, 2021, 56(5):1190-1202.ZHOU Feng, ZHANG Zhiyong, CHEN Huang, et al. Analysis of 3D finite-element forward modeling of CSEM data using three different formulas and unstructured grids[J]. Oil Geophysical Prospecting, 2021, 56(5):1190-1202.
[22]
汤文武, 邓居智, 黄清华. 非结构化网格与散度校正技术结合的三维CSEM矢量有限元正演[J]. 石油地球物理勘探, 2021, 56(4):891-901.TANG Wenwu, DENG Juzhi, HUANG Qinghua. Three-dimensional CSEM forward modeling using edge-based finite element method based on unstructured meshes and divergence correction[J]. Oil Geophysical Prospecting, 2021, 56(4):891-901.
[23]
赵宁, 黄明卫, 申亚行, 等. 高阶自适应有限元三维直流电阻率正演方法及其在沁水盆地煤气层压裂监测中的应用[J]. 石油地球物理勘探, 2021, 56(1):209-216.ZHAO Ning, HUANG Mingwei, SHEN Yahang, et al. Forward modeling of 3D DC resistivity based on high-order adaptive finite element and its application in Qinshui Basin[J]. Oil Geophysical Prospecting, 2021, 56(1):209-216.
[24]
TANG W, HUANG Q, DENG J, et al. Joint application of secondary field and coupled potential formulations to unstructured meshes for 3-D CSEM forward modeling[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60:1-9.
[25]
LIU G R, GU Y T. An Introduction to Meshfree Methods and Their Programming[M]. Springer, Dordrecht, 2005.
[26]
刘欣. 无网格方法[M]. 北京:科学出版社, 2011.
[27]
冯德山, 王洪华, 戴前伟. 基于无单元Galerkin法探地雷达正演模拟[J]. 地球物理学报, 2013, 56(1):298-308.FENG Deshan, WANG Honghua, DAI Qianwei. Forward simulation of Ground Penetrating Radar based on the element-free Galerkin method[J]. Chinese Journal of Geophysics, 2013, 56(1):298-308.
[28]
卢杰, 李予国. 无网格局部Petrov-Galerkin法大地电磁场二维正演模拟[J]. 地球物理学报, 2017, 60(3):1189-1200.LU Jie, LI Yuguo. Two-dimensional magnetotelluric modeling using the Mesh-free Local Petrov-Galerkin method[J]. Chinese Journal of Geophysics, 2017, 60(3):1189-1200.
[29]
BELYTSCHKO T, LU Y Y, GU L. Element-free Galerkin methods[J]. International Journal for Numerical Methods in Engineering, 1994, 37(2):229-256.
[30]
PANT M, SINGH I V, MISHRA B K. Numerical simulation of thermo-elastic fracture problems using element free Galerkin method[J]. International Journal of Mechanical Sciences, 2010, 52(12):1745-1755.
[31]
LIU F, CHENG Y. The improved element-free Galerkin method based on the nonsingular weight functions for inhomogeneous swelling of polymer gels[J]. International Journal of Applied Mechanics, 2018, 10(4):1850047.
[32]
CHEN L, LI X. A complex variable boundary elementfree method for the Helmholtz equation using regularized combined field integral equations[J]. Applied Mathematics Letters, 2020, 101:106067.
[33]
WU Q, PENG M, CHENG Y. The interpolating dimension splitting element-free Galerkin method for 3D potential problems[J]. Engineering With Computers, 2022, 38(4):2703-2717.
[34]
JIA X, HU T. Element-free precise integration method and its applications in seismic modelling and imaging[J]. Geophysical Journal International, 2006, 166(1):349-372.
[35]
王月英. 地震波正演模拟中无单元Galerkin法初探[J]. 地球物理学进展, 2007, 22(5):1539-1544.WANG Yueying. Study of element-free Galerkin method in the seismic forward modeling[J]. Progress in Geophysics, 2007, 22(5):1539-1544.
[36]
JIE Y, TANG X, LUAN M, et al. Adaptive element-free Galerkin method applied to analysis of earthquake induced liquefaction[J]. Earthquake Engineering and Engineering Vibration, 2008, 7(2):217-224.
[37]
ZHOU Z, JIA X, QIANG X. GPU-accelerated elementfree reverse-time migration with Gauss points partition[J]. Journal of Geophysics and Engineering, 2018, 15(3):718-728.
[38]
JIA X, QIANG X. A frequency-domain element-free method for seismic modeling and reverse-time migration[J]. Journal of Geophysics and Engineering, 2018, 15(4):1719-1728.
[39]
TAKEKAWA J, MIKADA H. A mesh-free finite-difference method for elastic wave propagation in the frequency-domain[J]. Computers & Geosciences, 2018, 118:65-78.
[40]
刘立彬, 段沛然, 张云银, 等. 基于无网格的地震波场数值模拟方法综述[J]. 地球物理学进展, 2020, 35(5):1815-1825.LIU Libin, DUAN Peiran, ZHANG Yunyin, et al. Overview of mesh-free method of seismic forward numerical simulation[J]. Progress in Geophysics, 2020, 35(5):1815-1825.
[41]
戴前伟, 王洪华. 基于随机介质模型的GPR无单元法正演模拟[J]. 中国有色金属学报, 2013, 23(9):2436-2443.DAI Qianwei, WANG Honghua. Element free method forward modeling of GPR based on random medium model[J]. The Chinese Journal of Nonferrous Metals, 2013, 23(9):2436-2443.
[42]
FENG D, GUO R, WANG H. An element-free Galerkin method for ground penetrating radar numerical simulation[J]. Journal of Central South University, 2015, 22(1):261-269.
[43]
苏洲, 胡文宝, 朱毅. 二维大地电磁正演中无网格算法研究[J]. 石油天然气学报, 2012, 34(5):87-90.SU Zhou, HU Wenbao, ZHU Yi. Meshfree method used in two-dimensional magnetotelluric forwarding[J]. Journal of Oil and Gas Technology, 2012, 34(5):87-90.
[44]
严家斌, 李俊杰. 无网格法在大地电磁正演计算中的应用[J]. 中南大学学报(自然科学版), 2014, 45(10):3513-3520.YAN Jiabin, LI Junjie. Magnetotelluric forward calculation by meshless method[J]. Journal of Central South University (Science and Technology), 2014, 45(10):3513-3520.
[45]
嵇艳鞠, 黄廷哲, 黄婉玉, 等. 起伏地形下各向异性的2D大地电磁无网格法数值模拟[J]. 地球物理学报, 2016, 59(12):4483-4493.JI Yanju, HUANG Tingzhe, HUANG Wanyu, et al. 2D anisotropic magnetotelluric numerical simulation using meshfree method under undulating terrain[J]. Chinese Journal of Geophysics, 2016, 59(12):4483-4493.
[46]
闫海涛, 张龙, 张继锋, 等. 全局弱式无网格法在大地电磁二维正演中的应用研究[J]. 地球物理学进展, 2019, 34(2):658-667.YAN Haitao, ZHANG Long, ZHANG Jifeng, et al. Application of global weak-form mesh-free methods in two-dimensional magnetotelluric forward[J]. Progress in Geophysics, 2019, 34(2):658-667.
[47]
LONG J, FARQUHARSON C G. On the forward modelling of three-dimensional magnetotelluric data using a radial-basis-function-based mesh-free method[J]. Geophysical Journal International, 2019, 219(1):394-416.
[48]
李俊杰, 严家斌. 重力异常二维正演中的无网格方法[J]. 煤田地质与勘探, 2018, 46(6):181-186.LI Junjie, YAN Jiabin. Meshfree method for 2-D forward modeling of gravity anomaly[J]. Coal Geology & Exploration, 2018, 46(6):181-186.
[49]
LONG J, FARQUHARSON C G. Three-dimensional forward modelling of gravity data using mesh-free methods with radial basis functions and unstructured nodes[J]. Geophysical Journal International, 2019, 217(3):1577-1601.
[50]
麻昌英, 柳建新, 刘海飞, 等. 基于全局弱式无单元法直流电阻率正演模拟[J]. 地球物理学报, 2017, 60(2):801-809.MA Changying, LIU Jianxin, LIU Haifei, et al. A global weak form element free method for direct current resistivity forward simulation[J]. Chinese Journal of Geophysics, 2017, 60(2):801-809.
[51]
麻昌英, 柳建新, 郭荣文, 等. 耦合有限单元法扩边的直流电阻率勘探无单元Galerkin法正演[J]. 地球物理学报, 2018, 61(6):2578-2588.MA Changying, LIU Jianxin, GUO Rongwen, et al. Element-free Galerkin forward modeling of DC resistivity using a coupled finite element method with extended the boundaries[J]. Chinese Journal of Geophysics, 2018, 61(6):2578-2588.
[52]
MA C, LIU J, LIU H, et al. 2.5D electric resistivity forward modeling with element-free Galerkin method[J]. Journal of Applied Geophysics, 2019, 162:47-57.
[53]
ATLURI S N, ZHU T. A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics[J]. Computational Mechanics, 1998, 22(2):117-127.
[54]
李俊杰, 严家斌. 无网格局部Petrov-Galerkin法大地电磁二维正演[J]. 煤田地质与勘探, 2014, 42(6):101-104, 109.LI Junjie, YAN iabin. Magnetotelluric two-dimensional forward modeling using meshless local Petrov-Galerkin method[J]. Coal Geology & Exploration, 2014, 42(6):101-104, 109.
[55]
MELENK J M, BABUŠKA I. The partition of unity finite element method:basic theory and applications[J]. Computer Methods in Applied Mechanics and Engineering, 1996, 139(1-4):289-314.
[56]
BABUŠKA I. MELENK J M. The partition of unity method[J]. International Journal for Numerical Methods in Engineering, 1997, 40(4):727-758.
[57]
CARPINTERI A, FERRO G, VENTURA G. The partition of unity quadrature in meshless methods[J]. International Journal for Numerical Methods in Engineering, 2002, 54(7):987-1006.
[58]
CARPINTERI A, FERRO G, VENTURA G. The partition of unity quadrature in element-free crack modelling[J]. Computers & Structures, 2003, 81(18-19):1783-1794.
[59]
RABCZUK T, ZI G. A meshfree method based on the local partition of unity for cohesive cracks[J]. Computational Mechanics, 2007, 39(6):743-760.
[60]
ZOU W, ZHOU J X, ZHANG Z Q, et al. A truly meshless method based on partition of unity quadrature for shape optimization of continua[J]. Computational Mechanics, 2007, 39(4):357-365.
[61]
RABCZUK T, BORDAS S, ZI G. On three-dimensional modelling of crack growth using partition of unity methods[J]. Computers & Structures, 2010, 88(23-24):1391-1411.
[62]
CAI Y, HAN L, TIAN L, et al. Meshless method based on Shepard function and partition of unity for two-dimensional crack problems[J]. Engineering Analysis with Boundary Elements, 2016, 65:126-135.