of fractures. The transmissivity decreased with the increase in confining pressure, while it increased with the increase in shear displacement. In addition, during the shearing process, the change of transmissivity with Reynolds number can be characterized by a fifth-order polynomial function. The Forchheimer coefficient (β) increased by 71.05-153.15% with the applied confining pressure, while it decreased by 42.72-69.52% with the increase in shear displacement. Variations in the equivalent permeability coefficient with the hydraulic gradient can be divided into three stages of being constant, gentle decrease and significant decrease, and corresponded to linear Darcy flow, weak inertia and strong flow nonlinearity, respectively.
Article highlights• An engraving technique was adopted to copy 3D self-affine fracture surfaces. • Hydro tests with perpendicular flow to shear direction were conducted. • Nonlinear flow responses of rough fractures during shearing were revealed.