首页 理论教育非稳定流抽水试验资料整理与降深比值法分析

非稳定流抽水试验资料整理与降深比值法分析

【摘要】:利用非稳定抽水试验资料求水文地质参数,可用配线法、雅柯布半对数直线图解法、降深比值法等。导水系数:如含水层为潜水,则给水度:如含水层为承压水,则储水系数:最后取平均值:3.计算实例现有非稳定流抽水试验资料,潜水井的抽水量为2718m3/d。附表1-2试算过程表附表1-3抽水试验资料因此得:u2=0.0019。

利用非稳定抽水试验资料求水文地质参数,可用(泰斯标准曲线)配线法、雅柯布半对数直线图解法、降深比值法等。其中前两者基本上都是利用作图法来求解参数,其方法在相关教材中(地下水动力学或地下水文学)已作详细的介绍,不再赘述。这里仅介绍降深比值法。该方法简单易行,在资料严重不足时(仅需两个时刻的降深数据即可)就能估求参数。

1.降深方法的原理

在做非稳定流抽水试验时,同一观测孔中两个不同时刻t1、t2所观测的水位降深分别为S1、S2,相应时刻的积分参数为u1、u2

根据泰斯式可分别列出以下两式[1]:

上两式中:μ为给水度;μ*为储水系数。

设以上两式的比值为α,则

若规定S1<S2,则

α<1,W(u1)<W(u2)

由井函数的性质可推知

u1>u2

假设u1/u2=A(A>1),则

所以

如利用两个观测孔的资料,则

将式(附1-4)代入式(附1-3),得:

2.资料整理的步骤[1]

(1)根据选定的两个时刻t1、t2的抽水试验资料,利用式(附1-5)或式(附1-6)、式(附1-3)计算出相应的A值和α值。

(2)根据算出的A、α值,利用α-u2的关系曲线(或试算法)反查出u2值。再按式(附1-4)求出u1值。

(3)将u1、u2值分别代入泰斯公式[式(附1-1)、式(附1-2)],即可计算

参数。

导水系数: 

如含水层为潜水,则给水度:

如含水层为承压水,则储水系数:

最后取平均值:

3.计算实例

现有非稳定流抽水试验资料(附表1-3),潜水井的抽水量为2718m3/d。观测孔2个,编号分别为B1、B2。试根据该抽水资料求水文地质参数。

解:先用同一孔不同时刻的降深求解。

(1)选用B1孔资料计算[2]。现选用t1=60min,t2=120min。据附表1-3,与这两个时刻对应的B1孔的降深分别为:

S1=0.878m,S2=1.000m

由式(附1-5)得: A=t2/t1=120/60=2

由式(附1-3)得: α=S1/S2=0.878/1.00=0.878

根据A、α值,由试算表(附表1-2)求得u2值。

在试算中,先假定u2=0.01和0.001,查附表1-1分别得到W(Au2)和W(u2)的值,发现其比值α=W(Au2)/W(u2)随着u2的增大而变小。由于真实的α值(0.878)在这两个比值之间,可以断定,所要求的u2值应在0.01~0.001之间。为此,再次假定u2=0.005,经试算后再次断定,真实的u2应在0.001~0.005之间。如此经过反复试算发现,当u2=0.0019时,α=W(Au2)/W(u2)的比值很接近真实值(0.878)。

附表1-1 与u对应的W(u)的值

注 本表范围以外的数据可据文献[3],利用表达式W(u)=18.109u0.256-31.370u0.16+13.482近似求出。

附表1-2 试算过程表

附表1-3 抽水试验资料(摘自 李炳森,2005)

因此得:u2=0.0019。由式(附1-4)求得:

u1=Au2=2×0.0019=0.0038

由试算表(附表1-2)得:

W(u1)=4.9993;W(u2)=5.6906

由式(附1-4)求出:

u1=Au2=2×0.0019=0.0038

于是,由式(附1-8)可求出导水系数为[2]:

所以

T=(T1+T2)/2=(1231.56+1230.83)/2=1231.20(m2/d)

由式(附1-9)又求得给水度[2]:

μ=(μ12)/2=2.098×10-4

(2)若用B2孔的资料求参。选t1=60min,t2=120min;则相应的降深为:

S1=0.402m,S2=0.518m

仿照以上步骤分别求得[2]:

A=2,α=0.7761,u2=0.03,u1=0.06;

W(u1)=2.2953,W(u2)=2.9591

T1=1234.96m2/d,T2=1234.58m2/d

μ1=2.077×10-42=2.078×10-4

(3)如同时用B1和B2孔的资料求参。方法如下:

选B1孔:t1=1min,S1=0.201m

选B2孔:t2=80min,S2=0.454m

由式(附1-6b)得[2]:

由式(附1-3)得:

利用两个观测孔资料时的计算过程见附表1-4。

附表1-4 利用两个观测孔资料时的试算过程表

由附表1-4可见,当u2=0.044时,可使α=W(Au2)/W(u2)=0.442,故取u2=0.044。再由式(附1-4)求出:

u1=Au2=5×0.044=0.22

据试算表查得,当u2=0.044时,W(u1)=1.1454,W(u2)=2.5899。由式(附1-8)得导水系数T为:T1=1234.54m2/d,T2=1233.86m2/d。由式(附1-9)得:μ1=2.027×10-42=2.029×10-4

根据以上三次计算结果[2]:

导水系数T=1230.83~1235.58m2/d,平均值为1233.21m2/d。

给水度μ=2.027×10-4~2.098×10-4,平均值为2.063×10-4

此结果与用泰斯标准曲线配线法得计算的结果(T=1272.77m2/d,μ=2.063×10-4)基本一致[2]

练习题

从附表1-3中随机挑出另两个时段数据,按照文中介绍的步骤和公式,再次计算该含水层的导水系数T和给水度μ。