裂隙网络非达西渗流REV及非达西系数张量研究  

Non-Darcy seepage rev and non-Darcy coefficient tensor in fracture network

作　　者：牛玉龙 王媛[1,3] 于可 冯迪[1,3] NIU Yulong;WANG Yuan;YU Ke;FENG Di(Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering,Hohai University,Nanjing 210098,China;China Three Gorges Corporation,Beijing 100038,China;College of Civil and Transportation Engineering,Hohai University,Nanjing 210098,China)

机构地区：[1]河海大学岩土力学与堤坝工程教育部重点实验室,江苏南京210098 [2]中国长江三峡集团公司,北京100038 [3]河海大学隧道与城市轨道工程研究所,江苏南京210098

出　　处：《水利学报》2020年第4期468-478,共11页Journal of Hydraulic Engineering

基　　金：国家自然科学基金项目(U1765204,41772340);浙江省水利厅科技项目(RA1503)。

摘　　要：裂隙网络渗流理论最初建立在裂隙渗流服从立方定律的基础上,而裂隙中非达西渗流更为广泛存在,因此,裂隙网络非达西渗流特性研究更具工程应用前景。本文首先对基于单裂隙立方定律的线素法模型进行改进,用以实现裂隙网络非达西渗流的计算,进而将其应用于实测裂隙数据生成的裂隙网络非达西渗流特性的研究。研究结果表明:根据改进线素法模型和Forchheimer定律计算得到的等效渗透系数张量与线素法模型和达西定律计算得到的具有一致性。对于同一裂隙网络而言,等效非达西系数张量确定的渗流REV要比等效渗透系数张量确定的大,且等效渗透系数张量与等效非达西系数张量的主轴方向一致,但两者之间的长轴互相垂直。The theory of fracture network seepage is initially established on the basis of the cubic law obeyed by fracture seepage,while non-Darcy seepage in fractures exists more widely.Therefore,the research on non-Darcy seepage characteristics of fracture network has more engineering application prospect.In this paper,the linear element method based on the cubic law of single fracture is improved to calculate non-Darcy seepage in fracture network.Then it is applied to the study of non-Darcy seepage characteristics of fracture network generated by measured fracture data.The results show that the equivalent permeability coefficient tensor calculated by improved linear element method model and Forchheimer's law is consistent with those calculated by linear element method model and Darcy's law.For the same fracture network,the seepage REV determined by the equivalent non-Darcy coefficient tensor is larger than that by the equivalent permeability coefficient tensor.The principal axes of the equivalent permeability coefficient tensor and the equivalent non-Darcy coefficient tensor are the same,but the long axes between them are perpendicular to each other.

关 键 词：裂隙网络 线素法模型 非达西渗流 表征单元体积 非达西系数张量 

分 类 号：TV223.6[水利工程—水工结构工程]