(2)按无偏估值公式计算统计参数
这里我们直接利用所编好程序嵌入实测径流数据得出偏差系数Cv=0.24, 偏态系数Cs=0.36;多年平均径流量由资料给出Q(平
均)=381.85.
(3)选配理论频率曲线
由上步求出的Q(平均),和Cv,Cs值先假设Cv=0.2,Cs=2Cv=0.4由程序可得出其对应的φp值,由Xp=Q(平均)×(1+Cv*φp)由已知年平均径流过程推得对应Xp,并绘制理论频率曲线;绘制完成后观察理论频率曲线发现其与经验点距相比较平缓上部靠下,下部靠上根据偏差系数Cv值对理论频率曲线的影响Cv增大会使曲线变陡,固第二步假设加大Cv值;
取Cv=0.26,Cs=2Cv=0.52,绘制理论频率曲线,绘制完成后较第一次假设拟合程度高,但还有偏差理论频率曲线较经验点距中间部分靠右,由偏态系数对理论频率曲线的影响,Cs越大曲线中部越向左偏,且上端越陡,下端越缓则应加大Cs值;经过调制最终得Cv=0.26,Cs=0.56;
(4)有所得理论频率曲线查得对应设计频率的年径流值:Q5%=559.56,Q50%=372.52,Q95%=236.5。计算过程见下表: 表(2) 设计年径流理论频率曲线试算表
均值X=∑Xi/n 计算值 381.85 0.24 一 381.85 二 381.85 三 381.85 取值 变差系数Cv 0.20 0.26 0.26 取值 偏态系数Cs 0.36 0.40 0.52 0.56 取值 P-Ⅲ理论曲线 频率% 至P=50%水平距离 X 0.4 Cs φp 0.6 一 Cv Cs 二 Cv Cs 三 0.2 0.4 0.26 0.52 φp Kp=φp*Cv+1 Xp=Kp*X φp Kp=φp*Cv+1 Xp=Kp*X φp 0.01 0.10 0.20 0.33 0.50 1.00 2.00 5.00 10.00 20.00 50.00 75.00 90.00 95.00 99.00 -3.72 -3.09 -2.88 -2.72 -2.58 -2.33 -2.05 -1.64 -1.28 -0.84 0.00 0.67 1.28 1.64 2.33 0.00 0.63 0.84 1.00 1.14 1.39 1.67 2.07 2.44 2.88 3.72 4.39 5.00 5.36 6.05 4.61 3.67 3.36 3.14 2.95 2.62 2.26 1.75 1.32 0.82 -0.07 -0.71 -1.23 -1.52 -2.03 5.10 3.96 3.60 3.35 3.13 2.75 2.35 1.80 1.33 0.80 -0.10 -0.72 -1.12 -1.45 -1.88 4.61 3.67 3.36 3.14 2.95 2.62 2.26 1.75 1.32 0.82 -0.07 -0.71 -1.23 -1.52 -2.03 1.92 1.73 1.67 1.63 1.59 1.52 1.45 1.35 1.26 1.16 0.99 0.86 0.75 0.70 0.59 733.92 662.13 638.45 621.65 607.14 581.94 554.45 515.50 482.66 444.47 376.50 327.63 287.91 265.77 226.82 4.90 3.84 3.50 3.27 3.06 2.70 2.31 1.78 1.33 0.81 -0.09 -0.72 -1.16 -1.48 -1.94 2.28 2.00 1.91 1.85 1.80 1.70 1.60 1.46 1.34 1.21 0.98 0.81 0.70 0.62 0.50 868.72 763.49 729.73 706.10 685.45 649.71 611.59 558.57 513.50 462.07 373.11 310.76 266.29 235.11 189.24 5.00 3.90 3.55 3.31 3.09 2.72 2.33 1.79 1.33 0.80 -0.09 -0.72 -1.14 -1.46 -1.91 Cv Cs 0.26 0.56 Kp=φp*Cv+1 Xp=Kp*X 2.30 2.01 1.92 1.86 1.80 1.71 1.61 1.47 1.35 1.21 0.98 0.81 0.70 0.62 0.50 878.45 769.24 734.50 710.27 689.03 652.29 613.37 559.56 513.70 461.67 372.52 310.57 268.47 236.50 192.22 表(3) 设计年径流经验频率计算表
经验频率计算表 年(月/日)份 水(流)量Y 序号n Xi大→小 至模比系数P=50%Ki-1 (Ki-1)^2 (Ki-1)^3 P=m/(n+1)*100% Ki=Xi/(X/n) 处的水平距离 28 1.58 1.35 1.30 1.28 1.26 1.24 1.18 1.17 1.13 1.11 1.10 1.06 1.06 1.05 1.02 0.95 0.92 0.88 0.84 0.79 0.77 0.76 0.75 0.74 0.71 0.71 0.65 0.64 1.6136 0.1229 0.1924 0.0445 0.0276 0.0214 0.0183 0.0147 0.0055 0.0050 0.0020 0.0012 0.0009 0.0002 0.0002 0.0001 0.0000 0.01 3.45 6.90 10.34 13.79 17.24 20.69 24.14 27.59 31.03 34.48 37.93 41.38 44.83 48.28 51.72 55.17 58.62 62.07 65.52 68.97 72.41 75.86 79.31 82.76 86.21 89.66 93.10 96.55 X 合计 10691.8 28 10691.8 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 270.1 405.9 245.9 602.3 270 390.3 303.4 517.2 287.8 1 2 3 4 5 6 7 8 9 602.3 517.2 497.2 487.8 482.4 475.4 449.3 447.2 429.6 422.2 419 405.9 403.7 401 390.3 361.5 350.9 335.2 320.5 303.4 295.2 289.3 287.8 282.1 270.1 270 249.4 245.9 -3.719 0.00 -1.819 1.90 -1.484 2.24 -1.262 2.46 -1.090 2.63 -0.945 2.77 -0.817 2.90 -0.702 3.02 -0.595 3.12 -0.495 3.22 -0.399 3.32 -0.307 3.41 -0.218 3.50 -0.130 3.59 -0.043 3.68 0.043 3.76 0.130 3.85 0.218 3.94 0.307 4.03 0.399 4.12 0.495 4.21 0.595 4.31 0.702 4.42 0.817 4.54 0.945 4.66 1.090 4.81 1.262 4.98 1.484 5.20 1.819 5.54 0.58 0.3333 0.35 0.1256 0.30 0.0913 0.28 0.0770 0.26 0.0693 0.24 0.0600 0.18 0.0312 0.17 0.0293 0.13 0.0156 0.11 0.0112 0.10 0.0095 0.06 0.0040 0.06 0.0033 0.05 0.0025 0.02 0.0005 350.9 10 295.2 11 422.2 12 482.4 13 429.6 14 320.5 15 249.4 16 447.2 17 419 18 449.3 19 401 20 475.4 21 282.1 22 487.8 23 361.5 24 497.2 25 335.2 26 403.7 27 289.3 28 -0.05 0.0028 -0.0002 -0.08 0.0066 -0.0005 -0.12 0.0149 -0.0018 -0.16 0.0258 -0.0041 -0.21 0.0422 -0.0087 -0.23 0.0515 -0.0117 -0.24 0.0587 -0.0142 -0.25 0.0607 -0.0149 -0.26 0.0682 -0.0178 -0.29 0.0856 -0.0251 -0.29 0.0858 -0.0251 -0.35 0.1203 -0.0417 -0.36 0.1268 -0.0451 四,推求个设计代表年的径流过程
(1)选取代表年:1,选取年径流量接近于设计年径流量的年份; 2,选取对工程较不利的年径流过程线;
根据上述原则初步选择54—55,58—59;和53—54,66—67为5%丰水年和95%枯水年的代表年;并根据年内分配和年平均流量接近多年平均情况的原则,选择56—57,74—75年为50%平水年的代表年。
根据上述初步选定的代表年绘制峰值、流量对比图如下: