瞬变电磁法对含地形的三维异常体的正演计算.pdfVIP

力跚，但尚需进一步深化研究。本文采用体积分方程法H对一系列含地形模型进行瞬变电磁响应的正演 模拟。 1理论，算法及实现 层状地层中三维异常体中的总电场和总磁场满足如下积分方程： 顾力=q(，)+(吒一巳)脬(，，，’皿◇’)w (1) ㈣2缉(，)+瓴一乃彬(，’，J剐圳 (2) 这里的易、一是由外加人工源如接地或不接地回线源产生的一次场，或层状地层产生的场，而不是由平 面波源产生。吒和乃分别为异常体和第／层的电导率，张量格林函数群和G尹分别与电场和磁场有关， 且在第1层和第．，层(包括1=歹)与r’处的电流源有关【4】。 计算体积分方程法时，需要对异常体进行离散化处理。为此将异常体划分为N个小元素，每一小 元素的电场值近似取其中心点的值，则第辨个小块中心的总电场可由下式近似： 互∞珥¨+乏k一巧蜉b‘皿吨) ，一。 (3) 然后将N个小块的响应叠加起来，用于解所有小元素内的总电场： 善【p．-一口，)晔(。；，1)一丸一谬-(r1)=一昱，(，-’ (4) 为了在计算中考虑精度同时又兼顾速度，在划分异常体时采用变步长法，即在测线附近网格要密一 些，远离测线对网格可以稀疏一些；考虑到趋肤效应，在靠近地面处(埋深较浅时)，网格也要小一些，在 较深处网格则可以大一些。计算时，根据实际情况需要，长方体元素还可以进一步细分为若干立方体的 小元素，然后沿这些小元素对张量格林函数进行积分。张量格林函数在长方体元素上的积分可以看成是 由其所细分的立方体小元素积分的累加嗍。 当求出了频率域响应后，数字线性滤波法提供了一种由频率域向时间域转换的方法。Newman等以 数字滤波法近似计算正弦和余弦变换为基础，提供了一种三维响应的傅立叶变换法嘲： h陬训sin(cot)d国 (5) 掣=吾f Jl(f)：一吾j必co《彩f)露彩 (6) 刀i 珊 其中k4域嘲为磁场分量的虚部。^(o)场在t=O时的初值。以(6)式为基础的余弦变换近似算法是相对合 适的，尤其在需要精确计算晚期响应时该法更是如此。 基于以上原理，编制了用体积分方程法计算三维异常体瞬态正演响应的软件。通过与均匀半空间、 一维层状介质等理论公式结果的比较，证明该算法是正确的。 2模型的建立。 为了研究地形对瞬变电磁法勘探的影响，共设计了四个模型。地形在本文研究中主要指山脊和山谷。 分别在山脊和山谷的右边及正下方设计异常体，通过体积分离散三维异常体，并绘制出异常体电场和磁 383 场分量分布异，分析比较地形瞬变电磁法勘探效果的影响。本文中设计的异常体模型的电阻率为lQ·111， 模型一和二，分别是目标体位于山脊水平方向的右侧和正下方的情况，即是图l左边图中所示，当取小 方框1和2时的情形。 同样山谷也建立了两个模型：模型三和四，分别是目标体位于山谷水平方向的右侧和正下方的情况 (也即是图l右边图中所示，分别取小方框3和4的情况)。需要注意的是山谷作为一个异常体，其剖 面高为200米，长也为400米，电阻率为50欧姆·米。另外，所有模型中，接收器均沿着Y方向排列， 异常体埋在第一层地层中。 L=J 知／№2000Q删U 3000m500Qn一 45∞札嘲0m_ 图1理论模型，有标号的小方框分别表示四种不同的情况 3数值计算和结果分析 在运用体积分对目标体和山谷进行剖分的时候，设计目标体为三维薄板形状，山谷为三角形，实际 为带有一定梯度的台阶状。经过建立三维地电模型，绘制出了四个电场分量等值线图和四个磁场分量等 值线图，下面分别进行比较和分析。 Th。¨啪蛔．fV-．．．¨-山即-，咖md-扔t，-20‰_协蛐of300m The麟咖．fhduo协shIRtopony-0md·3D“嗍with却山of300m 图2模型一的电场等值线图以及磁场等值线图，此时