Research Article  Open Access
Tianjun Zhang, Mingkun Pang, Xiufeng Zhang, Hongyu Pan, "Determining the Seepage Stability of Fractured Coal Rock in the Karst Collapse Pillar", Advances in Civil Engineering, vol. 2020, Article ID 1909564, 11 pages, 2020. https://doi.org/10.1155/2020/1909564
Determining the Seepage Stability of Fractured Coal Rock in the Karst Collapse Pillar
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 highpressure water. The destruction or instability of the cementation structure between the rocks can directly induce coalmine waterinrush accidents. To study the seepage stability of cemented and fractured coal rock under triaxial pressures, a selfdesigned triaxial seepage testing system was used and the permeability k and nonDarcy factor β of the cemented and fractured coal rock were tested. Furthermore, the 1D nonDarcy 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 longterm 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 waterinrush 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 waterinrush 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 highpressure 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 highpressure 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 damagebased model for simulating mininginduced waterinrush 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 fluidsolid coupling model of KCPs based on the seepage theory of porous media, using different methods to study waterinrush 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 waterinrush system is shown in Figure 1. When a KCP occurs, the fractured coalrock body will form a new cementation structure under the longterm 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 steadystate 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 waterinrush accidents.
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 longterm 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 01, so this test takesin order, corresponding to the sample number M1 and M2 to M10. 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.

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.
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 highpressure 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 coalrock 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, threeaxis percolation test cylinder, lateral pressure pump, and seepage pressure pump, with a computer, data acquisition device, voltage regulator, and so on. The threeaxis seepage system of cemented coal rock is shown in Figure 3.
Considering the influence of the gradation structure on the seepage characteristics of the fractured coalrock 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 M_{d} is the particle mass with particle size less than or equal to D in the sample, kg; M_{t} 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 coalrock specimens of four different grades were subjected to seepage tests under gradual loading by applying fivestage 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 steadystate 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.
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 coalrock threeaxis seepage test system and measuring the initial height h_{0} 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 σ_{1} 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.
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 heatshrinkable 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.
 
h_{0} 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 threedimensional force of coal rock, which is equivalent to the effective stress of coalrock 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 aswhere is the hoop stress, equal to the confining pressure σ_{3} applied in the test, is the axial stress, equal to σ_{1} loaded by the test machine, p_{1} is the inlet pressure, p_{2} 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_{2} = 0, the average effective stress under the triaxial stress condition of the fractured sample iswhere σ_{3} is the test lateral pressure, in MPa, and σ_{1} 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^{3}, and σ_{c} is the effective stress, in MPa. The curve correlation coefficient R^{2} is greater than 0.9.
(a)
(b)
(c)
(d)
(e)
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 pressurebearing skeleton structure changes and enters the skeleton readjustment stage.
3.3. NonDarcy 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 aswhere m is the quality of the cemented fractured sample, in kg; h is the initial height of the rock sample, in m; ρ_{1} is the density of the fractured coal particles, in kg/m^{3}; is the axial displacement, in m; and A is the crosssectional area of the sample, in m^{2}.
When liquid seeps through the rock sample, p_{1} and p_{2} are the pore pressures of the inlet and outlet ends of the seepage with respect to the atmosphere. p_{2} 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 nonDarcy 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 nonDarcy 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^{2}; is the percolation velocity, in m/s; ρ is the density of water, in kg/m^{3}; and β is a nonDarcy factor, indicating the degree of deviation from the linear law, in m^{−1}. When β = 0, the seepage obeys Darcy's law.
The distribution of β in nonDarcy 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^{−12} m^{2}, β is most likely to have a negative value and to appear in region I (the instability zone); when β ≥ 10^{13} m^{−1}, 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.
3.4. Estimate of Seepage Loss Stability
A large number of experiments have shown that waterinrush 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 nonDarcy seepage differential dynamic system of fractured coal rock [34–36]. For onedimensional (1D) passive flow, the 1D nonDarcy 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 KadomtsevPetviashvili equation [38] (such as the solitary wave solution) are obtained by this method. So, we can set
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 nonDarcy 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 abovementioned 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^{3}; 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^{3} and dynamic viscosity μ = 1.96×10^{−2} Pa·s. Taking the first set of samples as an example, H = 51.2 mm, the permeability k_{0} obtained by the test is 7.5 × 10^{−11}, β_{0} = 2.17 × 10^{11}, and △p = 0.5 × 10^{6} 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)
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 coalrock 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^{2} = 3.90 × 10^{−16} m^{3}.
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.

4. Conclusions
To study the evolution mechanism of KCPinduced mine waterinrush 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 nonDarcy 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^{−16} 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 WaterGas Coupling Fracture Expansion and Ultrasonic Characteristics of Coal Rock in Drilling Holes), 51874234 (Study on Mechanism and Parameter Optimization of Carbon Dioxide Deep Hole PreCracking Blasting), and 51604214.
References
 H. Keqiang, Y. Guangming, and L. Yaoru, “Palaeokarst collapse pillars in northern China and their damage to the geological environments,” Environmental Geology, vol. 58, no. 5, pp. 1029–1040, 2009. View at: Publisher Site  Google Scholar
 D. Ma, H. Duan, X. Li, Z. Li, Z. Zhou, and T. Li, “Effects of seepageinduced erosion on nonlinear hydraulic properties of broken red sandstones,” Tunnelling and Underground Space Technology, vol. 91, Article ID 102993, 2019. View at: Publisher Site  Google Scholar
 H. Bai, D. Ma, and Z. Chen, “Mechanical behavior of groundwater seepage in karst collapse pillars,” Engineering Geology, vol. 164, pp. 101–106, 2013. View at: Publisher Site  Google Scholar
 J. Wu, M. Feng, B. Yu, W. Zhang, X. Ni, and G. Han, “Experimental investigation on dilatancy behavior of watersaturated sandstone,” International Journal of Mining Science and Technology, vol. 28, no. 2, pp. 323–329, 2018. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 D. Ma, J. Wang, and Z. Li, “Effect of particle erosion on mininginduced water inrush hazard of karst collapse pillar,” Environmental Science and Pollution Research, vol. 26, no. 19, pp. 19719–19728, 2019. View at: Publisher Site  Google Scholar
 F. Jinyan, L. Xuhang, and Y. Zhiquan, “Numerical simulation study on the mininginduced fracture evolution of steep coal seam,” Journal of China Coal Society, vol. 42, no. 8, pp. 1971–1978, 2017. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 W. C. Zhu and C. H. Wei, “Numerical simulation on mininginduced water inrushes related to geologic structures using a damagebased hydromechanical model,” Environmental Earth Sciences, vol. 62, no. 1, pp. 43–54, 2011. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 Y. Xiang, “The mechanics model on flooding process of collapse pillar,” Coal Geology Exploration, vol. 21, no. 5, pp. 36–39, 1993. View at: Google Scholar
 H. Li, “Research on seepage properties and pore structure of the roof and floor strata in confined waterrich coal seams: taking the xiaojihan coal mine as an example,” Advances in Civil Engineering, vol. 2018, Article ID 9483637, pp. 1–8, 2018. View at: Publisher Site  Google Scholar
 T. Zhang, L. Zhang, S. Li, J.L. Liu, H.Y. Pan, and X. Ji, “Wave velocity and power spectral density of holecontaining specimens with different moisture content under uniaxial compression,” Energies, vol. 11, no. 11, p. 3166, 2018. View at: Publisher Site  Google Scholar
 G. Wang, P. Wang, Y. Guo, and W. Li, “A novel true triaxial apparatus for testing shear seepage in gassolid coupling coal,” Geofluids, vol. 2018, Article ID 2608435, pp. 1–9, 2018. View at: Publisher Site  Google Scholar
 B. Yu, Z. Chen, Q. Ding, and L. Wang, “Nondarcy 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 A. V. Garcia, R. M. Rached, and J. C. Santamarina, “Largescale true triaxial apparatus for geophysical studies in fractured rock,” Geotechnical Testing Journal, vol. 41, no. 4, Article ID 20170144, pp. 821–829, 2018. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 W. Gang, W. Pengfei, G. Yangyang, and W. Li, “A novel true triaxial apparatus for testing shear seepage in gassolid coupling coal,” Geofluids, vol. 2018, Article ID 2608435, pp. 1–9, 2018. View at: Google Scholar
 W. Du, Y. Zhang, X. Meng, W. Zhang, and W. Li, “Deformation and seepage characteristics of gascontaining coal under true triaxial stress,” Arabian Journal of Geosciences, vol. 11, no. 9, 2018. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 M. Li, J. Zhang, W. Zhang, A. Li, and W. Yin, “Experimental investigation of waterinrush risk based on permeability evolution in coal mine and backfill prevention discussion,” Geofluids, vol. 2019, Article ID 3920414, pp. 1–9, 2019. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 C. Yong and L. Biao, “New exact travelling wave solutions for generalized zakharovkuzentsov equations using general projective riccati equation method,” Communications in Theoretical Physics, vol. 41, no. 1, pp. 1–6, 2004. View at: Publisher Site  Google Scholar
 H. Hua, J. Yudang, and Z. Jinmei, “Some special types of multisolution of the modified kadomtsevpetviashvili equation,” Communication in Theoretical Dhysics, vol. 40, no. 3, pp. 262–264, 2003. View at: Google Scholar
Copyright
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.