融入坐标—时间函数(CT-PIM)的矿区时序InSAR形变预计——以淮安盐矿为例  

Incorporation of Coordinate−Time Function(CT−PIM)time-series InSAR deformation prediction for salt mining areas:Case study of the Huaian Salt Mine

作　　者：张腾飞 邢学敏[1,2] 彭葳 朱珺 刘祥彬 葛家旺 雷敏超 ZHANG Tengfei;XING Xuemin;PENG Wei;ZHU Jun;LIU Xiangbin;GE Jiawang;LEI Minchao(School of Traffic and Transportation Engineering,Changsha University of Science&Technology,Changsha 410114,China;Institute of Radar Remote Sensing Applications for Traffic Surveying and Mapping,Changsha University of Science&Technology,Changsha 410114,China)

机构地区：[1]长沙理工大学交通运输工程学院,长沙410114 [2]长沙理工大学交通测绘雷达遥感应用研究所,长沙410114

出　　处：《遥感学报》2024年第6期1615-1631,共17页NATIONAL REMOTE SENSING BULLETIN

基　　金：国家自然科学基金(编号:42074033,51878078,41701536,41904003,61701047);湖南省自然科学基金青年基金(编号:2017JJ3322,2019JJ50639,2020JJ5571);湖南省教育厅重点项目(编号:18A148,19C0042);洞庭湖生态环境遥感监测湖南省重点实验室开放项目(编号:2021-011);中国水利水电第八工程局有限公司科研项目(编号:2023060);长沙市杰出创新青年培养计划项目(编号:kq2209011)。

摘　　要：对矿区开展长期形变监测和后续形变预计对于预防矿山开采所引发的潜在安全问题和矿山生态环境保护具有十分重要的意义。利用合成孔径雷达干涉测量技术InSAR(Interferometric Synthetic Aperture Radar)与概率积分法PIM(Probability Integral Method)结合是较为有效的矿区形变预计方法。然而,目前矿区时序InSAR形变建模环节中大多数采用纯经验数学模型,未顾及开采沉陷机制,严重影响获取的形变观测值的精度;不准确的InSAR形变观测值将误差传递至后续形变预计,进而导致误差累积。本文将坐标—时间函数CT(Coordinate Time)引入时序InSAR建模环节,构建了坐标—时间函数预计模型(CT-PIM)以取代传统纯经验数学模型;再基于InSAR相位观测值求解未知参数,将构建的CT-PIM直接用于矿区形变预计,避免了误差的传递,以改善形变预计精度。分别开展了模拟实验和真实实验对提出预计方法进行验证。模拟实验显示:模型预计获取的时序形变与模拟真值的均方根误差(RMSE)为±4.6 mm。真实实验利用覆盖淮安市一盐矿区的35景Sentinel-1A SAR卫星影像开展,获取了研究区域2019年3月30日至2019年7月28日的预计结果,结果显示测区最大预计沉降量为152 mm。本文所提方法相比传统多速率模型,建模精度提升了38.2%;相比传统静态概率积分预计方法获取的形变结果,精度提升了39.1%。本文方法为预防盐矿区采矿引起的地质灾害问题提供了有力的工具,为盐矿区生产安全和生态环境保护提供参考。Long-term monitoring and the subsequential prediction of deformation for salt mining areas is essential to the safety prevention and environmental protection of mining areas.The combination of the interferometric synthetic aperture radar(InSAR)technique with the Probability Integral Method(PIM)has proven to be powerful in predicting the deformation of mining areas.However,single multitemporal InSAR(MT-InSAR)is limited because it can only obtain the deformation sequences during SAR acquisition dates,and the subsequent future displacement beyond the span of the SAR observations cannot be acquired.In addition,traditional mathematical empirical models are mostly used in the time-series modeling of mining areas,ignoring the underground mining mechanisms,which seriously affect the accuracy of the observations.Inaccurate InSAR deformation monitoring results transmit errors to forward predicted subsidence,which may induce considerable errors.In this study,the Coordinate-Time(CT)function is introduced into time-series InSAR deformation modeling,and a CT function prediction model(CT-PIM),which can well describe the dynamic evolution disciplines of the underground mining subsidence in InSAR deformation modeling,is constructed to replace the traditional mathematical empirical models.The unknown CT-PIM parameters can be estimated directly via InSAR time-series phase observations,and the constructed CT-PIM is directly used in the deformation prediction of the mining area,which can avoid the error propagation from the InSAR-generated deformations and improve deformation prediction accuracy.The new approach is tested by simulation and real data experiments.The simulation results show that the root mean square error between the time-series deformation prediction of the model and the simulated true value is estimated to be±4.6 mm,which implies that the proposed method is of promising accuracy.The real experiment was conducted using a total of 35 Sentinel-1A SAR images covering the salt mining area in Huaian City,and the deformati

关 键 词：遥感 INSAR 矿区 坐标—时间函数 地表沉陷 形变预计 

分 类 号：P237[天文地球—摄影测量与遥感] P2[天文地球—测绘科学与技术]