Influence of particle size on non-Darcy seepage of water and sediment in fractured rock
Yu Liu1,2 &
Shuncai Li1
SpringerPlus volume 5, Article number: 2099 (2016) Cite this article
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.
In China, water and sand inrush is very serious safety problem for coal mining in 20 years, there were many accidents which gave more damage to coal mining (Limin et al. 2015, 2016). The coal reserves are located at shallow depths and the thin bedrock and thick sand overburdens the strata layers, inducing connected cracks. Surface water, groundwater and sand can flow into the mine goaf through the fractured rock and lead to inrush of water and collapsing of sand, which can be seen in Fig. 1.
Water and sand inrush in the fracture
From the mechanical perspective, the result of water and sand erupting, permeating fractured rock reflects the instability of the strata layers. Therefore, studying the seepage properties of fractured rocks plays an important role in coal mining engineering. The inrush of water and sand compromises mine safety by causing instability in stress block beams, which creates surface subsidence and water resource run off.
Field tests that are conducted in order to replicate water and sand inrush are difficult; therefore, many scholars suggested conducting experimental simulations of inrushing 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:
NDJ-8S viscometer
Shear strain rate γ of water–sand is defined as
$$\gamma = \frac{{\pi n_{rot}^{{}} }}{30} \times \frac{d}{(D - d)}$$
where D is the diameter of outer cylinder, d is the diameter of the rotor.
Apparent viscosity \(\mu_{a}^{{}}\) of water–sand mixture were obtained from the NDJ-8S viscosimeter and the shear stress was calculated as follows
$$\tau = \mu_{a}^{{}} \gamma$$
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 \(\gamma - \tau\) scatter diagram. According to the shape of the \(\gamma - \tau\) diagram, the water–sand mixture was identified as non-Newton fluid, and the viscous parameter of water–sand was obtained through linear regression.
During the experiment, the diameters of sand particle are 0.038–0.044, 0.061–0.080, 0.090–0.109 and 0.120–0.180 mm. Firstly, Sand particles 0.061–0.080 mm with 20 kg/m3 sand at 20 °C was measured; the shear strain rate \(\gamma\), apparent viscosity \(\mu_{a}^{{}}\) and stress \(\tau\) of the water–sand mixture was gotten at various rotation rates (Table 1; Fig. 3).
Table 1 Angle strain rate, apparent viscosity and shear stresses at different rotating speed
Scatter plot of angle strain rate—shear stress (0.061–0.080 mm)
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
$$\tau = C\gamma_{{}}^{n}$$
where C is the consistency coefficient, n is the power exponent.
Combining Eqs. 2 and 3 yields expression of apparent viscosity as follows
$$\mu_{a}^{{}} = C\gamma^{n - 1}$$
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 \(\rho_{s}\).
Table 2 Relation of surface viscosity and shearing rate
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/m3 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; power exponent n increases along with the increase of mass concentration in water, and decreases along with the increase of sand particle.
Seepage test of water and sand in a fracture
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.
Seepage in parallel fracture
According to Fig. 4, we can get Eq. 5.
$$V = \frac{Q}{bh}$$
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).
$$R_{e} = \frac{\rho Q}{\mu b}$$
where Re is Reynolds number, ρ is the density, Q is the flow of seepage, μ is the fluid viscosity.
In the paper, Q is \(6.00 \times 10^{ - 4} {-}3.10 \times 10^{ - 3}\) m3/s, \(\rho = 1.02{-}1.08 \times 10^{3}\) kg/m3, \(\mu = 1.005\;{\text{mp}}_{\text{a}} \;{\text{s}}\).
So, \(R_e = \frac{\rho Q}{\mu b} = 76.5{-}421.2\) in case of higher Reynolds numbers (\(R_e \gg 1\)), the pressure losses pass from a weak inertial to a strong inertial regime, described by the Forchheimer equation (Forchheimer 1901; Chin et al. 2009; Cherubini et al. 2012, 2013; Javadi et al. 2010; Li et al. 2008), given by:
$$\rho c_{a} \frac{\partial V}{\partial t} = - \frac{\partial p}{\partial l} - \frac{\mu }{k}V - \rho \beta V^{2}$$
where \(\mu\) is fluid viscosity, \(\beta\) is non-Darcy factor, the pressure is \(p\), \(\frac{\partial p}{\partial 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 \(\mu_{e}\), \(k_{e}^{{}}\) to describe the water and sand of effective viscosity \(\mu_{e}\), effective permeability \(k_{e}^{{}}\), as shown in Eq. 8 (Liu 2014).
$$\rho c_{a}^{{}} \frac{\partial V}{\partial t} = - \frac{\partial p}{\partial l} - \frac{{\mu_{e}^{{}} }}{{k_{e}^{{}} }}V_{{}}^{n} - \rho \beta V_{{}}^{2}$$
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.
$$I_{e}^{{}} = \frac{{k_{e} }}{{\mu_{e} }}$$
The Eq. 8 can be changed into
$$\rho c_{a}^{{}} \frac{\partial V}{\partial t} = - \frac{\partial p}{\partial l} - \frac{1}{{I_{e}^{{}} }}V_{{}}^{n} - \beta \rho V_{{}}^{2}$$
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,
$$\frac{1}{{I_{e} }}V^{n} + \beta \rho V^{2} = - \frac{\partial p}{\partial l}$$
Substituting Eq. 5 into Eq. 11 yields Eq. 12
$$- dp = \frac{1}{{I_{e} }}\left( {\frac{Q}{bh}} \right)^{n} dl + \beta \rho \left( {\frac{Q}{bh}} \right)^{2} dl$$
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:
$$\left\{ {\begin{array}{l} {\left. p \right|_{x = 0}^{{}} = p_{0}^{{}} } \\ {\left. p \right|_{x = L}^{{}} = 0} \\ \end{array} } \right.$$
The definite integral of Eq. 12 on the interval [0, L] was
$$p = \frac{L}{{I_{e} }}\left( {\frac{Q}{bh}} \right)^{n} + \beta mL\left( {\frac{Q}{bh}} \right)^{2}$$
Introducing the sign \(\lambda_{1} = \frac{1}{{I_{e} }}\left( {\frac{1}{bh}} \right)^{n}\), \(\lambda_{2} = \frac{m\beta }{{(hb)^{2} }}\),
Therefore, Eq. 14 was obtained by using
$$\lambda_{1} Q^{n} + \lambda_{2} Q^{2} - p_{0} = 0$$
In the test, 5 flows were set as \(Q_{i}^{{}} ,i = 1,2, \ldots ,5\). Steady state values of inlet pressures were tested, and coefficients \(\lambda_{1}^{{}}\) and \(\lambda_{2}^{{}}\) were fitted. The specific process was as follows:
Equation 15 was obtained
$$\varPi = \sum\limits_{i = 1}^{5} {\left( {\lambda_{1}^{{}} Q_{i}^{n} + \lambda_{2}^{{}} Q_{i}^{2} - p_{0}^{i} } \right)_{{}}^{2} } = 0$$
In order to get the least value of the flow Q, Eq. 16 can be set as Eq. 17.
$$\left\{ \begin{aligned} \left( {\sum\limits_{i = 1}^{5} {Q_{i}^{n} Q_{i}^{n} } } \right)\lambda_{1}^{{}} + \left( {\sum\limits_{i = 1}^{5} {Q_{i}^{2} Q_{i}^{n} } } \right)\lambda_{2}^{{}} = \left( {\sum\limits_{i = 1}^{5} {Q_{i}^{n} p_{0}^{i} } } \right) \hfill \\ \left( {\sum\limits_{i = 1}^{5} {Q_{i}^{2} Q_{i}^{n} } } \right)\lambda_{1}^{{}} + \left( {\sum\limits_{i = 1}^{5} {Q_{i}^{2} Q_{i}^{2} } } \right)\lambda_{2}^{{}} = \left( {\sum\limits_{i = 1}^{5} {Q_{i}^{2} p_{0}^{i} } } \right) \hfill \\ \end{aligned} \right.$$
\(\lambda_{1}\) and \(\lambda_{2}\) were solved by Eq. 16, effective mobility \(I_{e}\) and non-Darcy \(\beta\) 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.
Scheme of system principles. 1 Sidebend; 2 pressure transmitter; 3 injection pipe; 4 vane pump; 5 agitator tank; 6 VVVF; 7 screw pump; 8 flow sensor; 9 piezometer; 10 filter box; 11 scheme of permeameter
Figure 6 illustrates the entire experimental procedure. The test steps were as follows:
Flow chart of the test
The test system was assembled according to Fig. 6 and the sample was loaded. The leakage of the experiment system was tested.
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/m3 in water.
To control the motor speed, flow and pressure under different rotational speeds were recorded while the fracture aperture 0.75 mm; the motor speeds, 200, 400, 600, 800, 1000 r/min were changed separately. Different pressures and seepage velocities of the fracture were obtained using a paperless recorder. The sand concentration \(\rho_{s}\) in water was 40, 60, 80 kg/m3 respectively.
The flow and pressure under different grain diameters (0.038–0.044, 0.061–0.080, 0.090–0.109 and 0.120–0.180 mm)were recorded during the different rotational speeds. In order to easily calculate the data, we choose the arithmetic mean of each range of the grain diameter, e.g. 0.041, 0.071, 0.100 and 0.150 mm.
According to Eqs. 15 and 16, \(I_{e}\) and \(\beta\) were calculated.
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 and Table 3. 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.
Relationship between pressure gradient and velocity
Table 3 Relationship between pressure gradient and velocity under different sand concentration
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/m3 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.
Curves of permeability parameters changing with \(d_{s}\). a Curve of \(I_{e} - d_{s}\) under 20 kg/m3, b curve of \(\beta - d_{s}\) under 20 kg/m3, c curve of \(I_{e} - d_{s}\) under 40 kg/m3, d curve of \(\beta - d_{s}\) under 40 kg/m3, e curve of \(I_{e} - d_{s}\) under 60 kg/m3, f curve of \(\beta - d_{s}\) under 60 kg/m3, g curve of \(I_{e} - d_{s}\) under 80 kg/m3, h curve of \(\beta - d_{s}\) under 80 kg/m3
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 \(\beta\) 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.
Table 4 permeability parameters of water and sand under different sand concentration
Fitting the curves of Fig. 8, the functional relationship between seepage parameters and sand concentration in water was used, as shown in Table 5.
Table 5 Fitted equations of permeability parameters changing with \(d_{s}\) at JRC 4–6
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 \(\beta\) and Sand concentration.
From Fig. 8 and Table 5, the following results were obtained:
The seepage of water and sand in a fracture is nonlinear.
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.
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.
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.
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.
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.
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.
Along with the change of the grain size of sediment, the relationship between effective fluidity \(I_{e}\) and mass concentration of sediment \(\rho_{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.
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.
b :
aperture of fracture
consistency coefficient
C a :
acceleration coefficient
d :
diameter of the rotor
diameter of outer cylinder
d s :
particle size of sand
h :
height of the fracture
I e :
effective fluidity
k e :
effective permeability
sample length
m :
mass of sand and water
power exponent
n rot :
rotate speed of rotor
β :
non-Darcy coefficient
Q :
flow of seepage
\(\tau\) :
shear stress
\(\mu\) :
fluid viscosity
\(\mu_{a}\) :
apparent viscosity of water and sand
\(\mu_{e}\) :
effective viscosity
V :
velocity of seepage
\(\gamma\) :
apparent viscosity
\(\rho\) :
\(\rho_{s}\) :
mass concentration of sand
\(\frac{\partial p}{\partial l}\) :
pressure gradient
Bear J (1972) Dynamics of fluids in porous media. New York, Dover
Boutt DF, Grasselli G, Fredrich JT, Cook BK, Williams JR (2006) Trapping zones: the effect of fracture roughnesson the directional anisotropy of fluid flow and colloid transport in a single fracture. Geophys Res Lett 33:L21402. doi:10.1029/2006GL027275
Cherubini C, Giasi CI, Pastore N (2012) Bench scale laboratory tests to analyze non-linear flow in fractured media. Hydrol Earth Syst Sci 16:2511–2522. doi:10.5194/hess-16-2511-2012
Cherubini C, Giasi CI, Pastore N (2013) Evidence of non-Darcy flow and non-Fickian transport in fractured media at laboratory scale[J]. Hydrol Earth Syst Sci 17:2599–2611. doi:10.5194/hess-17-2599-2013
Chin DA, Price RM, DiFrenna VJ (2009) Nonlinear Flow in Karst Formations[J]. Ground Water 47:669–674. doi:10.1111/j.1745-6584.2009.00574.x
Du F (2014) Experimental investigation of two phase water–sand flow characteristics in crushed rock mass. China University of Mining & Technology, Xuzhou
Elsworth D, Doe TW (1986) Application of non-linear flow laws in determining rock fissure geometry from single borehole pumping tests. Int J Rock Mech Min Sci 23(3):245–254. doi:10.1016/0148-9062(86)90970-8
Forchheimer P (1901) Wasserbewegung durch Boden. Z Ver Dtsch Ing 45:1781–1788
Javadi M, Sharifzadeh M, Shahriar K (2010) A new geometrical model for nonlinear fluid flow through rough fractures. J Hydrol 389:18–30
Javadi M, Sharifzadeh M, Shahriar K, Mitani Y (2014) Critical Reynolds number for nonlinear flow through rough-walled fractures: the role of shear process. Water Resour Res 50(2):1789–1804
Li SC, Miao XX, Chen ZQ, Mao XB (2008) Experimental study on Seepage properties of non-Darcy flow in confined broken rocks. Eng Mech 25(4):85–92
Limin F, Xiongde M, Ruijun J (2015) Progress in engineering practice of water-preserved coal mining in western eco-environment frangible area. J China Coal Soc 40(8):1711–1717. doi:10.13225/j.cnki.jccs.2015.0223
Limin F, Xiongde M, Hui J et al (2016) Risk evaluation on water and sand inrush in ecologically fragile coal mine. J China Coal Soc 41(3):531–536. doi:10.13225/j.cnki.jccs.2015
Liu Y (2014) Test study of non-Darcy flow of water–sand mixture in fractured rock. China University of Mining and Technology, Xuzhou
Lomize GM (1951) Filtratsia v treshchinovatykh porodakh, Seepage in Jointed Rocks. Gosudarstvennoe Energeticheskoe Izdatel'stvo, Moskva-Leningrad
Louis C (1969) A study of groundwater flow in jointed rock and its influence on the stability of rock masses. Rock mechanics research report. Imperial College, London, UK
Qian JZ, Zhan H, Chen Z, Ye H (2011) Experimental study of solute transport under non-Darcian flow in a single gracture. J Hydrol 399:246–254
Sui WH, Cai GT, Dong QH (2007) Experimental research on critical percolation gradient of quicksand across overburden fissures due to coal mining near unconsolidated soil layers. Chin J Rock Mech Eng 26(10):084–2091
Sui WH, Dong QH, Cai GT, Yang WF (2008) Mechanism and prevention of sand inrush during mining. Geological House
Tenchine S, Gouze P (2005) Density contrast effects on tracer dispersion in variable aperture fractures. Adv Water Resour 28:273–289
Tsang YW (1984) The effect of tortuosity on fluid-flow through a single fracture. Water Resour Res 20(9):1209–1215. doi:10.1029/WR020i009p01209
Tsang YW, Tsang CF (1987) Channel model of flow through fractured media. Water Resour Res 23(3):467–479. doi:10.1029/WR023i003p00467
Watson SJ, Barry DA, Schotting RJ, Hassanizadeh SM (2002) On the validity of Darcy's law for stable high-concentration displacements in granular porous media. Transp Porous Med 47:149–167
Wen Z, Huang G, Zhan H (2006) Non-Darcian flow in a single confined vertical fracture toward a well. J Hydrol 330(3–4):698–708. doi:10.1016/j.jhydrol.2006.05.001
Wu YP (2004) Condition analysis on sand inrush in shallow stope. Mine Press Support 4:57–59
Xu YC, Wang BS, You SW (2012) Mechanism and criteria of crushing sand near loosening sand stone aquifer. J Xi'an Univ Sci Technol 32(1):63–69
Yang WF (2009) Overburden failure in thin bedrock and characteristics of mixed water and sand flow induced by mining. China University of Mining Technology, Xuzhou
Yang WF, Sui WH, Ji YB, Zhao GR (2012) Experimental research on the movement process of mixed water and sand flow across overburden fissures in thin bedrock by mining. J China Coal Soc 37(1):141–146
Yeo IW, Ge S (2001) Solute dispersion in rock fractures by non-Darcian flow. Geophys Res Lett 28(20):3983–3986. doi:10.1029/2001GL013274
Zhang YJ, Kang YH, Liu X (2006) Predicting on inrush of sand of mining under loosening sandstone aquifer. J China Coal Soc 31(4):429–432
Zhang, GM, Zhang K, Wang LJ, Wu Y (2015a) Mechanism of water inrush and quicksand movement induced by a borehole and measures for prevention and remediation. Bull Eng Geol Environ 74:1395–1405
Zhang JX, Zhang Q, Sun Q, Gao R, Deon G, Sami A (2015b) Surface subsidence control theory and application to backfill coal mining technology. Environ Earth Sci 74(2):1439–1448
The work presented here was carried out in collaboration between all authors. Yu Liu defined the research theme, designed experiments methods and wrote the paper. Shuncai Li did the experiments, analyzed the data and explained the results. Both authors read and approved the final manuscript.
This research was supported by Natural science fund for colleges and universities in Jiangsu Province (14KJB440001), Jiangsu Normal University PhD Start Fund (14XLR032), Jiangsu Planned Projects for Postdoctoral Research Funds (1402055B), and National Natural Science Foundation of China (51574228), All the supports are gratefully acknowledged.
Both authors declare that they have no competing interests.
School of Mechanical and Electrical Engineering, Jiangsu Normal University, Xuzhou, 221116, Jiangsu, China
Yu Liu
& Shuncai Li
School of Mines, China University of Mining and Technology, Xuzhou, China
Search for Yu Liu in:
Search for Shuncai Li in:
Correspondence to Yu Liu.
Liu, Y., Li, S. Influence of particle size on non-Darcy seepage of water and sediment in fractured rock. SpringerPlus 5, 2099 (2016) doi:10.1186/s40064-016-3778-9
Accepted: 01 December 2016
Water–sand
Non-Darcy seepage
Non-Newton fluid