论文部分内容阅读
瞬变电磁法以其独特的优点被广大地球物理工作者所熟知,它被认为是一种非常具有发展前景的地球物理勘探方法。然而由于瞬变电磁理论的复杂性,真正熟悉和掌握该方法的人很少,严重制约了它的发展。为了推广和普及瞬变电磁法,应该让更多的人掌握电磁理论的必要知识。本文从电磁场所满足的麦克斯韦方程出发,推导出了大回线源装置下水平层状介质的时间域电磁场的表达式。该式是一个双重积分,内层积分为汉克尔型积分,外层积分为正弦积分或余弦积分。对于这两类高振荡积分,常规的数值方法可以求解,但是耗时较长。本文主要采用当前比较流行的线性数字滤波法求解这两类积分,并通过与均匀半空间情况下的解析解对比,验证了线性数字滤波法不仅计算速度快,而且精度很高,是一种较为理想的求解振荡积分的方法。本人编写了大回线源中心点一维正演程序,计算出了两层和三层介质共六种地电断面的频率域和时间域响应,细致分析了各响应曲线的形态,得到了一些有价值的结论,这对瞬变电磁资料的处理和解释具有重要的指导意义。目前,各种电磁资料拟地震处理解释研究日渐成为电磁勘探领域的热点课题,其中瞬变电磁波场变换成像技术是一种非常具有发展潜力的新方法。本文从麦克斯韦方程组出发,推导出了满足波动方程的“虚拟波场”和满足扩散方程的电磁场之间的积分关系式。波场变换式是典型的第一类Fredholm型算子方程,它的一个重要特征就是“不适定性”,离散化后所得的线性方程组是严重病态的,其病态程度往往随着方程组阶数的增加而更加严重。本文采用超松弛预条件处理有效地降低了方程组系数矩阵的条件数,然后采用正则化共轭梯度法解方程,并选用高斯脉冲作为虚拟波,采用符号积分法进行波场正变换得到时域电磁场,然后由时域电磁场和离散化后的系数矩阵反求“虚拟波场”,通过与理论波场曲线的对比验证了预条件正则化共轭梯度法求解此类严重病态方程组的有效性。文中还细致研究了子波的宽度、峰值位置和相邻子波距离等因素与波场反变换精度之间的关系,发现随着虚拟时间的延迟,反演子波具有幅值降低和波形展宽的现象。最后,进行了时间域响应加噪后的波场反变换试验。这些为瞬变电磁波场变换理论应用于实际资料的处理做了较好的铺垫。