Influence of particle size on non-Darcy seepage of water and sediment in fractured rock

Surface water, groundwater and sand can flow into mine goaf through the fractured rock, which often leads to water inrush and quicksand movement. It is important to study the mechanical properties of water and sand in excavations sites under different conditions and the influencing factors of the water and sand seepage system. The viscosity of water–sand mixtures under different particle sizes, different concentration was tested based on the relationship between the shear strain rate and the surface viscosity. Using the self-designed seepage circuit, we tested permeability of water and sand in fractured rock. The results showed that (1) effective fluidity is in 10−8–10−5 mn+2 s2−n/kg, while the non-Darcy coefficient ranges from 105 to 108 m−1 with the change of particle size of sand; (2) effective fluidity decreases as the particle size of sand increased; (3) the non-Darcy coefficient ranges from 105 to 108 m−1 depending on particle size and showed contrary results. Moreover, the relationship between effective fluidity and the particle size of sand is fitted by the exponential function. The relationship between the non-Darcy coefficient and the particle size of sand is also fitted by the exponential function.

water and sand. Yang (2009), Yang et al. (2012) and Sui et al. (2007) analyzed the angle of fluid using cemented sand to analyze the mechanisms supporting the inrushing of water and sand. The flow law was examined during various conditions and critical hydraulic gradients of sand inrush currents. Sui et al. (2008) and Xu et al. (2012) analyzed the initial position of inrushing sand based on the structure of water inrush.
Based on underground water dynamic theories, Zhang et al. (2006) created the critical condition and forecasting formula for the prevention of sand inrush by calculating the hydraulic head. Wu (2004) designed a mechanical model of sand inrush pseudo structures, and discussed the force during sand inrush and described the theory of expression of sand inrush. Zhang et al. (2015a) used a case study to discuss drills resulting in sand inrush based on the funnel model. Zhang et al. (2015b) studied the relationship between backfill and water through conducting crack zone. Moreover, river sediment engineering, the theory of sediment transmission and sediment transport mechanics are excellent subject matters to aid in studying the start and movement of sand in mines. Furthermore, the study of sediment engineering, sediment transport theory and practice, and sediment kinematics can aid in understanding the commencement, flow and inrushing sand problem (Du 2014). But others discussed water and sand form the pressure, water and sand flow in tunnel or broken rock (Limin et al. 2016;Du 2014), but the important is seepage in the fracture, which has not been discussed. The concentration and particle's influence on water and sand inrush.
In this work, permeability attributes of water-sand mixtures are obtained through testing by replicating the design system of water-sand seepage in fractures. The influence of mass concentration in water and particle size of sand on the seepage parameters are tested using specially designed instruments.

Viscosity test of water and sediment
Viscous parameters of water-sand mixture in stress-strain relationships were tested using a NDJ-8S viscosimeter in Fig. 2: Shear strain rate γ of water-sand is defined as where D is the diameter of outer cylinder, d is the diameter of the rotor. Apparent viscosity µ a of water-sand mixture were obtained from the NDJ-8S viscosimeter and the shear stress was calculated as follows By changing the rotational speed of the NDJ-8S viscosimeter, several values of shear strain rate γ and shear stress were obtained and plotted on a γ − τ scatter diagram. According to the shape of the γ − τ diagram, the water-sand mixture was identified as non-Newton fluid, and the viscous parameter of water-sand was obtained through linear regression.
It can be seen that the shear strain rate γ increases monotonously along with the shear strain of the water-sand mixture. Therefore, we assume that the water-sand mixture is a power law fluid as follows where C is the consistency coefficient, n is the power exponent.
Combining Eqs. 2 and 3 yields expression of apparent viscosity as follows  Through linear regression, viscous parameters of water-sand (consistency coefficient C and power exponent n) were obtained, as shown in Table 2. It was deduced that the water-sand mixture was a pseudo-plastic fluid, whose viscous parameters changed with sand particle d s and mass concentration of sand ρ s . Different consistencies were tested of coefficient C and power exponent n with the diameters of sand particle sizes 0.038-0.044, 0.061-0.080, 0.090-0.109 and 0.120-0.180 mm; and sand 20, 40, 60 and 80 kg/m 3 in the water. The testing results of consistency coefficient C and power exponent n are shown in Table 2.
From Table 2, consistency coefficient C increases with mass concentration in water as exponential relationship, and decreases along with the increase of sand particle; (4) µ a = Cγ n−1 power exponent n increases along with the increase of mass concentration in water, and decreases along with the increase of sand particle.

Test principle
Figure 4 demonstrates a model of seepage in a fracture. In this paper, the aperture of fracture is 0.75 mm, the length is 12.5 mm, the height is 75 mm. b is the aperture of fracture, h is the height of the fracture, and L is the sample length. According to Fig. 4, we can get Eq. 5.
where V is the velocity of seepage, Q is the flow of seepage. For the fracture, Re is defined as Eq. 6 (Javadi et al. 2014).
where Re is Reynolds number, ρ is the density, Q is the flow of seepage, μ is the fluid viscosity.

Fig. 4 Seepage in parallel fracture
where µ is fluid viscosity, β is non-Darcy factor, the pressure is p, ∂p ∂l is the pressure gradient, c a is the acceleration of water and sand, b 1 is two term coefficient.
Because of water and sand permeability parameter's particularity (permeability parameter is relevant to liquid and fracture), we use µ e , k e to describe the water and sand of effective viscosity µ e , effective permeability k e , as shown in Eq. 8 (Liu 2014).
As for one kind of non-Newton fluid, liquid viscosity and permeability in fracture of water-sand mixture were related to fluid properties and fracture aperture. Therefore, liquid viscosity and permeability were not obtained separately, and the effective fluidity I e was introduced to simplify the expression.
The Eq. 8 can be changed into Equation 10 calculated the momentum conservation of water-sand seepage in the fracture. For the seepage in Fig. 4, the steady-flow method was selected to measure water-sand seepage in the fracture. Equation 10 can be deduced into Eq. 11, Substituting Eq. 5 into Eq. 11 yields Eq. 12 b is the aperture of the fracture, m is the mass of sand and water.
For the length, the integrating range is [0, L]; the mass is m, the pressure of water and sand at the entrance wall were: The definite integral of Eq. 12 on the interval [0, L] was Introducing the sign 1 = 1 I e 1 bh n , 2 = mβ (hb) 2 , Therefore, Eq. 14 was obtained by using In the test, 5 flows were set as Q i , i = 1, 2, . . . , 5. Steady state values of inlet pressures were tested, and coefficients 1 and 2 were fitted. The specific process was as follows: Equation 15 was obtained In order to get the least value of the flow Q, Eq. 16 can be set as Eq. 17.
1 and 2 were solved by Eq. 16, effective mobility I e and non-Darcy β were obtained.

Experimental equipment and steps
Based on testing principles, a set of experimental system was designed and manufactured as shown in Fig. 5. Sand comes from the surface of the mine in northwest of China. The rock sample is the sandstone under −265 m from the Luan mine in Shanxi, China.
There are five specimens of rock fracture with Joint Roughness Coefficient (JRC) 4-6, the velocity of seepage was obtained. Figure 6 illustrates the entire experimental procedure. The test steps were as follows: 1. The test system was assembled according to Fig. 6 and the sample was loaded. The leakage of the experiment system was tested. 2. The sand grain with a diameter of 0.038-0.044 mm was placed into the mixing pool and the sand concentration was 20 kg/m 3 in water.

Pressure graduate
According to the pressure and velocity measured in the tests, the pressure gradient and velocity under different sand concentration in water were shown in Fig. 7 Table 3 lists the test result of pressure gradient and seepage velocity of the water and sand, the polynomial fitting formula and its coefficient, the power fitting formula and its coefficient. Figure 7 and Table 3 been presented above, we can obtain: the seepage velocity of water and sand increases with pressure gradient increasing, Moreover, the greater the sand concentration in water is, the lower the seepage velocity is.

Permeability of water and sand in the fracture
Keeping the fracture aperture 0.75 mm, the permeability parameters of water and sand seepage in the fracture under particle sizes 0.041, 0.071, 0.100 and 0.150 mm are tested at 20, 40 and 60 kg/m 3 sand concentration in water, as shown in Fig. 8, liquid viscosity and permeability were not obtained separately, and the effective fluidity I e was introduced to simplify the expression. Because of the permeability parameters of water and sand seepage in fracture are connected with water and sand, at the same time, the structure of fracture; so the permeability k is not enough to describe permeability parameters, the effective fluidity and non-Darcy factor β are used.
The 5 samples were used to obtain the permeability parameters in test, and we adopted the arithmetic mean values, as shown in Table 4.
Fitting the curves of Fig. 8, the functional relationship between seepage parameters and sand concentration in water was used, as shown in Table 5.
The exponential function was used to fit the relationship between effective fluidity, non-Darcy coefficient and particle sizes of sand. The power exponent equations are used to fit the relationship between effective fluidity I e , the non-Darcy factor β and Sand concentration.
From Fig. 8 and Table 5, the following results were obtained: 1. The seepage of water and sand in a fracture is nonlinear. 2. Along with the change of grain size of sediment, the relationship between effective fluidity I e and mass concentration of sand d s was the negative exponential relationship; the absolute value of the exponent increased along with the increase of sand particle in the water. 3. Non-Darcy factor β and sand concentration in water had a positive exponential relationship; the absolute value of the exponent increased along with the decrease of sand particle in water.

Discussion
It is non-Darcy flow in the paper, which was influenced by roughness, flow velocity, aperture of fracture, and so on. Roughness has a large influence on fracture flow, where non-Darcy also happened (Boutt et al. 2006;Lomize 1951;Louis 1969;Qian et al. 2011).
During the flow, Reynolds number and Forchheimer's number are important parameters to judge (Bear 1972): when Re > 100 or Re < 1, it will be nonlinear flow and does not conform to Darcy flow. What's more, the velocity of water and sand, the aperture of fracture and the tortuosity of fracture also have much influence on flow parameters (Tsang 1984;Tsang and Tsang 1987). The concentration and density also have influence on flow character in fracture (Watson et al. 2002;Tenchine and Gouze 2005). Here I e has relationship with the structure of fracture, and the character of mixture or water and sand. With the pressure drop increasing, the nonlinear flow became obvious (Elsworth and Doe 1986;Wen et al. 2006;Yeo and Ge 2001) the Forchheimer's law is well known classical approach to describe the nonlinear flow in fracture. Non-Darcy factor β is the parameter which reflected the deviation of Darcy of the seepage. Along with sand particle in water, the non-Darcy character became more obvious.

Conclusion
In this paper, the viscosity of water and sand mixture was discussed and the seepage of water and sand mixture in rude fracture was analyzed.  1. The seepage velocity of water and sand in a fracture increases along with the pressure of the fracture, but the relationship between them is nonlinear. 2. Consistency coefficient C becomes larger in conjunction with the mass concentration in water, but decreases along with the particle size of sand. The lower exponent n becomes enlarger along with mass concentration in water, but decreases along with particle size of sand. 3. Along with the change of the grain size of sediment, the relationship between effective fluidity I e and mass concentration of sediment ρ s in water is exponential. The absolute value of the exponent increases along with the increase of sand concentration in water. The non-Darcy factor β and sand concentration in water has a positive exponential relationship and the absolute value of the exponent increases along with the decrease of sand concentration in water. 4. For the future work, we will work for the different concentration, for particle and concentration both has influence to the flow character, but we should do some experiments to make sure which one is more influence. And acceleration, low velocity of water and sand how to change into water and sand inrush.