2. InSAR形变监测基本原理

InSAR（合成孔径雷达干涉测量）之所以能从太空“量出”矿区地表毫米级的形变，核心在于对两景SAR影像的相位差进行精密解读。但干涉相位并非只包含形变信息，它是一组物理量叠加的结果。理解这套相位的组成与分离逻辑，是判读矿区形变数据的第一道门槛。

干涉相位的组成

对同一地区在t₁、t₂两个时刻获取的两景单视复数影像逐像元共轭相乘，得到干涉相位Δφ。在忽略噪声与失相干的前提下，它可分解为四项：

• 平地相位（参考相位）：由地球椭球面几何与成像几何决定，即使地表完全平坦也会存在，量级可达数千弧度，必须去除。
• 地形相位：由地表高程起伏引起，与垂直基线成正比。通常借助外部DEM（如SRTM 30 m）或干涉图自身生成DEM来模拟并扣除。
• 形变相位：视线向（LOS）位移在雷达波传播方向上的投影，正是我们关心的目标信号。
• 大气延迟相位：电磁波穿过对流层、电离层时因折射率变化产生的附加延迟，是矿区短基线监测中最主要的误差源之一。

此外还有轨道误差相位与热噪声。工程上常把“平地+地形”合称几何相位，把“大气+轨道+噪声”合称误差相位，于是有：

Δφ = φ_flat + φ_topo + φ_def + φ_atm + φ_orb + φ_noise
φ_def ≈ Δφ − φ_flat − φ_topo − φ_atm − φ_orb

需要强调，干涉相位是缠绕在(−π, π]区间内的，上述各项都是“缠绕后”的相位，直接相减并不严格成立，这正是后续需要相位解缠的原因。

形变相位提取与相位解缠

去平地、去地形后得到的差分干涉图，其相位仍被包裹在2π周期内。相位解缠的任务，就是为每个像元确定一个整数2π模糊度，把缠绕相位还原为连续相位。解缠质量直接决定形变结果的可靠性，矿区常见的解缠失败往往表现为大范围“跳变”或整片区域被错误抬高/压低。

解缠后的相位φ_def与LOS形变d_LOS的关系为：

d_LOS = −(λ / 4π) · φ_def

式中λ为雷达波长（C波段约5.6 cm，L波段约23.6 cm）。负号源于相位随距离增大而减小的约定：目标朝向卫星运动（抬升或向视线方向靠近）时，斜距减小，相位增大，d_LOS取负值。判读时务必先确认所用软件的正负号约定，全文与图示保持一致。

视线向形变与地表位移的关系

SAR只能测到沿视线向的一维分量，而矿区地表位移通常是三维的。设地表位移为(d_E, d_N, d_U)，入射角θ，卫星飞行方位角与东向夹角为α（ ascending/descending 不同），则：

d_LOS = d_U·cosθ − sinθ·(d_E·sinα − d_N·cosα)

由于θ通常为30°~45°，垂直分量权重最大，水平分量被压缩。单一轨道难以分离三维位移，工程上常用升降轨联合观测，或结合GNSS、水准点约束反演。对矿区而言，沉陷盆地以垂直下沉为主，LOS形变与下沉量近似满足 d_LOS ≈ d_U·cosθ，这也是把LOS形变近似当作下沉量估算时需注明的前提。

矿区常见相位失真来源与判读要点

矿区场景有其特殊性，以下几类失真在判读时必须逐项排查：

• 失相干：植被覆盖、农田耕作、采动引起的地表剧烈破坏，会使干涉相干性骤降，表现为干涉图上“雪花状”噪声。判读时应结合相干图掩膜，相干性低于约0.3的区域不宜直接采信形变值。
• 大气延迟：山谷地形与昼夜温差会引入数厘米量级的假形变，空间上呈平滑长波。可通过多景时序平均、GACOS等大气改正产品或与GNSS比对识别。
• 大梯度解缠错误：沉陷盆地边缘形变梯度大，相邻像元相位差超过π时易发生解缠跳变，表现为形变场出现不自然的“台阶”或整片偏移。应检查解缠残差与闭合相位。
• 轨道与DEM残余：轨道误差呈线性斜坡，DEM误差与地形相关，二者都会在差分图上留下规则条纹，可通过拟合残余斜坡或使用更高精度DEM抑制。

综上，InSAR形变监测的本质是一条“相位分解—误差剔除—解缠还原—几何投影”的信号链。只有清楚每一项相位从何而来、如何被扣除、残余误差如何识别，才能把干涉图上的条纹读成矿区地表真实的沉陷信息，而不是被大气、失相干或解缠错误误导的假象。