《文章编号1004-4140(二零一六)01-0047-08》.docVIP

文章编号：1004-4140（2009）01-0047-08 复杂三维大地电磁的有限单元法正演模拟( 张三四，李四五?，赵六七 （XX大学XX信息物理工程学院，长沙123456） 摘要：复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性。对于复杂的二维和三维大地电磁正演问题，我们采用有限单元法进行求解。有限单元法最后形成一个线性方程组，系数矩阵是大型稀疏的带状对称复系数矩阵，并且其条件数远大于1，为严重病态矩阵，求解其对应方程组会遇到很多困难。不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解，并且具有速度快、精度高和稳定性好等优点；为了模拟无穷远边界及满足计算机的内存需求，在保证计算精度的情况下设计了非均匀网格剖分；在程序编制中，只存储有限元系数矩阵的非零元素，大大减少了正演计算的时间。通过对二维和三维模型电磁响应的计算，验证了算法的正确性。 关键词：大地电磁；有限单元法；正演模拟；Bi-CGSTAB算法 中图分类号：P631.3　　　文献标识码：A 二维/三维地电模型的大地电磁响应计算必须依靠数值模拟方法，数值解通常通过微分或积分方程的逼近和线性方程的求解来得到。对于大地电磁测深的正演模拟，最常用的3种数值模拟方法有：积分方程法、有限差分法和有限元法。总的来说，3种数值模拟方法在大地电磁测深正演计算方面都有一定的优势和不足：有限差分法的优点是方法简便易