Abstract
To study the effects of compression rate and weight ratio of sandstone to mudstone on permeability, water flow tests were performed on crushed rocks using an in-house designed and built water flow apparatus. Seepage properties of sandstone and mudstone were tested on specimens under six axial displacement levels. The weight ratios of the six specimens of sandstone to mudstone were set to 1 : 0, 4 : 1, 3 : 2, 2 : 3, 1 : 4, and 0 : 1. A non-Darcy law was adopted to calculate the permeability in the experiments based on verification of the Reynolds number. Non-Darcy flow is found to be most significant when the axial displacement or mudstone content increases. The permeability in the experiments is found to be in the range of 9.1 × 10−15–9.492 × 10−13 m2, which is calculated from a quadratic polynomial fit between the flow velocity and pressure gradient. A compression rate of 28.6% and weight ratio of sandstone to mudstone of 60% are found to be key values affecting permeability. When the axial displacement of specimen A (1 : 0) increases from 10 mm to 30 mm, the compression rate (ratio of axial displacement to original specimen height) increases from 9.5% to 28.6%, and the permeability decreases by 83.8% to 1.534 × 10−13 m2. When the axial displacement is 10 mm and the content of mudstone is increased from 0% to 60%, the permeability decreases by 77.1% to 2.172 × 10−13 m2.
1. Introduction
Underground mining can lead to strata movement and loss of equilibrium of the stress state, causing deterioration in the mechanical behavior of surrounding rocks [1–3]. Schatzel et al. [4] tested water flow properties of overburden strata at a mine site and demonstrated that the permeability increased by two to three orders of magnitude as a result of mining disturbance. Advancement of the working face may result in a series of hazards, such as groundwater inrush into a coal mine from an aquifer, which is one of the most frequent forms of mining accidents [5–7]. This type of accident seriously threatens safe production of the mine, and mine workers themselves [8, 9]. Since 2000, about 500 events of mine water inrush have occurred in China and more than 3000 people have lost their lives. For example, 36 people died in a groundwater inrush accident at the Wangjialing Coal Mine in Shanxi Province, China, on 28 March 2010 [8]. In recent years, and with continuing extension to more complex coal seams containing faults, karst collapse pillars and other weak tectonic structures, the hydrogeological conditions in such mines are deteriorating [10, 11].
Many scholars have studied the mechanics and the prediction and the control measures of water inrush [12–15]. Bukowski [16] proposed a risk assessment method to classify mine shafts according to the risk of water inrush. The method considers factors such as flow intensity, condition of water-bearing formations, suspended materials contained in the water, and protective coal pillars. Wu et al. [11] adopted the “three maps-two predictions” method to construct an index model of water abundance in overlying aquifers, with results shown on a map indicating potential water flow channels to aquifers.
A more accurate and in-depth understanding of seepage properties is important for prediction and control of mine water inrush accidents [17], and the study of permeability of rocks is the foundation of such. Tanikawa and Shimamoto [18] carried out water permeability tests on sedimentary rocks from the western foothills of Taiwan and found that the permeability increased slightly when the pressure gradient increased. Benavente et al. [19] studied the water permeability of 15 sedimentary rock types and proposed a model for calculating permeability that considers factors such as flow velocity, water properties (density and viscosity), porosity, interfacial tension, and contact angle.
Research into the permeability of rocks mainly focuses on intact rocks that have not been broken or damaged, while crushed rocks are widely distributed in fault zones, karst collapse pillars, and caved zones. These crushed rocks, of which the permeability is much higher than intact rock [20–25], may act as preferential migration pathways for groundwater into a mine [26].
In this paper, the effects of compression rate and weight ratio of sandstone to mudstone on permeability are investigated using an in-house designed and built water flow apparatus. Water flow tests were performed on crushed rocks containing sandstone and mudstone. The experimental results provide an improved understanding of the mechanism of water inrush through such preferential pathway materials and serve as a reference for prediction and control of water inrush accidents.
2. Experimental Materials and Testing Procedures
2.1. Experimental Materials
The sandstone and mudstone used in the experiments were collected from the 21105 working face of the Pingshuo Coal Mine in Shanxi Province, China. The mineral composition of the tested sandstone and mudstone was composed of feldspar, quartz, kaolinite, chlorite, siderite, and small amounts of other minerals. The sandstone and mudstone blocks were crushed using a hammer, piercer, and rock crusher in the laboratory. A sieve was used to separate the particles in the 10–15 mm size range, as shown in Figure 1. This size range closely correlates with site conditions, as 10–15 mm size sandstone and mudstone grains are typically widely distributed in fault zones, collapse pillars, and caved zones. Six specimens were made, each with a different weight ratio of sandstone to mudstone. These specimen were A (1 : 0), B (4 : 1), C (3 : 2), D (2 : 3), E (1 : 4), and F (0 : 1), the ratios of sandstone to mudstone shown in parentheses.

2.2. Experimental Equipment and Testing Procedures
The experimental apparatus is illustrated in Figure 2. The details of the apparatus setup are provided in the figure caption.

The testing procedures are conducted in two main steps, as described in detail as follows.(1)The axial displacement (S) is applied. Six axial displacement levels were set: 10 mm, 15 mm, 20 mm, 25 mm, 30 mm, and 35 mm. The target value was applied to the specimen and maintained until step (2) was completed, as shown in Figure 3. The porosity (φ) of the crushed rocks was calculated as follows:where A is the cross sectional area of the inner wall of the cylindrical barrel (1.27 × 10−2 m2), m is the mass (1.7 kg), ρ is the density, and Hn is the height of specimen at a given axial displacement and is calculated as follows:where H is the height of the cylindrical barrel (180 mm), hn is the length of part of piston that is outside of the cylindrical barrel (variable according to the six different axial displacement levels), h2 is the thickness of the perforated plate (9 mm), and h3 is the thickness of the felt pad (3 mm).(2)Water flow testing is undertaken. When the axial displacement becomes stable, a pore pressure control system is used to produce four different piston velocities: 10, 20, 30, and 35 mm/min. The relation between water flow velocity (Vs) and piston velocity (Vp) is written as follows:where ds and dp are the respective diameters of the specimen and the piston (12.7 cm and 5.5 cm, respectively). Thus, it is possible to obtain the pressure gradient and corresponding water flow velocities.

3. Permeability Calculation Based on a Non-Darcy Flow
The Forchheimer equation has been proven to accurately quantify water flow in crushed rock when water flow is stable [27, 28]. The value of the pressure gradient relative to flow velocity can be expressed as follows:where μ is the dynamic viscosity (tap water was used in the experiment: μ = 1.01 × 10−9 MPa·s), k is the permeability, ρw is the mass density of the fluid medium (1.0 × 10−6 kg/mm3), and β is a factor of non-Darcy seepage.
The non-Darcy equation (4) adopted in the experiments is verified by the Reynolds number Re, which is calculated as follows:where d is the average grain size of the crushed rocks (10–15 mm). Earlier studies have shown that the upper limit of Re is 5 when using Darcy law, and non-Darcy flow can be adopted when Re > 5. In our experiments, the minimum and maximum values of Re were 0.72 and 10.09, respectively. Thus, non-Darcy flow is more applicable.
There is a difference between the base pressure Pb and top pressure Pt of the specimen, and the pressure gradient, shown in (6), is approximately equal to the ratio of pressure difference to specimen height h. The downstream end is open to the atmosphere in the testing system, thus Pt = 0:
The permeability k is calculated from the quadratic polynomial fit between the flow velocity and the corresponding pressure gradient as follows:
4. Test Results and Discussion
As the water flow tests began, the pore pressure gradually increased, reaching a stable state after about 20 seconds (Figure 4). It was found that a piston velocity of 10 mm/min corresponds to the first axial displacement level in the water flow test, and the 10 mm/min curves show more significant changes than other curves as time increases.

(a)

(b)

(c)

(d)

(e)

(f)
The positions and distributions of the constituent grains changed as the tests continued, leading to the unstable state of the curves. The 10 mm/min curve in Figure 4(a) and the 20 mm/min curve in Figure 4(c) both show a significant decrease before the water flow stabilized. This resulted from the evolution of water flow channels. When the water flow tests were started, the channels through which water flows had not formed completely, and pore pressure was increasing sharply. As water flow continued, the channels became smoother, leading to a decrease in pore water pressure.
Values of water flow velocity (Vs), the corresponding pressure gradient, and the calculated permeability k are provided in Table 1, and the relationship between the pressure gradient and the flow velocity is shown in Figures 5 and 6. As the axial displacement and the content of mudstone increase, the curves of pressure gradient vs. flow velocity change from straight lines to curves. That is, the phenomenon of non-Darcy flow is most significant when the porosity decreases. The relationship between the pressure gradient and the axial displacement is shown in Figure 7. The permeability k was calculated from the quadratic polynomial fitting between the pressure gradient and flow velocity (Vs).



Table 2 shows that the permeability of crushed rocks in the tests was in the range of 9.1 × 10−15–9.492 × 10−13 m2, which is an approximate two orders of magnitude increase compared to fractured sandstone [29].
Figure 8 demonstrates that the permeability is closely correlated with the compression state and, when the axial displacement increases, the porosity decreases, resulting in a decrease in permeability. The decreasing rate is most significant when the axial displacement increased from 10 mm to 30 mm. As the axial displacement of specimen A (1 : 0) increased from 10 mm to 30 mm, the compression rate (ratio of axial displacement to original specimen height) increased from 9.5% to 28.6%, the porosity decreased from 0.427 to 0.275, and the permeability decreased by 83.8% from 9.492 × 10−13 to 1.534 × 10−13 m2. When the axial displacement increased to 35 mm, the permeability decreased slightly and the difference in permeability among different specimens decreased when the axial displacement increased.

For specimens with the same axial displacement, a higher content of mudstone led to a decrease in porosity and thus a decrease in permeability (Figure 9). When the content of mudstone increased from 0% to 60%, the changing rate of permeability was most significant. When the axial displacement was 10 mm and the content of mudstone increased from 0% to 60%, the porosity decreased from 0.427 to 0.400 and the permeability decreased by 77.1% from 9.492 × 10−13 to 2.172 × 10−13 m2. The decreasing rate of permeability declined when the content of mudstone was >60%. As the mudstone content increased, the difference in permeability among different axial displacement levels decreased.

5. Conclusions
Based on the work presented, the following conclusions can be drawn regarding the experiment and the results thereof:(1)Stable water flow in the flow tests is achieved after about 20 seconds, and the pore water pressure showed most significant changes at the first velocity levels for each axial displacement.(2)Non-Darcy flow was adopted for calculations of the permeability, based on the verification of the Reynolds number. Non-Darcy flow was most significant when the axial displacement or the mudstone content increased.(3)Permeability of the crushed rock samples was in the 9.1 × 10−15–9.492 × 10−13 m2 range, which correlates closely with the level of axial displacement and mudstone content.(4)The permeability decreases most significantly when the axial displacement increased from 10 mm to 30 mm. When the axial displacement of specimen A (1 : 0) increased from 10 mm to 30 mm, the compression rate increased from 9.5% to 28.6%, the porosity decreased from 0.427 to 0.275, and the permeability decreased by 83.8% from 9.492 × 10−13 to 1.534 × 10−13 m2. When the mudstone content was increased from 0% to 60%, the changing rate of permeability was most significant. When the axial displacement was 10 mm and the mudstone content increased from 0% to 60%, the porosity decreased from 0.427 to 0.400, and the permeability decreased by 77.1% from 9.492 × 10−13 to 2.172 × 10−13 m2.
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 work was supported by the Initial Scientific Research Fund of Young Teachers in Dongguan University of Technology (grant no. GC300501-069), the Social Science and Technology Development Program of Dongguan (major, grant nos. 2017507140429 and 20185071401603), the Guangdong Natural Science Foundation (grant no. 2018A030313864), the Research on Key Problems of Rail Flaw Detection Based on Ultrasound Guided Waves (grant no. 2016KZDXM053), and the Undergraduate Training Program for Innovation and Entrepreneurship (grant nos. 201811819109, 201811819382, and 201911819028).