三维离散裂隙网络渗流模型与实验模拟

《三维离散裂隙网络渗流模型与实验模拟》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、三维离散裂隙网络渗流模型与实验模拟

2、第1内容加载中...　摘要：本文根据天然裂隙系统发育规律及其渗透机制，将复杂的裂隙系统划分成带状断层、面状裂缝和管状孔洞三大类型，在忽略岩块渗透性的前提下，建立了由管状线单元、缝状面单元和带状体单元组合而成的三维裂隙网络渗流数值模型，为检验模型的合理性，在室内利用混凝土试块制成的离散体来模拟三维裂隙网络结构，设定不同的边界水位和降雨入渗对裂隙中渗流进行了实验，利用数值模型对实验结果进行了模拟，初步验证模型是合理而可靠的。关键词：裂隙网络；渗流；实验；数值模拟基金项目：国家自然科学基金资助项目（59479010；40172084）自

3、Wittke提出线素模型[1、2]以来，二维裂隙网络渗流模型有所发展[3、4]，进入20世纪80年代后，Long提出了圆盘裂隙网络三维渗流模型[5]，Dershowitz、万力等提出了多边形裂隙网络模型[6，7]，这些模型可描述面状裂缝相互切割所构成的裂隙网络渗流。然而，实际岩体中的裂隙系统并非都是由单一的面状裂缝所构成，为此，根据天然裂隙系统发育规律及其渗透机制，将复杂的裂隙系统划分成三大类型，在忽略岩块渗透性的前提下，建立了由管状线单元、缝状面单元、带状体单元组合而成的三维裂隙网络渗流数值模型[8]。为检验模型的合理性，在室内利用混凝土试块制成的离散体来模拟三维

4、裂隙网络结构，设定不同的边界水位和降雨入渗对裂缝中渗流进行了实验。本文就三维裂隙网络渗流实验结果和理论模型的数值模拟对比进行了分析，其结果初步表明理论模型是合理可行的。1三维裂隙网络渗流数值模型复杂岩体中的裂隙网络可定义为由不同成因、不同力学属性、不同规模的管状、面状和带状结构体相互交错切割所构成的网络状系统[8]，由管状结构的一维流、面状结构的二维流和带状结构的三维流在空间上组合，形成了岩体中复杂的三维裂隙网络渗流。因此，可按三种类型结构单元来建立三维裂隙网络渗流数值模型。图1裂隙单元示意1.1管状结构一维渗流方程管状结构被其它裂隙切割后为一条线段，以该线段作为线

5、单元，则一条管状结构可划分成数条线单元的组合。沿裂隙段方向建立局部的一维坐标系(图1(a)).设裂隙段中的水流为层流，线单元的渗流连续方程为：(1)1.2面状裂隙二维渗流方程面状裂缝类结构为典型的裂隙单元，被其它裂隙切割后，成为数个多边形裂缝（图1(b)），每个多边形裂缝在空间上表现为多边形面单元，按该单元所在的平面建立局部坐标系x′o′y′，并取局部坐标轴方向与裂缝的渗透主轴方向一致，设裂隙中水流为层流时，则局部坐标系下的二维渗流连续方程为：(2)1.3带状结构三维渗流方程带状断层类结构可分为构造岩带和影响带，被其它裂隙切割后，每个带就成为空间中的数个扁多面体（图

6、1（c）），其内结构为碎裂岩块的集合体（挤压面属面状结构），从渗透介质的特征上属于多孔介质。则根据带状断层的结构特点，对每一个扁多面体可用连续介质方法来描述其中的渗流。因此，当选取的局部坐标系与带状断层的3个渗透主轴一致时，其渗流连续性方程为：(3)式中：x′、y′、z′为局部三维坐标系；Kx′、Ky′、Kz′为3个主渗方向上的渗透系数；H为水头函数；μ为单位贮水率；Q为源汇项；t为时间变量。1.4裂隙网络三维渗流数值方程将上述3种类型结构单元体进行离散化，对断层单元的任一扁多面体划分成数个四面体、五面体或六面体的组合，对裂缝单元的任一多边形划分成数个三角形或四边形

7、的组合，对管状单元又可划分为数条线单元的组合。这样根据线单元、面单元、体单元在空间的组构关系，对各类单元统一编号，由于各种单元相交存在公共节点，则节点也实行统一编号，由上面的一维、二维和三维渗流方程，分别求得各类单元上的传导矩阵，然后，按单元的相互关系，对各节点上的单元传导矩阵进行迭加组合，形成总体传导矩阵，得到由多种裂隙单元构成的三维裂隙网络渗流数值模型：{T}{H}+{Q}={E}{dH/dt}(4)式中：{H}为节点水头列向量；{T}为总体传导矩阵；{E}为裂隙贮水矩阵；{Ｑ}为补给和排泄列向量。式（4）反映的介质结构为线单元、面单元和体单元的组合，其渗流是由

8、一维流、二维流和三维流在空间上的组合，是由局部的连续介质渗流组合成的整体上的非连续介质渗流。断层类结构按其性质将有不同的划分，张性断层、劈理带、岩脉等可划分成一个条带，压性断层可划分数个条带，中部挤压面则作为面状单元嵌入在体单元中。