梯级水库连续溃坝计算.docx

梯级水库连续溃坝计算陈春燕1，李远强2，吴春耀1，刘锐1，许征3(1.珠江水文水资源研究所勘测中心，广东广州510370;广东省水文局梅州水文分局，广东梅州514021;3.中水珠江规划勘测设计有限公司，广东广州510370)2.摘要:为了合理地规划水库群的防洪设计标准，应给出水库溃坝流量过程甚至梯级水库连续溃坝的流量过程。该文选取2个梯级水库(上游水库的库容远大于下游水库库容)，计算上游水库溃坝洪水，并研究上游水库溃决引起下游水库连续溃坝的可能性，采用简化水力计算法、经验公式等方法计算下游水库溃坝可能发生的时间以及连续溃坝流量过程，为类似梯级水库工程的设计和现有梯级水库的调度管理提供参考。关键词:混凝土坝;连续溃坝;溃坝最大流量;溃坝洪水演进;调洪演算中图分类号:TV642文献标识码:B文章编号:1008－0112(2011)09－0040－03随着水库群的发展，为了合理地规划水库群的防洪设计标准，必须给出水库溃坝流量过程;若梯级水库中上游水库的库容较大，远大于下游的水库库容，上游坝体溃决后，坝下游形成溃坝涌波和洪流巨浪，可能危及到下游一定距离范围内的水库安全，则需要考虑梯级水库的连续溃坝，以便合理地确定水库工程等级、防洪设计标准和下游安全防护措施，或为改善现有重力坝安全性以及制定避险措施的重要依据［1－3］。1工程概况某河道梯级水库，上游水库A距离下游水库B为630.83km。水库A为碾压混凝土重力坝，水库正常蓄水位为400m，对应坝前水深为200m，最大坝高为216.5m，正常蓄水位库容为272.7亿m3，枯水期河面高程为219m，水面宽为90～100m，水深为1～19.5m，两岸岸坡为32°～42°，坝下恒定流水深为19m;水库B为碾压混凝土重力坝，坝顶高程为64.5m，对应坝前水深为44.5m，库容为38.47亿m3，死水位为47.6m，正常蓄水位为61.0m，坝顶长为1047.98m。坝址河谷呈“U”形河谷，谷底宽为400～500m，两岸坡度为25°～40°，局部有陡立悬崖，坝下恒定流水深为13m。河道断面接近为抛物线型，则河谷断面形状指数约为1.5。可见，上游水库A的库容大大于下游水库B的库容，一旦上游水库发生溃决，有必要进行下游水库连续溃坝的分析或计算。22.1计算方法坝址断面溃坝流量过程线1)水库坝址溃坝最大流量对于瞬时溃坝，其水流流态是一种非恒定的不连续波的运动。综合连续波与间断波的解法，用抛物线概化河谷断面，去掉对下游水深和流速的控制条件得:溃坝最大流量公式:gH3/2q=λB0槡0m2m+1v0??=mm－1?2槡m+式中λ为流量参数，λ?;mgH??0?槡?1+2m为河谷断面形状指数，矩形、三角形、二次抛物线形、组合抛物线型分别为1、2、1.5、3;B为坝址处河谷宽度，0m;g为重力加速度，9.8m/s2;H为上游水深，m［4］。0坝址断面溃坝流量过程线2)溃坝流量过程线线型可概化为4次抛物线型或2.5次抛物线型，即溃坝初瞬流量突增到Q，在很短时间内max流量迅速降为溃坝前入库流量Q0形成下凹曲线，在工程上，多选择4次抛物线。先拟定过程线总历时T，验算过程线与Q=Q0直线间的水量是否等于溃坝库容，如不相等，则须调整初步定的T值，直到两者相等为止［4］。溃坝波的河道演进在下游坝址处，溃坝洪水波到达时，往往是以立波方式传递的，因此采用“线性河道”假定，但在峰顶以后2.2收稿日期:2011－06－24;修回日期:2011－08－15作者简介:陈春燕(1985－)，女，硕士，助理工程师，主要从事建设项目防洪影响评价，建设项目水资源论证、水文水资源调查评价，核电厂等重要工程可能最大暴雨、可能最大洪水计算研究。·40·陈春燕，等:梯级水库连续溃坝计算2011年9月第9期No．9Sep．2011的洪水水深过程线，线性化的结果带来的误差就较大，因此线性河道法主要用来计算下游最大水深及对应的时间。由于峰顶流量的能坡比降近似于恒定流的水面比降，而峰前较陡，峰后较缓。故在获得最大水深后，若峰前峰后比降皆用恒定流比降，将会使峰前峰后流量过程线带来一些误差，本文采用无因次概化过程线法。2)概化过程线法从概念上说，越往下游水深越小，但总历时越长，w总w系由坝址过程线可知总历时T==，其qmλ槡gH1.50中为与河槽断面指数有关之数，对于棱柱体河槽可视为常数，w总为洪水总量，也是常数，故T与水深的1.5次1)线性河道法方成反比，根据诺谟图作出两条珔T～珚Hmx曲线，即:［5］溃坝波以立波方式向下游传递，波速较大，摩擦损失与之相比，占次要地位，故忽略初始流速和摩擦损失，?h0h0珔T1=2珚Hmx－1.5=0.05H0(5)gHMgHt)。1(得t时间内平均流速为:C珔=珔T=0.1珚H－5.65+=0.42槡m槡m2mxH0式中m