洪水调节设计模板-带试算c语言程序

1、洪水调节课程设计任务书1、 设计目的1. 洪水调节目的：定量地找出入库洪水、下泄洪水、拦蓄洪水的库容、水库水位的变化、泄洪建筑物型式和尺寸间的关系，为确定水库的有关参数和泄洪建筑型式选择、尺寸确定提供依据；2. 掌握列表试算法和半图解法的基本原理、方法、步骤及各自的特点；3. 了解工程设计所需洪水调节计算要解决的课题；培养学生分析问题、解决问题的能力。2、 设计基本资料1. 某水利枢纽工程以发电为主，兼有防洪、供水、养殖等综合效益，电站装机为5000KW年发电量1372X104kwh,水库库容亿m3。挡水建筑物为混凝土面板坝，最大坝高。溢洪道堰顶高程，采用2孔8rnK6m（宽X高）的弧形门控制

2、。水库正常蓄水位。电站发电引用流量为10m3/s。2. 本工程采用2孔溢洪道泄洪。在洪水期间洪水来临时，先用闸门控制下泄流量q并使其等于洪水来水量Q，使水库水位保持在防洪限制水位不变；当洪水来水量Q继续增大时，闸门逐渐打开；当闸门达到全开后，就不再用闸门控制，下泄流量q随水库水位z的升高而增大，流态为自由流态，情况与无闸门控制一样。3. 上游防洪限制水位（注：X=”号最后1位/10,即），下游无防汛要求。3、 设计任务及步骤分别对设计洪水标准、校核洪水标准，按照上述拟定的泄洪建筑物的类型、尺寸和水库运用方式，分别采用列表试算法和半图解法推求水库下泄流量过程，以及相应的库容、水位变化过程。具体步

3、骤：1 .根据工程规模和建筑物的等级，确定相应的洪水标准；2 .用列表试算法进行调洪演算：根据已知水库水位容积关系曲线VZ和泄洪建筑物方案，用水力学公式求出下泄流量与库容关系曲线qZ,并将VZ,qZ绘制在图上；决定开始计算时刻和此时的qi、Vi,然后列表试算，试算过程中，对每一时段的q2、V2进行试算；将计算结果绘成曲线：CHt、qt在一张图上，Zt曲线绘制在下方。3 .用半图解法进行调洪计算:绘制三条曲线：fizL_q,f,zL9,qt2t2进行图解计算，将结果列成表格。4 .比较分析试算法和平图解法调洪计算的成果。4、 时间安排和要求1 .设计时间为1周;2 .成果要求：设计说明书编写要求

4、条理清楚、附图绘制标准； 列表试算法要求采用手工计算，熟悉过程后可编程计算，如采用编程计算需提供程序清单及相应说明； 设计成果请独立完成，如有雷同则二者皆取消成绩，另提交成果时抽查质询。5、 参考书3 .水利水电工程等级划分及洪水标准(SL252-2000)4 .水利水能规划附录：一、堰顶溢流公式3/2qnbm2gH0式中：q通过溢流孔口的下泄流量，m3/s；n溢流孔孔口数；b溢流孔单孔净宽，m;g重力加速度，s2;闸墩侧收缩系数，与墩头形式有关，初步计算可假设为;m流量系数，与堰顶形式有关，可查表，本工程取；H0堰顶水头，m.设计洪水过程时刻(h)、一，3Q实测(m/s).,一、3各频率Q(

5、m/s)%1%2%5%0503529201136296196162121231268052443235733491300727602524496020001220104073951670230013901130806612902100129010907757919175011901010698854311808537065419402895647505387103248174834003271129470943736227012264606398326243132345493482892161420447729425119515191440283230176161774142632191621

6、7164385245204151181503512241871391913732020417012520123286183152113211102571711421062210224015412796239722614411989249021213511183水位-库容曲线和库容表540530520)/.510位500见490小480470某水利枢纽水位与库容曲线图一,H1一，一一fr46045001000200030004000容积（万立方米）5000600070008000库容表高程（m）450460470480490500505库容（104m3）018高程（m）5105155205255

7、30535540库容（104m?）6670四、工程分等分级规范和洪水标准Z1J水利水电工程分等指标，工程等则工程嫌检水阵总璋容阴洪浩潜灌氨供水发电你安城他及工ir隹业的重:共性保护农田CIO'亩）泊涝而视南）据双向限（101供水前世/案性裳机摩室I大（烟10特丽弟整3sjOO学四门>150忡冽重要力120IK型IC-1.0*找5。5110口1505。桌安120-301空L.qQ.IC中等4。一155。舟中等3。5W小门.型力.ICCt.61一般30当-JJ5-0.*5-限5-1小刚也口口.CXjIY5<3<0<5L一<1表3+2一山瓯r丘隆区水利水电工程永

8、包性水工逢填物洪水标准重现期L年】域目m-L建馍物级别1£345他核设计土行取lo(X)fc51X1可他配人演火CPMK)成101)00-51X)(SWLHT0UMJ制卬心100-59Ldoti5'U-3011门口小tXJJO-20中j。24N)浪凝土坝.望勘石城53】。呼EUVOIJQU1QL>g-50(J5OUN。"2UU15调洪计算成果表项目频率一一一设计洪水校核洪水列表试算法最大泄量（m/s）sm3/s水库最高水位（m1半图解法最大泄量（m/s）ss水库最高水位（m）洪水调节演算过程洪水标准的确定1 .工程等别的确定：由设计对象的基本资料可知，该水利枢

9、纽工程以发电为主，兼有防洪、供水、养殖等其他综合效益，电站装机为5000kW水库库容。若仅由装机容量5000kW的指标，根据“水利水电工程分等指标”，可将工程等别定为V;若仅以水库总库容为指标，则可将工程等别定为田。综合两种指标，取等级最高的田等为工程最终等别。2 .洪水标准的确定：该水利工程的挡水建筑物为混凝土面板坝，由已确定的为田等的工程等别，根据“山区，丘陵区水利水电工程永久性水工建筑物洪水标准”，可查得，该工程设计洪水标准为10050年，校核标准为1000500年，不妨取设计标准为100年，校核洪水标准为1000年。试算法洪水调节计算3/21 .计算并绘制水库的q=f(V)关系曲线：应

10、用式qnbm：2gH0,根据不同水库水位计算H与q,再由HV关系曲线查得V,并计算于下表，绘制q=f(V)关系曲线图如下。2 .10324817483400327112947094373622701226460639832624313234549348289216142044772942511951519144028323017616177414263219162171643852452041511815035122418713919137320204170125201232861831521132111025717114210622102240154127962397226144119892

11、490212135111833.）米Z440某水利枢纽水位与库容曲线图（VZ）2000400060008000库容V （万立方米）55054053052051050049048047046045010000位水4.水库q=f（V）关系曲线计算表水位/m堰顶水头/mq溢(m3/s)q电(m3/s)总泄流量q库容V/10A4m351900101052011052121052231052341052451052561052671052781052891052910105301110531121013115321310533141053415105351610201366705361710537181

12、05381910260253920105402110库水位Z（m）)s/3ca.量泄总水库水位谷量关系曲线00 45524529库水位Z（R）5345390058007500560055容库0053519下泄流量与水位关系曲线（qZ量泄总ooo ooo 5 0 53 3 2o o o o0 52 1o o50水库q=f (V)关系曲线水库q=f (V)关系曲 线00 5345005500650075008500库容V(万m35.确定调洪的起始条件：起调水位也是防洪限制水位，Z=。相应库容x104m3在洪水期间洪水来临时，先用闸门控制下泄流量q并使其等于洪水来水量Q,使水库水位保持在防洪限制水位

13、不变；当洪水来水量Q继续增大时,闸门逐渐打开；当闸门达到全开后，就不再用闸门控制，下泄流量q随水库水位z的升高而增大，流态为自由流态，情况与无闸门控制一样。由公式：3qQnbmJ2gH0210=V29.816.2.m/s得调洪开始时的下泄流量为m3/s。所以在第一时段，以闸门控制入库流量等于下泄流量；以后时段闸门全开不再控制，下泄流量由试算计算。6.列表试算泄流量q,本过程米用C语言编程试算。基本原理：根据水库容积曲线V=f(Z)和堰顶溢流公式q=f(H),得出蓄泄方程q=f(V)。联立水量平衡方程V2Vi-t(Q2Qiq2q1)qf(V)2可得q=f(V)=g(q),即q=g(q)。编程公式

14、的主要过程a)已知的电站发电引用流量为10n3/s,结合堰顶溢流公式，得出下泄流q=nbm.,2gH02 +100000000000000000000000000000000000000000(1)b)水位高程Z与堰顶水头H的关系。基本材料可知溢洪道堰顶高程为519mWJH=Z-519mc)水库容积曲线V=f(Z)的近似化。根据该设计的蓄泄情况，水位高程的变化范围在525m535nt问，又由于水库容积曲线在水位高程属525m535r间的变化率较小,： 化、近似计算。由水位一库容表(HI),易算出 H=g(V)=V 3591 .68 V182 .02V 3226 .48 V215 .22联立(1

15、)、(2)式得3q = nb m、2gg2(V) 10。,d)将(3)式与水量平衡方程联立。/日，t _行V2 V1 一 (Q2 Q1 q233q = nb mv'2gg2(V)e) C语言程序源代码如下：9方便计算，故可将其分段直线化以简V=f (Z)及上式 H=Z 519rn 可得 V=f4683 .8,5593 .95593 .9,6670 .0 00 000 00(2)OOOOOOOOOOOOOOOOOOOOOOOOOOOOO (3)。)10。 (4)#include<>#include<>voidmain()floatV1,V2,Q1,Q2,q1,q2

16、,q3,t=;printf("V1=");scanf("%f",&V1);printf("Q1=");scanf("%f",&Q1);printf("Q2=");scanf("%f",&Q2);printf("q1=");scanf("%f",&q1);printf("q2=");scanf("%f",&q2);printf("nn");

17、loop:V2=V1+(Q1+Q2-q2-q1)*t/2;if(V2>=&&V2<=q3=(pow(/,)*+10;elseif(V2>=&&V2<=q3=(pow(/,)*+10;if(fabs(q3-q2)>q2=q3;gotoloop;printf("q2=%fn",q3);printf("V2=%fnnn",V2);7. 对设计洪水计算时段平均入库流量和时段入库水量。 将洪水过程表中P=1%勺洪水过程线划分计算时段，初选时段At=1h=3600填入下表第一栏，表中第二栏为按计算时段摘录

18、的入库洪水流量，计算的时段平均入库流量和时段入库水量分别填入第三栏和第四栏。泄流量的计算见第五，六，七栏。从表中第一，五栏可绘制下泄流量过程线。第一，十栏可绘制水位过程线； 为了枯水期能保证兴利部门的用水需求，当水位再次下降到调洪水位时，又需要用闸门控制下泄流量q并使其等于洪水来水量Q使水库水位保持在防洪限制水位不变。见第15时段q=f(V)的程序计算截图； 绘制Qt与qt曲线，如图所示。最大下泄流量qmax=m3/s发生在t=8h时，正好是qt曲线与Qt曲线的交点，即为所求的最大下泄流量；推求设计调洪库容V设和设计洪水位Z设。qmax=对应的库容和水位分别为万m3和，减去堰顶以下的库容万m3

19、即可得V设=万m3,Z设=mo第2时段试算法程序计算截图第时段试算法的程序计算截图回SLL 'E:cg=SifDebugCppl.exeB1=54*651=1021.592=1880.4375ql=996.53q2=996q2=99&.781677U2=S408.92187SPfrssanykevtocontinue设计洪水调节计算表时间t(h)入库洪水流量Q（m3/s）时段平均入库流量Q（平均）（m3/s）时段入库水量Q（平均）t（万m3）下泄流量q（m3/s）时段平均下泄流量（m3/s）时段下泄水量q（平均）t（万m3时段内水库存水量变化V（万m3）水库存水量V（万m3水库

20、水位Z（m）123456789100353501196196252436037274122051390130561290134071190124088539647750270104835651143746075012398133483731429432115283162632731724525424518224224019204214204214020183183021171177171177022154154023144149144149024135135Qt、qt曲线图0061入库流量下泄流量004002OOOO08641)0OO224681012141618202224时间t(h)水库水

21、位时间关系曲线（Zt）2422O28161412O18642O530987622225555m幺位水时间t(h)水旱存亡量之臣夫娄史或行Wt<i;)8.对校核洪水计算时段平均入库流量和时段入库水量。将洪水过程表中P=%勺洪水过程线划分计算时段，初选时段At=1h=3600填入下表第一栏，表中第二栏为按计算时段摘录的入库洪水流量，计算的时段平均入库流量和时段入库水量分别填入第三栏和第四栏。泄流量的计算见第五，六，七栏。从表中第一，五栏可绘制下泄流量过程线。第一，十栏可绘制水位过程线。 为了枯水期能保证兴利部门的用水需求，当水位再次下降到调洪水位时，又需要用闸门控制下泄流量q并使其等于洪水来

22、水量Q使水库水位保持在防洪限制水位不变。见第20时段q=f(V)的程序计算截图。 绘制Qt与qt曲线，如图所示。最大下泄流量qmax=s发生在t=8h时，正好是qt曲线与Qt曲线的交点，即为所求的最大下泄流量。 推求校核调洪库容V校和设计洪水位Z校。qmax=对应的库容和水位分别为万m3和，减去堰顶以下的库容万m3即可得V校二万m3,Z校=m。第2时段试算法程序计算截图I旧拿语言空档存宿'Debug'CppL®h三2回第|ui=4652.27601=296O2=S801=29£2-4502=448.806024U2=4654.876953Pressanyke

23、ytocontinue第时段试算法程序计算截图"J"E；ci§ugCpp1-回隋W1-6115.4301=1571.87502=1554.0625Ql=1548.392=15482-1548.521362112=6UE.593262Pressanykeytocentinlie校核洪水调节计算表时间t(h)入库洪水流量Q（m3/s）时段平均入库流量Q（平均）（m3/s）时段入库水量Q（平均）t（万m3下泄流量q（m3/s）时段平均下泄流量（m3/s）时段下泄水量q（平均）t（万m3时段内水库存水量变化V（万m3）水库存水量V（万m3水库水位Z（m）123456789

24、100501735017301296296268048831300990420001650594523002150774530621002200792717501925693154814658118014659895108178561335117097631260613549144775131544016414427173851835136819320202863032125722240240282232262332262330242122192122190某水库校核调洪计算表Qt、qt过程曲线2500r024681012141618202224时间t(h)水位时间关系曲线（Zt）024681

25、012141618202224时间h(t)水库存水量时间关系曲线水库存水量时间关系曲线时间t(h)三、半图解法洪水调节计算（以设计洪水标准进行调洪演算为例子）Vq1 .计算并绘制q楙辅助线。计算时段取tiho计算过程见下表。利用表中第五，七栏的相应数值绘制的辅助线如图。qlq辅助曲线计算表t2q=f（V/At+q/2）辅助曲线计算表水库水位Z（m）总库容V（万m3堰顶以上库容（万m3）V/At(m3/s)q(m3/s)q/2(m3/s)V/At+q/2(m3/s)V/A+q/2ZZ（D）q=f(V/ At+q/2)辅助曲线3500.00 3000.00& 2500.00sm 2000.

26、00(1500.00q 1000.00500.00 0.00q=f(V/ At+q/2)辅助曲0.005000.0010000.0015000.00V/ A t+q/2Z(m)Z（m）Nv/A2 .调洪计算求qt过程和库水位过程：由于作辅助线时t需取固定值，且Vqqf（2）是由蓄泄曲线q=f（V）转换而来的，故该法只适用于自由泄流（无闸或闸门全开）和t固定的情况。当有闸门控制泄流时，应按控制的流量调洪。所以对于第一时段Q<Q艮时，起调水位也是防洪限制水位，Z=。相应库容x104m3o在洪水期间洪水来临时，先用闸门控制下泄流量q并使其等于洪水来水量Q,使水库水位保持在防洪限制水位不变，即q

27、=Q可直接求出（V!q!）=m3/s和（"曳）=m3/s。从第二时段开始闸门全开，Q=196t2t2m3/s,Q3=524m3/s,q2=196m3/s,将（V2的）=m3/s代入3=t2t2Viqi”3q3Vq、一、（X-）+Qq-求得（匚V=m/s,由式qf（匚；）行，q=m/st2t2t23同法依次类推，可求出其他时段的泄量。其成果如表第四，六栏所示。3 .绘制Qt与qt过程线以及q乙过程线，求qmax。利用表中第一，二，四栏数值，可绘出Qt与qt过程线。取Qt与qt两曲线的交点的纵坐标数值，作为qmax=s,t=8h时。利用表中的库水位Z与泄流能力q,可绘制Zq关系曲线，如图所示半图解法调洪计算表(设计P=1%时间t(h)入库流量Q(m3/s)平均入库流量Q'(m3/s)V/t+q/2(m3/s)q(m3/s)水库水位Z(m)035351196196252436037274