Numerical simulations for the effects of normal loading on particle transport in rock fractures during shear

Fluid flow and tracer transport in a single rock fracture during shear is investigated using the finite element method (FEM) and streamline particle tracking, considering evolutions of aperture and transmissivity with shear displacement histories under different normal stresses, based on laboratory tests. The distributions of fracture aperture and its evolution during shear were calculated from the initial aperture fields, based on the laser-scanned surface roughness of feature replicas of rock fracture specimens, and shear dilations measured during the coupled shear-flow tests in laboratory. The coupled shear-flow tests were performed under two levels of constant normal loading (CNL). A special algorithm for treating the contact areas as zero-aperture elements was used to produce more accurate flow field simulations using FEM. The simulation results agreed well with the flow rate data obtained from the laboratory tests, showing complex histories of fracture aperture and tortuous flow channels with changing normal stresses and increasing shear displacements for the flow parallel with the shear direction. A greater increase was observed for flow in the direction perpendicular to the shear direction, due to the significant flow channels created by the shearing process.From the obtained flow velocity fields, particle transport was predicted using a streamline particle tracking method with the flow velocity fields (vectors) taken from the flow simulations, yielding particle travel times, breakthrough curves, and the Péclet number, Pe. The transport behavior in the fracture is also anisotropic, and advective transport is greater in the direction parallel with the shear direction. The effect of normal stress on the particle transport is significant, and dispersion becomes larger with increasing normal stress.