越流含水层抽水井附近非达西流与达西流区界面位置变化规律研究.PDF

水 利 学 报 2012年10月 SHUILI XUEBAO 第43卷 第10期 文章编号：0559-9350（2012）10-1171-08 越流含水层抽水井附近非达西流与达西流区 界面位置变化规律研究 王全荣 ，唐仲华 ，文 章 ，詹红兵1 1 1 2，3 （1. 中国地质大学（武汉）环境学院，湖北 武汉 430074；2. 中国地质大学（武汉）工程学院，湖北 武汉 430074； 3. 德克萨斯农工大学 地质与地球物理系，美国） 摘要：为了研究越流含水层系统中大口径完整井附近地下水流动规律，根据Reynolds数将渗流场分为非达西流与 达西流两个区域，提出了一种基于有限差分原理的迭代法，来模拟两区界面位置（RcD ）随时间的变化规律，并采 用实测数据验证本文解的实用性。结果表明：抽水初期，RcD 很小，本文解与全达西流解析解几乎重合；抽水中 期，RcD 逐渐增大，解析解与数值解的降深曲线差别很明显，而在非达西流区域，本文解与全非达西流解析解的 降深曲线斜率差别不大；在抽水后期，地下水流场达到稳定，RcD 达到最大值，同时，在非达西流区域中，本文 解与全非达西解吻合较好，在达西流区域，本文解与全达西流解吻合较好。RcD 对井筒半径和滤水管半径比较敏 感。当RcD 趋于无穷大时本文解转化为全非达西流解；当RcD 趋于0时本文解转化为全达西流解。 关键词：非达西流；边界条件；有限差分法；Izbash定律 中图分类号：TV139.14 文献标识码：A 1 研究背景 对于渗透性较弱的含水层，为了增强抽水井的储水效应，通常将距地表较近的井筒半径扩大(称 ［1］ 大口径井），这种井已经被广泛地应用于地下水开采中，尤其是在中东与亚洲地区 。由于大口径井 ［1］ 筒结构与普通竖井存在一定的差异，其水动力特征也有所不同，Papadopulos和Cooper 于1967年首 ［2］ 次采用Laplace变化法推导了承压含水层大口径井附近水头降深解析解。类似地，Lai和Su 采用同样 的方法推导了越流含水层系统中大口径井附近水头降深的解析解。由于解的复杂性，Laplace变化法 ［3］ 很难得出实空间下的解析表达式，于是Cimen 针对这个问题给出一种新的求解方法。然而，上述研 究的一个重要假设是地下水流动满足达西定律。 ［4-6］ ［7］ 事实上，由于井附近的水力梯度较大，水流可能不满足达西定律 ，文章等 推导了承压含水层 中大口径井附近非达西流近似解析解，其中假设整个含水层中的地下水流都满足非达西流定律。对于 ［8］ 径向流，地下水流速随距井中心距离的增大逐渐减小，因此，远离井的水流往往是达西流。Sen 将整 个渗流区域分为两个子区域（距离井中心某一距离的区域为非达西流区，此外为达西流区），并采用 ［9］ Boltzmann变换法推导了不同条件下抽水井附近非达西流的解析解。常安定等 指出Sen的解析公式适 用于充分紊动的非线性流情况，于是推导非线性指数小于等于2情况下的解析解。同时，刘元会等［10］ ［11］ 和王鹏举 也对非达西流问题进