重力异常及其数据处理答题.ppt

设趋势面为 为六个待定系数。在测区中按一定网格共选取m个测点，相应点的布格异常值为gi（i=1,2,…,m）。则要使二次曲面能与重力异常的变化在最小二乘意义下得到最佳拟合，系数aj应满足 （9-84） 根据多元函数求极值法，则式（9-84）成立的条件为各阶导数等于零。见式(9-85) 根据m个点的已知重力异常值,用相应方法解出各待定系数aj,然后代入式 (9-83),便可按每个网格点的坐标(xi,yi) ,计算出该点的趋势值。 最后绘出趋势异常 (区域异常)图和剩余异常图。 趋势分析法与平滑法的区别: 1. 求区域异常的趋势分析值时,利用全区测点上的数据，而平滑时只利用计算点附近的一个区间内若干个测点的值。 2. 趋势分析时坐标原点是固定的，因而所有待定系数均应求出，才能获得各测点上的趋势值。而平滑时，原点是滑动的，即逐点计算，故只需求取系数a0即可。 利用趋势分析方法应注意: 1.效果好坏取决于所选数学模型与实际区域异常的逼近程度，因而经常用不同阶次方程作试验，从中选取最能代表区域异常相剩余异常成果的阶次。 当测区范围较大，地质情况又比较复杂时，一般不宜应用此法，因为很难用一个多项式来表示该区的区域异常，或者应分区用不同阶次的多项式来进行划分。 2.多项式阶次的选择， 原则上应视区域异常的复杂程度来决定。阶次太高，会造成趋势值受局部异常的影响较大，因而会削弱局部异常的成分。同时，也使趋势面畸变，因为阶次增加时，方程组解的误差会急剧增加，致使趋势面面目全非。 一般说选用2-3阶即可，较复杂地区也只取4~5阶。 3．利用趋势分析也同样会在分场时出现虚假异常问题。必要时可采用多次迭代的办法予以消除。 第八节 频率域滤波法 （图9-26） 为了便于讨论，把各计算公式用一般形式来表示为 一般来说，计算精度较高的公式则重力异常传递误差较大，因而受局部干扰的影响较大，反之，计算精度低的公式传递误差较小。 两者是互相矛盾的，这是因为在重力异常中有用成分与干扰成分的“频谱”并不是截然分开的，我们只能权衡利弊，在满足一定精度下，尽可能地发挥gzz 计算公式的特长。由于gzz 对于叠加重力异常的分辨率较高，因而具有较好的突出被区域场掩盖、甚至被歪曲了的浅部地质体引起的次级异常的能力。 (四)高阶导数gzz的应用 图9-27是江苏某铁矿区gzz 异常实例。我们从?g平面等值线(a)图上很难发现次级断裂F1 ，尽管对F2 有些显示，但位置也难确定;但是从gzz 的(b)图上进行解释就容易多了，从中还可以分析出这两条断层的性质并不相同，F1 主要为岩层的水平错动;而F2 则主要为两侧岩层的相对升降。 另外在应用gzz 时，如果量板的基本半径R小于地质体埋深h时，则不同半径R计算的的gzz 曲线的两侧会出现交点，在剖面上交点的水平距离与矿体宽度相当。 图9-28是江苏某矿体上的实例。此外根据（R/h）与gzz计算值近似程度之间的关系，可以大致判别不同深度地质体的分布。 这是因为埋深浅的局部异常源如矿体等，由于其h不大，R增大时，(R/h)的值明显增大，致使算得的gzz幅值显著减弱，甚至消失。而埋深较大的异常源，如基岩隆起， 由于h很大，同样加大半径R时，(R/h)的值相对变化较小 所以在gzz异常上深部地质因素减弱并不明显。图9-29为某区的布格重力异常及其高阶导数图。 其中(a)是布格异常图， (b)、(c)、(d)是用艾勒金斯公式以半径R分别为400m,200m和l00m 进行计算的结果。 从中看出R=400m时，即半径较大时，gzz异常总体为北东向异常，可 反映矿区北东向基底构造形态。 R=200m时，gzz结果表明，在原来北东向异常上叠加了一 个北东向封闭异常，可认为它反映了基岩起伏。 R=100m时，封闭异常等值线更为精细，并显示有东西向的次级异常，对比地质资料发现，该次级异常恰与已知矿体吻合。这表明计算gzz时，选用不同大小的半径，则能反映不同深度上的信息，且半径越大反映的越深。 图9-27 江苏某矿区的?g与gzz平面图 图9-28 利用不同R计算gzz 以确定矿体边界的实例 图9-28 利用不同R计算gzz以确定矿体边界的实例 图9-29 利用不同R计算gzz以分析 不同深度异常特征的实例 第七节 解析延拓法 人们把由观测平面或剖面上的已知重力异常?g值换算出高于它的平面或剖面上的异常值的过程称为向上延拓，反之则称为向下延拓。 由于重力场值是与场源到测点距离的平方成反比，因此对于深度相差较大的两个场源体来说， 进行同一