重力勘探数据处理PPT课件

1、重力仪器观测数据的初步整理的内容是通过改正（消除）重力场以外的一些因素引起的重力仪观测数据的变换， 这种变化是由仪器本身的特性和工作方法引起的。例如 零点变化及海洋、航空重力测量工作中航速，航向的影响和垂直加速度干扰等异常。一. 观测资料的整理和质量评价第1页/共66页一. 观测资料的整理和质量评价1)1)普通点观测资料的整理和质量评价 消除仪器的零点漂移, ,求取各测点相对与基点的相对重力值; ; 仪器的零点漂移: :(a)(a)纯零点漂移( (弹性元器件疲劳等造成的读数变化) )(b)(b)重力日变( (重力固体潮的变化) )和温度影响未消除部分等叠加在一起的混合零点漂移; ;对中小比例尺

2、测量, ,这一整理称为混合零点校正; ;对大比例尺测量, ,要先进行固体潮校正, ,对余下的纯零点校正和温度影响做零点校正; ;第2页/共66页一. 观测资料的整理和质量评价2) 2) 基点网观测资料的整理和质量评价在消除仪器零点变化后, ,求得各相邻两基点的重力差值( (一个边段的段差).).常用三重小循环法, ,由于测量时间差小, ,可以较好地监测重力仪的零点漂移, ,且视为线性变化. .第3页/共66页二. 重力基点网平差0ig 第4页/共66页二. 重力基点网平差12:1.520 1.2611.361 1.5400.060:1.370 1.6811.751 1.5200.080iigg

3、mGalggmGal第一环第二环第5页/共66页0iigg为平差后各相邻点的重力差；二. 重力基点网平差第6页/共66页重力闭合差的来源：a)重力仪的观测精度；b)偶然或随机误差；通常它们与仪器的零点位移关系密切；因此，闭合差的分配应与各点的闭合时间成正比，即将闭合差按时间分配到各点。二. 重力基点网平差第7页/共66页 a) 对单独个闭环计 0 iiiiiitgttKtt若若闭闭合合环环路路测测量量总总时时间间为为： t t= =平平差差系系数数为为 K K = =若若某某一一边边的的闭闭合合时时间间为为 ，则则改改边边的的平平差差值值： g g各各边边平平差差数数之之和和必必须须满满足足：

4、g gK K（一一合合路路的的平平差差算算：二. 重力基点网平差第8页/共66页11222111121212,(,)mmiFCiFCiiFCFCKtK tKtK ttFCKKK Kt(b) 对上例的两公共边FC必须公用条件：gggg 分别是两个环路的闭合差；是边的闭合时间；，是待求的系数；公共边的平差个闭合环路的平为差计算：二. 重力基点网平差第9页/共66页(1)边对应的重力差标记在圆括号内，闭合时间标记在方括号内；(2)按顺时针方向计算各环路的闭合差；(3)计算各环路的总c) 对多闭合环路的波波夫逐闭合时间；(4)从闭合差最大次逼近法平差计的环路开算：始平差；二. 重力基点网平差第10页/

5、共66页121453251 3316142365142171tt下边还是以上边讲的两个环路为例讲解方法的实现：环的总时间：分2环的总时间：分二. 重力基点网平差第11页/共66页222220.080/171420.020360.016510.024420.020BACBFCAFKKKK 2环的平差系数为：K各边的平差值分别为：B-A:C-B:F-C:A-F:将各边的平差值写在各边的对应位置上二. 重力基点网平差第12页/共66页11120.0600.0240.036FCg考虑 环的平差：在 环的公共边上，将 环的平差值变号加到该边上，这样1环的闭合差为：二. 重力基点网平差第13页/共66页1

6、111110.036/161:510.011:320.007:450.010:330.008CFDCEDFEKCFgKDCgKEDgKFEgK 环的平差系数：各边的平差值为：再将公共边平差值变化分配转入2环，按上述方法持续下去，直到所有环路完全闭合。二. 重力基点网平差第14页/共66页222220.011/171,420.003360.007510.010420.008BACBFCAFKKKK 此时，2环又有0.011的闭合差，平差系数为：K各边的平差值分别为：B-A:C-B:F-C:A-F:将各边的平差值写在各边的对应位置上将公共边上的平差值变号转入1环，1环又有了0.003的闭合差，继续

7、对1环平差。二. 重力基点网平差第15页/共66页二. 重力基点网平差第16页/共66页3） 求各基点相对于总基点的重力差 设总基点为A则各基点相对于总基点的重力差为0 1.6641.6641.664 1.7650.1011.664 1.765 1.2891.1881.664 1.765 1.289 1.3500.1621.664 1.765 1.289 1.350 1.5491.387BCDEF B点: gC点: gD点: gE点: gF点: g二. 重力基点网平差第17页/共66页三. 重力异常计算和图示 1 1. . 重力观测结果的内部校正重力观测结果的内部校正 由重力测量经过零点改正后

8、所得到的结果是测点相对于某点的重力差值，通过基点的重力值可以推算出所有的测点绝对重力值，也称为重力观测值。基点的重力值可以通过将基点网中任一基点与国家基点联测得到。各测点的重力观测值中包含着地质因素以外额外部条件变化的影响。因此，因此，由重力观测值得到反映地质因素的重力异常，需要经过一系列的校正。 下边介绍各项校正的方法和由此得到的各种异常的地球物理意义。第18页/共66页1) 内部校正概述 野外测量通常采用多台仪器,其测量值包含实际的重力值,及前边提到的各种影响因素. 其中由仪器的结构不完善造成的干扰称为仪器的内部影响,对内部影响的校正称为内部校正. 内部校正包括:温度校正,零点校正,气压校

9、正,地磁校正,目前主要作前两项校正校正后值=观测值+校正值校正值=-影响值1. 重力观测结果的内部校正第19页/共66页(1) 温度校正0020000()()TiTiiTTTTTT 线性校正：非线性校正:仪器读数时的内部温度；仪器标定的工作温度；仪器标定温度下的校正系数；1. 重力观测结果的内部校正第20页/共66页(2) 零点校正用于勘探的重力仪都有零点变化的现象，在一定精度要求下，一定时间间隔内重力仪的零点变化可以看成是与时间成线性关系，可以相应的通过改正来消除零点变化的影响。进行零点改正的方法是：重力测量工作前，在全测区均匀布置一些基点，这些基点间的重力差值采用高一级精度重力仪器或者较为

10、复杂的工作方法首先测出。在进行测点重力测量时，首先在基点上观测，然后对测点进行观测，在规定的重力仪器线性零位变化的时间间隔内再到基点上观测，当然也可以回到工作开始的那个基点观测。1. 重力观测结果的内部校正第21页/共66页2) 零点校正单次观测的零点校正方法(1)同一个基点上的两次观测的时间为(2)校对两个不同基点A和B,而且两个基点的重力差已知为由校对基点读数算出的结果是,ABABt点读数的时间为t点读数的时间为则仪器的零点变化率为：ABgABg00000, ()iittggKgK tttt 和 ，测量结果分别为 g和 g ,则单位时间内零点的变化率为：1. 重力观测结果的内部校正第22页

11、/共66页,()ABABABBAABiiAttggKtttttK tt 因此，从 至 的时间内 仪器的零漂变化率为:从 至 的任一时间内零漂的校正量为:2) 零点校正1. 重力观测结果的内部校正第23页/共66页当重力仪的零点变化只是近似线性变化,而不是线性变化,而且,在一天的工作中校对的三个基点进行了重复观测，这时零点校正应按如下方法进行：1)在方格纸上,以读数或 g为纵坐标,时间t为横坐标将重复观测点按比例点在图上;2)将相同点用直线连接,形成大致平行的直线;3)以最长的一条线的中点为标准,过中点作所有连线的平行线,形成一组线束;2) 零点校正1. 重力观测结果的内部校正第24页/共66页

12、21()0,0miiiyKxaQyKxaQQKa1. 重力观测结果的内部校正4)用这些线束的端点数据进行最小二乘拟合得到一条拟合直线,拟合直线的斜率即为该时间段内的零点变化率.第25页/共66页重力的组成和影响因素第26页/共66页第27页/共66页1) 中间层校正2 2. . 重力观测的外部校正图中B为总基点,A为A的投影,与B处于同一大地水准面上,A与B的高差为h(A的海拔高程).A与B相比,多了一个密度为 的,厚度为h的物质层,称为中间层的引力作用.由于各测点的高度不同,中间层引力的垂直分量的影响也不同,因此,必须进行中间层校正.第28页/共66页2222222222 3/2000( ,

13、 , )(cos ,sin , )000,()()VVVVzRhVzdvdvFGGf x y z dxdydzf RRz Rd d dzRzgFdvzRd dRdzRz 222在我们的问题中：=(x- )(y- )(z- )且，重力为引力的垂直分量：1) 中间层校正2 2. . 重力观测的外部校正第29页/共66页222 3/2000222222()2-1(1)221-/(2 )0.419()RhzRgGd dRdzRzgGRhRhRhhRhRRgGhRhhR 对于原盘片的引力垂直分量：当 远大于 时：1) 中间层校正2 2. . 重力观测的外部校正第30页/共66页1) 中间层校正11120

14、.4190.4190,0;0,0ghghhghg 当地形起伏不大，下图中（ ）部分的影响可以忽略不计，中间层的影响与一个水平层的影响相似：水平板的重力异常：则中间层的校正为：2 2. . 重力观测的外部校正第31页/共66页2) 地形较正2 2. . 重力观测的外部校正测点A周围的地形起伏对A点观测值的影响可用右图说明.与地面平坦相比,高于A点的 质量对A点的引力使A点的重力值减小;低于A点的地形,由于缺少物质,也使A点重力减小.所以,地形影响的校正总是正值.但在大范围测量和海洋测量中,地形校正有正有负.A第32页/共66页2) 地形较正（Terrain Corrections）2 2. .

15、重力观测的外部校正第33页/共66页a) 扇形分区地形校正计算方法2) 地形较正（Terrain Corrections）2 2. . 重力观测的外部校正0202cosPzzdmgGrdvGrr左图所示的地形，高出测点的部分以P表示，低于测点的部分以Q表示，下边分析这两部分的影响：第34页/共66页2220022220002,cos,coscosPzzzzQzrdmdvrdmdvgGGrrrQdmdvdmdvgGGrrrdvGrr 类似地，对 部分，2) 地形较正（Terrain Corrections）2 2. . 重力观测的外部校正第35页/共66页02zdvPQgGrr 可见，对于 和

16、两部分都可以表示为：实际工作中采用近似方法计算，下边介绍近似方法的实现过程.2) 地形较正（Terrain Corrections）2 2. . 重力观测的外部校正第36页/共66页1103/222,iiiiHnoHd d dgG 以测点为圆心，作许多同心圆，并从圆心作一系列等夹角的射线，将同心圆分成许多小扇形域，以该扇形域为底，地形高程差为高的区域构成梯形柱ABCD:若射线数为扇形域相对于 点的平均高度内外半径分别为 和，两侧射线与x轴的交角分别为和，则每个小扇形柱在o点产生的引力的垂直分量（重力）2) 地形较正（Terrain Corrections）2 2. . 重力观测的外部校正第37

17、页/共66页11103/22203/222013/222()iiiiiiHHiiHd d dgGod dgGdd dG 设扇形柱的密度为常数 ，则每个扇形柱在o点产生的引力的垂直分量为：若扇形柱的密度 为常数，则每个扇形柱对 点产生引力的垂直分量（重力）2) 地形较正（Terrain Corrections）2 2. . 重力观测的外部校正第38页/共66页1103/222222211222()2 /22iiiiHiiiinnGdgdnGHHngggg 若射线将测点上的同心圆作 等分，对每个小柱形域的重力校正值为：对所有扇形柱的影响的校正为:实际工作中用细剖分计算和制表来实现。2) 地形较正（

18、Terrain Corrections）2 2. . 重力观测的外部校正第39页/共66页3) 自由空间校正(高程校正,法衣校正)2222222227234,:,()11()111 2 /3/1(1/)233.087.2 10( , )AAhAAAAhAAMMgGgGRhRgggGMRhRGMGMh Rh RRh RRGMhGMhhhRRG M R 经过中间层校正和地形校正后， 点与 相比 只剩下高度 的影响因为 比 距地心更远 正常重力值随之减小 其减小值为2 2. . 重力观测的外部校正第40页/共66页3) 自由空间校正(高程校正,法衣校正)6233,(6.37 10,980000)22

19、3.08(/)3.08 ()3.08 ()CgRm gmGalRdgCgmGal mdRRRdghghh mGaldRgh mGal 为消除高程变化造成的影响需要进行自由空间校正：当高程变化 时，则：因此，自由高程校正为：2 2. . 重力观测的外部校正第41页/共66页4) 纬度校正22020978.049(1 0.0052884sin0.0000059sin 2 )978.049(1 0.0053singg在进行大区域的重力测量时，需考虑由于测点纬度不同而引起的重力的规律性变化。这种变化与地质构造无关，只是由地球三维椭圆形和自转引起的：由正常重力公式：可以近似采用）2 2. . 重力观测的

20、外部校正第42页/共66页30345.17 10sin2 ()5.17 10sin2 ()0.812sin20.812sin2)ggmGalDgmGalDRDgDmGal当地面两点的纬度差为时，重力差为：为测点与总基点的纬间距离(km)纬度校正值为：（4) 纬度校正2 2. . 重力观测的外部校正第43页/共66页5) 区域校正根据研究目标，需要将重力异常中大范围的区域地质构造引起的异常与局部构造引起的局部异常分离开，在寻找局部构造和矿体时，需要进行区域异常校正。2 2. . 重力观测的外部校正第44页/共66页2 2. . 重力观测的外部校正中间层校正和自由空间校正合起来,称为布格校正由地形

21、校正,纬度校正,中间层校正和自由空间校正得到的重力异常称为布格重力异常. 前边的四项校正中,自由空间校正的值最大,中间层校正的值次之,约为自由空间校正值的1/3.第45页/共66页6)几种重力异常的物理意义3布测量布布校纬校地校剩剩布区校法测量法纬校自校（1）布格异常 gg g g g g主要反映地壳内各种密度不均匀体引起的重力异常，是重力勘探的主要内容；（2）剩余异常 gg g g（ ）法依异常 gg g g g2 2. . 重力观测的外部校正第46页/共66页6)几种重力异常的物理意义2 2. . 重力观测的外部校正第47页/共66页6)几种重力异常的物理意义2 2. . 重力观测的外部校

22、正第48页/共66页均均测量0中地校均校自均校（4）均衡异常 gg gg g g g gg是均衡校正，后边将详细介绍，主要用于补偿海平面以上物质对测量重力的影响。 2 2. . 重力观测的外部校正第49页/共66页6). 几种重力异常的物理意义1、绝对布格重力异常： 主要反映地壳内部各种密度不均匀体所引起的重力异常2、相对布格重力异常： 主要反映地壳内部各种密度不均匀体所引起的重力异常两者的区别：两者相差一个常数。2 2. . 重力观测的外部校正第50页/共66页3、自由空间异常： 自由空间异常只是将重力观测结果归算到同一大地水 准面上，所以它的校正项只有高度改正和正常场校正偏离正常地球密度分

23、布产生的重力异常 负异常：质量亏 正异常：质量盈4、均衡异常： 仅仅反映地壳内部各种密度不均匀体所引起的重力异常均衡异常存在的地区，反映当今地壳未达均衡德地区2 2. . 重力观测的外部校正第51页/共66页3. 重力异常图的表示 布格重力异常等值图 重力异常剖面图 重力异常剖面平面图1) 等值图第52页/共66页 2) 重力异常剖面图3) 重力异常平面剖面图3. 重力异常图的表示第53页/共66页 为了使重力异常的解择推断更方便，必须用各种图件将 重力异常的分布和变化形象地表示出来，这种图仲称为 重力异常图。传统的异常图件有三种：4) 重力异常平面图3. 重力异常图的表示第54页/共66页5

24、) 用颜色和等值线表示重力异常平面特征3. 重力异常图的表示第55页/共66页 6) 用三维可视化技术表示重力异常平面特征3. 重力异常图的表示第56页/共66页1) 1735 1745,(),Bouguer均衡及均衡异常均衡的发现年间 法国地球物理学家布格 领导法国大地测量队 在秘鲁测量子午线的弧长时发现了均衡现象. 他们发现Ades(安地斯山)对铅直线有一个水平引力,该水平引力使铅直线方向发生局部变化.四. 均衡及均衡重力异常第57页/共66页 P点的铅直线方向为G1 安第斯山的引力为g1 G1和g1的合力为G2 实际观测到的重力为g3这样就发现在山下存在质量亏损,即存在均衡现象. Airy设想莫霍面在山下加深,形成一个质量亏损.1854年,英