Next Article in Journal
Three General Double-Series Identities and Associated Reduction Formulas and Fractional Calculus
Next Article in Special Issue
Study on Uniaxial Compression Deformation and Fracture Development Characteristics of Weak Interlayer Coal–Rock Combination
Previous Article in Journal
Abstract Fractional Cauchy Problem: Existence of Propagators and Inhomogeneous Solution Representation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractal Characteristics of Overburden Rock Fractures and Their Impact on Ground Fissures in Longwall Coal Mining

1
School of Emergency Management and Safety Engineering, China University of Mining and Technology, Beijing 100083, China
2
Inner Mongolia Research Institute, China University of Mining and Technology (Beijing), Ordos 017001, China
3
State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines, Anhui University of Science and Technology, Huainan 232001, China
4
Information Institute of the Ministry of Emergency Management of PRC, Beijing 100029, China
5
Coal Mining and Designing Branch, China Coal Research Institute, Beijing 100013, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(10), 699; https://doi.org/10.3390/fractalfract7100699
Submission received: 7 August 2023 / Revised: 8 September 2023 / Accepted: 15 September 2023 / Published: 23 September 2023
(This article belongs to the Special Issue Applications of Fractal Analysis in Underground Engineering)

Abstract

:
Ground fissures are generated during the coal mining process due to overlying strata migration, which provides gas and water seepage channels and usually contributes to coal mining accidents in shallow buried coal seams with larger mining height working faces. Thus, the evolution features of ground fissures in large mining height working faces and shallow buried coal seams were explored by considering field observation data, similar simulation, and numerical simulations. The results show that the weathered rock layer above the thin bedrock inhibits fracture growth caused by coal mining. Overlying strata fracture expansion can be quantitatively divided into three stages based on fractal dimension: the overburden fracture formation stage, the overburden deformation fracture expansion stage, and the overburden fracture stable development stage. The movement deformation region of the ground fissure can be described by three characteristic zones: the boundary tension zone, the central compression zone, and the central dynamic tension and compression zone. Central dynamic fissures usually have 3–5 fissure counts with 2–5 m intervals between each group of fissures; central fissures develop twice from open to closed widths; the time period for two central fissures to reach their with maximum widths is 11–20 days; and border fissures expand rapidly with maximum values in a time period of 5–6 days.

1. Introduction

Coal is one of the basic energy sources. About 80% of coal mining in China is underground mining [1,2]. However, with the rapid development of the energy economy in China, the demand for coal energy has been increasing in the past 20 years, and the ecological and environmental damage caused by underground coal mining poses a threat to the sustainable development of society in mining areas [3,4]. In western mining areas of China, about 70% of all the raw coal is produced where the depth of the coal seam is generally between 50 and 250 m, the bedrock is relatively thin, the loose layer and coal seam are both thick, and the mining intensity is higher [1,5,6,7,8]. It is difficult to maintain the stability of the main key layer for a long time after the roofs of the coal seam are broken [2,7,9,10,11,12,13]. The uncoordinated movement of the upper strata affects the discontinuous settlement of the ground, which can cause collapse grooves, permanent ground fissures, and boundary ground fissures to appear on the ground. Some ground fissures induced by mining run through to the goafs, which may cause water bursting, sand bursting, and coal spontaneous combustion disasters in the subsequent mining process [1,8,11,14,15,16,17,18].
Underground coal seam mining causes the top plate to break and the overlying strata to subside. Eventually, it causes continuous and discontinuous ground deformation, resulting in ground cracks, step settlement, and collapse pits. At present, the analysis methods used are mainly empirical formula methods, theoretical analysis methods, simulation methods, and field measurements.
At present, the methods used to predict mining subsidence are based on the analysis of ground observation data, which cannot explain the damage caused by overlying strata migration. However, mining subsidence is caused by overlying strata breaking from the coal seam to the ground. The height of the broken roof and backfilling affects the bearing capacity transfer rule of the overlying strata, and finally causes the movement law of ground deformation [1,5,7,12,19,20]. Based on field observation data combined with physical similarity simulation experiments, fractal theory can be applied to explain the movement characteristics and fracture extension of the overlying strata [21,22]. The dynamic evolution of the overlying strata fractures can be analyzed and calculated by numerical simulation software. However, the use of numerical simulations to simply interpret the characteristics in the actual mining process needs to be verified by physical modeling or field engineering measurements [23,24]. The movement of the overlying strata and the damage to and deformation of the strata are discontinuous deformation problems in the longwall high-intensity mining of shallow coal seams. Based on continuous medium finite element numerical simulation analysis, the dynamic development of fractures cannot be accurately described [25,26]. In contrast, the most important advantage of PFC2D is the discrete element particle flow program, which can simulate the overlying strata migration processes of fracture, collapse, accumulation, and compaction [27,28,29].
In summary, this paper studies the evolution law of overlying rock fractures caused by longwall mining, by combining field examples and numerical analysis. Based on the field examples, the monitoring of surface fracture generation and development caused by overlying rock migration in the process of mining is analyzed. Using the PFC2D (Particle Flow code 2D) numerical analysis method, the characteristics of surface fractures and stress change induced by overburden fractures during coal excavation are simulated and analyzed. The research results have theoretical and practical significance for coal mining and environmental and ecological management in coal mining.

2. Geological Conditions

The No. 12401 panel of the No. 1–2 coal seam has a designed advanced length of 5254.8 m, width of 299 m, and height of 8.8 m. The ground elevation varies from 1188 m to 1300 m, with a coal seam floor elevation range of 1043–1066 m. The loose overburden is 0 to 27 m thick and consists mainly of loess and fine sand. The thickness of the overlying bedrock is 120–220 m. The false roof of the 12401 working face is mudstone with a compressive strength of 11.3–13.2 MPa and consolidation coefficient of 1.32. The immediate roof is gray fine sandstone with a compressive strength of 13.3–15.2 MPa and consolidation coefficient of 1.35. The main roof is grayish-white siltstone with a compressive strength of 14.5–36.6 MPa and consolidation coefficient of 2.32. The main roof is silty sandstone with a compressive strength of 14.5–36.6 MPa. As shown in Figure 1, there are three types of ground fissures in the Shendong coal field.

3. Overlying Strata Characteristics of the Working Face

Figure 2 shows the borehole column diagram of the 12401 working face. Ground boreholes R113, b279, and R106 on the outer side of the transport roadway and R110 on the outer side of the return laneway in working face No. 12401 are counted. The lithology statistics based on the surface boreholes around the 12401 working face in the Shangwan Mine are shown in Table 1. The hard rock strata above the main mining seam mainly consists of sandstone (siltstone, medium-grained sandstone, and fine sandstone) and soft rocks (sandy mudstone).
Based on the bedrock lithology analysis of the borehole column diagram near the work area, the bedrock has mudstone at the bottom and top. The rock strength is generally reduced by the oxidation zone. The clay content of the rocks increases relatively, and plasticity generally increases. This inhibits the upward development of overlying strata fractures to a certain extent. Therefore, the upper overlying strata oxide layer has an inhibitory effect on the overlying strata fracture in the shallow bedrock of working face No. 12401.
As shown in Figure 3, Figure 4 and Figure 5, the thicknesses of the bedrock and loose layers were counted by drilling around the working face of the 12401 working face. Contour plots of the bedrock and loose layers and the bedrock thickness ratio were plotted using Surfer 10.0 software. The analysis shows that the thickness of the loose layer above the working face gradually decreases from the air return lane to near the transport lane. The thickness of the bedrock gradually decreases along the working face advance direction. The ratio of the thickness of the loose layer above the working face to that of the bedrock gradually decreases from the air return lane to near the transport lane.

4. Simulation Experimental Study

4.1. Simulation Experimental Analysis of Dynamic Fractures

As shown in Figure 6, the size of the model frame was chosen as 1800 mm × 200 mm × 1200 mm (length × width × height). According to similarity theory [11], the following several types of similarities should be considered.
For geometric similarity,
L p L m = C L
where Lp is the size of the on-site prototype, Lm is the size of similar models, and CL is the size similarity ratio.
For dynamic similarity,
γ p γ m = C r
where γp is the bulk density of the on-site prototype, γm is the bulk density of similar models, and is the similarity ratio of bulk density.
For time similarity,
T p T m = C T
where Tp is the on-site prototype time, Tm is the similarity model time, and CT is the time similarity ratio.
The geometric similarity ratio is 200:1, the volumetric similarity ratio is 1.5:1, and the model time ratio is 96:1. The simulated materials are river sand as aggregate, calcium carbonate and gypsum as cement, mica as an auxiliary material, and a certain proportion of fly ash that is added to the coal seam.
To analyze the dynamic development law of the overlying strata transport and fractures, the fracture extension height and rock failure area of the overlying strata were determined by similar simulation. Figure 7 shows different advance distances of the rock failure area.

4.2. Simulation of the Caving Process of the 12401 Working Face Based on PFC

In terms of numerical simulation, we simplified the complex geological conditions without too much addition to the temperature field and seepage flow field. Based on the geometry of the Shangwan No. 12401 working face, a two-dimensional model based on PFC2D 5.0 software was established in order to further validate the simulated fracture evolution characteristics and laws of similar simulations, revealing the evolution characteristics of the overlying rock fractures at different mining stages. As shown in Figure 8, the model is 300 m long and 200 m high. The upper boundary is free, the left and right boundaries fix the x-displacement, and the bottom boundary fixes the y-displacement. The particles contact each other through parallel bonding, the particle diameter ranges from 0.3 to 0.5 m, and a total of 102,814 balls are generated. The model is used to simulate the longwall mining process, with advancement from left to right and a mining height of 8.8 m. Compression and tensile tests were carried out using PFC2D 5.0 software. The relationship between rock macroscopic parameters and rock fine-grained parameters was inverted to find the desired fine-grained parameters by inversion of the stratigraphic macroscopic parameters [27,30], as listed in Table 2.
As the working face advances, the overlying strata collapse and migrate, causing a new fracture extension, which eventually leads to overlying strata step settlement. As shown in Figure 9, the migration laws of the overlying strata from similar simulation and numerical simulations are consistent. From the force chain diagram, the overlying strata fracture extension height was 65 m when the working face advanced to 150 m. When the working face advanced to 200 m, the overlying strata fracture extension height was 160 m. When the working face advanced to 250 m, the overlying strata fracture height extended to the ground surface. The equilibrium structure of the entire overlying strata broke, and the overlying strata broke up into blocks. The settlement of the blocks eventually leads to ground deformation. From the horizontal displacement-x contour, it can be concluded that the horizontal displacement of the ground surface varies unevenly as the working face advances. Eventually, ground tension fissures, compression fissures, and step subsidence fissures form.

5. Results and Discussion

5.1. Overlying Strata Fracture Evolution Characteristics

As shown in Figure 10, a high-speed camera was used to record and photograph the transport characteristics of the overlying strata fracture development characteristics after each advance of the working face. Then, Photoshop and ImageJ image processing software were used to extract the range of fracture development characteristics after each advance of the working face. Finally, the fractal dimensions of the fracture network map were calculated using the box dimension method. First, an image of the main study area is intercepted. Then, the appropriate grayscale and contrast are adjusted to identify the difference between the fracture distribution and nonfracture pixels. The graphics of the study area are imported into ImageJ 1.0 software to derive the fractal dimensions of the fractal network, as shown in Table 3.
The characteristic of fractal dimension is scale invariance, which is universal to all kinds of scales. Fractal dimension can quantitatively describe the complex physical processes in mining engineering and reveal the internal mechanism of these physical processes. Prior research shows that the evolution of overburden fractures is a complex physical and mechanical process, during which many physical and mechanical parameters such as roof stress and deformation will change, and the magnitude and even magnitude of different production geological conditions are different, which make it difficult to use as an accurate prediction index. Therefore, fractal dimension can explain the complex fracture network and quantitatively characterize the evolution law of overcast rock fractures at each stage in the mining process, providing a reference for safe production and ecological restoration and management after mining. As shown in Figure 11, an overlying strata fracture was formed when the working face advanced to 20–60 cm, the fractal dimension was 0–1.325, the extension height of the overlying strata fracture was 20 cm, the rock damage area was 1485 cm2, and the height of the overlying strata fracture was 0.9 times the advance distance. When the working face advanced to 60–135 cm, the fractal dimension was 1.325–1.38, the extension height of the crack was 88 cm, the rock damage area was 7949 cm2, and the height of the crack in the overlying strata was 0.07 times the advance distance. This is the stage of deformation and fracture extension of the overlying strata. When the working face advanced to 142.5 cm, the rock damage area was 10,169 cm2. With the advance of the working face, the three stages of overburden fracture development were divided mainly through the upward–stable–downward trend of the fractal dimension of the overburden fracture. The overlying strata fracture height extension pattern was divided into three stages: the overlying strata fracture formation stage, the overlying strata deformation fracture extension stage, and the overlying strata fracture extension steady-stage. These three stages can better describe the disturbance of the overlying strata and the initial formation of fractures in the early stage of workface advancement. This method can quantitatively explain the evolution law of overburden fractures in the actual mining process of the coal mine.

5.2. Ground Fissure Movement and Deformation Zone

In the early stage of mining, the main roof stratum is bent and sliding, and the ground has no obvious deformation characteristics. With the continuous advance of the working face, the overlying strata will eventually have stepped fault settlement, and the overlying strata fault will penetrate the ground from bottom to bottom [31]. Eventually, boundary fissure zones and central fissure zones form. In the study of ground fissure dynamics, ground fissures begin to crack, expand, contract, stabilize, or close. Ground fissures close due to damage caused by gravity within the soil, resulting in central dynamic ground fissures experiencing a double cycle extension [2,32]. In order to better simulate the morphological changes of surface cracks caused by fracture expansion in the overlying strata, the development and expansion laws of overlying strata cracks in different mining stages were simulated by PFC. As can be seen from the force chain diagram, the entire equilibrium structure of the overlying strata has been disrupted. The settlement of the block ultimately leads to ground deformation. It can be seen from the horizontal displacement-x contour that the horizontal displacement of the surface changes unevenly. Finally, ground cracks, extrusion cracks, and step settlement cracks are formed. As shown in Figure 12, through physical simulation experimental and numerical simulation studies of the overlying strata migration, the main influencing factor is the key strata in the overlying strata. As the working face progressed, key strata were broken, leading to massive settlement displacement of the overlying strata and the development of ground fissures. The rock blocks in the goaf are irregular, the arrangement is also very irregular, and the loose coefficient is relatively large. However, the overlying rock strata are broken and sink again, causing surface cracks to develop again, ultimately forming stable ground fissures. Eventually, from the horizontal displacement-x contour, it can be concluded that the horizontal displacement of the ground surface varies unevenly as the working face advances. The ground fissure compressed zone, tensile zone, compression zone, and central fissure zone form.

5.3. Field Ground Fissure Monitoring Analysis

In previous studies, there were few continuous settlement observations, and the time difference was more than three days. However, the selected ground fissure was continuously observed for monitoring. The Shangwan mine studied in this paper is different from other mines. The characteristics of the Shangwan mine are a large mining height, fast advancing speed, shallow burial depth, and large area of overlying rock caving during mining, so that cracks develop in the overlying rock and penetrate the surface. However, differently from the observed phenomena of fractures in other mining areas, the ground fissures in the middle of the working face are a secondary periodic development in the Shangwan mining area. Due to the complex environment of the field, we could collect only limited data, which provided basic research for this paper. The ground fissure development monitoring data of the working face were counted and listed in Table 4, as shown in Figure 13a. With the advancement of the working face, the f14 ground fissure forms when the working face advances 360 m, and the f08 ground fissure forms when the working face advances 450 m. From 3 June to 7 June, the f14 ground fissure width expansion is 55 mm, which then decreases to 31 mm. From 13 June to 18 June, the f14 ground fissure width again expands to 39 mm, and then gradually decreases to 20 mm to maintain a stable width after 18 June. From 11 June to 15 June, the f08 ground fissure width expansion is 45 mm, and then the fissure width decreases to 24 mm from 15 June to 19 June. From 19 June to 25 June, the f08 ground fissure width again expands to 35 mm, and then gradually decreases to 18 mm to maintain a stable width after 25 June. As shown in Figure 13b, the hf4 ground fissure forms when the working face advances 360 m, and the hf5 ground fissure forms when the working face advances 375 m. From 3 June to 7 June, the ground fissure width expansion of hf4 is 2.8 mm, which gradually decreases to 1.5 mm from 7 June to 17 June, and maintains a stable width after 18 June. From 5 June to 9 June, the ground fissure width expansion of hf5 is 3.5 mm, which gradually decreases to 1.8 mm from 9 June to 15 June, and a stable width is maintained after 16 June.
As shown in Figure 13c, with the advancement of the working face, the hf3 ground fissure forms when the working face advances 360 m, and the hf2 ground fissure forms when the working face advances 465 m. From 3 June to June, the hf3 ground fissure width expansion is 42 mm, which gradually decreases to 25 mm from 9 June to 15 June, and maintains a stable width after 16 June. From 11 June to 16 June, the ground fissure width expansion of hf2 is 38 mm, which gradually decreases to 22 mm from 16 July to 23 July, and maintains a stable width after 24 June. The available findings indicate that periodic fractures in the overlying strata cause the periodic development of central ground fissures, which is confirmed by observations made under the same geological and mining conditions. With the advancement of the working face, the central ground fissures f14 and f08 are formed by a double cycle development process, while the boundary ground fissures hf2, hf4, hf3, and hf5 are formed by single cycle development processes.
Once the boundary fissure is formed, the opening and drop also expand rapidly, reaching a maximum within 3–4 days and finally close slowly to a certain width. The central ground fissures form in groups of 3–5, and the ground fissures are spaced 2–5 m apart. This indicates that the instability of the main key layer is the direct cause of ground settlement. The central ground fissure developed twice from opening to narrowing its width, and the time interval between when the maximum widths of the two central ground fissures developed was 11–20 days. The width of the central ground fissure that developed first is greater than the width of the second ground fissure. Boundary fissures develop only once, the fissures expand rapidly to the maximum width in 5–6 days, and finally the width of the fissures remains stable. In conclusion, once a boundary fissure forms, the opening and fall also expand rapidly, reaching the maximum value within 3–5 days, and the boundary fissure finally shrinks slowly to a certain width and remains stable. The central ground fissures form in groups of 3–5. The spacing between the ground fissures is 2–5 m. This indicates that the instability of the main key strata is the direct cause of ground settlement.

6. Conclusions

According to the engineering background of the No. 12401 working face of the Shangwan mine, similar simulation, numerical simulation, and fractal theory are used to study the evolution law of overlying rock fractures caused by longwall mining. This analysis and research method can be extended to other mining areas to study the evolution law of overlying rock fractures.
(1)
Laboratory experiments show the overlying strata fracture expansion of the working face is divided into three stages: the overlying strata fracture formation stage, the overlying strata deformation fracture expansion stage, and the overlying strata fracture stable extension stage. The overlying strata are broken into blocks, resulting in ground aeolian sand movement and deformation with the blocks. Finally, ground closure fissures, compression fissures, and open fissures formed;
(2)
Through physical simulation experimental and numerical simulation studies of the overlying strata migration, as the working face progressed, key strata were broken, leading to massive settlement displacement of the overlying strata and the development of ground fissures. From the horizontal displacement-x contour, it can be concluded that the horizontal displacement of the ground surface varies unevenly as the working face advances. Eventually, the ground fissure compression zone, tension zone, compression zone, and central fissure zone formed;
(3)
From the fieldwork, 3–5 central dynamic ground fissures formed in each group, and their spacing is 2–5 m. The central ground fissure developed twice from opening to narrowing its width, and the time interval between when the two central ground fissures developed maximum widths was 11–20 days. The width of the first central ground fissure to develop is greater than the width of the second ground fissure. Boundary fissures develop only once and expand rapidly to the maximum width in 5–6 days, and finally the width of the fissures remains stable. This indicates that in a western shallow buried mine, after mining a working face with a large mining height, the main key stratum cycle fracture leads to two central ground fissure development cycles.

Author Contributions

Conceptualization, C.L.; methodology, B.L.; formal analysis, C.L.; investigation, C.Z.; resources, T.T.; data curation, K.Z. and B.S.; visualization, J.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Nos. U1910206, 51874312, 51861145403), the major scientific and technological innovation project of Shandong Province (No. 2019SDZY01), and the National Science Fund for Distinguished Young Scholars (No. 52225402). These sources of support are gratefully acknowledged.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, L.; Wu, K.; Hu, Z.; Xu, Y.; Zhou, D. Analysis of Developmental Features and Causes of the Ground Cracks Induced by Oversized Working Face Mining in an Aeolian Sand Area. Environ. Earth Sci. 2017, 76, 135. [Google Scholar] [CrossRef]
  2. He, X.; Zhao, Y.; Zhang, C.; Han, P. A Model to Estimate the Height of the Water-Conducting Fracture Zone for Longwall Panels in Western China. Mine Water Environ. 2020, 39, 823–838. [Google Scholar] [CrossRef]
  3. Palchik, V. Formation of Fractured Zones in Overburden Due to Longwall Mining. Environ. Geol. 2003, 44, 28–38. [Google Scholar] [CrossRef]
  4. Sellers, E.J.; Kataka, M.O.; Linzer, L.M. Source Parameters of Acoustic Emission Events and Scaling with Mining-Induced Seismicity. J. Geophys. Res. Solid Earth 2003, 108, 2418. [Google Scholar] [CrossRef]
  5. Gao, R.; Kuang, T.; Zhang, Y.; Zhang, W.; Quan, C. Controlling Mine Pressure by Subjecting High-Level Hard Rock Strata to Ground Fracturing. Int. J. Coal Sci. Technol. 2021, 8, 1336–1350. [Google Scholar] [CrossRef]
  6. Lechner, A.M.; Baumgartl, T.; Matthew, P.; Glenn, V. The Impact of Underground Longwall Mining on Prime Agricultural Land: A Review and Research Agenda. Land Degrad. Dev. 2016, 27, 1650–1663. [Google Scholar] [CrossRef]
  7. Zhang, D.; Fan, G.; Liu, Y.; Ma, L. Field Trials of Aquifer Protection in Longwall Mining of Shallow Coal Seams in China. Int. J. Rock Mech. Min. Sci. 2010, 47, 908–914. [Google Scholar] [CrossRef]
  8. Xiao, W.; Hu, Z.; Fu, Y. Zoning of Land Reclamation in Coal Mining Area and New Progresses for the Past 10 Years. Int. J. Coal Sci. Technol. 2014, 1, 177–183. [Google Scholar] [CrossRef]
  9. Zhang, C.; Mitra, R.; Oh, J.; Hebblewhite, B. Analysis of Mining-Induced Valley Closure Movements. Rock Mech. Rock Eng. 2016, 49, 1923–1941. [Google Scholar] [CrossRef]
  10. Hu, Z.; Chen, C.; Xiao, W.; Wang, X.; Gao, M. Surface Movement and Deformation Characteristics Due to High-Intensive Coal Mining in the Windy and Sandy Region. Int. J. Coal Sci. Technol. 2016, 3, 339–348. [Google Scholar] [CrossRef]
  11. Yan, W.; Dai, H.; Chen, J. Surface Crack and Sand Inrush Disaster Induced by High-Strength Mining: Example from the Shendong Coal Field, China. Geosci. J. 2018, 22, 347–357. [Google Scholar] [CrossRef]
  12. Chen, C.; Hu, Z.; Wang, J.; Jia, J. Dynamic Surface Subsidence Characteristics Due to Super-Large Working Face in Fragile-Ecological Mining Areas: A Case Study in Shendong Coalfield, China. Adv. Civ. Eng. 2019, 2019, 12–17. [Google Scholar] [CrossRef]
  13. Zhuo, H.; Qin, B.; Shi, Q.; Li, L. Development Law of Air Leakage Fractures in Shallow Coal Seams: A Case Study in the Shendong Coalfield of China. Environ. Earth Sci. 2018, 77, 772. [Google Scholar] [CrossRef]
  14. Liu, C.; Li, H.; Mitri, H. Effect of Strata Conditions on Shield Pressure and Surface Subsidence at a Longwall Top Coal Caving Working Face. Rock Mech. Rock Eng. 2019, 52, 1523–1537. [Google Scholar] [CrossRef]
  15. Wang, J.; Park, H.D. Coal Mining above a Confined Aquifer. Int. J. Rock Mech. Min. Sci. 2003, 40, 537–551. [Google Scholar] [CrossRef]
  16. Guney, A.; Gul, M. Analysis of Surface Subsidence Due to Longwall Mining under Weak Geological Conditions: Turgut Basin of Yatağan-Muğla (Turkey) Case Study. Int. J. Min. Reclam. Environ. 2019, 33, 445–461. [Google Scholar] [CrossRef]
  17. Yang, H.; Liu, Z.; Zhu, D.; Yang, W.; Zhao, D.; Wang, W. Study on the Fractal Characteristics of Coal Body Fissure Development and the Law of Coalbed Methane Migration of around the Stope. Geofluids 2020, 2020, 9856904. [Google Scholar] [CrossRef]
  18. Xu, Y.; Wu, K.; Bai, Z.; Hu, Z. Theoretical Analysis of the Secondary Development of Mining-Induced Surface Cracks in the Ordos Region. Environ. Earth Sci. 2017, 76, 703. [Google Scholar] [CrossRef]
  19. Lian, X.; Hu, H.; Li, T.; Hu, D. Main Geological and Mining Factors Affecting Ground Cracks Induced by Underground Coal Mining in Shanxi Province, China. Int. J. Coal Sci. Technol. 2020, 7, 362–370. [Google Scholar] [CrossRef]
  20. Yang, Z.L. Stability of Nearly Horizontal Roof Strata in Shallow Seam Longwall Mining. Int. J. Rock Mech. Min. Sci. 2010, 47, 672–677. [Google Scholar] [CrossRef]
  21. Yang, B.; Liu, Y. Application of Fractals to Evaluate Fractures of Rock Due to Mining. Fractal Fract. 2022, 6, 96. [Google Scholar] [CrossRef]
  22. Wang, C.; Zhang, N.; Han, Y.; Xiong, Z.; Qian, D. Experiment Research on Overburden Mining-Induced Fracture Evolution and Its Fractal Characteristics in Ascending Mining. Arab. J. Geosci. 2015, 8, 13–21. [Google Scholar] [CrossRef]
  23. Lei, Q. Rock Mechanics Bulletin More Is Different: On the Emergence of Collective Phenomena in Fractured Rocks. Rock Mech. Bull. 2023, 2, 100080. [Google Scholar] [CrossRef]
  24. Zhang, C.; Bai, Q.; Han, P.; Wang, L.; Wang, X.; Wang, F. Strength Weakening and Its Micromechanism in Water–Rock Interaction, a Short Review in Laboratory Tests. Int. J. Coal Sci. Technol. 2023, 10, 10. [Google Scholar] [CrossRef]
  25. Chen, J.; Shi, K.; Pu, Y.; Apel, D.B.; Zhang, C.; Zuo, Y.; Chen, J.; Xu, L.; Gui, Z.; Song, L. Study on Instability Fracture and Simulation of Surrounding Rock Induced by Fault Activation under Mining Influence. Rock Mech. Bull. 2023, 2, 100037. [Google Scholar] [CrossRef]
  26. Zhang, C.; Bai, Q.; Han, P. A Review of Water Rock Interaction in Underground Coal Mining: Problems and Analysis. Bull. Eng. Geol. Environ. 2023, 82, 157. [Google Scholar] [CrossRef]
  27. Wang, G.; Wu, M.; Wang, R.; Xu, H.; Song, X. Height of the Mining-Induced Fractured Zone above a Coal Face. Eng. Geol. 2017, 216, 140–152. [Google Scholar] [CrossRef]
  28. Lv, X.F.; Zhou, H.Y.; Wang, A.W.; Feng, C.; Xiao, X.C. Characteristics of Stress Transfer and Progressive Fracture in Overlying Strata Due to Mining-Induced Disturbances. Adv. Civ. Eng. 2018, 2018, 8967010. [Google Scholar] [CrossRef]
  29. Deng, G.; Xie, H.; Gao, M.; Li, C.; He, Z. Numerical Simulation on the Evolution of Mining-Induced Fracture Network in a Coal Seam and Its Overburden under the Top Coal Caving Method. Adv. Civ. Eng. 2020, 2020, 8833193. [Google Scholar] [CrossRef]
  30. Lian, X.; Zhang, Y.; Yuan, H.; Wang, C.; Guo, J.; Liu, J. Law of Movement of Discontinuous Deformation of Strata and Ground with a Thick Loess Layer and Thin Bedrock in Long Wall Mining. Appl. Sci. 2020, 10, 2874. [Google Scholar] [CrossRef]
  31. Wang, H.; Jiang, Y.; Zhao, Y. Numerical Investigation of the Dynamic Mechanical State of a Coal Pillar During Longwall Mining Panel Extraction. Rock Mech. Rock Eng. 2013, 46, 1211–1221. [Google Scholar] [CrossRef]
  32. Sun, Y.; Zuo, J.; Karakus, M.; Wang, J. Investigation of Movement and Damage of Integral Overburden during Shallow Coal Seam Mining. Int. J. Rock Mech. Min. Sci. 2019, 117, 63–75. [Google Scholar] [CrossRef]
Figure 1. Location of the Shendong coalfield: (a) 12401 Panel (b) central dynamic fissure, (c) boundary fissure, (d) step fissure, (e) R113 borehole column.
Figure 1. Location of the Shendong coalfield: (a) 12401 Panel (b) central dynamic fissure, (c) boundary fissure, (d) step fissure, (e) R113 borehole column.
Fractalfract 07 00699 g001
Figure 2. Borehole column diagram.
Figure 2. Borehole column diagram.
Fractalfract 07 00699 g002
Figure 3. Contour map of loose layer thickness.
Figure 3. Contour map of loose layer thickness.
Fractalfract 07 00699 g003
Figure 4. Contour map of bedrock thickness.
Figure 4. Contour map of bedrock thickness.
Fractalfract 07 00699 g004
Figure 5. Contour map of thickness ratio of the loose layer to bedrock.
Figure 5. Contour map of thickness ratio of the loose layer to bedrock.
Fractalfract 07 00699 g005
Figure 6. Test model: (a) related equipment, (b) similar material model of design.
Figure 6. Test model: (a) related equipment, (b) similar material model of design.
Fractalfract 07 00699 g006
Figure 7. Different advance distances of rock failure area.
Figure 7. Different advance distances of rock failure area.
Fractalfract 07 00699 g007
Figure 8. Overlying strata model of the 12401 working face.
Figure 8. Overlying strata model of the 12401 working face.
Fractalfract 07 00699 g008
Figure 9. Evolution of horizontal displacement clouds and force chains in overlying strata at different advance distances.
Figure 9. Evolution of horizontal displacement clouds and force chains in overlying strata at different advance distances.
Fractalfract 07 00699 g009
Figure 10. Simulation experiment: (a) advance distance 67.5 cm, fracture network evolution, schematic diagram for calculating Df, (b) advance distance 120 cm, fracture network evolution, schematic diagram for calculating Df.
Figure 10. Simulation experiment: (a) advance distance 67.5 cm, fracture network evolution, schematic diagram for calculating Df, (b) advance distance 120 cm, fracture network evolution, schematic diagram for calculating Df.
Fractalfract 07 00699 g010
Figure 11. The characteristic curve of overlying strata failure.
Figure 11. The characteristic curve of overlying strata failure.
Fractalfract 07 00699 g011
Figure 12. Schematic diagram of overlying strata movement and ground fissure development and deformation.
Figure 12. Schematic diagram of overlying strata movement and ground fissure development and deformation.
Fractalfract 07 00699 g012
Figure 13. Ground fissure width variation curve: (a) f08, f14 central ground fissure, (b) hf4, hf5 central ground fissure, (c) hf2, hf3 central ground fissure, (d) location of ground fissure.
Figure 13. Ground fissure width variation curve: (a) f08, f14 central ground fissure, (b) hf4, hf5 central ground fissure, (c) hf2, hf3 central ground fissure, (d) location of ground fissure.
Fractalfract 07 00699 g013
Table 1. Statistics of borehole lithology around the No. 12401 working face.
Table 1. Statistics of borehole lithology around the No. 12401 working face.
No.Loose Sandy Layer Thickness (m)Bedrock Thickness (m)Thickness Ratio of the Loose Layer to BedrockThickness Ratio of SandstoneThickness Ratio of Mudstone
R11019.60165.9011.81%74.98%14.25%
R11322.10219.1010.09%68.08%10.86%
R10618.20172.2010.56%90.14%0.29%
b2799.00210.004.28%91.78%4.28%
Table 2. Physical and mechanical parameters of rocks in PFC model.
Table 2. Physical and mechanical parameters of rocks in PFC model.
LithologyDensity (KN/m3)Tensile Strength
(MPa)
Cohesion
(MPa)
Friction Angle (°)KratEmod
(GPa)
Pb-kn
(GPa)
Pb-ks
(GPa)
aeolian sand17--181.00.37.83.9
fine sandstone243.221.30371.013.311.610.4
siltstone254.052.55371.012.715.89.0
fine sandstone243.221.30361.013.311.610.4
siltstone254.052.55371.010.713.89.0
fine sandstone243.221.30361.013.311.610.4
medium sandstone242.584.35381.09.111.69.0
coal161.071.21391.05.15.82.6
siltstone254.052.55381.010.713.89.0
Strata macroscopic and mesoscopic physic and mechanical parameters. Krat is the normal-to-shear stiffness ratio, Emod the effective modulus, Pb-kn the normal stiffness, Pb-ks the shear stiffness.
Table 3. Overlying strata failure results of different advance distances.
Table 3. Overlying strata failure results of different advance distances.
Mining Length (cm)Fractal Dimension DCorrelation CoefficientOverlying Strata Failure Height
(cm)
Rock Failure Area
(cm2)
Failure Stage of Overlying Strata
60.01.3250.990114.41485immediate roof caving
67.51.3390.990121.63354first key stratum broken
120.01.3930.997668.86285primary key stratum broken
135.01.3950.997584.87949overlying strata bending
142.51.3630.9971110.010,169integral caving
Table 4. Dynamic fissure width.
Table 4. Dynamic fissure width.
Measurement DateAdvanced Distance
(m)
Central Fissure Width (m)Boundary Fissure Width (m)
f08f14hf4hf5hf2hf3
4 June3453.00.201.003.0
5 June3608.00.501.0017.0
6 June37510.00.601.6023.0
7 June39038.01.202.7038.0
8 June39553.02.303.1041.0
9 June41051.02.703.1042.0
10 June42040.02.602.0039.0
11 June43535.02.402.1038.5
12 June4504.032.02.002.004.036.0
13 June45715.031.51.801.809.035.5
14 June46036.032.51.701.9021.027.0
15 June47042.032.51.601.8532.025.0
16 June48036.036.01.501.8035.024.0
17 June49128.038.01.401.8030.023.2
18 June51026.039.01.401.7029.523.0
19 June51525.036.01.431.8029.023.3
20 June52325.535.01.451.8526.023.4
21 June53625.034.51.391.9026.523.6
22 June54526.033.01.381.8023.023.0
23 June56029.030.51.391.8021.523.0
24 June57232.030.01.4020.523.0
25 June58134.029.01.4020.023.0
26 June59333.025.01.4219.823.0
27 June62031.022.01.4519.5
28 June63329.021.51.4319.6
29 June64526.022.01.4219.5
30 June66925.022.51.4619.8
1 July67622.022.01.4019.3
2 July68120.022.01.3819.0
3 July69619.51.3719.5
4 July71018.01.3819.0
5 July71817.51.40
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ling, C.; Liu, B.; Zhang, C.; Teng, T.; Zhang, K.; Sun, B.; Zhou, J. Fractal Characteristics of Overburden Rock Fractures and Their Impact on Ground Fissures in Longwall Coal Mining. Fractal Fract. 2023, 7, 699. https://doi.org/10.3390/fractalfract7100699

AMA Style

Ling C, Liu B, Zhang C, Teng T, Zhang K, Sun B, Zhou J. Fractal Characteristics of Overburden Rock Fractures and Their Impact on Ground Fissures in Longwall Coal Mining. Fractal and Fractional. 2023; 7(10):699. https://doi.org/10.3390/fractalfract7100699

Chicago/Turabian Style

Ling, Chunwei, Bin Liu, Cun Zhang, Teng Teng, Kangning Zhang, Bo Sun, and Jinlong Zhou. 2023. "Fractal Characteristics of Overburden Rock Fractures and Their Impact on Ground Fissures in Longwall Coal Mining" Fractal and Fractional 7, no. 10: 699. https://doi.org/10.3390/fractalfract7100699

APA Style

Ling, C., Liu, B., Zhang, C., Teng, T., Zhang, K., Sun, B., & Zhou, J. (2023). Fractal Characteristics of Overburden Rock Fractures and Their Impact on Ground Fissures in Longwall Coal Mining. Fractal and Fractional, 7(10), 699. https://doi.org/10.3390/fractalfract7100699

Article Metrics

Back to TopTop