抽水试验确定渗透系数的方法及步骤要点

1、抽水试验确定渗透系数的方法及步骤1 .抽水试验资料整理试验期间，对原始资料和表格应及时进行整理。试验结束后，应进行资料分析、整理，提交 抽水试验报告。单孔抽水试验应提交抽水试验综合成果表，其内容包括：水位和流量过程曲线、水位和流量关系曲线、水位和时间（单对数及双对数）关系曲线、恢复水位与时间关系曲线、抽水成果、 水质化验成果、水文地质计算成果、施工技术柱状图、钻孔平面位置图等。并利用单孔抽水试验资料编绘导水系数分区图。多孔抽水试验尚应提交抽水试验地下水水位下降漏斗平面图、剖面图。群孔干扰抽水试验和试验性开采抽水试验还应提交抽水孔和观测孔平面位置图（以水文地质图为底图卜勘察区初始水位等水位线图、

2、水位下降漏斗发展趋势图（编制等水位线图系列 卜水位下降漏斗剖面图、水位恢复后的等水位线图、观测孔的S t、S-lg t曲线注、各抽水孔单孔流量和孔组总流量过程曲线等。注意：（1）要消除区域水位下降值；（2）在基岩地区要消除固体潮的影响；3）傍河抽水要消除河水位变化对抽水孔水位变化的影响。多孔抽水试验、群孔干扰抽水试验和试验性开采抽水试验均应编写试验小结，其内容包括：试验目的、要求、方法、获得的主要成果及其质量评述和结论。2 .稳定流抽水试验求参方法求参方法可以采用 Dupuit公式法和Thiem公式法。（1）只有抽水孔观测资料时的 Dupuit公式承压完整井：2开力M q式中K含水层渗透系数

3、（m/d）;Q抽水井流量 （m3/d）;sw抽水井中水位降深 （m）;M承压含水层厚度（m）;R 影响半径（m）;H 潜水含水层厚度（m）;h潜水含水层抽水后的厚度（m）;rw抽水井半径 （m）。（2）当有抽水井和观测孔的观测资料时的Dupuit或Thiem公式承压完整并、hL-hv =132 戒M rThiem 公式：h二 0 1口/2ttKM A潜水完整井工Thiem公式：煨一卜；=-；;1口工 戒Mq式中hw抽水井中水柱高度（m）;hl、h2与抽水井距离为ri和r2处观测孔（井）中水柱高度（m）,分别等于初始水位H0与井中水位降深 s之差，h1= H0 -si; h2= H0 s2。其余

4、符号意义同前。当前水井中的降深较大时，可采用修正降深。修正降深s'与实际降深 s之间的关系为：s'=s-s2/2H。3.非稳定流抽水试验求参方法3.1 承压水非稳定流抽水试验求参方法（1）Theis配线法在两张相同刻度的双对数坐标纸上，分别绘制Theis标准曲线W（u）-1/u和抽水试验数据曲线s-t,保持坐标轴平行，使两条曲线配合，得到配合点M的水位降深s、时间t、Theis井函数w（u）及1/u的数值，按下列公式计算参数（r为抽水井半径或观测孔至抽水井的距 离）：7的YI 网4”订以上为降深 一一时间法（s-t）。也可以采用降深 一时间距离法（s-t/r2）、降深距离法（s

5、-r） 进行参数计算。（2） Jacob直线图解法当抽水试验时间较长，u= r2/（4at）0.01时，在半对数坐标纸上抽水试验数据曲线s-t为一直线（延长后交时间轴于t0,此时s=0.00m）,在直线段上任取两点t1、si、t2、s2,则有225%S - 7-a -225名(3) Hantush拐点半对数法对半承压完整井的非稳定流抽水试验(存在越流量,K' /b为越流系数)，当抽水试验时间较长，u= r2/(4at)<0.1时，在半对数坐标纸上抽水试验数据曲线s-t,外推确定最大水位降深Smax,在s-lgt线上确定拐点Si = Smax/2 ,拐点处的斜率 mi及时间ti,则

6、有。（力求"yBD D2Tls- B rb' B(4)水位恢复法当抽水试验水位恢复时间较长，u= r2/(4at)<0.01数据曲线s-t,在直线段上任取两点t1, s1, t2时，在半对数坐标纸上绘制停抽后水位恢复 s2,则有 S2产，入a=10225tlTs = a(5)水位恢复的直线斜率法当抽水试验水位恢复时间较长，u= r2/(4at)<0.1数据曲线s-t,直线段的斜率为B,则有时，在半对数坐标纸上绘制停抽后水位恢复23。4霏B1%3.2 潜水非稳定流抽水试验求参方法 潜水参数计算可采用仿泰斯公式法、Boulton法和Numan法。（1）仿泰斯公式法Wu

7、22rr iiM A-at47Y式中H0、hw -初始水头及抽水后井中水头；W（u） 泰斯井函数；Q抽水井的流量（m3/d）;r到抽水井的距离（m）;t自抽水开始起算的时间（d）;T含水层的导水系数（m2/d）; T=Khm ;hm 潜水含水层的平均厚度 （m）;K含水层的渗透系数（m/d）;A 含水层的导压系数（1/d）;m 潜水含水层的给水度。具体计算时可采用配线法、直线图解法、水位恢复法等。（2）潜水完整井考虑迟后疏干的Boulton公式抽水中期：占=昌女式。）2万7DQ . r抽水晚期3券避,可根据抽水早期、中期、晚期的观测资料，采用相应的方法计算参数。3 3） Numan 法对于潜水

8、含水层完整井非稳定流抽水试验，也可以采用Numan模型求参，具体求参过程可参阅地下水动力学等教科书。4 .参数计算新技术新方法的应用采用AQUIFERYTEST 软件（图1）、数值模拟法（可采用 GMS、MODFLOW、FEFLOW 等软件）以及肖长来教授提出白全称曲线拟合法（图2）等一些新的软件、方法确定水文地质参数，效果非常好。CT井第一次抽水（HantushUmEb）1EYr1E+21E+7Conductivity:9.38E-2 m/d图1 AQUIFERYTEST 软件求参图示抽水廿间tffdiQ图2全称曲线拟合法求参图示5 .参数计算结果的验证上述参数计算结果的精度如何，取决于试验

9、场地水文地质条件的概化，也取决于观测数据的精度。对于所求得的参数，应将其代入相应的公式，通过对比计算降深与实测降深的差值， 分析所求参数的精度及其可靠性和代表性，最终确定抽水试验场地的有代表性意义的参数 值。方法（二）单孔稳定流抽水试验，当利用抽水孔的水位下降资料计算渗透系数时，可采用下列公式：1当Qs （或A h2）关系曲线呈直线时，1）承压水完整孔：K - - Jn 二 R"(8.2.1-1)2)承压水非完整孔：当 M>150r , l/M>0.1 时：K = - (In + jIn )二,：:(8.2.1-2)或当过滤器位于含水层的顶部或底部时：K 二 1 In +

10、 ;In 11+0,2 J2782M L尸/厂(821-3)3)潜水完整孔：3 Q (K 二八M六1口7T(H2 - A2)(821-4)4)潜水非完整孔： 当>150r, l>0.1 时：(821-5)7T (H2 一 A2)或当过滤器位于含水层的顶部或底部时: In 1 +0.2（821-6）式中 K渗透系数（m/d）;s-M- 出水量 （m3/d）;水位下降值（m）;- 承压水含水层的厚度（m）;- 自然情况下潜水含水层的厚度（m）;- 潜水含水层在自然情况下和抽水试验时的厚度的平均值(m)；rR潜水含水层在抽水试验时的厚度过滤器的长度（m）;-抽水孔过滤器的半径（m）;影响

11、半径(m)。(m)；2当Qs （或Ah2）关系曲线呈曲线时，可采用插值法得出Qs代数多项式，即:s= a1Q+a2Q2+ anQn式中 a1、a2an 待定系数。注：a1宜按均差表求得后，可相应地将公式(821-7)(821-1)、(821-2)、(821-3)中的 Q/s和公式（821-4）、（821-5）、（821-6）中的Q,以1/a1代换，分别进行计算。3当s/Q （或Ah2/Q）Q关系曲线呈直线时，可采用作图截距法求出a1后，按本条第二款代换,算。单孔稳定流抽水试验，当利用观测孔中的水位下降资料计算渗透系数时，s （或A h2）在s （或A h2）lgr关系曲线上能连成直线，可采用下

12、列公式: 1承压水完整孔：若观测孔中的值2 7T j'V/ 5 J 5 2（822-1 ）2潜水完整孔:（822-2）式中 si、s2在slgr关系曲线的直线段上任意两点的纵坐标值（m）;，* I1 IA/在A h2Igr关系曲线的直线段上任意两点的纵坐标值（m2）;门、r2在s （或Ah2）lgr关系曲线上纵坐标为 si、s2 （或） 的两点至抽水孔的距离（m）。单孔非稳定流抽水试验，在没有补给的条件下，利用抽水孔或观测孔的水位下降资料计算 渗透系数时，可采用下列公式：1配线法：1）承压水完整孔：W (uMsS r2=4KM * 72）潜水完整孔:(L08Q_ W(“)Istw、If

13、4Kh(8 2 3-4)（8.2.2-1）、（822-2）或下列公式:(823-5)Q2n(Ahl - 婕)式中 W （u）井函数；S承压水含水层的释水系数;科一一潜水含水层的给水度。2直线法：啪（或心）当4 KMl 4 K<0.01时，可采用公式1）承压水完整孔:_QK =-47rM（” 、*1）水完整孔：（823-6）式中 s1、s2观测孔或抽水孔在slgt关系曲线的直线段上任意两点的纵坐标值（m）；观测孔或抽水孔在Ah2lgt关系曲线的直线段上任意两点的纵坐标值（m2）;t1、t2在s （或A h2）lgt关系曲线上纵坐标为 s1、s2 （或） 两点的相应时间（min）。8.2.4

14、单孔非稳定流抽水试验，在有越流补给（不考虑弱透水层水的释放）的条件下，利用slgt关系曲线上拐点处的斜率计算渗透系数时，可采用下式：2.3Q（8.2.4）式中r观测孔至抽水孔的距离（m）;B 越流参数；mislgt关系曲线上拐点处的斜率。1/2确定曲线的拐注：1拐点处的斜率，应根据抽水孔或观测孔中的稳定最大下降值的 点位置及拐点处的水位下降值，再通过拐点作切线计算得出。炉22越流参数，应根据1 ,从函数表中查出相应的r/B ,然后确定越流参数 B。8.2.5稳定流抽水试验或非稳定流抽水试验，下列公式：1停止抽水前，若动水位已稳定，可采用公式当利用水位恢复资料计算渗透系数时，可采用（824）计算

15、，式中的mi值应采用恢复水位的曲线上拐点的斜率。2停止抽水前，若动水位没有稳定，仍呈直线下降时，可采用下列公式:1）承压水完整孔：QMs1（825-1 ）2）潜水完整孔:1nli十- 上丁(825-2)式中 tk-tTs-h一抽水开始到停止的时间-抽水停止时算起的恢复时间 水位恢复时的剩余下降值-水位恢复时的潜水含水层厚度(min)；(min)；(m)；(m)。注：1当利用观测孔资料时，应符合（或上g42n <0.01时的要求。2如恢复水位曲线直线段的延长线不通过原点时，应分析其原因，必要时应进行修正。利用同位素示踪测井资料计算渗透系数时，可采用下列公式：(826-1 )（826-2）式

16、中 Vf 测点的渗透速度（m/d）;I 测试孔附近的地下水水力坡度；r测试孔滤水管内半径（m）;r0探头半径（m）;t示踪剂浓度从 N0变化到Nt所需的时间（d）;N0 同位素在孔中的初始计数率；Nt同位素t时的计数率；Nb 放射性本底计数率；a流场畸变校正系数。0. 366QK=方法（三）insR=10SM在单孔抽水试验中，由于没有观测孔，只能根据抽水试验未稳定前的水位，做出BI深-半对数 时间图，以图解法来求渗透系数或根据水位恢复数据，以图解法来求岩石渗透系数.像这种联解方程，想用数学推导法来求解，是非常困难的.涉及募函数和指数涵数.如一矿山的抽水试验，涌水量为Q=1053吨/天，含水层厚

17、度为 m=241.3米，降深s=9.40米，抽 水管径r=0.084米.经过化简和代入后为：k-0.085lgk=1.41113可以用逼进法，在excel里计算.如k=1时，左边的式子，其得数是小于1的，显然不符合方程.如k=3时,左边的式子，其得数是大于2的,显然也是不符合方程.如k=2时,左边的式子，其得数是介于1-2之间的，这样就界定了 k值的大致范围然后再分别计算k=1.1 左边的式子k-0.085lgk=1.0965k=1.2,k-0.085lgk=1.19327k=1.3,k-0.085lgk=1.2903k=1.4,k-0.085lgk=1.38758k=1.5,k-0.085l

18、gk=1.4850显然 k 值小于 1.5,大于 1.4.然后再这个区间继续计算k-0.085k 的值 ,使之趋近于 1.41113.在 excel 里 ,这种计算非常快捷.很快就可以得出k=1.4244R=112.12水文地质参数确定方法确定水文地质参数的方法一般分为经验数据法、 经验公式法、 室内试验法和野外试验法四种。供水水文地质勘察中主要采用野外试验法，因为野外试验法求得的参数精确度较高。经验数据法根据长期的经验积累的数据， 列成表格供需要时选用。 渗透系数、 压力传导系数、释水系数、越流系数、弥散系数、降水入渗系数、给水度和影响半径等都有经验数据表可查。在评估地下水资源时，水文地质参

19、数常采用经验数据。经验公式法考虑到某些基本规律列出的公式， 并加上经验的修正。 渗透系数、 给水度等都可按经验公式计算，其值比选用经验数据的精确度要高。室内试验法在野外采取试件， 利用试验室的仪器和设备求取参数。渗透系数、给水度、降水入渗系数等水文地质参数，可通过室内试验法求得。野外试验法利用野外抽水试验取得有关数据， 再代入公式计算水文地质参数。 其计算公式分稳定流公式和非稳定流公式。 计算时根据含水层的状态（潜水或承压水）、 井的完整性（完整井或非完整井） 、边界条件（傍河或其他边界）、抽水孔状态（单孔抽水或带观测孔抽水）等条件选择。渗透系数、导水系数、压力传导系数、释水系数、越流系数、给

20、水度和影响半径等，都可用野外抽水试验法求得较精确的数据。 野外确定降水入渗系数还可采用地下水均衡试验场的实测数据，一般精度较高。8.2 渗透系数8.2.1 单孔稳定流抽水试验， 当利用抽水孔的水位下降资料计算渗透系数时， 可采用下列公式：1当Qs （或A h2）关系曲线呈直线时，1 ）承压水完整孔：2(8.2.1-1)27rtM f2）承压水非完整孔:当 M>150r , l/M>0.1 时:(凡x r ITtr(8.2.1-2)或当过滤器位于含水层的顶部或底部时:ln(l + 0.2 力(821-3)3）潜水完整孔:(821-4)4）潜水非完整孔:当>150r, l>

21、0.1 时:Q 八区卜h tt(H2 - A2)( 0 尸(821-5)或当过滤器位于含水层的顶部或底部时:Q jr(-11- M)(821-6)式中 K渗透系数（m/d）;Q出水量（m3/d）;s水位下降值（m）;M 承压水含水层的厚度（m）;H 自然情况下潜水含水层的厚度（m）;h潜水含水层在自然情况下和抽水试验时的厚度的平均值（m）;h潜水含水层在抽水试验时的厚度（m）;l过滤器的长度（m）;r抽水孔过滤器的半径 （m）;R影响半径（m）。2当Qs （或Ah2）关系曲线呈曲线时，可采用插值法得出Qs代数多项式，即:s= aiQ+a2Q2+anQn（821-7）式中ai、a2an待定系数。

22、注：ai宜按均差表求得后，可相应地将公式（821-1）、 （821-2）、（821-3）中的QQ/s和公式（821-4）、（821-5）、（821-6）中的 “一 " 以1/a1代换，分别进行计省一 歼。3当s/Q （或Ah2/Q）Q关系曲线呈直线时，可采用作图截距法求出a1后，按本条第二款代换，并计算。8.2.2 单孔稳定流抽水试验，当利用观测孔中的水位下降资料计算渗透系数时，若观测孔中的值s （或Ah2）在s （或A h2）lgr关系曲线上能连成直线，可采用下列公式：9In 二2tzM（5 1 S2） r I7T（人彳-Aftp厂1承压水完整孔:（822-1 ）2潜水完整孔：（8

23、.2.2-2）式中 &、s在slgr关系曲线的直线段上任意两点的纵坐标值（m）;N屁、-fli在Ah2lgr关系曲线的直线段上任意两点的纵坐标值（m2）;门、r2在s （或A h2）lgr关系曲线上纵坐标为 s1、s2 （或） 的两点至抽水孔的距离（m）。8.2.3 单孔非稳定流抽水试验，在没有补给的条件下，利用抽水孔或观测孔的水位下降资料 计算渗透系数时，可采用下列公式：1配线法:增2.3.1)阳 23-2)1）承压水完整孔:K - W yit )MsS r2"4 KM t2)潜水完整孔：Kii4 KH 7(823-3)<8.2.3-4)式中 W （u）井函数;S承压水含水层的释水系数;科一一潜水含水层的给水度。2直线法：当上£一 （或匕 <0.01时，可采用公式（822-1）4 KM八 4KhJ（822-2）或下列公式:1）承压水完整孔:兽1447rM(52 ' 1(823-5)水完整孔:(823-6)式中S1、S2(m2)；K =-Q In 丝2n(Ahl - 舒)/观测孔或抽水孔在Slgt关系曲线的直线段上任意两点的纵坐标值（m）;观测孔或抽水孔在A mlgt关系曲线的直线段上任意两点的纵坐标值tl、t2在s （或A h2）lgt关系曲线上纵坐标为 Si、S2 （或）两点的相应时间(min)