用重力异常分析地壳构造.ppt

全球空间重力异常图 量板法 任意形体的正反演：量板法 求任意形状地质体的重力异常的量板法，是以地质及岩（矿）石物质的资料为依据，给定地质体或地质构造的形状及剩余密度，然后用专门制作的量板求算出所引起的重力异常。量板法求重力异常的方法有很多，在本文中我们仅对二维体给出推导。 说明： ①正演：知形状及剩余密度求算重力异常 ②条件：形状为任意，剩余密度为已知 ③方法：量板法 任意形体的正反演：量板法 二维体正演量板法　 右图为二沿OY轴无限延伸的横截面积为S的二维异常密度体。选取以观测点为原点的地面坐标系（如图：称为容克夫量板又称矩形量板或异常扇形量板），由重力计算公式可知，该二维体在原点O的重力异常为 1 任意形体的正反演：量板法 容克夫量板　如图所示，作间隔为的等间隔同心圆，再以为中心作射线，使得它们满足　 则 扇形量板 因此，即每一小扇形对点的重力异常都等于a。若横截面所包含的扇形数为N，则整个二维体在O点的重力异常为： 问题：若记　　　　　　，此时 2 解析法 任意形体的正反演：解析法 利用电子计算机解任意形体重力异常正问题的解析法，实际上就是计算 所表示的三重积分的近似值。 其方法通常是将地质体划分为一系列的形状简单的单元，每个单元的重力异常可用解析法（积分或数值积分）求得。 于是整个地质体在各测点上的异常值就是这些单元分别在各测点上的异常之和（这个和也可通过数值积分求得）。 任意形体的正反演：解析法 体元法（长方体元法）　将三度体用三组平行于直角坐标面的平面进行分割，使物体分成许多直立长方体元，用解析法解出这些直立长方体元在计算点产生的异常，然后求和得到整个物体引起的异常值。 任意形体的正反演：解析法 体元法 首先看长方体下重力异常的表达式。取坐标原点位于直立长方体中心投影点，且观测点位于原点。再设长方体顶面为方形，边长为2a，顶深为h，高为H，即 把这些参数代入积分，经运算得 式中　　　　　　　　　　　　　　。此外还可求定各个高阶导数项。 34 任意形体的正反演：解析法 面元法（多边形截面法）　用一系列相互平行的铅垂面或水平面切割物体，使其分成许多直立或水平薄片。然后，用适当边数的多边形来代替每个薄片的轮廓，用解析法计算各多边形薄片在计算点产生的重力值。最后，对所有多边形薄片重力值进行数值积分，即得三度体在计算点产生的重力异常值。 任意形体的正反演：解析法 面元法 假设，计算点位于坐标原点处。当用水平面进行切割时，有 其中　　对应于铅垂薄片或水平薄片的重力效应。可用数值积分求得它的值。 可类似求得沿铅垂面进行切割时的结果。 5 6 任意形体的正反演：解析法 线元法　 789 任意形体的正反演：解析法 正反演理论的新进展 前面论述了正反演中的一些方法，它们是理论完整．效果明显．应用广泛的一些线性方法。然而．由于大多数地球物理问题部是非线性的，用线性反演方法解非线性问题，其结果必然是使解非线性方程中本来就很严重的非惟一性，更是雪上加霜，也容易陷入局部极小，闲此。一遇到非线性问题。特别是那些高维数，多极值的非线性反演问题时，非线性间题线性化的方法总感到“力不从心”． 随着地球物理学，特别是勘探地球物理学的发展，近年来，非线性反演方法发展突飞猛进。一些传统的非线性反演方法，如梯度法，牛顿法和蒙特卡洛法。 目前的一些现代方法。启发式的反演方法，如退火法，遗传算法，人工神经网络法，小波分析法等也应运而生。发展比较完善的非线性反演方法有遗传算法等。 5.重力异常与重力构造 重力异常与地球构造 研究地壳构造主要是研究地壳的内部物理特性、地壳深度、断裂情况以及地壳活动等。这些构造都会造成一定的质量或者是密度的差异。而重力异常是与地下物质质量的异常，尤其是存在物质密度差异时，重力异常就能测量并反映出来。地壳内部的属性我们无法直接研究，我们必须借助其他手段才能得到内部的信息。 中国大陆与邻区海域地壳、岩石圈厚度图（黄怀曾等，1994） 70～80地壳厚度/km；80～100岩石圈厚度 横向 纵向 岩石圈厚度分布　 重力异常与地球构造 上地幔密度分布　 中国上地幔密度分布示意图（单位：8/cm3）（冯锐，1985，转引自王相等，1996） 上地幔结构 重力是地球对物体的引力与地球自转所产生的离心力之和。小比例尺的重力测量所采集的数据，主要与地球形状和地壳深部构造有关，而大比例尺重力测量所采集的数据则主要与地壳物质的不均匀性有关。 纵向 横向 布格重力异常 用重力异常分析地壳构造 布格重力异常与地壳深部构造 　 中国10×10平均布格重力异常图 （殷秀华，史志宏，刘占坡和张玉梅编，引自丁国瑜等，1991） 利用重力资料研究地壳深部构造，主要是