Special Issue: Fluid Flow in Fractured Porous Media

The fluid flow in fractured porous media plays a significant role in the characteristic/assessment of deep underground reservoirs such as CO2 sequestration [1–3], enhanced oil recovery [4,5] and geothermal energy development [...]

The very first accepted paper of the Special Issue on "Fluid flow in fractured porous media" presented a series of flow tests on granite after high temperature treatment of up to 800 • C. The results showed that higher temperatures induce a larger number of micro-fractures and the larger decrement ratios of both the P-wave velocity and the density of granite. The tested equivalent permeability increases with increasing temperature, following an exponential function.
(2) Jin, Y. et al., Experimental Investigation of the Mechanical Behaviors of Grouted Sand with UF-OA Grouts. [35] This paper investigated the effect of a newly developed grouting material-which is a chemical material-ureaformaldehyde resin mixed with oxalate curing agent (UF-OA)-on the mechanical behaviors of sand specimens through uniaxial compression tests. The results showed that the presence of initial water can decrease both the uniaxial compressive strength (UCS) and the elastic moduli (E) of the grouted specimens; however, the peak strain (ε c ) does not change significantly after 14 days of curing. This paper analyzed the gas concentration and air pressure in the gob and proposed a roof-cutting mechanical model of gob-side entry retaining (GER) with a roadside backfill body (RBB) to determine the key parameters of the GER-TCS. Here, TCS denotes the thin coal seam. They reported that a greater RBB width results in greater roof-cutting resistance and the floor heave of GER increases with increasing RBB width. Finally, the validity of the proposed model is verified by comparing the calculated results with in situ measurements. This paper studied the effect of water soaked height in coal mines on the deformation and crushing characteristics of loose gangue backfill material (LGBM). They concluded that, with the increment of water socked height, both the crushing ratio and the maximum axial strain of LGBM specimens increase. The crushing ratio has a strong linear correlation with the maximum axial strain and a higher water soaked height leads to a lower deformation resistance and a more significant influence on the quality of backfilling. River Embankment. [38] This paper carried out a series of laboratory and in situ experiments to characterize the reinforcement effect of fracture grouting on the Yellow River Embankment. The results showed that the optimum mixing proportion of the compound cement-silicate grout is 70% cement, 25% fly ash and 5% bentonite. The interval distance of the injection hole is suggested to be 1.2 m, in which the lapping effect of fracture grouting veins play an import role. Roadway below a Contiguous Seam Goaf. [39] This paper studied the deformation of rocks surrounding the roadway using numerical simulation and proposed a mechanical model of water-dripping rock using a bolt support. They found that the erosion by water that is dripped on the roadway is substantial and changes significantly over time during roadway excavation. The suggested bolt line spacing is 0.7~0.9 m and the suggested tightening force should exceed 40 kN. This paper investigated the performance of pervious concrete, such as its porosity, pore distribution, permeability coefficient and mechanical properties, with laboratory tests. The main conclusions are that the standard deviation of permeability is less than 0.03 under different hydraulic gradients. The permeability is closely correlated with the volumetric porosity (or effective pore size d 50 ), following polynomial functions. This paper carried out orthogonal and single-factor tests to estimate the critical hydraulic gradient of internal erosion on soil-structure interfaces. They reported that the critical hydraulic gradient is mostly sensitive to degree of compaction, followed by roughness and clay content. The internal erosion is found to develop inside first and then extend to the outside along the interface.

(9) Cui, Y. et al., A New Pseudo Steady-State Constant for a Vertical Well with Finite-Conductivity
Fracture. [42] The authors proposed a new semi-analytical solution for the Pseudo Steady-State constant, which is defined as the difference between the dimensionless wellbore pressure and dimensionless average pressure and is an important parameter for describing Type Curves. They derived a new conductivity-influence function and a normalized conductivity-influence function, which are then utilized to derive the approximate solute of Pseudo Steady-State constant.
(10) Li, Z. et al., The Fracturing Behavior of Tight Glutenites Subjected to Hydraulic Pressure. [43] Numerical simulations were performed to model the fracturing behaviors of tight glutenites, in which the glutenite heterogeneity has been taken into account. The main conclusions are that the hydraulic fractures propagate along the gravels, influenced by stress anisotropy and gravel strength. When the gravel distribution orientation, size and axial ratio are considered, the fractures propagate past the gravel, forming a bypass fracture. Weakening Mechanism under the Effect of Thermal-Hydrological-Chemical Coupling. [44] This paper focuses on the microsctructure evolution of glaubeite under different leaching conditions. Cylindrical specimens with a size of 3 mm in diameter and 10 mm in length were cored from a bedded salt deposit buried at over 1000 m depth in Yunying salt formation, China. The results show that those specimens immersed in a solution at a low temperature generate cracks at the center of the disc. With incremental increase of the temperature of the solution, the crack nucleation appears and the dissolution rate increases significantly. Heating and Cooling with Circulating Water. [45] This paper experimentally studied the damage of granite using acoustic emissions after cyclic heating and cooling with circulating water. They found that the physical and mechanical parameter of granite exhibits exponential relationships with increasing number of cycles. The cumulative number of acoustic emission decreases with the increment of heating and cooling number and the temperature play a key role in the failure mode of samples.
(13) Chen, Y. et al., Coal Anisotropic Sorption and Permeability: An Experimental Study. [46] This paper presented an experimental study on the anisotropic gas adsorption-desorption and permeability of coal. The results show that an increase in the bedding plane angle of the specimen expands the length and area of the contact surface, thereby increasing the speed and quantity of adsorption and desorption. With an increase in the bedding angle, the number of pores and cracks is found to increase together with the volumetric strain. Laboratory experiments were carried out to estimate the relationship between pore fluid pressure and normal deformation of joints. They reported that during fluid injection, the joint opens and host rock dilates. The normal deformation of joints can be well described using the proposed exponential functions and the relative normal deformation of host rock during has a linear relationship with the pore fluid pressure. This paper investigated the creep characteristics of coal measures sandstone influenced by the effect of seepage. A seepage-creep coupling test system was established and scanning electron microscopy techniques were used. The results show that the maximum creep deformation gradually decreases with the increase in confining pressure. For a fixed confining pressure, the creep deformation, creep deformation rate, and accelerated creep deformation rate, increase dramatically with the increment of infiltration pressure. The authors adopted a high-order numerical manifold method to improve the simulation accuracy of fluid flow in heterogeneous porous media. They developed a regular mathematical mesh, which is independent of the physical domain and is very efficient and flexible. The proposed new numerical manifold method model is confirmed to be more accurate and efficient, compared with the traditional numerical manifold method model. This paper investigated the mechanism of hydraulic fracture propagation in coal seams with discontinuous natural fracture networks by using an innovative cohesive element method. The hydraulic fracture network characteristics, as well as the growth of secondary hydraulic fractures and the pore pressure distribution, were analyzed. The results depict that the hydraulic fracture network is spindle-shaped at a high stress difference condition. The ratio of total length of secondary fracture to that of main fracture is 2.11~3.62, which indicates that the secondary fracture plays a dominate role in coal seams. The process of hydraulic fracture propagation in deep coal seams shows very different behaviors with that in hand and/or tight rocks.

(19) Liu, W. et al., A Coupled Thermal-Hydraulic-Mechanical Nonlinear Model for Fault Water
Inrush. [52] A coupled thermal-nonlinear hydraulic-mechanical (THM) model for fault water inrush was proposed in this paper to study water-rock-temperature interactions and to predict the fault water inrush. The governing equations of the THM coupling model coupled the particle transport equation, the nonlinear flow equation, the mechanical equation and the heat transfer equation. Results show that the water velocity increases with the increasing water pressure. Temperature change in the fault zone is subjected to the interaction of the water pressure and the working face advanced distance. The temperature change in the fault zone can reflect the change in seepage field in the fault and confined aquifer, which can be used to predict the fault water inrush induced by coal mining activities. This paper examined the effect of geogrid reinforcement on the deformation and conductivity of bentonite-sand mixtures subjected to differential settlement. When the hydraulic conductivity is 5 × 10 −7 cm/s, the settlement for a 100 mm width bentonite-sand mixture with geogrid reinforcement is 100 mm and increases to 15 mm for the case without geogrid reinforcement. The obtained results are helpful for predicting the deformation and conductivity of a landfill cover system. A strain-based percolation model was proposed to model the coupled process of strain and CH 4 seepage under triaxial compressions. Using a volume-covering method, it is found that the multiscale fracture network is fractal, which has a linear correlation with the porosity of matrix. The proposed percolation model, which links the axial strain to the permeability, is used to interpret the transitions of volumetric strain and CH 4 seepage. A series of experiments were carried out to address the reinforcement by segmented split grouting in soft filling media. They establish a quantitative relationship between mechanical behaviors and hydraulic behaviors of soft filling media both before and after grouting. The results exhibit that three types of veins that are skeleton vein, cross-grid grouting vein and parallel dispersed grouting vein are formed with the increase in the distance from the grouting hole. After the skeleton support and compaction, the segmental split grouting shows a significant reinforcement effect on the soft filling media. This paper performed numerical simulations on stress distribution and fracture development above a protective coal seam in longwall mining in Jincheng, China. Five zones are divided for the rock mass above the protective coal seam, including: pre-mining zone, compression zone, expansion zone, recovery zone and re-compacted zone. The boreholes for gas drainage in the coal seam should be drilled in the volume expansion or the dilation zone, in which the gas concentration is high. This paper analyzed the overlying strata movement and studied the mechanical behaviors of breaking induced disaster of hard and thick magmatic rocks. The main results are that the magmatic rocks control the development of the bed separation. The distribution of the fracture zone alters with the advance of the working face. The gas and water are accumulated below the bed separation space during the breaking of magmatic rocks and tend to rush towards the working face. This paper performed regular shear tests and shear-flow coupled tests on artificial joints that have sawtooth asperities. The influences of joint roughness, normal stress and seepage pressure were systematically analyzed. The results indicate that the shear strength is weakened after a seepage process due to the mutual reinforcing of offsetting and softening effect. The evolution of hydraulic aperture during the tests shows a four-stage characteristic: shear contraction, shear dilation, re-contraction and stability. Fractures. [59] The authors simulated the temporal mixing behavior of solute transport through self-affine fractures. The influences of Hurst exponent and Peclet number were investigated and the mixing was studied by scalar dissipation rate in fractures. For the parallel plate fracture, the Hurst exponent played a negligible role on the temporal evolution of the scalar dissipation rate when the Peclet number was small. The longitudinal scalar dissipation rate overestimated the global scalar dissipation rate for a large Peclet number. This paper numerically studied the failure mechanisms of transversely isotropic rocks using a particle flow modeling method. The results show that three modes are observed in the transversely isotropic rock models, including tensile failure across interfaces, shear failure along interfaces, tensile failure along interfaces. The bounded joint cohesion and joint friction angle significantly influence the uniaxial compressive strengths, whereas the joint coefficient of friction and joint tensile strength slightly influence the uniaxial compressive strengths.
Finally, we would like to thank all the authors for their valuable contributions, all the reviewers for their comments that greatly improved the quality of the papers and all the editors from the editorial office for their substantial support for us to handle this Special Issue.