Advances in Civil Engineering
Volume 2020 (2020), Article ID 1909564, 11 pages
https://doi.org/10.1155/2020/1909564

Research Article

Determining the Seepage Stability of Fractured Coal Rock in the Karst Collapse Pillar

Tianjun Zhang,¹ Mingkun Pang,¹ Xiufeng Zhang,² and Hongyu Pan¹

¹College of Safety Science and Engineering, Xi’an University of Science and Technology, Xi’an 710054, China
²College of Science, Xi’an University of Science and Technology, Xi’an 710054, China

Correspondence should be addressed to Mingkun Pang; 1063862886@qq.com

Received 2 December 2019; Revised 21 May 2020; Accepted 2 July 2020; Published 25 July 2020

Academic Editor: Xiaodong Hu

Copyright © 2020 Tianjun Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The karst collapse pillar (KCP) is a common geological structure in the coal mines of northern China. KCPs contain many fractured coal rocks, which can easily migrate under the action of high-pressure water. The destruction or instability of the cementation structure between the rocks can directly induce coalmine water-inrush accidents. To study the seepage stability of cemented and fractured coal rock under triaxial pressures, a self-designed triaxial seepage testing system was used and the permeability k and non-Darcy factor β of the cemented and fractured coal rock were tested. Furthermore, the 1D non-Darcy seepage equations were used to calculate the evolution criteria of the seepage loss stability. The results show the following: (1) The cemented structure in the KCP under the triaxial pressures can be easily destroyed. The damaged coal and rock body mainly exists in bulk form, and the permeability depends mainly on the effective stress of the particles. (2) The seepage process in the KCP structure is a combination of pore flow, fracture flow, and pipe flow, and the transition of the seepage state is closely related to the change in the magnitude of β. (3) Under the long-term effect of confined underground water, the migration of small fractured particles in the KCP will increase the structural porosity. If the parameter reaches the threshold value, the seepage system will evolve into a pipeline flow state, eventually causing a water-inrush accident.

1. Introduction

The karst collapse pillar (KCP) is a common geological phenomenon in coal mining. It is a special geological structure formed by sedimentation and cementation of fractured coal rock [1, 2]. KCPs commonly occur above a natural aquifer containing a large amount of groundwater [3] and the advancement of the working face in the mine causes the destruction of the bottom layer, which provides conditions for the formation of the KCP water channel. The instability of the seepage state can directly lead to water-inrush accidents [4–6] where groundwater inrush can cause severe damage to the coal mining operation. Therefore, the study of KCP permeability and its stability has attracted worldwide attention [7].

In recent years, research on KCP water inrush focused mainly on the formation cause, evolution mechanism, the KCP mechanical model, and seepage characteristics. Extensive statistical data indicate that KCP water inrush is caused by karst sedimentation in the carbonate rock distribution area [1, 8, 9], and the destruction of the limestone structure enables the formation of a channel for the inflow of high-pressure groundwater from the aquifer into the mine [10, 11]. Furthermore, as the fractured rock undergoes dissolution and erosion, the resulting fine particles will migrate and disperse under the action of high-pressure water, and the pores and cracks in the rock around the KCP will be under mining conditions. As the inflow channel develops and expands, the seepage field and stress field inside and outside the karst collapse column change [12]. Zhu and Wei [13] established a fluid and damage-based model for simulating mining-induced water-inrush behavior considering the influence of a fault and a KCP. Bai et al. [3] established a plug model to describe the seepage flow characteristics of the coal seam floor containing a KCP. Yao et al. [14] established a fluid-solid coupling model of KCPs based on the seepage theory of porous media, using different methods to study water-inrush behavior. Xiang [15] conducted numerical experiments to simulate the groundwater influx of KCP. However, theoretical and experimental studies of KCP seepage characteristics are still lacking. Therefore, we carried out triaxial seepage tests to investigate the seepage state in KCPs, analyzed the seepage characteristics, and developed a new method for estimating seepage loss stability.

The evolution of the KPC instability water-inrush system is shown in Figure 1. When a KCP occurs, the fractured coal-rock body will form a new cementation structure under the long-term dissolution and erosion of the water flow. Considering the gradation structure, the influence of cementation on the seepage stability is not negligible. Based on the above research, we used the Talbot theory to create cemented and fractured coal samples with different grading structures, and considering the influence of multilevel stress on the pore structure of coal, we used a triaxial permeameter to conduct steady-state tests of cemented fractured coal samples. The seepage test explores the effects of the cementation, gradation structure, and effective stress on the seepage stability of the fractured coal rock, and our analysis provides an important theoretical basis for the prediction and early warning of coalmine water-inrush accidents.

Figure 1: Illustration of KCP instability water-inrush system evolution.

2. Materials and Methods

2.1. Sample Preparation

In order to study the seepage stability of the cemented structure between the fractured coal and rock, a series of cemented samples were prepared to test the seepage stability of the structure. The sample preparation process was as follows: First, to produce the specimens, fractured coal rock with different particle sizes was screened, and four particle sizes were selected for gradation. The particle sizes were 0–2.5 mm, 2.5–5.0 mm, 5.0–7.5 mm, and 7.5–10 mm. According to the rock test requirements [16, 17], D of the sample to the maximum particle diameter d should satisfy ; therefore, the maximum particle size in the test was 10 mm. Then, considering the long-term effect of KCP evolution, Wang et al. [18] carried out the experiment by preparing a weak cementation structure; according to the cement content requirements for the fracture testing of the fractured mudstone, the cementing material was selected as standard 425 cement, and the optimum amount of cement per specimen was selected as 50 g. The selection of the particle size is based on the Talbot index theory [18] and the continuous grading theory [19], which defines the Talbot index n to range from 0-1, so this test takesin order, corresponding to the sample number M-1 and M-2 to M-10. The coal particles were mixed with the cementing materials and an appropriate amount of water was added and the mixture stirred until evenly distributed. The ratio of the sample is shown in Table 1.

Table 1: The mass ratio of the sample.

Particle sizes (mm)
Talbot index n

0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0

0∼2.5
87.1
75.8
66.0
57.4
50.0
43.5
37.9
33.0
28.7
25.0

2.5∼5
6.2
11.3
15.3
18.4
20.7
22.5
23.7
24.4
24.9
25.0

5∼7.5
3.9
7.4
10.5
13.3
15.9
18.2
20.2
20.0
23.6
25.0

7.5∼10
2.8
5.5
8.2
10.9
13.4
15.8
18.2
22.6
22.8
25.0

Finally, applying the initial pressure of 2 kN and maintaining, each sample of φ 50  mm × 50  mm was prepared in the mold and stored in the mold in a cool place for 15 days and then dried in room air for 15 days, until the cement had completely solidified. The sample was then removed from the mold, sanded, and dried. A DDL600 rock mechanics test system was combined with a triaxial permeameter to carry out the seepage test on the cemented and fractured coal and coal rock. The triaxial seepage test system of the cemented coal rock is shown in Figure 2.

Figure 2: Cemented samples made at different values of n.

2.2. Test Equipment

From the perspective of rock mechanics, the KCP belongs to a special geological structure, which is mainly due to the dissolution of underground high-pressure water. The fractured coal rock is dissolved in the long term and then cemented. The stability of the KCP depends on the integrity of the cemented structure. In order to further understand the seepage failure mechanism of the cemented structure in the KCP, a set of cemented coal-rock seepage test device was independently developed to carry out the triaxial seepage test on the cemented and fractured coal rock. The main structure of the device includes DDL600 electronic universal testing machine, three-axis percolation test cylinder, lateral pressure pump, and seepage pressure pump, with a computer, data acquisition device, voltage regulator, and so on. The three-axis seepage system of cemented coal rock is shown in Figure 3.

Figure 3: Triaxial seepage test system for cemented fractured coal rock.

Considering the influence of the gradation structure on the seepage characteristics of the fractured coal-rock mass, the sample aggregate prepared in this paper includes four size ranges, and the distribution in the particle size range of each level obeys the Talbot theory [14, 20, 21]; namely,

where Md is the particle mass with particle size less than or equal to D in the sample, kg; Mt is the total mass of the fractured coal rock in the sample, mm; d is the particle size of the fractured coal rock, mm; D is the maximum particle size in the ratio, mm; and n is the Talbot index; the higher the value of n, the higher the content of large particles in the aggregate.

2.3. Test Methods

Structural coal-rock specimens of four different grades were subjected to seepage tests under gradual loading by applying five-stage loading (1, 5, 10, 20, and 30 kN) in the axial direction. The specimens were tested at five levels of seepage pressure (0.5, 1.0, 1.5, 2.0, and 2.5 MPa) under each axial load (Figure 4). The axial load F, the axial displacement , and the loading time t were recorded by the computer system. Then, the steady-state seepage method is used to carry out the test and record the data [22–24]. Previous research on the characteristics of cemented fractured rock showed that when the seepage time reached 3000 s, the decrease in the permeability was about 73%–96% of the total reduction at 6000 s, with the permeability decrease slowing in the later stages of penetration [19]. Therefore, in this test, to eliminate the test error caused by the creep effect of coal rock, the stress holding time at each stage was set at 3000 s.

Figure 4: Schematic diagram of seepage pressure loading.

In order to ensure the reliability of the test data during the test, the percolation test under the same test conditions was repeated 3 times, and the average value was taken as the final test result. The specific test procedures mainly include (1) assembling the coal-rock three-axis seepage test system and measuring the initial height h₀ of the sample; (2) starting the press, setting the initial pressure to 0.2 kN, and recording the initial displacement of the test machine as zero; (3) starting lateral pressure pump, injecting liquid saturation, loading axial pressure 10 kN, and maintaining; (4) keeping the lateral pressure constant, opening the seepage pressure pump, and starting the seepage test; (5) changing the lateral pressure through the hydraulic pump overflow valve, completing different seepage pressure tests, and recording data; (6) changing the value of lateral pressure, completing pressure tests of different lateral pressures, and recording data; (7) pressure relief, unloading, and end of the test.

In the triaxial seepage test, the schematic diagram of loading at each stage is shown in Figure 5. It includes (1) applying axial load, that is, providing axial stress σ₁ by loading of the press; (2) applying lateral pressure; that is, the circumferential pressure is provided by the hydraulic pump; (3) the seepage pressure applied; that is, the seepage pressure p of the stable pressure is supplied through another hydraulic oil pump. Subsequently, the line pressure is adjusted by opening the overflow valve to complete the penetration test under different axial pressures, different lateral pressures, and different seepage pressures.

Figure 5: Schematic diagram of different stages of triaxial seepage.

3. Results and Discussion

3.1. Basic Parameter Test

Each group of coal samples with different cracks is set to 5 levels of lateral pressure, and each stage of lateral pressure is set to 4 levels of seepage pressure. The lateral pressure is selected from the range of 0.5–2.0 MPa. In the test, it is usually ensured that the lateral pressure is 0.2–0.5 MPa larger than the pore pressure [25–27]. If the pore pressure is higher than the lateral pressure, the plastic insulating tape and the heat-shrinkable plastic sleeve will be fractured. Therefore, the pressure of the lateral pressure during the test was set to a total of five grades of 3.0–5.0 MPa, and the seepage flow at each stage of the test was recorded. Through analysis, the results of the percolation parameters are shown in Table 2.

Table 2: Coal sample test parameters.

Sample
M-1
M-2
M-3
M-4
M-5
M-6
M-7
M-8
M-9
M-10

Talbot n
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0

h₀ (mm)
51.2
49.3
51.8
49.8
52.0
52.1
51.9
51.2
50.9
50.3

σb (MPa)
2.043
1.086
1.042
1.069
2.951
3.113
3.447
2.694
3.662
2.017

Critical value
−15.0
−14.9
−13.1
−13.4
−15.1
−15.3
−14.1
−13.2
−13.8
−15.4

h₀ represents the initial height of the sample and σb represents the uniaxial peak compressive strength of the sample.

3.2. Seepage Parameter Analysis

Effective stress is the fundamental factor leading to particle deformation of fractured coal rock. The principle is to simplify the macroscopic three-dimensional force of coal rock, which is equivalent to the effective stress of coal-rock matrix, so as to simplify the complex stress conditions. Therefore, the effective stress can better reflect the stress of the cemented and fractured coal rock under the triaxial stress. The effective stress reflects the stress of the cemented fractured coal rock under triaxial stress. According to the effective stress principle of Terzaghi [28, 29], the effective stress can be written as

where is the hoop stress, equal to the confining pressure σ₃ applied in the test, is the axial stress, equal to σ₁ loaded by the test machine, p₁ is the inlet pressure, p₂ is the outlet pressure (approximately 0 here), the average fracture pressure is , and δ is the Kronecker delta. One has

The effective stress in the test is described by the mean effective stress:

When p₂ = 0, the average effective stress under the triaxial stress condition of the fractured sample is

where σ₃ is the test lateral pressure, in MPa, and σ₁ is the axial stress, in MPa.

The relationship between effective stress and permeability is shown in Figure 6. For all the tested grading structures, the permeability of the coal decreased with increasing effective stress as the permeate seeped into the cemented coal. A negative exponential curve can be fitted to the permeability results corresponding to the following regression equation:

where k represents the permeability, in kg/m³, and σc is the effective stress, in MPa. The curve correlation coefficient R² is greater than 0.9.

(a)

(b)

(c)

(d)

(e)

Figure 6: Relationship between effective stress and permeability of samples. (a) n = 0.1; n = 0.2. (b) n = 0.3; n = 0.4. (c) n = 0.5; n = 0.6. (d). n = 0.7; n = 0.8. (e) n = 0.9; n = 1.0.

The results can be explained as follows: when the lateral pressure is constant, the permeability of the fractured coal rock will decrease with increasing axial pressure. At this stage, σc = 0.26–2.05 MPa and the cemented fractured rock is quickly compacted and enters the fast compaction stage. The maximum decrease of k occurs when n = 0.9. When the axial pressure is further increased, σc > 2.05 MPa, and the decrease of k with increasing axial pressure is less sharp because the cementation between the particles is substantially destroyed, the particles of the coal rock are fractured, and fine particles are generated. With the migration, the pressure-bearing skeleton structure changes and enters the skeleton readjustment stage.

3.3. Non-Darcy Effect Analysis

For porous media, pore channels refer to microchannels that penetrate each other in porous media and are also the main parameters for fluid transport properties of porous media [30]. The void ratio is related to the porous medium skeleton structure and is also susceptible to external loads. For the pore structure of cemented and fractured coal rock, the change of porosity is an indispensable part of studying the permeability characteristics. The porosity φ of the specimen at any stage of the test can be calculated as

where m is the quality of the cemented fractured sample, in kg; h is the initial height of the rock sample, in m; ρ₁ is the density of the fractured coal particles, in kg/m³; is the axial displacement, in m; and A is the cross-sectional area of the sample, in m².

When liquid seeps through the rock sample, p₁ and p₂ are the pore pressures of the inlet and outlet ends of the seepage with respect to the atmosphere. p₂ is the outlet pressure(approximate 0 here), and the pressure gradient after each stage of compression is

is the height of the compressed rock sample, in m. The theoretical and experimental studies on the non-Darcy seepage characteristics of a fractured coal body are based on the Forchheimer empirical formula [31–33], which takes into account the inertia and turbulence effects of the permeate and is a method that accurately reflects the non-Darcy seepage state; the following relationship is satisfied between the pressure gradient and the percolation velocity :

where is the pressure gradient across the sample; μ is the dynamic viscosity of water, in Pa·s; k is the permeability of broken rock under a certain pressure, in m²; is the percolation velocity, in m/s; ρ is the density of water, in kg/m³; and β is a non-Darcy factor, indicating the degree of deviation from the linear law, in m⁻¹. When β = 0, the seepage obeys Darcy's law.

The distribution of β in non-Darcy flow in the tested specimens is shown in Figure 7. The data show that β may be either positive or negative. In the initial stage of the test, the axial pressure and seepage pressure were small and the porosity was large; β is positive and appears in zone II (the stable zone). As the pressure on the sample increases, the deformation and fracturing of the skeleton particles lead to a decrease in permeability. When k ≤ 10⁻¹² m², β is most likely to have a negative value and to appear in region I (the instability zone); when β ≥ 10¹³ m⁻¹, it is also likely to become negative and shift to zone III (the easy mutation zone). In both conditions, the fractured rock sample will have a stable seepage loss. For n = 0.3, 0.4 or 0.5, the β distribution of the skeleton structure is the most dispersed, the stability of the seepage state is poor, and the deviation from the stable zone is the most obvious, indicating that seepage loss is most likely to occur.

Figure 7: Distribution of β-factors in non-Darcy flow in the tested specimens.

3.4. Estimate of Seepage Loss Stability

A large number of experiments have shown that water-inrush or coal and gas outburst accidents caused by seepage loss in mining projects are often caused by sudden changes in the seepage parameters. Therefore, a dynamic method should be used to study the non-Darcy seepage differential dynamic system of fractured coal rock [34–36]. For one-dimensional (1D) passive flow, the 1D non-Darcy seepage dynamics equations of fractured coal rock can be expressed as

For an incompressible fluid and ignoring the acceleration and the volume force F, (10) can be simplified as

The first and second formulas of (11) are separately biased to x and t and then added to obtain

We define the variable Z as

Then, (12) can be simplified as

Among them, the traveling wave method is an important way to solve this kind of nonlinear wave equation. Zakharov equation [37] and the analytical solution of the Kadomtsev-Petviashvili equation [38] (such as the solitary wave solution) are obtained by this method. So, we can set

Substituting (15) into (14),

The stability criterion is then

When , the solution to (17) is

When , the solution to (17) is

When , the solution to (17) is

So, the dynamics of the system depend on the initial conditions; the stability conditions of the non-Darcy seepage differential dynamic system are as follows:

If , the system has two flat steady states and , and the equilibrium state is stable.

If , the system is unstable and it is impossible to reach a steady state regardless of the conditions of the percolation system.The meaning of the formula involved in the above-mentioned formula derivation is shown in Table 3.

In this test, the skeleton structure is cemented and fractured coal rock, and seepage stability was tested under triaxial pressure. The boundary conditions of the equation are ; , and the initial moment velocity distribution can be solved:

Then, , , and the conditions for instability are

Thus, the system instability conditions with the permeability, hydrodynamics, and pressure boundary conditions are obtained:

Here, μ is the kinematic viscosity of the permeate, in Pa·s; ρ is the permeate density, in kg/m³; is the pressure difference across the sample, in MPa; and H is the height of the coal sample, in m.

DTE22 hydraulic oil was used as the permeate in the tests, with a density of ρ = 874 kg/m³ and dynamic viscosity μ = 1.96×10⁻² Pa·s. Taking the first set of samples as an example, H = 51.2 mm, the permeability k₀ obtained by the test is 7.5 × 10⁻¹¹, β₀ = 2.17 × 10¹¹, and △p = 0.5 × 10⁶ Pa. The threshold of instability is then

Taking the logarithm on both sides of (24) yields the critical curve expression:

Similarly, using the above calculation method, the threshold value of the seepage loss for can be calculated. With lgβ and lgk being the coordinate axes, a straight line can be drawn on the direct coordinate system. This line is the dividing line between the stable zone and the unstable zone. The seepage states of the fractured coal body under different gradation structures are shown in Figure 8.

(a)

(b)

(c)

(d)

(e)

Figure 8: The seepage states of stability of fractured coal. (a) n = 0.1; n = 0.2. (b) n = 0.3; n = 0.4. (c) n = 0.5; n = 0.6. (d) n = 0.7; n = 0.8. (e) n = 0.9; n = 1.0.

For the stability analysis in Figure 8, the specimens were subjected to triaxial stress, axial pressure of 1–30 kN was applied to the fractured cemented coal specimens, and the seepage pressure was adjusted within 0.5–2.5 MPa. It was estimated that seepage occurred through 24% of the sample. There were some stages of instability. To ensure the stability of the seepage flow through the fractured coal-rock body, the seepage pressure and lateral pressure need to be maintained within certain ranges to satisfy the relationship:

where , when and Li is the sample that reached the critical value of the instability at the initial state.

The instability curve of the test is

At this point, the threshold of instability is βk² = 3.90 × 10⁻¹⁶ m³.

The analysis results can be explained as follows: (1) In cemented fractured coal rock with different grading structures, when there is a high content of small particles or the rock is finely fractured, blockage of the seepage channel is likely to cause the seepage state to change. (2) Under external loading, the pressure increases and damage to the skeletal structure can lead to the collapse of the pore channels, resulting in a sudden change in gas permeability, which ultimately leads to a stable runoff.

Table 3: Symbol comment list.

Symbol
Physical meaning and interpretation
Symbol
Physical meaning and interpretation

,
Density and initial density, kg/m³

New functions defined,

,
Porosity and initial porosity, %
, ,
Representation function and its various derivatives

Comprehensive compression factor
,
New variables defined,

Acceleration coefficient
,
Permeability and initial permeability, m²

Seepage pressure, MPa
,
Non-Darcy factor and initial factor, m⁻¹

Time, s

Facial strength, N

Seepage velocity, m/s

Integral constant

Seepage displacement, m

Constant

Dynamic viscosity, Pa·s

Height, m

4. Conclusions

To study the evolution mechanism of KCP-induced mine water-inrush events, a triaxial seepage test of fractured coal rock with different grading structures was carried out. The important factors affecting the permeability of fractured coal rock were obtained, and the formulation for estimating the seepage flow was developed. Two criteria were found for instability assessment. The following conclusions were drawn:

(1)The effective stress is a key factor in determining the permeability of fractured coal rock. The permeability of coal with different grades of rock will decrease exponentially as the effective stress increases. When the effective stress is less than 2.05 MPa, the skeleton structure of the coal rock is in the rapid compaction stage which increases the pressure in the skeleton structure; then, the skeleton pressure damage enters the structural readjustment stage.(2)The change in the order of magnitude of the non-Darcy flow factor β is the first criterion for estimating whether the seepage loss is stable. Thus, in the process of particle compaction in the skeleton structure, if the permeability k decreases by an order of 10 or the magnitude of β increases to 6 or more, the percolation system will lose its stability.(3)The value of the seepage parameter of the fractured coal body is the second criterion for judging the stability of the system. Based on our analysis of the mechanism of seepage loss of fractured coal rock, steady flow occurs when is less than the critical value of 3.9 × 10⁻¹⁶ under the initial pressure gradient.

Data Availability

The data used to support the findings of this study have not been made available because of confidentiality reasons.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant nos. 51774234 (Study on the Mechanism of Water-Gas Coupling Fracture Expansion and Ultrasonic Characteristics of Coal Rock in Drilling Holes), 51874234 (Study on Mechanism and Parameter Optimization of Carbon Dioxide Deep Hole Pre-Cracking Blasting), and 51604214.

References

H. Keqiang, Y. Guangming, and L. Yaoru, “Palaeo-karst collapse pillars in northern China and their damage to the geological environments,” Environmental Geology, vol. 58, no. 5, pp. 1029–1040, 2009.

D. Ma, H. Duan, X. Li, Z. Li, Z. Zhou, and T. Li, “Effects of seepage-induced erosion on nonlinear hydraulic properties of broken red sandstones,” Tunnelling and Underground Space Technology, vol. 91, Article ID 102993, 2019.

H. Bai, D. Ma, and Z. Chen, “Mechanical behavior of groundwater seepage in karst collapse pillars,” Engineering Geology, vol. 164, pp. 101–106, 2013.

J. Wu, M. Feng, B. Yu, W. Zhang, X. Ni, and G. Han, “Experimental investigation on dilatancy behavior of water-saturated sandstone,” International Journal of Mining Science and Technology, vol. 28, no. 2, pp. 323–329, 2018.

P. Bukowski, “Water hazard assessment in active shafts in upper silesian coal basin mines,” Mine Water Environment, vol. 30, no. 4, pp. 302–311, 2010.

C. Ping, W. Hua, J. Jin et al., “Experimental study of the fracture failure of sandstone containing hole and fissure under seepage water pressure,” China University of Mining and Technology, vol. 47, pp. 240–246, 2008.

K. He, S. Zhang, F. Wang, and W. Du, “The karst collapses induced by environmental changes of the groundwater and their distribution rules in north China,” Environmental Earth Sciences, vol. 61, no. 5, pp. 1075–1084, 2010.

X. Fang and Y. Fu, “Impact of coal mining on karst water system in north China,” Procedia Earth & Planetary Science, vol. 3, pp. 293–302, 2011.

H. Keqiang, G. Dong, D. Wen, and W. Ronglu, “The effects of karst collapse on the environments in north China,” Environmental Geology, vol. 52, no. 3, pp. 449–455, 2007.

D. Ma, J. Wang, and Z. Li, “Effect of particle erosion on mining-induced water inrush hazard of karst collapse pillar,” Environmental Science and Pollution Research, vol. 26, no. 19, pp. 19719–19728, 2019.

F. Jinyan, L. Xuhang, and Y. Zhiquan, “Numerical simulation study on the mining-induced fracture evolution of steep coal seam,” Journal of China Coal Society, vol. 42, no. 8, pp. 1971–1978, 2017.

H. Kong and L. Wang, “Seepage problems on fractured rock accompanying with mass loss during excavation in coal mines with karst collapse columns,” Arabian Journal of Geosciences, vol. 11, no. 19, 2018.

W. C. Zhu and C. H. Wei, “Numerical simulation on mining-induced water inrushes related to geologic structures using a damage-based hydromechanical model,” Environmental Earth Sciences, vol. 62, no. 1, pp. 43–54, 2011.

B. H. Yao, H. B. Bai, J. P. Wei, and D. Ma, “Numerical study on seepage property of karst collapse pillars under particle migration,” Computer Modeling in Engineering and Sciences, vol. 91, no. 2, pp. 81–100, 2013.

Y. Xiang, “The mechanics model on flooding process of collapse pillar,” Coal Geology Exploration, vol. 21, no. 5, pp. 36–39, 1993.

H. Li, “Research on seepage properties and pore structure of the roof and floor strata in confined water-rich coal seams: taking the xiaojihan coal mine as an example,” Advances in Civil Engineering, vol. 2018, Article ID 9483637, pp. 1–8, 2018.

T. Zhang, L. Zhang, S. Li, J.-L. Liu, H.-Y. Pan, and X. Ji, “Wave velocity and power spectral density of hole-containing specimens with different moisture content under uniaxial compression,” Energies, vol. 11, no. 11, p. 3166, 2018.

G. Wang, P. Wang, Y. Guo, and W. Li, “A novel true triaxial apparatus for testing shear seepage in gas-solid coupling coal,” Geofluids, vol. 2018, Article ID 2608435, pp. 1–9, 2018.

B. Yu, Z. Chen, Q. Ding, and L. Wang, “Non-darcy flow seepage characteristics of saturated broken rocks under compression with lateral constraint,” International Journal of Mining Science and Technology, vol. 26, no. 6, pp. 1145–1151, 2016.

S. Ou, L. Wang, P. Wang, Z. Wang, J. Huang, and D. Zhou, “Numerical analysis of seepage flow characteristic of collapse pillar under the influence of mining,” International Journal of Mining Science and Technology, vol. 23, no. 2, pp. 237–244, 2013.

Ji Xiong, Yu Zhang, N. Zhou, and M. Li, “Pore pressure evolution and mass loss of broken gangue during the seepage,” Royal Society open science, vol. 5, no. 10, Article ID 180307, 2018.

L. Wang, Z. Chen, H. Kong, and X. Ni, “An experimental study of the influence of seepage pressure and initial porosity on variable mass seepage for broken mudstone,” Journal of Mining & Safety Engineering, vol. 31, no. 3, pp. 462–468, 2014.

A. V. Garcia, R. M. Rached, and J. C. Santamarina, “Large-scale true triaxial apparatus for geophysical studies in fractured rock,” Geotechnical Testing Journal, vol. 41, no. 4, Article ID 20170144, pp. 821–829, 2018.

G. Zhou, L. Qiu, W. Zhang, and J. Xue, “Simulation analysis on water’s micro seepage laws under different pressure gradients using computed tomography method,” Advances in Civil Engineering, vol. 2018, Article ID 3401341, pp. 1–26, 2018.

L. Shi and R. N. Singh, “Study of mine water inrush from floor strata through faults,” Mine Water and Environment, vol. 20, no. 3, pp. 140–147, 2001.

W. Gang, W. Pengfei, G. Yangyang, and W. Li, “A novel true triaxial apparatus for testing shear seepage in gas-solid coupling coal,” Geofluids, vol. 2018, Article ID 2608435, pp. 1–9, 2018.

W. Du, Y. Zhang, X. Meng, W. Zhang, and W. Li, “Deformation and seepage characteristics of gas-containing coal under true tri-axial stress,” Arabian Journal of Geosciences, vol. 11, no. 9, 2018.

D. Ma, M. Rezania, H.-S. Yu, and H.-B. Bai, “Variations of hydraulic properties of granular sandstones during water inrush: effect of small particle migration,” Engineering Geology, vol. 217, pp. 61–70, 2017.

L Bo, Y. Liang, L. Zhang, and Q. Zou, “Breakage law and fractal characteristics of broken coal and rock masses with different mixing ratios during compaction,” Energy Science & Engineering, vol. 7, no. 3, pp. 1000–1015, 2019.

Z. Kai, Z. Boyang, L. Jiangfeng, D. Ma, and H. Bai, “Experiment on seepage property and sand inrush criterion for granular rock mass,” Geofluids, vol. 2017, pp. 1–10, 2017.

E. Celik and L. Hoang, “Maximum estimates for generalized forchheimer flows in heterogeneous porous media,” Journal of Differential Equations, vol. 262, no. 3, Article ID 9352618, 2017.

B. Zhang and Z. Lin, “A computing method for sand inrush quantity through a borehole in Longde coal mine,” Advance of Civil Engineering, vol. 2018, Article ID 4842939, 11 pages, 2018.

S. J. Peng, J. Xu, H. W. Yang, and D. Liu, “Experimental study on the influence mechanism of gas seepage on coal and gas outburst disaster,” Safety Science, vol. 50, no. 4, pp. 816–821, 2012.

W. Zhu, Q. Qi, Q. Ma, J. Deng, M. Yue, and Y. Liu, “Unstable seepage modeling and pressure propagation of shale gas reservoirs,” Petroleum Exploration and Development, vol. 43, no. 2, pp. 285–292, 2016.

M. Li, J. Zhang, W. Zhang, A. Li, and W. Yin, “Experimental investigation of water-inrush risk based on permeability evolution in coal mine and backfill prevention discussion,” Geofluids, vol. 2019, Article ID 3920414, pp. 1–9, 2019.

Z. Zhou, X. Cai, D. Ma et al., “Water saturation effects on dynamic fracture behavior of sandstone,” International Journal of Rock Mechanics and Mining Sciences, vol. 114, pp. 46–61, 2019.

C. Yong and L. Biao, “New exact travelling wave solutions for generalized zakharov-kuzentsov equations using general projective riccati equation method,” Communications in Theoretical Physics, vol. 41, no. 1, pp. 1–6, 2004.

H. Hua, J. Yudang, and Z. Jinmei, “Some special types of multisolution of the modified kadomtsev-petviashvili equation,” Communication in Theoretical Dhysics, vol. 40, no. 3, pp. 262–264, 2003.