重力异常的数据处理.ppt

第九章 重力异常的分离 布格重力异常包含了从深部到地表所有密度不均匀体的影响，不同地质因素引起的异常无论从幅度、分布范围，变化快慢等特征看均有所不同， 第三节 重力异常的平滑 通过野外实测所获得的观测数据，以及在室内进行各项校正中总是或多或少地存在误差，从而使所得到的异常不可能如理论曲线那样光滑;更重要的是，实测异常往往是由浅到深多种地质因素产生的叠加异常。因此，在对重力异常进行解释之前，首先要对实测异常进行数据处理，其目的是: 重力异常的数据处理 1．消除因重力测量和对测量结果进行各项校正时引进的一些偶然误差或与勘探目的无关的某些近地表小型密度不均匀体的干扰; 2．从叠加的异常中划分出与勘探目标有关的异常 3．进行位场转换以满足解异常反问题的需要，例如将?g转换成vzz、 vxz、vzzz等。 常。 第三节 重力异常的平滑 对原始重力异常在解释之前作的平滑处理是为了去掉数据中某些偶然误差，及由地表密度分布不均匀体引起的杂乱无章的重力效应，获得有意义的异常。 一、剖面异常的平滑法 (二)最小二乘平滑法 若xi 以剖面上的点距为单位，即?x=1, 取点方式如图9-10所示,则式9-3式中的xi=0，?1，?2 ……?m,把它们代入式（9-3）可解出 同理可得5点、7点、9点等平滑公式。 实际工作中究竟采用几点平均最合适，这需要根据乎滑的目的而定。一般说参加平滑的点越多，得出的曲线越平缓。 图9-11就是线性平滑效果的例子。图9-11中，参加平滑的点数越多，高频信息逐渐减弱。即短周期开始消失。 2. 二次曲线平滑法 若重力异常剖面曲线在一定范围内可视为二次曲线时，则在这个范围内，平滑公式可用下面的二次曲线方程来表示;即 二次曲线平滑公式 应用导数求极值的方法，将式 (9-6)分别对a0 、a1 和a2 求偏导数，并令其等于零，得 可由上述方程组解出a0 ，若取m=2，点距?x=1,选取被平滑的点做坐标原点，求得 图9-12 为各次曲线平滑的例子 平滑处理 （一）线性平滑公式 在重力异常平面图的一定范围内，若异常形态呈简单线性变化时，可对某一点（x,y）的异常值用下面方程来拟合表示 九点平滑公式 其中g(i,j) 是流动坐标中x=i ,y=j点的原始异常值。线性平滑取点的分布如图9-13所示。 (二) 二次曲面平滑公式 在平面图上，如果重力异常的分布在一定范围内可以用二次曲面拟合时，则平滑后的 异常值g(x,y)可用下面方程来表示，即 当x=0,y=0时，a0 值便是相应点的平滑值。a0 也是利用最小二乘法来确定， 下面直接给出常用的几个二次曲面平滑公式的系数。 1．当点数一定，阶次越低结果越平滑； 2．阶次一定，点数越多结果越平滑 3，不同阶次和不同点数的结合有时可能得到相似的平滑效果; 所以实际工作中在能达到目的的前提下，尽量利用较少的点参加平滑。这样既能节省计算工作量，又可减少周围点的损失。 上面介绍的平滑法是利用有限点的异常值计算出某一点的平滑值。若想平滑一条剖面或一个平面上各点的值，可以依次在所有点上进行滑动计算而求得。 平滑本意是为了消除研究点的偶然误差，但本着数据处理的目的，平滑法是大点距平滑的结果可以用来研究区域场形态，起到压制浅部干扰的作用（接第四节） (二)地壳深部的因素 根据天然地震及地壳测深资料，地壳结构的模式大体如图9-3所示。 在大陆区从地表直至前震旦系结晶基底的顶面，是厚度从零到十几公里的沉积岩层，结晶基底以下几十公里的范围内，是花岗岩类和玄武岩类的物质层，再往下则是橄榄岩类，在不同岩类的各分界面上，上下两侧地震波传播速度有明显的差异。莫霍洛维奇(简称莫霍面)作为地壳下玄武岩类与橄榄岩类 之间的界面，它在全球范围内基本上可连续追踪；花岗岩与玄武岩类之间也是一个密度分界面，被命名为康腊德界面.但该面在大陆区不能连续追踪，在大洋区，随花岗岩类的消失而消失。 地壳厚度的变化 (即莫霍面的起伏)、壳内各层物质密度和上地幔物质密度的横向变化对地表重力分布的影响，被称为地壳深部因素的影响。上地幔密度横向不均匀的影响是十分缓慢，大范围的、平均的布格异常特征主要是对应着莫霍面起伏(即地壳厚度变化)的。 图9-4为横贯我国东西向、重力异常和莫霍面深度对照图。可见，其异常幅值大、异常范围大，异常变化单调、平缓，因而较易识别和区分。 (二)结晶基岩内部的密度变化 由于