Two-dimensional magnetotelluric anisotropic finite element modeling with undulating terrains based on unstructured extrapolation multigrid method
PAN Kejia1, HAN Xu1, WANG Pengde1, WANG Jinxuan1, XIAO Xiao2, ZHANG Lincheng3
1. School of Mathematics and Statistics, Central South University, Changsha, Hunan 410083, China; 2. School of Geosciences and Info-Physics, Central South University, Changsha, Hunan 410083, China; 3. School of Information and Electronic Engineering, Hunan City University, Yiyang, Hunan 413002, China
Abstract:Magnetotelluric anisotropic finite element forward modeling based on unstructured grids has a slow convergence rate if the classical iterative method is used, and it is not suitable for the forward calculation of complex geoelectric models. The multigrid method, as an improvement to the classical iterative method, is an effective algorithm for solving discrete linear systems of elliptic equations. However, the classical multigrid method relies on nested orthogonal grids and cannot be directly applied to solve unstructured grid problems. Therefore, the authors propose the concept of semi-structured grids. In other words, an initial unstructured grid receives binary encryption level by level, and an extrapolation cascadic multigrid method (EXCMG) is used to quickly solve the large-scale complex linear system derived from two-dimensional magnetotelluric finite element forward modeling in anisotropic media. The EXCMG uses extrapolation and high-order interpolation techniques on a triangular grid to construct a new multigrid prolongation operator and construct the high-order approximation to the finite element solution on refined grids with numerical solutions on two coarse meshes, which is then used as the iterative initial value of the multigrid smoothing operator, namely, BiCGStab, so as to accelerate its convergence. The method is tested with international standard test models (COMMEMI-2D1 and COMMEMI-2D4), and the relative differences of the apparent resistivity and phase are within 1%. In addition, the solution time of the method is significantly reduced compared with that of BiCGStab and aggregation algebraic multigrid method (AGMG). The EXCMG shows excellent scalability and adaptability and can handle complex geoelectric models, arbitrary undulating terrains, and anisotropic problems.
COGGON J H.Electromagnetic and electrical mode-ling by the finite element method[J].Geophysics,1971,36(1):132-155.
[2]
AWDEEV D,KNIZHNIK S.3D integral equation modeling with a linear dependence on dimensions[J].Geophysics,2009,74(5):F89-F94.
[3]
陈桂波,汪宏年,姚敬金,等.各向异性海底地层海洋可控源电磁响应三维积分方程法数值模拟[J].物理学报,2009,58(6):3848-3857.CHEN Guibo,WANG Hongnian,YAO Jingjin,et al.Three-dimensional numerical modeling of marine controlled-source electromagnetic responses in a layered anisotropic seabed using integral equation method[J].Acta Physica Sinica,2009,58(6):3848-3857.
[4]
李建慧,胡祥云,曾思红,等.基于电场Helmholtz方程的回线源瞬变电磁法三维正演[J].地球物理学报,2013,56(12):4256-4267.LI Jianhui,HU Xiangyun,ZENG Sihong,et al.Three-dimensional forward calculation for loop source transient electromagnetic method based on electric field Helmholtz equation[J].Chinese Journal of Geophy-sics,2013,56(12):4256-4267.
[5]
ZHOU J,HU X,CAI H.Three-dimensional finite-element analysis of magnetotelluric data using Coulomb-gauged potentials in general anisotropic media[J].Pure and Applied Geophysics,2021,178(11):4561-4581.
[6]
王宁,汤井田,任政勇,等.基于有限体积法的二维大地电磁各向异性数值模拟[J].地球物理学报,2019,62(10):3912-3922.WANG Ning,TANG Jingtian,REN Zhengyong,et al.Two-dimensional magnetotelluric anisotropic forward modeling using finite-volume method[J].Chinese Journal of Geophysics,2019,62(10):3912-3922.
[7]
WANNAMAKER P E,STODT J A,RIJO L.A stable finite element solution for two-dimensional magnetotelluric modelling[J].Geophysical Journal International,1987,88(1):277-296.
[8]
FRANKE A,BÖRNER R U,SPITZER K.Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography[J].Geophysical Journal International,2007,171(1):71-86.
徐世浙,赵生凯.二维各向异性地电断面大地电磁场的有限元法解法[J].地震学报,1985(1):80-90.XU Shizhe,ZHAO Shengkai.Solution of magnetotelluric field equations for a two-dimensional,anisotropic geoelectric section by the finite element method[J].Acta Seismologica Sinica,1985(1):80-90.
[11]
惠鑫,吴小平.电阻率各向异性介质大地电磁二维非结构有限元数值模拟[J].物探化探计算技术,2018,40(4):468-478.HUI Xin,WU Xiaoping.Two dimensional modeling of magnetotelluric in anisotropic media using unstructured finite element method[J].Computing Techniques for Geophysical and Geochemical Exploration,2018,40(4):468-478.
[12]
LI Y G.A finite-element algorithm for electromagnetic induction in two-dimensional anisotropic conductivity structures[J].Geophysical Journal International,2002,148(3):389-401.
[13]
SHEWCHUK J R.Triangle engineering:a 2D quality mesh generator and Delaunay triangulator[C].Applied Computational Geometry Towards Geometric Engineering,1996,203-222.
[14]
LI Y G and PEK J.Adaptive finite element modelling of two-dimensional magnetotelluric fields in general anisotropic media[J].Geophysical Journal International,2008,175(3):942-954.
[15]
冯凯,秦策,李论,等.基于三次插值的大地电磁自适应有限元二维正演模拟[J].物探化探计算技术,2019,41(4):456-461.FENG Kai,QIN Ce,LI Lun,et al.Magnetotelluric two-dimensional forward modeling based on adaptive finite element of cubic interpolation[J].Computing Techniques for Geophysical and Geochemical Exploration,2019,41(4):456-461.
[16]
冯凯,秦策.大地电磁(MT)自适应有限元各向异性正演[J].吉林大学学报(地球科学版),2020,50(6):1887-1896.FENG Kai,QIN Ce.Anisotropy forward modeling of magnetotelluric (MT) adaptive finite element[J].Journal of Jilin University (Earth Science Edition),2020,50(6):1887-1896.
[17]
LIU Y,XU Z H,LI Y G.Adaptive finite element mo-delling of three-dimensional magnetotelluric fields in general anisotropic media[J].Journal of Applied Geo-physics,2018,151:113-124.
[18]
LI Y G,KEY K.2D marine controlled-source electromagnetic modeling:part 1-an adaptive finite-element algorithm[J].Geophysics,2007,72(2):WA51-WA62.
[19]
REN Z Y,TANG J T.3D direct current resistivity modeling with unstructured mesh by adaptive finite-element method[J].Geophysics,2010,75(1):H7-H17.
[20]
KEY K,OVALL J.A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic mo-delling[J].Geophysical Journal International,2011,186(1):137-154.
[21]
REN Z Y,KALSCHEUER T,GREENHALGH S,et al.A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling[J].Geophysical Journal International,2013,194(2):700-718.
[22]
YIN C C,ZHANG B,LIU Y H,et al.A goal-oriented adaptive finite-element method for 3D scattered airborne electromagnetic method modeling[J].Geophysics,2016,81(5):E337-E346.
[23]
YE Y X,LI Y G,LI G,et al.3-D adaptive finite-element modeling of marine controlled-source electromagnetics with seafloor topography based on secondary potentials[J].Pure and Applied Geophysics,2018,175(12):4449-4463.
[24]
OTHMAN M,ABDULLAH A R.An efficient multigrid Poisson solver[J].International Journal of Computer Mathematics,1999,71(4):541-553.
[25]
WANG Y,ZHANG J.Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D Poisson equation[J].Journal of Computational Physics,2009,228(1):137-146.
[26]
BORNEMANN F A,DEUFLHARD P.The cascadic multigrid method for elliptic problems[J].Nume-rische Mathematik,1996,75(2):135-152.
[27]
CHEN C M,HU H L,XIE Z Q,et al.Analysis of ex-trapolation cascadic multigrid method (EXCMG)[J].Science in China Series A:Mathematics,2008,51(8):1349-1360.
[28]
胡宏伶,陈传淼,谢资清.外推瀑布多网格法(EXCMG)——大规模求解椭圆问题的新算法[J].计算数学,2009,31(3):261-274.HU Hongling,CHEN Chuanmiao,XIE Ziqing.Extrapolation cascadic multigrid method (EXCMG)-a new algorithm for solving large scale elliptic problems[J].Mathematica Numerica Sinica,2009,31(3):261-274.
[29]
潘克家,汤井田,胡宏伶,等.直流电阻率法2.5维正演的外推瀑布式多重网格法[J].地球物理学报,2012,55(8):2769-2778.PAN Kejia,TANG Jingtian,HU Hongling,et al.Extrapolation cascadic multigrid method for 2.5D direct current resistivity modeling[J].Chinese Journal of Geophysics,2012,55(8):2769-2778.
[30]
PAN K J,TANG J T.2.5-D and 3-D DC resistivity modelling using an extrapolation cascadic multigrid method[J].Geophysical Journal International,2014,197(3):1459-1470.
[31]
PAN K J,WANG J X,HU S G,et al.An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials[J].Geophysical Journal International,2022,230(3):1834-1851.
[32]
ZHDANOV M S,VARENTSOV I M,WEAVER J T,et al.Methods for modelling electromagnetic fields results from COMMEMI-the international project on the comparison of modelling methods for electromagnetic induction[J].Journal of Applied Geophysics,1997,37(3-4):133-271.
[33]
NOTAY Y.An aggregation-based algebraic multigrid method[J].Electronic Transactions on Numerical Analysis,2010,37:123-146.
[34]
徐世浙,王庆乙,王军.用边界单元法模拟二维地形对大地电磁场的影响[J].地球物理学报,1992,35(3):380-388.XU Shizhe,WANG Qingyi,WANG Jun.Modeling 2-D terrain effect on MT by the boundary element method[J].Chinese Journal of Geophysics,1992,35(3):380-388.
[35]
霍光谱,胡祥云,黄一凡,等.带地形的大地电磁各向异性二维模拟及实例对比分析[J].地球物理学报,2015,58(12):4696-4708.HUO Guangpu,HU Xiangyun,HUANG Yifan,et al.MT modeling for two-dimensional anisotropic conductivity structure with topography and examples of comparative analyses[J].Chinese Journal of Geophysics,2015,58(12):4696-4708.
[36]
胡祥云,霍光谱,高锐,等.大地电磁各向异性二维模拟及实例分析[J].地球物理学报,2013,56(12):4268-4277.HU Xiangyun,HUO Guangpu,GAO Rui,et al.The magnetotelluric anisotropic two-dimensional simulation and case analysis[J].Chinese Journal of Geophysics,2013,56(12):4268-4277.