2) Устанавливаем возможность сезонного регулирования стока
Так как у нас Wг = млн. м. куб.>Uг = млн. м. куб., то сезонное регулировании возможно.
3) Устанавливаем режим работы водохранилища.
В моём случае будет наблюдаться однократный режим работы так как сумма недостатка Δ d=-5 млн. м3 с октября по март и сумма избытка ΔV=27 млн. м3 наблюдается с апреля по сентябрь. Условия Δd< ΔV выполняются.
4) Определяем начало водохозяйственного года. За начало водохозяйственного года берем первый месяц периодов избытка стока «Апрель».
5) Определяем величину полезного объема (V пл) при однотактном режиме.
Величина полезного объема будет равна дефициту стока
V плз = Δ d=5 млн. м3
2.2 Таблично-цифровой расчет водохранилища без учета потерь
Расчет сезонного регулирования без учета потерь.
Расчеты объемов наполнения и сброса водохранилища ведем в таблице 5.
Расчет по I – варианту правил регулирования стока.
Наполнение водохранилища до Vнпу=11 млн. м3 происходит за счет первых избытков стока.
Расчет ведем в следующем порядке:
а) Объем наполнения водохранилища на начало первого месяца принимаем равному мертвому объему Vмо =6 млн. м. куб.
б) Затем определяем фиктивное наполнение водохранилища на конец месяца
Vф = Vн + (Wр-U)
Определяем наполнение водохранилища на конец месяца исходя из того, что Vф находится между объемами Vмо иVнпу, то есть должно выполняться условие
Vнпу>Vф>Vмо
По этому если полученное значение Vф соотвецтвует этому условию, то конечный объем Vк=Vф.
в) Если Vф>V нпу, то Vк = Vнпу, а остальной объем уходит на сброс
Vсб = Vф – Vнпу.
Полученное значение Vк и Vсб записываем в 6 и 7 колонку.
г) Выполняем расчеты последующих месяцев в таком же порядке.
Расчет по II – варианту правилу регулирования стока.
Первые избытки стока сбрасываются, а за счет последующих наполняется водохранилище.
Расчет ведется против хода времени, то есть с последнего месяца:
а) Объемы наполнения водохранилища на конец последнего месяца принимаем Vк = Vмо.
б) Определяем фиктивное наполнение на начало последнего месяца
Vф = Vк + (Wр – U)
Проверяем условие Vнпу >Vф >Vмо
1) если это условие выполняется Vн = Vф и сбросы Vсб = 0
2) если Vф <Vмо, то объем на Vн = Vмо, а сброс равен Vсб = Vмо – Vф полученные значение заносим в колонку 8 и 9.
2.3 Таблично-цифровой расчет водохранилища с учета потерь
Расчет будем проводить по II – варианту.
Потери в этом случае определяю по ранее построенном графику потерь (рис. 2) в следующем порядке:
а) Определяем средний объем воды водохранилища за каждый месяц (10 колонка)
Vсб = (Vн + Vк)/ 2
б) По батеграфическим кривым (рис. 1) Vн = f (Н) находим площадь зеркала водохранилища по величине среднего объема ω =f (h) (11 колонка)
в) по соответствующему графику потерь по величине ωср величину потерь Vп и записываем в 12 колонку.
г) Определяем дефицит стока (13–14 колонка) с учетом потерь
Wр – U – Vг
д) Уточняем режим работы водохранилища величину Vплз и полного объема
Vплз = 5 млн. м3
Vнпу = 11 млн. м3
е) Производим перерасчет объемов наполнения и сброса с учетом потерь. В порядке зависящих от варианта правил регулирования стока.
2.4 Графические расчеты водохранилища сезонно-годового регулирования стока
Графические способы применяются для предварительных расчетов регулирования стока без учета потерь.
Различают три способа таких расчетов:
1) С помощью полных интегральных кривых стока и отдачи.
2) С помощью сокращенных интегральных кривых стока и отдачи.
3) с помощью разностных кривых стока и отдачи.
Для выполнения этих расчетов выполняем предварительный расчет в таблице 6.
Таблица 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: