注冊(cè)結(jié)構(gòu)工程師程序說(shuō)明c

1、C-3兩種供水保證的水庫(kù)徑流調(diào)節(jié)計(jì)算程序作者(水電部勘測(cè)設(shè)計(jì)院)勘測(cè)設(shè)計(jì)院)校核 潘東海(水電部能本程序可求出滿足兩種供水保證率的水庫(kù)興利及限制的數(shù)值范圍，并對(duì)在以上范圍內(nèi)選取的幾組興利及限制進(jìn)行徑流調(diào)節(jié)計(jì)算和供水破壞深度的比較。通過(guò)以上計(jì)算成果的比較，可幫助選取水庫(kù)的興利及限制，作為水庫(kù)設(shè)計(jì)的依據(jù)。二、計(jì)算方法概要本程序采取分別給興利賦值從輸入的允許最大興利及限制賦值，進(jìn)行試算，求其兩種供水的保證率。的開(kāi)始，逐次減少，直至計(jì)算出的保證率小于允許保證率為止。對(duì)于每一給定的，采用優(yōu)選法找出能達(dá)到允許供水保證率的限制的上、下限。再在以上上下限范圍內(nèi)選擇幾組及計(jì)算供水破壞深度。通過(guò)以上選擇和比較，最

2、后選出合適的水庫(kù)興利及限制。三、計(jì)算公式水庫(kù)徑流調(diào)節(jié)計(jì)算的主要公式為：式中：V-月末水庫(kù)蓄量；V0-上月末水庫(kù)蓄量； R-月徑流量；A 及 B 分別為兩種供水的月水量；E-水庫(kù)蒸發(fā)滲漏朋損失量。若 V0 表示保證率以年統(tǒng)計(jì)；K0 表示保證率以月統(tǒng)計(jì)；G0 表示需打印輸入數(shù)據(jù)(以便校對(duì)使用) G0 表示不需打印輸入數(shù)據(jù)L0 表示需繼續(xù)運(yùn)算；L0 表示結(jié)束運(yùn)算；Y(99)-存放供水破壞年份數(shù)組； M(99)-存放供水破壞月份數(shù)組； V(99)-存放水庫(kù)蓄量的數(shù)組； YE(11)-存放年份的數(shù)組； MO(11)-存放月份的數(shù)組； R(11)-存放月徑流量的數(shù)組；A(11)-存放高保證率供水月水量的數(shù)

3、組； B(11)-存放低保證率供水月水量的數(shù)組； E(11)-存放水庫(kù)蒸發(fā)滲漏月?lián)p失量的數(shù)組。數(shù)據(jù)文件順序A E年號(hào)R，BR，B年號(hào) R，B高保證率供水月水量(12 個(gè)數(shù))水庫(kù)蒸發(fā)滲漏月?lián)p失量起始月 終止月月徑流量和低保證率供水月水量月徑流量和低保證率供水月水量同上五、操作說(shuō)明將數(shù)據(jù)文件與運(yùn)行程序放在同一目錄下，啟動(dòng)程序，按照提問(wèn)輸入數(shù)據(jù)即可計(jì)算完畢。六、算例某水庫(kù)有 1958 年 6 月至 1978 年 5 月共 20 年的逐月徑流系列，其逐月工業(yè)用水(A)、灌溉用水(B)及水庫(kù)損失量(E)如表 1(為萬(wàn)立米)。工業(yè)用水的保證率為 90%，灌溉用水的保證率為 70%(保證率均以年統(tǒng)計(jì))，試用

4、時(shí)歷法進(jìn)行徑流調(diào)節(jié)計(jì)算并選取合適的興利及灌溉限制。啟動(dòng)后輸入 N20、VM7000、DX50、PA0.9、PB=0.7、K=1、MO=6、G=1。經(jīng)過(guò)較長(zhǎng)時(shí)間的運(yùn)算后打印出 VX 及 VY 及的上下限范圍及供水破壞的次數(shù)(見(jiàn)表 2)。在以上范圍內(nèi)不同的及值其供水保證率是相同的，但其供水破壞深度是不同的。一般情況下 VX 大則供水破壞濃度淺；在 VX 相同的情況下，VY 愈大，則工業(yè)供水破壞愈淺，灌溉用水破壞愈深。本次選四組 VX 及 VY，其分別為：(一)VX=6900、VY=1420；(二)VX =6800、VY=1700；(三)VX =6800、VY=1320；(四)VX =6500、VY

5、=1420。依次輸入進(jìn)行計(jì)算并打印出計(jì)算結(jié)果(見(jiàn)表 3)。通過(guò)對(duì)照比較選出合適的 VX 及 VY，作為水庫(kù)設(shè)計(jì)的依據(jù)。如果保證率是以月統(tǒng)計(jì)，工業(yè)用水及灌溉用水的保證率分別為 98.5%及 90%。啟動(dòng)后輸入相應(yīng)的數(shù)據(jù)、(其中 K=0、A=0.985、PB=0.90，其余數(shù)據(jù)同上)。數(shù)據(jù)文件為 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計(jì)算結(jié)果文件為 C-3.OUT:*兩種供水保證率的水庫(kù)徑流調(diào)節(jié)書*輸入數(shù)據(jù):徑流系列年數(shù) N = 20允許最大的興利VM =7000.0興利遞減的幅度 DX = 50高保證率(例如工業(yè)用水)的供水保證率 PA = 0.900低保證率(例如農(nóng)業(yè)用水)的供水保證率 PB = 0.700K0 表示保證率以年統(tǒng)計(jì),K=0 保證率以月統(tǒng)計(jì). K = 1每年起

12、始月份 (如日M0=1,水文年一般為 6) M0 = 6計(jì)算結(jié)果:年月 月徑流量高保證率供水 低保證率供水 蒸發(fā)滲漏月?lián)p失量5515719587224419581958101110675918231130205757057195961440182522573192695757195919591959101112307182564311000000575757575757935368545057575757196019601960196091011122248182225225311311269002690341047866061002690057575757575757575757575757

13、575757575757575719612383111961115133111962196219622346141493111821822692693414284071962196219621011126742053405182311311000575757575757571963268311269001963545031963891084242197497182182026957571963575757575719641031100019644430342253871826336635800269057575757575757575757575757196419641964196419648

14、910111212742788368649423618218218231131100026986319657875501300182204571966196612063113110341057575719664562851822705767657057575757005236825757570341000002700109457575757575757575757196719671967101112365173835218231131119681968196819689101112273523061821820057571213110000575757571969419475351969196

15、989215848771821823140575705757575757575757570341005671215781970882182027057571970103551820270042439057575757571971197167297614191821821971197110112186671823110005757575757575757573020302699770005757575757571972197210111331671823110027027057575757571973110311197410311057575702701974510040671821824754

16、0657571974719741974197410111211811822755713331100275575757197534381470131197519751975101112643179442182311311269000000057575757575757576762705757575757575757052356702702700000270575757575757575757高保證率供水破壞月數(shù)1977197719771011124536031231823113111978465182高保證率供水破壞年數(shù)低保證率供水破壞年數(shù)低保證率供水破壞月數(shù)興利限制6423興利限制VX= VY

17、=69001420破壞年數(shù) AY=破壞月數(shù) AM=破壞年數(shù) BY=高保證率供水高保證率供水低保證率供水152高保證率供水年供水保證率 PAY= 0.905低保證率供水年供水保證率 PAM= 0.714高保證率供水月供水保證率 PBY= 0.988低保證率供水破壞月數(shù) BM=水庫(kù)蓄量1541209101182165553841114201199855550201-7818低保證率供水月供水保證率 PBM= 0.921年1958196219621962196219621962196619661966196619671967月62345679101112124565VX= VY=6316657031

18、41868001700破壞年數(shù) AY=破壞月數(shù) AM=破壞年數(shù) BY=破壞月數(shù) BM=水庫(kù)蓄量15414161109911721555438311170014981277933628279196719671978興利限制高保證率供水高保證率供水低保證率供水低保證率供水15121高保證率供水年供水保證率 PAY= 0.905低保證率供水年供水保證率 PAM= 0.714高保證率供水月供水保證率 PBY= 0.992低保證率供水月供水保證率 PBM= 0.909年19581962196219621962196219621962196619661966196619661967月61234567891011121126316657031700149968001320破壞年數(shù) AY=破壞月數(shù) AM=破壞年數(shù) BY=破壞月數(shù) BM=水庫(kù)蓄量15411099117215554383111320109975545010