Movement Boundary Shape of Overburden Strata and Its Inﬂuencing Factors

: Strata movement boundary is not only a parameter for the prediction of overburden strata movement and deformation but also a key index of setting shafts, roadways and protective coal pillars. Based on physical and mechanical properties of rock mass, the overburden strata are divided into bedrock and unconsolidated stratum. By means of theoretical analysis, physical simulation and numerical simulation, this paper studies the movement boundary shapes of bedrock and unconsolidated stratum, builds ﬁtting equations of movement boundary of the two, analyzes the inﬂuence of key strata (KS) on the shape of strata movement boundary, and determines the principle of setting protective coal pillars. The results show that the movement boundaries of bedrock and unconsolidated strata are located at the outside of coal mining boundary. They are concave-upward power function curves that cannot be merged into a smooth one due to their different mechanisms of movement and deformation. The movement boundary of bedrock can approximate a straight line when lithology of the overburden is relatively uniform with thin strata in different positions; the surface movement boundary extends when the overburden has thick and stiff KS that are common in deeply buried coal seam. Therefore, the width of protective coal pillar is small if the movement boundary is regarded as a straight line. According to the curve movement boundary, the protective coal pillar for the passenger roadway of Panel 31010 of Pingdingshan No.1 mine is at least 99.4 m in width, larger than the designed one, which is the actual reason for its deformation and breakage.


Introduction
Energy resource is the basis for human survival and energy exploitation has made an important contribution to the development of national economy. However, most underground energy exploitation will lead to roof caving, overburden fracturing and surface subsidence, resulting in environmental damage and roadway convergence etc. [1][2][3]. Therefore, it is especially important to study the strata movement caused by coal mining. Nevertheless, in many cases, the protection of ground buildings and roadways only requires the determination of overburden movement boundary, namely the delimitation of the risk of collapse, rather than the movement and deformation values of the whole areas affected.
reveal that the features of unconsolidated stratum are similar to that of granules, while bedrock is hard and compact, close to continuum. Therefore, the study of the overburden movement boundary should cover two parts: bedrock and unconsolidated stratum, and particles of the latter move after the deformation of the interface of bedrock top. The simplest methods of forecasting deformations are methods based on empirical formulas [4], defined for different ways of exploitation and the position of the excavation and the selected field of a given size. Models-prognostic theories can be divided into four groups: geometric-integral, stochastic, based on rock mass (continuous medium) and numerical models based on discretization of the medium. Geometric-integral theories referring to asymptotic deformations assume the existence of the so-called the influence function, which determines the distribution of the decreases caused by the elemental (about infinitesimal volume) depending on the horizontal distance of the point of the object from it. The reduction due to the selection of a finite field is the sum of the decreases from all elementary exploitation of this field. Horizontal dislocations are described on the basis of additional assumptions defining mutual relations with depressions. This group of methods includes the theories of Knothe-Budryk [25] and Ehrhardt-Sauer [26]. They differ mainly in the shape of the influence function. Stochastic models treat a rock mass as a center subject to random laws (Litwiniszyn theory). So far, they have not found practical applications due to the complexity of the description and a large number of parameters. Models based on rock mass mechanics-the mechanics of continuous medium (on differential equations) are most often used to locally describe the impact of mining excavations in the rock mass in their small surroundings, which results from the need to simplify the rock mass structure and interactions in order to obtain analytical solutions (Sałustowicz's work [27]).

Movement Boundary Shape of Bedrock
Any stratum of the overburden is treated as a homogeneous isotropic plane strain beam, and meanwhile the height, as well as mechanical properties of it, is regarded to be constant axially. The top and bottom surfaces of it are influenced by friction of opposite direction. The friction coefficient, supporting force and loads of both surfaces are approximately equal, so the effect of friction together with axial force can be ignored. The subsidence of the beam is free from the impact of the uniformly distributed loads at its upper part because they have already existed in the process of sediment. Therefore, the beam is only affected by end bending moment and shear force, as shown in Figure 2.
Assuming that the deformation of foundation conforms to Winkler's assumption, the differential equation of the beam is as follow: where E i and I i are the elastic modulus and the section moment of inertia of the ith stratum, respectively; E if is the comprehensive elastic modulus of the foundation; H i is the distance from the ith stratum to the coal floor. (1) is: The general solution of the fourth order ordinary differential Equation (2) is: There is a boundary condition w i (∞) = 0 at infinity due to the limited influence of coal mining. From Equation (3), According to the relationship between the internal force and the beam subsidence: The boundary condition corresponding to the origin of coordinates is M i = −E i I i y 0 , F i = −E i I i y 0 . They are introduced into Equation (5) to obtain the equation of the beam subsidence expressed by boundary conditions: Trigonometrical transform is done on Equation (6) to get: where ϕ i is the initial phase of the displacement wave, cos . Based on Equation (7), under the action of bending moment and shear force, vertical displacement of strata is an attenuation cosine function in accordance with the law of negative exponent function e −β i x , and it can be neglected after half a cycle. Equation (7) is the subsidence of the ith stratum under its own loads without the combination of lower strata subsidence. When combining, the subsidence of the ith stratum is: The subsidence range of strata is calculated by viewing subsidence of the first solution as the movement boundary of the ith stratum. Then According to the summation principle of trigonometric functions, the angular frequency of a trigonometric function after the addition of trigonometric functions with different angular frequency is that of the smallest one before the sum, i.e., β = min β k (k = 1 . . . i). With the rise of strata level, β is decreasing, so β i < β i−1, and then β = β i . At this point, Equation (9) is: According to the geometric relationship between the angle of draw γ i and the depth H i of the ith stratum (see Figure 2), the difference between angles of draw of two adjacent strata is as follows: where ϕ i is the initial phase of ith stratum with the combination of lower strata subsidence.
The relationship between β i and H i can be simplified as i . β i is brought into Equation (11) to gain: Equation (12) can be simplified to get: where max k i H 3 4 i , k i−1 H 3 4 i−1 is the maximum value of the two. The subsidence of the overburden on the mining boundary is in line with the cosine function of initial phase, and it keeps increasing as the calculated strata approach the surface due to the combination of subsidence. As a result, ϕ i keeps falling, so ϕ i < ϕ i−1 . From Equation (13), ∆γ i < 0, so γ i < γ i−1 .
The relationship between the angle of draw and the curve pattern is shown in Figure 3. If CD is the strata movement boundary, the angle of draw will not change with different strata level when points C and O coincide, so the movement boundary is not a straight line. If the curve DAC is the strata movement boundary, it can be seen from Figure 3 that γ 1 < γ 2 . For the convex-upward curve DAC, the closer the strata to the surface after a certain stratum, the greater the angle of draw, so the movement boundary is not a convex-upward curve. If the curve DBC is the strata movement boundary, it can be found that γ 3 < γ 2 . For the concave-upward curve DBC, the closer the strata to the surface, the smaller the angle of draw, so the movement boundary is a concave-upward curve.

Movement Boundary Shape of Topsoil
The unconsolidated stratum is loose deposits formed by the rock that has not been hardened after the process of erosion-transport-deposition. They can be regarded as granular media. The velocity field and stress field of granular media have been investigated for long due to the wide application of its storage and usage to industry [28][29][30][31]. For incompressible granular media without cohesion, there are two main theoretical models to analyze their velocity field, stress field and movement boundary in silos, namely, elastic-perfectly plastic model and pure kinematic model [29,[32][33][34][35][36]. The boundary shapes of granular media in the process of discharge are found to be concave-upward curves [35][36][37][38][39][40][41][42], like parabolic [35], hyperbolic [41] or elliptic [40]. However, all of these results are based on a percentage of moving particles with the maximum speed as the movement boundary, which is not the real static boundary. In addition, the boundary conditions of unconsolidated stratum in coal seam differ from that of silos in industry.
According to the abundant data from surface observation stations at mining face, it can be speculated that the subsidence scope of unconsolidated stratum is larger than the mining area, but the movement boundary shape of them cannot be determined. Studies in laboratory and on site show that different from bedrock, unconsolidated stratum possesses mechanical properties and movement and deformation laws which are close to those of stochastic media. Based on the stochastic medium (SM) movement theory [43] by Litwiniszyn, a Polish scholar, the theoretical model of the probability integral method [44] is improved, as presented in Figure 4. The differential equation of any particle's subsidence in two-dimensional SM theoretical model is as follow: Equation (14) is a second order parabolic partial differential equation. In order to obtain its solution, the following boundary conditions are listed according to practice: where δ(x) is the Dirac number which is defined as: Based on the solution to the second order parabolic partial differential equation, Equation (17) can be derived: Equation (17) indicates that subsidence varies with strata, and the area affected by mining is theoretically infinite. To obtain the movement boundary shape of unconsolidated stratum, the size of movement boundary of different stratum can be obtained by taking a percentage of maximum subsidence as the boundary: B is a positive constant. Based on second derivative of Equation (18) together with the criteria theorem of concavity and convexity [45], it can be known that the movement boundary of unconsolidated stratum is a concave-upward curve, but it may not be the parabola presented by Equaiton (18) due to its viscosity more or less.

Parameter Study
The overburden above coal seam deposits in layers with different heights and mechanical properties of each layer. The thick and hard strata controlling the movement of the overburden are KS, among which the uppermost one is the primary key strata (PKS) [46]. Field measurements and simulated studies have illustrated that KS is capable of controlling the dynamic subsidence of surface [47,48], and it is synchronous with the subsidence and subsidence rate of the controlled overburden strata. Besides, it also exerts a significant influence on strata movement boundary.
The influence of KS on strata movement boundary can be compared through the movement of strata in the case of either KS or soft rock in the same layer, respectively. Based on Equation (9), strata movement boundary is related to the elastic modulus E, height h, initial phase ϕ of the calculated strata, the height H from it to the coal seam floor, as well as the elastic modulus of foundation E f . The ratio of elastic modulus of KS to that of foundation is set to be m. As can be found in Figure 5, the larger the values of m and h, the greater the size of movement boundary. In other words, the KS will expand the scope of strata movement.
For strata with same geological conditions, yet different lithology at the same layer, the initial phase remains unchanged while the elastic modulus and the height of KS increase, thus contributing to the expansion of movement compared with soft rock. The distance from KS movement boundary to the mining boundary rises with the increase of elastic modulus and height of KS. According to KS theory, the elastic modulus and the height of KS are greater than that of other strata. Therefore, the size of movement boundary of strata expands when there is typical KS in the overburden. The enlargement of subsidence area of KS inevitably leads to the increase of that of surface because the movement characteristics of strata above KS are determined by their physical and mechanical properties. This is consistent with the statistical data of observation of surface subsidence in numerous mines in the literature [18], offering a theoretical explanation for the expansion of surface subsidence when the overburden contains hard and stiff strata. The gradient of angle of draw is defined as the ratio of the difference between the angles of draw of two adjacent layers to the vertical distance between them. When the points C and O in Figure 3 coincide, the corresponding gradient of angle of draw is zero for straight movement boundary, so is the difference between angles of draw of two adjacent layers. From Equations (11) and (12), a smaller height of strata and difference between elastic modulus of two adjacent layers is conducive to eliminating the difference between angles of draw of them. In other words, strata movement boundary can be approximate to a straight line when the lithology is homogeneous and the height of different strata are small in the overburden (that is, there is no typical KS in it).
Although mining factors such as mining height and mining width have an impact on strata movement boundary, their influence is closely related to the overburden geological conditions. Surface subsidence scope can be little with small mining width in spite of large mining height, demonstrating that mining factors are inseparable from geological factors. Hence, the effect of mining factors on strata movement boundary is not studied in this paper.

Movement Boundary Shape and Empirical Model
It is of practical significance to determine the movement boundary equation of rock strata for protecting surface structures, arranging underground roadway and designing staggered distance between the mining face of the upper and lower coal seam. As the concept of green mining gets well accepted, the requirements for ground structures protection become higher, so greater importance is attached to the prediction of surface subsidence. Nevertheless, in many cases, the protection of ground buildings only requires the determination of surface movement boundary, rather than the prediction of the movement and deformation value of the ground affected by coal mining. The movement and deformation of bedrock should be aware before underground water protection and structures construction. The prediction of strata movement and deformation demands the position of strata movement boundary, whatever method is adopted. Owing to various characteristics of the overburden, strata movement boundaries are also distinguishable in different mining areas.
According to the movement boundary shapes of bedrock and unconsolidated stratum, it is assumed that the distance from the movement boundary of the two to the mining boundary and the vertical distance to the coal seam satisfy the power function relationship. The mining boundary is defined as coordinate origin to establish the coordinate system, as shown in Figure 6. The movement boundary equations of bedrock and unconsolidated stratum are Equations (19) and (20), respectively: [r(y) − r 0 ] n = By y < H r (19) [r(y) where B, B 1 , n and n 1 are coefficients, and n with n 1 are greater than 1. r is the main influence radius of surface; r 1 is that of the interface of bedrock top; and r 0 is that of coal seam. H r is the bedrock height. The whole overburden strata can be expressed by Equation (20) when unconsolidated stratum is thin or the movement boundaries of bedrock and unconsolidated stratum are continuously smooth. The exact value of unknown parameters can be derived according to back analysis of measured data, and their approximate value can be calculated based on simulation results or Equation (9). The physical significance of parameters in Equations (19) and (20) is shown in Figure 7 where r 0 = 0. As can be seen from Figure 7a, with the increase of n, the movement boundary becomes steeper and steeper, and the area of subsidence shrinks. From Figure 7b, the larger the value of B, the smoother the movement boundary, and the larger the area of subsidence.

Physical Simulation
In order to verify the shape of the overburden movement boundary, the following experimental scheme was designed with physical and mechanical parameters of bedrock and unconsolidated stratum conforming to practical situation. The model in the physical simulation was 250 cm in length, 20 cm in width and 123 cm in height. According to similar material proportioning test of unconsolidated stratum, the geometric similarity ratio was 100:1 with density similarity ratio of 1.67:1 and stress similarity ratio of 167:1 in the simulation. In this experiment, river sand was taken as the basic aggregate, while gypsum, calcium carbonate and mica powder as auxiliary materials. The physical and mechanical parameters of each stratum in the model were determined on the basis of the similarity theory, as listed in Table 1. In the model, 25 horizontal displacement measurement lines were arranged, of which one is in the coal seam, 12 in bedrock and unconsolidated strata respectively. There are 49 displacement measurement points on each line with the horizontal and vertical distance of 5 cm between two adjacent points. A total of 1225 displacement measurement points were arranged, as exhibited in Figure 8. In the process of coal stoping, DigiMetric 3d (Beijing TenYoun 3D Technology Co., Ltd., Beijing, China) photogrammetry system was adopted to monitor displacement with theoretical measurement accuracy 0.1 mm/4 m.
On the right side of the model, a 30-cm-wide protective coal pillar for boundary was reserved. The mining was carried out from right to left with single excavation space of 5 cm Displacement measurement and data processing were conducted every 10 cm of stoping. Further excavation proceeded when all points cease to move.   Figure 9 shows strata movement boundary with mining width of 125 m and 135 m, respectively. The blue and red curves represent the 20 mm and 30 mm strata movement boundary, respectively (The subsidence value has been converted according to the similarity ratio.). It is evident that the movement boundary shapes of bedrock and unconsolidated stratum are concave-upward curves, no matter which subsidence isoline is regarded as the movement boundary. The shapes of two subsidence isolines are consistent in the same figure. To illustrate whether Equation (16) can be considered as the equation of strata movement boundary, the movement boundaries of 20 mm bedrock and unconsolidated stratum are fitted, respectively. The result of bedrock (a correlation coefficient of 0.96) is shown in Figure 10a, while that of unconsolidated stratum (a correlation coefficient of 0.98) is in Figure 10b. As can be observed from Figure 9, the strata movement boundary is different in bedrock and unconsolidated stratum.
In Figure 9a, the change of movement boundary of unconsolidated stratum with layers is not as distinct as that of bedrock due to the arch structure in it. In Figure 9b, the arch structure has been destroyed and the movement of unconsolidated stratum approximately conforms to probability integral principle. The strata movement boundary turns into concave and convex shape at PKS where bedrock and unconsolidated stratum obey different movement and deformation laws. Therefore, there are differences in the movement boundary equations of the two (Figure 10a,b).

Numerical Simulation
The finite-difference software FLAC3D (5.0 64-bit, Itasca Consulting China, Ltd., Wuhan, China) was introduced in this section to simulate the movement boundary of bedrock for further analysis of strata movement boundary shape and the influence of KS on strata movement boundary. Six models were designed in two groups. Models 1-3 in Group 1 were used to study the effect of elastic modulus of KS on the boundary. Among them, Model 1 was a basic model with parameters shown in Table 2; Models 2 and 3 were the same with Model 1 in size, but their elastic moduli of KS were 3 and 5 times that of Model 1, respectively. Models 4-6 in Group 2 were adopted to research the effect of KS height on the boundary. Among them, Model 4 was a basic model with parameters shown in Table 2; Models 5 and 6 were the same with Model 1 in size, but their KS heights were 3 and 5 times that of Model 4, respectively. The calculated model was 1000 m long, 4 m wide and 233 m high, and the coal seam was 3 m in horizontal height. Mohr-Coulomb criterion was adopted for constitutive relationship of the model. At the bottom of the model, the vertical displacement was fixed, and the horizontal displacement was restricted by boundaries in four directions with the upper boundary a free surface.
The model was divided into five strata from the bottom up, namely coal floor, coal seam, immediate roof, KS and soft rock. The physical and mechanical parameters of each stratum in basic model are given in Table 2. After the completion of model, balancing calculation was performed under the effect of gravity, and then the monitoring variables were cleared. The working face was mined for 300 m in width at one time to obtain balanced default. Regulations of coal pillar design and extraction for buildings, water bodies, railways, main shafts and roadways (hereafter referred to the Regulations) specify that the point with subsidence value of 10 mm is the boundary point [2]. Hence, all these points are linked to form the movement boundary of overburden strata. The subsidence isolines formed by these points in six groups are shown in Figure 11a,b and only a part of strata in the length direction is drawn for the need of demonstration. It can be seen from Figure 11a,b that regardless of the elastic modulus and height of KS in the overburden, the movement boundary is a concave-upward curve. In addition, the shape of strata movement boundary remains unchanged while the movement scope increases when elastic modulus and height of KS become larger. Figure 11c exhibits the strata movement boundary for the overburden with same physical and mechanical parameters (no KS). As can be observed in Figure 11c, the boundary is a concave-upward curve, but it is more like a straight line compared with that of the overburden with KS. The numerical simulation results of the influence of KS on strata movement boundary are consistent with the theoretical analysis in Section 3. In fact, there is little difference between the elastic modulus and height of KS compared with that of other strata in practices, so the n of bedrock movement boundary can take a small value in general geological condition.

Design of Protective Coal Pillars
Strata movement boundary is the line to judge whether the overburden is affected by coal stoping. Inside of it, the displacement of the overburden changes under the impact of stoping, whereas outside of it, the overburden remains basically intact. As illustrated in Figure 12, the movement boundaries of bedrock and unconsolidated strata extend in a concave way from the coal seam to surface. The higher the coal seam, the larger the horizontal distance between movement boundary and mining boundary, yet the smaller the increasing trend of the distance. Therefore, the width of protective coal pillars should be designed according to strata movement boundary shape and the position of protectors to ensure that the roof roadway or underground structures are not affected by stoping. Specific procedures are as follows: calculating and drawing movement boundary based on actual situation after stoping; determining the intersection of the level of protected objects and strata movement boundary, then making its projection point on the coal seam; the distance between mining face boundary and projection point is the width of the protective coal pillar.
The overburden movement boundary after stoping is demonstrated in Figure 12. γ is both the included angle between traditional movement boundary and the vertical line and the angle between the line that links the intersection of curve movement boundary and surface with mining boundary and the vertical line. represents the roadway. After stoping, the roadway is located at A according to traditional method, and the width of the protective coal pillar should be a; the location of roadway is B based on curve movement boundary. It can be known from Figure 10 that a < b, and the roadway A is within the impact range of mining. Therefore, design requirements cannot be satisfied because the roadway is still affected by stoping and produces deformation when the pillar for roof roadway are designed based on traditional movement boundary. However, the pillars designed according to curve movement boundary can exclude roadway and underground structures from impact range of mining, thus better meeting engineering practice. In addition, the level of protected objects should also be taken into consideration in the design of the width of the pillar. In Figure 12, at different height, the roadway C is affected by mining while the roadway B is intact despite of the same width of pillars. Therefore, both the shape of strata movement boundary and the position of protected objects are important factors that should be considered in designing protective coal pillars.
As for the design of protective coal pillars in practices, it should be noted that different buildings and structures can resist varied hazards. And the basic factor in assessing the danger to surface construction objects or to linear objects is the horizontal specific deformation [49]. Consequently, the size of protective coal pillars can be optimized by the criterion of horizontal deformation, which can be divided into two groups. Criterion for the minimization of the resources trapped in the pillar is on the one hand and the criteria for facility protection (excavation) is on the other hand. Criteria for the protection of objects affect the size of the pillar in the direction of their growth, while the criterion of resource minimization acts in the opposite direction. Both criteria taken together optimize the size of the pillar in individual decks. It should be noted, however, that the designation of protective coal pillars, including optimization, is justified when geological and mining conditions as well as known parameters of rock mass theory are recognized.

Mining Conditions
The DingWu-3 passenger roadway of No.1 mine of Pingdingshan Tianan Coal Co., Ltd., Pingdingshan, China (hereafter referred to as passenger roadway) is 440 m in depth with loess of some 5.1 m. With a section of 4.2 m wide and 3 m high, the roadway is supported by the combination of bolt and anchor. The bolt is 20 mm in diameter and 2.2 m in length with raw and line space of 700 mm × 700 mm. The anchor is 20 mm in diameter and 7 m in length. Wu-10 coal seam contains Panels Wu10-21210 and Wu10-31010. Plan and cross-sectional views of passenger roadway and panels are shown in Figure 13.
As shown in Figure 13a, 1#~20# is the number of measurement points arranged in the passenger roadway. I~V on Wu10-31010 is the stoping position from January to May 2012, respectively, and the face mining is finished on 2 July 2012. VIII~XII on Wu-10-21210 is that from August to December 2011, respectively, and it is finished in January 2012. The passenger roadway is adjacent to DingWu-3 Eastern entry (hereafter referred to as Eastern entry) and DingWu-3 main return airway (hereafter referred to as return airway). Figure 13b shows that there is a 256-m-long protective coal pillar between Panels 21210 and 31010. To protect the passenger roadway, the pillar on the return roadway side of Panel 31010 is 90 m wide, and that on the railway roadway side of Panel 21210 is 165 m wide to protect the Eastern entry. The vertical distance between passenger roadway and Wu 10 coal seam is 85 m.

Field Measurement and Theoretical Analysis
In order to observe the deformation characteristics of passenger roadway under the impact of mining, a total of 20 monitoring points (spacing: 50 m) are arranged in it, as presented in Figure 13. In the later observation, the data of roof-floor are measured by the cross-distributing point method and its initial data are obtained on 21 January 2012.
Field measurement has been carried out on 12 February, 10 March, 6 April, 30 April, 7 July and 9 June 2012, respectively. The results of roof-floor subsidence are shown in Figure 14. After the mining of Panel 21210 stops, there is no convergence of roof-floor before the stoping of Panel 31010, showing that the mining of Panel 21210 has almost no effect on the roadway. During the coal mining of Panel 31010 until July 2012, deformation of passenger roadway occurs, so it is caused by mining of Panel 31010. It is obvious in Figure 15 that the passenger roadway is beyond the influence range of mining if it is designed according to traditional movement boundary of AC, that is, the 90-m-wide pillar designed based on AC can protect the passenger roadway from deformation. However, measured results show that there is great convergence of passenger roadway. The width of the protective coal pillar for passenger roadway should be 99.4 m which is greater than the former width of 90 m, according to the curve boundary movement proposed in this paper. To avoid deformation of the passenger roadway, the width of the pillar should be at least 99.4 m. The reason for the deformation of passenger roadway is that the traditional movement boundary is not the real strata movement boundary.
The calculation procedure with equation of curve movement boundary is given below. For lack of the observation data of surface subsidence after coal mining of Panels 21210 and 31010, the surface boundary angle approximately equals to main influence angle whose tangent is calculated according to average of 11 mining faces near the coal area based on the Regulations. The average is 1.85. The buried depth of the return airway of Panel 31010 is 538 m, so the surface impact scope of Panel 31010 is 290.8 m according to the relationship between the tangent of main influence angle and buried depth after stoping. Due to a thin unconsolidated stratum of the overburden on Panel 31010, the strata movement boundary can be indicated by a curve. In order to make a better comparison, r 0 in Equation (19) is set to zero because an endpoint of traditional movement boundary is on the mining boundary. To determine the two parameters n and B in the equation of a strata movement boundary, at least two points on the boundary are required, while only one point of C is fixed based on existing data. When n is 1.2, the movement boundary tends to be a straight line. However, the movement boundary cannot be approximated to straight owing to different lithology of each stratum in the overburden of Panel 31010 with a large depth. For better explanation of the problem, n is set to 1.4 in this paper. Put n = 1.4, r 0 = 0 and r = 290.8 into Equation (19), obtaining B = 5.23, thus the equation of movement boundary of Panel 31010 is determined, as shown in Figure 15.

Conclusions
Mining activities is inevitable in many countries. Therefore, the prediction of strata and surface movement caused by mining practices is quite important to social constructions and roadways setting. In most situations, however, only the delimitation of the risk of collapse is required. For this reason, the overburden movement boundary shape and its influencing factors was studied in this paper. The research works indicated that the movement boundary was not a traditionally regarded straight line. Based on this, an empirical model was presented to guide the design of protective pillars.
(1) Strata movement boundary is divided into two parts, namely bedrock and unconsolidated stratum. The movement boundaries of the two are both concave-upward curves. The horizontal projection distance from the movement boundaries to the mining boundary and the vertical distance to the coal seam satisfy the power function relationship, but parameters of the two movement boundaries are different. Due to the difference of sedimentary time, physical and mechanical properties of bedrock and unconsolidated stratum are quite different. The latter is similar to granules, while the former is similar to the continuum, so the same mechanical model cannot be applied to movement boundaries of the two. The movement boundary of bedrock can be analyzed with continuum mechanics, while that of unconsolidated stratum can be studied by the influence function method.
(2) KS is an important geological factor affecting the size of strata movement boundary. With typical KS in the overburden, the surface subsidence area increases with the expansion of KS subsidence range. When the lithology of the overburden is relatively uniform and the strata is thin, strata movement boundary can be approximated as a straight line.
(3) Based on strata movement boundary, the width of the protective coal pillar on Panel 31010 in Pingdingshan No.1 mine is at least 99.4 m. The traditional movement boundary is not the real strata movement boundary, and the width of the pillar designed according to it is smaller than actually needed width, which is the reason for the convergence of passenger roadway.
The proposed model can be applicable to stratified overburden or soil in underground mining, whether in coal mining or tunnel excavation. Regarding the strata inclination, however, only horizontal layering is considered in the current modeling. If the overburden strata are inclined, the movement boundary would be different from that in a horizontal overburden. Nevertheless, if the proposed model is extended to the nearly-horizontal strata (with the angle of <8 • ), it may still be applicable within the tolerable error accepted by practical engineering.