承压含水层中非达西径向稳定流动抽水模拟试验.pdfVIP

水 利 学 报 SHUILI XUEBAO2012年 1月 第 43卷 第 1期 文章编号：0559-9350（2012）01-0076-08 承压含水层中非达西径向稳定流动抽水模拟试验 李 健 1，2，3，黄冠华 1，3，文 章 4，詹红兵 4，5，6 （1. 中国农业大学 水利与土木工程学院，北京 100083；2. 农业部规划设计研究院，北京 100125； 3. 中国-以色列国际农业研究培训中心，北京 100083；4. 中国地质大学 环境学院，湖北 武汉 430074； 5. 中国地质大学 工程学院，湖北 武汉 430074；6. 德克萨斯农工大学 地质与地球物理系，美国） 摘要：为分析抽水井附近水流的非达西流动，采用装填均质河砂的扇形密闭水槽模拟承压含水层，开展了外边界 水头恒定条件下的定流量抽水试验，获得到了水流稳定流动时的水头降深曲线。通过比较试测水头降深与达西流 动理论水头降深，以及比较相邻两压力传感器之间的平均水力梯度与相同流速下的达西水力梯度，进行了非达西 流动的判别，比较结果发现本试验条件下存在明显的非达西流动现象。同时，应用基于 Izbash方程的非达西稳定 流动水头降深解析解，对不同流量下的试验数据进行了拟合分析，发现抽水流量对非达西稳定流动的 Izbash方程 系数值影响较小。通过不同平均方法获得了描述本试验中非达西流动的 Izbash方程系数的适宜取值为 k=7.17× 10-7，n=1.65，此时不同抽水流量下拟合水头降深曲线与试测水头降深曲线吻合良好，表明基于 Izbash方程的非 达西稳定流动水头降深解析解可较好地描述抽水井附近的非达西径向稳定流动。 关键词：非达西流动；水头降深；抽水模拟试验；Izbash方程 中图分类号：TV139.16 文献标识码：A 1 研究背景 抽水井附近承压含水层中水流流动的水力梯度和水流流速较大且变化剧烈，是发生非达西流动 的典型情形之一［1-2］，非达西流动的存在可能影响抽水井的抽水效率和水头降深曲线［3］。目前，国内 外学者多采用非线性方程——Forchheimer方程［4］或 Izbash方程［5］——描述非达西流动，通过推导特定 边界条件下的井流函数，对非达西流动及其影响进行相应的理论研究［6-10］。 室内渗流试验也是研究非达西流动的主要方法之一。这种试验主要包括一维渗流试验和径向渗 流试验两类［3，11-14］，而在日常生产试践中后一种流动问题更为多见，因此该问题已受到越来越多研究 者的重视。Wright［15］开展了径向渗流试验，分析了雷诺数与摩擦系数之间的关系，结果表明随雷诺数 增大，径向水流存在非达西流动现象。考虑到径向流动时水力梯度与流速随半径不断变化，Thiru? vengadam等［16］研究得到了描述径向非达西流动的 Forchheimer方程。Venkataraman等［17］结合径向流动 的边界条件，发现描述径向非达西流动的 Forchheimer方程的系数仅受径向边界条件的影响。然而， 上述试验研究均假设径向水流存在非达西流动，但未见利用试验资料对此假设进行验证，且非达西 流动方程主要根据水力梯度与流速之间的关系确定，这种方法并不适用于水力梯度与流速随位置不 断变化的径向流动问题。此外，不同抽水流量对径向非达西流动影响的试验研究还鲜有报道，同时 抽水井附近非达西流动理论也有待利用试测资料进行验证。 收稿日期：2009-12-09 基金项目：国家自然科学基金项目50979106；教育部新世纪优秀人才支持计划（NCET-05-0125），长江 学者创新团队计划与北京市重点学科资助计划 作者简介：李健（1983-），男，河北易县人，博士，主要从事非达西流动研究。E-mail：lijcau@ 通讯作者：黄冠华（1966-），江西武宁人，博士，教授，主要从事多孔介质中水分与溶质运移模拟研究。 E-mail：ghuang@ —— 76 图 2 压力传感器沿扇形槽侧壁布置（单位：cm） 为此，本文开展了承压含水层中外边界水头恒定条件下的定流量抽水模拟试验（承压含水层中圆 岛模型［18］），通过分析承压含水层中非达西稳定流动时的试测水头降深，对径向非达西流动进行判 别，并采用非达西稳定流动的解析解对试测降深资料进行拟合，反求描述承压含水层中非达西稳定 流动的 Izbash方程系数。 2 材料与方法 2.1 试验装置 本文设计了一套模拟抽水井附近承压含水层中水流非达西流动的试验装置，整个装 置由扇形水槽、数据采集装置、进水装置和出水装置 4部分组成，如图 1所示。 扇形水槽由钢板焊接而成，半径长 110cm，圆心角 60°，高 30cm。水槽顶盖与槽身由法兰连接可 保证密封，装填介质后用于模拟承压含水层。水