利用SBAS-InSAR技术对辽河平原进行沉降监测.pdf

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技术对辽河平原进行沉降技术对辽河平原进行沉降 监测监测 庞庞 宇宇 辽宁省自然资源卫星应用技术中心，辽宁 沈阳 收稿日期：2023年2月28日；录用日期：2023年4月5日；发布日期：2023年4月12日 摘摘 要要 针对大面积区域沉降监测难、耗时耗力成本高的问题，本文提出利用针对大面积区域沉降监测难、耗时耗力成本高的问题，本文提出利用SBAS-InSAR技术对面积大的区域进技术对面积大的区域进行沉降监测，该方法不仅效率高，同时相对于传统沉降监测方法能够有效降低成本，可有效获取监测区行沉降监测，该方法不仅效率高，同时相对于传统沉降监测方法能够有效降低成本，可有效获取监测区域全域的沉降信息。域全域的沉降信息。关键词关键词 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 problems 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 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:/creativecommons.org/licenses/by/4.0/1.引言引言 辽河平原地域辽阔，面积约 3 万多平方公里，长期的油气资源、矿产资源(煤炭)、地下水和地热资源的持续开采，使得辽河平原区域内地质环境问题逐步显现，地面沉降等地质灾害现象呈加剧发展趋势。地面沉降由于具备影响范围大，持续时间长，形成缓慢，不易察觉的特性，一旦形成，就会对地质环境系统造成不可逆的损害，破坏地面稳定性，造成建筑物倾斜、开裂等，严重威胁城市健康安全。传统沉降监测手段很难满足辽河平原监测需求1。合成孔径雷达干涉测量(InSAR)技术逐渐成熟并得到了工程化应用，已经成为地表形变监测的主要技术手段，在全球及区域性地形测图、大尺度地表形变监测中得到广泛应用，如地震前后的地表形变监测，由于地下资源过度开采等因素造成的城市地面沉降。因此本项目利用 SBAS-InSAR 技术，采用两种 SAR 数据源，对辽河平原整体地区进行沉降监测2 3 4 5。2.SBAS-InSAR 技术原理技术原理 SBAS-InSAR 方法可以通过限制空间基线和时间基线对满足条件的 SAR 影像进行干涉组合，能够有效降低因时间间隔较长而引起的差分干涉图去相干现象的影响。该方法除了将高质量的点目标作为地表形变信息的表达载体外，还能够充分利用 SAR 影像场景中具有一定相干性的分布式目标参与解译，进一步提高离散点位的监测密度，更加完整得反映形变区域的空间形态和演化规律。假设存在 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 数据源选择数据源选择 由于辽河平原地区面积大，且整个区域分为重点监测区域与非重点监测区域，因此为节约成本并同时满足监测需求，本项目采用 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)，兼具高分辨率、大幅宽以及复杂地物环境下的地面沉降监测能力。根据项目监测需求查询监测区的卫星影像覆盖情况，综合考虑研究区大小、数据空间覆盖、时间序列完整性、空间基线等因素，本次处理选择 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处理，时序分析部分由 StaMPS 处理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 进行配准，获取雷达坐标系下高程文件，为地理编码和差分干涉图做准备。干涉对组合与差分干涉图生成：设置一定时间与空间基线阈值，将配准后的 SLC 影像进行干涉对组合，然后进行差分干涉处理，获取短基线集差分干涉图。高相干目标点提取：将生成的差分干涉图转换格式，导入 StaMPS 中进行后续处理。迭代计算每个差分干涉图中每个候选点的噪声相位，结合 SLC 影像幅度离差对高相干目标点进行初选。计算初选点的噪声标准差，剔除噪点及质量不好的点，获取高相干目标点。相位解缠：利用三维相位解缠方法对高相干目标点进行相位解缠。初步获取形变相位，此时形变相位包含轨道误差、DEM 残余误差和大气相位等，需要分步进行剔除。轨道、大气相位去除：对解缠后的形变相位进行最小二乘估计，估算 DEM 残余误差和主影像大气分量，利用曲面拟合方法估计轨道残余相位，然后通过时空滤波方法去除大气影响，最终获取地表形变信息。3.3.Sentinel-1 数据处理数据处理 Figure 2.Sentinel-1 data processing flow chart 图图 2.Sentinel-1 数据处理流程图 数据导入与连接图生成：将原始卫星影像和精密轨道星历数据导入进行 TOPSAR 处理。生成 SLC数据以及索引文件、SLC 影像强度图、谷歌地球坐标文件、SARscape 的索引识别文件、图形矢量文件。选择 3D 解缠模式，平台会根据导入的数据的相对关系自动选取最佳的超级主影像，并根据影像数据选取合适的时间基线和空间基线阈值进行配对，其他辅影像会以主影像作为参考影像，以此来生成辅影像与主影像的时间基线、空间基线信息，生成 SAR 的数据连接图。流程如图 2 所示。干涉工作流及连接图编辑：因为使用的是哨兵数据，所以多视比设为 4:1。SBAS-InSAR 干涉流参数选取很大程度上决定了干涉结果的好坏，设置参数要根据影像数量及测区地表环境相关。利用覆盖研究区的 DEM 数据对影像进行辅助配准。本项目测区地质条件复杂，相干性较低，因此选择 Goldstein 算法进行干涉图滤波。由于研究区域内含有多处积水区，同时草木植被覆盖丰富，因此选取 Delaunay 三角网方法进行解缠。为了提高相位解缠质量，大幅降低相位解缠期间的错误率，可将分解等级调整为 2。庞宇 DOI:10.12677/gst.2023.112012 112 测绘科学技术 轨道精炼和重去平：这步工作主要是去除上一步成果中的轨道、相位坡道等误差。这一步最重要的步骤就是选取地面控制点(GCP)。选取 GCP 点的依据：要远离形变区，在相位平稳处选点；在相位连续处选点，要远离高频相位；要避开个别解缠错误的相位点。SBAS-InSAR 第一次反演：基于线性模型计算出所有像对的形变(日期、速度、加速度和加速度变化)和高程(校正值和新的 DEM)。SBAS-InSAR 第二次反演：计算时间序列上的位移，得到研究区地表的时序形变量。在第一步得到的形变速率基础上，进行大气滤波，估算和去除大气相位，得到最终位移结果。地理编码：以上成果是在斜距坐标系下，此步可将成果的斜距坐标转换为地理坐标。本项目选择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 以上的地区，基本都在地下资源开采区。因此，原沉降成因主要受三方面影响，一是地下煤炭资源开采。比如小康矿、大平矿的煤炭开采在三台子水库附近形成大面积的沉降区；调兵山市多个煤矿在全市大范围内引发整体沉降；西马煤矿在西马峰镇的西大堡村和野老村等地形成了严重的沉降区。二是地下油气资源开采。辽河油田有大量采油厂，多年以来的持续开采，致使周边区域地下油气储层等液体资源储层大面积采空，且开采范围逐步扩大，地面沉降范围和沉降量值也将逐渐随时间扩大。同时由于油气储层埋藏深，采空区面积巨大，因而一般以缓慢、累进性的整体沉陷为主。三是地下水资源开采。根据资料实地调查取证，沉降区域内大面积种植水稻、玉米等农作物以及草莓等经济作物，在农作物生长周期内需水量巨大。调查区内水位埋深浅，平均在 0.55 m 左右，因其灌溉水源和饮用水源往往来自于深层地下承压水的开采，且其开采之后并无回灌，造成沉降区承压水位持续下降，土体有效应力增加产生压缩变形，从而表现为以开采中心为原点的地面沉降漏斗。Figure 3.Overall settlement rate of Liaohe Plain 图图 3.辽河平原整体沉降速率图 庞宇 DOI:10.12677/gst.2023.112012 113 测绘科学技术 Figure 4.Monitoring results of settlement area I 图图 4.沉降区域一监测成果 Figure 5.Monitoring results of settlement area II 图图 5.沉降区域二监测成果 Figure 6.Monitoring results of settlement area III 图图 6.沉降区域三监测成果 庞宇 DOI:10.12677/gst.2023.112012 114 测绘科学技术 5.结束语结束语 通过本次项目作业，得出以下结论：1)SBAS-InSAR 技术对于地表沉降变化的监测是有效的，其结果的精度基本可以达到毫米级。2)SBAS-InSAR 技术特别适合对面积大的区域进行沉降监测。3)SBAS-InSAR 技术对于地表发生沉降变化尤其是长时间的缓慢变化非常敏感。4)SBAS-InSAR 技术成本低，目前国际上有哨兵免费的 SAR 数据源。5)本次作业发现辽河平原地区确实存在多处沉降区，其中沉降最大的位置在辽河油田，其次是煤矿区。对于地下资源开采的同时，采取科学的方法减少对我们生活环境的破坏是必要的。参考文献参考文献 1 汪钰晴,唐伶俐,王新鸿,周增光,李子扬,李传荣.基于时序 InSAR 技术的中巴走廊盖孜河谷段地表形变分析J.中国科学院大学学报,2022,39(3):343-351.2 刘广,郭华东,Ramon Hanssen,等.InSAR 技术在矿区沉降监测中的应用研究J.国土资源遥感,2008(2):51-55.3 李珊珊,李志伟,胡俊,等.SBAS-InSAR 技术监测青藏高原季节性冻土形变J.地球物理学报,2013,56(5):1476-1486.4 王桂杰,谢谟文,邱骋,等.D-INSAR 技术在大范围滑坡监测中的应用J.岩土力学,2010,31(4):1337-1344.5 邓辉,黄润秋.InSAR 技术在地形测量和地质灾害研究中的应用J.山地学报,2003,21(3):373-377.6 刘国祥,刘文熙,黄丁发.InSAR 技术及其应用中的若干问题J.学术动态,2001(8):10-12.7 葛大庆,王艳,郭小方,等.基于相干点目标的多基线 D-InSAR 技术与地表形变监测J.环境遥感学术年会.遥感学报,2007,11(4):144-150.8 赵超英,张勤,丁晓利,等.利用 InSAR 技术定位西安活动地裂缝J.武汉大学学报:信息科学版,2009,34(7):5.9 陈琦,刘江,鲁大江.InSAR 技术在城市更新中的应用研究J.北京测绘,2022,36(5):619-623.10 张训虎,章磊,郝树宾,等.合成孔径雷达干涉(InSAR)测量技术应用及展望J.北京测绘,2014(2):32-35.11 管力,唐伟,戴华阳,等.采用 SBAS-InSAR 技术监测郑州城区地面沉降J.北京测绘,2019,33(4):462-467.12 杨红磊,彭军还,张丁轩,等.轨道误差对 InSAR 数据处理的影响J.测绘科学技术学报,2012,29(2):118-121,126.13 何秀凤,何敏.InSAR 对地观测数据处理方法与综合测量M.北京:科学出版社,2012.14 吴斌,高旭芝,李森,等.InSAR 沉降监测量检验方法的研究简J.北京测绘,2019,33(2):212-215.15 王彦军.关于 InSAR 沉降监测结果精度比对的方法研究J.北京测绘,2017(5):109-111.