автореферат диссертации по строительству, 05.23.16, диссертация на тему:Противопаводковое управление водохранилищем многоцелевого назначения

СЮ340 Ю72

На правах рукописи

Светлов Евгений Александрович

Противопаводковое управление водохранилищем многоцелевого назначения

Специальность 05.23.16 - Гидравлика и инженерная гидрология

Автореферат диссертации на соискание ученой степени кандидата технических наук

Москва 2009

003481072

Работа выполнена в Московском государственном университете природообустройства на кафедре «Гидрология, метеорология и регулирование стока»

Научный руководитель: - кандидат технических наук, профессор,

Ильинич Виталий Витальевич

Официальные оппоненты: - доктор технических наук,

заслуженный мелиоратор РФ Антипова Татьяна Николаевна.

- кандидат технических наук, профессор Раткович Лев Данилович

Ведущая организация: ОАО «Институт Гидропроект им. СЛ. Жука»

Защита состоится 16 ноября 2009 г. в 1500 ч. на заседании диссертационного совета Д 220.045.02 в Московском государственном университете природообустройства по адресу: 127550, Москва, ул. Прянишникова, д. 19, аудитория 201/1.

С диссертацией можно ознакомиться в научной библиотеке Московского государственного университета природообустройства.

Автореферат разослан «15» октября 2009 г.

Ученый секретарь диссертационного совета, доцент,

кандидат технических наук - И.М. Евдокимова

Общая характеристика работы

Актуальность темы. Во многих странах мира и в большинстве субъектов Российской Федерации безопасность жизнедеятельности населения в значительной степени зависит от катастрофических дождевых паводков. Ежегодно от таких паводков гибнут люди, и наносится существенный ущерб экономике. В России ярким примером этого являлся дождевой паводок, произошедший в 2002 году в бассейне реки Кубань. В результате паводка было подтоплено 130 населенных пунктов, разрушено около 70 тыс. домов, пострадало более 130 тыс. человек, в том числе 96 человек - с летальным исходом. Общий ущерб от паводка составил 8,65 млрд. руб. Одной из эффективных мер противодействия разрушительной силе паводков является противопаводковое регулирование стока водохранилищем. В настоящее время при создании и эксплуатации водохранилищ выделяется специальная противопаводковая емкость, предназначенная для задержания избытков речного стока, таким образом, любое водохранилище должно являться либо комплексным, либо специально противопаводковым.

Задача установления рационального режима работы водохранилища комплексного назначения, весьма сложна из-за противоречивого и неоднозначного характера его использования. В частности, при определенных режимах речного стока и водопользования, характерных для юга России, в различные по водности годы в одноименные фазовые периоды могут происходить как засухи, так и катастрофические паводки. В этих случаях, для целей орошения и водоснабжения в маловодные периоды водохранилище должно находиться в наполненном состоянии, а для задержания максимальных расходов катастрофических паводков должно быть достаточно опорожнено. Указанное противоречие может быть отчасти устранено при заблаговременных и точных прогнозах о сроках прохождения паводков, их объемов и характере гидрографа. Однако современный уровень прогнозирования стока не может в полной мере обеспечить водное хозяйство соответствующими прогнозами. Этим обуславливается целесообразность разработки максимально безопасных правил управления водохранилищами при использовании детальной оценки ограниченной информации о гидрологическом режиме водоисточника. Соответственно, в настоящей работе разрабатывается и обосновывается методика противопаводкового регулирования стока на основе анализа формирования притока речных вод к водохранилищу, а также режима требований водопользователей.

Цель и задачи исследований. Основной целью диссертационной работы является разработка методики определения оперативных правил управления водохранилищем комплексного назначения в условиях паводкового режима, позволяющих достичь понижения вероятности аварии при регулировании максимального стока. Для достижения поставленной цели потребовалось решение комплекса взаимосвязанных задач:

• Анализ и оценка режима речного притока к водохранилищу;

• Разработка математической модели притока речных вод к водохранилищу;

• Проверка разработанной модели относительно точности исходной информации;

• Разработка модели функционирования водохранилища комплексного назначения;

• Разработка методики определения оперативных правил регулирования максимального стока водохранилищем комплексного назначения.

Методика исследований. Основой методологии являются современные подходы к математическому моделированию режимов стока рек и функционирования водохранилища комплексного назначения. В качестве математического аппарата используется имитационные стохастические модели и методы линейного программирования.

Объектом исследования является Краснодарское водохранилище и сток питающих его рек, характеризующихся прохождением катастрофических паводков в период интенсивного водопользования.

Научная новизна работы состоит в следующем:

• разработана методика формирования безопасных и рациональных правил оперативного управления водохранилищем на основе имитационного моделирования функционирования водохранилища с использованием методов линейного программирования;

• разработана вероятностная модель процесса стока по типу Марковской цепи с пентадной дискретностью;

• получены новые режимные характеристики Краснодарского водохранилища по выделенным внутригодовым интервалам.

Научная гипотеза, проверяемая в диссертационной работе, заключается в следующем: возможны такие правила регулирования паводкового стока, отличные от традиционных и ранее применявшихся, при которых: вероятность аварии на гидроузле будет уменьшена, а фактическая водоотдача и обеспеченность ее плановой величины будут стабильны.

Практическая ценность работы. На основе полученных результатов разработана методика оперативного управления трансформацией паводка водохранилищем с учетом краткосрочных гидрологических прогнозов. Данная методика представляет собой инструмент для моделирования режимов функционирования водохранилищ в различных сценариях изменяющихся природно-хозяйственных условий и обеспечивает рентабельное и безопасное регулирование стока. Предложенная методика может быть использована для построения диспетчерских правил управления водными ресурсами водохранилища. Доступный инструментарий методики позволяет применение ее непосредственно на объекте без использования специального программного обеспечения. Разработанная методика позволяет уменьшить вероятность

аварийной ситуации при прохождении катастрофических паводков. Результаты исследований приняты в Кубанском бассейном водном управлении для рассмотрения о внедрении в работу оперативных служб.

Достоверность полученных результатов подтверждается сравнением точности результатов моделирования и исходных материалов наблюдений, а также использованием современных подходов имитационного моделирования.

Апробация работы. Основные результаты диссертационной работы докладывались:

1) на заседании Кубанского бассейнового водного управлении и были приняты для рассмотрения о возможности внедрения в работу оперативных служб КБВУ (август 2009 г.),

2) на международной научно-практической конференции «Агротехнологии 21 века» (МСХА им. КА. Тимирязева, 2007 г.),

3) на научно-технических конференциях МГУП (2006-2009 гг.),

4) на заседаниях кафедры гидрологии, метеорологии и регулирования стока МГУП (2006-2009 гг.),

5) на занятиях с магистрами на кафедре гидрологии, метеорологии и регулирования стока МГУП (2006-2008 гг.)

Публикации. По теме диссертации опубликованы 6 печатных работ, две из которых опубликованы в журналах, входящим в список ВАК.

Структура и объём диссертации. Диссертация состоит из введения, четырех глав, заключения и списка используемой литературы. Работа изложена на 145 страницах машинописного текста, иллюстрированного 25 рисунками, и содержит 16 таблиц. Список используемой литературы состоит из 125 наименований.

Содержание работы.

Во введении обоснована актуальность темы диссертационной работы, сформулированы цель и задачи исследования, методы исследований, научная новизна, практическая ценность, достоверность полученных результатов и структура работы.

Первая глава посвящена рассмотрению современного состояния проблемы трансформации паводков водохранилищем комплексного назначения. В ней анализируются существующие методы определения противопаводковой емкости водохранилища комплексного назначения и методы управления водохранилищем в период катастрофических паводков.

С целью предотвращения ущербов от наводнений для водохранилищ комплексного назначения предусматривается специальная резервная емкость, предназначенная для временной аккумуляции излишков стока в период паводков. Такую емкость называют противопаводковой. Вопросам ее расчета и использования при трансформации паводков посвятили свои работы

Алексеев Г.А., Асарин А.Е., Бабурин БЛ., Железняк И.А., Жиркевич А.Н., Исмайылов Г.Х., Крицкий С.Н., Менкель М.Ф., Подольский Е.М., Потапов М.В., Раткович ДЯ., Ржаницын Н.А., Сванидзе Г.Г., Цингер В.Н. и многие другие исследователи. Проведенный анализ работ позволил выяснить, что большинство моделей регулирования максимального паводкового стока не используют в полной мере современные возможности оперативного управления режимом работы водохранилища. В частности, не используются хорошо оправдывающиеся краткосрочные прогнозы стока, практически не используются также и условные стохастические зависимости распределения вероятностей величин стока в будущем временном интервале от конкретных величин стока в текущем интервале.

Вторая глава посвящена рассмотрению гидрологических условий речного бассейна Краснодарского водохранилища и включает в себя анализ и оценку многолетних и внутригодовых колебаний притока речного стока к водохранилищу. Рассмотрены основные выдающиеся наводнения в бассейне реки Кубань, причины их возникновения, а также условия формирования стока рек, питающих Краснодарское водохранилище.

Краснодарское водохранилище контролирует около 96% стока реки Кубань и является самым крупным искусственным водоемом Краснодарского края. Основным источником питания водохранилища является река Кубань. Менее значительными донорами являются реки Белая, Лаба, Пшиш, Псекупс и др.

Анализ многолетних колебаний годового стока р. Кубань в створе г. Краснодара выполнен по материалам наблюдений за период 1911 - 2006 гг. (п = 95 лет). Годовой приток к Краснодарскому водохранилищу не отличается большой водностью и изменчивостью и достаточно четко отражается представленной в диссертационной работе теоретической кривой обеспеченности. Коэффициент автокорреляции г между величинами стока в смежные годы не превышает 0,3. Такая величина соразмерна со среднеквадратической ошибкой этого параметра.

Анализ разностной интегральной кривой годового стока показал наличие двух четко выраженных периодов изменения водности. С 1911 по 1972 гг. водность р. Кубань имела незначительные колебания различной направленности и с 1973 по 2006 гг. четко прослеживалась тенденция к увеличению значений годового стока. В связи с этим для более объективной оценки цикличности и тренда многолетних колебаний стока р. Кубань исходный 95-летний гидрологический ряд был разбит на два более коротких. Первый ряд п1 состоял из последовательных значений годового стока с начала наблюдений - 1911 г. по 1972 г. (п1 = 61 год). Второй ряд п2 состоял из оставшихся наблюдений с 1973 г. по 2006 г. (п2 = 33 года). Проведенный анализ включает оценку квазипериодичности и автокорреляции, тренда и статистической однородности рядов.

При анализе тенденции к группировке лет повышенной и пониженной водности (цикличности многолетних колебаний речного стока) была

рассмотрена автокорреляция исследуемых радов между стоком смежных лет (г1) и стоком лет, удаленных на один (г2) и два года (гЗ), а также использован непараметрический критерий серий (1(и)), основанный на подсчете числа серий лет (и) повышенной или пониженной водности. Для выяснения наличия в исследуемых рядах монотонного возрастающего или убывающего тренда использован непараметричсский критерий тренда Спирмсна (ге).

Табл. 1

Данные анализа цикличности и тренда но р. Кубань

п г! г2 гЗ и 1(и) ГБ

Ряд 1 (1911-1972) 61 0,26 0,13 -0,09 25 -0,002 -0,07

Ряд 2 (1973-2006) 33 0,30 0,09 -0,01 12 -0,007 0,51

Весь ряд (1911-2006) 95 0,30 0,15 -0,10 37 -0,005 0,19

Данные анализа цикличности и тренда р. Кубань (табл.1) показывают, что для обоих гидрологических рядов имеет место небольшая, но статистически достоверная корреляция между стоком смежных лет. Корреляция же для двухзвенной и трехзвенной Марковской цепи (г2 и гЗ) оказалась статистически недостоверной. Исходя из этого, можно принять гипотезу об отсутствии корреляции в несмежные годы для реки Кубань. Анализ критерия серий показал отсутствие статистически достоверной тенденции к группировке лет повышенной и пониженной водности для реки Кубань. Согласно оценке непараметрического критерия тренда Спирмена можно сделать вывод, что ряд многолетних колебаний стока реки Кубань приобрел явно возрастающий тренд в последнее тридцатилетие.

В соответствии с вышеизложенным, можно заключить, что в последнее тридцатилетие водность реки Кубань имеет тенденцию к незначительному увеличению за счет появления отдельных значительно отличающихся от среднемноголетнего стока его годовых величин.

Внутригодовое распределение притока речных вод к Краснодарскому водохранилищу отличается большой изменчивостью однофазных величин. Годовые гидрографы имеют пилообразный вид. Половодье начинается резким подъемом уровней, как правило, в мае и проходит в течение вссснне-летнего периода с апреля по август, вследствие наложения паводков ледникового и дождевого происхождения. Из анализа гидрографов внутригодового распределения притока воды к водохранилищу (рис. 1) для выборочных лет (1986 и 1989 гг.) видно, что практически невозможно построить достаточно достоверный типовой гидрограф расходов воды или получить типовое внутригодовое распределение стока.

1400 1200 1000 800 600 400 200 0

1 я —

А Ас \ /1 А ~

/ ч^ АГ \ /\ 1 * 11Л/ \ лЛ

Дут4« \ \ 1 1\ ) А , 1

1989 год

; -! 2 ? ; о - 2 л 2 « з ^ А А з 2 « о ¡з ¡4 ¡в <1 в й <• « Ч! *< и » <4 ** " м ** « ®

Рис. 1 Гидрографы внутригодового распределения притока речных вод к Краснодарскому водохранилищу для 1986 и 1989 гг..

Для статической характеристики внутригодового распределения стока по данным наблюдений были определены средние значения расходов притока к водохранилищу и их коэффициенты вариации по выделенным 72 внутригодовым интервалам. Из их анализа можно сделать вывод, что наибольшую опасность появления аварийной ситуации при пропуске катастрофических паводков Краснодарским водохранилищем представляет период с июня по август месяц.

Третья глава посвящена построению стохастической модели притока речных вод к водохранилищу комплексного назначения и стохастической модели функционирования Краснодарского водохранилища.

В качестве исходных данных для моделирования был принят гидрологический ряд суточных объемов суммарного притока воды к Краснодарскому водохранилищу с 1973 г. по 2005 г. (всего 32 водохозяйственных года), полученный в Кубанском Бассейновом Водном Управлении и, непосредственно, в ФГУ «Краснодарское водохранилище».

Для всех современных моделей стока и управления водохранилищами приходится решать вопрос выбора расчетных интервалов - дискретности. Сущность дискретизации состоит в выборе наиболее оптимального временного интервала, в пределах которого сток можно характеризовать осредненным значением за выбранный период. В проведенных исследованиях для объективного выбора дискретности моделирования были рассмотрены следующие критерии: а) время добегания паводка от истока реки Кубань до Краснодарского водохранилища; б) степень превышения максимального значения расхода в пределах расчетных временных интервалов над средним значением за эти интервалы; в) статистические характеристики по расчетным временным интервалам; г) стандартные краткосрочные прогнозы гидрометслужбы по осадкам и температурам.

Было установлено, что время добегания паводка от истока р. Кубань до Краснодарского водохранилища составляет примерно 5 суток. Значение коэффициента отношения максимальных значений суточных расходов, наблюдавшихся в пределах конкретного временного интервала, к

соответствующему среднему пентадпому (пятисуточному) расходу находится в пределах допустимой точности измерения самих этих величин. Анализ статистических характеристик стока (0ср, Су, Сб, г) для различных внутригодовых периодов также показал предпочтительность выбора пятисуточного временного интервала для дальнейших исследоватй. Также при дискретизации модели в 5 суток (пентада) имеется возможность учитывать краткосрочные прогнозы стока при оперативном управлении водными ресурсами Краснодарского водохранилища. Достоверность таких прогнозов достаточно высока, и вероятность их оправдываемости приближается к единице. Соответственно, достоверность среднесрочных прогнозов стока (заблаговрсменность до 10 суток) и долгосрочных прогнозов (заблаговременпость свыше 10 суток) гораздо меньше.

Таким образом, проведенный комплексный анализ выбора дискретности модели управления Краснодарским водохранилищем позволяет достаточно объективно сделать вывод, что наиболее подходящим расчетным интервалом в данном случае является пятидневный интервал (пентада).

Стохастическая модель притока речного стока к водохранилищу комплексного использования разрабатывалась на основе метода Монте-Карло в реализации способа фрагментов Г.Г. Сванидзе. Суть этого способа в двойном моделировании: искусственного гидрологического ряда из среднегодовых расходов воды, а затем модели внутригодового распределения стока, так называемого фрагмента.

В результате получались длительные искусственные ряды колебаний пятисуточных расходов воды. Далее решался вопрос о необходимой длительности искусственного ряда с точки зрения стабилизации статистических параметров искусственного гидрологического ряда при различной его длительности. С этой целью были построены графики изменения статистических характеристик в зависимости от числа лет смоделированного гидрологического ряда объемов притока к Краснодарскому водохранилищу. По таким графикам можно судить о том, когда стабилизируются все статистические характеристики. Анализировались средние пентадные значения стока, их коэффициенты вариации и ассиметрии Су и Сб. Анализ показал, что уже через 600-700 лет моделирования все характеристики принимают стабильное значение с небольшим колебанием около их среднего значения. При дальнейшем моделировании уменьшение степени таких колебаний не происходит. Основываясь на таком факте, можно сделать вывод, что для получения стабильных статистических характеристик искусственных гидрологических рядов расходов достаточно моделировать ряды длительностью 1000 лет.

Следующим этапом проводился анализ степени совпадения выборочных режимных статистических характеристик исходного и смоделированного гидрологических рядов. Несмотря на то, что способ фрагментов неоднократно проверялся различными исследователями на многочисленных экспериментах, которые подтвердили правомерность его

применения в водохозяйственных задачах, анализ степени совпадения исходного и смоделированного гидрологических рядов дает каждый раз возможность сделать конкретные выводы по отношению к конкретной задаче и конкретному объекту. На первом этапе использования способа фрагментов моделируются годовые расходы воды. Поэтому в первую очередь анализировалось совпадение исходного ряда годовых расходов и искусственного. Для решения этой задачи проверялось соответствие теоретической и эмпирической кривых обеспеченностей. Максимальное расхождение вероятностей превышения достигает 0,02, что не превышает критического значения 0,043, определенно по таблице Колмогорова А.Н.. Следовательно, принятая функция обеспеченности годового стока с параметрами С\'=0,24, Св=2Су не отвергается.

Норма стока по 32-летнему ряду равна 397 мЗ/с. Среднее многолетнее значение стока по 1000 летнему ряду равно 395 мЗ/с. Различие между этими двумя величинами составляет 0,5%, что вполне допустимо. Коэффициенты вариации и ассиметрии для величин стока искусственного и исходного рядов также идентичны.

Чтобы оценить соответствие внутригодового распределения стока исходного и искусственного гидрологических рядов, были сопоставлены и статистические параметры стока по 72 выделенным внутригодовым пентадным интервалам (рис 2, 3). Максимальное расхождение между коэффициентами вариации при внутригодовом распределении составляет 13%, что вполне допустимо. Максимальное различие между средними значениями расхода смоделированного и исходного рядов практически не превышает 10% (только за два внутригодовых интервала с 21 по 25 марта и с 26 по 31 марта составляет 15%), что тоже вполне допустимо, так как сами измерения максимальных расходов проводятся с точностью 10-15%, а относительные среднеквадратические ошибки этих параметров при имеющихся рядах наблюдений и принятых законах распределения составляют 10 -15%.

1-1-111ё з з § § Ш 1I Ж I ж Ж §•

2 2 ^ ~ « * 2 * 2.2 =2 2 # гч 2; « г- м -г- «м -г-

Рис. 2 Соотношение коэффициентов вариации по внутригодового распределения для естественного и смоделированного рядов. 10

юао аоо 800 го о

Оср исх

АОО ['-

Щ зоо 4

2 200 -1-— я

юо +

о -А

Рис. 3 Соотношение средних расходов притока по внутригодового распределению для естественного и смоделированного рядов.

Кроме того, была проведена статистическая оценка однородности наблюденного и искусственного гидрологических рядов. В качестве критериев однородности средних расходов использовался критерий Стьюдента I, а для оценки однородности дисперсий критерий Фишера Б. Для рядов годовых расходов критерий Фишера получился равным Р=012/022=1.<002, что значительно меньше его критического значения Ба =1,4. Значение критерия Стьюдента 1=0,14 также намного меньше критического (2,58). При статистической оценке однородности рядов для внутригодового распределения стока практически во всех случаях получились результаты, при которых нуль-гипотеза однородности расходов и дисперсий не может быть отвергнута.

Также был проведен анализ различия коэффициентов корреляции между величинами стока смежных и несмежных внутригодовых интервалов по исходному и смоделированному рядам пятисуточных расходов воды. Коэффициенты корреляции определялись между величинами стока за расчетный интервал ш и предыдущий интервал ш-1, а также между расчетным и предшествующий предыдущему ш-2. Результаты сопоставления коэффициентов корреляции г(ш-1) и г(ш-2) исходного и искусственного гидрологических рядов по внутригодовым интервалам также не дают повода ' отвергнуть гипотезу о принадлежности исходного наблюденного и искусственных рядов к одной генеральной совокупности случайных величин. При дальнейших водохозяйственных расчетах достаточно учитывать коэффициенты корреляции только между стоком только в смежные внутригодовые интервалы.

Таким образом, был сделан вывод о том, что принятый способ моделирования гидрологических рядов достаточно точен. Величины стока моделируются с достаточной точностью, и имеющиеся отклонения не могут существенно повлиять на характеристики регулирования стока Краснодарским водохранилищем. Более того, можно говорить о том, что

исходный ряд и смоделированные ряды относятся к одной генеральной совокупности случайных величин. Выбранный способ стохастического моделирования стока является достаточно корректным для условий Краснодарского водохранилища.

Кроме модели притока речных вод к водохранилищу, основанной на статистических испытаниях способом фрагментов, в представляемой диссертационной работе разработана еще одна вспомогательная модель, предназначенная для выбора рациональных правил регулирования паводкового стока в целях уменьшения вероятности аварии и рентабельности водоотдачи из водохранилища. В этой модели процесс стока представляется по типу простой Марковской цепи - совокупностью условных функций обеспеченности объемов стока, поступающих в водохранилище в пределах выбранных 72 внутригодовых интервалов. Для выявления рациональных правил регулирования паводкового стока целесообразно характеризовать условными функциями обеспеченности величины стока в интервалах, характеризующих период прохождения паводков и период интенсивной водоотдачи. Параметры условных функций обеспеченности определялись следующим образом. Условные средние значения объемов притока по формуле:

WCpм(Wм-l) = \УсрМ + гМ-1 См/ом-1 • - \VcpM.0 (1)

Где: \УСрм и \УсрМ_1 - соответствующие фактические средние значения притока в М-ый и (М-1)-ый интервалы; ом и ом.1 - фактические среднеквадратические отклонения в М-ый и (М-1)-ый интервалы; \VcpmCWm.i) -условные средние объемы притока

Второй параметр условных функций обеспеченности притока -коэффициент вариации вычисляется по формуле:

Сум/М., = (СУм • А^/ ^ГсрмО^м-О) • О-^млм)"2 (2)

где, помимо прежних обозначений: №срМ, С\'м — ч соответственно среднее значение притока и коэффициент вариации за рассматриваемый М-ый интервал.

Третий параметр условных функций обеспеченности - коэффициент асимметрии, согласно принятой гипотезе, жестко связан с принятым типом распределения случайных величин:

Сэм/мп = К • Сум/М-1 (3)

Где: К - коэффициент, принимаемый равным значению пропорциональности исходных (безусловных) функций.

Большинство водохозяйственных задач при эксплуатации водохранилища, невозможно решить, основываясь исключительно на аналитических методах, из-за ряда весомых причин. Альтернативой этим методам служит математическое моделирование, применимое для решения задач проектирования и эксплуатации, как водохозяйственных систем в целом, так и водохранилищ в частности. Как известно, основные параметры

водохранилища комплексного назначения и правила регулирования стока устанавливаются из сопоставления режимов стока и требуемого водопотребления, а также из условий пропуска паводков. При этом желательно получить такие правила регулирования, при которых расчетные объемы водохранилища (в том числе и используемая противопаводковая емкость) и ущерб от дефицита воды были как можно меньше, а величина обеспеченности плановой водоотдачи как можно больше. Такие задачи достаточно корректно можно решить с помощью математических моделей работы водохранилища. В представляемой работе способом решения поставленных задач принято стохастическое моделирование для анализа долговременной работы водохранилища, и оптимизационное для выбора наилучшего решения внутри каждого отдельно взятого интервала. Разработанная модель функционирования водохранилища комплексного назначения является имитационной в том смысле, что она предназначена для анализа долговременного функционирования водохранилища комплексного назначения. Внутри же каждого расчетного периода решается задача оптимизации, которая является вспомогательной для выбора наилучшего решения.

На основе моделирования процесса притока воды к Краснодарскому водохранилищу и его проектных характеристик была разработана математическая модель функционирования гидроузла. Модель составлялась исходя из особенностей режима притока к водохранилищу и соответствующего возможного прогноза притока. Заблаговременность краткосрочного прогноза стока по реке Кубань составляет 5 суток. Такой прогноз можно считать достаточно реальным. Поэтому пятисуточный прогноз притока в модели представлялся, как достоверный. Приток за последующую пентаду определялся с помощью условных функций обеспеченности притока в зависимости от притока в предыдущую пентаду, так как коэффициент корреляции между значениями притока в смежные интервалы составляет 0,5-0,9. Поскольку изменчивость условных функций обеспеченности притока незначительна, в качестве прогнозной величины принималась значение 50% обеспеченности.

В основе модели функционирования водохранилища было положено обычное балансовое уравнение для наполнений водохранилища по пентадным внутригодовым интервалам.

Ум = V; + -XI, -Х2; - Б; (4)

Где: - наполнение водохранилища к концу расчетного периода (пентады); V; - наполнение водохранилища к началу расчетного периода (пентады); - объем притока воды к водохранилищу; ХГ - объем фактической водоподачи из верхнего бьефа водохранилища включая потери на испарение и фильтрацию; X2i - объем фактической водоподачи из нижнего бьефа водохранилища; - объем холостого сброса в нижний бьеф.

Для условий возможного дефицита воды начальные правила управления водохранилищем выражались традиционным диспетчерскими

13

графиками, рассчитанными методом , Крицкого С.Н., Менкеля М.Ф.. Применительно к Краснодарскому водохранилищу правила регулирования стока сводятся к следующим граничным условиям:

1) Если Vi < Vmo , TOUl'i = 0 и U2*i = О

2) Если Vmo < Vi < УДЙСП3-,, toX1; = k3*U1; и Х2; = кЗ* U2;

3) Если Удасп3( < Vi < Удиспя , то Xli = к2* Uli и Х2,=к2* U2i

4) Если V№C,,2i < Vi < Vwcnn , то Xli = Kl* Uli и Х2; = к1* U2;

5) Если Vi > Vmc\ , то Xli = Ulj и X2i = U2i

Vmo < VOTCn3i < Vwcn2i < Vjwal,i < VHny ; 0 < КЗ < K2 < Kl < 1

где Ulj - объем требуемого водозабора из верхнего бьефа водохранилища; U2i - объем требуемого водозабора из нижнего бьефа водохранилища включая санитарно-экологический попуск;

Объемы фактической водоотдачи Xli и X2j определяются в зависимости от начального наполнения (Vi) в соответствии с заданным диспетчерским графиком. Потери воды из водохранилища на испарение и фильтрацию учитываются в величинах водоотдачи.

Для противопаводкового регулирования использовались следующие ограничения:

1) Если Vi+1 > Vnny ,то Xi2 + S;<W"nyic6;

2) Если Vi+i> VMny , то Xi2 + Si < \4vnyicS;

3) Если Vi+1 > Vфпy ,то X^ + Si^W^ice;

Где Wic6 ~ максимальный разрешенный сброс в нижний бьеф;

Во время прохождения катастрофических паводков при достижении НПУ VHiry=Vmaxb в нижний бьеф пропускаются расходы не превышающие значения Qmaxi до заполнения части противопаводковой емкости VMny=Vmax2, далее разрешено сбрасываются расходы не превышающие Qmax2. После заполнения всей основной противопаводковой емкости Vфпy=Vmaxз разрешено сбрасывать расходы не превышающие Qmax3. Только при невозможности пропустить паводок с такими условиями регулирования используются все пропускные устройства. Однако такие действия являются аварийными, поскольку в нижнем бьефе создается катастрофическая ситуация. В течение зимнего периода (декабрь, январь и февраль) условия пропуска паводков такие же, меняются только численные значения максимально допустимых сбросных расходов Qmaxi,Qmax2,Qmax3.

Изложенные правила регулирования стока были исходными для построения модели работы водохранилища и определения его параметров и вероятностных характеристик. Далее, при нормальных условиях стока и заданных начальных правилах эксплуатации водохранилища проводился обычный балансовый расчет пятисуточных наполнений водохранилища согласно уравнению 4. Если в процессе работы водохранилища общий период (за год) пониженного водопотребления был длительнее 5 суток, то такой год считался перебойным. Частота появления перебойных лет по 14

отношению к общему числу лет смоделированной работы водохранилища принималась за вероятность перебоев в работе.

Ра' = Ап/Ы *100 % (5)

где N - общее число лет моделирования; Ап - количество перебойных

лет.

Обеспеченность плановой отдачи воды потребителю вычислялось вычитанием вероятности перебоев из единицы.

Ра=100-Ра' (6)

В представляемой модели фиксировались и случаи полного перебоя. Если в процессе работы водохранилища подача воды в нижний бьеф прекращалась (за исключением санитарного расхода), и, если общий период (за год) такого дефицита был длительнее 10 суток, то фиксировался случай полного перебоя. Частота появления таких лет принималась за вероятность полного перебоя.

Ра" = Апп/Ы (7)

где Апп - количество лет с полными перебоями.

Частота появления аварийных лет по отношению к общему числу лет смоделированной работы водохранилища принималась за вероятность превышения максимального уровня ФПУ.

Величины вероятности появления аварии в зимний (Рз) и летний (Рл) периоды вычислялись отдельно.

Рз=Аз/Ы (8)

Рл=Ал/К (9)

где Рз и Рл - количество лет с зимними и летними авариями соответственно; N - общее число смоделированных лет.

В описываемой модели возможно варьировать величинами (2тахь (Зтах2, (Зтахз, С?тах21, (Зтахг2, (¿тахгз, Углах I, Утах2, Утах3, величиной полезной емкости - Уплз, а также величинами плановой отдачи и координатами трех противоперебойных линий диспетчерского графика.

В четвертой главе анализируется точность расчетов противопаводковой емкости, и излагается разработанная методика оперативного противопаводкового управления водохранилищем комплексного назначения.

В настоящее время традиционный способ расчетов трансформации паводков имеет широкое применение в проектировании объектов водного хозяйства в силу своей узаконенное™ документами СНиП 33-01-2003 и СНиП СП-33-101-03 и сравнительно несложной схемы расчетов. Для анализа возможных ошибок в определении величины противопаводковой емкости при расчетах трансформации паводка традиционным способом использовались длительные искусственные гидрологические ряды пятисуточных значений расхода притока к Краснодарскому водохранилищу.

15

Длительность каждого ряда составляла 1000 лет. Каждый смоделированный 1000-летний ряд разбивался на тридцать одну выборку максимальных среднепентадных расходов воды длительностью по тридцать три года, то есть аналогичной длительности исходного ряда наблюдений. Проведенные расчеты показали, что отклонения расчетных значений противопаводковой емкости от их математического ожидания достаточно значимы. А сами значения противопаводковой емкости могут различаться между собой более чем в два раза при различных выборках с одинаковыми исходными статистическими параметрами.

Таким образом, можно заключить, что при определении противопаводковой емкости традиционным способом целесообразно неоднократную проводить проверку полученных решений, в том числе при возможности - с помощью имитационного стохастического моделирования функционирования водохранилища по длительным искусственным гидрологическим рядам.

Главными особенностями задачи управления водохранилищем комплексного назначения являются стохастичность и многокритериальность. Цели управления водохранилищем часто носят противоречивый характер. В последнее время для усовершенствования методов управления водохранилищем комплексного назначения широко используется модели, использующие как имитационные, так и оптимизационные методы. В диссертационной работе разрабатывалась и использовалась имитационная стохастическая модель функционирования водохранилища с пентадной дискретностью на примере Краснодарского гидроузла, которая включала в себя решение оптимизационной задачи для каждого расчетного периода.

Методика оперативного противопаводкового управления водохранилищем комплексного назначения опирается на оптимизацию сбросов в нижний бьеф путем решение задачи линейного программирования в пределах заданного временного интервала (10 суток). Выбор такого временного интервала обосновывается тем фактом, что генетически обоснованный (метод «осадки-сток») гидрологический прогноз по притоку речных вод к Краснодарскому водохранилищу можно считать практически достоверным только в пределах пяти дней. При большей заблаговременности такого прогноза его вероятность свершения резко понижается. Поэтому в качестве расчетного значения притока за последующую пеитаду принималось прогнозное стохастическое значение, определенное при помощи вспомогательной стохастической модели стока.

Опасными условиями с точки зрения затопления территории были выбраны следующие параметры:

1) 950 млн.м3 < V; < Уфпу

2) №¡+1 > 250 млн.м3

Где: V, - наполнение водохранилища к началу расчетного периода (пентады); №¡+1 - прогнозируемый объем притока воды к водохранилищу.

При выполнении этих условий для каждых последующих двух интервалов, решалась задача линейного программирования, нацеленная на минимизацию затопления территорий верхнего бьефа Краснодарского водохранилища. Целевая функция, предназначенная для безопасного пропуска максимального стока через водохранилище, выражалась следующим образом:

г = (V, -У„пу) + (V*, -Унпу) —> гшп (12)

При этом использовались следующие ограничения:

1) Vi=,( = const

2) V,+i = V; + W; - XI; - X2| - S|

3) Xli = Uli

4) X2i + S; > U2|

5) X2i + Si <Wic6

6) vi+i >VyMO

7) Vi+1 < Унпу + bi

В) Vi+2 = Vi+1 + Wi+I-Xli+l-X2i+1

9) X2i+i + Sj+i > U2j+i

Ю) X2j+i + S+Ü <Wic6

П) Vi+2 >Уумо

12) Vi+2 < Vfmy + bi+I

Где: Vi+1 - наполнение водохранилища к концу расчетного периода (пентады); Vj - наполнение водохранилища к началу расчетного периода (пентады); Wj - объем притока воды к водохранилищу; X1; - объем фактической водоподачи из верхнего бьефа водохранилища включая потери на испарение и фильтрацию; Х2\ - объем фактической водоподачи из нижнего бьефа водохранилища; Sj - объем холостого сброса в нижний бьеф; Wic6 - максимальный разрешенный сброс в нижний бьеф; Uli - объем требуемого водозабора из верхнего бьефа водохранилища; - объем требуемого водозабора из нижнего бьефа водохранилища включая санитарно-экологический попуск; bi - превышение объема водохранилища над объемом НПУ.

Разработанная модель включает в себя 12 ограничений и 4 переменных.

Переменные - объемы холостого заблаговременного сброса воды из водохранилища за два расчетных временных интервала (S, и Sl+i) и объемы водохранилища в концах двух расчетных временных интервалов (V;+i и Vi+2).

Изложенный подход к управлению водохранилищем требует использования методов, приспособленных для практического их применения лицами, принимающими решения (ЛПР), работающими с ограниченным программным обеспечением. Поэтому способы решения поставленных задач опирались на инструментарии информационных технологий, имеющиеся в традиционном офисном обслуживании компьютеров. Разработанная

методика опирается на программное обеспечение MS Excel, имеющееся практически на каждом ПК, и не требующее специализированной подготовки специалистов.

Методика оперативного противопаводкового управления водохранилищем комплексного назначения включает в себя следующие основные этапы:

1) Формирование базы данных значений ежедневного притока к водохранилищу за длительный период.

2) Проведение анализа гидрологического режима притока речных вод к водохранилищу с целью объективного выбора внутригодового интервала для оценки возможности появления аварийной ситуации на гидроузле вследствие прохождения катастрофического паводка, а также для проведения дальнейших водохозяйственных расчетов.

3) Моделирование длительных искусственных гидрологических рядов методом Монте-Карло, включающим способ фрагментов Г.Г. Сванидзе по выбранным внутригодовым интервалам.

4) Разработка вспомогательной модели стока для получения значений стохастических прогнозных величин стока.

5) Разработка имитационной модели функционирования водохранилища комплексного назначения, основанной на стохастическом моделировании стока в совокупности с применением оптимизационных методов.

6) Анализ результатов, полученных с помощью имитационной модели, в которой учитываются множество сценариев, связанных с изменением входных переменных, условиями взаимодействия системы и ограничений. При этом принятие окончательных решений остается за специалистом (ЛПР).

На основе представленной методики разработана соответствующая технология управления водохранилищем, которая проверялась на наблюденном и смоделированных стохастических гидрологических рядах пятисуточных значений притока речных вод к Краснодарскому гидроузлу. При этом были получены следующие результаты: вероятность превышения установленного максимального уровня водохранилища ФПУ была понижена с 0,1% до 0,04%. Для самого многоводного наблюденного года (2002) величина форсировки объема была снижена более, чем на 25% с 635 млн. м3 до 442 млн. м3. На рис. 4 и 5 показаны фактический и смоделированный (на основе имитационной модели с учетом оптимизационной оперативной методики управления водохранилищем) режимы работы Краснодарского водохранилища во время прохождения катастрофического паводка 2002 года.

2000 1500 JOOO 500 О

2500

-фпу

нпу

11pilTOW, млн •МЗ

мла

1да. млн.мЗ

1 3

S u

3

Рис. 4 Режим работы Краснодарского водохранилища во время прохождения катастрофического паводка 2002 года.

Рис. 5 Режим работы Краснодарского водохранилища с учетом методики оперативного противопаводкового управления во время прохождения катастрофического паводка 2002 года.

На рисунках 6 и 7 показаны режимы работы Краснодарского гидроузла при пропуске паводка 0,01% обеспеченности, полученные при испытаниях на имитационной модели функционирования водохранилища комплексного назначения, без учета и с учетом методики оптимизационного оперативного управления. На рисунках видно, что максимальное превышение уровня ФПУ при прохождении паводка 0,01% обеспеченности при режиме работы водохранилища с учетом оперативной методики управления сократилось с 520 млн. м3 до 155 млн. м3.

Рис. 6 Режим работы Краснодарского водохранилища при пропуске паводка 0,01% .обеспеченности при помощи традиционных правил регулирования работы водохранилища.

Рис. 7 Режим работы Краснодарского водохранилища при пропуске паводка 0,01% обеспеченности при учете методики оперативного управления водохранилищем.

Также был проанализирован режим водопользования Краснодарского водохранилища, и было выявлено, что применение методики оперативного противопаводкового управления не приводит к снижению обеспеченности плановой отдачи водопотребителям.

На рис. 8 показана частота превышения НПУ на Краснодарском водохранилище за 1000 смоделированных лет при применении методики оперативного управления и без ее применения соответственно. Из анализа рисунка видно, что наибольшую опасность представляет период с 10 июня по

10 июля, так как именно в этот внутригодовой интервал существует значительная вероятность ущерба от аварийной ситуации на гидроузле.

Количество

. ___

НПУ 1 ООО лет

использования

методики сзПри использовании

"методика----

-Шт. в—I т Аз

И

I

В Ё = « В

2 2 £ £ £

§ Ё I & ¡2 2 8 1 5 2

-Л

2 3-

«Л -н

~ 3

Рис. 8 Количество превышений НПУ на Краснодарском гидроузле за 1000 лет.

ЗАКЛЮЧЕНИЕ

1. Анализ режима речного притока к Краснодарскому водохранилищу проводился по многолетним колебаниям, как годового, так и внутригодового распределения стока. Оценка разностной интегральной кривой годового стока основного донора Краснодарского водохранилища - реки Кубань показала наличие двух четко выраженных периодов изменения водности. С 1911 по 1972 гг. наблюдались незначительные колебания различной направленности и с 1973 по 2006 гг. четко прослеживалась тенденция к увеличению значений годового стока. Оценка цикличности и тренда показала, что водность реки Кубань имеет тенденцию к незначительному увеличению за счет появления отдельных значений, значительно превышающих среднемноголетний сток.

2. Внутригодовое распределение стока характеризуется значительной изменчивостью. В течение различных календарных лет в одни и те же фазовые периоды наблюдались как катастрофические паводки, так и ярко

выраженная межень. Например, 8 марта 1980 года приток к водохранилищу составил всего 110 м7с, а 8 марта 2004 года - 1310 м3/с. Этот факт говорит о том, что практически невозможно построить достаточно достоверный типовой гидрограф расходов воды или получить типовое внутригодовое распределение стока для реки Кубань применительно водохозяйственных расчетов.

3. Для разработки математической модели притока речных вод к Краснодарскому водохранилищу было обосновано моделирование длительных искусственных гидрологических рядов методом Монте-Карло, включающего в себя способ фрагментов Г.Г. Сванидзе.

4. Анализ соответствия искусственных смоделированных гидрологических рядов с исходным рядом пяти суточных значений объемов притока к Краснодарскому водохранилищу показал, что расхождения квантилей аналитической и эмпирической кривых обеспеченностей годового стока, а также различия выборочных статистических характеристик внутригодового распределения стока (я„. С„, С5) и коэффициентов корреляции между величинами смежных и несмежных интервалов находились в пределах точности исходной информации. Статистическая оценка однородности, проведенная по критериям Стьюдента и Фишера, подтверждает нуль-гипотезу однородности рассматриваемых рядов. Таким образом, была принята гипотеза о принадлежности исходного наблюденного ряда и смоделированных рядов к одной генеральной совокупности случайных величин.

5. Разработанная стохастическая модель работы водохранилища комплексного назначения, построенная на основе моделирования длительных искусственных гидрологических рядов расходов воды способом фрагментов, дает возможность рассматривать различные сценарии соотношения стока и водоотдачи при различных заданных параметрах водохранилища и определять величину противопаводковой емкости с оценкой возможной погрешности.

6. Разработанная вспомогательная стохастическая модель стока позволяет получать значения стохастического прогноза притока и представляет собой комплекс условных функций обеспеченности пятидневного притока к Краснодарскому водохранилищу с учетом корреляционных связей между смежными величинами.

7. Разработанная методика регулирования паводкового стока позволяет рассчитывать наполнения водохранилища и объемы необходимых сбросов в нижний бьеф, дает возможность снизить вероятность появления аварийной ситуации на гидроузле, уменьшить используемую часть противопаводковой емкости водохранилища и при этом не допустить увеличение дефицитов воды и частоту их появления.

8. Проведенные исследования подтвердили проверяемую научную гипотезу: возможны такие правила регулирования паводкового стока, отличные от традиционных и ранее применявшихся, при которых вероятность аварии на гидроузле уменьшается, а фактическая водоотдача и обеспеченность ее плановой величины остаются стабильны.

Основное положения диссертации опубликованы в следующих работах:

1. Светлов Е.А. Обоснование выбора дискретности при разработке математической модели водохранилища [Текст] / Е.А. Светлов // Материалы международной научно-практической конференции «Роль Природообустройства сельских территорий в обеспечении устойчивого развития АПК». Часть 1. -М.: МГУП, 2007. - С. 277-281.

2. Светлов Е.А. Технология повышения обеспеченности орошения из водохранилищ комплексного назначения [Текст] / В.В. Ильинич, Е.А. Светлов, H.JI. Жморщук, Д. Бутутау // Сборник трудов международной научно-практической конференции «Агротехнологии 21 века». -М.: МСХА им. К.А. Тимирязева, 2007. - С. 96-98.

3. Светлов Е.А. Стохастическая модель водохранилища [Текст] / В.В. Ильинич, Е.А. Светлов // Природообустройство. - 2008. - №4. - С. 59-64.

4. Светлов Е.А. Оценка однородности годового стока реки Кубань [Текст] / Е.А. Светлов // Мелиорация и водное хозяйство. - 2009. - №1. -С. 19-20. (ВАК)

5. Светлов Е.А. Анализ методики расчета противопаводковой емкости водохранилища [Текст] / В.В. Ильинич, Е.А. Светлов, НЛ. Жморщук // Гидротехническое строительство. - 2009. - №9. - С. 46-48. (ВАК)

6. Светлов Е.А. Управление водохранилищем комплексного назначения в период паводков [Текст] / В.В. Ильинич, Е.А. Светлов, A.B. Ильинко // Природообустройство. -2009. - №3. - С. 58-60.

Подписано в печать:

14.10.2009

Заказ № 2728 Тираж -100 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

Введение.

Глава 1. Современное состояние вопроса определения противопаводковой емкости (трансформации паводка) водохранилища ирригационно-противопаводкового назначения.

1.1 Существующие методы расчета противопаводковой емкости.

1.2. Область применения разрабатываемой методики и выбор объекта исследования.1,

Глава 2. Анализ и оценка гидрологических условий речного бассейна Краснодарского водохранилища.

2.1. Основные характеристики объекта исследования.

2.2. Условия формирования стока рек, питающих

Краснодарское водохранилище.

2.3. Анализ и оценка многолетних колебаний притока речных вод в Краснодарское водохранилище

2.4. Анализ и оценка внутригодового распределения притока речных вод в Краснодарское водохранилище.

Глава 3. Стохастическое моделирование процесса речного стока и режима работы водохранилища.

3.1. Выбор дискретности моделирования.

3.2. Стохастическая модель притока речного стока в водохранилище комплексного использования.

3.3. Верификация модели притока относительно точности исходной информации.ч.

3.4. Вспомогательная модель паводкового стока.

3.5. Стохастическая модель функционирования водохранилища комплексного использования.

Глава 4. Управление Краснодарским водохранилищем с учетом стохастического моделирования.

4.1 . Исследование традиционного способа регулирования паводкового стока водохранилищем комплексного использования.

4. 2. Методика оперативного противопаводкового управления водохранилищем комплексного использования.

4.3 Результаты применения методики оперативного противопаводкового управления водохранилищем комплексного использования на имитационной модели.

Актуальность темы диссертации. Во многих странах мира и в большинстве субъектов Российской Федерации безопасность жизнедеятельности населения в значительной степени зависит от катастрофических дождевых паводков. Ежегодно от таких паводков гибнут люди, и наносится существенный ущерб экономике. В России ярким примером этого являлся дождевой паводок, произошедший в 2002 году в бассейне реки Кубань. В результате паводка было подтоплено 130 населенных пунктов, разрушено около 70 тыс. домов, пострадало более 130 тыс. человек, в том числе 96 человек - с летальным исходом. Общий ущерб от паводка составил 8,65 млрд. руб. Одной из эффективных мер противодействия разрушительной силе паводков является противопаводковое регулирование стока водохранилищем. В настоящее время при создании и эксплуатации водохранилищ выделяется специальная противопаводковая емкость, предназначенная для задержания избытков речного стока, таким образом, любое водохранилище должно являться либо комплексным, либо специально противопаводковым.

Задача установления рационального режима работы водохранилища комплексного назначения, весьма сложна из-за противоречивого и неоднозначного характера его использования. В частности, при определенных режимах речного стока и водопользования, характерных для юга России, в различные по водности годы в одноименные фазовые периоды могут происходить как засухи, так и катастрофические паводки. В этих случаях, для целей орошения и водоснабжения в маловодные периоды водохранилище должно находиться в наполненном состоянии, а для задержания максимальных расходов катастрофических паводков должно быть достаточно опорожнено. Указанное противоречие может быть отчасти устранено при заблаговременных и точных прогнозах о сроках прохождения паводков, их объемов и характере гидрографа. Однако современный уровень прогнозирования стока не может в полной мере обеспечить водное хозяйство соответствующими прогнозами.

Этим обуславливается целесообразность разработки максимально безопасных правил управления водохранилищами при использовании детальной оценки имеющейся информации о гидрологическом режиме водоисточника. Соответственно, в настоящей работе формируется методика противопаводкового регулирования стока на основе анализа формирования притока речных вод к водохранилищу, а также режима требований водопользователей.

Цель и задачи исследования. Основной целью диссертационной работы является разработка методики определения оперативных правил управления водохранилищем комплексного назначения в условиях паводкового режима, позволяющих достичь понижения вероятности аварии при регулировании максимального стока. Для достижения поставленной цели потребовалось решение комплекса взаимосвязанных задач:

• Анализ и оценка режима речного притока к водохранилищу;

• Разработка математической модели притока речных вод к водохранилищу;

• Разработка модели функционирования водохранилища комплексного назначения;

• Разработка методики определения оперативных правил регулирования максимального стока водохранилищем комплексного назначения с учетом краткосрочных гидрологических прогнозов;

Методика и объект исследования. Общим методическим положением являются современные подходы к математическому моделированию режимов стока рек и функционирования водохранилища комплексного назначения. В качестве математического аппарата используется принцип имитационного стохастического моделирования с применением методов линейного программирования. В работе используются основные положения статистического анализа временных рядов.

Объектом исследования является Краснодарское водохранилище и сток питающих его рек, характеризующихся прохождением катастрофических паводков в период интенсивного водопользования.

Научная новизна работы состоит в следующем:

• представлена методика формирования безопасных и рациональных правил оперативного управления водохранилищем на основе имитационного моделирования функционирования водохранилища с использованием методов линейного программирования;

• представлен вероятностный процесс стока по типу Марковской цепи с пентадной дискретностью;

• получены новые режимные характеристики Краснодарского водохранилища по выделенным внутригодовым интервалам; Научная гипотеза, проверяемая в диссертационной работе, заключается в следующем: возможны такие правила регулирования паводкового стока, отличные от традиционных и ранее применявшихся, при которых: а) вероятность аварии на гидроузле будет уменьшена; б) фактическая водоотдача и обеспеченность ее плановой величины будут стабильны при неизменных параметрах водохранилища;

Практическая ценность работы. На основе полученных результатов разработана методика оперативного управления трансформации паводка водохранилищем с учетом краткосрочных гидрологических прогнозов. Данная методика представляет собой инструмент для моделирования режимов функционирования водохранилищ, в различных сценариях изменяющихся лриродно-хозяйственных условиях. Предложенная методика может быть использована для построения диспетчерских правил управления водными ресурсами водохранилища Доступный инструментарий методики позволяет применение ее непосредственно на объекте без использования специального программного обеспечения. Применение разработанной методики позволяет уменьшить вероятность аварийной ситуации при прохождении катастрофических паводков. Результаты исследований приняты в Кубанском бассейном водном управлении для рассмотрения о внедрении в работу оперативных служб.

Достоверность полученных результатов подтверждается сравнением результатов моделирования и материалов наблюдений, а также использованием современных подходов имитационного моделирования.

Апробация работы. Основные результаты диссертационной работы докладывались:

1) на заседании Кубанского бассейнового водного управлении и были приняты для рассмотрения о возможности внедрении в работу оперативных служб КБВУ (август 2009 г.),

2) на международной научно-практической конференции «Агротехнологии 21 века» (2007 г.),

3) на научно-технических конференциях МГУП (2006-2009 гг.),

3) на заседаниях кафедры гидрологии, метеорологии и регулирования стока (2006-2009 гг.).

Публикации. По теме диссертации опубликованы 6 печатных работ, две из которых опубликованы в журналах, входящим в список ВАК.

Структура и объём диссертации. Диссертация состоит из введения, четырех глав, заключения и списка используемой литературы. Работа изложена на 142 страницах машинописного текста, иллюстрированного 25 рисунками, и содержит 16 таблиц. Список используемой литературы состоит из 125 наименований.

ЗАКЛЮЧЕНИЕ

1. Водность основного донора Краснодарского водохранилища — реки Кубань имеет тенденцию к незначительному увеличению за счет появления отдельных значительно отличающихся от среднемноголетнего стока его годовых величин.

2. Внутригодовое распределение стока характеризуется значительной изменчивостью. В течение различных календарных лет в одни и те же фазовые периоды наблюдались как катастрофические паводки, так и ярко выраженная межень.

3. Метод Монте-Карло в реализации его способом» фрагментов. Г.Г. Сванидзе достаточно корректен, позволяет с большой- точностью моделировать искусственные гидрологические ряды расходов и может быть использован при расчетах трансформации- максимального стока водохранилищем комплексного назначения.

4. Стохастическая модель работы водохранилища комплексного назначения, построенная на основе моделирования длительных искусственных гидрологических рядов расходов воды способом фрагментов, дает возможность рассматривать различные сценарии соотношения стока и водоотдачи при различных заданных параметрах водохранилища и определять величину противопаводковой емкости с оценкой возможной погрешности.

5. При выборе рациональных правил оперативного регулирования стока в паводковый период следует учитывать корреляционные зависимости между смежными величинами стока внутригодовых интервалов.

6. Традиционный способ определения, противопаводковой емкости водохранилищ на реках с длительным паводковым режимом; основанный на трансформации построенного; гидрографа согласно требованиям СНиП. 3301-2003, не является достаточно точным и корректным, исключает объективную оценку возможной погрешности точности полученных результатов с точностью исходной информации.

7. Проектирование основных параметров водохранилища комплексного назначения целесообразно проводить с учетом будущих правил его эксплуатации.

8. При выборе рациональных правил оперативного регулирования стока в паводковый период следует учитывать корреляционные зависимости между смежными величинами стока внутригодовых интервалов.

9. Разработанная методика регулирования паводкового стока дает возможность уменьшить вероятность появления аварийной ситуации на гидроузле и уменьшить величину противопаводковой емкости водохранилища.

10. Предложенная методика регулирования стока водохранилищем комплексного назначения не приводит к увеличению дефицитов воды и частоты их появления.

1. Абакян А.Б. Влияние наводнений на жизнь общества и защита от них // Абакян А.Б., Плюшкин А.А. - Изв. АН СССР. - 1989. - №2. - С. 41-54.

2. Авакян, А.Б. Водохранилища // А.Б. Авакян, В.П. Салтанкин, В .А. Шарапов. М.: Мысль, 1987. - 320 с.

3. Алексеев Г.А: Расчеты паводочного стока рек СССР // Г.А. Алексеев. -Л.: Гидрометеоиздат, 1955.

4. Алексеев Г.А. Методы определения случайных погрешностей гидрометеорологической информации // Г. А. Алексеев. — Л.: Гидрометеоиздат, 1975. 94 с.

5. Андерсон Т. Статистический анализ временных рядов // Т. Андерсон. — М.: Мир, 1976. 755 с.

6. Андреянов B.F. Внутригодовое распределение стока // В.Г. Андреянов. Л.': Гидрометеоиздат, 1960. - 327 с.

7. Апполов Б.А. Гидрологические прогнозы,// Б.А. Апполов, Г.П. Калинин, В.Д. Комаров. Л.: Гидрометеоиздат, 1974. - 419 с.

8. Асарин А.Е. Методика определения максимального половодного стока к Среднеенисейской ГЭС // А.Е. Асарин, М.И. Джейбо, М.И. Львова. -Гидротехническое строительство. 1985. — №7. — С. 29-31.

9. Асарин А.Е. Водно-энергетические расчеты // А.Е. Асарин, К.М. Бестужева. М.: Энергия, 1986. - 223 с.

10. Бабурин Б.Л. Общий случай расчета- противопаводочной емкости в водохранилище гидроэлектростанции // Б.Л. Бабурин. Тр. Гидропроекта. -1964.-т. 12.-С. 200-210.

11. Бабурин Б.Л. Расчет резервной емкости в водохранилищах ГЭС // Б.Л. Бабурин. Гидротехническое строительство. - 1963. - №1. - С. 40-43.

12. Бек Н.Н. Статистические методы оптимизации в экономических исследованиях // Н.Н. Бек, Д.И. Голенко. -М.: Статистика, 1971. 135 с.

13. Беллиндер Е.Н. Основные положения и методика определения критериев безопасности ГТС // Е.Н. Беллиндер, И.Н. Иващенко, В:Д. Новоженин. — Безопасность энергетических сооружений. — 2000. — Вып.7.

14. Блохинов Е.Г. Распределение вероятностей речного стока // Е.Г. Блохинов. -М.: Наука, 1974. 169 с.

15. Бусалаев И.В. Сложные водохозяйственные системы // И.В. Бусалаев. — Алма-Ата: Наука, 1980.-231 с.

16. Бусленко Н.П. Моделирование сложных систем // Н.П. Бусленко. — М.: Наука, 1978.-339 с.

17. Важное А.Н. Анализ. и прогноз рек-Кавказа // А.Н. Важнов М.: Гидрометеоиздат, 1966. — 274 с.

18. Вилкас Э.И. Векторные критерии оптимизации. // Э.И. Вилкас // Математика и кибернетика в экономике. — М.: Экономика, 1971. — С. 20-221

19. Водноэнергетические расчеты методом Монте-Карло // Под ред. А.Ш. Резниковского. -М;: Энергия, 1969. — 304 с.

20. Воропаев Г.В. Моделирование водохозяйственных систем аридной зоны СССР // Г.В. Воропаев, F.X. Исмайылов, В.М. Федоров. М.: Наука, 1984.-312 с.

21. Вуглшский B.C. Водные ресурсы и водный баланс крупных водохранилищ СССР // B.C. Вуглинский. JL: Гидрометеоиздат, 1991. -223 с.

22. Гаврилов В.М. О влиянии разрезки годового гидрографа речного стока на результаты его регулирования // В.М. Гаврилов. — Сб. Труды МГМИ. — вып. №2.

23. Гаврилов В.М. Учет вероятностной сущности речного стока при диспетчеризации работы водохранилища // В.М. Гаврилов. — Сб. Известия АН Армянской СССР. Ереван: Технические науки.

24. Галямин Е.П. Оптимизация оперативного распределения водных ресурсов в орошении // Е.П. Галямин. JL: Гидрометеоиздат, 1981. — 271 с.

25. Гидрогеологическое прогнозирование // Под ред. М.Г. Андерсона. — М.: Мир, 1988.

26. Гриневич Г.А. Композиционное моделирование гидрографов // Г. А. Гриневич, Н.А. Петелина, А.Г. Гриневич. М.: Наука, 1972. - 181 с.

27. Дружинин И.П. Динамика многолетних колебаний речного стока // И.П. Дружинин, В.Р. Смага, А.Н. Шевнин. -М.: Наука, 1991. 175 с.

28. Дюбин Г.Н. Введение в прикладную теорию игр // Г.Н. Дюбин, В.Г. Суздаль М.: Наука, 1981. - 336 с.

29. Ефимович П. А. Вопросы водохозяйственных расчётов и гидрологии // П.А. Ефимович М.: Изд-во ОНТИ, 1935. - 320 с.

30. Железняк И.А. Регулирование паводочного стока // И.А. Железняк — JL: Гидрометеоиздат, 1965. 326 с.

31. Железняков Г.В. Точность гидрометрических измерений и расчетов // Г.В. Железняков, Б.Б. Данилевич. Л.: Гидрометеоиздат, 1966. -240 с.

32. Железняков Г.В. Теория гидрометрии // Г.В. Железняков. JL: Гидрометеоиздат, 1976. - 344 с.

33. Железняков Г.В. Гидрология, гидрометрия и регулирование стока // Г.В. Железняков, Т.А. Неговская, Е.Е. Овчаров. -М.: Колос, 1984. 439 с.

34. Жиркевич А.Н. Современные подходы к определению характеристик максимального стока (на примере проекта гидроузла Шонла на р. Да во

35. Вьетнаме) // А.Н. Жиркевич // Юбилейный сборник научных трудов Гидропроекта. -2000:-Выпуск 159.

36. Зайков Б.Д. Средний сток и его распределение в году на территории Кавказа // Б.Д. Зайков. Л.': Гидрометеоиздат, 1946. - 64 с.

37. Ильинич В.В. Имитационная модель работы водохранилища // В.В. Ильинич. Экспресс информация, 1982,— С. 1-4.

38. Ильинич В.В. Стохастическая модель водохранилища // В:В. Ильинич, Е.А. Светлов. — Природообустройство. 2008. - №4. - С. 59 -64.

39. Исмайылов Г.Х. Расчет паводков водохранилищем // Г.Х. Исмайылов. // Материалы докладов научно-технического семинара-совещания. Паводки-и защита от них. Хабаровск. - 1968.

40. Исмайылов Г.Х. Вероятностные расчеты водохранилищ комплексного (ирригационно-противопаводочного) назначения // Г.Х. Исмайылов // -Водные ресурсы. 1973. - №2. - С. 81-94.

41. Исмайылов» Г.Х. Формализация1 гидрологических особенностей в моделях управления водохозяйственными' системами // Г.Х. Исмайылов; К.Ю. Шаталова // Сб. Водные проблемы на рубеже веков. — М.: Наука, 1999. С. 279-290.

42. Казакевич Д.И. Основы теории случайных, функций, в задачах гидрометеорологии // Д.И. Казакевич. Л.: Гидрометеоиздат, 1989: - 230 с.

43. Калинин Г.П. Основы методики краткосрочных прогнозов водного режима // Г.П. Калинин. Тр. ЦИП, 1952: - Вып. 28 (55).

44. Картвелишвили Н.А. Теория вероятностных процессов в гидрологии и регулировании паводочного стока, // Н.А. Картвелишвили. — Л.: Гидрометеоиздат, 1967. — 292 с.

45. Картвелишвили Н.А. Регулирование речного стока // Н;А. Картвелишвили. Л.: Гидрометеоиздат, 1970.-218 с.

46. Картвелишвили Н.А. Стохастическая гидрология // Н.А. Картвелишвили. Л.: Гидрометеоиздат, 1981. -167 с.

47. Кучмент Л.С. Математическое, моделирование речного стока // Л.С. Кучмент. — Л.: Гидрометеоиздат, 1972. -191 с.

48. Кучмент Л.С. Модели формирования речного стока // JI.C. Кучмент. -J1.: Гидрометеоиздат, 1980.- 142 с.

49. Кучмент Л.С. Гидрологическое прогнозирование для управления водноресурсными системами// Л.С. Кучмент. М.: ВИНИТИ, 1981. — Т.4.53; Кучмент Л.С. Формирование речного стока // JI:C. Кучмент, ВЖДемидов, Ю.Г. Мотовилов^ — Mi: Наука, 19831 216 с.

50. Лурье ШМ: Водные ресурсы и водный баланс Кавказа// П;М; Лурье. -СП;: Гидрометеоиздат, 2002.-500с.55: Лурье П.М. Река Кубань. Гидрография и режим стока // П.М. Лурье, ВД. Панов, Ю.Ю. Ткаченко. СП1: Гидрометеоиздат, 2005;

51. Математическое моделирование в управлении водными ресурсами. -М.: Наука, 1988.-246 с.

52. Методика определения критериев безопасности гидротехнических сооружений И РД 153-34.2-21.342-00 РАО «ЕЭС России». -2001.

53. Митропольский А.К. Техника статистических вычислений // А.К. Митропольский — М.: Наука, 1971. 576 с.

54. Моисеев Н.Н. Математические задачи системного анализа // Н.Н. Моисеев. -М.: Наука, 1981. -487 с.

55. Музылев С.В. Стохастические модели в инженерной гидрологии // С.В. Музылев, В.Е. Привальский, Д.Я. Раткович. М.: Наука, 1962. - 184 с.

56. Обоснование стратегий управления водными ресурсами— М;: Научный мир, 2006. — 336 с.

57. Плешков Я.Ф. Регулирование речного стока // Я.Ф. Плешков — Л.: Гидрометеоиздат, 1972. — 560 с.

58. Подольский Е.М. Определение противопаводочной емкости с учетом сезонного регулирование стока // Е.М. Подольский, Б.Л. Бабурин. -Гидротехническое строительство. — 1960. —№6. — С. 30—35.

59. Потапов М.В. Сочинения // М.В. Потапов. М.: Сельхозгис, 1951. -т.З.

60. Пособие по определению расчетных гидрологических характеристик. -Л.: Гидрометеоиздат, 1984. — 447 с.

61. Правила эксплуатации Краснодарского водохранилища на р. Кубани. Книга 1. Гидрологическое обоснование. — Краснодар. —1990;

62. Правила эксплуатации- Краснодарского водохранилища. — Рабочий проект. Гидрологическое обоснование. — Рукопись. Краснодар. — 2001.

63. Пряжинская В.Г. Математическое моделирование в водном хозяйстве // В.Г. Пряжинская. М.: Наука, 1985. - 113 с.

64. Раткович Д.Я. Стохастические модели колебаний составляющих водного баланса речного бассейна // Д.Я. Раткович, М.В. Болгов. М. - 1997.

65. Раткович Д.Я: Многолетние колебания речного стока // Д.Я. Раткович. — Л.: Гидрометеоиздат, 1976. 255 с.

66. Раткович Д.Я. Математические модели колебаний речного стока // Д.Я. Раткович // Моделирование водохозяйственных систем. М.: ИВП РАН, 1992.-350 с.

67. Раткович Д.Я. Гидрологические основы водообеспечения // Д.Я. Раткович М.: ИВП РАН, 1993. - 428 с.

68. Расчеты паводочного стока. Методы расчетов на основе мирового опыта // Под ред. А.А. Соколова. Л:: Гидрометеоиздат, 1978. - 304 с.

69. Регулирование и комплексное использование водных ресурсов // Под ред. Г.В. Воропаева. М:: Наука, 1987. - 157 с.

70. Резниковский. А.Ш. Гидрологические основы гидроэнергетики // А.Ш. Резниковский. М.: Энергия, 1978. — 147 с.

71. Резниковский А.Ш: Диспетчерские правила управления режимами водохранилищ // А.Ш. Резниковский, М.И. Рубинштейн. — М.: Энергоатомиздат, 1984. 157 с.

72. Руководство пользователя' и техническое описание программного продукта «Mike 11» — 2008. — 58 с.

73. Ржаницин Н:А: Речная гидравлика и гидротехника // Н.А. Ржаницин. -М.: Наука, 1973.

74. Рождественский А.В. Статистические методы в гидрологии // А.В. Рождественский, А.И. Чеботарёв. Л.: Гидрометеоиздат, 1974. — 424 с.

75. Рождественский А.В. Оценка точности кривых распределения гидрологических характеристик // А.В. Рождественский. — Л.: Гидрометеоиздат, 1977. — 268 с.

76. Рождественский А.В. Оценка тонности гидрологических расчетов // А.В. Рождественский, А.В, Ежов, А.В, Сахарюк. Л:: Гидрометеоиздат, 1990. -275 с.

77. Саваренский А.Д. Регулирование речного стока водохранилищами // А.Д. Саваренский. М.: Изд-во АН СССР, 1951. - 251 с.

78. Самборский Т.В. Многолетние колебания стока рек с паводочным режимом: Автореферат диссертации на соискание ученой степени кандидата географических наук / Т.В. Самборский. М., 1997. - 23 с.

79. Светлов Е.А. Оценка однородности-годового стока реки Кубань // Е.А. Светлов. Мелиорация и водное хозяйство. — 2009. -№1. - С. 19-20;

80. Соколовский Д:Л: Речной сток // ДШ: Соколовский' Л.: Гидрометеоиздат, 1968. — 539 с.93; Статистические методы в гидрологии. Л.: Гидрометеоиздат, 1970. -270 с.

81. Строительные Нормы и> Правила 33-01-2003 РФ. Гидротехнические сооружения. — М., 2003.95; Строительные Нормы и-Правила СП-33-101-03. — Л.: Гидрометеоиздат, 2003.

82. Схема комплексного использования и охраны водных ресурсов реки Кубань // Гидрологическое обоснование. ОАО Проектно-изыскательный институт «КубаньВодПроект». Краснодар, 2003.

83. Федеральные указания по безопасности плотин. США, 1972.

84. Хранович И.Л. Управление водными ресурсами, потоковые модели // И.Л. Хранович. М.: Научный мир, 2001. - 295с.

85. Христофоров А.В. Надежность расчётов речного стока // А.В. Христофоров. — М.: Изд-во МГУ, 1993. 168 с.

86. Христофоров А.В. Теория случайных процессов в гидрологии //

87. A.В. Христофоров. М.: Изд-во МГУ, 1994. - 139 с.

88. Христофоров А.В. Стохастическая модель колебаний речного стока в паводочный период // А.В. Христофоров Г.В. Круглов, Т.В. Самборский. М.: Изд-во МГУ, 1998. - 146 с.

89. Христофоров А.В. Стохастическая модель колебаний паводочного стока рек Северного Кавказа // А.В. Христофоров, Н.М. Юмина // Водное хозяйство России. Проблемы, технологии, управление. 2005. -№4. - Т.7 — С. 359-377.

90. Цингер В.Н. Трансформация максимальных расходов водохранилищами // В.Н. Цингер. Л.: Гидрометеоиздат, 1960. - 121 с.

91. Чокин Ш.Ч. Методика расчета регулирования стока при комплексном использовании водных ресурсов // Ш.Ч. Чокин, В.А.Григорьев,

92. B.К. Редышн. Алма-Ата: Наука, 1977. - 300 с.

93. Шелутко В.А. Статистические модели и методы исследования многолетних колебаний стока // В.А. Шелутко. — JL: Гидрометеоиздат, 1984. 159 с.

94. Шенгелия Г.П. К вопросу о регулировании паводочных гидрографов методом фрагментов. // Г.П. Шенгелия. Сообщения АН ГрузССР. - 1972. - №2. - т.68. - С. 365-368.

95. Шикломанов И. А. Оценка изменений стока р. Кубани под влиянием хозяйственной деятельности // И.А. Шикломанов, А.С. Хорецкая. Труды ГГИ. - 1975. - вып. 229: - С. 55-71.

96. Юмина Н.М. Оценка стока рек Кавказа: Автореферат диссертации на соискание ученой степени кандидата географических наук / Н.М. Юмина. М., 2008. - 28 с.

97. Dawdy D.K. Application of systems analysis to network design //

98. D.K. Dawdy, M.E. Moss, N.C. Matalas. Cacebook of Hydrilogical.

99. Fahlbuch E.F. Spillway design floods and dam safety //

100. E.F. Fahlbuch. Hydropower & Dam. - USA. - Vol. Six. - Issue Four. - 1999.

101. Ilinich V.V. Search of Anti-Accident Function for Flood Flow Management by Water Reservoir // V.V. Ilinich. Proceedings of the 7-th International Conference on Hydroinformatics. — Nica. - 2006. — p.1025 — 1031.

102. Kovacs Z.P. Problems in the estimation of maximum flood peaks // Z.P. Kovacs. South African National Hydrological Simposium, - Pretoria. -1983.

103. Mialhe F. Synoptic assessment of water resource variability in reservoirs by remote sensing: General- approach and application to the runoff harvesting systems of south India // F. Mialhe, Y. Gunnell, C. Mering. — India. — 2008.

104. Montanari. A. Estimating the uncertainty of hydrological forecasts: A statistical approach // A. Montanari, G. Grossi. 2008.

105. Perelman. L. Cross Entropy multiobjective optimization for water distribution systems design // L. Perelman, A. Ostfeld, E. Salomons. 2008.

106. Smith, Т. J. Bayesian methods in hydrologic modeling: A study of recent advancements in Markov chain Monte Carlo techniques // T. J. Smith, L. A. Marshall.-2008.

107. Texas Water Develop. Board. A Completion Report on System Simulation for Management of Total Water Resources. Texas. Austin. — 1969. — P. 430.

108. Federal Energy Regulatory Commission Office of Hydropower Licensing. USA. - 1993.

109. Vrugt J. A. Treatment of input uncertainty in hydrologic modeling: Doing hydrology backward with Markov chain Monte Carlo simulation // J.A. Vrugt, J.F. ter Braak, M.P.Clark, J.M. Hyman, B.A. Robinson. 2008.

110. Wei. C. Derived operating rules for a reservoir operation system: Comparison of decision trees, neural decision trees and fuzzy decision trees // C. Wei, N. Hsu. 2008.

111. You. J. Hedging rule for reservoir operations I I J. You, X. Cai. 2008.

112. Young G. K. Fielding reservoir operating rules. Hydraulics Div. ASCE. - 1967. - Vol. 93. - P. 41-56.

113. Федеральное агентство водных ресурсов (Росводресурсы)

114. Кубанское бассейновое водное управление350020, г. Краснодар, ул. Красная, 180-а Тел.(8б1) 253-73-07: Факс(861) 253-73-05 E-mail: kbvu@mail. kuban. г и

115. От Ш ОУ- 0$ № Of-bj^^ на от

116. В специализированный ученый совет Московского Государственного Университета Природообустройства1. СПРАВКА

117. Заместитель руководителя Кубанского бассейнового вод, управления1. Белан Е.В.