Computational Methods and Applications to Simulate WaterRelated Natural Hazards
View this Special IssueResearch Article  Open Access
Mechanism of Fracturing in Shaft Lining Caused by HighPressure Pore Water in Stable Rock Strata
Abstract
With the increase in shaft depth, the problem of cracks and leakage in singlelayer concrete lining in porous waterrich stable rock strata has become increasingly clear, in which case the mechanism of fracturing in shaft lining remains unclear. Considering that the increase in pore water pressure can cause rock mass expansion, this paper presents the concept of hydraulic expansion coefficient. First, a cubic model containing spherical pores is established for studying hydraulic expansion, and the ANSYS numerical simulation, a finite element numerical method, was used for calculating the volume change of the model under the pore water pressure. By means of the multivariate nonlinear regression method, the regression equation of the hydraulic expansion coefficient is obtained. Second, based on the hydraulic expansion effect on the rock mass, an interaction model of pore water pressure–porous rock–shaft lining is established and further solved. Consequently, the mechanism of fracturing in shaft lining caused by highpressure pore water is revealed. The results show that the hydraulic expansion effect on the surrounding rock increases with its porosity and decreases with its elastic modulus and Poisson’s ratio; the surrounding rock expansion caused by the change in pore water pressure can result in the outer edge of the lining peeling off from the surrounding rock and tensile fracturing at the inner edge. Therefore, the results have a considerable guiding significance for designing shaft lining through porous waterrich rock strata.
1. Introduction
In recent years, with the implementation of China’s western development strategy, the abundant coal resources in some waterrich bedrock areas (Shaanxi, Ningxia, Gansu, Inner Mongolia, etc.) have been mined on a large scale [1]. Since the shaft is the passageway of mine production, the shaft lining should be designed with sufficient strength, rigidity, and waterproof performance to meet the requirements of safety, sealing, and durability. However, with the increase in coal mining depth, the nonminingrelated deformation and fracture of the shaft lining under the action of high water pressure is still very serious. At present, scholars at home and abroad have carried out a lot of considerable research on shaft fracture caused by groundwater activities. Regarding the study topic on groundwater seepage, Farmer [2] studied the hydrostatic stresses acting on shaft lining based on the theory of surrounding rock permeability. Bear [3] carried out a microscopic study on the seepage mechanism, indicating that the change in pore shape is a sensitive factor that causes permeability coefficient changes. Bruno [4] studied the effect of pore pressure on the tensile fracture of rock on the basis of Biot theory and noted the importance of fluidsolid coupling in hydraulic fracturing. Zimmerman [5] briefly deduced the linear porouselasticity and thermoelastic equation by using the threedimensional fluidsolid coupling seepage model. Jiang [6] used a highpressure permeability test to simulate the deformation law of rock mass fracturing under the action of pore water pressure. Regarding the study topic on hydrophobic settlement, Lou [7] derived a general formula for calculating the additional stress in the shaft lining based on the drainage consolidation theory. According to the principle of superposition and strain compatibility, Yang [8] derived the theoretical solution of vertical additional stress on the shaft lining by using a numerical method. Wang [9] obtained the additional stress value at complex alluvium strata by the negative friction coefficient between the shaft lining and surrounding rock. Zhang [10] studied the fluidsolid coupling numerical simulation of sharp severely aquifer drainage in underground mining and analyzed the influence of drainage location and drainage rate on the stability of shaft lining. Regarding the study on unstable rock strata, Sun [11] simulated and analyzed the deformation and stress characteristics of surrounding rock and shaft lining under the interactive geological conditions of soft and hard rock strata, and the conclusion that the shear failure occurred in shaft lining under the inhomogeneous pressure of surrounding rock was obtained. In addition, Yang [12] deduced the elastic approximate stress and displacement solution of irregular inclined shaft lining subjected to water pressure and proposed the optimal design for inclined shaft lining. Meng [13] classified and studied the constructing techniques for the inclined shaft penetrating the drift sand stratum, and the problems such as the stability of inclined shaft structure and the sealing of water and sand were solved. It is well known that shaft lining in porous waterrich rock strata is constructed by the temporary drainage method or the frozen water method. When the drainage or freezing is stopped, the pore water pressure of the surrounding rock will rise, and at this time, cracking and leakage of lining often occurs. However, no leakage occurs before the shaft lining fractures, inhibiting seepage in the surrounding rock. Additionally, the theory of vertical additional stress caused by mining hydrophobic settlement is applicable to mainly alluvium. Clearly, although some existing theories have matured, the fracture mechanism of shaft lining in highpressure waterrich stable rock strata has not been clarified.
Based on the compressive and expansive deformation properties of the porous media under the action of stress, in analogy with the coefficient of thermal expansion, this paper defines the linear expansion coefficient of the rock caused by unit pore water pressure as the hydrostatic expansion coefficient, which is expressed as . With the advantage of solving complex problems that cannot be solved by theoretical and experimental research and having the characteristics suitable for any problem geometry and boundary conditions [14], the finite element numerical method has become one of the most effective methods to solve scientific and engineering problems. Therefore, the ANSYS numerical simulation is adopted to obtain the volume change of the cubic model under the pore water pressure. Based on the calculation results, the approximate analytical solution of the hydraulic expansion coefficient is derived by means of the multivariate nonlinear regression method. By comprehensively analyzing the influence of various factors on the shaft lining stress, the mechanism of fracturing in shaft lining due to highpressure pore water is clarified, providing an effective and scientific basis for lining safety.
2. Hydraulic Expansion Coefficient Solution
2.1. Basic Assumptions
The rock mass is composed of an infinite number of identical cubic microunits containing spherical pores of equal diameter, and the rock is a homogeneous and isotropic linear elastic medium.
The effect of pore water pressure on a microunit can be regarded as free expansion, and the effect on the rock is the superposition of the expansive effect of each microunit [15].
The elastic modulus, Poisson’s ratio, and porosity of the cubic model containing spherical pores are expressed by E_{0}, μ_{0,} and n_{0}, respectively. Since the elastic modulus reflects the degree of rock elastic deformation, E_{0} is clearly inversely proportional to α. According to the basic assumptions, the hydraulic expansion coefficient of the entire porous waterrich rock is equal to that of any single cubic microunit containing spherical pores under conditions of a free boundary. Assuming is the volumetric hydraulic expansion coefficient of the rock, when the infinitesimal of higher order is ignored. Considering that n_{0} and μ_{0} are both dimensionless parameters, then αE_{0} can be studied as dimensionless. When the compressive stress is specified as positive, the following equations can be derived from the theory of elastic mechanics.where is the total volume of the cubic model, is the pore water pressure, r_{0} is the pore radius, and is the half length of the microunit.
2.2. Solution Result
According to the characteristics of the solution model, a 1/8 cubic model containing spherical pore is established by ANSYS numerical simulation, and the volume change in the model with different values of n_{0} and μ_{0} under the action of pore water pressure is calculated. When , taking n_{0}=0.1, μ_{0}=0.2 as an example, the model after mesh generation is shown in Figure 1, and the equivalent displacement cloud of the model after solution is shown in Figure 2.
By means of the same method, can be calculated with different parameters of cubic model, as shown in Figure 3.
Based on the changing laws of the data in Figure 3, the multivariate nonlinear regression method is adopted to obtain the expression of for n_{0} and μ_{0}, which satisfies that when , . Toward this goal, the objective function can be expressed aswhere n_{0} satisfies the general rock porosity value; i.e., , and are the regression coefficients, the values of which are shown in Table 1.
 
Note: R^{2} is the goodness of fit. 
Substitute (4) into (1). Then,
3. The Variation Law of the Coefficient of Hydraulic Expansion
3.1. Regular Analysis
Analysis of the relationship between the hydraulic expansion coefficient and the various parameters is helpful to better understand the influence law of hydraulic expansion on underground structures in different waterrich rock strata and has an important guiding significance for mitigation of the shaft lining fracture problem. Therefore, the partial derivatives of n_{0} and μ_{0} for (5) are solved, and (6) and (7) are always satisfied.
The results show that the hydraulic expansion coefficient increases with n_{0} and linearly decreases with μ_{0}.
3.2. Influence of Pore Shape on Hydraulic Expansion Effect
Rock is a typical porous medium, in which the internal pores are not a sphere but an irregular space with various forms and a fractal dimension [16]. Therefore, an equalporosity cubic model with a regular Npolyhedron is adopted to analyze the influence of pore shape on the hydrostatic expansion effect, where N=4, 6, 8, 12, 20, and ∞ (sphere). If the hydraulic expansion coefficient of the model with a regular Npolyhedron pore is expressed as α_{N}, then . Taking μ_{0}=0.2 as an example, the relationship between and 1/N at different n_{0} values can be obtained by the ANSYS numerical simulation (see Figure 4). The expressions of the fitting function in Figure 4 can be expressed in the form of (8), and the corresponding fitting coefficient values are shown in Table 2.
 
Note: Adj. R^{2} is the goodness of fit after correction. 
As shown in Figure 4, increases with n_{0} and decreases with N. The hydraulic expansion coefficient of the model containing the tetrahedral pore is the largest, while that of the model containing the spherical pore is the smallest.where η_{1} and η_{2} are both functions in terms of n_{0}.
Based on the data in Table 2, the polynomial regression formulae of η_{1} and η_{2} can be calculated, as shown in Figures 5 and 6, respectively.
Since both and are functions of regarding n_{0} and N, the function f(n_{0}, N) can be used to represent the value of . When μ_{0}=0.2, the expression of f(n_{0}, N) can be derived from (5) and (8) as follows:
The value of f(n_{0}, N) for different rock porosities can be solved by (9). Therefore, when μ_{0} is between 0.1 and 0.3 and n_{0} is between 0.02 and 0.4, the range of f(n_{0}, N) can be obtained, as seen in Table 3.

As seen in Table 3, the hydraulic expansion coefficient of the model containing regular polyhedral pores is approximately (1.0~7.8) α and decreases with the number of pore faces.
4. Mechanism of Shaft Lining Fracture Analysis
4.1. Stress and Displacement Solution
4.1.1. Basic Assumptions
In the absence of water, the elastic moduli of the surrounding rock and the shaft lining are E_{1} and E_{2}, the corresponding Poisson’s ratios are μ_{1} and μ_{2}, the corresponding porosities are n_{1} and , respectively, and the porosity of the contact surface between the surrounding rock and the shaft lining is .
According to the nature of the problem, the following assumptions are made.
The elastic modulus and Poisson’s ratio of the surrounding rock mass and the rock matrix approximately satisfy the linear relations as follows:where k_{1} and k_{2} are constants related to pore shape, for spherical pores [17]: k_{1}=2.08, and k_{2}=0.345.
The porosity of contact surface n_{3} can be calculated according to the probability statistical method, namely, n_{3}=n_{1}+n_{2}n_{1}n_{2}.
The shaft lining is a homogeneous, continuous, and isotropic linear elastic medium.
The expansion effect of the pore water pressure on the surrounding rock can be equivalent to the volume expansion of elastomer, which is similar to the thermal expansion.
The inner and outer radii of the shaft lining are r_{1} and r_{2}, respectively, and the inner and outer radii of the surrounding rock are r_{2} and r_{3}, respectively. The outer boundary of the surrounding rock is fixed, as shown in Figure 7.
Assume that the interaction force between the surrounding rock and the shaft lining under the action of pore water pressure is a compressive stress, and the stress direction is defined as positive. The radial displacement is expressed by u_{r}, and the radial stress and circumferential stress are expressed by and , respectively, and the radial strain and circumferential strain are expressed by and , respectively. The total stress at the interface is p_{0}, and the effective stress is p_{1}. The bonding strength is f_{1}, the uniaxial compressive design strength and the uniaxial tensile design strength of the concrete are and , respectively, and the compressive design strength and the tensile design strength of the reinforcement are and , respectively. To facilitate the analysis, the relevant parameters are dimensionless as follows:
4.1.2. Stress and Displacement Solution in the Surrounding Rock Zone
The equations of equilibrium in the surrounding rock zone can be expressed as [18]
The equations of geometry for the axisymmetric plane strain problem can be written as
The equations of constitution in the surrounding rock zone under the action of water pressure expansion can be expressed as
The boundary conditions of the surrounding rock can be expressed as
According to (13)~(18), the analytical solutions of the stress and displacement in the surrounding rock zone can be obtained as
4.1.3. Stress and Displacement Solution in the Shaft Lining Zone
The equations of equilibrium and geometry in the shaft lining zone are the same as those in the surrounding rock zone, and the equations of constitution can be expressed as
The boundary conditions of the shaft lining can be expressed as
The analytical solutions of the stress and displacement in the shaft lining zone can be derived as
4.2. Stress Analysis at the Interface between the Surrounding Rock and Shaft Lining
4.2.1. Total Stress Analysis
Considering the conditions of singlevalued displacement at the interface between the surrounding rock and the shaft lining, can be solved by (21) and (28) as follows:
The partial derivative of for is
According to (30), the smaller the value of is, the larger the value of is. Therefore, when the outer boundary of the surrounding rock is farther from the interface between the lining and the surrounding rock, the interaction force at the interface is greater under the unit pore water pressure. When , the interaction force is a constant value equal to
By analyzing (31), it can be found that(1)when , ;(2)when , ;(3)when , .
The results show that when the elastic modulus of the surrounding rock is too large compared with the shaft lining, p_{0} is close to 0; otherwise, p_{0} is constant. Additionally, when the lining thickness is very small, p_{0} approaches 0; when the pore water pressure increases (), p_{0} is a compressive stress, and a compressive fracture forms at the inner edge of the shaft lining. Finally, when the pore water pressure decreases (), p_{0} is a tensile stress, and a tensile fracture forms at the inner edge of the shaft lining.
4.2.2. Effective Stress Analysis
According to the composition of the stress at the contact surface between the surrounding rock and the shaft lining, the effective stress (or skeleton stress) at the contact surface can be expressed as
By analyzing (32), it can be found that
when , then ; that is, the contact skeleton stress is equal to 0;
when the pore water pressure increases, if , then , and the skeleton stress at contact surface is tensile stress; if , then , and the skeleton stress at contact surface is compressive stress;
when the pore water pressure decreases, if , then , and the skeleton stress at contact surface is compressive stress; if , then , and the skeleton stress at contact surface is tensile stress.
The results show that when the skeleton stress at the contact surface is subjected to tensile stress, the shaft lining has a risk of peeling off the surrounding rock.
4.3. Shaft Lining Fracture Condition Analysis
4.3.1. Fracture Condition at Inner Edge of Shaft Lining
When p_{0} is the compressive stress, the inner edge of the shaft lining is in compression. According to the Code [19], the compressive failure condition at the inner edge of the shaft lining iswhere is the improvement coefficient of the concrete strength under multiaxial stress, is related to and (in general, .2), and is the minimum steel content.
When p_{0} is a tensile stress, the inner edge of the shaft lining is in tension. According to the Code [20], the tensile failure condition at the inner edge of the shaft lining is
4.3.2. Stripping Condition at Outer Edge of Shaft Lining
To avoid the hydrostatic action of groundwater on the entire outside surface of the shaft lining, there must be a certain bonding strength between the shaft lining and surrounding rock to ensure that the two are not separated. Therefore, when the contact surface skeleton is subjected to tension, the condition for peeling of the outer edge of the shaft lining away from the surrounding rock is
The Code [20] stipulates that the bonding strength of shotcrete should be no less than 0.8 MPa with Class I and Class II surrounding rock and should be no less than 0.5 MPa with Class III surrounding rock, respectively. However, the bonding strength between the pouring concrete and the surrounding rock is less than these values. In fact, the maximum bonding strength between the shaft lining and the surrounding rock in engineering is less than 1 MPa.
5. Engineering Example
A shaft passes through waterrich bedrock strata with Class III surrounding rock. The radius of the shaft is r_{1}=3 m, and the shaft lining is made of C60 concrete and HRB335 steel, that is, f_{c}=27.5 MPa, f_{t}=2.04 MPa, f_{y}=300 MPa, and MPa. The other parameters are as follows: E_{0}=15 GPa, E_{2}=36 GPa, μ_{0}=0.16, μ_{2}=0.2, n_{1}=0.15, n_{2}=0.12, ν=0.2%, f_{1}=0.8 MPa. For a change in groundwater level in the waterrich bedrock strata of 500 m, the fracturing in the lining with different thicknesses is shown in Figures 8 and 9.
The calculation results show that . Figures 8 and 9 show that when the water level rises, there is no compressive failure at the inner edge of the shaft lining, while stripping failure occurs at the outer edge of the shaft lining when ; when the water level drops, only the inner edge of the shaft lining undergoes tensile fracturing.
To further study the relationship between the stress and the deformation at the outer edge of the shaft lining, the lateral stiffness coefficient K_{c} of the lining can be calculated from (28) and (32), which can be derived as
According to (36), the values of K_{c} with different thicknesses are shown in Figure 10.
As shown in Figure 10, K_{c} decreases approximately linearly with , indicating that the increase in shaft lining thickness can effectively improve the lateral stiffness coefficient.
6. Conclusions
Based on the cubic model of hydraulic expansion containing a spherical pore, its volume change is calculated by the ANSYS numerical simulation, and the approximate analytical solution of the hydraulic expansion coefficient is obtained by the multivariate nonlinear regression method. Finally, the mechanism of fracturing in singlelayer concrete lining caused by highpressure pore water in stable rock strata is revealed by analyzing the interaction between the pore water pressure and the porous rock and shaft lining, and the following conclusions are obtained.
The hydraulic expansion effect increases with the porosity of the surrounding rock and decreases with the elastic modulus and Poisson’s ratio of the surrounding rock. Therefore, by reducing the porosity and improving the strength of the surrounding rock, the expansion effect of the surrounding rock on shaft lining can be weakened.
The mechanism of fracturing in shaft lining in waterrich bedrock strata is as follows: the hydraulic expansion effect on the surrounding rock may lead to stripping failure at the outer edge or tensile failure at the inner edge of the shaft lining.
When designing the shaft lining of waterrich bedrock section through the stable rock strata, the influence of pore water pressure change should be fully considered, which is to ensure that the shaft lining not only can take full advantage of deformation and compression effect but also has high bearing capacity.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Key Research and Development Program of China (2016YFC0600904). The authors would like to thank the editor and reviewers for their contributions on the paper.
References
 Z. S. Yao, H. Cheng, and C. X. Rong, “Shaft structural design and optimization of deep freezing,” Journal of China Coal Society, vol. 35, no. 5, pp. 760–764, 2010. View at: Google Scholar
 I. W. Farmer and D. H. Jennings, ““Effect of strata permeability on the radial hydrostatic pressures on mine linings,” International Journal of Mine Water, vol. 2, no. 3, pp. 17–24, 1983. View at: Publisher Site  Google Scholar
 J. Bear, “Dynamics of fluids in porous media,” Engineering Geology, vol. 7, no. 2, pp. 174175, 1972. View at: Google Scholar
 M. S. Bruno and F. M. Nakagawa, “Pore pressure influence on tensile fracture propagation in sedimentary rock,” International Journal of Mining Science and Technology, vol. 28, no. 4, pp. 261–273, 1991. View at: Publisher Site  Google Scholar
 R. W. Zimmerman, “Coupling in poroelasticity and thermoelasticity,” International Journal of Rock Mechanics and Mining Sciences, vol. 37, no. 12, pp. 79–87, 2000. View at: Publisher Site  Google Scholar
 Z. Jiang, S. Feng, and S. Fu, “Coupled hydromechanical effect of a fractured rock mass under high water pressure,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 4, no. 1, pp. 88–96, 2012. View at: Publisher Site  Google Scholar
 G. D. Lou and L. F. Su, “Analysis of loading on shaft lining subjected to alluvium settlement due to waterdrainage,” Journal of China Coal Society, vol. 16, no. 4, pp. 54–62, 1991. View at: Google Scholar
 W. H. Yang and H. L. Fu, “Theoretical investigation on the vertical additional force of shaft lining in special stratum,” Journal of China University of Mining and Technology, vol. 9, no. 2, pp. 130–136, 1999. View at: Google Scholar
 Y. Wang, C. Zhang, L. Xue, and X. Huang, “Prediction and safety analysis of additional vertical stress within a shaft wall in an extrathick alluvium,” Mining Science and Technology, vol. 20, no. 3, pp. 350–356, 2010. View at: Google Scholar
 Z. Zhang, X. Tang, and Q. Yu, “Impact of aquifer drainage on shaft lining in underground mining of Eastern China,” Electronic Journal of Geotechnical Engineering, vol. 21, no. 13, pp. 4721–4729, 2016. View at: Google Scholar
 X. Sun, G. Li, C. Zhao, Y. Liu, and C. Miao, “Investigation of deep mine shaft stability in alternating hard and soft rock strata using threedimensional numerical modeling,” Processes, vol. 7, no. 1, p. 17, 2019. View at: Publisher Site  Google Scholar
 R. S. Yang, Q. X. Wang, and S. Z. Chen, “Elactic analysis of irregular inclined shaft lining subjected to water pressure,” Journal of China University of Mining and Technology, vol. 46, no. 1, pp. 48–57, 2017. View at: Google Scholar
 Q. B. Meng, L. J. Han, R. J. Shi et al., “Study and application of construction technology for inclined shafts penetrating drift sand strata in coal mine,” Chinese Journal of Geotechnical Engineering, vol. 37, no. 5, pp. 900–910, 2015. View at: Google Scholar
 O. C. Zienkiewicz, R. L. Taylor, and J. Z. Zhu, The Finite Element Method: Its Basis and Fundamentals, Elsevier/Butterworth Heinemann, Amsterdam, Netherlands, 7th edition, 2013. View at: Publisher Site  MathSciNet
 J. Rutqvist, L. Börgesson, M. Chijimatsu et al., “Thermohydromechanics of partially saturated geological media:govening equations and formulation of four finite element models,” International Journal of Rock Mechanics and Mining Sciences, vol. 38, no. 1, pp. 105–127, 2001. View at: Publisher Site  Google Scholar
 C. H. Zhang, F. Jin, Y. L. Hou et al., DiscreteContactFracture Analysis of Rock and Concrete, Tsinghua University Press, Beijing, China, 2008.
 D. L. Bo, Hydraulic Load of Underground Structure in Pore Water Bearing Rock, China University of Mining and Technology, 2015.
 S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, McGrawHill, New York, NY, USA, 1970. View at: MathSciNet
 M. O. H. U.R. D. P. R. China, Code for Design of Concrete Structures (GB 500102010), China Architecture & Building Press, Beijing, China, 2015.
 C. M. C. A., “Technical Code for Engineering of Ground Anchorages and Shotcrete Support (GB 50086),” Tech. Rep., China Planning Press, Beijing, China, 2015. View at: Google Scholar
Copyright
Copyright © 2019 Jihuan Han et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.