论文部分内容阅读
随着勘探精细化和三维探测的发展,大规模和海量数据的电磁法三维正反演计算成为当前研究难点和热点问题之一。近年来,电磁法三维正反演在计算速度提升、算法优化和并行加速等方面取得重大进步;但受内存空间、计算效率以及正演精度等限制,仍难于在实际勘探中得到广泛应用。区域分解算法(DDM)将大规模问题转化成若干小问题进行独立求解,极大缩小计算规模并节省内存空间,已成为解决大规模复杂数值问题的有力手段。据此,本文以大地电磁(MT)二维正演为例,开展多种区域分解算法研究及对比分析,为大规模电磁法三维快速正反演研究提供参考和依据。本文首先简要介绍MT二维有限差分正演的基本理论,然后介绍区域分解算法的子域剖分方式及分解模式。在此基础上,提出四种不同区域分解算法(基于元素剖分的Schur补偿算法(E_Schur)、分级区域分解算法(HDD)、基于顶点剖分的Schur补偿算法(V_Schur)和Schwarz交替法(Schwarz))实现MT二维正演模拟。最后,利用典型地电模型试算验证算法的准确性和可行性,并对子域分解方法、子域剖分形状、子域剖分组合、子域数目以及子域重叠大小等对区域分解算法的计算效率和影响因素进行综合分析,综合评价各种不同区域分解算法的计算效率和所需内存空间大小。通过多个典型地电模型大地电磁二维正演模拟,及与传统全域正演算法相比表明,四种区域分解算法的CPU计算时间均随子域数目增多而增加,所需内存随子域数目增多而降低;其中E_Schur、HDD和V_Schur三种方法所需内存随子域数目增多先急剧下降,后呈缓慢上升,但总体仍比全域求解所需内存少。Schwarz交替法所需内存随子域数目增加而降低,重叠子域组合方式和重叠规模大小对计算效率有一定影响,需进行合理优化。四种区域分解算法中,HDD算法的计算效率相对较好,E_Schur补偿算法次之,Schwarz交替法最慢。综上所诉,区域分解算法能大大降低计算所需内存,求解大规模问题具有极大优势,为电磁法多维大规模问题的正反演计算提供了一种新的思路。