Таблица 6. Ордината интегральных кривых
Месяц |
W р |
U |
∑W р |
∑W р -U |
∑W р -К∑U |
1 |
2 |
3 |
4 |
5 |
8 |
IV |
8 |
4 |
0 |
0 |
0 |
8 |
4 |
4 |
|||
V |
12 |
4 |
|||
20 |
8 |
12 |
|||
VI |
10 |
4 |
|||
30 |
12 |
18 |
|||
VII |
8 |
4 |
|||
38 |
16 |
22 |
|||
VIII |
7 |
4 |
|||
45 |
20 |
25 |
|||
IX |
6 |
4 |
|||
51 |
24 |
27 |
|||
X |
4 |
4 |
|||
55 |
28 |
27 |
|||
XI |
4 |
4 |
|||
59 |
32 |
27 |
|||
XII |
3 |
4 |
|||
62 |
36 |
26 |
|||
I |
3 |
4 |
|||
65 |
40 |
25 |
|||
II |
2 |
4 |
|||
67 |
44 |
23 |
|||
III |
3 |
4 |
|||
70 |
48 |
22 |
Графические способы расчета сезонного регулирования стока
Режим работы водохранилища по I Варианту
1) С 1го апреля по 6е мая период заполнения полезного объёма.
2) С 6го мая по 31е октября идёт сброс лишней воды при полном заполнении полезного объёма.
3) С 1го ноября по 31е марта проходит период сработки полезного объёма.
Режим работы водохранилища по II Варианту
1) С 1го апреля по 26е июня наблюдается период сброса.
2) С 27го июня по 31е октября происходит процесс наполнения полезного объёма.
3) С 1го ноября по 31е марта проходит период сработки водохранилища.
3. Расчет методом Крицкого – Менкеля
В этом методе сезонная составляющая полезного объема так же, как в балансовом методе:
Из исходных данных:
Коэффициент за регулирование стока ά=0,7;
Длительность межени в долях года tМ = 7/12 =0,58;
Доля меженного периода mМ =0,3
βсез =ά*(tм – mм )=0,9*(0,58–0,3)=0,25
W – среднемноголетний объем годового стока = 180 млн. м³.
Vсез = βсез *Wг =180*0,25=8,05 млн. м³.
Многолетняя составляющая βмн в этом методе определяется с помощью графиков «Сванидзе».
В зависимости от коэффициентов вариации Сv и асимметрии Сs речного стока.
Коэффициент корреляции стока смежных лет r, расчетной обеспеченности Р=90% и коэффициент регулирования стока ά =0,7.
В начале выбирают расчетный график в зависимости от Сv /Сs , r, обеспеченности р%.
Затем по этому графику в зависимости от Cv и ά определяют β мн.
Vмн = β мн * Wг = 0,56*180=100,8 млн. м³.
βмн = 0,56
Vплз = Vсез + Vмн = 8,05+10= 18,05 млн. м³.
Расчет методом Монте – Карлом.
Он основан на моделирование искусственных гидрологических рядов большой продолжительности (1000 – и более лет).
При этом расчеты обычно выполняют способом обратной задачи.
Зная величину βмн устанавливаем вероятность отдачи Рu%.
Порядок расчетов:
1) По величинам Cs и Cv по заданной реке рассчитываем теоретическую кривую обеспеченности таблица 7.
Таблица 7. Расчет теоретической кривой обеспеченности
Р% |
0,1 |
1 |
3 |
5 |
10 |
20 |
30 |
50 |
70 |
80 |
90 |
95 |
97 |
99 |
К |
2,49 |
2,02 |
1,79 |
1,66 |
1,48 |
1,28 |
1,51 |
0,96 |
0,78 |
0,69 |
0,58 |
0,49 |
0,44 |
0,36 |
И строим теоретическую кривую обеспеченности
1) Моделируем ряд значений обеспеченности годового стока Р%.
Методом генерации случайных чисел.
Дальнейший расчет ведем в табличной форме таблица 8.
2) колонку записываем случайные числа. Они равны обеспеченности.
3) колонка – по кривой обеспеченности определяем модульные коэффициенты Кi и записываем.
4) колонка – для каждого года определяется величина βpi = βнi – ά + Кi
βнi – наполнение водохранилища в начале года. Для первого расчетного года будет равен 0.
Таблица 8. Вычисляем вероятностные характеристики методом Монте-Карло βмн = 0,2 ά = 0,7
№ п.п. |
Р i % |
К i |
βр i |
βк i |
di |
Βсб |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
1 |
66 |
0,82 |
0,12 |
0,12 |
0 |
0 |
2 |
44 |
1,02 |
0,52 |
0,2 |
0 |
0,32 |
3 |
72 |
0,76 |
0,58 |
0,2 |
0 |
0,38 |
4 |
57 |
0,89 |
0,39 |
0,2 |
0 |
0,19 |
5 |
31 |
1,14 |
0,64 |
0,2 |
0 |
0,44 |
6 |
37 |
1,06 |
0,56 |
0,2 |
0 |
0,36 |
7 |
52 |
0,93 |
0,43 |
0,2 |
0 |
0,23 |
8 |
51 |
0,92 |
0,42 |
0,2 |
0 |
0,22 |
9 |
62 |
0,86 |
0,36 |
0,2 |
0 |
0,16 |
10 |
87 |
0,62 |
0,12 |
0,12 |
0 |
0 |
11 |
66 |
0,82 |
0,24 |
0,2 |
0 |
0,04 |
12 |
18 |
1,3 |
0,8 |
0,2 |
0 |
0,6 |
13 |
10 |
1,5 |
1,0 |
0,2 |
0 |
0,8 |
14 |
48 |
1,01 |
0,51 |
0,2 |
0 |
0,31 |
15 |
25 |
1,22 |
0,72 |
0,2 |
0 |
0,52 |
16 |
82 |
0,68 |
0,18 |
0,18 |
0 |
0 |
17 |
13 |
1,4 |
0,9 |
0,2 |
0 |
0,7 |
18 |
48 |
1,01 |
0,51 |
0,2 |
0 |
0,31 |
19 |
3 |
1,79 |
1,29 |
0,2 |
0 |
1,09 |
20 |
17 |
1,31 |
0,81 |
0,2 |
0 |
0,61 |
21 |
75 |
0,74 |
0,24 |
0,2 |
0 |
0,04 |
22 |
12 |
1,42 |
0,92 |
0,2 |
0 |
0,72 |
23 |
4 |
1,7 |
1,2 |
0,2 |
0 |
1,0 |
24 |
59 |
0,88 |
0,38 |
0,2 |
0 |
0,18 |
25 |
4 |
1,7 |
1,2 |
0,2 |
0 |
1,0 |
26 |
75 |
0,74 |
0,24 |
0,2 |
0 |
0,04 |
27 |
60 |
0,87 |
0,37 |
0,2 |
0 |
0,17 |
28 |
17 |
1,31 |
0,81 |
0,2 |
0 |
0,61 |
29 |
4 |
1,7 |
1,2 |
0,2 |
0 |
1,0 |
30 |
92 |
0,52 |
0,02 |
0,02 |
0 |
0 |
Определяем 5 колонку Объем наполнения на конец года βкi, в зависимости от полученного значения βрi:
а) βр > βмн в этом случае наблюдаются избытки воды и принимаем βкi = βмн.
Дефицит стока di = 0, а величина сброса равна βсб = βр – βмн.
б) βмнi > βрi > 0 в этом случае βсб = 0, βdi =0, βкi = βр.
в) βрi < 0 d, в этом случае βкi =0, βdi = βрi, βсб = 0.
Выполнив расчет за первый год, преступаем к расчетам для последующего года.
βр= βн + Кi – ά
βн – принимается βк – за предыдущий год.
Рассчитав полностью таблицу. Определяем обеспеченность отдачи.
Для этого определяем вероятность перебоев по формуле:
Аd =nd/N
где nd – число лет с дефицитом стока. Определяем по 6 колонке
nd = 0
N – Продолжительность ряда равное 30.
Аd = 0/30=0
Находим фактическую обеспеченность по формуле:
Рu = (1 – Аd)* 100%=(1 – 0)*100% = 100%.
Это значение сравнивается с расчетным:
Фактическая обеспеченность Рu =100% больше расчетной обеспеченности Р = 90%, то βмн рассчитано верно.
4. Расчет трансформации паводка водохранилищем способом Качерина
Для предварительных расчетов допускается применять упрасченный способ Качерина.
В этом способе сделано два допущения:
1) Рассчитанный гидрограф max стока представлен в виде треугольника или трапеции.
2) Считается, что увеличение сбросных расходов происходит по линейному закону.
Этот способ применим в случае водослива без затворов с отметкой гребня на уровне НПУ, точность этого способа 5 – 10%, что в полнее достаточно.
Порядок расчетов.
1) Чертится расчетная таблица 9.
водохранилище паводок трансформация сток
Таблица 9. Расчет сбросных расходов
№ п.п |
H ф, м |
ФПУ, м |
V фпу |
V ф |
Qc |
q с |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
1 |
0,5 |
64,7 |
13,5 |
2,5 |
295,9 |
58,29 |
2 |
1,0 |
65,2 |
14,4 |
3,4 |
294,4 |
166,55 |
3 |
1,5 |
65,7 |
17,3 |
6,3 |
289,5 |
305,97 |
4 |
2,0 |
66,2 |
18,0 |
7,0 |
288,4 |
471,07 |
5 |
2,5 |
66,7 |
21,5 |
10,5 |
282,6 |
658,35 |
6 |
3,0 |
67,2 |
22,6 |
11,6 |
280,8 |
865,42 |
2) заполняем 2 колонку задаваясь значением слоем форсировки hф.
Вычисляем отметку ФПУ по формуле:
ФПУ = НПУ + hф
НПУ – находим по графику батеграфических и объемных кривых
3) По батеграфической кривой объемов, по значениям ФПУ определяем объем водохранилища Vфпу и записываем в 4 колонку.
4) Определяем объем форсировки по формуле:
Vф = Vфпу – Vнпу млн. м.³
Объем Vнпу = 11 млн.м.³ И записываем в 5 колонку.
Определяем слой форсировки и объем форсировки.
Рассчитываем объем паводка по формуле:
Wп = 0,5* Т *Qmax.
Из исходных данных:
Максимальный расход паводка Qmax= 300 м/с;
Продолжительность паводка Т = 14 суток.
Wп = 0,5* Qmax *Т* 86400=0,5*300*14*86400=181440000≈181 млн. м³
Рассчитываем max сбросной расход по формуле:
Qс = Qmax*(1 – Vф/Wп)
Полученный результат записываем в 6 колонку.
Рассчитываем сбросные расходы, пропускаемые водосливом.
qc = mb *B√2g* hФ м³/с
Из исходных данных:
Коэффициент расхода водослива mb = 0,47;
Ширина водослива B =80
Строится график
По колонке 2 и 6 строится график Qс = f (hф)
По колонке 2 и 7 строится график qc =f (hф)
Точка пересечения графиков дает нам расчетное значение форсировки, и сброс.
По расчетной величине hф определяем отметку ФПУ и Vф.
Слой форсировки hф = 109 м.;
Объем форсировки Vф = 8,8 млн. м³
Литература
1. Практикум по гидрологии, гидрометрии и регулированию стока / Е.Е. Овчаров, Н.
2. Н. Захаровская, И.В. Прошляков и др.; Под ред. Е.Е. Овчарова. – М.: Агропромиздат, 1988. – 224 с.
29-04-2015, 01:00