4. 矿区InSAR干涉图判读基础

干涉图是InSAR形变监测中最直观、信息密度最高的中间产品。它把两景SAR影像的相位差以条纹形式呈现出来,一条条纹对应一个完整的2π相位循环。对矿区监测而言,能否从干涉图上快速判断“哪里在动、动得多快、哪里不可信”,直接决定了后续解缠与形变反演的成败。本节围绕干涉条纹含义、条纹密度与形变梯度关系、失相干识别、矿区典型特征和干涉图质量评价五个方面展开。

一、干涉条纹到底表示什么

干涉相位由平地相位、地形相位、形变相位、大气延迟相位和噪声相位叠加而成。经过差分处理(去除参考DEM地形相位和平地相位)后,残余条纹主要反映视线向形变。基本关系为:

Δφ_def = -4π/λ · Δd_LOS
一条条纹(2π)对应 Δd_LOS = λ/2

以C波段(λ≈5.6 cm)为例,一条完整条纹约对应2.8 cm的视线向形变量;L波段(λ≈23.6 cm)一条条纹约对应11.8 cm。这意味着同一形变量在C波段干涉图上条纹更密、更易失相干,而L波段条纹更稀疏、穿透植被能力更强。判读时先确认波长,才能把“条纹数”换算成“形变量”。

条纹的疏密、走向和连续性都携带信息:条纹密集表示形变梯度大;条纹走向大致垂直于形变梯度方向;条纹突然中断或呈碎斑状,往往意味着失相干或相位突变。需要注意,相位存在2π模糊,单条条纹无法区分“动了2.8 cm”还是“动了2.8 cm加一个整周”,这正是后续解缠要解决的问题。

二、条纹密度与形变梯度的定量关系

条纹密度通常用“每公里条纹数”或相邻条纹间距描述。设视线向形变梯度为 g(cm/km),波长为 λ(cm),则条纹密度近似为:

N ≈ 2g / λ   (条/km)

举例:C波段(λ=5.6 cm)下,若某沉陷盆地边缘梯度为 5 cm/km,则 N≈2×5/5.6≈1.8 条/km,条纹间距约550 m;若梯度升至 20 cm/km,N≈7.1 条/km,间距缩至约140 m。当条纹间距接近或小于像元尺寸(Sentinel-1约20 m×5 m)时,相位在相邻像元间发生混叠,解缠必然失败。因此判读时要先估算梯度,再判断该干涉对是否“可解缠”。矿区沉陷盆地边缘正是梯度最大的区域,条纹最密,也是最容易出问题的地方。

干涉条纹密度与形变梯度关系示意(C波段 λ=5.6 cm) 缓梯度区(盆地内部) 梯度 ≈ 2 cm/km,条纹间距 ≈ 1.4 km 条纹稀疏,解缠可靠 陡梯度区(沉陷盆地边缘) 梯度 ≈ 20 cm/km,条纹间距 ≈ 140 m 条纹密集,接近混叠,解缠风险高 梯度增大

三、失相干区域的识别

失相干表现为干涉图上相位随机、呈颗粒状噪声、无规则条纹的区域,相干系数通常低于0.3。矿区常见的失相干来源有四类:一是植被覆盖与季节变化,采空区上方复垦植被生长快,两景间隔过长时相干性骤降;二是地表大形变本身造成的失相干,梯度超过临界值后相位完全随机;三是水体与积水塌陷坑,镜面反射导致几乎完全失相干;四是开采活动引起的快速地表改变,如排土场、运输道路的频繁扰动。

识别方法上,先看相干图(Coherence Map),把相干系数低于阈值的区域掩膜掉;再对照光学影像或历史开采边界,判断失相干是“自然植被”还是“形变过快”。前者可通过缩短时间基线、选择冬春季影像改善;后者则说明该区域形变已超出该干涉对的测量能力,需要改用短基线集或L波段数据。判读时切忌把失相干区误判为“无形变”,它恰恰可能是形变最剧烈的地方。

四、矿区典型干涉图特征

矿区干涉图有几个可辨识的“签名”。第一,沉陷盆地呈同心椭圆或“牛眼”状条纹,中心条纹最密、向外逐渐稀疏,条纹闭合圈数反映累计沉陷量。第二,工作面推进方向上的条纹呈不对称分布,超前影响范围条纹稀疏、滞后压实区条纹相对密集。第三,多工作面叠加时出现“8”字形或串珠状条纹组合,需结合开采时序区分。第四,断层或地裂缝附近条纹出现错断、相位跳变,是判断地表破裂的重要线索。第五,排土场、边坡区域条纹杂乱且相干性低,常与植被和持续蠕变叠加。

实际判读时,建议把干涉图与开采平面图、工作面推进时间表叠加对照:条纹中心是否落在采空区上方、条纹闭合圈数与预计沉陷量是否量级一致、边缘条纹密度是否与预计梯度吻合。三者一致,判读结论才可靠。

五、干涉图质量评价指标

评价一张干涉图能否用于形变解缠,常用以下指标综合判断:

指标含义参考判据(示例)
平均相干系数研究区相干均值>0.5 良好;0.3~0.5 可用;<0.3 慎用
有效像元比例相干高于阈值的像元占比>70% 良好
条纹可见度条纹连续、无碎斑目视连续可追踪
残余平地/地形相位差分后残余趋势无明显线性斜坡
解缠误差闭合环/残差点数量残差点少、无大面积跳变
时间/空间基线干涉对几何条件时间基线尽量短、垂直基线小

这些指标应结合使用:平均相干高但局部失相干面积大,仍可能在解缠时产生区域误差;条纹可见但残差点密集,说明噪声水平偏高。工程上通常先做相干掩膜,再统计有效像元比例,最后用残差点密度和闭合环检查解缠可靠性。只有指标与目视判读相互印证,干涉图才具备进入下一步形变反演的条件。

判读口诀:先看波长定标尺,再看条纹估梯度,相干掩膜辨真假,对照开采验结论。失相干不等于无形变,条纹密不等于形变大,必须结合波长与梯度换算后再下判断。