基于格点型有限体积法的大地电磁二维正演

摘要：

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

关键词：大地电磁法、正演模拟、有限体积法、格点型、二维模型

所属期刊栏目：33

分类号：P631

资助基金：国家自然科学基金;湖南省自然科学基金资助项目

在线出版日期：2023-07-12（万方平台首次上网日期，不代表论文的发表时间）

页数：共13页

页码：1682-1694