The effect of model size on fluid flow through fractal rough fractures under shearing is investigated using a numerical simulation method. The shear behavior of rough fractures with selfaffine properties was described using the analytical model, and the aperture fields with sizes ¶ Corresponding author. This is an Open Access article published by World Scientific Publishing Company. It is distributed under the terms of the Creative Commons Attribution 4.0 (CC-BY) License. Further distribution of this work is permitted, provided the original work is properly cited. varying from 25 to 200 mm were extracted under shear displacements up to 20 mm. Fluid flow through fractures in the directions both parallel and perpendicular to the shear directions was simulated by solving the Reynolds equation using a finite element code. The results show that fluid flow tends to converge into a few main flow channels as shear displacement increases, while the shapes of flow channels change significantly as the fracture size increases. As the model size increases, the permeability in the directions both parallel and perpendicular to the shear direction changes significantly first and then tends to move to a stable state. The size effects on the permeability in the direction parallel to the shear direction are more obvious than that in the direction perpendicular to the shear direction, due to the formation of contact ridges and connected channels perpendicular to the shear direction. The variances of the ratio between permeability in both directions become smaller as the model size increases and then this ratio tends to maintain constant after a certain size, with the value mainly ranging from 1.0 to 3.0.
1840001-1