An improved hydro-mechanical model for 3D rough-walled rock joints considering contact area during shearing

The presence of rock joints has a huge impact on the strength, deformability, and permeability of fractured rock masses by providing surfaces of weakness on which further deformation are more prone to occur, and by serving as the major channels of water flow due to low matrix permeability. In earlier studies, water flow-through rock joints is often assumed as laminar flow between two ideal smooth parallel plates, and the well-known cubic law was derived and widely used to characterise the hydraulic behaviour of rock joints. However, natural rock fractures are irregular and of spatially varied aperture. Thus, the concept of equivalent hydraulic aperture was introduced to modify the cubic law. Numerous works have been conducted on the relation between hydraulic aperture and mechanical aperture incorporating various geometrical parameters. The surface roughness is one of the most sensitive factors for water flow-through rock joints. The aperture distributions are irregular in natural roughed fractures relative to the parallel plates. Barton et al (1985) and Olsson and Barton (2001) incorporated the joint roughness coefficient JRC in evaluating hydraulic aperture and proposed an empirical equation based on fluid flow tests. However, their equations do not conform to the dimensional consistency. Many other researchers (Zimmerman and Bodvarsson, 1996; Xiong et al, 2011; Li et al, 2019) considered aperture distributions in their equations. But the aperture distributions change during shear was always neglected. The contact of rock joints also has a significant impact on the hydro-mechanical behaviour of rock joint, since fluid tends to flow-through a rough joint along connected channels while bypassing the contact obstacles with tortuosity. Nonetheless, previous research mostly considered joint roughness or aperture changes based on 2D joint profiles, while the contact and tortuosity based on 3D surfaces were often overlooked. Only a few works have been made on quantifying the effect of contact area on flow behaviour (Walsh, 1981; Zimmerman et al, 1992; Zimmerman and Bodvarsson, 1996; Yeo, 2001), and mostly ignored contact area change induced by shearing and the associated impact on flow behaviour. Summarily, the quantification of fracture geometry effects on water flow in single rock joints still needs further investigation, especially when considering coupled shear-flow processes. This paper, therefore, aims to study the impacts of contact area and aperture distributions of the single rock joint on the flow behaviour. The evolutions of contact area and aperture distributions induced by shear are considered. Then, an improved model for estimating the hydraulic aperture with the mechanical aperture incorporating the above factors is developed. Finally, a series of shear and shear-flow tests are conducted on artificial joint samples and the model performance is verified by comparing against the results of experimental tests.