2024年5月6日发(作者：查夏萱)
南京
第2期
2000年6月
水利科学研究院
No.2
Jun.2000
水利水运科学研究
JOURNALOFNANJINGHYDRAULICRESEARCHINSTITUTE
混凝土面板坝渗流分析方法初探
盛金保,李 雷
(南京水利科学研究院,水利部大坝安全管理中心,江苏南京 210029)
a
摘要:
对混凝土面板坝的渗流分析方法,重点是面板与高挡墙底板间水平接缝止水破坏后的渗流分析方法进行了
探讨,可供面板坝设计参考.
关键词:
混凝土面板坝;渗流;分析方法
中图分类号:TV641.43 文献标识码:A 文章编号:1001-3962(2000)02-0039-05
混凝土面板堆石坝是近二三十年内发展起来的一种新坝型,因其具有断面小(与土坝相比)、就地取材、
施工导流简便、工期短和造价低等优点,已在世界范围内得到迅速推广.我国目前已建成或正在修建的面板
坝有50多座.最大坝高达180m.
由于混凝土面板坝的发展历史尚短,而且堆石坝体具有较高的强度,可以维持较陡的边坡,粗颗粒形成
的骨架能很好地保持其渗透稳定性,因此在1993年青海省沟后水库钢筋混凝土面板砂砾石坝溃决之前,国
际上尚无一例面板坝因整体失稳或渗漏破坏而导致大坝溃决的实例.迄今为止,面板坝的设计基本上是经验
性的,采用工程类比的经验方法为多,很少进行分析计算,即使做了计算分析,也仅作为参考.这主要是由于
现有的面板坝计算理论还不够成熟和完善.
面板坝的计算理论包括稳定、变形、应力变及渗流分析等四个部分,沟后坝的溃决与这几方面都有直接
联系.沟后面板坝失事的惨痛教训表明,面板坝(特别是用弱透水材料填筑坝体的面板坝)的建设仅停留在经
验水平上是不够的,应在前人工作的基础上,继续发展和完善面板坝的设计计算方法,使其更加合理.
对于堆石或卵、砾石面板坝,因其坝体透水性极强,以及对于坝内有竖向有效排水系统的面板坝,因渗水
进入坝体后排水十分通畅,均不会形成高浸润面,大部分坝体无孔隙水压力和渗流作用.本文主要针对用弱
透水材料(砂砾,含泥石渣、卵砾石,含细粒多的风化料等)填筑的、且未设竖向排水的面板坝渗流计算方法进
行探讨,重点对面板与高挡墙底板水平接缝止水破坏后的渗流进行分析.
1 正常运行条件下的面板坝渗流分析
正常运行情况下的面板坝的渗流分析,主要是对面板缝隙产生的渗漏计算.
研究表明
[1]
,当面板厚度
t
p
(即入渗渗径)与缝隙宽度
b
之比(
t
p
/
b
)<3～4时,通过缝隙的水头损失为
零,即可忽略缝隙水头损失.这时,由于垫层具有半透水性,几乎90%以上的水头损失都集中在垫层区,垫层
区控制着渗流量.计算结果表明,此时缝隙的宽度对渗流量影响很小.当(t
p
/b)>3～4时,通过缝隙的水头损
失不可忽略,此时,缝隙的几何尺寸及缝面糙率将对渗流起控制作用.
钢筋混凝土面板的厚度一般大于30cm,而大量的观测资料表明
[1]
,面板的裂缝宽度为毫米量级.如国
内多座面板坝中,以小干沟坝的面板裂缝宽度最大,也仅为2mm.面板的接缝,特别是周边缝的张开宽度一
般较大,最大可达10cm左右.即使如此,在止水失效情况下,仍可基本满足(t
p
/b)>3～4的条件,故假定:¹
a
收稿日期:1999-06-24
作者简介:盛金保(1966-),男,安徽青阳人,工程师,硕士,主要从事大坝安全评价、监测与管理工作.
水 利 水 运 科 学 研 究 2000年6月
40
缝隙的几何尺寸及缝面糙率对渗流起控制作用;º坝体内渗流属层流运动,即满足达西定律;»堆石体各料
区均质;¼基岩相对不透水;½面板除缝隙外相对不透水.因此,可以近似地认为,通过面板缝隙的渗流量就
是整个坝体的渗流量.
根据C.路易斯的研究,当水流通过粗糙裂缝且为层流时,其流速
1
v=
gb
i
12T
1+8.8(K/D
n
)
1.5
2
[1]
(1)
式中:K为裂缝表面的绝对粗糙度(即粗糙高度);裂缝的水力直径D
n
=2b;b为裂缝宽度;k/D
n
为相对粗糙
度(对于光滑缝面为0,当粗糙高度等于裂缝宽度时为0.5);
g
为重力加速度;
T
为水的运动粘滞系数;
i
为裂
缝渗流水力比降,近似为作用水头与面板厚度之比.故通过单位长度裂缝的渗透量为
gb
3
1
q
f
=i
12T1+8.8(K/D
n
)
1.5
若相对粗糙度
K
/
D
n
=0.5,水温为20℃,则式(2)可简化为
q
f
=2×10
5
b
3
i
现取单位长度坝体作为分析对象,分段建立各坝区的流量方程,计算简图见图1.
(3)
(2)
图1 渗透流量计算简图
Fig
.1
Calculatingchartforseepagedischarge
通过面板的渗流量
q=
∑
lbi
3
nnn
×10
5
(4)
式中:l
n
为相应于缝宽为b
n
的缝隙总长;i
n
为裂缝渗流水力比降,并按下式计算:
H
1
-h
1
(x≤L
1
)
t
px
i
n
=
H
1
-xtanA
1
(x>L
1
)
t
px
式中:
H
1
,
h
1
分别为面板上、下游水深;
A
1
为面板坡角;
t
px
为(
x
,
xtanA
1
)处的面板厚度.
通过垫层区的渗流量按下式计算
[1]
:
q=
k
1
22222
1+m
2
1
(h
1
-h
2
-b
s
cos
A
1
sinA
1
2b
s
sinA
1
)
(5)
式中:m
1
,A
1
分别为垫层区的坡度系数和坡角;h
1
,h
2
分别为垫层区上、下游水深;k
1
为垫层区渗透系数;b
s
为
垫层厚度.
通过坝身的渗流量按下式计算
[2]
:
q=k
L-m
2
H
2
+
(h
2
-H
2
)
2
(L-m
2
H
2
)-m
2
(h
2
-H
2
)
222
(h
2
-H
2
)H
2
(h
2
-H
2
)T
+k
L-0.5m
2
H
2
+k
0
L+m
1
h
2
+0.88T
(6)
式中:H
2
为下游水深;T为坝基覆盖层厚度;m
2
为下游坝坡坡度系数;k为主堆石区渗透系数;k
0
为坝基覆
第2期 盛金保,等:混凝土面板坝渗流分析方法初探
41
盖层渗透系数.
联立解式(4)～(6),便可求得
q
,
h
1
和
h
2
.
渗透逸出高度h
0
可按下式计算
[2]
:
q
k1+
m
2
+0.5
(m
2
+0.5)(h
0
-H
2
)+
q
(m
2
+0.5)H
2
(m
2
+0.5)(h
0
-H
2
)+0.5H
2
[2]
(m
2
+0.5)H
2
m
2
H
2
k>k
0
+
k
0
T
(m
2
+0.5)(h
0
-H
2
)+m
2
H
2
h
0
-H
2
=
2(m
2
+0.5)
k
0
T
+0.44T
(7)
k
m
2
k≤k
0
1++
m
2
k
0
+0.44T
坝体内的浸润线计算式为
y-h
0
y
2
-h
2
0
x=-k
0
T+k
q′2q′
2
h
2
2
-h
0
h
1
-h
0
q′=k+k
0
T
m
1
L+m
2
h
1
-m
1
h
0
2(L+h
2
-m
2
h
0
)
2m
1
+1
(8)
其中:
也可采用有限单元法对面板坝进行渗流分析.此时,只需将缝隙简化为水头边界或流量边界:当(t
p
/b)
<3～4时,忽略通过缝隙的水头损失,将缝隙简化为一条线,水头通过该线直接作用在垫层上;当(
t
p
/
b
)>3
～4时,将缝隙处理成流量边界,作用在垫层上的流量按式(2)或式(3)计算.
有限元分析方法的支配方程为:
55h55h55h
(k
x
)+(k
y
)+(k
z
)=0
5x5x5y5y5z5z
边界条件为:
水头边界:hû
S
1
=f
1
(x,y,z);流量边界:hû
S
2
=f
2
(x,y,z)
式中:
h
(
x
,
y
,
z
)为水头分布函数;
k
x
,
k
y
,
k
z
分别为以
x
,
y
,
z
轴为主轴方向的渗透系数;
f
1
(
x
,
y
,
z
)为渗流域
边界上的水头分布函数;f
2
(x,y,z)为渗流域边界上的流量分布函数.
(9)
2 高挡墙底板与面板间水平接缝止水破坏后的坝体渗流分析
高挡墙底板与面板间的水平接缝止水遭到破坏后,当库水超过该缝高程时,将成为重要的渗流通道.
面板与高挡墙水平接缝位于坝顶,张开变形较大,同时缝两边的相对沉降量也较大,致使过缝渗径变短.
例如,沟后水库大坝面板与高挡墙水平接缝张开宽度b=6cm,缝两边的相对沉降差$s=14cm,而顶部面板
与高挡墙底板厚t
p
=30cm,故过缝渗径l=t
p
-$s=16cm,l/b=2.67<3～4,接缝水头损失可忽略.顶部面
板与高挡墙底板间的相对沉降量过大,还会造成接缝下的局部区域垫层受剪切作用而破坏,起不到控制渗流
的作用.所以,计算分析中可将该缝视为一条等水头线性边界,库水通过该缝直接作用于主堆石区.
由于接缝下的局部区域垫层因剪切变形过大遭到破坏,能承受的作用水头可能很小,使得其下游部分范
围内的主堆石区内渗流不再满足达西定律.但对主堆石区为弱透水材料的面板坝来说,这些超过达西定律上
限的局部区域与整个流场相比是很小的,仍可按达西定律来分析整个流场的渗流,即渗流支配方程仍为拉普
拉斯方程.
3 高挡墙与面板接缝止水防渗失效时的坝体浸湿过程
大量原观资料表明
[1]
,通过面板与垫层后的渗流剩余水头一般占总渗流水头的10%左右,因此,在面板
与高挡墙水平接缝开始进水之前,坝体内的浸润面总是很低的.如果坝体内设有竖直排水设施并与面板与高
挡墙水平接缝相接,因排水通畅,过缝渗水也不会形成高浸润面.否则有可能抬高浸润面,影响坝体的渗透稳
定和下游坝坡的抗滑稳定性.
水 利 水 运 科 学 研 究 2000年6月
42
库水自面板与高挡墙水平接缝渗入坝内起,将经历以下两个阶段:
(1)非稳定自由渗漏阶段 过缝渗水从缝底垫层逐步湿润至坝面和起始浸润线,并不断抬高坝体浸润线.
(2)顶托渗漏阶段 渗水不断补给,坝体浸润线将逐步上升而达缝底.当坝体浸润线与库水连成一体时,
过缝渗漏将受到坝内水体的顶托,在经过一定的时间以后,达到稳定状态,即可按前述方法进行渗流计算.
现分析过缝渗水饱和浸湿坝体的过程,并推导渗水前锋线在非饱和坝体中的推进速度和距离.计算简图
见图2.
图2 渗水前锋线推进计算简图
Fig
.2
Calculatingchartforfrontlineofseepagewater
由于非饱和土体中存在毛细管吸力,渗流水力比降按下式计算:
i=(h+h
c
)/l
式中:h为作用水头;h
c
为吸压力水头;l为渗径.
当不能保证在接缝下面局部区域内的渗流为层流运动时,可根据雷诺数来判断其流态:
Re=
vd
T
式中:v为坝体中的渗流速度;d为坝体料的平均颗粒直径;T为水的运动粘滞系数.达西定律上限的临界雷
诺数Re约在1～10之间,或确切一些说Re=5.因此可约定,当Re<5时,为层流运动;当Re>5时,开始
进入紊流状态,且达西定律不再成立.流速与渗流水力比降之间的关系可表示为指数的形式:
v=ki
m
约从1降到0.7.
对于层流运动,
t
时刻渗水前锋线上任一水质点
M
(见图2)的推进速度:
)
h(t)+lcosH+h
c
v
m
(t)=
dl
=
K(H
dtn
e
l
(12)
(11)
式中的指数m决定于坝体填筑料的密实程度与有效粒径.根据试验,对于粗砂和砾(d
10
>0.4mm)的指数m
[2]
[2]
(10)
式中:H为l与铅垂线夹角;K(H)为沿H方向的渗透系数;n
e
为非饱和土的气孔体积孔隙率;h(t)为t时刻缝
上水深;
l
为渗径,即缝至
M
点的距离;
h
c
为非饱和坝体的吸压力水头.
解式(12)就可确定任一时刻过缝渗水前锋面的位置,进而可计算渗水自坝背水坡出逸的最短时间,也可
计算过缝渗水何时到达坝体下部的饱和区.
4 算 例
对沟后水库大坝面板与高挡墙水平接缝渗水浸湿坝体的过程进行计算.为简单起见,只考虑沿水平和铅
直向的流动,并假定渗流为层流运动.
(1)沿高挡墙底板底的水平向流动(H=90°):K(90)=1×10cm/s=36m/h;h(t)=0.3+0.033t,
其中,高挡墙底板厚为0.3
m
,
t
的单位为
h
;上部坝体的平均孔隙率
n
=0.265,含水量
X
=3.81%,干重度
C
d
=19.70kN/m,算得n
e
=0.190;因高挡墙底板下存在水平向卵砾石透水层,故取吸压力水头h
c
=0.将所
有变量代入式(12),并分离变量后得
3[3]
00[3][4]
第2期 盛金保,等:混凝土面板坝渗流分析方法初探
43
(13)
ldl=(113.68+6.25t)dt
lû
t=0
=0
对上式进行积分后,得
l=3.125t
2
+113.68t(14)
将水平渗径长度l=13.45m代入上式,可求得t=1.47h.也就是说,如果高挡墙底板下存在水平向强
透水层(
k
=1×10
0
cm
/
s
),过缝渗水只需约1.5
h
就会从坝背水坡出逸.
(2)沿垂直向的流动(H=0°):K(0)=1×10
-2
cm/s=0.36m/h;n
e
=0.153
[3]
;h
c
=0.3m;为简化计算,取
缝上平均水深,即h(t)=0.45m.将所有变量代入式(12),并简化后得
ldl
=dt
1.76+2.35l
lû
t=0
=0
积分后,得
1
2
[1.76+2.35l-1.76ln(1.76+2.35l)]-0.14(16)
2.35
假定在接缝开始渗水前,垫层后剩余水头为10%,则接缝到达起始浸润线的垂直距离
l
=50
m
,将其代
入式(16)得
t
=19.93
h
,即过缝渗水需要约20
h
才能到达坝体下部的饱和区.
t=
(15)
5 结 语
本文的混凝土面板坝渗流分析方法,可以计算面板坝接缝止水全部或部分失效情况下的渗漏总量及坝
体浸润线,从而能核算在不利运行条件下的坝体渗透稳定和整体稳定性.此外,还可估算在面板与高挡墙水
平接缝进水后,渗流自坝体背水坡出逸的时间,进而估计自接缝进水到发生坝坡滑动的时间,为防洪减灾决
策提供参考.
参考文献:
[1] 傅志安,风家骥.混凝土面板堆石坝[M].武汉:华中理工大学出版社,1993.
[2] 毛昶熙.渗流计算分析与控制[
M
].北京:水利电力出版社,1990.
[3] 李 雷,盛金保.沟后坝砂砾料的工程特性分析[R].南京:水利部大坝安全管理中心.1996.
[4] 盛金保.沟后坝溃坝渗流初步分析[J].大坝观测与土工测试,1996,(3).
Preliminarystudyonseepageanalysismethodsofconcrete
-
facedams
SHENGJin
-
bao
,
LILei
(NanjingHydraulicResearchInstitute,DamSafetyManagementCenteroftheMinistryofWater
Resources,Nanjing 210029,China)
Abstract:Seepagecharacteristicsofconcrete-facedams,especiallyafterthejointsealbetweenthefaceand
thebaseplateofthehighwavewallistorn,hodstoanalyze
theseepageofconcretefacedamsarealsopresented
,
whichmaybehelpfulfordesigningconcrete
-
face
damsinfuture.
Keywords:concrete-facedam;seepage;analysismethod
2024年5月6日发(作者：查夏萱)
南京
第2期
2000年6月
水利科学研究院
No.2
Jun.2000
水利水运科学研究
JOURNALOFNANJINGHYDRAULICRESEARCHINSTITUTE
混凝土面板坝渗流分析方法初探
盛金保,李 雷
(南京水利科学研究院,水利部大坝安全管理中心,江苏南京 210029)
a
摘要:
对混凝土面板坝的渗流分析方法,重点是面板与高挡墙底板间水平接缝止水破坏后的渗流分析方法进行了
探讨,可供面板坝设计参考.
关键词:
混凝土面板坝;渗流;分析方法
中图分类号:TV641.43 文献标识码:A 文章编号:1001-3962(2000)02-0039-05
混凝土面板堆石坝是近二三十年内发展起来的一种新坝型,因其具有断面小(与土坝相比)、就地取材、
施工导流简便、工期短和造价低等优点,已在世界范围内得到迅速推广.我国目前已建成或正在修建的面板
坝有50多座.最大坝高达180m.
由于混凝土面板坝的发展历史尚短,而且堆石坝体具有较高的强度,可以维持较陡的边坡,粗颗粒形成
的骨架能很好地保持其渗透稳定性,因此在1993年青海省沟后水库钢筋混凝土面板砂砾石坝溃决之前,国
际上尚无一例面板坝因整体失稳或渗漏破坏而导致大坝溃决的实例.迄今为止,面板坝的设计基本上是经验
性的,采用工程类比的经验方法为多,很少进行分析计算,即使做了计算分析,也仅作为参考.这主要是由于
现有的面板坝计算理论还不够成熟和完善.
面板坝的计算理论包括稳定、变形、应力变及渗流分析等四个部分,沟后坝的溃决与这几方面都有直接
联系.沟后面板坝失事的惨痛教训表明,面板坝(特别是用弱透水材料填筑坝体的面板坝)的建设仅停留在经
验水平上是不够的,应在前人工作的基础上,继续发展和完善面板坝的设计计算方法,使其更加合理.
对于堆石或卵、砾石面板坝,因其坝体透水性极强,以及对于坝内有竖向有效排水系统的面板坝,因渗水
进入坝体后排水十分通畅,均不会形成高浸润面,大部分坝体无孔隙水压力和渗流作用.本文主要针对用弱
透水材料(砂砾,含泥石渣、卵砾石,含细粒多的风化料等)填筑的、且未设竖向排水的面板坝渗流计算方法进
行探讨,重点对面板与高挡墙底板水平接缝止水破坏后的渗流进行分析.
1 正常运行条件下的面板坝渗流分析
正常运行情况下的面板坝的渗流分析,主要是对面板缝隙产生的渗漏计算.
研究表明
[1]
,当面板厚度
t
p
(即入渗渗径)与缝隙宽度
b
之比(
t
p
/
b
)<3～4时,通过缝隙的水头损失为
零,即可忽略缝隙水头损失.这时,由于垫层具有半透水性,几乎90%以上的水头损失都集中在垫层区,垫层
区控制着渗流量.计算结果表明,此时缝隙的宽度对渗流量影响很小.当(t
p
/b)>3～4时,通过缝隙的水头损
失不可忽略,此时,缝隙的几何尺寸及缝面糙率将对渗流起控制作用.
钢筋混凝土面板的厚度一般大于30cm,而大量的观测资料表明
[1]
,面板的裂缝宽度为毫米量级.如国
内多座面板坝中,以小干沟坝的面板裂缝宽度最大,也仅为2mm.面板的接缝,特别是周边缝的张开宽度一
般较大,最大可达10cm左右.即使如此,在止水失效情况下,仍可基本满足(t
p
/b)>3～4的条件,故假定:¹
a
收稿日期:1999-06-24
作者简介:盛金保(1966-),男,安徽青阳人,工程师,硕士,主要从事大坝安全评价、监测与管理工作.
水 利 水 运 科 学 研 究 2000年6月
40
缝隙的几何尺寸及缝面糙率对渗流起控制作用;º坝体内渗流属层流运动,即满足达西定律;»堆石体各料
区均质;¼基岩相对不透水;½面板除缝隙外相对不透水.因此,可以近似地认为,通过面板缝隙的渗流量就
是整个坝体的渗流量.
根据C.路易斯的研究,当水流通过粗糙裂缝且为层流时,其流速
1
v=
gb
i
12T
1+8.8(K/D
n
)
1.5
2
[1]
(1)
式中:K为裂缝表面的绝对粗糙度(即粗糙高度);裂缝的水力直径D
n
=2b;b为裂缝宽度;k/D
n
为相对粗糙
度(对于光滑缝面为0,当粗糙高度等于裂缝宽度时为0.5);
g
为重力加速度;
T
为水的运动粘滞系数;
i
为裂
缝渗流水力比降,近似为作用水头与面板厚度之比.故通过单位长度裂缝的渗透量为
gb
3
1
q
f
=i
12T1+8.8(K/D
n
)
1.5
若相对粗糙度
K
/
D
n
=0.5,水温为20℃,则式(2)可简化为
q
f
=2×10
5
b
3
i
现取单位长度坝体作为分析对象,分段建立各坝区的流量方程,计算简图见图1.
(3)
(2)
图1 渗透流量计算简图
Fig
.1
Calculatingchartforseepagedischarge
通过面板的渗流量
q=
∑
lbi
3
nnn
×10
5
(4)
式中:l
n
为相应于缝宽为b
n
的缝隙总长;i
n
为裂缝渗流水力比降,并按下式计算:
H
1
-h
1
(x≤L
1
)
t
px
i
n
=
H
1
-xtanA
1
(x>L
1
)
t
px
式中:
H
1
,
h
1
分别为面板上、下游水深;
A
1
为面板坡角;
t
px
为(
x
,
xtanA
1
)处的面板厚度.
通过垫层区的渗流量按下式计算
[1]
:
q=
k
1
22222
1+m
2
1
(h
1
-h
2
-b
s
cos
A
1
sinA
1
2b
s
sinA
1
)
(5)
式中:m
1
,A
1
分别为垫层区的坡度系数和坡角;h
1
,h
2
分别为垫层区上、下游水深;k
1
为垫层区渗透系数;b
s
为
垫层厚度.
通过坝身的渗流量按下式计算
[2]
:
q=k
L-m
2
H
2
+
(h
2
-H
2
)
2
(L-m
2
H
2
)-m
2
(h
2
-H
2
)
222
(h
2
-H
2
)H
2
(h
2
-H
2
)T
+k
L-0.5m
2
H
2
+k
0
L+m
1
h
2
+0.88T
(6)
式中:H
2
为下游水深;T为坝基覆盖层厚度;m
2
为下游坝坡坡度系数;k为主堆石区渗透系数;k
0
为坝基覆
第2期 盛金保,等:混凝土面板坝渗流分析方法初探
41
盖层渗透系数.
联立解式(4)～(6),便可求得
q
,
h
1
和
h
2
.
渗透逸出高度h
0
可按下式计算
[2]
:
q
k1+
m
2
+0.5
(m
2
+0.5)(h
0
-H
2
)+
q
(m
2
+0.5)H
2
(m
2
+0.5)(h
0
-H
2
)+0.5H
2
[2]
(m
2
+0.5)H
2
m
2
H
2
k>k
0
+
k
0
T
(m
2
+0.5)(h
0
-H
2
)+m
2
H
2
h
0
-H
2
=
2(m
2
+0.5)
k
0
T
+0.44T
(7)
k
m
2
k≤k
0
1++
m
2
k
0
+0.44T
坝体内的浸润线计算式为
y-h
0
y
2
-h
2
0
x=-k
0
T+k
q′2q′
2
h
2
2
-h
0
h
1
-h
0
q′=k+k
0
T
m
1
L+m
2
h
1
-m
1
h
0
2(L+h
2
-m
2
h
0
)
2m
1
+1
(8)
其中:
也可采用有限单元法对面板坝进行渗流分析.此时,只需将缝隙简化为水头边界或流量边界:当(t
p
/b)
<3～4时,忽略通过缝隙的水头损失,将缝隙简化为一条线,水头通过该线直接作用在垫层上;当(
t
p
/
b
)>3
～4时,将缝隙处理成流量边界,作用在垫层上的流量按式(2)或式(3)计算.
有限元分析方法的支配方程为:
55h55h55h
(k
x
)+(k
y
)+(k
z
)=0
5x5x5y5y5z5z
边界条件为:
水头边界:hû
S
1
=f
1
(x,y,z);流量边界:hû
S
2
=f
2
(x,y,z)
式中:
h
(
x
,
y
,
z
)为水头分布函数;
k
x
,
k
y
,
k
z
分别为以
x
,
y
,
z
轴为主轴方向的渗透系数;
f
1
(
x
,
y
,
z
)为渗流域
边界上的水头分布函数;f
2
(x,y,z)为渗流域边界上的流量分布函数.
(9)
2 高挡墙底板与面板间水平接缝止水破坏后的坝体渗流分析
高挡墙底板与面板间的水平接缝止水遭到破坏后,当库水超过该缝高程时,将成为重要的渗流通道.
面板与高挡墙水平接缝位于坝顶,张开变形较大,同时缝两边的相对沉降量也较大,致使过缝渗径变短.
例如,沟后水库大坝面板与高挡墙水平接缝张开宽度b=6cm,缝两边的相对沉降差$s=14cm,而顶部面板
与高挡墙底板厚t
p
=30cm,故过缝渗径l=t
p
-$s=16cm,l/b=2.67<3～4,接缝水头损失可忽略.顶部面
板与高挡墙底板间的相对沉降量过大,还会造成接缝下的局部区域垫层受剪切作用而破坏,起不到控制渗流
的作用.所以,计算分析中可将该缝视为一条等水头线性边界,库水通过该缝直接作用于主堆石区.
由于接缝下的局部区域垫层因剪切变形过大遭到破坏,能承受的作用水头可能很小,使得其下游部分范
围内的主堆石区内渗流不再满足达西定律.但对主堆石区为弱透水材料的面板坝来说,这些超过达西定律上
限的局部区域与整个流场相比是很小的,仍可按达西定律来分析整个流场的渗流,即渗流支配方程仍为拉普
拉斯方程.
3 高挡墙与面板接缝止水防渗失效时的坝体浸湿过程
大量原观资料表明
[1]
,通过面板与垫层后的渗流剩余水头一般占总渗流水头的10%左右,因此,在面板
与高挡墙水平接缝开始进水之前,坝体内的浸润面总是很低的.如果坝体内设有竖直排水设施并与面板与高
挡墙水平接缝相接,因排水通畅,过缝渗水也不会形成高浸润面.否则有可能抬高浸润面,影响坝体的渗透稳
定和下游坝坡的抗滑稳定性.
水 利 水 运 科 学 研 究 2000年6月
42
库水自面板与高挡墙水平接缝渗入坝内起,将经历以下两个阶段:
(1)非稳定自由渗漏阶段 过缝渗水从缝底垫层逐步湿润至坝面和起始浸润线,并不断抬高坝体浸润线.
(2)顶托渗漏阶段 渗水不断补给,坝体浸润线将逐步上升而达缝底.当坝体浸润线与库水连成一体时,
过缝渗漏将受到坝内水体的顶托,在经过一定的时间以后,达到稳定状态,即可按前述方法进行渗流计算.
现分析过缝渗水饱和浸湿坝体的过程,并推导渗水前锋线在非饱和坝体中的推进速度和距离.计算简图
见图2.
图2 渗水前锋线推进计算简图
Fig
.2
Calculatingchartforfrontlineofseepagewater
由于非饱和土体中存在毛细管吸力,渗流水力比降按下式计算:
i=(h+h
c
)/l
式中:h为作用水头;h
c
为吸压力水头;l为渗径.
当不能保证在接缝下面局部区域内的渗流为层流运动时,可根据雷诺数来判断其流态:
Re=
vd
T
式中:v为坝体中的渗流速度;d为坝体料的平均颗粒直径;T为水的运动粘滞系数.达西定律上限的临界雷
诺数Re约在1～10之间,或确切一些说Re=5.因此可约定,当Re<5时,为层流运动;当Re>5时,开始
进入紊流状态,且达西定律不再成立.流速与渗流水力比降之间的关系可表示为指数的形式:
v=ki
m
约从1降到0.7.
对于层流运动,
t
时刻渗水前锋线上任一水质点
M
(见图2)的推进速度:
)
h(t)+lcosH+h
c
v
m
(t)=
dl
=
K(H
dtn
e
l
(12)
(11)
式中的指数m决定于坝体填筑料的密实程度与有效粒径.根据试验,对于粗砂和砾(d
10
>0.4mm)的指数m
[2]
[2]
(10)
式中:H为l与铅垂线夹角;K(H)为沿H方向的渗透系数;n
e
为非饱和土的气孔体积孔隙率;h(t)为t时刻缝
上水深;
l
为渗径,即缝至
M
点的距离;
h
c
为非饱和坝体的吸压力水头.
解式(12)就可确定任一时刻过缝渗水前锋面的位置,进而可计算渗水自坝背水坡出逸的最短时间,也可
计算过缝渗水何时到达坝体下部的饱和区.
4 算 例
对沟后水库大坝面板与高挡墙水平接缝渗水浸湿坝体的过程进行计算.为简单起见,只考虑沿水平和铅
直向的流动,并假定渗流为层流运动.
(1)沿高挡墙底板底的水平向流动(H=90°):K(90)=1×10cm/s=36m/h;h(t)=0.3+0.033t,
其中,高挡墙底板厚为0.3
m
,
t
的单位为
h
;上部坝体的平均孔隙率
n
=0.265,含水量
X
=3.81%,干重度
C
d
=19.70kN/m,算得n
e
=0.190;因高挡墙底板下存在水平向卵砾石透水层,故取吸压力水头h
c
=0.将所
有变量代入式(12),并分离变量后得
3[3]
00[3][4]
第2期 盛金保,等:混凝土面板坝渗流分析方法初探
43
(13)
ldl=(113.68+6.25t)dt
lû
t=0
=0
对上式进行积分后,得
l=3.125t
2
+113.68t(14)
将水平渗径长度l=13.45m代入上式,可求得t=1.47h.也就是说,如果高挡墙底板下存在水平向强
透水层(
k
=1×10
0
cm
/
s
),过缝渗水只需约1.5
h
就会从坝背水坡出逸.
(2)沿垂直向的流动(H=0°):K(0)=1×10
-2
cm/s=0.36m/h;n
e
=0.153
[3]
;h
c
=0.3m;为简化计算,取
缝上平均水深,即h(t)=0.45m.将所有变量代入式(12),并简化后得
ldl
=dt
1.76+2.35l
lû
t=0
=0
积分后,得
1
2
[1.76+2.35l-1.76ln(1.76+2.35l)]-0.14(16)
2.35
假定在接缝开始渗水前,垫层后剩余水头为10%,则接缝到达起始浸润线的垂直距离
l
=50
m
,将其代
入式(16)得
t
=19.93
h
,即过缝渗水需要约20
h
才能到达坝体下部的饱和区.
t=
(15)
5 结 语
本文的混凝土面板坝渗流分析方法,可以计算面板坝接缝止水全部或部分失效情况下的渗漏总量及坝
体浸润线,从而能核算在不利运行条件下的坝体渗透稳定和整体稳定性.此外,还可估算在面板与高挡墙水
平接缝进水后,渗流自坝体背水坡出逸的时间,进而估计自接缝进水到发生坝坡滑动的时间,为防洪减灾决
策提供参考.
参考文献:
[1] 傅志安,风家骥.混凝土面板堆石坝[M].武汉:华中理工大学出版社,1993.
[2] 毛昶熙.渗流计算分析与控制[
M
].北京:水利电力出版社,1990.
[3] 李 雷,盛金保.沟后坝砂砾料的工程特性分析[R].南京:水利部大坝安全管理中心.1996.
[4] 盛金保.沟后坝溃坝渗流初步分析[J].大坝观测与土工测试,1996,(3).
Preliminarystudyonseepageanalysismethodsofconcrete
-
facedams
SHENGJin
-
bao
,
LILei
(NanjingHydraulicResearchInstitute,DamSafetyManagementCenteroftheMinistryofWater
Resources,Nanjing 210029,China)
Abstract:Seepagecharacteristicsofconcrete-facedams,especiallyafterthejointsealbetweenthefaceand
thebaseplateofthehighwavewallistorn,hodstoanalyze
theseepageofconcretefacedamsarealsopresented
,
whichmaybehelpfulfordesigningconcrete
-
face
damsinfuture.
Keywords:concrete-facedam;seepage;analysismethod