岩体粗糙单裂隙非达西流动各向异性特性  

Anisotropy in non-Darcy flow in individual rough-walled rock fractures

作　　者：朱寅斌 廖震 李长冬[1,3] 刘洪斌 姜茜慧 ZHU Yinbin;LIAO Zhen;LI Changdong;LIU Hongbin;JIANG Xihui(School of Engineering,China University of Geosciences,Wuhan 430074,China;Yajiang Clean Energy Science and Technology Research(Beijing)Corporation Limited,Beijing 100038,China;Badong National Observation and Research Station of Geohazards,China University of Geosciences,Wuhan 430074,China)

机构地区：[1]中国地质大学(武汉)工程学院,湖北武汉430074 [2]雅江清洁能源科学技术研究(北京)有限公司,北京100038 [3]中国地质大学(武汉)湖北巴东地质灾害国家野外科学观测研究站,湖北武汉430074

出　　处：《煤田地质与勘探》2025年第2期130-146,共17页Coal Geology & Exploration

基　　金：国家自然科学基金重大项目(42090054);湖北省自然科学基金创新群体项目(2022CFA002)。

摘　　要：【目的】研究岩石粗糙单裂隙中发生的流体流动行为是全面了解裂隙网络渗流特性的基础。【方法】基于分形理论生成人工合成岩石裂隙,采用COMSOL Multiphysics软件开展渗流数值模拟研究;利用3D打印技术制备与数值模型具有相同几何特征的裂隙试样,并结合自主搭建的渗流实验装置开展一系列不同注入流量下的渗流试验,研究了不同注入方向对三维粗糙单裂隙非达西流动微观和宏观流动行为的影响。【结果和结论】结果表明:(1)不同注入方向条件下岩石裂隙表面粗糙度的差异会显著影响裂隙的渗流特性,在表面粗糙度更大的方向上流体流动阻力、流动路径迂曲度和非达西效应更大。(2)裂隙中涡旋区发育程度与表面粗糙度呈正相关,在表面粗糙度更大的注入方向上涡旋区的发育程度比表面粗糙度较小的注入方向上涡旋区发育程度更高。(3)Forchheimer方程在描述裂隙中的非达西流动时不能适用于整个雷诺数范围,只有在非达西流动充分发生时才能有足够的准确性。(4)裂隙表面粗糙度各向异性程度越大则裂隙在粗糙度更大的方向上越容易发生非达西流动,两个注入方向上的流动行为差异也越明显。(5)基于三维裂隙渗流模拟研究建立了裂隙临界雷诺数与表面粗糙度各向异性因子以及平均开度的定量化表征模型,并通过开展渗流试验验证了所建立模型的有效性。研究成果可为更全面研究粗糙裂隙中非达西流动行为的各向异性提供借鉴。[Objective]Investigating the flow behavior of fluids in individual rough-walled rock fractures is fundamental for fully understanding the seepage characteristics of fluids in rock fracture networks.[Methods]Rock fractures were artificially synthesized based on the fractal theory,and numerical simulations of seepage along these fractures were performed using the COMSOL Multiphysics software.Then,fracture specimens with the same geometric characteristics as numerical models of fractures were prepared using 3D printing technology.In combination with the self-designed equipment for seepage experiments,this study conducted seepage experiments on these fracture specimens under different injection flow rates to investigate the impacts of varying injection directions on the microscopic and macroscopic behavior of non-Darcy flow in 3D individual rough-walled fractures.[Results and Conclusions]The results indicate that the differences in the surface roughness of rock fractures under different injection directions significantly influenced the seepage characteristics of fractures,with greater flow resistance,flow path tortuosity,and non-Darcy flow effect occurring in the injection direction with higher fracture surface roughness.The developmental degree of eddies in fractures was posit-ively correlated with the fracture surface roughness.Specifically,eddies in the injection direction with higher fracture surface roughness manifested a higher developmental degree compared to those in the direction with lower surface roughness.The Forchheimer equation did not apply to the entire range of Reynolds numbers when describing non-Darcy flow in fractures.Instead,it was sufficiently accurate only in the case of the sufficient occurrence of non-Darcy flow.The higher anisotropy in the surface roughness,non-Darcy flow was more prone to occur along fractures in the direction with higher surface roughness,accompanied by more pronounced differences in flow behavior between both injection directions.Quantitative characterization models for d

关 键 词：岩石单裂隙 各向异性 非达西流动 数值模拟 渗流试验 

分 类 号：P618.11[天文地球—矿床学]