基于SBAS-InSAR技术的老采空区注浆充填地表变形时空分布特征分析  

Analysis of spatial-temporal distribution characteristics of surface deformation during grouting filling in goaf based on SBAS-InSAR technology

作　　者：李幸丽 戴华阳[2] 方军 张豪磊 LI Xingli;DAI Huayang;FANG Jun;ZHANG Haolei(CCTEG Ecological Environment Technology Co.,Ltd.,Tianjin 300450,China;College of Geoscience and Surveying Engineering,China University of Mining and Technology(Beijing),Beijing 100083,China)

机构地区：[1]中煤科工生态环境科技有限公司,天津300450 [2]中国矿业大学(北京)地球科学与测绘工程学院,北京100083

出　　处：《中国矿业》2024年第11期86-94,共9页China Mining Magazine

基　　金：中煤科工生态环境科技有限公司科技创新基金项目资助(编号:0206KGST0027);国家重点研发计划项目资助(编号:2023YFF1306003)。

摘　　要：InSAR技术已广泛应用于矿区地表沉降监测,但老采空区注浆治理区域地表变形时空分布特征研究较少。为获取老采空区注浆充填地表变形时空特征,本文提出了一种基于SBAS-InSAR技术的数据处理方法,该方法联合使用了D-InSAR技术、PS-InSAR技术、SBAS-InSAR技术等多种技术,以SBAS图像对选取策略选择时空基线较优的图像对进行组合,以D-InSAR技术为基础,生成干涉图集,通过PS方法选择稳定的散射相干目标点,对相干目标点相位信息进行解算,获取了相干目标点地表沉降速率及时间序列等信息。通过与部分水准数据的对比,表明该方法与传统水准测量具有一致性。基于39期RadarSAT-2数据及40期Sentinel数据,获取了研究区地表在注浆前、注浆中及注浆后的变形特征:研究区在注浆前(2015年7月—2016年12月)区域最大沉降速率达-98mm/a,最大抬升速率达134mm/a;在注浆中(2019年8月—2022年9月),区域最大抬升速率达18.7mm/a,最大沉降速率为-11.2mm/a;在注浆后(2023年1月—2023年12月)最大抬升速率达16.0mm/a,最大沉降速率达-10.0mm/a。区域地表变形速率呈现下降趋势。通过区域沉降速率、时序累计沉降值及典型相干点分析,获得了研究区从注浆前、注浆中到注浆后的地表变形分布特征,揭示了研究区地表变形从注浆前开采引起地表下沉及残余变形到注浆中地表缓慢下沉并抬升到注浆后地表变形趋于稳定的变化规律。InSAR technology has been widely used in monitoring of surface subsidence in mining areas,but it is poor to research on the spatial-temporal distribution characteristics of surface deformation in grouting areas.In order to obtain the spatial-temporal distribution characteristics of surface deformation during grouting filling in goaf,this paper proposes a data processing method based on SBAS-InSAR technology,which jointly uses various technologies such as D-InSAR technology,PSInSAR technology,SBAS-InSAR technology,etc.,and selects the image pairs with better spatialtemporal baselines to be combined by the SBAS image pair selection strategy.The D-InSAR technology is used to generate the interferogram images,select the stable scattering coherent target points by the PS method,solve the phase information of the coherent target points,and obtain the information of the surface subsidence rate and time series of the coherent target points.The comparison with some leveling data shows that the method is consistent with traditional leveling survey.Based on the 39-phase RadarSAT-2 data and 40-phase Sentinel data,the deformation characteristics of surface during grouting filling in goaf have been obtained:the max deformation rate before grouting filling(from July 2015 to December 2016)reaches−98 mm/a,and the max uplift rate reaches 134 mm/a;the max uplift rate during grouting filling(from August 2019 to September 2022)reaches 18.7 mm/a,and the max subsidence rate is−11.2 mm/a;after grouting filling(from January 2023 to December 2023),the max uplift rate reaches 16.0 mm/a,and the max subsidence rate reaches−10.0 mm/a.The regional surface deformation rate shows a decreasing trend.Through the analysis of the average deformation rate,the series accumulated subsidence and typical coherent points,obtains the distribution characteristics of surface deformation in the research area from before grouting filling,during grouting filling to after grouting filling,the rule of the surface deformation in the study area is discovered

关 键 词：注浆充填 SBAS-InSAR技术 地表变形 老采空区 开采沉陷 

分 类 号：TD327[矿业工程—矿井建设]