抽水试验确定渗透系数的方法与步骤

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

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

3、水井流量(m3/d);sw抽水井中水位降深(m);M承压含水层厚度(m);R影响半径(m);H潜水含水层厚度(m);h潜水含水层抽水后的厚度(m);rw抽水井半径(m).4) )当有抽水井和观测孔的观测资料时的Dupuit或Thiem公式承压完整井,-h=一过一山-“2md八Thiem.公式士潜水完整井,Thiem.公式士式中hw抽水井中水柱高度m;hl、h2与抽水井距离为ri和r2处观测孔井中水柱高度m,分别等于初始水位H0与井中水位降深s之差,h1=H0-si;h2=H0s2.其余符号意义同前.当前水井中的降深较大时,可采用修正降深.修正降深s'与实际降深s之间的关系为：s'

4、;=s-s2/2H.3.非稳定流抽水试验求参方法3.1承压水非稳定流抽水试验求参方法1Theis配线法在两相同刻度的双对数坐标纸上,分别绘制Theis标准曲线Wu-1/u和抽水试验数据曲线s-t,保持坐标轴平行,使两条曲线配合,得到配合点M的水位降深s、时间t、Theis井函数Mu及1/u的数值,按以下公式计算参数r为抽水井半径或观测孔至抽水井的距离丁二一4M动国TK=MATt*二r21rlir-以上为降深一一时间法s-t.也可以采用降深一时间距离法s-t/r2、降深一距离法s-r进行参数计算.2Jacob直线图解法当抽水试验时间较长,u=r2/4at0.01时,在半对数坐标纸上抽水试验数据曲

5、线s-t为一直线延长后交时间轴于t0,此时s=0.00m,在直线段上任取两点t1、si、t2、s2,那么有（3） Hantush拐点半对数法对半承压完整井的非稳定流抽水试验存在越流量,K'/b为越流系数,当抽水试验时间较长,u=r2/4at0.1时,在半对数坐标纸上抽水试验数据曲线s-t,外推确定最大水位降深Smax,在s-lgt线上确定拐点Si=Smax/2,拐点处的斜率mi及时间ti,那么有国4一国占等=的.5,求白了!G183R产叫盲二迎Brkt4水位恢复法当抽水试验水位恢复时间较长,u=r2/4at0.01时,在半对数坐标纸上绘制停抽后水位恢复数据曲线s-t,在直线段上任取两点

6、t1,s1,t2,s2,那么有5水位恢复的直线斜率法当抽水试验水位恢复时间较长,u=r2/4at0.1时,在半对数坐标纸上绘制停抽后水位恢复数据曲线s-t,直线段的斜率为B,那么有3.2潜水非稳定流抽水试验求参方法潜水参数计算可采用仿泰斯公式法、Boulton法和Numan法.1仿泰斯公式法式中H0、hw-初始水头及抽水后井中水头;Wu泰斯井函数；Q抽水井的流量m3/d;r到抽水井的距离m;t自抽水开始起算的时间d;T含水层的导水系数m2/d;T=Khm;hm潜水含水层的平均厚度m;K含水层的渗透系数m/d;A含水层的导压系数1/d;m潜水含水层的给水度.具体计算时可采用配线法、直线图解法、水

7、位恢复法等.2潜水完整井考虑迟后疏干的Boulton公式1M句1日11E隹1E-31E+41E+51E/1E+71E+21E*2抽水早期:即山.,口,4牙D0尸抽水中朗：±=上反?一抽水晚期：3=力.可根据抽水早期、中期、晚期的观测资料,采用相应的方法计算参数.3 3Numan法对于潜水含水层完整井非稳定流抽水试验,也可以采用Numan模型求参,具体求参过程可参阅?地下水动力学?等教科书.4 .参数计算新技术新方法的应用采用AQUIFERYTEST软件图1、数值模拟法可采用GMS、MODFLOW、FEFLOW等软件以及肖长来教授提出的全称曲线拟合法图2等一些新的软件、方法确定水文地质

8、参数,效果非常好.C1井第一次抽水(Hantcish-Jacob)Conductivity:9.38E-2m/d图1AQUIFERYTEST软件求参图示憎术叶Mtbniui图2全称曲线拟合法求参图示5 .参数计算结果的验证上述参数计算结果的精度如何,取决于试验场地水文地质条件的概化,也取决于观测数据的精度.对于所求得的参数,应将其代入相应的公式,通过比照计算降深与实测降深的差值,分析所求参数的精度及其可靠性和代表性,最终确定抽水试验场地的有代表性意义的参数值.方法二单孔稳定流抽水试验,当利用抽水孔的水位下降资料计算渗透系数时,可采用以下公式：1当Qs或Ah2关系曲线呈直线时,1承压水完整孔：8

9、.2.1-12承压水非完整孔：当M>150r,l/M>0.1时：RM-lyA2MyIn+7In27rxM(8.2.1-2)1丁I,或当过滤器位于含水层的顶部或底部时:QK二27rsM口存竽缶1+0.2第(821-3)3潜水完整孔:821-44潜水非完整孔：当>150r,l>0.1时:Q(.Rh-I.1.12Ar八(111+f*111)(821-5)水广力2)'厂I兀厂'或当过滤器位于含水层的顶部或底部时:821-6式中K渗透系数m/d;Q出水量m3/d;s水位下降值m;M承压水含水层的厚度m;H自然情况下潜水含水层的厚度m;h潜水含水层在自然情况下和抽水

10、试验时的厚度的平均值m;h潜水含水层在抽水试验时的厚度m;l过滤器的长度m;r抽水孔过滤器的半径m;R影响半径m.2当Qs或Ah2关系曲线呈曲线时,可采用插值法得出Qs代数多项式,即：s=a1Q+a2Q2+anQn821-7式中al、a2an待定系数.注：al宜按均差表求得后,可相应地将公式821-1、821-2、821-3中的Q/s以1/a1代换,分别进行计算.和公式821-4、821-5、821-6中的a1后,按本条第3当s/Q或Ah2/QQ关系曲线呈直线时,可采用作图截距法求出二款代换,并计算.单孔稳定流抽水试验,当利用观测孔中的水位下降资料计算渗透系数时,假设观测孔中的值s或Ah2在s

11、或Ah2lgr关系曲线上能连成直线,可采用以下公式：1承压水完整孔：822-12潜水完整孔：8.2.2-2式中s1、s2在slgr关系曲线的直线段上任意两点的纵坐标值m;在Ah2lgr关系曲线的直线段上任意两点的纵坐标值m2;r1、r2在s或Ah2lgr关系曲线上纵坐标为s1、s2或的两点至抽水孔的距离m.单孔非稳定流抽水试验,在没有补给的条件下,利用抽水孔或观测孔的水位下降资料计算渗透系数时,可采用以下公式：1配线法：1承压水完整孔：2潜水完整孔:Sil2出.£4KH/=呼叫)Ax(8233)军22式中Wu井函数；S承压水含水层的释水系数；潜水含水层的给水度.2直线法：当0.01时

12、,可采用公式822-1、822-2或以下公式:1承压水完整孔：823-5水完整孔：823-6式中si、s2观测孔或抽水孔在slgt关系曲线的直线段上任意两点的纵坐标值m；观测孔或抽水孔在Ah2lgt关系曲线的直线段上任意两点的纵坐标值m2;t1、t2在s或Ah2lgt关系曲线上纵坐标为si、s2或两点的相应时间min.8.2.4单孔非稳定流抽水试验,在有越流补给不考虑弱透水层水的释放的条件下,利用slgt关系曲线上拐点处的斜率计算渗透系数时,可采用下式：824式中r观测孔至抽水孔的距离m;B越流参数；mislgt关系曲线上拐点处的斜率.注：1拐点处的斜率,应根据抽水孔或观测孔中的稳定最大下降值

13、的1/2确定曲线的拐点位置及拐点处的水位下降值,再通过拐点作切线计算得出.2越流参数,应根据,从函数表中查出相应的r/B,然后确定越流参数Bo8.2.5稳定流抽水试验或非稳定流抽水试验,当利用水位恢复资料计算渗透系数时,可采用以下公式：1停止抽水前,假设动水位已稳定,可采用公式824计算,式中的mi值应采用恢复水位的曲线上拐点的斜率.2停止抽水前,假设动水位没有稳定,仍呈直线下降时,可采用以下公式：1承压水完整孔：825-12潜水完整孔：825-2式中tk抽水开始到停止的时间min;tT抽水停止时算起的恢复时间min;s水位恢复时的剩余下降值m;h水位恢复时的潜水含水层厚度m.注：1当利用观测

14、孔资料时,应符合当0.01时的要求.2如恢复水位曲线直线段的延长线不通过原点时,应分析其原因,必要时应进行修正.利用同位素示踪测井资料计算渗透系数时,可采用以下公式：826-1826-2式中Vf测点的渗透速度m/d;I测试孔附近的地下水水力坡度；r测试孔滤水管半径m;r0探头半径m;t示踪剂浓度从N0变化到Nt所需的时间d;N0同位素在孔中的初始计数率；Nt同位素t时的计数率；Nb放射性本底计数率；a流场畸变校正系数.方法三K=ms在单孔抽水试验中,由于没有观测孔,只能根据抽水试验未稳定前的水位,做出BI深-半对数时间图,以图解法来求渗透系数或根据水位恢复数据,以图解法来求岩石渗透系数.像这种

15、联解方程,想用数学推导法来求解,是非常困难的.涉及募函数和指数涵数如一矿山的抽水试验,涌水量为Q=1053吨/天,含水层厚度为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

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

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

18、系数、压力传导系数、释水系数、越流系数、给水度和影响半径等,都可用野外抽水试验法求得较精确的数据.野外确定降水入渗系数还可采用地下水均衡试验场的实测数据,一般精度较高.8.2渗透系数8.2.1单孔稳定流抽水试验,当利用抽水孔的水位下降资料计算渗透系数时,可采用以下公1当Qs或Ah2关系曲线呈直线时,1承压水完整孔：8.2.1-12承压水非完整孔：当M>150r,l/M>0.1时:厂QR-I,M2MK=-lln+7In(8.2.1-2)(821-3)2ksM/几厂或当过滤器位于含水层的顶部或底部时：父QnRM-xnoM1K=-In-H：In(1+0.22心MLfI'厂'

19、;3)潜水完整孔：821-44潜水非完整孔：当>150r,l>0.1时:或当过滤器位于含水层的顶部或底部时:7T1口,+红产山(1+0.2表)(821-6)式中K渗透系数m/d;Q出水量m/d；s水位下降值m;M承压水含水层的厚度m;H自然情况下潜水含水层的厚度m;h潜水含水层在自然情况下和抽水试验时的厚度的平均值m;h潜水含水层在抽水试验时的厚度m;l过滤器的长度m;r抽水孔过滤器的半径m;R影响半径m.2当Qs或Ah2关系曲线呈曲线时,可采用插值法得出Qs代数多项式,即:式中ai、a2ans=aiQ+a2Q2+anQn待定系数.(821-7)注：a1宜按均差表求得后,可相应地将

20、公式8.2.1-1、8.2.1-2、8.2.1-3中的Q/s和公式821-4、821-5、821-6中的以1/ai代换,分别进行计算.3当s/Q或Ah2/QQ关系曲线呈直线时,可采用作图截距法求出ai后,按本条第二款代换,并计算.8.2.2单孔稳定流抽水试验,当利用观测孔中的水位下降资料计算渗透系数时,假设观测孔中的值s或Ah2在s或Ah2lgr关系曲线上能连成直线,可采用以下公式：1 承压水完整孔：8.2.2-12 潜水完整孔：8.2.2-2式中si、s2在slgr关系曲线的直线段上任意两点的纵坐标值m;在Ah2lgr关系曲线的直线段上任意两点的纵坐标值m2;门、r2在s或Ah2lgr关系曲

21、线上纵坐标为si、s2或的两点至抽水孔的距离m.8.2.3 单孔非稳定流抽水试验,在没有补给的条件下,利用抽水孔或观测孔的水位下降资料计算渗透系数时,可采用以下公式：i配线法：i承压水完整孔：U.=2潜水完整孔:i-;-W或0.08Q、-I)垣.2324KHr4K式中Wu井函数;S承压水含水层的释水系数;it潜水含水层的给水度.2直线法:当0.01时,可采用公式822-1、822-2或以下公式:1承压水完整孔:(823-5)水完整孔:(823-6)式中&、S2观测孔或抽水孔在slgt关系曲线的直线段上任意两点的纵坐标值(m)；观测孔或抽水孔在Ah2lgt关系曲线的直线段上任意两点的纵坐标值m2;tl、t2在s或Ah2