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

: Ground ﬁssures 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 ﬁssures in large mining height working faces and shallow buried coal seams were explored by considering ﬁeld 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 ﬁssure 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 ﬁssures usually have 3–5 ﬁssure counts with 2–5 m intervals between each group of ﬁssures; central ﬁssures develop twice from open to closed widths; the time period for two central ﬁssures to reach their with maximum widths is 11–20 days; and border ﬁssures expand rapidly with maximum values in a time period of 5–6 days.


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.

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.2MPa and consolidation coefficient of 1.32.The immediate roof is gray fine sandstone with a compressive strength of 13.3-15.2MPa and consolidation coefficient of 1.35.The main roof is grayish-white siltstone with a compressive strength of 14.5-36.6MPa and consolidation coefficient of 2.32.The main roof is silty sandstone with a compressive strength of 14.5-36.6MPa.As shown in Figure 1, there are three types of ground fissures in the Shendong coal field.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).

Overlying Strata Characteristics of the Working Face
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.

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 Figures 3-5, the thicknesses of the bedrock and loose layers were counted by drilling around the working face of the 12,401 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.As shown in Figures 3-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.As shown in Figures 3-5, the thicknesses of the bedrock and loose layers were counted by drilling around the working face of the 12,401 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.

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.
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.For geometric similarity, where L p is the size of the on-site prototype, L m is the size of similar models, and C L is the size similarity ratio.
For dynamic similarity, where γ p is the bulk density of the on-site prototype, γ m is the bulk density of similar models, and Cγ is the similarity ratio of bulk density.For time similarity, where T p is the on-site prototype time, T m is the similarity model time, and C T 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.

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.
advance distance 120 cm advance distance 135 cm advance distance 142.5 cm

Simulation of the Caving Process of the 12,401 Working Face Based on PFC
In terms of numerical simulation, we simplified the complex geological condi without too much addition to the temperature field and seepage flow field.Based o geometry of the Shangwan No. 12401 working face, a two-dimensional model base PFC2D 5.0 software was established in order to further validate the simulated fra evolution characteristics and laws of similar simulations, revealing the evolution ch teristics of the overlying rock fractures at different mining stages.As shown in Figu the model is 300 m long and 200 m high.The upper boundary is free, the left and boundaries fix the x-displacement, and the bottom boundary fixes the y-displacem The particles contact each other through parallel bonding, the particle diameter ra from 0.3 to 0.5 m, and a total of 102,814 balls are generated.The model is used to sim the longwall mining process, with advancement from left to right and a mining heig 8.8 m.Compression and tensile tests were carried out using PFC2D 5.0 software.Th lationship between rock macroscopic parameters and rock fine-grained parameters inverted to find the desired fine-grained parameters by inversion of the stratigraphic roscopic 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.

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.

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  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 cm 2 , 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 cm 2 , 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 cm 2 .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.
times the advance distance.When the working face advanced to 60-135 cm, th dimension was 1.325-1.38,the extension height of the crack was 88 cm, the rock area was 7949 cm 2 , and the height of the crack in the overlying strata was 0.07 ti advance distance.This is the stage of deformation and fracture extension of the ov strata.When the working face advanced to 142.5 cm, the rock damage area wa cm 2 .With the advance of the working face, the three stages of overburden fractur opment were divided mainly through the upward-stable-downward trend of th dimension of the overburden fracture.The overlying strata fracture height extens tern was divided into three stages: the overlying strata fracture formation stage, t lying strata deformation fracture extension stage, and the overlying strata fractur sion steady-stage.These three stages can better describe the disturbance of the ov strata and the initial formation of fractures in the early stage of workface advan This method can quantitatively explain the evolution law of overburden fracture actual mining process of the coal mine.

Ground Fissure Movement and Deformation Zone
In the early stage of mining, the main roof stratum is bent and sliding, and the has no obvious deformation characteristics.With the continuous advance of the w face, the overlying strata will eventually have stepped fault settlement, and the ov strata fault will penetrate the ground from bottom to bottom [31].Eventually, bo fissure zones and central fissure zones form.In the study of ground fissure dy ground fissures begin to crack, expand, contract, stabilize, or close.Ground fissur due to damage caused by gravity within the soil, resulting in central dynamic fissures experiencing a double cycle extension [2,32].In order to better simulate t phological changes of surface cracks caused by fracture expansion in the overlyin the development and expansion laws of overlying strata cracks in different minin were simulated by PFC.As can be seen from the force chain diagram, the entire rium structure of the overlying strata has been disrupted.The settlement of the b timately leads to ground deformation.It can be seen from the horizontal displac contour that the horizontal displacement of the surface changes unevenly.Finally, cracks, extrusion cracks, and step settlement cracks are formed.As shown in Fi through physical simulation experimental and numerical simulation studies of th lying strata migration, the main influencing factor is the key strata in the overlyin As the working face progressed, key strata were broken, leading to massive set displacement of the overlying strata and the development of ground fissures.T

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

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

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

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.

Figure 2
Figure 2 shows the borehole column diagram of the 12,401 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 12,401 working face in the Shangwan Mine are shown in Table1.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.

Figure 3 .
Figure 3. Contour map of loose layer thickness.

Figure 3 .
Figure 3. Contour map of loose layer thickness.

Figure 3 .
Figure 3. Contour map of loose layer thickness.

Figure 5 .
Figure 5. Contour map of thickness ratio of the loose layer to bedrock.

Figure 5 .
Figure 5. Contour map of thickness ratio of the loose layer to bedrock.

Figure 5 .
Figure 5. Contour map of thickness ratio of the loose layer to bedrock.

Figure 6 .
Figure 6.Test model: (a) related equipment, (b) similar material model of design.

Figure 6 .
Figure 6.Test model: (a) related equipment, (b) similar material model of design.

FractalFigure 7 .
Figure 7. Different advance distances of rock failure area.4.2.Simulation of the Caving Process of the 12,401 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

Figure 7 .
Figure 7. Different advance distances of rock failure area.

Figure 7 .
Figure 7. Different advance distances of rock failure area.

Figure 8 .
Figure 8. Overlying strata model of the 12,401 working face.

Figure 8 .
Figure 8. Overlying strata model of the 12401 working face.

Fractal
Fract.2023, 7, x FOR PEER REVIEW 9 of 16 150 m displacement-x contour 150 m contact force chain 200 m displacement-x contour 200 m contact force chain 250 m displacement-x contour 250 m contact force chain

Figure 9 .
Figure 9. Evolution of horizontal displacement clouds and force chains in overlying strata at different advance distances.

Figure 9 .
Figure 9. Evolution of horizontal displacement clouds and force chains in overlying strata at different advance distances.

Figure 10 .
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 .
Figure 10.Simulation experiment: (a) advance distance 67.5 cm, fracture network evolution, schematic diagram for calculating D f , (b) advance distance 120 cm, fracture network evolution, schematic diagram for calculating D f .

Figure 11 .
Figure 11.The characteristic curve of overlying strata failure.

Figure 12 .
Figure 12.Schematic diagram of overlying strata movement and ground fissure development and deformation.

Figure 12 .
Figure 12.Schematic diagram of overlying strata movement and ground fissure development and deformation.

Table 1 .
Statistics of borehole lithology around the No. 12401 working face.

Table 1 .
Statistics of borehole lithology around the No. 12401 working face.

Table 2 .
Physical and mechanical parameters of rocks in PFC model.

Table 3 .
Overlying strata failure results of different advance distances.