注册结构工程师程序说明c

1、C-3两种供水保证的水库径流调节计算程序作者(水电部勘测设计院)勘测设计院)校核 潘东海(水电部能本程序可求出满足两种供水保证率的水库兴利及限制的数值范围，并对在以上范围内选取的几组兴利及限制进行径流调节计算和供水破坏深度的比较。通过以上计算成果的比较，可帮助选取水库的兴利及限制，作为水库设计的依据。二、计算方法概要本程序采取分别给兴利赋值从输入的允许最大兴利及限制赋值，进行试算，求其两种供水的保证率。的开始，逐次减少，直至计算出的保证率小于允许保证率为止。对于每一给定的，采用优选法找出能达到允许供水保证率的限制的上、下限。再在以上上下限范围内选择几组及计算供水破坏深度。通过以上选择和比较，最

2、后选出合适的水库兴利及限制。三、计算公式水库径流调节计算的主要公式为：式中：V-月末水库蓄量；V0-上月末水库蓄量； R-月径流量；A 及 B 分别为两种供水的月水量；E-水库蒸发渗漏朋损失量。若 V0 表示保证率以年统计；K0 表示保证率以月统计；G0 表示需打印输入数据(以便校对使用) G0 表示不需打印输入数据L0 表示需继续运算；L0 表示结束运算；Y(99)-存放供水破坏年份数组； M(99)-存放供水破坏月份数组； V(99)-存放水库蓄量的数组； YE(11)-存放年份的数组； MO(11)-存放月份的数组； R(11)-存放月径流量的数组；A(11)-存放高保证率供水月水量的数

3、组； B(11)-存放低保证率供水月水量的数组； E(11)-存放水库蒸发渗漏月损失量的数组。数据文件顺序A E年号R，BR，B年号 R，B高保证率供水月水量(12 个数)水库蒸发渗漏月损失量起始月 终止月月径流量和低保证率供水月水量月径流量和低保证率供水月水量同上五、操作说明将数据文件与运行程序放在同一目录下，启动程序，按照提问输入数据即可计算完毕。六、算例某水库有 1958 年 6 月至 1978 年 5 月共 20 年的逐月径流系列，其逐月工业用水(A)、灌溉用水(B)及水库损失量(E)如表 1(为万立米)。工业用水的保证率为 90%，灌溉用水的保证率为 70%(保证率均以年统计)，试用

4、时历法进行径流调节计算并选取合适的兴利及灌溉限制。启动后输入 N20、VM7000、DX50、PA0.9、PB=0.7、K=1、MO=6、G=1。经过较长时间的运算后打印出 VX 及 VY 及的上下限范围及供水破坏的次数(见表 2)。在以上范围内不同的及值其供水保证率是相同的，但其供水破坏深度是不同的。一般情况下 VX 大则供水破坏浓度浅；在 VX 相同的情况下，VY 愈大，则工业供水破坏愈浅，灌溉用水破坏愈深。本次选四组 VX 及 VY，其分别为：(一)VX=6900、VY=1420；(二)VX =6800、VY=1700；(三)VX =6800、VY=1320；(四)VX =6500、VY

5、=1420。依次输入进行计算并打印出计算结果(见表 3)。通过对照比较选出合适的 VX 及 VY，作为水库设计的依据。如果保证率是以月统计，工业用水及灌溉用水的保证率分别为 98.5%及 90%。启动后输入相应的数据、(其中 K=0、A=0.985、PB=0.90，其余数据同上)。数据文件为 C-3.:182,182,182,182,182,311,311,311,311,182,182,18257,57,57,57,57,57,57,57,57,57,57,571958,393,551,2244,66,2483,319,187,204,106,302,759,0,130,086,0,454,

6、0,967,0,1949,269,578,01959,1440,522,1117,572,271,528,333,319,307,269,144,269564,0,237,0,209,0,725,0,711,0,528,01960,787,935,780,368,213,545,2248,0,191,269,225,269225,0,157,0,38,269,211,0,133,341,120,01961,132,478,109,660,123,610,132,0,48,269,513,092,0,59,0,61,269,41,269,49,341,73,01962,122,428,112,4

7、07,2230,269,405,0,674,0,2053,0405,0,55,0,68,269,142,0,229,0,4503,01963,1118,791,128,335,8424,66,2197,0,497,269,0,00,0,0,0,134,0,260,0,4303,0,4225,01964,387,633,554,66,1274,358,2788,0,3686,0,494,269236,0,61,0,114,0,327,0,366,0,140,2691965,333,863,7550,83,1300,204,72,66,24,341,13,00,0,0,0,6,341,290,0,

8、62,66,85,2701966,93,615,289,506,111,676,37,270,18,0,24,27063,0,19,0,89,0,251,0,858,0,273,5231967,277,682,4627,137,277,374,944,0,365,341,1738,0352,0,135,0,102,0,117,0,104,270,453,01968,160,1094,2013,121,835,440,2735,0,2306,0,142,270121,0,157,0,852,0,754,0,1947,0,535,01969,443,776,320,154,2158,314,487

9、7,0,571,341,67,37483,0,70,0,198,0,348,341,614,0,475,01970,917,567,271,121,82,578,414,0,355,270,142,27090,0,98,0,121,0,238,0,248,270,227,01971,2976,424,1419,390,2127,154,183,270,218,0,667,095,0,98,0,258,302,821,0,481,302,86,01972,455,699,4186,121,174,852,134,77,133,0,167,022,0,75,0,118,0,74,302,898,0

10、,1634,01973,45,567,1695,121,212,578,143,0,94,270,0,2700,0,0,0,14,0,102,0,121,270,1004,01974,0,475,67,406,2848,380,127,0,1181,275,162,275133,0,147,0,266,275,38,69,1470,0,131,01975,122,791,654,335,19600,66,640,0,643,269,179,0442,0,102,0,336,0,126,0,171,0,123,01976,57,615,774,506,611,676,332,270,105,0,

11、75,27046,0,27,0,23,0,14,0,192,0,387,5231977,238,567,1249,121,584,578,319,0,453,270,603,270123,0,104,0,363,0,150,0,65,270,38,0计算结果文件为 C-3.OUT:*两种供水保证率的水库径流调节书*输入数据:径流系列年数 N = 20允许最大的兴利VM =7000.0兴利递减的幅度 DX = 50高保证率(例如工业用水)的供水保证率 PA = 0.900低保证率(例如农业用水)的供水保证率 PB = 0.700K0 表示保证率以年统计,K=0 保证率以月统计. K = 1每年起

12、始月份 (如日M0=1,水文年一般为 6) M0 = 6计算结果:年月 月径流量高保证率供水 低保证率供水 蒸发渗漏月损失量5515719587224419581958101110675918231130205757057195961440182522573192695757195919591959101112307182564311000000575757575757935368545057575757196019601960196091011122248182225225311311269002690341047866061002690057575757575757575757575757

13、575757575757575719612383111961115133111962196219622346141493111821822692693414284071962196219621011126742053405182311311000575757575757571963268311269001963545031963891084242197497182182026957571963575757575719641031100019644430342253871826336635800269057575757575757575757575757196419641964196419648

14、910111212742788368649423618218218231131100026986319657875501300182204571966196612063113110341057575719664562851822705767657057575757005236825757570341000002700109457575757575757575757196719671967101112365173835218231131119681968196819689101112273523061821820057571213110000575757571969419475351969196

15、989215848771821823140575705757575757575757570341005671215781970882182027057571970103551820270042439057575757571971197167297614191821821971197110112186671823110005757575757575757573020302699770005757575757571972197210111331671823110027027057575757571973110311197410311057575702701974510040671821824754

16、0657571974719741974197410111211811822755713331100275575757197534381470131197519751975101112643179442182311311269000000057575757575757576762705757575757575757052356702702700000270575757575757575757高保证率供水破坏月数1977197719771011124536031231823113111978465182高保证率供水破坏年数低保证率供水破坏年数低保证率供水破坏月数兴利限制6423兴利限制VX= VY

17、=69001420破坏年数 AY=破坏月数 AM=破坏年数 BY=高保证率供水高保证率供水低保证率供水152高保证率供水年供水保证率 PAY= 0.905低保证率供水年供水保证率 PAM= 0.714高保证率供水月供水保证率 PBY= 0.988低保证率供水破坏月数 BM=水库蓄量1541209101182165553841114201199855550201-7818低保证率供水月供水保证率 PBM= 0.921年1958196219621962196219621962196619661966196619671967月62345679101112124565VX= VY=6316657031

18、41868001700破坏年数 AY=破坏月数 AM=破坏年数 BY=破坏月数 BM=水库蓄量15414161109911721555438311170014981277933628279196719671978兴利限制高保证率供水高保证率供水低保证率供水低保证率供水15121高保证率供水年供水保证率 PAY= 0.905低保证率供水年供水保证率 PAM= 0.714高保证率供水月供水保证率 PBY= 0.992低保证率供水月供水保证率 PBM= 0.909年19581962196219621962196219621962196619661966196619661967月61234567891011121126316657031700149968001320破坏年数 AY=破坏月数 AM=破坏年数 BY=破坏月数 BM=水库蓄量15411099117215554383111320109975545010