Three-dimensional CSEM forward modeling using edge-based finite element method based on unstructured meshes and divergence correction
TANG Wenwu1,2, DENG Juzhi2, HUANG Qinghua1
1. School of Earth and Space Sciences, Peking University, Beijing 100871, China; 2. State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang, Jiangxi 330013, China
Abstract:The iterative computation with a system of linear equations derived from the three-dimensional (3D) controlled-source electromagnetic(CSEM) forward modeling of the electric field equation suffers from slow convergence. Moreover, unstructured meshes can make the system of linear equations more ill-posed. In view of this, we propose an algorithm for finite element forward modeling based on unstructured tetrahedral meshes and divergence correction. Starting from the divergence equation of current density, we derive the corrected divergence equation of the potential on geo-electrical interfaces. Solving the system of linear equations is accelerated with the preconditioned quasi-minimal residual (QMR) method and the alternate divergence correction during the iteration. A three-layer medium model is subjected to the forward modeling under two conditions (with/without divergence correction) to verify the reliability of the proposed algorithm. The iterative convergence and the accuracy of numerical solution of the system of linear equations indicate that the divergence correction is effective to accelerate the iteration and improve the forward modeling accuracy. On this basis, a 3D geo-electric model is built, the electromagnetic response of which is employed for the comparison of numerical solutions between the proposed algorithm and the forward modeling based on the quadratic coupling potential equation. It further confirms the high accuracy of the algorithm in this study. The modeling of a complex oil and gas monitoring system demonstrates the application potential of the CSEM method in oil and gas monitoring.
Nabighian M,Corbett J. Electromagnetic Methods in Applied Geophysics:Theory[M].SEG,1988.
[2]
He Z,Liu X,Qiu W,et al. Mapping reservoir boundary by borehole-surface TFEM:Two case studies[J].The Leading Edge,2005,24(9):896-900.
[3]
邓居智,陈辉,殷长春,等.九瑞矿集区三维电性结构研究及找矿意义[J].地球物理学报,2015,58(12):4465-4477.DENG Juzhi,CHEN Hui,YIN Changchun,et al. Three-dimensional electrical structures and significance for mineral exploration in the Jiujiang-Ruichang District[J].Chinese Journal of Geophysics,2015,58(12):4465-4477.
[4]
汤井田,何继善.可控源音频大地电磁法及其应用[M].湖南长沙:中南大学出版社,2005.
[5]
Wannamaker P E,Hohmann G W,and San Filipo W A. Electromagnetic modeling of three-dimensional bodies in layered earths using integral equations[J].Geophysics,1984,49(1):60-74.
[6]
Wannamaker P E.Advances in three-dimensional magnetotelluric modeling using integral equations[J].Geophysics,1991,56(11):1716-1728.
[7]
李静和,何展翔,孟淑君,等.三维地形频率域井筒电磁场区域积分方程法模拟[J].物理学报,2019,68(14):202-211.LI Jinghe,HE Zhanxiang,MENG Shujun,et al. Domain decomposition based integral equation modeling of 3D topography in frequency domain for well electromagnetic[J].Acta Physica Sinica,2019,68(14):202-211.
[8]
Newman G A,Alumbaugh D L. Frequency-domain modelling of airborne electromagnetic responses using staggered finite differences[J].Geophysical Prospecting,1995,43(8):1021-1042.
[9]
沈金松.用交错网格有限差分法计算三维频率域电磁响应[J].地球物理学报,2003,46(2):281-288.SHEN Jinsong. Modeling of 3D electromagnetic responses in frequency domain by using staggered-grid finite difference method[J].Chinese Journal of Geophysics,2003,46(2):281-288.
[10]
Streich R. 3D finite-difference frequency-domain mode-ling of controlled-source electromagnetic data:Direct solution and optimization for high accuracy[J].Geophysics,2009,74(5):F95-F105.
[11]
邓居智,谭捍东,陈辉,等.CSAMT三维交错采样有限差分数值模拟[J].地球物理学进展,2011,26(6):2026-2032.DENG Juzhi,TAN Handong,CHEN Hui,et al. CSAMT 3D modeling using staggered grid finite difference method[J].Progress in Geophysics,2011,26(6):2026-2032.
[12]
Haber E,Ascher U,Aruliah D,et al. Fast simulation of 3D electromagnetics problems using potentials[J].Journal of Computational Physics,2000,163:150-171.
[13]
Jahandari H and Farquharson C. Finite-volume model-ling of geophysical electromagnetic data on unstructured grids using potentials[J].Geophysical Journal International,2015,202:1859-1876.
[14]
陈辉,殷长春,邓居智.基于Lorenz规范条件下磁矢势和标势耦合方程的频率域电磁法三维正演[J].地球物理学报,2016,59(8):3087-3097.CHEN Hui,YIN Changchun,DENG Juzhi. A finite volume solution to 3D frequency-domain electromagnetic modelling using Lorenz-gauged magnetic vector and scalar potentials[J].Chinese Journal of Geophy-sics,2016,59(8):3087-3097.
[15]
Badea E A,Everett M E,Newman G A,et al. Finite-element analysis of controlled-source electromagnetic induction using Coulomb-gauged potentials[J].Geophysics,2001,66(3):786-799.
[16]
徐志锋,吴小平.可控源电磁三维频率域有限元模拟[J].地球物理学报,2010,53(8):1931-1939.XU Zhifeng,WU Xiaoping. Controlled source electromagnetic 3-D modeling in frequency domain for finite element method[J].Chinese Journal of Geophysics,2010,53(8):1931-1939.
[17]
Puzyrev V,Koldan J,Puente J,et al. A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling[J].Geophysical Journal International,2013,193(2):678-693.
[18]
Ansari S,Farquharson C G. 3D finite-element forward modeling of electromagnetic data using vector and scalar potentials and unstructured grids[J].Geophy-sics,2014,79(4):E149-E165.
[19]
Cai H,Hu X,Li J,et al. Parallelized 3D CSEM mode-ling using edge-based finite element with total field formulation and unstructured mesh[J].Computers & Geosciences,2017,99:125-134.
[20]
刘长生,汤井田,任政勇,等.基于非结构化网格的三维大地电磁自适应矢量有限元模拟[J].中南大学学报(自然科学版),2010,41(5):1855-1859.LIU Changsheng,TANG Jingtian,REN Zhengyong,et al. Three-dimension magnetotellurics modeling by adaptive edge finite-element using unstructured meshes[J].Journal of Central South University (Science and Technology),2010,41(5):1855-1859.
[21]
杨军,刘颖,吴小平.海洋可控源电磁三维非结构矢量有限元数值模拟[J].地球物理学报,2015,58(8):2827-2838.YANG Jun,LIU Ying,WU Xiaoping. 3D simulation of marine CSEM using vector finite element method on unstructured grids[J].Chinese Journal of Geophysics,2015,58(8):2827-2838.
[22]
殷长春,张博,刘云鹤,等.面向目标自适应三维大地电磁正演模拟[J].地球物理学报,2017,60(1):327-336.YIN Changchun,ZHANG Bo,LIU Yunhe,et al. A goal-oriented adaptive algorithm for 3D magnetotelluric forward modeling[J].Chinese Journal of Geophy-sics,2017,60(1):327-336.
[23]
刘颖,李予国,韩波.可控源电磁场三维自适应矢量有限元正演模拟[J].地球物理学报,2017,60(12):4874-4886.LIU Ying,LI Yuguo,HAN Bo. Adaptive edge finite element modeling of the CSEM field on unstructured grids[J].Chinese Journal of Geophysics,2017,60(12):4874-4886.
[24]
叶益信,李予国,刘颖,等.基于局部加密非结构网格的海洋可控源电磁法三维有限元正演[J].地球物理学报,2016,59(12):4747-4758.YE Yixin,LI Yuguo,LIU Ying,et al. 3D finite element modeling of marine controlled-source electromagnetic fields using locally refined unstructured meshes[J].Chinese Journal of Geophysics,2016,59(12):4747-4758.
[25]
赵宁,王绪本,余刚,等.面向目标自适应海洋可控源电磁三维矢量有限元正演[J].地球物理学报,2019,62(2):779-788.ZHAO Ning,WANG Xuben,YU Gang,et al. 3D MCSEM parallel goal-oriented adaptive vector finite element modeling[J].Chinese Journal of Geophysics,2019,62(2):779-788.
[26]
Smith J T. Conservative modeling of 3D electromagnetic fields,part Ⅱ:Biconjugate gradient solution and an accelerator[J].Geophysics,1996,61(5):1319-1324.
[27]
陈辉,邓居智,谭捍东,等.大地电磁三维交错网格有限差分数值模拟中的散度校正方法研究[J].地球物理学报,2011,54(6):1649-1659.CHEN Hui,DENG Juzhi,TAN Handong,et al. Study on divergence correction method in three dimensional magnetotelluric modeling with staggered-grid finite difference method[J].Chinese Journal of Geophysics,2011,54(6):1649-1659.
[28]
Farquharson C G,Miensopust M P. Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction[J].Journal of Applied Geophysics,2011,75:699-710.
[29]
Liu Y H,Yin C C. Electromagnetic divergence correction for 3D anisotropic EM modeling[J].Journal of Applied Geophysics,2013,77:19-27.
[30]
秦策,王绪本,赵宁,等.频率域电磁法三维有限元正演线性方程组迭代算法[J].地球物理学报,2020,63(8):3180-3190.QIN Ce,WANG Xuben,ZHAO Ning,et al. Research on the iterative solver of linear equations in three-dimensional finite element forward modeling for frequency domain electromagnetic method[J].Chinese Journal of Geophysics,2020,63(8):3180-3190.
[31]
汤文武,李耀国,柳建新,等.基于二次电场的可控源电磁法三维矢量有限元正演模拟[J].石油物探,2015,54(6):665-673.TANG Wenwu,LI Yaoguo,LIU Jianxin,et al. Three-dimensional controlled-source electromagnetic forward modeling by edge-based finite element using secondary electrical field[J].Geophysical Prospecting for Petroleum,2015,54(6):665-673.
[32]
Jin J. The Finite Element Method in Electromagne-tics[M].Wiley Interscience,2002.
[33]
Saad Y. Iterative Methods for Sparse Linear Systems[M].SIAM,2003.
[34]
Hang Si. TetGen:a delaunay-based quality tetrahedral mesh generator[J].ACM Transactions on Mathematical Software,2015,41(2):1-36.
[35]
Tang W W,Li Y G,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].Geo-physical Prospecting,2015,63(4):1491-1507.