《地下水动力学渗流理论基础3专课件》由会员分享,可在线阅读,更多相关《地下水动力学渗流理论基础3专课件(63页珍藏版)》请在金锄头文库上搜索。
1、16 渗流的连续性方程 在渗流区内以P点取一无限小的平行六面体,其边长分别为x、 y、 z,并且和坐标轴平行,设P点沿坐标轴的渗透速度分量为vx、vy、vz,液体密度为,则P点处,单位时间内通过垂直于坐标轴方向单位面积的水流质量分别为 vx、 vy、 vz。地下水动力学 渗流理论基础3专那么,通过abcd面中点 的单位时间单位面积的水流质量为:用Taylor级数展开:略去二阶导数以上的高次项,得t时间内由abcd面流入单元体的质量为:同理,通过abcd面流出单元体的质量为:地下水动力学 渗流理论基础3专沿x轴方向流入和流出单元体的质量差为:同理,可得到沿y轴和z轴方向流入和流出这个单元体的液体
2、质量差,分别为:地下水动力学 渗流理论基础3专在t时间内,流入与流出这个单元体的总质量差为:在均衡单元体中,孔隙体积为n x y z,其内液体质量为 n x y z, t时间内,单元体内液体质量的变化为:根据质量守恒定律,上二式应相等,因此,消去t得此式为渗流的连续性方程(研究地下水运动的基本方程)(研究地下水运动的基本方程)。地下水动力学 渗流理论基础3专17 承压水运动的基本微分方程 假设条件: (1) 水流服从Darcy定律; (2) K不随= (p)的变化而变化; (3) s和K也不受n变化的影响; (4) 含水层侧向无压缩,即x、 y为常量,只有垂直方向z的压缩。在连续性方程的右端项
3、中,有三个变量,随压力p的变化而变化。 三个变量随时间的变化转化成压力随时间的变化。地下水动力学 渗流理论基础3专 液体压缩后,质量不变。即密度和体积V变化,二者乘积不变。 d( V)= dV +Vd=0 得:由水的压缩系数: 得:所以,d=dp 前面给出了含水层厚度z和孔隙度n随压力p的变化关系: d(z)= zdp ;dn=(1-n) dp 式中:为多孔介质压缩系数。 将三式代入连续方程右端项得:地下水动力学 渗流理论基础3专 于是连续性方程变为:地下水动力学 渗流理论基础3专将 化为 :因为 ,故有:p=(H-z)=g(H-z)或:将d=dp代入,得:即,因为水的压缩性很小, p忽略不计
4、,地下水动力学 渗流理论基础3专代入前式,得第二项非常小,忽略不计,于是上式变为:地下水动力学 渗流理论基础3专根据Darcy定律: 1. 在各向同性介质中,有:代入上式,得因为s=g(+n)所以上式变为:地下水动力学 渗流理论基础3专两边消去单元体体积x yz,得:此式为非均质各向同性介质承压水流微分方程此式为非均质各向同性介质承压水流微分方程。 2. 对于各向异性介质:非均质各向异性介质承压水流微分方程为:地下水动力学 渗流理论基础3专 3. 对于均质各各向同性介质,K为常数,承压水流微分方程为: 4. 地下水流为二维流时,非均质各向同性介质承压水流微分方程为:两边乘含水层厚度M,得或地下
5、水动力学 渗流理论基础3专 5. 柱坐标:如果能用柱坐标表示,则x = rcos、y = rsin,代入可化成式 6. 有源汇项,用W表示。 源:在垂向上有水流入含水层称源。W为正。 汇:在垂向上有水流出含水层称汇。W为负。有源汇项时,只需在上述方程中左边加W即可。如各向同性介质:地下水动力学 渗流理论基础3专7. 稳定流:水位H不随时间变化,即 ,上述微分方程的右端项等于零,即非均质各向同性:非均质各向异性:均质各向同性:二维流,去掉上式中左边第三项。地下水动力学 渗流理论基础3专小结: 承压水三维非稳定流非均质各向同性:非均质各向异性:均质各向同性: 承压水二维非稳定流非均质各向同性:非均
6、质各向异性:均质各向同性:地下水动力学 渗流理论基础3专承压水三维稳定流非均质各向同性:非均质各向异性:均质各向同性: 承压水二维稳定流非均质各向同性:非均质各向异性:均质各向同性:地下水动力学 渗流理论基础3专柱坐标表示的三维流方程:井流方程: 非稳定流: 稳定流: 思考题地下水动力学 渗流理论基础3专18 越流含水层中地下水非稳定运动的基本微分方程 越流含水层(半承压含水层):越流含水层(半承压含水层):当承压含水层的上、下岩层(或一层)为弱透水层时,承压含水层可通过弱透水层与上、下含水层发生水力联系,该承压含水层为越流含水层。 越流:越流:当承压含水层与相邻含水层之间存在水头差时,地下水
7、便会从高水头含水层通过弱透水层流向低水头含水层,这种现象称越流。 假设条件:假设条件: (1) 水流服从Darcy定律; (2) K不随= (p)的变化而变化;地下水动力学 渗流理论基础3专 (3) s和K也不受n变化的影响; (4) 含水层侧向无压缩,即x、 y为常量,只有垂直方向z的压缩。 (5)当弱透水层的渗透系数K1比主含水层的渗透系数K小很多时,近似认为水基本上是垂直地通过弱透水层,折射90后在主含水层中基本上是水平流动的。(如K1与K相差较小时,用等效渗透系数,非越流)。 (6)弱透水层和主含水层释放的水及相邻含水层的越流量相比,弱透水层本身释放的水量小到可以忽略不计。地下水动力学
8、 渗流理论基础3专越流含水层的微分方程: 如图,主含水层的厚度为M,水头为H;上、下各有一厚度为m1和m2、渗透系数为K1和K2的弱透水层。上覆潜水含水层的水位为H1,下伏承压含水层的水位为H2。在含水层全厚度上取一单元体。水平长,宽为x,y。地下水动力学 渗流理论基础3专单元体的中点为P,在P点沿x方向单位时间,通过P点,面积yM的流量为Qx,沿x轴流量的变化率为 。则沿x轴流入单元体的水量为:沿x轴流出单元体的水量为:沿x轴单位时间流入流出单元体的水量差为:同理,可得沿y轴单位时间流入流出单元体的水量差为:地下水动力学 渗流理论基础3专在Z轴方向: 由下部承压含水层单位时间流入越流含水层单
9、元体的水量为:Qz2=v2 x y 向上部含水层单位时间流出越流含水层单元体的水量为:Qz1=v1 x y单位时间内沿Z轴方向流入和流出单元体的水量为:(v2- v1 )x y 地下水动力学 渗流理论基础3专由质量守恒定律有:因为:代入上式,得地下水动力学 渗流理论基础3专消去x y,得:此式为非均质各向同性越流含水层中地下水非稳定运动的基本微分方程。地下水动力学 渗流理论基础3专对于非均质各向异性介质,有: 对于均质各向同性介质,有:地下水动力学 渗流理论基础3专有源、汇项的情况下:非均质各向同性非均质各向异性均质各向同性地下水动力学 渗流理论基础3专 上式中的K/m(用表示)为越流系数。
10、越流系数:当主含水层和供给越流的含水层间水头差为一个长度单位时,通过主含水层和弱透水层间单位面积界面上的水流量。 称为越流因素,量纲为L。 地下水动力学 渗流理论基础3专19 研究潜水运动的基本微分方程一、Dupuit假设 潜水面是弯曲的,等水头面也是弯曲的,潜水流的运动不是水平的。实际上潜水面的坡角很小。 Dupuit假设:假设潜水面比较平缓,等水头面铅直,水流基本上水平,可忽略速度的垂直分量。同一剖面各点的渗透速度相等。地下水动力学 渗流理论基础3专 Dupuit假设是忽略了渗流速度的垂直分量vz,但是,有些地段垂向分速度较大,不能采用Dupuit假设。也就是说,垂向分速度不能忽略。 此外
11、,还有降深较大的抽水井附近,地下水动力学 渗流理论基础3专二、Boussinesq方程 在Dupuit假设和不考虑水的压缩性的条件下。 考虑二维问题(含水层不水平),在渗流场内取一土体。如图。水平宽为x,y。 单元体的中点为P,在P点沿x方向单位时间,通过面积yh的流量为Qx,沿x轴流量的变化率为 。则沿x轴流入单元体的水量为:地下水动力学 渗流理论基础3专沿x轴流出单元体的水量为:沿x轴单位时间流入流出单元体的水量差为:同理,可得沿y轴单位时间流入流出单元体的水量差为: 单位时间内,垂直方向的补给量为: W xy t时间流入流出单元体的水量差为:地下水动力学 渗流理论基础3专 土体内的水量变
12、化引起潜水面的升降。假设潜水面的变化速率为: ,则t时间内,土体内水的增量为: 据质量守恒原理,两个增量应相等。即将:代入上式,得地下水动力学 渗流理论基础3专消去xy t,得:此式为非均质各向同性潜水二维运动的微分方程。 非均质各向异性: 均质各向同性:地下水动力学 渗流理论基础3专 隔水底板水平时 非均质各向同性: 非均质各向异性: 均质各向同性:地下水动力学 渗流理论基础3专 隔水底板水平时,潜水一维流: 非均质各向同性: 均质各向同性: 从微分方程中可知,潜水运动的微分方程是非线性的。地下水动力学 渗流理论基础3专 对于三维流,这是考虑垂向分速度,其微分方程同承压水流微分方程。 非均质
13、各向同性: 非均质各向异性: 地下水动力学 渗流理论基础3专 对于稳定流: 二维流: 非均质各向同性: 非均质各向异性: 三维流: 非均质各向同性: 非均质各向异性:地下水动力学 渗流理论基础3专110 定解条件 前面我们给出了,不同类型地下水的微分方程。分类有:地下水动力学 渗流理论基础3专地下水动力学 渗流理论基础3专地下水动力学 渗流理论基础3专 微分方程是描述地下水运动的方程,不同类型的地下水有相应的微分方程。也就是说,如果只有微分方程,每一类微分方程描述的是这一类地下水运动的所有情况。如:承压水一维稳定流。 或者 因此,仅用微分方程是不能够确切的描述某一地区的地下水运动的。必须同时考
14、虑其它条件的限制。如,补给、径流、排泄条件,以及边界性质、边界形状等。地下水动力学 渗流理论基础3专 因此,除微分方程外,我们还应给出下列条件: (1) 方程中有关参数的值。主要有:渗透系数、导水系数、给水度、贮水率、贮水系数、越流系数,还有源、汇项等。 (2) 渗流区的范围和形状。 (3) 边界条件:渗流区边界所处的条件。(三维流是一个封闭的曲面;二维流是一条封闭的曲线) (4) 初始条件:在某一选定的初始时刻(t=0)渗流区内水头H的分布情况。 定解条件定解条件:就是指边界条件和初始条件。 非稳定流既有边界条件又有初始条件;稳定流只有边界条件。 数学模型数学模型:一个或一组数学方程与其定解
15、条件加在一起,构成一个描述某实际问题的数学模型。地下水动力学 渗流理论基础3专一、边界条件 边界条件分为三类: 1. 第一类边界(也叫给定水头边界)条件:在某一部分边界上,各点在每一时刻的水头都是已知的。用下列式子表示: 三维流: 二维流:式中:H表示S1或1某点(x,y,z)或(x,y)在t时刻的水头;1、2为已知函数。 可作为第一类边界边界的有:(1) 与含水层水力联系好的河流;(2) 泉和抽水井;(3)含水层无限延伸时,影响范围以外的边界。 注意:区别给定水头边界和定水头边界。地下水动力学 渗流理论基础3专 2. 第二类边界(也叫给定流量边界)条件:某一部分边界单位面积上流入(或流出)的
16、流量是已知的。用下列式子表示: 三维流: 二维流:式中:n边界外法线方向;q1、q2为已知函数。 常见的二类边界有: 零流量边界: ,有隔水边界、分水岭和流线。 抽水井和注水井: 泉水。地下水动力学 渗流理论基础3专 3. 第三类边界(混合边界)条件:某边界上水位和法向水力坡度的线性组合是已知的。即 如图,当含水层的边界上有一弱透水层与其它水体或含水层相隔时,常为三类边界。 越过边界单位面积的流量为单位面积越流量为: 二者相等,并化简得地下水动力学 渗流理论基础3专 4. 潜水面边界 在三维流中,潜水面作为边界处理。 稳定流 浸润曲线上压强等于大气压,任意一点的水头等于该点的纵坐标。为一类边界
17、。 非稳定流 二类边界: 非稳定流中,地下水位是随时间变化的,如图。当浸润曲线下降时,从浸润曲线边界流入渗流区的单位面积流量q为:地下水动力学 渗流理论基础3专二、初始条件 初始条件:指某一选定时刻渗流区内各点的水头值。用下列式子表示: 三维流: 二维流:地下水动力学 渗流理论基础3专111 描述地下水运动的数学模型及其解法 一、地下水流问题的数学模型 1. 数学模型 数学模型分随机模型和确定性模型。 随机模型:数学关系式中含有一个或多个随机变量的模型。 确定性模型:数学模型中各变量之间有严格确定的关系。 本书讨论的模型为确定性模型。 描述实际地下水流的确定性模型的必备条件: 有一个或一组描述
18、地下水运动规律的偏微分方程; 相应的定解条件; 识别模型(或校正模型)。地下水动力学 渗流理论基础3专 2. 建立数学模型的步骤 查明地质、水文地质条件; 确定研究区的范围和形状; 对研究区的水文地质条件进行概化; 给出微分方程和定解条件,以及模型中需要的参数; 模型识别。 3. 数学模型必须满足的三个条件: 解的存在性; 解的唯一性; 解的稳定性,指当参数和定解条件发生微小变化时,所引起的解的变化也是微小的。地下水动力学 渗流理论基础3专 4. 例题 例1,研究区地质情况如图,W(x,y,t)代表单位时间、单位面积上的垂向补给量,P(x,y,t)为计划开采区单位面积上的抽水流量,试写出它的数
19、学模型。地下水动力学 渗流理论基础3专 查明地质、水文地质条件; 确定研究区的范围和形状; 对研究区的水文地质条件进行概化; 给出微分方程和定解条件; 地下水动力学 渗流理论基础3专 数学模型为:地下水动力学 渗流理论基础3专 例2,其它条件同例1,在河底有一弱透水层,河水与地下水无直接的水力联系,只是通过弱透水层越流补给地下水。地下水动力学 渗流理论基础3专其数学模型为: 另外,此时的河流可作为二类边界处理(略)。 地下水动力学 渗流理论基础3专 5. 正问题和逆问题 正问题:在数学模拟中,给定含水层的水文地质参数和定解条件,求解水头H,这类问题叫正问题或水头预报问题。 逆问题:根据动态观测
20、资料或抽水试验资料反过来确定水文地质参数,为逆问题或反求参数问题。二、地下水流问题的解法 三种解法:解析法、数值法、模拟法。地下水动力学 渗流理论基础3专本章总结一、渗流的基本概念 1 多孔介质多孔介质,实际应用中能够区分:孔隙介质、裂隙介质、岩溶(Karst)介质。 地下水的流动类型可归纳为两类:渗流和流动。 2 地下水和多孔介质的性质 (1)地下水的状态方程 (2)多孔介质的某些性质 多孔介质的孔隙性:孔孔隙隙度度、有有效效孔孔隙隙、有效孔隙度有效孔隙度。 多孔介质的压缩性。地下水动力学 渗流理论基础3专 3 贮水率和贮水系数 贮水率、贮水系数、二者关系、说明贮水率、贮水系数、二者关系、说
21、明(4点点) 4 渗流 (1)渗流 假假想想水水流流应应有有以以下下特特点点(4点点) ,渗渗流流区区或或渗渗流流场场 (2)典型单元体 5 渗流速度 过水断面、过水断面、渗流速度、渗流速度与实际流速的关系。渗流速度、渗流速度与实际流速的关系。 6 地下水的水头和水头坡度 地下水的水头和水位、等水头面地下水的水头和水位、等水头面 、等水头线、水等水头线、水力坡度。力坡度。地下水动力学 渗流理论基础3专 7 地下水运动特征的分类 稳定流、非稳定流、一维运动、二维运动和三维稳定流、非稳定流、一维运动、二维运动和三维运动。运动。 8 地下水流态的判断 层流和紊流。二、渗流基本定律 1 Darcy定律
22、及其适用范围。定律及其适用范围。 2 渗透系数、渗透率、导水系数、渗透系数与导水渗透系数、渗透率、导水系数、渗透系数与导水系数的关系、渗透系数和渗透率说明(系数的关系、渗透系数和渗透率说明(3点)、渗点)、渗透系数与导水系数说明(透系数与导水系数说明(2点)。点)。 3 非线性运动方程地下水动力学 渗流理论基础3专三、岩层透水特征分类和渗透系数张量 1 岩层透水特征分类 均质岩层、非均质岩层、各向同性、各向异性均质岩层、非均质岩层、各向同性、各向异性 2 渗透系数张量四、突变界面的水流折射和等效渗透系数 1 折射定律、折射定律、4点结论点结论 2 层状岩层的等效渗透系数(水平和垂直)层状岩层的
23、等效渗透系数(水平和垂直)五、流网 1. 流线和迹线流线和迹线 2. 流线的方程 3. 流函数: 流函数有下列特性流函数有下列特性地下水动力学 渗流理论基础3专 4 流网及其性质 流网、流网的性质流网、流网的性质 5 流网的应用流网的应用(做到会作)六、渗流的连续性方程七、承压水运动的基本微分方程 假设条件、基本方程假设条件、基本方程。八、越流含水层中地下水非稳定运动的基本微分方程 越流含水层、越流、越流系数、越流因素越流含水层、越流、越流系数、越流因素九、研究潜水运动的基本微分方程 Dupuit假设、基本微分方程假设、基本微分方程地下水动力学 渗流理论基础3专十、定解条件 边界条件、初始条件
24、、数学模型、第一类边界边界条件、初始条件、数学模型、第一类边界(也叫给定水头边界也叫给定水头边界)条件、第二类边界条件、第二类边界(也叫给定也叫给定流量边界流量边界)条件、第三类边界条件、第三类边界(混合边界混合边界)条件条件十一、数学模型及其解法 1 建立数学模型的步骤 2 正问题和逆问题。 3 地下水流问题的三种解法:解析法、数值法、模拟法。地下水动力学 渗流理论基础3专课堂作业: 1. 一口井位于无限分布的均质、各向同性潜水含水层中,初始时刻潜水位在水平不透水底板以上高度为H0(x,y),试写出下列两种情况下流向井的非稳定流数学模型。已知水流为二维非稳定流。(1)井的抽水量保持不变;(2)井中水位Hw保持不变。地下水动力学 渗流理论基础3专2. 如图,为非均质各向同性含水层。开采后,采区仍属承压含水层,且未影响到两侧的流线边界。已知开采强度为,上游区的补给强度为W,水流为非稳定二微流,试写出两条流线间渗流区的数学模型。地下水动力学 渗流理论基础3专