重力勘探中场分离怎么做？趋势面法与三维正演的标定实践

1. 场分离这件事，先想清楚你想在重力资料里找什么

1.1 布格重力异常是“多源信号叠在一起”的结果

做重力资料解释的人，第一眼看到布格异常图，往往容易误判：一块高值区摆在那，好像直接就能画个密度体上去。实际上，布格重力异常几乎是地下所有密度不均匀体的综合响应，地表到地幔都贡献信号，你看到的那个“大高值”，很可能是几十公里尺度上的密度体造成的区域性背景，你的目标矿体只是背景场里一个不起眼的弯曲。

我处理某次矿区老资料时就有这种强烈感受。新做的1:1万重力测量，布格异常图上一眼望去只有个北西走向的隆起，像一个巨大平台，但把平台压平后，两侧确实有微弱的重力高。后来钻探证实在这个微弱高值下方400米深处有一个隐伏大理岩化和矿化体。问题在于，如果直接拿布格异常做反演或圈定剩余密度，模型会跑偏到深部区域场源上去；必须先把那些“不好看但真实存在”的大尺度区域场剥掉，再去解释浅部局部异常，而“剥掉”这一步就是位场分离。

区域场和局部异常没有严格物理分界，它只取决于研究尺度。你要找浅部矿体，那么埋深几公里、横向延伸十几公里以上的背景，就可以当作区域场；你要研究盆地基底起伏，可能就得把更深的地壳构造当作区域场。所以在我理解里，场分离本质是一种尺度滤波——把“对当前目标来说太宽缓的那些分量”从数据里拆出去。

1.2 趋势场法为什么能承担区域场提取？它是靠什么拟合的

趋势场法的核心思想很朴素：既然区域场形态宽缓、波长长，那它就适合用一个在平面上光滑变化的多项式曲面来逼近。想象拿出一张弹性薄膜，从布格异常顶部盖下去，膜不会跟着每一个小凹坑起伏，只反映大背景形状，膜的形态就是趋势面，而原始异常减去趋势面，剩下的就是剩余异常。

数学上，平面趋势面写成：

Z(x, y) = a₀ + a₁x + a₂y + a₃x² + a₄xy + a₅y² + …

x、y是测点平面坐标，Z是拟合值。一次多项式是倾斜平面，二次是抛物面，三次往后可以描述更复杂的宽缓形状。用最小二乘求系数，使所有测点上(Z - Δg)²的总和最小，这就是趋势场法的全部本质。

说“全部本质”可能让人觉得太简单，但它确实不复杂。为什么那么多人还在用？因为它稳定、快速、可解释。区域场到底长什么样没有先验信息时，用低阶光滑曲面去近似，比频率域高通滤波更直观，也不会像向上延拓那样需要一个合理延拓高度，更不会出现明显振铃。这套方法唯一风险在于——当你选的阶次太高，趋势面开始“刻画”局部异常细节，那剩余场就会被吃掉一部分甚至被反转。这个环节在第四部分我会用模型实验仔细展示。

2. 三维正演：先把“已知答案”做出来，再谈怎么分离

2.1 三维正演在整个流程里是什么角色

很多人会把重力正演当成独立的验证工具，用来“算算某个地质体能不能产生这么强的异常”。但我在做趋势场法研究时，反而更喜欢先用正演造一套“剧本模型”——地下密度结构是我设定好的，理论布格异常由计算机算出来，接着再对这套理论数据做场分离，看分离结果能不能还原当初设定的浅部局部异常。

这一步的价值很容易被忽略。真实重力资料是“未知谜底”，你做完分离只能凭地质经验判断效果好坏；而正演模型是所有解的“标准答案”，分离得好不好，直接拿结果和理论局部异常一对比就知道。所以正演不是和场分离并列的两件事，而是后者的“测试平台”和“标定工具”。三维正演在这里不是跑一个复杂软件，而是把地下简化成若干密度单元，算出地表每个测点的重力响应。

2.2 模型剖分与单元密度设计

三维正演常用做法是长方体单元离散。把地下空间按X、Y、Z方向切成规则网格，每个小长方体块赋予一个剩余密度值。剩下就是沿所有测点叠加每个块的垂向重力效应，输出理论布格异常。