承压含水层中抽水井附近非达西流的近似解析解.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 承压含水层中抽水井附近非达西流的近似解析解 文 章 黄冠华 李 健 詹红兵 中国农业大学水利与土木工程学院 北京 中国 以色列国际农业研究培训中心 北京 德克萨斯农工大学地质与地球物理系 美国 摘要 本文将线性化方法和 变换相结合研究承压含水层中单一抽水井附近的非达西流问题 得到了水位降 深在抽水后期和抽水稳定阶段的近似解析解 并对抽水后期近似解析解的适用性进行了讨论 通过 数值 逆变换得到了任意时刻任意位置水位降深的半解析解 并采用数值解对线性化方法所得到的近似解进行了 验证 研究结果表明 在抽水后期 水位降深随着 定律中的两个常数的增大而减小 在抽水后期 水位降深近 似为时间的幂函数 在抽水稳定阶段 水位降深可以近似为距离的幂函数 在抽水后期 线性化方法所得到近似解 与数值解吻合很好 而在抽水初期线性化方法则存在一定误差 会低估水位降深 关键词 非达西流 定律 变换 线性化 中图分类号 文献标识码 研究背景 达西定律己被广泛应用于解决地下水渗流问题 然而研究发现 当水流流速过快或过慢时 达西定 律所描述的水流通量与水力梯度之间的线性关系将不复存在 所有偏离达西定律的流动问题都可 以归结为非达西流 非达西流可以分为两大类 前线性流和后线性流 前线性流是指由于流速较慢而 导致偏离达西定律的流动 主要发生在石油工程或弱透水层中的水流中 后线性流是指由于流速过快 而偏离达西定律的流动 比如抽水井附近 有关研究 建立了一系列描述非达西流的经验公式 但应 用最为广泛的是 定律 和 定律 常安定等 在 定律的基础上采用 变换方法得到了两种流态区域 非达西与达西 条件下井流的解析解 刘元会等 和王鹏举 也对非达西流问题进行了相关研究 分别对不同条件下抽水井附近的非达西流问题进行了研 究 采用 变换方法得到了相应的解析解 等 在实验室发现单一裂隙介质中很容易发生 非达西流 等 也采用 变换方法得到了单一裂隙介质中抽水井附近非达西流情况下的解 析解 采用数值方法得到了非达西流情况下的数值解 由于非达西流的高度非线性 通常情况下关于非达西流的解析解较难获得 尽管很多学者 采 用 变换得到了非达西流情况下的解析解 但 变换在求解非达西径向流 存在一 定的问题 主要是因为边界条件不能表示为 变量 ? 的完全函数 导致了积分常数是变 量 或时间 的函数 并进而导致 变换的解析解产生一定的误差 本文基于 定律描述 非达西流的条件下 采用线性化方法和 变换来求解非达西流问题 得到了相应的近似解析解 并 将该结果与 变换解进行比较 收稿日期 基金项目 国家自然科学基金项目 教育部新世纪优秀人才支持计划 教育部创新团队计划 作者简介 文章 男 湖北公安人 博士生 主要从事地下水非达西流动的研究 问题的提出与求解 考虑承压含水层中抽水井附近的水流运动 如图 所示 假设含水层均质各向同性 抽水前整个含 水层中水位降深相等且都等于零 在这种情况下 定解问题可以如 下描述 式中 为离抽水井中心距离 为抽水时间 为含水层的储水系 图