布格重力异常地形校正方法在微重力勘探中的缺点及其纠正方法.docx

布格重力异常地形校正方法在微重力勘探中的缺点及其纠正方法高庆余(中国煤田地质总局物探研究院)【摘要】多年来,不论重力勘探程度如何,在布格重力异常计算中都必须经过地形改正和中间层改正。本文通过分析在计算布格重力异常时地改和中间层改正对测点的重力补偿,提出了取消中间层改正以适应微重力勘探精细解释需要的地形校正方法。该方法建立在对实际地形(岩性)的正演基础上,可以根据施工地区的地质条件合理选择重力基准面进行可变密度地形校正。使用该方法可以比较好地消除地形起伏和不均匀岩性对测点产生的重力影响,从而得到比较可信的重力异常数据。【关键词】微重力勘探地形校正重力补偿【中图法分类号】P631.14近年来,随着重力仪观测精度的大幅提高和普及使用,利用微重力勘探研究弱异常地质构造成为可能,因此引起了人们的高度重视。但是,由于野外地质(表)条件的复杂性,实际对重力数据的地改和中间层改正误差常常大于地质体引起的重力异常。所以,要收到满意的勘探效果,还必须辅之以高质量的地形校正。多年来,我们一直使用?g=?g观+?g正+?g中+?g高+?g地计算出的数据,即布格重力异常来衡量地下地质体在地面引起的重力变化。这种计算方法在地改和中间层改正中只能使用平均(或综合)密度,所以在自然条件下(岩性不均匀)布格重力异常数据的可靠程度较低,不利于对弱异常做出合理的推断和解释。本文在分析布格重力异常地改和中间层改正基础上提出了取消中间层改正的地形校正方法。该方法可以在最大程度上克服密度使用不当造成的地形校正误差,在相对不利的地质(表)条件下提高勘探质量,得到容易评价而又比较可靠的重力异常数据。地形校正的理论基础12从牛顿的万有引力定律F=GM13M2/R可知,任一质量单元dm对相距r的某点单位质量产生的引力与物质质量dm成正比,与物质到该点距离r的平方成反比dF=Gdm/r2。收稿日期:1998-06-24;改稿日期:1998-11-16地址:河北涿州中国煤田地质总局物探研究院(邮政编码:072750)4期高庆余:布格重力异常地形校正方法在微重力勘探中的缺点及其纠正方法·321·按此定律,一个质量为M的较大物质对该点单位质量产生的引力垂直分量(即重力或重力场强度)为质量M中所有微小质量单元dm在该点产生的引力垂直分量之和dm?g=G∫r2cosΗ其中Η为r与计算点铅垂方向间的夹角如果质量为M的物质密度Θ均匀,则上式可写成(1)dv?g=GΘ∫r2cosΗ(2)为了消除由于所处环境不同周围物质给各测点重力读数造成的不同重力效应,布格重力异常在地改和中间层改正中使用平均密度,即利用式(2)作为地改和中间层改正的基本公式。2布格重力异常地形校正方法的若干缺点2.1地形校正方法按照通常计算布格重力异常时的地形校正方法,为了把各点重力值归算到统一的基准面上,必须首先进行地形改正,把高于测点地形部分实际存在的物质“铲平”,同时把低于测点的物质缺失部分即自由空间用同一密度物质“填满”,使测点相当于处在上部无物质(引力影响),下部又完全被同一密度物质填充满的无限延伸平面上。由于地改后各测点所处平面比重力基准面多了一层(测点比基准面高)或少了一层(测点比基准面低)无限平板物质对各测点的不同引力影响,所以地改后必须进行中间层引力改正(由于在中间层中包含地形起伏部分物质,所以在本文中我们把中间层改正纳入地形校正范畴)。2.2岩性均匀条件下重力基准面对地形校正精度的影响通常,在某一地区开展的重力勘探中,为了计算方便,一般把重力总基点水平面作为进行中间层改正(或布格重力异常)计算的重力参考基准面。在重力勘探理论(和实际工作)中,对重力总基点的所在高程没有特殊限制。我们假设图1中S为某一自然起伏地形表面,A为此表面上任一测点,L1、L2分别为A点水平面和重力参考基准面。图1测点比基准面高的地改和中间层改正重力补偿示意图(1)测点比基准面高对于某一固定高程的基准面来说,当测点A比基准面高(如图1(a)所示)时,A点重力值经过地改后,A点高度以上S地表面与A点水平面L1包围的物质V0部分和A点水平面L1以下自由空间(图中V1+V4+V3+V5)人为“填满”物质都使A点重力值增加(图中用“+”表示)。根据式(2),A测点得到的重力补偿为dvdv?g=?gV0+?gV1+?gV4+?gV3+?gV5;?gV0=GΘ∫Vr2cosΗ;?gV1=GΘ∫V1r2cosΗ0以下类推·322·物探化探计算技术20卷由于A点重力值经过地改后比基准面多了V1+V2+V3中间层物质引力贡献,因此,在进行中间层改正时应该减去这部分物质引力影响(图1(b)中用“-改正后,A点得到的实际重力补偿为?g=?g地+?g中=”表示)。经过地改和中间层(?gV0+?gV1+?gV4+?gV3+?gV5)-(?gV1+?gV2+?gV