大地电磁资料静态校正方法研究.docVIP

大地电磁（MT）静态校正方法研究及应用 龙 宇1 汪 洋2 姚 强2 作者单位：1、中石化石油工程地球物理有限公司华东分公司 江苏华东八一四地球物理勘查有限公司 联系方式： 邮编： 摘要：大地电磁（MT）资料静态效应的校正是大地电磁测深法资料处理中的关键问题，如果校正不当，会使后续的反演解释得出错误的结果。本文讨论了静态效应对MT参数资料的影响，并从实用性和可操作性的角度，提出了采用阻抗张量方法和利用非电场观测的电磁法进行静态校正，并在实际生产中得到了理想的应用效果。 关键词：大地电磁（MT）、静态效应、阻抗张量、瞬变电磁法(TEM) 引言：在大地电磁测深法（MT）中，通常地表存在电性不均匀体，在电流流过不均匀体界面时，形成界面积累电荷，由此产生一个与外电场成正比的附加电场，给大地电磁测深响应带来畸变，这种影响称为静态效应。实践证明这种由地表不均匀体及复杂地形和地貌引起的静态效应严重影响到MT资料的解释和应用，所以必须通过静态校正来消除静态效应对资料的影响。 对静态校正国内外大地电磁工作者曾先后提出了一系列校正方法和技术，这些方法技术在一定条件下对静态效应起到了抑制作用。但有些方法由于条件的限制，只能在特定情况下使用。有些方法由于一些不确定因素较多，处理后的结果可信度无法把握，所以难以推广。 目前国际通用的静态校正方法有：电磁阵列剖面法(EMAP)、阻抗张量分析法、联合反演方法、利用相位资料的静校正方法和空间滤波法等方法。 在仅有MT资料的情况下，比较实用的静态校正方法是阻抗张量分解法和曲线平移法。阻抗张量分解法有个最小二乘拟合过程，所获得结果不稳定，加之假设的2D/3D模式，可能与实际情况不符，所以对这种方法的应用总是持谨慎的态度。曲线平移法采用不受静态效应影响的电磁测深数据对MT资料进行静态校正，在实际应用中效果更佳，但相应增加了施工成本，两种方法各有所长。 阻抗张量静态校正法原理及方法：在大地电磁测深法中，某一频率的电、磁场分量（Ex，Ey，Hx，Hy）满足下列方程组： Ex=ZxxHx+ZxyHy Ey=ZyxHx+ZyyHy (1-1) 写成矩阵的形式有： E=ZH (1-2) 其中： E= H= Z= (1-3) Z为阻抗张量，当测点近地表存在有局部不均匀体时，其产生的静态影响使电场发生畸变，畸变后的电场为： E?=SE （1-3）＇ 其中S为电场畸变张量，与频率无关，电场的畸变使得由其计算出的阻抗张量也发生了畸变，由式（1-2）可知畸变的阻抗张量Z?与未受静态影响的真实阻抗张量Z具有如下关系: Z?=SZ （1-4） 令C=S-1，则由上式解得Z为: Z=CZ? （1-5） C称为静态校正张量，与频率无关，显然如果能求得C，则就可获得正常的阻抗张量Z，实现静态校正。 下面来讨论C的求法，在MT勘探中，由于趋肤效应的作用，高频的勘探范围很小，所以Z在理论上应满足一维地电模型条件下的阻抗张量特征，则： Zh= （1-6） 下标h表示高频时的阻抗张量。 由式（1-5）和式（1-6）解得： C=Z0 （Zh?）-1 （1-7） 由于 （1-8） 则： （1-9） （1-10） D是一旋转不变量，其指数形式可写为，它的平方根定义为有效阻抗。式（1-6）和（1-7）中Z0随后的讨论可知，其大小不会影响视电阻率和相位的曲线形态以及电性主轴、二维偏离度等参数，可暂时令其等于有效阻抗，则（1-7）式写为： C= （1-11） 上式给出了静态校正张量的计算方法，将C乘以由实测资料计算出的Z?，则得到静态校正后的阻抗张量，尔后对Z进行常规的阻抗张量旋转等处理，求得电性主轴方位的视电阻率和相位以及其它相关参数。 对由Z计算出的视电阻率曲线再采用各种平移方法进一步进行校正，以获得合理的测点增益，完成MT的静态校正工作。 图1 某测点阻抗张量静态校正前结果 图2 某测点阻抗张量静态校正后结果 视电阻率曲线 b.相位曲线 c.电性主轴 d.二维偏离度(SKEW) 图1是某测点阻抗张量静态校正前的结果，ρxy和ρyx两条曲线严重分离，说明资料中有很强的静态干扰，特别是相位资料，φyx出现了超相位现象(