复杂二维/三维大地电磁的有限单元法正演模拟策略

摘要：

复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性.对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解.有限单元法最后形成一个线性方程组,系数矩阵足大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难.不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.通过对二维和三维模型电磁响应的讣算,验证了算法的正确性.

关键词：大地电磁、有限单元法、正演模拟、Bi-CGSTAB算法

所属期刊栏目：18

分类号：P631.3

资助基金：国家自然科学基金60672042;教育部高等学校博士学科点专项科研基金20070533102

在线出版日期：2009-04-30（万方平台首次上网日期，不代表论文的发表时间）

页数：共8页

页码：47-54