论文部分内容阅读
高密度电阻率法(High-density resistivity method)是以岩、土等介质的导电性的差异作为基础,通过人工施加稳定电流场,以达到研究地下传导电流分布规律的一种电法勘探方法。在解决水文、煤田、工程地质、以及地壳构造的研究中都得到较为广泛的应用。高密度电阻率法的正演问题是在已知地球模型和初始边界条件的情况下,求得地下电场的分布。而高密度电阻率法反演的前提之一是提出准确和完备的正演问题,并找出较为快速的算法。所以,在研究反演问题之前必须进行正演问题的研究。由于电场的复杂性,除了个别十分简单的规则形体具有解析表达式外,而其它大多数地电模型一般很难得到它们的解析解,只能通过数值计算法求得它们的近似数值解。为了解决复杂地电模型的模拟,在阅读了大量国内外的电阻率法正演数值模拟的文章基础上,采用了常用的有限差分数值模拟方法,编写了高密度电阻率法二维的正演程序,并对不同模型进行了正演分析和总结。对于点源的二维地电断面的情况,对电位沿y方向进行傅氏变换,从而引入波数域中的变换电位,进而可以得到简化了的电位的微分方程,同时采用差分方程来替代微分方程。采用基于传统的混合边界条件作为约束条件,用LLT优化措施解线性方程,得到波数域中的变换电位,再通过反傅氏变换中波数选取的最优化方法进行反傅氏变换,便可得到空间域的电位,从而根据不同的装置形式,得到相应装置的视电阻率值。采用垂直接触带模型,与解析解进行对比,算法运行良好,精度满足要求,初步验证了算法的可行性和正确性。分别计算水平低阻体与水平高阻体模型,并与国外的Pilisecky A and KnightR正演程序的计算结果进行对比,进一步验证了正演程序的正确性。对直立低阻体和直立高阻板状体模型,分别采用不同的装置进行有限差分正演模拟。对水平方向的低阻和高阻排列、双低阻排列以及双高阻排列电性异常进行正演模拟,对垂直方向的上高下低,上低下高,双高阻,双低阻不同异常体分别进行正演模拟,并对它们的结果进行对比和分析,同时对其视电阻率进行反演,来认识不同电性异常体的分辨规律。