Vector operator-based description of evolutional velocity field for water flow between rough fractures in granite under large compression-shear displacements
LIU Yintong1，CHENG Jianchao1，ZHANG Liao1，HOU Mengdong1，MAO Tingting1，TIAN Jinzhu1，GAO Weiqi1，LIU Shenggui1，XUE Dongjie1，2
(1. School of Mechanics and Architectural Engineering，China University of Mining and Technology(Beijing)，Beijing 100083，China；2. State Key Laboratory of Water Resource Protection and Utilization in Coal Mining，National Institute of Clean-and-Low-Carbon Energy，Beijing 102211，China)
Abstract：An accurate field description of fracture flow still presents a huge obstacle in flow mechanics. Unlike some classic flow models，i.e. parallel plate flow，focusing on the velocity distribution along the longitudinal section of the fracture，very little attention is paid to the velocity distribution within the fracture cross-section. The spatial distribution of fracture aperture indicates a geometric field，which significantly impacts the flow velocity distribution within the cross-section，while lacking proper correlation models induces a difficult understanding of statistical analysis of geometric field?s effect on fracture flow distribution. To solve the challenge of carrying out a seepage test with relative compression or shear displacement exceeding 5% in lab，the Reynolds equation-based discrete solution of fracture flow considering the large shear-compression displacements for Beishan granite is made，and through the gradient of seepage rate and velocity divergence，the differential characterization of the seepage field's field function was achieved. For the same group of cracks，three fields of flow velocity，gradient and divergence presents a consistent spatial geometric distribution under different shear or compression displacements and nearly all fields conform to the normal distribution. In addition，with the increase in closed aperture and contact area，the full-scale fracture flow turns into local-scale flow and two vector operators of gradient and divergence describing the evolutional distribution of fracture flow are still effective. The results show that in the statistical sense，under low-permeability conditions，a rough fracture flow，characterized by the average gradient and the average divergence of the flow velocity field tend towards 0，can be equivalently treated as uniform flat plate flow with low velocity.
刘殷彤1，程建超1，张 辽1，侯孟冬1，毛婷婷1，田金柱1，高伟琪1，刘升贵1，薛东杰1，2. 压剪大位移下花岗岩粗糙裂隙渗流场演化矢量算子描述[J]. 岩石力学与工程学报, 2025, 44(3): 651-677.
LIU Yintong1，CHENG Jianchao1，ZHANG Liao1，HOU Mengdong1，MAO Tingting1，TIAN Jinzhu1，GAO Weiqi1，LIU Shenggui1，XUE Dongjie1，2. Vector operator-based description of evolutional velocity field for water flow between rough fractures in granite under large compression-shear displacements. , 2025, 44(3): 651-677.
[1] TALON L，AURADOU H，HANSEN A. Permeability estimates of self-affine fracture faults based on generalization of the bottleneck concept[J]. Water Resources Research，2010，46(7)：1–5.
[2] JIN Y，DONG J B，ZHANG X Y，et al. Scale and size effects on fluid flow through self-affine rough fractures[J]. International Journal of Heat and Mass Transfer，2017，105(2)：443–451.
[3] WANG Z，XU C，DOWD P. A modified cubic Law for single-phase saturated laminar flow in rough rock fractures[J]. International Journal of Rock Mechanics and Mining Sciences，2018，103：107–115.
[4] PERATTA A，POPOV V. A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media[J]. Advances in Water Resources，2006，29(1)：42–61.
[5] 王恩志，张 东，刘晓丽，等. 裂隙岩体多结构多流态渗流模型与模拟[J]. 地球科学与环境学报，2022，44(6)：894–902.(WANG Enzhi，ZHANG Dong，LIU Xiaoli，et al. Model and simulation of multi structure and multi fluid flow in fractured rock mass[J]. Journal of Earch Sciences and Environment，2022，44(6)：894–902.(in Chinese))
[6] YU P，LI Y M，ZHAO J Z. A novel approach to simulate the stress and displacement fields induced by hydraulic fractures under arbitrarily distributed inner pressure[J]. Journal of Natural Gas Science and Engineering，2016，35(1)：1 079–1 087.
[7] SUN L，TANG X H，ABOAYANAH K R，et al. Coupled hydro-mechanical two-phase flow model in fractured porous medium with the combined finite-discrete element method[J]. Engineering with Computers，2024，40：2 513–2 535.
[8] HU X D，DONG E J，ZHOU F J，et al. A new numerical model for simulation of flow on rough fracture[J]. Geoenergy Science and Engineering，2024，234：212 540.
[9] BISDOM K，BERTOTTI G，NICK H M. The impact of different aperture distribution models and critical stress criteria on equivalent permeability in fractured rocks[J]. Journal of Geophysical Research：Solid Earth，2016，121(5)：4 045–4 063.
[10] WANG S D，ZHANG Q，ZHAO L，et al. Seepage characteristics study of single rough fracture based on numerical simulation[J]. Applied Sciences，2022，12(14)：7 328–7 345.
[11] SHI Z L，YAO Q L，WANG W N，et al. Size effects of rough fracture seepage in rocks of different scales[J]. Water，2023，15(10)：1 912– 1 934.
[12] ZHANG B，WANG L，LIU J. Effect of fracture geometry parameters on the permeability of a random three-dimensional fracture network[J]. Processes，2023，11(8)：2 237–2 261.
[13] ZIMMERMAN R W，BODVARSSON G S. Hydraulic conductivity of rock fractures[J]. Transport in Porous Media，1996，23(1)：1–30.
[14] 常 远，张晓虎，何满潮. 不饱和十字交叉型裂隙流动的尺寸效应[J]. 矿业科学学报，2018，3(4)：342–348.(CHANG Yuan，ZHANG Xiaohu，HE Manchao. Size effect of unsaturated cross shaped fracture flow[J]. Journal of Mining Science and Technology，2018，3(4)：342–348.(in Chinese))
[15] XIONG F，WEI W，XU C S，et al. Experimental and numerical investigation on nonlinear flow behaviour through three dimensional fracture intersections and fracture networks[J]. Computers and Geotechnics，2020，121(3)：129–145.
[16] XU D，LIU J F，LIANG C，et al. Fracture surface morphology effect on radial seepage flow in a horizontal single granite fracture[J]. Bulletin of Engineering Geology and the Environment，2024，83(1)：55–68.
[17] WANG M，CHEN Y F，MA G W，et al. Influence of surface roughness on nonlinear flow behaviors in 3D self-affine rough fractures：Lattice Boltzmann simulations[J]. Advances in Water Resources，2016，96：373–388.
[18] ZHANG Q，LUO S H，MA H C，et al. Simulation on the water flow affected by the shape and density of roughness elements in a single rough fracture[J]. Journal of Hydrology，2019，573(8/12)：456–468.
[19] BABADAGLI T，REN X J，DEVELI K. Effects of fractal surface roughness and lithology on single and multiphase flow in a single fracture：An experimental investigation[J]. International Journal of Multiphase Flow，2015，68：40–58.
[20] 甘 磊，刘 玉，张宗亮. 岩体裂隙粗糙度表征及其对裂隙渗流特性的影响[J]. 岩土力学，2023，44(6)：1 585–1 592.(GAN Lei，LIU Yu，ZHANG Zongliang. Characterization of rock fracture roughness and its influence on fracture seepage characteristics[J]. Rock and Soil Mechanics，2023，44(6)：1 585–1 592.(in Chinese))
[21] 钟 振，孟 醒，胡云进，等. 考虑基质渗透性的岩体单裂隙渗流及影响因素的室内和数值试验研究[J]. 岩石力学与工程学报，2023，42(9)：2 148–2 163.(ZHONG Zhen，MENG Xing，HU Yunjin，et al. Interior and numerical experimental research on single fracture seepage and influencing factors in rock mass considering matrix permeability[J]. Chinese Journal of Rock Mechanics and Engineering，2023，42(9)：2 148–2 163.(in Chinese))
[22] FINENKO M，KONIETZKY H. Numerical simulation of turbulent fluid flow in rough rock fracture：2D Case[J]. Rock Mechanics and Rock Engineering，2024，57(1)：451–479.
[23] LIU R C，LI B，JIANG Y J. Critical hydraulic gradient for nonlinear flow through rock fracture networks：The roles of aperture，surface roughness，and number of intersections[J]. Advances in Water Resources，2016，88：53–65.
[24] 刘才华，陈从新，付少兰. 二维应力作用下岩石单裂隙渗流规律的实验研究[J]. 岩石力学与工程学报，2002，21(8)：1 194–1 198.(LIU Caihua，CHEN Congxin，FU Shaolan. Experimental study on the seepage law of a single rock fracture under two-dimensional stress[J]. Chinese Journal of Rock Mechanics and Engineering，2002，21(8)：1 194–1 198.(in Chinese))
[25] 李 博，黄嘉伦，钟 振，等. 三维交叉裂隙渗流传质特性数值模拟[J]. 岩土力学，2019，40(9)：3 670–3 678.(LI Bo，HUANG Jialun，ZHONG Zhen，et al. Numerical simulation of mass transfer characteristics in three-dimensional cross fracture seepage[J]. Rock and Soil Mechanics，2019，40(9)：3 670–3 678.(in Chinese))
[26] KOU M M，LIU X R，WANG Z Q，et al. Laboratory investigations on failure，energy and permeability evolution of fissured rock-like materials under seepage pressures[J]. Engineering Fracture Mechanics，2021，247(1/4)：107 694.
[27] 薛东杰，周宏伟，任伟光，等. 北山花岗岩深部节理间距分布多重分形研究[J]. 岩土力学，2016，37(10)：2 937–2 944.(XUE Dongjie，ZHOU Hongwei，REN Weiguang，et al. Multi-fractal characteristics of joint geometric distribution of granite in Beishan[J]. Rock and Soil Mechanics，2016，37(10)：2 937–2 944.(in Chinese))
[28] BERKOWITZ B. Characterizing flow and transport in fractured geological media：A review[J]. Water Resource，2002，25(8/12)：861–884.
[29] TAUCARE M，VIGUIER B，DANIELE L，et al. Connectivity of fractures and groundwater flows analyses into the western andean front by means of a topological approach(Aconcagua Basin，Central Chile)[J]. Hydrogeol，2020，28(7)：2 429–2 438.
[30] LIU X Y，ZHANG X W，KONG L W，et al. Multiscale structural characterizations of anisotropic natural granite residual soil[J]. Canadian Geotechnical Journal，2023，60(9)：1 383–1 400.
[31] LIU X Y，ZHANG X W，KONG L W，et al. Inherent and stress-induced stiffness anisotropy of natural granite residual soil[J]. Acta Geotech，2023，18(11)：5 681–5 699.
[32] WEN J，TANG L，ZHANG S，et al. Qualitative and quantitative investigations on the failure effect of critical fissures in rock specimens under plane strain compression[J]. Materials，2023，16(2)：611–628.
[33] JONES T A，DETWILER R L. Fracture sealing by mineral precipitation：The role of small-scale mineral heterogeneity[J]. Geophysical Research Letters，2016，43(14)：7 564–7 571.
[34] YU X，VAYSSADE B. Joint profiles and their roughness parameters[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts，1991，28(4)：333–336.
[35] BARBATI A C，DESROCHES J，ROBISSON A，et al. Complex fluids and hydraulic fracturing[J]. Annual Review of Chemical and Biomolecular Engineering，2016，7(1)：415–453.
[36] IM K，ELSWORTH D，FANG Y. The influence of preslip sealing on the permeability evolution of fractures and faults[J]. Geophysical Research Letters，2018，45(1)：166–175.
[37] MADADI M，SAHIMI M. Lattice Boltzmann simulation of fluid flow in fracture networks with rough，self-affine surfaces[J]. Physical Review E，2003，67(2)：026309.
[38] BRUSH D J，THOMSON N R. Fluid flow in synthetic rough-walled fractures：Navier-Stokes，Stokes，and local cubic law simulations[J]. Water Resources Research，2003，39(4)：1 037–1 041.
[39] CHAKRABORTY G. A note on methods for analysis of flow through microchannels[J]. International Journal of Heat and Mass Transfer，2008，51(17/18)：4 583–4 588.
[40] WATANABE N，HIRANO N，TSUCHIYA N. Determination of aperture structure and fluid flow in a rock fracture by high-resolution numerical modeling on the basis of a flow-through experiment under confining pressure[J]. Water Resources Research，2008，44(6)：W06 412.
[41] CRANDALL D，BROMHAL G，KARPYN Z T. Numerical simulations examining the relationship between wall-roughness and fluid flow in rock fractures[J]. International Journal of Rock Mechanics and Mining Sciences，2010，47(5)：784–796.
[42] XUE D J，LIU Y T，ZHOU H W，et al. Fractal characterization on anisotropy and fractal reconstruction of rough surface of granite under orthogonal shear[J]. Rock Mechanics and Rock Engineering，2020，53(7)：1 225–1 242.
[43] GONG H L，LUO Y，ZHOU J R，et al. Fracture behaviors and damage evolution anisotropy of granite under coupling of multiaxial confinement and dynamic loading[J]. Rock Mechanics and Rock Engineering，2023，56(4)：2 515–2 534.
[44] DU H B，DAI F，WEI M D，et al. Dynamic compression-shear response and failure criterion of rocks with hydrostatic confining pressure：an experimental investigation[J]. Rock Mechanics and Rock Engineering，2021，54(2/3)：955–971.
[45] 胡盛斌，杜国平，徐国元，等. 基于能量测量的声呐渗流矢量法及其应用[J]. 岩土力学，2020，41(6)：2 143–2 154.(HU Shengbin，DU Guoping，XU Guoyuan，et al. Sonar seepage vector method based on energy measurement and its application[J]. Rock and Soil Mechanics，2020，41(6)：2 143–2 154.(in Chinese))
[46] XUE D J，GAO L，LU L，et al. An acoustic emission-based cluster damage model for simulating triaxial compression behaviors of granite[J]. Rock Mechanics and Rock Engineering，2020，53(4)：4 201–4 220.
[47] HE M，YU L Y，LIU R C，et al. Experimental investigation on mechanical behaviors of granites after high-temperature exposure[J]. Journal of Central South University，2022，29(4)：1 332–1 344.
[48] LI Y Z，YUAN L，ZHANG Q H. Brittle rock mass failure in deep tunnels：The role of infilled structural plane with varying dip angles[J]. International Journal of Rock Mechanics and Mining Sciences，2024，176：105 721.
[49] XUE D J，ZHOU H W，ZHAO Y，et al. Real-time SEM observation of mesoscale failures under thermal-mechanical coupling sequences in granite[J]. International Journal of Rock Mechanics and Mining Sciences，2018，112：35–46.
[50] LEE H P，SCHULTZ R A，OLSON J E. Inferring friction and stress states from the angle between wing fractures and solution surfaces[J]. Journal of Structural Geology，2022，162(2)：104 706.
[51] HUANG J Z，TANG H Z，PAN X，et al. Investigation of the compression-shear fracture propagation for rocks accounting for confining pressure and crack surface friction[J]. Theoretical and Applied Fracture Mechanics，2022，119(3)：103346.
[52] 窦子豪. 岩石裂隙剪切力学行为的温度与水岩作用效应[博士学位论文][D]. 北京：清华大学，2022.(DOU Zihao. Effects of temperature and water rock interaction on shear behaviors of rock fractures[Ph. D. Thesis][D]. Beijing：Tsinghua University，2022.(in Chinese))
[53] BROWN S R，STOCKMAN H W，REEVES S J. Applicability of the Reynolds Equation for modeling fluid flow between rough surfaces[J]. Geophysical Research Letters，1995，22(18)：2 537–2 540.
[54] STEPANIUK V P，TARAU C，ÖTÜGEN M V，et al. Evaluation of a laser velocity gradient probe and measurements in a boundary layer[J]. Experiments in Fluids，2004，36(3)：510–514.
[55] KUZNETSOV E A，MIKHAILOV E A. Slipping flows and their breaking[J]. Annals of Physics，2022，447(2)：169 088.
2025, Vol. 44 
