承压水的井损问题
6.1.2.1 承压水井损的计算公式
承压水井损的计算采用如下公式:
供水水文地质计算
式中:
可根据多落程单井抽水试验资料用相关分析法求出a1、a2,然后代入(6-1-2)式计算出不同涌水量Qi时的水头损失值
根据多落程单井抽水试验资料用相关分析法求a1、a2的计算程序如下。
6.1.2.2 承压水井损的计算程序
c zuocheng.for左乘转置矩阵法建立多元随机模型,然后求承压井损(样本容量n,自变量个数m,待定系数m1个),样本容量n,
c 自变量个数m,待定系数m1可在第1句任意给定,后边语句皆以n,m,m1出现;程序中还输出了DPP绘图的数据文件sq-q,可在GRAPHER中绘曲线并找出直线方程
parameter(n=3,m=1,m1=2)
dimension s1(n),x(n,m),y(n),b(m1),ab(m1,m1),yj(n)
c 从’qs’数据文件读入x(n,m):m个自变量的n组观测值,y(n):因变量的n组观测值
open(1,file=’qs’,status=’old’)
read(1,*)((x(i,j),j=1,m),s1(i),i=1,n)
close(1)
do 110 i=1,n
do 100 j=1,m
100 x(i,j)=x(i,j)*60.0*60.0*24.0/1000.0
y(i)=s1(i)/x(i,1)
110 continue
c 上面的y(i)即井损公式中的s(i)/q(i),x(i,1)即井损公式中的q(i)
open(1,file=’sq-q’,status=’new’)
do 140 i=1,n
140 write(1,150)x(i,1),y(i),
close(1)
150 format(2f15.8)
c 006-008:x阵的第1列全为零
do 40 k=1,n
x(k,0)=1.0
40 continue
c 009-014:因变量左乘x阵的转置矩阵得线性代数方程组的常数阵
do 60 i=1,m1
b(i)=0
do 50 k=1,n
b(i)=b(i)+x(k,i-1)*y(k)
50 continue
60 continue
c 015-022:自变量阵左乘x阵的转置矩阵得线性代数方程组的系数阵
do 70 i=1,m1
do 71 j=1,m1
ab(i,j)=0
do 72 k=1,n
ab(i,j)=ab(i,j)+x(k,i-1)*x(k,j-1)
72 continue
71 continue
70 continue
c 调用gas2全主元高斯消去法子程序解线性代数方程组得随机模型的待定系数b(i),i=1.m1 call gas2(m1,ab,b,0.0000000001,kwji)
yp=0
c 025-027:求因变量实测值的均值
do 80 i=1,n
yp=yp+y(i)/n
yj(i)=0
c 028-030:用新建立的多元随机模型求因变量的回归值
do 75 j=1,m1
yj(i)=yj(i)+b(j)*x(i,j-1)
75 continue
80 continue
huigui=0
piancha=0
c 035:回归平方和 036:偏差平方和
do 90 i=1,n
huigui=huigui+(yj(i)-y(i))*(yj(i)-y(i))
piancha=piancha+(y(i)-yp)*(y(i)-yp)
90 continue
huigui=huigui/(n-2)
piancha=piancha/(n-1)
r=sqrt(1-huigui/piancha)
write(*,10)(i,b(i),i=1,m1)
10 format(’b(’,i2,’)=’,f16.8)
write(*,20)r
20 format(’R=’,f6.4)
do 120 i=1,n
ds=b(2)*x(i,1)**2
ss=ds/s1(i)
120 write(*,130)i,ds,ss
130 format(i3,’ds’,f12.8,’ds%=’,f12.4)
stop
end
c 全主元高斯消去法子程序(见参考文献[25]P127,见例题)
subroutine gas2(m,ab,b,ep,kwji)
dimension ab(m,m),b(m),l(3)
do 20 i=1,m
20 l(i)=i
do 10 k=1,m
c=0
do 11 i=k,m
do 11 j=k,m
q=ab(i,j)
if(abs(q).le.abs(c))goto 11
c=q
io=i
jo=j
11 continue
if(abs(c).le.ep)return
if(jo.eq.k)goto 19
do 21 i=1,m
t=ab(i,jo)
ab(i,jo)=ab(i,k)
21 ab(i,k)=t
j=l(k)
l(k)=l(jo)
l(jo)=j
19 if(io.eq.k)goto 13
do 12 j=k,m
t=ab(io,j)
ab(io,j)=ab(k,j)
12 ab(k,j)=t
t=b(io)
b(io)=b(k)
b(k)=t
13 c=1/c
b(k)=b(k)*c
if(k.eq.m)goto 10
k1=k+1
do 23 j=k1,m
23 ab(k,j)=ab(k,j)*c
do 15 i=k1,m
do 14 j=k1,m
14 ab(i,j)=ab(i,j)-ab(i,k)*ab(k,j)
15 b(i)=b(i)-b(k)*ab(i,k)
10 continue
do 17 ii=2,m
i=m+1-ii
i1=i+1
do 17 j=i1,m
17 b(i)=b(i)-ab(i,j)*b(j)
do 24 k=1,m
lk=l(k)
24 ab(1,lk)=b(k)
do 25 k=1,m
25 b(k)=ab(1,k)
return
end
6.1.2.3 承压水井损的计算实例
以KC1单井抽水试验三落程为例,其三落程si分别为2.16m、4.40m、6.69m,所对应的涌水量分别为915.84m3/d、1581.12m3/d、1944.00m3/d,经用该程序计算,相关方程为:
si/Qi=0.00137763+0.00000100Qi。
经用该相关方程计算,KC1单井抽水试验三落程所对应的水头损失值
中文英文对照
井损well loss 井位well location 静水压力 hydrostatic pressure 局域( 部) 流动系统 local flow system 决策变量 decision variable 均质homogeneity K 开采...再生水 regenerated water 张力tension 张量tensor 蒸发evaporation 蒸发岩( 盐) evaporite 蒸发蒸腾 evapotranspiration 正定positive definiteness 正演问题 forward...
井附近的水位降深
在以后各节建立的公式中,都不考虑井损和有效井半径问题。在本章最后一节中,将介绍有关井损和有效井半径的确定方法。图3—2 承压含水层中的水位降深和有效井径(据Jacob)(a)—裸井;(b)—下过滤器的井;(c)—填砾的井 在以后几节中,除了特别提到的以外,一般都采用了以下假设:(1)...
井附近的水位降深
2)有过滤器的井(图3-2(b)),一般情况下过滤器的半径即井半径(rw),由于水流通过过滤孔眼以及进水泵时,产生水头损失(井损),使井壁内外水位不一致,通常井管外的水头高于井管内的水头。3)在过滤器外面充填石料时,如图3-2(c)所示,此时的半径采用有效半径,有效半径等于井轴到井管外壁...
Dupuit公式的应用
前面导出的Dupuit公式,可以解决下列两类问题:(1)求含水层参数。这时可将Dupuit公式写成便于求参数的形式。对于承压水井有:地下水动力学(第二版)对于潜水井有:地下水动力学(第二版)对其它公式,也可作类似变换,不再一一列举了。将抽水试验趋近稳定时测得的流量Q及抽水井或观测孔的水位降深s...
思念贾福海先生
报到以后,贾福海主任工程师要我担任坝址工程地质队的第一地质队长,全面负责坝址的工程地质填图和钻探编录、压水试验等工作,负责初步设计阶段坝址工程地质勘察...以贾总为首的专家组,肯定了我推导的“消除井损计算渗透系数的公式”(日后,该公式纳入了供水水文地质勘察国家标准、专著和水文地质手册,并获得了地质矿产部...
戊弯健儿: 必须有
临武县13053252475: 地下水动力学中的井损和渗出面有什么区别? - ?
戊弯健儿: 井损:井内水位与含水层中靠近井壁的水位不相同.井损的成因:1)井管外的水经过过滤器的孔眼进入井内的损失;2)井内的水流紊动损失;3)井管内的水向上运动的沿程损失等.渗出面:当潜水流入井中时,井壁水位hs高于井中水位hw,称其为渗出面(水跃).
临武县13053252475: 为什么井里的水不会满? - ?
戊弯健儿: 井水都是地下水,地下水分为两种,一种是潜水,一种是承压水,一般情况下的井都地在高于潜水位的地方挖的,所以水不会流到地面上来. 在承压区挖的井才会自动流到地面上来,那种井叫做承压井.
临武县13053252475: 自流井,承压井,潜水井的区别是什么? - ?
戊弯健儿: 自流井:井所处的位置低于承压水位时,承压水将涌出地表,形成自流井.自流井只针对承压水位而言. 承压井:透过隔水层,钻到承压水中的井叫承压井,承压井中的水受到压力的影响...
临武县13053252475: 向斜盆地承压水的成因 - ?
戊弯健儿: 承压水(confined groundwater ),充满两个隔水层之间的含水层中的地下水.承压水由于顶部有隔水层,它的补给区小于分布区,动态变化不大,不容易受污染.它承受静水压力.在适宜的地形条件下,当钻孔打到含水层时,水便喷出地表,形成...
临武县13053252475: 地下承压水是怎样受冲击形成自流井的? - ?
戊弯健儿: 水源于大气降水,湖底有粘土组成的隔水层;而井水却源于深部承压含水层,此含水层延伸较远,且补给区位置较高,当地下水以径流形式运移时,逐渐转为承压,形成地下水承压区,含水层被切穿时,承压水涌出地面就形成了上升泉,挖井含水层被切穿时,承压水涌出地面就形成了上升泉,且水位高出地面,就是自流井.
临武县13053252475: 地下水 传我们大都饮承压水(因潜水层水不干净),请问掘井的时候先出来脏水而后干净水最终不还是脏水? - ?
戊弯健儿: 打井的时候,在不想要的水层,会在外面用胶泥填一圈,不让哪一层水(也就是你说的脏水)渗到井里来,所以最终的井水是干净的,符合我们的生活用水标准.
临武县13053252475: 承压水出露地表形成 自流井 和 泉,二者有何不同?一个是人为打出的井,一个是自己找缝喷出的? - ?
戊弯健儿: 地下水分为浅层地下水和承压水,浅层地下水埋藏浅,可以抽出来,需要人为影响.承压水埋藏深,水受岩层压力大,或水位高出地表.当承压水沿着地层裂缝喷涌而出时,就形成了喷泉.这是天然泉,不受人为影响的.自流井就是人工打的井,不用人抽,水可自流.承压井分布最广最有名的要数大自流盆地(澳大利亚大盆地)
临武县13053252475: 潜水与承压水的根本区别在于? - ?
戊弯健儿: 地下水有两种不同的埋藏类型,即埋藏在第一个稳定隔水层之上的潜水和埋藏在上下两个稳定隔水层之间的承压水.潜水和承压水除了埋藏条件不同外,还有一定的区别. 潜水的补给主要是当地的大气降水和部分河湖水.承压水则是依靠大气降...
临武县13053252475: 承压水的包含水层 - ?
戊弯健儿: 充满于两个隔水层之间的含水层中的地下水.典型的承压含水层可分为补给区、承压区及排泄区三部分.补给区含水层裸露,具有自由水面,实际上分布着潜水,可接受外界补给. 排泄区承压水通过上升泉和向浅部含水层越流等方式排泄. 承压区含水层所处的位置高程较小,充满着承受压力的水.当井或钻孔揭穿隔水顶板时,井孔中的水上涌到含水层顶面以上一定高度才停止下来.静止水面的高程即承压含水层的测压水位.测压水位高出含水层顶面的距离为承压水的压力水头.在一定的地形、地质条件下,测压水位高出地面,井孔喷发自流水,成为自流井.