{'fraction_non_alphanumeric': 0.02493601854457263, 'fraction_numerical': 0.02324771207675095, 'mean_word_length': 4.910041178333861, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 9, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mechanical Properties and Failure Mechanism of the Weakly Cemented Overburden in Deep Mining

: With increases in the mining depth and area in the Ordos coal ﬁeld, the failure law of the super thick sandstone in the Zhidan group leads to frequent disasters, such as rock bursts and mine earthquakes, which have become a signiﬁcant issue, restricting large-scale continuous mining. To adequately understand the movement mechanism of the super-thick and weakly cemented overburden, and to promote the large-scale mining of the coal resources under it, this study analyzes the physical and mechanical properties, along with the microstructural characteristics, of the weakly cemented overburden of the Yingpanhao Coal Mine through mechanics tests, scanning electron microscope tests (SEM) and hydrolysis experiments. A two-dimensional discrete element model of the survey region is then built to explore the temporal and spatial evolution laws of the overburden failure. The results show that, even though poorly cemented strata such as the Cretaceous Zhidan group sandstone and the Zhiluo group sandstone are weak in lithology, their unique mineral composition and microstructural characteristics give them a greater rigidity when their thickness reaches a certain value. The surface subsidence exhibits a sudden increase, and the dynamic disaster range of the overlying strata is wide when deep multi-face mining was carried out under the super-thick and weakly cemented overburden. The temporal and spatial evolution laws of the strata subsidence and inﬂuence boundary are closely related to their depth, and their relationships evolve into the Boltzmann function and Boltzmann–parabolic function, respectively. The failure mode of the super-thick and weakly cemented overburden is ‘beam–arch shell–half arch shell’, and the failure boundary exhibits arch fractures.


Introduction
The Ordos Basin is one of the five energy bases of China, and its deep coal resources are of strategic significance for China's national energy development. However, the failure mechanism of its weakly cemented overburden differs from that of the weak rock and hard strata in other regions, which significantly limits the large-scale mining of deep coal resources in the area.
At present, most research focuses on the physical and mechanical properties of weakly cemented sandstone. For example, by carrying out SEM and mechanical tests, Li et al. found that both dissolution and cementation have a significant influence on a rock's strength [1]. By studying the influence of hydration on the physical and mechanical properties of weakly cemented sandstone in the Ordos Basin, Zhao et al. discovered that rock with increased stratification or joints will form pores and crevices owing to mineral shedding during water injection. The porosity can increase by up to 2.2 times, suggesting that the influence of hydration on the physical and mechanical properties of weakly cemented sandstone cannot be ignored [2]. After collecting samples of Zhidan group sandstone from the Xutuan and Zhaoji Coal Mines, Jin et al. observed their microstructure using an electron microscope. It was determined that this kind of rock has a lower compaction degree, weaker cementation, unconsolidated structure with poor conformity, and a higher volume of highly linked micron-sized porosity [3]. Park et al. used the discrete element method to assess the fracture mechanics behavior of weakly cemented sandstone and discussed variations in the fracture toughness with particle size, notch length, and specimen size [4]. Li et al. also collected the weakly cemented sandstone in Bultai Coal Mines and carried out uniaxial compression, triaxial compression, and Brazilian split tests. They discovered that the weakly cemented sandstone is characterized by a lower cement content and larger pores between the cemented sandstone particles. Further, weakly cemented sandstone is mainly fractured along the grain boundary, which is I-shaped from a macroscopic perspective [1]. Li et al. used X-ray diffraction instruments and environmental scanning electron microscopes to study the composition and microstructure of weakly cemented sandstone, finding that the compressive strength and the modulus of elasticity are much smaller than those of conventional sandstone. However, the peak strain was found to be 1.6 times that of conventional sandstone [5]. Zhang et al. carried out triaxial compression tests and creep tests on weakly cemented sandstone, finding that it has significant creep properties. Next, based on the test data, an improved Burgers creep model was built to describe the creep properties [6].
Significant advancements have been made in research on the failure mechanism of weakly cemented strata in western China. For example, Zhen et al. used the Wanli No. 1 Coal Mine as an example to analyze the surrounding rock stability of argillaceous cementation sandstone during large-span caving mining. In addition, according to the roadway supporting conditions, several reasonable and effective supporting countermeasures and key parameters have been proposed [7]. Sun et al. carried out crushing experiments and disintegration tests to study the physical and mechanical properties of weakly cemented overburdens, conducting a difference analysis of mid-east rocks. Then, the discrete element numerical model of the research area was formulated to study the failure mechanism of the weakly cemented overburden, ascertaining the distribution characteristics of the caving zone and its collaborative evolution law with the fracture zone of the large mining height work face [8]. Using the measured data, Lin found a smaller surface subsidence in the Yingpanhao Coal Mine and then explored the structural characteristics of the overburden based on the key stratum theory. They believed that the combined action of multi-key layers is one of the main reasons for the special strata movement [9]. By assessing the physical and mechanical parameters along with SEM tests, Gao found that the Tuha coalfield overburden has typical weak cementation characteristics such as low mechanical strength, poor cementation (primarily particle contact cementation), high porosity, and a loose structure. Furthermore, it was found that the development of water-conducting fracture zones in the weakly cemented overburden is still half-space shaped, and its development height conforms to a general empirical formula [10]. Zhang et al. took 21,103 fully mechanized caving workfaces in the No.4 Mine in Yili, Xinjiang, as an example to study the movement rule of the weakly cemented overburden. They found that when there is no obvious key strata structure in the weakly cemented overburden, the mining-induced fractures in the overburden develop upward synchronously with the advancement of the working face, which manifests as the temporal and spatial evolution laws of fracture development-main roof breakdown-crack closure [11]. Zhang et al. used digital photography to monitor the similarity material model and study the movement law of a super-thick and weakly cemented overburden, finding that the first breaking span of the immediate roof reached 240 m and that the cyclic fracturing length was about 60 m [12,13]. Lu et al. studied the mining thickness effect of the water-conducting fractured zone in thick coal seam mining through field coring and laboratory mechanical tests, finding that the fracture zone exhibited linear growth and the collapse zone presented a stepwise growth pattern [14]. Li et al. explored Minerals 2022, 12, 1276 3 of 30 a non-hydrophilic simulation material and carried out physical simulation tests to study the spatial distribution and dynamic evolution of the fractured zone in a weakly cemented overburden [15]. Du et al. used small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technology to study the surface subsidence law of deep mining under a weakly cemented overburden and found that the surface subsidence boundary followed a circular-parallelogram-trapezoidal shape [16]. Gong et al. established a refined numerical model to study the surface subsidence law of deep mining under the weakly cemented overburden and suggested that the surface was fully mined when the goaf length and width both reached about 3.3 h (where h is the average mining depth) [17]. Gong et al. also studied the setting method for vertical joint spacing in Universal Distinct Element Code (UDEC) and suggested that with increases in the vertical joint spacing, the maximum surface subsidence first drops, then rises, and reaches a minimum when the vertical joint spacing is set at approximately 7× the bed thickness [18]. Zhang et al. explored the surface movement rule and the influence of overburden characteristics on strata movement through field measurements and numerical simulation, suggesting that the Zhidan group sandstone has a strong control effect, whose first breakage results in the surface sinking in a fractured manner [19].
Based on the aforementioned research, it can be concluded that the physical and mechanical properties, along with the mineral composition of the weakly cemented sandstone, have been fully understood, consisting of the failure characteristics of the weakly cemented overburden in shallow coal seams. Some researchers have studied the overburden movement law in deep mining through field monitoring and numerical simulation, but little research has focused on the failure mechanism of high-level, super-thick, and weakly cemented overburdens of deep mining.
Therefore, based on our research on the physical and mechanical properties of weakly cemented sandstone in the Yingpanhao Coal Mine, this study builds a two-dimensional discrete element model of 22 mining areas to elucidate the movement rule and failure characteristics of super-thick and weakly cemented overburdens in deep mining. Finally, a theoretical basis is provided for strata movement control.

Mechanical Properties of Weakly Cemented Sandstone
To gain insights into the failure mechanism of super-thick and weakly cemented overburdens, we collected sandstone from the Zhidan group, Zhiluo group, and Yan'an formation from rock holes 1, 2, and 3 in the Yingpanhao Coal Mine. They were processed into multiple standard test blocks for mechanics experiments. The mechanical parameters of the partial sandstone are shown in Table 1, and pictures of partial sandstone samples are shown in Figure 1.
Minerals 2022, 12, x FOR PEER REVIEW 3 of 30 zone in thick coal seam mining through field coring and laboratory mechanical tests, finding that the fracture zone exhibited linear growth and the collapse zone presented a stepwise growth pattern [14]. Li et al. explored a non-hydrophilic simulation material and carried out physical simulation tests to study the spatial distribution and dynamic evolution of the fractured zone in a weakly cemented overburden [15]. Du et al. used small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technology to study the surface subsidence law of deep mining under a weakly cemented overburden and found that the surface subsidence boundary followed a circular-parallelogram-trapezoidal shape [16]. Gong et al. established a refined numerical model to study the surface subsidence law of deep mining under the weakly cemented overburden and suggested that the surface was fully mined when the goaf length and width both reached about 3.3 h (where h is the average mining depth) [17]. Gong et al. also studied the setting method for vertical joint spacing in Universal Distinct Element Code (UDEC) and suggested that with increases in the vertical joint spacing, the maximum surface subsidence first drops, then rises, and reaches a minimum when the vertical joint spacing is set at approximately 7× the bed thickness [18]. Zhang et al. explored the surface movement rule and the influence of overburden characteristics on strata movement through field measurements and numerical simulation, suggesting that the Zhidan group sandstone has a strong control effect, whose first breakage results in the surface sinking in a fractured manner [19]. Based on the aforementioned research, it can be concluded that the physical and mechanical properties, along with the mineral composition of the weakly cemented sandstone, have been fully understood, consisting of the failure characteristics of the weakly cemented overburden in shallow coal seams. Some researchers have studied the overburden movement law in deep mining through field monitoring and numerical simulation, but little research has focused on the failure mechanism of high-level, super-thick, and weakly cemented overburdens of deep mining.
Therefore, based on our research on the physical and mechanical properties of weakly cemented sandstone in the Yingpanhao Coal Mine, this study builds a two-dimensional discrete element model of 22 mining areas to elucidate the movement rule and failure characteristics of super-thick and weakly cemented overburdens in deep mining. Finally, a theoretical basis is provided for strata movement control.

Mechanical Properties of Weakly Cemented Sandstone
To gain insights into the failure mechanism of super-thick and weakly cemented overburdens, we collected sandstone from the Zhidan group, Zhiluo group, and Yan'an formation from rock holes 1, 2, and 3 in the Yingpanhao Coal Mine. They were processed into multiple standard test blocks for mechanics experiments. The mechanical parameters of the partial sandstone are shown in Table 1, and pictures of partial sandstone samples are shown in Figure 1.    As can be seen from Table 1 and the corresponding data regarding other types of weak cementation rock, the compressive strength of the weak cementation rock was determined to be between 10.10 and 44.02 MPa, the elastic modulus was between 777.77-9467.59 MPa, the Poisson's ratio was between 0.21 and 0.34, the cohesion was between 1.31 and 11.15 MPa, the internal friction angle was between 23 and 29 degree, and the density was between 1.93 and 2.55 g/cm 3 . Figure 2 helps analyze the lithologies of different strata. It can be seen that the compressive strength, elastic modulus, and cohesion of the coarse sandstone, medium sandstone, fine sandstone, and silt sandstone of the Cretaceous Zhidan group were all found to be smaller than those of the corresponding strata of the Jurassic Zhiluo group. In the Cretaceous Zhidan and Jurassic Zhiluo groups, the compressive strength and elastic modulus of coarse sandstone, medium sandstone, fine sandstone, and silt sandstone were reduced accordingly. The cohesion of coarse sandstone, medium sandstone, fine sandstone, and siltstone decreased in the Jurassic Zhiluo group but increased in the Cretaceous Zhidan group. As can be seen from Table 1 and the corresponding data regarding other types weak cementation rock, the compressive strength of the weak cementation rock was d termined to be between 10.10 and 44.02 MPa, the elastic modulus was between 777.7 9467.59 MPa, the Poisson's ratio was between 0.21 and 0.34, the cohesion was betwe 1.31 and 11.15 MPa, the internal friction angle was between 23 and 29 degree, and t density was between 1.93 and 2.55 g/cm 3 . Figure 2 helps analyze the lithologies of different strata. It can be seen that the co pressive strength, elastic modulus, and cohesion of the coarse sandstone, medium san stone, fine sandstone, and silt sandstone of the Cretaceous Zhidan group were all fou to be smaller than those of the corresponding strata of the Jurassic Zhiluo group. In t Cretaceous Zhidan and Jurassic Zhiluo groups, the compressive strength and elastic mo ulus of coarse sandstone, medium sandstone, fine sandstone, and silt sandstone were duced accordingly. The cohesion of coarse sandstone, medium sandstone, fine sandston and siltstone decreased in the Jurassic Zhiluo group but increased in the Cretaceous Z dan group. 14  The rock microstructure is the reason why the rock mechanics properties of the abo strata differ from what is conventionally known. In the following sections, the microstru tural characteristics of weakly cemented sandstone are studied through detailed expe ments.

Physical Structural Characteristics of Weakly Cemented Sandstone
In this paper, some sandstone samples of the Zhidan group, Anding group, Zhil group, and Yan'an group in the Yingpanhao Coal Mine were collected. Then, scanni electron microscope (SEM) tests and mineral composition identification tests were carri out.

SEM Tests of Weakly Cemented Sandstone
To clearly elucidate the crack development in rock samples, they were polished usi an Argon-ion Polisher-697. To obtain clear SEM images, samples were carbon coated. F ure 3 shows the equipment used (SIGMA, Carl Zeiss AG, Oberkochen, Germany).  The rock microstructure is the reason why the rock mechanics properties of the above strata differ from what is conventionally known. In the following sections, the microstructural characteristics of weakly cemented sandstone are studied through detailed experiments.

Physical Structural Characteristics of Weakly Cemented Sandstone
In this paper, some sandstone samples of the Zhidan group, Anding group, Zhiluo group, and Yan'an group in the Yingpanhao Coal Mine were collected. Then, scanning electron microscope (SEM) tests and mineral composition identification tests were carried out.

SEM Tests of Weakly Cemented Sandstone
To clearly elucidate the crack development in rock samples, they were polished using an Argon-ion Polisher-697. To obtain clear SEM images, samples were carbon coated. Figure 3 shows the equipment used (SIGMA, Carl Zeiss AG, Oberkochen, Germany).  The rock microstructure is the reason why the rock mechanics properties of the above strata differ from what is conventionally known. In the following sections, the microstructural characteristics of weakly cemented sandstone are studied through detailed experiments.

Physical Structural Characteristics of Weakly Cemented Sandstone
In this paper, some sandstone samples of the Zhidan group, Anding group, Zhiluo group, and Yan'an group in the Yingpanhao Coal Mine were collected. Then, scanning electron microscope (SEM) tests and mineral composition identification tests were carried out.

SEM Tests of Weakly Cemented Sandstone
To clearly elucidate the crack development in rock samples, they were polished using an Argon-ion Polisher-697. To obtain clear SEM images, samples were carbon coated.  However, in the argon-ion polishing experiment, the Zhidan group sandstone sample was always uneven and even broken locally. This phenomenon is closely related to its particle size, void, compactness and cement. Therefore, SEM proceeded with untreated sandstone samples (Figure 3d).  However, in the argon-ion polishing experiment, the Zhidan group sandstone sample was always uneven and even broken locally. This phenomenon is closely related to its particle size, void, compactness and cement. Therefore, SEM proceeded with untreated sandstone samples (Figure 3d). show the geometrical and microscopic characteristics of Anding formation sandstone, Zhiluo formation sandstone, and Yan'an formation sandstone, respectively.

Zoom
As shown in Figure 4, the grain size of Zhidan group sandstone was generally large, and could be divided into coarse (500-1000 μm) and medium sandstone (250-500 μm), according to the different levels of grain size. The pores were relatively large, some of which had a width of 603 μm. Moreover, the rocks were relatively loose, and some pores were filled with fine sandstone. Further, no fractures developed. Note that the test sample of Zhidan group sandstone was not polished, and its conductivity was not good. Under the irradiation of the electron beam, a discharge was produced, and white bright spots and patches appeared on the image.  As shown in Figure 5, slight, gently inclined and lateral stratification developed in the Zhidan group sandstone, while vertical stratification did not develop. Therefore, it was susceptible to transverse shear failure along the stratification, forming a flaky structure. As can be seen in Figure 6, the Anding formation sandstone had a small grain size, comprising mostly siltstone (5-50 μm), according to the grain size, mixed with a little fine sand (50-250 μm). Moreover, obvious fractures and voids developed.  As can be seen in Figure 6, the Anding formation sandstone had a small grain size comprising mostly siltstone (5-50 μm), according to the grain size, mixed with a little fine sand (50-250 μm). Moreover, obvious fractures and voids developed. There are three kinds of pores in rock: pores, holes and cracks. Pores and holes are spaces that develop in three dimensions and vary in size. A hole diameter is greater than 2 mm while that of a pore is less than 2 mm. A crack is a gap which develops in two dimensions and the difference between a crack and a hole is morphological. The ratio of the major and minor axes in a crack is larger than 10, while that of a hole is less than 10.
As shown in Figure 7, the Zhiluo group sandstone had a moderate grain size, com prising mostly fine sandstone, with some medium sandstone. The rock was relatively There are three kinds of pores in rock: pores, holes and cracks. Pores and holes are spaces that develop in three dimensions and vary in size. A hole diameter is greater than 2 mm while that of a pore is less than 2 mm. A crack is a gap which develops in two dimensions, and the difference between a crack and a hole is morphological. The ratio of the major and minor axes in a crack is larger than 10, while that of a hole is less than 10.
As shown in Figure 4, the grain size of Zhidan group sandstone was generally large, and could be divided into coarse (500-1000 µm) and medium sandstone (250-500 µm), according to the different levels of grain size. The pores were relatively large, some of which had a width of 603 µm. Moreover, the rocks were relatively loose, and some pores were filled with fine sandstone. Further, no fractures developed. Note that the test sample of Zhidan group sandstone was not polished, and its conductivity was not good. Under the irradiation of the electron beam, a discharge was produced, and white bright spots and patches appeared on the image.
As shown in Figure 5, slight, gently inclined and lateral stratification developed in the Zhidan group sandstone, while vertical stratification did not develop. Therefore, it was susceptible to transverse shear failure along the stratification, forming a flaky structure. dense, with no large fracture, stratification, or joint development; holes only developed rarely. As shown in Figure 7e, the cementation contained relatively high amounts of Si and Al, as well as a quantifiable amount of Mg and a small amount of Fe, K, and Ca. This mineral may be plagioclase. As can be seen in Figure 8, the rock particles were closely bound to the cements in the Yan'an formation sandstone. The irregular material may have been a mineral matrix or cementation agent, mostly comprising fine sandstone, according to the grain size classification. No large fractures, stratification, and joints developed, and holes developed rarely.

Cementation Types and Mineral Analysis of Weakly Cemented Sandstone
Before mineral composition analysis, the larger block rock sample was crushed, ground into powder by a crusher, and filtered using a 300-mesh filter screen to obtain more than 10 g of powdered rock, which was the amount used for the two mineral analysis experiments. The experimental equipment and process are shown in Figure 9. As can be seen in Figure 6, the Anding formation sandstone had a small grain size, comprising mostly siltstone (5-50 µm), according to the grain size, mixed with a little fine sand (50-250 µm). Moreover, obvious fractures and voids developed.
As shown in Figure 7, the Zhiluo group sandstone had a moderate grain size, comprising mostly fine sandstone, with some medium sandstone. The rock was relatively dense, with no large fracture, stratification, or joint development; holes only developed rarely. As shown in Figure 7e, the cementation contained relatively high amounts of Si and Al, as well as a quantifiable amount of Mg and a small amount of Fe, K, and Ca. This mineral may be plagioclase.
As can be seen in Figure 8, the rock particles were closely bound to the cements in the Yan'an formation sandstone. The irregular material may have been a mineral matrix or cementation agent, mostly comprising fine sandstone, according to the grain size classification. No large fractures, stratification, and joints developed, and holes developed rarely.

Cementation Types and Mineral Analysis of Weakly Cemented Sandstone
Before mineral composition analysis, the larger block rock sample was crushed, ground into powder by a crusher, and filtered using a 300-mesh filter screen to obtain more than 10 g of powdered rock, which was the amount used for the two mineral analysis experiments. The experimental equipment and process are shown in Figure 9.

Cementation Types and Mineral Analysis of Weakly Cemented Sandstone
Before mineral composition analysis, the larger block rock sample was crushed, ground into powder by a crusher, and filtered using a 300-mesh filter screen to obtain more than 10 g of powdered rock, which was the amount used for the two mineral analysis experiments. The experimental equipment and process are shown in Figure 9. Then, the mineral composition was determined using a mineral composition analyzer (AXiOX MAX, PANalytical B.V., Almelo, Netherlands), the identification results of which are shown in Table 2. As can be seen in Table 2, the clay mineral content of the Zhidan group and Yan'an group sandstones was less than that of the Anding group and Zhiluo group sandstones. The Zhidan group sandstone was dominated by quartz, accounting for 54.1%. The Yan'an formation sandstone was dominated by dolomite, accounting for 42.9%. Subsequently, the mineral compositions of the energy spectrum of each rock sample mineral and cementation were analyzed, and a quantitative analysis of rock sample mineral compositions was performed.
As can be seen in Figure 10, the mineral contained large amounts of Si and Al, with a certain amount of K. Based on the mineral composition percentage of the Anding Formation sandstone in Table 2 and the element proportion contained in each mineral, it could be determined that this mineral was potash feldspar.
As shown in Figure 11, there were large amounts of Fe and Si, with smaller amounts of K, in the cementation component. According to the mineral composition percentage in Table 2 and the fuchsia appearance of the Anding Formation sandstone, it could be confirmed that this mineral was hematite.
In Figures 12 and 13, large amounts of Ca can be identified, with a certain amount of Mn or Fe in the cementation. However, the cement in Figure 12 does not contain Mg, while the cementation in Figure 13 does. By comparing the mineral composition percentage of the Zhidan sandstone group in Table 2 and the element proportion contained in each mineral, it can be determined that the cementation in Figure 12 is calcite and the cementation in Figure 13 is dolomite. the mineral compositions of the energy spectrum of each rock sample mineral and cementation were analyzed, and a quantitative analysis of rock sample mineral compositions was performed.
As can be seen in Figure 10, the mineral contained large amounts of Si and Al, with a certain amount of K. Based on the mineral composition percentage of the Anding Formation sandstone in Table 2 and the element proportion contained in each mineral, it could be determined that this mineral was potash feldspar. As shown in Figure 11, there were large amounts of Fe and Si, with smaller amounts of K, in the cementation component. According to the mineral composition percentage in Table 2 and the fuchsia appearance of the Anding Formation sandstone, it could be confirmed that this mineral was hematite. In Figures 12 and 13, large amounts of Ca can be identified, with a certain amount of Mn or Fe in the cementation. However, the cement in Figure 12 does not contain Mg, while the cementation in Figure 13 does. By comparing the mineral composition percentage of the Zhidan sandstone group in Table 2 and the element proportion contained in each mineral, it can be determined that the cementation in Figure 12 is calcite and the cementation in Figure 13 is dolomite.  In Figures 12 and 13, large amounts of Ca can be identified, with a certain amount of Mn or Fe in the cementation. However, the cement in Figure 12 does not contain Mg, while the cementation in Figure 13 does. By comparing the mineral composition percentage of the Zhidan sandstone group in Table 2 and the element proportion contained in each mineral, it can be determined that the cementation in Figure 12 is calcite and the cementation in Figure 13 is dolomite.  As shown in Figure 14, the rock particles contained large amounts of Si and a small amount of Al. Using the mineral composition percentage of the Zhidan group sandstone in Table 2 and the element proportion contained in each mineral, the mineral could be identified as quartz. As shown in Figure 14, the rock particles contained large amounts of Si and a small amount of Al. Using the mineral composition percentage of the Zhidan group sandstone in Table 2 and the element proportion contained in each mineral, the mineral could be identified as quartz. As shown in Figure 15, there was a high Si content in the rock particles. This mineral could be recognized as quartz according to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral. As shown in Figure 16, the rock particles contained relatively large amounts of Si and Al, as well as a quantifiable amount of Mg and small amounts of Fe, K, and Ca. According to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral, this mineral could be determined to be plagioclase. As shown in Figure 15, there was a high Si content in the rock particles. This mineral could be recognized as quartz according to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral. As shown in Figure 15, there was a high Si content in the rock particles. This mineral could be recognized as quartz according to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral. As shown in Figure 16, the rock particles contained relatively large amounts of Si and Al, as well as a quantifiable amount of Mg and small amounts of Fe, K, and Ca. According to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral, this mineral could be determined to be plagioclase. As shown in Figure 16, the rock particles contained relatively large amounts of Si and Al, as well as a quantifiable amount of Mg and small amounts of Fe, K, and Ca. According to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral, this mineral could be determined to be plagioclase.
As shown in Figure 16, the rock particles contained relatively large amounts of Si and Al, as well as a quantifiable amount of Mg and small amounts of Fe, K, and Ca. According to the mineral composition percentage of the Zhiluo group sandstone in Table 2 and the element proportion contained in each mineral, this mineral could be determined to be plagioclase.   Figure 17 shows that the rock particles contained large amounts of Si and Al and a certain amount of K. According to the mineral composition percentage of the Yan'an formation sandstone in Table 2 and the element proportion contained in each mineral, it could be determined that this mineral was potash feldspar.
Minerals 2022, 12, x FOR PEER REVIEW 13 of 30 Figure 17 shows that the rock particles contained large amounts of Si and Al and a certain amount of K. According to the mineral composition percentage of the Yan'an formation sandstone in Table 2 and the element proportion contained in each mineral, it could be determined that this mineral was potash feldspar. As shown in Figure 18, the rock particles had high Si contents. According to the mineral composition percentage of Yan'an Formation sandstone in Table 2 and the element proportion contained in each mineral, it was determined to be quartz.

Hydrolytic Properties of Weakly Cemented Sandstone
Some of the Zhidan group sandstone, Anding group sandstone, Zhiluo group sandstone, and Yan'an group sandstone were placed in a 5000 mL container. The container was filled with water until the experimental samples were submerged, and sealed with plastic wrap. In the hydrolysis experiment, there was no artificial shaking of the container, As shown in Figure 18, the rock particles had high Si contents. According to the mineral composition percentage of Yan'an Formation sandstone in Table 2 and the element proportion contained in each mineral, it was determined to be quartz.  Figure 17 shows that the rock particles contained large amounts of Si and Al and a certain amount of K. According to the mineral composition percentage of the Yan'an formation sandstone in Table 2 and the element proportion contained in each mineral, it could be determined that this mineral was potash feldspar. As shown in Figure 18, the rock particles had high Si contents. According to the mineral composition percentage of Yan'an Formation sandstone in Table 2 and the element proportion contained in each mineral, it was determined to be quartz.

Hydrolytic Properties of Weakly Cemented Sandstone
Some of the Zhidan group sandstone, Anding group sandstone, Zhiluo group sandstone, and Yan'an group sandstone were placed in a 5000 mL container. The container was filled with water until the experimental samples were submerged, and sealed with

Hydrolytic Properties of Weakly Cemented Sandstone
Some of the Zhidan group sandstone, Anding group sandstone, Zhiluo group sandstone, and Yan'an group sandstone were placed in a 5000 mL container. The container was filled with water until the experimental samples were submerged, and sealed with plastic wrap. In the hydrolysis experiment, there was no artificial shaking of the container, which means that the whole experiment process occurred in a stable state. Figures 19-22 show comparative diagrams of the corresponding rock samples before and after hydrolysis.  Figure 20 is a hydrolysis experiment comparison diagram of the Zhidan group sandstone. When the hydrolysis experiment was carried out, it was found that the high porosity of the immediately became saturated with the fluid (Note that the Zhidan group sandstones have been on site for some time and are not in their natural state), and some phenomena-such as bubbles and creaking noises-were transiently observed (Figure 20b). After one month of the hydrolysis experiment, although a layer of tiny particles was deposited at the bottom of the vessel, the overall structure of the Zhidan group sandstone was relatively complete (Figure 20c). This was because there were large gaps between particles in the rock mass; therefore, water absorption occurred. In addition, the Zhidan group sandstone contained few clay minerals, accounting for only 8.5% of the total sample amount. The cements were mostly potash feldspar and plagioclase, which are poorly soluble in water. Thus, the Zhidan group sandstone will not collapse.   (Figure 21b). Several days later, large fractures developed in the rock sample; however, disintegration did not occur, and the rock sample almost maintained its original geometry (Figure 21c). This was because the Zhiluo group sandstone contains clay minerals, which are mostly insoluble. After encountering water, the Zhiluo group sandstone was only destroyed along the tiny fractures and cavities in the rock mass. Thus, the original geometry of the rock sample was maintained. interlayer positions), which provided the pressure for tensile failure of the individual components. Then, the rock sample expanded and failed along the internal pores and fractures, leading to the disintegration of the rock sample.  Figure 20 is a hydrolysis experiment comparison diagram of the Zhidan group sandstone. When the hydrolysis experiment was carried out, it was found that the high porosity of the immediately became saturated with the fluid (Note that the Zhidan group sandstones have been on site for some time and are not in their natural state), and some phenomena-such as bubbles and creaking noises-were transiently observed (Figure 20b). After one month of the hydrolysis experiment, although a layer of tiny particles was deposited at the bottom of the vessel, the overall structure of the Zhidan group sandstone was relatively complete (Figure 20c). This was because there were large gaps between particles in the rock mass; therefore, water absorption occurred. In addition, the Zhidan group sandstone contained few clay minerals, accounting for only 8.5% of the total sample amount. The cements were mostly potash feldspar and plagioclase, which are poorly soluble in water. Thus, the Zhidan group sandstone will not collapse.   (Figure 21b). Several days later, large fractures developed in the rock sample; however, disintegration did not occur, and the rock sample almost maintained its original geometry (Figure 21c). This was because the Zhiluo group sandstone contains clay minerals, which are mostly insoluble. After encountering water, the Zhiluo group sandstone was only destroyed along the tiny fractures and cavities in the rock mass. Thus, the original geometry of the rock sample was maintained.   (Figure 22c). This was because the clay mineral content in Yan'an Formation sandstone is low, and it is poorly soluble in water. In addition, holes and fractures almost did not develop, and the rock samples maintained a good overall structure. hydrolysis. (c) After hydrolysis. Figure 22 is a hydrolysis experiment comparison diagram of the Yan'an formation sandstone. No abnormality was detected during the hydrolysis experiments on this sandstone (Figure 22c). This was because the clay mineral content in Yan'an Formation sandstone is low, and it is poorly soluble in water. In addition, holes and fractures almost did not develop, and the rock samples maintained a good overall structure. By analyzing the physical and mechanical properties of weakly cemented sandstones, it was found that most sandstones in the Zhidan group were medium sandstone or coarse sandstone, with a relatively large particle size, large internal voids, calcareous cementation, non-disintegration in the presence of water, and relatively weak lithology. Most Anding Formation sandstones were sandy mudstones of small particle size with internal fissures and holes. Additionally, they easily disintegrated and had a weak lithology. Moreover, the cementation was clay minerals or hematite. The Zhiluo group sandstone was mostly sandy mudstone, with a moderate particle size, relatively dense rocks, clay mineral cementation, non-disintegration in the presence of water, and weak mediumhard lithology. The Yan'an Formation sandstones were mostly fine sandstones or siltstones, with dense rocks, ferruginous cementation, non-disintegrating in the presence of water, and a medium-hard lithology.

Strata Failure Rule of Deep Mining with the Super-Thick and Weakly Cemented Overburden through Numerical Simulation
According to the geological and mining conditions of Yingpanhao Coal Mine, a twodimensional discrete element model of the study area was constructed. According to the physical and mechanical properties and micro-geometric characteristics of weakly cemented sandstone, the mechanical parameters of rock joints and stratification were adjusted until the simulation results were consistent with the measured data. The model size  Figure 19 is a hydrolysis experiment comparison diagram of the Anding Formation sandstone. When the hydrolysis experiment of the Anding Formation sandstone was being conducted, the water near the rock sample gradually became turbid, and fine particles separated from the overall structure of the rock sample and continuously fell to the bottom of the container (Figure 19b). After several hours, the overall structure of the rock sample broke down and formed flake-or block-like structures (Figure 19c). This was because the clay mineral content of the Anding Formation sandstone was up to 41.3%, and the pores and fractures developed internally. Some of the clays may have expanded (water in interlayer positions), which provided the pressure for tensile failure of the individual components. Then, the rock sample expanded and failed along the internal pores and fractures, leading to the disintegration of the rock sample. Figure 20 is a hydrolysis experiment comparison diagram of the Zhidan group sandstone. When the hydrolysis experiment was carried out, it was found that the high porosity of the immediately became saturated with the fluid (Note that the Zhidan group sandstones have been on site for some time and are not in their natural state), and some phenomena-such as bubbles and creaking noises-were transiently observed (Figure 20b). After one month of the hydrolysis experiment, although a layer of tiny particles was deposited at the bottom of the vessel, the overall structure of the Zhidan group sandstone was relatively complete (Figure 20c). This was because there were large gaps between particles in the rock mass; therefore, water absorption occurred. In addition, the Zhidan group sandstone contained few clay minerals, accounting for only 8.5% of the total sample amount. The cements were mostly potash feldspar and plagioclase, which are poorly soluble in water. Thus, the Zhidan group sandstone will not collapse. Figure 21 is a hydrolysis experiment comparison diagram of the Zhiluo group sandstone. No abnormal phenomena occurred during the hydrolysis experiment of the Zhiluo group sandstone (Figure 21b). Several days later, large fractures developed in the rock sample; however, disintegration did not occur, and the rock sample almost maintained its original geometry (Figure 21c). This was because the Zhiluo group sandstone contains clay minerals, which are mostly insoluble. After encountering water, the Zhiluo group sandstone was only destroyed along the tiny fractures and cavities in the rock mass. Thus, the original geometry of the rock sample was maintained. Figure 22 is a hydrolysis experiment comparison diagram of the Yan'an formation sandstone. No abnormality was detected during the hydrolysis experiments on this sandstone (Figure 22c). This was because the clay mineral content in Yan'an Formation sandstone is low, and it is poorly soluble in water. In addition, holes and fractures almost did not develop, and the rock samples maintained a good overall structure.
By analyzing the physical and mechanical properties of weakly cemented sandstones, it was found that most sandstones in the Zhidan group were medium sandstone or coarse sandstone, with a relatively large particle size, large internal voids, calcareous cementation, non-disintegration in the presence of water, and relatively weak lithology. Most Anding Formation sandstones were sandy mudstones of small particle size with internal fissures and holes. Additionally, they easily disintegrated and had a weak lithology. Moreover, the cementation was clay minerals or hematite. The Zhiluo group sandstone was mostly sandy mudstone, with a moderate particle size, relatively dense rocks, clay mineral cementation, non-disintegration in the presence of water, and weak medium-hard lithology. The Yan'an Formation sandstones were mostly fine sandstones or siltstones, with dense rocks, ferruginous cementation, non-disintegrating in the presence of water, and a medium-hard lithology.

Strata Failure Rule of Deep Mining with the Super-Thick and Weakly Cemented Overburden through Numerical Simulation
According to the geological and mining conditions of Yingpanhao Coal Mine, a twodimensional discrete element model of the study area was constructed. According to the physical and mechanical properties and micro-geometric characteristics of weakly cemented sandstone, the mechanical parameters of rock joints and stratification were adjusted until the simulation results were consistent with the measured data. The model size and parameters are detailed in Tables 3 and 4. In this paper, we studied the temporal and spatial evolution laws of the overburden movement when workfaces 2201, 2202, 2203, and 2204 were mined.  Table 4. Partial parameters of the numerical simulation model. As can be seen from Figure 23, when the width of working face was 300 m, the height of the water-conducting fracture zone developed to 116 m above the coal seam, which was consistent with the measured results (the height of the water-conducting fracture zone develops to 115 m above the coal seam). At this time, the surface subsidence of 407 mm was larger than the actual measured value (326 mm). This was because the model built by UDEC was two-dimensional, the advance direction could be regarded as infinite mining, and the phenomenon of the inter-embedding among blocks occurred, so that the simulation results were larger than the measured results under the same mining width. This shows that the numerical model established in this paper is reasonable.

Displacement Evolution of Overburden Strata
In this section, the authors only discuss the displacement nephogram and the movement curves of the overlying strata at different buried depths when workfaces 2201, 2202, 2203, and 2204 were 300 m in length, as shown in Figures 24-27.

Displacement Evolution of Overburden Strata
In this section, the authors only discuss the displacement nephogram and the movement curves of the overlying strata at different buried depths when workfaces 2201, 2202, 2203, and 2204 were 300 m in length, as shown in Figures 24-27.

Displacement Evolution of Overburden Strata
In this section, the authors only discuss the displacement nephogram and the movement curves of the overlying strata at different buried depths when workfaces 2201, 2202, 2203, and 2204 were 300 m in length, as shown in Figures 24-27.    It can be seen from these figures that after the coal seam had been mined, the subsidence was more significant from the immediate roof to 66 m above the coal seam. The separation layers of workforce 2201, 2202, 2203, and 2204 developed to 96 m, 66 m, 66 m, and 86 m above the coal seam, respectively.
With the expansion of the mining range, the compression of coal pillar 1, coal pillar 2 and coal pillar 3 exhibited a decreasing trend, which was consistent with the simulation results of similar materials. In Figure 28, Pillar 1 fractured and pressed into the floor strata, whereas Pillar 2 and 3 fractured the floor strata. It can be seen from these figures that after the coal seam had been mined, the subsidence was more significant from the immediate roof to 66 m above the coal seam. The separation layers of workforce 2201, 2202, 2203, and 2204 developed to 96 m, 66 m, 66 m, and 86 m above the coal seam, respectively.
With the expansion of the mining range, the compression of coal pillar 1, coal pillar 2 and coal pillar 3 exhibited a decreasing trend, which was consistent with the simulation results of similar materials. In Figure 28, Pillar 1 fractured and pressed into the floor strata, whereas Pillar 2 and 3 fractured the floor strata.
With increases in the mining range, the compression of coal pillar 1, coal pillar 2, and coal pillar 3 exhibited a decreasing trend, which was consistent with the simulation results for similar materials. Influenced by the coal pillar, the subsidence curve of the overburden rock was wavy, which was weakened when it passed through the thick Zhiluo group sandstone and, thus, gradually became a single gentle subsidence basin. To demonstrate the relationship between the surface subsidence and mining degree more intuitively, the surface subsidence values of different width-depth ratios were found and the corresponding surface subsidence coefficients were calculated, as shown in Table 5. In addition, the corresponding surface subsidence curves and the relationship curves between the width-depth ratio and surface subsidence coefficient were drawn, as shown in Figure 29.  With increases in the mining range, the compression of coal pillar 1, coal pillar 2, and coal pillar 3 exhibited a decreasing trend, which was consistent with the simulation results for similar materials. Influenced by the coal pillar, the subsidence curve of the overburden rock was wavy, which was weakened when it passed through the thick Zhiluo group sandstone and, thus, gradually became a single gentle subsidence basin.
To demonstrate the relationship between the surface subsidence and mining degree more intuitively, the surface subsidence values of different width-depth ratios were found and the corresponding surface subsidence coefficients were calculated, as shown in Table 5. In addition, the corresponding surface subsidence curves and the relationship curves between the width-depth ratio and surface subsidence coefficient were drawn, as shown in Figure 29.  As shown in Figure 29a, after mining workfaces 2201, 2202, 2203, and 2204, the corresponding maximum surface subsidence values were 407 mm, 1910 mm, 4750 mm, and 5450 mm, respectively, exhibiting obvious jumps. When the length of workface 2203 was 120 m, the surface subsidence was the most significant, and the maximum surface subsidence value increased by 1090 mm. The numerical simulation results were consistent with the simulation results of similar materials, further verifying the reliability of the results found in this paper.
From Figure 29b, it can be seen that with increases in the width-depth ratio, the surface subsidence coefficient increased as a Boltzmann function, and the correlation coefficient R 2 = 0.998. The mathematical expression is as follows: In addition, the subsidence values of the overlying strata at different depths and the distance of the subsidence boundary to the goaf boundary with increases in the mining space were counted, and the corresponding relationship curves were drawn and fitted, as shown in Tables 6 and 7 and Figures 30 and 31: As shown in Figure 29a, after mining workfaces 2201, 2202, 2203, and 2204, the corresponding maximum surface subsidence values were 407 mm, 1910 mm, 4750 mm, and 5450 mm, respectively, exhibiting obvious jumps. When the length of workface 2203 was 120 m, the surface subsidence was the most significant, and the maximum surface subsidence value increased by 1090 mm. The numerical simulation results were consistent with the simulation results of similar materials, further verifying the reliability of the results found in this paper.
From Figure 29b, it can be seen that with increases in the width-depth ratio, the surface subsidence coefficient increased as a Boltzmann function, and the correlation coefficient R 2 = 0.998. The mathematical expression is as follows: In addition, the subsidence values of the overlying strata at different depths and the distance of the subsidence boundary to the goaf boundary with increases in the mining space were counted, and the corresponding relationship curves were drawn and fitted, as shown in Tables 6 and 7 and Figures 30 and 31:   3  340  370  340  340  33  400  410  370  370  43  410  410  380  370  66  430  420  380  380  186  450  430  390  380  236  470  430  380  380  276  470  440  380  380  312  480  450  390  390  612  550  470  410  410  725  600  520  420  420 Minerals 2022, 12, x FOR PEER REVIEW  20 of 30   66  430  420  380  380  186  450  430  390  380  236  470  430  380  380  276  470  440  380  380  312  480  450  390  390  612  550  470  410  410  725  600  520  420  420 It can be seen from Figure 30 that, when workfaces 2201, 2202, 2203, and 2204 were mined, the maximum subsidence value of the overlying strata at different buried depths and the height above the coal seam exhibited a Boltzmann function relationship, and the corresponding correlation coefficient R 2 = 0.999, 0.994, 0.956, and 0.944. With mining range expansion, the correlation coefficient of the Boltzmann function tended to decrease. As shown in Figure 31, when workfaces 2201, 2202, 2203, and 2204 were mined, the distance of the subsidence influence boundary to the goaf boundary at different buried depths and the height above the coal seam (∆H) had a Boltzmann function, with a correlation coefficient R 2 = 0.992, 0.967, 0.872, and 0.918. The relevant mathematical expression is as follows:  By analyzing Figures 30 and 31, it can be seen that when the goaf width was small, there was a Boltzmann function relationship between the Wmax at different buried depths and the height above the coal seam; the same was true for the SW at different buried depths and the height above the coal seam. With the continuous expansion of the goaf, the functional relationship between the Wmax at different buried depths and the height above the coal seam was constant, while the functional relationship between the SW at different buried depths gradually changed from a Boltzmann function to a parabolic function. Furthermore, due to the continuously expanding scope of the goaf, the tangent of the main influence angle increased gradually, and its variation range was between 1.21 and 1.71, which was within the common value range.

Evolution Law of the Stress Fields of the Overlying Strata
This study extracted the stress evolution data of overlying strata. Figure 32 shows a schematic of the evolution law of the surrounding rock stress field, and Table 8 reflects the corresponding vertical stress value statistics. The specific analysis was as follows: (1) When the lengths of workface 2201 were 180 m, 240 m, and 300 m, the vertical stress of the coal wall at coal pillar 1 gradually increased, reaching 28.6 MPa. It can be seen from Figure 30 that, when workfaces 2201, 2202, 2203, and 2204 were mined, the maximum subsidence value of the overlying strata at different buried depths and the height above the coal seam exhibited a Boltzmann function relationship, and the corresponding correlation coefficient R 2 = 0.999, 0.994, 0.956, and 0.944. With mining range expansion, the correlation coefficient of the Boltzmann function tended to decrease.
As shown in Figure 31, when workfaces 2201, 2202, 2203, and 2204 were mined, the distance of the subsidence influence boundary to the goaf boundary at different buried depths and the height above the coal seam (∆H) had a Boltzmann function, with a correlation coefficient R 2 = 0.992, 0.967, 0.872, and 0.918. The relevant mathematical expression is as follows: By analyzing Figures 30 and 31, it can be seen that when the goaf width was small, there was a Boltzmann function relationship between the W max at different buried depths and the height above the coal seam; the same was true for the S W at different buried depths and the height above the coal seam. With the continuous expansion of the goaf, the functional relationship between the W max at different buried depths and the height above the coal seam was constant, while the functional relationship between the S W at different buried depths gradually changed from a Boltzmann function to a parabolic function. Furthermore, due to the continuously expanding scope of the goaf, the tangent of the main influence angle increased gradually, and its variation range was between 1.21 and 1.71, which was within the common value range.

Evolution Law of the Stress Fields of the Overlying Strata
This study extracted the stress evolution data of overlying strata. Figure 32 shows a schematic of the evolution law of the surrounding rock stress field, and Table 8 reflects the corresponding vertical stress value statistics. The specific analysis was as follows: (3) When mining workface 2203, the vertical stress of coal pillar 1 gradually d creased to 6.1 MPa, and the phenomenon of strata behavior appeared at the immedi roof of workface 2202; the vertical stress increased to 27.7 MPa. Moreover, the verti stress of coal pillar 2 increased gradually, and the vertical stress values were 39.4 M 39.9 MPa, 46.9 MPa, and 50.80 MPa, in sequence. In addition, the vertical stress values the coal wall at pillar 3 first increased and then decreased, at 25.7 MPa, 33.1 MPa, 4 MPa, and 30.3 MPa, in sequence.
(4) When mining workface 2204, the vertical stress value of coal pillar 1 increas rapidly, from 7.01 MPa to 101 MPa, 102 MPa, and 103 MPa. Moreover, the vertical str of coal pillar 2 and coal pillar 3 increased gradually, and the vertical stress at the imme ate roof of the working face 2201 remained almost unchanged. In addition, slight str behavior manifested at workfaces 2202 and 2203, and the vertical stress gradually creased.
By analyzing the stress evolution law of coal pillar 1, it can be seen that the dynam range of the super-thick and weakly cemented overburden was wide during deep minin For example, the stress concentration degree of the segment coal pillar at 600 m away fro workface 2204 was further increased.    (1) When the lengths of workface 2201 were 180 m, 240 m, and 300 m, the vertical stress of the coal wall at coal pillar 1 gradually increased, reaching 28.6 MPa.
(2) When mining workface 2202, the vertical stress of coal pillar 1 continued to increase, and the stress value reached 67.8 MPa. The stress of the coal wall at pillar 2 on both sides of the goaf increased, and the vertical stress value increased to 42.5 MPa.
(3) When mining workface 2203, the vertical stress of coal pillar 1 gradually decreased to 6.1 MPa, and the phenomenon of strata behavior appeared at the immediate roof of workface 2202; the vertical stress increased to 27.7 MPa. Moreover, the vertical stress of coal pillar 2 increased gradually, and the vertical stress values were 39.4 MPa, 39.9 MPa, 46.9 MPa, and 50.80 MPa, in sequence. In addition, the vertical stress values of the coal wall at pillar 3 first increased and then decreased, at 25.7 MPa, 33.1 MPa, 42.9 MPa, and 30.3 MPa, in sequence.
(4) When mining workface 2204, the vertical stress value of coal pillar 1 increased rapidly, from 7.01 MPa to 101 MPa, 102 MPa, and 103 MPa. Moreover, the vertical stress of coal pillar 2 and coal pillar 3 increased gradually, and the vertical stress at the immediate roof of the working face 2201 remained almost unchanged. In addition, slight strata behavior manifested at workfaces 2202 and 2203, and the vertical stress gradually increased.
By analyzing the stress evolution law of coal pillar 1, it can be seen that the dynamic range of the super-thick and weakly cemented overburden was wide during deep mining. For example, the stress concentration degree of the segment coal pillar at 600 m away from workface 2204 was further increased.
To further study the failure law of the super-thick and weakly cementation overburden in deep mining, we analyzed the evolution law of the pressure arch in the overburden during the mining process: It can be seen from Figure 33a-c that when the width of workface 2201 was 180 m, a beam fracture occurred in thin sandstone (immediate roof), and a double-pressure arch developed above the goaf, forming a stable arch shell structure to support the overburden load. When the width of workface 2201 was 240 m, three pressure arches developed above the goaf, and arch-shell failure occurred in thick sandstone. When the width of workface 2201 was 300 m, the three-pressure arch above the goaf continued to develop upwards, and the arch-shell failure range of thick sandstone was further expanded.
It can be seen from Figure 33d-g that when the width of workface 2202 was 120 m, beam fracture occurred in thin sandstone (immediate roof), and a double-pressure arch developed above the goaf, forming a stable arch shell structure to support the rock load above the immediate roof. At this time, four pressure arches developed above workface 2201, and double arch-shell failure occurred in the thick sandstone. When the width of workface 2202 was 180 m, there were two large pressure arches above the goaf of workfaces 2201 and 2202, producing a stable arch shell structure to support the overburden load. Moreover, beam failure continued to occur in the thin sandstone (immediate roof) above workface 2202, forming a masonry beam structure. When the width of workface 2202 reached 240 m, two large pressure arches above the goaf of workfaces 2201 and 2202 developed slightly upwards to support the overlying strata load sequentially, and archshell failure occurred in thick sandstone above workface 2202. When the width of workface 2202 reached 300 m, the range of the two large pressure arches above the goaf of workfaces 2201 and 2202 expanded slightly to support the overlying strata load, and the arch-shell failure range of the thick sandstone above workface 2202 is slightly increased.
It can be seen from Figure 33h,i that, when the width of workface 2203 reached 120 m, the scope of the double large pressure arch above the goaf of workfaces 2201, 2202, and 2203 slightly expanded, and the arch-shell failure height continued to develop upwards. When the width of workface 2203 reached 300 m, the Zhidan group sandstone became severely damaged, losing its bearing capacity, and the stable arch-shell structure disappeared. Multiple pressure arches and semi-pressure arches developed at the bottom of the Zhidan group sandstone, with arch feet concentrating in coal pillar 1, coal pillar 2, and the caving zone of workface 2202. Moreover, small stress arches developed in the fracture zone. Beam failure occurred in the thin sandstone above workface 2203, and arch-shell type failure occurred in the thick sandstone.
It can be seen from Figure 33j-l that, when the width of workface 2204 reached 120 m and 180 m, the development pattern of pressure arches in the overburden was similar to that described in Figure 33i. The thin sandstone above workface 2204 experienced beamtype damage, forming a beam structure, and a double-pressure arch developed above the goaf of workface 2204, forming a stable arch-shell structure to support the above-rock load. When the width of workface 2204 reached 300 m, the small pressure arches disappeared in the caving fracture zone of workfaces 2201 and 2202, the half arch-shell failure structure developed in the caving fracture zone of workface 2203, three pressure arches developed in the rock layer above workface 2204, and arch-shell failure occurred in thick sandstone. (k) (l) To deeply analyze the failure law of the super-thick and weakly cemented overburden in deep mining, the study simulated and analyzed the influence of the mining degree of the advancing direction and incline direction on the evolution law of the overburden pressure arch, without considering the influence of the section coal pillars on the overburden failure law.
It can be seen from Figure 34 that, with the continuous advancement of the working face, the pressure arch in the overlying strata developed forward. When the main key layer structure was broken, a stress concentration area formed in the broken place and the caving fracture zone below it, and the small pressure arches developed forwards with the advancement of the working face. At this time, the thick sandstone structure inside the small pressure arch was damaged through semi-arch shell type failure. To deeply analyze the failure law of the super-thick and weakly cemented overburden in deep mining, the study simulated and analyzed the influence of the mining degree of the advancing direction and incline direction on the evolution law of the overburden pressure arch, without considering the influence of the section coal pillars on the overburden failure law.
It can be seen from Figure 34 that, with the continuous advancement of the working face, the pressure arch in the overlying strata developed forward. When the main key layer structure was broken, a stress concentration area formed in the broken place and the caving fracture zone below it, and the small pressure arches developed forwards with the advancement of the working face. At this time, the thick sandstone structure inside the small pressure arch was damaged through semi-arch shell type failure.
By analyzing the evolution law of the overburden pressure arch in Figures 32 and 33, it can be seen that the failure law of the super-thick and weakly cemented overburden in deep mining presented a combined 'beam-arch shell' type failure, which transitioned into 'beam-semi-arch shell' type failure with the continuous expansion of the mining space.  By analyzing the evolution law of the overburden pressure arch in Figures 32 and 33, it can be seen that the failure law of the super-thick and weakly cemented overburden in deep mining presented a combined 'beam-arch shell' type failure, which transitioned into 'beam-semi-arch shell' type failure with the continuous expansion of the mining space.

Discussions
In order to verify the reliability of the numerical simulation results, we also carried out similar material model tests by simplifying the strata. In the test plan, the geometric similarity ratio and the similarity ratio of the unit weight of this physical experiment were 1:400 and 1:1.6, respectively. Before similar material model tests, test pieces had to be produced to carry out the compression test. Figure 35 illustrates the test equipment and experiment process. From the results of the compression test, the final ratio of the aggregate and the binder for similar material model was ensured, as shown in Table 9.  [13,20].

Discussions
In order to verify the reliability of the numerical simulation results, we also carried out similar material model tests by simplifying the strata. In the test plan, the geometric similarity ratio and the similarity ratio of the unit weight of this physical experiment were 1:400 and 1:1.6, respectively. Before similar material model tests, test pieces had to be produced to carry out the compression test. Figure 35 illustrates the test equipment and experiment process. From the results of the compression test, the final ratio of the aggregate and the binder for similar material model was ensured, as shown in Table 9. By analyzing the evolution law of the overburden pressure arch in Figures 32 and 33, it can be seen that the failure law of the super-thick and weakly cemented overburden in deep mining presented a combined 'beam-arch shell' type failure, which transitioned into 'beam-semi-arch shell' type failure with the continuous expansion of the mining space.

Discussions
In order to verify the reliability of the numerical simulation results, we also carried out similar material model tests by simplifying the strata. In the test plan, the geometric similarity ratio and the similarity ratio of the unit weight of this physical experiment were 1:400 and 1:1.6, respectively. Before similar material model tests, test pieces had to be produced to carry out the compression test. Figure 35 illustrates the test equipment and experiment process. From the results of the compression test, the final ratio of the aggregate and the binder for similar material model was ensured, as shown in Table 9.  [13,20].  The test results of similar material model were highly consistent with the numerical simulation results. In Figure 36, "beam-arch shell" failure occurred on the overburden in the similar material model, and the overburden failure boundary was arched in the UDEC model. This phenomenon can be explained by masonry beam and pressure arch theory. The test results of similar material model were highly consistent with the numerical simulation results. In Figure 36, "beam-arch shell" failure occurred on the overburden in the similar material model, and the overburden failure boundary was arched in the UDEC model. This phenomenon can be explained by masonry beam and pressure arch theory.
In this paper, the span of the pressure arch above the goaf gradually increased when the width of the working face 2202 increased from 60 m to 180 m, and the height of the pressure arch gradually developed upwards. The stress in the pressure arch reached the sandstone ultimate compression strength, "arch shell" failure occur when the width of the working face 2202 was 240 m. When the width of the working face 2202 was 300 m, the height of the fractured arch no longer retained to develop upward, the fractured arch continued to go forward, and the overburden was then destroyed in the form of a half arch. From the view of the mechanics, the failure mode in the semi-arch mechanical model was similar to that of the full arch pressure mechanics model. The bending moment, shearing force and axial force of the arbitrary section D of the half arch were the same as that of an arch. The "semi-arched shell" structure in Figure 36e was pressure shear damage, because the arch waist position bore the vertical downward load of the overburden and upward sloping forces of the lower strata and reached the shear ultimate strength of ultra-thick and weak cementation overburden.  [13,20].
In this paper, the span of the pressure arch above the goaf gradually increased when the width of the working face 2202 increased from 60 m to 180 m, and the height of the pressure arch gradually developed upwards. The stress in the pressure arch reached the sandstone ultimate compression strength, "arch shell" failure occur when the width of the working face 2202 was 240 m. When the width of the working face 2202 was 300 m, the height of the fractured arch no longer retained to develop upward, the fractured arch continued to go forward, and the overburden was then destroyed in the form of a half arch. From the view of the mechanics, the failure mode in the semi-arch mechanical model was similar to that of the full arch pressure mechanics model. The bending moment, shearing force and axial force of the arbitrary section D of the half arch were the same as that of an arch. The "semi-arched shell" structure in Figure 36e was pressure shear damage, because the arch waist position bore the vertical downward load of the overburden and upward sloping forces of the lower strata and reached the shear ultimate strength of ultra-thick and weak cementation overburden.

Conclusions
In this paper, the Yingpanhao Coal Mine was taken as the research object. Rock samples of the Zhidan group sandstone and Zhiluo group sandstone were obtained through drilling and coring. The physical properties and microstructural characteristics of weakly cemented rocks were further studied using an electronic microscope and other equipment. Based on the experimental results, a two-dimensional discrete element model of the study area was reasonably constructed to analyze the movement law and failure characteristics of the super-thick and weakly cemented overburden in deep mining, and the following conclusions were obtained: (1) The uniaxial compressive strength of the Cretaceous Zhidan group sandstone was found to be between 10 and 20 MPa, and the uniaxial compressive strength of the Jurassic sandstone was between 20 and 40 MPa. The mechanical parameters of different lithologic sandstones in the Cretaceous Zhidan group sandstone were generally smaller than those of the corresponding Jurassic Zhiluo group sandstones. The mechanical parameters of different lithologic sandstones in the Cretaceous Zhidan group decreased with increases in particle uniformity. The compressive strength and elastic modulus of different lithologic sandstones in the Jurassic Zhiluo group decreased with increases in particle uniformity, while their cohesion increases.
(2) The Zhidan group sandstone was characterized by a large grain size, large internal voids, and calcareous cementation. Moreover, it did not disintegrate in water. Further, the Anding group sandstone was characterized by a small grain size, along with clay mineral and hematite cementation. Moreover, internal fractures and pores developed, and the Anding group sandstone disintegrated easily. The Zhiluo group sandstone was relatively compact, cemented by clay minerals, and did not easily disintegrate in water. The Yan'an group sandstone was compact, cemented with iron, and did not disintegrate in water.
(3) In deep multi-workface mining with super-thick and weakly cemented overburdens, the dynamic behavior range of the overburden was wide, and the influence range extended to 600 m beyond the goaf boundary. The compressive deformation of coal pillars at different positions differed significantly. The compressive deformation of the coal pillar between the first mining face and the adjacent workface was the largest, and that of other coal pillars decreased sequentially. The surface subsidence developed significantly, and the surface subsidence coefficient exhibited a Boltzmann functional relationship with the width-depth ratio.
(4) With the continuous expansion of the goaf, the tangent of the main influence angle increased gradually; its change ranged is between 1.21 and 1.71, falling within the common value range. In addition, the correlation degree of the Boltzmann functional relationship between the W max of different buried depths and the height above the coal gradually weakened, and the functional relationship between the S W values at different buried depths and the heights above the coal gradually changed from a Boltzmann function to a parabolic function.
(5) The failure mode of the super-thick and weakly cemented overburden was 'beam-arch shell' failure, and the failure boundary was an arch fracture. The failure mode of thin weakly cemented sandstone was 'beam' failure, which finally formed the masonry beam structure. The failure mode of thick weakly cemented sandstone was 'arch-shell' failure, and with the continuous expansion of the mining scope, it gradually changed from complete 'arch-shell' failure to 'half-arch-shell' failure.

Data Availability Statement:
The data supporting reported results can be found by contacting the corresponding author (g_j_zhang@cumt.edu.cn).