裂隙岩体网络注浆渗流规律.docVIP

裂隙岩体网络注浆渗流规律 ) 摘 要：本文在裂隙结构模拟的基础上，建立了裂隙网络非牛顿流体的渗流模型，并应用所建立的分析模型进行注浆渗透过程模拟，建立了确定注浆参数的一般方法。在此基础上，本文还讨论了裂隙岩体注浆过程中的渗流和应力耦合问题，分析了应力和渗流耦合对注浆过程的影响。 关键词：裂隙模拟；非牛顿流体；渗流；应力；注浆参数 收稿日期：2000-03-09作者简介：杨米加(1971-)，男，安徽和县人，中国矿业大学博士，现为河海大学渗流室在站博士后，主要从事裂隙渗流方面的研究。 　　天然岩体内分布着大量的裂隙和孔隙。一般地，岩体中的孔隙尽寸较小但数量较多，而裂隙数量少但其渗透性好，因此在岩体渗流中对渗流起控制作用的主要是裂隙。特别是对于注浆工程，由于注浆浆液都为比较粘稠的悬浮液，在孔隙中几乎是不可渗，因此我们可假定在岩体注浆中，浆液只在岩体裂隙网络中渗流。 1 裂隙岩体网络模型 　　对岩体裂隙网络模型的研究是从研究地下水活动规律开始的[1～3]。早期的研究一般认为岩体中的裂隙规则分布，或把岩体作为一连续体，把裂隙均匀到岩体的各个质点上去作多孔介质处理。而实际上，裂隙是在长期地质历史作用下形成的，而且由于地质组成的不确定性以及多次地应力的作用，裂隙在岩体中的分布，既体现了一定的规律性又表现出一定的随机性。为了形成即便于计算又反映裂隙赋存状态和裂隙渗流规律的离散裂隙网络模型，采用蒙特卡洛法模拟生成岩体裂隙网络是最合适不过的。蒙特卡洛模拟(或称蒙特卡洛法)是一种具有独特风格的数值计算方法，可求解具有随机性的不确定性问题。利用蒙特卡洛法求解问题时，首先要建立一个概率模型，把需要求解的问题同概率模型结合起来，然后通过随机抽样试验得到某些所求问题解的概率分析，这就是蒙特卡洛法的基本思想。裂隙网络结构模型的蒙特卡洛法模拟就是根据现场裂隙的统计参数及规律，应用蒙特卡洛法产生裂隙参数的随机抽样，并通过计算机图形功能，实现裂隙结构的再现[4]。 2 岩体裂隙网络注浆渗流规律 2.1 岩体裂隙网络注浆 对于注浆浆液这样有较大粘度的液体，它只能在岩体裂隙中流动[2、5]。对于裂隙发育的岩体可采用等效连续介质模型来模拟它在岩体中的流动，但对于裂隙不太发育的岩体，其单个裂隙的几何参数及位置将对浆液流动过程产生重要影响，若应用连续介质模型，可能会造成较大的误差，因此有必要建立岩体裂隙网络注浆渗流模型。 　　岩体裂隙网络是指岩体内不同成因、不同力学性质、不同规模和不同方向的裂隙个体在空间上相互交叉构成的网络状系统，称为裂隙网络系统。裂隙网络系统因受个体裂隙非均质特性的影响，其整体性质可表现为均质性、非均质性、各向异性或各向同性等多种形式，从而也使得裂隙网络岩体的性质更趋于复杂[5]。 2.2 裂隙岩体注浆渗流控制方程 裂隙网络渗流模型是根据裂隙交叉的几何特征，把裂隙交叉点作为节点，节点与节点之间的裂隙作为线单元，并根据各线单元流向共同节点的流量等于零(稳定流)或等于贮存量的变化(非稳定流)，建立裂隙岩体网络的渗流方程，然后再结合初始条件和边界条件就构成了裂隙网络渗流模型(见图1). 图1 应用蒙特卡洛模拟生成的裂隙网络 根据各线单元流向共同节点的流量等于零(稳定流)或等于贮量的变化(非稳定流)，有 (1) 式中：Hfi为i节点上的压力；,Si为裂隙以i为中心的表征体积单元的弹性贮存系数，wj为i点的垂直入渗量，bj、lj为裂隙段的宽度和长度，j为连接i的各裂隙段；N′为特征单元i点的线单元数。 考虑到注浆过程中无垂直补给，且注浆浆液具有非牛顿流体特性，大多属于宾汉姆体，有 q=b=γb3/12ηJ+τ30/3ηbγ2J2-τ0b/4η (2) 结合式(1)有 ()i+Qi=-didHfi/dt(i=1,2,…N) (3) 式中：q为裂隙单元单宽流量；为流速；η为流体动力粘度；τ0为流体的剪切强度；βj=γb3/12η+τ30/3ηbγ2J3-τ0b/4ηJ;Jj为i点表征单元第j个裂隙的渗透梯度。 则式(3)即为裂隙岩体注浆渗流控制方程。 2.3 非牛顿体非稳定渗流场分析及注浆参数的确定 根据上述的结论，作者自己编制了程序，对裂隙岩体的渗流厨行了分析，并比较了牛顿流体和非牛顿流体渗流过程的差异。计算选取模型如图2所示。计算中选取粘度η=0.034Pa·s，注入点压头为5m，渗出点压头为0m，流体为牛顿流体，服从达西渗流定律。第一种方案，以19，20点作为注入点，边界点1，2，3，4，5，10，14，23，28，32，37，38，39，40，41作为渗出点，则计算所得的结果见表1.第二种方案以21点作为注入点，注入点压头为5m，以边界点1，2，3，4，5，10，14，19，23，28，32，37，38，39，40，41为渗出点，渗出点压头为0m，