大地电磁三维标量有限元正演及MPI并行计算

摘要

本文从麦克斯韦方程组出发,用加权余量法推导了磁场源情况下的变分方程,用六面体对研究区域进行剖分,然后在单元内进行线性插值,最后得到总体刚度矩阵.用Fortran编程实现了三维大地电磁正演程序,最后通过模型计算结果对比,验证了程序的正确性;用前人模型进行计算,并与其结果对比,二者一致.基于MPI实现了频点间的并行,对一个三维模型进行计算,并行后开启16进程时加速比达到了11,大大减少了所需时间.