1、重力勘探重力异常的数据处理第四章 重力异常的数据处理布格重力异常 反映了地壳内部物质密度的不均匀性，即从地表到地下 几十公里的地壳深部，只要物质密度横向发生变化，在地下不同 的空间和范田内形成剩余质量，就可以引起地表的重力异常。定性解释 侧重于判断引起异常的地质原因，并粗略估计产生异常的地 质体的形状、产状及埋深等。定量解释 则是通过理论计算对地质体的规模、形状、产状及埋深等 作出具体解答。重力异常的推断解释的步骤： 阐明引起异常的地质因素 具体地说， 就是确定异常是浅部因素还是深部因 素引起，是矿体还是构造或其它密度不均匀体 （岩性变化、侵入体等 ）的反映。定性解释 划分和处理实测异常 重力
2、异常图往往是地表到地球深处所有密度不均匀体 产生的异常的叠加图象。为了获取探测对象产生的异常，需要将它们进行划分。不同的研究 目的提取的异常信息不同，例如，矿产调查要提取队是矿体或没部构造产生的局部异常；而 深部重力研究的目标正好相反，需要划分出的是反映地壳深部及上地幔的区域异常。 确定地质体或地质构造的赋存形态 一是根据已知地质体或地质构造的 形状、产状及埋深等研究它们引起的异常的特征，包括异常的形状、幅度、梯度及变化规 律等。二是根据异常的形态及变化规律等，确定地质体或地质构造的形状、产状、埋深及规 模等。前者足由源求场，称为止 （演 ）问题；后者是由场求源，称为反 （演 ）问题。正问题是
3、反问 题的基础，而求解反问题则是定量解择的最终目的。4.1 重力异常的主要地质原因一地壳深部因素莫霍洛维奇面：地壳与上地馒之间存在着一个界西地壳厚度各地不同，大陆平原地区大约 2030km，高山区为4060km，西藏高原达60km以上，海洋区为1020km，最薄处仅数公里。这一界面上下物质 密度差达0.3g/cm3以上，界面以上的硅镁层密度为 2.83.0g/cm3,硅侣层为2.52.7g/ cm3,界面以下物质密度为3.33.4g/ cm3。该界面的起伏引起地表重力变化的特点是导常分布植围广，幅度变化大。地 形海拔越高，地壳越厚，布格重力异常就越低，而海洋地区一般显示重力异常高。地壳厚度变化
4、与布格力异常存在近似的线性关系。内比可见，地壳深部对重 力异常的影响主要来源于莫氏界面的起伏。重力异常形态与地形起伏呈镜像关系。图4-1青藏高原布格重力异常，地壳厚度和地形剖面图的异常（图4-1）;大陆平原地区地壳厚度减小，沉积层仅数公里，重力异常强度一 般为正负几百重力单位；海洋地区显示重力正异常，强度最大可达30004000g.u.2.结晶基岩引起的重力异常在沉积岩不厚、结晶基底较浅的地台区，其基底一般由强烈变质的各种类型的结晶片岩、 片麻岩、大理岩等组成。同时还有各种类型的侵入岩，成分变化较大。此外，由于受到构造 运动的影响，结晶基岩形成褶皱和断裂，这就造成结晶基岩内部的物质密度发生横向
5、变化而 引起地表重力异常。在结晶基底深度较大地区，覆盖层有巨厚的沉积岩系，在这种情况下， 基岩内部的岩性变化和构造对重力异常影响相对减弱，基岩顶面的起伏变化则相对的对重力 异常有较大影响。结晶基岩的平均密度为 2.72.8g / cm3,上覆沉积岩平均密度为 2.52.6g/cm3，所以它们之间往往形成剩余密度为 0.10.3g/ cm3的密度分界面。在结晶基底密度较均匀的情况下，重力异常主要反映了结晶基岩表面的起伏，它的隆起和坳陷是划分地质构造 单元的重要依据。由于我国大多地区下古生界地层（以灰岩为主）和结晶基底无明显的密度变化， 而下古生界地层与上覆沉积岩之间却有明显的密度差别。这时，根据
6、重力异常研究基底起伏时，基底就 可能是包含了结晶基底和下古生界的沉积基底，人们往往把两者统称为基底。3.沉积岩层的构造和岩相变化沉积岩系不同时代和不同岩性的地层往往存在着密度差异。明显的密度界面除前面所说 的下古少代基底的顶界瓦之外，还有州 L个界面上下存在着密度捏 Sri，并且这些界面又往往与地质界而相吻合，这是利用重力勘探研究沉积岩层区域性和局部性构造的依据。这些外面 往往也是地震勘探的波阻抗界面、屯法勘探的导电性界简，例如煤系地层的底部与央陶系灰 岩的界面。当沉积岩内部有裕皱、尖灭、断裂构造使密度界面陨之起伏和紫断时，就会引起 巫力异常变化，这类异常的幅度一般为 100g.u 左右或者更
7、小，分布范围也较小，一般数平方公里到数百平方公里。沉积岩内部岩相变化、砾石层局部堆积、岩浆岩入侵体等，也能引起局部重力异常的变化。4.金属矿体和其它地质原因一般金属矿体的密度比围岩大的多，通常有 13g/cm3的密度差别，但因其范围较小，如埋藏较深时，引起的重力异常很微弱，所以勘探金属矿应使用高精度重力仪。在地表附近的浅层，如果物质密度不均匀也会引起重力异常，如地下洞穴、煤矿区的采 空区或巷道、陷落柱等都能引起重力异常。上述各种因素所发现的重力导常特征是不相同的，在研究某种地质原因引起的重力异常 时，必须根据其特征加以区分。数据处理的目的消除重力测量和对测量结果进行校正时引进的一些偶然误差或
8、与勘探目的无关的某些近地表小型密度不均匀体的干扰； 从叠加的异常中划分出与勘探目的有关的异常；进行位场转换，以满足解异常反问题的需要。 4.2 重力异常的平滑、剖面异常的平滑方法1、徒手平滑法要求： 圆滑偏差小于实测异常的均方差; 曲线围成的面积不变，中心不变,2、线性平滑公式一 1三点： g（） 3【g（ i）g（）g）一 i五点： g(0) -g( 2) g( 1) g(0) g(1) g(2)53、二次曲线圆滑处理g(x) a。ax a?x2_ 1五点：go 品17g 12(g 1 gj 3(g 2 g?) 1七点：go 肩7g 6(g 1 gj 3(g 2 g?) 2(g 3 g3)二
9、、平面异常的平滑方法1、 线性2、 二次曲面平滑公式三、多次线性内插 4.3 重力异常的识别、区域异常和局部异常区域异常：分布较广的中、深部地质因素所引起的重力异常特点：幅度和范围都比较大，梯度小局部异常：相对区域因素而言范围有限的研究对象（构造、矿产）引起的异常特点：幅度和范围都比较小，梯度大剩余异常：从布格异常中去掉区域异常后的剩余部分。局部异常和区域异常是一对相对的概念，如研究浅部地质构造问题，研究对象引起的异 常称为局部异常，把较深部的地质原因引起的异常作为区域性异常。当研究结晶基底时，就 把结晶基底引起的异常作为局部异常，地壳界面引起的异常称为区域异常。除区域异常和局部异常外，凡是比
10、勘探对象更浅层或地表的原因引起的异常统称为浅层 干扰异常。此外，干扰异常也包括偶然的测量误差、各项校正的误差等。这些干扰因素往往 使重力异常出现某些跳变，使异常曲线呈锯齿状。在对异常进行解释前，必须消除这些无意 义的跳变，使曲线光滑。叠加在区域背景场上的局部异常可使实测异常产生某些畸变，但从布格异常等值线平面 团上仍能识别局部异常的某些特征。阳1应一从巾略*帑薯直帰平面总中说的应力高 中疋m 一艸尽下尹宜幢的申Qi，制茁址* *牺童力育4：31二、重力异常的划分划分重力异常一般采用两类方法：一类是将总异常分解为几个不同的部分，把需要的保 留下来，不需要的消除掉；另一类是位场转换方法。至于从复杂
11、的干扰背景上划分有用异常, 常用的方法有数字滤波及相关分析等。1、图解法1） 直线法2） 圆滑曲线法在一定范围内（剖面、平面上）的区域异常可视为线性变化1）偏差值法图9-2! KI周法最佳半径的怙计弓 - f .区域异常变化平缓，在量板圆周范围内可以认为是呈线性变化的，因此圆周各点上区域异常的平均值应近似等于量板中心点的区域异常使。但局部异常的分布范围却十分有限，在 量板圆周上的局部异常呈现出正、负交替或幅值迅速衰减的情况，它们的平均值几乎趋于零。 因此，只要量板的圆周半径取得合适， 圆周上各点的布格异常平均值 g就不包含局部异常成分，而只含区域异常，且与量板中心处的区域异常值相近，即g g区
12、域所以，从量板中心的布格异常值中减去平均重力值，就近似得到了局部异常值。g区域 g g下是江苏某铁矿区：布格异常等值线平面图。是用六边形量板得到的局部重力异常等值线平面图，图中出现了明显的异常封闭圈。在穿过局部异常中心的 1-1面上，272号钻孔于156m深处见到厚达144m的铁矿层。阳呑江拣莫蝶旷庞利常上嵐用Kill早均摆奕稠 4.4 位场转换、位场转换的主要内容 由观测平面上 同一平面上的 Vxz Vzz由观测平面上 角 f任意高度上的在无质量分布区域，引力位 V满足拉普拉斯方程重力异常厶g也满足拉氏方程第一类边界条件:f (x,y,z)解得:f(x,y，z) 于 )2 f(，0)2 d
13、d2 (x ) (y ) z 向上延拓:g(0,0, h) ( 2 g(2 Wd d2 ( h )位函数的计算：把观测面看作为面密度为V(xy, z) 2 g( 0) dd2 (x ) (y ) z 二度体:g(o, h)二、重力异常的解析延拓1、二度体向上延拓公式：g(0, h) 0.2951 g(0,0)0.1653 g(h,0) g( h,0)0.0660 g(2h,0) g( 2h,0)0.0325 g(3h,0) g( 3h,0)0.0190 g(4h,0) g( 4h,0)0.0124 g(5h,0) g( 5h,0)0.0087.2、向下延拓公式:g(O,h) 3.7048 g(
14、0,0)1.1653 g(h,0)g( h,0)0.066q g(2h,0)g( 2h,0)0.0325 g(3h,0)g( 3h,0)0.0190 g(4h,0)g( 4h,0)0.0124 g(5h,0)g( 5h,0)0.0087由上两式可见，随着|x|增大,式中的系数迅速减少。换算时，应根据异常的大小及精度要 求，合理地取计算点数。三、重力位高次导数法2重力位二次垂向导数 Wzz（即 一g ）和三次垂向导数 Wzzz（即一/ ）可改变异常所含成分的z z比例，有利于异常的划分。丨Bi丄/一收 不同埋泵鞫悴在爲次吓向辱捡异常申的反眺上图为两个大小不一、中心埋深也不相同的球体。浅部小球的异常在实测异常中并不显 著，但对该异常求取二次垂向导数（即Wzzz）后深部大球引起的“区域异常”受到压抑而浅 部小球引起的“局部异常”得到充分显示，因此，高次导数异常就是局部异常。4* 5理 甲 Jlfl nMKS)恂利用虫力垂向导数区分拄加并#重力位高次导数还可以区分多个地质体的叠加异常。征两个平行排列的水平圆柱体上方,它们的 g异常已经叠加在一起，完全反映不出下面的地质体。重力二次垂向导数清晰地显观出了它们各自的异常。
