洪水调节设计试算法和半图解法带试算C语言程序

Documentserialnumber【UU89WT-UU98YT-UU8CB-UUUT-UUT108】Documentserialnumber【UU89WT-UU98YT-UU8CB-UUUT-UUT108】洪水调节设计试算法和半图解法带试算C语言程序《洪水调节课程设计》任务书设计目的洪水调节目的：定量地找出入库洪水、下泄洪水、拦蓄洪水的库容、水库水位的变化、泄洪建筑物型式和尺寸间的关系，为确定水库的有关参数和泄洪建筑型式选择、尺寸确定提供依据；掌握列表试算法和半图解法的基本原理、方法、步骤及各自的特点；了解工程设计所需洪水调节计算要解决的课题；培养学生分析问题、解决问题的能力。设计基本资料某水利枢纽工程以发电为主，兼有防洪、供水、养殖等综合效益，电站装机为5000KW，年发电量1372×104kw·h，水库库容亿m3。挡水建筑物为混凝土面板坝，最大坝高。溢洪道堰顶高程，采用2孔8m×6m（宽×高）的弧形门控制。水库正常蓄水位。电站发电引用流量为10m3/s。本工程采用2孔溢洪道泄洪。在洪水期间洪水来临时，先用闸门控制下泄流量q并使其等于洪水来水量Q，使水库水位保持在防洪限制水位不变；当洪水来水量Q继续增大时，闸门逐渐打开；当闸门达到全开后，就不再用闸门控制，下泄流量q随水库水位z的升高而增大，流态为自由流态，情况与无闸门控制一样。上游防洪限制水位（注：X=+学号最后1位/10，即），下游无防汛要求。设计任务及步骤分别对设计洪水标准、校核洪水标准，按照上述拟定的泄洪建筑物的类型、尺寸和水库运用方式，分别采用列表试算法和半图解法推求水库下泄流量过程，以及相应的库容、水位变化过程。具体步骤：根据工程规模和建筑物的等级，确定相应的洪水标准；用列表试算法进行调洪演算：根据已知水库水位容积关系曲线V～Z和泄洪建筑物方案，用水力学公式求出下泄流量与库容关系曲线q～Z，并将V～Z，q～Z绘制在图上；决定开始计算时刻和此时的q1、V1，然后列表试算，试算过程中，对每一时段的q2、V2进行试算；将计算结果绘成曲线：Q～t、q～t在一张图上，Z～t曲线绘制在下方。用半图解法进行调洪计算：绘制三条曲线：，，；进行图解计算，将结果列成表格。比较分析试算法和半图解法调洪计算的成果。时间安排和要求设计时间为1周；成果要求：设计说明书编写要求条理清楚、附图绘制标准；列表试算法要求采用手工计算，熟悉过程后可编程计算，如采用编程计算需提供程序清单及相应说明；设计成果请独立完成，如有雷同则二者皆取消成绩，另提交成果时抽查质询。参考书《水利水电工程等级划分及洪水标准》(SL252-2000)《水利水能规划》附录：堰顶溢流公式式中：q——通过溢流孔口的下泄流量，m3/s；n——溢流孔孔口数；b——溢流孔单孔净宽，m；g——重力加速度，s2；ε——闸墩侧收缩系数，与墩头形式有关，初步计算可假设为；m——流量系数，与堰顶形式有关，可查表，本工程取；H0——堰顶水头，m。设计洪水过程时刻（h）Q实测（m3/s）各频率Q（m3/s）%1%2%5%05035292011362961961621212312680524432357334913007276025244960200012201040739516702300139011308066129021001290109077579191750119010106988543118085370654194028956475053871032481748340032711294709437362270122646063983262431323454934828921614204477294251195151914402832301761617741426321916217164385245204151181503512241871391913732020417012520123286183152113211102571711421062210224015412796239722614411989249021213511183水位－库容曲线和库容表库容表高程（m）450460470480490500505库容（104m3）018高程（m）510515520525530535540库容（104m3）6670工程分等分级规范和洪水标准调洪计算成果表频率项目设计洪水校核洪水列表试算法最大泄量（m3/s）sm3/s水库最高水位（m）半图解法最大泄量（m3/s）ss水库最高水位（m）洪水调节演算过程洪水标准的确定工程等别的确定：由设计对象的基本资料可知，该水利枢纽工程以发电为主，兼有防洪、供水、养殖等其他综合效益，电站装机为5000kW，水库库容。若仅由装机容量5000kW为指标，根据“水利水电工程分等指标”，可将工程等别定为Ⅴ；若仅以水库总库容为指标，则可将工程等别定为Ⅲ。综合两种指标，取等级最高的Ⅲ等为工程最终等别。洪水标准的确定：该水利工程的挡水建筑物为混凝土面板坝，由已确定的为Ⅲ等的工程等别，根据“山区，丘陵区水利水电工程永久性水工建筑物洪水标准”，可查得，该工程设计洪水标准为100~50年，校核标准为1000~500年，不妨取设计标准为100年，校核洪水标准为1000年。试算法洪水调节计算计算并绘制水库的q=f(V)关系曲线：应用式，根据不同水库水位计算H与q,再由H~V关系曲线查得V，并计算于下表，绘制q=f(V)关系曲线图如下。1032481748340032711294709437362270122646063983262431323454934828921614204477294251195151914402832301761617741426321916217164385245204151181503512241871391913732020417012520123286183152113211102571711421062210224015412796239722614411989249021213511183水库q=f(V)关系曲线计算表水位/m堰顶水头/mq溢(m3/s)q电(m3/s)总泄流量q库容V/10^4m35190010105201105212105223105234105245105256105267105278105289105291010530111053112101311532131053314105341510535161020136670536171053718105381910260253920105402110确定调洪的起始条件：起调水位也是防洪限制水位，Z=。相应库容×104m3。在洪水期间洪水来临时，先用闸门控制下泄流量q并使其等于洪水来水量Q，使水库水位保持在防洪限制水位不变；当洪水来水量Q继续增大时，闸门逐渐打开；当闸门达到全开后，就不再用闸门控制，下泄流量q随水库水位z的升高而增大，流态为自由流态，情况与无闸门控制一样。由公式：=m3/s得调洪开始时的下泄流量为m3/s。所以在第一时段，以闸门控制入库流量等于下泄流量；以后时段闸门全开不再控制，下泄流量由试算计算。列表试算泄流量q，本过程采用C语言编程试算。基本原理：根据水库容积曲线V=f（Z）和堰顶溢流公式q=f（H），得出蓄泄方程q=f（V）。联立水量平衡方程可得q=f(V)=g（q），即q=g（q）。编程公式的主要过程已知的电站发电引用流量为10m3/s，结合堰顶溢流公式，得出下泄流q=nbm+10。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。（1）水位高程Z与堰顶水头H的关系。基本材料可知溢洪道堰顶高程为519m则H=Z-519m；水库容积曲线V=f（Z）的近似化。根据该设计的蓄泄情况，水位高程的变化范围在525m~535m之间，又由于水库容积曲线在水位高程属525m~535m之间的变化率较小，为方便计算，故可将其分段直线化以简化、近似计算。由水位—库容表V=f（Z）及上式H=Z519m，可得V=f（H），易算出H=g（V）=。。。。。。。（2）联立（1）、（2）式得。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。（3）d）将（3）式与水量平衡方程联立。得。。。。。。。。。。。。。。。。。。。。。。。。。。。（4）e）C语言程序源代码如下：#include<>#include<>voidmain(){ floatV1,V2,Q1,Q2,q1,q2,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("\n\n");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=%f\n",q3); printf("V2=%f\n\n\n",V2);}对设计洪水计算时段平均入库流量和时段入库水量。将洪水过程表中P=1%的洪水过程线划分计算时段，初选时段Δt=1h=3600填入下表第一栏，表中第二栏为按计算时段摘录的入库洪水流量，计算的时段平均入库流量和时段入库水量分别填入第三栏和第四栏。泄流量的计算见第五，六，七栏。从表中第一，五栏可绘制下泄流量过程线。第一，十栏可绘制水位过程线；为了枯水期能保证兴利部门的用水需求，当水位再次下降到调洪水位时，又需要用闸门控制下泄流量q并使其等于洪水来水量Q，使水库水位保持在防洪限制水位不变。见第15时段q=f（V）的程序计算截图；绘制Q~t与q~t曲线，如图所示。最大下泄流量=m3/s发生在t=8h时，正好是q~t曲线与Q~t曲线的交点，即为所求的最大下泄流量；推求设计调洪库容和设计洪水位。=对应的库容和水位分别为万m3和，减去堰顶以下的库容万m3即可得=万m3，=m。第2时段试算法程序计算截图第时段试算法的程序计算截图设计洪水调节计算表时间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)1234567891003535011961962524360372741220513901305612901340711901240885396477502701048356511437460750123981334837314294321152831626327317245254245182242240192042142042140201831830211711771711770221541540231441491441490241351350对校核洪水计算时段平均入库流量和时段入库水量。将洪水过程表中P=%的洪水过程线划分计算时段，初选时段Δt=1h=3600填入下表第一栏，表中第二栏为按计算时段摘录的入库洪水流量，计算的时段平均入库流量和时段入库水量分别填入第三栏和第四栏。泄流量的计算见第五，六，七栏。从表中第一，五栏可绘制下泄流量过程线。第一，十栏可绘制水位过程线。为了枯水期能保证兴利部门的用水需求，当水位再次下降到调洪水位时，又需要用闸门控制下泄流量q并使其等于洪水来水量Q，使水库水位保持在防洪限制水位不变。见第20时段q=f（V）的程序计算截图。绘制Q~t与q~t曲线，如图所示。最大下泄流量=s发生在t=8h时，正好是q~t曲线与Q~t曲线的交点，即为所求的最大下泄流量。推求校核调洪库容和设计洪水位。=对应的库容和水位分别为万m3和，减去堰顶以下的库容万m3即可得=万m3，=m。第2时段试算法程序计算截图第时段试算法程序计算截图校核洪水调节计算表时间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)123456789100501735017301296296268048831300990420001650594523002150774530621002200792717501925693154814658118014659895108178561335117097631260613549144775131544016414427173851835136819320202863032125722240240282232262332262330242122192122190某水库校核调洪计算表半图解法洪水调节计算（以设计洪水标准进行调洪演算为例子）计算并绘制辅助线。计算时段取。计算过程见下表。利用表中第五，七栏的相应数值绘制的辅助线如图。辅助曲线计算表q=f(V/Δt+q/2)辅助曲线计算表水库水位Z（m）总库容V（万m3）堰顶以上库容（万m3）V/Δt(m3/s)q(m3/s)q/2(m3/s)V/Δt+q/2(m3/s)调洪计算求q~t过程和库水位过程：由于作辅助线时需取固定值，且是由蓄泄曲线q=f（V）转换而来的，故该法只适用于自由泄流（无闸或闸门全开）和固定的情况。当有闸门控制泄流时，应按控制的流量调洪。所以对于第一时段Q<Q限时，起调水位也是防洪限制水位，Z=。相应库容×104m3。在洪水期间洪水来临时，先用闸门控制下泄流量q并使其等于洪水来水量Q，使水库水位保持在防洪限制水位不变，即q=Q；可直接求出=m3/s和=m3/s。从第二时段开始闸门全开，Q2=196m3/s，=524m3/s，=196m3/s，将=m3/s代入=（）+，求得=m3/s，由式得，=m3/s。同法依次类推，可求出其他时段的泄量。其成果如表第四，六栏所示。绘制Q~t与q~t过程线以及q~Z，过程线，求。利用表中第一，二，四栏数值，可绘出Q~t与q~t过程线。取Q~t与q~t两曲线的交点的纵坐标数值，作为=s，t=8h时。利用表中的库水位Z与泄流能力q，可绘制Z~q关系曲线，如图所示。半图解法调洪计算表（设计P=1%）时间t(h)入库流量Q（m3/s)平均入库流量Q'(m3/s)V/Δt+q/2(m3/s)q(m3/s)水库水位Z（m）035351196196252436037274122051390130561290