Developing an in situ, hydromechanical coupling, true triaxial rock compression tester and investigating the deformation patterns of reservoir bank slopes

Interactions between water and rocks are the main factors affecting the deformation of rock masses on sloped banks by reservoir impoundment. The technology used in laboratory tests of water-rock interaction mechanisms cannot simulate the coupling of water, the rock structure and the initial stress environment. In this work, we develop an in situ hydromechanical true triaxial rock compression tester and apply it to investigate the coupling response of reservoir bank rocks to changing groundwater levels. The tester is composed of a sealed chamber, loader, reactor, and device for measuring deformation, which are all capable of withstanding high water pressures, and a high-precision servo controller. The maximum axial load, lateral load and water pressure are 12 000 kN, 3 000 kN and 3 MPa, respectively. The dimensions of the test specimens are 310 mm×310 mm×620 mm. The test specimens are grey-black basalts with well-developed cracks from the Xiluodu reservoir area. The results show that increasing water pressure promotes axial compression and lateral expansion, while decreasing water pressure causes axial expansion and lateral compression. A water pressure coefficient, K, is introduced as a measure of the hydromechanical coupling effect (expansion or compression) with changing groundwater level. A mechanical tester can be used to perform accurate field tests of the response of wet rocks to hydromechanical coupling. The test results provide new information about the deformation patterns of rock slopes in areas surrounding high dams and reservoirs.Thematic collection: This article is part of the Role of water in destabilizing slopes collection available at: https://www.lyellcollection.org/cc/Role-of-water-in-destabilizing-slopes

An Open-Source Code for Fluid Flow Simulations in Unconventional Fractured Reservoirs

In this article, an open-source code for the simulation of fluid flow, including adsorption, transport, and indirect hydromechanical coupling in unconventional fractured reservoirs is described. The code leverages cutting-edge numerical modeling capabilities like automatic differentiation, stochastic fracture modeling, multicontinuum modeling, and discrete fracture models. In the fluid mass balance equation, specific physical mechanisms, unique to organic-rich source rocks, are included, like an adsorption isotherm, a dynamic permeability-correction function, and an Embedded Discrete Fracture Model (EDFM) with fracture-to-well connectivity. The code is validated against an industrial simulator and applied for a study of the performance of the Barnett shale reservoir, where adsorption, gas slippage, diffusion, indirect hydromechanical coupling, and propped fractures are considered. It is the first open-source code available to facilitate the modeling and production optimization of fractured shale-gas reservoirs. The modular design also facilitates rapid prototyping and demonstration of new models. This article also contains a quantitative analysis of the accuracy and limitations of EDFM for gas production simulation in unconventional fractured reservoirs.

New Leakage Model for Naturally Fractured Formations Considering the Effects of Dual-System Hydro-Mechanical Coupling

Within the existing leakage model accounting for drilling mud loss in naturally-fractured formations, the leak-off velocity is assigned to a fixed value or described by the Cater model, which does not consider the influence of dual-system hydromechanical coupling effects between fracture-wall and fracture-inner systems. The dual-system between the formation and fracture is controlled by the flowing net pressure inside the fracture, which determines the dynamic width of the natural fracture and leak-off velocity. In this study, first, the leak-off velocity under the hydromechanical coupling of the fracture-wall system was obtained based on the coupled governing equations of the solid and liquid phases of the natural fracture-wall, as well as Darcy’s law. Second, the leakage-front invasion velocity, leakage rate, and leakage volume under the hydromechanical coupling of the fracture-inner system were clarified according to the geometric governing of the natural fracture morphology. Finally, the dual-system coupling leakage model was developed considering the continuous equation, while the numerical solution was obtained through a time-step deduction. Results show that at a given time, a greater formation permeability leads to a greater leakage rate and volume, with a smaller leakage front distance. The leakage rate increases with an increase in formation permeability, well bottom differential pressure, and initial width of the natural fracture, while it decreases with an increase in the fracture normal stiffness, yield stress, and plastic viscosity. The new leakage model and numerical method concerning time-step deduction are assessed by solving the issues of fully coupled fracture-wall and fracture-inner systems considering drilling fluid leak-off. The new model may be utilized to simulate various problems related to the invasion of drilling fluids into the fractures, including predicting the dynamic width of natural fracture and borehole ballooning/breathing phenomena.

An Experimental and Numerical Investigation on the Initiation and Interaction of Double Cracks in Rocks under Hydromechanical Coupling

The growth of double cracks is the main factor leading to progressive rock failure under hydromechanical coupling. The initiation modes and interaction behaviors of double cracks were investigated by using laboratory tests, and the influences of water pressure were analyzed. The maximum energy release rate criterion was modified to determine the crack growth characteristics. A numerical model was established and then verified by the test results. Based on the simulation, the distribution of stress fields and key fracture parameters of double cracks was investigated. Then, initiation characteristics and interaction behaviors of parallel and nonparallel cracks were quantitatively analyzed. The results indicate that the increase in water pressure leads to the crack initiation being inclined to the original surfaces and the growth length along the crack fronts tending to be uniform; the small tensile stress zones are formed close to the crack tips, and significant compressive stress zones are formed at both sides of the crack surfaces; stress superposition and interaction occur when crack spacing is less than 2.5a; the interactive weakening effect is mainly present in the inner side (rock bridge zone) of cracks, while a certain degree of interactive enhancement effect exhibits in the outer sides; the cracks are much easier to initiate at the outer wing cracks when the spacing is less than the critical length (0.5a); and cracks with a dip angle of 45° are much easier to initiate at the endpoints of long axis. The research results provide certain theoretical guidance for the safety assessment of underground engineering.