Influence of Fracture Roughness and Void Space Morphology on Nonlinear Fluid Flow Through Rock Fractures

Yong Luo,Zhenyu Zhang,Lei Zhang,Kangsheng Xue,Kun Long
DOI: https://doi.org/10.1140/epjp/s13360-022-03499-5
2022-01-01
The European Physical Journal Plus
Abstract:The effects of fracture roughness and geometric morphology of void space between two fracture walls on nonlinear fluid flow through rock fractures were investigated by performing fluid dynamic computation on mated and non-mated rock fractures. The fractal dimension D was used to characterize to the morphology of fracture void space, and it shows a positive correlation with either the root mean square of the height of the fracture void space morphology or the standard deviation of roughness angle. Forchheimer equation describes the nonlinear flow behavior through rock fractures well. Compared to mated rock fractures, the unmatched morphology of fracture void space of non-mated rock fractures increased the flow heterogeneities, producing prominent preferential flow and obvious eddy flow in non-mated rock fractures. This renders the nonlinear coefficient in the Forchheimer equation of non-mated rock fractures is generally greater than that of mated rock fractures of identical fracture aperture and roughness. For mated rock fractures, a power-law relationship was proposed to quantify the nonlinear coefficient b in terms of fracture peak asperity R z , the first derivative of the profile Z 2 and fracture aperture e h , and then, the critical Reynolds number for the onset of nonlinear fluid flow was derived. To further describe the influence of fracture void space morphology on the nonlinear fluid flow through non-mated rock fractures, an extended power-law model was proposed by quantifying b in terms of fracture surface roughness parameters R z , Z 2 , aperture e h and fractal dimension D , and the critical Reynolds number to demark the onset of nonlinear flow was subsequently derived. The predicted critical Reynolds number agrees well with that of fluid dynamic computation for both mated and non-mated rock fractures, validating the proposed power-law and extended power-law relationships. Our research also shows that the critical Reynolds number generally decreased with the increase in fractal dimension.
What problem does this paper attempt to address?