# A Computing Method for Sand Inrush Quantity through a Borehole in Longde Coal Mine

Authors: Boyang Zhang; Zhibin Lin
Journal: Advances in Civil Engineering
 (2018)

---

## Abstract

A quicksand disaster through a borehole occurred in Longde coal mine. A lot of aeolian sand, the volume of which is between 310,000 m3 and 380,000 m3, has submerged into the underground space in about 70.5 h. The volume flux of quicksand cannot be calculated accurately by the empirical method. Based on the method of fluid mechanics, an all-purpose computing method for quicksand disaster through a borehole was proposed. The result shows that the inrush volume of sand into underground space was between 310,000 m3 and 350,000 m3, which was consistent with the actual result. To apply and popularize this method, the impact laws of water yield properties of an aquifer on the volume flux were discussed. The all-purpose computing method can be suitably used for the volume flux calculation of quicksand disaster through the borehole.


---

## Body

## 1. Introduction

Water inrush is one of the most severe disasters during mining in China [1–4]. There are vast deserts [5] in the ground and abundant coal resources underground [6] in Shanxi, Shaanxi, Inner Mongolia, Sinkiang, and Qinghai provinces in China. The hidden dangers of quicksand exist extensively in Western China during mining [7]. Quicksand disaster happened occasionally with the water inrush, which have caused serious economic losses [8–11].There are two types of quicksand disaster. Firstly, the quicksand disaster caused by the aeolian sand submerging into underground space through caved zone as shown in Figure1(a). When the mining depth is shallow, the mining-induced fractures may go through the rock layers to the sand bed, or even to the ground surface. In this situation, the water and sand can move to the goaf through the fractures and granular rock mass, which results in the quicksand disaster [14]. Secondly, the quicksand disaster caused by a borehole from ground surface to underground. For example, a disaster of water inrush and quicksand caused by a borehole occurred in Longde coal mine. The borehole directly transfixed the underground space under the condition of no casing pipe, which resulted in the disaster of water inrush and quicksand, as shown in Figure 1(b). A lot of aeolian sand submerged into the underground space. The sand cleaning work lasted for almost one year, which caused serious economic losses and mining difficulties. The quicksand disasters caused by the borehole have been reported for many times in papers [15, 16], but this phenomenon has not been explained reasonably until now.Figure 1
Quicksand disaster schematic diagram. (a) Quicksand trip through overburden strata above a longwall panel [12, 13]. (b) Quicksand trip through drill.
(a)
(b)Flow in the borehole (i.e., in the circular pipe) has been researched by researchers, and a lot of empirical formulas have been put forward.Manning improved the Manning coefficient based on Chery’s research [17]. The Chery formula is(1)v=1nR1/6RJ,where v is the velocity, R is the hydraulic radius, J is the hydraulic slope, and n is the Manning coefficient.Hazen and Williams [18] proposed the Hazen–Williams formula for transition region:(2)v=0.355Chd0.63J0.54,where d is the diameter of circular pipe and Ch is the empirical coefficient of roughness height.Based on the research of wood-stave pipe, Scobey [19] put forward the Scobey formula:(3)hf=KsLQ1.9d4.9,where hf is the frictional head loss, L is the length of circular pipe, Q is the volume flux, and Ks is the empirical coefficient.In addition, there are many other empirical formulas such as Blasius formula [20] and Churchill formula [21]. Blasius formula is suitable for the Reynolds number below 105 based on a number of experiment results. Churchill formula can be used in both laminar flow and turbulence flow.The volume flux of Longde quicksand disaster was calculated by using the empirical formula with the common value of empirical parameters, and the result is shown as Table1. The actual volume flux of quicksand is 4,400–5,400 m3/h (Section 2.3). As shown in Table 1, the calculation results of Chezy formula and Scobey formula are far from the actual volume flux, and that of the Hazen–William formula contains the actual flux. The simple form of these empirical formulas contributes to the convenient usage in the engineering. However, there are two shortages: on the one hand, the empirical formulas are hard to select under different conditions of disasters. The Hazen–William formula can be used to calculate the volume flux in Longde quicksand disaster. However, if the geological conditions of pressure, water head, or diameter of borehole differ from that of Longde disaster, the calculation results of the Hazen–William formula can also appear with errors. On the other hand, each formula has an empirical parameter which is hard to determine in the calculation process.Table 1
Empirical formula and calculation result [17–19].
Formula
Equation
Common value of empirical parameters
v (m/s)
Q (m3/h)
Chezy formula
v
=
1
/
n
R
1
/
6
RJ
0.010–0.030
5.44–16.32
1537–4611
Hazen–William formula
v
=
0.355
C
h
d
0.63
J
0.54
100–140
15.40–21.56
4352–6093
Scobey formula
h
f
=
K
s
L
Q
1.9
/
d
4.9
0.3–0.6
0.48–0.70
137–1977
Actual result
—
—
—
4400–5400Therefore, an all-purpose computing method was proposed for quicksand through the borehole in this paper. The Bernoulli formula [22], Darcy–Weisbach formula [17], semiempirical Nikuradse formula [23, 24], and Colebrook–White formula [25, 26] were combined to calculate the volume flux of quicksand. Besides, the volume concentration, density, and viscosity of sand-water mixture were tested for the method used in similar engineering. The all-purpose computing method for the calculation of volume flux can be used in all quicksand disasters though the borehole, which has a certain guiding significance in engineering.


## 2. Illustration of Water Inrush and Quicksand in Longde Coal Mine

### 2.1. The Process of the Disaster

Longde coal mine is located in southwest of Yulin in Shaanxi province in China, and the shape of this field is similar to a pistol, as shown in Figure2. This area is vulnerable to long-term effects of polar continent. As a result, the continental climate is remarkable, and the ground is dominated by desert and hills in this area.Figure 2
Landscapes of Longde mine field.A water inrush and quicksand disaster through a borehole happened in Longde coal mine at 14:30, September 17, 2012. A geological team undertook the construction of cable borehole for the central water pump room. The wrong statistics of drilling depth and the cable borehole directly getting into the underground space without the protection measures such as steel casing and seal ring caused immediate disappearance of the drill tower in sand seam and submergence of the water pump room by sand. Then, the central substation could not deliver electricity normally. As a result, the whole coal mine was submerged by sand.From September 17 to September 20, although emergency measures had been taken by Longde coal mine, a lot of sand flowed into the railway, even main shaft, and auxiliary shaft, as shown in Figure3. The sandline elevation of quicksand position, main shaft, and auxiliary shaft was 1,227.1 m, 1,081 m, and 1,100 m, respectively. However, the waterline was only 1,029 m less than the sandline elevation of main and auxiliary shaft. It was mainly because that the aeolian sand has become a fluid with the lubrication action of water.Figure 3
The profile after water inrush and quicksand disaster.


### 2.2. Sand Volume Flux Analysis

The total volume of the underground space was 488,000 m3 in Longde coal mine. The statistical data of sand cleaning presented that the actual volume of sand was 380,000 m3. The 77.9% volume of the underground space was submerged by sand. After the disaster, in order to protect the infrastructure such as the air shaft square, protective measures were taken; for instance, the gravity dam was set up with stones and concrete. And based on geological conditions, the sand subsidence area did not present a cone type but an irregular cone type (prismoid), as shown in Figure 4. The inflection points of sand subsidence area were monitored after the disaster happened 70.5 h. Figure 4 is the aerial view of the aeolian sand subsidence area. The inflection points were recorded by Xi’an 80 coordinate [27]. For the sake of illustration, the origin point was transformed to C1 point, as shown in Figure 4.Figure 4
Aerial view in subsidence area of aeolian sand.


#### 2.2.1. Disaster Duration

The coordinate of 7 inflection points is shown in Figure4. It is assumed that two adjacent points are in a straight line, and the area of aerial view is calculated by the following equation:(4)A=12∑i=17xi⋅yi+1−xi+1⋅yi,where A is the area of aerial view, x is the longitude (east) coordinate, and y is the latitude (north) coordinate.The calculation result shows that the area of aerial view is 22,019.75 m2 until 70.5 h after the disaster, but it only increases 40 m2 from 70.5 h to 238.5 h. The area percent change only increases 0.18%, which is almost negligible. Therefore, it can be confirmed that the disaster mainly happens within 70.5 h.


#### 2.2.2. The Volume of Sand Submerged in Underground Space

It is necessary to calculate the volume of sand subsidence area for the obtainment of quicksand volume. Consequently, the irregular cone type or prismoid type can be assumed for presenting the sand subsidence area based on the method in [15]. Considering that the friction coefficient of the sand is constant, the friction coefficient is tested, as shown in Figure 5. The friction force is generated by the gravity and contact force among particles. The equation of the friction coefficient is k = tanθ, where θ can be obtained as shown in Figure 5, and the average value of friction coefficient k is 0.617.Figure 5
Test of friction coefficient.To calculate the volume of sand subsidence area, the profile plane equation of each plane is established firstly. The uniform form of plane equation is as follows [28]:(5)ax+by+cx+e=0,where a, b, c, and e are calculating parameters of plane equation.The plane equation parameters of each profile plane are given in Table2. The volume is calculated by the method [29]. The total volume of sand subsidence area is 310,000 m3, which is less than the statistical data of 380,000 m3 from the sand cleaning. It can be explained that the sand volume varies easily under different confining pressure. The confining pressure is nearly zero when the cleaning sand work starts. But the confining pressure of the ground sand increases with depth. The compaction part of sand causes the volume statistical data of subsidence area bigger than the sand cleaning. In conclusion, the aeolian sand is submerged into the underground space of Longde coal mine in 70.5 h. The volume of sand is between 310,000 m3 and 380,000 m3, and the volume flux is between 4,400 m3/h and 5,400 m3/h.Table 2
Parameters of plane equation.
Parameter
C1-C2
C2-C3
C3-C4
C4-C5
C5-C6
C6-C7
C7-C1
a
0.8750
0.0334
−0.0057
−0.0019
−0.0164
−0.0071
−0.3080
b
1.0000
0.0071
0.0099
0.0081
−0.0167
−0.0017
1.0000
c
2.1501
0.0552
−0.0185
0.0135
−0.0379
0.0118
−2.7304
e
0.0000
1.0000
1.0000
1.0000
1.0000
1.0000
0.0000


### 2.3. Water Volume Flux Analysis

The water level of auxiliary shaft was monitored after the disaster happened 12 h. Figure6 is the water volume flux change curve with time in auxiliary shaft. From the figure, it can be known that the water volume flux was only 12.291 m3/h after the disaster happened 12 h. It decreased to 0.500 m3/h gradually when the disaster happened 30 h. Assuming that water flowed into main shaft, auxiliary shaft, and air shaft at the same rate, then the sectional area of main shaft, auxiliary shaft, and air shaft were 17.813 m2, 22.873 m2, and 28.260 m2, respectively, and the total volume flux was 3.014 times of auxiliary shaft. After the disaster happened 12 h, the total volume flux was 37.05 m3/h, while it was only 1.51 m3/h after the disaster happened 30 h. There would also be a small amount of water flow near to 1.51 m3/h in the mine without disaster. Therefore, the water volume flux variation caused by the disaster could be negligible.Figure 6
Volume flow rate change curve of auxiliary shaft with time in Longde coal mine.In conclusion, the volume flux of sand was between 4,400 and 5,400 m3/h, while water flow was only 37.05 m3/h. The volume flux of sand was much bigger than that of water. Consequently, it can be confirmed that the disaster of water inrush and quicksand is nearly the flow of solid particle in Longde coal mine. Jaeger et al. [30] and MiDi [31] have also presented that tiny solid particle can show the characteristics of fluid.


## 2.1. The Process of the Disaster

Longde coal mine is located in southwest of Yulin in Shaanxi province in China, and the shape of this field is similar to a pistol, as shown in Figure2. This area is vulnerable to long-term effects of polar continent. As a result, the continental climate is remarkable, and the ground is dominated by desert and hills in this area.Figure 2
Landscapes of Longde mine field.A water inrush and quicksand disaster through a borehole happened in Longde coal mine at 14:30, September 17, 2012. A geological team undertook the construction of cable borehole for the central water pump room. The wrong statistics of drilling depth and the cable borehole directly getting into the underground space without the protection measures such as steel casing and seal ring caused immediate disappearance of the drill tower in sand seam and submergence of the water pump room by sand. Then, the central substation could not deliver electricity normally. As a result, the whole coal mine was submerged by sand.From September 17 to September 20, although emergency measures had been taken by Longde coal mine, a lot of sand flowed into the railway, even main shaft, and auxiliary shaft, as shown in Figure3. The sandline elevation of quicksand position, main shaft, and auxiliary shaft was 1,227.1 m, 1,081 m, and 1,100 m, respectively. However, the waterline was only 1,029 m less than the sandline elevation of main and auxiliary shaft. It was mainly because that the aeolian sand has become a fluid with the lubrication action of water.Figure 3
The profile after water inrush and quicksand disaster.


## 2.2. Sand Volume Flux Analysis

The total volume of the underground space was 488,000 m3 in Longde coal mine. The statistical data of sand cleaning presented that the actual volume of sand was 380,000 m3. The 77.9% volume of the underground space was submerged by sand. After the disaster, in order to protect the infrastructure such as the air shaft square, protective measures were taken; for instance, the gravity dam was set up with stones and concrete. And based on geological conditions, the sand subsidence area did not present a cone type but an irregular cone type (prismoid), as shown in Figure 4. The inflection points of sand subsidence area were monitored after the disaster happened 70.5 h. Figure 4 is the aerial view of the aeolian sand subsidence area. The inflection points were recorded by Xi’an 80 coordinate [27]. For the sake of illustration, the origin point was transformed to C1 point, as shown in Figure 4.Figure 4
Aerial view in subsidence area of aeolian sand.


### 2.2.1. Disaster Duration

The coordinate of 7 inflection points is shown in Figure4. It is assumed that two adjacent points are in a straight line, and the area of aerial view is calculated by the following equation:(4)A=12∑i=17xi⋅yi+1−xi+1⋅yi,where A is the area of aerial view, x is the longitude (east) coordinate, and y is the latitude (north) coordinate.The calculation result shows that the area of aerial view is 22,019.75 m2 until 70.5 h after the disaster, but it only increases 40 m2 from 70.5 h to 238.5 h. The area percent change only increases 0.18%, which is almost negligible. Therefore, it can be confirmed that the disaster mainly happens within 70.5 h.


### 2.2.2. The Volume of Sand Submerged in Underground Space

It is necessary to calculate the volume of sand subsidence area for the obtainment of quicksand volume. Consequently, the irregular cone type or prismoid type can be assumed for presenting the sand subsidence area based on the method in [15]. Considering that the friction coefficient of the sand is constant, the friction coefficient is tested, as shown in Figure 5. The friction force is generated by the gravity and contact force among particles. The equation of the friction coefficient is k = tanθ, where θ can be obtained as shown in Figure 5, and the average value of friction coefficient k is 0.617.Figure 5
Test of friction coefficient.To calculate the volume of sand subsidence area, the profile plane equation of each plane is established firstly. The uniform form of plane equation is as follows [28]:(5)ax+by+cx+e=0,where a, b, c, and e are calculating parameters of plane equation.The plane equation parameters of each profile plane are given in Table2. The volume is calculated by the method [29]. The total volume of sand subsidence area is 310,000 m3, which is less than the statistical data of 380,000 m3 from the sand cleaning. It can be explained that the sand volume varies easily under different confining pressure. The confining pressure is nearly zero when the cleaning sand work starts. But the confining pressure of the ground sand increases with depth. The compaction part of sand causes the volume statistical data of subsidence area bigger than the sand cleaning. In conclusion, the aeolian sand is submerged into the underground space of Longde coal mine in 70.5 h. The volume of sand is between 310,000 m3 and 380,000 m3, and the volume flux is between 4,400 m3/h and 5,400 m3/h.Table 2
Parameters of plane equation.
Parameter
C1-C2
C2-C3
C3-C4
C4-C5
C5-C6
C6-C7
C7-C1
a
0.8750
0.0334
−0.0057
−0.0019
−0.0164
−0.0071
−0.3080
b
1.0000
0.0071
0.0099
0.0081
−0.0167
−0.0017
1.0000
c
2.1501
0.0552
−0.0185
0.0135
−0.0379
0.0118
−2.7304
e
0.0000
1.0000
1.0000
1.0000
1.0000
1.0000
0.0000


## 2.2.1. Disaster Duration

The coordinate of 7 inflection points is shown in Figure4. It is assumed that two adjacent points are in a straight line, and the area of aerial view is calculated by the following equation:(4)A=12∑i=17xi⋅yi+1−xi+1⋅yi,where A is the area of aerial view, x is the longitude (east) coordinate, and y is the latitude (north) coordinate.The calculation result shows that the area of aerial view is 22,019.75 m2 until 70.5 h after the disaster, but it only increases 40 m2 from 70.5 h to 238.5 h. The area percent change only increases 0.18%, which is almost negligible. Therefore, it can be confirmed that the disaster mainly happens within 70.5 h.


## 2.2.2. The Volume of Sand Submerged in Underground Space

It is necessary to calculate the volume of sand subsidence area for the obtainment of quicksand volume. Consequently, the irregular cone type or prismoid type can be assumed for presenting the sand subsidence area based on the method in [15]. Considering that the friction coefficient of the sand is constant, the friction coefficient is tested, as shown in Figure 5. The friction force is generated by the gravity and contact force among particles. The equation of the friction coefficient is k = tanθ, where θ can be obtained as shown in Figure 5, and the average value of friction coefficient k is 0.617.Figure 5
Test of friction coefficient.To calculate the volume of sand subsidence area, the profile plane equation of each plane is established firstly. The uniform form of plane equation is as follows [28]:(5)ax+by+cx+e=0,where a, b, c, and e are calculating parameters of plane equation.The plane equation parameters of each profile plane are given in Table2. The volume is calculated by the method [29]. The total volume of sand subsidence area is 310,000 m3, which is less than the statistical data of 380,000 m3 from the sand cleaning. It can be explained that the sand volume varies easily under different confining pressure. The confining pressure is nearly zero when the cleaning sand work starts. But the confining pressure of the ground sand increases with depth. The compaction part of sand causes the volume statistical data of subsidence area bigger than the sand cleaning. In conclusion, the aeolian sand is submerged into the underground space of Longde coal mine in 70.5 h. The volume of sand is between 310,000 m3 and 380,000 m3, and the volume flux is between 4,400 m3/h and 5,400 m3/h.Table 2
Parameters of plane equation.
Parameter
C1-C2
C2-C3
C3-C4
C4-C5
C5-C6
C6-C7
C7-C1
a
0.8750
0.0334
−0.0057
−0.0019
−0.0164
−0.0071
−0.3080
b
1.0000
0.0071
0.0099
0.0081
−0.0167
−0.0017
1.0000
c
2.1501
0.0552
−0.0185
0.0135
−0.0379
0.0118
−2.7304
e
0.0000
1.0000
1.0000
1.0000
1.0000
1.0000
0.0000


## 2.3. Water Volume Flux Analysis

The water level of auxiliary shaft was monitored after the disaster happened 12 h. Figure6 is the water volume flux change curve with time in auxiliary shaft. From the figure, it can be known that the water volume flux was only 12.291 m3/h after the disaster happened 12 h. It decreased to 0.500 m3/h gradually when the disaster happened 30 h. Assuming that water flowed into main shaft, auxiliary shaft, and air shaft at the same rate, then the sectional area of main shaft, auxiliary shaft, and air shaft were 17.813 m2, 22.873 m2, and 28.260 m2, respectively, and the total volume flux was 3.014 times of auxiliary shaft. After the disaster happened 12 h, the total volume flux was 37.05 m3/h, while it was only 1.51 m3/h after the disaster happened 30 h. There would also be a small amount of water flow near to 1.51 m3/h in the mine without disaster. Therefore, the water volume flux variation caused by the disaster could be negligible.Figure 6
Volume flow rate change curve of auxiliary shaft with time in Longde coal mine.In conclusion, the volume flux of sand was between 4,400 and 5,400 m3/h, while water flow was only 37.05 m3/h. The volume flux of sand was much bigger than that of water. Consequently, it can be confirmed that the disaster of water inrush and quicksand is nearly the flow of solid particle in Longde coal mine. Jaeger et al. [30] and MiDi [31] have also presented that tiny solid particle can show the characteristics of fluid.


## 3. Semiempirical Method for Quicksand in Longde Coal Mine

### 3.1. The Generalization of Quicksand Disaster

The photograph of the borehole and chamber after cleaning is shown in Figure7. It can be clearly seen that the borehole was not damaged in the whole process of quicksand disaster, and the drill pipe remained in the borehole all the time. The profile of the quicksand process is shown in Figure 8:(1)
There is aeolian sand layer with a thickness of 30 m in the ground. The lithology is silty fine sand with a distinct thickness variation. The sand seam constitutes a unified aquifer together with the Sara Wusu group aquifer below.(2)
The Sara Wusu group aquifer is the main aquifer of the coal roof. The main lithology is based on silty fine sand, medium-coarse sandstone, and loam. The units-inflow is between 0.0441 and 0.0569 L·s−1·m−1. The aquifer belongs to low water abundance aquifer, as shown in Table 3.(3)
The borehole diameter is 325 mm in quicksand disaster, including a built-in drill pipe with the diameter of 75 mm. And the borehole and drill pipe have not been damaged in the disaster. Therefore, the cross-sectional areaS, wetted perimeter χ, and hydraulic radius R are 0.0785 m2, 1.256 m, and 0.0625 m, respectively.Figure 7
Borehole and chamber photo after quicksand disaster.Figure 8
Cutaway view of sediment-water flow model.Table 3
Division of water abundance level [32].
Water abundance level
Higher
High
Middle
Low
Units-inflow L·s−1·m−1
>5
1∼5
0.5∼1
<0.5


### 3.2. The Building of Computing Method

Before finding the methods used for quicksand in Longde coal mine, the fluid condition of sand-water mixture flow has to be judged firstly, as different formulas need to be applied for the calculation of the sand-water mixture flow in different fluid conditions. For instance, the flow is easy to obtain with an accurate analytic solution in laminar flow, while the turbulence flow is complicated, which has to be calculated by semiempirical method. The fluid condition can be judged by Reynolds number:(6)Re=ρvdμ,where Re is the Reynolds number, ρ is the density, v is the velocity, d is the diameter, equivalent diameter d = 4R (in this paper), and μ is the dynamic viscosity.In general,Re of laminar flow is below 2,300 and Re of turbulence flow is above 2,300. Nikuradse [23] divided the turbulence flow into effectively smooth flow, transition region, and rough flow.When the flow is in effectively smooth flow, the range of Re is(7)4000<Re<26.98dks8/7,where ks is the roughness height.When the flow is in the transition region, the range of Re is(8)6.98d2ks8/7<Re<4160d2ks8/7.When the flow is in rough flow, the range of Re is(9)Re>4160d2ks8/7.The velocity of sand-water mixture flow is calculated by viscous Bernoulli formula [22]. The equation is(10)z1+p1ρg+α1v122g=z2+p2ρg+α2v222g+hw,where hw is the head loss, z is the third Cartesian coordinate of cross-sectional area, as shown in Figure 9, p is the relative pressure (that the absolute pressure subtracts the atmosphere pressure) of cross-sectional area, and α1 and α2 are correction factors, which are related to the velocity of cross-sectional area. The 1-1 area and 2-2 area (shown as Figure 9) are the same, so α1 is approximately equivalent to α2.Figure 9
Computing model for sediment-water flow.From the continuity equation of fluid mechanics [22], it can be known that there is a relationship between the volume flux Q of 1-1 and 2-2 cross-sectional area as follows:(11)Q1=Q2.For the borehole drilled with steel casing all the while, the sectional areaS of borehole remains unchanged, and v1 is same to v2. The sand-water mixture flow contacts with air in 2-2 area. The relative pressure of 2-2 area is near to zero. The 2-2 area is taken as the base level, and the third Cartesian coordinate z2 is zero. Equation (10) is simplified to(12)hw=z1+p1ρsg.The head losshw is constituted by frictional head loss hf and local head loss hj. From Figure 9, it can be seen that the borehole is a long straight pipe with a constant sectional area. There are only low local head loss hj in both inlet and outlet. Due to the value of hj is much less than hf, hj is negligible. So, hw is approximately equal to hf. The hf can be expressed as [22](13)hf=λldv22g,where l is the distance between 1-1 area and 2-2 area and λ is the flow resistance ratio and it can be expressed with different equations in different flow conditions.In laminar flow, the equation is [33](14)λ=64Re.In effectively smooth flow, the equation is [24](15)1λ=2lgReλ2.51.In the transition region, the equation is [25](16)1λ=−2lgReλ2.51+ks3.7d.In rough flow, the equation is [23](17)1λ=2lg3.7dks.From Equations (15)–(17), it can be known that λ is only related to Re in effectively smooth flow. λ is only related to ks in rough flow. λ is related to both λ and ks in the transition region.Put Equation (13) into Equation (12), and the velocity v can be expressed as follows:(18)v=2gdz1+p1/ρsgλl.In conclusion, the computational model diagram is drawn as Figure10. The velocity will be calculated under different flow conditions after all flow parameters have been determined. Then, Re will be obtained by the calculating result of velocity. Finally, the flow condition will be verified with Re one by one until the right velocity is worked out.Figure 10
Computational model diagram.


### 3.3. The Calculation of Computing Method

The calculating parameters are shown in Table4. The density and viscosity of aeolian sand are obtained by experiments.Table 4
Parameters of semiempirical calculation.
Parameters
Density (kg/m3)
Dynamic viscosity (Pa·s)
Roughness height (mm)
Water pressure (MPa)
Value
1652.89
9.25
0.1–0.4
0.187The obtaining way of the roughness height in Table4 is as follows: researchers have studied on different roughness heights of different rocks, as shown in Table 5. The roughness heights of all rocks are between 0.9 and 16.2 mm. The rock in Longde coal mine is mainly claystone and sandstone. The heights of them are within the range from 1.0 to 5.3 mm without quicksand. But the changing of roughness height is difficult to obtain in the process of quicksand. It can be known that the volume flux of quicksand is about 4,400 to 5,400 m3/h, and the velocity is about 15 to 19 m/s. The roughness height of borehole will decrease gradually with the high-velocity friction. The phenomenon is just like the grinding wheel sanding samples, and the common speed of grinding wheel is 35 m/s (diameter is 200 mm and rpm is 3,150 r/min) [38]. Therefore, the phenomenon in quicksand disaster can be described as follows: the velocity is slow and the roughness height is high at the beginning. Then, the roughness height decreases gradually under the continuous impact of aeolian sand while the velocity increases with the decrease of roughness height. The interaction between velocity and roughness height repeats until they reach a stable value finally. As the velocity is near to the speed of the grinding wheel, the value of roughness height is near to the minimum roughness height of technical standard, as shown in Table 6. In this way, the value of 0.1 mm varies to 0.4 mm, which is the minimum roughness height value taken as the calculation basis of the semiempirical formula in this paper.Table 5
Roughness height and angle of borehole [34–37].
Rock type
Roughness height (mm)
Roughness angle (°)
Seol et al.
Gneiss
1–7
1.1–8.0
Seidel and Collingwood
Claystone, sandstone
1.7–16.2
1.9–18.9
Shale, limestone, etc.
0.9–6.6
1.0–7.6
Lee et al.
Granite, Gneiss
1–4
1.1–4.6
Sandstone, andesite
1–3.5
1.1–4.0
Nam
Claystone
3.6–5.3
4.1–6.1
Limestone
3.2–3.7
3.7–4.2Table 6
Max roughness of technical standard in part country.
Country
Max roughness of technical standard (mm)
Technical standard
China
0.4
GB1031-83
USA
0.1
BS1134-61/ASAB46.1-62
Switzerland
0.2
VSM10321-62
Italy
0.1
UNI13963-60
Japan
0.4
JISB0601-70
Poland
0. 32
PN58/M042-51The obtaining of the water pressure is as follows: there are 4 hydrological boreholes near by the position of quicksand disaster. The water level elevation is recorded aperiodically, as shown in Figure11. From the figure, it can be known that the water level elevation decreases in the underground engineering but has not decreased violently with the quicksand disaster. The decreasing laws of aquifer are similar before and after the quicksand disaster. Therefore, the average value of the observation data of 4 hydrological boreholes on October 29 is selected as the water level elevation value of the computing method. And the water level elevation of aquifer h1 is 1,198.20 m, and the elevation of aquifer roof h2 is 1,179.15 m. The water pressure can be obtained according to the Pascal law:(19)p1=ρwgh1−h2.Figure 11
Water level elevation of Sarah Wusu group aquifer in hydrological borehole.From Equation (19), the water pressure of aquifer roof p1 is 0.187 MPa.Takingz1 equal to 157 m, p1 equal to 0.187 MPa, ρs equal to 1,652.89 kg/m3 in Equation (12), it can be obtained that hw and hf are equal to 168.54 m and the equivalent diameter d = 4R is equal to 0.25 m. Taking the above parameters into Equation (18),(20)v=5.26λ.The velocity is calculated by the simultaneous Equations (6), (14), and (20), Equations (6), (15), and (20), Equations (6), (16), and (20), and Equations (17) and (20) when the sand-water mixture flow is, respectively, in laminar flow, effectively smooth flow, transition region, and rough flow. The calculation results are shown in Table 7, from which it can be known that the sand-water mixture flow belongs to transition region, and the velocity and volume are 15.32 to 17.78 m/s and 310,000 to 350,000 m3, respectively. The volume result is contained in the actual result 310,000 to 380,000 m3.Table 7
Computing result.
Flow condition
Laminar flow
Effectively smooth flow
Transition region
Rough flow
Velocity (m/s)
5.87 × 104
21.14
15.32–17.78
15.43–18.19
Re
2.62 × 109
9.44 × 105
6.84 × 105–7.94 × 105
6.89 × 105–8.13 × 105
Re range
<2300
4000–4.23 × 104/2.01 × 105
4.23 × 104/2.01 × 105–2.95 × 106/1.44 × 107
>2.95 × 106/1.44 × 107
Verification result
Wrong
Wrong
Right
Wrong


## 3.1. The Generalization of Quicksand Disaster

The photograph of the borehole and chamber after cleaning is shown in Figure7. It can be clearly seen that the borehole was not damaged in the whole process of quicksand disaster, and the drill pipe remained in the borehole all the time. The profile of the quicksand process is shown in Figure 8:(1)
There is aeolian sand layer with a thickness of 30 m in the ground. The lithology is silty fine sand with a distinct thickness variation. The sand seam constitutes a unified aquifer together with the Sara Wusu group aquifer below.(2)
The Sara Wusu group aquifer is the main aquifer of the coal roof. The main lithology is based on silty fine sand, medium-coarse sandstone, and loam. The units-inflow is between 0.0441 and 0.0569 L·s−1·m−1. The aquifer belongs to low water abundance aquifer, as shown in Table 3.(3)
The borehole diameter is 325 mm in quicksand disaster, including a built-in drill pipe with the diameter of 75 mm. And the borehole and drill pipe have not been damaged in the disaster. Therefore, the cross-sectional areaS, wetted perimeter χ, and hydraulic radius R are 0.0785 m2, 1.256 m, and 0.0625 m, respectively.Figure 7
Borehole and chamber photo after quicksand disaster.Figure 8
Cutaway view of sediment-water flow model.Table 3
Division of water abundance level [32].
Water abundance level
Higher
High
Middle
Low
Units-inflow L·s−1·m−1
>5
1∼5
0.5∼1
<0.5


## 3.2. The Building of Computing Method

Before finding the methods used for quicksand in Longde coal mine, the fluid condition of sand-water mixture flow has to be judged firstly, as different formulas need to be applied for the calculation of the sand-water mixture flow in different fluid conditions. For instance, the flow is easy to obtain with an accurate analytic solution in laminar flow, while the turbulence flow is complicated, which has to be calculated by semiempirical method. The fluid condition can be judged by Reynolds number:(6)Re=ρvdμ,where Re is the Reynolds number, ρ is the density, v is the velocity, d is the diameter, equivalent diameter d = 4R (in this paper), and μ is the dynamic viscosity.In general,Re of laminar flow is below 2,300 and Re of turbulence flow is above 2,300. Nikuradse [23] divided the turbulence flow into effectively smooth flow, transition region, and rough flow.When the flow is in effectively smooth flow, the range of Re is(7)4000<Re<26.98dks8/7,where ks is the roughness height.When the flow is in the transition region, the range of Re is(8)6.98d2ks8/7<Re<4160d2ks8/7.When the flow is in rough flow, the range of Re is(9)Re>4160d2ks8/7.The velocity of sand-water mixture flow is calculated by viscous Bernoulli formula [22]. The equation is(10)z1+p1ρg+α1v122g=z2+p2ρg+α2v222g+hw,where hw is the head loss, z is the third Cartesian coordinate of cross-sectional area, as shown in Figure 9, p is the relative pressure (that the absolute pressure subtracts the atmosphere pressure) of cross-sectional area, and α1 and α2 are correction factors, which are related to the velocity of cross-sectional area. The 1-1 area and 2-2 area (shown as Figure 9) are the same, so α1 is approximately equivalent to α2.Figure 9
Computing model for sediment-water flow.From the continuity equation of fluid mechanics [22], it can be known that there is a relationship between the volume flux Q of 1-1 and 2-2 cross-sectional area as follows:(11)Q1=Q2.For the borehole drilled with steel casing all the while, the sectional areaS of borehole remains unchanged, and v1 is same to v2. The sand-water mixture flow contacts with air in 2-2 area. The relative pressure of 2-2 area is near to zero. The 2-2 area is taken as the base level, and the third Cartesian coordinate z2 is zero. Equation (10) is simplified to(12)hw=z1+p1ρsg.The head losshw is constituted by frictional head loss hf and local head loss hj. From Figure 9, it can be seen that the borehole is a long straight pipe with a constant sectional area. There are only low local head loss hj in both inlet and outlet. Due to the value of hj is much less than hf, hj is negligible. So, hw is approximately equal to hf. The hf can be expressed as [22](13)hf=λldv22g,where l is the distance between 1-1 area and 2-2 area and λ is the flow resistance ratio and it can be expressed with different equations in different flow conditions.In laminar flow, the equation is [33](14)λ=64Re.In effectively smooth flow, the equation is [24](15)1λ=2lgReλ2.51.In the transition region, the equation is [25](16)1λ=−2lgReλ2.51+ks3.7d.In rough flow, the equation is [23](17)1λ=2lg3.7dks.From Equations (15)–(17), it can be known that λ is only related to Re in effectively smooth flow. λ is only related to ks in rough flow. λ is related to both λ and ks in the transition region.Put Equation (13) into Equation (12), and the velocity v can be expressed as follows:(18)v=2gdz1+p1/ρsgλl.In conclusion, the computational model diagram is drawn as Figure10. The velocity will be calculated under different flow conditions after all flow parameters have been determined. Then, Re will be obtained by the calculating result of velocity. Finally, the flow condition will be verified with Re one by one until the right velocity is worked out.Figure 10
Computational model diagram.


## 3.3. The Calculation of Computing Method

The calculating parameters are shown in Table4. The density and viscosity of aeolian sand are obtained by experiments.Table 4
Parameters of semiempirical calculation.
Parameters
Density (kg/m3)
Dynamic viscosity (Pa·s)
Roughness height (mm)
Water pressure (MPa)
Value
1652.89
9.25
0.1–0.4
0.187The obtaining way of the roughness height in Table4 is as follows: researchers have studied on different roughness heights of different rocks, as shown in Table 5. The roughness heights of all rocks are between 0.9 and 16.2 mm. The rock in Longde coal mine is mainly claystone and sandstone. The heights of them are within the range from 1.0 to 5.3 mm without quicksand. But the changing of roughness height is difficult to obtain in the process of quicksand. It can be known that the volume flux of quicksand is about 4,400 to 5,400 m3/h, and the velocity is about 15 to 19 m/s. The roughness height of borehole will decrease gradually with the high-velocity friction. The phenomenon is just like the grinding wheel sanding samples, and the common speed of grinding wheel is 35 m/s (diameter is 200 mm and rpm is 3,150 r/min) [38]. Therefore, the phenomenon in quicksand disaster can be described as follows: the velocity is slow and the roughness height is high at the beginning. Then, the roughness height decreases gradually under the continuous impact of aeolian sand while the velocity increases with the decrease of roughness height. The interaction between velocity and roughness height repeats until they reach a stable value finally. As the velocity is near to the speed of the grinding wheel, the value of roughness height is near to the minimum roughness height of technical standard, as shown in Table 6. In this way, the value of 0.1 mm varies to 0.4 mm, which is the minimum roughness height value taken as the calculation basis of the semiempirical formula in this paper.Table 5
Roughness height and angle of borehole [34–37].
Rock type
Roughness height (mm)
Roughness angle (°)
Seol et al.
Gneiss
1–7
1.1–8.0
Seidel and Collingwood
Claystone, sandstone
1.7–16.2
1.9–18.9
Shale, limestone, etc.
0.9–6.6
1.0–7.6
Lee et al.
Granite, Gneiss
1–4
1.1–4.6
Sandstone, andesite
1–3.5
1.1–4.0
Nam
Claystone
3.6–5.3
4.1–6.1
Limestone
3.2–3.7
3.7–4.2Table 6
Max roughness of technical standard in part country.
Country
Max roughness of technical standard (mm)
Technical standard
China
0.4
GB1031-83
USA
0.1
BS1134-61/ASAB46.1-62
Switzerland
0.2
VSM10321-62
Italy
0.1
UNI13963-60
Japan
0.4
JISB0601-70
Poland
0. 32
PN58/M042-51The obtaining of the water pressure is as follows: there are 4 hydrological boreholes near by the position of quicksand disaster. The water level elevation is recorded aperiodically, as shown in Figure11. From the figure, it can be known that the water level elevation decreases in the underground engineering but has not decreased violently with the quicksand disaster. The decreasing laws of aquifer are similar before and after the quicksand disaster. Therefore, the average value of the observation data of 4 hydrological boreholes on October 29 is selected as the water level elevation value of the computing method. And the water level elevation of aquifer h1 is 1,198.20 m, and the elevation of aquifer roof h2 is 1,179.15 m. The water pressure can be obtained according to the Pascal law:(19)p1=ρwgh1−h2.Figure 11
Water level elevation of Sarah Wusu group aquifer in hydrological borehole.From Equation (19), the water pressure of aquifer roof p1 is 0.187 MPa.Takingz1 equal to 157 m, p1 equal to 0.187 MPa, ρs equal to 1,652.89 kg/m3 in Equation (12), it can be obtained that hw and hf are equal to 168.54 m and the equivalent diameter d = 4R is equal to 0.25 m. Taking the above parameters into Equation (18),(20)v=5.26λ.The velocity is calculated by the simultaneous Equations (6), (14), and (20), Equations (6), (15), and (20), Equations (6), (16), and (20), and Equations (17) and (20) when the sand-water mixture flow is, respectively, in laminar flow, effectively smooth flow, transition region, and rough flow. The calculation results are shown in Table 7, from which it can be known that the sand-water mixture flow belongs to transition region, and the velocity and volume are 15.32 to 17.78 m/s and 310,000 to 350,000 m3, respectively. The volume result is contained in the actual result 310,000 to 380,000 m3.Table 7
Computing result.
Flow condition
Laminar flow
Effectively smooth flow
Transition region
Rough flow
Velocity (m/s)
5.87 × 104
21.14
15.32–17.78
15.43–18.19
Re
2.62 × 109
9.44 × 105
6.84 × 105–7.94 × 105
6.89 × 105–8.13 × 105
Re range
<2300
4000–4.23 × 104/2.01 × 105
4.23 × 104/2.01 × 105–2.95 × 106/1.44 × 107
>2.95 × 106/1.44 × 107
Verification result
Wrong
Wrong
Right
Wrong


## 4. Discussion

As the coal roof aquifer belongs to low water abundance aquifer in Longde coal field, the sand-water mixture flow can be seen as the flow of solid particle of aeolian sand. However, there may be high or higher water abundance aquifers (Table3) even in the desert area. At this point, the volume concentration of sand C in the sand-water mixture will decrease, and the density ρ and viscosity μ in Equation (3) will vary with the change of volume concentration of sand. In order to apply the computing method in this paper expediently, the density ρ and viscosity μ of sand-water mixture with different volume concentration of sand are tested.


### 4.1. Density of Sand-Water Mixture Flow

The density test results of differentC are obtained, as shown in Figure 12. It can be seen that the density of sand-water mixture firstly increases and then decreases with the increase of volume concentration of sand. The max value is 2,095.81 kg/m3 when the volume concentration of sand is 0.78. When the volume concentration of sand is greater than 0.78, the density of sand-water mixture decreases gradually, while the volume concentration of sand tends to 1, the density tends to 1,652.89 kg/m3 which is the dry density of aeolian sand.Figure 12
Density of sand-water mixture change curve with volume ratioC.


### 4.2. The Dynamic Viscosityμ of Sand-Water Mixture Flow

The dynamic viscosityμ of sand-water mixture is tested, as shown in Figure 13. The max viscosity ratio of sand-water mixture is 9.70, and the viscosity ratio of the dry sand is 9.23. The test results are fitted. When S is below 0.49, the fitting equation is(21)μμ0=1−0.9S−0.25,where μ0 is the viscosity of water.Figure 13
Dynamic viscosity ratio change curve with sand concentration of sand-water mixture.WhenS is between 0.49 and 0.78, the fitting equation is(22)μμ0=11.7−1−0.9S−0.25.The aeolian sand is assumed sufficient in the ground. The roughness height is 0.4 mm. The values of water head heightz, pressure p, and drilling diameter d are evaluated according to the approach in Section 3. The calculation results of velocity under the condition of different volume concentrations of sand are obtained, as shown in Figure 14. It can be known that the velocity decreases with the increase of volume concentration of sand, and the overall trend is close to a straight line. It can be explained that both the density in the denominator of Equation (6) and the viscosity in the numerator of Equation (6) increase with the increase of volume concentration of sand. But the increased rate of the viscosity is greater than that of the density; thus, both the Re and velocity tend to decrease with the increase of volume concentration of sand. To sum up, with the same roughness height, the higher the water abundance aquifer is, the smaller the volume concentration of sand will be (Table 3). Therefore, the velocity of sand-water mixture will increase with higher water abundance aquifers, and then larger scale disaster will be caused.Figure 14
Velocity of sand-water flow change curves with volume concentration of sand.


## 4.1. Density of Sand-Water Mixture Flow

The density test results of differentC are obtained, as shown in Figure 12. It can be seen that the density of sand-water mixture firstly increases and then decreases with the increase of volume concentration of sand. The max value is 2,095.81 kg/m3 when the volume concentration of sand is 0.78. When the volume concentration of sand is greater than 0.78, the density of sand-water mixture decreases gradually, while the volume concentration of sand tends to 1, the density tends to 1,652.89 kg/m3 which is the dry density of aeolian sand.Figure 12
Density of sand-water mixture change curve with volume ratioC.


## 4.2. The Dynamic Viscosityμ of Sand-Water Mixture Flow

The dynamic viscosityμ of sand-water mixture is tested, as shown in Figure 13. The max viscosity ratio of sand-water mixture is 9.70, and the viscosity ratio of the dry sand is 9.23. The test results are fitted. When S is below 0.49, the fitting equation is(21)μμ0=1−0.9S−0.25,where μ0 is the viscosity of water.Figure 13
Dynamic viscosity ratio change curve with sand concentration of sand-water mixture.WhenS is between 0.49 and 0.78, the fitting equation is(22)μμ0=11.7−1−0.9S−0.25.The aeolian sand is assumed sufficient in the ground. The roughness height is 0.4 mm. The values of water head heightz, pressure p, and drilling diameter d are evaluated according to the approach in Section 3. The calculation results of velocity under the condition of different volume concentrations of sand are obtained, as shown in Figure 14. It can be known that the velocity decreases with the increase of volume concentration of sand, and the overall trend is close to a straight line. It can be explained that both the density in the denominator of Equation (6) and the viscosity in the numerator of Equation (6) increase with the increase of volume concentration of sand. But the increased rate of the viscosity is greater than that of the density; thus, both the Re and velocity tend to decrease with the increase of volume concentration of sand. To sum up, with the same roughness height, the higher the water abundance aquifer is, the smaller the volume concentration of sand will be (Table 3). Therefore, the velocity of sand-water mixture will increase with higher water abundance aquifers, and then larger scale disaster will be caused.Figure 14
Velocity of sand-water flow change curves with volume concentration of sand.


## 5. Conclusion

An all-purpose computing method was constructed with Bernoulli formula, Darcy–Weisbach formula, semiempirical Nikuradse formula, and Colebrook–White formula to calculate the volume of quicksand. The variation laws of density and dynamic viscosity under different water yield properties were analyzed based on the experimental results. The influence of water yield properties on volume flux was discussed:(1)
It can be seen that the total volume of quicksand is between 310,000 m3 and 380,000 m3, and the quicksand disaster duration is about 70.5 h from the statistical results of sand cleaning and the variation laws of sand subsidence area with time.(2)
The value of water volume flux is much less than that of sand volume flux. The sand-water mixture can be regarded as dry sand in Longde water inrush and quicksand disaster.(3)
The volume of sand between 310,000 m3 and 350,000 m3, i.e., calculated by all-purpose computing method, is contained with the actual result between 310,000 m3 and 380,000 m3.(4)
Both the density and dynamic viscosity of sand-water mixture increase with the increase of sand volume concentration. With the fixed value of water head, roughness height, and diameter of borehole, the velocity decreases with the increase of sand volume concentration, which means that the sand volume concentration will be smaller and the velocity will be faster under a higher abundance level aquifer.



---

*Source: 4842939-2018-10-24.xml*