大地电磁静态效应的二维模拟研究.pdfVIP

增刊 杜华坤等：大地电磁静态效应的二维模拟研究 135 别命名为TE模式和1M模式，可统一表示为 V·(下Vu)+Au=0 (I) 其中 V是二维哈密顿算子；V2荽+寺c，。 对于TE模式： 2∥一面占 (2) Ⅱ=Z=，孑2瓦1，)c 对于TM模式： 比=璁，芋=-七，)c=卸’ ’ (3) 盯一‘埘占 把边值问题转为下列等价变分问题 m)=厶【争…)2-争Ⅱ2】吲 +L争肭r } (4) 图1视电阻率曲线对比 of curve Ⅱ≥B_t J。 Fig．1Comparisonapparentresistivity 8F(n)=0 公式(4)即为有限元正演所依赖的变分公式， 其中AB是上边界；CD是下边界，侧边界为齐次边 界条件．不予以考虑。 2程序正确性检验 建立三层模型，其参数如下： 00012·m Pl=10012·m，P2=10Q·m，P3=l 000m hl=600m，h2=2 其中P为层电阻率；^为示层厚度。计算结果如 图1和图2所示。结果表明：视电阻率误差TE模 式0．38％，TM模式0．24％，相位误差TE模式为 0．24％，TE模式0．12％。由此可证明该程序正确 图2视相位曲线对比 且计算精度高。 of curve Fig．2Comparisonapparentphase 3静态效应的计算模拟与分析 所引起。从图6可以明显的看出由于地表不均匀 体的影响，导致不同测点的测深曲线纵向整体抬 首先，建立一个二维地电模型，如图3所示，共 升，而曲线的形态不变。在不均匀体上面的测点， X 有三层地层，在地表有一个4 40m2的不均匀体。静态效应比较明显；但随着测点远离不均匀体，静 图中不均匀体中心为坐标原点，下页图4至图lO 态效应逐渐变弱，最后直至消失。而相位曲线几乎 中d代表测点到原点的距离，p代表视电阻率，^ 不受静态效应的影响，只是高频段曲线有所变化， Solution指没有不均匀 代表异常体厚度，Analytical 这是地表不均匀体异常响应所致。总之，在二维地 体的层状解析解。 电条件下，静态效应对TE模式的影响很小，TM模 3．1 TE模式和TM模式下的静态效应 式中视电阻率受静态效应影响很大，从图7可以看 从下页图4和图5可以看出，静态效应对TE出，静态效应对相位的影响很小。 模式下的视电阻率和相位影响很小，图中曲线在高 3．2不均匀体电阻率变化时静态效应 频段有所变化，这已经不是静态位移，静态位移主 见下页图8，当不均匀体为低阻时，静态效应 要产生在低频段。因为高频段探测深度较浅，地表 使整个视电阻率曲线抬升，而且当不均匀体电阻率 不均匀体为低阻体，曲线的变化是低阻异常的响应 与围岩电阻率对比度越高，抬升的幅度越大；当不 136 物探化探计算技术 33卷 4m 600m 2000m 图3模型示意图