基于SBAS-InSAR的采煤区地表形变监测与分析.pdf

第46 卷第9期2023年9月地测与水害防治基于SBAS-InSAR的采煤区地表形变监测与分析煤炭与化工Coal and Chemical IndustryVol.46 No.9Sep.2023宋旭斌，路鑫（潞安化工集团有限公司古城煤矿，山西长治0 46 0 0 0）摘要：山西长治地区煤炭资源丰富，因地下开采引起的地质灾害频发，针对基于“点”的传统形变测量技术难以胜任大范围地表的监测任务问题，本文以2 9景Sentinel-1A雷达影像为数据源，利用小基线集（SBAS-InSAR）技术获取了2 0 2 1年10 月12 日至2 0 2 2 年10 月7日期间古城煤矿S1306工作面开采引起的地表沉降结果，并结合特征点的时序累积沉降信息对研究区地表形变特征和时空演化过程做了细致分析。结果表明，S1306工作面周边地表自2022年1月16 日开始发生明显形变，并逐渐在始采线附近形成一个沉陷盆地，盆地中心可探测到的最大下沉值约-17 8.6 mm，且盆地有逐渐向西北方向发育的趋势；由于工作面开采尚未结束，截止2 0 2 2 年10 月7 日，形变影响沿采动方向仅传递至工作面中部位置，与实际开采位置基本吻合，说明SBAS在矿区形变监测中具有良好的适用性。研究成果对矿区合理开采和地表建（构）筑物的保护具有重要的现实意义。关键词：小基线集；开采沉陷；形变监测；Sentinel-1A；沉降漏斗中图分类号：TD163文献标识码：A文章编号：2 0 9 5-59 7 9（2 0 2 3）0 9-0 0 51-0 6Monitoring and analysis of surface deformation incoal mining area based on SBAS-InSARSong Xubin,Lu Xin(Gucheng Mine,Luan Chemical Group Co.,Ltd,Changzhi 046000,China)Abstract:The coal resources in Changzhi area of Shanxi Province are abundant,and the geological disasters caused byunderground mining are frequent.Aiming at the problem that the traditional deformation measurement technology based onpoint was dificult to be competent for the monitoring task of large-scale surface,this paper took No.29 field Sentinel-1 Aradar images as the data source,and used the small baseline set（SBA S-In SA R )t e c h n o l o g y t o o b t a i n t h e s u r f a c esubsidence caused by the mining of No.S1306 Face in Gucheng Mine from October 12,2021 to October 7,2022.Combinedwith the time-series cumulative setlement information of feature points,the surface deformation characteristics andspatial-temporal evolution process of the study area were analyzed in detail.The results showed that the surface around No.S1306 Face had undergone obvious deformation since January 16,2022,and a subsidence basin had gradually formed nearthe initial mining line.The maximum subsidence value detected in the center of the basin was about-178.6 mm,and thebasin had a tendency to gradually develop to the northwest.Because that the mining of the working face had not beencompleted,as of October 7,2022,the deformation influence was only transmitted to the middle position of the working facealong the mining direction,which was basically consistent with the actual mining position,indicating that SBAS had goodapplicability in the deformation monitoring of the mining area.The research results had important practical significance forthe rational mining of mining areas and the protection of surface buildings.Key words:small baseline set;mining subsidence;deformation monitoring;sentinel-1A;settlement funnel责任编辑：张彤D0I:10.19286/ki.cci.2023.09.013作者简介：宋旭斌（197 2 一），男，山西长治人，高级工程师。引用格式：宋旭斌，路鑫.基于SBAS-InSAR的采煤区地表形变监测与分析 J煤炭与化工，2 0 2 3，46（9）：51-56.512023年第9期0引言煤炭作为我国的主体能源，长期担负着国家能源安全和经济高速发展的重任。但大规模、高强度的地下开采往往会破坏地层的原始应力平衡，而且随着采空区范围扩大，直接顶在自身重力和上覆岩层荷载作用下破断，最终引起地表的移动变形。因开采导致的地表沉陷易造成建筑物墙体开裂、交通道路损坏、农田积水等地质灾害问题，严重制约了矿区城市的可持续发展。因此，对采空区地表进行实时精密形变监测，获取动态沉降信息势在必行。常规开采监测手段主要有精密水准测量和全球导 航系统（Global Navigation Satellite System,GNSS)，但在实际工程中往往面临观测点易丢失、观测精度易受干扰、成本投人巨大、时空分辨率低等问题，无法准确获取矿区整体的形变特征。合成孔径雷达干涉测量(Interferometry Synthetic ApertureRadar,InSAR)是一种新型的对地观测技术，它仅对同一区域的多景SAR影像进行复共轭相乘即可提取大范围、毫米级的地面形变信息，具有全天时、高精度、自动化等优势。而小基线集（SmallBaseline Subsets InSAR,SBAS-InSAR）技术在差分InSAR（D i f f e r e n t i a l In SA R,D-In SA R）基础上提出，有效克服了后者易受大气效应、时空失相干影响的缺陷3，它通过对SAR影像上的高相干点进行相位模型建立和相位变化分析，即使在永久散射体PS(Persistent Scatterer）密度较低的矿区也能获取更为精确的长时序形变结果。许多学者都曾将 SBAS技术运用于矿区的形变监测，并取得了一系列丰硕的成果。李达等 4利用SBAS技术对13景TerraSAR-X进行处理，获取了榆林某矿的地表沉降信息，发现地表下沉值与时间呈线性关系；栾元重等 借助19景Sentinel-1A影像和SBAS技术对沉陷漏斗区域内的工作面进行时序监测，并利用水准数据验证了SBAS技术的可靠性；柴华彬等回提出利用反距离加权插值方法(Inverse Distance Weighted,IDW）将实测水准数据与SBAS解译结果相融合，以解决部分区域因失相干导致InSAR形变值缺失的问题；张香凝等以蒲河矿区为例，结合SBAS技术和数值模拟共同探究了煤矿区地面沉降的时空变形规律和机制。山西长治地区煤炭资源丰富，但因地下开采引起的地质灾害频发，呕需开展相关的地表沉陷监测工作。本文在前人研究的基础上，借助52煤炭与化工SBAS-InSAR技术对2 9 景Sentinel-1A影像进行解译处理，以期获取古城煤矿S1306工作面在2 0 2 1年10 月12 日至2 0 2 2 年10 月7 日期间的时序沉降结果和开采影响范围。研究成果对矿区开采规划和地表塌陷的预警防治均具有重要的现实意义。1研究区与数据源1.1研究区概况本文研究区位于长治市古城煤矿东南方向一采区，整体地势呈西高东低，地貌类型以河谷平原为主。研究区分布着大面积农田、村庄房屋以及电力线等重要的基础设施，且覆盖有较厚的第四系黄土沉积；地下山西组3号煤层目前有多个工作面正在回采，其中S1306工作面走向长约198 2 m，倾向宽约36 3m，平均采厚6.4m，煤层倾角约为0 6，底板标高-352 一-39 7 m。开采方法以倾斜长壁式综采法为主，全部落法管理顶板。目前，该采区周边已出现大量明显的地裂缝及路面损坏情况，呕需对该区域进行地表沉陷监测工作。研究区的地理位置如图1所示。1125550E112570E1125810EN.0.S1.96N.0.ST.96N.0.t190图例1125550E112570E1125810E图1研究区地理位置Fig.1 Geographical location of the study area1.2数据源本文从欧空局(European Space Agency,ESA)哥白尼开放数据访问中心（https:/scihub.copernicus.eu/dhus/#/home）收集了2 9 景Sentinel-1A升轨单视复数影像，时间跨度为2 0 2 1年10 月12 日至2022年10 月7 日，工作模式均为干涉宽幅(Interferometric Wide-swath，IW）模式，该模式基于递进式地形扫描(Terrain Observation withProgressive Scans SAR，T O PSA R)技术，能够将 3个子条带合并成1景高质量的SAR雷达数据叫，其地面分辨率可达5m20m。实验影像的其他参数见表1。另外，在数据处理过程中，引人了POD第46 卷42.80098550宋旭斌等：基于SBAS-InSAR的采煤区地表形变监测与分析精密轨道文件用来去除卫星轨道的系统性误差，并利用30 m分辨率的STRM1DEM数据消除了干涉SAR成像时间极化方式幅宽/km人射角/）车20211012、2 0 2 110 2 4、2 0 2 1110 5、2 0 2 11117、20211129、2 0 2 112 11、2 0 2 112 2 3、2 0 2 2 0 10 4、20220116、2 0 2 2 0 12 8、2 0 2 2 0 2 0 9、2 0 2 2 0 2 2 1、20220305、2 0 2 2 0 317、2 0 2 2 0 32 9、2 0 2 2 0 410、20220422、2 0 2 2 0 50 4、2 0 2 2 0 516、2 0 2 2 0 52 8、20220609、2 0 2 2 0 6 2 1、2 0 2 2 0 7 0 3、2 0 2 2 0 7 15、20220727、2 0 2 2 0 8 0 8、2 0 2 2 0 8 2 0、2 0 2 2 0 90 1、2技术原理与数据处理流程2.1SSBAS-InSAR原理小基线集技术最早由国外学者Berardino等人在2 0 0 2 年提出，其基本原理如下 12-15：假设选取（t o,t i,t 2,t n）时段内覆盖同一区域的N+1幅SAR影像，基于设定的时空基线阈值，组合成M幅多视差分干涉对，其中M的取值范围满足：d,(c,r)=0(ts x,)-d(s,x,4 d式中：入为雷达波长；中(tB，x,r)和(ta,x,r)分别表示2 幅影像上的相位值；d(tB,x,r)和d(t,x,r)分别为 ta和 ts时刻相对于初始时刻to在雷达视线向(Line of Sight,LOS)上的形变，一般假设 to 时刻的形变相位值为O；U;为ts至ts时段内的平均形变速率。将用于干涉处理的M对影像数据按照时序进行排列，则所有干涉图相位可表示为如下矩阵形式：8 d=BV式中：B为1个MN的系数矩阵；V为速度矢量。由于矩阵B的秩小于N，根据最小二乘法无法得到唯一的结果，因此需要采用奇异值分解法(Singular Value Decomposition,SVD）联合多个基线集，来解决法方程的秩亏问题，并得到速度V的广义逆解。之后将各个时间段内的形变速率在时间域上进行积分，即可获取整个监测时段内地表在LOS向上的累积时序形变结果。2.2SBAS-InSAR数据处理流程SBAS-InSAR数据处理关键步骤和参数设置如下：选取2 0 2 1年11月2 9日的影像作为主影像，其余从影像与之配准，并设置空间基线百分比阈值2023年第9期图中地形相位的影响。表1Sentinel-1A基本参数Table 1 Basic parameters of Sentinel-1AVV(垂直极化)20221007干涉处理，并生成了第j幅已解缠的差分干涉图，在去除地形相位、大气延迟相位和噪声误差之后，任意像元位置（x,r)（和r分别为方位向和距离向坐标）处的差分干涉相位可表示为：入为45%，时间基线阈值为6 0 d，多视比为4：1(距离向：方位向），最后共生成10 4个干涉对；利用最小费用流(Minimum Cost Flow,MCF)方法进行相位解缠，Goldstein滤波技术去除干涉图噪声，并逐个检查调整不理想的像对；参考干涉处理步骤生成的相干图和解缠相位图，选取2 0 30 个稳定且分布均匀的地面控制点（Ground Control Point,GCP）完成轨道精炼，去除残余轨道误差相位；选择线性模型（LinearModel）进行二次解缠，并(3)估算平均位移速率和残余地形相位；根据大气延迟相位、形变相位和噪声相位表现出的不同的时空特性，利用时间域的高通滤波（36 5d）和空间域的低通滤波（12 0 0 m）分离并剔除大气延迟误差，得到覆盖整个观测时间的地表形变时间序列；将最终解译结果从斜距坐标系转换至WGS-84坐标系，即完成地理编码。本文的数据处理流程如图2所示。3解译结果与分析3.1室空间演化分析利用SBAS-InSAR技术对覆盖S1306工作面的29景影像进行解译处理，获取了监测时段内研究区LOS向的时序累积形变结果。之后忽略水平位53轨道号(Path)25038.97N+12MN(W+1)2若在ta 和 ts（t b t a）时刻对2 幅SAR 影像进行(2)图幅号(Frame)113111(1)2023年第9期移的影响，利用公式d，=d i o s /c o s o 将解译结果从视线向投影至垂直向，其中为卫星人射角。另N+1.景SLC数据煤炭与化工外，解译结果中的正值表示地面抬升，负值表示地面下沉。第46 卷主辅影像选取时空滤波估算位移速率和残余地形生成干涉对STRMDEM滤波和去平相位解缠和辐射定标去除残余轨道和平地相位S1306工作面自2 0 2 2 年1月开始回采，目前仅开采到距离始采线约8 50 m的位置。但图3显示，在监测时段内，该工作面周边地表已经发生了明显沉陷。在形变监测的初期阶段，即2 0 2 1年10 月12日至2 0 2 1年12 月2 3日（图3a图3c），研究区地表总体较为稳定，零星区域发生小幅沉降，量值约-13mm，产生该现象的原因可归结为工作面所处的农田区域相干性较差，导致 SBAS在反演形变的过程中发生解缠误差，影响了结果的准确性；从2022年1月16 日起（图3d），工作面的始采线位置开始出现明显下沉，最大沉降量约-32.6 mm，112560E112570ECn20211024N.0.S19%112560E112570WE1120560WE20570E(d 2022016N.0.S109%奇异值分解(SVD)去除大气相位和残余地形形变时间序列诺干个GCP地理编码图2 SBAS数据处理流程Fig.2 Data dealing process of SBAS符合工作面开始回采的时间；随着工作面持续推进，地表沉降范围逐渐朝西南方向（即工作面开采方向）扩张；到了2 0 2 2 年5月16 日（图3h），始采线附近地表已经形成了一个较为明显的沉降盆地，盆地中心可探测到的最大下沉值约-17 8.6mm；随着采动影响逐渐传递至工作面中部位置（图3i图3k），始采线区域的下沉速度逐渐趋于平缓，但仍处在下沉状态；直至2 0 2 2 年10 月7 日(图3m），此时S1306工作面的开采影响边界稍超过实际开采进度位置，但总体与实际情况相符，说明SBAS技术在矿区形变监测中具有一定的可靠性和准确性。112560E112570*Eby2021117N.0.S1.90N.0.S198N.0.tL.96N.0.t109(e)20220221112560E112570E(e2021123N.0.ST.9N.0.t1N.O.tI112560E112570E112056:0WE1120N.0.S196o98112560E112571OWE112056:0WE112002020317N.0.S.90570E570WEN.0.S19N.0.t9654N.0.tIN.0.t1112560*E112570E9N.0.tIN.O.tI112560E112570E98N.O.tI112560E112570E宋旭斌等：基于SBAS-InSAR的采煤区地表形变监测与分析2023年第9期112560E112570*E(220220410N.0.S190112560E112570E02020516N.0.S1096N.0.S196112560E112570Ei20220621N.0.S1.96N.0.S1096N.0.S196N.O.t1N.0.t19619100128560E112112.360E112。3 7:0 E20221715N.0.ST.98N.0.t1N.0.t10981.00118560570E560E112570E1ii25704202080N.0.S1090N.0.S1096N.0.t19002022090N.0.S196N.0.S1096112560E112570E112560E112570EN()2021007人图例S1306工作面N.0.S1096累积形变量/mm-314-188-187-129-128-86.786.655.2-55.1-33.1N.0.t90-33-16.4-16.3-3.483.477.6117.62-24.224.3-157112560E112570E图3 S1306工作面时序形变Fig.3 Time series deformation of No.S1306 Face3.2日时间演化分析为了进一步探究研究区地表在采动影响下的时序形变过程，在S1306工作面的走向和倾向上共选取8 个特征点（编号点A点H)，并以2 0 2 1年10月12 日影像为参考基准，生成各个点的时序累积沉降曲线。具体点位分布如图4所示。由图4可知，在2 0 2 1年10 月12 日至2 0 2 2 年1月16 日期间，走向线上4个特征点没有发生明显形变，由于点A和点B分别位于始采线和沉陷盆地中心，易受开采扰动，这一时段内两点均出现小幅度下沉，最大沉降值约-10 mm；而点C和点D分布在工作面中部和停采线的位置，点位表现出反复微弱抬升的迹象，可将这种情况其视为SAR图像噪声带来的解算误差。2022年1月16 日之后，点A和点B开始持续加速下沉，且分别于2 0 2 2 年3 月17 日和2 0 2 2 年5月16 日，其形变速率才逐渐放缓，但形变量仍在持续增加，这是由于采动作用对始采线区域地表112560E112570E的影响即将达到“临界点”，该区域的形变经历了“开始下沉一剧烈下沉一缓慢停止下沉”的过程，一定程度上符合常规开采沉陷规律。点C和点D在2 0 2 2 年8 月8 日之前，下沉曲线吻合度较高，但之后点C的沉降量陡增至-8 7.6 mm，说明工作面中部区域此时开始受采动的强烈影响，与实际工程进度一致。同时，发生迅速沉降的原因还可能是由于点处在土质松软的农田，在夏季强降水的渗透下，土层的抗剪强度会大大降低，从而引起地表的湿化崩解和垮落12 。倾向线上的4个特征点均匀分布在工作面的两侧，在2 0 2 2 年1月16 日之前，由于点E和点H远离开采中心，点F和点G位于沉陷盆地边缘，因此后者呈现小幅持续下沉的态势，前者沉降与起伏交替，受采动影响较小；之后点F和点G开始加速下沉，但点F的沉降幅度与速率要大于点G。同样，点E和点H于2 0 2 2 年4月10 日也开始经历一个缓慢下沉的过程，但总体上点E的沉降量552023年第9期级要高于点H。因此结合图3 的时序形变结果，可以发现沿开采方向右侧的地表形变值要大于左侧区域，即112560E112570EN人N.0.S1.96112560E煤炭与化工始采线附近的沉陷盆地实际上有朝西北方向发育的趋势，后期需要对东李高村附近地表进行重点监测。0wu/鲁250100-150N.0.ST-200(b)21/101 21/2272/0324 2/061922/0914 2/1/109830N.0.t1wu/鲁305001.000112570E第46 卷AB时间098-60-90E时间-369Fig.4 The location of feature points and their time-series cumulative subsidence values4结 论（1）在采动作用下，S1306工作面周边地表自2 0 2 2 年1月16 日开始出现下沉，之后在始采线附近迅速产生了一个沉陷盆地，且沉陷盆地有朝西北方向发育的趋势。（2）目前采动影响仅波及至工作面的中部区域，与实际开采到达位置基本吻合，说明SBAS-InSAR技术在矿区形变监测中具有良好的适用性。（3）除了利用时序InSAR技术获取不同时段的地表沉陷情况，仍可以考虑利用UAV激光雷达扫描技术或建立地面移动观测站来对局部重点区域进行详细监测，真正发挥“空天地”协同监测技术在开采沉陷领域的重要作用。参考文献：1何国清，杨伦，凌娣，等.矿山开采沉陷学M.徐州：中国矿业大学出版社，1991.2 朱建军，李志伟，胡俊.InSAR变形监测方法与研究进展J.测绘学报，2 0 17，46（10)：17 17-17 3 3.3朱建军，杨泽发，李志伟.InSAR矿区地表三维形变监测与年平均形变速率/(mm*yl)图4特征点位及其时序累积沉降值4李达，邓喀中，高晓雄，等.基于SBAS-InSAR的矿区地表沉降监测与分析J.武汉大学学报（信息科学版），2018,43(10):1531 1537.5栾元重，梁耀东，纪赵磊，等.基于SBAS-InSAR技术采动地表沉降监测与分析J.煤炭科学与技术，2 0 2 0，48（10）：198 204.6柴华彬，胡吉彪，耿思佳.融合实测数据的地表沉降SBAS-InSAR监测方法J.煤炭学报，2 0 2 1，46（S1）：17-2 4.7张香凝，贺黎明，刘翠芝，等.基于SBAS-InSAR技术的煤矿开采沉降监测与分析J1遥感技术与应用，2 0 2 2，3 7（4）：1 021-1 028.8欧阳伦曦，李新情，惠风鸣，等.哨兵卫星Sentinel-1A数据特性及应用潜力分析J.极地研究，2 0 17，2 9（2）：2 8 6-295.9谭志祥，杨嘉威，邓喀中.基于SBAS-InSAR的矿区全盆地开采沉陷求参方法研究J.煤炭科学技术，2 0 2 1，49（1）：312-318.10 李珊珊，李志伟，胡俊，等.SBAS-InSAR技术监测青藏高原季节性冻土形变J.地球物理学报，2 0 13，56（5）：1 476 1 486.11 曹发伟，廖维谷.SBAS技术在矿区地面沉降监测中的应用J.测绘通报，2 0 2 1（3）：156-158，16 3.12 周亚涛.襄垣县采煤诱发地质灾害的工程地质研究D.太原：太原理工大学，2 0 0 2.175120(c)21/1001 21/2/2722/03/24 2/06/1922/09/1422/1/10预计研究进展J.测绘学报，2 0 19，48（2）：13 5144.56