#### Abstract

In this paper, the transient pulse test is used to study the permeability and hydromechanical coupling effect of the fractured limestone. The permeability parameters (permeability, *β* factor of non-Darcy flow, and acceleration coefficient) of non-Darcy flow in fractured limestone are obtained by experimental data. The experimental results show that, in the process of transient seepage test of fractured Maokou limestone, the relationship between hydraulic pressure gradient and seepage velocity does not conform to Darcy’s law but meets Forchimer relationship. The relationship between hydraulic pressure difference and time can be fitted by quartic polynomial. The larger the confining pressure is, the more obvious the non-Darcy seepage effect of fractured rock seepage is. The seepage of rock fracture under high confining pressure is a highly nonlinear time-varying seepage. The permeability coefficient of rock decreases with the increase of volume stress. Under the action of low volume stress, the relationship between permeability coefficient and stress is more sensitive, while under the action of high volume stress, the relationship between permeability coefficient and volume stress is not significant. In the process of volume stress increasing, the *β* factor of non-Darcy flow appears negative. Under the action of low volume stress, the acceleration coefficient and *β* factor of non-Darcy flow increase, while under the action of high volume stress, the acceleration coefficient and *β* factor of non-Darcy flow decrease.

#### 1. Introduction

Fractured rock mass, which is potentially dangerous, is widely existing in nature. The seepage instability of fractured rock mass is an important reason for rock engineering destruction and even large-scale geological disasters. Therefore, the hydraulic coupling and seepage characteristics of fractured rock masses are generally concerned [1–6]. Consequently, experts and scholars have proposed a variety of mathematical models of seepage. For example, Dverstorp and Andersson [7] use a discrete fracture network model to conduct seepage analysis of fractured rock mass, and Zhao [8] has proposed a dual-medium multifield coupling model. And many experiments of fluid-solid coupling have been done [9–15].

The rock permeability is usually studied by the steady-state method or the transient method. Miao et al. used steady-state method to test the seepage performance of broken sandstone under different porosity conditions [16–19]. Wang et al. measured the rock seepage rate by transient method and analyzed the hydrodynamic characteristics of transient method in rock permeability test by using the hydrodynamic theory in porous media [20, 21]. Hsieh et al. found that the transient method is suitable for low permeable rocks and tested the permeability of a variety of different geological materials with the transient method [22–27].

It is found that, under the condition of low water pressure, the coarse fractured rock mass has Darcy flow characteristics, and under the condition of high hydraulic pressure, the seepage of rock mass has Darcy flow characteristics, and the hydraulic pressure plays an important role in the development of rock fracture [28–32]. Li et al. [33] discovered the variation of rock permeability with rock porosity by studying rocks with different particle sizes. Zhao et al. [34] found that the seepage curve of low permeability rock has obvious non-Darcy seepage characteristics through transient pulse method and used a new data analysis method to study the relationship between velocity and hydraulic gradient. Nguyen et al. [35] studied the effect of stress path on the mechanical properties and coupled permeability evolution of quartz sand filling and soft sandstone. The seepage of the fractured rock mass is controlled by the shape of the fracture, so many models and parameters have been proposed to characterize the fracture [36–43]. In addition, factors affecting rock permeability include rock type, loading method, and damage history [44–50].

It can be seen from overview of the references that although there are many researches on the seepage of fractured rock mass, these researches mainly focus on the use of steady-state and transient methods to study the factors affecting the seepage of rock, the extraction of non-Darcy flow permeability coefficient mainly focuses on permeability, and the research on the extraction of acceleration coefficient and non-Darcy flow *β* factor is less. In this paper, the transient pulse test is carried out on the fractured limestone produced by splitting method. The test data of the hydraulic pressure varying with time of fractured limestone under the action of constant axial pressure and variable confining pressure is obtained by MTS815 [34, 51], the permeability coefficient of non-Darcy flow is extracted, and its permeability characteristics and hydromechanical coupling effect are studied, and the seepage law of fractured Maokou limestone under different confining pressures is revealed. According to the data about permeability coefficient, volume stress, and hydraulic pressure difference, the law of change among them can be obtained. It has a certain guiding effect on the seepage problem of fractured rock mass in practical engineering.

#### 2. Test Method

##### 2.1. Specimen Preparation

In this paper, Maokou limestone from Ningxiang coal mine of Hunan Province was used to make a sample of Φ 50 × 100 mm. The limestone in Maokou is very dense, so it is necessary to induce some fractures artificially for seepage experiment. The artificial fractures were obtained by Brazil splitting method. Firstly, the standard rock sample of Φ 50 × 100 mm with no damage on the edge and smooth end face was selected. Then a layer of glue on the outer surface of the rock sample was evenly applied, and a layer of thermoplastic tube with the same height as the rock sample was covered. Next the heat gun was used to close the thermoplastic tube from top to bottom. Finally, the Brazilian splitting method was used to conduct the splitting test on the rock sample. In order to prevent the instant brittle failure of the test sample, which leads to the failure of rock sample production, the loading rate of the test is strictly controlled to 10 N/s until the rock fracture occurs; the test is stopped immediately.

Figure 1 is the induced fracture diagram of Brazil splitting method, and Figure 2 is the fracture end face photo and fracture path sketch of the test piece.

The sample fracture was marked with a pen at different parts, especially at the places where the gap width is obviously different, so that the measured value can represent the gap width of the fracture more accurately. The fracture of the end face of the rock sample at the serial number was measured with a fracture microscope, each reading was recorded, and the average value of all the measured values is the gap width of the fracture. Table 1 is the measurement value of initial crack width of rock sample terminal face.

##### 2.2. Seepage Experiments

The test is carried out on MTS815 rock mechanics multifunctional testing machine, which is produced by MTS company of the United States and is specially used for testing rock and concrete and other materials. It has four independent control systems, which are axial force, confining pressure, hydraulic pressure, and temperature, and it is at international leading level in testing rock and concrete materials. It is mainly used for various dynamic tests and static tests of rock and concrete materials, including uniaxial compression test, uniaxial tensile test, fracture strength, triaxial pressure, osmotic pressure, hydraulic pressure, temperature simulation, and other rock-related physical property analysis applications, which provides the test means for the research and analysis of dynamic characteristics of these materials and structures. It is the most advanced indoor rock mechanics test equipment (Figure 3). The main technical parameters of MTS815 are as follows: the axial force system: the overall rigidity of testing machine frame is 10.5 × 109 N/m, maximum axial force can be loaded to 4600 kN, and maximum axial tension can be loaded to 2300 kN; the confining pressure system: the maximum confining pressure can be loaded to 140 MP; the hydraulic pressure system: the maximum hydraulic pressure can be loaded to 140 MPa, and the pressure difference and seepage flow can be measured; the temperature system: it can provide a temperature between 20°C and 200°C.

The seepage test of rock fracture under cyclic loading and unloading was carried out in the triaxial pressure chamber of MTS815. During the test, when the stress reaches the preset value, the hydraulic pressure loading system is used to apply the same pressure on both ends of the rock sample and then reduce the hydraulic pressure at one end. At this way, the hydraulic pressure difference between the two ends of the rock sample is formed, collecting the series of the hydraulic pressure difference of the fracture varying with time and calculating the permeability characteristics. Table 2 shows the stress and hydraulic pressure used in the test. Figure 4 is the flowchart of seepage test of rock fracture under cyclic loading and unloading.

#### 3. Permeability Parameter Calculation

There are two methods to determine the permeability of rock samples: steady-state method and transient method. The transient method is to record the change value of the hydraulic pressure gradient with time in the process of rock seepage and calculate the change rate of the hydraulic pressure gradient. The permeability characteristics of rock can be obtained through the fitting curve of the hydraulic pressure gradient and its rate of change.

Figure 5 shows the test principle of permeability characteristics. The volume of both chambers in MTS hydraulic pressure system is B, and their hydraulic pressures are *P*1 and *P*2, respectively. The cross-sectional area and height of rock sample are A and H, respectively. Since the hydraulic pressure at both ends of rock sample is different at the initial time (), there will be a hydraulic pressure gradient . The liquid in chamber 1 will enter into chamber 2 through the rock sample, which will cause the hydraulic pressure in chamber 1 to decrease continuously, while the hydraulic pressure in chamber 2 will rise continuously until the hydraulic pressure in the two chambers is equal, so as to reach the equilibrium state. Suppose that the mass flow of the water penetrating the rock sample from chamber 1 is *q*, and the rock sample is saturated with water. Then the mass flow of water from the rock sample into the chamber 2 is *q* as well, and the seepage velocity of water in the rock sample is . From the compressibility of the fluid, we can obtain

In (1), is the compression coefficient and *ρ* is the density.

According to the formulaswe can obtain

From (4) and (5), we can getwhere is the hydraulic pressure gradient of rock sample, .

##### 3.1. Calculation of Darcy Flow Permeability Parameters

For Darcy flow, the permeability velocity and pressure gradient follow the following formula:where is the dynamic viscosity of the permeable liquid and is the permeability of Darcy flow of the rock sample. Substituting (7) into (8),

In the experiment, sampling at equal interval *t* in the test, the total number of samplings is *n*. The end of sampling is . The hydraulic pressure gradient at the end of sampling is . Integrating (9), we can getwhere the hydraulic pressure gradients and are negative and is positive. Therefore, is meaningful. The permeability of rock sample can be calculated from (10). One has

Equation (11) is the calculation formula of rock sample permeability obtained by transient method which is using MTS815 at present.

##### 3.2. Calculation of Non-Darcy Permeability Parameters

The experimental study shows that, in many cases, the seepage flow of rock sample is non-Darcy flow, and its seepage law satisfies the Forchimer relation; that is,where is acceleration coefficient, is dynamic viscosity, is permeability, and is non-Darcy flow factor.

From (7), we get

It is assumed that the time series of fracture hydraulic pressure difference collected by equal time interval *t* during the test is (*i* = 1, 2, 3, …, *n*). In this way, the sequence of hydraulic pressure gradient (*i* = 1, 2, 3, …, *n*) can be calculated. The formula of time series for calculating the seepage velocity *V* and change rate can be obtained by the difference between (7) and (13):

From (12),

Constructive functional is as follows:

The extremum conditions of the functional are as follows:

By combining (18) ∼ (20), we can get , , and .

Sometimes, the single time series of permeability velocity and its change rate calculated by difference is not smooth enough. In this case, polynomial can be used to fit the curve of the hydraulic pressure gradient varying with time scatter diagram, (7) and (13) can be derived to obtain a single time series of permeability velocity *V*, and its change rate in (18) ∼ (20) becomes

#### 4. Analysis of Test Results

##### 4.1. Non-Darcy Effect of Rock Fracture Seepage

The transient pulse penetration test was carried out on Maokou limestone. The percolation liquid was distilled water, its mass density is 1000 kg/m^{3}, its dynamic viscosity is 1.01 × 10–3 Pa·s, and its compression coefficient is 0.556 × 10^{–9} Pa^{−1}. The volume of two chambers in MTS hydraulic pressure system is 0.336 × 10^{–6} m^{3}.

Figure 6 shows the decay curve of hydraulic pressure difference of M07 with time under an axial pressure of 6 MPa and different confining pressures. From Figure 6, we can find that the hydraulic pressure difference decreases rapidly in the early stage, slows down in the later stage, and finally approaches plateau. This is due to the decrease in the hydraulic pressure difference, which reduces the water seepage velocity. With the increase of confining pressure, the time required for the hydraulic pressure difference to reach equilibrium state increases significantly. This is due to the increase of confining pressure, which closes the rock fissures and creates greater resistance to water seepage. Using (11) to fit the test curve, the fitting curve is also shown in each figure. From each figure, it is found that the curve drawn by exponential formula (11) derived from Darcy’s law is quite different from the test data, and with the increase of confining pressure, the deviation between fitting curve and test data is greater.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

The hydraulic pressure difference decay-time curve can be fitted by the quartic polynomial, and the fitting has a very high fitting precision. Table 3 is the quartic polynomial fitting equation of hydraulic pressure difference varying with time of M07 under different confining pressures.

This shows that the flow of fractured rock is a non-Darcy flow, which is time-varying. With the increase of confining pressure, the non-Darcy seepage effect of fractured rock is more obvious. The non-Darcy effect of fractured rock seepage is affected by confining pressure. When the confining pressure is low, the seepage of rock fracture is close to Darcy flow, and the seepage of rock fracture under high confining pressure is significant non-Darcy flow.

The volume pressure can be expressed as .

The fitting equation under axial pressure of 6 MPa and confining pressure of 8 MPa is

Hydraulic pressure gradient can be expressed as

The seepage velocity yields

According to the fitting equation, there is no linear relationship between the hydraulic pressure gradient and the seepage velocity, indicating that the seepage of Maokou limestone does not obey Darcy flow.

According to the data collected from the experiment, the permeability characteristic parameters of non-Darcy flow are calculated.

##### 4.2. Seepage Characteristics under the Action of Hydromechanical Coupling

From Figure 7, it can be observed that the permeability of rock sample M07 decreases with the increase of volume stress. Under the action of low volume stress, the permeability change is more sensitive, and under the action of high volume stress, the permeability change is slow, because, at the initial stage of volume stress loading, the rock fracture has closed to a large extent, which has great resistance to the seepage of water and reduces the permeability. With the increase of volume stress, the rock fracture is further closed. When the fracture opening is close to the residual gap width, further increase of volume stress will maintain the residual gap width rather than close the rock fracture further, resulting in the fact that the permeability of fractured rock has little correlation with stress under the condition of high volume stress.

The *β* factor of non-Darcy is a parameter to indicate the degree of seepage deviating from Darcy flow, which is related to seepage coefficient and porosity. When *β* factor of non-Darcy is 0, the seepage flow is in Darcy flow state. When *β* factor of non-Darcy is negative, the seepage is in a state of instability. From Figure 8, it is shown that when the volume stress increases from 14 MPa to 34 MPa, the *β* factor of non-Darcy flow fluctuates greatly. When the volume stress is 14 MPa, the value of *β* factor is very small, which means that the seepage is close to Darcy flow. When the volume stress is 26 MPa, the *β* factor is at the maximum value, and the seepage deviates from Darcy flow to the maximum extent. When the volume stress is 30 MPa, the *β* factor is in a negative value, indicating that the seepage is in a state of instability. When the volume stress is 34 MPa, the value of *β* factor decreases to a lower value, and the degree of seepage deviation from Darcy flow is lower. It is noted that, in the early stage of the volume stress loading, the *β* factor generally increases, which is due to the closure of the rock fracture under the volume stress, making the seepage deviate from the Darcy flow. In the later stage of volume stress loading, the fracture and pore passage in the rock have complex changes, and the degree of seepage deviating from Darcy flow is reduced.

Acceleration coefficient *c*_{a} is an index reflecting the seepage inertia. The greater the acceleration coefficient is, the more difficult it is to change the flow state. The process of crack and fracture expansion and porosity increase is the process of inertia weakening, that is, the process of acceleration coefficient decreasing, as well as the process of transformation from steady flow to unsteady flow. From Figure 9, it can be observed that when the acceleration coefficient *ca* is always increasing, the volume stress is loaded from 14 MPa to 22 MPa. When the volume stress is from 26 MPa to 34 MPa, the value of acceleration coefficient *ca* is small, with small variation. It is noted that, under the action of low volume stress, rock fracture is gradually compacted and acceleration coefficient *ca* increases gradually, while under the action of high volume stress, complex changes have taken place in the internal fractures and pore channels of rocks and the acceleration coefficient *ca* decrease.

##### 4.3. Seepage Law under the Combined Action of Volume Stress and Hydraulic Pressure Difference

By analyzing the test data of non-Darcy flow obtained from the seepage test of fractured rock, the fitting relationship between the permeability coefficient of fractured rock, volume stress, and hydraulic pressure difference can be obtained based on the transient pulse methodwhere *a*, *b*, and c are the fitting parameters; is volume stress, unit: MPa; is the permeability coefficient, unit: m/s; and is the hydraulic pressure difference, unit: MPa.

Equation (28) is used to fit the test data, and the equation of fluid-solid coupling of fractured rock permeability coefficient is obtained. The fitting curve is shown in Figure 10.

From the analysis of the test data, we can draw the conclusions that, under the same axial pressure, the permeability coefficient of rock sample decreases continuously with the volume stress increases. Under the action of low volume stress, the permeability coefficient change is more sensitive, and under the action of high volume stress, the permeability coefficient change is slow. This is because, at the initial stage of confining pressure loading, the rock fracture has closed to a large extent, which has great resistance to the seepage of water and reduces the permeability coefficient. With the increase of confining pressure, the rock fracture is further closed. When the fracture opening is close to the residual gap width, further increase of confining pressure will maintain the residual gap width rather than close the rock fracture further, resulting in the fact that the permeability coefficient of rock fracture has little correlation with stress under high volume stress.

#### 5. Conclusions

The transient seepage test of fractured Maokou limestone is carried out based on MTS815 rock mechanics test system. The non-Darcy seepage characteristics of fractured rock under confining pressure are extracted by a single time series of hydraulic pressure difference to analyze the relationship among seepage characteristics, stress. Through the experimental study, the main conclusions of this paper are as follows:(1)Under the action of confining pressure, the hydraulic pressure difference of fractured rock is not in accordance with the specific exponential decay law. There is a significant difference between fractured rock seepage and Darcy flow. In the process of fractured rock transient pulse seepage test, the hydraulic pressure gradient and seepage velocity meet the Forchimer relationship rather than Darcy law. The equation of hydraulic pressure difference varying with time can be fitted by quartic polynomial.(2)With the confining pressure increase, the non-Darcy seepage effect of rock fracture seepage is more obvious. And the seepage of rock fracture which under high confining pressure is a kind of highly nonlinear time-varying seepage.(3)The permeability coefficient of fractured rock is affected by volume stress; the permeability coefficient decreases with the increase of volume stress. The relationship between permeability coefficient and stress is more sensitive while being under the action of low volume stress and it is not significant while being under the action of high volume stress.(4)During the process of volume stress increasing, the *β* factor of non-Darcy flow appears negative. Under the action of low volume stress, the acceleration coefficient and *β* factor of non-Darcy flow increase, while under the action of high volume stress, the acceleration coefficient and *β* factor of non-Darcy flow decrease.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research was supported by the National Natural Science Foundation of China (nos. 51774131, 51274097, and 51434006), the CRSRI Open Research Program (CKWV2017508/KY), and the Open Projects of State Key Laboratory of Coal Resources and Safe Mining, CUMT (no. SKLCRSM16KF12).