采矿地球物理学概论第八章微重力课件.ppt

采矿地球物理学概论八、开采引起的重力场,8 开采引起的重力场,重力法是一种地球物理方法。它是根据地层中岩石介质质量分布的不均匀性来测量重力的异常变化。重力的变化取决于地层的尺寸、形式、埋藏深度以及该地层与周围岩层之间的密度差异。 测量仪器的高精密度（210-9ms-2）可使其广泛用于小结构，岩相变化的测量以及由于采矿引起的密度分布变化。 目前，采矿重力法主要应用于开采引起的岩体体积变化；地层震动的预测；小范围内煤层构造的变化；局部空洞的定位。,根据球体公式V=4r3/3和密度公式=M/V，通过简单的数学变换，可以将由牛顿定律所求出的地心处重力表示为： 对均匀球体而言，地表的理论重力值应该只与地理纬度有关。但实际上，不仅地球的地面起伏甚大，内部的物质密度分布也极不均匀，在结构上还存在着显著差异。这些都使得实测的重力值与理论值之间有明显的偏离，在地学上称之为重力异常。,对某地的实测重力值，通过高程及地形校正后，再减去理论重力值，差值称作重力异常值。如为正值，称正异常；如为负值，则称为负异常。前者反映该区地下的物质密度偏大，后者则说明该区地下物质密度偏小。地球物理勘探中的重力勘探方法，就是利用这一原理，通过发现各地的局部重力异常来进行找矿和勘查地下地质构造的。,重力正异常,重力负异常,重力异常在找矿中的应用,8.1 物理基础 根据牛顿定律，物体之间的引力是密度、大小、形状及距离的函数：,（8-1),式中 dF单元的引力G 引力常数G=6.6710-11 Nm2kg2dm单元质量r 单元距测量点的距离,由于岩石介质密度分布的不均匀性，可以测量重力的变化。重力（法向）与地理宽度的变化可用下式表示a(1+Sin2Sin2) (8-2)式中法向重力值a 考虑相对平面的法向重力值、系数地理宽度,8.2 井下测量重力的递减法 为了同一平面观测到的重力进行相互比较，采用Bouger递减法，如图8-1所示。假设巷道的观测点处于所取平面的最低高度。,图8-1 测量重力变化示意图,gBB= (hd+hs)Ghd+Ghg+gt+gg=0.3086(hd+hs)+0.04187(ghg-dhd)+gt+gg,（8-3）,式中：g0地球表面的平均重力值R地球的平均半径（6370m）引力常数hd测量点距约定水平的高度hs测量仪器的高度，精确到0.01mhg从地表平均高度起距测量点的深度g，d岩层hg 和hd的平均密度gt地表重力修正值gg井下井巷中观测点的重力修正值,在重力差法的测量中，地表的重力修正值gt是个常数，而实际岩石的平均密度值为2.6103kg/m2，因此，所测点的Bouger递减的微重力异常为：gg0gBB0 （8-4） 重力微异常为差的异常，计算第i次测量与第一次测量之差(gig1)或相邻二次之差(gigi-1)或第i次与所选的某次之差（gigk），式中ki。,8.3 测量方法与观测仪器 重力法进行测量时，重力仪布置在危险区域附近，测量巷道的倾角不大于10，测量巷道最好在观测区域的上部或下部，最好是石门穿过2个以上的煤层。测量点的间距一般为1050m，测点的高度精确到0.01m。,测量仪器有美国公司Worden，LaCwta-Romberg和加拿大公司Scintex生产的重力仪。重力仪主要测量重力值的变化，其灵敏度很高（10-7到10-8重力值）。其原理如图8-2所示，这种情况下，重力的微小变化，将使产生较大的倾斜度。,系统臂长，2回转轴，质量，弹簧，弹性作用臂，平衡状态时与平面的夹角，mg重力图8-2 重力仪原理示意图,8.4 多水平重力剖面 对于地层中的上下两层岩层，由于其密度不同，而且在各分层边界面不平的情况下，则在该边界的上下形成重力异常，如图8-3所示。,1复合层上部的重力异常曲线；2复合层下部的重力异常曲线图8-3 密度分布不均边界21的重力异常g,实际中，对于两水平来说，可用重力异常来： 1）揭露测量水平之间的地质情况； 2）研究开采引起的岩层密度的变化； 3）揭露岩体引力升高的影响下，膨胀过程的发展。,8.5 冲击地压的预测预报 重力法广泛用于矿山震动与冲击矿压的预测预报之中。其主要是确定在有冲击地压危险的区域，微重力异常曲线RBGA趋向性的变化特征。 一般情况下，在发生震动与冲击矿压前，岩体的体积将会增加，从而使岩体的密度降低，微重力异常值RBGA将发生变化。其变化趋势为T，该趋势的斜率估计量为A，见图8-4所示。,AMPi1定义为i1的振幅度。它是RBGA曲线的梯度值，其变化对应着岩石密度梯度的变化，与岩层开采引起的应力变化具有简单的关系。AMPi1增加时对应着压力的集中，可以提前预计观测区域内岩层释放的能量。AMPi1减小时，对应岩层应力释放，卸压。,图8-4 微重力异常的趋势及变化幅度,（8-5）,波兰与德国的研究结果表明，AMP(t)和(t)存在着如下关系： AMPs(E)D(E)Rez(E) (8-6) 式中 AMPs(E)研究期间内能量总和对应的AMPi1增加总和。 D(E)经验趋势，实际中D()按最小二乘法确定 D=AE+B Rez(E)正态分布的随机变量。,采用观测站的形式，可以用以前的AMP为基础，根据当前最大的AMPs(t)来预测将来岩体内释放的能量值E。 根据测量重力变化，可以确定应力集中区域，对冲击地压及震动危险区域进行定位。图8-5为STASZIC矿观测的重力异常RBGAi-1与趋势Ti-1的分布图。,图8-5 STASZIC矿观测的重力异常RBGAi-1与趋势i-1的分布图,由图可见，在3337之间存在AMPi-1的异常变化。说明该部分是卸压区。在第8系列，AMP出现了较大的增长，说明压力迅速增长。发生能量为E8107J的震动。而AMP9-1和AMP10-1表明，震动发生后产生的压力降。而后在AMP11-1表明，压力又回到原来的水平。而AMP14-1和AMP16-1的压力跳动是两次破坏顶板和放卸压振动炮的结果。其后高应力保持3周到第17循环的AMP17-1。,8.6 揭露老巷及其充填程度 地下的老巷和采空区对地表建筑物及工业区有很大的威胁。根据老巷与采空区上方重力场中的密度与周围岩体的密度之差，可以采用重力垂直梯度的分布确定其位置。 重力垂直梯度的分布定义为：,（8-7）,式中： z 垂直方向h重力测量水平之高zzw重力垂直梯度的导数，单位为etwesz（）,图8-6为两条巷道上方重力垂直梯度的分布，根据曲线zzw的极小值，可以很好的确定巷道的位置。 图8-7表示了在地表测量的巷道上方重力垂直梯度zzw的分布。这些巷道距地表13.5m。从图上可以清楚的看到，巷道上方的zzw出现最小值。 图8-8表示了保护煤柱上方重力垂直梯度的变化情况。zzw出现的最大值的情况与压力分布以及煤柱的破裂区域联系在一起。,图8-6 巷道上方12.0m处的重力垂直梯度分布,图8-7 地下巷道的定位（其中1为重力异常曲线，2为测点）,图8-8 煤柱与巷道上方zzw分布,根据Fajlclewicz的研究，zzw的负振幅与巷道中非充填部分的体积成比例关系，这样就可以计算开采后充填的程度。,8.7 揭露井筒周围岩层条件及变形 确定井筒周围岩层密度的方法称之为重力垂直剖面法PPGR（某个岩层的重力是五倍厚上覆岩层重力的叠加） 图8-9介绍了密度变化的五层岩层模型在有井筒的情况下，密度的变化。其中C部分中有空洞P环绕井筒周围，井筒半径为3.5m。S为该模型的垂直重力曲线。在这种情况下，空洞形成了曲线S的异常。井筒周围的空洞是岩石压出的根源，因此可以借助于PPGR的异常关系来确定和研究井筒周围的空洞的分布。经验表明压出过程中PPGR的负异常将增加其幅度，因此PPGRi+1PPGRi（1,2,3 ）是研究的基础。在同一点测量PPGR值。PPGR曲线之差,可以给出井筒附近压凸出的具体数量的程度。PPGRi-1(i1)曲线及垂直坐标轴的面积与岩石压出量成比例。,（8-8）,式中引力常数积分段（垂直方向）出现压出过程的区域宽度（未知）i-1=iYi次测量与第一次测量之间压出的岩石 量，i为线密度或用两个量相比，求出其压出量的相对增长率。,图8-9 井筒的模拟研究,。,图8-10为Maciej井的研究情况，井深199m，剖面为第四纪，三叠纪和石炭纪。第四纪，三叠纪为近水平，石炭纪为75SE，在3049和5778m有水的影响 测了四个系列的PPGR，时间间隔为69个月。从地表到81m，每隔3m设一个测点，深部每915m设一个例点，采用Worden Master重力仪，均方差控制在0.16ms-2。从PPGR2-1的曲线分布可以注意到，岩石的压出集中在第四纪和三叠纪，深度从地表到50m。PPGR3-1表明，压出很可能既在垂直也在水平方向上出现。考虑压出量的增加，则：,上式表明，在第二、三测量段压出量比上一段增加了30%。同样，在第四纪和三叠纪一、二测量段，其值为400%。 第四纪的压出强度为50%，而PPGR3-1表示，该动力过程只出现在下部分。但PPGR4-1表明。在第四次测量时，压出活动已降低。,图8-10 Maciej 井垂直重力剖面测量结果a井壁周围表面影响的地形重力修正值b井筒支护、邻近巷道影响的采矿重力修正值c地质剖面d各层密度与平均密度eBouger重力异常分布fPPGR随时间异常差的变化,8.8 微重力方法在地质灾害调查中的应用实例,8.8.1 概述 微重力测量是国际上世纪70年代形成的新技术，它具有精度高、分辨率高、快速、经济等优点，已在诸多领域获得广泛应用。 主要特点有两个： 一是探测对象小，测量范围小； 二是测量精度高，可达微伽级（10-8ms-2或10-9g）。,微重力测量和传统重力测量的主要差别: 最基本的差别在于测量精度提高了三个数量级，传统重力测量是毫伽级（10-5ms-2），微重力测量是微伽级（10-8ms-2）。测量方法和途径有很大发展。传统重力测量测点稀疏，而且只在地表进行，获得二维数据。微重力测量测点密集，可深入到竖井和多层坑道，甚至在物体内部进行，获得丰富的三维数据。 现今的数字模拟和反演技术是传统的重力勘探时期不可比拟的。如今几乎可以任意修正模型去逼近观测值，因而能得到更真实的解。,用微重力测量进行重力勘探之所以能找矿、发现空穴和解决地质问题，是基于各种岩矿、洞穴、构造与周围物质存在密度差异。在解释重力异常时岩矿密度又是解释的基础。因此在微重力测量中，岩矿密度的测定是一项重要工作，因此，在测量工作中，专门对岩矿密度测定的取样、测定和数据整理给以详细明确地规定。,重力加速度g是最基本的物理量之一，重力加速度g及其微小变化量g的精确测量，在生产、国防和科学实践中有极其重大的作用。 重力垂直梯度是地球重力沿铅垂方向的变化率，它在重力测量学和应用地球物理学的重力勘探中的地位是仅次于重力（加速度）的一个参数。研究这个参数的分布与变化，可以提供关于地球形状与地下密度异常体及其活动性的一些信息。,8.8.2 岩溶塌陷探测 广西桂林市某研究所生活区位于桂林市东郊、桂林孤峰平原、漓江河谷东侧的坡洪积台地上，原地貌呈垅岗缓坡地形，全区均为第四系沉积物粗盖。土层厚度一般大于15m，土层中除填土层外，主要为红黄、棕黄色粘土或亚粘砾（卵）石层；下伏基岩为上泥盆统桂林组（D3t）深灰色石灰岩，岩溶较为发育。该所处在一近南北向溶槽的东侧，属易发育土洞和塌陷的典型粗盖岩溶区。,该所1979年底土建工程基本结束。近年来由于大量抽吸岩溶地下水，改变了地下水的赋存状态，加速了土洞的发育。1991年5月的一深夜，一场大雨过后，在距该所锅炉房墙基不足lm处发生了直径约3m、坑深约3m的岩溶塌陷。为查明隐伏溶（土）洞的发育情况，搞清隐患以进行综合治理，在院内开展了点距为2m的微重力详查工作。,根据实测重力异常的分布规律，查出了四条近南北向延伸与隐伏岩溶发育有关的串珠状局部重力低异常带，这些串珠状的重力低异常正好反映了隐伏岩溶发育的不均匀性。在I号异常带上实施钻探工程以验证物探异常，在位于塌坑南约65m处（即37.5/5号测点）的ZK4孔，于深度23.433.92m深处见到10.52m厚土层中的软塑和流塑层（见图5）。此软层上覆有23m厚的粘土、粘土砾卵石持力层对大型建筑物的威胁是不可忽视的。,1粘土；2.粘土卵砾石；3软塑状粘土；4.石灰岩；5.溶洞；6.钻孔及编号 图5 桂林市某研究所第5测线物探地质综合剖面图,物探第10测线的g布剖面曲线在50/10号测点附近（即发生塌陷地段）出现了明显的重力低异常（约-40微伽），推断此地段应存在与溶（土）洞发育有关的低密度体。后在距48/10测点南约1.5m处（即塌坑西侧）进行钻探（ZK1孔），在地下33.1541.99m见到了8.84m高的半充填溶洞（图6）。看来物探异常与钻探结果吻合，物探先行，为准确定孔位、查明隐患提供了可靠的依据。,1.素填土；2.枯土；3.枯土卵砾石；4.石灰岩；5.溶洞；6塌坑；7.钻孔及编号图6桂林市某研究所第10测线物探地质综合剖面图,8.8.3 土洞探测 湖南省某化工厂1957年建成投产，由于工业、生活废水排放不合理，造成土洞发育，自八十年代初期以来，厂区地面塌陷时有发生，其中影响最大的一次是1983年12月突然发生的地陷，陷坑深约6.5m，面积2.8 4.5m2为了查明原漏勘和新发育的土洞，排除隐患，在区内进行了地面微重力工作。,微重力工作主要布置在原有重要建筑物附近和拟建厂房的建筑场地上。重力点距采用25m不等。物探结果不但在已知塌陷区发现明显的重力低异常，而且还发现了一批重力异常。后立即实施钻探工作，验证物探异常的属性，几乎在所有的重力异常区均分别打到了软弱、流塑层或土洞、溶洞。其中异常最强的是44/24 50/24号测点（约-35微伽）和56/50- 62/50号测点（约-30微伽）的两个重力低异常。见到了3.5m高和2.3m高的土洞（图8）。为排除隐患，提供了依据。,1.素填土；2.粘土卵砾石；3.土洞；4.石灰岩；5.钻孔及编号图7 湖南省某化工厂第24测线物探地质综合剖面图,1.粘土；2.软塑状粘土含砾石；3.流塑状粘土含砾石；4.土洞；5.石灰岩，6.溶隙图8 湖南省某化工厂第50测线物探地质综合剖面图,8.8.4 四川省攀枝花市某厂滑坡探测 该厂位于我国金沙江畔的斜坡地带，是我国山体滑坡易发区之一。1991年一场大雨过后，厂区旁的公路大范围山体滑坡。采用微重力方法，用以探测昔格达组地层与花岗岩、辉长岩的密度界面，从而达到间接寻找滑坡界面的目的。野外工作中顺地形等高线方向布置了三条重力剖面，点距为510m不等。实侧布格重力剖面曲线在山体滑动地段出现-320微伽的重力低异常，依据这些重力资料，采用微机反演昔格达组地层与火成岩的密度界面深度。,这一解释结果后来被钻探所证实，坡体外的K16孔（185/1号测点）花岗岩埋深仅有lm（图9）。通过综合分析密度分层断面图判断火成岩沟槽构造和地表、地下水对昔格达组地层的软化及润滑是该区山体产生滑坡的主要原因。认为滑坡体位于高密度、坚硬火成岩的沟槽部位，此沟槽便于地表、地下水的汇集，当上覆的昔格达组岩层被水软化、润滑到一定程度时，便形成滑动面而产生大面积的山体滑波。,1.粘土岩；2.滑坡堆积物；3.花岗岩；4.微重力异常推断的密度界面图9 四川省攀枝花市某厂第1测线物探地质综合创面图,docin/sanshengshiyuandoc88/sanshenglu,更多精品资源请访问,