Next Article in Journal
The Multiple Frequency Conversion Sinusoidal Chaotic Neural Network and Its Application
Previous Article in Journal
Dynamic Compressive Mechanical Property Characteristics and Fractal Dimension Applications of Coal-Bearing Mudstone at Real-Time Temperatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractal Study of the Development Law of Mining Cracks

1
School of Safety Engineering, North China Institute of Science and Technology, Langfang 065201, China
2
Hebei State Key Laboratory of Mine Disaster Prevention, Langfang 065201, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(9), 696; https://doi.org/10.3390/fractalfract7090696
Submission received: 15 August 2023 / Revised: 6 September 2023 / Accepted: 8 September 2023 / Published: 19 September 2023

Abstract

:
Studying mining fracture development is vital for geotechnical and mining engineering and geological disaster prevention. This research assesses crack effects on rock mass stress equilibrium during coal mining, potentially causing geological disasters such as land subsidence and landslides. Using fractal geometry theory, the present study investigates the development of horizontal and vertical mining cracks, revealing their propagation patterns. The fractal dimension generally increases as the propulsion distance increases; however, fluctuations vary from 250 to 287.5 m, forming a wavering line chart. The proportion of mining fracture area relative to mining space area increases with greater propulsion distance, indicating expanded upward mining space due to separation layers. The horizontal distribution of mining cracks persists, while the vertical distribution decreases, suggesting ground subsidence results from upward transmission. The fastest increase in fractal dimension occurs at 87.5–100 m. At 250 m, it peaks at 1.4136, indicating complex crack structures. During propulsion, the fractal dimension decreases due to upward mining space expansion through overlying rock layer collapse, forming new cracks. The proportion of mining crack area to mining space area increases gradually throughout the mining process. The present study presents a simulation model for crack identification, noting limitations in identifying tiny cracks.

1. Introduction

During coal mining activities, the stress within the rock body is redistributed, causing a disruption in the stress equilibrium of the overlying rock layer above the coal seam. This disturbance results in the subsidence and displacement of the overburden layer, ultimately leading to the formation of a complex network of staggered fractures. The evolution of mining fracture networks directly affects the permeability of coal and rock masses, as well as the transport of gases. When the evolution of mining cracks is large, geological disasters such as land subsidence and mine water inrush can occur. Therefore, studying the development of mining fractures plays a crucial role in geotechnical engineering, mining engineering, and geological disaster prevention [1,2,3,4].
In recent years, scholars have used similar simulation experiments, numerical simulations, physical simulation models, and other methods to study the deformation of overlying strata and the evolution, shape, and distribution of fractures under mining conditions. Wang et al. [5] and Liu [6] used similar simulation experiments to simulate the formation process and distribution of overburden fractures caused by mining, revealing their evolution. Studies have found that fissures develop highly discontinuous jumps. Xu et al. [7] and Huang [8] used physical simulation and field measurements to study the evolution of overlying formations and fractures caused by mining. The results show that rock separation fractures and vertical fractures occur in the overlying strata, and the overburden damage evolution is trapezoidal. Zhou et al. [9] and Wei [10] established a physical simulation model based on geological conditions, studied the distribution and evolution of mining-induced fractures in overburden strata, and extracted three typical crack networks from digital photos using image analysis. Zhou et al. [11] and Wei [12] used theoretical analysis and numerical simulation, combined with actual survey data, to study the formation and development mechanisms of ground cracks caused by coal mining. With the increase in the number of coal seams, the fracture rate and development height of the overburden increase.
Scholars have performed substantial research on cracks and microcracks using fractal theory. Numerous experiments have demonstrated that cracks in nature are fractals. The fractal property is one of the major properties of complex networks. Current research is concentrated on single-layer static networks and mesoscale applications. The fractal dimension can show the fractal properties of the network and the central nodes. Fractal properties can be used to address real problems involving network vulnerability metrics and node importance identification [13]. For example, Nayak et al. [14] proposed a generalized box count (grayscale invariant DBC) based on the fractal dimension (FD) in image analysis. The proposed model exhibits high accuracy in terms of a small fitting error when detecting surface roughness from a given dataset.
Simulations of the specimens show that fractal dimensions affect the crack growth process [15]. Liu et al. [16] studied the evolution of mining fractures and quantitatively described their fractal laws. The results demonstrate that the fractal dimension of the fractures generated by coal mining in the overlying rock undergoes two periodic changes as the coal surface advances. This fractal dimension serves as an indicator of fracture evolution and can be used for stability assessment. Wei [10] found that overburden formation pressure plays an important role in the evolution of mining-induced fractures, and the fractal dimension of fractures increases with the increase in overburden pressure. Gao et al. [17,18] used the box dimension method to quantitatively analyze the changes of coal sample cracks under different load conditions; with the increase in mining width, the fractal dimension of the crack network showed a staggered upward trend, which was divided into the rapid fracture dimension rising, slow fracture size rising, and fracture dimension stability stages. Cai [19] used the improved box counting method to measure the dynamic characterization of the extension of the fracture network in porous media, and the fractal dimensions of the fracture network were related to drying time, average pore size, moisture content, and fracture porosity, respectively. Zheng [20] combined digital image processing technology and a rock failure process analysis system to study the influence of fracture geometric distribution and microscopic heterogeneity on material structure and proposed a damage variable index based on fractal theory. Chen et al. [21,22] introduced fractal dimension to describe fracture surfaces with complex geometries and established a functional relationship between rock fracture bulk density and fractal dimension. Deng [23] proposes a three-dimensional modeling method based on discrete fracture networks and fractal theory. After mining, the lower cracks increase linearly. In higher formations, most of the cracks occur behind mining faces, and lower formations produce more cracks than the upper layers. The results of Feng [24] show that the coal crack propagation and internal crack surface show significant fractal characteristics. In addition, during dynamic loading, the multifractal feature becomes increasingly prominent. Gou [25] studied the evolution and fractal characteristics of buried rock fractures under mining disturbances. The results show that the correlation coefficient of the fractal dimension of rock fractures at different propulsion distances is greater than 0.90, and the rock fracture network under karst geomorphology has high self-similarity.
Under certain static rock stress conditions, the existence of microcomponents in coal will affect fracture abundance in similar coal reservoirs and then affect the initial permeability of the reservoir. The change in coal composition and rupture strength influence the rheological characteristics of coal system reservoirs [26]. Bandyopadhyay et al. [27] proposed the deployment of a calibrated discrete fracture network (DFN) modeling technology based on mesoscale coal images. Furthermore, a sensitivity analysis of the key fracture properties was performed to quantify their impact on coal permeability. Widely accepted predicted crack features include crack width, spacing, and number of cracks. Theodoros et al. [28] studied the mechanical factors of tensile deformation from a cracking behavior perspective, established empirical equations to predict the number of cracks, crack width and crack spacing, and predicted the average width of cracks. In addition, Jiang et al. [29] conducted hydraulic fracturing experiments using coal blocks and showed that fracture behavior at different interfaces mainly depends on the vertical stress and the interfacial friction coefficient. In the present study, a hydraulic fracture prediction model that predicts the propagation behavior of hydraulic cracks during hydraulic fracturing with high accuracy is developed.
In summary, fractal theory can be used to study the evolution and distribution characteristics of fractures, but the influence of fracture height, area, fractal dimension, and other parameters on mining activities and the propagation of fractures requires further research. In this study, the changes in shallow-buried coal seams obtained by simulating coal seam mining through similar simulation experiments were studied, and Image-Pro Plus image processing software 6.0 and Photoshop 2019 were used to identify and extract mining cracks. This paper analyzed the relationship between propulsion distance, mining crack area, coal mining space, and fractal dimension. The fracture changes of the overlying rock layer created during coal seam mining were studied based on fractal theory. The degree of rock mass failure was characterized by fractal dimensions. The development of mining fractures and their impact on the surface were predicted.

2. Materials and Methods

2.1. Experimental Method

Similarity simulation experiments are a common approach for understanding and analyzing traffic and roadway problems. It overcomes numerous inconveniences in field investigations and experiments and makes up for deficiencies in theoretical analysis.
Based on the principle of self-similarity theory, a similar material model is built on the two-dimensional similarity simulation test bench based on the ramming filling method, and the coal seam in the model is used to simulate mining. The experiment process is recorded, and the displacement, deformation, collapse, and destruction of overlying strata of coal seams after mining are captured.

2.2. Model Mining and Observation

After a similar material model is built, the mining process is simulated by tapping small wood blocks. After each tap, the overburden is waited for to halt collapsing in the mining space before advancing to the subsequent block. A coherent light source is set up for image acquisition. For instance, halogen lamps are placed around a similar analog material model, and shooting is performed after stopping the collapse after each tap of the small wood block to obtain image data. Therefore, the deformation of the overlying strata of the coal seam can be observed. Relevant image recognition software and programs are used to extract data.

3. Simulation Experiment of Mining Crack

The outline of the roadway system in the Shangwan Coal Mine is approximately rectangular, and the floor of the coal seam working face is black-gray mudstone. The overlying strata of the coal seam are mainly yellow medium-grained eolian sand, fine-grained eolian sand, loose soil, and unconsolidated soil. The simulation experiment was based on the Shangwan wellfield, and a similar engineering background material model was built based on the regional geological overview of the mine and the mechanical parameters of the rock mass. The similar material model adopts a two-dimensional similarity simulation test bench, and the ratio of the model to the actual mine is 1:250. The dimensions of the model are: length × width × height = 3200 mm × 250 mm × 870 mm. The stratification of the model is shown in Table 1. The mining model diagram is shown in Figure 1. Twenty-seven small wooden blocks were used instead of coal seam mining space, each with a width of 25 mm, corresponding to an actual distance of 12.5 m.
After modeling similar materials, this study extracted small wooden blocks to simulate the mining process. Each small wooden block removed represents a 12.5 m propulsion distance. After each extraction, wait for the overlying rock layer to stop collapsing before continuing to extract the next wooden block until all the small wooden blocks are removed, which is equivalent to 337.5 m of mining coal seam advance. When the propulsion distance was equal to 87.5 m (the 7th small wooden block), the overlying rock layer of the coal seam began to collapse, and mining cracks appeared in the observation area. Crack image data were obtained by shooting after each twitch of the small wooden block after collapse. Thus, deformation of the overlying rock layer on the coal seam was observed. The process from the beginning of the collapse of the rock formation (the 7th block) to the last small block is shown in Figure 2.

4. Image Processing and Data Analysis Simulation

4.1. Fractal Theory

In fractal geometry, the box dimension is a method of measuring the fractal dimension in distance spaces, such as Euclidean space R. In practical applications, the box dimension can be understood as constructing squares (bodies, called boxes) with sides of r to cover the set (which can be a line, surface, or volume), calculating the number of “boxes” with different r values and F intersecting N(r), and box dimension D can be obtained from the slope value of the straight line in the double logarithmic coordinate system of N(r) and r. D is the box dimension of traditional scale r and characterizes the efficiency of a small set of the same shape covering a universal set. Box dimensions can be defined as:
dim box ( s ) = lim r 0 log N ( r ) log ( 1 / r )
Box counting is commonly used to determine D, originally defined by [30].
N ( r ) ~ r D f
Equation (2) is called the fractal scale law, where N(r) represents the number of boxes occupied by pixels, and r is the size of the box or scale used to measure stellar islands/clouds/galaxies.

4.2. Fractal Dimension

We used Photoshop software to color correct images, denoise, reduce background interference, and reduce calculation errors. Image-Pro Plus was used to binarize the cropped pictures, and then the fractal dimension of the mining crack was calculated based on the box dimension method using MATLAB 2021a. The image binarization effect and the MATLAB calculation result for an advancement of 175 m (Block 14) are shown in Figure 3. The full calculation results are shown in Table 2.
From Table 3 and the percentage increase in fractal dimension with propulsion distance, the fractal dimension is shown to increase with the increase in propulsion distance. This suggests that the complexity of mining-induced fractures grows as the propulsion distance increases. At a propulsion distance of 87.5 m, the overlying rock layer of the coal seam began to collapse, and mining cracks appeared in the observation area. The mining crack structure was the simplest at this time. At a propulsion distance of 250 m, the fractal dimension reached its maximum, indicating that the mining crack structure was the most complex at this time.
When advancing from 87.5 to 100 m, the fractal dimension increased from 1.2559 to 1.2950, a total of 3.11%; the percentage increase in this stage was the largest in the selected observation images, indicating that the damage and failure transmission of mining cracks was the fastest. When advancing from 312.5 to 325 m, the fractal dimension increased from 1.3730 to 1.3746, a total of only 0.12%; the percentage increase in fractal dimension at this stage was the smallest in the selected observation images (Figure 4), indicating that the damage and failure conduction of mining cracks was the slowest.
Although the fractal dimension exhibits an overall upward trend with the increase in the advancing distance, the change in fractal dimension is significant during the advance from 250 to 287.5 m, and the line chart of fractal dimension changes with the advancing distance is shown as an approximate wavy line.

4.3. Mining Crack Area

The image parameters were extracted using Image-Pro Plus, and mining space area data, mining fracture area data, mining crack area data of overlying coal seams, and mining fracture area data in mining space are listed in Table 4. Details of these areas are shown in Figure 5.
In the propulsion phase from 250 to 287.5 m, the fractal dimension changes with the propulsion distance, and the line chart appears as an approximate wavy line. The advancement process at this stage was analyzed; from the data of the mining crack area (Table 4), the mining crack areas at 250, 262.5, 275, and 287.5 m were found to be 260.07662, 285.79540, 339.95989, and 296.40685 cm2, respectively. The relationship between the fractal dimension and the mining crack area cannot be seen from the area change.
Using a section selection tool, the binarized images at 250, 262.5, and 275 m were divided into ten aliquots of landscape orientation. The slice results are shown in Figure 6. The area data identified and extracted are shown in Table 5. According to Table 5, When advancing from 250 to 275 m, new large-scale mining cracks appeared in slices 1–3, collapsed in slices 4, did not change significantly in slices 5–9, and increased greatly in area in slice 10 owing to the increase in propulsion distance. Combined with the fractal dimension, the area change in slices 1–4 was studied. The following conclusions were drawn: when the fractal dimension of the mining crack area increases, a large area of the formed mining crack collapses and compacts, and a large area expansion occurs vertically and horizontally to form a new mining crack, which is the main factor for the decrease in the fractal dimension of the area increase. After the mining crack is compacted, as advancement continues, the mining crack area increases. The compacted part does not change significantly, and a new collapse and mining crack form in the mining space; the fractal dimension also increases.

4.4. Connection between Mining Cracks and Mining Spaces

Table 6 compares mining fracture area S1 of the overlying rock layer of the coal seam for each propulsion distance and mining fracture area S2 in the mining space with the total area. A line chart of the percentage of mining fracture area in the overlying rock layer as a function of the mining space area is plotted in Figure 7.
When advancing 87.5 m, the mining crack area of the overlying rock mass of the coal seam accounted for 25.04% of the mining space, and the intermediate process increased as a whole, although there were slight undulations. When it advanced to 337.5 m, it finally conducted to the surface to form a subsidence basin, at which time the mining cracks of the overlying rock mass accounted for 50.4% of the mining space. With the increase in the propulsion distance, the percentage of mining fracture area for the overlying rock mass increases in the mining space; that is, the overall upward propagation of the mining space is enhanced as the cracks widen.
Figure 7 shows that the change rate of the mining crack area of the overlying rock mass for each propulsion distance to the mining space area has two slower parts, which are from 175 to 262.5 m and 287.5 to 337.5 m. In these two areas, the change in the percentage of the mining fracture area of the overlying rock layer was relatively slow, but both first decreased and then increased. This is because the mining space propagates upwards in a separatory manner, and the overlying rock layer on the coal seam collapses, filling in the old mining cracks to form new ones and gradually compacting in the process.
From Table 6, the ratio of mining fracture area to mining space area is shown to be relatively stable, and the average value is 22.59%. Pores left by the collapse of the overlying rock layer create the new “cracks”. However, these are not true mining cracks; therefore, the transmission of the mining cracks in the mining space is uniform.
In the mining process, the cracks were distributed in a trapezoidal shape, and a certain height formed during the upward propagation. The fracture height of the image at different propulsion distances was measured, and then the statistics were calculated. The results are shown in Table 7.
The average value of fracture height was selected to study the mining process. In Figure 8, the correlation coefficient between the fractal dimension and the advancing distance is 0.56, indicating that the correlation between the fractal dimension and the advancing distance is not high, which is due to the subsequent compaction process during the advancing process.
Figure 8 and Figure 9 show that the propulsion distance is positively correlated with the fracture development height and mining fracture area of the overlying rock layer. The correlation coefficient reaches 94%. The height of fissure development varies significantly with propulsion distance.

4.5. Conduction Law of Mining Cracks

To study the conduction law of mining fractures, we must study the spatial distribution of mining cracks in rock formations during the final transmission to the surface, forming a surface subsidence basin. During the final propagation of the mining crack to the surface, the surface subsidence basin advanced to 337.5 m. This is depicted in Figure 10 and Figure 11, which is divided into ten parts evenly in the horizontal and vertical directions after the cropping operation.
Number of slices: Horizontal slices are numbered from top to bottom; vertical slices are numbered from left to right, 1–10. We used Photoshop to count the number of pixels and calculate the area of the measured area based on the conversion relationship between pixels and area. The data are listed in Table 8. The percentage of mining cracks in the horizontal and vertical slices is plotted, and the results are shown in Figure 12.
The horizontal slices reflect the propagation of the mining crack in the vertical direction, and the vertical slice reflects that in the horizontal direction.
Mining cracks propagate from the bottom up; thus, when studying the propagation law of mining cracks, it is logical to look at the transverse sections starting with slice 10. Slices 9 and 10 were divided into the first part, and the percentage of the mining crack area to mining space area exhibited a sharp and significant decrease, indicating that mining cracks were denser in these two slices. Slices 1–8 were divided into the second part, and the percentage value of the mining crack area to the mining space area was evenly distributed; the average value was 1.22%. Figure 12 shows that the mining crack gradually weakens when propagating vertically because it has been repeatedly compacted and closed, forming new mining cracks until ground subsidence occurs. The new mining cracks are formed in the upper layer of the compacted rock layer.
Similarly, in the vertical slices, the slice is divided into two parts. Slices 1–2 and 9–10 are the first part, and slices 3–8 are the second part. The average value of the first part is 8.93%, and that of the second part is 3.85%. Within these two parts, the distribution of mining cracks is relatively uniform. In the second part, the percentage of mining crack area to mining space area is considerably smaller than that of the first part, because the second part is located in the middle of the mining space, and its comprehensive rock and soil stress is less than that at both ends. Thus, it is more prone to collapse, compaction, and closure.
Based on the data in Table 8, the total percentage of mining cracks in the rock formation in the mining space is:
S 1 S = 271.35 461.22 × 100 % = 58.83 %
which is significantly smaller than the mining space area and is primarily caused by the propagation of the mining space from bottom to top through separation. Theoretically, the subsidence area of the earth can be calculated.
A s 1 = ( 1 S 1 S ) × S = 189.8   cm 2
The actual surface subsidence area was calculated by Photoshop and was A s 2 = 75.33   cm 2 .
In Equation (3), A is the subsidence area, S1 is the mining fracture area, and S is the mining space area.
The swelling of rock is mainly characterized by soft rock, and the volume of roof rock in the goaf area increases. The broken expansion coefficient is usually described by the fracture height.
k = H 1 H 1 H 2 = 193.623 193.623 8.73 = 1.05
In Equation (4), k is the broken expansion coefficient, H1 is the fissure development height maximum, and H2 is the thickness of the coal seam.
These data show that the broken expansion coefficient reached 1.05 during mining, indicating that the fragmentation expansion was created and the area increased.

5. Conclusions

  • The fractal dimension increases the fastest when advancing from 87.5 to 100 m; therefore, the mining crack conducts the fastest at this stage. When advancing to 250 m, the fractal dimension reached a maximum value of 1.4136, indicating that the mining crack structure was the most complex. During propulsion, the fractal dimension has a negative growth because the mining space propagates upward in a separational manner and the overlying rock layer on the coal seam collapses, filling in the old mining cracks to form new mining cracks and gradually compacting.
  • During the entire mining process, the percentage of the mining crack area of the overlying rock layer in relation to the mining space area gradually increases; that is, the mining crack gradually increases when the mining crack propagates upward. The mining crack is formed by the collapse of the overlying rock layer in the mining space. Thus, the ratio of the mining crack area to the mining space area is relatively stable. However, the mining crack area in the middle of the mining space is slightly smaller than the two ends.
  • The propulsion distance is positively correlated with the fracture development height and mining fracture area of the overlying rock layer. The correlation coefficient reaches 94%. The height of fissure development varies significantly with propulsion distance.
  • The percentage of mining fracture area in the mining space area during the final transmission of the mining crack to the surface was analyzed to study the propagation of the mining fracture and its final impact on the surface. The distribution of mining cracks is relatively uniform in the horizontal direction, and the distribution in the vertical direction is gradually weakened; therefore, the propagation of mining cracks is gradually reduced. Slices 9 and 10 are concentrated areas of mining cracks, and the structure of mining cracks is also more complex. The area of mining cracks is significantly smaller than the area of mining space.
  • In conclusion, future research in this field should aim to address the limitations of current methods, enhance our understanding of mining crack behavior, and develop practical solutions to minimize the risks associated with mining cracks in both underground and surface mining operations.

Author Contributions

Methodology, S.Y. (Sihai Yi) and J.W.; formal analysis, Y.L.; investigation, Y.W.; resources, D.H.; data curation, S.Y. (Shangxian Yin); writing—original draft, D.X.; writing—review and editing, J.W. and D.X.; funding acquisition S.Y. (Sihai Yi) and S.Y. (Shangxian Yin). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (52274165), National Natural Science Foundation of China (51974126), National Natural Science Foundation of China (51204069) and Funding Projects for Basic Scientific Research Business Fees in Central Universities (3142015093).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Z.G.; Zhou, H.W.; Xie, H.P. Research on fractal characterization of mined crack network evolution in overburden rock stratum under deep mining. Rock Soil Mech. 2009, 30, 2403–2408. [Google Scholar]
  2. Cheng, Y.P.; Yu, Q.X.; Yuang, L. Experimental research of safe and high-efficient exploitation of coal and pressure relief gas in long distance. J. China Univ. Min. Technol. 2004, 33, 132–136. [Google Scholar]
  3. Liu, Z.G.; Liang, Y. Study on coal seam roof gas drainage from the strike of annular fracture areas by the long drill method. Eng. Sci. 2004, 6, 32–38. [Google Scholar]
  4. Tu, M. Study on the growth height of separation fracture of mining rock in Panxie area. Mei T’an Hsueh Pao 2004, 29, 641–645. [Google Scholar]
  5. Wang, C.; Zhang, N.C.; Han, Y.F.; Xiong, Z.Q.; Qian, D.Y. 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]
  6. Liu, X.Y. Experimental Study of Dynamic Development and Fractal law of Mining Rock Fracture. Appl. Mech. Mater. 2011, 71, 3428–3432. [Google Scholar] [CrossRef]
  7. Xu, X.L.; Zhang, N.; Tian, S.C. Mining-induced movement properties and fissure time-space evolution law in overlying strata. Int. J. Min. Sci. Technol. 2012, 22, 817–820. [Google Scholar] [CrossRef]
  8. Huang, Q.X.; He, Y.P.; Cao, J. Experimental investigation on crack development characteristics in shallow coal seam mining in China. Energies 2019, 12, 1302–1318. [Google Scholar] [CrossRef]
  9. Zhou, H.W.; Liu, J.F.; Xue, D.J.; Yi, H.Y.; Xue, J.H. Numerical simulation of gas flow process in mining-induced crack network. Int. J. Min. Sci. Technol. 2012, 22, 793–799. [Google Scholar] [CrossRef]
  10. Wei, X.J.; Gao, M.Z.; Lv, Y.C.; Shi, X.C.; Gao, H.L.; Zhou, H.W. Evolution of a mining induced fracture network in the overburden strata of an inclined coal seam. Int. J. Min. Sci. Technol. 2012, 22, 779–783. [Google Scholar] [CrossRef]
  11. Zhou, D.W.; Wu, K.; Bai, Z.H.; Hu, Z.Q.; Li, L.; Xu, Y.K.; Diao, X.P. Formation and development mechanism of ground crack caused by coal mining: Effects of overlying key strata. Bull. Eng. Geol. Environ. 2019, 78, 1025–1044. [Google Scholar]
  12. Yang, W.F.; Ji, Y.B.; Shen, D.Y.; Xia, X.H. Discrete element numerical simulation of crack evolution in multi-coal seam mining. In Proceedings of the International Conference on Computer Application and System Modeling, Taiyuan, China, 22–24 October 2010; pp. 383–385. [Google Scholar]
  13. Wen, T.; Cheong, K.H. The fractal dimension of complex networks: A review. Inf. Fusion 2021, 73, 87–102. [Google Scholar] [CrossRef]
  14. Nayak, S.R.; Mishra, J. Fractal dimension-based generalized box-counting technique with application to grayscale images. Fractals 2021, 29, 17. [Google Scholar] [CrossRef]
  15. Chang, X.; Ren, M.F.; Guo, X. A cohesive zone model for self-similar fractal crack propagation. J. Mech. Sci. Technol. 2017, 31, 4763–4769. [Google Scholar]
  16. Liu, C.; Xue, J.H.; Yu, G.F.; Cheng, X.Y. Fractal characterization for the mining crack evolution process of overlying strata based on microseismic monitoring technology. Int. J. Min. Sci. Technol. 2016, 26, 295–299. [Google Scholar] [CrossRef]
  17. Gao, M.Z.; Xie, J.; Guo, J.; Lu, Y.Q.; He, Z.Q.; Li, C. Fractal evolution and connectivity characteristics of mining-induced crack networks in coal masses at different depths. Géoméch. Geophys. Geo-Energy Geo-Resour. 2021, 7, 1–15. [Google Scholar] [CrossRef]
  18. Gao, M.Z.; Zhang, J.G.; Li, S.W.; Wang, M.; Wang, Y.W.; Cui, P.F. Calculating changes in fractal dimension of surface cracks to quantify how the dynamic loading rate affects rock failure in deep mining. J. Cent. South Univ. 2020, 27, 3013–3024. [Google Scholar] [CrossRef]
  19. Cai, J.C.; Wei, W.; Hu, X.G.; Liu, R.C.; Wang, J.J. Fractal characterization of dynamic fracture network extension in porous media. Fractals 2017, 25, 1750023. [Google Scholar] [CrossRef]
  20. Zheng, L.L.; Liu, H.; Zuo, Y.J.; Zhang, Q.P.; Lin, W.; Qiu, Q.; Liu, X.R.; Liu, Z.Q. Fractal study on the failure evolution of concrete material with single flaw based on DIP technique. Adv. Mater. Sci. Eng. 2022, 2022, 1–15. [Google Scholar] [CrossRef]
  21. Chen, Y.; Jin, Y.; Chen, M.; Yi, Z.C.; Zheng, X.J. Quantitative evaluation of rock brittleness based on the energy dissipation principle, an application to type II mode crack. J. Nat. Gas Sci. Eng. 2017, 45, 527–536. [Google Scholar] [CrossRef]
  22. Ma, D.; Duan, H.G.; Liu, J.F.; Li, X.B.; Zhou, Z.L. The role of gangue on the mitigation of mining-induced hazards and environmental pollution: An experimental investigation. Sci. Total Environ. 2019, 664, 436–448. [Google Scholar] [CrossRef] [PubMed]
  23. Deng, G.D.; Xie, H.P.; Gao, M.Z.; Li, C.; He, Z.Q. 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, 6077187. [Google Scholar] [CrossRef]
  24. Feng, J.J.; Wang, E.Y.; Huang, Q.S.; Ding, H.C.; Dang, L. Time-varying multifractal analysis of crack propagation and internal fracture process of coal under dynamic loading. Fractals 2021, 29, 2150089. [Google Scholar] [CrossRef]
  25. Gou, R.T.; Jiang, C.Y.; Liu, Y.; Wang, C.; Li, Y.L. Study on fractal characteristics of evolution of mining-induced fissures in Karst landform. Energies 2022, 15, 5372. [Google Scholar] [CrossRef]
  26. Mukherjee, S.; Rajabi, M.; Esterle, J. Relationship between coal composition, fracture abundance and initial reservoir permeability: A case study in the Walloon Coal Measures, Surat Basin, Australia. Int. J. Coal Geol. 2021, 240, 103726. [Google Scholar] [CrossRef]
  27. Bandyopadhyay, K.; Mallik, J. An optimum solution for coal permeability estimation from mesoscopic scale calibrated stochastic and deterministic discrete fracture network models. Fuel 2023, 331, 125626. [Google Scholar] [CrossRef]
  28. Theodoros, A.C.; Vassilis, P.P. Experimental investigation on cracking behavior of reinforced concrete tension ties. Case Stud. Constr. Mater. 2022, 16, e00810. [Google Scholar]
  29. Jiang, Y.; Lian, H.; Nguyen, V.P.; Liang, W. Propagation Behavior of Hydraulic Fracture across the Coal-rock Interface under Different Interfacial Friction Coefficients and A New Prediction Model. J. Nat. Gas Sci. Eng. 2019, 68, 102894. [Google Scholar] [CrossRef]
  30. Yu, B.M. Analysis of flow in fractal porous media. Appl. Mech. Rev. 2008, 61, 050801. [Google Scholar] [CrossRef]
Figure 1. Mining model diagram.
Figure 1. Mining model diagram.
Fractalfract 07 00696 g001
Figure 2. (a) Propulsion 87.5 m original image (Block 7); (b) Propulsion 100 m original image (Block 8); (c) Propulsion 300 m original image (Block 24); (d) Propulsion 312.5 m original image (Block 25); (e) Propulsion 325 m original image (Block 26); (f)Propulsion 337.5 m original image (Block 27).
Figure 2. (a) Propulsion 87.5 m original image (Block 7); (b) Propulsion 100 m original image (Block 8); (c) Propulsion 300 m original image (Block 24); (d) Propulsion 312.5 m original image (Block 25); (e) Propulsion 325 m original image (Block 26); (f)Propulsion 337.5 m original image (Block 27).
Fractalfract 07 00696 g002
Figure 3. (a) Advance (Block 14) 175 m binarized image; (b) fractal dimension calculation at 175 m in advance (Block 14).
Figure 3. (a) Advance (Block 14) 175 m binarized image; (b) fractal dimension calculation at 175 m in advance (Block 14).
Fractalfract 07 00696 g003
Figure 4. Fractal dimension and a propulsion distance plot.
Figure 4. Fractal dimension and a propulsion distance plot.
Fractalfract 07 00696 g004
Figure 5. Area schematic diagram.
Figure 5. Area schematic diagram.
Fractalfract 07 00696 g005
Figure 6. (a) Horizontal section when advancing 250 m (Block 20); (b) Horizontal section when advancing 262.5 m (Block 21); (c). Horizontal section when advancing 275 m (Block 22).
Figure 6. (a) Horizontal section when advancing 250 m (Block 20); (b) Horizontal section when advancing 262.5 m (Block 21); (c). Horizontal section when advancing 275 m (Block 22).
Fractalfract 07 00696 g006
Figure 7. Percentage of mining fractures in mining space area in overlying rock layer line.
Figure 7. Percentage of mining fractures in mining space area in overlying rock layer line.
Fractalfract 07 00696 g007
Figure 8. Fractal dimension, fracture development height average, and propulsion distance.
Figure 8. Fractal dimension, fracture development height average, and propulsion distance.
Fractalfract 07 00696 g008
Figure 9. Relationship between mining fracture area, fracture development height, and propulsion distance in the overlying rock layer.
Figure 9. Relationship between mining fracture area, fracture development height, and propulsion distance in the overlying rock layer.
Fractalfract 07 00696 g009
Figure 10. Cutting of mining fracture area.
Figure 10. Cutting of mining fracture area.
Fractalfract 07 00696 g010
Figure 11. (a) Horizontal slice; (b) Vertical slice.
Figure 11. (a) Horizontal slice; (b) Vertical slice.
Fractalfract 07 00696 g011
Figure 12. Broken-line diagram of horizontal slice and vertical slice mining fractures.
Figure 12. Broken-line diagram of horizontal slice and vertical slice mining fractures.
Fractalfract 07 00696 g012
Table 1. Model layering scheme.
Table 1. Model layering scheme.
Rock LayerOriginal Layer Thickness
(m)
Model Layer Thickness (cm)The Layer NumberLayer Thickness (cm)Bulk
Density
Ratio
Number
Aeolian sand6.502.3312.601.79982
Fine grained sandstone26.9510.48410.781.79655
Sandy mudstone20.217.8048.081.73673
Medium grained sandstone8.893.4913.561.77646
Sandy mudstone36.7014.21814.681.73673
Coarse grained sandstone39.5715.26715.831.77537
Sandy mudstone21.908.1568.761.73673
Fine grained sandstone33.3712.931013.351.79655
Coal seam9.263.4913.700.93682
Fine grained sandstone23.008.8539.201.79655
Table 2. Table of fractal dimension calculation results.
Table 2. Table of fractal dimension calculation results.
The Number of Small Wooden Blocks That Are PropelledPropulsion Distance (m)Fractal DimensionError of the Fractal Dimensions
787.51.25590.0000
8100.01.2950−0.0015
9112.51.33450.0049
13162.51.36440.0023
14175.01.3507−0.0051
17212.51.37170.0001
18225.01.3885−0.0069
20250.01.41360.0242
21262.51.3780−0.0461
22275.01.40700.0454
23287.51.3737−0.0250
24300.01.36970.0026
25312.51.37300.0018
26325.01.3746−0.0067
27337.51.39050.0000
Table 3. Percentage growth of fractal dimension.
Table 3. Percentage growth of fractal dimension.
Propulsion Process
(m)
Fractal Dimension ChangesPercentage Growth in Fractal
Dimension
87.5–100.01.2559–1.29503.11%
100.0–112.51.2950–1.33453.05%
112.5–162.51.3345–1.36442.24%
162.5–175.01.3644–1.3507−1.00%
175.0–212.51.3507–1.37171.55%
212.5–225.01.3717–1.38851.22%
225.0–250.01.3885–1.41361.81%
250.0–262.51.4136–1.3780−2.52%
262.5–275.01.3780–1.40702.10%
275.0–287.51.4070–1.3737−2.37%
287.5–300.01.3737–1.3697−0.29%
300.0–312.51.3697–1.37300.24%
312.5–325.01.3730–1.37460.12%
325.0–337.51.3746–1.39051.16%
Table 4. Mining area and mining crack area.
Table 4. Mining area and mining crack area.
Propulsion Distance
(m)
Fractal
Dimension
Mining Crack Area Data
(cm2)
Mining Crack Area of Overlying Coal Seam S1
(cm2)
Mining Fracture Area in Mining Space S2
(cm2)
Mining Area S
(cm2)
87.51.255969.1987331.8859737.31276127.33084
100.01.295078.0410530.9465347.09452146.20998
112.51.334595.4190261.8478733.57115162.68079
162.51.3644163.81067109.5557854.25489236.17639
175.01.3507185.72440129.7907955.93361254.90379
225.01.3885216.32550168.7241447.60136324.06072
250.01.4136260.07226190.3796169.69265357.71567
262.51.3780285.79540188.4200397.37537376.46887
275.01.4070339.95989225.80058114.15931393.85330
287.51.3737296.40685201.8433394.56352412.35750
300.01.3697291.98729199.5189592.46834437.66464
312.51.3730279.47761204.0353475.44226454.78760
325.01.3746308.02875213.3005994.72816470.89362
337.51.3905349.14444243.84036105.30408483.76484
Table 5. Section mining fracture area comparison.
Table 5. Section mining fracture area comparison.
Serial NumberPropulsion 250 mPropulsion 262.5 mPropulsion 275 m
PixelAreaPixelAreaPixelArea
1000060.01844
20015,28046.9725214,15143.50184
30025507.8390030819.47136
414,01443.0806912543.8549413794.23921
512543.8549411813.6305316184.97392
616445.0538517305.3182220066.16668
729168.9641326888.26323330210.15074
831709.7449523787.31025362211.13446
916,97352.1770115,96249.0690718,03255.43250
1042,715131.3109549,115150.9853162,966193.56492
aggregate82,686254.1865292,138283.24309110,163338.65407
Table 6. The proportion of mining fracture area S1 and mining fracture area S2 in the overlying rock layer of the coal seam.
Table 6. The proportion of mining fracture area S1 and mining fracture area S2 in the overlying rock layer of the coal seam.
Propel Distance (m)S1/Mining Area SS2/Mining Area S
87.525.04%29.30%
100.021.17%32.21%
112.538.02%20.64%
162.546.39%22.97%
175.050.92%21.94%
212.549.54%20.29%
225.052.07%14.69%
250.053.22%19.48%
262.550.05%25.87%
275.057.33%28.99%
287.548.95%22.93%
300.045.59%21.13%
312.544.86%16.59%
325.045.30%20.12%
337.550.40%21.77%
Table 7. Statistical table of fracture development height.
Table 7. Statistical table of fracture development height.
Propel Distance (m)Mean Value of Fracture Development Height (m)Minimum Fracture Development Height (m)Fissure Development Height Maximum (m)
87.511.522188.9495712.8692
100.011.574148.9433512.8609
112.521.4055219.152622.5522
162.530.2200227.482733.1722
175.060.0007558.00962.47
212.559.5528957.724361.8436
225.0120.4479118.586121.996
250.0120.001119.033121.627
262.5160.48158.378162.33
275.0167.6193166.038169.451
287.5190.8538188.655193.279
300.0190.3442187.77193.721
312.5190.5738187.824193.892
325.0189.9195187.499193.587
337.5190.2236187.56193.623
Table 8. Horizontal slice and vertical slice mining fracture area and proportion of mining space area. Area and proportion of mining space area.
Table 8. Horizontal slice and vertical slice mining fracture area and proportion of mining space area. Area and proportion of mining space area.
Serial NumberThe Area of Cracks Mined in Each Horizontal Section S1The Area of Cracks in Each Vertical Section S2Mining Space Area SS1/SS2/S
11.0267646.07181461.222230.22%9.99%
210.0062538.432622.17%8.33%
33.2708623.510850.71%5.10%
42.2810019.797320.49%4.29%
54.9923716.818501.08%3.65%
66.2189412.213471.35%2.65%
79.0471313.775121.96%2.99%
88.1771520.479771.77%4.44%
973.4622035.1586915.93%7.62%
10152.8697445.0942433.14%9.78%
aggregate271.35240271.3524058.83%58.83%
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

Wu, J.; Xie, D.; Yi, S.; Yin, S.; Hu, D.; Li, Y.; Wang, Y. Fractal Study of the Development Law of Mining Cracks. Fractal Fract. 2023, 7, 696. https://doi.org/10.3390/fractalfract7090696

AMA Style

Wu J, Xie D, Yi S, Yin S, Hu D, Li Y, Wang Y. Fractal Study of the Development Law of Mining Cracks. Fractal and Fractional. 2023; 7(9):696. https://doi.org/10.3390/fractalfract7090696

Chicago/Turabian Style

Wu, Jinsui, Dongyu Xie, Sihai Yi, Shangxian Yin, Dezhi Hu, Yuanyuan Li, and Yun Wang. 2023. "Fractal Study of the Development Law of Mining Cracks" Fractal and Fractional 7, no. 9: 696. https://doi.org/10.3390/fractalfract7090696

APA Style

Wu, J., Xie, D., Yi, S., Yin, S., Hu, D., Li, Y., & Wang, Y. (2023). Fractal Study of the Development Law of Mining Cracks. Fractal and Fractional, 7(9), 696. https://doi.org/10.3390/fractalfract7090696

Article Metrics

Back to TopTop