高密度电阻率法的正演数值模拟研究.pdfVIP

高密度电阻率法的正演数值模拟研究 摘 要 高密度电阻率法（High-density resistivity method ）是以岩、土等介质的导电性的差异作为 基础，通过人工施加稳定电流场，以达到研究地下传导电流分布规律的一种电法勘探方法。 在解决水文、煤田、工程地质、以及地壳构造的研究中都得到较为广泛的应用。 高密度电阻率法的正演问题是在已知地球模型和初始边界条件的情况下, 求得地下电场 的分布。而高密度电阻率法反演的前提之一是提出准确和完备的正演问题，并找出较为快速 的算法。所以,在研究反演问题之前必须进行正演问题的研究。 由于电场的复杂性，除了个别十分简单的规则形体具有解析表达式外, 而其它大多数地 电模型一般很难得到它们的解析解, 只能通过数值计算法求得它们的近似数值解。为了解决 复杂地电模型的模拟,在阅读了大量国内外的电阻率法正演数值模拟的文章基础上，采用了常 用的有限差分数值模拟方法，编写了高密度电阻率法二维的正演程序，并对不同模型进行了 正演分析和总结。 对于点源的二维地电断面的情况，对电位沿y方向进行傅氏变换，从而引入波数域中的 变换电位，进而可以得到简化了的电位的微分方程，同时采用差分方程来替代微分方程。采 用基于传统的混合边界条件作为约束条件，用LLT 优化措施解线性方程，得到波数域中的变 换电位，再通过反傅氏变换中波数选取的最优化方法进行反傅氏变换，便可得到空间域的电 位，从而根据不同的装置形式，得到相应装置的视电阻率值。 采用垂直接触带模型，与解析解进行对比，算法运行良好，精度满足要求，初步验证了 算法的可行性和正确性。分别计算水平低阻体与水平高阻体模型，并与国外的PiliseckyAand KnightR正演程序的计算结果进行对比，进一步验证了正演程序的正确性。对直立低阻体和 直立高阻板状体模型，分别采用不同的装置进行有限差分正演模拟。对水平方向的低阻和高 阻排列、双低阻排列以及双高阻排列电性异常进行正演模拟，对垂直方向的上高下低，上低 下高，双高阻，双低阻不同异常体分别进行正演模拟，并对它们的结果进行对比和分析，同 时对其视电阻率进行反演，来认识不同电性异常体的分辨规律。 关键词： 高密度电阻率法有限差分正演计算视电阻率值 I Forward modeling study of high-density resistivity method Abstract The high-density resistivity method (High-density resistivity method) based on the differences conductivity of rock, soil and other media a conduction current distribution as a basis, the by artificial current field is applied, in order to achieve the research underground electrical prospecting method. Have been more widely used in the study of hydrology, coalfield, engineering geology and tectonic. In the high-density resistivity method, the forward problem of high-density resistivity method is to obtain the underground high density electric field distribution by knowing the conditions of the earth model and the initial - boundary conditions. One of the prerequisites for high density inversion