基于非重叠区域分解算法的大地电磁法二维正演模拟

摘要

区域分解算法采用分而治之的思想,将大规模问题转化为若干个小问题进行求解,已成为大规模数值计算领域的常用算法之一。将非重叠区域分解算法引入到大地电磁法二维正演模拟中。首先将整个求解区域分解为多个互不重叠的子域,子域之间共享边界元素;然后对每个子域采用有限差分进行离散,采用Schur补偿算法解耦得到共享边界节点上的未知数,并作为子域问题的边界条件得到关于子域内部节点的线性方程组;最后,利用直接求解算法对上述方程组进行求解,实现了大地电磁法二维正演。该算法的准确性和可行性通过多个地电模型的对比试算得到了验证。此外,还统计分析了采用不同子域分区方式和分解个数时的计算耗时,结果表明子域分区的方式对计算效率影响不大,但子域分解个数的影响则较大,进行区域分解时需要选择合适的子域个数。