Monitoring the hydraulic properties within subsurface fractures is vitally important in the contexts of geoengineering developments and earthquakes. Geophysical observations are promising tools for remote determination of subsurface hydraulic properties; however, quantitative interpretations are hampered by the paucity of relevant geophysical data for fractured rock masses. This study explored simultaneous changes in hydraulic and geophysical properties of natural rock fractures with increasing normal stress and correlated these property changes through coupling experiments and digital fracture simulations. We show that electrical resistivity is linked with permeability and flow area regardless of fracture roughness, whereas elastic wave velocity is roughness dependent. We also are able to categorize fracture flow patterns as aperture-dependent, aperture-independent, or disconnected flows, with transitions at specific stress levels. Elastic wave velocity offers potential for detecting the transition between aperture-dependent flow and aperture-independent flow, and resistivity is sensitive to detect the connection/disconnection of the fracture flow.

In an effort to reveal the subsurface hydraulic changes in fractures by seismic monitoring, aperture-related velocity changes need to be investigated. We developed a numerical approach for calculating changes in elastic wave velocity with fracture aperture opening by determining the internal energy of a digitized fracture model based on natural rough surfaces. The simulated local elastic energy revealed that the interaction energy converged within 1.5 mm of the mean fracture position, and was insignificant unless the fractures intersected. This energetic approach clarified the aperture–velocity relationship and reproduced the experimental results. Further calculations using digital fractures with various sizes and density demonstrated that the velocity can be accounted for by the superposition of a linear function of fracture density and quadratic function of aperture, and is insensitive to the fracture size. Although the relationship between fracture permeability and elastic wave velocity (i.e., the k-V relationship) depends on the fracture density, the offset-normalized k-V relationship shows clear linearity with the fracture density. The proposed k-V relationship as a function of the aperture and fracture density indicates that laboratory-scale fracture properties of a single fracture can be applied to multiple fractures on a larger scale. Our findings can be used to interpret temporal changes in seismic observations and to monitor fluid flow in fractures.