第六章 磁异常的处理与转换 中国地质大学（武汉）(2).docxVIP

磁异常的转换处理是磁法勘探解释理论的一个重要组成 部分，目前磁异常的转换处理主要有圆滑、划分异常(如区 域场和局部场的分离，深源场与浅源场的分离等)、磁异常 的空间转换(由实测异常换算其他无源空间部分的磁场，也 称解析延拓);分量换算(由实测异常进行、、及之间的分 量换算)、导数换算(由实测异常计算垂向导数、水平方向 导数等)、不同磁化方向之间的换算(如化磁极等)以及曲 面上磁异常转换等等。 磁异常转换处理的方法包括空间域和频率域两类。频率 域方法由于速度快，方法简单等优点，已成为磁异常转换处 理的主要方法。 早在20世纪50年代，诸如导数异常的计算，磁场解析延 拓，化磁极等方法已相继提出。但是， 一直到20世纪60、70 年代以后，由于电子计算机的广泛应用，磁异常的转换处理 才变得容易实现，其理论和方法得到了迅速的发展，并不断 得到完善。 第二节 观测数据的圆滑、 插值 与网格化 一 、观测数据的圆滑 若 以剖面上的点距为单位，即 X, 取点的方式 如图所示： △xXi=1 △x Xi x,=0,±1,±2 -2 -1 -2 -1 0 1 把它们代则上式中 把它们代 入上式可 因此，当x=0 寸，g(0)=ao,即 由此可见，按上图取数平滑某一点的值，实际上是在剖面 上以该点为中心取奇数点的算术平均值。当m=±1 时，得 三点平滑公式为： 同理可得5点、7点、9点等平滑公式。 观 测 数 据 的 圆 滑 线性平滑效果 曲线D 形态包含了各种周期变 化的成分。 随着参加平均的点数增加， 高频成分逐渐减弱，即短 周期干扰逐渐消失。用7点平 均时， B、C 两种异常基本被 平滑掉了，只保留了原来的 低频成分A 了。在9点平 滑后，同样保留了低频成分， 只是更平滑了。 A—低频 B= 中 频 C一高频 D=A+B+C E—2 点平均 F—3 点平均 G—5 点平均 H— 7 点平均 I—9 点平均 不同点数平滑效果的对比示意图 2.二次曲线平滑公式 若异常剖面曲线在一定范围内 可视为二次曲线时，则在这个范围内，平滑公式可用下 面的二次曲线方程来表示。即 g(x,)=a?+a?x+a?x2 同样可以使用最小二乘法求出上面方程中的系数。即 2 应用导数求极值的方法，将上式分别对ao,q,a? 偏导 从而可由上述方程组解出a 若取m=2, 即采用五点圆滑时，当数据是等间距，点 距△x=1, 选取被平滑的点在坐标原点时，即取x;=0, 求得 同理可得七点二次平滑公式为 +3[g(2)+g(-2)]-2[g(3)+g(-3)]} 不同点数不同阶次平滑效果实例 平 面 异 常 的 平 滑 法 (一)线性平滑公式 在异常平面图的一定范围内，若异常形态呈简单 线性变化时，可对某一点(x,y)的异常值用下面方程来 拟合表示为 g(x,y)=ao+a?x+a?y 其中ao,a,a? 为待定系数，同样利用最小二乘法来确 定。 当x=0,y=0 时，可知g(x,y)=a? 平 面 异 常 的 平 滑 法 下面直接给出五点和九点平滑公式的结果。 五点平滑公式为 九点平滑公式为 +g(-2,0)+g(- 1,0)+g(0,- 1)+g(0,-2) 其中g(i,) 为流动坐标中x=i,y=j点的原始异常值。 ,平面异常的平滑法 y (-3,3) (-2,3) (-1,3) (3, 0) (1,3) (2,3) (3,3) (-3, 2) (-2. 3) (-1 2) (0, 2) (1, 2) (2.2) (3, 2) (-3, 1) (-2 1) (-1 1) (0, 1) (1.1) (2.1) (3.1) (-3, 0) (-2 0) (-1 0) (0, 0) (1, 0) (2.0) (3.0) X (-3, -1) (-2 -1) (-1, -1) (0, -1) (1, -1) (2. -1) (3.-1) (-3, -2) (-2 -2) (-1 -2) (0 -2) (1.-2) (2, -2) (3, -2) (-3,-3) (-2 -3) (-1, -3) (0, -3) (1.-3) (2.-3) (3, -3) 平面异常九点平滑时取点位置分布图 平 面 异 常 的 平 滑 法 圆滑前 2m+1=13 n=6 2m+1=19 n=10 2m+1=7 n=2 2m+1=19 n=2 不同点数不同阶次平滑效果对比 平 面 异 常 的 平 滑 法 (二)二次曲面平滑公式 在平面图上，如果重力异常的分布在一定范围内 g(x,y) 可用下面方程来表示，即 当x=0,y=0时， ao值便是相应点的平滑值，即 g(0,0)=a? 其中ao也是利用最小二乘法来确定，也可直接给出推 出的几个二次曲面