基于二次场算法的大地电磁二维有限单元法正演
来源：同济大学学报·自然科学版杂志2015年第8期北京时间：
作者：张志勇
单位：中南大学,同济大学 海洋与地球科学学院
摘要：推导了同时考虑电阻率与磁导率变化的大地电磁二维方程，并应用有限单元法进行数值模拟.为了提高计算精度与效率、简化计算节点生成，采用格林定理处理二次场方程源项，并设计实现了一种基于二叉树结构的收缩网格剖分算法；采用基于最少填入元思想的稀疏矩阵符号分析方法，实现了稀疏线性方程组的LDLT求解.利用二次场算法进行模型试算，结果表明所采用的新计算方法大大减少了计算单元数量，提高了计算精度与效率.
关键词：大地电磁, 有限单元法, 数值模拟, 稀疏矩阵
基金资助：国家自然科学基金项目（41304055，41304056，41264004）
