基于地基雷达干涉测量技术的大坝边坡形变监测及应用

李翔宇 雷添杰 陈文晋

摘要 为了解决传统大坝边坡形变监测方法仅能获取点变形信息，且监测周期长、监测精度低、无法进行实时监测等问题，以云南牛栏江红石岩边坡为例，提出了基于地基雷达干涉测量技术（GB-InSAR）的大坝边坡形变监测。为了有效监测边坡形变情况，在详细分析地基干涉雷达测量技术的原理及关键技术的基础上，将其应用于该边坡稳定性的实时监测中。结果表明，采用地基雷达干涉测量技术对该边坡稳定性进行实时监测，可以有效監测边坡形变情况，达到边坡灾害预警的目的。

关键词 大坝边坡;形变监测;地基雷达干涉测量技术;灾害预警

Abstract In order to solve the problems of traditional dam slope deformation monitoring methods， such as only obtaining point deformation information， long monitoring period， low monitoring accuracy and unable to carry out realtime monitoring， this paper took the Hongshiyan slope of Niulanjiang River in Yunnan Province as an example， and put forward the deformation monitoring of dam slope based on groundbased radar interferometry （GBInSAR） technology. In order to effectively monitor the deformation of slope， the principle and key technology of groundbased radar interferometry was analyzed in detail and applied in the realtime monitoring of slope stability. The results showed that the realtime monitoring of slope stability by groundbased radar interferometry could effectively monitor the deformation of slope and achieve the purpose of early warning of slope disasters.

Key words Dam slope;Deformation monitoring;GBInSAR;Disaster warning

我国是全世界名副其实的坝工大国之一。然而由于自然灾害、地貌因素、水文因素、气候因素以及人为因素等方面原因，导致我国每年都会面临边坡崩塌、边坡滑坡、泥石流等严重的地质灾害，给社会的生态环境和人民的生命财产安全造成了巨大的损失。因此，通过对边坡形变情况进行实时监测，并将监测到的形变数据进行实时分析与处理，对于准确掌握边坡的形变情况、形变特征和形变趋势以及制定安全的防范措施及应急预案有着十分重要的意义。

传统的大坝边坡形变监测手段有全站仪、水准仪、测量机器人等，虽然这些测量手段已经得到了广泛应用，但其经常无法进行实时监测，精度易受气候时间、气候、通视等条件限制。

合成孔径雷达[1-2]（synthetic aperture radar，SAR）是一种高分辨率雷达，其特点是分辨率高，能进行全天时、全天候（24 h、全天候条件）、大范围（数公里全面覆盖）、高精度（监测精度达亚毫米）的实时自动化监测，可以在十分恶劣的环境条件下识别监测物体并获取监测物体的高分辨率雷达影像。差分干涉雷达（D-InSAR）是SAR的一个重要应用，该技术已经在近十余年内得到了快速发展和应用[3-4]。合成孔径雷达干涉技术（InSAR）主要包括星载和地基合成孔径干涉技术，星载InSAR技术可以全天时、全天候工作，能够获取大范围区域的地表沉降变形监测数据，但是由于星载InSAR技术时常受大气干扰，监测精度低，且重返周期长，分辨率低，连续监测能力差，因此导致星载InSAR技术无法在大范围区域进行高精度的实时监测。地基干涉合成孔径雷达技术（GB-InSAR）不仅具有连续观测性，而且观测周期短，监测精度高（可达亚毫米），分辨率高，监测点可达百万级。该技术属于非接触性测量方法范畴，适用于对边坡进行实时变形监测[5]。笔者以云南牛栏江堰塞湖红石岩边坡为例，针对红石岩边坡的地貌、水文等特点，采用GB-InSAR技术对红石岩边坡进行形变监测，同时结合边坡的形变数据以及数据处理分析的结果，进一步对该技术在红石岩边坡变形监测方面的可行性进行分析与评价。

1 GB-InSAR的测量原理

GB-InSAR是采用一种微波干涉技术的创新雷达，集成合成孔径雷达（SAR）、干涉测量技术（D-InSAR）和步进频率连续波技术（SFCW）等多种先进技术。其基本原理是通过SAR技术来提高GB-InSAR系统的方位向分辨率，通过SFCW技术来提高GB-InSAR系统的距离向分辨率，通过干涉测量技术获取GB-InSAR系统的高精度视线向形变。GB-InSAR技术具有高精度、高空间分辨率、高采样频率和多角度观测等突出的技术优势[6]。

1.1 SAR技术的基本原理

合成孔径雷达成像是一种高分辨率微波遥感成像技术，该技术可以对边坡进行全天时、全天候监测、体积小、监测周期短、监测范围广、精度高、且受恶劣气候条件影响较小。合成孔径雷达（SAR）是利用小天线作为单个辐射单元，向某一固定方向移动，在不同的位置上接收同一监测物体返回的雷达信号并进行相关处理，进行成像[7]。通过小天线的运动形成一个等效的大天线，就可以获得高分辨率的星载合成孔径雷达图像（图1）。根据合成孔径雷达成像技术原理得，设GB-InSAR系统发射的信号带宽为，雷达信号波在空气中的传播速度为，则GB-InSAR系统的斜距向分辨率为c/2B，若GB-InSAR系统发射的信号波长为λ，系统在轨道上移动的最远距离为L，则GB-InSAR系统的方位向分辨率为λ/2L[8]。该研究采用的是MPDMR-05-LSA1701型号的GB-InSAR系统，该系统主要通过雷达信号接收器沿着滑动轨道进行移动，从而形成合成孔径效果，用于测量雷达小天线接收信号的幅度与相位信息，并通过差分干涉测量技术获取地基雷达监测区域的地形信息和相对形变信息，从而达到监测边坡形变的目的。GB-InSAR以固定的视角不断地发射和接收回波信号，经过聚焦处理后形成极坐标形式的二维SAR影像。在影像像元内，距离向分辨率是固定不变的，而方位向分辨率与像元夹角及目标距离有关，将距离向与方位向进行结合，监测区域被分为若干个二维小像元（图2），监测距离越远，方位向分辨率越低[9]。

1.2 D-InSAR原理

2 GB-InSAR监测流程及技术

2.1 GB-InSAR监测流程

红石岩边坡所处地势陡峭，地形复杂，外部山体形态变化不稳定，易受于人类活动和恶劣天气的影响，这些因素往往会对边坡监测产生极大的困难。因此，在监测前要整体规划，制定一个完整规范的监测流程。该工作流程主要包括：构建GB-InSAR监测体系、建立GB-InSAR监测系统、采集实验数据、数据处理、制作边坡形变图、建立边坡预警模型以及制定相应的安全措施和应急预案等。具体监测流程如图4所示。

2.2 GB-InSAR关键技术

GB-InSAR关键技术主要含有图像配准、生成干涉图、相位解缠、相位估计、地理编码等技术。首先图像配准技术是生成干涉图的基础，也是GB-InSAR形变监测的关键步骤之一，通过GB-InSAR系统在轨道上的2次滑动将所获取的SAR影像中同一监测物体的像元匹配到同一位置的点位上，并利用该点位在2幅不同的时间获取的SAR影像重叠的相位信息计算SAR影像的相干值即可完成SAR影像的图像配准。在完成图像配准后，通过提取由二维SAR影像数据生成高质量干涉图的干涉相位为相位解缠做准备。由于噪声以及平地效应的存在，在完成相位解缠工作前不仅需对干涉图进行去噪处理，而且还要还需要去除平地效應，从而可以获取更加直观的高程信息和稀疏的干涉条纹，为相位解馋工作提供便利。干涉图噪声及平地效应会严重影响GB-InSAR数据处理中SAR影像的图像质量，导致相位解缠工作无法实现。只有有效减少噪声及平地效应对干涉图的干扰，才能保证获取高质量的GB-InSAR形变图。相位解缠主要是将GB-InSAR获取的干涉相位的主值还原为真值，也是GB-InSAR形变监测的关键技术之一，相位解缠的准确性直接影响到GB-InSAR监测结果的精度[11]。大气相位校正是指当雷达电磁波信号经过大气层时，大气会对电磁波信号产生折射和散射影响，导致其传播路径延迟，从而形成大气效应，对干涉结果造成严重影响。因此，在数据处理过程中，需要考虑到大气相位对监测结果的影响，以提高监测精度。最后完成以上工作的同时，采用多源数据融合技术，将GB-InSAR技术与其他监测技术结合并分析出边坡综合形变信息，再对获取的综合形变监测结果进行地理编码，并将雷达坐标系中的识别到的边坡的点、线、面等属性映射到同一个地理坐标系下。

3 红石岩边坡监测试验与结果分析

3.1 试验区概况

云南牛栏江红石岩堰塞湖工程位于牛栏江下游河段，该工程左岸部分位于云南省昭通市巧家县境内，右岸部分位于鲁甸县境内，中心位置位于103°36′E、26°99′N。坝址距上游小岩头水电站厂房约23 km，距离下游天花板水电站取水坝约17 km，堰塞体长度约为910 m，高约103 m，总方量约为1 000万m3。该工程兼有防洪、供水、发电、灌溉、旅游等综合利用效益。枢纽工程上主要由堰塞体坝、右岸溢洪洞、右岸泄洪冲沙洞、右岸引水系统以及右岸岸边主副厂房等建筑物组成，设计灌溉面积2 400 hm2，设计供水114.5万人，装机容量20.1万kW，总投资35.9亿元。

该工程边坡为典型岩质边坡，无植被覆盖，边坡高度落差大，地势陡峭，该边坡是由2014年8月3日鲁甸6.5级地震所导致山体大面积崩塌从而形成的。目前边坡下游部位坡面岩体为地震拉裂、损伤岩体、卸荷拉张裂隙发育、岩体破碎。该边坡右岸大部为厚层、巨厚层栖霞、茅口组灰岩，在高程1 380 m附近分布有厚度20～40 m的梁山组页岩、泥岩夹砂岩。根据现场地质人员对整个边坡下游侧进行实地勘察，对开挖外围的地质条件做出了简单的评价及危岩体的圈定。边坡现状如图5所示。

该边坡危岩体主要分为Ⅰ类和Ⅱ类，Ⅰ类危岩区主要表现为平行边坡和垂直边坡的2组陡倾角裂隙切割，且宽度大，与坡面基本分离，呈孤立块体，稳定性极差;Ⅱ类危岩区主要表现为为平行边坡和垂直边坡的2组陡倾角裂隙切割，裂隙短小密集，但宽度不大，稳定系较差。

由于Ⅰ类危岩区岩体稳定性极差，将边坡右岸下游侧（1 765 m高程以下）危岩体细分分为7块，即Ⅰ-6、Ⅰ-7、Ⅰ-8、Ⅰ-9、Ⅰ-10、Ⅰ-11、Ⅰ-12。Ⅱ类危岩区稳定性较差，将边坡右岸下游侧（1 765 m高程以下）危岩体细分分为4块，即Ⅱ-1、Ⅱ-2、Ⅱ-3、Ⅱ-4。该边坡危岩体等级具体划分如表1所示。

3.2 GB-InSAR的作业条件

为了更好地监测目标区域，在选择观测房时，应考虑到以下条件[12]：

①持续供电。为了保证设备可以连续24 h进行监测，GB-InSAR系统需要持续供电。因此，在安装GB-InSAR设备时，应保证对该设备的持续供电。若现场出现短时间断电情况，应采用UPS继续供电。

②监测距离及范围。在施工现场布设过程中，应根据施工现场所处位置的地貌、水文等条件将雷达设备的监测距离及监测范围控制在合理的范围内，雷达设备最远监测距离应小于5 km，最大覆盖范围应小于10 km2，在此基础上，监测距离越大，雷达接收回波信号越弱，其误差效果越大，监测效果越差，因此应根据现场实际情况选择合适的监测距离。

③通视条件。GB-InSAR系统应该在能见度高的条件下工作。若能见度较低，则会对GB-InSAR系统的监测结果造成一定的误差，影响其监测精度。同时，在雷达设备与监测目标之间不能有障碍物的存在，比如岩石、树木、设备、无关人员等。如果存在障碍物，监测目标的反射强度会减弱，最终会影响监测设备的数据处理结果。

④设备安放点稳定性。在监测过程中雷达监测设备应平稳放置，不能受到振动。因此，在选址过程中要考虑将监测房安置于水平地面上。

⑤仪器架设位置。观测边坡时，应对设备视线方向和主滑移方向的夹角进行权衡，夹角越小，监测雷达对形变信号的强度越敏感，但不利于接收回波信号。

⑥滑坡植被。GB-InSAR系统的雷达波频率较高，波长较短，其雷达信号对目标物体的形变比较敏感，但如果在监测区域内存在大量的植被，由于空气流动的影响，植被会发生摇摆，雷达无法很好地接收监测目标的形变信号，最终会影响到GB-InSAR系统的监测效果。因此，监测区域应选择植被稀疏的区域。

3.3 试验方案与数据分析

根据该边坡的实际条件，结合被监测边坡的地质条件及危险区域的分布范围，将GB-InSAR设备安置于被监测边坡正南方向的地基较稳定处。为了确保GB-InSAR设备在监测过程中不受降水、大风等恶劣天气条件的影响，并满足其在恶劣天气条件下仍然正常连续监测的要求，在红石岩边坡正南方向相对稳定的基岩上建立监测房。将GB-InSAR设备安置于监测房中对边坡进行形变监测，监测房的建立可以有效地避免设备因不利天气条件而对监测结果产生的影响，保证设备24 h连续监测。同时，将监测边坡的系统服务器放置于监测房中，监测人员可以通过客户端对GB-InSAR设备的监测结果进行数据分析以及实时预警。GB-InSAR设备采用24 h连续监测，测量频率7 min/轨，水平方向观测角度-34.6°～34.6°，监测雷达的中心视线与水平夹角约为25°，距离向分辨率0.2 m，方位向分辨率为5.5 mrad。最后，对边坡监测区域的形变情况进行实时监测及实时数据分析与处理[13]。

该监测站于2018年12月20日开始部署，并于2018年12月21日开始试运行且对红石岩边坡進行监测，试运行时间段为2018年12月21日至12月31日。该设备于2019年1月1日正式运行，该研究选取监测时间段为2019年1月1日至1月31日所采集到的数据并对其进行数据分析。地基监测雷达在监测期间共获取3 410幅SAR影像，该GB-InSAR设备观测距离为300～500 m，观测范围覆盖整个边坡区域。该设备主要针对监测区域的1～10号点位进行数据分析，如图6所示。

依据该边坡的形变结果，对该边坡形变区域进行了现场实地考察，该边坡右岸上方由于施工方正在进行清渣作业，导致部分石头和尘土下滑，导致7号与10号点位累计形变量增加。该边坡左岸下方由于施工方进行爆破作业，同时加上近期雨雪天气的干扰，导致8号与9号点位累计形变量增加。总体来看，GB-InSAR监测结果为施工现场的安全保障提供了有效的数据支持。红石岩工程项目部技术人员根据现场监测数据及时地进行了施工人员撤离并对边坡左岸上方崖壁进行了防护网加固工作，有效避免了边坡灾害的发生。

4 结论

该试验结果表明，该系统可以不受暴风雨、雾霾等恶劣天气条件的影响，能够全天候地进行低边坡实时监测，该系统在雷达视线上监测精度可达毫米级别，且具有距离远、范围大、连续空间覆盖、全自动等优点。笔者通过MPDMR-05-LSA1701型号的GB-InSAR系统对云南牛栏江红石岩边坡进行监测试验，并通过试验分析发现边坡7～10号点位累计形变量较大，并根据试验数据及时对边坡左岸上方崖壁进行了防护网加固措施，有效避免了边坡灾害的发生。

参考文献

[1] BERARDINO P，COSTANTINI M，FRANCESCHETTI G，et al.Use of differential SAR interferometry in monitoring and modelling large slope instability at Maratea（Basilicata，Italy）[J].Engineering geology，2003，68（1/2）：31-51.

[2] 王桂杰，谢谟文，邱骋，等.D-InSAR技术在大范围滑坡监测中的应用[J].岩土力学，2010，31（4）：1337-1344.

[3] ZHANG P，ZHAO Z.Evaluation of data applicability for DINSAR in areas covered by abundant vegetation[C]//The ISPRS Technical Commission III Midterm Symposium on "Developments，Technologies and Applications in Remote Sensing".北京：国家基础地理信息中心天地图工作部，2018.

[4] LI X，HUANG G M，KONG Q L.Atmospheric phase delay correction of DINSAR based on sentinel-1A[C]//The ISPRS Technical Commission III Midterm Symposium on "Developments，Technologies and Applications in Remote Sensing".北京：国家基础地理信息中心天地图工作部，2018.

[5] 杨红磊，彭军还，崔洪曜.GB-InSAR监测大型露天矿边坡形变[J].地球物理学进展，2012，27（4）：1804-1811.

[6] 麻源源，左小清，麻卫峰，等.地基雷达在水电站边坡形变监测中的应用[J].工程勘察，2018（12）：52-57，62.

[7] 张毅.基于InSAR技术的地表变形监测与滑坡早期识别研究：以白龙江流域中游为例[D].兰州：兰州大学，2018.

[8] 赵东寅，申其鸿，马海涛，等.国产地基合成孔径雷达监测预警系统在紫金山金铜矿露天采场边坡位移监测的应用[J].中国安全生产科学技术，2015（4）：54-58.

[9] 黄其欢，岳建平，贡建兵.GBInSAR隔河岩大坝变形监测试验[J].水利水电科技进展，2016，36（3）：47-51.

[10] 刘作利，刘景玉，申修强，等.唐山马兰庄铁矿露天开采边坡变形监测的GB-InSAR技术[J].现代矿业，2018（4）：165-170.

[11] 麻德明，刘焱雄，徐文学，等.GB-InSAR的监测流程和关键技术[J].海洋开发与管理，2018（8）：81-85.

[12] 刘斌，葛大庆，李蔓，等.地基InSAR技术及其典型边坡监测应用[J].中国地质调查，2018，5（1）：73-81.

[13] 于怀昌，王春磊，王硕楠，等.基于物联网技术的栾川魏家沟滑坡监测预警系统研究与实践[J].华北水利水电大学学报（自然科学版），2018，39（6）：46-50.