InSAR在滑坡监测中的应用

1、InSAR在滑坡监测中的应用一、概述合成孔径雷达（SAR）最早发展于20世纪50年代，作为有源系统，SAR具 有全天候、全天时的工作能力。可在不同的微波频段、不同极化状态下得到地面 目标的高分辨率图像为人们提供各种有用的信息，广泛应用于地形测绘、地质研 究、防灾减灾等诸多方面。20世纪60年代末合成孔径雷达干涉测量技术（INSAR） 发展起来，它是利用合成孔径雷达的相位信息提取地表的三维信息和高程变化信 息的一项技术。它的出现大大扩展了 SAR的应用领域，能够获取高精度的地形信 息，同时还可以检测地球陆地表面和冰雪表面的微小变化，监测时间从几天到几 年，可获得全球高精度的（毫米级）高可靠的地表

2、变化信息。近年来地震、火山、滑坡和地面沉降等地质灾害越来越严重地威胁着人类的 生存空间，针对这种灾害而发展起来的地表形变监测和测量技术就显得尤为重 要。合成孔径雷达干涉测量可以大面积的采样、测量时间短，同时成本也比较低， 而且合成孔径雷达干涉测量，具有全天侯、全天时、高分辨率和连续空间覆盖的 特征,能够提供短周期内空间连续曲面形变信息，可以满足沉陷监测要求,弥补地 面常规测量离散点的不足，特别是能够提供高水平和高垂直分辨率的三维数字模 型。并且雷达数据下载快捷，时间延误少，加之越来越成熟的配套处理软件，使得 地面沉降数据的提取十分迅速，可接近准实时动态监测，是现有滑坡监测手段的 有益补充，用于

3、地表及山体滑坡监测具有很大的优势和很好的前景。二、INSAR干涉基本原理从历史上来看，现代合成孔径雷达干涉测量技术（InSAR技术）的发展源于 Thomas Young于1801年所完成的“杨氏双狭缝光干涉实验”。从两条狭缝射出 的波束具有不同的传播距离，也就是波的相位不一致，这样引起波的叠加或削减， 从而在白板上可观察到明暗相间的条纹，也就是所谓的“干涉条纹”InSAR正 是受这一实验启发而发展起来。合成孔径雷达干涉测量技术是以合成孔径雷达复数据提取的相位信息为信息源 获取地表的三维信息和变化信息的一项技术。该方法基于时间测距的成像机理， 充分利用了雷达回波信号所携带的相位信息来提取地表高程

4、信息InSAR技术原 理是通过两副天线同时观测（单轨双天线模式）或两次近平行的观测（重复轨道 模式），获取地面同一景观的复数图像对（重复观测数据）。由于目标与两天线位 置的几何关系，在复图像上产生了相位差（同一目标对应的两个回波信号之间存 在相位差），形成干涉纹图（interferogram）。干涉条纹包含了斜距向上的点与 两天线位置之差的精确信息。因此，利用传感器高度、雷达波长、波束视向及天 线基线距之间的几何关系，可以精确地测量出图像上每一点的三维位置和变化信 息。图1为InSAR技术的基本原理图所示图1 InSAR技术的基本原理假设飞行平台上同时架设了两部天线，S2，若由S1发射信号，和

5、S2同 时接收从目标T返回的信号，由式u = e j = URe -（e押）+巾 （可知，S1或S2接收到的从目标丁返回的信号之相位可以表示为甲=七（P t + P，）=甲 t + 甲 r（2） 式中，下标t，r分别表示发射和接收信号的相关参数；入为波长；p为从天线到 目标的斜距。两幅天线所接收到的信号之相位差为（3）式中，P为系数，若是一幅天线用于发射信号，在干涉图上只反映出单程（信号 的返程）的相位差（单轨双天线模式）；若是两幅天线都发射和接收信号，则P=2, 即反映出往返双程的相位差（单天线重复轨道模式）。u当两幅天线所形成的复数影像精确配准后，对应像素值共轭相乘就可以在每个像 素上得到

6、如式（4）所示的相位差，形成通常我们所说的干涉图或干涉条纹图（Interferogram），艮口复数影像1：u = u |e j 91复数影像2：u = ue j 9 2干涉图：_ _ * u = u - uInt12Int一 9 2）相位差为中=arctan、(u 或)(5)Rm (/ )e Int从（5）式中可以看出，由复数影像七川2求得的相位差，实际上只是0, 2n 区间的主值，即对2n取模的值甲w。确切地说式（5）应写成中=w B，可 以得出0 =以一 arcsin抑2兀PB（9）(10 a)(10 b)h = H - p. - cos0y = p - sin0式（9）和（10）揭示了

7、干涉相位差6和高程h之间的关系。也就是说，如果天 线位置（参数H、B、a ）和雷达成像系统参数（。）等，就可以从6计算出地 表的高程值h。三、D-INSAR干涉测量原理INSAR技术进一步发展为差分合成孔径雷达干涉技术（differential synthetic aperture radar interferometry， D-INSAR），它可以用来监测地表 目标的水平和垂直运动，目标相对位移量的测量精度可达到厘米级甚至毫米级。 常用的D-INSAR数据处理方法有三轨道法、四轨道法和已知DEM的双轨道法。如前面提到，实际上两次重复观测（或成像）在技术上不可能达到完全重复 的轨道和其他参数。

8、通过比较两幅SAR影像的相位变化（干涉图）还不足以求出 地表的形变，但是可以有干涉图导出视差的变化，从而反推出地形信息。本文以 三轨道法为例，其成像几何示意图如图2所示。P3图2 D-INSAR的成像几何示意图设天线和L分别接收到的信号的传播路径P 1，P 2的差为第，这里设定由S1s 一 获得的影像为主影像，由2获得的影像为从影像。由几何关系可得21p2 = p2 + B2 2p - B - cos(9 0 )2PB - sin(0 一以)+P1 + P 2B2兰一B - sin(0以)+ 2 P18p = p - p =B 2P1 + P2(11)由于B p 1，通常为了简单起见果的影响在

9、1cm之内。这样就可以得到8p 兰-B - sin（ 0 - a ） = B H让等式右边的第二项近似为零，此时对最终结（12）式中B|l为基线B在视线方向的分量，即B | 兰 B - sin( 0 a )（13）于是得到综合式（12）和式（13）以及式(14)表明：用干涉方式所量测得到的相位差（或干涉相位）与视线方向基线分量成正 比。在重复轨道干涉成像模式下，式（14）中的P=2即有2 n（15）9=9 1 f 2 = T P( P 2考虑在同一区域先后形成的两幅干涉图的情形。如果在同一地区再一次成像得到 第三幅SLC影像。如上图所示，传感器的位置为S3,信号的传播路径为P3。与 第一对干涉

10、影响中的主影像形成第二幅干涉条纹图，那么两幅干涉条纹图的相位 就可以相互比较。形成第二幅干涉图的基线距B（ Si和S3之间的基线），基线与水平方向的夹角为a，相位差为9，0不变。由式（15）可得(16)由式（15）和式（16）可以得虹=9 B II（17）即相位的比值等于基线距水平分量的比值，与地形本身无关。注意这里的相位是 指未缠绕的实际相位值。如果在第三次成像之前，地表发生形变，并假设这种形变与雷达分辨率单元 相比很小，可以认为雷达信号仍是相关的。那么这时相位信息除了与地形信息有关，还包含了雷达视线方向的形变量AP。第二幅干涉图的相位可以表示为9= 丁 (B +Ap) 人(18)同理考虑到

11、式（15）也可以将视线方向的形变量Ap的表达式写为:AP9- 9(19)根据式（13），基线的比值还可以表示为B sin(。- a )(20)B nsin(。一 a ) 由式(19)和式(20)就可算出地面形变。D-InSAR技术可以达到毫米级别，相对InSAR技术来说具有更高的观测精度。 但是D-InSAR技术实施过程中，受到时间、空间去相关和大气效应等因素的影响， 存在测量范围减小和精度减弱等问题。四、PS INSAR技术PS ( Permanent Scatters,又称永久散射体)技术是20世纪末由意大利学者A. Fer2retti，C. Prati和F. Rocca首先提出的，以解决

12、常规干涉中大气影响、 失相干、DEM误差等问题，极大地拓展了 InSAR技术的应用前景，导致对地观 测技术发生了又一次革命性的飞跃，为精确研究地壳形变提供了强有力工具。所谓永久散射体，即是指在相当长的时间内仍然保持稳定反射特性的散射 体。PS技术是D-InSAR技术不断创新的结果，其主要目标就是研究某一区域内 地表长期缓慢变形的过程及变化规律。我们知道，InSAR技术通过比较相位差测 量地形和形变，但形变测量值对大气造成的相位延迟也相当敏感,大气误差极大 影响了 InSAR测量的精度。同时，当两次成像时间间隔比较长时，或是成像季节 差异比较大时，地表的反射特性可能会发生很大的变化，这势必会造成

13、图像的严 重时间失相干，导致干涉处理无法进行。而PS技术作为对D-InSAR技术的创新， 在解决时间失相关的同时，计算并消除大气影响，保证干涉处理的正常进行。PS InSAR的基本原理是：在给定的一组雷达干涉图像中，按照某种准则选 择相位稳定的一系列点作为永久散射体，也就是PS点，并且根据给定的相位模 型，去除这些PS点的大气附加相位，DEM误差和其他噪声，进而得到准确的地 表升降变化值。通常，PS技术需要处理同一地区30景以上的干涉图，以选取一 组反射特性强、时间上稳定的地物目标作为PS点，保证计算结果的精确性。最早利用PS InSAR技术进行的研究有，法国学者对巴黎市区因地铁建设造 成的沉

14、降研究，美国洛杉矶盆地的沉降研究等，研究甚至观测到博物馆建筑上金 属结构的季节性热胀冷缩现象。F. Rocca于2002年将PS InSAR技术用于监测 城市大面积长时间小幅度沉降(可监测到的年沉降速度为几毫米)。德国GFZ的 夏耶最早将人工角反射器应用于研究三峡滑坡监测，虽然在三峡地区由于各种 原因一时还无法得到最终结果，但该方法已显示出强大的生命力，国内许多应用 研究部门对此表现出浓厚兴趣。由中国科技部、欧空局等单位合作的“龙计划”， 在三峡库区安装了角反射器，用于监测三峡地区泥石流、滑坡等地质灾害；中国 地震局地壳应力研究所张景发、英国伦敦大学学院Peter等研究人员在西藏当雄 活动断裂

15、带区域安装了角反射器，用于监测地壳运动形变。与常规InSAR技术相比，PS技术具有以下几个优势：长期监测。精确基线计算。大气相位估计。单个散射体高度和运动量计算。PS INSAR技术的特点是可以长时间对一个地区进行地面沉降的监测，并且测量 的精度非常高，缺点是数据处理特别复杂。该技术主要用于城市的地面沉降，机 场等大建筑物的沉降监测、采矿区的开采沉陷监测。五、INSAR数据与处理流程图3是从接收雷达信号到产生数字高程模型的典型处理流程图。干涉数据的处理一般包括如下工作：(1)干涉SAR数据的特殊要求基线的长度要满足相干的要求，基线的长度超过一定的值则两幅图像不相干。 要获得好的干涉效果，两幅图

16、像间的垂直基线(基线是两遥感平台间的矢量距离) 最好在100 m到300 m之间；而监测地表变化，垂直基线最好在50 m到100 m 之间；检测如海浪运动等的微小变化，最好小于5 m21。相干图像获取期间成像区域下垫面变化要小，这样才能使图像相干并满足精确 几何配准时对光谱相似性的要求。处理成SLC(single look complex)格式数据。这要求专门的处理软件来完成。图3 SAR干涉测量流程六、INSAR在滑坡监测应用中的优势与不足通过各个国家的学者开展的实证研究表明，INSAR技术在监测滑坡具有如 下其它方法所无可比拟的优势：(1)全天候、大范围卫星雷达干涉测量一次能覆盖几百至上千

17、平方公里范围，可获得整个覆盖 范围内与此相应的形变位移数据。另外，由于卫星雷达成像能穿透云层且没有昼 夜之分，雷达数据下载快捷时间延误少，加之越来越成熟的配套处理软件，使得 数据提取十分迅速，可以接近准实时动态监测。(2)高分辨率、高精度雷达差分干涉测量技术能以厘米量级甚至更小尺度测量地表形变，它不仅 能够提供宏观的静态信息，而且能够给出定量的动态信息。另外，INSAR技术了提供时间跨度较大的SAR影像数据，例如ERS可提供 从1991年开始的十年间的数据。而且利用这项技术可以花费与其他方法相比较 少的费用来获取观测数据，这使得长期对滑坡进行监测成为可能。因此，总的来 说利用INSAR技术能得到十分详尽的地形高程图，而且取样率远远大于最详尽的 调查。滑坡常常发生在地形比较陡峭的地区，因此进行滑坡监测必须要有精确地 DEM数据，但是在数据的采集过程中，会有很多因素带来误差，影响干涉图的质 量。而且滑坡常常是在很短的时间内造成地表的下降，从而能导致干涉图像对失 相干，产生去相关问题。另外，在雷达数据采集的时期内，要想获得具有好的相 关性的干涉图有赖于季节和天气条件。通常，潮湿的天气和高植被