论文部分内容阅读
本文采用易于拟合起伏地形和倾斜界面的非结构三角单元进行网格剖分和网格自动细化加密技术,并利用对偶加权后验误差估计算法指导网格自动细化过程,实现了2.5D直流电阻率自适应有限元数值模拟正演计算。在模型分析中,我们利用所研发的算法及软件模拟了垂直岩脉模型的直流电阻率响应,并将数值计算结果与解析解进行对比。结果表明,本文算法的数值结果精度很高,除了离电源点最近测点的数值解误差较大外,其它测点处的相对误差都在1%以内。在对本算法的有效性和计算程序的正确性进行验证后,我们模拟计算了起伏地形2.5D地电模型的视电阻率,并利用比较法对其进行了地形改正。地形改正后的视电阻率与水平地形时的结果很相似,这说明比较法可以较好地消除地形影响,突出局部地质体异常。通过比较各向同性模型直流视电阻率响应与各向异性模型视电阻率响应,我们发现电阻率各向异性对直流视电阻率响应有较大影响,在资料解释时各向异性的影响不应被忽略,考虑地下介质各向异性的影响是很有必要的。考虑了各向异性的直流电阻率数值模拟方法既能够模拟更符合实际地质情况的地电模型,也能使得异常解释更加准确。在各向异性介质正演计算的基础上,我们绘制了地层倾斜各向异性情况下,倾斜角不同时的视电阻率响应曲线和电流密度矢量图,电流密度主要沿着低电阻率地层方向流动,电流密度的大致流动方向即为低电阻率地层倾斜的方向,角度约等于倾斜角。