Surface Subsidence Prediction Method for Coal Mines with Ultrathick and Hard Stratum

School of Mines, China University of Mining and Technology, Xuzhou 221116, China State Key Laboratory of Coal Resource and Safe Mining, China University of Mining and Technology, Xuzhou 221116, China State Key Laboratory of Groundwater Protection and Utilization in Coal Mining, China Energy Group, Beijing 100011, China Disaster Prevention Research Institute, Shaanxi Changwu Tingnan Coal Co., Ltd, Xianyang 713602, China

e profile function and influence function methods are now widely employed in all methods due to their practicality and ease of use [23].For conditions in which the topsoil is thick and the overburden strata include no ultrathick and hard stratum (UTHS), the aforementioned methods have relatively high precision for predicting surface subsidence.However, for overburden strata with UTHS, the prediction results often differ greatly from measured results.For example, when mining under the ultrathick igneous rock of the Haizi Coal Mine, the predicted maximum surface subsidence value obtained using influence function methods was significantly larger than the measured value (1270 mm versus 457 mm, respectively) [24].
is occurs because most of the above prediction methods consider the overburden strata as a homogenous medium [25].When the composition of the overburden strata is fairly uniform, this simplification is reasonable [25].However, surface subsidence is the result of a gradual development of overburden strata from bottom to top after mining, and different compositions of overburden strata have strong impacts on surface subsidence.Although the existing prediction methods take into account the influence of lithology of overburden strata on surface subsidence when choosing the prediction parameters, consideration of overburden remains insufficient [26].For instance, in China, one of the influence function methods, the probability integral method [14], has been used widely.To study surface subsidence, it regards overburden strata as a granular random medium and classifies overburden strata into three types, hard, medium, and soft, when choosing prediction parameters [14].When the overburden strata includes UTHS, considering the overburden strata as a homogenous medium is obviously an oversimplification of the overburden condition, which ignores the control function of UTHS on surface subsidence.It does not relate surface subsidence with the movement of UTHS and does not consider the strong influence of the characteristics of movement of UTHS on surface subsidence.Many case studies, such as the movement of thick and hard sandstone in the Datong mine area, the thick glutenite in the Huafeng Coal Mine in Xinwen, and the thick and hard sandstone in the Lutugin Coal Mine in Ukraine, have shown that the control function of UTHS on surface subsidence leads to overall movement of all of the overlying strata up to the ground surface [26].Considering the previously mentioned shortcomings, the existing prediction methods of subsidence are no longer applicable under the condition of UTHS.us, a surface subsidence prediction method for coal mining with the presence of UTHS should be established.
In China, an UTHS is generally defined as a rock stratum that is about hundreds of meters thick [24], with a uniaxial compressive strength of more than 40 MPa [27].
e overburden strata condition of UTHS exists widely in China and other countries.Some examples include the ultrathick conglomerate strata of the Huafeng Coal Mine in China [28], the ultrathick igneous rock of the Haizi Coal Mine in China [24], the ultrathick sandstones of the Tingnan Coal Mine in China, and the ultrathick dolerite sill of the Churcha West Colliery in India [29].A surface subsidence prediction method for coal mining with the presence of UTHS must consider the control function of UTHS on surface subsidence.To improve the precision of surface subsidence prediction for coal mining with the presence of UTHS, we must know the movement characteristics of UTHS accurately.It should be pointed out that when the thickness of an alluvium layer is more than 100 m, the thick alluvium layer will also have a strong impact on the surface subsidence [30].erefore, this paper only considers surface subsidence prediction when the overburden strata have UTHS without a thick alluvium layer.
e existing surface subsidence prediction methods ignore the strong impact and control function of UTHS on surface subsidence.ey also lack correct understanding of the movement characteristics of UTHS and do not relate surface subsidence with UTHS movement.us, they are not applicable for surface subsidence prediction for coal mining with the presence of UTHS.We conducted actual measurements of surface subsidence attributable to coal mining with the presence of UTHS and established UTHS internal rock movement monitoring drill holes to monitor the UTHS subsidence in the Tingnan Coal Mine.We discovered the relationship between surface subsidence and UTHS movement.We studied the UTHS movement characteristics using Winkler's theory of beams on an elastic foundation, developed a surface subsidence prediction method for coal mining with the presence of UTHS, and evaluated the influences of UTHS thickness, strength, and layer position on surface subsidence.We took the Tingnan Coal Mine as an example and calculated the subsidence curves of strike and dip main sections with different mining dimensions.en, we predicted the surface subsidence after the mining of the 204, 204 + 205, 204 + 205 + 206, and 204 + 205 + 206 + 207 working surfaces and verified the accuracy of the predicted results using the actual measurement results of surface subsidence of the 204, 204 + 205, and 204 + 205 + 206 working surfaces.

Actual Measurement of Surface and UTHS Subsidence for Coal Mines with UTHS
Ignoring the influence of UTHS often leads to distorted prediction results of surface subsidence, as for the previously mentioned Haizi Coal Mine.To study the influence of UTHS on surface subsidence, we conducted longterm surface subsidence measurements at the Tingnan Coal Mine and actual subsidence measurements of the UTHS.

Surface Subsidence Measurement.
e Tingnan Coal Mine is located in the Changwu County, Xianyang, Shaanxi Province, and it has an average mining height of 7.0 m and a coal seam dip angle of 0 °in the second panel.
ere is an ultrathick (222.14 m) and hard sandstone stratum, which is a moderately water-rich aquifer with an average uniaxial compressive strength of 42.2 MPa. e recovery periods of working faces 204, 205, and 206 were November 2011-November 2012, June 2013-November 2014, and January 2015-April 2016, respectively.Currently, the panel is recovering working face 207. Figure 1 shows the layout of the working area in the second panel 2 Advances in Civil Engineering and the positions of surface subsidence measurement points, and the integrated stratigraphic column in the second panel can be seen in Figure 2.
Considering the potential influences of the surface landform and streams, a surface subsidence observation line was deployed along the dip of the second panel to observe surface subsidence, as shown in Figure 1.After the mining of the 204, 204 + 205, and 204 + 205 + 206 working surfaces, the measured maximum surface subsidence was 0.248, 1.194, and 1.925 m, respectively.Because of the influence of surface reservoir construction, most of the measurement points could no longer be used to monitor surface subsidence after the recovery of working face 207. Figure 3 shows the observed surface subsidence curves in the dip direction (measurement point no.34-A1, shown in Figure 1) after the mining of working faces 204, 205, and 206, respectively.Based on the mining experience [14], when the dimensions of the gob in both the strike and dip directions reach 1.2-1.4times the mining depth of the seam, the surface reaches sufficient mining, which means that when the second panel's mining dimensions in both the strike and dip directions reach 657-766 m, the surface should reach sufficient mining.After the mining of working faces 204, 205, and 206, the strike and dip mining dimensions of the second panel gob were 2200 m and 660 m, respectively.erefore, the surface of the second panel should have approached sufficient mining.Based on the empirical subsidence coefficient of 0.65 of the Tingnan Coal Mine under sufficient surface mining, after the mining of working faces 204, 205, and 206, the maximum empirical subsidence value should have approached 4.55 m, but the measured maximum surface subsidence value was only 1.925 m, much smaller than 4.55 m, indicating that the surface was far from reaching sufficient mining.Under the UTHS condition, the required dimensions of the gob in the strike and dip directions are much greater than 1.2-1.4times the mining depth in the seam.In addition, the predicted surface subsidence curves of the observation line were also plotted using the influence function method [14] with empirical parameters collated from the Tingnan Coal Mine.
e maximum values of the predicted surface subsidence are 2.4, 3.71, and 4.07 m, respectively (Figure 3).e influence function method calculations produced results significantly larger than those obtained by measurements.Based on the above results, we could be certain that UTHS has a strong impact on surface subsidence.

Measurement of UTHS Subsidence.
As the UTHS is inside the overburden, it is impossible to measure its exact amount of deflection.After the mining of working face 206, to monitor the relative movement between the UTHS and the surface during the mining of working face 207, we drilled a strata displacement-monitoring borehole Y1-1 inside the overburden strata of working face 206.e position of the borehole is shown in Figure 1, and no separation spaces were detected during drill work.e measuring point was placed 130 m below the UTHS top boundary and connected to a surface displacement recorder using a wire rope to monitor the relative movement between the measuring point and the surface, as shown in Figure 2 and Figures 4(a) and 4(b).As of September 3, 2018, working face 207 had been mined 916.8 m past the monitoring borehole.e surface of Y1-1 had subsided another 1.352 m since the mining of working face 207 and was still subsiding at 0.41 mm per day without showing signs of becoming stable.However, the measuring point inside the UTHS on the other side subsided only 0.077 m and had already stabilized.
e monitoring results are shown in Figure 4(c).Evidently, the UTHS and the overburden strata up to the surface have approximately made an overall movement, and the UTHS subsidence can approximately represent the surface subsidence.Advances in Civil Engineering e measured surface and UTHS subsidence at the Tingnan Coal Mine show that UTHS has a strong impact on the surface subsidence.us, the lack of a correct understanding of the surface subsidence rules caused by UTHS movement in different mining dimensions could lead to a distorted prediction result of surface subsidence.As a possible solution, we could take the subsidence of UTHS as an approximate surface subsidence and then predict the surface subsidence by calculating the UTHS subsidence.In this paper, the UTHS subsidence represents the surface subsidence.

Mechanical Model for Subsidence Prediction for Coal
Mining with UTHS.After a coal seam is mined, a breaking movement and a flexural subsidence are expected in the overburden strata of the mining site, which is composed of the caved zone, fractured zone, and continuous deformation zone [12,26].e UTHS, which is generally in the continuous deformation zone, is subject to a curve subsidence rather than a breaking movement [12,26].To simplify the process, we consider the goaf to be filled up with bulked rocks from the caved and fractured zones [31]; the rocks underlying the UTHS are considered as an elastic foundation complying with the Winkler elastic foundation hypothesis [31][32][33][34][35][36]. e UTHS is then thought of as an elastic foundation beam on which both the dead weight of the UTHS itself and its load act as uniform loads [31][32][33][34][35], and it has a dip angle of 0 °.In Chinese coal mines, the most commonly used mining method is longwall mining in which the stoping space created by the stoping of one or more working faces is rectangular in shape, and the mining area contains two main sections: the dip main section I and the strike main section II.For each of the two main sections, based on the aforementioned assumptions, we established a mechanical model for subsidence prediction for the coal mining with the presence of UTHS, as shown in Figure 5.
e symbol q (MPa) represents the weight of the UTHS and the strata on it.Furthermore, H (m) is the thickness of the UTHS, E (GPa) is the elastic modulus of the UTHS, I (m 4 ) is the inertia moment of the UTHS, L 1 (m) is the mining dimension of main section I, L 2 (m) is the mining dimension of main section II, k (N/m 3 ) is the foundation coefficient of the UTHS at the coal side, k c (N/m 3 ) is the foundation coefficient of the UTHS at the goaf side, p s (MPa) is the foundation reaction of the UTHS at the coal side, and p c (MPa) is the foundation reaction at the goaf side.
If the calculated main section in Figure 5(b) is main section I, we assume that main section II has reached sufficient mining.As the UTHS in Figure 5(b) is symmetric about x � L 1 /2, only half of the UTHS was used for our analysis, as shown in Figure 5(c).e differential equations of the deflection curves can be established for the UTHS as follows: In the equation, w Is (x) (m) and w Ic (x) (m) are the UTHS deflection curve equations of main section I at the coal side and the goaf side, respectively, I (m 4 ) is the moment of inertia of the UTHS, where I � bH 3 /12, and b (m) is the width of the beam, considered to be equal to unity, i.e., b � 1.
en, the general solution for equation ( 1) is ( In the equation, a I1 � ����� k/4EI 4 √ .As x ⟶ −∞ and w Is (x) ⟶ q/k, these conditions can only be satisfied when C Ι3 and C Ι4 are equal to zero.Additionally, because q/k in equation ( 3) is caused by the premining compression of the seam, the UTHS deflection Advances in Civil Engineering curve equation at the coal side induced by coal mining, w Iw (x), is equal to en, the general solution for equation ( 2) is  .e UTHS deflection curve equation at the goaf side, w Ic (x), is shown in equation (5).From the relationship among the rotation angle, bending moment, and shear force of the UTHS with its deflection curve equation in equation (6) [32,34,36], by substituting the deflection curve equations (4) and ( 5) into (6), we can derive the rotation angle θ Ιw (x), bending moment M Ιw (x), and shear force Q Ιw (x) of the UTHS at the coal side as well as the rotation angle θ Ιc (x), bending moment M Ιc (x), and shear force Q Ιc (x) of the UTHS at the goaf side.
e deflection curve equations ( 4) and ( 5) contain six unknown parameters: C Ι1 , C Ι2 , C Ι5 , C Ι6 , C Ι7 , and C Ι8 .e UTHS is continuous at x � 0, where the subsidence, rotation angle, bending moment, and shear force of the UTHS are equal.e model is symmetric about x � L 1 /2, where the shear force and rotation angle of the UTHS are equal to zero [34,36].Using these continuity and boundary conditions for the UTHS in equation ( 7), we can derive the six unknown parameters.By substituting those values into equations ( 4) and ( 5), we obtain the UTHS deflection curve equations at the coal side and the goaf side.
When w Ic (L 0 /2) � W max , i.e., when the mining width L 1 is equal to L 0 , the UTHS will be under sufficient mining and the maximum subsidence [14,27] of the UTHS reaches the maximum W max under this geological mining condition.
e symbol δ represents the subsidence coefficient [14,27], which can be found in the literatures [14,27], whereas M is the mining thickness of the seam.
When L 1 ≤ L 0 , the unknown values of C Ι1 , C Ι2 , C Ι5 , C Ι6 , C Ι7 , and C Ι8 in w Iw (x) and w Ic (x) are determined by equation (7); as L 1 increases, these parameters will also change, i.e., the maximum subsidence of the UTHS will change with the increase of L 1 .en, when main section ΙΙ reaches sufficient mining and main section I does not, the UTHS deflection curve equation of main section I can be written as follows: Advances in Civil Engineering When L 1 > L 0 , to determine the six unknown parameters C Ι1 , C Ι2 , C Ι5 , C Ι6 , C Ι7 , and C Ι8 in w Iw (x) and w Ic (x), the continuity conditions in equation (7) remain the same at x � 0, but the boundary conditions change to θ c (L 0 /2) � 0 and Q c (L 0 /2) � 0. ese unknown values will no longer vary with the increase of L 1 .us, the maximum subsidence of the UTHS will not vary with the increase of L 1 .en, when both main section I and main section II reach sufficient mining, the UTHS deflection curve equation of main section I is If the calculated main section in Figure 5(b) is main section II, we assume that main section I has reached sufficient mining.In the same way, if we replace the parameters of main section I (x, L 1 , w Is (x), w Ic (x), , and Q Ιc (x)) with the parameters of main section II (y, L 2 , w IIs (y), w IIc (y), w IIw (y), a ΙΙ1 , a ΙΙ2 , , and Q ΙΙc (y)) in the calculation, we can obtain the deflection curve equation of section II with UTHS when main section I reaches sufficient mining and main section II does not: w ΙΙw (y) � C ΙΙ1 e a ΙΙ1 y cos a ΙΙ1 y  + C ΙΙ2 e a ΙΙ1 y sin a ΙΙ1 y , y ≤ 0, When both main section I and main section II reach sufficient mining, the UTHS deflection curve equation of main section II is where w IIw (y) (m) is the UTHS deflection curve equation of main section II at the coal side, and w IIc (y) (m) is the UTHS deflection equation of main section II at the gob side.
To solve the above equations, it is necessary to determine the parameters E, H, M, L 1 , L 2 , q, k, and k c .e first five of these (E, H, M, L 1 , and L 2 ) can be determined in relation to 8 Advances in Civil Engineering the borehole columns, the physical-mechanical parameter test results, and the actual mining conditions.e uniform load on the UTHS, q, can be calculated as [26] q �  z j�1 c j h j .(13) Here, h j is the thickness of the UTHS and any of the overburden load strata, c j is the bulk density of the UTHS and any of the overburden load strata, and z is the total number of the UTHS and the overburden load strata.
e foundation coefficient k of the UTHS at the coal side is calculated as [21,35,36] e variable H i represents the thickness of any of the foundation strata, E i is the elastic modulus of any of the foundation strata, and n is the total number of the foundation strata, often taken as the total number of strata between the UTHS and the thick and hard rock stratum in the coal floor.
Because the model assumes that the fractured rock in the gob region satisfies Winkler's elastic foundation hypothesis, we can calculate the ground coefficient k c of UTHS at the gob side if we know the subsidence value [36] of UTHS and the corresponding load bearing of the foundation [36] in the gob region.When the surface reaches sufficient mining, UTHS subsidence reaches the maximum value W max under this mining condition [14], and the corresponding load bearing of the foundation at the gob side is q [26].en, the ground coefficient in the gob region k c can be written as

Subsidence Prediction Equation of Main Sections.
According to the actual mining dimensions of the two main sections, the surface has four mining extents [14]: main section I does not reach sufficient mining, but main section II does; main section I reaches sufficient mining, but main section II does not; both main section I and main section II reach sufficient mining; and neither of the main sections reaches sufficient mining.
(1) When main section I does not reach sufficient mining but main section II does, equation ( 9) is the UTHS deflection curve equation of main section I. e UTHS deflection curve equation of main section II can be calculated through equation ( 12), but the equation needs to be multiplied by a coefficient, λ Ι , which is named the mining extent coefficient [14] of main section I and represents the times of reduction of the main section II subsidence value when main section I is at different levels of insufficient mining rather than sufficient mining.
w Ic (L 1 /2) is the subsidence value of the point with the coordinate of L 1 /2 on main section I when main section II is under sufficient mining, and it is calculated through equation ( 9).(2) When main section II does not reach sufficient mining but main section I does, equation ( 11) is the UTHS deflection curve equation of main section II.e UTHS deflection curve equation of main section I can be calculated through equation (10), but the equation needs to be multiplied by a coefficient, λ ΙΙ , which is named the mining extent coefficient [14] of main section II and represents the times of reduction of the main section I subsidence value when main section II is at different levels of insufficient mining rather than sufficient mining.
w IIc (L 2 /2) is the subsidence value of the point with the coordinate of L 2 /2 on main section II when main section I is under sufficient mining, and it is calculated through equation ( 11).(3) When both main section I and main section II reach sufficient mining, the UTHS deflection curve equations of main section I and main section II are equations ( 10) and ( 12), respectively.(4) When neither of them reaches sufficient mining, the UTHS deflection curve equation of main section I can still be calculated through equation ( 9), but the equation needs to be multiplied by a coefficient, λ ΙΙ , which is shown in equation ( 17). e subsidence of main section II can still be calculated through equation (11), but this equation also needs to be multiplied by a coefficient, λ Ι , which is shown in equation ( 16).

Subsidence Prediction Equation on Any Arbitrary Cross
Section Parallel to the Main Sections.Section 3.2 solves the subsidence prediction equations of the two main sections under four mining extents in the mining area.e subsidence distributions parallel to the cross sections of main section I and main section II are similar to those of the corresponding main sections [14], but the values are smaller than the predicted subsidence values of the main sections obtained in Section 3.2.e magnitudes are ζ ΙΙ times the deflection curve equation of main section I and ζ Ι times the deflection curve equation of main section II in Section 3.2, respectively.erefore, if the cross section parallel to main section I is y � y 0 and the cross section parallel to main section II is x � x 0 , we have the following expression:

Advances in Civil Engineering
In the equation, w IIw (y 0 ) and w IIc (y 0 ) are the subsidence values of point y 0 on main section II when main section I has reached sufficient mining.When main section II has not reached sufficient mining, w IIw (y 0 ) and w IIc (y 0 ) can be calculated through equation (11); when main section II has reached sufficient mining, w IIw (y 0 ) and w IIc (y 0 ) can be calculated through equation (12).
In the equation, w Iw (x 0 ) and w Ic (x 0 ) are the subsidence values of point x 0 on main section I when main section II has reached sufficient mining.When main section I has not reached sufficient mining, w Iw (x 0 ) and w Ic (x 0 ) can be calculated through equation ( 9); when main section I has reached sufficient mining, w Iw (x 0 ) and w Ic (x 0 ) can be calculated through equation (10).
After solving the subsidence prediction equation on any arbitrary cross section parallel to the main sections, we can predict the subsidence of the whole mining area.

Example of Calculation and Verification
Using the Proposed Surface Subsidence Prediction Method

Example of Surface Subsidence Calculation for the Tingnan
Coal Mine.e model calculation parameters, as shown in Table 1, were chosen according to the geological mining conditions of the second panel in Tingnan Coal Mine, the rock mechanical parameter test results, and the model parameter determination method described in Section 3.1.
To calculate the surface subsidence under different mining dimensions L 1 of main section I, we considered that with the increase of mining number on the working face of the second panel, the mining dimension L 1 of main section I increases.e average propelling distance of main section II is approximately 2200 m.We took L 2 � 2200 m and calculated according to the calculation method in Section 3.According to the calculation method in Section 3 and the symmetric nature of the mechanical model, the subsidence curves of main section I when taking different mining dimensions L 1 are shown in Figure 6.
To calculate the surface subsidence of main section I when the mining dimension L 1 of main section I is 200 m, 430 m, 660 m, 890 m, 1120 m, 1350 m, and 1580 m, with increase of L 1 , main section I changes from insufficient mining to sufficient mining.According to the calculation methods in Section 3 and the symmetric nature of the mechanical model, the subsidence curves of main section II when taking different mining dimensions L 1 are shown in Figure 7.
From Figures 6 and 7, when L 1 � L 2 � L 0 � 1350 m, both main section I and main section II reach sufficient mining, and the maximum subsidence value reaches the maximum value W max of 4.55 m.When L 1 ≤ L 0 , with increase of L 1 , the maximum subsidence value increases gradually.When L 1 > L 0 , with increase of L 1 , the maximum subsidence value does not change further.
us, under the UTHS mining condition in the Tingnan Coal Mine, the required dimension (1350 m) in the strike and dip directions in the gob region when the surface reaches sufficient mining is much larger than the empirical value (1.2-1.4 times the seam mining depth of 657-766 m), which verifies the actual measurement results in Section 2.1.
e calculated subsidence curves of main section I and main section II under different mining dimensions L 1 of main section I can be used in the calculation method in Section 3.3 to solve the subsidence prediction equation for any arbitrary cross section parallel to the main sections, and then the subsidence of the whole mining area can be predicted.Figure 8 shows the subsidence contour maps after the mining of working faces 204, 205, 206, and 207 calculated by the proposed prediction method, respectively.

Verification of Actual Measurement of Surface Subsidence in the Tingnan Coal Mine.
e subsidence at each measurement point in Figure 1 after the recovery of working faces 204, 205, and 206, respectively, can be obtained by the proposed prediction method, and comparing these values with the measured subsidence values can verify the correctness of the proposed prediction method.e predicted subsidence at each point can be obtained from Figure 8, and the comparison between the predicted and measured results is shown in Figure 9.As indicated by the figure, the predicted results are consistent with the measured results, which proves the correctness of the proposed prediction method.

Discussion
We investigated the influences of UTHS thickness, strength, and layer position changes on the surface subsidence of main section I when main section II of the Tingnan Coal Mine reaches sufficient mining with the mining dimension of main section I of L 1 � 660 m using the proposed prediction method.In the calculation parameters in Table 1, by prescribing UTHS thickness H as 100 m, 160 m, 222 m, and 300 m, respectively, or strength E as 11.5 GPa, 19.5 GPa, 27.5 GPa, and 35.5 GPa, respectively, we calculated surface subsidence with different H or E to show the influences of UTHS thickness and strength on surface subsidence.e

Conclusions
UTHS has a strong impact on the distribution of surface subsidence, but current surface subsidence prediction methods lack consideration of the influence of UTHS on surface subsidence.
ese methods ignore the control function on surface subsidence of UTHS.Furthermore, they do not appropriately consider the movement characteristics of UTHS and do not relate surface subsidence attributable to coal mining with the presence of UTHS.erefore, our paper describes a surface subsidence prediction method for coal mining with the presence of UTHS.e results obtained are summarized below.
In our study, actual measurements of surface subsidence for coal mining with UTHS were conducted, and UTHS internal rock movement monitoring drill holes were established to monitor the subsidence of UTHS in the 0.38 1.17 x (m) Left mining boundary     Advances in Civil Engineering Tingnan Coal Mine.It was found that there was a relationship between surface subsidence and UTHS movement.e results showed that under the UTHS mining condition, the required dimension in the strike and dip directions in the gob region is much larger than the empirical value when the surface reaches sufficient mining and that the actual measured value of surface subsidence is much smaller than the empirical value.e subsidence of UTHS is approximately equal to the surface subsidence.e movement of UTHS has a strong impact on surface subsidence and has a controlling function for it.It was proposed that surface subsidence could be approximately predicted by calculating the UTHS subsidence.
In this work, the UTHS movement characteristics were studied using Winkler's theory of beams on an elastic foundation, the subsidence prediction equation of the main sections in the strike and dip directions was obtained under different mining dimensions, and the subsidence prediction equation of any arbitrary cross section parallel to the two main sections was established.en, the surface subsidence prediction method for coal mining with the presence of UTHS was developed.
e Tingnan Coal Mine was taken as an example, and the subsidence curves of the strike and dip main sections were calculated using different mining dimensions.e surface reached sufficient mining when the mining region reached  Advances in Civil Engineering 13 1350 m in both strike and dip directions, which is much larger than the empirical range of 657-766 m (1.2-1.4 times the mining depth of the seam).Subsequently, the surface subsidence after the mining of working faces 204, 205, 206, and 207, respectively, was predicted, and the prediction method was verified by comparing the results with the measured surface subsidence results of working faces 204, 205, and 206.e influences of UTHS thickness, strength, and layer position on the surface subsidence were discussed when main section II of the Tingnan Coal Mine reaches sufficient mining and main section I has a mining dimension of L 1 � 660 m. ere are significant differences in the surface subsidence curves under different UTHS thicknesses, strengths, and layer positions, and variations in overburden strata compositions have a strong impact on surface subsidence.e greater the UTHS thickness or strength, the greater the range of surface subsidence and the lower the maximum surface subsidence value.When the distance between UTHS and seam increases, the maximum surface subsidence value decreases significantly and the range of surface subsidence increases slightly.

Figure 1 :
Figure 1: Map of the layout of the working area in the second panel and the positions of surface subsidence measurement points.

Figure 2 :
Figure 2: e integrated stratigraphic column in the second panel.

Figure 3 :
Figure 3: Comparison between the measured subsidence value and the predicted subsidence value by the probability integral method at the surface subsidence measurement point (measurement point no.34-A1).

Figure 4 :
Figure 4: Subsidence monitoring of surface and UTHS after recovery of the 207 working face.(a) Schematic diagram of UTHS subsidence monitoring drill hole.(b) UTHS subsidence monitoring equipment.(c) UTHS and surface subsidence curves.

Figure 5 :
Figure 5: Mechanical model for subsidence prediction for coal mining with UTHS.(a) Simplified floor plan of working faces of second panel in Tingnan Coal Mine.(b) Cross section plot of overburden strata with UTHS in the main section I (II).(c) Subsidence prediction mechanical model for coal mining with presence of UTHS in the main section I (II).
) and 10(b), with increase of UTHS thickness or strength, the range of surface subsidence increases significantly and the maximum surface subsidence value decreases significantly.Figure10(c)shows the surface subsidence prediction results after moving the UTHS by +131.55 m, +67.2 m, 0, and −51.55 m, respectively.It represents the influences of different layer positions on the surface subsidence.From Figure10(c), we can determine that when the distance between UTHS and the seam increases, the maximum surface subsidence value decreases significantly and the range of surface subsidence increases slightly.Comparing Figures10(a)-10(c), it can be seen that there are significant differences in surface subsidence curves under different UTHS thicknesses, strengths, and layer positions.us, different overburden strata compositions have a strong impact on the distribution of surface subsidence.

Figure 6 :
Figure 6: Subsidence curves of main section I when taking different mining dimensions L 1 .

Figure 7 :
Figure 7: Subsidence curves of main section II when taking different mining dimensions of main section I.

Figure 8 :
Figure 8: Subsidence contour maps after the mining of working faces 204, 205, 206, and 207, respectively, in the second panel.(a) After the mining of working face 204.(b) After the mining of working faces 204 and 205.(c) After the mining of working faces 204, 205, and 206.(d) After the mining of working faces 204, 205, 206, and 207.

Figure 9 :Figure 10 :
Figure 9: Comparison between predicted results obtained by the proposed method and measured results at the measurement points of surface subsidence.

Table 1 :
e model calculation parameters in the Tingnan Coal Mine.