复杂二维/三维大地电磁的有限单元法正演模拟策略
童孝忠
柳建新
郭荣文
1.中南大学信息物理工程学院,长沙,4100832.中南大学信息物理工程学院,长沙,4100833.中南大学信息物理工程学院,长沙,410083
摘要：复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性.对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解.有限单元法最后形成一个线性方程组,系数矩阵足大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难.不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.通过对二维和三维模型电磁响应的讣算,验证了算法的正确性.
关键词：大地电磁有限单元法正演模拟Bi-CGSTAB算法
分类号：P631.3(天文学、地球科学)
资助基金：国家自然科学基金(60672042)高等学校博士学科点专项科研基金(20070533102)
论文发表日期：2009-01-01
在线出版日期：2025-08-15（本平台首次上网日期，不代表文献的发表时间）
页数：8( 47-54 )
英文信息