单裂隙高速非达西渗流特征及参数研究.docx

单裂隙高速非达西渗流特征及参数研究 摘 要 随着我国水利、矿山、石油、核废料储存等事业的发展，地下空间的开发与 利用过程中常常遇到裂隙岩体这类复杂介质，其渗透特性研究成为当前岩体力学 研究领域的热点之一，在传统的研究中，仍然采用基于达西定律的方法。然而理 论研究及实践工程表明，裂隙中的渗流往往不再遵循线性的达西定律，故本文开 展了高速非达西渗流的研究，并基于纳维斯托克斯方程采用有限体积法对单裂隙 中的水流进行数值模拟，探讨了单裂隙中非达西渗流的非线性参数及特征，主要 工作如下： 1）系统总结了非达西渗流的研究成果，探讨了非线性机理和达西渗流与非达 西渗流的判别准则，将描述高速非达西渗流的两类运动方程进行对比，结果表明 Forchheimer 方程优于指数方程；利用有限元软件实现了二维非达西渗流的模拟， 通过对比分析达西渗流与非达西渗流，得出两者的计算结果随着水力梯度的增高 相差越来越大，揭示了非达西渗流研究的必要性； 2）建立了与室内实验对应的光滑裂隙模型，利用流体数值分析软件模拟了光 滑裂隙中的高速非达西渗流，通过数值模拟与室内实验的对比，证实了 Fluent 软 件模拟单裂隙中高速非达西流的可靠性；在数值模拟结果的基础上，提出用量纲 分析的方法来研究非达西影响系数，得到了光滑裂隙中的非达西影响系数与隙宽 之间的关系式； 3）根据粗糙裂隙渗流的影响因素建立了粗糙裂隙模型，从而开展了粗糙裂隙 中的高速非达西渗流研究，结果表明可以用 Forchheimer 方程来描述此间的水流运 动；对比分析不同粗糙裂隙模型，得到粗糙裂隙的非线性程度取决于粗糙度的形 状和凸起度，与间隙比无关； 4）揭示了单裂隙中非线性渗流的根本原因是有效渗透率的降低，围绕非达西 效应研究了单裂隙中非达西渗流的非线性特征，提出了用非达西效应来作为非线 性的临界值。 关键词：裂隙渗流、高速非达西渗流、有效渗透率、非达西效应、非达西影响系 数、Forchheimer 数。 I 目录 目录 目 录 摘 要 I ABSTRACT III 目 录 V 第一章 绪论1 1.1 课题研究背景和意义 1 1.2 研究历史及现状 1 1.2.1 理论及室内实验研究 2 1.2.2 数值方法研究8 1.3 本文研究内容及方法 12 第二章 非达西渗流基本特征 15 2.1 达西定律适用范围 15 2.2 高速非达西渗流产生机理 16 2.3 描述高速非达西流的方程 17 2.3.1 平方多项式形式的非达西流运动方程 18 2.3.2 指数形式的非达西流运动方程20 2.4 非达西渗流的判别准则21 2.4.1 临界值 21 2.4.2 非达西效应 22 2.4.3 非达西影响系数24 2.5 非达西流数值模拟 26 2.5.1Comsol Multiphysics 软件 26 2.5.2 达西与非达西渗流对比 26 2.6 本章小结28 第三章 裂隙水流及数值模拟基本理论29 3.1 裂隙岩体地下水流动模型 29 3.2 单裂隙水流 30 3.2.1 基于达西定律的立方定律 30 3.2.2 非达西渗流基本方程 31 3.2.3 定解条件 33 3.3Fluent 软件模拟介绍 33 3.3.1 软件简介 34 3.3.2 流动基本方程34 3.3.3 求解步骤 35 3.3.4 求解算法 36 3.3.5 流动模型 37 3.4 本章小结38 第四章 光滑单裂隙高速非达西流数值模拟 41 V 河海大学硕士学位论文 河海大学硕士学位论文 VI VI 4.1 光滑单裂隙水流数值模拟 41 4.1.1 数值模型的建立与模拟 41 4.1.2 数值模拟结果43 4.2 单裂隙非达西渗流特征参数分析 48 4.2.1 水流运动方程48 4.2.2 非达西影响系数50 4.2.3 有效渗透率 56 4.2.4 非线性特征值探讨 57 4.3 本章小结60 第五章 粗糙单裂隙高速非达西流数值模拟 63 5.1 粗糙裂隙模型建立 63 5.2 数值模拟结果 64 5.2.1 固有渗透率 64 5.2.2Forchheimer 拟合 66 5.2.3 非达西影响系数68 5.3 数值模拟结果分析 69 5.3.1 有效渗透率 69 5.3.2 粗糙因素对非线性渗流的影响70 5.4