论文部分内容阅读
可控源音频大地电磁法(Controlled Source Audio-frequency Magnetotellurics,简称CSAMT法)是一种频率域电磁测深方法。其原理是采用人工源激发交变电磁场,通过观测电磁响应来计算波阻抗及视电阻率响应。频率域可控源大地电磁法沿用了AMT的测量方式,人工源采用磁性源和电性源两种方式,由于电性源CSAMT法探测深度较大,多用于煤田探测、地热、油气和固体矿产深部找矿等,本文只研究电偶极源。本文在前人的基础之上推导了水平层状介质电偶源频率域电磁场各个分量公式,正演计算过程中引入前人的两种线性滤波系数对汉克尔积分进行了求解,自己编制了FORTRAN程序实现了正演计算。在反演方法的选择上吸收采纳了前人的建议并对其进行了适当的修改,通过对多种模型的数值模拟实现了可控源音频大地电磁法的一维全资料反演。正演计算:首先假设地层为水平层状,对于均匀各向异性介质中电偶源形成的电磁场而言,将麦克斯韦方程组变换成频率域的形式;在符合物理意义的情况下通过数学方法将其转换为直角坐标系下场分量,最终得到层状介质下的正演计算公式,其中最为关键的问题是解决汉克尔变换。对于汉克尔变换,本文作者采用两套线性数字滤波系数,对比分析均匀半空间下了两套线性数字滤波系数的计算精度和计算速度,得到了一些有价值的结论,其中第二套滤波方案是首次应用在CSAMT一维正演模拟中,经过正演模拟结果显示效果较好,也算是作者的一点点创新。在推导了正演公式后,作者编写了可控源音频大地电磁一维正演程序,并对大量模型进行了数值模拟,包括H、K、A、Q等典型地电模型,同时分析了正演响应曲线的形态,做了大量的对比试验,对于可控源电磁法的处理解释有指导意义。反演成像:可控源音频大地电磁法目前多用线性反演,本文利用博斯蒂克反演结果作为初始模型,最优化反演采用广义逆线性反演方法。最优化反演中初始模型选取的好坏直接影响了反演过程,博斯蒂克反演作为一种快速成像方法在大地电磁理论中具有广泛的应用,尽管结果不够精确,但是运算简便,能直观地给出地下电阻率随深度的变化形式,利用其结果作为初始模型后可以得到较合理的初始值。后期使用广义逆反演方法对几种典型地电模型进行反演后,得到较好的结果。经过对比分析,综合运用以上两种方法,以博斯蒂克快速反演做定性解释;再以博斯蒂克反演的结果作为广义逆反演的初始模型,做精细反演,试验证明两种方法相结合不仅可以提高反演的计算效率,计算结果也更接近真实模型。