摘要
提高正演计算的精度和速度,对数值模拟的有效应用具有重要的意义。基于谱元法(SEM)的大地电磁正演,从电磁场满足的2阶偏微分方程出发,采用高阶基函数对电磁场进行插值,用bicgstab求解线性方程组后得到电磁场和视电阻率的分布。为进一步提高计算速度,采用OpenMP编程模式实现了多个频点的并行计算。为了验证算法的正确性,采用国际模型作为测试算例,对计算结果进行对比分析。结果表明,谱元法相较于有限元,使用较少的网格即能达到同等的精度,有效减少了计算时间;采用并行计算实现频点之间的并行,有效减少了正演计算的时间。
-
单位桂林理工大学; 广西隐伏金属矿产勘查重点实验室