基于非结构外推多重网格法的带地形二维大地电磁各向异性有限元正演

基于非结构化网格的大地电磁各向异性有限元正演模拟若采用经典迭代法,收敛速度慢,不适用于复杂地电模型的正演计算.多重网格法作为经典迭代法的改进算法,是求解椭圆方程离散线性系统的一种有效算法.然而,经典多重网格法依赖于嵌套的正交网格,无法直接应用于非结构化网格问题的求解.笔者提出一种基于半结构化网格,即通过对一个初始的非结构网格进行逐层二分加密,并利用外推瀑布式多重网格算法(EX-CMG)快速求解各向异性介质中二维大地电磁有限元正演大规模复线性系统.EXCMG算法运用三角形网格上的外推和高次插值技术,构造新的多网格延拓算子,通过两层粗网格上的数值解构造下层密网上有限元解的高阶逼近,作为多网格磨光算子BiCGStab的迭代初值,加速迭代收敛.对国际标准测试模型(COMMEMI-2D1和COMMEMI-2D4)运用该方法进行测试,得到的视电阻率和相位相对误差均在1％以内,求解时间较BiCG-Stab、聚合型代数多网格法(AGMG)大幅缩短.EXCMG算法具有较好的拓展性和适应性,可处理复杂地电模型、任意起伏地形和各向异性问题.