Seep(渗流计算)V30使用手册.docVIP

?????????????????????????????????? ?????????????????????????? ??? 二??四年一月.1水工渗流的危害及渗流分析计算的任务 流体在多孔介质中的运动称为渗流。水是最为常见的流体，水利水电工程中由于广泛建造堤、坝、围堰、水闸等挡水建筑物形成了水头差，这些建筑物或其地基通常是透水的多孔介质，因此水工渗流现象十分普遍。 水工渗流造成多方面的危害。渗流造成水库、渠道水量损失；渗流使堤坝、围堰土体饱和，降低坝体的有效容重和抗剪强度，可能导致坝坡失稳；建筑物地基渗流对建筑物底部产生扬压力，也不利于建筑物的稳定；渗流流速过大时，还可能造成坝体或建筑物地基的土体颗粒流失，发生渗透变形，从而使堤坝崩塌或建筑物滑移、倾覆；水库渗流还可能引起下游地下水位升高，导致农田冷浸渍害、盐碱化，使作物减产；拦污坝渗流造成地下水环境污染。 水工渗流分析计算的任务就是要研究水在渗流区域的渗流流速、流量、水头分布及浸润线等，从而为采取合理的渗流控制措施提供依据，以避免或减缓渗流危害。 1.2达西定律 19世纪50年代，法国工程师亨利·达西（H.Darcy）h，反比于渗径长度L，且与土粒结构及流体性质有关。这就是著名的达西定律，可用公式表达为： (1.2.1) 式中h—测压管水头，，z为位置高度，p为压强，γ为水的容重。因为渗流的流速一般很小，流速水头可忽略，故。 k—反映土粒结构及流体性质的系数，即渗透系数，对于某一具体的流体（比如水）而言，k值仅与土粒结构有关。 J—渗透坡降，。 式中的负号“-”表示水总是流向水头减小的方向。 应当注意，达西定律中的流速是全断面上的平均流速v，而不是土体孔隙中的流速，这两种流速存在以下关系： (1.2.2) 式中n为体积孔隙率，可见达西流速小于土体孔隙中的流速。 还应注意，达西定律只能适用于层流状态的渗流运动。在水利工程中，除了堆石坝、堆石排水体等大孔隙介质中的渗流为 流之外，绝大多数渗流都属于层流，达西定律都可适用。对于非层流渗流，其流动规律可用以下公式形式表达: （1.2.3） 上式中当m=1时，为层流渗流；当m=2时，为完全紊流渗流；当1＜m＜2时，为层流到紊流的过渡区。 将式（1.2.1）等号两边向x、y、z轴投影，便得到空间直角坐标系中的达西公式： （1.2.4） 1.2.3渗流运动连续性方程 图1.2.1 微元体进出流量示意 从渗流场中取某一微分单元体（如图1.2.1），其体积，流入左侧面的水体质量的速率为，流出右侧面的水体质量的速率为，则左右面进出流量之差为 同样，对于前后面和上下面作流进流出的流入量计算，最后累加各净有流入量，得到微元体内总得流入量为 将上式展开为 式中后一项与前一项相比，小得可以忽略，故改写成 （1.2.5） 上式即为水体质量在微元体内积累的速率，根据质量守恒原理，它应等于微元体内水体质量M随时间的变化速率： （1.2.6） 式中，n为土体的孔隙率，ρ为水的密度，V为微元体的体积。式（1.2.6）右边三项分别代表土体骨架、孔隙体积及流体密度的改变速率。 引入弹性压缩理论，可导出（1.2.6）式中得等价表达为： （1.2.7） 式（1.2.7）中 α——土颗粒骨架的压缩性（即压缩模量） β——水的压缩性 h——渗流测压管水头（m） 根据质量守恒原理，式（1.2.5）与式（1.2.7）相等，则得到 （1.2.8） 假定水和土不可压缩时，上式变为： （1.2.9） 式（1.2.9）为不可压缩流体在刚体介质中流动的连续性方程。 1.2.4渗流微分方程 将达西公式代入渗流连续性方程（1.2.8）可得： （1.2.10） 式（1.2.10）为非稳定渗流微分方程的一般形式，既适合于承压含水层，也适合于无压渗流。式中称为单位储存量（尺度为1/L），其含义是：单位体积得饱和土体在水头下降1m时，由于土体压缩（）和水体膨胀（）所释放出来得储存水量。 当介质为均质各向同性（即）时，式（1.2.10）变为 （1.2.11） 当假