第四章-地下水向完整井的非稳定运动-2.pptVIP

§2 有越流补给的完整井流 一、基本方程 越流系统：有越流含水系统，弱透水层和相邻的含水层共同组成的系统； 假设条件： （1）越流系统中每一层都是均质各向同性，无限延伸的；含水层底部水平，含水层和弱透水层都是等厚的； （2）含水层中水流服从Darvy定律； （3）虽然发生越流，但是相邻含水层在抽水过程中水头保持不变； （4）弱透水层本身的弹性释水可以忽略，通过弱透水层的水流可视为垂向一维流； （5）抽水含水层天然水力坡度为零，抽水后为平面径向流； （6）抽水井为完整井；井径无限小，定流量抽水。 教学模型及其解： 其解为： 二、公式讨论 1. 降深—时间曲线的形状 上述公式变换为： 利用表中的数据，可作的W(u)—1/u 曲线，如下图。图中最外曲线是无越流的Theis曲线，其他曲线为有越流的曲线，而且s—t曲线与其中的某一曲线相对应。 因此，有越流的曲线可分为三段： （1）抽水早期：降深曲线同Theis曲线一致。表明越流尚未进入主含水层，抽水量几乎全部来自主含水层的弹性释水。 另外，在其它条件一定时,r/B的数值不同，曲线与Theis曲线吻合的时间也不一样，r/B小时（即越流系数σ′=K1/m1越小），吻合的时间长，也就是说，越流系数越小，越流开始进入含水层的水流越晚。 （2）抽水中期：降深曲线开始偏离Theis曲线，说明越流已经开始进入抽水含水层，这时，抽水量由两部分组成：一是抽水含水层的弹性释水，二是越流补给。 （3）抽水后期：曲线趋于水平，即降深不随时间变化，由非稳定流转化为稳定流。这时抽水量全部由越流补给，这时，上方程变为稳定流方程： 2. 水头下降速度 前面推出无越流时的公式： 相比，越流含水层水位下降速度比无越流含水层慢。 当t足够大时， 在一定范围内，水位下降速度是相同的。 三、 利用抽水试验资料确定越流系统的参数 1. 配线法： 两边取对数，得： 步骤： （1）在双对数坐标纸上绘制W(u,r/B)—1/u标准曲线； （2）在另一张模数相同的透明双对数纸上，点绘s—t实测数据曲线； （3）在保持对应坐标轴彼此平行的前提下，相对移动两坐标纸，在标准曲线中找最优重合曲线； （4）两曲线重合后，任选一配点，记下对应的四个坐标值W(u,r/B)、1/u、s、t，代入上式求得： （5）从曲线上都得r/B，可求得： 2. 拐点法 （1）原理 ① 拐点的斜率 前面，水头下降速度中，给出 则 由上式可看出，同一观测孔的s—lgt曲线的斜率变化规律是由小到大，又由大变到小，存在着拐点。 拐点位置确定： 设拐点为P ，上式对lgt求二阶偏导，得 故在拐点有： 解的拐点处的时间tp为： 相应的u值为： 将此二式代入式： 的 即拐点处的斜率为： ② 拐点处降深 ③ 拐点P处降深sp与斜率ip的关系 前二式相除，得 这三个公式，就是求参用到的公式，即 （a）有一个观测孔 步骤： ① 在单对数坐标纸上绘制s—lgt曲线，如图。用外推法确定最大降深Smax，并用下式计算拐点处的降深Sp。 ② 根据Sp确定拐点位置，并从图上读出拐点出现的时间tp。 ③ 作拐点p处曲线的切线，并从图上确定拐点p处切线的斜率ip。 ④ 根据下式 求 ，求得后查表4-7，确定r/B和 值； ⑤ 根据r/B值，求B值： 按下二式求T和μ*值： ⑥ 验证，因为图解出的Smax和Sp常有较大随意性而引起误差，所以进行验证是必要的。将求得参数，代入下式： 给出不同的t值，计算理论降深，然后把它同实测降深比较，如果不吻合，则应该重新图解计算。 （b） 有多个观测孔时 抽水时间较短时，抽水未趋于稳定，这时难以确定Smax，应用此法，如图。 根据拐点斜率公式， 两边取对数： 化为： r与lgip呈线性关系，可利用r和ip资料求参， 条件：必须由三个以上观测孔。 步骤： （1）汇每个观测孔s—lgt的曲线，并从图上确定每条曲线直线段的斜率ip，近似地代替拐点处的斜率。 （2）根据各孔的斜率作r—lgip曲线，应为一直线。从图上确定直线的斜率ip，代入下式求B。 （3）将r—lgip直线延长交横轴于一点，读得这一点的(ip)值，代入下式，求T、σ?： ? （4）利用所求得B值，得r/B，查表得K0(r/B)，和所求得T一起代入式： 计算出不同观测孔的拐点降深sp。利用s