技术交流 SMALL HYDRO POWER 2011No2,Total No158 中小型水库溃坝洪水估算 李庆华(云南省水文水资源局临沧分局 云南临沧677000) 【摘要】对水库溃坝洪水进行估算,可预测洪水对下游造成的淹没范围和程度,以便做出防范措施和编制应急计划。 介绍了溃坝洪水估算的方法,对于在溃坝数据不足及计算工具不完善情况下的溃坝洪水计算有一定的帮助。图2幅,表4个。 【关键词】中小型水库渍坝 洪水估算 表1 瞬时全溃情况下几种河槽溃坝最大流量计算公式 修建大坝、利用水库蓄水调节天然流量是综合 利用水资源的主要手段,也是优化水资源配置的主 要方法。但是如果因某些偶然因素导致坝体突遭破 坏溃决,水库蓄水和上游入库洪水经溃口下泄就会 形成骤发性的洪水,给下游人民的生命财产带来严 重危害。因此,有必要分析计算溃坝洪水,预估大 坝一旦失事后溃坝洪水对下游造成的危害,以便做 出防范措施和编制应急计划。 注:式中,曰为坝址断面的平均宽度(m);H为溃坝 前的坝前水深(m);Q 为溃坝时坝址处最大流量(m3/s); h为坝下游恒定流水深(m)。 1计算方法 1.1坝址溃坝最大流量计算 水库大坝可分为混凝土坝和土石坝,混凝土坝 又可分为重力坝、拱坝和支墩坝3种类型;土石坝 1.1.2横向局部溃决情况 在瞬时全溃情况计算公式的基础上考虑大坝溃 口的宽度,即为横向局部溃决情况下的计算公式。 根据不同断面情况其计算式也不同(见表2)。 表2横向局部溃决情况下几种河槽溃坝最大流量计算公式 则包括土坝、堆石坝、土石混合坝等。根据大坝的 坝型、筑坝材料、基础条件、溃坝原因等因素,大 坝溃决可分为瞬时全部溃决、横向部分溃决到底及 横向和纵向均为部分溃决3种类型。不同类型有不 同的计算方式,在实际应用中应选取最不利情况进 行分析计算。 1.1.1瞬时全溃情况 瞬时全溃情况主要是根据不稳定流基本方程 式,在假定河床底坡为0且忽略阻力的基础上,针 对河槽的各种不同断面形状又有不同的计算方式 (见表1)。 注:式中,b为溃口宽度(m),小型水库取b为 , 中型水库b为0.6B一0.7B;其余各符号意义同表1。 1949年,肖克利奇又根据试验提出了1种经 验公式: Qm= (鲁) /2 1.1.3横向和纵向均为部分溃决情况 1)这种情况下多采用此经验公式 Qm= ( )6 (日 ) 收稿日期:2011一Ol一26 作者简介:李庆华(1980一),男,助理工程师,主要从事 式中,h 为坝残留高度(m);其余各符号意 义同表1。 水文水资源规划设计、水情预报等工作。 Email:lqh198236@126.tom ・ 2)铁道科学院提出的经验公式 8 ・ 小水电2011年第2期(总第158期) qm=0.27 L) ( ) Bh(It一技术交流 ) /2 水库水量泄空时间(s); 为溃坝时水库人流量 式中, 为水库长度(km); 为系数。 3)美国水道实验站公式 (m3/s);K为系数,与溃坝流量过程线线型有关, 对4次抛物线来说,K一般为4~5,对2.5次抛物 线来说, 为3.5。 1.3渍坝流量过程线推求 qm= 丽BH)。_28Bh /2 1.2水库泄空时间计算 T=KW/(Q—qo) 溃坝流量过程采用表3所列无因次概化过程 线,根据已计算出的Q 和71便可求得概化溃坝流 量过程线的纵横坐标为:t=aT;q= (见表3)。 式中, 为水库溃坝后下泄水量(m3);T为 表3无因次概化溃坝泄流过程线 注:式中,t为水库溃坝后泄流时间(s);Q 为水库泄流时间t时的流量(m3/s);Qo为人库流量(m]/s)。 1.4坝下游断面最大流量计算 后,将断面流量过程线概化为三角形,从而求得 。 T3= 当溃坝址处最大流量向下游演进至距坝址£ (m)流程时的最大流量,采用下式估算: nVL—W/Q +L/(一一一 . ) 2实例计算 某拟建水库位于云南山区,工程区河谷呈v 字型。根据设计,水库为1座以发电为主兼顾防洪 式中,pI 为距坝址 (in)控制断面溃坝最大 流量(m3/s);V为洪水期河道断面的最大流速 (m/s);K为经验系数。 的中型水库,大坝为混凝土双曲拱坝,顶拱最大中 心角101.6。,大坝顶拱弧长274.8 m,最大坝高为 67 in。大坝洪水设计标准为100年一遇,校核洪水 有流量与断面资料时, 可采用特大洪水的最 大流速。无资料时,山区一般可采用3~5 m/s;半 山区一般可采用2 3 m/s;平原区可采用1—2 m/s。 一标准为500年一遇,水库校核洪水位下对应库容为 7 400万m3。由于水库下游14 km为1个重要集镇, 对溃坝洪水进行估算对下游集镇的防洪减灾具有重 要意义。 般情况下,山区 可采用1.1—1.5,半山 区采用1.0,平原区采用0.8~0.9。 1.5洪水传播时间计算 1)洪水起涨时间 。 2.1坝址溃坝最大流量计算 该水库大坝型式为混凝土双曲拱坝,为刚性 坝,考虑为瞬时全部溃决,按圣维南公式计算。由 为起涨时间(s);h0为下游断面初 为系数,取(0.5~7.5)×10~。 ,1.4 式中, 于坝下游恒定流水深未知,只能从工程区河道河槽 断面形状来选择计算公式。工程区河道河槽呈V 字型,故从表1中选择断面形状为等腰三角形的公 始水深(m); 2)最大流量到达时间 712= 2—W—02 H o.5—h—0.25 .式计算,通过计算,坝址溃坝最大流量为85 400 m3/s。 m2.2水库泄空时间计算 式中, (m)。 为最大流量到达时间(s); 2为系 与溃坝流量相比,溃坝时入库流量显得很小, 故忽略溃坝时入库流量,根据1.2节公式,得水库 泄空时间为3 900 s。 2.3溃坝流量过程线推求 数,取0.8~2.0;h 为最大流量时的平均水深 3)溃坝下游断面流量过程线 一般说来,溃坝洪水起涨陡,峰值到达快,峰 忽略溃坝时入库流量,根据1.3节公式,得水 库溃坝流量过程线(见图1、表4)。 ・后流量下降比较慢,根据上述公式求得 l及 9・ 技术交流 SMALL HYDRO POWER 2011No2,Total No158 图2溃坝洪水达到下游集镇处的概化洪水过程线 图1 水库溃坝泄流过程线 2.4溃坝洪水到达下游集镇最大流量计算 该工程位于山区,根据1.4节公式,相关系数 3结语 溃坝研究的目的是确定危害范围和洪峰的传播 时间,以便做出防范措施和编制应急计划。 溃坝洪水的计算由于情况复杂,资料又往往不 足,各种方法的计算成果可能相差较大。目前对于 按山区河道情况选取,计算得溃坝洪水到达下游集 镇最大流量为19 600 m3/s。 2.5溃坝洪水到达下游集镇时间计算 根据1.5节公式计算。由于无实测资料,溃坝 十分重要的水库,都倾向于用物理模型结合详细算 法进行,而简化公式仅作为粗估之用。 ■ 洪水到达集镇时对应的平均水深h 假定为9 m, 汁算得洪水起涨时间 为988 S,最大流量到达时 间 为1 200 s,退水时间 为8 555 S(见图2)。 责任编辑—吴昊 S屯.;|L . L —喜屯 ; .喜屯 舢 .喜 .5止.址.S .址.S屯 .工IL.j L.S^L (上接第26页)较差。 裂隙作为滑裂体的底裂面 (摩擦系数 :0.3,凝聚力C =0.05 MPa),在拱 不得超过1%e。基岩钻 ̄L:fL深5.2 m,加长杆(声26 钢管)长5 m,底环采用j580法兰接头。底环以上 0.4 m及底环以下0.2 m用膨胀水泥砂浆填塞。 加长杆顶部连接CF_12测缝针(规格为拉 12 Inln,压1 mln),外套保护管(声18O钢管,长 269 1TIITI)与基岩锚固。保护管顶端加盖板(马蹄 形钢板),引出观测电缆进入设置在坝后桥的集线 箱内,通过数字比例电桥进行变形观测。 坝平面布置时尽量使拱端上移,使I_2裂隙在拱槽 出露点位于坝体上部;同时,在满足坝肩稳定的前 提下,确定拱端开挖深度(开挖深度达15 m)。在 采取上述措施后,还有必要在 裂隙附近上下 30 m范围内进行基岩变形观测。 2.2坝肩变形观测设计 沿左坝肩442 m、457 m、472 m高程等高线距 3结语 左坝肩拱端下游10 m、30 m、50 m基岩表面布置测 点,共计9个测点,用以观测左坝肩的位移。采用 坝体位移结合库前水位、坝内温度的实测值, 相对封拱灌浆时的测时,通过逐步回归分析,可建 立回归方程。根据水位、温度的变化以及观测的时 三角控制网(与坝体表面测点采用同一控制网)进 行左坝肩位移观测。为了更加直观、精确观测左坝 肩的位移,在左坝肩拱槽442 m,457 m,472 m高 程布置3只基岩变形仪。 基岩变形仪埋设在1/2拱端坝厚处,与水平面 夹角为300。钻机造孔,孔径 100 1Tlln,最大孑L斜 ・ 间,即可计算坝体的位移,可以长期监视大坝的安 全运行。 ■ 责任编辑吴昊 1 0・