Open access peer-reviewed chapter

Numerically and Analytically Forecasting the Coal Burst Using Energy Based Approach Methods

By Faham Tahmasebinia, Chengguo Zhang, Ismet Canbulat, Onur Vardar and Serkan Saydam

Submitted: May 22nd 2017Reviewed: October 23rd 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71879

Downloaded: 311

Abstract

Coal burst is referred to as the violent failure of overstressed coal, which has been recognised as one of the most critical dynamic failures in coal mines. This chapter aims to analytically and numerically evaluate the energy transformation between the different strata and coal layers. An accurate closed-form solution is developed. Due to the complexity of the causes and mechanisms contributing to the coal burst occurrence, 3D finite element modelling as well as discrete element models will be developed to validate the suggested analytical assessments of rock/coal burst occurrence. The energy concept is emphasised in order to improve the understanding of the underlying mechanisms of coal burst. Only with enhanced understanding of the driving mechanisms, a reliable coal burst risk assessment can be achieved.

Keywords

  • analytical modelling
  • numerical modelling
  • released energy
  • coal burst
  • failure mechanism

1. Introduction

One of the critical engineering problems faced by the coal mining industry is coal burst. It is caused by a dynamic release of energy within the overstressed rock mass/coal during the mining process. It occurs under the effects of complex environments of geology, stress and mining conditions. It has been recognised that the unstable releases of potential energy of the rock around the excavations, mainly in the form of kinetic energy, contributes to the coal burst occurrence. Interactions between the coal and rock interface, as well as the confinement, can completely determine the failure mode and the ultimate bearing capacity of coal pillars, influencing the amount of stored energy within a pillar. Many authors define rock/coal burst as a sudden, rapid rupture of the rock mass with a violent, explosive release of elastic/strain energy from the surface of an excavation, which is generally associated with a seismic event and produces rock particle ejections [1, 2, 3, 4, 5]. The coal burst source is the mechanism that triggers or induces the damage mechanism visible on the excavation surface. The coal burst source is generally associated with a seismic event that can be performed at a wide range of local magnitudes, normally ranging from undetectable up to 5 [6]. Indeed, mining-induced seismicity can reach moderate values of ground velocity and acceleration, and in some cases its effects on the surface can be compared with low-intensity earthquakes [7]. The mechanism that produces the seismic event is a sudden release of the strain energy that has been stored above a critical level within the rock/coal mass. Some portion of this energy is demolished by crack development, and the rest of the energy is converted into the kinetic energy [8, 9]. When the energy source is located near the roadway, the released energy may lead to coal fragmentation. At the place of the source of the energy, where it is located in a plane of weakness inside the coal mass, the released energy provokes shear displacement along the plane, which in revolve generate vibrations that persuade coal ejections when they are situated near the excavation boundaries [7]. Tarasov and Randolph [6] have explained a number of special and inconsistent behaviours of hard rock at the significant depth that are directly related to rock failure mechanisms in deep excavations. They determined that the procedures of the shear failure, with respect to the significant low friction, can be classified as the main reason to release energy. Based on the suggested frictionless mechanism, the level of the brittleness of the confined rock/coal masses might be increased under high stress conditions. This may result in reducing the overall ductility which would in line with the abrupt fracture failure. Under an energy-balance approach, the methods to predict coal burst risk are based on energy indexes such as energy release rate (ERR) [8, 9, 10], energy storage rate (ESR), strain energy storage index (WET) [11], potential energy of elastic strain (PES) or strain energy density (SED) (i.e., the elastic strain energy in a unit volume of the coal mass, which can be computed by the uni-axial compressive strength of the coal and the relevant unloading tangential modulus), and burst potential index (BPI). A combination of both analytical as well as numerical methods, where they can comprehensively evaluate the structural performance of the mine scale, would be broadly addressed in the current research. Thus, the following aims explicitly will be addressed.

  1. Develop a full 2D and 3D finite element as well as discrete element models to compute the inducted energies in a single pillar with different high to width ratios. In this approach, different loading conditions varying from the static, quasi-static as well as dynamic loading will be exclusively examined.

  2. Considering the effect of the energy transformations between the rock/coal layers due to the different contact/joint properties.

  3. Suggest empirical equations to predict the amount of the released strain energy due to the mining activities.

The main novelty of this research is to simulate the effect of the failure and post-failure of the engaged material as well as joint/contact properties on the energy transformation.

2. Numerical modelling strategy

Numerical simulations can be considered as an individual tool to predict possible failure modes and the actual capacity of the mining setting. It is mostly useful to undertake parametric and sensitivity analyses to gain better understanding the nature and level of indecision, or remaining hazard, associated with design process.

First, a finite element model is developed by taking advantage from the commercial software package ABAQUS/Explicit. All the geotechnical components, including the rock and coal, were modelled by the eight-node linear brick element (C3D8R) available in the ABAQUS library. Element C3D8R relies on reducing integration and hourglass control. The assigned meshes were established by using the structured technique available in ABAQUS. The solution to the nonlinear problem was sought using the explicit dynamic analysis procedure available in ABAQUS. In the current study, Figure 1 presents a quarter of a single pillar.

Figure 1.

Illustration of a typical single pillar model using ABAQUS/Explicit.

Thus, by taking advantage from the symmetrical boundary conditions, a finer mesh was assigned to the model. Finding the right input material properties would be a significant assumption, which has not been appropriately studied in the available literature. Modelling of mechanical behaviour of the coal under both compression and shear stresses would be very complicated, since there are no articulated reports which might be concerned with the uni-axial and tri-axial behaviour of coal under both static and dynamic loading conditions. According to the elastic analysis, the stress analysis and energy computations were organised in line with the linear relationship between the stress and the strain in coal and overburden properties. The peak and post-peak behaviour of coal and surrounding rock masses will be ignored. Therefore, in the current literature, the computed stress, strain and kinetic energy have been noticeably overestimated. At the second stage, a combination of the 2D and 3D discrete element models using UDEC and 3DEC was developed. Figure 2 illustrates the pillar model incorporating half of coal, roof and floor along the symmetrical centre-line of the pillar. The height of the roof and floor was 20 m and the mining height was fixed at 3 m, while the pillar widths varied in order to simulate the pillars with width to height (w/h) ratios from 1 to 5.

Figure 2.

Geometry and zoning of coal pillar model using UDEC.

Figure 3.

Geometry and zoning of coal pillar model using 3DEC.

A Mohr-Coulomb (MC) material that presents a constant strength after failure and a Mohr-Coulomb strain-softening material that can reach the peak strength and then decrease to a residual strength have been considered. A quasi-static loading condition as a velocity was applied on the top and bottom of the model. The applied velocity was started with a very small, constant velocity to represent a relative loading system to promote a model of a coal failure that progresses slowly. Simulating a proper loading/displacement condition is significantly crucial, specifically, gaining a sound understanding of the structural reaction of a single coal sample under dynamic or quasi-static loading conditions. Consideration was also given to defining a joint interface. A Coulomb Slip (CS) joint interface property, where it is represented by displacement softening parameters, was taken into account to simulate the interface properties between the different joints.

The uniform zone size of 0.1 m was applied to the coal, and a smooth variation of zoning from the coal to the boundaries was used for roof and floor with appropriate aspect ratios to avoid numerical instability. Roller boundaries were applied along the side of the roof and floor, the bottom of the floor and the vertical line. The same trend was applied to develop the three-dimensional discrete element using 3DEC (see Figure 3).

3. Analytical approach

An analytical method is developed to evaluate shear stress and strain distributions between the engaged surfaces throughout different joint layers by considering the beam theory method in different directions with respect to the different planes, where it can independently calculate shear forces between the different layers and shear strain as well as the curvature distribution along the different layers that have been extracted. The main concept to derive the following equations was extracted from [12, 13]. The cross-sectional analysis is based on the assumption of the Euler-Bernoulli beam model. The strain distribution across the section can be calculated by ε=εry×κ, where εris the strain at the reference point (which can be determined at any point), yis the distance between the selected point and location of the neutral axis of the cross-section and κis the curvature across the section in different strata layers. A vector can be introduced by KDwhich will be included in the internal action N(axial forces) and M(internal moment). External loads, which might be due to the effect of the self-weight of the strata layers as well as the possible applied forces due to the vertical or horizontal displacement in the different layers, can induce the external axial force Neand external moment Me. The relationship between the internal and external actions can be presented by:

ε=εrκE1
rε=NME2
re=NeMeE3
rε=reThis is the vector for strainE4

By considering the nonlinear interactions, the presented equations can be re-written by:

rεi+1=rεi+rtεi×Δεi=reE5
rtεi×Δεi=rRiE6
rRi=rerεiE7
Nεriκiεr×Δεri+Nεriκiκ×Δκ=NRiE8
NRi=NeNεriκiE9
MRi=MeMεriκiE10

All the equations can be re-presented in matrix format:

rtεi=NεriκiεrNεriκiκMεriκiεrMεriκiκE11
Δεi=ΔεriΔκiChanging strain and curvatureE12
rRi=NRiMRiE13

The partial derivatives of Nand Mwith respect to εrand κcan be re-arranged in a more practical form, recalling the definitions of internal actions as:

Nεriκiεr=σεrdAE14
Nεriκiκ=σκdAE15
Mεriκiεr=yσεrdAE16
Mεriκiκ=yσκdAE17

where the values of the stress depend on the constitutive models adopted for the materials and on the magnitude of the strain

Nεriκiεr=σεrdA=σε×εεrdA=σε×εry×κεrdA=σεdAE18
Nεriκiκ=σκdA=σε×εry×κκdA=y×σεdAE19
Mεriκiεr=y×σεrdA=y×σε×εry×κεrdA=y×σεdAE20
Mεriκiκ=y×σεrdA=y×σε×εry×κκdA=y2×σεdAE21
σ=E×εforεεpelastic strainE22
σ=fpforε>εpplastic strainE23
σε=E×εε=Eforεεpelastic strainE24
σε=fpε=0forε>εpplastic strainE25
Nεriκi=σdA=j=1njσyjεriκi×AjE26
Mεriκi=yσdA=j=1njyj×σyjεriκi×AjE27
Nεriκiεr=σεdA=j=1njσyjεriκiε×AjE28
Nεriκiκ=y×σεrdA=j=1njyj×σyjεriκiε×AjE29
Mεriκiεr=y×σεrdA=j=1njyj×σyjεriκiε×AjE30
Mεriκiκ=y2×σεdA=j=1njyj2×σyjεriκiε×AjE31
rxde=NxdeMxde=AσxydedAAyσxydedA=j=1njσxydeAjj=1njyjσxydeAjE32
I=abfxdx=ba2×11fa+b2+ba2×x¯=ba2×k=1nGwk×fa+b2+ba2×xk¯E33
kexde=LBTxrxdedx=L2k=1nGwkBTxkrxkdeE34
qe=LNeTxpxdx=L2k=1nGwkNeTxkpxkxk=L2xk¯+1E35
uxvx=Nu1x00Nv1x0Nv2xNu2x0Nu3x00Nv3x0Nv4x×uLvLθLuMuRvRθR=NexdeE36
Nu1x=13xL+2x2L2E37
Nu2x=4xL+4x2L2E38
Nu3x=xL+2x2L2E39
Nv1x=13x2L2+2x2L3E40
Nv2x=x2x2L+x3L2E41
Nv3x=3x2L22x3L3E42
Nv4x=x2L+x3L2E43
ux=Nu1xuL+Nu2xuM+Nu3xuRvx=Nv1xvL+Nv2xθL+Nv3xvR+Nv4xθRE44
Bx=Nu1x00Nv1x0Nv2xNu2x0Nu3x00Nv3x0Nv4xE45
Nu1x=3L+4xL2E46
Nu2x=4L8xL2E47
Nu3x=1L+4xL2E48
Nv1x=12xL36L2E49
Nv2x=6xL24LE50
Nv3x=6L212xL3E51
Nv4x=6xL22LE52
kexdei=L2k=1nGwkBTxkrixkdei=L2k=1nGwk×Nu1xk00Nv1xk0Nv2xkNu2xk0Nu3xk00Nv3xk0Nv4xk×NxkdeiMxkdeiE53
qe=L2k=1nGwk×Nu1xk00Nv1xk0Nv2xkNu2xk0Nu3xk00Nv3xk0Nv4xk×nxkwxkE54
Nxkdei=j=1njσxkyideiAjE55
Mxkdei=j=1njyjσxkyideiAjE56

Thus, by calculating stress and strain at the different points in the different layers of the overburden, the internal axial forces as well as internal moments can be calculated. It was assumed that the strain energy can be calculated by:

A=12×σxx×εxx+σyy×εyy+σzz×εzz+σxy×εxy+σxz×εxz+σyz×εyzdxdydzE57

where σxx×εxx,,σyz×εyzcan be calculated, according to the principal of the virtual work and virtual deformation δA=δR1+δR2, when the induced stresses and strains cannot be directly extracted from the simulated model.

4. Energy calculation based on the numerical approach

According to Xie et al. [14], the coal burst proneness of a coal can be determined by the coal burst proneness assessments. Special attentions were devoted by the number of researchers to develop coal burst proneness indexes, which are broadly utilised, such as elastic energy, impact energy, dynamic failure time as well as elastic deformation and stiffness ratio indexes. The elastic energy index WET is defined as the ratio of the elastic strain energy and the strain energy dissipation at point E (75–85% of the peak strength). As shown in Figure 4, the ratio of the area SEAC (between the unloaded line EA and the strain axis) and the area SEOA (between the load and unload line) is the elastic energy index

WET=SEACSEOAE58

Figure 6.

Failure mode of a single pillar with the different w/h ratios.

The impact energy index KE is defined as the ratio of the pre-peak area and the post-peak area, KE namely, the ratio of energy in the pre-peak stage and the energy released in the post-peak stage.

5. Energy calculation based on the analytical approach

According to Xie et al. [14], dissipated and released energy can play a significant role which may result in coal deformation and failure. Based on the failure mechanism, the fracture procedure of a coal mass might be started from a partial fracture which would be followed by local damage. This procedure will be finally resulting in collapsing the mining structures. The failure process is thermodynamically permanent, which includes released and dissipated energy. The dissipated energy can cause damage as well as a permanent deformation of the coal mass, which is followed by weakening of strength. A sudden release of the strain energy may lead to a catastrophic failure, which clearly indicates a certain condition where the coal mass collapses. The released and dissipated energy from the coal mass, individually, plays an essential role in the relevant sudden failure, which would be one of the major requirements to investigate the procedure of the deformation and failure of a coal mass. Figure 5 is a typical compression curve of coal under a constant displacement.

Figure 4.

A typical stress-strain curve.

Figure 5 explicitly demonstrates the calculation of the dissipated, released and residual energies. With respect to the constant development of the inner micro-defects, energy dissipation is an indispensable characteristic of the deformation and failure of the coal mass. The evolution declines the strength of the coal, which may result in total failure. In this content, the dissipated energy is directly concerned with the damage as well as mitigating strength of the coal.

ud1=0ε1σa=0ε1ε1iε1i1×σai+σai14Dissipated energy before peakE59
ud2=ε3ε4σc=ε3ε4εjεj1×σcj+σcj14Dissipated energy after peakE60
ur=ε1ε2σm=ε1ε2εkεk1×σmk+σmk14Released elastic energyE61
ued=ε2ε3σp=ε2ε3εlεl1×σpl+σpl14Residual elastic energyE62

Tables 1 and 2 presents a comparison between the different elastic and post-failure energy components using UDEC and 3DEC output results as well as semi-close form solutions. As it can be found, there is a good agreement between the suggested semi-analytical methods as well as the calculated key energy components which were extracted from the UDEC and 3DEC output results.

Elastic strain energy (kJ/m3)Dissipated elastic strain energy (kJ/m3)The amount of the energy in the pre-peak stage (kJ/m3)The energy released in the post-peak stage (kJ/m3)
Ratio (w/h)UDECAnalyticUDECAnalyticUDECAnalyticUDECAnalytic
11.561.630.780.773.613.7312.4710.59
1.51.921.350.80.787.897.4412.4811.03
22.06212.0040.9910.9210.229.8318.3117.09
2.54.704.821.161.1114.4713.4321.1723.14
311.1310.582.512.4135.82532.8711.7310.6
414.7213.274.074.6060.2655.0056.3470.02
516.6316.245.3345.3775.8373.6791.1984.04

Table 1.

A comparison between the different energy components (UDEC and the analytical solution).

Elastic strain energy (kJ/m3)Dissipated elastic strain energy (kJ/m3)The amount of the energy in the pre-peak stage (kJ/m3)The energy released in the post-peak stage (kJ/m3)
Ratio (w/h)3DECAnalytic3DECAnalytic3DECAnalytic3DECAnalytic
12.582.651.7671.554.914.7714.8714.77
1.52.942.371.880.788.888.4616.5615.23
24.244.012.8912.9212.5511.9821.2320.14
2.56.726.843.183.1115.3715.4124.3523.99
313.1512.625.515.4437.1536.8728.4527.68
416.7616.277.077.2060.2662.3359.1157.88
519.8319.288.3348.2276.2275.1296.5492.66

Table 2.

A comparison between the different energy components (3DEC and the analytical solution).

6. Progress of the failure in different pillar ratios

Different loading conditions varying from the quasi-static to dynamic loading has been applied to the coal pillar with the different width to height (w/h) ratio to determine the pillar capability as well as the possible observed failure modes. A strain-based criterion, as a major failure criterion, was implemented in the ABAQUS/Explicit to predict of the cracking path due to the different applied loadings as well as different pillar geometrical properties. A quarter coal pillar model based on the symmetrical boundary conditions with respect to the different width by height (w/h) ratios of 1–10 were developed. It was observed that when the w/h ratios are less than 4, the failure mode of pillar can be either a double or a single diagonal shear failure in which the trajectory cracking starts from the edges and progresses towards the centre of the pillar. While the w/h ratios are greater than 4, the obtained possible failure mode would be a combination of the shear and compression failure modes. Thus, the trajectory of the cracking due to the pure compression failure would be propagated from the centre to the corners where a pillar gradually starts towards fully squashed (see Figure 6).

Figure 5.

Analytically calculation of dissipated energy and released energy (Xie et al. [14]).

7. Remarkable conclusions

Analytical method is an important part of coal burst evaluation and forecasting. Analytical forecasting methods, either alone or combined with numerical simulations, can be used to estimate both in situ stress and induced stress, which leads to the prediction of failure-prone areas and calculation of critical values of the energies. The behaviour of a single pillar under different applied loads ranging from the quasi-static towards the dynamic loading conditions was simulated using commercial finite element package ABAQUS/Explicit. A strain-based failure condition was evaluated to determine the failure modes in a single pillar by respecting to the different w/h ratios. The observed numerical failure modes can be classified by shear and compression failures as well as a combination of both shear and compression were comprehensively illustrated. The released energy or residual energy is either transferred into kinetic energy or dissipated energy in non-elastic behaviour such as joint shear and fracturing. The unstable release of potential energy of the coal around the excavations, mainly in the form of kinetic energy, causes coal burst.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Faham Tahmasebinia, Chengguo Zhang, Ismet Canbulat, Onur Vardar and Serkan Saydam (December 20th 2017). Numerically and Analytically Forecasting the Coal Burst Using Energy Based Approach Methods, Finite Element Method - Simulation, Numerical Analysis and Solution Techniques, Răzvan Păcurar, IntechOpen, DOI: 10.5772/intechopen.71879. Available from:

chapter statistics

311total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Linear Thermo-Poroelasticity and Geomechanics

By Horacio Florez

Related Book

First chapter

Modeling the Physical Phenomena Involved by Laser Beam – Substance Interaction

By Marian Pearsica, Stefan Nedelcu, Cristian-George Constantinescu, Constantin Strimbu, Marius Benta and Catalin Mihai

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us