Numerical Research of Fluid Flow and Solute Transport in Rough Fractures under Different Normal Stress

The effects of roughness and normal stress on hydraulic properties of fractures are significant during the coupled shear flow test. Knowing the laws of fluid flow and solute transport in fractures is essential to ensure the nature and safety of geological projects. Although many experiments and numerical simulations of coupled shear flow test have been conducted, there is still a lack of research on using the full Navier-Stokes (N-S) equation to solve the real flow characteristics of fluid in three-dimensional rough fractures. The main purpose of this paper is to study the influence of roughness and normal stress on the fluid flow and solute transport through fractures under the constant normal stiffness boundary condition. Based on the corrected successive random addition (SRA) algorithm, fracture surfaces with different roughness expressed by the Hurst coefficient ( H ) were generated. By applying a shear displacement of 5 mm, the sheared fracture models with normal stresses of 1 MPa, 3 MPa, and 5 MPa were obtained, respectively. The hydraulic characteristics of three-dimensional fractures were analyzed by solving the full N-S equation. The particle tracking method was employed to obtain the breakthrough curves based on the calculated flow field. The numerical method was verified with experimental results. It has been found that, for the same normal stress, the smaller the fracture H value is (i.e., more tough the fracture is), the larger the mechanical aperture is. The ratio of hydraulic aperture to mechanical aperture ( e h / e m ) decreases with the increasing of normal stress. The smaller the H value, the effect of the normal stress on the ratio e h / e m is more significant. The variation of transmissivity of fractures with the flow rate exhibits similar manner with that of e h / e m . With the normal stress and H value increasing, the mean velocity of particles becomes higher and more particles move to the outlet boundary. The dispersive transport behavior becomes obvious when normal stress is larger.

Evolution of Fracture Aperture in Quartz Sandstone under Hydrothermal Conditions: Mechanical and Chemical Effects

Fractures efficiently affect fluid flow in geological formations, and thereby determine mass and energy transport in reservoirs, which are not least exploited for economic resources. In this context, their response to mechanical and thermal changes, as well as fluid–rock interactions, is of paramount importance. In this study, a two-stage flow-through experiment was conducted on a pure quartz sandstone core of low matrix permeability, containing one single macroscopic tensile fracture. In the first short-term stage, the effects of mechanical and hydraulic aperture on pressure and temperature cycles were investigated. The purpose of the subsequent intermittent-flow long-term (140 days) stage was to constrain the evolution of the geometrical and hydraulic fracture properties resulting from pressure solution. Deionized water was used as the pore fluid, and permeability, as well as the effluent Si concentrations, were systematically measured. Overall, hydraulic aperture was shown to be significantly less affected by pressure, temperature and time, in comparison to mechanical aperture. During the long-term part of the experiment at 140 °C, the effluent Si concentrations likely reached a chemical equilibrium state within less than 8 days of stagnant flow, and exceeded the corresponding hydrostatic quartz solubility at this temperature. This implies that the pressure solution was active at the contacting fracture asperities, both at 140 °C and after cooling to 33 °C. The higher temperature yielded a higher dissolution rate and, consequently, a faster attainment of chemical equilibrium within the contact fluid. X-ray µCT observations evidenced a noticeable increase in fracture contact area ratio, which, in combination with theoretical considerations, implies a significant decrease in mechanical aperture. In contrast, the sample permeability, and thus the hydraulic fracture aperture, virtually did not vary. In conclusion, pressure solution-induced fracture aperture changes are affected by the degree of time-dependent variations in pore fluid composition. In contrast to the present case of a quasi-closed system with mostly stagnant flow, in an open system with continuous once-through fluid flow, the activity of the pressure solution may be amplified due to the persistent fluid-chemical nonequilibrium state, thus possibly enhancing aperture and fracture permeability changes.

Prediction of contaminant transport in fractured carbonate aquifer types: a case study of the Permian Magnesian Limestone Group (NE England, UK)

Viruses and bacteria which are characterized by finite lives in the subsurface are rapidly transported via fractures and cavities in fractured and karst aquifers. Here, we demonstrate how the coupling of a robust outcrop characterization and hydrogeophysical borehole testing is essential for prediction of contaminant velocities and hence wellhead protection areas. To show this, we use the dolostones of the Permian Magnesian Limestone aquifer in NE England, where we incorporated such information in a groundwater flow and particle tracking model. Within this aquifer, flow in relatively narrow (mechanical aperture of ~ 10−1–1 mm) fractures is coupled with that in pipe cavities (~ 0.20-m diameter) following normal faults. Karstic cavities and narrow fractures are hydraulically very different. Thus, the solutional features are represented within the model by a pipe network (which accounts for turbulence) embedded within an equivalent porous medium representing Darcian flowing fractures. Incorporation of fault conduits in a groundwater model shows that they strongly influence particle tracking results. Despite this, away from faulted areas, the effective flow porosity of the equivalent porous medium remains a crucial parameter. Here, we recommend as most appropriate a relatively low value of effective porosity (of 2.8 × 10−4) based on borehole hydrogeophysical testing. This contrasts with earlier studies using particle tracking analyses on analogous carbonate aquifers, which used much higher values of effective porosity, typically ~ 102 times higher than our value, resulting in highly non-conservative estimates of aquifer vulnerability. Low values of effective flow porosities yield modelled flow velocities ranging from ~ 100 up to ~ 500 m/day in un-faulted areas. However, the high fracturing density and presence of karstic cavities yield modelled flow velocities up to ~ 9000 m/day in fault zones. The combination of such flow velocities along particle traces results in 400-day particle traces up to 8-km length, implying the need for large well protection areas and high aquifer vulnerability to slowly degrading contaminants.

Development of a New Method for the Quantitative Generation of an Artificial Joint Specimen with Specific Geometric Properties

A rock joint is a planar discontinuity that has significant influence on the mechanical and hydraulic characteristics of rock mass. Laboratory experiments are often conducted on a joint to investigate and provide fundamental information for rock mass analysis. Although joint roughness and mechanical aperture exert great effects on the experimental results, controlling them in quantitative manner is quite complicated and consumptive in terms of specimen preparation. A new and simple method for the quantitative generation of the joint specimen was proposed in this study. Based on random midpoint displacement method, a joint specimen with a void space inside was generated. Parametric studies for the roughness and mechanical aperture were carried out, and as a result, the two joint properties could be controlled by manipulating input parameters of random midpoint displacement method. In order to validate the proposed method, two joint specimens, which had different levels of roughness and aperture, were generated and printed. Surface coordinates of the specimens were obtained by a 3D laser scanner, and calculated to make a comparison between the target values and the estimated values. Results showed that the method was capable of generating joint specimens with satisfactory precision.

Nonlinear Flow Characteristics of a System of Two Intersecting Fractures with Different Apertures

The nonlinear flow regimes of a crossed fracture model consisting of two fractures have been investigated, in which the influences of hydraulic gradient, surface roughness, intersecting angle, and scale effect have been taken into account. However, in these attempts, the aperture of the two crossed fractures is the same and effects of aperture ratio have not been considered. This study aims to extend their works, characterizing nonlinear flow through a system of two intersecting fractures with different apertures. First, three experiment models with two fractures having different apertures were established and flow tests were carried out. Then, numerical simulations by solving the Navier-Stokes equations were performed and the results compared with the experiment results. Finally, the effects of fracture aperture on the critical pressure difference and the ratio of hydraulic aperture to mechanical aperture were systematically analyzed. The results show that the numerical simulation results agree well with those of the fluid flow tests, which indicates that the visualization techniques and the numerical simulation code are reliable. With the increment of flow rate, the pressure difference increases first linearly and then nonlinearly, which can be best fitted using Forchheimer’s law. The two coefficients in Forchheimer’s law decrease with the increasing number of outlets. When increasing fracture aperture from 3 mm to 5 mm, the critical pressure difference increases significantly. However, when continuously increasing fracture aperture from 5 mm to 7 mm, the critical pressure difference changes are negligibly small. The ratio of hydraulic aperture to mechanical aperture decreases more significantly for a fracture that has a larger aperture. Increasing fracture aperture from 5 mm to 7 mm, that has a negligibly small effect on the critical pressure difference will however significantly influence the ratio of hydraulic aperture to mechanical aperture.

Numerical Investigation on Hydraulic Properties of Artificial-Splitting Granite Fractures during Normal and Shear Deformations

This study explores the effects of normal loading and shearing on hydraulic properties in roughness-walled rock fractures. The geometries of five fractures were measured by the 3D scanning technology. The flow simulation was performed for rough rock fractures with large displacements during normal loading and shearing by finite volume method (FVM). The results demonstrate that the deformation of fracture with increasing normal stress and shear causes nonuniform changes in void space geometry and further influences fracture permeability. Associated with normal displacement are an increase in contact area and a decrease in mechanical aperture. The transmissivity is decreasing by 3 orders of magnitude response to applied normal displacement values of 0.0 mm to 1.8 mm. In contrast, an increase in mechanical aperture and contact ratio that occurs with increasing shear displacement values of 0.0 mm to 4.0 mm is associated with decreasing distinctly transmissivity by 1.5–2 orders of magnitude. Based on the numerical results, an empirical equation is proposed to evaluate the effects of contact area and roughness of fracture on the hydraulic aperture. The good agreement between numerical results and the predicted results by the new model indicates that the proposed model is capable of estimating the hydraulic aperture of rock fractures through parametric analyses, compared with other published models from available literature. In addition, the new model succeeds in predicting the transmissivity in Develi and Babadagli (2014) water flow experiments.

Effect of Shear Stress and Displacement on the Hydraulic Properties of a Marble Fracture with Sand

A test apparatus was developed to study shear-flow coupling properties for a marble fracture with sand. Experimental tests include shear stress-flow coupling process and shear displacement-flow coupling process. The shear stress-flow coupling test results, revealed that mechanical aperture keeps constant during shear, and that hydraulic aperture and conductivity decrease almost linearly with increasing shear stress. The shear displacement-flow coupling test results showed that hydraulic behavior under shear displacement includes two phases. During the first about 0.05mm shear displacement, hydraulic aperture and conductivity decrease with increasing shear displacement, and mechanical aperture remains constant or increases only by several microns of magnitude. After shear displacement is up to about 0.05mm, the fracture dilates remarkably, and hydraulic aperture and conductivity increase with shear displacement increasing.