Trans. Nonferrous Met. Soc. China 24(2014) 1248-1254
Numerical simulation of dynamic surface deformation based on DInSAR monitoring
Yue-guan YAN1,2, Hua-yang DAI1, Lin-lin GE2, Jun-ting GUO1, Alex Hay-Man NG2, Xiao-jing LI2
1. College of Geoscience and Surveying Engineering, China University of Mining and Technology, Beijing 100083, China;
2. Geoscience and Earth Observing Systems Group, School of Civil & Environmental Engineering, University of New South Wales, Sydney 2052, Australia
Received 10 July 2013; accepted 24 December 2013
Abstract: Differential interferometric synthetic aperture radar (DInSAR) technology is a new method to monitor the dynamic surface subsidence. It can monitor the large scope of dynamic deformation process of surface subsidence basin and better reflect the surface subsidence form in different stages. But under the influence of factors such as noise and other factors, the tilt and horizontal deformation curves regularity calculated by DInSAR data are poorer and the actual deviation is larger. The tilt and horizontal deformations are the important indices for the safety of surface objects protection. Numerical simulation method was used to study the dynamic deformation of LW32 of West Cliff colliery in Australia based on the DInSAR monitoring data. The result indicates that the subsidence curves of two methods fit well and the correlation coefficient is more than 95%. The other deformations calculated by numerical simulation results are close to the theory form. Therefore, considering the influence, the surface and its subsidiary structures and buildings due to mining, the numerical simulation method based on the DInSAR data can reveal the distribution rules of the surface dynamic deformation values and supply the shortcomings of DInSAR technology. The research shows that the method has good applicability and can provide reference for similar situation.
Key words: DInSAR monitoring; numerical simulation; dynamic surface deformation; supplement
1 Introduction
In Australia, most underground coal mines that operate at a depth greater than 300 m mostly employ the longwall mining technique [1]. However, compared with room and pillar mining and other mining ways, longwall mining brings the serious surface deformation and threats the environment, the structures and buildings, railways and ancillary facilities. It affects the normal production of people life. Surface movement and deformation is a dynamic process. Therefore, it is very important to understand the rules of dynamic surface deformation due to longwall mining.
DInSAR is a new deformation monitoring technology in recent years. It is able to provide a large scope of deformation information of continuous surface in a short period. The accuracy can reach centimeter, even millimeter [2-5]. There have been many studies showing that the subsidence curves measured by DInSAR are consistent with the actual condition and the subsidence basin can be obtained. But in fact, the other deformation values (such as tilt, curvature and horizontal deformation) are the key to the most important protected objects. Usually, DInSAR data are influenced by noise and other factors and their deformation and actual curve have great deflection. So it is hard to actual application. Under this circumstance of no other actual measurement data, to reveal the surface dynamic deformation rules to protect the surface objects, other ways must be found to supply the shortcomings of DInSAR. The scholars both at home and abroad use Knothe function [6-9] and numerical simulation method [10-13] to study dynamic rules. But the parameter choice of these methods is mostly relying on the experience. No measured data could be referenced. So, based on DInSAR actual monitoring, dynamic mining of LW32 of West Cliff colliery in Australia is taken to study numerical simulation analysis of dynamic surface deformation.
2 Geology and coal seam condition of West Cliff colliery
The test area extracted by West Cliff colliery belongs to the southern coalfield in Australia, one part of the Sydney-Gunnadah Basin. There are several Permian period coal seams in this area and the top coal seam is Bulli seam. From the typical stratigraphic section of Southern Coalfield [14], the overlying rocks are major sandstone interbedded with other rocks and though shales and claystones are quite extensive in places, the sandstone predominates whose properties usually belong to the medium hard rock.
The longwall plan of West Cliff colliery is extracted from Bulli seam [1] shown in Fig. 1 [15]. LW31 was completed in February 2007. LW32 was mined from 12 February 2007 to 08 June 2008. The mining speed was slower before the middle May 2007 (around 23-30 m per week). The speed was accelerated from middle May 2007 (around 32-45 m per week). LW33 was started from the end of July 2008. The length of LW32 was 3222 m and the width was 305 m. The depth was 470-540 m; the thickness was 2.2-2.8 m; the dip angle is 2°. The mining method was the longwall mining.
Fig.1 Longwall plan of West Cliff colliery (unit: mm)
3 DInSAR results about studying area
In order to study the ground movement of LW32 in West Cliff colliery, ten ALOS PALSAR data acquired from 29 June 2007 to 01 October 2008 were used. Nine differential interferograms were generated with temporal baseline of one revisiting cycle (46 d) except the last pair (92 d). The 2-pass method [5] was used for data processing. The interferograms including deformation phase were generated by two repeat orbit SAR images. Then an external high precision digital elevation model (DEM) was subtracted from the interferogram. The track parameters are shown in Table 1 and the deformation value is calculated.
Table 1 ALOS PALSAR data condition (Ascending)
25 m resolution DEM provided by the New South Wales land bureau in Australia was used and the different stages of DInSAR results in line-of-sight (LOS) direction are shown in Fig. 2 [15]. The negative value indicates the subsidence.
4 Numerical simulation analyses
4.1 Model description
In order to ensure rapid and accurate simulation calculation, considering the typical stratigraphic section and the dimensions of the LW31 and LW32 (Table 2), the model with dimensions of 4000 m in longitudinal (x) direction, 2000 m in transverse (y) direction and 500 m in vertical (z) direction was selected as the initial geometry (Fig. 3). The entire model was divided by brick element, and the whole calculation modeling had a total of 332800 zones and 391245 grid points. Four sides of model have displacement boundary conditions, four vertical boundary faces and undersides were fixed, and the surface was free. The Mohr-Coulomb yield criterion of the rock mass destruction judgment was adopted [12]. The calculation parameters of the lithology and coal characteristics used in the model were chosen according to the geological investigation and similar projects.
To compare the subsidence with the DInSAR measured results, the mining distances need to be equal to the actual mining distance at different time stages. The residual subsidence influence of LW31 could not be ignored when LW32 was extracted. However, the information on the extraction time for LW31 was not available, therefore the calculation extracted speed was uniform velocity which was 200 m per month.
4.2 Modelling results
Using the software ArcGIS, the dynamic surface deformation process extracted in LW32 could be achieved. The vertical subsidence and the surface horizontal movement along the longitudinal and transverse directions at different time stages can be calculated as shown in Fig. 4. The black box is the relevant mining area.
Fig. 2 Subsidence in LOS direction measured by DInSAR computing time started from 29/06/2007 to 14/08/2007 (a), 29/09/2007 (b), 14/11/2007 (c), 30/12/2007 (d), 14/02/2008 (e), 31/03/2008 (f), 16/05/2008 (g), 01/07/2008 (h) and 01/10/2008 (i)
Table 2 Dimensions of longwalls 31 and 32
Fig. 3 Initial geometry model and position of workface of West Cliff colliery
5 Comparison between DInSAR and modeling results
The displacement measured from DInSAR is in LOS direction. Therefore, the 3D displacements calculated from modeling need to project into the same direction to compare with the DInSAR measurements. Figure 5 shows the modeling displacement along the LOS direction of LW32 at different time stages calculated using the following equation [16]:
[con θ -sin θ cos α sin θ sin α]
where θ is the incidence angle at the target point; α is the satellite heading vector angle (positive clockwise from north); DU is the vertical displacement; DE is the easting displacement; DN is the northing displacement; △R is the ground displacement between two acquisitions along LOS direction.
Fig. 4 Vertical subsidence of LW32 of modeling results computing time started from 29/06/2007 to 14/08/2007 (a), 29/09/2007 (b), 14/11/2007 (c), 30/12/2007 (d), 14/02/2008 (e), 31/03/2008 (f), 16/05/2008 (g), 01/07/2008 (h) and 01/10/2008 (i) (unit: mm)
Fig. 5 Modelling displacement along LOS direction computing time started from 29/06/2007 to 14/08/2007 (a), 29/09/2007 (b), 14/11/2007 (c), 30/12/2007 (d), 14/02/2008 (e), 31/03/2008 (f), 16/05/2008 (g), 01/07/2008 (h) and 01/10/2008 (i) (unit: mm)
Fig. 6 Comparison of DInSAR and modeling subsidence
The subsidence curves of DInSAR and numerical modeling along the longitudinal and transverse directions at different mining stages are shown in Fig. 6 [15]. As can be seen from the graph, good correlations between displacements along the longitudinal and transverse directions estimated from both methods were observed. The correlation coefficient is more than 95%. The boundary angle is 45°. It is also indicated that the results of numerical modeling are the same with actual condition.
Tilt is an first order derivative of subsidence. The tilt contrast of DInSAR and numerical simulation along the longitudinal and transverse directions are shown in Fig. 7. With the working face advancing, the maximum tilt value of open-off cut located at the same position increases gradually, and finally tends to be stable. The results of DInSAR and numerical simulation show performance for the deformation value gradually increasing with face of advance, and the maximum deformation value and the effects range between results are closer along two directions. The tilt deformation curves of numerical simulation are similar with theory distribution. But tilt deformation curves distribution of DInSAR have large fluctuations. The maximum tilt is always at the middle of the workface. This cannot reflect the actual surface deformation, and regularity is not apparent. Measured and theoretical tilt deformation and distributions should be consistent if there are no special geological and geomorphologic conditions in the study area. It is due to mainly noise and other factors that cause the phenomenon of inconsistent.
Curvature is second derivative of subsidence. Figure 8 shows curvature deformation comparison chart of DInSAR and numerical simulation along the transverse directions. The curvature values of DInSAR deviated from the theoretical ones are large along the transverse main sections. The maximum curvature value is not stable and there are no regular rules to be found. But the results of numerical simulation are closer to the theoretical form. From Fig. 8(b), there are three areas on the surface. The tensile deformation zones are on both sides and the compressive deformation zone is in the middle of mining-out area. Being affected by LW31 mining, the maximum curvature value of this side is less than the other side, but the influence range is increased.
The numerical simulation with the appropriate parameters can well reflect the dynamic surface subsidence based on DInSAR monitoring data. But due to the noise and other factors impact on DInSAR data, its analysis of dynamic deformation of tilt and curvature is a little poor than numerical simulation. It is difficult to evaluate the surface deformation and the buildings of the safety implications only with DInSAR data. Therefore, when researching on the surface and its subsidiary structures and buildings and calculating other surface deformation parameters due to longwall mining, the numerical simulation method can supply the shortcomings of DInSAR technology. Both unions can better describe the dynamic surface movement deformation and provide the basis for safe mining under buildings.
Fig. 7 Tilt comparison of DInSAR and modeling
Fig. 8 Curvature comparison of DInSAR and modeling
6 Conclusions
1) DInSAR technology can reflect the actual dynamic surface subsidence with all weather and 24-hours. But due to the noise and other factors, the DInSAR monitoring results are difficult to reflect the actual surface tilt and curvature deformation.
2) The numerical simulation method is used for modeling dynamic surface deformation based on DInSAR monitoring data. The subsidence correlation between simulation and DInSAR results is high, and the tilt and curvature deformations calculated from modeling are consistent with theory form.
3) When the tilt and curvature deformations need to be calculated based on DInSAR data, the numerical simulation method can supply the shortcomings of DInSAR technology.
Acknowledgements
This research work has been supported by the Cooperative Research Centre for Spatial Information through Project 4.09 funded by the Australian Commonwealth’s Cooperative Research Centers Program. The Australian Coal Association Research Program has also supported research for ground subsidence monitoring using DInSAR.
The authors wish to thank the Earth Remote Sensing Data Analysis Center (ERSDAC) for providing the ALOS PALSAR data. The PALSAR L-1.1 products were produced and distributed to the IAG Consortium for the Mining Subsidence Monitoring.
References
[1] NSWDOP. Impacts of underground coal mining on natural features in the southern coalfield: Strategic review [R]. Australia: New South Wales Department of Planning, 2008.
[2] COLESANTI C, FERRETTI A, NOVALI F, PRATI C, ROCCA F. SAR monitoring of progressive and seasonal ground deformation using the permanent scatterers technique[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(7): 1685-1701.
[3] GE L L, CHANG H C, RIZOS C. Mine subsidence monitoring using multi-source satellite SAR images [J]. Photogrammetric Engineering and Remote Sensing, 2007, 73(3): 259-266.
[4] STOW R J, WRIGHT P. Mining subsidence land surveying by SAR interferometry [C]//Proceedings of the 3rd ERS Symposium. Noordwijk, Pays-Bas: ESA Publications Division, 1997: 525-530.
[5] MASSONNET M, ROSSI D, CARMONA C, ADRAGNA F, PELTZER G, FEIGL K, RABAUTE T. The displacement field of the landers earthquake mapped by radar interferometry [J]. Nature, 1993, 364: 138-142.
[6] TENG Yong-hai. Computation study of ground movement and deformation during mining process [J]. Mine Survey, 1997(4): 17-20, 50. (in Chinese)
[7] MONGE L R A, OYANGUREN P R, FUENTE R D L. Subsidence development with time-experiences from longwall operations in the Appalachian coalfield [J]. International Journal of Mining and Geological Engineering, 1990, 8(2): 61-73.
[8] LIAN X G, ANDREW J, JOSE S R, DAI H Y. Extending dynamic mining subsidence models [J]. Transactions of Nonferrous Metals Society of China, 2011, 21: s536-s542.
[9] WU Kan, HUANG Zhen-zhen, WANG Xin. Self-contained forecast model of mining subsidence [J]. Coal Mining Technology, 2006, 11(4): 4-6. (in Chinese)
[10] MONGE L R A, OYANGUREN P R, de la FUENTE R. A method to predict final subsidence basins by means of a finite difference computer code [J]. Journal of Rock Mechanics and Mining Sciencces and Geomechanics Abstrats, 1996, 33(8): 251-258.
[11] FALON S, GAVETE L, RUIZ A. A model to simulate the mining subsidence problem: Coding and implementation of the algorithm [J]. Computers & Geosciences, 1996, 22(8): 897-906.
[12] SHAHRIAR K, AMOUSHAHI S, ARABZADEH M. Prediction of surface subsidence due to inclined very shallow coal seam mining using FDM [C]//Underground Coal Operators’ Conference. Wollongong: University of Wollongong and the Australsian Institute of Mining and Metallurgy, 2009: 130-139.
[13] DAI Hua-yang, LI Wen-chang, LIU Yi-xin, JIANG Yao-dong. Numerical simulation of surface displacement laws under different unconsolidated layers thickness [J]. Transactions of Nonferrous Metals Society of China, 2011, 21: s599-s603.
[14] MSEC. Report on the prediction of subsidence parameters and the assessment of mine subsidence impacts on natural features and surface infrastructure resulting from the extraction of proposed longwalls 34 to 36 in Area 5 at West Cliff Colliery in support of the SMP application [R]. Sydney: Mine Subsidence Engineering Consultants, 2007: 230.
[15] NG Alex Hay-Man, GE Lin-lin, YAN Yue-guan, LI Xiao-jing, CHANG Hsing-Chung, ZHANG Kui, RIZOS Chris. Mapping accumulated mine subsidence using small stack of SAR differential interferograms in the Southern coalfield of New South Wales, Australia [J]. Engineering Geology, 2010, 115: 1-15.
[16] FIALKO Y, SIMONS M, AGNEW D. The complete (3-D) surface displacement field in the epicentral area of the 1999 Mw 7.1 Hector Mine Earthquake, California, from space geodetic observations [J]. Geophysical Research Letters, 2001, 28(16): 3063-3066.
基于DInSAR监测的动态地表移动变形数值模拟分析
阎跃观1,2,戴华阳1,葛林林2,郭俊廷1,Alex Hay-Man NG2,李晓静2
1. 中国矿业大学(北京) 地球科学与测绘工程学院,北京 100083;
2. Geoscience and Earth Observing Systems Group, School of Civil & Environmental Engineering, University of New South Wales, Sydney 2052, Australia
摘 要:合成孔径雷达干涉测量(DInSAR)技术是监测动态地表沉陷的新方法,能够大范围地监测地表下沉盆地的动态演化过程,可以较好地反映各阶段的地表下沉形态。但是受到噪声等因素的影响,基于DInSAR数据计算的倾斜、水平变形等曲线规律性较差,与实际偏差较大,而倾斜、水平变形又是研究地表受护对象是否安全的重要指标。基于DInSAR监测数据,采用数值模拟方法,对澳大利亚West Cliff煤矿长壁工作面32开采动态地表移动变形进行分析研究。结果表明,数值模拟得到的下沉曲线与DInSAR实测数据比较吻合,相关系数可达95%以上,且数值模拟方法求算的地表倾斜、曲率和水平变形与理论形态一致。因此,进行矿区地表及其附属建(构)筑物受采动影响分析时,应以DInSAR实测数据为基础确定模型参数,结合数值模拟方法揭示地表动态其他移动变形值的分布规律。研究表明,该方法具有良好的适用性,可为类似的情况提供参考。
关键词:DInSAR监测;数值模拟;动态地表移动变形;补充
(Edited by Hua YANG)
Foundation item: Project (20110023110014) supported by Specialized Research Fund for the Doctoral Program of Higher Education of China; Project (2010QD01) supported by Fundamental Research Funds for the Central Universities, China
Corresponding author: Yue-guan YAN; Tel: +86-13811548084; E-mail: yyg_0720@163.com
DOI: 10.1016/S1003-6326(14)63186-1