Abstract

Research on the damage and seepage characteristics of unloading rock with high water pressure can help to further understand the mechanism of water inrush in deep mine floor and prevent water inrush. This paper used the RFPA2D-flow finite element software to study the failure and seepage characteristics of unloading rock with high water pressure and high stress and comparatively analyzed the failure modes and seepage characteristics of unloading rock with and without water pressure. The effects of different water pressure differences on the failure of unloading rock and the law of seepage were investigated by analyzing the change of acoustic emission and permeability coefficient with stress. The results showed that the unloading rock without water pressure was brittle failure, and the initial damage of the unloading model with water pressure was earlier than that of the model without water pressure and showed greater brittleness, and its cracks first break through at the bottom of the sample with higher osmotic pressure. With the increase in unloading, the permeability of rock increased gradually until it appeared an abrupt change. The failure mode and permeability law of the rock with different water pressure differences were basically the same, but the greater the pressure difference, the smaller the effective unloading capacity when the permeability coefficient changes suddenly, and the greater the possibility of water inrush in the rock.

1. Introduction

From the perspective of the nature of mechanics, the excavation and support of underground engineering belong to the category of alternating loading and unloading mechanics. In many cases, the unloading rock is simultaneously affected by groundwater seepage pressure. For example, most roof and floor water inrush failures in underground mining are caused by water pressure and unloading.

In [1], the unloading test of rock under high confining pressure and high water pressure was conducted. The analysis found that the degree of unloading damage after the rock reached the peak strength was greater than that before the peak. The existence of water pressure weakened the unloading strength of the rock in the process of unloading. In [2], the mechanical properties and failure mechanism of rock under loading and unloading were compared and analyzed by RFPA2D. It was found that the unloading failure was mainly tensile failure, and the angle between the failure surface and the direction of the maximum principal stress was less than that of the loading failure surface. In paper [3], the fluid-solid coupling of rock mass with complex fractures was studied, and it was found that the hydraulic gap width of the fractures was negatively correlated with the stress. In the article [4], a triaxial compression test of fractured rock mass under seepage action was carried out, the fluid-solid coupling equation of a single fracture under triaxial conditions was proposed, and the mechanism of stress affecting fracture permeability was revealed. In [5], based on the principle of effective stress, the mechanical mechanism of crack propagation under water pressure was studied. In [6], the UDEC numerical analysis software was used to study the permeability characteristics of fractured rock masses under various stress states, and it was obtained that the law of equivalent permeability coefficient was influenced by the normal closure and the dilatation effect of fractures. In [7], three coupling evaluation models (elastic and plastic strain-seepage coupling, rheology-seepage coupling, and variable parameter rheology-seepage coupling) were proposed, and FLAC3D software was used to analyze the delayed water inrush and its seepage mechanism in the fractured structural zone of the mine. In [8], RFPA2D-flow calculation program was applied to simulate and analyze the whole process of fracture formation and transfixion of the floor (including concealed small faults) under the joint action of mining and high confined water and revealed the water inrush evolution mechanism of the whole floor under mining and water pressure. In [9], the FLUENT module of ANSYS was used to study the nonlinear characteristics of fluid flow and the variation of critical hydraulic gradient with the opening, roughness coefficient, and the number of intersection points of the crack grid. In [10], the Navier–Stokes equation was solved using the COMSOL metaphysics program, and the seepage change characteristics of a single fracture with different shear displacements were calculated, and the change rules of equivalent hydraulic gap width, mechanical gap width, and volume flow velocity associated with the shear displacement of the rock were obtained. It can be seen from the above that previous studies have deepened the understanding of rock permeability and deformation characteristics under the coupled action of the rock seepage field and stress field. However, the influence of water pressure difference on the confining pressure unloading rock still needs to be explored to further understand its regularity.

In this paper, considering the influence of water pressure difference on confining pressure unloading rock, RFPA2D-flow software was used to simulate the confining pressure unloading rock with and without water pressure and the confining pressure unloading rock under different water pressure differences. According to the change of acoustic emission and permeability coefficient with stress, the initiation of new fractures, the expansion of fractures, and the evolution of permeability of seepage rock were studied in the process of unloading. It has important practical significance for the safe construction and technical guidance of underground engineering.

2. Principle of RFPA2D-Flow

The seepage of unloaded rock mass essentially means that the rock is subjected to stress and water pressure. And, the cracks in the rock stratum are generated, expanded, and penetrated until the rock stratum is unstable and damaged. It is a typical fluid-solid coupling process. The porous media deform under the fluid load, which in turn affects the distribution and magnitude of the flow field load.

RFPA2D-flow uses elastic mechanics as a stress analysis tool, elastic damage theory, and modified Mohr–Coulomb failure criterion as a medium deformation and failure analysis method to conduct a coupled analysis of stress, deformation, and seepage evolution during rock failure [11]. In the elastic state, the mesoelement satisfies the relationship between the permeability-stress-strain function, and the permeability of the element increases sharply after the element is damaged and ruptured. Then, through the coupling loop iteration, the coupling calculation between the elementʼs permeability and stress is realized. The classical equation of seepage coupling theory is as follows.

Equilibrium equation is as follows:

Geometric equation is as follows:

Constitutive equation is as follows:

Seepage equation is as follows:

Seepage and stress coupling equation is as follows:

In the formula, and are Lame coefficient and shear modulus, respectively; and are the total stress and effective stress, respectively, and their unit is Pa; xj, uj, and fj are coordinates, displacement, and volume forces in the direction of j; and are total strain and volumetric strain; is the Kronecker constant; k is the permeability coefficient; S is the water storage coefficient; and , , and a are the permeability coefficient jump ratio, coupling coefficient, and pore water pressure coefficient, respectively. The p is the pore water pressure, and the unit is Pa. t is time, and the unit is s.

In the process of stress-seepage coupling calculation, in order to make the rock unit better realize the coupling of the seepage field and stress field, the equivalent permeability coefficient was calculated according to the stress and pore water pressure of the model unit at the end of each cycle, and then the permeability coefficient was reassigned to each element of the model. It should connect the evolution of pores and fractures in rock with the permeability coefficient so that the seepage field and stress field interact with each other [12].

3. Numerical Experiment Model

3.1. Parameters of the Model

According to the principles of the RFPA2D-flow simulation software, a numerical simulation model was established to carry out loading and unloading seepage tests on rock with water pressure. Figure 1 shows the two-dimensional plane stress thin plate model. The width of the model was 50 mm, and the height was 80 mm. The model was divided into 100 × 160 = 16,000 units. The physical and mechanical parameters of the model are shown in Table 1.

The bottom boundary of the model was imposed with a displacement constraint in the vertical direction; an axial pressure p1 was applied on the top, and a confining pressure p2 was applied on all sides, and upper water pressure p3 and lower water pressure p4 are applied on the top and bottom boundaries of the model, respectively. The step-by-step loading method is adopted, the seepage time step is set as 1 day, and the maximum iteration number is set as 300.

3.2. Process of Simulation

Considering that if the confining pressure of the rock is released after the peak, the stress has exceeded the yield strength of the rock, and the specimen is in a state of failure. The Mohr–Coulomb strength criterion was used to simulate the unloading of seepage rock before the peak [13]. The simulation steps were as follows:(1)First of all, we applied axial pressure p1 and confining pressure p2 at the same time and kept them equal at all times. After loading to a predetermined pressure step by step, we applied certain pore water pressures p4 and p3 at the bottom and upper boundaries of the model to realize the rock seepage simulation.(2)The confining pressure p2 was kept unchanged, and displacement control was used to carry out the axial loading on the rock model. The loading displacement increment was ΔS = 0.005 mm until the axial pressure value p1 reached a certain state before rock failure. That was the critical unloading point. The critical value of the unloading stress should be much greater than the compressive strength under uniaxial compression and less than the compressive strength of the three-axis test under the corresponding confining pressure [14]. The starting point of unloading in this simulation was taken as 75% of the maximum strength of three-axis compression.(3)The axial pressure p1 was kept unchanged, and the confining pressure p2 was gradually removed Each unloading stage p = 0.2 MPa until the model was damaged.

3.3. Validation of the Model

According to the above process, the numerical simulation model of the test was established. The model needed to be tested to verify its reliability. In the process of verification, the models with different grid sizes were used for calculation. The material parameters are shown in Table 1. It can be seen from Figure 2 that the larger the grid size of the model, the higher the peak intensity; the smaller the grid size, the lower the peak intensity. When the grid size was less than 0.5 mm, the peak strength tended to be stable. However, the smaller the grid size, the more the number of elements and the longer the time of calculation. Therefore, 0.5 mm was selected as the optimal grid size after comprehensive consideration.

In [15], a cylindrical marble specimen with a diameter of 50 mm and a height of 100 mm was used to carry out the prepeak unloading tests with high confining pressure (20 MPa) under nonwater pressure and high water pressure (10 MPa) conditions. It is found that the existence of high water pressure reduces the initial confining pressure, that is, reduces the effective normal stress on the fracture surface and failure surface of rock, accelerates the fracture of rock, and reduces the strength of rock. Under the condition of high water pressure, the unloading rock shows more brittle failure, and its failure is more severe than that without water pressure, and its strength is significantly lower than that without water pressure.

RFPA-flow software was used to simulate the above experiment. A two-dimensional stress model with a width of 50 mm and a height of 100 mm was established. The grid was divided into 100 × 200 = 20,000 square elements. The physical and mechanical parameters of the model were adopted in reference [15].

The numerical simulation results of the test models with and without water pressure were compared with the results of the experiment in reference [15]. It can be found from Figure 3 that the stress-strain curve of the results of the simulation was approximately consistent with the results of experimental. The peak strength of the prepeak unloading test without water pressure was 148 MPa, and the peak strength of the prepeak unloading test with water pressure was 133 MPa. The peak strength of the simulation results without water pressure was 152 MPa, and the peak strength of the simulation results with water pressure was 128 MPa. There was little difference between the simulation results and the test results. Therefore, according to the verification in this section, the mesh and calculation of this numerical model are feasible.

4. Results of Numerical Simulation

4.1. Simulation Results with and without Water Pressure

First, the prepeak unloading of rocks with or without water pressure was simulated. The rock at a depth of 800 m was selected, and its gravity stress was approximately 20 MPa. The confining pressure was 0.6–0.8 times of the gravity stress. In this simulation, the confining pressure was 15 MPa. The simulation process of unloading rock with water pressure was similar to that without water pressure. On the basis of the model without water pressure, pore water pressure p3 and p4 were applied to the upper and lower boundaries of the model, which were, respectively, 8 MPa and 10 MPa. According to the change of acoustic emission and permeability coefficient with stress, the initiation of new fractures, the expansion of fractures, and the evolution of permeability of seepage rock were studied in the process of unloading.

4.1.1. Changes of Axial Stress and Acoustic Emission and Fracture Mode

By comparing and analyzing Figures 4 and 5, it could be seen that the number of acoustic emissions reached the maximum 289 when the unloading model without water pressure was at step 80, and almost no acoustic emission occurred before this. The ultimate fracture mode of the rock was brittle failure. The fracture started from the lower left corner to the upper right. The model with water pressure is different from the model without water pressure. Acoustic emission occurred before the maximum number of acoustic emissions in the unloading model with water pressure. The initial damage was earlier than the model without water pressure. The number of acoustic emissions reached the maximum 482 when the unloading model with water pressure was at step 115, and it showed more brittleness. Because the water pressure at the bottom boundary of the model was greater than that at the top boundary, the failure of the model mainly occurred in the lower part of the model.

4.1.2. Permeability Coefficient and Seepage Vector

Figures 6 and 7 are the permeability coefficient change curve and seepage vector diagram of the unloading rock with water pressure, respectively. During the loading process, the primary pores and microfractures gradually closed and the permeability coefficient of rock decreased. During the unloading stage, the permeability coefficient increased and even abruptly changed because of the continuous initiation, expansion, and acceleration of rock microfractures. The seepage path became increasingly obvious, and the seepage vectors were concentrated around the fracture zone.

4.2. Simulation Results of Different Water Pressure Differences

In order to explore the failure process of confining pressure unloading rocks with different water pressure differences, the following simulation schemes were made, as shown in Table 2. The physical and mechanical parameters and simulation steps of the rock were consistent with the previous simulation process with or without water pressure.

According to the numerical simulation calculation results of confining pressure unloading rock with different water pressure differences, the strength characteristics and failure evolution of confining pressure unloading rocks with different water pressure differences are analyzed from three perspectives, including the models of rock fracture, the changes of axial stress and acoustic emission and permeability coefficient during the experiment, and relationship between unloading and permeability coefficient.

4.2.1. Models of Rock Fracture

Figures 8 and 9 are the final fracture and seepage vector of rock specimens with different water pressure differences. It can be seen from Figure 8 that with the increase in water pressure difference, the range and degree of rock fracture increase continuously. The seepage vector distribution of Δp = 2 MPa in Figure 9 is relatively uniform, indicating that the fracture is not too large. With the increase in water pressure difference, the more seepage vectors are generated near the fracture zone. This is in accordance with the law of optimizing the seepage path, that is, the velocity of macrofractures is larger than that of microfractures [16].

4.2.2. Changes of Axial Stress and Acoustic Emission and Permeability Coefficient

Figures 10 and 11 are the AE and axial stress and permeability coefficient of rock specimens with different water pressure differences. It can be seen from Figure 10 that during the initial loading stage, the rock with different water pressure difference has no acoustic emission, and the rock has not been damaged; as the calculation proceeds, small failure points and acoustic emission appear. When the acoustic emission reaches the maximum, it indicates that the rock has been greatly damaged. When the water pressure difference is 2 MPa, 4 MPa, and 6 MPa and the corresponding load step is 117, 97, and 95, the acoustic emission number reaches its maximum value. It shows that the number of running steps to make rock failure decreases with the increase in water pressure difference. It can be seen from the change process of permeability coefficient in Figure 11 that rock failure modes with different pressure differences are basically same; during the loading process, the primary pores and microcracks gradually closed and the permeability decreases; in the unloading stage, the fractures constantly sprout, expand, and penetrate and the permeability increases continuously and even jumps.

4.2.3. Relationship between Unloading and Permeability Coefficient

To further describe the relationship between unloading volume and permeability coefficient, the unloading ratio of effective stress is introduced. It is the ratio of the difference between the initial stress and the effective stress corresponding to a load step in the unloading process and the initial stress, which is given as follows:

In the formula, is the unloading ratio of effective stress, is the initial stress, and is the stress value corresponding to a load step in the unloading process.

Figure 12 shows the relationship between the permeability coefficient and the unloading when the water pressure difference is 2 MPa, 4 MPa, and 6 MPa. In the unloading elastic stage, the permeability coefficient increases little and changes little. When entering the plastic stage, the reduction of confining pressure increases the rock fracture and the brittleness of the rock. The deformation of rock changes from compacting deformation to dilating deformation. Percolation transitioned from the percolation of compressed pores to the percolation of shear fractures. The permeability coefficient increased greatly. The penetration of unloading rock with different water pressure differences is also quite different. When the water pressure difference is 2 MPa, 4 MPa, and 6 MPa and the corresponding unloading capacity is 68.15%, 50.66%, and 45.58%, the permeability has obviously increased. It shows that with the increase in water pressure difference, the unloading quantity to make rock failure becomes smaller and smaller. In addition, when p4 is 6 MPa, 8 MPa, and 10 MPa, the permeability coefficient is basically 0.0178 m/d at the initial stage of unloading, and the permeability coefficient after failure is 0.149 m/d, 0.1512 m/d, and 0.2343 m/d, respectively. It is 8.37 times, 8.51 times, and 13.16 times the original. It indicates that with the increase in water pressure difference, the change of permeability coefficient increases with the failure of rock. It shows that the permeability has changed greatly after rock failure. The water flow channel is pores composed of compacted particles, which are transformed into continuous cracks formed by shear failure.

5. Conclusions

According to the comparative analysis of the simulation results, the following conclusions can be obtained:(1)The unloading rock without water pressure is a brittle failure. The initial damage of the unloading model with water pressure is earlier than that of the model without water pressure. The water pressure at the bottom boundary of the model is greater than that at the top boundary, so the failure of the model mainly occurs at the bottom.(2)With the increase in water pressure difference, the extent and degree of rock fracture increase. The increase in water pressure difference makes the seepage path gradually obvious, and a larger seepage vector is generated near the fracture zone.(3)The rock failure modes under different water pressure differences are basically identical; in the loading process, the primary pores and microfractures gradually close and the permeability of rock decreases, and in the unloading stage, the fractures constantly sprout, expand, and penetrate and the permeability increases and even jumps.(4)In the elastic stage of unloading, the permeability coefficient increases little. When entering the plastic stage, the reduction of confining pressure causes the rock deformation to change from compacted deformation to dilated deformation, the opening of the fracture increases, and the corresponding seepage transition from compressed pore seepage to shear fracture seepage and the permeability coefficient increase greatly. With the increase in water pressure difference, the amount of unloading when the failure is reached becomes smaller and smaller.

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 they have no conflicts of interest.

Acknowledgments

This research was supported by the First-Rate Discipline Construction Program of Mining Engineering (01LX03303).