牤牛河设计洪水计算书

1、牤牛河设计洪水计算书1 概述牤牛河发源于东营燕子沟，在大屯乡大屯村汇入兴洲河，流域面积344.76km2，河流长度33.5km。窟窿山水库坝址在滦平县城南7.5km处的安匠屯乡铧子炉村，水库以上控制面积142.2km2，占全河流域面积的41.2%，窟窿山水库大坝以上主河道长22km，河道平均坡降20.3。因坡陡流急、汛期洪水迅猛而得名“牤牛河”。牤牛河由南至北纵贯滦平县城，是滦平县城区主要行洪河道。线路于大屯乡营房村南侧约400m处跨越牤牛河，位于牤牛河入汇兴洲河河口上游约700m，根据河北省滦平县牤牛河县城下游段河道综合治理工程设计变更报告，河北省滦平县牤牛河县城下游段河道综合治理工程终止点

2、为国道112公路桥，该桥所处位置位于线路跨越牤牛河断面下游约500m,该段河道治理在原有河堤基础上修建防洪堤和平整河道，迎水面采用浆砌石仰斜式挡土墙，迎水坡边坡1：0.25，背水坡边坡1：0.15，顶宽1.1m，高5m，底宽1.85m；外边坡1：2。防洪堤内填筑砂砾石，采用浆砌石护岸。该段河道治理后，防洪标准为10年一遇。2 设计洪峰流量2.1设计暴雨根据承德地区水文实用图集，由“河北省承德地区年最大24小时降雨量多年平均值等值线图”和“河北省承德地区年最大3日降雨量多年平均值等值线图”，断面以上流域年最大24小时降雨量多年平均值位于70mm与80mm等值线之间，年最大3日降雨量多年平均值位于

3、80mm与85mm等值线之间，由此可见该流域的暴雨集中在1日，采用1日设计暴雨进行洪水计算符合该地区暴雨洪水特点。根据年最大24小时降雨量多年平均值等值线查得牤牛河流域最大24小时降雨量平均值为：H1=75（mm）由“河北省承德地区年最大24小时降雨量变差系数（Cv）等值线图”查得：Cv1=0.40，24小时Cv值的折算系数为1.11。Cv1 = Cv1×1.11=0.44；Cs1/Cv1 = 3.5；由“P-型曲线Kp值表(Cv=3.5)”查得最大24小时降雨量Kp值：P = 1，Kp = 2.48；P = 2，Kp = 2.21；P = 3.33，Kp = 2.02；P = 20

4、, Kp = 1.301；牤牛河流域各频率24小时设计暴雨量为：P = 1，H1% = 80×2.48 = 198.4（mm）；P = 2，H3.33% = 80×2.21 = 176.8（mm）；P = 3.33，H3.33% = 80×2.02 = 161.6（mm）；P = 20，H20% = 80×1.301 = 104.1（mm）；牤牛河流域位于燕山迎风区，计算区间流域面积202km2，查“各区不同历时点面折减系数表”，最大24小时点面折减系数为0.92：点面折算后，流域设计面暴雨量H为：频率流域24小时设计面暴雨量HP = 1198.4

5、15;0.92=182.5mmP = 2176.8×0.92=162.7mmP = 3.33161.6×0.92148.7mmP = 20104.1×0.92=95.8mm2.2河道长度/流域面积根据河北省滦平县牤牛河县城下游段河道综合治理工程设计变更报告，自窟窿山水库至国道112公路桥区间流域面积为202km2。根据1:50000地形图量算，主河道长度为14.306km，按等高线量算流域主河道比降，区间流域图见图2.2-1，比降量算表见2.2-1。图2.2-1 窟窿山水库至线路跨越断面区间流域图表2.2-1 窟窿山水库至牤牛河区间流域主河道比降量算表序号间距（m

6、）高程（m）1811.744021110.445031351.24604882.347051587.84806993.449071507.850082279.551091295.653010351.4540111387.755012637.460013109.765014690总长度14306比降9.5992.3区间设计洪峰流量分别采用承德水文图集中的合理化公式法和水利部科研院法计算线路跨越水文断面100年一遇、50年一遇、30年一遇和5年一遇的设计洪峰流量。 合理公式法根据承德水文手册，计算洪峰流量的合理公式如下：Q1%=0.278(i-)F （1）i= Sp /n （2）v=0.278L/

7、 （3）上式中：Q1%百年一遇设计洪峰流量（m3/s）；0.278单位换标系数；i暴雨强度（mm/h）；产流历时内平均入渗率（mm/h）；F流域面积（km2）；Sp雨力（mm/h）；汇流时间（h）；n暴雨递减指数；L河流长度（km）；v汇流速度（m/s）全流域的河道坡度为：J = 9.599。牤牛河流域属兴洲河的一级支流，由“河北省承德地区汇流参数分区图”可知，该流域属于区，查询“洪峰流量推理公式分区参数表”得：H=1.412，B=0.52，x=0.0036，y=1.88；查“河北省承德地区短历时暴雨递减指数等值线图”100年一遇暴雨衰减指数n1=0.65×1.03=0.67，n2=

8、0.75，n3=0.65×1.060.69；=H（L/J0.5）B = 1.412×（14.306/9.5990.5）0.52 = 3.13（h）；100年一遇设计洪水计算如下：Sp=H24P/241-n=182.5/241-0.69= 68.1（mm/h）；i= Sp /n=68.1/3.130.69=31.0（mm/h）；=xSpy=0.0036×68.11.88=10.1（mm/h）；Q1%=0.278(i-)F=0.278×（31.0-10.1）×202=1174（m3/s）；查承德水文图集“设计洪峰流量换算系数表”，潮河流域设计洪峰流

9、量换算系数及各频率洪峰流量见表2.3-1（30年一遇通过Cv、Cs/Cv查P-型曲线Kp值表，采用比值计算）。表2.3-1牤牛河设计洪水计算成果表频率CvCs/Cv换算系数洪峰流量（m3/s）1%0.842.51.011742%0.842.50.861010 3.3%0.842.50.74869 20%0.842.50.37434 水科院采用EHP软件中的水科院法计算区间设计洪水，由于流域内主要为山丘区，地面植被为零星杂草及低矮的树木，山区土层薄，多处岩石出露，土质多为黄土，汇流参数取m=1.2；损失系数取10.1mm/h，水科院法计算过程如附录1所示，计算结果如表2.3-2。表2.3-2 牤

10、牛河水科院法计算成果频率洪峰流量（m3/s）1%16002%13003.3%111020%512 区间设计洪水成果分析将合理化公式法与水科院法计算结果进行比较，水科院法各频率设计洪峰流量均明显高于合理化公式法，。合理化公式法中的参数取值均来自承德水文图集，能很好体现牤牛河流域特征，符合该流域洪水特点，合理化公式法计算成果相比于水科院法更为合理、可靠，因此选用合理化公式法计算成果，100年一遇、30年一遇和5年一遇洪水位设计洪峰流量分别为1134m3/s、839m3/s、420 m3/s。2.4 窟窿山水库洪水分析计算线路跨越牤牛河断面位于窟窿山水库下游14km处，所布设的水文断面来水为水库下泄

11、洪水与区间洪水的叠加组合。 窟窿山水库概况窟窿山水库于1958年3月开工兴建，位于滦平县城以南7.5km处的铧子炉村，控制流域面积142.2km2，总库容1430万立方米，调洪库容898万立方米，是一座以防洪为主结合供水的中型水利工程。水库设计标准为100年一遇，校核标准为1000年一遇。窟窿山中型水库于2005年完成了除险加固。窟窿山水库历史最大入库流量为948 m3/s，历史最大出库流量为334 m3/s，发生于1976年7月23日。 水库入库设计洪水根据河北省滦平县牤牛河县城下游段河道综合治理工程设计变更报告中已有成果，不同重现期设计洪水过程线见表2.4-1。表2.4-1 窟窿山水库不同

12、重现期设计洪水过程线成果表时间(h)不同重现期设计流量(m3/s)50年20年10年12.62.62.628.93.93.4321.212.310.1454.843.135.35283.4141.7115.96500.4310.3253.97656.5380.4311.28850.0660.0480.09576.6453.5371.010405.9371.1303.611195.8154.0126.01247.037.030.3窟窿山水库100年一遇入库洪峰流量为1110 m3/s，根据各频率入库洪峰流量点绘P-频率曲线，通过频率曲线查得30年一遇入库洪峰流量为771 m3/s，5年一遇入库洪

13、峰流量为326 m3/s，详细计算过程见附录2。100年一遇洪水过程线按洪峰流量同倍比法对50年一遇洪水过程线进行缩放，缩放系数为1110/850=1.31；30年一遇洪水过程线按洪峰流量同倍比法对20年一遇洪水过程线进行缩放，缩放系数为771/660=1.17；5年一遇洪水过程线按洪峰流量同倍比法对10年一遇洪水过程线进行缩放，缩放系数为326/480=0.68。各频率设计洪水过程线见表2.4-2。表2.4-2 窟窿山水库各频率入库洪水过程线时间(h)不同重现期设计流量(m3/s)100年50年30年20年10年5年13.4 2.62.62.62.61.8211.6 8.94.63.93.4

14、2.3327.7 21.214.412.310.16.9471.6 54.850.343.135.3245370.1 283.4165.5141.7115.978.76653.5 500.4362.5310.3253.9172.47857.3 656.5444.4380.4311.2211.481110.0 8507716604803269753.0 576.6529.8453.537125210530.1 405.9433.5371.1303.6206.211255.7 195.8179.915412685.61261.4 4743.23730.320.6 水库下泄洪水（本节所采用的高程系统

15、为1956黄海高程）窟窿山水库调度运用方式为：主汛限水位采用566.36m，后汛限水位567.36m，水位超过566.36m时由溢洪道泄洪，20年一遇洪水限泄300m3/s，10年一遇洪水限泄294 m3/s。窟窿山水库水位容积泄量关系见表2.4-3。表2.4-3 窟窿山水库水位容积泄量关系水位（m）容积（万m3）面积（km2）泄量（m3/s）泄洪洞正常溢洪道非常溢洪道合计564.3600566.362830.856369.0369.035673390.0889104.68104.685684260.9302169.48169.485695250.9691243.920243.92570625

16、1.0108326.88108.45435.325717341.0493417.56373.93791.495728421.0819515.36721.301236.665739591.1435619.781217.491837.27573.510171.1745674.351482.862157.21水库已有的调蓄成果见表2.4-4。2.4-4 窟窿山水库不同重现期调洪成果表洪水重现期50年20年10年设计入库洪峰流量(m3/s)850660480出库洪峰流量(m3/s)695300294本次计算需对水库100年一遇、30年一遇、5年一遇下泄洪水进行计算，因此在已有的窟窿山水库水位容积泄量关

17、系的基础上对各频率洪水进行调洪演算，以推求各频率下泄洪水，水库起调水位为566.36m。经调洪演算，窟窿山水库100年一遇下泄洪水为915 m3/s，30年一遇下泄洪水为497 m3/s，5年一遇下泄洪水为174 m3/s。各频率调洪演算过程详见附录3。2.5 考虑水库调蓄后线路跨越断面设计洪水线路跨越断面洪水来源为水库下泄洪水与区间洪水的组合，参考河北省滦平县牤牛河县城下游段河道综合治理工程设计变更报告，水库下泄洪水与区间洪水错峰叠加，错峰时间间隔为1h，考虑到水库大坝距离线路跨越断面约14km，而区间洪水来源为短历时高强度的局部暴雨的山区洪水，且洪水特点为暴涨暴落，因此断面洪水组合以区间洪

18、水为主，水库下泄洪水为辅，即叠加计算时水库下泄洪峰比区间洪峰晚1h出现。洪水错峰叠加计算见表2.4-5。表2-7 线路跨越断面设计洪水成果表组合形式项目不同重现期设计洪峰流量(m3/s)100年50年30年 5年区间设计水库相应区间设计1134975839420窟窿山水库相应入流1110850771326窟窿山水库相应出库915695493174下泄洪峰前1h流量817361297157线路跨越断面1951133611365773 设计洪水位（1）水文断面为推求线路跨越牤牛河处100年一遇、30年一遇设计洪水位，在穿越附近布设垂直于河道的水文断面。实测河道横断面高程间距见表3-1。表3-1跨越

19、断面河道高程间距表序号起点距高程备注10465.672公路路面27.58461.75河岸外苗圃、耕地354.933461.7494111.309461.1365188.492461.126197.779460.4857229.481460.858左岸河堤8231.389461.7099233676461.1210234.21460.04211238.205458.659河槽12242.165458.59913242.165458.0414273.021457.74215279.882457.14116287.052457.05817296.8457.66518309.703457.427193

20、15.322461.397右岸河堤20324.059462.57221325.987462.33622329.672461.69523407.914461.245河岸外耕地24469.055461.19125583.01461.27326583.322462.311山脚下路面根据向滦平县水务局收资了解，线路跨越河段规划修建河堤和清理河道，在原有河堤的基础上加高加固河堤，河堤高度为5m，迎水面采用浆砌石仰斜式挡土墙，迎水坡边坡1：0.25，背水坡边坡1：0.15，顶宽1.1m，高5m，底宽1.85m；外边坡1：2。防洪堤内填筑砂砾石，其相对密度不小于0.65。浆砌石护岸采用M10水泥砂浆砌筑并勾

21、平缝，每10m设一道伸缩缝，采用2cm厚闭孔泡沫板填缝。在距河底底部0.3m处设PVC排水孔，排水孔间距3m，孔径10cm。排水孔在墙内侧端口采用土工布包裹砂砾料封堵。（2）河段比降根据现场实测资料，利用EHP软件计算河段比降为1.38，从结果来看，实测比降成果偏小，主要原因是冬季水面结冰导致水流流动不明显。而根据河北省滦平县牤牛河县城下游段河道综合治理工程设计变更报告，该河段位于王家北沟至国道112之间，设计河道纵坡为4.2，经过分析比较，采用设计河段比降J = 4.2。表3-2 河段比降计算表点号间距（m）高程（m）146.4457.511268.3457.359362457.279464

22、.2457.183550.3456.7026104456.6977122.8456.6868456.565平均坡降1.38（3）河道糙率线路跨越处河底为砂砾或卵石组成，底坡较均匀，床面基本平整，上游顺直,下游接缓弯, 水流较通畅，河道为单式河槽，据“天然河道糙率表”，主槽糙率取0.03。左、右岸均有简易堆石堤坝，右岸岸外为农田和苗圃，种植玉米和密集的树苗，左岸岸外为大片地势平坦的农田，种植玉米等农作物，并长有稀疏的小树。根据滩地糙率表，左岸滩地糙率取0.075，右岸堤外滩地糙率均取0.06。洪水主要在河道中行洪，当大洪水漫溢出堤岸，将在左右岸外的农田中形成堤外行洪。 （4）设计洪水位根据跨越水

23、文断面实测资料，采用曼宁公式推求水位流量关系曲线，详见附录三。根据水位流量关系，现状河道堤内过水能力1608m3/s（相应水位为右岸堤顶高程461.7m），小于100年一遇设计洪峰流量1335m3/s，当发生100年一遇洪水时，河道洪水将漫溢出岸，造成堤外行洪。通过布设的水文断面水位流量关系，各频率设计洪水及洪水位见表3-3。表3-3 不同重现期下的设计洪水位频率（%）设计洪峰流量（m3/s）设计洪水位（m）11951461.913.31136460.8820577459.73根据现场调查，常年洪水位高于河道护槽顶部约1m，实测常年洪水位为459.60m。附录1 牤牛河设计洪水水科院法计算过程

24、附录2 窟窿山水库调洪演算过程附录3 曼宁公式推求牤牛河水位流量关系附录1 牤牛河设计洪水水科院法计算过程计算原理根据电力工程水文勘测计算手册，水利电力部科学研究院推理公式法是以推理成因为基础，并结合了经验性的参数定量方法。分析中采用了各种概化假定条件，使公式带有一定的条件性。本法只适用于汇水面积在500km2以下的、全面产流条件下的全面和部分汇流条件。 设计洪峰流量计算公式为： 当tc时 当tc时 式中： Qm 设计洪峰流量（m3/s） 洪峰径流系数Sp 设计频率的雨力（mm/h），计算公式为： 汇流时间（h） =1时的汇流时间Tc 产流历时（h） 损失参数（mm/h）m 汇流参数n 暴雨衰

25、减指数F 流域面积（km2）L 河长（km）J 沿L的加权平均纵比降（km）1 推理公式计算1.1 流域基本参数 本流域面积(F)=202km2，河长(L)=14.306km，平均坡降(J)=9.599（计算见下表），（L/J1/3）=67，汇流参数m=1.2。表1.1-1 平均坡降计算表811.74401110.44501351.2460882.34701587.8480993.44901507.85002279.55101295.6530351.45401387.7550637.4600109.7650690平均坡降9.5991.2 暴雨与损失参数表1.2-1 暴雨与损失参数表频率(%)H

26、24(mm)(mm/h)N1N2N31182.510.10.670.750.692162.710.10.670.750.693.33148.710.10.670.750.692095.810.10.670.750.69注：H2424小时最大降水（mm） 损失参数（mm/h） N暴雨衰减指数1.3 计算结果表1.3-1 推理公式计算结果表频率(%)洪峰流量(m3/s)汇流时间(h)洪水历时(h)洪水总量(万m3)116002.59.91420213002.610.412203.3311102.710.81080205123.313601表1.3-2 5点法过程线表频率(%)点1点2点3点4点51

27、时间01.22.54.99.9流量0256160018102时间01.32.65.210.4流量0209130014703.33时间01.42.75.410.8流量01781110126020时间01.63.36.513流量081.851257.80时段流量图附录2 窟窿山水库调洪演算过程一 100年一遇来水调洪演算1 水库调洪计算1.1 水库基本参数 起调水位为:566.35m，水位(H)、库容(V)、泄流量(q)关系见下表。表1.1-1 水位-库容-泄流曲线表水位(H):m 库容(V):万m3 出流(q):m3/sHVq564.3600566.3628369.03567339104.685

28、68426169.48569525243.92570625435.32571734791.495728421236.665739591837.27573.510172157.211.2 来水过程表1.2-1 来水过程表时间(T):小时 来水流量(Q):m3/sNoTQ113.42211.63327.74471.655370.166653.577857.3881110.099753.01010530.11111255.7121261.41.3 调洪结果 出库流量过程见下表。表1.3-1 调洪出流过程表时间序列(小时)设计洪水(m3/s)出库流量(m3/s)库水位(m)库容(万m3)13.468.

29、7566.35282211.663.5566.20260327.759.8566.09245471.659.0566.072425370.177.8566.522976653.5176568.094357857.3396569.7960481110817571.067409753915571.2876410530.1696570.7370511255.7472570.106361261.4300569.2955413206568.4947414157567.8040915120567.233591693.2566.793211774.0566.452911864.8566.242661959.

30、4566.082432054.4565.942232149.8565.802042245.6565.681872341.8565.571712438.3565.471572535.1565.381442632.1565.291322729.4565.211212826.9565.141102924.7565.071013022.6565.0192.63120.7564.9684.83219.0564.9177.73317.4564.8671.23415.9564.8265.23514.6564.7859.73613.3564.7554.73712.2564.7150.13811.2564.68

31、45.93910.2564.6642.0409.38564.6338.5418.59564.6135.2427.87564.5932.3437.21564.5729.6446.60564.5527.1456.05564.5424.8465.54564.5222.7475.07564.5120.8484.65564.4919.0494.25564.4817.4503.90564.4716.0513.57564.4614.6523.27564.4513.4532.99564.4512.3542.74564.4411.2552.51564.4310.3562.30564.439.43572.1156

32、4.428.64581.93564.427.91591.77564.417.25601.62564.416.63611.48564.406.07621.36564.405.56631.24564.405.10641.14564.394.66651.04564.394.28660.955564.393.91最大泄量(m3/s)915最高水位(m)571.28相应库容(万m3)764二 30年一遇来水调洪演算1 水库调洪计算1.1 水库基本参数 起调水位为:566.35m，水位(H)、库容(V)、泄流量(q)关系见下表。表1.1-1 水位-库容-泄流曲线表水位(H):m 库容(V):万m3 出流(

33、q):m3/sHVq564.3600566.3628369.03567339104.68568426169.48569525243.92570625435.32571734791.495728421236.665739591837.27573.510172157.211.2 来水过程表1.2-1 来水过程表 时间(T):小时 来水流量(Q):m3/sNoTQ112.6224.63314.44450.355165.566362.577444.48877199529.81010433.51111179.9121243.21.3 调洪结果 出库流量过程见下表。表1.3-1 调洪出流过程表时间序列(小

34、时)设计洪水(m3/s)出库流量(m3/s)库水位(m)库容(万m3)12.668.7566.3528224.663.2566.19259314.458.7566.06241450.356.5566.002325165.560.8566.122496362.590.7566.753177444.4163567.894178771297569.285539529.8497570.1764410433.5486570.1464011179.9378569.705951243.2243568.9952413190568.2845314145567.6239315111567.093471687.25

35、66.693111769.2566.362831863.3566.192591958.0566.042382053.1565.902182148.6565.771992244.5565.651832340.8565.541672437.4565.441532534.2565.351402631.3565.271282728.7565.191182826.3565.121082924.1565.0698.73022.0565.0090.43120.2564.9582.83218.5564.9075.83316.9564.8569.53415.5564.8163.63514.2564.7758.3

36、3613.0564.7453.43711.9564.7148.93810.9564.6844.83910.0564.6541.0409.16564.6337.6418.39564.6034.4427.68564.5831.5437.04564.5628.8446.44564.5526.4455.90564.5324.2465.41564.5222.2474.95564.5020.3484.53564.4918.6494.15564.4817.0503.80564.4715.6513.48564.4614.3523.19564.4513.1532.92564.4412.0542.68564.44

37、11.0552.45564.4310.0562.24564.429.20572.05564.428.42581.88564.417.71591.72564.417.06601.58564.416.47611.45564.405.92621.32564.405.42631.21564.404.97641.11564.394.55651.02564.394.17660.932564.393.82最大泄量(m3/s)497最高水位(m)570.17相应库容(万m3)644二 5年一遇来水调洪演算1 水库调洪计算1.1 水库基本参数 起调水位为:566.35m，水位(H)、库容(V)、泄流量(q)关系

38、见下表。表1.1-1 水位-库容-泄流曲线表水位(H):m 库容(V):万m3 出流(q):m3/sHVq564.3600566.3628369.03567339104.68568426169.48569525243.92570625435.32571734791.495728421236.665739591837.27573.510172157.211.2 来水过程表1.2-1 来水过程表 时间(T):小时 来水流量(Q):m3/sNoTQ111.8222.3336.944245578.766172.477211.488326992521010206.2111185.6121220.61.3

39、 调洪结果 出库流量过程见下表。表1.3-1 调洪出流过程表时间序列(h)设计洪水(m3/s)出库流量(m3/s)库水位(m)库容(万m3)11.868.7566.3528222.363.1566.1925936.958.2566.0523842454.6565.94224578.754.3565.932236172.460.3566.112477211.474.7566.462928326116567.183549252157567.8140910206.2174568.064321185.6167567.974231220.6140567.5538713110567.083461486.4

40、566.673101568.9566.362821663.1566.192591757.8566.032371852.9565.892171948.5565.761992044.4565.651822140.7565.541672237.2565.441532334.1565.351402431.2565.261282528.6565.191172626.2565.121072724.0565.0698.42822.0565.0090.12920.1564.9482.53018.4564.8975.63116.9564.8569.23215.5564.8163.43314.2564.7758.

41、13413.0564.7453.23511.9564.7048.73610.9564.6844.6379.97564.6540.9389.13564.6237.4398.36564.6034.3407.66564.5831.4417.01564.5628.8426.42564.5526.3435.88564.5324.1445.39564.5222.1454.94564.5020.2464.52564.4918.5474.14564.4817.0483.79564.4715.5493.47564.4614.2503.18564.4513.0512.91564.4411.9522.67564.4410.9532.44564.4310.0542.24564.429.17552.05564.428.40561.88564.417.70571.72564.417.05581.57564.416.45591.44564.405.91601.32564.405.41611.21564.394.95621.11564.394.53631.01564.394.15