一种改进的二维大地电磁场有限元数值模拟方法

摘要

目的 提出一种改进的有限元算法,解决计算机处理地质资料时编程复杂、存储量大、运算速度慢、模型设计复杂等问题.方法 采用有限元法对二维地电模型进行剖分,在计算网格的基本结构上,通过对泛函极值求解过程的分析,得出完全类似于有限差分法中的显式迭代格式.结果 实验结果表明利用该方法计算的视电阻率值误差较小,在电场和磁场两种极化模式下有较高的分辨率.结论 将迭代法应用于泛函求解过程中,避免了生成刚度矩阵.克服了传统有限元法在二维地电模型应用中存储量大,计算速度慢的缺点.