1. 课程导入与任务总览

重磁勘探通过观测地表的重力场或磁场变化,推断地下密度或磁化率的不均匀分布。然而,地表实测异常是地下所有深度场源叠加的结果——浅部小尺度地质体产生的高频、窄幅异常与深部大尺度构造产生的低频、宽幅异常混叠在一起,直接解释往往难以区分。本课程要解决的核心问题就是:如何通过向上延拓压制浅部干扰、突出深部场源,并通过边界增强技术锐化地质体边界,使重磁异常从"混叠信号"变为"可解释图像"。

向上延拓的物理本质是位场在无源空间中的解析延拓。根据拉普拉斯方程,位场在无源区域满足调和函数性质,因此可以通过频率域滤波实现延拓:将实测异常做傅里叶变换,乘以延拓算子 \(e^{-kz}\)(\(k\) 为波数,\(z\) 为延拓高度),再反变换回空间域。延拓高度越大,高频成分衰减越剧烈,浅部小尺度场源的影响被压制,深部区域场得以凸显。边界增强则走的是另一条路径——通过计算水平导数、垂直导数或总水平导数(THD)、倾斜角(Tilt)等,放大场源边缘处的梯度变化,使密度或磁性界面的位置更加清晰。

向上延拓与边界增强并非互相替代,而是互补关系:延拓用于"去浅取深"、分离区域场与局部场;边界增强用于在延拓后的场上进一步"锐化边界"。实际处理中通常先延拓、后增强,顺序颠倒会导致浅部边界被过度放大。

本课程面向的实际应用场景包括:区域地质调查中分离沉积盆地基底与上覆盖层异常;矿产勘查中识别隐伏岩体边界与断裂构造;工程物探中压制浅表不均匀体干扰、突出目标层位;以及深部构造研究中通过多高度延拓分析场源的垂向分布特征。这些场景的共同需求是:从叠加场中提取有意义的地质信息,而向上延拓与边界增强正是实现这一目标的基础工具链。

完整的工作流程可以概括为五个环节:数据准备与网格化 → 延拓高度选择与计算 → 延拓结果质量检查 → 边界增强算子应用 → 综合解释与成果输出。数据准备阶段需要将离散测点数据网格化为规则网格,检查网格间距、边界效应和噪声水平;延拓高度选择需要结合目标场源深度与数据分辨率,通常从较小高度(如 1~2 个网格间距)开始逐级试验;质量检查关注延拓后异常是否出现边界振荡、幅值是否合理衰减;边界增强阶段根据目标类型选择 THD、Tilt 或 θ 图等算子;最终将延拓场与增强结果叠合分析,圈定地质体边界与构造走向。

重磁异常向上延拓与边界增强 · 任务地图 原始观测数据 重力/磁力测点 · 坐标 · 高程 噪声水平 · 测网密度 网格化与预处理 规则网格化 · 网格间距选择 边界扩边 · 去趋势/去噪 频谱分析与参数设定 径向功率谱 · 场源深度估计 延拓高度初选 · 算子选择 关键输入 网格间距 Δ 延拓高度 z 向上延拓计算(频率域) F[ΔT_up] = F[ΔT] · e^(−k·z) k = 2π√(kx²+ky²) 为径向波数,z 为延拓高度 高频衰减快 → 压制浅部小尺度场源 低频保留 → 突出深部区域场与大型构造 延拓高度分级试验 z₁ = 1Δ(浅层压制) z₂ = 3Δ(中深分离) z₃ = 5Δ(区域场提取) Δ = 网格间距 延拓质量检查 ✓ 幅值合理衰减 ✓ 无边界振荡 ✗ 过延拓 → 信息丢失 ✗ 欠延拓 → 浅部残留 边界增强算子 总水平导数 THD = √(∂f/∂x)²+(∂f/∂y)² 倾斜角 Tilt = arctan(VDR / THD) θ 图 · 垂直导数 · 水平梯度 目标:锐化密度/磁性界面边界 综合解释与成果输出 延拓场 + 增强结果叠合分析 圈定地质体边界 · 断裂走向 划分区域场与局部异常 输出:边界图 · 构造纲要图 · 解释剖面 处理顺序原则 ① 先延拓 → 分离深浅场源 ② 后增强 → 锐化目标边界 ③ 顺序颠倒 → 浅部边界过度放大 延拓高度与增强算子需匹配目标深度

课程任务地图按"数据准备—延拓计算—质量检查—边界增强—综合解释"五段组织,但各段之间并非简单串联:延拓高度选择需要参考频谱分析得到的场源深度估计,边界增强算子的选择又取决于延拓后场的特征。例如,若延拓后区域场占主导、局部异常已很微弱,THD 可能放大噪声,此时 Tilt 或 θ 图因具有归一化特性而更稳健。这种"参数联动"关系是后续各章节反复强调的重点。

为了让学习过程有明确的检验标准,本课程在每个实操环节都设置了可观察的结果判据。以延拓为例:合理的延拓结果应表现为异常幅值随高度增加而单调衰减、异常形态趋于圆滑、边界处不出现明显振荡;若延拓后异常出现"环状"假异常或幅值反常增大,通常意味着扩边不足或网格化参数不当。边界增强的判据则是:目标地质体边界处梯度极大值清晰、不同走向构造的边界响应可区分、背景噪声未被过度放大。这些判据将在后续章节结合具体算例逐一展开。

学习建议:本课程涉及频率域滤波与空间域求导两类运算,建议在进入实操章节前,先确认自己熟悉傅里叶变换的基本性质(尤其是卷积定理与频率域滤波的对应关系)以及方向导数的几何含义。这两项预备知识直接决定延拓算子与边界增强算子的理解深度。

从任务分工看,向上延拓解决的是"看多深"的问题——通过选择不同延拓高度,实现场源垂向分层;边界增强解决的是"看多清"的问题——通过梯度算子放大界面响应。两者结合,才能从重磁异常中同时获得深度信息与边界信息。本课程后续将依次展开:位场延拓的数学基础与频率域实现、延拓高度的定量选择方法、各类边界增强算子的原理与适用条件、噪声与边界效应对处理结果的影响机制,以及基于实际工区数据的完整处理案例。每一章都围绕"可操作、可检验"展开,最终使学习者能够独立完成从原始网格到边界解释图的全流程处理。