复杂三维介质的大地电磁正演模拟

摘要

存在多个三维且不规则异常体的复杂模型，与实际地质模型更为接近，用有限单元法模拟无疑有着明显优势。但求解过程中产生的超大型稀疏矩阵占用存储空间大，且线性方程组为严重病态的，用直接法或普通迭代法求解困难。在引入不完全LU分解改善其条件数后，用稳定双共轭梯度法求解具有迭代速度快、收敛稳定等特点，适合求解大地电磁三维有限元正演问题。通过与一维模型解析计算结果的对比和COMMEMI-3D模型正演结果分析，验证了编制程序的有效性。为后续构建大地电磁三维反演提供基础。