# Deformation failure characteristics of coal body and mining induced stress evolution law.

1. IntroductionWith the deterioration of occurrence condition of coal resources, mechanical environment for coal mining, organization structure of coal and rock mass, its mechanical behavior and failure characteristics become complicated, so the characteristics of engineering response change [1-3]. Meanwhile, the large space mining model complicates time-space relationship and kinetics features of mining stress field rock burst and coal and gas outburst, which causes serious damage and a number of casualties. The main reasons for these problems are lack of knowledge of the space-time evolution law of mining stress field with different mining conditions, roadway excavation and maintenance, and the advance of heading face at the wrong time and space.

The domestic and foreign scholars usually use classical mechanics, (e.g., elastic-plastic mechanics) to study the damage development process of coal body and the mining stress field that loads on the coal body during the process under the condition of mining activity [4, 5]. In this analytical method, ideal coal body is under the condition of static stress. But coal body in the stope is under dynamic loading induced by the bending and breaking of the overlaying strata, which weakens the mechanical properties of coal body and continuously transfer the mining stress inside the coal body [6]. So the quantitative study on the space-time evolution mechanism and its law of mining stress field with the deterioration development of mechanical properties of coal body is necessary [7]. It is helpful to the roadway excavation, maintenance, and the advance of the heading face, which has very important practical significance.

2. Macromechanical Structure of the Stope

The overlying strata of coal seam can be divided into two parts: the outrange of overlying strata spatial structure of stope that consists of no significant move rock beams in the stress arch that does not move significantly and the overlying strata spatial structure that is composed of kinetic rock strata in the breaking arch that affects the ground pressure of stope directly.

With the advance of heading face, the hanging space of stope continuously increases, to a certain degree, the overlying strata start to break from bottom to top, which forms the breaking arch [8-10]. At the same time, the stress in the spatial structure strata is redistributed and the gravity of overlying strata, originally supported by the mining coal body, is loaded on the two sides of the coal body and rock mass at the same level, as shown in Figure 1. The total stress in the two sides of the coal and rock mass within a specific range consists of two parts: (1) the stress caused by the gravity of overlying strata within the outrange of breaking arch and (2) the stress transferred from broken rock beam within the breaking arch of stope. If the total stress is greater than its strength, the coal and rock mass will be broken, and then the peak value of abutment pressure will be transferred inside the coal and rock body [11]. In the breaking process of every rock beam, it forms stress arch that is composed of abutment pressure peak value of every rock strata outside of the breaking arch. Its shape is parabolic in the vertical plane, as shown in Figure 1.

2.1. Characteristics of Mechanical Structure of Breaking Arch and Its Evolution Law. The rock strata in the breaking arch play a leading role on the mining pressure, so the research on the characteristics of mechanical structure of breaking arch has a significant meaning for the structural stability of mining stope.

The mining process of heading face can be divided into two stages: (1) no-full mining stage, which means the advancing length of heading face [L.sub.x] is less than that of the heading face width [L.sub.0] and (2) full mining stage, which means the advancing length of heading face [L.sub.x] is greater than that of the heading face width [L.sub.0].

In the no-full mining stage, the height of the overlying spatial structure of stope is linear to the advancing distance of the heading face, which develops forward along the coal mine strike and upward in the space. The height is about half of the short side width of goaf. But this development law is limited. From the analysis above, the height of the overlying strata spatial structure of stope is mainly determined by the width of heading face, which means the maximum height of the overlying strata spatial structure of stope is constant when the width of heading face is fixed. When the advancing length of heading face is less than that of the heading face width, the height of the stope structure is relevant to the advancing length. When the advancing length is greater than the heading face width, the height of the stope structure is half of the heading faces width.

2.2. Mechanical Characteristics of Stress Arch and Its Evolution Law. The rock stratum in the stress arch is the main bearing body, which bear and transfer the load of overlying rock strata. The breaking arch structure is located in the distressed zone of the stress arch [12, 13]. When the structure of overlying rock strata in the stress arch fails, the rock burst may occur. So clearly the dynamic evolution process of stress arch must be done.

Assuming the total amount of rock strata is K in the overlying rock structure of stope, the n + 1 rock stratum which lies on the breaking arch is called supporting stratum. When the ith rock stratum breaks, the strata load, originally supported by the ith rockstratum, is transferred to the outside of the breaking arch, as shown in Figure 2. So the load, supported by the unit length of ith rock stratum outside of the breaking arch, is

[q.sub.i] = [q.sub.1(k-i)] + [q.sub.2(k-i)] = [gamma][H.sub.i][L.sub.i] + [gamma][H.sub.i] = [gamma][H.sub.i](1 + [L.sub.i]), (1)

where [q.sub.1(k-i)] and [q.sub.2(k-i)] are the loads of ith rockstratum inner and outer field of breaking arch, respectively, [H.sub.i] is the depth of ith rock beam, and [L.sub.i] is the breaking length of ith rock beam.

The tensile failure of rock beam occurs; lots of tensile fracture may appear in the breaking edge and the load increases quickly, which decrease the strength of rock beam at the breaking position and transfer the abutment pressure peak outside of the breaking arch, as shown in Figure 3(a). It usually happens in the deep coal mines. If the strength of the rockbeam is big enough to support the load transferred from the breaking arch with no damage, the abutment pressure peak is at the breaking position, as shown in Figure 3(b). It usually happens in the shallow coal mines.

The width of the stress arch is [L.sup.stress] = [L.sub.0] +2[S.sub.e]. The height of the stress arch is [H.sub.stress] = [L.sub.0]/2 + [H.sub.n+1].

3. Mechanical Characteristics of Coal Body Destruction

On the microlevel, the coal body has lots of randomly distributed defects (e.g., pores and microcracks). With the increase of the abutment pressure, the coal body has cumulative damage. If the abutment pressure is greater than the yielding strength of the coal body, the defects propagate quickly; the coal body will become block shaped, as shown in Figure 4.

The broken coal body has irregular shape, so it is really hard to model the breaking process of coal body through mathematical method [14, 15]. Fractal method provides an idea to solve fractal dynamic behavior of geotechnical material by using fractal theory [16]. But the mechanical mechanisms of geotechnical material deformation and failure are not set up, and the relationship between fractal evolution and its physical and mechanics parameters is not built up.

3.1. Mechanical Parameters of Macroscopic Coal Body

3.1.1. Mechanical Parameter of the External Stress Field. It should be noted that the abutment pressure range in front of the stope coal wall is up to its maximum when the advancing distance of heading face equals to the width of heading face. With the advance of heading face, the range is basically constant. Its model is shown in Figure 5.

When the advancing distance reaches the width of heading face, the overlying rock strata form an increasing pressure zone with the width of due to its own gravity. Without considering the gangue bearing capacity of goaf, the following expression is built up:

(2[L.sub.0][S.sub.x] + 2[C.sub.x][S.sub.x] + 2[S.sup.2.sub.x])([K.sub.a] - 1)[gamma]H = [L.sub.0][C.sub.x][gamma]H-1/2[pi](1/2 [L.sub.0])([H.sub.g])[[gamma].sub.g][C.sub.x], (2)

where [K.sub.a] is the average value of stress concentration factor; [H.sub.g] is the height of stope stress arch with the unit of m; [C.sub.x] is equal to mining face width that is [L.sub.0].

Then,

[S.sub.x]

= 2[L.sub.0] [+ or -] [square root of (4[L.sup.2.sub.0] + 2 (([L.sup.2.sub.0] - (1/4) [pi][L.sup.2.sub.0][H.sub.g]) / ([k.sub.a] - 1) H))]/2. (3)

Simplification,

Sx [approximately equal to] ([square root (1 + 4H - [pi][H.sub.g]/8([k.sub.a] - 1)H)] - 1)[L.sub.0]. (4)

3.1.2. Mechanical Parameters of the Internal Stress Field

(1) Stress Intensity.

Based on the assumption that the vertical abutment pressure distributed within the internal stress field around the goaf equals the gravity of the basic roof strata prior to the coming of the initial pressure on the heading face, after formation of the breaking arch, the dynamic structural mechanical model of stope becomes steady, and the distribution range of internal and external stress fields remains stable. As follows in Figure 5, the internal stress field is given by

[[sigma].sub.y] max/[K.sub.max][gamma]H = [S.sub.0]/2/[S.sub.1]. (5)

Then,

[[sigma].sub.ymax] = [S.sub.0][K.sub.max][gamma]H/2[S.sub.1], (6)

where [S.sub.1] is the distance between peak abutment pressure location and coal wall with the unit of m, [S.sub.0] is the range of the internal stress field with the unit of m, H is the mining depth with the unit of m, [gamma] is the stratum unit weight of kN/[m.sup.3], and [K.sub.max] is the stress concentration factor.

(2) Distribution Range.

As can be seen in Figure 6,

1/2 [[sigma].sub.ymax] [S.sub.0] = 1/2 [H.sub.g][C.sub.i][gamma]/2, (7)

where [C.sub.i] is the cycle stress step of basic roof strata with the unit of m and [h.sub.i] is the strata thickness of ith layer with the unit of m.

Substituting (7) into the above equation, the range of the internal stress field is given by

S0 = [square root of (2[C.sub.i][H.sub.g][S.sub.1]/[K.sub.max]H)]. (8)

3.2. Mechanical Parameters of Microscopic Coal Body. Coal body damages under the condition of abutment pressure and its structural forms are very complicated and have different size. The quantitative analysis for the distribution characteristics of coal blocks is essential for the development of abutment pressure. But the stope movement and the development of abutment pressure is a dynamic process, which makes the development and change of pore size have the characteristics of space-time effect during the dynamic process. So the reasonable determination of the relationship between the pore size of the coal body and the dynamic movement of overlying rock strata is helpful to build up the transition between the microscopic quantity and macromechanics parameters.

In the breaking process of rock body, it produces lots of rock blocks with different size, but these blocks are characterized by self-similarity. Assuming the pore size of a rock is r, the distribution function is F(r,) when the compressed volume of coal wall is [epsilon]; namely, accumulative percentage of pore size of r is

F(r, [epsilon]) = N (r)/[N.sub.0], (9)

where N(r) is the amount of pores when the compressed volume of coal wall is [epsilon] and [N.sub.0] is the initial amount of the pores.

For the pores in the internal stress field, assuming segments stand for pore size, according to the length of segments, they can be divided into n rows [17]. The ith row has the amount of pores [N.sub.i] (i = 1,2, 3,..., n), as shown in Figure 7.

In Figure 7, [[summation].sup.n.sub.i=1] [N.sub.i] = [N.sub.0]; then,

N (x) = [N.sub.0]/[(r/[r.sub.max]).sup.-D], (10)

where [N.sub.0] is constant and D is the fractal dimension of pore size.

Substituting (10) into (9), in accordance with F([r.sub.max]) = 1, the expression for the distribution fractal dimension of pore size of coal body under different compressed volume is given by

F (r, [epsilon]) = [(r/[r.sub.max]).sup.-D], (11)

where r is accumulative percentage of pore size and D is the fractal dimension of pore size.

3.3. Analysis of Extension Mechanism of Internal Stress Field. On the microlevel, last section analyses the outburst mechanism of coal wall, but the parameters are hard to obtain. Based on the field application, this paper provides a brief analysis of the extension mechanism of internal stress field from the macroperspectives.

As can be seen from Figure 8, the stress states of coal body at the coal wall are biaxial compressive stress, and the external coal body is triaxial compressive stress.

Assuming that the porosity of coal body is n, the vertical stress is [sigma], then the lateral pressure is p = f(n, [sigma]). The vertical stress [sigma] is basically constant, so the change of lateral pressure can be analyzed by using the porosity change indirectly.

As can be seen in Figure 9, the volume of coal body with in the unit length of internal stress field is [V.sub.0] = h[S.sub.1], after compression, the volume becomes

[V.sub.1] = (2h - [epsilon])([S.sub.1] + [xi])/2. (12)

Then the volume difference between the uncompressed volume and compressed volume is given by

[DELTA]V = [V.sub.1] - [V.sub.0] = h[xi] - [epsilon]/2 ([S.sub.1] - [xi]), [epsilon] = [xi] tan [theta]. (13)

Then,

[DELTA]V = [xi](h - [S.sub.1]/2 tan [theta]). (14)

(1) If h - ([S.sub.1]/2) tan [theta] > 0, then [DELTA]V > 0. The volume of coal body increases due to the increase of the porosity. At this moment, the lateral stress decreases at the interface between the external and the internal of stress field and the coal body becomes imbalance, which extends the range of internal stress field.

(2) If h - ([S.sub.1] /2) tan [theta] < 0, then [DELTA]V < 0. The volume of coal body decreases due to the compression reducing the porosity. At this moment, the lateral stress increases at the interface between the external and the internal of stress field and the coal body becomes balance, so the range of internal stress field is basically constant.

(3) If h- ([S.sub.1]/2) tan [theta] = 0,then AV = 0. The volume and the porosity of coal body do not change, so the range of internal stress field is basically constant.

3.4. Solution of the Range of Coal Deformation. The Sierpinski gasket, which has the characteristics of self-similarity, is used to analyse the fractal particle size distribution of real unit system.

In Figure 10, the initial model is an equilateral triangle. Then it is cut into quarters, removing the middle one and keeping the three sides; this is called one generation [18, 19]. Cutting the other three small equilateral triangles into quarters, with removing the middle parts and keeping their sides, is called the second generation. Repeating the operation until infinity, the Sierpinski gasket will be done.

Assuming that the gasket size is limited, the pore size is also limited. The pore size in the kth generation is

[L.sub.k], (k = 0,1,2,...) with the maximum size of [L.sub.c] and minimum size of [L.sub.f]. N stands for the amount of pores whose size is bigger than [L.sub.k]. For the initial model, [L.sub.0] = [L.sub.c]. So the relationship between the amount of pores and the pore size is shown in Table 1.

From Table 1, the amount of pores N is given by

N = 1/2([3.sup.k] - 1), [3.sup.k] = [2.sup.k(ln3/ln2)] = [2.sup.KD] = [L.sup.-D.sub.k][L.sup.D.sub.c]. (15)

Then,

N = 1/2 [L.sup.D.sub.c] ([L.sup.-D.sub.k] - [L.sup.-D.sub.c]). (16)

Equation (16) is the expression for fractal distribution of pore during the breaking process of coal body.

In order to qualitatively analyse the mechanical feature during the compressive process of coal body, we assume that the elements and pore body are generated by using the Sierpinski gasket. The area of initial element is 1, and the relationship between the amount and the size is shown in Table 2.

The area of pore body of coal body is given by

[S.sub.0] = 3k/[2.sub.2k], (k = 1,2,3,...). (17)

Then, the total volume of pore in coal body is given by

[V.sub.pore] = N[S.sub.0] = (1-(3/4)k), (k=1,2,3,...). (18)

With the assumption of the area of deformed coal body of [V.sub.coal] = [S.sub.1]h, the range of coal outburst is given by

[xi] = [S.sub.1] [1-[(3/4).sup.k]], (k=1,2,3,...), (19)

where [S.sub.1] is the width of the deformed coal body with the unit of m and h is the height of the deformed coal body with the unit of m.

4. Conclusion

(1) The spatial form of stope model can be divided into two parts: the outrange of overlying strata spatial structure of stope that consists of no significant move rock beams in the stress arch that does not move significantly and the overlying strata spatial structure that is composed of kinetic rock strata in the breaking arch that affects the ground pressure of stope directly.

(2) According to mechanical balance model of stope, it solves the mechanical parameters of stress field and modifies the expression for the distribution range of internal and external stress field. The range of internal stress field is 1.4 times greater than that of itself calculated by the traditional method.

(3) The solution for the distribution range of abutment pressure can be divided into two stages. When the advancing length of heading face is less than that of the heading face width, the width is expressed by [S.sup.x] [approximately equal to] ([square root of (1 + ((4H - [pi][H.sub.g])/8([k.sub.a] - 1)H))] - 1)[L.sub.0]. When the advancing length of heading face is greater than that of the heading face width, the width is [S.sub.xmax] [less than or equal to] [L.sub.0](4H - [pi]Hg)/8([k.sub.a'] - 1)H. The expression for the range of internal stress field is [S.sub.0] = [square root of (2[C.sub.i][H.sub.g][S.sub.1]/[K.sub.max]H)].

(4) The broken coal body has the characteristic of self-similarity. This paper studies the relationship between block degree of coal body and damage variables and extension mechanism of internal stress field. According to the Sierpinski gasket method, the pore area in coal body is given by [V.sub.pore] = N[S.sub.0] = (1-[(3/4).sup.k]) and the range of coal outburst is given by [xi] = [S.sub.1][1 - [(3/4).sup.k]],(k = 1,2,3,...).

http://dx.doi.org/ 10.1155/2014/714507

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by National Basic Research Program of China under Grant no. 2012CB72310402, National Natural Science Foundation of China under Grant no. 51244010, no. 51304126 and no. 51374046, Science Research Innovative Group of College of Resources and Environmental Engineering of SDUST no. 2012ZHTD06, New Teachers' Fund for Doctor Stations of Ministry of Education under Grant no. 20123718120009, Open Project of State Key Laboratory Breeding Base for Mining Disaster Prevention and Control (Shandong University of Science and Technology) no. MDPC2012ZR01), The Research Fund for Excellent Young and Middle-Aged Scientists of Shandong Province under Grant no. BS2013NJ007, Fok Ying Tung Education Foundation under Grant no. 141046, China Postdoctoral Science Foundation no. 2013M541918, and High Schools' Outstanding Young and Middle-Aged Backbone Teachers of International Cooperation Training Project Funding in Shandong Province.

References

[1] X. X. Miao and M. G. Qian, "Broken feature of key strata and its influence on rock pressure in super-length fully-mechanized coal face," Chinese Journal of Rock Mechanics and Engineering, vol. 22, no. 1, pp. 45-47, 2003.

[2] F. X. Jiang, X. M. Zhang, S. H. Yang, L. Xun, Q. Ma, and H. Wang, "Discussion on overlying strata spatial structures of longwall in coal mine," Chinese Journal of Rock Mechanics and Engineering, vol. 25, no. 5, pp. 979-984, 2006.

[3] H. Kanamori and M. Kikuchi, "The 1992 Nicaragua earthquake: a slow tsunami earthquake associated with subducted sediments," Nature, vol. 361, no. 6414, pp. 714-716, 1993.

[4] R. Mitra and E. C. Westman, "Investigation of the stress imaging in rock samples using numerical modeling and laboratory tomography," International Journal of Geotechnical Engineering, vol. 3, no. 4, pp. 517-525, 2009.

[5] Y. L. Tan, T B. Zhao, and Y. X. Xiao, "Researches on floor stratum fracturing induced by antiprocedure mining underneath close-distance goaf," Journal of Mining Science, vol. 46, no. 3, pp. 250-259, 2010.

[6] C. J. Hou and X. H. Li, "Stability principle of big and small structures of rock surrounding roadway driven along goaf in fully mechanized top coal caving face," Journal of China Coal Society, vol. 26, no. 1, pp. 1-7, 2001.

[7] G. Zhang, Y. Du, Y. Zhang et al., "Analyses and the study of syngas production in dual-gas resources polygene ration," Asian Journal of Chemistry, vol. 25, no. 13, pp. 7287-7289, 2013.

[8] F. X. Jiang, "Viewpoint of spatial structures of overlying strata and its application in coal mine," Journal of Mining & Safety Engineering, vol. 23, no. 1, pp. 30-33, 2006.

[9] K. Yang and G. X. Xie, "Stress field of surrounding rocks of fully mechanized top-coal caving gateway with small pillars," Journal of Mining Safety Engineering, vol. 24, no. 3, pp. 311-315, 2007

[10] Z. Q. Song, Practical Ground Pressure Control, China University of Mining and Technology Press, Xuzhou, China, 1988.

[11] A. Hasson and M. Jweeg, "Soil organic carbon sequestration under pastures in arid region," Nature Environment & Pollution Technology, vol. 12, no. 1, pp. 57-62, 2013.

[12] G. Zhang, Y. Du, Y. Zhang, and Y. Xu, "Desulfurization reaction model and experimental analysis of high sulfur coal under hydrogen atmosphere," Journal of Industrial and Engineering Chemistry, vol. 20, no. 2, pp. 487-493, 2014.

[13] J. L. Sheng and R. G. Zhu, "The fractal evaluation of joint roughness coefficient," Journal of Wuhan University of Science and Technology (Natural Science Edition), vol. 1, no. 23, pp. 1-3, 2000.

[14] B. B. Mandelbrot, The Fractal Geometry of Nature, San Francisco Freeman, San Francisco, Calif, USA, 1982.

[15] Q. H. Qin Huang and J. Angelier, "Fracture spacing and its relation to bed thickness," Geological Magazine, vol. 126, no. 4, pp. 355-362, 1989.

[16] F. L. Ladeira and N. J. Price, "Relationship between fracture spacing and bed thickness," Journal of Structural Geology, vol. 3, no. 2, pp. 179-183, 1981.

[17] B. M. Yu, "Sorne fraetal charaeters of porous media," Fractals, vol. 9, no. 3, pp. 365-372, 2001.

[18] J.-W. Xiao and Y. Yi, "Coupled-adaptive synchronization for Chen chaotic systems with different parameters," Chaos, Solitons and Fractals, vol. 33, no. 3, pp. 908-913, 2007

[19] P. P. Prochazka, "Application of discrete element methods to fracture mechanics of rock bursts," Engineering Fracture Mechanics, vol. 71, no. 4, pp. 601-618, 2004.

Zhijie Wen, (1,2,3) Guanglong Qu, (1) Jinhao Wen, (1) Yongkui Shi, (1) and Chuanyang Jia (1)

(1) State Key Laboratory of Mining Disaster Prevention and Control, Shandong University of Science and Technology, Qingdao, Shandong 266590, China

(2) Research Center of Geotechnical and Structural Engineering, Shandong University, Jinan, Shandong 250061, China

(3) Aalto University, 00076 Aalto, Finland

Correspondence should be addressed to Zhijie Wen; sdust0532@gmail.com

Received 23 March 2014; Accepted 29 April 2014; Published 19 May 2014

Academic Editor: Guojie Zhang

TABLE 1: The void distribution of Sierpinski gasket. k 0 1 2 [L.sub.k] [L.sub.c] [2.sup.-1][L.sub.c] [2.sup.-2][L.sub.c] N 0 1 4 k 3 4 ... [L.sub.k] [2.sup.-3][L.sub.c] [2.sup.-4][L.sub.c] ... N 13 40 ... k k ... [L.sub.k] [2.sup.-5][L.sub.c] ... N ([3.sup.k] -1)/2 ... TABLE 2: The area distribution of Sierpinski gasket. k 0 1 2 [L.sub.k] [L.sub.c] [2.sup.-1][L.sub.c] [2.sup.-2][L.sub.c] [S.sub.pore] 0 1/4 7/16 k 3 4 ... [L.sub.k] [2.sup.-3][L.sub.c] [2.sup.-4][L.sub.c] ... [S.sub.pore] 37/64 91/256 ... k k ... [L.sub.k] [2.sup.-5][L.sub.c] ... [S.sub.pore] 1 - [(3/4).sup.k] ...

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Wen, Zhijie; Qu, Guanglong; Wen, Jinhao; Shi, Yongkui; Jia, Chuanyang |

Publication: | The Scientific World Journal |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 4326 |

Previous Article: | Decision support model for introduction of gamification solution using AHP. |

Next Article: | Acaricidal activity of petroleum ether extract of leaves of tetrastigma leucostaphylum (Dennst.) Alston against rhipicephalus (Boophilus) annulatus. |

Topics: |