Field Measurement and Numerical Modelling Study on Mining-Induced Subsidence in a Typical Underground Mining Area of Northwestern China

Mining-induced subsidence is a great concern for environmental protection in underground mining areas in China and all over the world. In view of the fact that the research on land degradation above underground coal mines are completely or partially independent of coal mining activity and the fact that the mechanism behind mining-induced subsidence has not been well understood, this study presents a field measurement and numerical study of mining-induced subsidence with respect to mining activity of three adjacent longwall panels in a coal mine in Northwest China. -is study shows that surface subsidence lags far behind panel extraction or mining activity. -e profiles of ground surface are dominated and manifested by the subsurface strata structures. -e subsidence influence throughout the whole length of a longwall panel varies. Stability of strata structures within overburden before the final subsidence controls the stability of ground surface land. Chain pillars of 20 m between panels of 240 m wide with cover depth of 600 m have been crushed in the gob and do not have any function in supporting the overburden strata. -e final subsidence of the three adjacent panels is far to come in the future and the land reuse above underground coal mines should be carefully planned by making sure that the gob is completely compacted or no potential secondary subsidence occurs in the future.


Introduction
ere are many environmental issues associated with mining including air pollution due to dust and gas emission, water contamination, loss of groundwater resources, destruction to the surface structures, and degradation of ecosystems [1,2]. Underground mining, in particular, has a significant influence on ground movement or is known as subsidence, which induces serious environmental disturbances on land [3]. Many studies have been conducted on land degradation above underground coal mine workings such as water content change within the soil, evolution of physical and chemical properties of soil, and vegetation deterioration as shown in Figure 1.
In fact, back to the early days of the 1980s and 1990s in the US, the impact of mine subsidence on aquifer, farmland, agriculture, and other environmental issues had been studied [4][5][6]. Booth et al. [7] discussed the impact of mining on groundwater including decline and recovery of water levels. Utilizing a probability (frequency ratio) model, a statistical (logistic regression) model, and a Geographic Information System (GIS), [8] constructed a subsidence hazard map around abandoned underground coal mines at Samcheok City, Korea. Reference [9] conducted a comparison study on groundwater system before and after longwall mining with the condition of shallow cover depth.
Using subsurface subsidence prediction model, Cheng et al. [10] carried out a study on strata deformation over longwall coal pillars. Marino et al. [11] presented an empirical approach for evaluating ground movement due to mining, which allows for realistic understanding of the historical and adverse subsidence-induced movements. Shi et al. [12] carried out a study on the influence of mining-induced subsidence on soil bacterial communities in aeolian sand area of western China. Marino et al. [13] proposed a highaccuracy method to identify the location of mining disturbances and the time when they occurred by analyzing Landsat data. Guney and Gal [14] reported a case study analysis of surface subsidence of a Turkey longwall mine under weak geological conditions. Ma et al. [15] presented a study of soil fertility decrease due to mining-induced subsidence in arid and semiarid regions. A finite element analysis was performed by Marian et al. [16] to study the state of stresses on the structures of buildings subjected to the impact of underground mining of hard coal seams in the Jiu Valley Basin (Romania).
As for the land destruction due to mining-induced subsidence, a number of studies on mitigation and reclamation of land have been carried out. Dawkins [17] discussed potential management and rehabilitation requirements of environmental effects from longwall subsidence on streams, lakes, and groundwater systems. Haibach et al. [18] provided some mitigation techniques to cope with the impacts of subsidence on streams. Hu et al. [19] presented a new reclamation technology, which is called concurrent mining and reclamation plans for yet-to-be-stable subsiding land against a case study mining area in northern Anhui province, China. Li et al. [20] carried out a comparative study of bacterial diversity of soils from mining subsidence and reclamation areas through pyrosequencing analysis. Based on replicated analysis across 18 mitigation sites, a study of the restoration of macroinvertebrates, fish, and habitats in streams following mining subsidence is presented by Tim et al. (2017). Wang et al. [21] carried out a study on separation and fracturing in overlying strata disturbed by longwall mining in a mineral deposit seam. Bi et al. [22] studied the alleviation of root damage stress induced by simulated coal mining subsidence ground fissures by arbuscular mycorrhizal fungi. In view of the significance of reclamation time in coal mining subsidence areas after mining, Wang et al. [23] revealed that fertilization can shorten the restoration time and reclamation time is crucial for the restoration of soil bacterial communities. e mechanism behind land subsidence due to mining is in urgent need to take positive measures for environmental protection in the future mining which is good for homonymous development of mining industry and local environment. More appropriate measures are only able to be taken when the mechanism behind it is cleared. However, the above findings of surface subsidence or physical or chemical changes of soil and mining activity are completely or partially independent of each other; the relationship between these findings and mining activity has not been studied deeply or completely. In another word, the mechanism behind these findings has not been well understood. Instead of taking some positive measures, some passive measures are usually taken to protect associated hazards. And most importantly, those studies were mostly on large scale or in macrolevel. As a result, the accuracy and dimension of the above studies are hard to be adopted to guide land restoration and reclamation specifically and accordingly. As is known, the degree of damage of mining on land or soil cannot be evaluated precisely, because the damage degrees are different from region to region, mine to mine, panel to panel, and even patch to patch. For instance, subsidence changes the structure of the overlying strata; deformation of the surface above the panel center is different from that above the barrier pillars or chain pillars. e damage that topsoil suffered must be different; then, the vegetation restoration plan must be adjusted to those conditions, so as to reach better land protection or reclamation results. erefore, this study presents a study of mining subsidence with respect to mining activity of three adjacent longwall panels in a coal mine in Northwest China in order to provide some insights which we hope can provide a reference for the land subsidence prediction and land rehabilitation measure selection.

Characteristics of the Study Area
e case study mine is located in the southernmost region of Wushen County, Ordos, Inner Mongolia, China, as shown in Figure 1. e mine is 63 km away from Wushen County government. e annual output of the mine is 8 million tons per year and the service life is 71 years. Elevation of the surface of the region is +1138.0～+1125.1 m above sea level. e surface is mainly aeolian sand dune field dispersed with many stable and quasiflow crescent-shaped sand dunes. e natural vegetation consists of typically poplar trees. No important or large buildings or houses are there. e mine is using a fully mechanized longwall mining method. e case study longwall panels are located in the southern-eastern part of the first miming district which is shown in Figure 2.
e case study longwall panels are

Field Measurement Instrumentation Plan
ere are 4 monitoring lines above the case study longwall panels with one line across the three panels in the dip direction and two-line along the strike direction as shown in Figure 4.
Line N is in the middle of 31102 panel across the setup room and along the strike with a length of 1275 m. ere are 50 points in total, which are numbered from N01 to N45 with a spacing of 25 m. Monitoring of this line started from Nov 21, 2017, to Oct 24, 2018. Line A is in the middle of 31103 panel. 62 monitoring points are numbered AK1-AK3 and A1-A61; the data were collected from July 27, 2019, to July 5, 2020, with 22 groups of data. e invalid points include AK1, AK2, A1, A3, A7, A10, A11, A14, A21, A24, A54, and A56 due to weather or other reasons.
For each point, a survey peg was mounted using C20 cement with steel reinforcement. e peg dimension is 0.5 m × 0.5 m × 2 m as shown in Figure 5. e buried depth of each peg is no less than 0.6 m. e number of each point is marked on the top of the peg. en total stations are used to measure the displacement of each point.
During the early stage of panel extraction, monitoring was carried out comprehensively every month; at this stage, the surface subsidence rate was less than 50 mm/month. During the fast subsidence stage, 2-3-times of survey had to be carried out each month. For the later stage when the subsidence nearly ceased (less than 30 mm subsidence in a period of 6 months), 1-2 times of complete survey were required.

Field Measurement Data Analysis and Discussion
Survey for Line N was conducted 15 times from 11/16/2017 to 5/14/2019 as shown in Table 2. e profile of this line is shown in Figure 6.

Advances in Civil Engineering
As we can see, the subsidence did not cease yet until May 14, 2019, although Point N 45 was 684 m away from the setup room and the face advance was 3099 m. It can be found from the plot that subsidence lagged far behind the panel extraction. Although the thickness of the coal seam is 5.5 m, after the final extraction of 3099 m of the panel on May 14, 2019, the maximum subsidence was still less than 1.5 m,    However, for a period of time of over 8 months from Sep 12, 2018, to May 14, 2019, the subsidence increment of each point was less than 20 mm, that is, less than 2.5 mm per month. is seems to indicate that the subsidence ceased. e abrupt adjustment of subsidence occurred between Dec 2, 2017, when the face advance was 932 m, and Feb 5, 2018, when the face was 1268 m. e face had already passed Point N 45 by almost 250 m on Dec 2, 2017. During this period of time of 2 months, the subsidence increased over 1 m. erefore, it is deemed that this was due to the dramatic structural adjustment of masonry broken strata that bridged with the overburden.
Survey for Line Q was conducted 15 times from 11/16/ 2017 to 5/14/2019 as shown in Table 3. e profile of this line is shown in Figure 7. For comparison, the final subsidence of panel 31101 is presented in Figure 8. As we can see, the maximum subsidence point for 31101 panel is Q 26 with a value of 1.486 m.
Due to the influence of the mined-out 31101 panel, the maximum subsidence point on Line Q was not located in the middle of 31102 panel (see Figure 4, Q33 or Q34); rather, it was Point Q 29 which was about 120 m away from the middle of the panel and located at the chain pillar edge between 31101 and 31102 panels. e maximum subsidence value of Point Q 29 was 2.27 m. It is also noted that the profiles of all the curves did not present a flat-bottomed characteristic; that is, supercritical subsidence basin did not  NO3  NO2  NO1  N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12  N13  N14  N15  N16  N17  N18  N19  N20  N21  N22  N23  N24  N25  N26  N27  N28  N29  N30  N31  N32  N33  N34  N35  N36  N37  N38  N39  N40  N41  N42  N43  N44  N45 Point no.  Survey of Line A was conducted for a long period from 7/ 27/2019 to 7/5/2020. e details of survey times and corresponding face advance distances are given in Table 4. 22 groups of data were collected. e profile of this line based on the data is shown in Figure 9. KQ1  KQ2  KQ3  Q1  Q2  Q3  Q4  Q5  Q6  Q7  Q8  Q9  Q10  Q11  Q12  Q13  Q14  Q15  Q16  Q17  Q18  Q20  Q21  Q22  Q23  Q24  Q25  Q26  Q27  Q28  Q29  Q30  Q31  Q32  Q33  Q34  Q35  Q36  Q37  Q38  Q39  Q40  Q41  Q42  Q43  Q44  Q45  Q46  Q47  Q48  W1  W2  W3  W4  W5  W6  W7 2017.11.21 (extraction of 2600 m)  Advances in Civil Engineering It shows that the profiles after 08/16/2019 stayed similar. is indicates that the structures were formed within the overburden strata. e profiles of ground surface are dominated and manifested by the subsurface strata structures [24]. In other words, the profile of the ground surface is representative of the unseen structures within the subsurface strata. e profiles after 08/16/2019 only slightly changed. e subsidence value increased with time but the profiles stayed still similar, and even most parts of profiles after 04/ 26/2020 were overlapped from what time the subsidence had slowed down significantly after a big adjustment. We can deduce from the change of the profiles that the structures of the overburden strata hardly changed or may be only adjusted slightly. e subsidence magnitude of each point increased while still remaining the profile that first appeared on 08/16/2019. Extraction of 31102 panel started from 4/1/2020. In order to eliminate the influence of extraction of 31102 panel on subsidence analysis of single 31103 panel, a group of data on 3/ 17/2020 was chosen for curve fitting as shown in Figure 10. It shows that the maximum subsidence point was A43. e subsidence value and rate were given in Table 5. e subsidence rate of the point is presented in Figure 11(a), and the subsidence value of the point is presented in Figure 11(b).
Point no.  Ak1  Ak3  A-2  A-4  A-6  A-8  A-10  A_12  A_14  A_16  A_18  A_20  A_22  A_24  A_26  A_28  A_30  A_32  A_34  A_36  A_38  A_40  A_42  A_45  A_47  A_49  A_51  A_54  A_56  A_58     Advances in Civil Engineering e data and figures show that, despite the large face advance distance of 1657 m, however, the ground movement still did not cease, because the coal seam averages about 5.5 m, but the maximum subsidence now was only 1.52 m. From the subsidence rate of point A43, it is demonstrated that subsidence was becoming slower and slower as time went by. Nevertheless, subsidence is going to last for a long time and in some cases, over ten years before the final subsidence is achieved. is phenomenon indicates that the land reuse above underground coal mines should be carefully planned making sure that no potential secondary subsidence occurs in the future, especially sudden and violent future subsidence.
Monitoring for Line Q started from 7/27/2019; 22 groups of data were collected by 7/5/2020. e detail of the survey time and the corresponding face advance distance are shown in Table 6. e profile of the line is presented in Figure 12. Figure 12 is the profile of Line Q in the dip direction, while Figure 8 is the profile of Line A in the strike direction. Comparing the two figures, it is noticed that the maximum subsidence values of the two monitoring lines for one panel are of great difference. e maximum subsidence of Line Q is 2.74 m while that of Line A is 1.52 m. e reason for this is that the structures formed by the overburden strata are different along both the strike direction and dip direction. For example, the subsidence profile of Line Q must be different from another line parallel to but away from it. Figure 12 also indicates that there must be twice significant strata structure adjustments within the overburden. One is around 04/29/2020 near point Q44; the other one is around 06/17/2020 near point Q 40. But please note that the location of the strata structure adjustment is not vertically corresponding to the surface point. It is known that there is an area of influence on the surface subsidence due to mining. e plot also demonstrates that the maximum subsidence point is located in the middle of 31102 panel from Point Q28 to Point Q34 meaning that the chain pillars between the three panels did not have a symptomatic impact on surface subsidence. is may indicate that the 20 m of chain pillars had already been crushed and did not have any function in supporting the overburden strata.
In order to determine the maximum subsidence when the ground surface hardly moved, a group of data on 07/03/ 2020 was chosen for curve fitting as shown in Figure 12. It shows that the maximum subsidence point is Q32, with a little increase from Points Q28 to Q34. is more or less indicates a characteristic of supercritical subsidence basin. However, the mining height is 5.5 m; the maximum subsidence value is only 2.74. And 2.74/5.5 � 49.8%. In other words, only 49.8% of the maximum possible subsidence was achieved.
is indicates that the combined width of the adjacent three panels was far from critical width. But considering the flat-bottomed depression from the profile, it is more likely that the final subsidence is still to come in the future. e subsidence velocities of the point at different times are presented in Table 7 and the corresponding figure with fitted curve is shown in Figure 13. e subsidence magnitude with respect to time is presented in Figure 14. e inference that can be drawn for explaining the subsidence is given in Figure 15.
It is revealed by the data and figures that the stability of strata structures within overburden before the final subsidence controls the profile of the ground surface. As caving of roof strata is bulking controlled process, the compaction of the caved zone is critical for stability evaluation of the ground surface. If the subsidence is complete which means that there is no chance for potential future subsidence, then it is safe to reuse the land. For the case study coal mine, the surface land above the panels is not safe for reuse at the moment and one or two years ahead.

Numerical Modelling Study and Discussion
Due to the limitation of the number and the density of monitoring points on the ground surface in the field study, it was impossible to take a more detailed look at the strata movement.
erefore, a numerical modelling study was carried out in which any point in the model can be monitored. As the panel length or advance distance, to be more precise, is significantly larger than panel width, the width is the controlling factor on surface subsidence profile. erefore, the extraction of the three panels, that is, 31101, 31102, and 31103 panels, in the dip direction was simulated. e constructed model is given in Figure 16. Finer zone resolution in the vicinity of the material boundaries was used. e side boundaries were roller constrained, and the bottom boundary was fixed both horizontally and vertically. Mohr-Coulomb criterion was used to determine the failure of the rock groups. e rock mass engineering parameters were developed based on laboratory tests, supplemented by geologic data and field observation. en by trial and error, the original parameters were modified slightly by comparing the field subsidence data with the simulated data. e calibrated parameters used for numerical modelling are given in Table 8.
As shown in Figure 8, the monitored maximum subsidence of 31101 panel is 0.15 m. e numerical simulation result is 0.21 m (see Figure 17). e simulated result and the monitored result are in good agreement with each other. After extraction of the 31101 and 31102 panels, as shown in Figure 18, the monitored maximum subsidence in Figure 7 is 2.27 m, and simulated result in Figure 18 is 2.58 m. is result is believed to be reliable. In numerical modelling, the result is obtained when the calculation is totally completed; that is, the state of equilibrium of the whole model is reached. However, in the field, the subsidence is still not settled completely. erefore, the slight deviation is normal and the result from the numerical modelling is expected to be slightly larger than the monitored value. e numerical modelling result demonstrates that the pillar between the two panels has no influence on surface subsidence. e field measurement also demonstrates no pillar influence on the surface subsidence. Please also note that both field measurement and numerical modelling do not see the maximum possible subsidence at the center of the trough meaning that the critical width of extraction is not reached even though the adjacent two longwall panels are extracted with a 20 m of pillar left in between. is may be due to the large panel depth of over 600 m. Figure 19 shows that the simulated maximum subsidence after extraction of the adjacent three panels is 3.92 m. e field monitored result is only 2.72 m (see Figure 12). is result is also believed to be reasonable although the simulated result is much larger than field data because the field data was updated on July 03, 2020. But the subsidence is expected to be totally settled in 1 or 2 or more years. Also, the pillars do not have an obvious effect on the surface subsidence profile. It is noticed that a flat-bottomed depression  -3  W20  W19  W18  W17  W15  W14  W13  W12  W11  W10  W9  W8  W7  W6  W5  W4  W3  W1  Q48  Q46  Q45  Q44  Q43  Q42  Q41  Q40  Q39  Q38  Q37  Q36  Q34  Q33  Q32  Q31  Q30  Q29  Q28  Q27  Q26  Z33  Q20  Q17   occurs in both numerical modelling and the field meaning that a supercritical width of extraction is reached. In the numerical modelling, the flat-bottomed trough is about 270 m; in the field, it is about 240 m from points Q26 to Q36 (see Figure 12). e results match well. Figures 20-22 show the strain increments after extraction of 31101, 31102, and 31103 panels sequentially. Strain increments ( Figure 20) indicate that the failure of rock strata penetrates from the edge of 31101 panel and propagates upwards in a shape of trumpet flower (two blue ovals). But the obvious propagation stops at about 245 m above the coal seam. is suggests that the extraction of the 31101 panel does not have a significant influence on the surface subsidence. However, after 31101 and 31102 panels were extracted as shown in Figure 21, the zones of strain increment extend both vertically and transversely within overburden strata and start from the gob edges. e vertical strain increments extend all the way to the surface. is finding is consistent with the field observation that some cracks were found on the boundary of the subsidence trough. After extraction of the three panels as shown in Figure 22, the area of the two strain increment zones highlighted in closed blue lines expands substantially compared with that after extraction of the panels. Many studies indicated that width-to-depth   Figure 19, it is also concluded that, after extraction of the three panels, the critical width of extraction is about to be reached. is is in good agreement with many other researchers' studies that panel width-to-depth ratio of 1.2 is the critical value for the occurrence of critical width of extraction supported by a wide range of surface subsidence data worldwide (Figure 23). If the influence of chain pillars between panels is ignored ( Figure 24 also shows that the pillars are completely crushed meaning that they have no significant supporting influence or function), then the width-to-depth ratio of extraction of the three panels is 760/ 600, equaling 1.27, which is slightly larger than 1.2. As we can see from the figures, the subsidence trough is more or less flat-bottomed although the proportion of the flat bottom     is rather small and the simulated maximum subsidence after extraction of 31103 panel is 3.92 out of 5.5 m (71.3%). All in all, the numerical modelling results are generally in good agreement with field monitoring results. Please note that FLAC3D which is a finite difference software was used for numerical modelling study. As a result, caving or surface cracks or structural formed by the overburden strata cannot be modelled explicitly. However, the advantage of the software is that the subsidence value of any point in the model can be monitored; it is more applicable for continuous and uniform subsidence with small subsidence values. For other geological conditions where large subsidence or cracks occur, UDEC or other types of similar software are better.

Conclusions
A field measurement and numerical study on mining-induced subsidence of three adjacent longwall panels in a coal mine in Northwest China is presented in this study. e panels are 241 m wide and over 3000 m long extracting a coal seam of 5.5 m thick.
It is found that surface subsidence lags far behind panel extraction or mining activity. e profiles of ground surface are dominated and manifested by the subsurface strata structures. Structures formed by the overburden strata are different along both strike direction and dip direction. As a result, the subsidence influence throughout the whole length of a longwall panel varies. Dramatic adjustment of strata structures may induce violent surface subsidence. Stability of strata structures within overburden before the final subsidence controls the stability of ground surface land. Chain pillars of 20 m between panels of 240 m wide with a cover depth of 600 m have been crushed in the gob and do not have any function in supporting the overburden strata.
e study also suggests that the final subsidence of the three adjacent panels is far to come in the future, and the land reuse above underground coal mines should be 14 Advances in Civil Engineering carefully planned by making sure that the gob is completely compacted or no potential secondary subsidence occurs in the future. It is hoped that the investigations of this study can provide some insight into taking positive measures for environmental protection and land reuse plan for postmining cities.

Data Availability
e data used to support the findings of this study are available upon request.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
Kai Zhang contributed to conceptualization and supervision; Lu Bai and Pengfei Wang performed formal analysis, investigation, and data curation; wrote the original draft; provided the software; and reviewed and edited the article; Zhuang Zhu checked typos and grammatical errors; Kai Zhang contributed to funding acquisition.