基于迭代求解法的三维有限元瞬变电磁正演计算方法

本发明涉及地球物理勘探技术领域，具体提供了一种基于迭代求解法的三维有限元瞬变电磁正演计算方法，包括：基于矢量有限元法进行四面体空间离散，将电场赋予在空间四面体的六条棱边上；构建后推欧拉有限元方程组，并基于该方程组组装由大型稀疏矩阵构成的线性系统；基于FGMRES算法对线性系统进行迭代求解，获取空间四面体单元棱边对应电场值；基于二次嵌套方法优化空间四面体网格，并基于三阶后推欧拉差分格式优化时间步长；采用网格预加密方法，降低迭代求解的计算误差，并使用区域分解方法实现并行计算，大幅度提高计算效率。本发明可满足三维瞬变电磁正演计算中的高精度和高效率需求，为复杂地质环境下的资源勘探提供了先进的计算工具。