地下水向边界井及不完整井的运动.PPT

第五章 地下水向边界井及不完整井的运动 主要内容 §5.1 镜像原理及直线边界附近的井流 §5.2 扇形含水层中的井流 §5.3 条形含水层中的井流 §5.4 地下水向不完整井运动的特点 §5.5 地下水向不完整井的稳定运动 §5.6 地下水向承压不完整井的非稳定运动 §5.6地下水向承压不完整井的非稳定运动 1.基本方程 设不完整井所在的承压含水层符合推导Theis公式时用的假设条件，并有上覆潜水含水层通过弱透水层发生越流补给。忽略弱透水层的弹性释水，上下含水层初始水头一致。柱坐标系的取发如图5-17所示。 为了研究方便，将越流补给量等效地看作是随坐标变化的含水层内的垂向补给量，这样就可以把上覆透水层视为隔水层。在这些假设条件下，由（1-83）可知，地下水向不完整井的运动满足下列方程： 相应定解条件为： 式中，B为越流因素；L为过滤器长度； d为含水层顶板至 过滤器顶部的距离。 2) 定解问题的解 Hantush给出了上述定解问题，并导出了相应的解。对于过滤器与顶板相接d=0的情况，Hantush给出的几个解。 1）越流含水层中的非稳定流 式中， 式中含有z值，实用时常用沿过滤器长度积分，取平均值的方法消去z，即取平均降深作为降深的近似值，则有： 2) 承压含水层中的非稳定流 当B→∞，无越流补给时，由（5-57）式有 消去z，采用平均降深，由〈5-48〉式有： 3) 越流含水层中的稳定流 当抽水时间很长，趋于稳定流时，在理论上可设t→∞，因而u→0。这时，可有下列近似关系： 把上述结果代入（5-47）式，不难得出稳定流条件下的相 应表达式，代入（5-58）式，则得近似表达式： 其中， 上述各式表明，在非稳定流情况下，降深也由两部分组成， 前者代表相应的完整井降深，后者表示由抽水井不完整性引起 的由抽水井附近流线弯曲所造成的附加降深，它是z的函数。 附加降深的值除与（1）井流量、（2）导水系数有关外， 还与（3）过滤器长度L、（4）不完整程度 和（5）计算断 面到抽水井的相对距离 有关。 式(5-59)、(5-62)和(5-64)所代表的附加阻力称为附加阻力系数，其值随r的增大而减小。 如Hantush指出的，从实用角度看，当r≥1.5M时，这些阻力系数可忽略不计，简化为相应的完整井公式。对于各向异性含水层，这个范围则按 r≥1.5M 确定。 当承压含水层的厚度很大时，允许忽略隔水底板的影响。当M→∞时，由（5-61）式可导出半无限含水层中不完整井非稳定流的相应表达式。它与（5-45）式、(5-46)式很相近，只是对数项的系数略有不同，为1.47而已。 当抽水时间很长，趋于稳定时，只要 ≤0.01,就可利用第四章中提及的近似关系式： 式中， ，即影响半径。 在t很大，u→0时有 把这些关系代入(5-61)式、(5-62)式，便得承压含水层中不完 整井稳定流表达式： 式中， （如表5-1所示）。 表5-1 数值表 如取 则有 ，代入（5-65）可 得相应的有关的表达式。把它和(5-25)式相减，便得只有一个观测孔时的承压含水层中不完整井的稳定流表达式。 3.根据抽水实验资料确定水文地质参数 在不完整井抽水实验中，如在r1.5M区有观测孔，则可根 据相应的完整井公式计算。第四章介绍的各种求参数的方法在这里都有效。 在r1.5M区有观测孔（包括抽水井）时，必须按相应的不 完整井公式计算。因为不完整井的井函数包含变量较多，目前 还没能做出通用的标准曲线供求参数使用，所以下面仅以无 越流补给的不完整井为例，简单介绍据抽水试验资料确定参数 的方法。 在r1.5M区内，如有一个观测孔或多个观测孔的长时间观 测资料，均可用配线法求参数。为此，把(5-61)式简化为： 由此可见，在双对数纸上的 曲线和s-t曲线 或 曲线的形状是相似的。 表5-2给出了常用范围内井函数 值。根据这些资 料，可在双对数坐标纸上绘出三张 标准曲线 (图5-18)。然后，在同模数的透明双对数纸上绘出实测的 曲线或s-t曲线。把实测曲线叠置在相应 某一 标准曲线上，二者应能拟合。重合后，任选一匹配点，并记下该点坐标。如用s-t拟合，