摘要

伴随着大地电磁勘探方法的广泛应用,发展精确和快速的正反演算法成为应用地球物理学者和工程师们的研究重点。在大多数二维地电模型中,考虑到大地电磁响应不存在解析解,本文建立了格点型有限体积正演模拟算法。首先,从大地电磁场满足的变系数亥姆霍兹方程边值问题出发,导出了格点型有限体积法控制方程的离散表达式,对边界条件近似处理后求解线性方程组得到大地电磁场数值解。其次,通过对均匀半空间模型和一维介质模型的模拟计算,与解析法计算结果对比,验证了有限体积正演算法的正确性和稳定性,并分析了纵向单元剖分间距对大地电磁响应的影响。然后,利用格点型有限体积算法精确地计算了COMMEMI 2D-1模型的大地电磁响应值,该方法的计算消耗和精度都优于有限差分法。最后,通过对典型二维地电模型的大地电磁响应模拟,总结了异常响应规律,这能为实测数据的定性解释提供指导。