基于格点型有限体积法的直流电测深2.5维正演模拟

摘要

采用格点型有限体积法,建立直流电测深2.5维正演模拟算法。首先,从点电源电位所满足的边值问题出发,导出格点型有限体积法控制方程的离散表达式,对边界条件近似处理后求解线性方程组即得离散电位数值解。其次,通过模拟计算均匀半空间模型和一维层状介质模型,与解析法计算结果对比,验证格点型有限体积正演算法的正确性和稳定性。分别采用本文建立的直流电测深2.5维正演模拟算法、Res2dmod开源软件构建的有限差分算法对典型二维地电模型进行正演模拟,得到2种算法的视电阻率拟断面图。研究结果表明:直流电测深2.5维正演模拟算法在一维求解域内的计算结果与解析解之间的相对误差小于2%,满足数值计算精度要求;该算法能定性判断异常体电阻率与围岩电阻率相对大小和异常体的空间展布,与经过验证的成熟软件Res2dmod相比较,电阻率数值解的绝对误差不超过3.4Ω·m,证明该算法能正确分析直流电测深2.5维异常响应规律,可为实测数据的定性解释提供指导。