地下水动力学00-第03讲.ppt

3-2-15 3.2.5 等效厚度法 等效厚度的基本思路 以图3-2-3所示的双层结构潜水含水层为例介绍等效厚度法。在保持边界条件不变的前提下，将下层的渗透系数K2化为上层的K1，同时以厚度Md来代替原有的M，以保持其过水断面的过水能力不变。 图3-2-3 等效厚度法算例图 由此形成一个假想的的均质潜水含水层，以代替原有的双层结构的非均质潜水含水层。当然，新形成的假想的均质潜水含水层中的渗流问题必须有现成的解或其解比较容易获得。 等效厚度的计算 等效厚度Md为 K1Md=K2M 3-2-9 即 3-2-10 因此，假想均质含水层1、2端面的含水层厚度分别为 3-2-11 3-2-12 3-2-14 那么，任意断面x的含水层厚度为 已知 所以，水头线方程可写为 流量方程可写成 3-2-13 3.2.6 吉林斯基势函数法 概述 在层状非均质潜水含水层中，当潜水下降穿越层界面时，则难以利用分段法、等效厚度法求解。前苏联学者吉林斯基提出了吉林斯基势函数法，用于求解层状非均质含水层中的地下水流动问题。 吉林斯基势函数法与分段法、等效厚度法相比，它更具有广泛的应用价值。 在研究均质含水层中的地下水流动问题时，一些学者曾经引用“势”的概念来表述地下水的流动： 在均质、隔水底板水平的潜水含水层平面二维流中(引入裘布依假定)，势定义为 在均质、等厚的承压含水层平面二维流中，势定义为 那么单宽流量q的形式则可统一地表示为 3-2-17 3-2-18 若分别用式3-2-15和3-2-16代入式3-2-17、18，则可得到相应的潜水流和承压流的单宽流量公式3-1-10和3-1-7。 3-2-16 吉林基斯势函数的原理 吉林基斯1946年提出吉林基斯势函数Φg解决层状非均质含水层中的流量和水头计算问题。 吉林基斯首先研究透水性在垂向上渐变的含水层中的地下水运动问题(图3-2-4)。假定隔水底板水平，基准面取在隔水底板上z＝0，渗透系数沿垂直方向变化，沿水平方向不变。当平面上流线彼此平行时，h=h(x) 在含水层任一铅垂线上(图3-2-4)，取微分厚度dz，对应K=K(z)，水力坡度水平分量为 ，则通过dz断面的微分单宽流量 3-2-19 整个断面的单宽流量q为 b—含水层的顶面高度，对于承压：b=M；对于潜水：b＝h。 吉林基斯将渗透系数垂向变化的含水层中的流量势定义为 3-2-20 式中： φg—吉林基斯势函数； H—为水头值，从隔水底板算起； b—含水层的顶面高度，对于承压水：b=M，承压含水层厚度；对于潜水：b=h，潜水面高度(从隔水底板算起)。 由莱布尼兹微分法则，对式3-2-20求导得 由此可见，对层状非均质含水层引入吉林基斯势函数后，可以得到与均质含水层同样的单宽流量q的微分形式。 从吉林基斯势函数出发得到，当渗透系数K=常量时，对于均质潜水含水层， ；对于承压含水层， 。吉林基斯势函数应用广泛，既可用于均质，也可用于层状非均质含水层；既可用于潜水流、承压流，也可用于无压—承压流。 将吉林基斯势函数用于层状非均质含水层，则 得 其中：M=M1+M2+…+Mn，而 同理 3-2-21 3-2-22 式中： Ki—为第i层的渗透系数； Mi—为第i层的含水层厚度； zi—为第i层的含水层中点的高度(从隔水底板。 即 根据3-2-22式，先求出断面1的吉林基势函数 ，再求出断面2的 ；然后，利用流量公式3-2-18 求得流量。 图3-2-5 吉林斯基势函数法算例图 吉林基斯势函数的应用 有一水平层状含水层(图3-2-5)。已知：K1=3m/d，M1=6m；K2=5m/d，M2=2m；K3=7m/d，M3=4m；1断面的水头值H1=13m；2断面的水头值H2=4m；1、2断面之间的距离l=600m。求含水层的单宽流量及水头线。 求解单宽流量 根据吉林斯基势函数，第一断面的势为 其中，z1=3m，z2=7m，z3=10m，则 而第二断面的势为 则 求解任意断面的水头值 首先，确定承压转向无压的断面位置。此断面的水头值Hr=M1+M2+M3=12m。接着求出Hr=12m处的断面距离的r值为多少？ 因为Hr已知，所以r断面处的 可求得 ∴ 而 ∴ 水头线方程的求解 把式3-1-23 代入 ，可得