Numerical investigation on the non-Darcy flow in fractures under the joint influence of fracture dilation and inertial effect
-
摘要:
裂隙渗流是核废料、CO2地质储存等重大工程的关键因素,常涉及高压环境。随着水力梯度增加,惯性效应增强,导致非达西效应(表观渗透率减少)。然而,前人研究常忽略高压下水力耦合作用对裂隙的扩张及相应的表观渗透率上升效应,造成刻画和模拟渗流过程不准确。为此,基于二维粗糙单裂隙和裂隙-基质系统,采用直接数值模拟方法研究裂隙扩张和惯性效应2种机制共同作用下的非达西渗流发展规律,探讨不同岩石基质力学性质对裂隙非达西渗流的影响。结果表明:(1)当压力梯度较小时,渗流处于达西流态,2种机制对渗流的影响都可以忽略,随着压力梯度逐渐升高,渗流进入非达西流态,2种机制都起到重要作用;(2)在2种机制竞争作用下,非达西渗流由惯性效应主导转向裂隙扩张主导,与之对应,裂隙等效水力开度呈现出先减小后增大的规律;(3)基质岩石的抗变形能力越弱,裂隙扩张主导效应越显著,并且使得由惯性效应转向裂隙扩张作用主导的临界值也越小,反之惯性效应主导显著。研究成果能为准确研判高压环境下裂隙介质中后达西渗流机制提供科学依据。
Abstract:The fluid flow in fractures is crucial for major engineering projects such as nuclear waste disposal and CO2 geological storage, which often involve with high-pressure environments. Previous studies have revealed that as hydraulic gradients increase, inertial effects become stronger, leading to non-Darcy effects (i.e., reduced equivalent permeability). Unfortunately, past studies often neglected the influence of hydromechanical coupling on modifying fracture morphology. This oversight leads to a failure in capturing the effects of pressure-induced dilation and the consequent increase in equivalent permeability, resulting in inaccuracies in characterizing fluid flow in fractures. To this end, this study, based on two-dimensional rough single fractures and fracture-matrix systems, used direct numerical simulation methods to investigate the joint effects of two mechanisms: fracture dilation and inertial effects on the development of non-Darcy flow, and explored the influence of different mechanical properties of rock matrix on non-Darcy flow in fractures. The study leads to following key findings: (1) when the pressure gradient is small, the fluid flow is in Darcy flow regime, and the effects of both mechanisms on fluid flow can be ignored. As pressure gradient gradually increases, the non-Darcy flow regime emerges, where both mechanisms play important roles. (2) Under the joint effects of above two mechanisms, non-Darcy regime can be further divided into two phases: inertial effect domination and fracture dilation domination. Correspondingly, the equivalent hydraulic aperture decreases first and then increases. (3) Moreover, the weaker the deformation resistance of the matrix rock, the more significant the dominant effect of fracture dilation, and the smaller the critical value at which the transition from inertial effect to fracture dilation occurs. Conversely, inertial effects dominate significantly. The findings of this study can provide a scientific basis for accurately assessing fluid flow in fractured media in high-pressure environments, ultimately helping better design and manage major engineering projects.
-
-
图 1 裂隙中不同流态下的宏观渗流概念图 (修改自Zhou等 [10])
Figure 1.
表 1 基质岩性的力学性质
Table 1. List of the mechanical properties of five typical types of rock matrix
页岩 砂岩 石灰岩 花岗岩 玄武岩 杨氏模量/GPa 15 35 45 60 70 泊松比 0.15 0.22 0.28 0.25 0.25 表 2 不同岩性Forchheimer方程拟合参数
Table 2. Summary of Forchheimer equation fitting parameters for different types of rock matrix
参数 页岩 砂岩 灰岩 花岗岩 玄武岩 无基质 A/(1010 kg·m−5·s−1) 4.88 4.88 4.88 4.88 4.89 4.96 B/(1013 kg·m−8) 6.09 7.59 7.99 8.13 8.23 8.24 R2 0.99 0.99 0.99 0.98 0.96 0.98 -
[1] VISWANATHAN H S,AJO-FRANKLIN J,BIRKHOLZER J T,et al. From fluid flow to coupled processes in fractured rock:Recent advances and new frontiers[J]. Reviews of Geophysics,2022,60(1):e2021RG000744. doi: 10.1029/2021RG000744
[2] 尚慧,柳思航,甘智慧,等. 浅埋煤层群开采覆岩垮落及导水裂隙带发育规律研究[J/OL]. 水文地质工程地质,[2024-03-16]. [SHANG Hui,LIU Sihang,GAN Zhihui,et al. Development of overlying strata collapse and water-conducting fractured zone in shallow coal seams mining[J/OL]. Hydrogeology & Engineering Geology,[2024-03-16]. Https://doi.org/10.16030/j.cnki.issn.1000-3665.202309025. (in chinese with English abstract)]
SHANG Hui, LIU Sihang, GAN Zhihui, et al. Development of overlying strata collapse and water-conducting fractured zone in shallow coal seams mining[J/OL]. Hydrogeology & Engineering Geology, [2024-03-16]. Https://doi.org/10.16030/j.cnki.issn.1000-3665.202309025. (in chinese with English abstract)
[3] 焦婷婷,邓亚平,钱家忠,等. 基于集合平滑数据同化与直接采样法融合水文地球物理数据刻画裂隙网络分布[J]. 水文地质工程地质,2024,51(4):88 − 100. [JIAO Tingting,DENG Yaping,QIAN Jiazhong,et al. Characterizing fracture networks by integrating hydrogeophysical data based on the ESMDA-DS method[J]. Hydrogeology & Engineering Geology,2024,51(4):88 − 100. (in Chinese with English abstract)]
JIAO Tingting, DENG Yaping, QIAN Jiazhong, et al. Characterizing fracture networks by integrating hydrogeophysical data based on the ESMDA-DS method[J]. Hydrogeology & Engineering Geology, 2024, 51(4): 88 − 100. (in Chinese with English abstract)
[4] LEE S H,LEE K K,YEO I W. Assessment of the validity of Stokes and Reynolds equations for fluid flow through a rough-walled fracture with flow imaging[J]. Geophysical Research Letters,2014,41(13):4578 − 4585. doi: 10.1002/2014GL060481
[5] KIM D,YEO I W. Flow visualization of transition from linear to nonlinear flow regimes in rock fractures[J]. Water Resources Research,2022,58(11):e2022WR032088. doi: 10.1029/2022WR032088
[6] ZOU Liangchao,JING Lanru,CVETKOVIC V. Roughness decomposition and nonlinear fluid flow in a single rock fracture[J]. International Journal of Rock Mechanics and Mining Sciences,2015,75:102-118.
[7] CHAUDHARY K,CARDENAS M B,DENG W,et al. The role of eddies inside pores in the transition from Darcy to Forchheimer flows[J]. Geophysical Research Letters,2011,38(24):L24405.
[8] CHEN Yifeng,HU Shaohua,HU Ran,et al. Estimating hydraulic conductivity of fractured rocks from high-pressure packer tests with an Izbash’s law-based empirical model[J]. Water Resources Research,2015,51(4):2096 − 2118. doi: 10.1002/2014WR016458
[9] VOGLER D,SETTGAST R R,ANNAVARAPU C,et al. Experiments and simulations of fully hydro-mechanically coupled response of rough fractures exposed to high-pressure fluid injection[J]. Journal of Geophysical Research (Solid Earth),2018,123(2):1186 − 1200. doi: 10.1002/2017JB015057
[10] ZHOU Jiaqing,CHEN Yinfeng,TANG Huiming,et al. Disentangling the simultaneous effects of inertial losses and fracture dilation on permeability of pressurized fractured rocks[J]. Geophysical Research Letters,2019,46(15):8862 − 8871. doi: 10.1029/2019GL083355
[11] ZHOU Jiaqing,CHEN Yinfeng,WANG Lichun,et al. Universal relationship between viscous and inertial permeability of geologic porous media[J]. Geophysical Research Letters,2019,46(3):1441 − 1448. doi: 10.1029/2018GL081413
[12] XU Zengguang,ZHANG Jinsong,SHAO Lei,et al. Numerical simulation investigation on the seepage characteristics of aperture dilated rough single fracture under two-way hydro-mechanical coupling effect[J]. Gas Science and Engineering,2023,110:204865. doi: 10.1016/j.jgsce.2022.204865
[13] WITHERSPOON P A,WANG J S Y,IWAI K,et al. Validity of cubic law for fluid flow in a deformable rock fracture[J]. Water Resources Research,1980,16(6):1016 − 1024. doi: 10.1029/WR016i006p01016
[14] RONG Guan,YANG Jie,CHENG Long,et al. A forchheimer equation-based flow model for fluid flow through rock fracture during shear[J]. Rock Mechanics and Rock Engineering,2018,51(9):2777 − 2790. doi: 10.1007/s00603-018-1497-y
[15] 朱寅斌,李长冬,周佳庆,等. 粗糙岩石单裂隙非达西流动的试验和数值模拟研究[J]. 岩土工程学报,2023,45(6):1278 − 1284. [ZHU Yinbin,LI Changdong,ZHOU Jiaqing,et al. Experimental and numerical studies on non-Darcian flow in single rough-walled rock fracture[J]. Chinese Journal of Geotechnical Engineering,2023,45(6):1278 − 1284. (in Chinese with English abstract)] doi: 10.11779/CJGE20220307
ZHU Yinbin, LI Changdong, ZHOU Jiaqing, et al. Experimental and numerical studies on non-Darcian flow in single rough-walled rock fracture[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(6): 1278 − 1284. (in Chinese with English abstract) doi: 10.11779/CJGE20220307
[16] WANG Lichun,CARDENAS M B,SLOTTKE D T,et al. Modification of the Local Cubic Law of fracture flow for weak inertia,tortuosity,and roughness[J]. Water Resources Research,2015,51(4):2064 − 2080. doi: 10.1002/2014WR015815
[17] GUDMUNDSSON A. Rock Fractures in Geological Processes[M]. Cambridge:Cambridge University Press,2011.
[18] 马亚楠,杨志兵,熊小锋,等. 粗糙裂隙各向异性对非达西渗流特性的影响[J]. 人民长江,2023,54(1):233 − 239. [MA Yanan,YANG Zhibing,XIONG Xiaofeng,et al. Influence of anisotropy on non-Darcian seepage flow behavior in rough fractures[J]. Yangtze River,2023,54(1):233 − 239. (in Chinese with English abstract)]
MA Yanan, YANG Zhibing, XIONG Xiaofeng, et al. Influence of anisotropy on non-Darcian seepage flow behavior in rough fractures[J]. Yangtze River, 2023, 54(1): 233 − 239. (in Chinese with English abstract)
[19] CHEN Yifeng,ZHOU Jiaqing,HU Shaohua,et al. Evaluation of Forchheimer equation coefficients for non-Darcy flow in deformable rough-walled fractures[J]. Journal of Hydrology,2015,529:993 − 1006. doi: 10.1016/j.jhydrol.2015.09.021
[20] XING Kun,QIAN Jiazhong,MA Haichun,et al. Characterizing the relationship between non-darcy effect and hydraulic aperture in rough single fractures[J]. Water Resources Research,2021,57(9):e2021WR030451. doi: 10.1029/2021WR030451
[21] XING Kun,QIAN Jiazhong,ZHAO Weidong,et al. Experimental and numerical study for the inertial dependence of non-Darcy coefficient in rough single fractures[J]. Journal of Hydrology,2021,603:127148. doi: 10.1016/j.jhydrol.2021.127148
-