Reader TTS

ABSTRACT In this study, the change law of the surface subsidence coefficient under the condition of thick coal seam layered mining was investigated. The study is based on the measured

subsidence data of the 1210-working face of the Mengba mine surface mobile observation station after the first- and second-layer mining. UDEC numerical simulation software was used to

simulate the variation of surface subsidence coefficient after the first, second, third, fourth, fifth-, and sixth-layer mining when the thickness of slicing mining is 5 m. The maximum

relative error between the simulated result and the measured result of the subsidence coefficient _q_ is 2.7%, which further verifies the correctness of the established model. Moreover, the

simulation results show that with the increase of the cumulative mining thickness, the subsidence coefficient _q_ of the surface presents a segmented characteristic. When the cumulative

mining thickness does not reach 25 m, the subsidence coefficient of the surface gradually increases with the increase of the mining thickness. On the other hand, when the cumulative mining

thickness reaches 25 m, the subsidence coefficient of the surface will tend to a constant value and no longer change with the increase of the mining thickness. Finally, the calculation

formula between the surface subsidence coefficient and the cumulative mining thickness of layered mining under the condition of medium hard roof is fitted, which provides a parameter basis

for coal seam mining with similar geological conditions. SIMILAR CONTENT BEING VIEWED BY OTHERS ANALYSIS OF THE ROOF DAMAGE RANGE IN CLOSE-PROXIMITY GENTLY INCLINED COAL SEAMS MINING AND THE

FEASIBILITY OF UPWARD MINING Article Open access 13 February 2025 INVESTIGATION ON THE CHARACTERISTICS OF SURFACE MOVEMENT AND DEFORMATION IN MINING UNDER THICK UNCONSOLIDATED LAYER IN

JINING AREA, CHINA Article Open access 19 May 2025 SUBSIDENCE PREDICTION OF OVERBURDEN STRATA AND GROUND SURFACE IN SHALLOW COAL SEAM MINING Article Open access 23 September 2021

INTRODUCTION The stratum and ground surface will move and deform due to the destruction during the primary coal mining. The underground strata movement and surface movement deformation will

occur once more whenever mining its upper or lower coal seam again. This process is known as repeated mining1. Nowadays, given the large average recoverable thickness of thick coal seams, a

fully mechanized top-coal caving mining method has been mainly adopted. Meanwhile, the large mining intensity of the fully mechanized caving will inevitably aggravate the damage to the rock

seam above. The ground subsidence coefficient (_q)_, the main influence angle tangent (tanβ), and other angular parameters are different from those of the primary mining. At present, the

research on the law of repeated mining of thick coal seams mainly focuses on the mining depth ranging from 3.5 to 15 m2, while few studies investigated mining depth of more than 15 m3. The

slice mining of thick coal seams is popular, and the space above the mining area is larger than that of single-layer mining, and the movement range of the overlying rock layer is also

increased. Therefore, the study on surface subsidence and defamation law is of high significance. Yu et al.4 investigated the relationship between the development of the water suture zone in

the overlying strata and different mining thicknesses, which revealed the relationship between the development height of the water suture zone and the mining thickness of the coal seam,

with certain guiding significance. Li et al.5 applied the numerical simulation software UDEC to simulate and study the development height of “three zones” after the mining of thick coal

seams, and the results showed that the sinking value of the overlying strata roof gradually decreased with the increase of the distance between the roof and the coal seam to a constant

value. Wang et al.6 obtained the surface subsidence law and related rock movement parameters after mining of thick coal seams under comprehensive discharge conditions through the analysis of

the actual measurement study of the surface movement observatory in the Baodian mine, which provided a certain basis for the mining of coal seams under similar geological mining conditions

in China. The subsidence law of the surface caused by mining when the accumulative mining thickness is greater than 15 m is yet to be studied. Mastering the movement and deformation law of

the surface of thick coal seam under the condition of stratified mining can provide the basis for the mining of coal seams in other mining areas under similar conditions in China. In this

paper, taking the Barapukuli coal mine (referred to as Mengba mine) in Bangladesh as an example, a model was constructed by combining the measured subsidence data of surface station to

simulate the law of the surface subsidence coefficient in the layered mining of thick coal seams, and based on the simulation results of the UDEC numerical simulation software. The

simulation data and the actual measurement data were analyzed by regression to ensure the reliability and reasonableness of the simulation parameters, and the formula between the surface

subsidence coefficient and the accumulative mining thickness under the layered mining conditions was fitted, which is expected to provide a parameter basis for the mining of similar coal

mines in China. OVERVIEW OF THE MINING AREA The Mengba Mine is a modern large mine located in the northwest of Bangladesh with a design capacity of 1 million tons per year and a service life

of 60 years. The mine covers 4.9 km in length from north to south and 0.3–1.9 km in length from east to west. The topography of the mine area is mainly made up of a plain formed by the

alluvial deposits of the Ganges and Jamuna rivers, with a flat surface and a small slope, high in the north and low in the south. From the generalized geological section of the mine area, it

can be observed that the mine area contains 7 layers from top to bottom, namely I, II, III, IV, v, VI, and VII coal seams, with a total thickness of about 75 m and a coal-bearing factor of

20%. The recoverable thickness of VI coal reaches 36 m, which is distributed in the whole area and is the main coal seam for mining. The direct roof of the coal seam is a medium-hard

coarse-grained feldspar sandstone with a thickness of 95 m. Due to the large thickness of the coal seam, the method of slicing mining is adopted, and the mining thickness of a single layer

is about 5 m. The first stratum was mined from 2005 to 2012, and mining of the second stratum of VI coal is currently underway. THE ESTABLISHMENT OF THE SURFACE OBSERVATION STATION The

1210-mining face has an east–west full trend line of 1,200 m length and a north–south full incidence line of 500 m length. The former has 50 points and the latter 20 points. The thickness of

single-layer mining is 5 m. MATHEMATICAL MODELING GEOMETRIC MODELING This simulation of the movement changes of the surface after coal seam mining utilizes the numerical simulation software

UDEC4.0 to model the geological mining conditions of the 1210-mining face in the Mengba mine area with the thickness of a single seam mining being 5 m. The actual strata with thin thickness

and similar rock properties were combined together7,8,9,10,11 and treated as one layer to facilitate modeling during the numerical simulation. The overburden above the 1210-working face was

divided into 17 parts after the combination. Where part 16 represents the Lower Dupi Tila (LDT) and part 17 is the Upper Dupi Tila (UDT), as shown in Fig. 1. SIMULATION OF THE MINING

RECOVERY This article mainly simulates the subsidence of the surface after layered mining of thick coal seams, and the numerical simulation adopts a one-time full-mining method. According to

the mining plan of the 1210-working face, the first stratum was mined first, and the second, third, fourth, fifth, and sixth stratum will be mined successively after the various movement

and deformation of the ground surface being balanced. To obtain the subsidence coefficient of the ground surface after fully mining the coal seam, the mining length along the trend direction

is 600 m (1.5 times the average mining thickness), and a boundary of 400 m is taken on each side, so the total length of the model built is 1,400 m. PHYSICAL PARAMETERS OF EACH ROCK

FORMATION The Mohr–Coulomb yield criterion was used to determine the failure of the rock mass throughout the simulation, and plastic flow was not considered. The physical parameters of each

rock layer of the numerical simulation model were determined according to the generalized geological section of the Mengba mine. The physical parameters of the main rock formations used in

the simulation are shown in Table 1. NUMERICAL SIMULATION ANALYSIS ANALYSIS OF OVERBURDEN DISPLACEMENT AND SUBSIDENCE LAW Before analyzing the changes in surface subsidence coefficient and

subsidence activation coefficient with the accumulative mining thickness, the correctness of the established model should be verified first, and the specific comparison results are shown in

Fig. 2, Tables 2, and 3. The maximum sinking value of the ground was 3,550 mm after the recovery of the first stratum of the 1210-working face, and the numerical simulation results showed

that the maximum sinking value of the ground was 3550 mm, corresponding to the mining thickness of 5 m. The sinking coefficient q = 0.71 was then obtained, and the sinking coefficient q = 

0.70 was obtained by fitting the measured sinking data from the surface station using the referencing software. The maximum height of the hydraulic fracture zone after mining is 93.43 m, and

the maximum development height of the water-conducting fractured zone after mining the coal seam was 96 m. The maximum relative error between the simulation results and the measured results

q was 2.5%, and the maximum relative error of the water-conducting fractured zone was 2.7%. The maximum measured surface subsidence after the second stratification was 7,800 mm, and the

value obtained from the numerical simulation was 7,900 mm, which resulted in a subsidence coefficient of q = 0.79. The subsidence coefficient q = 0.78 obtained by fitting the measured

subsidence data from the surface station with the reference software was 1.2%, and the relative errors were less than the specified limit of 5%12. Therefore, the established simulation model

of coal seam mining under specific geological mining conditions of Mengba Mine is correct. To analyze the changes of surface subsidence law after the mining with the increase of cumulative

mining thickness and the reliability of the numerical simulation, the cumulative mining thickness of 5 m, 10 m, 15 m, 20 m, 25 m, and 30 m were simulated under the same conditions as the

actual mining of single seam, mining thickness, and mining method. The curves of surface subsidence, tilt, curvature, horizontal movement, and horizontal deformation with increasing

cumulative mining thickness are shown in Figs. 3, 4, 5, 6, 7. From the sinking curve simulating the change in mining thickness in Figs. 3, 4, 5, 6, 7, it can be seen that the surface sinking

curve had flattened out after the coal seam recovery, indicating that the strike direction had reached the full mining level. There are two maximum points on the slope curve, a positive

slope maximum point, and a negative slope maximum point, where the slope value of the mining center is zero. The curvature curve finally appears with three maximum curvatures; two positive

maximum curvatures appear outward from the protected coal pillar; a negative maximum curvature appears between the protected coal pillar and the mining area; two maximum points of horizontal

movement appear in the horizontal movement curve; three maximum points of deformation values appear in the horizontal deformation curve; a maximum negative point of horizontal deformation

value appears between the protected coal pillar and the center of the mining area; a pair of maximum positive points of horizontal deformation value appear outward from the protected coal

pillar, and the horizontal movement and horizontal deformation curves are similar to the tilt and curvature deformation curves. As the cumulative mining thickness increases, the values of

surface subsidence, tilt, curvature, horizontal movement, and horizontal deformation all increase accordingly. In the sinking curve, the sinking value of the ground eventually reaches the

maximum at the center of the mining area. The further away from the center of the mining area, the smaller the value of subsidence of the surface. After the mining, the coal seam can be

divided from the center of the mining area upwards in order: caving zone, fractured zone, and bending subsidence zone. In this paper, the subsidence law of the ground surface after coal seam

mining was investigated, so only the overburden deformation above the simulated excavation and the displacement of the ground surface in two directions, X and Y were extracted from the

simulation results. Figure 8 shows the displacement clouds of overlying strata corresponding to different thicknesses. Afturer mining, the direct top plate was the first to collapse with the

maximum movement and deformation value in the vertical direction. The further from the direct top plate, the smaller the deformation value of the overlying rock layer. Meanwhile, the

deformation of the overlying rock layer and the sinking value of the ground gradually increases with the increase of the accumulative mining thickness, and the range of the central part of

the ground surface sinking basin gradually increases as indicated by Fig. 8a,f. In the horizontal direction, the extent of the basin where the surface is affected by mining after coal seam

mining increases gradually with the increase of cumulative mining thickness. RESULT ANALYSIS OF SUBSIDENCE COEFFICIENT According to the numerical simulation results, the maximum surface

subsidence and sinking coefficient of different total mining thicknesses are shown in Table 4. For overlying rock formations with different properties, the activation coefficient for

subsidence under each repeated mining condition can be calculated according to the following equation: $$ q_{r1} = (1 + {\text{a}})q_{i} $$ (1) $$ q_{{r{2}}} = (1 + {\text{a}})q_{{r{1}}} $$

(2) where a is the repeated mining sinking activation coefficient,\(q_{i}\),\(q_{r1}\),\(q_{r2}\) correspond to the corresponding sinking coefficients of initial mining, primary repeated

mining, and secondary repeated mining, respectively. According to the above equation, changes in the activation coefficient of subsidence during repeated mining are shown in Table 5.

According to the numerical simulation results in Tables 3 and 4, the following conclusions could be drawn: (1) During the slicing mining of the thick coal seam, the sinking coefficient of

the first repeated mining increased by 20% when the accumulative mining thickness was less than 15 m. The sinking coefficient of the secondary repeated mining increased by 10%, and the

sinking coefficient of the subsequent repetitive mining no longer increased13. The cumulative mining thickness of the Mengba mine is 36 m, and the numerical simulation results showed that

the surface subsidence coefficient was 0.79 after the secondary repeated mining, an increase of 11% relative to the primary repeated mining; the surface subsidence coefficient was 0.84 after

the tertiary repeated mining, an increase of 7% relative to the secondary repeated mining; the surface subsidence coefficient was 0.89 after the quaternary repeated mining, an increase of

6% relative to the tertiary repeated mining; the surface subsidence coefficient was 0.92 after the fifth repeated mining, an increase of 3% relative to the quaternary repeated mining. After

the cumulative mining thickness reached 25 m, q is intended to be a fixed value and keep constant. Furthermore, the repeated mining coefficients until the subsidence coefficient no longer

increased after the fifth repeated mining. (2) The activation coefficient a1 of sinking for primary repeated mining was 0.11, a2 for the secondary repeated mining was 0.06, a3 for the

tertiary repeated mining was 0.06, and a4 for the quaternary repeated mining was 0.03. When the repeated mining time was more than five times, the activation coefficient of sinking no longer

changed. (3) The relationship between the surface subsidence coefficient and the cumulative thickness obtained by numerical analysis software is shown in Fig. 9: The relationship between

the surface subsidence coefficient and the cumulative thickness is as follows: $$ {\text{q}} = {0}{\text{.0426}}\sum {{\text{m}} + {0}{\text{.696 }}} {\text{m}} \le {25} $$ (3) $$ {\text{q}}

= {0}{\text{.692 m}} > {25} $$ (4) where \(q\) is the subsidence coefficient,\(\sum m\) is the cumulative thickness, m. THEORETICAL ANALYSIS After the mining, the direct roof plate is

subject to downward movement and deformation by gravity and the overburdened rock. Rock movement produces different forms of deformation in different areas, and according to the different

forms of damage, the damaged rock layer above the goaf can be divided into three parts, as shown in Fig. 10. In the caving zone where the rock layer collapses fully, the gap between the

collapsed rocks increases and the volume expands. The volume of the rock after the span is larger than that before the collapse, and the collapsed rock layer fills the mining area. The

overlying rock layer above the caving zone produces cracks, departures, and fractures of a certain width, but the zone still maintains the original laminated structure, which is called the

fractured zone. Water and sand can easily enter the mining area through the fracture zone, so when mining coal underwater, special attention should be paid to the development height of the

fracture zone. The bending zone refers to the area above the fracture zone to the surface, where the rock layers maintain a monolithic and laminated structure with good water barrier

performance. The surface movement process in the bending zone is continuous and regular, and there are few off-layer cracks, which make it difficult for water and sand to pass through. There

are two main characteristics of slicing mining: the cumulative thickness of mining is large, and the surface is subject to repeated mining. After the destruction of the direct roof of the

coal seam during the primary mining, the overlying rocks in the entire mining area are damaged and deformed, and the collapsed rocks of the overburdened rocks fill the mining void area, and

the surface is deformed to a certain extent. The ruptured rocks of the lower coal seam damaged by secondary mining will cause the overburdened rocks to continue to collapse and fill the

mining void area so that the gap between the rocks destroyed by the primary mining will be closed gradually, which increases the surface movement and deformation compared with the primary

mining14,15,16,17,18. However, when the accumulated mining thickness reaches a certain amount, the gap inside the rock formation will be gradually compacted and will no longer increase with

the mining operations. The number of repeated mining is related to the lithology of the overburdened rock on the roof. When the coal seam is repeatedly mined, there is a rock group with

thickness separating the different layers between the mined coal seams, without disturbance during the primary mining, and the undisturbed rock layer plays a certain protective role for its

upper layer during the repeated mining, thus making the repeated mining sinking coefficient gradually decrease with the repeated mining. When the coal seam is repeatedly mined, there is a

rock group with thickness separating the different layers between the mined coal seams, without disturbance during the primary mining, and the undisturbed rock layer plays a certain

protective role for its upper layer during the repeated mining, thus making the repeated mining sinking coefficient gradually decrease with the repeated mining. When the mining thickness

reaches a certain value, all layers encounter destructive deformation, the protective effect gradually weakens or disappears19,20,21,22,23, and the repeated mining sink coefficient will

gradually tend to zero and no longer change. According to the numerical simulation results, the subsidence coefficient of the ground gradually increased after the second, third, and fourth

strata mining compared with the previous layer, and when the accumulated mining thickness reached 25 m, the subsidence coefficient of the ground surface tended to a constant value of 0.92.

The experimental results are the same as the theoretical analysis. CONCLUSION * (1) The subsidence coefficient q of the surface exhibits segmented characteristics as the cumulative mining

thickness increases. Numerical simulation results show that the surface subsidence coefficient was 0.79 after the secondary slicing repeated mining, an increase of 11% than that after the

primary slicing repeated mining; the surface subsidence coefficient was 0.84 after the tertiary slicing repeated mining, an increase of 7% than that after the secondary slicing repeated

mining; the subsidence coefficient was 0.89 after the quaternary slicing repeated mining, an increase of 6% than that after the tertiary slicing repeated mining; the surface subsidence

coefficient was 0.92 after the fifth slicing repeated mining, an increase of 3% than that after the quaternary slicing repeated mining. After the cumulative mining thickness reached 25 m, q

is intended to be a fixed value and keep constant. And the repeated mining coefficients until the subsidence coefficient no longer increased after the fifth repeated mining. The activation

coefficient a1 of sinking for primary repeated mining was 0.11, a2 for the secondary repeated mining was 0.06, a3 for the tertiary repeated mining was 0.06, and a4 for the quaternary

repeated mining was 0.03. When the repeated mining time was more than five times, the activation coefficient of sinking no longer changed. * (2) The calculation formula between the ground

subsidence coefficient and cumulative mining thickness of layered mining under the condition of the medium-hard roof was fitted, which is expected to provide a parameter basis for the mining

of coal seams in similar mines. DATA AVAILABILITY The datasets used and/or analysed during the current study available from the corresponding author on reasonable request. REFERENCES *

Salmi, E. F., Nazem, M. & Karakus, M. Numerical analysis of a large landslide induced by coal mining subsidence. _Eng. Geol._ 217, 141–152 (2017). Article  Google Scholar  * Erarslan, N.

Microstructural investigation of subcritical crack propagation and Fracture Process Zone (FPZ) by the reduction of rock fracture toughness under cyclic loading. _Eng. Geol._ 208, 181–190

(2016). Article  Google Scholar  * Yja, B. _et al._ A new approach to estimate fugitive methane emissions from coal mining in China. _Sci. Total Environ._ 543, 514–523.

https://doi.org/10.1016/j.scitotenv.2015.11.024 (2016). Article  ADS  CAS  Google Scholar  * Harun-Or-Rashid, S. M. _et al._ Impact of coal mining on soil, water and agricultural crop

production: A cross-sectional study on Barapukuria coal mine industry, Dinajpur, Bangladesh. _J. Environ. Sci. Res_ 1(1), 1–6 (2014). Google Scholar  * Zhenfeng, L. I. & Jin, X.

Numerical simulation research on scope division of “three zones” in roof with UDEC. _Min. Saf. Environ. Protection_ 42, 21–24 (2015). Google Scholar  * Tan, Z. X., Wang, Z. S., Li, Y. J.,

Sun, X. M. & Deng, K. Z. Field research on ground subsidence rules of intensive fully-mechanized mining by sublevel caving. _J. Min. Saf. Eng._ 1, 59–62 (2008). Google Scholar  * Ren,

G., Li, J. & Buckeridge, J. Calculation of mining subsidence and ground principal strains using generalised influence function method. _Trans. Instit. Min. Metall. Sect. A: Min.

Technol._ 119(1), 34–41 (2010). Google Scholar  * Liu, Y. K., Zhou, F. B., Liu, L., Liu, C. & Hu, S. Y. An experimental and numerical investigation on the deformation of overlying coal

seams above double-seam extraction for controlling coal mine methane emissions. _Int. J. Coal Geol._ 87(2), 139–149 (2011). Article  CAS  Google Scholar  * Ren, G., Li, G. & Kulessa, M.

Application of a generalised influence function method for subsidence prediction in multi-seam longwall extraction. _Geotech. Geol. Eng._ 32(4), 1123–1131 (2014). Article  Google Scholar  *

Sui, W., Hang, Y., Ma, L., Wu, Z. & Wei, L. Interactions of overburden failure zones due to multiple-seam mining using longwall caving. _Bull. Eng. Geol. Environ._ 74, 1–17 (2014).

Google Scholar  * Salmi, E., Nazem, M. & Karakus, M. The effect of rock mass gradual deterioration on the mechanism of post-mining subsidence over shallow abandoned coal mines. _Int. J.

Rock Mech. Min. Sci._ 91, 59–71 (2017). Article  Google Scholar  * Helm, P. R., Davie, C. T. & Glendinning, S. Numerical modelling of shallow abandoned mine working subsidence affecting

transport infrastructure. _Eng. Geol._ 154, 6–19 (2013). Article  Google Scholar  * Suchowerska, A. M., Merifield, R. S. & Carter, J. P. Vertical stress changes in multi-seam mining

under supercritical longwall panels - ScienceDirect. _Int. J. Rock Mech. Min. Sci._ 61, 306–320 (2013). Article  Google Scholar  * Coggan, J., Gao, F., Stead, D. & Elmo, D. Numerical

modelling of the effects of weak immediate roof lithology on coal mine roadway stability. _Int. J. Coal Geol._ 90–91, 100–109 (2012). Article  Google Scholar  * Xu, N. _et al._ Surface

subsidence prediction for the WUTONG mine using a 3-D finite difference method. _Comput. Geotech._ 48, 134–145 (2013). Article  Google Scholar  * Yu, B., Zhao, J., Kuang, T. & Meng, X.

In situ investigations into overburden failures of a super-thick coal seam for longwall top coal caving. _Int. J. Rock Mech. Min. Sci._ 78, 155–162 (2015). Article  Google Scholar  * Yang,

L., Kong, D., Wu, G. & Zhen, Q. M. Development and application of a physical model for longwall coal face failure simulation. _Int. J. Min. Mineral Eng._ 8(2), 131–145 (2017). Article 

CAS  Google Scholar  * Mt, A., Lh, B., Hx, A. & Qm, B. Experimental study of deformations and failures of the coal wall in a longwall working face. _Eng. Fail. Anal._ 125, 105428 (2021).

Article  Google Scholar  * Kong, D. Z., Cheng, Z. B. & Zheng, S. S. Study on the failure mechanism and stability control measures in a large-cutting-height coal mining face with a

deep-buried seam. _Bull. Eng. Geol. Environ._ 78(2), 6143 (2019). Article  Google Scholar  * Behera, B., Yadav, A., Singh, G. & Sharma, S. A numerical modeling approach for evaluation of

spalling associated face instability in longwall workings under massive sandstone roof. _Eng. Fail. Anal._ 117, 104927 (2020). Article  Google Scholar  * Li, L., Sun, S., Wang, J., Yang, W.

& Fang, Z. Experimental study of the precursor information of the water inrush in shield tunnels due to the proximity of a water-filled cave. _Int. J. Rock Mech. Min. Sci._ 130, 104320

(2020). Article  Google Scholar  * Xue, Y. _et al._ Effect of water-cooling shock on fracture initiation and morphology of high-temperature granite: Application of hydraulic fracturing to

enhanced geothermal systems. _Appl. Energy_ 337, 120858 (2023). Article  Google Scholar  * Xue, Y. _et al._ Influence mechanism of brine-gas two-phase flow on sealing property of anisotropic

caprock for hydrogen and carbon energy underground storage. _Int. J. Hydrogen Energy_ 48(30), 11287–11302 (2023). Article  CAS  Google Scholar  Download references ACKNOWLEDGEMENTS The

financial support by the Found of the National Key Research and Development Program and Key Special Fund Project (No.2018YFC0808202) are gratefully acknowledged. AUTHOR INFORMATION AUTHORS

AND AFFILIATIONS * China University of Mining and Technology, Beijing, 100083, China Xiugang Liu & Zaibing Jiang * China Coal Research Institute, Beijing, 100013, China Xiugang Liu &

 Zaibing Jiang * CCTEG Xi’an Research Institute (Group)Co, Ltd, Xi’an, 710077, China Xiugang Liu, Zaibing Jiang, Yi Wang & Jie Zhang * National Time Service Center, Chinese Academy of

Science, Xi’an, 710600, China Fei Wei * University of Chinese Academy of Science, Beijing, 100049, China Fei Wei * JiangSu Key Laboratory of Resource and Environment Information Engineering,

China University of Mining and Technology, Xuzhou, 221116, China Zhixiang Tan Authors * Xiugang Liu View author publications You can also search for this author inPubMed Google Scholar *

Fei Wei View author publications You can also search for this author inPubMed Google Scholar * Zhixiang Tan View author publications You can also search for this author inPubMed Google

Scholar * Zaibing Jiang View author publications You can also search for this author inPubMed Google Scholar * Yi Wang View author publications You can also search for this author inPubMed 

Google Scholar * Jie Zhang View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS X.L. conceived the study and, together with F.W., Z.T. Z.J.,

Y.W., and J.Z. did the literature search, selected the studies. X.L. and F.W. extracted the relevant information. X.L. synthesised the data. F.W.conducted numerical simulation.X.L.and

F.W.wrote the first drafts of the paper. Z.T., Z.J., Y.W., and J.Z. critically revised successive drafs of the paper. All authors approved the final drafts of the manuscript. X.L. is the

guarantor of the study. CORRESPONDING AUTHOR Correspondence to Xiugang Liu. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION

PUBLISHER'S NOTE Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. RIGHTS AND PERMISSIONS OPEN ACCESS This article

is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give

appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in

this article are included in the article's Creatife Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's

Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Liu, X., Wei, F., Tan, Z. _et al._ Study on

the law of surface subsidence in layered mining of thick coal seam with medium hard roof. _Sci Rep_ 13, 14735 (2023). https://doi.org/10.1038/s41598-023-42012-5 Download citation * Received:

20 April 2023 * Accepted: 04 September 2023 * Published: 07 September 2023 * DOI: https://doi.org/10.1038/s41598-023-42012-5 SHARE THIS ARTICLE Anyone you share the following link with will

be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt

content-sharing initiative