论文部分内容阅读
谱激电法(SIP)也称作复电阻率法(CR),它是70年代国际上发展起来的一种新的激电分支方法。近年来,复电阻率法已成为国内外电法勘探领域的研究热点之一。时间域谱激电法(TSIP)则是由频率域的复电阻率法发展而来的。与逐个频率进行频谱测量的复电阻率法相比,观测时间谱能大大提高生产效率并降低成本,因而使时间域谱激电法更具应用前景。本论文是对时间域谱激电法的基础理论研究,其中包括一维正演和2.5维有限元正演模拟及其算法的研究工作。在一维正演中本文很好的解决了汉克尔变换的数值积分问题,采用连分式算法得到了高精度的结果。2.5维有限元正演是本文的重点。与孟永良,罗延钟(2000)的算法不同,本文从频率域麦克斯韦方程出发,导出了更为简化的两个独立的电场和磁场方程。本文采用能模拟复杂介质的三角形网格剖分,基于成熟的有限单元法离散化求解方程,并采用二次场算法完成了2.5维的模拟过程。其中,为了解决二次场算法,本文特别推导了波数域一次场的直接积分公式。对于频率域波数域中的结果,依次采用反傅氏变换和逆拉氏变换变换到了时间空间域。此外,本文的一维及2.5维结果均引入了Cole-Cole模型,考虑了模型的激电效应。基于上述理论及算法,本人利用FORTRAN编写了一维正演代码,而使用MATLAB开发了2.5维有限元正演模拟程序。一维数值模拟结果与解析解的吻合说明了一维正演代码的正确性,而同一模型下2.5维有限元结果与一维正演结果的对比说明了2.5维程序的可靠性。本文一共建立了七个地电模型进行数值模拟研究。其中模型一、二、五的数值模拟结果验证了一维和2.5维正演代码的可靠性;其他模型的模拟结果则表明,在同时考虑电磁效应和激电效应的条件下,激电效应均能在时间晚期被区分开来,而且存在一定的规律性。因此,通过观测地表总场瞬变响应能够在时间晚期获取地下低阻极化体的激电信息。