AbstractThe non‐Darcy flow behavior in unfilled and fully filled rough fractures has been investigated thoroughly for decades. Natural fractures usually be partially filled with porous media due to long‐term internal and external dynamic disturbances and water flow erosion. However, how to evaluate the nonlinear seepage characteristics of the partially filled rough fracture (PFRF) is never scrutinized. In the present study, 2D direct numerical simulations are conducted to investigate the non‐Darcy behavior for flow through 6 PFRF models under different filling conditions (i.e., filling degree, roughness, and hydraulic conductivity capability of the filling medium), which included coupling between the free and seepage flow in the unfilled and filled region, respectively. The results show that the critical non‐Darcy factor E = 0.02 can be used to judge the onset of the non‐Darcy regime in PFRF, which is also the critical point of eddy region (ER) volume mutation. The higher the filling degree and the weaker the hydraulic conductivity of the filled medium, the more significant the non‐Darcy behavior. Furthermore, the variation of ER volume can be divided into three stages with the increasing inertia effect, that is, stable stage (E < 0.02), slow nonlinear growth stage (0.02 ≤ E ≤ 0.1), and rapid linear growth stage (E > 0.1). In addition, a new model for quantitatively evaluating the non‐Darcy behavior of PFRF is developed via gene expression programming based on 735 simulation data set, which is further employed to establish the relationship between the friction factor.
Read full abstract