Fractal Study of the Development Law of Mining Cracks

: 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, ﬂuctuations vary from 250 to 287.5 m, forming a wa-vering 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 identiﬁcation, noting limitations in identifying tiny cracks.


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.

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.

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.

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.seven small wooden blocks were used instead of coal seam mining space, width of 25 mm, corresponding to an actual distance of 12.5 m.After modeling similar materials, this study extracted small wooden blo late the mining process.Each small wooden block removed represents a 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.

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: Box counting is commonly used to determine D, originally defined by [30].

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: Box counting is commonly used to determine D, originally defined by [30].
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.

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.
Fractal Fract.2023, 7, x FOR PEER REVIEW 6 of 17 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.

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.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  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.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.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.

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.

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 cm 2 , respectively.The relationship between the fractal dimension and the mining crack area cannot be seen from the area change.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 cm 2 , 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.
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.

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

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 seam accounted for 25.04% of the mining space, and the intermediate process in as a whole, although there were slight undulations.When it advanced to 337.5 m, i conducted to the surface to form a subsidence basin, at which time the mining cr the overlying rock mass accounted for 50.4% of the mining space.With the increas propulsion distance, the percentage of mining fracture area for the overlying roc increases in the mining space; that is, the overall upward propagation of the minin is enhanced as the cracks widen.
Figure 7 shows that the change rate of the mining crack area of the overlyi mass for each propulsion distance to the mining space area has two slower parts are from 175 to 262.5 m and 287.5 to 337.5 m.In these two areas, the change in the p age of the mining fracture area of the overlying rock layer was relatively slow, b first decreased and then increased.This is because the mining space propagates u in a separatory manner, and the overlying rock layer on the coal seam collapses, fi the old mining cracks to form new ones and gradually compacting in the process.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.The average value of fracture height was selected to study the mining proces Figure 8, the correlation coefficient between the fractal dimension and the advancing tance is 0.56, indicating that the correlation between the fractal dimension and the adv ing distance is not high, which is due to the subsequent compaction process durin advancing process.
Figures 8 and 9 show that the propulsion distance is positively correlated wit fracture development height and mining fracture area of the overlying rock layer correlation coefficient reaches 94%.The height of fissure development varies signific with propulsion distance.

Conduction Law of Mining Cracks
To study the conduction law of mining fractures, we must study the tion of mining cracks in rock formations during the final transmission to t ing a surface subsidence basin.During the final propagation of the min surface, the surface subsidence basin advanced to 337.5 m.This is depic and Figure 11, which is divided into ten parts evenly in the horizontal an tions after the cropping operation.

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 Figures 10 and 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: which is significantly smaller than the mining space area and is primarily caused by the 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 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.
The actual surface subsidence area was calculated by Photoshop and was A s2 = 75.33cm 2 .
In Equation (3), A is the subsidence area, S 1 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.
In Equation ( 4), k is the broken expansion coefficient, H 1 is the fissure development height maximum, and H 2 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.

1.
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.

2.
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.

3.
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.
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.5.
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.

Figure 4 .
Figure 4. Fractal dimension and a propulsion distance plot.

Figure 4 .
Figure 4. Fractal dimension and a propulsion distance plot.

Figure 7 .
Figure 7. Percentage of mining fractures in mining space area in overlying rock layer line.

Figure 7 .
Figure 7. Percentage of mining fractures in mining space area in overlying rock layer line.

Figures 8
Figures 8 and 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.

Figure 9 .
Figure 9. Relationship between mining fracture area, fracture development heigh distance in the overlying rock layer.

Figure 9 .
Figure 9. Relationship between mining fracture area, fracture development height, and propulsion distance in the overlying rock layer.

Figure 10 .
Figure 10.Cutting of mining fracture area.

Figure 10 .
Figure 10.Cutting of mining fracture area.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.

Figure 12 .
Figure 12.Broken-line diagram of horizontal slice and vertical slice mining fractures.

Figure 12 .
Figure 12.Broken-line diagram of horizontal slice and vertical slice mining fractures.

Table 2 .
Table of fractal dimension calculation results.

Table 2 .
Table of fractal dimension calculation results.

Table 3 .
Percentage growth of fractal dimension.

Table 3 .
Percentage growth of fractal dimension.

Table 4 .
Mining area and mining crack area.

Table 4 .
Mining area and mining crack area.

Table 5 .
Section mining fracture area comparison

Table 5 .
Section mining fracture area comparison.

Table 6 .
The proportion of mining fracture area S1 and mining fracture area S2 in the overlying rock layer of the coal seam.

Table 7 .
Statistical table of fracture development height.

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.

Table 8
, the total percentage of mining cracks in the rock formation in the mining space is: