1、Geomatics Science and Technology 测绘科学技术测绘科学技术,2023,11(2),109-114 Published Online April 2023 in Hans.https:/www.hanspub.org/journal/gst https:/doi.org/10.12677/gst.2023.112012 文章引用文章引用:庞宇.利用 SBAS-InSAR 技术对辽河平原进行沉降监测J.测绘科学技术,2023,11(2):109-114.DOI:10.12677/gst.2023.112012 利用利用SBAS-InSAR技术对辽河平原进行沉降技术对
2、辽河平原进行沉降 监测监测 庞庞 宇宇 辽宁省自然资源卫星应用技术中心，辽宁 沈阳 收稿日期：2023年2月28日；录用日期：2023年4月5日；发布日期：2023年4月12日 摘摘 要要 针对大面积区域沉降监测难、耗时耗力成本高的问题，本文提出利用针对大面积区域沉降监测难、耗时耗力成本高的问题，本文提出利用SBAS-InSAR技术对面积大的区域进技术对面积大的区域进行沉降监测，该方法不仅效率高，同时相对于传统沉降监测方法能够有效降低成本，可有效获取监测区行沉降监测，该方法不仅效率高，同时相对于传统沉降监测方法能够有效降低成本，可有效获取监测区域全域的沉降信息。域全域的沉降信息。关键词关键词
3、SBAS-InSAR，辽河平原，沉降监测辽河平原，沉降监测 SBAS-InSAR Technology is Used to Monitor the Subsidence of Liaohe Plain Yu Pang Liaoning Natural Resources Satellite Application Technology Center,Shenyang Liaoning Received:Feb.28th,2023;accepted:Apr.5th,2023;published:Apr.12th,2023 Abstract In order to solve the probl
4、ems of difficult and costly settlement monitoring in large area,this pa-per proposes to use SBAS-InSAR technology for settlement monitoring in large area,at the same time,compared with the traditional settlement monitoring method can effectively reduce the cost,can effectively obtain the monitoring
5、area of the entire settlement information.Keywords SBAS-InSAR,Liaohe Plain,Settlement Monitoring 庞宇 DOI:10.12677/gst.2023.112012 110 测绘科学技术 Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativec
6、ommons.org/licenses/by/4.0/1.引言引言 辽河平原地域辽阔，面积约 3 万多平方公里，长期的油气资源、矿产资源(煤炭)、地下水和地热资源的持续开采，使得辽河平原区域内地质环境问题逐步显现，地面沉降等地质灾害现象呈加剧发展趋势。地面沉降由于具备影响范围大，持续时间长，形成缓慢，不易察觉的特性，一旦形成，就会对地质环境系统造成不可逆的损害，破坏地面稳定性，造成建筑物倾斜、开裂等，严重威胁城市健康安全。传统沉降监测手段很难满足辽河平原监测需求1。合成孔径雷达干涉测量(InSAR)技术逐渐成熟并得到了工程化应用，已经成为地表形变监测的主要技术手段，在全球及区域性地形测图、大尺
7、度地表形变监测中得到广泛应用，如地震前后的地表形变监测，由于地下资源过度开采等因素造成的城市地面沉降。因此本项目利用 SBAS-InSAR 技术，采用两种 SAR 数据源，对辽河平原整体地区进行沉降监测2 3 4 5。2.SBAS-InSAR 技术原理技术原理 SBAS-InSAR 方法可以通过限制空间基线和时间基线对满足条件的 SAR 影像进行干涉组合，能够有效降低因时间间隔较长而引起的差分干涉图去相干现象的影响。该方法除了将高质量的点目标作为地表形变信息的表达载体外，还能够充分利用 SAR 影像场景中具有一定相干性的分布式目标参与解译，进一步提高离散点位的监测密度，更加完整得反映形变区域的
8、空间形态和演化规律。假设存在 N+1 幅经配准的时序 SAR 影像，经基线组合选取后得到 M 个差分干涉对，且 M 满足条件()()11 2NMN N+，第i个差分干涉对的第x个像素的相位模型为6 7 8：,x idef x ix ix iorb x in x i=+(1)其中，,def x i为雷达视向上的形变相位；,.E x i为 DEM 误差造成的残余地形相位；,x i为大气影响相位；,orb x i为轨道信息不精确引入的相位；,n x i为其它噪声相位。对每个干涉图的模型相位进行相位解缠，通过奇异值分解(SVD)便能获取到地表形变信息9 10 11。3.数据处理数据处理 3.1.SAR
9、 数据源选择数据源选择 由于辽河平原地区面积大，且整个区域分为重点监测区域与非重点监测区域，因此为节约成本并同时满足监测需求，本项目采用 RADARSAT-2 及 Sentinel-1 两种 SAR 数据，其中 RADARSAT-2 属于商业采购部分，应用于重点监测区域，其他区域采用 Sentinel-1 数据，其是免费的。3.2.RADARSAT-2 数据简介及处理方法数据简介及处理方法 3.2.1.RADARSAT-2 数据简介数据简介 RADARSAT-2 卫星超宽精细(Extra Fine,XF)模式影像空间分辨率为 5 m，幅宽 125 km，采用 C 波段成像(波长 5.6 cm)
10、，兼具高分辨率、大幅宽以及复杂地物环境下的地面沉降监测能力。根据项目监测需求查询监测区的卫星影像覆盖情况，综合考虑研究区大小、数据空间覆盖、时间序列完整性、空间基线等因素，本次处理选择 2021 年 10 月影像至 2022 年 10 月期间两组共 16 景 SAR 数据。Open AccessOpen Access庞宇 DOI:10.12677/gst.2023.112012 111 测绘科学技术 3.2.2.RADARSAT-2 数据处理数据处理 针对 RADARSAT-2 数据，采用 GAMMA 结合 StaMPS 软件进行处理。数据预处理部分由 GAMMA处理，时序分析部分由 StaM
11、PS 处理12 13 14 15。处理流程见图 1。Figure 1.RADARSAT-2 data processing flow chart 图图 1.RADARSAT-2 数据处理流程图 成像处理：导入 RADARSAT-2 数据，将原始影像转换成软件可以处理的文件类型，获取每景影像对应的 slc 及 slc.par 参数文件。影像配准及外部 DEM 配准：根据时间基线、空间基线选取配准主影像，将其他影像与主影像进行配准。SAR 影像配准精度要求小于 1/8 像元。同时，将主影像与外部 DEM 进行配准，获取雷达坐标系下高程文件，为地理编码和差分干涉图做准备。干涉对组合与差分干涉图生成：
12、设置一定时间与空间基线阈值，将配准后的 SLC 影像进行干涉对组合，然后进行差分干涉处理，获取短基线集差分干涉图。高相干目标点提取：将生成的差分干涉图转换格式，导入 StaMPS 中进行后续处理。迭代计算每个差分干涉图中每个候选点的噪声相位，结合 SLC 影像幅度离差对高相干目标点进行初选。计算初选点的噪声标准差，剔除噪点及质量不好的点，获取高相干目标点。相位解缠：利用三维相位解缠方法对高相干目标点进行相位解缠。初步获取形变相位，此时形变相位包含轨道误差、DEM 残余误差和大气相位等，需要分步进行剔除。轨道、大气相位去除：对解缠后的形变相位进行最小二乘估计，估算 DEM 残余误差和主影像大气分
13、量，利用曲面拟合方法估计轨道残余相位，然后通过时空滤波方法去除大气影响，最终获取地表形变信息。3.3.Sentinel-1 数据处理数据处理 Figure 2.Sentinel-1 data processing flow chart 图图 2.Sentinel-1 数据处理流程图 数据导入与连接图生成：将原始卫星影像和精密轨道星历数据导入进行 TOPSAR 处理。生成 SLC数据以及索引文件、SLC 影像强度图、谷歌地球坐标文件、SARscape 的索引识别文件、图形矢量文件。选择 3D 解缠模式，平台会根据导入的数据的相对关系自动选取最佳的超级主影像，并根据影像数据选取合适的时间基线和空间
14、基线阈值进行配对，其他辅影像会以主影像作为参考影像，以此来生成辅影像与主影像的时间基线、空间基线信息，生成 SAR 的数据连接图。流程如图 2 所示。干涉工作流及连接图编辑：因为使用的是哨兵数据，所以多视比设为 4:1。SBAS-InSAR 干涉流参数选取很大程度上决定了干涉结果的好坏，设置参数要根据影像数量及测区地表环境相关。利用覆盖研究区的 DEM 数据对影像进行辅助配准。本项目测区地质条件复杂，相干性较低，因此选择 Goldstein 算法进行干涉图滤波。由于研究区域内含有多处积水区，同时草木植被覆盖丰富，因此选取 Delaunay 三角网方法进行解缠。为了提高相位解缠质量，大幅降低相位
15、解缠期间的错误率，可将分解等级调整为 2。庞宇 DOI:10.12677/gst.2023.112012 112 测绘科学技术 轨道精炼和重去平：这步工作主要是去除上一步成果中的轨道、相位坡道等误差。这一步最重要的步骤就是选取地面控制点(GCP)。选取 GCP 点的依据：要远离形变区，在相位平稳处选点；在相位连续处选点，要远离高频相位；要避开个别解缠错误的相位点。SBAS-InSAR 第一次反演：基于线性模型计算出所有像对的形变(日期、速度、加速度和加速度变化)和高程(校正值和新的 DEM)。SBAS-InSAR 第二次反演：计算时间序列上的位移，得到研究区地表的时序形变量。在第一步得到的形变
16、速率基础上，进行大气滤波，估算和去除大气相位，得到最终位移结果。地理编码：以上成果是在斜距坐标系下，此步可将成果的斜距坐标转换为地理坐标。本项目选择WGS-84 坐标系，并得到 LOS 方向的形变结果。输出的结果还包括对应时期的累积形变量、最终的平均形变速率以及平均形变速率精度、平均形变加速度等。4.成果分析成果分析 本项目共监测出沉降区域十三处，沉降最严重地区沉降速率可达到 90 mm/年。辽河平原整体沉降速率图见图 3。部分沉降严重区域沉降速率图见图 46。从图 4 中可看出，该区域大部分地区地表稳定，沉降值在 10 mm 以内，但中心区域的村庄发生了较大的沉降，最大可达到 60 mm；图 5 中大部分地区沉降值在 10 mm 以内，但有多个村庄沉降值达到 40 mm；图 6 中中心区域沉降达到 60 mm，且周围有多处沉降达到 30。正常地表沉降在 10 mm 以内我们都认为是正常现象，20 mm 以上是我们的关注区。通过外业踏勘我们发现，沉降达到 60 mm 以上的地区，基本都在地下资源开采区。因此，原沉降成因主要受三方面影响，一是地下煤炭资源开采。比如小康矿、大平矿的煤炭开采