基于SBAS-InSAR技术的煤矿开采沉陷监测与分析

地下煤炭资源大量开采导致的地表形变,引发严重的安全和环境隐患,雷达干涉测量技术是高精度、大范围地表形变监测的重要手段之一。以辽宁省沈阳市沈北新区蒲河煤矿为例,采用SBAS-InSAR技术探测2018—2019年矿区地表形变结果,获取了采煤引起地表形变的时空分布特征,结合采场所在区域的地质条件和变形诱发因素,利用数值模拟技术对观测形变结果进行模拟分析,进而讨论了蒲河煤矿地面沉降在时间和空间上的变形规律和机制。InSAR形变监测结果显示,开采区域内存在两处沉降漏斗,且数值模拟结果与InSAR形变观测值分布规律一致,反演结果接近实际情况,可为相关部门制定地面沉降防治措施提供科学依据。     

Fuzzy probability measures (FPM) based non-symmetric membership function: Engineering examples of ground subsidence due to underground mining

This paper introduces a new FPM method for prediction of surface subsidence due to inclined coal seam mining. Based on the non-symmetric membership function and the definition of the fuzzy probability measure, the mathematical model for the two-dimensional problem is developed and applied to the analysis of ground subsidence due to underground mining of inclined coal seam. The results of calculation with the new model indicate that the predicted value has good agreement with the measured data.     

Geomorphological changes associated with underground coal mining in the Fushun area,northeast China revealed by multitemporal satellite remote sensing data

Fushun is a famous coal-mining city in northeastern China with more than 100 years of history. Long-term underground coal mining has caused serious surface subsidence in the eastern part of the city. In this study, multitemporal and multisource satellite remote sensing data were used to detect subsidence and geomorphological changes associated with underground coal mining over a 10-year period (1996–2006). A digital elevation model (DEM) was generated through Synthetic Aperture Radar (SAR) interferometry processing using data from a pair of European Remote Sensing Satellite (ERS) SAR images acquired in 1996. In addition, a Shuttle Radar Topography Mission (SRTM) DEM obtained from data in 2000 and an Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM from 2006 were used for this study. The multitemporal DEMs indicated that the maximum vertical displacement due to subsidence was around 13 m from 1996 to 2006. Multitemporal ASTER images showed that the flooded water area associated with subsidence had increased by 1.73 km² over the same time period. Field investigations and ground level measurements confirmed that the results obtained from the multitemporal remote sensing data agreed well with ground truth data. This study demonstrates that DEMs derived from multisource satellite remote sensing data can provide a powerful tool to map geomorphological changes associated with underground mining activities.     

The present status of subsiding land vulnerable to roof collapse in the Jharia Coalfield,India, as obtained from shorter temporal baseline C-band DInSAR by smaller spatial subset unwrapped phase profiling

In this paper, we identified recently subsiding areas in Jharia Coalfield, Jharkhand, India from the shorter temporal baseline Radarsat-2 C-band interferometric synthetic aperture radar (InSAR) data pairs of 2012. Although shorter wavelength C-band differential InSAR (DInSAR) is more sensitive to slow deformation and better suited for higher precision land subsidence measurement, the dynamic and adverse land cover in mining areas and resulting temporal decorrelation problem poses a serious problem for DInSAR observation in mining areas. We used smaller temporal baseline data pairs and adopted InSAR coherence-guided incremental filtering with smaller moving windows to highlight the deformation fringes over temporal decorrelation noise. We identified the deformation fringes and validated them based on ground information to prepare the land subsidence map of the coalfield in 2012. Several new, previously unreported subsidence areas were detected in the present study with a total subsiding area of 6.9 km². The recent incidence of roof collapse on 15 November 2014 at Angar Patra village in Katras region of the coalfield where 45 houses collapsed and 10 people were injured is situated in a highly subsiding vulnerable area as obtained from the present study. Due to spatial discontinuities of InSAR coherence, DInSAR phase unwrapping for the entire study area in one go did not appear feasible. To avoid this problem, we performed DInSAR processing in smaller spatial subsets and unwrapping of the subset interferograms by a ‘minimum cost flow’ algorithm. Subsequently, we plotted unwrapped phase profiles across the deformation fringes and retrieved the maximum deformation phase with respect to background phase and translated them into radar line of sight (LOS) displacement rates. For obtaining the average subsidence rates, we adopted InSAR coherence-weighted LOS displacement rates taking into account the contribution of each data pair as a function of DInSAR phase quality of the fringe areas. Ground-based subsidence measurements by precision levelling were conducted in four test sites that had been undergoing active underground mining during the observation period. We compared space-borne DInSAR-based subsidence rates obtained by the adopted technique with precision levelling measurements. Overall, the results are found to agree well. In the four test sites with gentle to flat topography, land subsidence occurs at slow to moderate rates due to compression of in-filled material (resulting from sand stowing in underground mining), without any evidence of roof collapse. In such cases, the horizontal displacement component is less significant, and overall surface displacement occurs essentially in the vertical direction. However, we assessed the nature of subtle horizontal strain to infer relative shrinkage or dilation of the land surface which could be additive or subtractive to vertical displacement in DInSAR-based LOS displacement.     

基于D-InSAR技术的淮南矿区地面沉陷监测

针对常规的通过大地水准测量、GPS测量监测矿区地面沉陷的技术存在监测周期长、成本高、无法全面监测等缺陷,提出了一种基于D-InSAR技术的矿区地面沉陷监测方法。以淮南矿区为试验区,采用两轨法D-InSAR技术,利用该地区2个时相的ALOS PALSAR数据获取了淮南矿区试验时间段内的地面形变图,分析了淮南矿区各矿的地面沉陷信息。结果表明,煤矿开采区存在5~25 cm不同程度的沉陷,与实际情况相符,因此,基于D-InSAR技术的监测方法可以作为一种获取矿区大范围的地表沉陷信息的有效方法。     

Combining differential SAR interferometry and the probability integral method for three-dimensional deformation monitoring of mining areas

With the exploitation of coal resources, ground surface subsidence continues to occur in mining areas, destroying the ecological environment and significantly affecting the daily productivity and life of humans. The differential synthetic aperture radar interferometry (D-InSAR) technique is widely used to monitor ground surface deformation because of its unique advantages such as high accuracy and wide coverage. However, conventional D-InSAR technology provides only one-dimensional (1D) displacement monitoring along the radar line of sight (LOS). This article proposes a method based on an analysis of the mining subsidence law for true three-dimensional (3D) displacement monitoring by combining D-InSAR and a subsidence prediction model based on the probability integral method. In this approach, 1D displacement, obtained using D-InSAR, is then combined with the prediction model to obtain the 3D displacement of ground surface target points. Here, 3D displacement curves were obtained for the Fengfeng mining area (China) using RadarSat-2 images obtained on 9 January and 2 February 2011. True ground surface displacement was measured simultaneously by levelling when the 152under31 s working face was being exploited in Jiulong mine. Vertical displacement and inclined deformation calculated using the proposed method were compared with levelling survey data and the results showed average differences of 3.2 mm and 0.1 mm m^?1, respectively; the calculated maximum displacement in the east–west and south–north directions were 106 and 73 mm, respectively. The spatial distribution of the displacements was in accordance with the mining subsidence law. Thus, the new method can retrieve highly accurate 3D displacements caused by mining subsidence.     

An approach for accurately retrieving the vertical deformation component from two-track InSAR measurements

As is well known, both conventional differential synthetic aperture radar interferometry (D-InSAR) and multi-temporal synthetic aperture radar interferometry (MT-InSAR) have a common limitation that they only can measure the deformation component along the radar line of sight (LOS) direction. However, in the majority of disaster investigations, there is more interest in the vertical deformation component than that of the horizontal direction, for example, in measuring ground subsidence in urban areas, or ground subsidence due to underground mining. To estimate the vertical deformation component accurately, it is in theory necessary to create at least three independent equations, and solve the vertical, the North–South, and the East–West deformation components by exploiting at least three-track InSAR LOS measurements or combining at least two-track InSAR LOS measurements with azimuth measurements. However, these methods are greatly limited and sometimes not even practical because there is typically little chance of obtaining a three-track SAR data set covering the same area in the same time span, and the accuracy of the azimuth measurements is far lower than that of InSAR LOS measurements. In this article, we found that it is possible to solve for the vertical deformation component from two-track InSAR LOS measurements in some circumstances. Then, an approach for accurately retrieving the vertical deformation component from two-track InSAR LOS measurements is proposed, and the analytical expression is presented. The approach is illustrated through an investigation of the ground subsidence in an area of Beijing, China. Unlike previous methods, this approach can accurately retrieve the vertical deformation component from two-track InSAR LOS measurements, and provide more reliable results for improving the interpretation of ground subsidence phenomena.     

Integrated space geodesy for mapping land deformation over Choushui River Fluvial Plain,Taiwan

Groundwater is an important part of the precious water resources. As the fresh surface water resources become scarcer because of climate change, population growth, and industrial activities, more and more groundwater has been extracted to meet the demands of various water uses (e.g. municipal, industrial, and agricultural). Excessive groundwater extraction leads to severe ground subsidence which compromises the safety of surface and underground infrastructures. Modelling the effects of groundwater extraction is vital to the management and sustainable use of groundwater. However, results of such modelling have to be validated with inputs such as the field survey of ground subsidence. Levelling and continuous global positioning system (GPS) receiver networks are routinely used to collect these field measurements. Unfortunately, these techniques have limitations in terms of areal coverage and density of survey marks and, as a result, subsidence hot spots can be easily missed out. In order to provide a comprehensive picture of subsidence to aid geotechnical modelling and to assess the effectiveness of measures used to mitigate ground subsidence, satellite imaging radar interferometry techniques (interferometric synthetic aperture radar (InSAR) can be used to complement other deformation monitoring techniques. In this study, 20 Advanced Land Observing Satellite (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) images acquired from 31 December 2006 to 26 February 2011 were used to map the land displacement over the Choushui River Fluvial Plain (CRFP), Taiwan. The GPS measurements acquired at 10 continuously operating reference stations (CORS) were used to refine the orbit error in the each differential interferogram obtained from each radar image pair. The displacement time series over the distributed scatterers and the persistent scatterers were analysed. Several subsidence bowls were identified in CRFP. A quantitative comparison was conducted to compare the radar measurements to the GPS measurements over 36 GPS CORS stations. Good agreement between both measurements was observed with coefficient of determination (R²) of 0.97, absolute mean difference of 3.2 mm year^?1, and standard deviation of 4 mm year^?1. The InSAR-measured Line-of-Sight displacement and GPS-measured horizontal displacement were integrated to derive the vertical displacement map. Two displacement maps were generated using two ALOS-2 PALSAR-2 pairs acquired between 2015 and 2016. Similar subsidence patterns were found in the two maps compared to the 2006–2011 displacement rate map, suggesting the land over the same region might have continued to fall.     

Monitoring of slow ground deformation by ERS radar interferometry on the Vauvert salt mine (France): Comparison with ground-based measurement

The differential SAR Interferometry (DInSAR) technique has been applied to a test site near Vauvert (France) to detect and monitor ground deformation. This site corresponds to the location of an industrial exploitation of underground salt using the solution mining technique. An area of subsidence has been observed using in situ measurements. Despite conditions unfavorable for InSAR because of the vegetal cover, we show that radar remote sensing observations provide valuable information which substantially improves our knowledge of the phenomenon. An adaptive phase filtering process has been used to improve the coherence level. In particular, our study shows that the geometry of the subsidence bowl is different to that previously assumed using ground-based techniques only. The size of the subsidence bowl (8 km) is larger than expected. This information will be useful for further modeling of the deformation and to improve the coverage of the in situ measurement networks. It also shows that radar interferometry can be used for the long-term monitoring of such sites and to predict potential environmental issues.     

基于NPSI方法的西安市地裂缝灾害链地表形变监测与演化态势分析

受构造背景、地下水采掘、活断层等因素的综合作用,西安市地裂缝于近年间不断加剧,引发多处不均匀地表沉降,并演化为对地表及地下建筑物均有强烈破坏作用的城市地质灾害链。为探明西安市地裂缝灾害链区域近年间的地表沉降态势、辨析不均匀沉降与地裂缝发育的联系,引入在城市地表形变监测精度和可靠性方面具有显著优势的网络化永久散射体时序雷达干涉测量（NPSI）方法,以2017年3月至2018年3月间成像的15期Sentinel-1A卫星SAR影像为数据源针对西安地区开展了时序监测分析,结合水准测量数据验证发现NPSI监测结果的精度达到 ±4.75 mm。实验结果表明：西安市地裂缝正向西南郊发育,地下水采掘及地上地下工程的建设在地裂缝发育趋势下加剧了不均匀沉降灾害,在形成地裂缝灾害链的鱼化寨、电子城、曲江新区以及地铁3号线等危害严重区域需要实时监测地裂缝发育趋势,并合理规划地下水开采及工程建设活动。相关研究结果可为路政、城建等部门的业务工作及相关研究提供参考信息。     

基于D-InSAR技术的唐山万达广场地区沉降监测及其安全性评价

唐山市万达广场南边为煤矿区,由于开采引起的地面沉降,直接影响到城市建筑物的安全,因此,采用D-InSAR形变监测手段获取煤矿区及周边地区的沉降时空分布特征显得尤为重要。收集了2004年～2008年唐山市万达广场及周边地区的3景ENVISAT ASAR数据,使用合成孔径雷达差分干涉测量(D-InSAR)技术来研究地面沉降对万达广场建设及周边地区的影响,获得了两期地面沉降结果,并将结果结合GIS技术对广场的建设进行安全性分析。实验结果表明:唐山煤矿地下开采引起的地面沉降从2004年开始没有再向万达广场的方向扩展,所以不会影响万达广场的建设,从而保证了广场的顺利建设和周边地区安全、有序的发展,也进一步验证了D-InSAR技术在煤矿区地面沉降监测中应用的可行性。     

Data fusion for investigating land subsidence and coal fire hazards in a coal mining area

Coal mining areas all over the world are often threatened by serious environmental hazards such as the occurrence of coal fires, land subsidence, etc. Coal fires burn away the natural non-renewable coal resources, locally raise the temperature of the area, emit polluting gases such as oxides of carbon, sulphur and nitrogen, and when present underground are even the cause of land subsidence. Mining-induced subsidences, on the other hand, cause horizontal and vertical movements in the land surface, and open cracks and fissures that serve as inlets for oxygen, which in turn aggravate the problem of coal fires. These inter-related phenomena often render the mining areas unfit for human inhabitation and the commercial exploitation of coal nearly impossible in some parts. In this study, satellite data acquired in three regions of the electromagnetic spectrum, namely optical, thermal and microwave, along with field data, are used to identify the areas affected by coal fires and land subsidence in a coalfield in north-west China. Data fusion techniques are used for an integrated analysis of this complex problem.     

基于SBAS-InSAR的关闭矿井地表形变规律研究

矿井关闭后,煤岩体在应力、地下水、氧气等多种因素的作用下发生风化劣化、强度降低,采动破裂岩体的应力和承载能力发生改变,易导致采空区地表发生二次形变。现有研究多是集中在矿区开采过程中的地表形变监测,而针对矿井关闭后的地表形变监测研究较少。为探明矿井关闭后地表的形变规律,利用2015-12-21—2019-12-24的62景Sentinel-1A影像,采用短基线集干涉测量(SBAS-InSAR)技术获取了徐州西部关闭矿井(包括庞庄矿、夹河矿和张小楼矿)4 a内时间序列的地表形变结果和形变规律。研究结果表明:3个矿区内地表最大沉降速率达-48 mm/a,4 a累计最大下沉178 mm;夹河矿地表呈现先下沉后抬升现象,而庞庄矿、张小楼矿地表持续下沉;4 a内3个矿区地表最大倾斜变形为1.70 mm/m,最大曲率变形为-0.039 mm/m 2;通过插值统计计算得到庞庄矿、夹河矿和张小楼矿地表形变大于10 mm的形变面积分别为10.5,13.7,11.6 km 2,且呈逐年上升趋势。     

Investigation on mining subsidence over Appin–West Cliff Colliery using time-series SAR interferometry

Advanced Time Series InSAR (ATS-InSAR) generally refers to those TS-InSAR methods with an external distributed scatterer selection module, e.g. SqueeSAR?, and GEOS-ATSA (Advanced Time-Series Analysis). It is being known as a very efficient tool for monitoring ground deformation over suburban or non-urban regions with great success. However, research conducted using C-band Envisat-based ATS-InSAR failed to produce reasonable outcome within Appin Colliery primarily due to the underground mining effect, which is located in the southeastern corner of the Southern Coalfield, New South Wales, Australia; thus, the general underground mining pattern cannot be formed. This work presents a modified ATS-InSAR method for mapping the ground deformation over underground mining region. More specifically, in order to achieve the best outcome, a modified measurement scatterer (MS) pixel selection method is introduced by including less reliable MS pixels through an Inverse Distance Weighted-based integration method. In addition, the proposed method is also applied to C-band Sentinel-1 image stacks for testing purpose, and the final result proved to be efficient to offer sufficient information to the mining industry and government for risk management purpose.     

Back analysis for rock model surrounding underground roadways in coal mine based on black hole algorithm

The huge plastic deformation is the characteristic of the underground roadways in coal mine. Therefore, to compute the stability of underground roadways, a elastic–plastic constitutive model of surrounding rock must be obtained. Many elastic–plastic constitutive models for rock mass have been proposed. In this study, a generalized constitutive law for an elastic–plastic constitutive model is applied. Using this generalized constitutive law, the problem of model identification is transformed to a problem of parameter back analysis, which is a typical and complicated optimization. To improve the efficiency of the traditional optimization method, a black hole algorithm is applied in this study. Combining the generalized constitutive law for an elastic–plastic constitutive model and black hole algorithm, a new back analysis method for model identification of rocks surrounding underground roadways in coal mine is proposed. Using this new method, the elastic–plastic constitutive models for two underground roadways in Huainan coal mine has been back-calculated. The results are compared with those of traditional genetic algorithm, fast genetic algorithm and immune continuous ant colony algorithm, that proposed in previous studies. The results show that the new model back analysis algorithm can significantly improve the computation efficiency and the computation effect, and is a very good method for back analysis the rock model surrounding underground roadways in coal mine.     

Detection of mining related ground instabilities using the Permanent Scatterers technique—a case study in the east of France

Ground stability is a major concern for land use planning and both natural and anthropogenic risk assessment, especially in urbanized areas. Space‐borne differential radar interferometry provides a unique tool able to give a synoptic view of ground deformation with centimetric to millimetric vertical precision. Approaches for combining a wide range of radar images such as the permanent scatterers (PS) technique allow the estimation of the deformation history of single buildings. The PS approach has been exploited to investigate a test site particularly exposed to ground deformation hazards, namely the iron mining basin in Lorraine (France). In this Letter, a specific focus was set on the case of Roncourt, where precursor signs of a collapse affecting an area of ～300×300?m² have been identified.     

PSI-HSR: a new approach for representing Persistent Scatterer Interferometry (PSI) point targets using the hue and saturation scale

Persistent Scatterer Interferometry (PSI) is a recently developed radar processing technique that enables the measurement of ground motion with millimetre accuracy. This technique makes use of an analysis of multi-temporal Synthetic Aperture Radar (SAR) images to detect stable PSI point targets, so-called persistent scatterers (PSs). These PSs are usually integrated in a geographic information system (GIS) as a standard point-based map, using colour coding on their velocities, with separation of ascending and descending data. In this paper we introduce the PSI-HSR (PSI Hue–Saturation Representation), a new method for representing PSI point targets using the hue–saturation scale. The aim was to render a unique colour for each PSI point target, based on a combination of displacements assessed along two different lines of sight. PSI-HSR provides a straightforward way to describe ground movement. To demonstrate the usefulness of this approach, an example of detecting ground subsidence in the Pistoia-Prato-Firenze basin is presented.     

“三下”采煤地表变形预计系统的开发及应用

文章介绍了一种用于任意形状采煤工作面开采后地面沉陷变形预测的、融数据可视化录入、数据管理、数据处理、可视化输出于一体的"三下"采煤地表变形预计系统。简述了该系统的组成及功能,详细介绍了该系统的计算模型、算法实现及数据的录入与输出方法,并给出了应用实例。应用实例表明,该系统的预测计算较准确可靠,系统功能达到了设计要求。     

Recent subsidence in Tianjin,China: observations from multi-looking TerraSAR-X InSAR from 2009 to 2013

Tianjin, China, has been suggested to have serious ground subsidence due to excessive extraction of groundwater. It is essential to monitor this subsidence, which has potential hazards and risks. Time series InSAR (TS-InSAR), such as small baselines subset (SBAS), is a powerful tool that can monitor ground deformation with high accuracy and at high spatial resolution over a long time interval. However, the high computational complexity may exceed computer memory limit when high-spatial resolution SAR (such as TerraSAR-X, TSX) images are used. In this article, the multi-look approach is introduced to the SBAS tool from StaMPS/MTI (Stanford method for persistent scatter/multi-temporal InSAR) in order to balance the spatial resolution and subsidence information in detection. The looks used for multi-looking are first fixed in terms of the accuracy of deformation and the density of coherent points. Then, the recent subsidence in Tianjin is extracted using multi-looking SBAS based on 48 TSX images acquired from 2009 to 2013. The results are validated by levelling measurements with a root mean square error (RMSE) of 4.7 mm year^–1, which demonstrates that SBAS analysis can effectively monitor deformation based on multi-looking TSX acquisitions in the area under investigation. Besides, the results also show that Tianjin has been suffering from subsidence during this period, and there were two separate large subsidence basins located in this study area with more than 500 mm cumulative subsidence. Moreover, the subsidence rate increased after December 2010 in Tianjin.     

开采厚度对沿空掘巷围岩稳定性的影响分析

目前大多煤矿根据平均开采厚度来确定煤柱宽度,进而指导沿空掘巷,然而煤层在形成过程中受各种因素影响,存在同一煤层厚度变化较大的情况。针对综放工作面煤层开采厚度变化大,导致沿空掘巷围岩变形差异大及破坏机理复杂等问题,采用FLAC 3D软件建立巷道模型,分析平均开采厚度下的围岩变形和破坏规律,并确定合理的煤柱宽度:平均开采厚度为18 m时,在实体煤帮侧,煤体内支承压力峰值与煤柱宽度呈正相关,且煤柱宽度大于8 m后,支承压力增长幅度变缓,因此合理的煤柱宽度应为8 m。在煤柱宽度确定的情况下,研究开采厚度对沿空掘巷围岩稳定性的影响,结果表明:煤柱宽度为8 m时,随着开采厚度的增加,顶板剪破坏面积增大,覆岩变形范围与顶板下沉量增大,但两帮剪破坏面积和两帮移近量减小;当煤层开采厚度小于18 m时,煤柱内支承压力峰值与煤层开采厚度呈负相关;当煤层开采厚度大于18 m时,煤柱内支承压力峰值与煤层开采厚度呈正相关,但增长幅度较小。根据数值模拟结果得出结论:开采厚度的增大对沿空巷道两帮的围岩控制有一定益处,但对顶板维护不利,对开采厚度较大的部位应及时补加锚杆进行强化支护。现场实际应用验证了本文研究的可靠性和有效性。