УДК 502/504:628.(1-21):628.113
В. В. ИЛЬИНИЧ, Е. А. СВЕТЛОВ, А. В. ИЛЬИНКО
Федеральное государственное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства»
ПРОТИВОПАВОДКОВОЕ РЕГУЛИРОВАНИЕ СТОКА ВОДОХРАНИЛИЩЕМ КОМПЛЕКСНОГО НАЗНАЧЕНИЯ
Предлагается технология противопаводкового управления водохранилищем комплексного назначения на основе использования методов линейного программирования с учетом краткосрочного прогноза стока. Технология проверяется с помощью стохастических моделей стока и функционирования водохранилища.
Водохранилище, гидрологический ряд, стохастическая модель, диспетчерский график, линейное программирование.
There is proposed a technology of flood control by the multipurpose reservoir on the basis of usage of the linear programming method taking into consideration a short-term flow forecast. The technology is checked by means of flow stochastic models and reservoir functioning.
Water reservoir, hydrological series, programming.
Целью настоящей работы является разработка методики оперативного управления водохранилищем комплексного назначения. Обозначенная цель обусловила необходимость решения следующих задач: а) разработка стохастической модели притока речного стока к Краснодарскому водохранилищу; б) разработка имитационной модели функционирования водохранилища; в) уточнение правил управления водохранилищем комплексного назначения на основе рассмотрения различных сценариев притока и управления.
Объект исследования — Краснодарское водохранилище на реке Кубань, за счет которого летом обеспечивается орошение 9 5 % рисовых систем Краснодарского края. Возможен дефицит водопользования, есть вероятность появления катастрофических паводков (например, в июле 2002 г. ущерб от паводка составил 8,65 млрд р).
Одним из наиболее удобных способов стохастического моделирования процесса стока является метод фрагментов, суть которого заключается в моделировании сначала искусственного гидрологического ряда из годовых расходов воды, а затем внутригодового рас-
stochastic model, operations schedule, linear
пределения стока [1]. В качестве исходных параметров для моделирования послужили статистические характеристики, определенные по имеющимся рядам водохозяйственных лет. Фрагменты — гидрографы из 72 пентадных (пятисуточных) значений (i = 1...72) [2]. Стохастическая модель такого типа применительно к Краснодарскому водохранилищу представлена в работе [3].
Моделировали искусственные 1000-летние ряды значений притока к водохранилищу, которые обеспечивают представительство различных сочетаний водности года и различных фазовых интервалов. В результате получили длительный 10-тысячелетний гидрологический ряд пятисуточного притока к Краснодарскому водохранилищу. Анализ соответствия исходного и смоделированных гидрологических рядов проводился на основе сравнения аналитической кривой обеспеченности искусственного ряда с эмпирической кривой обеспеченности исходного ряда наблюдений, оценки степени аналогии внутригодового распределения стока, а также оценки однородности рядов по статистическим критериям [3, 4]. Был сделан следующий вывод: используемый
ПРИРОДООБУСТРОЙСТВО
способ моделирования гидрологических рядов достаточно точен, и имеющиеся отклонения не могут оказать существенного влияния на параметры водохранилища при имитационных расчетах его функционирования. Исходный и смоделированный гидрологические ряды можно отнести к одной генеральной совокупности.
На основе моделирования процесса притока воды к Краснодарскому водохранилищу с учетом особенностей режима притока к водохранилищу и соответствующего возможного прогноза притока была составлена математическая модель его функционирования. Было выявлено, что в пределах пентады суточные расходы воды несущественно отличаются от среднего значения за пентаду. Продолжительность подъема частных паводков обычно не превышает одной декады, а их повторяемость в течение сезона существенно различается от года к году. Одни и те же фазы года могут быть как многоводными, так и маловодными.
Краткосрочный прогноз стока по реке Кубань выдается на 5 сут. Такой прогноз можно считать вполне реальным, так как по предварительным исследованиям время добегания паводковых вод к Краснодарскому водохранилищу с его водосбора составляет от 2 до 5 сут, при этом опасные максимумы формируются от выпадения осадков продолжительностью не менее чем за трое суток. Пяти-суточный прогноз притока в модели представлен как достоверный. Приток за последующую пентаду определен с помощью условных функций обеспеченности притока в зависимости от притока в предыдущую пентаду, так как коэффициент корреляции между значениями притока в смежные интервалы составил 0,5...0,9. Поскольку изменчивость условных функций обеспеченности притока незначительная, в качестве прогнозной величины принималось значение 50 % обеспеченности. И полученное таким образом прогнозное значение притока считалось также достоверным.
Максимальный возможный разрешенный сброс из водохранилища за 10 сут составляет 1382 млн м3, т. е. 138,2 млн м3/сут, или 1600 м3/с при величине полезного объема 1606 млн м3. Таким образом, принимая уровень наполнения водохранилища на отметке НПУ за начало расчета, можно опорожнить его практически на 70 % за две пентады, а если заблаговременно (исходя из прогноза) сбросить некоторый объем, то и на все 100 % .
Учитывая тот факт, что максимальный приток воды к водохранилищу, рассчитанный на обеспеченность, равную 0,1 % , составлял 2500 м3/с, т. е. 215 млн м3/сут, или 2150 млн м за декаду, можно говорить о возможности достаточно надежной защиты водохранилища от переполнения за счет заблаговременного (исходя из прогноза) опорожнения водохранилища на величину, равную 0.540 (2150 — 1606 = 540) млн м3.
Основа модели функционирования водохранилища — обычное балансовое уравнение для наполнения водохранилища по пентадным внутригодовым интервалам:
У.+ = V. + Ж — Х1. — Х2. — Б., (1)
1+1 1 1 1 1 I7 ^ '
где V — наполнение водохранилища к концу расчетного периода (пентады); V. — наполнение водохранилища к началу расчетного периода (пентады); Ж. — объем притока воды к водохранилищу; Х1. — объем фактической во-доподачи из верхнего бьефа водохранилища, включая потери на испарение и фильтрацию; Х2. — объем фактической водоподачи из нижнего бьефа водохранилища; Б. — объем холостого сброса в нижний бьеф.
Для условий возможного дефицита воды начальные правила управления водохранилищем выражались традиционным диспетчерскими графиками, рас-считаными методом С. Н. Крицкого, М. Ф. Менкеля:
если V. < V , то С1* = 0 и С2*. = 0;
1 мо7 1 1 7
= К3С1
если V < V. < Удисп„., то Х1
и Х2. = К3С2.;
1 3 1'
3.
если
Удисп3. < V. <
КЛ1. и Х2. = КЛ2
2 1 12
если
^дисп < 2.
V
клт. и Х2. = клт 1 . . 1
если V. > Удис" . 1.
= С/2.;
Vдис^., то Х1. =
2.
< Vдисп1., то Х1. =
1.
то Х1. = и1. и Х2
№ 3' 2009
(бэ|
V < удиоп < удиоп < удиоп < V ■
мо 3i 2i 1i нпу'
0 < К3 < К2 < К < 1,
где U1. — объем требуемого водозабора из верхнего бьефа водохранилища; U2. — объем требуемого водозабора из нижнего бьефа водохранилища, включающий оанитарно-экологичеокий попуок.
Объемы фактичеокой водоотдачи X1. и X2. определялиоь в завиоимооти от начального наполнения V в ооответотвии с
i
заданным диопетчероким графиком.
Для противопаводкового регулирования иопользовалиоь оледующие ограничения:
еоли V+ > V , то X.2 + S. < Жнпу. б;
i+1 нпу7 i i .об7
еоли V.+ > V , то X.2 + S. < Жмпу. б;
i+1 мп^ i i .об7
еоли V+ > V , то X.2 + S. < Жфпу. б,
i+1 фпу i i .об'
где W — макоимальный разрешенный оброо в нижний бьеф.
Опаоными о точки зрения затопления территории являютоя оледующие уоловия:
1400 м3 < V. < V, ;
i фпу'
W.+, > 250 млн м3.
i+1
При их выполнении для каждых пооледующих двух интервалов реша-лаоь задача линейного программирования, нацеленная на минимизацию затопления территорий. Целевая функция, предназначенная для безопао-ного пропуока макоимального отока через водохранилище, выражалаоь оледу-ющим образом:
min Z = Ь. + Ь.,,
. i+17
где Ь. — превышение объема водохранилища над объемом нормального подпорного уровня.
При этом иопользовалиоь оледую-щие ограничения: Ь. > 0; i = 1, ..., n;
V. л = const;
.=1 7
X1 — X2. — S.;
i i i i7
X1. = U1.;
ii
X2. + S. > U2.;
i i i'
X2. + S. < W. б;
i i .об7
V+ > V ;
i+1 умо7
V.+1 < V + b.;
i+1 нпу i7
Vi+2 = Vi+1 + Wi+1 — X1i+1 — X2i+1 — Si+1;
V.+, = V. + W.
i+1 i 1
X2i+1 + Si+1 > u^;
X2i+1 + s+1i < Wlсб; V.+2 > V ;
1+2 умо7
V.+2 < V + b .
1+2 нпу i+1
Технология управления водохранилищем проверялась на исходном и смоделированном стохастическом рядах. При использовании имитационной стохастической модели были получены следующие результаты: вероятность превышения максимального значения форсированного подпорного уровня была понижена с 0,1 до 0,04 %. Для самого многоводного года наблюдений величина форсировки объема была снижена с 709 до 350 млн м3.
Список литературы
1. Сванидзе, Г. Г. Основы расчета регулирования речного стока методом Монте—Карло [Текст] / Г. Г. Сванидзе. — Тбилиси : Мецниереба, 1964. — 272 с.
2. Светлов, Е. А. Обоснование выбора дискретности при разработке математической модели водохранилища [Текст] / Е. А. Светлов // Роль приро-дообустройства сельских территорий в обеспечении устойчивого развития АПК : материалы Международной научно-практической конференции МГУП,
2007. — С. 277-281.
3. Светлов, Е. А. Стохастическая модель водохранилища [Текст] / Е. А. Светлов, В. В. Ильинич // Природообустройсво. —
2008. — № 4. — С. 59-64.
4. Светлов, Е. А. Оценка однородности годового стока реки Кубань [Текст] / Е. А. Светлов // Мелиорация и водное хозяйство. — 2009. — № 1. — С. 19-20.
Материал поступил в редакцию 20.03.09. Ильинич Виталий Витальевич, кандидат технических наук, профессор
Тел.8 (495) 976-17-45
E-mail: [email protected]
Светлов Евгений Александрович, аспирант
Тел.: 8 (495) 907-73-69, 8 (926) 167-15-10
E-mail: [email protected]
Ильинко Ататолий Васильевич, доцент
Тел. 8 (495) 153-97-66