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

抽水试验确定渗透系数的方法及步骤1.抽水试验资料整理
试验期间，对原始资料和表格应及时进行整理。

试验结束后，应进行资料分析、整理，提交抽水试验报告。

单孔抽水试验应提交抽水试验综合成果表，其内容包括：水位和流量过程曲线、水位和流量关系曲线、水位和时间（单对数及双对数）关系曲线、恢复水位与时间关系曲线、抽水成果、水质化验成果、水文地质计算成果、施工技术柱状图、钻孔平面位置图等。

并利用单孔抽水试验资料编绘导水系数分区图。

多孔抽水试验尚应提交抽水试验地下水水位下降漏斗平面图、剖面图。

群孔干扰抽水试验和试验性开采抽水试验还应提交抽水孔和观测孔平面位置图（以水文地质图为底图）、勘察区初始水位等水位线图、水位下降漏斗发展趋势图
（编制等水位线图系列）、水位下降漏斗剖面图、水位恢复后的等水位线图、观测孔的S—t、S—lg t曲线［注卜各抽水孔单孔流量和孔组总流量过程曲线等。

注意：（1）要消除区域水位下降值；（2）在基
岩地区要消除固体潮的影响；3）傍河抽水要消除河水位变化对抽水孔水位变化的影响。

多孔抽水试验、群孔干扰抽水试验和试验性开采抽水试验均应编写试验小结，其内容包括：试验目的、要求、方法、获得的主要成果及其质量评述和结论。

2.稳定流抽水试验求参方法
求参方法可以采用Dupuit公式法和Thiem公式法。

⑴只有抽水孔观测资料时的Dupuit公式承压完整井：
4 2 In P
Ms. AT q
井:
潜水完整f
—ln^
7V（H2-h2）r v
2 2斗届7
式中K- -一含水层渗透系数（m/d）;
Q-抽水井流量（m3/d）;
sw抽水井中水位降深（m）;
M ----- 承压含水层厚度（m）;
R――影响半径（m）;
H ――潜水含水层厚度（m）;
h——潜水含水层抽水后的厚度（m）;
rw ——抽水井半径（m）。

（2）当有抽水井和观测孔的观测资料时的
Dupuit 或Thiem 公式
當压完整
井；松一饨-。

In'1
Thiem公
h -h = ° In 9
式；
潜水完荃井
£岸-比=In—
Thiem公式
陽圧三° hi乜
匕
式中hw ——抽水井中水柱高度（m）;
hl、h2——与抽水井距离为r1和r2处观测孔（井）中水柱高度（m），分别等于初始水位H0
与井中水位降深s之差，h1= H0 -1; h2= H0 -s2o
其余符号意义同前当前水井中的降深较大时，可采用修正降深
正降深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 为抽水井半径或观测孔至抽水井的距离)：
1>]
T
—
M
_ AT[t]
£=——
2 r I
尸[-]
«=—[-]
A[t]w
以上为降深时间法(s-t)。

也可以米用降
深---时间距离法(s-t/r2 )、降深---距离法(s-r) 进行参数计算。

(2) Jacob直线图解法
当抽水试验时间较长，u= r2/(4at)<0.01时，在半
对数坐标纸上抽水试验数据曲线s-t为一直线
(延长后交时间轴于to，此时s=0.00m),在直
线段上任取两点t1、si、t2、s2,则有
s2 ~ s i
2.257^
2.2510
(3) Hantush拐点半对数法
对半承压完整井的非稳定流抽水试验(存在越流量，K /b为越流系数)，当抽水试验时间较长，u= r2/(4at)<0.1时，在半对数坐标纸上抽水试验数据曲线s-t，外推确定最大水位降深Smax，在s-lgt线上确定拐点Si = Smax/2，拐点处的斜率mi及时间ti，则有
叫 B B B
一(H叫
叫
Br
K l T
(4) 水位恢复法
当抽水试验水位恢复时间较长，u= r2/(4at)<0.01
时，在半对数坐标纸上绘制停抽后水位恢复数据曲线s-t，在直线段上任取两点t1 , si, t2, s2,
则有
T二°-1835ln h
巧—f 1
2.25^
J
a
(5) 水位恢复的直线斜率法
当抽水试验水位恢复时间较长，u= r2/(4at)<0.1 时，在半对数坐标纸上绘制停抽后水位恢复数据曲线s-t，直线段的斜率为B，则有
Aft S
3.2潜水非稳定流抽水试验求参方法
潜水参数计算可采用仿泰斯公式法、Boulton 法和Numan法。

(1)仿泰斯公式法
式中HO 、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)潜水完整井考虑迟后疏干的 Boult on 公式 ■ 1-厂 匚曲、+------ : ------ sfm 、 J. (F 〕 1%
B- ■ ■■ —
■1
磊5
可根据抽水早期、中期、晚期的观测资料，采用捕水中朗: 1/
抽水早Q J 1
4皿 4R
抽水晩
4 at 477
相应的方法计算参数。

（3） Numan 法
对于潜水含水层完整井非稳定流抽水试验，也可 以采用Numan 模型求参，具体求参过程可参阅
《地下水动力学》等教科书。

4.参数计算新技术新方法的应用
采用AQUIFERYTEST 软件（图1）、数值模拟 法（可采用 GMS 、MODFLOW 、FEFLOW 等 软件）以及肖长来教授提出的全称曲线拟合法
（图2）等一些新的软件、方法确定水文地质参 数，效果非常好。

井第 1 取抽水（Hant L shJacob}
Con ductivity:9.38E-2 m/d
图1 AQUIFER YTEST 软件求参图示
+2+1+0-1-2 E E- E E ^rr- d
— HI- HI di rTb ‘ _________
: I
■—
一
•
，” I E
怕农H 间t 〔EjjP
图2全称曲线拟合法求参图示
5.参数计算结果的验证 上述参数计算结果的精度如何，取决于试验场地 水文地质条件的概化，也取决于观测数据的精 度。

对于所求得的参数，应将其代入相应的公式， 通过对比计算降深与实测降深的差值，分析所求 参数的精度及其可靠性和代表性，最终确定抽水 试验场地的有代表性意义的参数值。

方法（二）
单孔稳定流抽水试验，当利用抽水孔的水位下 降资料计算渗透系数时，可采用下列公式：
1当Q 〜s （或△ h2）关系曲线呈直线时，
1）承压水完整孔：
(8.2.1-1)
Q 卜严
2TTS M 厂
2）承压水非完整孔：当M>150r , l/M>0.1 时:
Q
2nsM
(
涉+气召皿'
v r I nr i
(821-2)
或当过滤器位于含水层的顶部或底部时:
与厶(1 +0.2弓)]
(8.2.1-3)
3）潜水完整孔:
(821-4)
4) 潜水非完整孔:
当>150r, l > 0.1 时:
Q
7T (FI- h2)
(8.2.1-5)
ln R
r• In
1 ・12h \
Trr) /心
(821-6)
式中K -------- 渗透系数(m/d);
Q——出水量(m3/d);
s――水位下降值(m);
M ----- 承压水含水层的厚度(m);
H ――自然情况下潜水含水层的厚度(m )；h——潜水含水层在自然情况下和抽水试验时的厚度的平均值(m)；
h――潜水含水层在抽水试验时的厚度(m)；
I 过滤器的长度(m)；
r――抽水孔过滤器的半径(m)；
R ――影响半径(m )。

2当Q〜s (或△ h2)关系曲线呈曲线时，可采用插值法得出Q〜s代数多项式，即：
s = a1Q+a2Q2+ ........... anQn (8.2.1-7)
式中al、a2 ........... an ----- 待定系数。

注：al宜按均差表求得后，可相应地将公式(821-1)、(821-2)、(8.2.1-3)中的Q/s 和公
q
式(821-4)、(821-5)、(8.2.1-6)中的厂以
1/a1代换，分别进行计算。

3当s/Q （或△ h2/Q）〜Q关系曲线呈直线时，可采用作图截距法求出al后，按本条第二款代换，并计算单孔稳定流抽水试验，当利用观测孔中的水位下降资料计算渗透系数时，若观测孔中的值s（或△ h2）在s （或△ h2）〜Igr关系曲线上能连成直线，可采用下列公式：
1 承压水完整孔：
K = 说
（8.2.2-1）
2潜水完整孔：
（822-2）
式中si、s2——在s〜Igr关系曲线的直线段上
任意两点的纵坐标值（m）;
——在厶h2〜Igr关系曲线的直线段
r1、r2 -------- 在s （或△ h2）〜Igr关系曲线上
纵坐标为si、s2 （或）的两点至抽水孔的距离
(m )。

单孔非稳定流抽水试验，在没有补给的条件下， 利用抽水孔或观测孔的水位下降资料计算渗透 系数时，可采用下列公式:
1配线法： 1)承压水完整孔:
2)潜水完整孔:
式中 W (u)——井函数；
S ——承压水含水层的释水系数; 卩 ——潜水
含水层的给水度。

2直线法:
-°川乞—
K = ----------:— W («)
AA 2
“ F 2
H - --------- ■—
AKH '
U ~ ------------- T
4Kh
4 KM
总仏忒）＜0.01时，可采用公式
（822-1）、（8.2.2-2）或下列公式：
1）承压水完整孔:
r _ M弋r ■ln
（823-5）
水完整孔:
K ------ ---- * 亠
2*△卮—/1
（823-6）
igt关式中si、s2——观测孔或抽水孔在s
系曲线的直线段上任意两点的纵坐标值（m）; W 时——观测孔或抽水孔在△ h2〜Igt关系曲线的直线段上任意两点的纵坐标值（m2）;
t1、t2 在s （或△ h2）〜Igt 关系曲
两点的相应时间
s 〜Igt 关系曲线上拐点处的斜率计算渗透系数 时，可米
用下式：
K = 2.3Q
4K * M • nij * e r/tl
（824）
式中r ――观测孔至抽水孔的距离
（m ）；
B ---- 越流参数；
mi ----- s 〜Igt 关系曲线上拐点处的斜率。

注：1拐点处的斜率，应根据抽水孔或观测 孔中的稳定最大下降值的1/2确定曲线的拐点位 置及拐点处的水位下降值，再通过拐点作切线计 算得出。

e T
^K^=2.3 —
2越流参数，应根据
"匕，从
函数表中查出相应的r/B ，然后确定越流参数B 。

8.2.5稳定流抽水试验或非稳定流抽水试验，当 利用水位恢复
线上纵坐标为si 、s2 （或） （min ）。

8.2.4单孔非稳定流抽水试验,
（不考虑弱透水层水的释放）
在有越流补给 的条件下，利用
资料计算渗透系数时，可采用下列公式：
1停止抽水前，若动水位已稳定，可采用公式（8.2.4）计算，式中的mi值应采用恢复水位的
曲线上拐点的斜率。

2停止抽水前，若动水位没有稳定，仍呈直线 下降时，
可采用下列公式：
1）承压水完整孔：
K = Q ln （l + 判
（8.2.5-1）
2）潜水完整孔:
K =
&
2?r （H 2 一 h 2y
（8.2.5-2）
式中tk --------- 抽水开始到停止的时间
（min ）；
tT ——抽水停止时算起的恢复时间 （min ）； s ――水位恢复时的剩余下降值
（m ）；
h ------ 水位恢复时的潜水含水层厚度
（m ）。

注：1
<0.01时的要求。

2如恢复水位曲线直线段的延长线不通过
lull + k
当利用观测孔资料时，应符合当
仏二 i
原点时，应分析其原因，必要时应进行修正
利用同位素示踪测井资料计算渗透系数时，可采用下列公式：
（826-1）
式中Vf――测点的渗透速度（m/d）；
I ――测试孔附近的地下水水力坡度；r ――测试孔滤水管内半径（m）；
r0 ――探头半径（m）；t 示踪剂浓度从NO变化到Nt所需的时间
（d）；
NO——同位素在孔中的初始计数率；Nt 同位素t时的计数率；
Nb ――放射性本底计数率；a――流场畸变校正系数。

方法（三）
().；3&6Q
111S
R-los,nr
在单孔抽水试验中，由于没有观测孔，只能根据抽水试验未稳定前的水位，做出降深-半对数时间图,以图解法来求渗透系数或根据水位恢复数据,以图解法来求岩石渗透系数.
像这种联解方程，想用数学推导法来求解，是非常困难的•涉及幂函数和指数涵数.
如一矿山的抽水试验，涌水量为Q=1053吨/天. 含水层厚度为m=241.3米，降深s=9.40米,抽水管径r=0.084 米.
经过化简和代入后为：
k-0.085lgk=1.41113
可以用逼进法,在exceI里计算.
如k=1时，左边的式子，其得数是小于1的,显然不
符合方程.
如k=3时，左边的式子，其得数是大于2的,显然也
是不符合方程•
如k=2时，左边的式子，其得数是介于1--2之间的, 这样就界定了k值的大致范围
然后再分别计算
k=1.1 左边的式子k-0.085lgk=1.0965
k=1.2,k-0.085lgk=1.19327
k=1.3,k-0.085lgk=1.2903
k=1.4,k-0.085lgk=1.38758
k=1.5,k-0.085lgk=1.4850
显然k值小于1.5,大于1.4.
然后再这个区间继续计算k-0.085k的值,使之趋
近于1.41113.
在excel里，这种计算非常快捷.很快就可以得出
k=1.4244
R=112.12
水文地质参数确定方法
确定水文地质参数的方法一般分为经验数据法、
经验公式法、室内试验法和野外试验法四种。

供水水文地质勘察中主要采用野外试验法，因为野外试验法求得的参数精确度较高。

经验数据法根据长期的经验积累的数据，列成表格供
需要时选用。

渗透系数、压力传导系数、释水系数、越流系数、弥散系数、降水入渗系数、给水度和影响半径等都有经验数据表可查。

在评估地下水资源时，水文地质参数常采用经验数据。

经验公式法考虑到某些基本规律列出的公式，并加上经验的修正。

渗透系数、给水度等都可按经验公式计算，其值比选用经验数据的精确度要高。

室内试验法在野外采取试件，利用试验室的仪器和设备求取参数。

渗透系数、给水度、降水入渗系数等水文地质参数，可通过室内试验法求得。

野外试验法利用野外抽水试验取得有关数据，再代入公式计算水文地质参数。

其计算公式分稳定流公式和非稳定流公式。

计算时根据含水层的状态（潜水或承压水）、井的完整性（完整井或非完整井）、边界条件（傍河或其他边界）、抽水孔状态（单孔抽水或带观测孔抽水）等条件选择。

渗透系数、导水系数、压力传导系数、释水系数、越流系数、给水度和影响半径等，都可用野外抽水试验法求得较精确的数据。

野外确定降水入渗系数还可采用地下水均衡试验场的实测数据，
一般精
度较高。

8.2渗透系数
8.2.1单孔稳定流抽水试验，当利用抽水孔的水位下降资料计算渗透系数时，可采用下列公式：
1当Q〜s （或△ h2）关系曲线呈直线时，
1）承压水完整孔：
瓦Q 1 R
k- 22山r
（8.2.1-1）
2）承压水非完整孔：
当M>150r , l/M>0.1 时：
k二毛於+XT
或当过滤器位于含水层的顶部或底部时:
K二為［详十忙吐屮普）］
（8.2.1-3）
3）潜水完整孔:
Q卞
7T}
(821-4)
4）潜水非完整孔:
当>150r, I > 0.1
K = 一（门+ 亙宀ln^）兀（H? —/八）'厂1兀厂'（8.2.1-5）
或当过滤器位于含水层的顶部或底部时:
式中K ------- 渗透系数（m/d）;
Q ――出水量（m3/d）;
s――水位下降值（m）;
M ----- 承压水含水层的厚度（m）;
H――自然情况下潜水含水层的厚度
（m）；
h——潜水含水层在自然情况下和抽水试验时的厚度的平均值（m）;
h――潜水含水层在抽水试验时的厚度
（m）；
I 过滤器的长度（m）;
r――抽水孔过滤器的半径（m）;
R――影响半径（m ）。

2当Q〜s （或△ h2）关系曲线呈曲线时，可采用插值法得出Q〜s代数多项式，即：
s=
a i Q+a2Q2+ a n Qn (821-7)
式中a i、a2 a n 待定系数。

注：a i宜按均差表求得后，可相应地将公式(8.2.1-1)、( 8.2.1-2)、( 8.2.1-3) 中的Q/s 和公式(8.2.1-4)、( 8.2.1-5)、( 8.2.1-6)中
Q |
的/亦以1/a i代换，分别进行计算。

3当s/Q (或△ h2/Q)〜Q关系曲线呈直线时，可采用作图截距法求出a i后，按本条第二款代换，并计算。

8.2.2单孔稳定流抽水试验，当利用观测孔中的水位下降资料计算渗透系数时，若观测孔中的值s (或△ h2)在s (或△ h2)〜Igr关系曲线上能连成直线，可采用下列公式：
1 承压水完整孔：
2潜水完整孔:
(822-2)
式中S i、S2——在s〜Igr关系曲线的直线段上任意两点的纵坐标值（m）;
uh歸------ 在厶h2〜Igr关系曲线的直线段上任
意两点的纵坐标值（m2）;
r i、r2 --------------- 在s （或△ h2）〜Igr关系曲线上纵坐标为s i、S2（或）的两点至抽水孔的距离（m）8.2.3单孔非稳定流抽水试验，在没有补给的条件下，利用抽水孔或观测孔的水位下降资料计算渗透系数时，可采用下列公式：
1配线法：
1）承压水完整孔:
2）潜水完整孔:
K 冒wq
!J r
it = ---- = * — 4KA $ 式中 W （u ）——井函数；
S ——承压水含水层的释水系数;
|1——潜水含水层的给水度。

2直线法：
当一-——<0.01时，可采用公式
4/vVf 八皿"
（8.2.2-1）、（ 822-2）或下列公式：
1）承压水完整孔：
J 1 il 4 KM 点23-2〉
4/<f/
（823-5） 水完整孔:
式中 s i 、s 2 观测孔或抽水孔在 s 〜Igt 关 系曲线的直线段上任意两点的纵坐标值
（m ）; M 、对 观测孔或抽水孔在 △ h 2〜Igt 关系 曲线的直线段上任意两点的纵坐标值（m 2）;
t i 、t 2 -------- 在s （或△ h 2）〜Igt 关系曲线
上纵坐标为S i 、s 2 （或）两点的相应时间
（min ）。

8.2.4单孔非稳定流抽水试验，在有越流补给
（不考虑弱透水层水的释放） 的条件下，利用 s 〜Igt 关系曲线上拐点处的斜率计算渗透系数 时，可米用下式： 4
TT ■ M • z ■
____ Q 4TT M - X])
・In
(823-6)
(824)
式中r――观测孔至抽水孔的距离(m);
B ---- 越流参数；
mi ----- s〜Igt关系曲线上拐点处的斜率。

注：1拐点处的斜率，应根据抽水孔或观测孔中的稳定最大下降值的1/2确定曲线的拐点位置及拐点处的水位下降值，再通过拐点作切线计算得出。

2越流参数，应根据丨、心，从函数表中查出相应的r/B，然后确定越流参数B。

8.2.5稳定流抽水试验或非稳定流抽水试验，当利用水位恢复资料计算渗透系数时，可采用下列公式：
1停止抽水前，若动水位已稳定，可采用公式
（824）计算，式中的m i值应采用恢复水位的
曲线上拐点的斜率
2停止抽水前，若动水位没有稳定，仍呈直线下降时，可采用下列公式：
1）承压水完整孔：
(825-1)
2）潜水完整孔:
(8.2.5-2)
式中t k――抽水开始到停止的时间
（min ）；
t T -------- 抽水停止时算起的恢复时间（min）；
s――水位恢复时的剩余下降值（m）；
h——水位恢复时的潜水含水层厚度
（m ）。

注：1当利用观测孔资料时，应符合当
点（读厶j V0.01时的要求。

八 4 to?
2如恢复水位曲线直线段的延长线不通过原点时，应分析其原因，必要时应进行修正。

826利用同位素示踪测井资料计算渗透系数
时，可采用下列公式:
(826-1)
亢（厂一一厂Q M
v
- 2总山、厂A；（826-2）
式中Vf ------ 测点的渗透速度（m/d）；
I ――测试孔附近的地下水水力坡度；r――测
试孔滤水管内半径（m）；
r o――探头半径 (m);
t 示踪剂浓度从N o变化到N t所需的时间(d); N o——同位素在孔中的初始计数率；
N t 同位素t时的计数率；
N b 放射性本底计数率；
a――流场畸变校正系数。