Research on Surface Failure Law of Working Faces in Large Mining Height and Shallow Buried Coal Seam

In the process of high-intensity and large-space mining in Shendong mining area, various surface cracks are generated on the surface, resulting in serious damage to the surface buildings and the local ecological environment. To study the influence of overlying rock movement on surface failure of near-field single key strata of near-shallow buried and large mining height working face, the relationship between overburden movement, strata pressure appearance, and surface failure at working face 52307 in Daliuta mining area was analyzed by fieldmeasurement and numerical simulation.+e results show the following: (1) there is only one thick and hard key stratum in the overburden of large mining height and near-shallow buried working face. Under the condition of presplitting roof blasting, the first weighting step is still as high as 95m, and the periodic breaking step of roof is 20–30m. During the weighting, the working resistance of support is still close to the rated resistance. (2) +e single key stratum plays an obvious role in controlling overburden movement. After the first weighting of the working face, a stepped subsidence crack appears on the surface within a short time, and the crack lags behind the working face for about 5m. (3) During each periodic weighting process, the breaking and subsidence of key blocks are accompanied by surface cracks.


Introduction
In the coal mining of working faces, the surface usually has a certain response. For example, surface soil strata can be seriously deformed, and mining cracks appear on the surface [1,2]. is kind of destruction is quite remarkable and poses a great threat to surface buildings, roads, and railways [3,4]. erefore, it is necessary to analyze and study the mechanism of surface cracks.
Many scholars have studied the development mechanism of mining-induced surface fracture. In 1981, Singh et al. [5] concluded that surface subsidence is caused by overburden fault movement. is is the beginning of a deep study on the surface and spatial-temporal evolution of overburden strata. Ghose [6], Greco [7], and Homoud [8] used fracture mechanics and plastic mechanics theory to study the mechanical mechanism of mining landslide. Luo [9] applied computer simulation technology to investigate the dynamic deformation calculation model of mining slope. e basic laws of surface cracks have also been studied. Herrera et al. [10] and Camec and Delacourt [11] monitored land subsidence through D-InSAR and detected tensile deformation area. Hahn et al. [12] used the equipment to investigate and classify surface cracks.
With the westward movement of China's coal mining, large-space and high-intensity mining of shallow and thin bedrock coal seam has caused irreversible impact on the surface and local ecology under the special geological conditions in Northwest China [13,14]. A lot of research on the mechanism and development of mining-induced cracks in western mining area has been performed [15,16]. Liu et al. [17] analyzed the mechanism and the dynamic development law of surface cracks and surface subsidence.
is study provides technical references for the placement of surface cracks. e movement of overlying strata in shallow seam mining is characterized by the full cutoff of roof strata and the large bedrock fracture angle; and the fracture directly develops to the surface, forming significant subsidence cracks [18]. In addition, under the condition of high-intensity mining, the surface cracks have secondary development characteristics, and there are two development cycles and corresponding two peaks [19,20]. At the same time, the development mechanism of surface cracks under the condition of wide and fast advance working face was analyzed in detail [21,22]. Hu et al. [23,24] focused on indepth study of surface subsidence and surface cracks in aeolian sand area and proposed suggestions on the restoration of ecological environment after the destruction of aeolian sand mining. is study provides a basis for ecological restoration at that time. Sun et al. [25] put forward a "quasi-hyperbolic" settlement model on the basis of key strata theory and basic mechanics. is study provides a new theoretical basis for the study of surface subsidence and damage. Ju and Xu. [26] studied the spatial and temporal effects of different failure structures in different locations of key strata on surface damage by changing the location of key strata through similar simulation methods. Yan et al. [27], Yang et al. [28], and Zhou et al. [29] performed a lot of research on the occurrence and disaster causing of surface cracks in shallow coal seam mining and provided favorable measures for the restoration of surface ecological environment and land reclamation. Other studies on the mechanism and development characteristics of surface cracks in shallow seam mining in Western China were conducted [30][31][32], contributing to local ecological protection and safe mining.
However, the above studies are focused on the shallow seam mining [33,34]. ere are few studies on the mechanism and development of surface cracks in the near-shallow seam mining (coal seam depth > 150 m) [35]. Based on the geological environment of working face 52307 with shallow burial and high mining height in the Daliuta mining area, the spatial-temporal relationship between surface damage and breakage of key strata is explored in this study. In this paper, the generation and development of surface cracks under the condition of near-shallow burial is explained. is study provides a theoretical basis for the research on mininginduced cracks.

Geological Conditions of Mining Working Face.
Daliuta Coal Mine is located in Wulanmulun Town, Shenmu County, Yulin City, Shaanxi Province. With the border to Ordos city of Inner Mongolia autonomous region, this mine belongs to Shendong Coal Group. With the subordinate of Daliuta Well and Huojitu, this mine is a supersize, highyield, and high-efficiency mine in China, with an annual output of more than 20 million tons. Figure 1 shows the geographic location of the Daliuta Coal Mine.
Working face 52307 is a fully mechanized mining face with large mining height in the Daliuta mine, the ground elevation is +1120.2-1217.1, and the coal seam floor elevation is +987.7-1025.5. Most of the surface overburden strata are Quaternary loose sediments, which is 32 m in working face 52307, and the thickness of overlying bedrock is 158 m. e strike length of the working face is 4462.6 m and the inclination length is 301 m. e roadway layout of the working face is shown in Figure 2(a). Coal 5-2 is the main working face with an average coal thickness of 7.2 m, an average coal seam depth of 190 m, and there is a stratum of giant fine-grained sandstone with a thickness of nearly 31 m at 4 m above the coal seam, as shown in Figure 2 In the advance process of working face, various types of surface cracks appear, as shown in Figure 3. Surface cracks in Figure 3 are placed in other coal mine nearby the Daliuta mining area. Since the geological and mining conditions in coal mining area nearby the Daliuta mine are similar, the characteristics and mechanism of surface cracks caused by underground mining are the same.

Discrimination for Key Strata Location.
e location of key strata is judged by two criteria of critical strata [36], namely, stiffness criterion and strength criterion. e criterion of stiffness is where q n+1 and q n are the loads on the key strata of the first strata when calculated to n + 1 strata and n strata, respectively. To judge whether n + 1 strata can be the key strata, the strength criterion should be satisfied, namely, the breakage distance of the lower strata should be less than that of the upper strata, shown as follows: where l n is the breakage distance of the n-th strata and k is the determined number of hard rock strata in equation (1). Based on the results of stiffness and strength discrimination of key strata and field settlement measurement of  overburden, fine-grained sandstone with a thickness of about 31 m in working face 52307 is determined as the main key strata in overburden strata of working face. Table 1 shows parameters of the roof and floor strata of working face.

Analysis of Weighting Law of Working Face
e mining pressure data of working face 52307 were collected by the monitored research group in the initial mining stage. Six Uloka pressure gauges are arranged at the head of the working face machine, which are located at 39#, 40#, 42#, 43#, and 44# brackets, respectively. Eight Uloka pressure gauges are arranged at the middle of the working face, which are located at 70#, 71#, 72# (the center of the working face), 73#, 74#, 75#, and 76# brackets, respectively. Five Uloka pressure gauges are arranged at the tail of the working face machine, located at 114#, 115#, 116#, 117#, and 118# brackets. Uloka data are collected once every 1 s. Figure 4 shows the layout of extension.   In this paper, the bracket in the middle of the working face is selected to analyze the pressure data of the bracket. Figure 5 shows the working resistance of the bracket. e formula of weighting criterion is where p t is the mean end-cycle resistance and σ pi is the mean square deviation. As shown in Figure 5, when the working face advances to 955 m, the support resistance increases rapidly, and the first weighting of the working face occurs, namely, the first weighting of the working face (when the cantilever of the main roof reaches its limit, the fracture of the main roof forms a three-hinged arch-like balance, while the broken rock blocks rotate unstably, resulting in a sharp drop of the working full roof. e load on the bracket is generally increased, and this phenomenon is called the initial weighting of the main roof ). e weighting step is about 95 m. When the working face advances to 120 m, the second weighting occurs, namely, the first periodic weighting of the working face (the roof weighting phenomenon caused by periodic instability of rock strata in fracture zone is called periodic weighting of the working face roof ), and the weighting step is 25 m. When the working face advances to 140 m, the second periodic weighting occurs, and the weighting step is 20 m; within the 204 m range of the working face advance, there are six types of obvious periodic weighting. e largest weighting step exceeds 30 m, averaging about 20 m. Under the special overburden condition of single thick and hard main key strata in the near field, presplitting blasting is carried out at the cutting end of the working face, but the distance between the initial mining weighting step and throwing step can reach 95 m, and the interval between periodic weighting step is also in the range of 20 m-30 m.
rough detailed analysis of data, it is concluded that the average end-cycle resistance of the bracket is 17111.04 kN during weighting and 10451.08 kN during nonweighting.

Relationship between Generation of Surface Cracks and Weighting of Working Face
In the initial mining stage of working face 52307 in the Daliuta Mining Area, there is a stratum of fine sandstone about 31 m thick above the working face. en this stratum is determined as the main key strata of the working face by the discriminant condition. e first weighting of the basic roof fracture occurs when the working face advances to 95 m. e step subsidence of 50 cm is generated near the direction of 70 m, and the subsidence of the surface measuring point is 746.7 mm. e relationship between underground periodic pressure and surface fracture is obtained, as shown in Table 2. e negative value in the table represents the distance between the ground fracture lagging behind the underground pressure point, while the positive value represents the distance between the ground fracture ahead of the underground pressure point. According to Table 2, when the pressure of the underground working face is applied, cracks will occur above the surface; if pressure is not applied on the working face, there are no cracks. It indicates that cracks can be generated in the working surface 52307 during nonweighting period. It further shows that the influence of overburden movement on surface deformation is intense.  ree different locations in the working face are selected to arrange the hydraulic support of the Euloka station to ensure the accuracy of data collection. e setup entry is the starting place of 52307 working face, as shown in Figure 6, during the initial weighting process of the working face, the "fixed supported beam" structure of the key bearing strata breaks in the middle, and the overburden subsides and forms a tensile force toward the middle. e overburden breaks and subsides on both sides of the corresponding surface and forms a tensile action to form the surface crack. e surface crack is shaped like a funnel in the "C" shape and the "C" direction (as shown in black box). In the latter weighting process, the periodic weighting step of the support in the middle of the working face is obviously smaller than that of the support at both ends, which results in unequal weighting step. e overburden rock in the middle of the working face has been broken, but the overburden rock at both ends has not been broken [37]. Besides, there is a "dislocation." e first fracture and subsidence of the middle roof overburden directly form a tensile force on the upper surface loess strata of this part, so that the surface cracks appear in the "C" shape during the period of periodic weighting.
In the face of the periodic weighting, when the key stratum of the middle support of the working face is broken, typical cracks on the surface are compared. According to the field observation, the surface is often accompanied by various kinds of cracks during each weighting, as shown in Figure 7.

Fracture Genesis Analysis.
With the advance of the working face, uneven settlement will occur on the surface. e horizontal displacement and subsidence of the surface soil in adjacent positions are different, which will inevitably lead to the deformation of the surface overburden soil. Because of small tensile and poor shear capacity of the Quaternary loose strata [38], surface cracks occur when the tensile deformation reaches the maximum tensile deformation energy of the surface soil.
To fully study the actual deformation of the surface, an observatory is set up on the surface above the open-off cut. e horizontal section of the observatory is shown in Figure 8. e actual length of the observation station is 626 m, the distance between the observation points is 20 m, and the number of working observation points is 23 m. By observing the surface movement and deformation of the working face and the calculation of observation data, the surface movement curve is obtained. e subsidence curve of the survey line in each observation period is shown in Figure 9. e average depth of working face 52307 is 190 m, and the length of working face is 301 m; the mobile basin with stable mining will be fully exploited. During the last observation, the working face 52307 is advanced to 637.4 m, the     According to Figures 9 and 10, in the advance of the working face, the surface subsides under the influence of mining, and the horizontal deformation of the surface soil changes with fluctuations. Near the measurement point A23, the maximum of surface deformation is accompanied by the maximum horizontal deformation. It indicates that the surface movement is intense, different forms of surface cracks will be formed in this area, and even multiple cracks will appear. In the actual observation, with the advance of working face, the speed of surface subsidence in the working face is fast, and the intense surface horizontal movement can be observed. e movement of the overburden subsidence  Advances in Civil Engineering results in a sharp subsidence of the topsoil strata. However, there exists the uneven subsidence because the subsidence speed between adjacent soils is not uniform in time. At this moment, the soil bears great tensile stress. When the tensile force is greater than the limit bearing force of the soil, dislocation is caused between the adjacent soils in the vertical direction, and surface cracks are generated.

Field Measurement of Surface Crack Morphology.
Surface cracks in the stations are recorded, and the characteristics of ground deformation and cracks are studied. e crack situation is shown in Figure 11. With the advance of the working face, surface cracks are generated. e distribution of cracks in the middle of the surface is basically parallel to the working face, some of which bend towards the goaf, and the overall shape of the cracks presents a "C" shape. At the beginning of mining in the working face, the surface dynamic cracks are mainly distributed in the middle of the working face. When the working face advances to about 450 m, fewer permanent cracks are found above the open-off cut, along the grooves on both sides of the working face and the adjacent areas. Meanwhile, surface subsidence is serious in the process of on-site investigation and mining, and there will be different forms of surface cracks in the field, as shown in Figure 12. Surface cracks caused by underground mining can be classified as tensile cracks, hysteretic extrusion crack, and stepped cracks. Tensile cracks are caused by tensile stress, hysteretic extrusion cracks are caused by compressive stress, and stepped cracks are caused by shearing stress in Quaternary loose sediments.
rough the analysis of the surface tensile cracks and the advance of the underground working face, it is known that the position of the cracks always lags behind the mining position of the working face for a certain distance. Stepped cracks are formed by mining-induced overburden breakage and surface subsidence, generally lagging behind the working face. eir mechanism can be divided into basic roof breaking principle based on thin plate theory and key strata theory. Based on the analysis of stepped cracks in working face 52307 of Daliuta, it is considered that the occurrence of stepped cracks is closely related to the lag distance of cracks, subsidence value, and horizontal deformation value of surface points of cracks. ese cracks are located in the tensile deformation area, rather than the extrusion area, and their deformation value and subsidence value are relatively large.
Traditionally, it is generally believed that dynamic cracks open first and then close gradually with the advance of the working face. e width of cracks generally presents a single peak period and increases from small to large and finally tends to close. Based on the statistics of all cracks in working face 52307, the development time of tensile, extrusion, and stepped cracks are analyzed, respectively. Tensile crack, which is located behind the working face in Daliuta, generally undergoes the process of "opening-large-closing." It has a small deformation value and short duration (averaging 4.8 days).  Advances in Civil Engineering crack, but the damage degree of extrusion crack is far less than the stepped crack. ey have small change range and rarely change after the generation, undergoing the process of "opening-unchanged." ese cracks generally lie behind the working face and last for an average of 5.3 days. e maximum distance pushed by the working face is 33.3 m, the minimum is 21.8 m, and the average is 29.5 m.

Model Design. CDEM software is the abbreviation for
Continuum-Based Discrete Element Method, which is the coupling of continuous medium algorithm (such as finite element) and discontinuous medium algorithm (such as discrete element), developed by Institute of Mechanics, Chinese Academy of Sciences. By controlling the block (finite element calculation in the block, discrete element calculation in the boundary), to simulate the deformation and movement characteristics of materials in discontinuous and continuous states, the failure process of materials can be realized. e numerical model of the hydraulic support is embedded in the simulation. According to the data of drilling strata in the initial mining stage of the working face, the model has 45 strata and the strike length is 600 m; Table 3 shows rock mechanics parameters in numerical simulation. In numerical simulation, the upper boundary is free face, and the bottom boundary is the vertical displacement restriction; the right and left boundary are both the horizontal displacement restriction. e constitutive model of hydraulic support adopts the numerical constitutive model which characterizes the relationship between the support shrinkage and the working resistance of the support. According to the calculation principle and characteristics of the software, to eliminate the influence of boundary effect, 100 m boundary pillars are left on both sides, the actual mining length is 400 m, and the simulated working face advances 400 m along the openoff cut. e establishment of surface morphology is based on the field data such as the contrast map, the layout map of surface survey lines, as shown in Figure 13. ere are 9065 block units in the model. According to the different lithology and thickness of each stratum, the block size of each stratum is different, and the characteristics of overburden breakage in situ can be restored to the greatest extent. And the plane numerical model is suitable for reflecting the relationship between underground mining and surface cracks. Because the main key stratum belongs to the near-field key strata, according to the theoretical formula, the stratified breakage supplementary caving zone is generated at the thickness of 16 m, and the remaining 15 m is the basic top of the key strata, that is, the key bearing strata of the whole overburden. When the working face advances 95 m, the "fixed beam" on the basic top of the key strata breaks the first weighting of the working face, and then surface cracks occur, as shown in Figure 14.
Due to the existence of tensile force, there are two obvious surface cracks on both sides above the open-off cut, the support and the coal wall, respectively. e cracks above the open-off cut are permanent cracks, which will not close with the advance of working face; the width of the cracks above the working face and the coal wall is between 0.5 m and 1 m. e penetration depth of the crack is about 5 m, and the crack hysteresis working face is 10 m. e first periodic weighting occurs when the working face advances to 115 m. e tensile cracks formed in the upper stage are closed due to the extrusion force behind them, while the breaking and subsidence of the key bearing blocks lead to the generation of new surface cracks.
According to Figure 15(a), two surface cracks are formed in this stage, both of which are hysteretic surface cracks.
e first crack lags the working face within the range of 6-9 m, the second surface crack lags the working face is about 25 m. Due to the different distance from the side of goaf, the farther the distance, the greater the tension effect. us, the penetration depth of the first crack is greater than that of the second one. ere is little difference in the lag distance of cracks at this stage in the simulation and in the field. e second periodic weighting occurs when the working face advances 130 m, as shown in Figure 15(b).
e key blocks of the key bearing strata break down and subside, and the surface cracks are very obvious. e crack lags behind the working face between 12 m and 15 m, and the width of the crack opening is about 1 m, which is larger than that of the first periodic weighting and deeper penetration into the bedrock.
Because there is only one key strata structure in the working face, and it belongs to the near-field key strata,    Advances in Civil Engineering when the key block of the key bearing strata breaks and has the subsidence displacement, the overburden bedrock will show the breakage and the subsidence displacement simultaneously. As a result, the surface also shows subsidence displacement, and the working face the first weighting. Within the initial weighting range of the working face, the key strata basic top break first subsides displacement, the overlying strata above follow, and the surface also has a certain subsidence at this time.
In the numerical simulation, the measured points are arranged along the surface of the model, and the surface subsidence at each step of excavation is monitored in real time. Figure 16 shows the surface subsidence curve.
As shown in Figure 16, when the initial weighting range of the working face is pushed to 100 m, the surface has already subsided. When the working face advances to 100, the maximum subsidence of surface survey points is 571 mm. When the working face advances 400 m, that is, when the excavation of the model is completed, the final subsidence of the surface survey line is 4,175 mm, and the surface subsidence value is 3,560 mm in the field measurement. Although there are some errors in the field measurement, the errors are unavoidable due to the difference between the numerical simulation and the rock strata change in the field.

Conclusions
(1) By distinguishing the overburden lithology and key strata of working face 52307 with shallow burial and high mining height in the Daliuta mining area, it is concluded that fine-grained sandstone with a thickness of 31 m or so above the coal seam is the main key strata of the overburden of the working face. According to the analysis of support pressure data of working face, after a roof prefracturing project, the first weighting step of working face still reaches about 95 m, and the periodic weighting step is in the range of 30 m-20 m, which indicates that the weighting step of working face is larger. (2) rough the analysis of statistical results of the observation stations on the surface above the working face, after fully mining, the surface survey lines are observed; the maximum subsidence of the surface is 3560 mm from the last observation, because the key strata can control the whole overburden movement from the key strata to the surface. At the same time, the horizontal deformation curve also shows that the surface deformation is serious along with the advance of the working face. According to the statistics of the cracks, with the  breaking movement of the overburden, various types of surface cracks appear on the surface, mainly tensile ground cracks, which lag behind the working face.
(3) rough the validation of numerical simulation, with the advance of the working face, the key bearing strata breaks and subsides, and the overburden bedrock moves as a whole, resulting in obvious lagging surface cracks.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.