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

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

□03054366 МУРАЩЕНКОВА Наталья Владимировна

ОЦЕНКА ХАРАКТЕРИСТИК РЕЧНОГО СТОКА В ИЗМЕНЯЮЩИХСЯ ПРИРОДНО-ХОЗЯЙСТВЕННЫХ УСЛОВИЯХ

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

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

МОСКВА 2007

003054966

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

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

заслуженный деятель науки РФ Исмайылов Г.Х.

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

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

кандидат технических наук, доцент Раткович Л.Д.

Ведущая организация - ЗАО ПО «СОВИНТЕРВОД»

Защита состоится 23 апреля 2007г. в на заседании диссертационного

совета Д 220 045 02 в Московском государственном университете природообустройства по адресу 127550, Москва, ул Прянишникова, д 19, аудитория zoi

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

Автореферат разослан «20» марта 2007 г

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

/ ИМ Евдокимова

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

- анализ многолетней изменчивости годовых величин речного стока и их учет при совершенствовании методов оценки количественных и качественных характеристик речного стока;

- исследование процесса трансформации речного стока по длине реки,

- оценка изменения показателей качества воды в устьевой области реки;

- построение зависимостей изменения речного стока и его минерализации с учётом стохастического характера данных процессов

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

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

- результаты анализа многолетних колебаний годового стока рек Волги и Дона, позволившие установить, что временным рядам годового стока этих рек в целом свойственна нестационарность на достаточно длительном отрезке времени (100 лет и более) и квазистационарность на отдельных отрезках рассматриваемых временных рядов (30-40 лет),

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

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

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

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

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

Апробация работы. Основные результаты диссертационной работы докладывались на 3й международной Российско-Иранской конференции «Сельское хозяйство и природные ресурсы» (2002 г), на научно-технических конференциях МГУП (2004-2006 гг), на заседаниях кафедры гидрологии, метеорологии и регулирования стока, а также изложены в ряде публикаций

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

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

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

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

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

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

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

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

О масштабах влияния хозяйственной деятельности на речной сток можно судить по объёмам общего и безвозвратного водопотребления. В 80-е годы XX века водопотребление достигло своего максимального значения, например, в бассейне Волги составило 35-37 км3/год, Дона - 20-22 км3/год, при чём безвозвратное водопотребление составило соответственно - 14 и 8 км3/год С начала 90-х годов ХХ-го столетия в связи с падением промышленного и сельскохозяйственного производства наметилась тенденция к снижению водопотребления, а, следовательно, и безвозвратного расхода По данным 2000 г водопотребление в бассейне Волги сократилось до 21,8 км3/год, Дона - до 6 км3/год, при этом безвозвратный расход уменьшился соответственно до 5 и 2,2 км3/год.

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

Кроме того, зарегулирование стока рек Волги и Дона изменило естественный водный режим пойм и устьевых областей В результате создания водохранилищ Волжско-Камского каскада речной сток уменьшился на 5-6 км3, что вызвало трансформацию экосистем Волго-Ахтубинской поймы и дельты Волги Снижение стока р Дон привело к снижению притока в Таганрогский залив и Азовское море, что неблагоприятно сказалось на изменении солености воды, а, следовательно, и на их продуктивности

Анализ многолетних колебаний годового стока р Волги в створе г Волгограда выполнен по материалам наблюдений за стоком воды за период 1881/82 - 1994/95 гг (п = 114 лет) При этом были использованы две версии исследуемого временного ряда фактическая (наблюденная) и условно-естественная (данные института Гидропроект) Для 1881/82 — 1933/34 гг объёмы годового стока для обеих версий совпадают, а для 1934/35 — 1994/95 гг объёмы стока в первой версии отражают интегральное влияние антропогенных факторов на годовой сток

см см см см см см см см см см см см

со О) о ^ со ю со С: со оз

т»

СО о о см со тг ю со со о

СО со о О) о> о о О) о о ст> О)

т- т— т- г— т*

Рис 1 Разностные интегральные кривые условно-естественного (1) и наблюденного (2) рядов годового стока р Волги

При сопоставлении разностных интегральных кривых условно-естественного и наблюденного годового стока проявляется влияние хозяйственной деятельности на сток (рис 1 ) Интегральное снижение годового стока р Волги за последние 60 лет составило около 1040 км3 или 17 км3/ год, те около 7% от среднемноголетнего условно-естественного стока Столь незначительное снижение годового стока свидетельствует о том, что основную роль в формировании стока Волги играют природные факторы и, прежде всего климатические

Важной особенностью процесса многолетних колебаний условно-естественного годового стока Волги является его цикличность, т е последовательное чередование лет (йли групп лет) с различной водностью В частности можно выделить период снижения годового стока р Волги, охватывающий 12 лет (1926/27 - 1937/38 гг), с изменением стока с 385 до 159 км3/год и период увеличения стока с 1938/39 по 1950/51 гг. (11 лет) со 159 до 333 км3/год, то есть, имеем полный цикл колебаний годового стока длительностью 22 года Аналогично можно выделить и другие циклы, что позволяет говорить о приближении характера колебаний к циклическому

В таблице 1 приведены статистические параметры годового стока р Волги за 1881/82 - 1994/95 гг, которые оценивались для различных отрезков исследуемых временных рядов Сравнение оценок выборочных средних и дисперсий исследуемых рядов с использованием критериев Стьюдента (г ) и Фишера (Р) показало, что исследуемые ряды можно считать однородными, так как Т =1,66 меньше ткр =1,96, при уровне значимости а =0,05 и Р=1,05 меньше Ркр=1,35 при а =0,05 Но при оценке степени однородности отрезков временных рядов для

периода 1934/35-1994/95 гг., оказалось, что ряд фактического годового стока в отличие от условно-естественного ряда, в результате воздействия антропогенных факторов приобретает черты нестационарности в части математического ожидания. Проверка гипотезы о независимости годового стока с помощью критерия Андерсона показала, что она должна быть отвергнута как для исходных рядов, так и для выделенных их отрезков и, таким образом, годовые объёмы стока р Волги не являются совокупностью независимых случайных величин

Таблица 1 Основные статистические параметры годового стока р Волги у г. Волгограда за 1881/82-1994/95гг

Период, число лет Статистические параметры

км3/год СГЖ ,км3/год Су Сэ ЖО

1 Наблюденный ряд

1881/82- 1994/95 гг. 247 44 0,18 2 Су 0,406

1881/82- 1933/34 гг 258 44 0,17 2 Су 0,345

1934/35 -1994/95 гг 239 42,6 0,18 2 Су 0,394

2 Условно-естественный ряд

1881/82- 1994/95 гг 257 45,1 0,18 2 Су 0,431

1881/82- 1933/34 гг 258 44,0 0,17 2 Су 0,345

1934/35 -1994/95 гг 256 46,5 0,18 2 Су 0,490

Кроме того, проведен анализ степени однородности выборочных оценок статистических параметров (среднего, стандарта, коэффициента автокорреляции) для скользящих 30-летних отрезков условно-естественного ряда и для периодов с различным типом атмосферной циркуляции (западной, восточной и меридиональной)

Проверка гипотезы об однородности средного и стандарта скользящих 30-летних отрезков ряда с помощью критериев однофакторного дисперсионного анализа и Бартлета показала, что при 5% уровне значимости расхождение в оценках среднего и стандарта может быть признано статистически недостоверным В тоже время существуют различия между оценками коэффициента автокорреляции Следовательно, анализируемый ряд является неоднородным с точки зрения корреляции между его смежными членами

Проверка гипотезы об однородности среднего и стандарта с помощью критериев однофакторного дисперсионного анализа и Бартлета для периодов, различающихся типом атмосферной циркуляции, показала, что условно-естественный ряд годового стока Волги за период 1881/82 - 1994/95 гг является неоднородным по математическому ожиданию и коэффициенту автокорреляции

Порядковый номер 30-летки

Рис. 2 Динамика скользящих 30-летних значений статистических параметров условно-естественного годового стока р Волги у г Волгограда за

1881/82-1994/95 гг

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

W{t) = Wcp + R(\)[W(t -1 )-Wcp] + m (1)

A

где fV(t) и Wit — 1) - объёмы годового стока в t и предшествующий ему ( / — 1) -ый годы, Wcp - среднемноголетний объём годового стока, R(l) - коэффициент автокорреляции, - гауссовский «белый шум» с нулевым средним

<Ts=<rw4\ -ад2

Для анализа многолетних колебаний годового стока р Дон использованы данные Государственного океанографического института (ГОИН) по наблюденному и условно-естественному стоку Дона за период с 1881/82 по 1999/00 гг в замыкающем отборе бассейна р Дон - ст-ца Раздорская

Для периода с 1881/82 по 1935/36 гг данные по наблюденному и условно-естественному стоку совпадают, а с 1936/37 по 1999/00 гг первый - отражает влияние антропогенных факторов, тогда как второй - влияние климатических факторов на сток

Рис.3 Разностные интегральные кривые условно-естественного (1) и наблюденного (2) рядов годового стока р Дон

Совместный анализ разностных интегральных кривых наблюденного и условно-естественного стока р Дон (рис 3) показывает, что интегральное снижение годового стока Дона за 60 последних лет составило 347 км3, или в среднем 5,8 км3/год, те 21% среднемноголетнего условно-естественного стока Столь явное снижение годового стока свидетельствует о том, что основную роль в формировании стока р Дон играет антропогенный фактор

Колебания годового стока р Дон имеют циклический характер -выделяется фаза пониженного стока с 1900/01 гг. по 1914/15 гг и повышенного - с 1915/16 гг. по 1932/33 гг, то есть имеется полный цикл годового стока длительностью 23 года Аналогично можно выделить и другие циклы

В таблице 2 приведены статистические параметры годового стока р Дон за 1881/82 — 1999/00 гг, которые оценивались для различных отрезков исследуемых временных рядов Оценка степени однородности отрезков исследуемых рядов (условно-естественного и наблюденного) для периода с 1952/53 по 1999/00 гг. показала, что ряды можно считать однородными с точки зрения дисперсии и неоднородными по математическому ожиданию

Аналогично выполнена оценка степени однородности отрезков наблюденного ряда с 1881/82 по 1951/52 гг и с 1952/53 по 1999/00 гг При сравнении оценок дисперсий по критерию Фишера получено, что фактическое значение Б =1,65 больше критического^ =1,58, при а-5%, а также при

сравнении оценок среднего по критерию Стьюдента г =3,32 больше Тк =1,98, следовательно, различие между ними статистически достоверно

Таблица 2 Основные статистические параметры годового стока р Дона у ст-цы Раздорская за 1881/82-1999/00 гг

Период, число лет Статистические параметры

\Уср, км3/год <УШ, км3/год Су Сэ 11(1)

1 Наблюденный ряд

1881/82-1999/00 гг 25,11 9,96 0,40 2 Су 0,16

1881/82- 1951/52 гг 27,50 11,0 0,40 2 Су 0,08

1952/53 -1999/00 гг 21,58 6,66 0,31 2 Су 0,17

2 Условно-естественный ряд

1881/82-1999/00 гг 27,32 9,44 0,35 2 Су 0,06

1881/82-1951/52 гг. 27,66 11,0 0,40 2 Су 0,09

1952/53 -1999/00 гг 26,81 6,35 0,24 2 Су 0,12

-0,2

1 6 11 16 21 26 31 36 41 46 51 56 61 66 71 76 81 86 91 Порядковый номер 30-летки

Рис 4 Динамика скользящих 30-летних значений статистических параметров условно-естественного годового стока р, Дон у ст-цы Раздорская за 1881/82-1999/00 гг

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

Проверка гипотезы об однородности среднего и стандарта скользящих 30-летних отрезков условно-естественного временного ряда годового стока р Дон позволила признать ряд однородным В тоже время различия между оценками коэффициента автокорреляции вполне очевидно (рис 4)

Анализируя автокорреляционную (АКФ) и частную автокорреляционную (ЧАКФ) функции условно-естественного ряда годового стока р Дон имеем, что первая имеет выброс при сдвиге по времени равном 1, а ЧАКФ экспоненциально колебательно затухает с изменением знака. Считаем, что рассматриваемый процесс условно-естественного стока может быть идентифицирован моделью авторегрессии проинтегрированного скользящего среднего АРПСС (0,1,1) с параметром ©i = 0,628.

Проведенная оценка степени однородности основных статистических характеристик годового стока реки Волги показала наличие изменений в их динамике с точки зрения математического ожидания и коэффициента корреляции между смежными членами ряда, и для годового стока реки Дона — в части математического ожидания и дисперсии

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

При исследовании многофакторных гидрологических процессов возникает необходимость установления зависимости между несколькими переменными Исходными данными для получения такой зависимости служат материалы наблюдений над переменой у и определяющими её переменными Xi,X2, ., Xk Результаты совместных наблюдений можно представить в виде {у,; xtl, х2„ , xkl}, i=l,2, , N

Уравнение множественной линейной регрессии имеет следующий вид

У, = а о + аххи + сс2х21 +... + акхь + Z,, (2)

где а0,ОС1,...,ОСк - параметры данного уравнения, которые оцениваются методом

наименьших квадратов, Z, - случайная величина, представляющая ошибку, причем M(Z,)=0. Ошибка связана как с неточностью измерения, так и с влиянием на нее неучтенных факторов

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

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

Определение кубической сплайн-функции

Пусть на отрезке [а, б] оси X имеется сетка узлов (Оы = {а = x{(x2(...(xN — b}, тогда кубическая сплайн-функция на каждом из отрезков [х(, х)+1 ] представляет собой полином третьей степени, сшитый в узлах X,, таким образом, чтобы в узлах сохранилась непрерывность функции S3(x) вплоть до второй производной

Представлением кубической сплайн-функции является выражение

S3(x) = a, +b,(x-xl) + cl(x-xl)2 +dl{x-xlf (3)

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

Ь N

/jCg)= Jg"wr^ = min и ^(gOc,)-^)2^, (4)

а '=1

где R > 0 - заданное число, называемое параметром сглаживания, R - параметр сглаживания, представляющий собой коэффициент автокорреляции первого порядка, подсчитанный по ряду остатков

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

синхронных измерений величин Уг 5-^2, построить функцию F. В

МОО, как и в методе множественной линейной регрессии функция ищется в виде

к

F(xl,...,xk) = YiafJ(xJ) + €. (5)

7=1

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

На первом шаге МОО в качестве первой переменной X, выбирается та, которая обладает наиболее тесной статистической связью с у, Далее по значениям

томощью кубического сплайна строится аппроксимирующая функция

влиянием на переменную у, оставшихся неучтенными к — 1 переменными На

вычисляются новые отклонения Ду2, = уи — /2Хх2^и т Д Процесс продолжается до тех пор, пока не будут учтены все рассматриваемые факторы, либо эти факторы перестанут быть значимыми Значимость факторов оценивается по коэффициенту парной линейной корреляции, если связь линейная В качестве меры тесноты статистической связи произвольной формы между двумя переменными служит коэффициент автокорреляции

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

Участок реки Волги ниже Волгоградского гидроузла занимает территория Волго-Ахтубинской поймы протяженностью более 450 км и шириной от 15 до 35 км На этом участке Волга протекает по засушливой зоне и не имеет притоков Широко развитый микрорельеф, наличие огромных пойменных пространств, обширной дельты с множеством рукавов в сочетании с климатическими условиями обуславливают на этом участке бассейна Волги существенное перераспределение стока и его потери К настоящему времени на данном участке сформировался сложный природно-хозяйственный комплекс, существование которого определяется попуском из водохранилища и его дальнейшей трансформацией при движении по руслу

Расчет трансформации речного стока основан на данных наблюдений за среднемесячными уровнями воды в нижнем бьефе Волгоградской ГЭС, створах Черный Яр и Верхнее Лебяжье (вершина дельты) за период с 1971 г по 1984 г, который включает в себя годы различной обеспеченности, в том числе исключительно маловодный 1975г. (р=97%) и исключительно многоводный 1979 г (р=4,3%)

этом полученные отклонения Ау1; — JV, ~ f\t (-^i,) обусловлены

следующем шаге МОО строится аппроксимирующая функция /2(^2)и

В результате была построена зависимость между среднемесячным уровнем воды в створе Волгоградской ГЭС, как главным определяющим фактором, и искомым среднемесячным уровнем воды в русле Волги - створе Черный Яр (рис 5) в виде кубического сплайна. Аналогично построена зависимость между среднемесячным уровнем воды в створе Волгоградской ГЭС, как главным определяющим фактором, и искомым уровнем воды в русле Волги - створе Верхнее Лебяжье

О 200 400 600 800 1000

100 300 500 700 900

Н фактический Н расчётный

Уровень воды в н б Волгоградской ГЭС, см

Рис 5 Зависимость между уровнем воды в н б Волгоградской ГЭС и уровнем воды в створе по руслу Волги — Верхнее Лебяжье

Полученный коэффициент корреляции между фактическими и расчётными значениями равен 0,942 (участок, нижний бьеф Волгоградской ГЭС - Чёрный Яр) и 0,845 (участок нижний бьеф Волгоградской ГЭС - Верхнее Лебяжье), показывает, что корреляционная связь достаточно тесная и свидетельствует об устойчивости и эффективности полученного уравнения

Используя полученную зависимость между уровнем воды в нижнем бьефе Волгоградской ГЭС и уровнем воды в характерных створах русла р Волги -Черный Яр и Верхнее Лебяжье, рассчитывались трансформированные гидрографы по руслу р Волги в характерных створах (рис 6) В качестве исходных приняты расчетные уровни воды в нижнем бьефе Волгоградского гидроузла, полученные в результате численной реализации имитационной модели Нижневолжской водноресурсной системы, разработанной на кафедре гидрологии, метеорологии и регулирования стока МГУП

На рисунке 6 представлены фактические и трансформированные гидрографы в створе нижнего бьефа Волгоградской ГЭС и в характерном створе по руслу р. Волги - Верхнее Лебяжье (вершина дельты), в разные по водности годы При сравнении фактических и трансформированных гидрографов видно неплохое соответствие.

1971 - 1972 г г , сут

Рис 6 Фактические (сплошные кривые) и расчетные (штриховые) гидрографы Волги по створам н б Волгоградской ГЭС (1) и Верхнее Лебяжье (2)

Акватория Куйбышес "кого водохранилища

Куйбышевское водохранилище - одно из крупнейших долинных водохранилищ и наиб», нее мощная ступень в Волжско-Камском каскаде водохранилищ Водохрани шще представляет собой ряд озеровидных расширений, соединенных между собой > зкостями, вытянуто в меридианальном направлении по р. Волге и в северо-восточно 1 направлении по р Каме и имеет длину по судовому ходу между Волжской, Чебоксарской и Нижнекамской ГЭС соответственно 480 и 540 км На Куйбышевском водохранилище осуществляется сезонное, недельное и суточное регулирование стока вследствие большого объёма (58 км3) оно является основным регулятором стока на Средней Волге

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

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

створах в Куйбышевское водохранилище (по руслу р Волги и р Камы) Для получения данной зависимости использован метод кубических сплайн-функций в сочетании с методом остаточных отклонений

В качестве входных на акватории Куйбышевского водохранилища выбраны участки г Чебоксары - с. Вязовые (по руслу Волги), г. Набережные Челны (по руслу Камы) - с Сокольи Горы, г. Вятские Поляны (по руслу Вятки) - с Сокольи Горы со среднемесячными уровнями воды соответственно Х1, Х3, Х2 Замыкающим - выходным является участок г Старая Майна — г Тольятти Этот участок охватывает 40% акватории водохранилища Расчёты для построения зависимости велись по среднемесячным уровням воды за период с 1958 г (после наполнения водохранилища) по 1980 г. (200 членов ряда), каждый четвертый элемент массива проверочный (76 членов)

В результате построена искомая зависимость в виде-

У = -0,96795 +1,01948*,! +1,076Ш12, (6)

/л =в| -хи)+с,(хх -хь)2 +алхх -хьу

112=а, +ЬХХЪ-Х3,) + с,{Х з-Х3,)2+<1,(Х3-Х31)3

где У - уровень воды у г Тольятти,

/, - преобразованные переменные, ? = 1, ..,п, ] = 1,2.

месяцы

Рис 7 Фактические (1) и расчетные (2) уровни воды р Волги у г Тольятти

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

Коэффициент множественной корреляции равен 11=0,875 с интервалом 0,842-0,907 По проверочному массиву построены фактические и расчётные уровни воды у г. Тольятти (рис 7) Анализ расчётов по проверочному массиву показывает, что данная методика дает относительно устойчивые результаты Коэффициент корреляции между расчётным и фактическим значением равен 0,77

Оценка изменения солености воды на акватории Таганрогского залива Азовского моря. Таганрогский залив представляет собой мелководный полузамкнутый водный объект (24,6 км3), сформированный на стыке двух крупных природных систем — реки Дон и Азовского моря Протяженность вытянутого по длине залива составляет 140 км. Большую часть площади залива (около 75%) занимают глубины менее 3 м. В результате мелководности залива, а, следовательно, хорошего прогрева водной массы и малой солености воды в нём, Таганрогский залив является благоприятным местом обитания и размножения ценных пород рыб

Соленость воды залива определяется как влиянием стока р Дон, так и притока соленых вод из Азовского моря. Если сток Дона в течение года значительно изменяется лишь в весенне-летний период (март-июнь), оставаясь относительно устойчивым в остальные сезоны года, то морской приток зависит от ветровой активности При нагонных ветрах юго-западного и западного направления морская вода поступает в залив, повышая соленость его вод При сильных сгонных ветрах и при стоке р Дон во время половодья опреснение вод распространяется до морской границы залива, а иногда выходит за его пределы Дон обеспечивает около 70 % от общего речного стока в море, поэтому именно воды Дона оказывают наибольшее влияние на солевой состав вод Таганрогского залива

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

Для расчета использованы среднемесячные данные наблюдений прибрежных морских гидрологических станций (МГ) Таганрог и МГ Ейск за период с 1966 г по 2000 г, характеризующие соответственно восточную и западную части залива и среднемесячные расходы воды р Дон в замыкающем створе - ст-ца Раздорская.

Расчёт среднемесячной солености воды Таганрогского залива для МГ Ейск.

В качестве входных (независимых характеристик) выбраны следующие факторы: среднемесячный расход воды р Дон за предшествующий месяц (Х0, температура воды (Х2), ветровой параметр, определяющий сгонно-нагонные явления залива (Х3), соленость воды предшествующего месяца, определяющая начальные условия (Х4) Среднемесячные расходы воды р Дон выбраны с учетом времени добегания распресняющего влияния стока, которое для МГ Ейска равно двум месяцам

В качестве выходной характеристики выступает соленость воды Таганрогского залива для МГ Ейск в текущем месяце ( у )

Расчёт осуществлялся по-шагово, начиная с наиболее значимой переменной определенной по коэффициенту парной линейной корреляции и коэффициенту автокорреляции

Уравнение искомой зависимости будет иметь вид-

^ = 0,44 + (0,941 ± 0,0458)/; (Х4) + (0,860 ± 0Д496ЩХ,) + (1,186 ± 0Д736)/з(^Гз) + +(1,208 ± 0,1787ЩХ2) (7)

где ^(Х4), Г2(Х]), ^(Хз), Г4(Х2) - кубические сплайн-функции

Переменная Х3, то есть ветровой параметр, определяющий сгонно-нагонные явления, взятый в интервале дискретности один месяц, для данного расчета оказался статистически не значим

По полученной зависимости на рис 8 приведены фактические и расчетные значения среднемесячной солености воды Таганрогского залива для МГ Ейск

Коэффициент корреляции между наблюденными и расчётными значениями равен 0,89 с доверительными интервалами 0,858 и 0,916 Точность полученных результатов неплохо согласуется с фактическими данными наблюдений по МГ Ейск (отклонение не превышает 25-30%) Полученное уравнение также проверено на независимом материале (ряд из 66 членов) Коэффициент корреляции между наблюденными и расчетными значениями оказался равным 0,86

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

Рис 8 Фактическая (1) и расчетная (2) среднемесячная соленость воды Таганрогского залива для МГ Ейск за период 1966 по 2000 гг.

Расчет среднемесячной солености воды Таганрогского залива для МГ Таганрога выполнен аналогично предыдущему расчету Установлено, что соленость воды имеет наиболее тесную положительную корреляцию с соленостью воды предыдущего периода (11=0,78), отрицательную - с речным стоком р. Дон (Ы = -0,34) Ветровой параметр, определяющий сгонно-нагонные явления, взятый в интервале дискретности один месяц для данного расчета оказался статистически не значим

Фазу понижения солености вйды, вызванную половодным стоком р Дон, характеризует среднемесячная соленость июня (у) В качестве предикторов для уравнения множественной регрессии выбраны суммарный сток р Дон за половодье ((Зпол) с учетом времени добегания распресняющего влияния стока, соленость воды, наблюдающаяся до влияния стока половодья (8П0Л) (среднемесячная соленость марта для МГ Таганрог) и ветровой параметр (Ур)

Общее уравнение искомой зависимости солености воды для МГ Таганрога в июне I - года будет иметь вид

уг = 1,89+0,3 88пол-0,000бдпол+0,0003Ур (8)

Фактические и расчетные значения среднемесячной солености июня для МГ Таганрог за период с 1966 гг по 2000 гг представлены на рис 9 Коэффициент множественной корреляции равен 0,74.

Рис 9 Фактическая (1) и расчетная (2) среднемесячная соленость воды для МГ Таганрог, июнь

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

Теоретической основой метода композиции является обобщенный метод, разработанный С Н Крицким, М Ф. Менкелем и А Д Саваренским

Если водохранилище работает в нормальном режиме, то есть без сбросов и дефицита, уравнение водного баланса водохранилища имеет вид

х = у + 2-ит-ий+а> (Уто<х<Ун/7У\* = 0 , ¿ = 0 (9)

где X и у - соответственно конечное и начальное наполнение водохранилища,

г - приток воды к водохранилищу; в - холостой сброс из водохранилища,

(1 - дефицит полезных отдач, 1/т - полезная водоотдача из водохранилища,

и. - потери воды из водохранилища (испарение + фильтрация), IУ - объём сточных вод поступающих в водохранилище при минерализации

сточных вод

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

г+ю-и,

Г тт тт . ,Л

^ г-ит-Цн+а>

г-и„-С1„+а>

-V у-£-, (10)

г+й)-С/л

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

При постоянстве параметров IIт, IIп , (О, ¡ла в расчётном отрезке времени Т минерализация воды водохранилища /их будет зависеть от начального наполнения,

расчётного притока к водохранилищу и соответствующей им минерализации

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

рыг = т,Р<.У,)'Р,;", (и)

то

гдеР(л:()г, />[//(,)]г - полная вероятность /-го интервала соответственно

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

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

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

Полезная емкость Кайраккумского водохранилища была разделена на семь

групп (интервалов) наполнения - от Ушо=1,2 км3 до Уипу =3,8 км3 Для каждой

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

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

После построения условных функций распределения конечных наполнений водохранилища и соответствующих им минерализаций вычисляются ординаты безусловных функций распределения их вероятностей по формулам (11) и (12) (рис 10,11)

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

Рис 10 Безусловные функции распределения вероятностей конечного наполнения

Кайраккумского водохранилища, июнь 1 - с учетом ограничений по минерализации, 2 - без учёта ограничений

Рис 11 Безусловные функции распределения вероятностей минерализации конечного наполнения Кайраккумского водохранилища, апрель-июнь (IV - VI)

Заключение

1. Анализ пространственно-временных закономерностей изменчивости характеристик речного стока Волги и Дона, особенно за последние 20-30 лет, показывает наличие статистически достоверных трендов в их динамике, обусловленных как изменением регионального климата, так и антропогенным воздействием на водные ресурсы. В связи с этим требуется усовершенствование методов анализа и оценки характеристик речного стока в изменяющихся природно-хозяйственных условиях

2 Анализ многолетних колебаний годового стока р Волги по степени однородности за период с 1881/82 по 1994/95 гг показал наличие изменений в динамике годового стока, обусловленных как климатическими, так и антропогенными факторами Анализируемый ряд годового стока р Волги можно считать неоднородным с точки зрения математического ожидания и коэффициента корреляции между смежными членами ряда Стационарность данного процесса имеет место лишь на тех отрезках исследуемого ряда, где преобладают климатические факторы над антропогенными Для описания многолетних колебаний условно-естественного стока р Волги принята модель авторегрессии первого и второго порядка

3. При анализе многолетних колебаний годового стока р Дон установлено, что исследуемый ряд условно-естественного стока является однородным с точки

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

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

5 Проведенные исследования по оценке трансформации речного стока по длине реки на отдельных участках русла реки Волги показали перспективность использования метода кубических сплайнов и метода остаточных отклонений в имитационных моделях функционирования водохозяйственных систем речных бассейнов Получены аналитические зависимости, описывающие трансформацию речного стока, как между узлами управления, так и в пределах водохранилищ применительно к водохозяйственной системе речного бассейна р Волги В частности, для Куйбышевского водохранилища построена зависимость, связывающая уровень воды в замыкающем створе с уровнями воды на входных створах по руслам рек Волги и Камы Данная зависимость показала, что среднемесячный уровень воды на приплотинном участке водохранилища достаточно точно определяется средними уровнями воды на входных русловых участках рек Волги и Камы Для участка р Волги ниже Волгоградского гидроузла получены трансформированные гидрографы в характерных створах р Волги -Черный Яр и Верхнее Лебяжье, определяющие рациональный водный режим Волго-Ахтубинской поймы в интересах рыбного и сельского хозяйства

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

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

8. Для описания изменения речного стока и его минерализации предложен метод композиции в сочетании с балансовым методом Данный подход реализован для Кайраккумского водохранилища на р Сырдарье. Построены условные и безусловные функции распределения вероятностей конечного наполнения водохранилища и его минерализации и соответственно установлена связь между минерализацией и водностью реки (притока к водохранилищу) и минерализацией и наполнением водохранилища

9 Проведенные исследования показали уникальность пространственно-временной изменчивости характеристик речного стока как Волги, так и Дона. Залогом дальнейшего развития научных исследований по данному направлению и совершенствования моделей водноресурсных систем речных бассейнов является обеспечение непрерывных рядов наблюдений по существующим постам и расширение наблюдательной сети на реках России

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

1 Account of transformation of a flood at movement it on the river bed on the basis of spline-functions Proceedings of The 3rd International Iran and Russia Conference «Agriculture and Natural Resources». Moscow, 2002. P. 252 (В соавторстве).

2. Расчет трансформации паводка при движении его по руслу реки на основе сплайн-функций //Проблемы научного обеспечения развития эколого-экономического потенциала России Сборник научных трудов. М ■ МГУП, 2004 С. 77-80

3. Анализ многолетних колебаний годового стока р. Дон //Природообустойство и рациональное природопользование - необходимые условия социально -экономического развития России Сборник научных трудов Часть 1. - М . МГУП, 2005 С. 131-135 (В соавторстве)

4. Применение композиционного метода для определения минерализации воды водохранилища //Материалы международной научно-практической конференции «Роль природообустройства в обеспечении устойчивого функционирования и развития экосистем» Часть II.-М. МГУП, 2006. С 133-138.

5 Вероятностный расчет минерализации воды при регулировании стока водохранилищем // Мелиорация и водное хозяйство, 2006. №5. С. 56-58.

г

Московский государственный университет природообустройства (МГУП)

Зак. № iZti Тираж 100

Введение.

Глава 1. Современное состояние вопроса определения характеристик речного стока.

1.1 Существующие методы оценки характеристик речного стока

1.1.1. Оценка многолетних колебаний речного стока.

1.1.2. Оценка трансформации руслового стока.

1.2. Влияние антропогенных факторов на изменение характеристик речного стока.

1.3. Влияние возможного изменения климата на характеристики речного стока.

Глава 2. Анализ и оценка многолетних колебаний речного стока.

2.1. Анализ и оценка многолетних колебаний годового стока р. Волги.

2.2. Анализ и оценка многолетних колебаний годового стока р. Дон.

Глава 3. Методика оценки основных характеристик речного стока.

3.1. Метод множественной линейной регрессии.

3.2. Метод сплайн-функций.

3.3. Метод остаточных отклонений (МОО).

3.4. Расчёт трансформации речного стока на участке р. Волги.

3.5. Расчёт солености воды Таганрогского залива.

3.5.1. Характеристика Таганрогского залива.

3.5.2. Изменчивость солености воды Таганрогского залива.

3.5.3. Расчёт солености воды Таганрогского залива методом множественной линейной регрессии.

3.5.4. Расчёт солености воды Таганрогского залива методом сплайн-функций и методом остаточных отклонений.

Глава 4. Композиционный метод определения основных характеристик речного стока.

4.1. Постановка задачи.

4.2. Объект исследований.

4.3. Расчёт и построение кривых обеспеченностей стока и минерализации воды в Кайраккумском водохранилище.

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

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

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

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

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

Для достижения поставленной цели потребовалось решение комплекса взаимосвязанных задач:

- анализ многолетней изменчивости годовых величин речного стока и их учёт при совершенствовании методов оценки количественных и качественных характеристик речного стока;

- исследование процесса трансформации речного стока по длине реки;

- оценка изменения показателей качества воды в устьевой области реки;

- построение зависимостей изменения речного стока и его минерализации с учётом стохастического характера данных процессов.

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

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

- результаты анализа многолетних колебаний годового стока рек Волги и Дона, позволившие установить, что временным рядам годового стока этих рек в целом свойственна нестационарность на достаточно длительном отрезке времени (100 лет и более) и квазистационарность на отдельных отрезках рассматриваемых временных рядов (30-40 лет);

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

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

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

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

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

Апробация работы. Основные результаты диссертационной работы докладывались на 3й международной Российско-Иранской конференции «Сельское хозяйство и природные ресурсы» (2002 г.), на научно-технических конференциях МГУП (2004-2006 гг.), на заседаниях кафедры гидрологии, метеорологии и регулирования стока, а также изложены в ряде публикаций.

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

ЗАКЛЮЧЕНИЕ

1. Анализ пространственно-временных закономерностей изменчивости характеристик речного стока Волги и Дона, особенно за последние 20-30 лет, показывает наличие статистически достоверных трендов в их динамике, обусловленных как изменением регионального климата, так и антропогенным воздействием на водные ресурсы. В связи с этим требуется усовершенствование методов анализа и оценки характеристик речного стока в изменяющихся природно-хозяйственных условиях.

2. Анализ многолетних колебаний годового стока р. Волги по степени однородности за период с 1881/82 по 1994/95 гг. показал наличие изменений в динамике годового стока, обусловленных как климатическими, так и антропогенными факторами. Анализируемый ряд годового стока р. Волги можно считать неоднородным с точки зрения математического ожидания и коэффициента корреляции между смежными членами ряда. Стационарность данного процесса имеет место лишь на тех отрезках исследуемого ряда, где преобладают климатические факторы над антропогенными. Для описания многолетних колебаний условно-естественного стока р. Волги принята модель авторегрессии первого и второго порядка.

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

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

5. Проведенные исследования по оценке трансформации речного стока по длине реки на отдельных участках русла реки Волги показали перспективность использования метода кубических сплайнов и метода остаточных отклонений в имитационных моделях функционирования водохозяйственных систем речных бассейнов. Получены аналитические зависимости, описывающие трансформацию речного стока, как между узлами управления, так и в пределах водохранилищ применительно к водохозяйственной системе речного бассейна р. Волги. В частности, для Куйбышевского водохранилища построена зависимость, связывающая уровень воды в замыкающем створе с уровнями воды на входных створах по руслам рек Волги и Камы. Данная зависимость показала, что среднемесячный уровень воды на приплотинном участке водохранилища достаточно точно определяется средними уровнями воды на входных русловых участках рек Волги и Камы. Для участка р. Волги ниже Волгоградского гидроузла получены трансформированные гидрографы в характерных створах р. Волги - Чёрный Яр и Верхнее Лебяжье, определяющие рациональный водный режим Волго-Ахтубинской поймы в интересах рыбного и сельского хозяйства.

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

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

8. Для описания изменения речного стока и его минерализации предложен метод композиции в сочетании с балансовым методом. Данный подход реализован для Кайраккумского водохранилища на р. Сырдарье. Построены условные и безусловные функции распределения вероятностей конечного наполнения водохранилища и его минерализации и соответственно установлена связь между минерализацией и водностью реки (притока к водохранилищу) и минерализацией и наполнением водохранилища.

9. Проведенные исследования показали уникальность пространственно-временной изменчивости характеристик речного стока как Волги, так и Дона. Залогом дальнейшего развития научных исследований по данному направлению и совершенствованию моделей водноресурсных систем речных бассейнов является обеспечение непрерывных рядов наблюдений по существующим постам и расширение наблюдательной сети на реках России.

1. Алберг Дж., Нильсон Э., Уолш Дж. Теория сплайнов и её применения. М.: Мир, 1972.316 с.

2. Анисимов О.А. Последствия изменения климата в регионах северного полушария // В кн.: Изменения климата и их последствия. СПб.: Наука, 2002. С.239-250.

3. Антропогенные воздействия на водные ресурсы России и сопредельных государств в конце XX столетия / Отв. ред.: Н.И. Коронкевич, И.С. Зайцева. М.: Наука, 2003. 367 с.

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

5. Бокс Дж., Дженкинс Г. Анализ временных рядов. Прогноз и управление. Вып. 1., 2., М.: Мир, 1974.

6. Боровкова Т.Н., Никулин П.И. Куйбышевское водохранилище (краткая физико-географическая характеристика). Куйбышевское книжное издательство, 1972.

7. Бронфман A.M., Хлебников Е.П. Азовское море: основы реконструкции. JL: Гидрометеоиздат, 1985. 271 с.

8. Вагер Б.Г., Серков Н.К. Сплайны при решении прикладных задач метеорологии и гидрологии. JI.: Гидрометеоиздат, 1987. 160 с.

9. Водогрецкий В.Е. Влияние агромелиораций на годовой сток. JL: Гидрометеоиздат, 184 с.

10. Вода России: Речные бассейны / Под ред. Черняева А. ФГУП РосНИИВХ. Екатеринбург: АКВА-ПРЕСС, 2000. 536 с.

11. Воды России (состояние, использование, охрана) 1996-2000 гг. Екатеринбург: РосНИИВХ, 2002. 254 с.

12. Воинов А.А., Комилов Ф.С. Имитационная модель Кайраккумского водохранилища. М.: Вычислительный центр АН СССР, 1986. 39 с.

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

14. Воропаев Г.В., Исмайылов Г.Х., Федоров В.М. Проблемы управления водными ресурсами Арало-Каспийского региона. М.: Наука, 2003.427 с.

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

16. Вучков И., Бояджиева Д., Солаков Е. Прикладной линейный регрессионный анализ. М.: Финансы и статистика. 1987. 155 с.

17. Гаргопа Ю.М. Современное распреснение Азовского моря и его связь с многолетними колебаниями атмосферной циркуляции // Водные ресурсы. 2002. Т.29. №6. С.747-754.

18. Георгиевский В.Ю., Цыценко К.В., Шалыгин A.JI. Оценка притока поверхностных вод в Каспийское море // В кн.: Гидрометеорологические аспекты проблемы Каспийского моря и его бассейна. СПб.: Гидрометеоиздат, 2003. С. 217-229.

19. Георгиевский В.Ю. Изменение стока рек России и водного баланса Каспийского моря под влиянием хозяйственной деятельности и глобального потепления: Автореферат диссертации на соискание ученой степени доктора географических наук. СПб., ГТИ, 2005. 55 с.

20. Геткер М.И., Куропатка JI.M., Рубинова Ф.Э. Об изменении общей минерализации воды р.Сырдарьи в связи с развитием орошения в её бассейне // Тр. САРНИГМИ. 1975. Вып. 23 (104).

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

22. Гидрометеорологический режим озер и водохранилищ СССР. Водохранилища Средней Азии / Под ред. Ю.Н. Иванова. Д.: Гидрометеоиздат, 1991. 163 с.

23. Голубев B.C., Цыценко К.В. Структура и динамика потерь речного стока в низовьях Волги // Метеорология и гидрология. 1995. №8. С. 85-93.

24. Груза Г.В., Ранькова Э.Я. Обнаружение изменений климата: состояние, изменчивость и экстремальность климата // Метеорология и гидрология. 2004. №4. С. 50-66.

25. Грушевский М.С. Неустановившееся движение воды в реках и каналах. Л.: Гидрометеоиздат, 1982. 288 с.

26. Дрейпер Н., Смит Г. Прикладной регрессионный анализ: В 2-х кн. М.: Финансы и статистика. 1986. Т.1. 366 е.; 1987. Т.2. 351 с.

27. Дружинин И.П., Коваленко З.П. Речной сток и геофизические процессы (связи, цикличность и предвидение). М.: Наука, 1966.

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

29. Езекиэл М., Фокс К. Методы анализа корреляций и регрессий линейных и криволинейных. М.: Статистика, 1966. 557 с.

30. Ефимова Н.А., Строкина JI.A., Байкова И.М. Изменение основных элементов климата на территории СССР в 1967- 1990 гг. // Метеорология и гидрология. 1996. №4. С. 34-41.

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

32. Завьялов Ю.С., Квасов Б.И., Мирошниченко B.JI. Методы сплайн-функций. М.: Наука, 1980. 352 с.

33. Зажимаров С.А. Сравнительный анализ и область применения упрощенных методов расчёта неустановившегося движения речных потоков // Моделирование и прогнозы гидрологических процессов. СПб.: РГГМУ, 1999. Вып. 120. С. 58-65.

34. Закономерности экосистемных процессов в Азовском море / Под ред. Г.Г. Матишова. М.: Наука, 2006. 304 с.

35. Знаменский В.А., Чигиринский П.Ф. Куйбышевское водохранилище. JL: Гидрометеоиздат, 1978.

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

37. Исмайылов Г.Х. Федоров В.М. Анализ многолетних колебаний годового стока Волги // Водные ресурсы. 2001. №5. С.517-525.

38. Исмайылов Г.Х., Шаталова К.Ю. Учет трансформации руслового стока в имитационной модели Волжской ВРС // В кн.: Г.В. Воропаев, Г.Х. Исмайылов, В.М. Федоров Проблемы управления водными ресурсами Арало-Каспийского региона. М.: Наука, 2003. С. 386-392.

39. Исмайылов Г.Х., Шаталова К.Ю. Исследование возможности применения сплайнов для формализации гидрологических процессов при решении водохозяйственных задач // Водные ресурсы. 2003. Т.30. №2. С. 245-249.

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

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

42. Калинин Г.П., Милюков П.И. Приближенный расчёт неустановившегося движения водных масс. Тр. ЦИП. 1958. Вып. 66.

43. Карасев И.Ф. Трансформация расходов воды и русловой водный балпнс при неустановившемся движении речных потоков // Метеорология и гидрология. 1991. №2. С. 84-91.

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

45. Константинов А.Р., Химин Н.М. Применение сплайнов и метода остаточных отклонений в гидрометеорологии. JI.: Гидрометеоиздат, 1983. 184 с.

46. Коронкевич Н.И. Водный баланс Российской равнины и его антропогенные изменения. М.: Наука, 1990. 205 с.

47. Коронкевич Н.И., Барабанова Е.А., Зайцева И.С., Малик JI.K. Водность года по речному стоку и хозяйственная деятельность. // Мелиорация и водное хозяйство. 2002. №3. С. 10-12.

48. Крицкий С.Н., Менкель М.Ф. Водохозяйственные расчёты. JL: Гидрометеоиздат, 1952. 392 с.

49. Кукса В.И. Южные моря (Аральское, Каспийское, Азовское и Чёрное) в условиях антропогенного стресса. СПб.: Гидрометеоиздат, 1994. 320 с.

50. Кукса В.И., Гаргопа Ю.М. Современная оценка гидрологических условий формирования биопродуктивности Азовского моря // Водные ресурсы. 2004. Т.31. №4. С. 489-497.

51. Кучмент JI.C. Математическое моделирование речного стока. JL: Гидрометеоиздат, 1972. 191 с.

52. Кучмент JI.C. Модели формирования речного стока. JL: Гидрометеоиздат, 1980. 142 с.

53. Кучмент JI.C. Гидрологическое прогнозирование для управления водноресурсными системами. М.: ВИНИТИ, 1981. Т.4.

54. Кучмент JI.C., Демидов В.Н., Мотовилов Ю.Г. Формирование речного стока. М.: Наука, 1983. 216 с.

55. Лисицына Л.В. Особенности солевого состава вод Азовского моря. Автореферат диссертации на соискание ученой степени кандидата географических наук. М, 1997. 21 с.

56. Лобанов В.А. Расчёты речного стока в неоднородных и нестационарных условиях: Автореферат диссертации на соискание ученой степени доктора технических наук. СПб., ГГИ, 1998. 40 с.

57. Международное руководство по методам расчёта основных гидрологических характеристик. Л.: Гидрометеоиздат, 1984. 248 с.

58. Митропольский А.К. Техника статистических вычислений. М.: Наука, 1971.576 с.

59. Нежиховский Р.А. Русловая сеть бассейна и процесс формирования стока воды. Д.: Гидрометеоиздат, 1971. 476 с.

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

61. Огиевский А.В. О соответственных уровнях и их использовании. М.: Геофизика, 1932. Т.2., вып. 2.

62. Парниковый эффект, изменение климата и экосистемы / Под ред. Б. Болина, Б.Р. Десса. Д.: Гидрометеоиздат, 1989. 557 с.

63. Полонский В.Ф., Лупачёв Ю.В., Скриптунов Н.А. Гидролого-морфометрические процессы в устьях рек и методы их расчёта. СПб.: Гидрометеоиздат, 1992. 383 с.

64. Поляков В.Н., Шлыгин И.А. Водный баланс и возможности изменения солености вод Таганрогского залива // Тр. ГОИН. 1979. Вып. 143. С. 2637.

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

66. Пространственно-временные колебания стока рек СССР / Под ред. Рождественского А.В. Л.: Гидрометеоиздат, 1988. 376 с.

67. Развитие антропогенной гидрологии в Средней Азии (опыт критического анализа) / Под ред. Ю.Н. Иванова. М.: Гидрометеоиздат, 1991. 55 с.

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

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

70. Раткович Д.Я. Об учёте нестационарности речного стока при оценке перспектив его использования // Метеорология и гидрология. 1992. №6. С. 101-105.

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

72. Регулирование и комплексное использование водных ресурсов. М.: Наука, 1987. 160 с.

73. Ремизова С.С. Водный баланс // Каспийское море. М., 1969. С. 107-138.

74. Родионов Н.А. Гидрология устьевой области Дона. М.: Гидрометеоиздат, 1958. 98 с.

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

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

77. Рубинова Ф.Э. Влияние водных мелиораций на сток и гидрологический режим рек бассейна Аральского моря // Тр. САРНИИ. 1987. Вып. 124 (205). 161 с.

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

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

80. Симов В.К. К вопросу о расчёте солёности вод Таганрогского залива // Тр. ГОИН. 1975. Вып. 125. С. 35-38.

81. Симов В.Г. Гидрология устьев рек Азовского моря. М.: Гидрометеоиздат (Московское отделение), 1989. 327 с.

82. Симонов А.И. Гидрология и гидрохимия устьевого взморья в морях без приливов // Тр. ГОИН. 1969. Вып. 92. 230 с.

83. Скриптунов Н.А. Течения на устьевом взморье Дона (Таганрогский залив) // Тр. ГОИН. 1978. Вып. 139. С. 43-47.

84. Соколовский Д.Л. О величине стока р. Волги в Каспийское море // Гидротехническое строительство. 1934. №7.

85. Статистические методы в гидрологии. Л.: Гидрометеоиздат. 1970. 123 с.

86. Тарасова М.Н., Бесчетнова Э.И. Гидрохимия Нижней Волги при регулировании стока // Гидрохимические материалы. 1987. Т. 101. 96 с.

87. Уланова Е.С., Забелин В.Н. Методы корреляционного и регрессионного анализа в агрометеорологии. Л.: Гидрометеоиздат, 1990. 150 с.

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

89. Христофоров А.В. Теория случайных процессов. М.: Изд-во МГУ, 1994. 139 с.

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

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

92. Шикломанов И.А., Кожевников В.П. Потери стока в Волго-Ахтубинской деятельности // Тр. ГТИ. 1974. Вып. 221. С.3-47.

93. Шикломанов И.А. Антропогенные изменения водности рек. Л.: Гидрометеоиздат, 1979. 302 с.

94. Шикломанов И.А. Влияние хозяйственной деятельности на речной сток. Л.: Гидрометеоиздат, 1989. 330 с.

95. Шикломанов И.А., Георгиевский В.Ю. Методы и результаты оценки влияния хозяйственной деятельности на водные ресурсы игидрологический режим рек // Тр. V Всесоюзного гидрологического съезда. Т.4. Л.: Гидрометеоиздат, 1990. С. 224-237.

96. Шикломанов И.А., Георгиевский В.Ю. Влияние антропогенных факторов на сток рек бывшего СССР // В кн.: Географическое направление в гидрологии. М.: РАН, 1995. С. 96-107.

97. Шикломанов И.А., Георгиевский В.Ю. Влияние антропогенных изменений климата на гидрологический и водные ресурсы // В кн.: Изменения климата и их последствия. СПб.: Наука, 2002. С.152-164.

98. Шульц В.А. Реки Средней Азии. Л.: Гидрометеоиздат, 1965. 690 с.

99. Schoenberg I. J. Contribution to the problem of approximation of equidistant data by analytic functions // Quart. Appl. Math., 1946. N.4. Pp. 45-99, 112141.