承压含水层中大口径井附近非达西流的线性化近似解与数值解 Approximate analytic and numerical solutions for non-Darcian flow towards large-diameter well in confined aquifers.pdfVIP

水 利 学 报 2009年7月 SHUIUXUEBAO 第40卷第7期 文章编号：0559．9350(2009)07．0863．07 承压含水层中大口径井附近非达西流的线性化近似解与数值解 文 章1’2，黄冠华1’2，李健1’2，詹红兵3 (1．中国农业大学水利与土木-I-程学院，北京100083；2．中国．以色列国际农业研究培iJiI中心，北京100083； 3．德克萨斯农工大学．美国) 摘要：研究了承压含水层中大1：3径井附近的非达西流．采用线性化方法和Laplace变换得到了无量纲水位降深在拉 氏空间下的解析解，采用Stehfest数值逆变换得到了实空间下的水位降深，同时得到了抽水后期和抽水稳定阶段水 位降深的解析解。定义了反映非达西流动特征的水力传导度e，采用相应的数值解对线性化近似解进行了验证。 研究结果表明。在抽水初期，不同亭值所对应的井筒中的水位降深基本相同，并在双对数坐标下表现为一直线；在 抽水后期，无论是在井简还是含水层中，e越大，水位降深越小；在抽水稳定时期，水位降深随距离增大并呈幂函数 递增，线性化方法在抽水初期会带来一定误差，而在抽水后期结果比较精确。 关键词：非达西流；lzbash定律；大121径井；含水层；水位降深；Laplace变换；线性化近似 中圈分类号：TVl39．1 文献标识码：A 1研究背景 在发展中国家，大口井(指抽水条件下井中的储水效应不能忽略时的抽水井)已被广泛应用于地下 水开采中，目前关于大口井的水力学问题的研究都是以达西流为基础的¨q1，但在孔隙较大的多孔介质 含水层或裂隙较为丰富的含水层中，大口井附近易发生非达西流¨。41，即水力梯度与水流通量之间的关 系为非线性的水流。关于抽水井附近的非达西流问题，国内外很多学者分别从实验¨“1或理论上口’7q] 达西流问题进行了研究，但其仅对井筒中的水位降深进行了分析，而含水层中的水位降深则没有进行分 析。Wen等归。分析了非达西流情况下大口井中和含水层中的水位降深规律，同时Wen等∽1也指出采用 Boltzmann变换求解非达西流问题可能会带来一定的误差。 本文采用线性化和Laplace变换相结合的方法求解大口井附近非达西流情况下的水位降深规律，为 利用大口井开采地下水提供相应的理论依据，抽水后期以及抽水稳定阶段的近似解析解可以用来估计 含水层参数，最后采用有限差分得到了相应的数值解，并对线性化近似解进行了验证。 2问题的提出与求解 致，如图1所示，假设含水层均质各向同性，抽水前整个含水层中水位降深相等且都等于零，并假设抽水 井为定流量Q。 考虑大IzI井井储的影响，定解问题可以如下描述： 牧稿日期：2008．03．30 市重点学科水文学与水资源建设项目 作者简介：文章(1982一)，男，湖北公安人，博士生，主要从事地下水非达西流动的研究。E-mail：m崦¨@印，ail．咖 —-——863‘_—— 万方数据 ㈩ 吗≯+掣：曼m芎≠dr r df s(r，0)=0 (2) s(∞，t)=0 (3) ＼． ．y 互 丑 2兀rw ，／，／／，／／／／／／／，， 。／／／，，／，，／，，，，，． mq(r洲一．一矶2掣：一Q(4) 式中：r为离抽水井中心距离；t为抽水时间；S为含 水层的储水系数；m为含水层的厚度；q(r，￡)为在距 E 杰