沙永莲，王晓文，刘国祥，张 瑞，张 波

(1.西南交通大学地球科学与环境工程学院,成都 611756；2.北京理工大学重庆创新中心,重庆 401120;3.高速铁路运营安全空间信息技术国家地方联合工程实验室,成都 611756)

0 引言

煤炭是我国能源结构的主体，对煤炭资源开采活动进行合理的规划和管理对我国经济建设和生态环境保障具有重要意义[1]。我国煤炭资源多是通过地下开采，而这种开采模式往往会破坏地层内部原有的力学平衡状态，导致采空区上方地表发生移动和变形[2]。煤田矿区地表不均匀变形可能导致地面形成裂缝或者塌陷，破坏原有地质环境条件并引发地质灾害，从而严重影响地表相关建筑物等基础设施的安全使用，给矿区当地经济可持续发展造成隐患[3]。

针对矿区不均匀地表形变的时空演化特征进行监测，可为当地地表沉陷隐患治理，安全生产保障和当地生态环境保护等提供关键支撑信息。目前，水准测量、全球导航卫星系统(GNSS)和合成孔径雷达干涉(InSAR)等多种大地测量方法已在矿区地表沉陷监测中得到广泛应用[4-5]。相比于传统大地测量手段，InSAR具有全天候观测、覆盖范围广和观测效率高等优势[6-7]。常规基于两景SAR影像进行干涉处理测量地表形变的两轨法DInSAR (Differential InSAR)，由于无法对形变观测模型中的误差项(如地形误差、轨道误差和大气延迟误差等)进行校正，常被用于瞬时、大幅度地表形变(如地震和火山)测量[8-9]。近十年来，众多学者基于多景SAR影像提出多种时序InSAR分析方法以克服常规DInSAR相位失相干、大气延迟等影响，并实现长时间序列地表形变监测[10-11]。……