
Задание 5.
Построить кривую обеспеченности минимальных среднемесячных (за многолетний период) расходов источника и определить минимальный среднемесячный расход 85% обеспеченности.
№ п/п |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
Год |
1985 |
1986 |
1987 |
1988 |
1989 |
1990 |
199130 |
Q л/с |
30 |
19 |
38 |
34 |
29 |
31 |
22 |
№ п/п |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
Год |
1992 |
1993 |
1994 |
1995 |
1996 |
1997 |
1998 |
Q л/с |
26 |
25 |
13 |
27 |
24 |
28 |
32 |
Для построения теоретической кривой обеспеченности расходов источника можно воспользоваться функцией распределения Шарлье, которая соответствует общему виду кривой распределения
1)по формуле Q0 = Σ Qi / n вычисляется норма расхода, т.е. среднее арифметическое из n наблюдений расхода.
Q0 = 378/14 = 27 л/с
2)определим модульный коэффициент к = Qi / Q0 , где Qi – значение расхода для этого года наблюдений
к1 = 1,11 к8 = 0,96
к2 = 0,70 к9 = 0,93
к3 = 1,41 к10 = 0,48
к4 = 1,26 к11 = 1
к5 = 1,07 к12 = 0,89
к6 = 1,15 к13 = 1,04
к7 = 0,81 к14 = 1,19
Σкi = 14
3)определим коэффициент вариации по формуле Сv =( Σ(кi – 1)2 /(n-1))1/2
Сv = 0.235
4) по формуле х=(к-1)/Сv вычисляется нормированная величина признака, в нашем случае расхода
х1 = 0,46 х8 = 0,17
х2 = 1,25 х9 = 0,29
х3 = 1,71 х10 = 2,17
х4 = 1,08 х11 = 0
х5 = 0,29 х12 = 0,46
х6 = 0,63 х13 = 0,17
х7 = 0,79 х14 = 0,79
5) По величине х определяется функция φ(х) и затем Σ φ(х)
φ(х1 ) = 0,3589 φ (х8) = 0,3932
φ (х2 ) = 0,1826 φ(х9) = 0,3825
φ (х3) = 0,0925 φ(х10) = 0,0379
φ (х4 ) = 0,2227 φ(х11) = 0,3938
φ(х5) = 0,3825 φ(х12) = 0,3589
φ(х6) = 0,3271 φ(х13) = 0,3932
φ(х7) = 0,2920 φ(х14) = 0,2920
Σ φ(хi ) = 4,1149
6)определяем значение хi * φ(хi ) и затем Σ хi * φ(хi )
х1 φ(х) = 0,165094
х2 φ(х) = -0,22625
х3 φ(х) = 0,158175
х4 φ(х) = 0,240576
х5 φ(х) = 0,110925
х6 φ(х) = 0,206073
х7 φ(х) = -0,23068
х8 φ(х) = -0,066844
х9 φ(х) = -0,110925
х10 φ(х) = -0,082243
х11 φ(х) = 0
х12 φ(х) = -0,165094
х13 φ(х) = 0,066844
х14 φ(х) = 0,23068
Σ хi * φ(хi ) = 0,294274
7)вычисляем взвешенные моменты по формуле b0 =2√π/n * Σ φ(хi )
b0 = 1,042
b1 = 0.053
8) Определяем эксцесс по формуле E = 4(b0 – 1)
E = 0.168
9) по формуле S = 4b1 вычисляется характеристика асимметрии
S = 4*0.053 = 0.2
10) по асимметрии S находим отклонение ординат кривых обеспеченности от середины для разных процентов обеспеченности ( при Е = 0 и Сv = 1)
11) по эксцессу Е находим поправку на эксцесс
12) Прибавив поправку на эксцесс к отклонению ординат кривой обеспеченности от середины ( при Е=0, Сv = 1), получаем отклонение х` с учетом эксцесса.
13) по формуле к1 = Сv * х` + 1 вычисляем модульный коэффициент
14) по формуле Q = Q0 * к1 вычисляем расходы разной обеспеченности и по ним строим кривую обеспеченности, по которой строим определяем расход 85% обеспеченности.
Q0 = 27 л/с Сv = 0,24