重力课件 三篇八章一节精要.pptxVIP

重力勘探 Gravitational Prospecting ;第八章：重力异常反问题;什么是正问题与反问题？;正问题：已知地质体的形状、产状和剩余密度等，通过理论计算来求得异常的分布和规律。;;重力异常反演问题 依据观测重力异常求取引起它的场源体，具体的讲是要获得地下剩余密度的空间分布，进而对其做出地质属性的判断，可分为定性解释和定量解释。;（一）叠加异常的相互影响 （二）场源物性参数的不均匀性问题 （三）测量误差与数值计算的影响 （四）反演的多解性问题 解决途径 充分搜集和利用已知的地质、钻探和物性资料，反演的多解性就会大大减小，甚至可以得到唯一的解；另外，利用多种地球物理资料的反演解释结果彼此补充，互相约束也能够有效地减少反演的多解性。 ;● 定量解释;1. 三度体剩余量的估计 位于任意封闭曲面S以内的质量M产生的引力场强，沿该曲面外法线方向的通量是 ; 在具体计算时，不可能遍及全平面，因而要分为下述两部分:一部分为实测异常在有限测区内的二重积分I‘;另一部分为积分的余项△I。 I’值可用二重增量求和的形式来计算;当观测面以等面积 划分面元时，便有;2. 三度体重心水平坐标的计算 设第三度体内某一质量单元的坐标为 ，剩余质量为m1， 则在地面P（x,y,0）点处引起的重力异常为: ; 由力学知识可知，若上述质量单元的重心在地面上的坐标为 ，应有 ;3. 二度体横截面积的求法 设在垂直二度体走向的X方向剖面上实测重力异常为△g。假定在二度体的横截面S上取一面元 ，其剩余密度为 ，dS的坐标为 。由基本公式可计算这个面元在X轴上任一点P（x,0）处引起的重力异常为 ，则下列积分: ; 如果近似积分只在有限区间（一x, x）内进行，则其误差可通过计算余项来估计。假定x值足够大，因而对积分区间以外的重力异常，可以近似地用无限长水平圆柱体的异常来代替，则积分余项可写成 ;4.二度体横截面重心水平坐标的求法; 特征点法(或任意点法)系根据异常曲线上的一些点或特征点(如极大值点、零值点、拐点)的异常值及相应的坐标求取场源体的几何或物性参数，仅适用于剩余密度为常数的几何形体。这个方法是将计算重力异常及其导数的解析式中自变量与因变量的函数关系反过来，即将异常值和相应的测点点位作为已知量，而将该异常体的产状要素和剩余密度作为未知量，并把异常体参数与异常值、点位的关系表示为显函数形式。一般情况下，这种显函数都是非线性函数，如三角函数、对数函数等，因此特征点法实质上是一种非线性反演。; 一 、 应用条件 1）基于简单规则形体的正演公式导出的，用实测曲线进行反演时，事前应对异常曲线作平滑处理，并要尽量准确地确定坐标原点的位置。 2)仅适用于单个孤立的几何形体所引起的单一异常，而实测异常往往是多种地质因引起的异常的叠加，故反演之前应进行相应的异常分离。 3) 不同几何形体其反演公式不同，使用前应明确该异常的场源体接近于何种可能的几何形体，然后选用相应的反演公式。;二、反演公式及反演方法;球体 1. △g的反演 利用主剖面上Δg正演公式;由上式可得求球体中心埋藏深度D通式为:; ; 由于实测异常并非严格对称，所以求D时最好用 的平均值，会取得更好的效果。若条件允许，还可以分别用量取的 等来求D值。当求得的数个D值较为接近，说明场源体确实近于球形，则取D的平均值会更为合理。当求得的D值彼此差别甚大或沿深度方向有规律的递增，则说明该异常不能当作球体引起的异常来解释，或者场源体沿深度还有一定的延伸。;2. 异常曲线零值点的位置指明了球心在地面投影处的位置，也即是坐标原点。通过求极值坐标的方法来解反问题，是针对这类曲线常用的方法。 ;3.;4. ;◇ 选择法;;;; 主要是由解释者凭借经验修改模型体的形状，以实现理论异常同观测异常的拟合。 具体步骤如下： 1）建立直角坐标系，横坐标沿测线方向（x），垂直向上坐标表示重力异常(△g),垂直向下坐标表示深度(z)。观测重力异常一般为一些离散分布的点值。 2) 根据重力异常的形态及已知的地质资料(例如地质体的露头位置)在x-z坐标面内，用光笔画一个初始模型，例如一个矩形，这时