1、稳定井流直线图解法确定含水层参数摘要：对于潜水井流，利用 Dupuit 公式计算的参数往往相差很大，难以直接选用。对 Dupuit 公式进行线性化，建立降深(s)与距离(r)的 s-lnr 直线关系，利用直线斜率、直线在 x 轴上的截距，可以直接求得唯一的渗透系数 K 和影响半径 R。利用内蒙古突泉县突泉镇周边河谷地区的承压水井抽水试验资料和内蒙古锡林浩特市 27 个潜水抽水井及其观测孔的抽水试验资料，应用直线图解法分析计算，得到唯一的含水层参数T、K 和 R。将计算结果与直接利用 Dupuit 公式所计算的结果相比较，前者的求参效果较好。 关键词：直线图解法；完整井稳定流；渗透系数 Abst
2、ract: For the diving well flow, using Dupuit formula to calculate the parameters often vary greatly, it is difficult to directly use. Dupuit equation is linearized, establish deep down ( s ) and distance ( r ) s-lnr linear relationship, using the slope of the line, intercept of line in the X line, c
3、an be directly calculated only permeability coefficient K and radius R. Using Inner Mongolia County dash forward fontal town surrounding valley region of artesian well pumping test data of Xilinhaote city and Inner Mongolia 27 submersible pumping well and observation well pumping test data, using li
4、near graphic method of calculation and analysis, is the only aquifer parameters T, K and R. The calculation results and the direct use of Dupuit formula calculation results are compared; the former parameters effect is good. Key words: linear graphic method; complete well steady flow; permeability c
5、oefficient 中图分类号：P641.72 献标识码： A 文章编号：2095-2104（2012）11-0020-02 引言 含水层参数的计算与确定是地下水资源计算与评价、开发利用与保护及科学管理的重要依据。确定含水层参数的方法，通常有抽水试验法1、压水试验法、注水试验法、示踪试验法、室内实验法、动态资料推求法、数值方法和同位素方法等。利用 Dupuit 公式或 Thiem 公式进行稳定井流含水层参数计算十分普遍。在进行参数计算时，是根据现有规范及计算方法的适用条件，采取相应方法进行计算，最常用的是野外稳定流抽水试验法2。利用 Dupuit 公式，将抽水试验达到稳定或趋近稳定时的流量
6、Q 及抽水井或观测孔的水位降深 s 代入公式中，可直接求出含水层渗透系数 K 和影响半径 R3-5。 同一试验场所资料相同、方法不同或资料不同、方法相同所得到的含水层参数值相差很大，这给含水层参数的确定和选用带来了极大的困难。这样确定的参数是否合理、是否符合实际水文地质条件以及是否具有较高的精度一直是专业技术人员非常关心的问题。含水层参数表征了含水层的内部结构特征，是刻画水文地质实体特征的基本参数，它的正确与否将直接关系到研究成果的精度。本文提出一种新的确定含水层参数的方法稳定井流的直线图解法，可唯一地确定一个试验场所的含水层参数，且所求参数能很好地反映实际水文地质条件。 1 基本原理 参数确
7、定的是否合理，直接影响着地下水资源计算与评价、地下水模拟预测的准确性。本文在含水层参数的传统计算基础上，尝试利用直线图解法对稳定流抽水资料进行水文地质参数的确定。 1.1 传统稳定井流求参方法 承压水完整井流的 Dupuit 公式为： (1)或(2) (3)潜水完整井流的 Dupuit 公式为： (4) 或(5) (6) 式中：为含水层的渗透系数(m/d)；为抽水井流量(m3/d)；为抽水井半径(m)； 为抽水井中水位降深（m） ；， ， ， 为观测井距抽水井的距离和井中的水位降深（m） ，其中 i，j= 1，2，3，4，ij；R 为影响半径（m） ；， 为观测孔水位（m） ；，为分别为潜水含
8、水层初始厚度和承压含水层厚度（m） 。 式（1） 、 （4）分别为用抽水井和任一个观测井计算参数的承压水完整井、潜水完整井的稳定流 Dupuit 公式；式（2） 、 （5）分别为用两个观测井计算参数的承压水完整井、潜水完整井的稳定流 Dupuit 公式，又称Thiem 公式6。 1.2 稳定井流的直线图解法 将 Dupuit 公式进行线性化，采用直线图解法确定含水层参数，基本思路如下。 承压水完整井的稳定流 Dupuit 公式中，s- lnr 之间呈线性关系： （6） 潜水完整井的稳定流 Dupuit 公式中，s- lnr 之间呈非线性关系，所以首先要将 Dupuit 公式线性化，令： （7）
9、 对 Dupuit 公式中的降深（任意观测孔的降深）进行修正，得一关于s- lnr 的线性方程： （8） 公式（6） 、 （7）表明在半对数坐标纸上，s-lnr 或 s-lnr 为一条直线，其斜率分别为， ；由此可导出如下 K 值计算公式。 承压水稳定流完整井： （9） 潜水稳定流完整井： （10） （s 或 s）- lnr 直线在横轴 r 上的截距为影响半径 R，即（s 或 s）= 0 处的 r 值为 R。 因此，利用抽水试验中多个观测孔的实测资料，建立（s 或 s）- lnr 直线关系，求出直线斜率，直接带入（9） 、 （10）中求出 K 值；直线在横轴 r 上的截距即影响半径 R，或令
10、s=0，根据方程（6） 、 （8）可直接求出 R。 利用上述的直线图解法能够充分利用所有观测孔的观测资料，所计算出的参数具有唯一性，参数值符合实际情况，从而提高了所用参数的精度。 1.3 公式适用条件 将 Dupuit 公式进行线性化，所得公式的适用条件与原 Dupuit 公式的适用条件相同。原 Dupuit 公式的假设条件为： （1）含水层均质，各向同性，产状水平，厚度不变，分布面积很大，可视为无限延伸； （2）抽水前的地下水面是水平的，并视为稳定的； （3）含水层中的水流服从 Darcy 定律，并在水头下降的瞬间水就释放出来。如有弱透水层，则忽略其弹性释水； （4）水流为水平径向流，即流线
11、为指向井轴的径向直线，等水头面为以井共轴的圆柱面，并和过水断面一致； （5）通过各过水断面的流量处处相等，并等于井的流量。 2 实例应用 2.1 承压完整井稳定流直线图解法求内蒙古突泉县突泉镇周边河谷地区水文地质参数，地层普遍分布有松散岩类孔隙潜水（上部）和孔隙承压水（下部） ，承压含水层分布较均匀，主要由下更新统砂砾石层组成。采用公式（1）（3）对该区 2 个抽水井的抽水试验资料进行求参，所求取的参数不是唯一的，且相差很大，不同资料可以得到不同的参数值，很难确定该区的水文地质参数，这给地下水资源的科学评价和合理开发、利用带来了很大的困难。而应用直线图解法所求的参数却是唯一的，且能很好地反映出
12、该区的水文地质条件（表 1） ，与实际情况相符合。 表 1 承压含水层参数计算结果 2.2 潜水完整井稳定流直线图解法求参 内蒙古锡林浩特市北 40km，面积 500km2，潜水含水层分布较为均匀，主要由中、上更新统砾石和卵石层组成。利用公式（4）（6）对该地区 27 个抽水井及其观测井的试验资料进行参数计算，每两个观测孔可以计算一个 K 和 R，所求取的参数并不是唯一的，同一眼抽水井周围的参数相差很大，所求参数与实际水文地质条件不符，很难从中选择有代表性的水文地质参数。 利用直线图解法求参，以抽水井 zk207 及其观测孔的抽水试验观测资料为例，绘出 s- lnr 线性关系曲线图（图 1）
13、，求取的参数是唯一的，并能真实地反映出该区的水文地质情况（表 2） 。 图 1 zk207 号井 s- lnr 线性关系图 表 2 潜水含水层参数计算结果 2.3 分析与讨论 据计算结果分析，采用承压完整井、潜水完整井的稳定流 Dupuit 公式或 Thiem 公式确定含水层参数时，利用抽水井与观测孔资料求参，渗透系数 K 与观测井到抽水井的距离 r 有关，随 r 的增大，K 和 R 也增大，出现利用不同距离的观测孔可以求得不同的参数的现象。每两个观测孔之间都可以计算出一个参数值，并且距主孔的距离不同，所求取的参数值相差很大。 利用本文提出的直线图解法计算参数，每个抽水井及其观测孔资料所求出的
14、参数不因观测孔的距离不同而不同，参数是相对稳定的且具有唯一性，所求参数基本是抽水井及其周围含水层的具有代表性的参数。因而具有实际意义。 3 结论 （1）利用原始的 Dupuit 公式和直线图解法分别对内蒙古突泉县突泉镇周边河谷地区承压含水层参数 T、K、R 以及内蒙古锡林浩特市潜水含水层参数 K、R 进行计算，计算结果表明，利用 Dupuit 公式计算所得的参数值相差很大，而利用直线图解法计算的参数值具有唯一性，并能很好地反映研究区的水文地质条件。 （2）采用 Dupuit 公式确定含水层参数时，随 r 的增大，K 和 R 也增大，出现利用不同距离的观测孔可以求得不同的参数的现象，并且距主孔的
15、距离不同，所求取的参数值相差很大。利用直线图解法计算参数，每个抽水井及其观测孔资料所求出的参数不因观测孔的距离不同而不同，参数是相对稳定的且具有唯一性，所求参数基本是抽水井及其周围含水层的具有代表性的参数。 参考文献： 1贝尔雅. 地下水水力学M. 许涓铭，等，译. 北京：地质出版社，1985. 2陈雨孙，颜明志. 抽水试验原理与参数测定M. 北京：水利电力出版社，1985. 3杜绍敏，姜英俊，王国春. 抽水试验中的参数计算问题J. 黑龙江水专学报，1994（4）. 4石中平. 水文地质求参实际问题探讨J. 西安理工大学学报，2002，18 （1）. 5张远明. 山前冲洪积扇水文地质参数的选择及计算在广汉市三星堆水源地的应用J. 四川地质学报，2001，21 （1）. 6薛禹群. 地下水动力学M.北京：地质出版社，1997.
