5.第二章 渗流理论.doc

第二章 地下水在岩土介质中的渗流 浙江大学硕士学位论文 2006年8月 第二章 地下水在岩土介质中的渗流 杭州地区粉砂土中基坑降水面的数值模拟 -  PAGE 26 - - PAGE 25- 第二章 地下水在岩土介质中的渗流 水在岩土等孔隙或裂隙介质中的流动就称为渗流，其流动性质取决于作为渗流骨架的岩土性质以及其中流动的水的性质。 2.1 地下水渗流模拟技术发展概况 19世纪中叶至2O世纪初，地下水流模拟研究刚刚起步，经历了从单层稳定流模型发展到多层越流非稳定流模型，但直到2O世纪5O年代末，对含水层内部结构和不规则边界条件的研究仍无明显进展。60年代以来，数值计算方法在水文地质学中的应用及电子计算机技术的广泛使用，使复杂含水层系统中地下水流动及溶质运移的模拟成为可能，此阶段是地下水流数值模型迅速发展的时期，先后出现了二维流平面(削面)模型、准三维流模型、三维流模型和耦合模型等。国外90年代做的几个模型，例如R．Bravo(美国)等做的休斯敦模型，A．RJvera(法国)等做的墨西哥城模型，G．Gambolati(意大利)做的拉温纳区域地下水流模型，都属于三维流模型或准三维模型”。我国地下水数值模型研究起步较晚，开始于2O世纪7O年代，但经过数学工作者和水文地质工作者以及科研院所的共同努力，现已取得了长足发展。 渗流计算是在已知定解条件下解渗流基本方程，以求得渗流场水头分布和渗流量等渗流要素，是工程设计的重要内容，本论文主要关注的是渗流场内水头的分布，也就是地下水自由水面的分布情况。由于无压渗流有渗流自由面（浸湿线），且非稳定渗流自由面随基坑地下水位的降低而变动，加之一般渗流场有不同程度的非均质性和各向异性，几何形状和边界条件较复杂，解析求解在数学上带来了不少困难，仅能对一些简单的流动情况获得解析解。虽然有人提出了很多近似解，但由于实际情况的特殊性，应用上仍受一定条件的限制，实际工程往往借助模拟试验求解。随着计算机的普及和数值计算方法的发展，特别是有限单元法提出后，推进了渗流数学模型的发展，为渗流计算提供了有效的方法。 2.2 渗流基本方程和定解条件 2.2.1 渗流运动方程 地下水运动方程可由作用到液体上各力的平衡求得，这些力包括有液体表面的水压力、重力、渗流受到的阻力、加速力，这些力可概括为表面力和体积力两个类别。推导过程与一般流体力学中的运动方程相同，可以直接引用一般流体运动方程，只要把水质点运动速度当作多孔介质中孔隙水流运动速度，再按照孔隙水流真实速度ν’与全断面上平均渗流速度ν的关系（）转换为ν即可。 对于不可压缩流体，纳维-司托克斯方程为： (2-1) (2-2) (2-3) 上式可写为向量式： (2-4) 上式表明了质量力、流体压力、流动阻力与加速度力的平衡关系，是描述能量守恒的运动方程。当粘滞性ν等于零时，最后一项消去，变为理想流体运动的欧拉方程。 对于多孔介质中的渗流，可把上式中的水质点真实流速ν’改换为全断面平均流速ν除以孔隙率n得到相应的运动方程： (2-5) 由于，将其求导并展开得到： (2-6) 渗流速度及其在各座标方向的导数很小，可以略去，则得到： (2-7) 上式中的单位质量力f，只有一个沿z方向向下的重力，且因，，上式变为： (2-8) 式（2-7）中的最后一项在流体力学中相当于牛顿粘滞性液体的内部摩擦力，而在渗流中则应为液体对土颗粒表面的摩擦力，液体质点之间的摩擦力相对很小却可以忽略。因此可以引用达西定律表示该项仅有的阻力，对于单位质量液体来说，渗流阻力应为沿流程S单位长度的能量损失，即： (2-9) 代入式（2-8）则可得不可压缩流体在不变形多孔介质中的纳维-司托克斯方程： (2-10) 上式一般被称为地下水运动方程。如果是不随时间改变的稳定渗流，上式就简化为重力和阻力控制的达西流动，即： (2-11) 运动方程式（2-10）中包括有两个未知数ν和h，故仍需另一个方程才能求解，它就是下面要介绍的连续方程式。应用这两个基本方程，结合边界和初始条件，就可以计算流速和水头。 2.2.2 连续性方程 图 2-1 微单元体各面上进出流量示意图 地下水运动的连续性方程，可从质量守恒原理出发来考虑可压缩土体的渗流加以引证，即渗流场中水在某一单元体内的增减速率等于进出该单元体流量速率之差。如图（2-1）的微分单元体体积dxdydz，通过左面流进水体质量的速率为，通过右面流出水体质量的速率为，则左右面流出质量之差，即净有的流入量或积累量为。同样对于前后面和上下面作流进流出的流入量计算，最后累加各净有流入量，则得单元体内总的进行流量为： (2-12) 将上式括号内的展开得