IBIS-L地基InSAR系统是地基InSAR系统的一种, 它将步进频率连续波技术、SAR技术以及InSAR技术融合在一起来获取目标的地形/形变信息。该系统的工作模式一般采用连续模式对工作区进行不间断观测。

地基InSAR形变测量示意如图1所示。图1中设置固定观测基站观测目标点 P$\begin{array}{c}P\end{array}$, P'$\begin{array}{c}\mathit{P\text{'}}\end{array}$为 P$\begin{array}{c}P\end{array}$移动了距离 d$\begin{array}{c}d\end{array}$后的位置。形变前后的主副图像相位分别为 φM$\begin{array}{c}{\phi }_M\end{array}$和 φS$\begin{array}{c}{\phi }_S\end{array}$, 两者的干涉相位 ΔφMS$\begin{array}{c}\Delta {\phi }_{\mathit{MS}}\end{array}$可表示为

ΔφMS=φS−φM=4π(SP'−MP)λ=4πλd$\begin{array}{c}\Delta {\phi }_{\mathit{MS}}={\phi }_S-{\phi }_M=\frac{4\pi (\mathit{SP\text{'}}-\mathit{MP})}{\lambda }=\frac{4\pi }{\lambda }d\end{array}$, (1)

式中: MP$\begin{array}{c}\mathit{MP}\end{array}$和 SP'$\begin{array}{c}\mathit{SP\text{'}}\end{array}$分别表示传感器到 P$\begin{array}{c}P\end{array}$和 P'$\begin{array}{c}\mathit{P\text{'}}\end{array}$的距离; λ$\begin{array}{c}\lambda \end{array}$为雷达波长, 由式(1)可得目标点 P$\begin{array}{c}P\end{array}$的形变量d为

d=λ4πΔφMS$\begin{array}{c}d=\frac{\lambda }{4\pi }\Delta {\phi }_{\mathit{MS}}\end{array}$, (2)

式(2)是在空间基线为0, 且不考虑干扰项的情况下得到的形变理论关系式。但在非连续观测模式下, 空间基线可能不等于0。因此, 干涉相位的完整表达为

ΔφMS=φGeom+φDefo+φAtmo+φNoise+2kπ$\begin{array}{c}\Delta {\phi }_{\mathit{MS}}={\phi }_{\mathit{Geom}}+{\phi }_{\mathit{Defo}}+{\phi }_{\mathit{Atmo}}+{\phi }_{\mathit{Noise}}+2\mathit{k\pi }\end{array}$, (3)

式中: φGeom$\begin{array}{c}{\phi }_{\mathit{Geom}}\end{array}$为地形相位; φDefo$\begin{array}{c}{\phi }_{\mathit{Defo}}\end{array}$为形变相位; φAtmo$\begin{array}{c}{\phi }_{\mathit{Atmo}}\end{array}$为大气延迟相位; φNoise$\begin{array}{c}{\phi }_{\mathit{Noise}}\end{array}$为噪声相位。

	Figure Option
	图1 地基InSAR形变测量示意图[2]Fig.1 Sketch of GB-InSAR for deformation measurement[2]

充分收集观测区域基础地理数据、地质数据和遥感数据等资料, 通过野外实地勘察, 掌握边坡位置、交通、地质环境和结构特征等。观测方案设计应充分考虑以下内容:

(1)持续供电。为了保证长时间的连续监测, 地基InSAR系统需要持续供电。该套设备的备用电池仅能维持10 h左右, 在长时间的连续监测过程中, 应采用交流电持续供电。所以, 选址时应选择离电源比较近的地方, 或采用发电机供电。

(2)交通便利。地基InSAR系统设备是精密而又贵重的仪器, 在搬运过程中要轻拿轻放。由于该套设备体积大且较重, 在选址时应该选择交通比较便利或者容易搬运的地方。

(3)监测距离和范围。在实际情况中, 应根据地形地貌条件将地基InSAR的监测距离控制在合理的范围内, 即最远监测距离不能超过4 km, 最大覆盖范围不能超过7 km²。在此前提下, 监测距离越大, 误差越大, 监测效果越差, 所以应该尽量缩短监测距离。

(4)通视条件。地基InSAR应该在通视条件下工作, 在工作点与监测区域之间不能有障碍物存在。如果存在障碍物, 监测对象的反射强度会受到影响, 影响最终的数据处理结果。

(5)仪器安放点稳定性。在监测过程中, 地基InSAR系统应该保持水平, 不能受到扰动。所以在选址时要考虑仪器安放点的稳定性, 最好安放在稳定的基岩上。

(6)观测墩的修建。根据观测区域实际情况和精度要求, 考虑是否采用观测墩。

(7)仪器架设位置。观测边坡时, 应对仪器视线方向和主滑移方向的夹角进行权衡, 夹角越小, 对形变信号越敏感, 但不利于接收回波信号。

(8)滑坡植被。地基InSAR雷达波频率高, 波长短, 其信号对目标体的变形比较敏感。但如果在监测范围内植被较多, 由于植被在风的影响下会发生摆动, 雷达将不能很好甚至不能接收到地表变形信号, 这就会严重影响测量结果。因此, 监测区应选在植被稀疏的区域。

(9)反射器的安装。野外监测区域难免出现植被覆盖情况, 以及考虑到后期数据的地理编码等问题, 根据地形地貌条件在一些稳定的、接近滑坡的地方做一些参考点, 可以在这些参考点上安置角反射器。另外, 在滑坡上一些感兴趣点也可以安置角反射器。

IBIS-L地基InSAR形变时序分析流程主要包括原始数据聚焦、影像质量检查、特征点选取和形变分析等4部分^[3, 4]。

1.3.1 原始数据聚焦处理

天线波束内单个目标的回波信号会出现在线性轨道的每一个位置以及所有频率的原始数据中, 此时的图像处于散焦状态。对原始数据的每一个像元来说, 理想点目标具有相同的强度, 并且相位是发射频率和传感器在轨道上位置的函数, 强度和相位观测值需要转换到具有空间分辨率的格网上。图像聚焦分为距离向和方位向聚焦, 聚焦使得图像具有距离向和方位向空间分辨率。

1.3.2 影像质量检查

由于SAR系统采用侧视成像, 将会在SAR影像中产生叠掩、阴影和透视收缩等几何失真现象, 可利用地基InSAR初步观测一组数据, 通过强度图、热信噪比、估计信噪比、时间相干性、空间相干性和相位稳定度等查看该组影像的观测质量, 以判断地基InSAR观测位置、方位和仰角等是否满足要求, 及时调整观测角度以满足监测要求。

1.3.3 特征点选取

地基InSAR形变时序分析是基于相干点目标进行的, 一般通过设置估计信噪比、相干系数和相位稳定度满足一定阈值的像元, 筛选出稳定的相干点目标。由于IBIS-L系统采用Ku波段雷达(波长 λ$\begin{array}{c}\lambda \end{array}$为17.4 mm), 其对形变信息非常敏感, 初始形变时序分析得到的形变信息不但包含实际形变, 还包括温度、湿度、气压和粉尘等环境变化引起的相位变化。因此, 需要在稳定区域内选取一定数量的估计信噪比、相干系数和相位稳定度高的点作为地面控制点(ground control point, GCP), 通过对所有GCP在视线向的形变时序曲线加权平均, 建立环境相位时序校正函数, 消除环境因素的影响。

1.3.4 形变分析

地基InSAR时序分析主要采用简单网络和冗余网络组合方式对聚焦后的影像对进行干涉处理, 为提高处理效率, 工作中可采用简单网络组合方式对相邻时刻影像两两干涉处理。IBIS-L系统一般采用连续观测模式, 且观测距离较短, 空间基线为0, 因此无需进行图像配准和地形相位补偿。由于连续模式观测数据采集时间短, 目标的形变速率 v$\begin{array}{c}v\end{array}$在相邻的SAR影像获取时间间隔 Δt$\begin{array}{c}\mathit{\Delta t}\end{array}$内, 一般能够满足|v|<λ/4Δt$\begin{array}{c}\left|v\right|<\lambda /4\mathit{\Delta t}\end{array}$, 因此避免了相位解缠。

2014年9月1日06:00至2日05:00重庆奉节大树场镇发生了大面积山体滑坡, 冲毁了3栋楼房并阻断了奉竹公路。滑坡区位于梅溪河左岸斜坡地段, 滑坡区上部为治理后的危岩带, 如图2(a)所示。

	Figure Option
	图2 大树场镇滑坡Fig.2 Landslide in Dashuchang Town

该滑坡区为构造剥蚀侵蚀中山斜坡地貌, 总体上东高西低, 最高点和最低点分别位于滑坡区东侧山顶和梅溪河支流崔家河河谷。主要地貌类型包括陡崖、斜坡和河谷堆积地貌, 陡崖分布在滑坡区斜坡中上部, 南高北低, 陡崖带长约440 m, 高度为30~65 m, 近乎直立, 2011年完成危岩带治理工程, 处于稳定状态; 斜坡为滑坡区的主要地貌类型, 分布在陡崖带的上部和下部, 上部陡斜坡总坡高约80 m, 平均坡角约43° , 下部斜坡总坡高约200 m, 呈折线形, 其上部较陡, 坡角约55° , 下部稍缓, 平均坡角约32° , 斜坡脚部为生活区; 河谷堆积地貌分布于斜坡坡脚外, 为梅溪河支流崔家河河谷堆积地貌区, 地势比较平坦, 该区域位于在建的大树场镇新区, 已修建了多栋住宅楼。

本次观测将地基InSAR系统架设于滑坡体正前方一个宽敞、稳定的场地上, 设备前方视野开阔, 无任何物体遮挡, 扫描区域能够完全覆盖治理后的危岩带和滑坡区, 2个观测区域植被覆盖稀少, 满足了监测的基本条件(图2(b))。观测时间自2014年9月23日17:06至26日8:53, 考虑到高角度斜坡观测, 天线倾角设置为15° , 有利于回波信号接收, 地基InSAR系统其他参数设置如表1第1类所示, 共获得654景SAR影像。观测期间, 高温、阴雨和强降雨等恶劣天气对IBIS-L地基InSAR系统几乎没有任何影响^[5]。

表1 IBIS-L地基InSAR系统观测参数配置Tab.1 Parameters configuration of IBIS-L GB-InSAR

本次观测数据处理选取了位于治理后危岩带上的4个GCP, 由于观测时间相对较短, 滑坡体上部治理后的危岩带在观测期间可认为基本保持稳定; 当地相关地质部门全球定位系统长期观测结果也表明危岩带基本稳定。IBIS-L系统测量精度能控制在[-0.5 mm, +0.5 mm]范围内。该观测区经过时空滤波的累积形变如图3所示。

	Figure Option
	图3 大树场镇滑坡地基InSAR观测累积形变Fig.3 GB-InSAR cumulative displacement of landslide in Dashuchang Town

图3表明, 观测期间治理后的危岩带基本处于稳定状态; 但滑坡区上存在2处显著活动区, 分别位于滑坡区左侧中上部(图2(b)区域Ⅰ 和图3区域Ⅰ )和右侧中上部(图2(b)区域Ⅱ 和图3区域Ⅱ ), 这2处形变区域主要是由堆积松散土的滑动引起的, 滑移区域Ⅰ 和Ⅱ 最大变形分别达到120 mm和75 mm。该观测区典型像元时序形变分析如图4所示。

	Figure Option
	图4 大树场镇滑坡典型像元时序形变分析Fig.4 Time-series displacement analysis of typical pixels of landslide in Dashuchang Town

从图4可以看出, 危岩带相干点P30[-37.2 m, 443.8 m]形变时序分析表明, 治理后的危岩带比较稳固。滑坡区2处显著活动区的形变时序特征大体一致: 9月23日17:00至24日10:00, 当地天气多云转阴, 空气湿度为70%~75% rH, 相干点都以较小的下滑速度滑移, Ⅰ 处(如P59[-6.2 m, 351.7 m])下滑速率最大约为1.5 mm/h; 24日10:00左右, 当地开始降小雨, 空气湿度陡升至90% rH左右, 区域Ⅰ 和Ⅱ 表面整体发生了陡滑现象; 14:00降雨突然变大并持续2 h, 滑移区域下滑速率加快; P59处至观测结束仍处于快速下滑状态, 而Ⅱ 处(如P62[46.5 m, 275.8 m]和P65[72.2 m, 304.3 m])在25日12:00基本处于稳定。根据区域Ⅰ 和Ⅱ 特征点的形变时序分析, 降雨冲刷松散土表层是造成2处活动区加速下滑的主要原因。

人工不稳定斜坡位于神农架林区某度假酒店建设场地的西北部, 狭长槽谷谷坡的下部。该区总体形态呈“ 舌” 型, 自然地形呈南高北低, 地形坡度为17° ~20° , 局部陡坎为40° ~60° 。坡面受人类活动的影响, 植被较少, 地形整体条件不利于坡体的稳定。

根据现场勘查, 后缘边界位于缓坡平台, 后缘高程为1 720 m; 前缘至神农架旅游公路以南, 因工程建设施工场地平整切坡, 其前缘为人工切坡坡脚, 前缘高程为1 640 m, 相对高差为80 m; 东侧以山脊的坡脚处为边界, 有基岩出露; 西侧以小型冲沟为边界, 出露神农架群(P_tsh)白云岩。不稳定斜坡纵长为230 m, 宽度为150~200 m, 不稳定斜坡上松散堆积体约46.1× 1 04$\begin{array}{c}{0^4}^{}\end{array}$m³。该区域堆积的第四系残坡积( Qel+dl4$\begin{array}{c}{Q}_4^{\mathit{el}+\mathit{dl}}\end{array}$)碎、砾石土厚度较大, 其危险性大, 危害亦大, 现场情况如图5(a)所示。

	Figure Option
	图5 神农架林区不稳定斜坡Fig.5 Unstable slope in Shennongjia forest region

在前期的勘查中, 拟将地基InSAR系统安置在神农架旅游公路北侧的冲沟内, 其沟口正对着不稳定斜坡, 扫描区域能够完全覆盖, 沟口右侧地表部分基岩出露, 较为稳固。但该冲沟坡降大, 且沟内残存有建设神农架飞机场的大量弃渣, 在暴雨工作情况下多次发生过泥石流。考虑到观测期间持续降雨(中到大雨), 存在爆发泥石流的可能性, 因此地基InSAR系统架设在公路北侧路边, 这种近距离的观测造成了观测区域只覆盖了不稳定斜坡部分区域^[6]。

本次观测时间自2014年9月27日11:48至19:06, 共获取74景影像, 观测参数如表1第2类所示。由于观测距离较短, 为了获得较大的视场, 天线倾角设置为20° , 可令视线方向与形变方向的夹角较小, 有利于形变信号的提取。近距离的观测保证了空间分辨率和观测精度, 但较大的天线倾角不仅减弱了雷达回波信号, 同时造成了影像遮掩现象, 使得不稳定斜坡的中部和上部平台的底部无法观测到。

扫描到的区域大部分点回波信号较好, 视场中部0+165桩附近出露基岩处选取了2个基准GCP, 用以校正形变时序结果。不稳定斜坡底部大部分位置为基岩, 并进行了防滑处理, 整体较为稳定。出露基岩上的相干点(如P20[-30.0 m, 92.9 m])的时序分析结果表明本次观测的测量精度基本能控制在[-0.5 mm, 0.5 mm]内。地基InSAR观测不稳定斜坡的累积形变如图6所示。

	Figure Option
	图6 神农架林区不稳定斜坡地基InSAR观测累积形变Fig.6 GB-InSAR cumulative displacement of unstable slope in Shennongjia forest region

由图6可知, 不稳定斜坡中上部有整体靠近仪器的趋势, 即整体下滑; 斜坡平台自下往上累积形变量逐渐增大; 较大形变出现在不稳定斜坡的第一级平台的流水沟附近(图5(b)区域Ⅰ 和图6区域Ⅰ )以及第三级平台中部的土堆处(图5(b)区域Ⅱ 和图6区域Ⅱ ), 最大累积形变达到25 mm, 这主要是暴雨积水导致的第三级平台中部堆积土层塌陷以及第一级平台表层覆土滑动。该区典型像元时序形变分析如图7所示。

	Figure Option
	图7 神农架林区不稳定斜坡典型像元时序形变分析Fig.7 Time-series displacement analysis of typical pixels of unstable slope in Shennongjia forest region

不稳定斜坡第二、三级平台(如P36[3.6 m, 116.6 m], P1[-20.0 m, 145.9 m])在观测期间前4 h内下滑速度较大, 随后下滑速度放缓, 这与降雨情况非常相关, 观测区前4 h内一直持续暴雨, 之后转为小雨。第一级平台流水沟处(如P34[10.9 m, 72.3 m])在前3 h下滑速度较小, 随着持续降雨造成土体松动, 随后呈现出加速下滑的趋势。

备战铁矿矿区位于天山主峰―博罗霍洛山北坡的中高山区, 山体走向近EW, 地势南高北低, 高程为3 862~3 416 m, 地形坡度为25° ~35° 之间, 沟深坡陡, 属高山深切地貌。矿体位于山谷地形, 东、南、西三面环山, 仅有北面为狭长山谷。目前采用“ 露天开采+公路开拓汽车运输” 方案, 自上而下水平分层台阶式采矿方法。全矿圈定了一个露天采矿场, 采场共设17个最终台段, 最终台段高度为24 m, 最终边坡角为50° ~51° , 设计露天开采标高为3 704~3 284 m。矿区属次不稳定区, 现状地质灾害不发育, 露天开采已形成了大采坑, 对地质环境有一定破坏和影响, 矿区地质环境质量中等, 现场情况如图8(a)所示。

	Figure Option
	图8 备战铁矿3号矿体露天采矿边坡Fig.8 Open-pit slope in No.3 ore body of Beizhan iron mine

2015年6月29日21:06至7月2日14:37, 对备战露天铁矿3号矿体开采面进行观测, 地基InSAR系统安置在矿体开采面对面的观景台旁边, 设备所处位置稳定, 观测视野开阔, 无障碍物遮挡, 适合架设观测设备。本次虽为高边坡观测, 由于设备所处高程与开采面中部相当, 考虑到能尽可能观测到开采面全貌, 故采用增高垫抬高仪器, 并设置天线倾角为5° , 雷达系统其他参数设置如表1第3类所示, 共获得雷达影像650景。观测期间, 系统经历了大风扬尘、放炮作业和雨雪冰雹等恶劣条件影响, 但仍圆满完成了边坡持续观测^[7]。

由于矿体所处山谷地形气象条件复杂多变, 观测视场内不断有矿车经过扬起漫天尘土, 因此数据处理过程中环境因素校正工作十分棘手。本文提出了在稳定区域均匀选取GCP的方法, 根据这些GCP利用加权最小二乘法估计校正系数, 准确提取了观测区域的形变信息(图9)。

	Figure Option
	图9 备战铁矿露天采矿边坡地基InSAR观测累积形变Fig.9 GB-InSAR cumulative displacement of open-pit slope in Beizhan iron mine

坚硬铁矿岩石上的相干点(如P3[72.8 m, 634.2 m])的时序分析结果表明, 本次观测的测量精度能控制在[-0.5 mm, 0.5 mm]内。地基InSAR累积形变结果表明, 虽然开采现场环境复杂(多变的天气、漫天粉尘、矿车和铲车的扰动等), 但开采面整体保持稳定; 形变较大的区域位于开采面底部(图8(b)区域Ⅰ 和图9区域Ⅰ ), 最大形变量超过15 mm, 主要是由爆炸点处的堆积铁矿石滑动以及铲车移动铁矿石等因素造成的; 开采面上方的“ V” 型松散土层凹谷(图8(b)区域Ⅱ 和图9区域Ⅱ )表层也有稍微的滑动, 最大滑动速率约2 mm。该区典型像元时序形变分析如图10所示。

	Figure Option
	图10 备战铁矿露天采矿边坡典型像元时序形变分析Fig.10 Time-series displacement analysis of typical pixels of open-pit slope in Beizhan iron mine

相干目标时序分析结果显示, 开采面上方的“ V” 型松散土层凹谷表层有轻微持续下滑趋势(如P47[-75.7 m, 612.1 m]), 在雨季需要进行关注。形变最大区域位置上方的相干点(如P5[0.0 m, 568.8 m])则呈现出较大的持续下滑趋势, 主要是在铲车和矿车等人为扰动下, 松散铁矿石层因自身重力下滑导致的。开采面底部形变最大区域位置(如P32[-7.1 m, 538.7 m])观测初期处于较为稳定的状态, 由于铲车移动铁矿石, 于7月2日00:00开始加速下滑, 最大下滑位移超过30 mm。

备战铁矿矿区东北部数百米处即为天山山脊, 山脊线多为尖棱状, 常年冰川覆盖, 具有典型冰川地貌特征, 现代冰川、冰蚀洼地、冰川“ U” 型谷和冰蚀崖较发育。选矿厂及生活区与矿区通过一条简易公路连接, 分别位于天山主脉南北, 为穿越天山主脉在海拔3 800 m处修建了一条隧道。本次冰川运动观测试验旨在研究地基InSAR探测冰雪边坡与冰川活动的可行性。所选的冰川位于隧道南口旁边(图11(a)), 由于近几年的铁矿开采、大量矿车扬起的路面粉尘以及全球气温升高等因素影响, 导致冰舌退缩加剧, 冰川面积缩小。

	Figure Option
	图11 备战隧道南口冰川Fig.11 Glacier in the south of Beizhan tunnel

地基InSAR系统安置在靠近冰舌末端的简易公路边平台处(图11(b)), 观测平台稳定, 视野开阔, 满足地基InSAR观测要求。数据获取时间跨度为7月3日12: 46— 17: 29, 共采集43景SAR影像, 具体参数如表1第4类所示。冰川表层滑动情况复杂, 确定稳定的GCP具有一定的难度, 同时由于冰舌属于低倾角边坡, 本次地基InSAR观测试验选取雷达视角5° 扫描, 这种情况下两侧山体雷达信号强度高, 便于在两侧山体选取GCP, 冰舌表碛以及正对的冰雪覆盖层亦具有较强的回波信号^[8]。

数据处理所选的GCP位于两侧山体的基岩上, 观测时间仅4.7 h, 认为GCP没有发生移动。稳定区域观测点(如P32[-156.5 m, 1 009.3 m])地基InSAR时序分析结果表明, 隧道南口冰川运动的观测精度控制在[-0.5 mm, 0.5 mm]内, 能够满足监测需要。该区累积形变及典型像元时序形变分析分别如图12和图13所示。

	Figure Option
	图12 备战隧道南口冰川地基InSAR观测累积形变Fig.12 GB-InSAR cumulative displacement of glacier in the south of Beizhan tunnel

	Figure Option
	图13 备战隧道南口冰川典型像元时序形变分析Fig.13 Time-series displacement analysis of typical pixels of glacier in the south of Beizhan tunnel

本次观测的前一晚该地区普降大雪, 观测期间冰雪融化使得冰舌表碛产生滑动, 远处的冰川积雪也产生了一定的滑动, 远处冰川和冰舌上段都发生了朝向雷达方向的移动(如P26[-308 m, 582.9 m]和P30[-232.1 m, 668.6 m]), 最大滑动位移超过15 mm。左侧山体底部松散土层个别区域发生了下滑(如P18[-172.6 m, 881.1 m]), 相对于雷达传感器来说, 发生了远离雷达传感器的运动, 在地基InSAR累积形变图上呈现正值。本次观测时间较短, IBIS-L地基InSAR系统采用的波长较短, 观测结果反映的一般是冰川表碛和积雪等的滑动情况。为了获得冰川的整体活动情况, 以后的工作中需要采取在冰川上架设角反射器、选取与冰川整体运动情况较为一致的天然散射体等措施, 采用波长较长(如X和C波段)的地基InSAR更利于冰川整体运动观测。