体积分方程法模拟井间三维异常体电磁响应.pdfVIP

体积分方程法模拟井间三维异常体电磁响应+ 杨曦1潘和平1 ． 1．中国地质大学(武汉)魄球物理s空闯信怠学院．武汉430074 摘要 基于井间电磁探测方式，应用体积分方程方法对磁偶极子激发时地下介质中三维异常体 频率域电磁响应进行了数值模拟。针对并矢格林函数的计算和奇异性问题，用解析形式计算 二次场体电流分量，而用差分方式对二次场面电荷分量进行逼近，从而简化了格林函数的计 算。从本文给出的数值模拟实例可以看出，体积分方程法能有效模拟井问地下介质中电磁波 的传播，并且精度及稳定性均优于微分方程法。由于体积分方程法能够解决复杂环境下非均 匀介质体电磁散射以及大尺寸复杂目标的电磁散射问题，因此该方法适合用于井间电磁探测 的正演研究。 1 引言 对于三维井间电磁场数值方法而言，用微分方程法(有限差分法和有限元法)的缺陷在 于：必须对整个研究区域进行离散，导致了巨大的计算量，不适用于反演成像；而截断边界 条件的引入，不但增加了计算量，也带来了截断边界误差；由于无法直接通过计算得到远场 响应，必须在近场的基础上利用等效原理外推得到远场。而利用积分方程法仅需对异常体进 行剖分，因此数据量相对较小；由于引入了格林函数，使积分方程可以避免微分方程方法中 场值传递造成的误差累积问题，同时也不需要施加吸收边界条件。因此积分方程方法在模拟 三维异常体电磁响应时更为有效“11。 在求解积分方程的方法中目前应用最为广泛的是矩量法(MOM)，并在此基础上发展出了 多种改进方法，包括共轭梯度迭代结合快速傅立叶变换方法(CGFFT)、混合积分方程方法 (VSIE)、多层快速多极子方法(MLFMA)等口1。但是，三维井问电磁探测尚处于发展阶段， 目前仅仅是方法研究和试验验证，方法理论体系还需要进一步完善。因此，为了能够反演出 高精度的图像，有必要对地下介质中电磁波的传播，即对正问题作详细研究。本文利用体积 分方程方法(VIE)，通过引入并矢格林函数，采用脉冲基函数和点选配方法计算了磁偶源激 发的，均匀地下导电介质中三维异常体的井间电磁响应，实现了三维井间电磁场数值模拟， 从而高精度地重现了地下介质中低频电磁波的传播过程并研究了其规律。 4基金：中国地质大学(武汉)2006年度研究生学术探索与创新基金项目。 +作者介绍：杨曦，硕士研究生，就读于中国地质大学(武汉)地球物理与空间信息学院。 Email：yxkml69@163．com 481 2积分方程理论及其矩量法解 将井间三维电磁响应总场分为一次场和二次场，一次场是磁偶极子在均匀地下介质中产 生的场，二次场是三维异常体散射场，于是可以得到地下介质中三维异常体电磁响应总场积 分方程为心1 (1) E(，．)=Ep(，．)+IG(，．，r’)·E(，．’)·(ra(r’)咖’ V o 其中E。为偶极子产生的一次场，假如围岩介质均匀，则一次场仅与场点r有关，E(r’)为异 常体内电场，G(厂，r’)为空间并矢格林函数，吒为异常体与地下介质电导率之差。根据矩量 法(MOM)的概念【4】，将地下介质中的三维异常体剖分为N个单位立方体单元，并设每个体 元上的电场为定值，并且都等于其中心点处电场值。选取有限个基函数之和对(1)式进行逼 近，设r为单个体元上的并矢格林函数，于是有 喜l詈。屯卜栌一Ep(rm) ㈣ 其中当m=n时如为单位并矢，m≠n时am。为零并矢。此时5m。相当于脉冲基函数，将 (2)式离散，得到矩阵方程 [z]·[E]_卜E，l (3) 其中[Z]为3N×3N阶的满阵。由矩阵方程求出异．常体内电场后，即可求出空间任意一点的 电场值，应用法拉第定律，同理也可以得到空间任意一点的磁场值。 3数值模拟实例及讨论 本文应用上述方法对井间异常体三维电磁响应进行了数值模拟，不考虑磁损耗并且忽略 地磁场影响。设有无限大均匀地下围岩，其电性参数为：￡=5．6