第2期 王桂杰等：差分干涉合成孔径雷达技术在广域滑坡动态辨识上的实验研究 ·137° 0.0667m月-1（位移值大于0.2m)的四个区域来 同时进行GPS监测的山R6号滑坡为例，将其 显示.在水平滑动位移形变状态图中，按水平位移 D-ISA监测结果进行详细分析.图8为山R-6号 与垂直位移所成比例的倍数，将抬升区域由水平位 滑坡上D-SAR监测的详细滑动形变图.在垂直滑 移绝对值大小按红色由浅到深划分为0~0.774m 动形变位移图中按滑动速率小于0.0167m月-1 和大于0.774两个区域.将沉降区的水平位移按 (位移值小于0.05m,滑动速率为0.0167~ 移动速率绝对值及位移绝对值大小也进行相同的划 0.0334m月-1（位移值为0.05~01m)和滑动速 分，按绿色由浅到深划分为滑动速率为0~0.0647 率大于0.0334m月1（位移值大于0.1m从慢到 m月-1（位移值为0~0194m)、滑动速率为 快分为三个区域，在水平形变位移图中按绝对值滑 0.0647-0.129m月（位移值为0.194~0.387 动位移速率小于0.0647m月-1（位移值小于0194 m)滑动速率为0.129-0.258m月-（位移值为 m)、滑动速率为0.0647~0.129m月-1（位移值为 0.387-0774m以及滑动速率大于0.258m月-1 0.194~0.387m和滑动速率大于0.129m月（位 (位移值为大于0.774m的四个区域.其中水平位 移值大于0.387m)也从慢到快分为三个区域分别 移是指南北向的水平位移，图中向南符号为负，向北 用黄、蓝和红三种颜色来显示.从形变图中能清晰 符号为正.形变图中沉降区域的垂直和水平位移分 地分辨出不同区域的滑动速率和位移情况.图中将 级图之所以按着一定的所成比例来划分，是因为得 山R-6号滑坡划分为四个区域2-，其中I区的部 到的研究区域的滑动位移是按着沿卫星视线向的形 分区域处在垂直向滑动速率小于0.0167m月- 变，根据卫星侧视倾角和轨道参数分解而得到的，所 水平向小于00647m月一)的滑动区间上，Ⅱ区 以相同区域的水平位移值与垂直位移值具有基本相 的部分区间和区处在垂直向滑动速率大于 同的比例系数，因而图中水平形变分级图和垂直形 0.0334m月-1（冰平向大于0.129m月-）的滑 变分级图的变化区域是基本相一致的，水平形变分 动区间上，其余区域都处在垂直向滑动速率为 级图仅反映了其水平滑动的速率和位移大小趋势. 0.0164-00334m·月-1（水平向0.0647~ 从研究区域的滑动形变位移图可知，其沉降区 0.129m月-1)的滑动区间上.图中红色线体为长 域最大滑动位移值为25.3四由于所用A1OS卫星 江三峡院现场监测划定的重点灾害区，可知处在 PALSAR传感器的波长为I波段23四其形变监测 D-InSR监测结果滑动位移值为5~13.3m的区 灵敏度为12波长即为12.5m(沿卫星视线向 间上. 的，再考虑到各种可能误差（后面分析，其区域内 a 观测时段92内的变动位移很小.可看到其红色抬 升区域的抬升位移值绝大多数都在小于20m的范 00.018-0050 ■0.050-0100 围内，因此抬升区域的抬升位移值大部分都处在误 m0.100-0.133 差影响范围内.因此对研究区域垂直沉降滑动位移 速率大于667m月-(0.0667m月-1位移值大 于20四水平滑动位移速率25.8m月-1(0.258 m月1)位移值大于77.4的区域划定为可能发 b 生滑坡和滑坡活动的风险区域，即为图中用紫色线 ■-0.5140387 圈圈出的两个区域.经与研究区域地质资料及光学 ■-0387-0194 ■-0.194-0.070 影像确认，此两区域一个区域位于金沙江右岸一个 较大泥石流沟的沟口附近，另一区域位于该泥石流 沟沟中汇水区的中部.其形变位移值较大的原因除 可能因为崩滑原因外，也可能是由于沟中物质迁移 图8L1R-6号滑坡的滑动形变图（单位：m,()92d的垂直 所引起的. 滑动形变图：(b)92的水平滑动形变图 F琴8 Depmation map of山R-6and时de(ni诉m:(a)92d 3与GPS对比研究 venical defmation map (b)92 d horgontal depmation map 3.1与GP结果的对比 将DSR获得的山R-6号滑坡的形变位移 为证明DSAR技术在滑坡辨识监测上的可 值与己获得的地面GPS监测结果相比较，图8中绿 行性及正确性，本文以研究区域内正处活动状态并 色字体编号的紫色点为现场GPS监测点，其中点名第 2期 王桂杰等:差分干涉合成孔径雷达技术在广域滑坡动态辨识上的实验研究 0.066 7 m·月 -1 (位移值大于 0.2 m)的四个区域来 显示.在水平滑动位移形变状态图中 , 按水平位移 与垂直位移所成比例的倍数 ,将抬升区域由水平位 移绝对值大小按红色由浅到深划分为 0 ～ 0.774 m 和大于 0.774 m两个区域.将沉降区的水平位移按 移动速率绝对值及位移绝对值大小也进行相同的划 分 ,按绿色由浅到深划分为滑动速率为 0 ～ 0.064 7 m·月 -1 (位移值为 0 ～ 0.194 m)、滑动 速率为 0.064 7 ～ 0.129 m·月 -1 (位移值为 0.194 ～ 0.387 m)、滑动速率为 0.129 ～ 0.258 m·月 -1 (位移值为 0.387 ～ 0.774 m)以及滑动速率大于 0.258 m·月 -1 (位移值为大于 0.774 m)的四个区域 .其中水平位 移是指南北向的水平位移 ,图中向南符号为负 ,向北 符号为正.形变图中沉降区域的垂直和水平位移分 级图之所以按着一定的所成比例来划分, 是因为得 到的研究区域的滑动位移是按着沿卫星视线向的形 变 ,根据卫星侧视倾角和轨道参数分解而得到的,所 以相同区域的水平位移值与垂直位移值具有基本相 同的比例系数,因而图中水平形变分级图和垂直形 变分级图的变化区域是基本相一致的, 水平形变分 级图仅反映了其水平滑动的速率和位移大小趋势. 从研究区域的滑动形变位移图可知, 其沉降区 域最大滑动位移值为 25.3cm,由于所用 ALOS卫星 PALSAR传感器的波长为 L波段 23 cm, 其形变监测 灵敏度为 1 /2 波长即为 12.5 cm(沿卫星视线向 的 ),再考虑到各种可能误差 (后面分析),其区域内 观测时段 92 d内的变动位移很小 .可看到其红色抬 升区域的抬升位移值绝大多数都在小于 20 cm的范 围内, 因此抬升区域的抬升位移值大部分都处在误 差影响范围内.因此对研究区域垂直沉降滑动位移 速率大于 6.67 cm·月 -1 (0.066 7 m·月 -1 )位移值大 于 20 cm、水平滑动位移速率 25.8 cm·月 -1 (0.258 m·月 -1 )位移值大于 77.4 cm的区域划定为可能发 生滑坡和滑坡活动的风险区域 ,即为图中用紫色线 圈圈出的两个区域.经与研究区域地质资料及光学 影像确认,此两区域一个区域位于金沙江右岸一个 较大泥石流沟的沟口附近 , 另一区域位于该泥石流 沟沟中汇水区的中部 .其形变位移值较大的原因除 可能因为崩滑原因外, 也可能是由于沟中物质迁移 所引起的. 3 与 GPS对比研究 3.1 与 GPS结果的对比 为证明 D-InSAR技术在滑坡辨识监测上的可 行性及正确性,本文以研究区域内正处活动状态并 同时进行 GPS监测的 L1R-6 号滑坡为例 , 将其 D-InSAR监测结果进行详细分析.图 8为 L1R--6号 滑坡上 D-InSAR监测的详细滑动形变图.在垂直滑 动形变位移图中按滑动速率小于 0.016 7 m·月 -1 (位 移 值 小 于 0.05 m)、滑 动 速 率 为 0.016 7 ～ 0.033 4 m·月 -1 (位移值为 0.05 ～ 0.1 m)和滑动速 率大于 0.033 4 m·月 -1 (位移值大于 0.1 m)从慢到 快分为三个区域, 在水平形变位移图中按绝对值滑 动位移速率小于 0.064 7 m·月 -1 (位移值小于 0.194 m)、滑动速率为 0.064 7 ～ 0.129 m·月 -1 (位移值为 0.194 ～ 0.387m)和滑动速率大于 0.129m·月 -1 (位 移值大于 0.387 m)也从慢到快分为三个区域, 分别 用黄 、蓝和红三种颜色来显示 .从形变图中能清晰 地分辨出不同区域的滑动速率和位移情况 .图中将 L1R--6号滑坡划分为四个区域 [ 12--13] , 其中 Ⅰ区的部 分区域处在垂直向滑动速率小于 0.016 7 m·月 -1 (水平向小于 0.064 7 m·月 -1 )的滑动区间上, Ⅱ区 的部分区 间和 Ⅲ区 处在垂直 向滑动速率 大于 0.033 4 m·月 -1 (水平向大于 0.129 m·月 -1 )的滑 动区间上, 其余区域都处在垂直向滑动速率为 0.016 4 ～ 0.033 4 m·月 -1 (水 平 向 0.064 7 ～ 0.129 m·月 -1 )的滑动区间上 .图中红色线体为长 江三峡院现场监测划定的重点灾害区, 可知处在 D-InSAR监测结果滑动位移值为 5 ～ 13.3 cm的区 间上 . 图 8 L1R--6号滑坡的滑动形变图(单位:m).(a)92d的垂直 滑动形变图;(b)92d的水平滑动形变图 Fig.8 DeformationmapofL1R-6 landslide(unite:m):(a)92d verticaldeformationmap;(b)92dhorizontaldeformationmap 将 D-InSAR获得的 L1R--6号滑坡的形变位移 值与已获得的地面 GPS监测结果相比较, 图 8中绿 色字体编号的紫色点为现场 GPS监测点, 其中点名 · 137·