基于时序InSAR的即墨市区地表形变监测

地面沉降是一种对地面及地下基础设施造成安全隐患,对经济可持续发展和环境保护产生破坏影响的地质灾害现象。本文使用2017年5月至2018年5月16景Sentinel-1A卫星SAR影像,根据D-InSAR的初步形变监测结果将即墨城区内沉降明显的区域作为研究区,基于PS和SBAS两种时序InSAR方法对该区域进行地面沉降监测,获得的沉降分布和形变时序结果吻合。地面沉降的分布与新建高层建筑区吻合,地面形变趋势与区域降水量和地下水水位变化有较高的相关性。研究结果有助于了解即墨城区的地表沉降状况以及原因,为地面沉降综合治理和地下水水资源开发利用提供参考依据。     

基于SAR干涉点目标分析技术的城市地表形变监测

通过深入研究干涉点目标的相位模型,提出基于空间搜索的邻近点目标干涉相位差解缠方法,用以计算点目标的地形残差和线性形变,以及分离点目标大气延迟相位和非线性形变相位的时空域滤波方法,解决干涉点目标分析中的关键问题.最后,以苏州地区地表沉降监测为应用试验,利用形成的SAR干涉点目标形变信息提取技术,获取苏州市区1992-2002年间的地表沉降信息.研究结果与已有文献记录保持比较好的一致性,证明SAR干涉点目标技术完全可以发展成为应用于城市地表形变监测的实用化技术.

Abstract：

Interferometric point target analysis (IPTA) is one of the latest developments in radar interferometric processing. It is achieved by analysis of the interferometric phases of some individual point targets, which are discrete and present temporarily stable backscattering characteristics, in long temporal series of interferometric SAR images. This paper analyzes the interferometric phase model of point targets, and then addresses two key issues within IPTA process. Firstly, a spatial searching method is proposed to unwrap the interferometric phase difference between two neighboring point targets. The height residual error and linear deformation rate of each point target can then be calculated, when a global reference point with known height correction and deformation history is chosen. Secondly, a spatial-temporal filtering scheme is proposed to further separate the atmosphere phase and nonlinear deformation phase from the residual interferometric phase. Finally, an experiment of the developed IPTA methodology is conducted over Suzhou urban area. Totally 38 ERS-1/2 SAR scenes are analyzed, and the deformation information over 3 546 point targets in the time span of 1992-2002 are generated. The IPTA-derivecl deformation shows very good agreement with the published result, which demonstrates that the IPTA technique can be developed into an operational tool to map the ground subsidence over urban area.     

Assessment of Land Surface Subsidence Due to Underground Metal Mining Using Integrated Spaceborne Repeat-Pass Differential Interferometric Synthetic Aperture Radar (DInSAR) Technique and Ground Based Observations

Differential interferometric synthetic aperture radar (DInSAR) is a novel remote sensing technique to measure earth surface deformation. It is capable of obtaining dense information related to the deformation of a large area efficiently, economically and effectively. Therefore, DInSAR is a promising technology for monitoring the earth surface deformation related to some natural hazardous events, such as earthquake, volcano eruption, land subsidence, landslide. In present study, Conventional DInSAR technique have been applied to a mineral rich zone, coming under the Khetri copper belt, a part of Northern Aravali range of hillocks in India, predominant with mining activities since late 1960’s to address the possibility of deformation phenomena due to hard rock underground metal mining. Four interferometric SAR data sets of Radarsat-2 was used for the study area to address the subsidence/uplift phenomena. Further, results obtained from conventional DInSAR technique using Radarsat-2 data sets compared with results obtained from ground based observation technique for its validity. In both the techniques, deformation results obtained in terms of average subsidence rate in mm (quarterly basis) of points under study within mining zone of Mine-A has well agreed to each other. Further, it has been observed that average subsidence rate in mm (quarterly basis) obtained from space based observation and ground based observation are 5.6 and 6.67, respectively over the points under study in mining zone of Mine-A.     

Deriving time-series three-dimensional displacements of mining areas from a single-geometry InSAR dataset

This paper presents a method for deriving time-series three-dimensional (3-D) displacements of mining areas from a single-geometry interferometric synthetic aperture radar (InSAR) dataset (hereafter referred to as the SGI-based method). This is mainly aimed at overcoming the limitation of the traditional multi-temporal InSAR techniques that require SAR data from at least three significantly different imaging geometries to fully retrieve time-series 3-D displacements of mining areas. The SGI-based method first generates the multi-temporal observations of the mining-induced vertical subsidence from the single-geometry InSAR data, using a previously developed method for retrieving 3-D mining-related displacements from a single InSAR pair (SIP). The weighted least-squares solutions of the time series of vertical subsidence are estimated from these generated multi-temporal observations of vertical subsidence. Finally, the time series of horizontal motions in the east and north directions are estimated using the proportional relationship between the horizontal motion and the subsidence gradient of the mining area, on the basis of the SGI-derived time series of vertical subsidence. Seven ascending ALOS PALSAR images from the Datong mining area of China were used to test the proposed SGI-based method. The results suggest that the SGI-based method is effective. The SGI-based method not only extends the SIP-based method to time-series 3-D displacement retrieval from a single-geometry InSAR dataset, but also limits the uncertainty propagation from InSAR-derived deformation to the estimated 3-D displacements.     

百年煤城地表沉降融合PS/SBAS InSAR监测——以徐州市为例

融合PS/SBAS算法成为当前InSAR技术的研究热点,本文通过融合PS/SBAS的优势,利用128景ALOS PALSAR影像和24景Radarsat-2影像获取了百年煤城徐州地区2007-2015年的地表形变场及变化趋势。试验结果表明：①2007-2011年,徐州主要有4个沉降区,分布在沛县、丰县、铜山区和贾汪区。②2012-2015年,沛县、丰县和铜山区地表沉降范围有所扩大;贾汪区地表沉降范围和速率明显减小,说明近年来贾汪区重点治理采矿塌陷工作已初步见效;另外,睢宁县县城出现明显沉降。③已有文献对沛县大屯中心区的沉降监测结果验证了试验结果是可靠的。④首次提取了徐州近8年间地表沉降的信息,可为该地区开展地表沉降调查监测与防治工作提供参考;试验方法可为全国开展地表沉降调查监测工作提供参考。     

京津冀地区1992—2014年三阶段地面沉降InSAR监测

京津冀地区是我国地面沉降灾害发育最严重的地区,几乎每年都造成巨大的经济损失。本文提出了改进的时间序列InSAR技术——多主影像相干目标小基线干涉技术（MCTSB-InSAR）,利用4颗卫星摄取的3个时段的时间序列SAR影像：ERS-1/2 SAR（1992—2000年）、ENVISAT ASAR（2003—2010年）、RADARSAT-2（2012—2014年）,获取了京津冀地区1992—2014年间3个时段的地面沉降信息。经与京津两地120个以上水准测量数据进行比较,3个时段的地表监测结果的精度分别为8.7、4.7、5.4mm/a。分析了北京和天津两市22年间地面沉降的时空变化特征,其中北京市地面沉降呈不断加重趋势;天津市地面沉降在1992—2010年间发展迅猛,在2010年以后有所减缓。同时,本文也表明MCTSB-InSAR技术是有效可靠的,在大区域地面沉降监测中具有推广应用价值。     

Radar interferometry time series analysis of Mashhad subsidence

A large agricultural area located in 20 km north of the city of Mashhad in the north-east of Iran is subject to land subsidence. The subsidence rate was achieved in a couple of sparse points by precise leveling between 1995 and 2005, and continuous GPS measurements obtained from 2005 to 2006. In order to study the temporal behavior of the deformation in high spatial resolution, the small baseline subset (SBAS) algorithm was used to generate the interferometric SAR time series analysis. Time series analysis was performed using 19 interferograms calculated from 12 ENVISAT ASAR data spanning between 2003 and 2006. The time series results exhibited that the area is subsiding continuously without a significant seasonal effect. Mean LOS deformation velocity map obtained from time series analysis demonstrated a considerable subsidence rate up to 24 (cm/yr). In order to evaluate the time series analysis results, continuous GPS measurements as a geodetic approach were applied. The comparisons showed a great agreement between interferometry results and geodetic technique. Moreover, the information of various piezometric wells distributed in the area corresponding to 1995 to 2005 showed a significant decline in water table up to 20 meters. The correlation between the piezometric information and the surface deformation at well’s locations showed that the subsidence occurrence in Mashhad is due to the excess groundwater withdrawal.     

基于差分雷达干涉测量的矿区地面沉降监测研究

收集了1992年12月至1998年6月JERS-1L波段雷达数据,利用差分雷达干涉测量技术进行地面沉降监测研究。通过"二轨法"进行差分处理,选择其中的6景差分雷达干涉图进行分析。在研究区内发现了4个沉降区域,并按照时间序列通过对地面形变图进行分析,获取其中两个沉降区最大的垂直沉降量。     

联合InSAR和SSA的膨胀土边坡形变特征分析——以南水北调工程为例

膨胀土边坡具有吸水剧烈膨胀、失水剧烈收缩及反复胀缩形变的特性,工程中与外界接触极易引发膨胀土滑坡。南水北调中线工程是一项缓解我国华北平原水资源短缺、优化水资源配置的国家重大战略工程。中线总干渠约有1/3段穿越膨胀土区域,增加了南水北调工程沿线边坡失稳的可能性,亟须开展膨胀土边坡形变监测和特征分析。基于此,本文收集了112景Sentinel-1A升轨SAR数据,采用时间序列InSAR技术获取了中线某渠道沿线2017年3月12日—2021年4月8日的时间序列形变,并基于奇异谱分析(SSA)对时间序列形变进行分解提取膨胀土边坡形变特征。研究结果表明:干渠挖方膨胀土边坡由于卸荷回弹作用表现为抬升形变,最大年抬升速率约为18 mm/a;填充膨胀土边坡由于土体固结作用表现为下沉形变,最大年下沉速率约为15 mm/a;膨胀土形变量大小与填挖深度正相关,挖方越深抬升形变越大,填充越高下沉形变越大;降雨和温度是影响干渠膨胀土边坡形变的主要外部因素,其造成的形变会时滞2~3个月。本文研究将为南水北调工程膨胀土边坡稳定性评估提供科学依据。     

PALSAR和ASARPSI显著地表沉降探测与分析

永久散射体雷达干涉(PSI)已被广泛应用于城市缓慢沉降监测,但针对显著性地表沉降监测及PSI适用性的探讨较少。本文选取天津市与廊坊市城郊结合部的典型工业化城镇为研究区,分别以2007至2010年间的23幅L波段(23.6 cm)ALOS PALSAR影像和2007至2009年间的23幅C波段(5.6 cm)ENVISAT ASAR影像为数据源,使用基于相位空间相关性分析的PSI方法进行时空沉降提取,并对两个平台的结果进行对比分析和交叉验证。在此基础上,结合长短波SAR对沉降敏感性不同的特征探讨该PSI方法应用于显著沉降解算的可靠性,并结合地质条件及地下水开采信息对实验区沉降的时空分布进行深入分析和解释。研究表明,两平台的沉降结果十分吻合,在沉降漏斗区二者均方根误差(RMSE)为6.5 mm/a;实验区地表沉降呈现出显著的非均匀性,最大沉降速率超过200 mm/a,造成严重沉降的主要原因为地下水的大量开采;证实了PSI方法在显著地表沉降监测方面的有效性和可靠性。     

Sentinel-1A TS-DInSAR京津冀地区沉降监测与分析

近年来,京津冀地区的不均匀地表沉降日趋严重,对公路、铁路等基础设施安全造成严重威胁,已引起国内外广泛关注。合成孔径雷达时序差分干涉TS-DIn SAR(Time Series Differential Interferometric Synthetic Aperture Radar)作为一种高效的广域形变测量手段,在地面沉降调查监测中已被广泛应用,但如何针对高现势性的影像数据高精准提取大范围区域地表形变是当前深化的热点。本文利用欧洲太空局发布的Sentinel-1A新型数据源,针对TOPS(Terrain Observation by Progressive Scans)模式影像间存在多普勒中心不一致问题,借助外部高精度POD(Precise Orbit Determination)轨道和DSM(Digital Surface Model)数据进行频移滤波迭代配准,并针对Sentinel-1A数据特征集成优化了基于点目标时序分析的大区域地表形变监测方法,以2015年—2016年期间的29景影像为实验数据开展了京津地区沉降监测研究,提取了北京东、廊坊及天津西等地区的沉降结果,并结合区域内典型区域人口密度、产业分布、地表覆盖和线路剖面等信息深入分析了沉降的时空分布特征和成因。结果表明,Sentinel-1A时序干涉结果在大范围地表沉降调查监测上具有可靠的应用精度。     

时间序列InSAR技术中的形变模型研究

时间序列InSAR技术,包括PS-InSAR技术和小基线集InSAR技术, 能有效克服传统D-InSAR的失相干限制, 逐步成为形变测量中的实用化技术。但是, 现有时间序列InSAR技术主要采用线性函数对真实形变进行模拟, 在真实形变呈现强烈的非线性时, 这种处理将不能得到正确的形变结果。本文就时间序列InSAR的形变模型问题展开研究, 首先从干涉相位模型解算的方法入手, 深入分析了线性形变模型的不足, 当干涉点目标的密度不够并且真实形变的非线性较强时, 干涉相位方程的解将会发散。根据魏尔斯特拉斯逼近定理, 提出以高阶多项式取代线性形变模型, 并给出了基于多项式形变模型的干涉相位方程解算方法。 利用太原市2003-2009年的23景ENVISAT ASAR影像, 分别采用线性形变模型和三次多项式形变模型, 利用小基线集技术进行了形变反演。将这两种方法得到的结果分别与水准测量结果进行了比较。结果表明, 采用多项式形变模型不仅能取得更高的形变测量精度,而且能提高点目标的密度。由于高阶多项式总能比低阶多项式更准确地拟合连续函数, 因此本文提出的多项式形变模型在时间序列InSAR形变监测中具有广泛的应用价值。     

改进的PS探测方法及其应用

提出一种改进的永久散射体PS(Persistent Scatterer)探测方法,首先基于时序SAR影像的振幅信息,利用振幅和振幅离差指数双阈值法获取PS候选点,然后根据干涉相位空间相关性,对所有候选点进行相位稳定性分析,筛选出既满足雷达反射强度高、又符合散射稳定性要求的真实PS点,最后基于所选PS点提取形变信息。选取上海市局部区域作为实验区,以16幅TerraSAR-X(TSX)影像为数据源,进行PS探测和形变信息提取。结果表明,相比于已有的PS探测方法,本文方法能够更准确、可靠地识别出时序SAR影像中的PS点。形变结果显示,研究区域内最大沉降速率为-51mm/a,平均沉降速率为-15.8mm/a,以水准实测数据进行检核,验证了基于改进的PS探测方法进行地表形变信息提取结果的可靠性。     

Ground deformation associated with post-mining activity at the French–German border revealed by novel InSAR time series method

We present a novel methodology for integration of multiple InSAR data sets for computation of two dimensional time series of ground deformation. The proposed approach allows combination of SAR data acquired with different acquisition parameters, temporal and spatial sampling and resolution, wavelength and polarization. Produced time series have combined coverage, improved temporal resolution and lower noise level. We apply this methodology for mapping coal mining related ground subsidence and uplift in the Greater Region of Luxembourg along the French–German border. For this we processed 167 Synthetic Aperture Radar ERS-1/2 and ENVISAT images acquired between 1995 and 2009 from one ascending (track 29) and one descending (track 337) tracks and created over five hundred interferograms that were used for time series analysis. Derived vertical and east–west linear deformation rates show with remarkable precision a region of localized ground deformation located above and caused by mining and post-mining activities. Time series of ground deformation display temporal variability: reversal from subsidence to uplift and acceleration of subsidence in the vertical component, and horizontal motion toward the center of the subsidence on the east–west component. InSAR results are validated by leveling measurements collected by the French Geological Survey (BRGM) during 2006–2008. We determined that deformation rate changes are mainly caused by water level variations in the mines. Due to higher temporal and spatial resolution the proposed space-borne method detected a larger number of subsidence and uplift areas in comparison to leveling measurements restricted to annual monitoring of benchmark points along roads. We also identified one deformation region that is not precisely located above the mining sites. Comparison of InSAR measurements with the water levels measured in the mining pits suggest that part of the water that filled the galleries after termination of the dewatering systems may come from this region. Providing that enough SAR data is available, this method opens new opportunities for detecting and locating man-made and natural ground deformation signals with high temporal resolution and precision.     

Monitoring of urban subsidence with SAR interferometric point target analysis: A case study in Suzhou,China

Ground subsidence, mainly caused by over exploitation of groundwater and other underground resources, such as oil, gas and coal, occurs in many cities in China. The annual direct loss associated with subsidence across the country is estimated to exceed 100 million US dollar. Interferometric SAR (InSAR) is a powerful tool to map ground deformation at an unprecedented level of spatial detail. It has been widely used to investigate the deformation resulting from earthquakes, volcanoes and subsidence. Repeat-pass InSAR, however, may fail due to impacts of spatial decorrelation, temporal decorrelation and heterogeneous refractivity of atmosphere. In urban areas, a large amount of natural stable radar reflectors exists, such as buildings and engineering structures, at which radar signals can remain coherent during a long time interval. Interferometric point target analysis (IPTA) technique, also known as persistent scatterers (PS) InSAR is based on these reflectors. It overcomes the shortfalls in conventional InSAR. This paper presents a procedure for urban subsidence monitoring with IPTA. Calculation of linear deformation rate and height residual, and the non-linear deformation estimate, respectively, are discussed in detail. Especially, the former is highlighted by a novel and easily implemented 2-dimensional spatial search algorithm. Practically useful solutions that can significantly improve the robustness of IPTA, are recommended. Finally, the proposed procedure is applied to mapping the ground subsidence in Suzhou city, Jiangsu province, China. Thirty-four ERS-1/2 SAR scenes are analyzed, and the deformation information over 38,881 point targets between 1992 and 2000 are generated. The IPTA-derived deformation estimates correspond well with leveling measurements, demonstrating the potential of the proposed subsidence monitoring procedure based on IPTA technique. Two shortcomings of the IPTA-based procedure, e.g., the requirement of large number of SAR images and assumed linear plus non-linear deformation model, are discussed as the topics of further research.     

Spatial-Correlation Based Persistent Scatterer Interferometric Study for Ground Deformation

Interferometric Synthetic Aperture Radar (InSAR), nowadays, is a precise technique for monitoring and detecting ground deformation at a millimetric level over large areas using multi-temporal SAR images. Persistent Scatterer Interferometric SAR (PSInSAR), an advanced version of InSAR, is an effective tool for measuring ground deformation using temporally stable reference points or persistent scatterers. We have applied both PSInSAR and Small Baseline Subset (SBAS) methods, based on the spatial correlation of interferometric phase, to estimate the ground deformation and time-series analysis. In this study, we select Las Vegas, Nevada, USA as our test area to detect the ground deformation along satellite line-of-sight (LOS) during November 1992–September 2000 using 44 C-band SAR images of the European Remote Sensing (ERS-1 and ERS-2) satellites. We observe the ground displacement rate of Las Vegas is in the range of ?19 to 8 mm/year in the same period. We also cross-compare PSInSAR and SBAS using mean LOS velocity and time-series. The comparison shows a correlation coefficient of 0.9467 in the case of mean LOS velocity. Along this study, we validate the ground deformation results from the satellite with the ground water depth of Las Vegas using time-series analysis, and the InSAR measurements show similar patterns with ground water data.     

Measurement of subsidence in the Yangbajing geothermal fields,Tibet, from TerraSAR-X InSAR time series analysis

Yangbajain contains the largest geothermal energy power station in China. Geothermal explorations in Yangbajain first started in 1976, and two plants were subsequently built in 1981 and 1986. A large amount of geothermal fluids have been extracted since then, leading to considerable surface subsidence around the geothermal fields. In this paper, InSAR time series analysis is applied to map the subsidence of the Yangbajain geothermal fields during the period from December 2011 to November 2012 using 16 senses of TerraSAR-X stripmap SAR images. In the case of the TerraSAR-X data, most orbital fringes were removed using precise orbits during the interferometric processing. However, residual orbital ramps remain in some interferograms due to the uncertainties in the TerraSAR-X orbits. To remove the residual orbital ramps, we estimated a best-fit ‘twisted plane’ for each epoch interferogram using quadratic polynomial models based on a network approach. This method removes most of the long-wavelength signals, including orbit ramps and atmospheric effects. The vertically strati?ed component (Topography Correlated Atmospheric Delay, TCAD) was also removed using a network approach. If the influence of seasonal frozen ground (SFG) is not taken into consideration, our results show that the subsidence rate around power plant I (the south plant) is approximately 20?mm/yr with a peak of 30?mm/yr. The subsidence rate around power plant II (the north plant) is approximately 10?mm/yr, when accounting for the influence of SFG on the power plant and its surrounding ground surface. Our results show that ground motion is caused by seasonal frozen ground and is strongly related to the temperature change.     

Interferometric SAR Time Series Analysis for Ground Subsidence of the Abandoned Mining Area in North Peixian Using Sentinel-1A TOPS Data

The North Peixian mining area of China has rich coal resources, with total proven reserves of 2.37 billion tons. However, the underground coal mining activities have resulted in ground collapse, which has caused serious harm to the environment and threatened the lives and properties of local residents. In this study, 12 Sentinel-1A terrain observation by progressive scans (TOPS) mode acquisitions between 30 July 2015 and 13 May 2016 over the abandoned mining area in North Peixian were analyzed using the interferometric synthetic aperture radar (InSAR) time series method to detect the ground subsidence, with the maximum ground subsidence reaching 83 mm/a and an average value of about 12.7 mm/a. The subsidence results derived from the Sentinel-1A TOPS mode dataset were proven to be effective in investigating and monitoring the ground subsidence in the North Peixian mining area. Compared to the rapid deformation during the ongoing period of mining excavation, the ground subsides slowly in abandoned mining areas and shows a linear relationship with time over a relatively long period of time. Spatial correlation between the subsidence distribution and land cover was found, in that the magnitude of the subsidence in urban areas was smaller than that in rural areas, which is associated with the controlled coal mining activities under buildings, railways, and water bodies. The results demonstrate that Sentinel-1A TOPS SAR images can be used to effectively and accurately detect and monitor ground subsidence in a mining area, which is critically important when investigating land subsidence in a large-scale mining area.     

基于PALSAR数据的青藏高原冻土形变检测方法研究

季节性冻胀和融沉导致的地面形变是青藏高原冻土区建设施工与维护的主要问题。对冻融造成的形变进行有效监测是青藏铁 路建设与维护的前提。差分干涉测量技术是地表形变监测的重要手段之一,PALSAR（L波段的合成孔径雷达）数据在非城市区域具 有较高的相关性,适合青藏高原冻土区的地表形变监测。本文选用4景覆盖研究区域的PALSAR数据,研究利用该数据进行冻土形变 检测的方法,并对其检测结果进行了分析。结果表明,该方法与水准测量方法有较好的一致性。     

From ERS-1/2 to Sentinel-1: two decades of subsidence monitored through A-DInSAR techniques in the Ravenna area (Italy)

Land subsidence due to underground resources exploitation is a well-known problem that affects many cities in the world, especially the ones located along the coastal areas where the combined effect of subsidence and sea level rise increases the flooding risk. In this study, 25 years of land subsidence affecting the Municipality of Ravenna (Italy) are monitored using Advanced Differential Interferometric Synthetic Aperture Radar (A-DInSAR) techniques. In particular, the exploitation of the new Sentinel-1A SAR data allowed us to extend the monitoring period till 2016, giving a better understanding of the temporal evolution of the phenomenon in the area. Two statistical approaches are applied to fully exploit the informative potential of the A-DInSAR results in a fast and systematic way. Thanks to the applied analyses, we described the behavior of the subsidence during the monitored period along with the relationship between the occurrence of the displacement and its main driving factors.