1. 课程导论与重力勘探基础
重力勘探是一种以地球重力场为观测对象的地球物理方法。它通过在地表或空中测量重力加速度的微小变化,反推地下岩石密度的横向不均匀性,从而圈定构造、识别岩体、寻找与密度差异有关的矿产资源与地下构造。要真正掌握这门方法,关键不在于记住几个公式,而在于理解一条完整的逻辑链:密度差异 → 重力场扰动 → 观测数据 → 各项改正 → 布格异常 → 反演解释。本课程围绕这条链条展开,最终目标是让学习者能够独立完成从原始重力数据到布格异常图、再到简单密度反演的完整流程。
一、重力勘探的基本原理
牛顿万有引力定律告诉我们,任何有质量的物体都会对单位质量产生引力。地球是一个近似旋转椭球体,其表面重力加速度约为 9.78~9.83 m/s²,随纬度升高而增大、随海拔升高而减小。实际观测到的重力值 g 是地球总质量引力、自转离心力以及地下局部密度不均匀体共同作用的结果。重力勘探关心的正是其中“局部”那一部分——通常只有毫伽(1 mGal = 10⁻⁵ m/s²)量级,甚至更小。现代相对重力仪的分辨率可达 0.01 mGal,足以捕捉埋深数百米至数公里、密度差 0.1~0.5 g/cm³ 的地质体所产生的微弱信号。
这里必须先明确一个贯穿全课程的符号约定:本课程采用国际大地测量与地球物理联合会(IUGG)约定,重力方向向下为正,重力异常定义为“观测重力值减去同一点的理论正常重力值”。因此,密度偏高的地质体(如基性岩体、金属矿体)会产生正异常,密度偏低的构造(如盐丘、松散沉积盆地、断裂破碎带)会产生负异常。全文与所有图示均遵循这一约定,后续章节不再重复说明。
二、重力异常的定义与分类
“重力异常”是一个层层剥离的概念。最原始的是绝对重力异常,即观测值与正常重力之差,但它混杂了纬度、高程、地形和中间层质量的影响,不能直接用于地质解释。于是我们按改正顺序逐级定义:
- 自由空气异常:仅做高度改正,把观测点归算到大地水准面,保留地形质量的影响;
- 布格异常:在自由空气异常基础上,再扣除观测点与大地水准面之间“无限大平板”状中间层的引力,并做地形改正;
- 均衡异常:进一步考虑地壳均衡补偿的影响,多用于地壳尺度研究。
此外,按异常的空间尺度还可分为区域异常与局部异常;按场源性质可分为位场异常与梯度异常。理解这些分类的关键,是明白每一种异常对应“扣掉了什么、保留了什么”,而不是死记名称。
三、布格异常在勘探中的核心作用
布格异常是本课程的主线。它之所以成为勘探中最常用的异常形式,是因为它把与地下密度分布无关的纬度、高程、中间层和地形效应尽可能剥离,使剩余场主要反映地下横向密度不均匀性。换句话说,布格异常图就是一张“地下密度差异的投影图”。在工程实践中,布格异常可用于:识别盆地基底起伏与沉积厚度、圈定隐伏岩体与断裂、评价盐丘与火山岩分布、为密度反演提供输入数据。下面的信号链图展示了从野外观测到布格异常、再到反演解释的完整数据流,也标出了每一步的关键改正项。
需要强调的是,布格异常并非“纯净”的地下密度信号。它仍包含区域场与局部场的叠加、深部与浅部场源的耦合,以及观测噪声。因此,布格异常计算只是解释工作的起点,后续还需进行异常分离、延拓、导数换算等处理,才能进入反演环节。
四、课程整体学习路径与实战目标
本课程共 30 章,按“基础—改正—计算—反演—实战”的逻辑递进组织。前段建立重力场与异常的基本概念,明确符号与坐标约定;中段集中讲解各项改正的原理与实现,重点是自由空气改正、布格改正、地形改正的计算细节与常见错误;后段进入布格异常图的编制、异常分离与密度反演,最终以一个完整的示例数据集串联全流程。
为了让学习目标可检验,这里给出一个贯穿课程的示例任务(数据为教学假设,非实测):某测区面积约 20 km×15 km,测网 500 m×500 m,共约 1200 个测点,相对重力仪观测精度 0.02 mGal,测点高程由 GNSS 测得、精度 ±0.1 m。中间层密度取 2.67 g/cm³。若某测点高程为 350 m、纬度为 30°,正常重力约 9.793 m/s²,观测值为 9.792 500 m/s²,则自由空气改正约 +108.0 mGal,布格改正约 −39.2 mGal,可据此得到该点的布格异常值。学员应能独立复现这一计算,并判断结果是否落在合理区间。
完成本课程后,学员应达到三个实战目标:其一,能根据测区条件正确选择改正项与参数,独立完成布格异常计算;其二,能识别并排查常见错误,如高程基准混用、地形改正符号错误、密度取值不当导致的系统偏差;其三,能对布格异常图作出基本地质判读,并完成一次简单的密度反演,给出合理的解释结论。这三项目标将贯穿后续每一章的训练。