scholarly journals Evaluation of the Non-Darcy Effect of Water Inrush from Karst Collapse Columns by Means of a Nonlinear Flow Model

Water ◽  
2018 ◽  
Vol 10 (9) ◽  
pp. 1234 ◽  
Author(s):  
Yi Xue ◽  
Teng Teng ◽  
Lin Zhu ◽  
Mingming He ◽  
Jie Ren ◽  
...  

Karst collapse columns (KCCs) are naturally formed geological structures that are widely observed in North China. Given their influence on normal mining operations and the progress of mining work, collapse columns pose a hidden danger in coal mining under the influence of manual mining. By communicating often with the aquifer, the water inrush from KCCs poses a serious threat to construction projects. This paper adopts three flow field models, namely, Darcy aquifer laminar flow, Forchheimer flow, and Navier–Stokes turbulent flow, based on the changes in the water inrush flow pattern in the aquifer and laneway, and uses COMSOL Multiphysics software to produce the numerical solutions of these models. As the water inrush flow velocity increases, the Forchheimer flow shows the effect of additional force (inertial resistance) on flow in KCCs, in addition to the effect of viscous resistance. After the joint action of viscous resistance and inertial resistance, the inertial resistance ultimately dominates and gradually changes the water inrush from the KCCs to fluid seepage. Forchheimer flow can comprehensively reflect the nonlinear flow process in the broken rock mass of KCCs, demonstrate the dynamic process from the Darcy aquifer to the final tunnel turbulence layer, and quantitatively show the changes in the flow patterns of the water inrush from KCCs.

Water ◽  
2019 ◽  
Vol 11 (4) ◽  
pp. 855 ◽  
Author(s):  
Yi Xue ◽  
Yang Liu ◽  
Faning Dang ◽  
Jia Liu ◽  
Zongyuan Ma ◽  
...  

Underground fault water inrush is a hydrogeological disaster that frequently occurs in underground mining and tunnel construction projects. Groundwater may pour from an aquifer when disasters occur, and aquifers are typically associated with fractured rock formations. Water inrush accidents are likely to occur when fractured rock masses are encountered during excavation. In this study, Comsol Multiphysics, cross-platform multiphysics field coupling software, was used to simulate the evolution characteristics of water flow in different flow fields of faults and aquifers when water inrush from underground faults occurs. First, the Darcy and Brinkman flow field nonlinear seepage models were used to model the seepage law of water flow in aquifers and faults. Second, the Forchheimer flow field was used to modify the seepage of fluid in fault-broken rocks in the Brinkman flow field. In general, this phenomenon does not meet the applicable conditions of Darcy’s formula. Therefore, the Darcy and Forchheimer flow models were coupled in this study. Simulation results show that flow behavior in an aquifer varies depending on fault permeability. An aquifer near a fault is likely to be affected by non-Darcy flow. That is, the non-Darcy effect zone will either increase or decrease as fault permeability increases or decreases. The fault rupture zone that connects the aquifer and upper roadway of the fault leads to fault water inrush due to the considerably improved permeability of the fractured rock mass.


Processes ◽  
2018 ◽  
Vol 6 (8) ◽  
pp. 120 ◽  
Author(s):  
Weitao Liu ◽  
Jiyuan Zhao ◽  
Ruiai Nie ◽  
Yuben Liu ◽  
Yanhui Du

A coupled thermal-nonlinear hydraulic-mechanical (THM) model for fault water inrush was carried out in this paper to study the water-rock-temperature interactions and predict the fault water inrush. First, the governing equations of the coupled THM model were established by coupling the particle transport equation, nonlinear flow equation, mechanical equation, and the heat transfer equation. Second, by setting different boundary conditions, the mechanical model, nonlinear hydraulic-mechanical (HM) coupling model, and the thermal-nonlinear hydraulic-mechanical (THM) coupling model were established, respectively. Finally, a numerical simulation of these models was established by using COMSOL Multiphysics. Results indicate that the nonlinear water flow equation could describe the nonlinear water flow process in the fractured zone of the fault. The mining stress and the water velocity had a great influence on the temperature of the fault zone. The temperature change of the fault zone can reflect the change of the seepage field in the fault and confined aquifer. This coupled THM model can provide a numerical simulation method to describe the coupled process of complex geological systems, which can be used to predict the fault water inrush induced by coal mining activities.


2021 ◽  
Vol 2021 ◽  
pp. 1-20
Author(s):  
Bin Yang ◽  
Tianhong Yang

For water inrush induced by fracture network flow, the critical velocity of the incipient motion of sand particles was obtained, and the functional relation between critical velocity and particle size was established through a series of tests on the nonlinear flow characteristics of a filling fracture network. The influence of the particle size distribution, hydrodynamic force, and geometric features of the fracture network on the characteristics of particle loss; distribution laws; and water-sand, two-phase migration was also explored. Moreover, the interactions amongst water, movable particles, the surface of the skeleton, and fracture walls, and the formation mechanism of the flow channel were qualitatively analyzed. In addition, the change rules of the mass loss characteristics and porosity of the samples with time were tested successfully. The calculation methods of the permeability and non-Darcy factor of the filling fracture network were also determined.


2014 ◽  
Vol 753 ◽  
pp. 49-79 ◽  
Author(s):  
Ory Schnitzer ◽  
Itzchak Frankel ◽  
Ehud Yariv

AbstractSmoluchowski’s celebrated electrophoresis formula is inapplicable to field-driven motion of drops and bubbles with mobile interfaces. We here analyse bubble electrophoresis in the thin-double-layer limit. To this end, we employ a systematic asymptotic procedure starting from the standard electrokinetic equations and a simple physicochemical interface model. This furnishes a coarse-grained macroscale description where the Debye-layer physics is embodied in effective boundary conditions. These conditions, in turn, represent a non-conventional driving mechanism for electrokinetic flows, where bulk concentration polarization, engendered by the interaction of the electric field and the Debye layer, results in a Marangoni-like shear stress. Remarkably, the electro-osmotic velocity jump at the macroscale level does not affect the electrophoretic velocity. Regular approximations are obtained in the respective cases of small zeta potentials, small ions, and weak applied fields. The nonlinear small-zeta-potential approximation rationalizes the paradoxical zero mobility predicted by the linearized scheme of Booth (J. Chem. Phys., vol. 19, 1951, pp. 1331–1336). For large (millimetre-size) bubbles the pertinent limit is actually that of strong fields. We have carried out a matched-asymptotic-expansion analysis of this singular limit, where salt polarization is confined to a narrow diffusive layer. This analysis establishes that the bubble velocity scales as the $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}2/3$-power of the applied-field magnitude and yields its explicit functional dependence upon a specific combination of the zeta potential and the ionic drag coefficient. The latter is provided to within an $O(1)$ numerical pre-factor which, in turn, is calculated via the solution of a universal (parameter-free) nonlinear flow problem. It is demonstrated that, with increasing field magnitude, all numerical solutions of the macroscale model indeed collapse on the analytic approximation thus obtained. Existing measurements of clean-bubble electrophoresis agree neither with present theory nor with previous models; we discuss this ongoing discrepancy.


RSC Advances ◽  
2018 ◽  
Vol 8 (3) ◽  
pp. 1656-1665 ◽  
Author(s):  
Xian'gang Hou ◽  
Wenhao Shi ◽  
Tianhong Yang

A non-linear flow model that couples three flow types is built based on flow transition to investigate the flow behavior of water inrush induced by KCC.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-15 ◽  
Author(s):  
Shichuan Zhang ◽  
Baotang Shen ◽  
Yangyang Li ◽  
Shengfan Zhou

Water inrush in underground mines is a major safety threat for mining personnel, and it can also cause major damage to mining equipment and result in severe production losses. Water inrush can be attributed to the coalescence of rock fractures and the formation of water channel in rock mass due to the interaction of fractures, hydraulic flow, and stress field. Hence, predicting the fracturing process is the key for investigating the water inrush mechanisms for safe mining. A new coupling method is designed in FRACOD to investigate the mechanisms of water inrush disaster (known as “Luotuoshan accident”) which occurred in China in 2010 in which 32 people died. In order to investigate the evolution processes and mechanisms of water inrush accident in Luotuoshan coal mine, this study applies the recently developed fracture-hydraulic (F-H) flow coupling function to FRACOD and focuses on the rock fracturing processes in a karst collapse column which is a geologically altered zone linking several rock strata vertically formed by the long-term dissolution of the flowing groundwater. The numerical simulation of water inrush is conducted based on the actual geological conditions of Luotuoshan mining area, and various materials with actual geological characteristics were used to simulate the rocks surrounding the coal seam. The influences of several key factors, such as in situ stresses, fractures on the formation, and development of water inrush channels, are investigated. The results indicate that the water inrush source is the Ordovician limestone aquifer, which is connected by the karst collapse column to No. 16 coal seam; the fracturing zone that led to a water inrush occurs in front of the roadway excavation face where new fractures coalesced with the main fractured zone in the karst collapse column.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Wenhao Shi ◽  
Tianhong Yang

A large number of statistics indicate that broken rock mass always transforms into a flowing channel and leads to water inrush disasters in mining engineering, such as fault, karst, and strongly weathered rock mass. During the process of water inrush, the structure of the broken rock mass is constantly changing due to seepage erosion under high-velocity flow. Therefore, it is of vital importance to quantitatively evaluate the flow behavior of the water inrush related to the seepage erosion in order to prevent or reduce the risks. This study described a coupled nonlinear flow model, which couples the high-velocity seepage, the small particle migration, and the evolution of the broken rock mass structure. The model was verified firstly for simulation of nonlinear flow behavior by comparing with the traditional one. Then, the proposed model was used to simulate the evolution of particle migration and seepage properties of the water inrush through broken rock mass by a numerical case. The simulation results generally agree well with the existing experimental results. The simulations indicate that small particle migration causes the unstable characteristics of the seepage and the heterogeneity properties of the broken rock mass, which lead to the nonlinear flow behavior of the water inrush in both time and space. From a different perspective, it also indicates that the proposed model is capable of simulating the interaction of high-velocity seepage, small particle migration, and evolution of broken rock mass structure in the process of water inrush.


2016 ◽  
Vol 2016 ◽  
pp. 1-13 ◽  
Author(s):  
Yue-bao Deng ◽  
Gan-bin Liu ◽  
Rong-yue Zheng ◽  
Kang-he Xie

A nonlinear flow relationship, which assumes that the fluid flow in the soil skeleton obeys the Hansbo non-Darcian flow and that the coefficient of permeability changes with void ratio, was incorporated into Biot’s general consolidation theory for a consolidation simulation of normally consolidated soft ground with or without vertical drains. The governing equations with the coupled nonlinear flow model were presented first for the force equilibrium condition and then for the continuity condition. Based on the weighted residual method, the finite element (FE) formulations were then derived, and an existing FE program was modified accordingly to take the nonlinear flow model into consideration. Comparative analyses using established theoretical solutions and numerical solutions were completed, and the results were satisfactory. On this basis, we investigated the effect of the coupled nonlinear flow on consolidation development.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-8
Author(s):  
Zhengzheng Cao ◽  
Yulou Ren ◽  
Qinting Wang ◽  
Banghua Yao ◽  
Xinchao Zhang

There are many karst collapse columns in coal seam roof in the southern coal field in China, which are different from those in coal seam floor in the northern coal field, due to the stratum characteristics. The karst collapse column in coal seam roof tends to reactivate and conduct water and induce the serious water inrush disaster, when the karst collapse column communicates with the overlying aquifer. In order to reveal the evolution mechanism of water-conducting channel of collapse column in karst mining area of southwest China, the aquifers and water inflow rule in 1908 working face in Qianjin coal mine are analyzed. Besides, the particle size distribution and mineral component of collapse column are researched by the X-ray diffraction test and the screening method, which are the basis for researching the water inrush mechanism in karst collapse column. On this basis, the water inrush of roof collapse column under the influence of mining is researched by establishing the numerical calculation model with the UDEC numerical software. The results show that the water flowing into the 1908 working face comes from the Changxing formation aquifer and Yulongshan formation aquifer above the coal seam, and the proportion of coarse particles and fine particles in collapse column is 89.86% and 10.14%, respectively. With the advance of working face, the water-conducting channel connected the working face with the aquifer, or the surface is formed by collapse pits, karst caves, and collapse column. The research results can be treated as an important basis for the water-preserved mining in southern coal field in China.


Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-18 ◽  
Author(s):  
Banghua Yao ◽  
Zhongwei Chen ◽  
Jianping Wei ◽  
Tianhang Bai ◽  
Shumin Liu

Although the impact of Karst Collapse Pillars (KCPs) on water inrush has been widely recognized and studied, few have investigated the fluid-solid interaction, the particles migration inside KCPs, and the evolution feature of water inrush channels. Moreover, an effective approach to reliably predict the water inrush time has yet to be developed. In this work, a suite of fully coupled governing equations considering the processes of water flow, fracture erosion, and the change of rock permeability due to erosion were presented. The inverse velocity theory was then introduced to predict the water inrush time under different geological and flow conditions. The impact of four different controlling factors on the fracture geometry change, water flow, and inrush time was discussed in detail. The results showed that the inverse velocity theory was capable of predicting the occurrences of water inrush under different conditions, and the time of water inrush had a power relationship with the rock heterogeneity, water pressure, and initial particle concentration and an exponential relationship with the initial fracture apertures. The general approach developed in this work can be extended to other engineering applications such as the tunneling and tailing dam erosion.


Sign in / Sign up

Export Citation Format

Share Document