3. 向上延拓的物理意义与数学基础

向上延拓是把观测面（或某一高度面）上的位场值，换算到该面之上更高一个平面的处理过程。它回答的是一个很朴素的问题：如果我把仪器抬到空中某个高度再测一次，磁场或重力场会变成什么样？从物理上看，向上延拓相当于把观测面“推远”，让深部、规模大、横向延伸广的场源在总场中占比上升，而浅部小尺度异常被迅速抹平。因此它既是一种位场变换手段，也是一种天然的“低通滤波”和“深源增强”工具，常被用于压制浅部干扰、突出区域背景场、为向下延拓和求导提供稳定的中间结果。

调和函数性质：延拓为什么可行

在无源空间（场源之外）中，重力位和磁位都满足拉普拉斯方程：

∇²U = ∂²U/∂x² + ∂²U/∂y² + ∂²U/∂z² = 0

满足该方程的函数称为调和函数。调和函数有一个关键性质：它在区域内部的值完全由边界上的值决定，这就是位场延拓的理论根基。给定一个平面上的场值，理论上可以唯一地确定无源半空间内任意一点的场值，向上延拓和向下延拓都是这一性质的具体应用。

调和函数还带来两个实用推论。其一，位场在无源空间中是无限次可微且解析的，不存在真正的“突变”，所以延拓结果总是光滑的，不会凭空产生新的高频细节。其二，调和函数满足极值原理：在无源区域内，场值不会出现新的局部极大或极小，极值只能出现在边界上。这解释了为什么向上延拓后异常形态趋于平缓、峰值降低、范围展宽——它是在“抹平”而不是“制造”结构。

频率域表达式：延拓的数学核心

在空间域，向上延拓是一个卷积运算，核函数随高度增大而展宽。但真正便于计算和理解的是频率域表达。对平面场做二维傅里叶变换后，设波数为 k = √(kx² + ky²)，延拓高度为 h（向上为正），则向上延拓的频率域响应为：

U(kx, ky, z = -h) = U(kx, ky, 0) · exp(-k·h)
其中 k = √(kx² + ky²)，h > 0 表示向上延拓

这个指数因子 exp(−kh) 就是向上延拓的全部“秘密”。它是一个随波数增大而迅速衰减的低通滤波器：低频（小 k，对应深部大尺度场源）衰减慢，几乎被保留；高频（大 k，对应浅部小尺度场源）衰减快，被强烈压制。延拓高度 h 越大，截止效应越明显，结果越平滑。

实际计算流程通常是：对网格化数据做二维 FFT → 乘以 exp(−kh) → 反变换回空间域。需要注意两个工程细节：一是要合理选择网格边长和补零（padding），抑制周期性和边缘效应；二是高频端噪声会被压制，但若原始数据含长波长误差（如区域背景倾斜），向上延拓无法去除，反而会保留甚至相对放大其占比。

延拓高度与分辨率损失的关系

从 exp(−kh) 可以直接读出分辨率损失的定量规律。设某异常的特征波长为 λ，对应波数 k = 2π/λ，则其幅值经高度 h 延拓后被乘以 exp(−2πh/λ)。当 h = λ/(2π) ≈ 0.16λ 时，幅值衰减到约 37%；当 h = λ/2 时，衰减到约 4%。换句话说，延拓高度接近异常波长的一半时，该异常基本被“延拓掉”。这给出一个非常实用的经验判据：

• 若希望保留某尺度异常，延拓高度应明显小于其波长的一半，例如 h ≤ 0.2λ；
• 若希望压制某浅部干扰，延拓高度应达到或超过其波长的一半，例如 h ≥ 0.5λ；
• 延拓高度每增加一个“特征波长量级”，可分辨的最小横向尺度大致同比例增大，分辨率损失与 h 近似成反比关系。

用一个教学示例说明（以下为合理假设的示例参数，非实测结果）：设网格点距 100 m，某浅部磁性体引起的异常半波长约 400 m（波长约 800 m）。若向上延拓 h = 200 m，该异常幅值保留约 exp(−2π×200/800) ≈ 0.21，即衰减近 80%，基本被压制；而一个波长约 4000 m 的深部区域异常，幅值保留约 exp(−2π×200/4000) ≈ 0.73，仍清晰可辨。若把 h 提高到 600 m，浅部异常保留不足 1%，区域场保留约 0.39，此时图面已非常平滑，适合提取区域背景，但不宜再用于刻画小构造。

实践中的判读与参数选择

实际作业时，建议先对原始网格做频谱分析，观察高频噪声与有效信号的分界波数，再据此估算合适的延拓高度。计算时对数据补零到 2 的整数次幂边长，减少 wraparound 效应；输出后重点检查三点：一是浅部窄小异常是否按预期消失，二是区域场走向是否与地质认识一致，三是测区边缘是否出现明显的人工梯度。若边缘畸变严重，可先做扩边或余弦镶边再延拓。向上延拓的结果通常作为中间产品，与向下延拓、水平导数、解析延拓组合使用，才能兼顾深部背景与浅部细节的刻画需求。