重磁电勘探简介

重磁电勘探简介重力勘探一、重力勘探的基本概念1．重力重力的实质是牛顿万有引力和离心力的合力。

万有引力是牛顿总结前人伽里略研究行星运动规律提出来的，认为任何物体相互之间都有吸引力，吸引力的大小与两物体的质量乘积成正比，与两物体之间的距离平方成反比，其相互之间量的关系为122m m R F f R R=⋅ (6—1) 式中 m 1，m 2——分别为任意两物体的质量；R ——两物体相互间的距离；f ——引力常数，其值在CGS 制中为6．67×10—8cm 3／g ·s 2。

上式即为牛顿万有引力定律，F 力的方向对m l 来说，是由m l 指向m 2，对m 2来说则相反。

地球是有质量的，对地球表面上任一物体来说，都有地球的吸引力。

设地球的质量为M ，地面上任一物体的质量为m ，则它们之间相互的吸引力F 可根据式(6—1)来确定，其方向如图6—1(a)所示。

由于地球近似一个球体，对地面的m 物体来说，其引力的方向指向地心。

由于地球在不断地自转，地球表面上任何物体都具有一个离心力P ，其大小由下式来决定 2P mr ω= (6—2)式中 r ——m 到地轴的垂向距离；ω——地球自转的角速度。

力P 的方向如图6—1(a)所示，径向指向外。

离心力P 随纬度的不同而变化，随着r 向两极减小而减小，从赤道的最大值减小到两极为零。

为了描述重力的空间分布，通常采取直角坐标系，以数学解析式表示，如图6—1(b)所示。

设地心为坐标原点，z 轴与地球的自转轴重合，x ，y 轴在赤道面上。

设任意点A 的坐标为(x ，y ，z)，地球内部某一质量单元dm 的坐标为(,,ξηζ)，A 点到dm 的距离为r ，则dm 对A 点单位质量的引力为 2dm r dF fr r = (6—3) ()()()12222r x y z ξηζ⎡⎤=-+-+-⎣⎦式中 r r——A 到dm 方向的单位矢量，其方向是从A 到dm 。

r 对x ，y ，z 三个坐标方向的余弦分别为：,,x y z r r rξηζ---，那么dF 在x ，y ，z 三个坐标方向的引力分量为图6—1 地球引力示意图()()()232323dm x x dF x ff dm r r r dm y y dF y f f dm r r rdm z z dF z f f dm r r r ξξηηζζ--=⋅=--=⋅=--=⋅= 地球的全部质量对A 点所产生的引力分量为()()()333V V V xF x f dm r y F y f dm r z F z f dm r ξηζ-=-=-=⎰⎰⎰ 积分号下的V 表示对整个地球进行积分。

关于离心力的三个分量分别为()2P x x ω= ()2P y y ω= ()0P z =这样重力g 在x ，y ，z 三个坐标方向的分量分别为()()()23233V V V xg x f dm xr y g y f dm y r z g z f dm r ξωηωζ-=+-=+-=⎰⎰⎰ 规定不同位置均用1g 质量所受到的重力来衡量受力的大小，这个单位质量所受到的重力通常称为该点的重力场强度。

根据牛顿第二定律G ma mg ==式中 G ——物体所受的重力；g ——重力加速度。

令m ＝1，则G=g采用单位质量所受的重力来衡量重力场强度，它在数值上和重力加速度相等，因此常用重力加速度代表重力场强度，单位同加速度，为cm/s 2，在重力勘探中称之为1伽，用Gal 表示。

实测时单位太大，常用1伽的千分之一为单位，称之为毫伽，用mGal 表示。

随着重力测量仪器精度的不断提高，取毫伽的千分之一做单位，称为微伽(μGal)。

即 lcm ／s 2二1Gal1Gal ／1000＝lmGal[又称1米盖(mGal))lmGal ／1000＝l μGal实际上地球不是圆形的，而是一个偏心率近似为1／291的椭圆体，并且地球在不断地自转着，从而使万有引力和离心力随着不同位置而变化，它的变化可以按国际正常重力公式计算()220978.031810.0053024sin 0.0000059sin 2γϕϕ=+- (6—4)式中 0γ——大地水准面上的正常值；ϕ——地球的纬度。

国际正常重力公式是为了在全球范围内预测地球表面重力值而规定出的几个关系式之一，还有其他的公式，只是常数不同而已。

2．重力异常由于地球是个椭圆球，在不断地自转，从而引起地球表面上重力值的变化。

对于石油勘探来说，主要研究的是地壳密度的横向不均匀性，即由于各种地质原因使得地壳密度不均匀引起重力的变化。

如图6—2所示，地下埋藏一个密度较大的地质体，设其密度为0σ，围岩的密度为1σ，0σ >1σ，那么在其地表上，把密度为1σ的围岩在地面上产生的重力值认为是正常重力值， 在图6—2(a)中以0γ值的一条平行x 轴的直线表示。

当地下存在密度为0σ的地质体，并且其密度大于围岩密度时，球形空间里的质量就会比完全为均匀密度1σ时的质量要大，即较原先的情况会有多余的质量，通常称之为剩余质量，用M 表示，M ＝V(0σ—1σ)，V 为地质体的体积。

按照万有引力定律，这个剩余质量就会使得其相应上方地表A 1，A 2，A 3，…处的重力值比正常重力值有所增大，如图6—2(a)所示,在地质体的正上方A l 处，增加的值用F l 示之，其方向是沿垂直向下，与正常重力方向重合。

图中的A 2，A 3，…处，它们离球体越远，其重力的增加亦愈小，以AF 2，AF 3，…示之，它们的方向离地质体愈远，偏离正常重力方向的角度越大，但它们指向地质体的中心。

重力勘探所能观测的是F l ，F 2，F 3，…的垂直分量g l ，g 2，g 3，…，而不是它们的本身。

将观测到的g l ，g 2，g 3，…标在其上方的图中，作成g 曲线，这个曲线称为重力异常曲线。

其重力值的变化，称为重力异常。

当0σ >1σ时，即地下埋藏一个密度较小的地质体，如图6—2(b)所示，那么其异常曲线与图6—2(a)相反，都比正常重力0γ小，为一负重力异常曲线。

如果地层是水平地层，尽管它们之间有密度差，但不会有重力异常，如图6—2(c)所示。

可见产生重力异常的关键是水平方向岩石密度要有变化。

如图6—2(d)所示，在围岩密度为1σ的岩石中存在着密度为2σ的岩石，在密度为2σ的岩石中又存在密度为3σ的岩石，1σ >2σ >3σ。

设密度为2σ的岩石，其分布范围较密度为3σ的岩石分布的范围广，则g 的异常曲线[如图6—2(d)所示]。

反映出岩石密度的横向变化。

岩石2所引起的异常，是相对正常重力(如图中的虚线1所示)来说的；而岩石3所引起的异常，其正常重力可以认为是虚线2。

这样，要研究岩石3所引起的异常，只需用虚线2作零线即可，不必用虚线1作零线。

这与读构造图时，为了确定构造高点位置和形态，只需知道其相对高差，而无需知道它的绝对标高的情况一样。

也就是说，利用重力异常图研究局部情况，只要求知道重力的相对值就可以了。

(c)水平层状介质的重力异常(d)介质横向密度变化的重力异常图6—2 不同密度地质体的重力异常3．岩石的密度重力异常是由于地壳内部岩石密度分布的不均匀所引起的，因而对于岩石的密度及其分布情况的了解是十分必要的。

岩石密度是指在自然蕴藏条件下，岩石单位体积的质量。

根据观测结果表明，不同种类的岩石有不同的密度值；同种类岩石，在不同的地质条件下，也会有不同的密度值。

影响岩石密度的主要因素有两个，即岩石中的矿物成分和孔隙度。

一般岩浆岩和变质岩比较致密。

表6—1列出一些岩石和矿石的密度。

由表中可知，岩浆岩、变质岩所含矿物的密度比较大，大约为2.2～5.3g／cm3；大多数沉积岩，其孔隙度较岩浆岩、变质岩大，最大可达30％～50％，而且一般沉积岩所含的重矿物也较岩浆岩、变质岩为少，所以沉积岩的密度在很大程度上取决于孔隙度。

不过在沉积岩中，水、化学沉积岩的密度和成分存在着明显关系，例如石膏为2.7，岩盐为2.1，通常沉积岩的密度大约在1.1～3.0之间，比岩浆岩、变质岩小。

同是沉积岩，其密度常随埋藏深度从浅到深而增加，起初增加很快，达到一定值后，增加越来越不明显，这种密度随深度增加的关系，是因为上覆岩层的巨大压力使孔隙度减小的结果。

表6—1 常见岩石和矿石的密度岩石名称密度，g／cm3岩石名称密度，g／cm3沉积岩类火成岩类土 壤1.1～1.3 花 岗 岩2.5～3.7 砂 岩1.8～2.8 安 山 岩 2.5～2.8 页 岩2.4～3.0 辉 长 岩 2.9～3.1 石 灰 岩2.3～3.0 玄 武 岩 2.7～3.2 石 膏2.7～3.0 橄 榄 岩 2.9～3.3 岩 盐2.1～2.2 矿 石 类 变质岩类赤 铁 矿 4.9～5.3 片 麻 岩2.4～2.9 磁 铁 矿 4.9～5.2 蛇 纹 岩2.6～3.2 黄 铁 矿4.9～5.2 石 英 岩2.6～2.9 铬 铁 矿 4.5～4.6 大 理 岩 2.6～2.9 重 晶 石 4.3～4.6二、重力值的测量与校正1．重力测量的基本原理从原则上说，凡是与重力有关的物理现象，如物体的自由降落、振摆的摆动、重荷使弹簧的伸长等，都可以用来测量重力值，把它们归结起来可以分两个方面，即重力绝对值的测定和重力相对值的测定。

重力勘探所采用的是相对值的测定，其基本原理如下：如图6—3所示，它是一个由弹簧悬挂着一个重荷m 的弹簧秤，当重力有变化时，重荷将发生相应的位移，其位移的大小正比于重力大小。

当弹簧秤位于测点A 时，则根据虎克定律有如下的关系()0A A mg l l τ=-式中 m 一—重荷的质量；τ——弹簧的弹性系数；A l ——弹簧在重荷作用下的长度；0l ——弹簧不受重荷作用时的原始长度。

当弹簧秤移到B 点时，得到()0B B mg l l τ=-以上两式相减后有()()AB A B A B AB m g m g g l l lg l C l m τττ=-=-===⋅上式中C 是仪器常数，它与弹簧的性能、重荷的质量有关。

它表示重荷移动单位长度时相应的重力值的变化，称之为重力仪的格值。

测定格值的方法是借已知重力变化g 来观测重荷移动后弹簧长度的相应变化l ，从而求得格值g C l= 由此可见，已知格值就可以通过测量l 来确定任意测点间的重力g 。

图6—3 弹簧秤的基本原理 图6—4 弹簧重力仪的原理2．重力仪的原理重力仪的基本原理可以用图6—4来说明。

图6—4示出的是一根可以绕水平轴、并在垂直面上自由转动的摆杆，摆杆的一端固定着一个质量为M 的重荷，并用两个不相同的弹簧将摆杆悬挂起来，构成一个弹簧秤。

同时有两个力作用在摆杆上，即重力和弹力，重荷在重力的作用下，带动摆杆以0点为轴心向上转动，用Mgl 表示重力产生的力矩，其中l 为摆杆的长度，g 为重力值。