广域电磁法三维积分方程法正演研究.pdfVIP

中国地球物理2011 ·235· (7)电磁方法研究与应用 广域电磁法三维积分方程法正演研究 3Dforward forwidefield methods modeling electromagnetic with method integralequation 邓锋华’李帝铨 Li Fenghua Deng Diquan 中南大学地球科学与信息物理学院长沙410083 广域电磁法继承了CSAMT使用人工场源克服场源随机性的优点，摒弃了CSAMT远区信号微弱的劣势， 扩展了观测适用的范围，使得测深能够在一定区域的非远区进行，在同等发收距上可勘探的深度增大，即勘探 同等深度可用较小的发收距。只测量一个场分量(电分量或磁分量)，提高了相对精度和工作效率。采用伪随 机编码信号，一次同时发送和接收强度相当的多个频率信号，用适合于全域的不进行简化的公式进行计算机编 程迭代反演，从而大大拓展了人工源电磁法的观测范围，提高了观测速度、精度和野外工作效率。可见，广域 电磁法具有普遍的可行性，同时对适合广域电磁法特征的电磁信号正反演理论进行研究也有着重要意义。因此， 广域电磁法三维正演研究，对广域电磁法资料解释水平和方法的应用效果具有理论和现实意义。 人工源三维正反演问题已成为地球电磁感应领域研究的前沿课题，人工源三维电磁法正演数值模拟研 究非常有意义，同时也为三维反演研究做准备。现阶段，电磁场三维数值模拟方法主要有有限差分法、有 限元法和积分方程法，与需要全空间离散的有限差分和有限元法相比，积分方程法只需对异常体内部进行 离散处理，不涉及微分方法中的吸收边界等复杂问题，当异常体仅分布在有限空间范围时，用积分方程法 求解电磁场可大大减少求解变量的数目，有利于提高计算效率及精度。由于现代计算机的迅速发展，对异 常体进行三维网格剖分和数值求积已越来越方便，同时也大大降低了计算时间，从而使三维电磁数值模拟 可以成为一种廉价、快速、能推广的解释技术。 本文基于麦克斯韦方程组，利用格林函数建立电磁场积分方程，求出已知源分布的电磁场，通过格林 函数可以将电场E或磁场Ⅳ所满足的偏微分方程加边界条件的问题，转化为积分问题进行求解。应用一种 差分近似的方法求解三重积分张量格林函数电荷源积分项一次部分，采用二次剖分的算法解决计算张量格 林函数的奇异值问题，再者，计算含有贝塞尔函数的积分项时，利用一种结合连分式展开的高斯求积代替 常规的快速汉克尔变换方法，提高了计算精度。最后根据广域电磁法定义的视电阻率公式： AT， 1 ．． 一 D。k。。二!』盟————二——一 其中， 疋，(ikx)一1—3sin。驴+e一咿(1+ikr) I u 2。1x ’‘4。 I fE。Ex(ikx) 进行计算求得广域电磁三维正演结果，验证该方法的可行性。然后通过分析典型三维地电模型的电磁响应， 给出4中不同装置的测量结果，对比广域视电阻率与Cagniard视电阻率对典型三维地质体的分辨能力，指 出广域电磁法对不同地质体的最优测量装置，完善现有的广域电磁法理论，为实际资料的正确解释提供必 要的理论依据和指导。积分方程法只需对异常区进行剖分，只需计算小体积异常区的场，计算速度快，但 只适合模拟简单模型，精确度也较低；有限差分法与有限元法虽然精确度较高，但要求对全部区域进行离 散化，占用的计算机容量较大，计算时间较长，因此本论文将对并行算法的三维电磁场正演进行研究，对 ～些复杂地形进行计算，并给出正演结果，为以后进一步对广域电磁法进行各向异