工程水文及水利计算%20.docx

工程水文及水利计算 武 汉 大 学 水利水电 学院 港口海岸及治河工程 专业课 程 设 计 说 明 书题 目： 工程水文及水力计算 设计期 限： 2012年1月 设 计 人： 班级及学号： 2009301580016 港航2班 指导 教 师： 艾学山 2012年 1月 29日目录工程水文及水利计算课程设计2水库设计及水电站水利水能计算2一设计任务3二 流域与资料概括3 1防洪限制水位的选取4 2出力系数4 3预留水头损失4 4正常蓄水位选取4三 设计年径流量及其年内分配的推求4 1.设计代表年径流量的计算4 2.设计代表年径流的年内分配4四.水库死水位的选择5 1. 绘制水库水位容积曲线和水电站下游水位流量关系曲线5 2. 根据泥沙资料计算水库的淤积体积和水库相应的淤积高程7 3.由水轮机实际情况可确定水库死水位。7五水库正常蓄水位的确定8 1.保证出力的计算8 2.确定水库正常蓄水位9六.装机容量与多年平均发电量计算10 1.保证出力的计算10 2.中水年、丰水年的水能计算10 3.出力历时曲线11 4.装机容量和平均发电量的计算11七设计洪水的推求13八.水库洪水特征水位的计算和坝顶高程的确定 21 1.防洪高水位的计算21 2.设计洪水位的计算24 3.校核洪水位的计算27 4.坝顶高程计算30一设计任务在某河流上，拟建一水库，因而要进行水库水文分析及水利水能计算，其具体任务是：1. 设计保证率下的年径流及其年内分配2. 水库死水位的选择3. 水库正常蓄水位的确定4. 装机容量与多年平均发电量计算5. 设计洪水的推求6. 水库洪水特征水位的计算和坝顶高程的确定二 流域与资料概括西苕溪为太湖流域一大水系，流域面积为2260km2，发源于浙江省安吉县天目山，干流全长150km，上游坡陡流急，安城以下堰塘遍布，河道曲折，排泄不畅，易遭洪涝灾害，又因流域拦蓄工程较少，灌溉水源不足，易受旱灾。根据解放后二十多年的统计，仅安吉县因洪涝旱灾每年平均损失稻谷1500万斤，严重的196163年，连续三年洪水损失稻谷9300万斤，冲毁耕地万余亩。该水库为根治西苕溪流域水旱灾害骨干工程之一，控制西苕溪主要支流西溪，坝址以上流域面积328km2。流域内气候温和、湿润、多年平均雨量孝丰站为1450mm，国民经济以农、林业为主，流域内大部为山区，小部为丘陵，平地较少。流域水系及测站分布见图1。水库以防洪为主，结合发电、灌溉、航运及水产，是一座综合利用水库。图1 西苕溪流域图1防洪限制水位的选取 起调水位（防洪限制水位）为78.4 m。2出力系数 A=7.53预留水头损失HF0.5 m。4正常蓄水位选取 初步确定正常蓄水位为80+(学号后两位/100)米；若通过水能计算后能满足保证出力要求就作为确定的正常蓄水位，否则，以0.2米为单位进行调整。三 设计年径流量及其年内分配的推求1.设计代表年径流量的计算对观测的流量资料进行年径流量频率计算，求出丰、中、枯三个代表年（频率分别为85%、50、15）的年径流量，选择典型年，计算缩放倍比，成果表1。表1 设计年径流及典型年径流量 单位：m3/s代表年设计频率设计年径流量典型年典型年径流量缩放倍比枯水年P=855.5619735.351.039中水年P=507.5019577.111.055丰水年P=1510.3719679.951.0422.设计代表年径流的年内分配在坝址下游1公里处设有潜渔水文站，自1954年开始有观测的流量资料。通过频率计算，得各设计频率的设计年径流量，选择典型年，计算缩放倍比如题目中给出。通过典型缩放得丰（P=15）、中(P=50)、枯(P=85)设计年径流过程，见表1。年平均流量为：5.56+7.50+10.373=7.81（m3）表1 潜渔站设计年径流过程 月份 设计枯水年 设计中水年 设计丰水年典型年Q设计年Q典型年Q设计年Q典型年Q设计年Q314.0014.566.476.8211.1011.58418.6019.346.937.3112.8013.35519.5020.2820.8021.9423.0023.9961.311.368.248.6927.6028.7975.215.4221.8023.003.974.148.06.065.505.801.241.299.04.044.324.5623.2024.2010.41.432.812.967.657.98111.701.771.421.501.361.4212.75.783.603.80.44.461.70.731.101.16.84.8821.881.962.332.466.116.3764.1666.7385.3290119.31124.45四.水库死水位的选择1. 绘制水库水位容积曲线和水电站下游水位流量关系曲线根据表2和表3，绘出水库水位容积曲线和水电站下游水位流量关系曲线，并将原始数据处理之后，插入AUTO CAD中，进行精确画图，以便进行后面数据的查阅，提高所得数据的准确度。表2 水位容积曲线水位（m）485052556065容积（106 m3）00.10.62.38.018.0水位（m）707580818283容积（106 m3）35.760.394.4102.8111.3120.0水位（m）848586878889容积（106 m3）129.0138.6148.3158.8170.0181.5水位（m）9091容积（106 m3）194.5207.0水位Z/mV/106m3水库水位库容曲线表3 水库下游水位流量关系曲线水位（m）46.046.246.446.646.847.047.2流量（m3/s）01.03.67.311.916.220.5水位/m流量/m3s-1 下游水位流量曲线2. 根据泥沙资料计算水库的淤积体积和水库相应的淤积高程根据实测泥沙资料得多年平均含沙量kg/m3，泥沙干容重kg/m3，泥沙沉积率m90，孔隙率p0.3。表1可得，多年平均径流总量为：=3652436007.81=246.30106m3,则悬移质年淤积体积为：V1=0.237246.30106m30.9（1-0.3）1650=0.0455106m3因为推移质与悬移质淤积量之比值15，水库的设计运行年限为50年，所以，总的淤积体积为：。可查得淤积高程为：55.4002m。由安全值为2+0.16=2.16m，所以死水位Z死55.4002+2.16=57.5602m3.由水轮机实际情况可确定水库死水位。水轮机最小水头为16m，则可计算得到最大发电流量为： Q=12509.8116=7.964m3/s；查水电站下游水位流量曲线，可确定下游水位为46.63m，则可确定上游水位为：46.63+16=62.63m。综合上述情况，可得水库的最低死水位为62.63m，查水库水位容积曲线可得水库的死库容为13.26106 m3。五水库正常蓄水位的确定根据本地区的兴利、发电等综合利用要求，保证出力不低于800千瓦，发电保证率为85，灌溉及航运任务不大，均可利用发电尾水得到满足，因此，初步确定正常蓄水位为：80+0.16=80.16m .V兴=V总-V死=95.5237-13.26=82.2637106m3由此可知水库的兴利库容为82.2637106m3.1.保证出力的计算根据本地区的兴利、发电等综合利用要求，保证出力不低于800千瓦，发电保证率为85，选取设计枯水年进行计算：假设3月到5月为蓄水期，计算蓄水期的发电流量为：Q调=W蓄-V兴T蓄=14.56+19.34+20.28-32.513=7.32m3/s；由于经计算得到的蓄水期的发电流量Q调Q来，因此假设的供水期不正确。重新假设6月，8月到次年2月为供水期，计算供水期的发电流量为：Q调=W供+V兴T供=1.36+0.06+0.04+0.43+1.77+0.78+0.73+1.96+32.228=4.92m3/s此时计算得到的供水期的发电流量满足Q调Q来，因此假设的供水期正确。因此可根据每月的天然来水和发电流量确定每月水库的水量，从而根据水库水位流量关系确定上游水位，再根据发电流量确定下游水位，利用出力公式即可计算得到每月的出力。出力公式： N=AQH-HF其中出力系数A7.5，Q为发电流量，预留水头损失HF0.5 m。表4 水库枯水年出力计算表月份天然来水（m3/s）发电用水（m3/s）余水量不足水量月末水库蓄水（106 m3）平均蓄水（106 m3）月平均水位（m）下游水位（m）水头（m）出力（kw）流量（m3/s）水量（106 m3）流量（m3/s）水量（106 m3）213.26314.567.327.24 19.03 32.29 22.77 66.35 46.60 19.75 1056.74 419.347.3212.02 31.59 63.88 48.08 72.52 46.60 25.92 1395.36 520.287.3212.96 34.06 97.93 80.90 78.02 46.60 31.42 1697.57 61.364.923.569.36 88.58 93.26 79.83 46.47 33.36 1212.62 75.425.4288.58 88.58 79.15 46.50 32.65 1306.75 80.064.924.8612.77 75.81 82.19 78.21 46.47 31.74 1152.76 90.044.914.8712.80 63.01 69.41 76.34 46.47 29.87 1081.38 100.434.924.4911.80 51.21 57.11 74.35 46.47 27.88 1010.37 111.774.923.158.28 42.93 47.07 72.31 46.47 25.84 935.08 120.784.924.1410.88 32.05 37.49 70.36 46.47 23.89 863.23 10.734.924.1911.01 21.04 26.54 67.41 46.47 20.94 754.37 21.964.922.967.78 13.26 17.15 64.57 46.47 18.10 649.61 2.确定水库正常蓄水位由于保证出力为供水期的平均出力，计算得：Np=1212.62+1152.76+1081.38 +1010.37 +935.08+863.23+754.37+649.618=956.38kw，根据本地区的兴利、发电等综合利用要求，保证出力不低于800千瓦，通过水能计算后能满足保证出力要求。因此水库的正常蓄水位为：h正=80.63m六.装机容量与多年平均发电量计算1.保证出力的计算先对丰、中、枯三个代表年以月为时段进行水能计算，计算出各月的水流出力。设计枯水年的水能已经在上述过程中进行计算，取其供水期的平均出力为保证出力，即为957kw。2.中水年、丰水年的水能计算应用相同的方法进行丰水年和中水年的水能计算，如表6、7所示：表65 水库丰水年出力计算表月份天然来水（m3/s）发电用水（m3/s）余水量不足水量月末水库蓄水（106 m3）平均蓄水（106 m3）月平均水位（m）下游水位（m）水头（m）出力（kw）流量（m3/s）水量（106 m3）流量（m3/s）水量（106 m3）213.26311.5811.5813.26 13.26 62.63 46.79 15.84 1332.28 413.3513.3513.26 13.26 62.63 46.87 15.76 1527.91 523.9914.929.07 23.84 37.10 25.18 67.03 46.94 20.09 2191.86 628.7914.9213.87 36.45 73.55 55.32 73.99 46.94 27.05 2970.73 74.147.83.669.62 63.93 68.74 76.24 46.62 29.62 1703.35 81.297.86.5117.11 46.82 55.37 74.00 46.62 27.38 1572.41 924.214.929.28 24.39 71.21 59.01 74.74 46.94 27.80 3054.70 107.987.9871.21 71.21 76.60 46.63 29.97 1763.74 111.427.86.3816.77 54.44 62.82 75.37 46.62 28.75 1652.63 120.467.797.3319.26 35.18 44.81 71.85 46.62 25.23 1444.94 10.887.796.9118.16 17.02 26.10 67.29 46.62 20.67 1178.29 26.377.81.433.76 13.26 15.14 63.57 46.62 16.95 962.30 表6 水库中水年出力计算表月份天然来水（m3/s)发电用水（m3/s余水量不足水量月末水库蓄水（106 m3）平均蓄水（106 m3）月平均水位（m）下游水位（m）水头（m）出力流量（m3/s水量（106 m3）流量（m3/s水量（106 m3）213.2636.826.8213.26 13.26 62.63 46.57 16.06 795.69 47.317.3113.26 13.26 62.63 46.60 16.03 851.41 521.947.1414.80 38.89 52.15 32.71 69.15 46.59 22.56 1181.49 68.697.131.56 4.10 56.25 54.20 73.76 46.59 27.17 1426.19 7237.1415.86 41.68 97.93 77.09 77.46 46.59 30.87 1626.37 85.87.781.985.20 92.73 95.33 80.11 46.62 33.49 1924.92 94.567.783.228.46 84.27 88.50 79.13 46.62 32.51 1867.96 102.967.784.8212.67 71.60 77.94 77.59 46.62 30.96 1777.57 111.57.786.2816.50 55.10 63.35 75.45 46.62 28.83 1652.78 123.87.783.9810.46 44.64 49.87 72.88 46.62 26.26 1502.97 11.167.786.6217.40 27.24 35.94 70.05 46.62 23.43 1337.78 22.467.785.3213.98 13.26 20.25 66.13 46.62 19.50 1108.85 3.出力历时曲线将上述求得的丰水年、中水年和枯水年共36个月的出力值按从大到小的顺序排列并分别求出相应的频率，作出出力历时曲线如图3。4.装机容量和平均发电量的计算应用装机年利用小时法确定多年水库的装机容量和平均发电量。具体方法是：利用水电站的多年平均发电量E年除以装机容量N装得到机组多年平均的年利用小时数。因此，推算装机容量，可利用公式：N装=E年h装由于水电站的E年和N装有关，因此利用上式计算装机容量时需采用试算法。 图3 出力历时曲线先假定装机容量大于36个月中的最大出力，即3054.7kw，此时计算水电站的多年平均发电量如表8。根据上述计算结果分别计算得到枯水年、中水年和丰水年的年发电量为：E枯=9574569.5（kwh）E中=12449408.5（kwh）E丰=15589237.9（kwh）则水电站的多年平均发电量为：E年=E枯+E中+E丰3=9574569.5+12449408.5+15589237.93=12537637.5（kwh）因此可求得水电站的装机容量为：N装=E年h装=12537637.53440=3630kw由于水电站的最大出力为3054.7kw，所求的装机容量大于此值，因此符合要求，即水电站的装机容量为：N装=3630kw，多年平均发电量为E年=12537637.5kwh. 表7 水电站发电量计算表月份枯水年中水年丰水年出力（kw）发电量（kwh）出力（kw）发电量（kwh）出力（kw）发电量（kwh）31056.74 771418.4 795.69 580851.2 1332.28 972563.7 41395.36 1018615.6 851.41 621527.4 1527.91 1115372.5 51697.57 1239229.7 1181.49 862484.6 2191.86 1600058.7 61212.62 885211.9 1426.19 1041118.5 2970.73 2168629.5 71306.75 953928.5 1626.37 1187253.0 1703.35 1243446.1 81152.76 841512.8 1924.92 1405190.9 1572.41 1147855.7 91081.38 789407.2 1867.96 1363610.9 3054.70 2229933.6 101010.37 737568.9 1777.57 1297628.1 1763.74 1287529.6 11935.08 682605.6 1652.78 1206532.4 1652.63 1206420.5 12863.23 630160.2 1502.97 1097169.1 1444.94 1054802.6 1754.37 550692.3 1337.78 976582.3 1178.29 860149.3 2649.61 474218.3 1108.85 809460.3 962.30 702476.3 七设计洪水的推求 确定洪本水库为大（）型水库，工程等别为等，永久性水工建筑级别为2级。下游防洪标准为5%，设计标准为1，校核标准为0.1%，需要推求5%、1、0.1%设计洪水过程线。 按年最大值选样方法在实测资料中选取最大洪峰流量及各历时洪量，根据洪水特性和防洪计算的要求，确定设计历时为7天，控制历时为1天和3天，因而可得洪峰和各历时的洪量系列。表8 潜渔站洪峰及定时段洪量统计表年份洪峰Qm（m3/s）24小时洪量W1d（106m3）三天洪量W3d（106m3）七天洪量W7d（106m3）195470227.9458.4019552848.1713.30195674829.8036.00195740222.8037.1952.4619582008.7215.85221319.8032.90196047815.7020.8033.20196165952.5079.1088.20196258543.7049.2053.101963116055.6086.6095.90196440914.3231.7040.70196551015.6224.4027.0019662329.5014.0025.40196724411.8219.0028.0019681679.9018.4035.50196938720.9032.8048.40197030517.2031.9035.60197150023.4031.8035.3019721085.3410.2012.23197348419.8742.85197428716.1639.05197516611.5822.0519761198.2919.9519772387.6120.45由于7天洪量只有19571972年，用相关分析方法延长插补：采用三天洪量为相关变量，取用线性回归的的方法建立W7d倚W3d 的回归方程进行延长插补系列。由表9及回归方程的定义及公式可知： Y=0.9967X+9.064125表9 回归方程计算表年份3D X（106m3）7D Y（106m3）KXKYKX-1KY-1(KX-1)2(KY-1)2(KX-1)(KY-1)195737.1952.461.138311.2602250.138310.2602250.019130.0677170.035992195815.8522.150.4851360.5321-0.51486-0.46790.2650850.218930.240905195919.832.90.6060370.790343-0.39396-0.209660.1552070.0439560.082597196020.833.20.6366450.79755-0.36335-0.202450.1320270.0409860.073561196179.188.22.4210892.1187921.4210891.1187922.0194941.2516951.589902196249.253.11.5059111.2755990.5059110.2755990.2559460.0759550.139429196386.695.92.6506492.3037661.6506491.3037662.724641.6998052.152059196431.740.70.9702720.977719-0.02973-0.022280.0008840.0004960.000662196524.4270.7468340.64861-0.25317-0.351390.0640930.1234750.0889619661425.40.4285110.610174-0.57149-0.389830.3265990.1519650.222781196719280.5815510.672632-0.41845-0.327370.17510.107170.136987196818.435.50.5631860.852802-0.43681-0.14720.1908060.0216670.064298196932.848.41.0039411.1626930.0039410.1626931.55E-050.0264690.000641197031.935.60.9763940.855204-0.02361-0.14480.0005570.0209660.003418197131.835.30.9733330.847997-0.02667-0.1520.0007110.0231050.004053197210.212.230.3122010.293796-0.6878-0.70620.4730670.4987240.485726由表格可得：x=X(Kx-1)2n-1=22.003 106m3 y=Y(Ky-1)2n-1=22.4765 106m3 r=0.975702由显著性检验当=0.01时，r=0.4869 。由于 r r 所以为总体正相关。由公式计算的结果可知空缺处的结果依次为：67.27141；22.32024；44.94533；51.77272；47.98526；31.04136；28.94829；29.44664。 根据调查1922年9月1日在坝址附近发生一场大洪水，推算得潜渔站洪峰流量为1350m3/s。这场洪水是发生后至今最大的一次洪水。缺测年份内，没有大于1160m3/s的洪水发生。对于洪峰加入历史特大洪水采用独立样本法计算经验频率，计算结果见表10。表10 洪峰经验频率计算表调查期或实测期系列年数洪水序列洪水年份洪峰Qm（m3/s）经验频率P统一样本法n(实测）N（调查）m（实测）M（调查）调查期N156（1922-1977）1192213500.0175442196311600.035088实测期n24（1954-1977）1196311600.035088219567480.075293319547020.115497419616590.155702519625850.195907619655100.236111719715000.276316819734840.316521919604780.3567251019644090.396931119574020.4371351219693870.4773391319703050.5175441419742870.5577491519552840.5979531619672440.6381581719772380.6783631819592370.7185671919662320.7587722019582000.7989772119681670.8391812219751660.8793862319761190.9195912419721080.959795对于1日、3日和7日洪量频率计算，按照公式进行计算。结果列于表11中。表11 1日、3日和7日洪量频率计算表年份24小时洪量W1d（106m3）经验频率P年份三天洪量W3d（106m3）经验频率P年份七天洪量W7d（106m3）经验频率P196355.60.04196386.60.04196395.90.04196152.50.08196179.10.08196188.20.08196243.70.12195458.40.12195467.271410.12195629.80.16196249.20.16196253.10.16195427.940.20197342.850.20195752.460.20197123.40.24197439.050.24197351.772720.24195722.80.28195737.190.28196948.40.28196920.90.321956360.32197447.985260.32197319.870.36196932.80.36195644.945330.36197017.20.40197031.90.40196440.70.40197416.160.44197131.80.44197035.60.44196015.70.48196431.70.48196835.50.48196515.620.52196524.40.52197135.30.52196414.320.56197522.050.56196033.20.56196711.820.60196020.80.60195932.90.60197511.580.64197720.450.64197531.041360.64195911.130.68197619.950.68197729.446640.6819689.90.72195919.80.72197628.948290.7219669.50.761967190.761967280.7619588.720.80196818.40.801965270.8019768.290.84195815.850.84196625.40.8419558.170.881966140.88195522.320240.8819777.610.92195513.30.92195822.150.9219725.340.96197210.20.96197212.230.96由以上两表可得频率曲线如下： 利用绘图软件画出频率曲线，并查询可得频率为5%，1%，0.1% 的洪峰值和各历时段的洪量值，列入表12.表12 相应设计频率的洪峰值和历时洪量值频率（%）洪峰Qm（m3/s）24小时洪量W1d（106m3）三天洪量W3d（106m3）七天洪量W7d（106m3）0.12278.52110.91167.35173.0211566.3674.11113.15117.2651067.2348.6777.5080.68由最大洪峰Qm=1160（m3/s）,由最大一日流量W1d=56.79106 m3 由最大三日流量W3d=89.5284106 m3 由最大七日流量W7d=92.93106 m3 求得放大倍比系数K列入表12.表13 典型年洪水过程的放大倍比系数频率（%）洪峰Qm（m3/s）24小时洪量W1d（106m3）三天洪量W3d（106m3）七天洪量W7d（106m3）0.11.9611.5321.7231.66911.3531.3111.1931.20950.9210.8620.8210.882根据所给的典型洪水过程及放大倍数，得到设计洪水过程所得结果列入表14.表14 设计洪水过程表时段（t1h）流量（m3/s）流量P=0.1%（m3/s）流量P=1%（m3/s）流量P=5%（m3/s）时段（t1h）流量（m3/s）流量P=0.1%（m3/s）流量P=1%（m3/s）流量P=5%（m3/s）0610.13227.25585.289650286493.0068341.055234.83461711.82098.46516.171251295508.521351.7875242.22452711.82098.46516.171252290499.902345.825238.1193813.50969.67447.052853263453.3594313.6275215.949341525.330518.139513.22454277477.4926330.3225227.444752338.840127.813920.276855265456.807316.0125217.591563457.415841.116229.974456252434.3976300.51206.917274372.614151.277535.307357245422.331292.1625201.169585796.255967.972546.802758235405.093280.2375192.9585973125.837487.052559.940359224386.1312267.12183.926410100172.38119.2582.1160212365.4456252.81174.073211130224.094155.025106.74361205353.379244.4625168.325512171294.7698203.9175140.408162192330.9696228.96157.651213230396.474274.275188.85363180310.284214.65147.79814260448.188310.05213.48664161277.5318191.9925132.197115317546.4446378.0225260.288765155267.189184.8375127.270516335577.473399.4875275.068566145249.951172.9125119.059517362624.0156431.685297.238267130224.094155.025106.74318373728.469486.765320.668168118203.4084140.71596.889819350683.55456.75300.89569118203.4084140.71596.889820330644.49430.65283.70170103177.5514122.827584.573321315615.195411.075270.80557197167.2086115.672579.646722305595.665398.025262.20857292158.5896109.7175.541223293572.229382.365251.89217385146.523101.362569.793524293572.229382.365251.89217482141.351697.78567.330225445869.08