Next Article in Journal
On Hybrid Energy Utilization in Wireless Sensor Networks
Next Article in Special Issue
Optimizing the Structure of the Straight Cone Nozzle and the Parameters of Borehole Hydraulic Mining for Huadian Oil Shale Based on Experimental Research
Previous Article in Journal
Energy Management and Control of Plug-In Hybrid Electric Vehicle Charging Stations in a Grid-Connected Hybrid Power System
Previous Article in Special Issue
Heat Transfer Characteristics and Prediction Model of Supercritical Carbon Dioxide (SC-CO2) in a Vertical Tube
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Cascade Disaster Caused by Geological and Coupled Hydro-Mechanical Factors—Water Inrush Mechanism from Karst Collapse Column under Confining Pressure

1
State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, Jiangsu, China
2
Shanxi China Coal Huajin Energy Limited Liability Company, Hejin 043300, Shanxi, China
3
School of Resources and Geosciences, China University of Mining and Technology, Xuzhou 221116, Jiangsu, China
*
Authors to whom correspondence should be addressed.
Energies 2017, 10(12), 1938; https://doi.org/10.3390/en10121938
Submission received: 23 September 2017 / Revised: 7 November 2017 / Accepted: 17 November 2017 / Published: 23 November 2017
(This article belongs to the Special Issue Mathematical and Computational Modeling in Geothermal Engineering)

Abstract

:
The water inrush from karst collapse column (KCC) is a cascading, vicious cycle disaster caused by geological and mining activities, that can cause serious casualties and property losses. The key to preventing this risk is to study the mechanism of water inrush under confining pressure. Aiming at the investigationg the characteristics of the KCC named X1 in Chensilou mine, a series of methods, including connectivity experiments, water pressure monitoring tests in two side-walls, and numerical simulations based on plastic damage-seepage (PD-S) theory have been developed. The methods are used to test the security of the 2519 mining area, the damage thickness, pore water pressure, and seepage vector in the X1. The results indicate that the X1 has a certain water blocking capacity. In addition, with the decrease of confining pressure and increase of shear stress, deviatoric stress could cause the increase of permeability, the reduction of strength, and the reduction of pore water pressure in KCC. Therefore the increased effective stress in the rock will force the rock to become more fractured. Conversely, the broken rock could cause the change of stress, and further initiate new plastic strains, damage and pore water pressure until a new equilibrium is reached. This cascading water inrush mechanism will contribute to the exploitation of deep coal resources in complex geological and hydrogeological conditions.

1. Introduction

Mine water inrush and consequent secondary disasters have caused huge losses in many countries, including USA [1], Russia [2], Hungary [3], Germany [4], Australia [5] and China [6]. As a unique geological phenomenon in north China, karst collapse column (KCC), which is caused by the karst subsidence that occurs in Ordovician limestone caves, forms cone, barrel, inclined tower or irregular shapes [7], connects coal seams and aquifers (Figure 1), and may cause serious groundwater inrush disasters under the influence of the confined water pressure and the ground pressure, including axial pressure and confining pressure [8].
Between the 1950s to 1990s, people mainly study the formation process, distribution law and exploration method of KCC by geological and hydrogeological theories [9,10,11]. In recent years, with the increase of mining depth, serious water inrush accidents from KCC have happened in many mines and attracted scholars' attention [12,13]. Different theories and methods have been used to study the water inrush mechanism of KCC. Li et al. [14] and Zhang et al. [15] assumed that the rock mass is an elastic material, and the stresses in the internal and external of KCC under the mining influence were calculated by solid mechanics theory. Yin and Zhang [16] and Ou Subei et al. [17] calculated the distribution of the stress and plastic zone of KCC under the influence of water pressure by numerical simulation method. According to the formation process of the KCC, Ma et al. [4] believed that the fluid-solid coupling theory can better reflect the geological condition of KCC, and the distribution of plastic zone and pore water pressure under the mining influence were calculated by numerical simulation. Zhang [18] simplified the floor as a clamped-clamped beam and calculated the tensile failure area of the beam under the influence of confined water pressure by material mechanics theory, the results were verified by water injection test. Zhu et al. [19] used the “strain method” to test the strain of floor rock mass during mining advance, monitoring the deformation of rock strata in different depths.
The axial stress or shear stress, deformation and range of plastic zone under the confined water pressure of Ordovician limestone are the focus of the above literature. However, they focus less on the effects of plastic strain, damage, seepage, especially the effect of confining pressure on the water inrush mechanism under complex hydrogeological conditions, including a KCC and multiple limestone aquifers in the Taiyuan formation. During the deep mining process, the variation of confining pressure could cause large plastic deformation of the KCC [20], and lead to the formation of damage zone. In this case, the pore water pressure decrease, the permeability coefficient and seepage velocity increase, and eventually, a cascade disaster is formed in the KCC.
In this paper, the 2519 mining face of Chensilou mine has been chosen as the study area, where the floor contains a KCC and three high confined aquifers, meanwhile the confining pressure is 19.3 MPa. In order to simulate the cascade effect of such complex geology and hydrogeology on the KCC, a new plastic damage-seepage coupled model that conformed to the characteristics of mechanical and hydraulics response of rock under confining pressure was established. The damage distribution and seepage vector in KCC during mining advance were obtained by programming. Instead of “strain method” and “single side-wall water injection test”, a new in-situ test, that is, water pressure monitoring test in two side-walls has been developed to test the damage thickness in the KCC’s side-walls where it has a high probability of water inrush. By these methods, we found that the damage thickness at the left side-wall was slightly larger than that at the right side-wall, resulting in the instantaneous seepage velocity to the maximum. But the total water inflow at the right side-wall is greater due to the reduced confining pressure. This conclusion provides a basis for coal mine safety production.

2. Study Area

Compared to the late Paleozoic coalfields exploited in other parts of the world, such as the Appalachian coalfield in the United States, the Yorkshire coalfield in Great Britain, the Ruhr coalfield in Germany and the Donetz coalfield in Ukraine, there is a great risk of water inrush in the North China coalfields due to the presence of karst collapse column (KCC) or other complex geological structures [21]. Furthermore, compared to other coal mine area in North China, geological conditions in the Chensilou mine are even worse: 2519 mining face contains a KCC filled with water, the KCC directly connected three high pressure aquifer and coal seam and there is a synclinal basin in 2519 mining area. Therefore, the study of the damage and seepage law of the KCC under the mining influence can provide a reference for other KCCs in North China.

2.1. Overview of the Chensilou Mine

Chensilou coal mine is located northeast of Yongcheng City in Henan Province (Figure 2). In the coal seam floor, the stratigraphic succession consists of two obvious rock formations from bottom to top: Ordovician (O), containing a thick limestone aquifer, and Carboniferous (C), containing a coal bed, two limestone aquifers and several sandstones and mudstones. In the horizontal direction, there is a synclinal basin to the southwest of the mine.

2.2. Engineering Situations at 2519 Mining Face

We selected the 2519 mining face as the study area for the following three reasons: (1) as the coal in the southeastern mine has been mined, coal will now be mined in the southwest (Figure 2); (2) water inrush accidents have occurred more than 20 times since 1997 (Figure 2). As the mining face has advanced towards the syncline basin, geological and hydrogeological conditions have become even worse. In the east of the south part, the mining depth is 390 m, the water pressure in the three limestone aquifers are 2.6 MPa, 1.5 MPa and 1.3 MPa from bottom to top, the distance between the coal seam and the limestone 1 is 60 m, and the confining pressure is 6.6 MPa. However, in the 2519 mining face, the mining depth increases to 890 m, the water pressures in the three limestone aquifers are as high as 7.0 MPa, 3.3 MPa and 3.0 MPa from bottom to top, and the confining pressure is 28.6 MPa. The distance between the coal seam and the upper limestone aquifer is only 24 m (Figure 3), while the sandstone in floor is more broken due to tectonic stress in the syncline structure; (3) more serious is the fact that the 2519 mining area has a KCC named X1. According to the drilling data, the coring recovery rate (the ratio of the total length of intact rock and broken rock to drilling length) is only 52%, 48% and 23% in the upper, middle and lower part of X1, respectively (Figure 3). Meanwhile, the X1 connects the three limestone aquifers and coal seam, making this weak geological structure have a high risk of water inrush.

3. In-Situ Test

Rock damage caused by mining activities will cause a significant change in the permeability of KCC. Thus the main purpose of in-situ test is to test the damage characteristics of KCC under the mining influence. According to the geological and hydrogeological conditions of the 2519 mining face as described in Section 2.2, it might be dangerous to directly test the damage thickness of X1 during mining advance. Therefore, a connectivity experiment was carried out firstly to test the initial permeability of the KCC.

3.1. Connectivity Experiment

The principle and the process of the connectivity test is a follows: three boreholes (i.e., P01, P02 and P03) were drilled into the middle of the KCC with the vertical depths of 28 m, 33 m and 38 m respectively (Figure 3 and Figure 4). We injected high pressure water continuously into water injection hole (P01 or P03) and observed the water flow from the observation hole (P02), the permeability of rock mass between the two holes can be determined [7]. A plugging operation is required after the test to avoid interfering with the results of the “water pressure monitoring test”. The results of water flow in the X1 at different water pressure conditions are shown in Table 1 and Table 2.
Table 1 and Table 2 indicate that the permeability of the X1 depends greatly on the water pressure under initial conditions. Although the water pressure in P01 continues to increase, however, the water flow in P02 increased slowly when water injection pressure was less than 5 MPa. Once the water pressure exceeded 5 MPa, water inrush will occur in P02. The connectivity experiment between P03 and P02 has similar change laws, while the water inrush phenomenon occurred when water pressure in P03 reached 5.5 MPa. Compared the two groups of connectivity test, the final water inflow in P02 (P01 test) was more than that of the P02 (P03 test), and the water pressure in P01 hole is less than P03.
The above results indicate that the permeability between the P03 and P02 is higher than that between the P01 and P02 in the initial state. However, the ultimate permeability between the P03 and P02 is less than that between the P01 and P02 due to the confining pressure in the former section higher than that in the latter section. In other words, the water blocking capacity of X1 increases with the increase of confining pressure.
In addition, as shown in Table 1 and Table 2, the rock pillar with 5 m long in the X1 cannot withstand a water pressure of more than 5 MPa, but the X1 is still safe under the action of three aquifers after analysis. The main reasons are as follows: (1) the water pressure is only 3 MPa in limestone 1, and the water pressure in limestone 1 will not exceed 5 MPa under the mining influence (Section 5.1). The water pressure in limestone 3 is 7 MPa, but the distance between the coal seam and the limestone 3 is 82 m (Figure 3). (2) According to the formation process of the KCC and the coring recovery rate (Figure 3), we know that the permeability of the X1 decreases from bottom to top. This internal structure will gradually consume the energy of confined water, and contributes to the water barrier performance of KCC. In a word, the X1 has a certain water blocking capacity, and it is estimated that the water inflow will not exceed the maximum drainage capacity (1100 m3/h) of the mine when mining face pass through the X1 (Section 5.1). Besides that, there is not much water in fractures in limestone and KCC due to the effect of high confining pressure. However, a certain amount of groundwater is contained in the dissolved pores. And fortunately, the amount and size of pores is small. The method of mining KCC can avoid tunnel re-excavation, pipelines and hydraulic supports re-arrangement, greatly reduce production costs. In summary, we judge that mining X1 is safe, and can improve production efficiency.

3.2. Water Pressure Monitoring Test in Two Side-Walls

The channel for groundwater inrush in the KCC is generally located in the side-wall [12]. In the traditional water injection test [4], engineers inject water into several boreholes at the middle of the KCC in a short time, and test the water flow loss of each borehole. However, for the X1 with certain permeability, the water flow in boreholes will interfere with each other during the water injection test. Considering this drawback, and in order to evaluate the damage thickness at different side-walls, we have developed a new in-situ experiment, that is, water pressure monitoring test in two side-walls.
The basic principle of the test is simple, the more broken the rock, the shorter the time of water pressure relief. We can use this “relief time” to determine the damage thickness of the KCC in two side-walls. The test process is as follows: the drilling rig was placed in the 2521 materials haulage roadway, 950 m away from the starting cut (Figure 4). We drilled six boreholes at the two side-walls of the X1, retained the 20 cm water injection section (P1i and P2i in Figure 3 and Figure 4, i = 1, 2, 3, represented vertical depths of 28 m, 33 m and 38 m, respectively) at the end of each borehole. High pressure water with 3.8 MPa was injected into a hole, and the relief time would be recorded. If the water pressure dissipated, the same steps would be performed in the next hole. The relief time were started to record if the distance between the mining face and the central of X1 is 50 m, and no longer recorded if the distance is over 30 m. This test lasted a month, and the results shown in Figure 5.
The data indicate that when the distance between the mining face and the P11 (left side of the X1) was more than 15 m, the water pressure would last about 15–18 min. With mining advance, there is a downward trend of the relief time. As the face passes through the X1, damage appeared at the left side-wall and quickly increased to 28–33 m, and the relief time decreased from 14.08 to 1.01 min (P11), which meant that the floor damage degree and range continued to increase; however, for P13, the relief time is over 3.67 min, which means that the maximum damage thickness does not exceed 38 m or the damage at the depth of 38 m is slight. The relief time of P11 and P12 achieved a stable value after 0.95 min and 1.21 min, respectively, until the mining face passes the center of X1 at 15 m. The comparison of test data from boreholes P2i and P1i showed that the relief time of the two have similar change laws. However, the time in P13 is less than that of P23, indicating that the damage thickness at the left side-wall is greater than that of the right side-wall, and the maximum damage depth is about 38 m.

4. Constitutive Model and Verification

Through the above in-situ test, it is possible to accurately obtain the damage thickness at some points during the mining process. However, numerical simulation is needed to obtain the damage thickness as well as stress, damage, seepage, pore water pressure at each point of the X1. The water inrush mechanism is easily inferred based on the above data.

4.1. Plastic Damage-Seepage Coupled Model

The accuracy of the numerical simulation results is closely related to the accuracy of the constitutive model. The complex hydrogeological conditions (Section 2.2) indicate the rationality of PD-S model to explain the mechanism of water inrush from KCC. The traditional constitutive models for studying KCC, such as plastic-seepage theory [4,17], and elastic damage-seepage theory [19], are not applicable to X1. The main reasons are as follows: (1) for the deep rock mass, the mining influence will produce plastic deformation, so the theory of elastic damage [19] is not appropriate; (2) although the plastic-seepage theory was used in some literature [4,17], they ignored the effect of confining pressure on rock plastic deformation, and the effect of plastic deformation on seepage. In addition, they considered that the permeability coefficient is a function of porosity, which means that it can only be used for the elastic stage of rock.
The three basic equations used in this model are identical to the literature [22], including the equilibrium equation (Equation (1)), the geometric equation (Equation (2)), and the continuous equation (Equation (3)):
σ t F = 0
ε = ( u + u ) / 2
q = ε v t + n K f p w t
where there are eight model parameters: F is the body force vector; u is the displacement vector; q is the seepage velocity vector; t is the time; εv is the volumetric strain; ε is the strain tensor; σt is the total stress tensor; and pw is the pore water pressure. There are two fitting parameters: n is porosity, it can be obtained by the water absorption under the vacuum experimental conditions; Kf is the bulk modulus of water.
However, the remaining equations in the PD-S model, including constitutive equation, strength criteria, plastic potential functions and permeability coefficient equation, are more consistent with the response of deep rock under high confining pressure.
According to the principle of damage mechanics [7], the Boit effective stress theory [8] can be expressed by σt, pw, ε and the effective stress σe, and the relation is shown in Equation (4):
σ t = σ e + α p w = ( 1 d ) D 0 ( ε ε p ) + α p w
where there are two fitting parameters: α is the effective-stress coefficient, determines the effect of pore pressure on rock and can be defined as α = 1 − K’/Ks, Ks is the effective bulk modulus of the solid constituent, and K’ is the drained bulk modulus of the porous medium [14]; D0 is the undamaged elastic stiffness matrix, we can get it by elasticity theory [23]. The remaining parameters belong to the model parameters: εp is the plastic strain vector; d is the damage value, which is the function of model parameter γ ^ p eq , and:
d = ( e γ ^ p eq / a 1 ) / ( e 1 / a 1 )
where γ ^ p eq is a function of confining pressure pc, and γ ^ p eq = d γ p eq / [ A 1 ( p c σ c ) + A 2 ] , its range is 0–1. γ p eq is the equivalent shear plastic strain; σc is uniaxial compressive strength; a, A1 and A2 are constant, and can be obtained as follows: through the cyclic loading and unloading experiment, we can get the plastic deformation and the elastic modulus during each cycle, so that the curves of d γ p eq and γ p eq p c / σ c can be obtained, and then fitting the two curves by linear functions, three material parameters can be obtained.
Based on the modified Mohr-Coulomb yield criterion (Equation (6)), using the non-associated flow rule expressed by Equation (7) [4], the internal variable of εp in the program (Figure 6) can be obtained:
F = ( σ n e tan φ + c ) + d ( p w tan φ c ) τ n e = 0
where, σ n e and τ n e are the normal and tangential effective stresses on the shear plane; and the strength parameters φ’ and c’ (fitting parameters) are the internal friction angle and cohesive force, respectively, both of them are functions of γ ^ p eq , and the relation between γ ^ p eq and φ’ and c’ can be obtained by the following method: (a) by cyclic loading and unloading test, the deviatoric stress and plastic strain under different confining pressure can be obtained. (b) Take equal interval value of γ ^ p eq , the deviatoric stress is linearly interpolated according to the corresponding γ ^ p eq , thus the deviatoric stress in the same γ ^ p eq under different confining pressures can be obtained. (c) according to the above data, the values of φ’ and c’ in Mohr-Coulomb criterion ( σ 1 = 2 c cos φ 1 sin φ + 1 + sin φ 1 sin φ σ 3 ) can be obtained.
The plastic potential function expressed by the effective hydrostatic stress pe and the effective deviatoric stress qe (model parameters) can better reflect the geometric characteristics of plastic potential surface (Figure 6), the expression is shown in Equation (7):
G = ( δ c 0 tan ψ ) 2 + ( R q e ) 2 p e tan ψ
where the parameters are fitting parameters except for pe and qe. δ is eccentricity, which defines the rate at which the function approaches the asymptote, δ = 0.1 [24]; c 0 is the cohesion without plastic deformation; ψ is the dilation angle measured in the peqe plane at high confining pressure; c 0 and ψ can be obtained by cyclic loading and unloading test. R controls the shape of the plastic potential surface in the deflection plane [24], and:
R = R c [ 4 ( 1 e 2 ) cos 2 θ + ( 2 e 1 ) 2 ] 2 ( 1 e 2 ) cos θ + ( 2 e 1 ) 4 ( 1 e 2 ) cos 2 θ + 5 e 2 4 e
where Rc is the polar diameter under the polar coordinate system, its value can be obtained by the conventional triaxial test; the θ is the polar angle; e is the eccentricity of the partial plane, which controls the shape of the plastic potential surface in the range of θ = 0 – π/3. Compared with a variety of strength criteria, such as Matsuoka-Nakai [25], Lade [26], Yao [27] and Yu [28], the criteria used in the DP-S model are relatively conservative (Figure 6). For the accident of water inrush from KCC, it is safer to choose this strength criterion.
In order to get the expression of the permeability coefficient, we need a simplified hydraulic calculation model. Figure 7a shows the results of scanning electron microscopy (SEM) of rock samples when they are loaded to failure (pc = 10 MPa). Figure 7b shows that the rock sample formed a macroscopic fracture surface after loading. Both the results of mesoscopic test and the macroscopic observations reflect the fact that the rock sample yielded due to an increase of axial pressure, creating a rupture surface (damaged part) and two rock blocks (non-damaged part). The angle between the rupture surface and the horizontal plane is different due to different confining pressure, and the length of the flow path is changed, however, the k of rock sample is mainly determined by the plastic deformation and damage factors (d) [29]. In other words, the model of representative volume element (RVE) shown in Figure 7c can basically reflect the change of the permeability when the plastic deformation and damage grow in the rock sample. For the non-damaged part, the permeability is marked as k1, and the permeability of damaged part is marked as k2. For undamaged rock without plastic deformation (εp = 0), the permeability coefficient k is mainly determined by porosity n and elastic volume strain εve [22]. According to the cube law of fracture seepage and the definition of permeability coefficient [30], we can deduce the permeability coefficient k of the RVE. The equation is as follows:
k = { k 1 { ( 1 + ε ve ) 3 / n ( 1 n ) / [ n ( 1 + ε ve ) 1 / 3 ] } ,   ε p = 0 ( 1 d ) k 1 + d k 2 ( 1 + d ε ^ vp ) 3 ,   ε p > 0
where, n is the porosity; εve is the elastic volume strain; 0 ε ^ vp 1 , is a function of confining pressure pc, and ε ^ vp = d ε vp / [ B 1 ( p c / σ c ) + B 2 ] , εvp is the absolute value of plastic volume strain, σc is uniaxial compressive strength, the method of obtaining B1 and B2 is the same as that of A1 and A2.
Based on above equations, the PD-S model is established, and the numerical simulation process is shown in Figure 8.

4.2. Validation by Laboratory Tests

Laboratory tests are an effective means of validating a new constitutive model. Taking sandstone as an example, the verification process of the PD-S model is illustrated. For numerical simulations, calculation model is established in the finite element software ABAQUS, the parameters of rock specimen (Table 3) can be obtained by cyclic loading and unloading tests and permeability tests. The results obtained by the two methods are shown in Figure 9. The results of calculation are close to that of test, which shows the rationality of the PD-S model.
As shown in Figure 9, the stress-strain curves calculated by PD-S model agree well with the experimental results, indicating that the strength criterion, the non-associated flow rule and the mechanical parameters adopted in the model are appropriate. Similarly, the calculated values of k−ε1 are close to the experimental results under high confining pressure (pc ≥ 20 MPa), which indicates that the PD-S model and the calculation parameters are reasonable. However, in the condition of low confining pressure (pc ≤ 10 MPa), there is a discrepancy between the predicted value and the measured value at the end of the curve, and the reason is as follows: In the process of triaxial compression test, the rock sample breaks down at the softening section after the peak stress, and cause the permeability coefficient (k) to increase sharply. The confining pressure (pc) has a significant impact on the k in high confining pressure. Although the rock failed, it is difficult for the dislocation of rock blocks on both sides of the rupture surface to occur under the action of high confining pressure, such as the curves of k (pc = 20 MPa) and k (pc = 30 MPa). However, under low confining pressure, pc has little effect on the crack width after rock fracture, and the k increases first, then decreases and then increases: (1) when the axial displacement is applied, the fracture surface is uneven, the crack width increases gradually, and the k increases rapidly (Figure 10b). (2) As the axial displacement continues to increase, the rock blocks on both sides of fracture surface have a relative displacement, if the two “fit each other”, that is, the peak of a rock block is close to the trough of another rock block, making the fissure width decrease, thus the k decreases (Figure 10c). (3) If the axial displacement continues to increase, the peaks of two rock blocks are closed, the crack width increases, the k increases (Figure 10d).

5. Results and Discussion

5.1. Material Parameters and Calculation Results

The material parameters of the other materials can be obtained by cyclic loading and unloading tests and seepage tests under different confining pressures as described in Section 4.2. However, these parameters need to be reduced to one-sixth to one-tenth due to the presence of cracks in the natural rock mass [31]. In the process of mining advance, the goaf area is filled by a soft elastic material to simulate the support capability of the fallen rock from the roof. Young's modulus for this material is set at 470 MPa, and Poisson’s ratio is 0.25 [32]. The physical properties (or coring recovery rate) at different sections of the X1 are different (Figure 3), leading to the differences of the mechanical and hydraulic properties. Thus we defined the calculation parameters of the different parts of the X1. As shown in Table 4.
As mentioned above, the confining pressure pc affects the plastic strain γ ^ p eq and ε ^ p (Equation (7)), thereby affecting the permeability coefficient k (Equation (9)), the damage variable d (Equation (5)), as well as the strength parameters φ’ and c’ (Equation (6)). The effect of pc on the γ ^ p eq and ε ^ p can be determined by A1, A2, B1, B2 in Table 4, while the relationship between the γ ^ p eq , ε ^ p and k, φ’ and c’ are shown in Figure 11. And the relation between φ’, c’ and γ ^ p eq can be obtained by the statements in Section 4.1. It should be noted that, the rock below the limestone 1 and the lower part of X1 does not undergo plastic deformation by trial, so just showing the relation of ε ^ p k , γ ^ p eq c and γ ^ p eq φ’ of sandstone, the upper and middle part of X1.
According to the geological and hydrogeology conditions of the 2519 mining face (Figure 3), a plane strain numerical model with length of 300 m and height of 200 m was established. This model assumed that the longitudinal profile of the X1 was trapezoidal, with a top side 12.6 m, a bottom side 36 m and a height of 140 m, as shown in Figure 12. To obtain the mechanical and hydraulic responses, reasonable boundary conditions are established for the numerical model. Pressures of 17.7 MPa and 28.6 MPa were applied on the top and left and right boundary to simulate the vertical stress (or confining pressure) from the overburden and tectonic stress in the syncline basin (Figure 2), respectively, while pore water pressures of 6.5 MPa, 3.3 MPa and 3.0 MPa were applied in the three limestone aquifers from bottom to top. The vertical displacement constraints were applied to the bottom boundary of the model. Based on measured data, 20 m length was added after the first roof caving of 32 m, and backfill of 12 m length was repeated added in periodic caving. In order to explain the mechanism of water inrush from KCC, we chose the distance of mining face at L0 = −15 m, −7 m, 0 m before the center of the X1 and L0 = 7 m, 15 m after the X1, respectively. The results of the shear stress, confining pressure (or vertical stress), damage zone, pore water pressure and seepage vector were analyzed, as shown in Figure 12.

5.1.1. Shear Stress and Vertical Stress (i.e., Confining Pressure)

After the coal is excavated, the shear stress concentrates and forms a concentration zone with bubble shape (the first row, blue region, Figure 12) at the front of the coal face. Meanwhile the KCC under the goaf experiences a confining pressure release and forms a relief area with a flat elliptical shape (the second row, red region, Figure 12). The distribution of shear stress and confining pressure changes as mining advances. When L0 < −38 m (Figure 4), the stress in the column varies little. However, in the range of −15 m < L0 < 15 m, the stress in the column changed dramatically. The maximum shear stress increased from 0.014 MPa (the initial state) to 6.28 MPa (L0 = −7.0 m) or 6.69 MPa (L0 = 7 m), and then the KCC located in the relief area, the confining pressure deduced from 14.2 MPa (the initial state) to 3.59 MPa (L0 = 0 m). After experiencing such a special cycle loading and unloading path, damage could appear in the column.

5.1.2. Damage Zone

The damage zone starts from the above shear stress concentration zone, which determines the damage thickness of the X1, while the relief area of confining pressure determines the ultimate damage degree. When the coal face is away from the X1, the stress within the column changes little, and no damage occurs. When the L0 = −15 m, the damage thickness of the X1 reaches 31.6 m suddenly and the corresponding damage zone is focused on the left side-wall. In the range of −7 m < L0 < 7 m, although the shear stress near the two side-walls is close, the damage thickness at the left wall (37.7 m) is slightly larger than that of the right wall (35.1 m) due to the influence of pore water pressure. When L0 ≥ 15 m, the maximum damage thickness of the X1 is maintained at about 37.7 m, this will create a channel for groundwater inrush between L1 and coal seam (the third row, Figure 11).
It is important to note that in the relief zone, the damage degree continues increasing as the confining pressure decreases. Although this change is evident only in the depth range within 14 m of the floor near the free face of the goaf, the damage increase can cause drastic changes of the permeability coefficient of the KCC according to Equations (5) and (9), leading to the changes in pore water pressure and seepage velocity in the X1.

5.1.3. Pore Water Pressure and Seepage Vector

As mining advances, the pore water pressure in the KCC strongly depends on the confining pressure. Therefore, the distribution law of the two is similar. When L0 < −26 m, the pore water pressure in the X1 is almost unchanged. The maximum water pressure at the left wall increased from 2.8 MPa to 4.9 MPa when L0 = −15m. Once the coal face passes through the left wall (L0 = −7 m), the water pressure drops rapidly to −13.1 MPa. However, the phenomenon that the pore pressure increases first and then decreases does not occur at the right wall position due to the good permeability of the X1. The minimum pore pressure is −14.2 MPa when L0 = 15 m, the negative pore water pressure combined with the reduced confining pressure exacerbates the damage degree of the rock at the right wall.
The seepage velocity vector is the result of the combined action of the above factors (Equation (9)). When −15 m < L0 < −7 m, the high water pressure at the front of the coal face pushes the water in the X1 into limestone 1, 2, and even 3. The groundwater flows to the lower part of goaf along the limestone stratum, rapidly enters the low water pressure area and causes the water inrush from sandstone floor. In this process, the X1 acts as a vertical aquifer. However, the X1 acts as a dangerous vertical water-conducting channel when −7 m < L0 < 15 m. At this moment, the high water pressure area is located outside the X1, and drives the water in the limestone aquifer to enter the X1. Meanwhile the water in sandstone also enters the X1 under the influence of filling body, and emit from the goaf with low water pressure. The seepage rate reached the maximum at the left wall, 8.4 × 10−7 m/s (L0 = −3 m). While the seepage velocity is 6.2 × 10−7 m/s (L0 = 10 m) at the right wall. Although the seepage velocity is small, but the X1 has a greater seepage area and a smaller confining pressure at L0 = 10 m, resulting in the water inflow up to 53.8 m3/h, and the water inflow in left side-wall is only 32.5 m3/h. Therefore, the right side-wall has a greater risk of water inrush.

5.2. Discussion

In this article, we take Chensilou mine as an example, and study the mechanism of water inrush from KCC in complicated geological conditions by a series of methods, including connectivity experiments, water pressure monitoring tests in two side-walls, and numerical simulations based on the PD-S model. Connectivity experiments are an important in-situ test to ensure the safety of mining faces under the influence of KCC. The results obtained by the other two methods show that the KCC’s damage thickness at the right side-wall is less than that at the left side-wall, where the instantaneous seepage velocity reaches the maximum. However, due to the influence of confining pressure, the most dangerous position is the right side-wall, and the water inflow is larger than that at the left wall (Figure 12). This finding is significant because it clearly indicates that the most dangerous position is not in the maximum damage thickness for an inhomogeneous KCC connected three aquifers. The KCC’s damage zone, distribution of pore water pressure and seepage vector were obtained by numerical simulation based on the PD-S model, and verified by water pressure monitoring tests in two side-walls. To our knowledge, the constitutive theory considering the effect of confining pressure on rock strength, plastic deformation, damage and seepage, and the method of monitoring the relief time of water pressure were used herein for the first time in searching water inrush mechanism of KCC. As is known, the KCC is subjected to the cyclic loading and unloading during the mining advance [33]. The change of confining pressure could lead to the change of rock strength, plastic strain and permeability coefficient. Elasticity theory [14,15], elastic-plastic theory [16,17] and traditional fluid-solid coupling theory [19] have been used previously for predicting the safety of the mining face. However, the water inrush mechanism in Chensilou mine cannot be concluded from those models because they were conducted by the lack of the influence of confining pressure on permeability coefficient and plastic deformation (Equations (5) and (9)), lack of the influence of pore water pressure on damage and seepage (Equations (4) and (7)), and lack of the influence of complex geological factors, including multiple limestone aquifer and an inhomogeneous KCC (Figure 3 and Table 4). The new method for in-situ test, and the new model verified by laboratory tests (Figure 9), makes the results more meaningful.
To date, most researchers have used “strain methods” for their in situ tests [19,33]. This method could partially explain the change of strain in floor with the mining advance, but the reason for water inrush from KCC is that the development of rock fissures under the influence of mining, and there seems to be no direct relationship between fissures and strain [34]. Although a method of directly testing KCC’s damage was reported in a previous study using water injection tests [4], the exact location of the most dangerous position in the KCC was still unclear because they set the test section in the middle of the KCC. In addition, the water inflow in boreholes could interfere with each other during the water injection tests. In order to obtain damage thickness in the two side-walls of the KCC, an improved method named “water pressure monitoring test in two side-walls” has been developed. Through the above in-situ test and numerical simulation method, we obtained that the water inrush from KCC is a cascading and vicious cycle accident (Figure 13).
The PD-S model and water pressure monitoring test could be valuable tools with applications in deep buried, and cyclic loading and unloading conditions, such as coal mining in the seabed [35], and the co-exploitation of uranium, oil, gas, and coal under the influence of watery faults [36].

6. Conclusions

The main conclusions drawn from the investigations are summarized as follows:
  • The reasonableness of the in-situ test method determines the accuracy of the test results. Connectivity experiment is an important method to ensure the safety of coal mining. The water pressure monitoring test in two side-walls, which is suitable for good permeability of the KCC, has been developed. According to the relief time of water pressure, the damage thickness of the side-wall far away from the starting cut is smaller than that near the starting cut, and the maximum depth of damage zone reaches limestone 1.
  • Shear stress is one of the important factors in determining the damage thickness. In the process of mining advance, the shear stress concentration area below the coal face leads to formation of a damage zone, and this phenomenon is more pronounced near the side-wall of KCC due to a strong shear dislocation between the two materials with different mechanical properties. Consistent with the results of in-situ tests, the maximum damage thickness is about 37.7 m, and it is located in the left side-wall.
  • Confining pressure plays a crucial role in the water inrush mechanism of KCC. With the decrease of confining pressure, deviatoric stress can cause further plastic strain increases, rock damage and permeability, and a reduction of strength and stiffness. Meanwhile the pore water pressure changes from positive to negative, and the effective stress in the rock will damage the deeper rock, or force the rock to become more fractured in the damage zone. This broken rock can cause stress changes, especially the confining pressure in the KCC. Thus a cascade of vicious water inrush disaster cycles forms.

Acknowledgments

The author sincerely thanks Wang Qiang at Southeast University for providing workstation and related data access. This work was supported by the Research Innovation Program for College Graduates of Jiangsu Province (KYLX16_0537) and National Natural Science Foundation of China (51323004). The authors would like to thank anonymous referees for their careful reading of this article and valuable suggestions.

Author Contributions

J.W. and Z.M. conceived and designed the experiments; K.M. and G.W. performed the experiments; Y.D. and S.H. analyzed the data; H.B. contributed analysis tools; H.L. wrote the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Booth, C.J. Groundwater as an environmental constraint of longwall coal mining. Environ. Geol. 2006, 49, 796–803. [Google Scholar] [CrossRef]
  2. Odintsev, V.N.; Miletenko, N.A. Water inrush in mines as a consequence of spontaneous hydrofracture. J. Min. Sci. 2015, 51, 423–434. [Google Scholar] [CrossRef]
  3. Istvan, B.; Lucien, D.; Ferenc, S. Bayesian analysis of underground flooding. Water Resour. Res. 1982, 18, 1110–1116. [Google Scholar]
  4. Ma, D.; Miao, X.; Bai, H.; Huang, J.; Pu, H.; Wu, Y.; Zhang, G.; Li, J. Effect of mining on shear sidewall groundwater inrush hazard caused by seepage instability of the penetrated karst collapse pillar. Nat. Hazards 2016, 82, 73–93. [Google Scholar] [CrossRef]
  5. Cook, A.C. Australian Black Coal; Australian Institute of Mining and Metallurgy: Illawarra, Australia, 1975. [Google Scholar]
  6. Huang, Z.; Jiang, Z.; Qian, Z.; Cao, D. Analytical and experimental study of water seepage propagation behavior in the fault. Acta Geodyn. Geomater. 2014, 176, 361–370. [Google Scholar] [CrossRef]
  7. Li, H.; Bai, H.B.; Wu, J.J.; Zhao, H.M.; Ma, K. A method for prevent water inrush from karst collapse column: A case study from sima mine, china. Environ. Earth Sci. 2017, 76, 493. [Google Scholar] [CrossRef]
  8. Bai, H.B.; Ma, D.; Chen, Z.Q. Mechanical behavior of groundwater seepage in karst collapse pillars. Eng. Geol. 2013, 164, 101–106. [Google Scholar] [CrossRef]
  9. Qian, X.P. The formation of gypsum karst collapse-collum and its hydrogeological significance. Carsologica Sin. 1988, 4, 344–348. [Google Scholar]
  10. Yuan, D. Karstology of China; Geological Publishing House: Beijing, China, 1994. [Google Scholar]
  11. Shi, J.; Lian, D.; Yang, S. A discussion of karst collapse. Geol. Min. Res. North China Mag. 1998, 13, 264–266. [Google Scholar]
  12. Yin, S.X.; Wu, Q.; Wang, S.X. Water-bearing characteristics and hydro-geological models of karstic collapse columns in north China. Chin. J. Rock. Mech. Eng. 2005, 24, 77–82. [Google Scholar]
  13. Zhang, J.F.; Zhang, H.L.; Meng, D.; Cao, J. Numerical simulation of rock deformation and seepage field with a fully-water karstic collapse column under mining influence. Chin. J. Rock. Mech. Eng. 2009, 28, 77–82. [Google Scholar]
  14. Li, Z.; Feng, G.; Zhai, C. Study on “triangle” water-inrush mode of strong water-guide collapse column. J. Cent. South. Univ. 2016, 23, 2402–2409. [Google Scholar] [CrossRef]
  15. Zhang, C.; Tu, S. Control technology of direct passing karstic collapse pillar in longwall top-coal caving mining. Nat. Hazards 2016, 84, 17–34. [Google Scholar] [CrossRef]
  16. Yin, S.X.; Zhang, J. Impacts of karst paleo-sinkholes on mining and environment in northern China. Environ. Geol. 2005, 48, 1077–1083. [Google Scholar] [CrossRef]
  17. Ou, S.; Wang, L.; Wang, P. Numerical analysis of seepage flow characteristic of collapse column under the influence of mining. Int. J. Min. Sci. Technol. 2013, 23, 237–244. [Google Scholar] [CrossRef]
  18. Zhang, J.C. Investigations of water inrushes from aquifers under coal seams. Int. J. Rock. Mech. Min. Sci. 2005, 42, 350–360. [Google Scholar] [CrossRef]
  19. Zhu, W.C.; Wei, C.H. Numerical simulation on mining-induced water inrushes related to geologic structures using a damage-based hydromechanical model. Environ. Earth Sci. 2011, 62, 43–54. [Google Scholar] [CrossRef]
  20. Liu, E.L.; He, S.M. Effects of cyclic dynamic loading on the mechanical properties of intact rock samples under confining pressure conditions. Eng. Geol. 2012, 125, 81–91. [Google Scholar] [CrossRef]
  21. Bai, H.B.; Miao, X.X. Hydrogeological characteristics and mine water inrush prevention of late Paleozoic coalfields. J. China Univ. Min. Technol. 2016, 45, 1–10. [Google Scholar]
  22. Li, L.C.; Tang, C.A.; Liang, Z.Z.; Ma, T.H.; Zhang, Y.B. Numerical simulation on water inrush process due to activation of collapse columns in coal seam floor. J. Min. Saf. Eng. 2009, 26, 158–162. [Google Scholar]
  23. Boresi, A.P.; Lynn, P.P.; Hung, Y.Y. Elasticity in Engineering Mechanics; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  24. Fei, K.; Zhang, J.W. Application of ABAQUS in Geotechnical Engineering; China Water & Power Press: Beijing, China, 2013. [Google Scholar]
  25. Matsuoka, H.; Nakai, T. Relationship among Tresca, Mises, Mohr-Coulomb and Matsuoka-Nakai failure criteria. Soils Found. 2008, 25, 123–128. [Google Scholar] [CrossRef]
  26. Woodward, P.K. Earth pressure coefficients based on the Lade-Duncan failure criterion. Eng. Struct. 1997, 19, 733–737. [Google Scholar] [CrossRef]
  27. Yao, Y.P.; Zhang, B.Y.; Zhu, J.G. Behaviours, constitutive models and numerical simulation of soils. China Civ. Eng. J. 2012, 45, 127–150. [Google Scholar]
  28. Zan, Y.W.; Yu, M.H. Generalized nonlinear unified strength theory of rock. J. Southwest Jiaotong Univ. 2013, 48, 616–624. [Google Scholar]
  29. Luo, Z.H.; Zhu, Z.D.; Ruan, H.N.; Shi, C. Extraction of microcracks in rock images based on heuristic graph searching and application. Comput. Geosci. 2015, 85, 22–35. [Google Scholar] [CrossRef]
  30. Yin, L.M.; Chen, J.T.; Sun, W.B. Experimental Study on three dimensional coupled stress-seepage law of single fracture. Appl. Mech. Mater. 2013, 353, 524–528. [Google Scholar] [CrossRef]
  31. Peng, S.S. Longwall Mining; Science Press: Beijing, China, 2011. [Google Scholar]
  32. Cheng, Y.M.; Wang, J.A.; Xie, G.X.; Wei, W.B. Three-dimensional analysis of coal barrier pillars in tailgate area adjacent to the fully mechanized top caving mining face. Int. J. Rock. Mech. Min. Sci. 2010, 47, 1372–1383. [Google Scholar] [CrossRef]
  33. Yin, H.; Lefticariu, L.; Wei, J. In situ dynamic monitoring of stress revolution with time and space under coal seam floor during longwall mining. Environ. Earth Sci. 2016, 75, 1249–1260. [Google Scholar] [CrossRef]
  34. Wang, Y.; Yang, W.; Li, M.; Liu, X. Risk assessment of floor water inrush in coal mines based on secondary fuzzy comprehensive evaluation. Int. J. Rock. Mech. Min. Sci. 2012, 52, 50–55. [Google Scholar] [CrossRef]
  35. Gao, Y.; Zhang, Q.S.; Xu, B.S.; Li, W. Study of mining roof abutment pressure distribution law. Rock Soil Mech. 2010, 31, 1309–1313. [Google Scholar]
  36. Wang, Y.; Yang, W.L.; Deng, J.; Wu, B.L.; Li, Z.Y.; Wang, M.Z. Accumulation system of cohabitating muti-energy minerals and their comprehensive exploration in sedimentary basin—A case study of Ordos basin, NW China. Acta Geol. Sin. 2014, 88, 815–824. [Google Scholar]
Figure 1. Formation process of KCC. (a) Connection of caves; (b) initial stage; (c) growth stage; (d) final stage.
Figure 1. Formation process of KCC. (a) Connection of caves; (b) initial stage; (c) growth stage; (d) final stage.
Energies 10 01938 g001
Figure 2. Location of Chensilou mine and study area.
Figure 2. Location of Chensilou mine and study area.
Energies 10 01938 g002
Figure 3. The profile and core of X1, and stratigraphic characteristics of surrounding rock.
Figure 3. The profile and core of X1, and stratigraphic characteristics of surrounding rock.
Energies 10 01938 g003
Figure 4. The location of the test points (top view).
Figure 4. The location of the test points (top view).
Energies 10 01938 g004
Figure 5. The results of water pressure monitoring test in different vertical depth.
Figure 5. The results of water pressure monitoring test in different vertical depth.
Energies 10 01938 g005
Figure 6. Shape of yield surface in deviatoric plane and meridian plane.
Figure 6. Shape of yield surface in deviatoric plane and meridian plane.
Energies 10 01938 g006
Figure 7. Simplified calculation model of k. (a) microscopic fracture surface; (b) macroscopic fracture surface; (c) RVE model.
Figure 7. Simplified calculation model of k. (a) microscopic fracture surface; (b) macroscopic fracture surface; (c) RVE model.
Energies 10 01938 g007
Figure 8. Numerical calculation process of PD-S model.
Figure 8. Numerical calculation process of PD-S model.
Energies 10 01938 g008
Figure 9. Comparison of experimental and calculated results under different confining pressures.
Figure 9. Comparison of experimental and calculated results under different confining pressures.
Energies 10 01938 g009
Figure 10. Fracture morphology and variation of permeability coefficient of rock under low confining pressure. (a) Intact rock sample; (b) Cracks appear in the rock sample and the k value increase; (c) The rock blocks on both sides of the rupture face slide and agree with each other, and the k value decreases; (d) The crack width increases under low confining pressure, and the k value increases.
Figure 10. Fracture morphology and variation of permeability coefficient of rock under low confining pressure. (a) Intact rock sample; (b) Cracks appear in the rock sample and the k value increase; (c) The rock blocks on both sides of the rupture face slide and agree with each other, and the k value decreases; (d) The crack width increases under low confining pressure, and the k value increases.
Energies 10 01938 g010
Figure 11. The evolution curve of c’, φ’ and k. (a) The relationship between c’, φ’ and γ ^ p eq ; (b) The relationship between k and ε ^ vp .
Figure 11. The evolution curve of c’, φ’ and k. (a) The relationship between c’, φ’ and γ ^ p eq ; (b) The relationship between k and ε ^ vp .
Energies 10 01938 g011
Figure 12. Numerical calculation model and results.
Figure 12. Numerical calculation model and results.
Energies 10 01938 g012
Figure 13. The process of cascade and vicious cycle for water inrush from KCC.
Figure 13. The process of cascade and vicious cycle for water inrush from KCC.
Energies 10 01938 g013
Table 1. Connectivity experiment between P01 and P02.
Table 1. Connectivity experiment between P01 and P02.
Water Pressure in P01/MpaInjection Time/minWater Flow in P02/m3
0–1150.101
1–2150.132
2–4150.227
4–4.5150.413
4.5–5150.682
Table 2. Connectivity experiment between P03 and P02.
Table 2. Connectivity experiment between P03 and P02.
Water Pressure in P03/MpaInjection Time/minWater Flow in P02/m3
0–2150.110
2–4150.152
4–4.5150.235
4.5–5150.416
5–5.2150.635
Table 3. Calculation parameters of sandstone (rock samples).
Table 3. Calculation parameters of sandstone (rock samples).
E0/GPaμeaψc’/MPaφ’/°pw/MPak0(k1)/ms−1k2/ms−1αn/%A1A2B1B2
20.400.260.6670.01062325.2323.51.98 × 10−124.56 × 10−90.763.17.620.36.170.54
Table 4. Material parameters (rock mass).
Table 4. Material parameters (rock mass).
E0/GPaμeaψc’/MPaφ’/°pw/MPak0(k1)/ms−1k2/ms−1αn/%A1A2B1B2
Coal1.520.250.6670.112010.22309.73 × 10−127.17 × 10−90.013.04.120.214.830.37
Mudstone1.310.320.6670.15174.72103.62 × 10−126.18 × 10−90.012.53.160.122.230.20
Sandstone8.730.260.6670.0102215.22809.37 × 10−101.04 × 10−70.8087.820.336.060.74
Limestone16.70.220.6670.0072818.7353.0/3.3/7.01.05 × 10−99.25 × 10−60.858.39.310.457.570.82
Upper X11.030.250.6670.12205.3262.04.76 × 10−93.24 × 10−60.887.22.860.142.160.30
Middle X10.950.300.6670.15183.0233.38.5 × 10−97.18 × 10−60.918.53.580.202.540.25
Lower X10.780.320.6670.22122.1187.04.36 × 10−88.23 × 10−50.979.64.640.373.660.28

Share and Cite

MDPI and ACS Style

Li, H.; Bai, H.; Wu, J.; Ma, Z.; Ma, K.; Wu, G.; Du, Y.; He, S. A Cascade Disaster Caused by Geological and Coupled Hydro-Mechanical Factors—Water Inrush Mechanism from Karst Collapse Column under Confining Pressure. Energies 2017, 10, 1938. https://doi.org/10.3390/en10121938

AMA Style

Li H, Bai H, Wu J, Ma Z, Ma K, Wu G, Du Y, He S. A Cascade Disaster Caused by Geological and Coupled Hydro-Mechanical Factors—Water Inrush Mechanism from Karst Collapse Column under Confining Pressure. Energies. 2017; 10(12):1938. https://doi.org/10.3390/en10121938

Chicago/Turabian Style

Li, Hao, Haibo Bai, Jianjun Wu, Zhanguo Ma, Kai Ma, Guangming Wu, Yabo Du, and Shixin He. 2017. "A Cascade Disaster Caused by Geological and Coupled Hydro-Mechanical Factors—Water Inrush Mechanism from Karst Collapse Column under Confining Pressure" Energies 10, no. 12: 1938. https://doi.org/10.3390/en10121938

APA Style

Li, H., Bai, H., Wu, J., Ma, Z., Ma, K., Wu, G., Du, Y., & He, S. (2017). A Cascade Disaster Caused by Geological and Coupled Hydro-Mechanical Factors—Water Inrush Mechanism from Karst Collapse Column under Confining Pressure. Energies, 10(12), 1938. https://doi.org/10.3390/en10121938

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop