Рациональное управление водными ресурсами является одной из первостепенных задач по жизнеобеспечению населения. В связи со значительным сокращением густоты речной сети за последние 200 лет [1, 7 и др.] безотлагательно необходимы научно обоснованные оценки происходящих и будущих изменений водного режима рек и водных ресурсов, а также выявление закономерностей этих изменений в современных социально-экономических и климатических условиях.
Малые реки особенно чутко реагируют на антропогенные воздействия и служат интегральным индикатором сложных природно-антропогенных процессов, происходящих на их водосборах. В современных исследованиях [3, 4, 6 и др.] бассейновые территориальные структуры рассматриваются как иерархические общности пространственных отношений, определяемых стоком воды, наносов и растворенных веществ. Интегратором эколого-хозяйственных условий на территории речных бассейнов служит величина и характеристики стока. Однако длинные ряды наблюдений на малых реках крайне редки.
На берегах малых рек проживает значительная часть населения России. Городские малые реки часто становятся объектом хищнического потребления ресурсов речной системы городскими структурами, зачастую не сильно заботящихся о сохранении и восстановлении водной экосистемы. Поэтому актуальным вопросом при ведении хозяйственной деятельности является анализ и прогноз водного режима городских малых рек.
Цель работы – используя результаты моделирования, получить прогноз изменения стока малой реки в условиях антропогенной трансформации и изменяющегося климата.
Материалы и методы
Объектом исследований выбран правый приток Северского Донца – Болховец (Везёлка). Это типичная малая река V-го порядка длиной 27 км, которая является основной рекой города Белгорода: 55 % площади городских земель приходится на ее бассейн (394 км2).
Наиболее интегральное представление о миграции и аккумуляции потоков вещества в пределах бассейна позволяет получить организованный в замыкающем створе водотоков гидроэкологический мониторинг [2, 3]. Исходными данными для анализа речного стока за 1947–2010 гг. послужили наблюдения по гидропосту на р. Везелка в черте Белгорода. Для определения влияния климатических факторов на динамику расходов воды использованы данные многолетних наблюдений метеорологических станций Харькова и Белгорода.
Учитывая высокую динамичность и сложность многих процессов гидрологического цикла, адекватность их моделирования могут обеспечить не менее сложные методы [8]. В работе применяли методы статистики (анализ временных рядов, нейронные сети), картографии и ГИС-технологии. Пространственное моделирование проводили в ArcGIS 9.3; анализ временных рядов – в STATISTIСA 6.0, модуль Neural Network. Решая проблему выбора метода для создания универсальной нелинейной модели в целях прогнозирования динамики водного режима малых рек в зависимости от климатических факторов, мы использовали возможности искусственных нейронных сетей. Принцип работы нейромодели основан на том, что нейроны первого слоя получают входные сигналы, преобразуют их и через точки ветвления передают нейронам второго слоя. Далее срабатывает второй слой и т.д. до j-го, который выдает выходные сигналы аппроксимации.
Результаты
Известно, что там, где в реки поступают наносы с урбанизированных частей водосборов, донные отложения загрязнены тяжелыми металлами и другими поллютантами. Но ближе к устью содержание загрязнителей в донных отложениях снижается. Однако в случае с Везелкой, у которой нижняя треть речной долины приходится на территорию Белгорода, эта закономерность не наблюдается.
Расход речного стока 50 %-ной обеспеченности оценивается для этой реки в 1,25 м3/с, для среднего маловодного года он составляет 0,91, а для маловодного года – 0,48 м3/с [5]. Анализ формирования водного режима можно разделить на два этапа (рис. 1): природное колебание стока реки (1947–1980 гг.) и стабильное снижение расхода воды (1981–2010 гг.). Максимальный расход воды в реке Везелка за весь период составил 2,18 м3/с (1981 г.).
Рис. 1. Изменение расхода воды в реке Везелка за период 1947–2010 гг.
Резкое снижение расходов во втором периоде объясняется сильным антропогенным преобразованием русла Везелки. В 1980–1982 гг. и 1999–2005 гг. проводили углубление и расширение русла (с 3–5 до 37–45 м) на протяжении 2,5 км, используя земснаряд на локальных участках. По сравнению с периодом 1947–1980 гг., когда в среднем расход воды составлял 1,19 м3/с, он уменьшился на 0,40 м3/с (т.е. на 34 %).
Климат является одним из основных факторов, определяющих общую величину и внутригодовое распределение стока. Для определения степени влияния климатических факторов на динамику стока, а также для обеспечения ввода климатических данных в разрабатываемую нейромодель, проведены статистические преобразования временных рядов. Вначале с помощью одномерного спектрального анализа Фурье были определены сезонные составляющие. Для динамики температуры воздуха цикл составил – 10 лет, а для динамики суммы годовых осадков и речного стока – 11 лет. Далее методом Censusu I (сезонная декомпозиция) проведена корректировка ряда (удаление сезонной составляющей) и пятиточечное центрирование для получения сглаженных тренд-циклических рядов. Это позволяет определить закономерности развития временных рядов в «чистом» виде (рис. 2).
Рис. 2. Статистически преобразованные среднегодовые значения изучаемых показателей
При комплексном анализе влияния климатических факторов на динамику среднегодового стока реки (рис. 3) наблюдаются следующие сценарии формирования стока:
– при уменьшении среднегодовой температуры воздуха и минимальных, но стабильных осадках сток реки возрастает;
– при увеличении среднегодовой температуры воздуха и увеличение суммы годовых осадков сток реки стремительно уменьшается;
– при увеличении суммы осадков и уменьшения температуры сток реки увеличивается.
Для оценки влияния климата на динамику речного стока была использована кусочно-регрессионная модель. Увеличение уровня аппроксимации такой модели и определения условий формирования стока реки Везелка в зависимости от климатических факторов было достигнуто путем разделения временного ряда на два явно выраженных периода (до и после 1981 г.). Комплексная уточненная кусочно-регрессионная модель имеет вид:
Рис. 3. Зависимость стока реки Везелка от климатических факторов
Рис. 4. Моделирование влияния климатических факторов на расходы воды реки Везелка по периодам их формирования
На первом этапе коэффициент корреляции влияния факторов на динамику стока реки составил – 0,86. Объясненная доля дисперсии влияния факторов равна 0,74. Коэффициент корреляции влияния факторов на динамику стока реки составил – 0,90.
Для прогнозирования динамики стока реки в зависимости от климатических факторов создана искусственная многослойная нейромодель вида (табл. 1):
где, y – сток реки, м3/с; х – сумма годовых осадков; z – среднегодовая температура воздуха.
Таблица 1. Характеристика нейронной модели для прогнозирования динамики стока малой реки в зависимости от климатических факторов
Архитектура |
Производительность |
Ошибка |
Обучение/эпохи |
Количество нейронов в скрытом слое |
|||||||||
обучения |
контрольная |
тестовая |
обучения |
контрольная |
тестовая |
1 слой |
2 слой |
||||||
МП s9 2:18-9-7-1:1 |
0,223 |
0,450 |
0,561 |
0,051 |
0,102 |
0,129 |
ОР100, |
9 |
7 |
||||
Функция коррекции весовых коэффициентов |
|
||||||||||||
Функция отклика сети |
|
||||||||||||
Алгоритм обучения сети |
|
||||||||||||
Архитектура нейросети |
Свойства нейромодели |
||||||||||||
|
Четырехслойный перцептрон с девятью нейронами в первом скрытом слое и семью нейронами во втором скрытом слое. Метод обучения: обратное распределение (100 эпох) и связанных градиентов (25 эпох). Матрица искусственной нейронной сети состоит из 232 весовых коэффициентов. Функция активации нейронов – синусоид-гиперболического тангенса |
||||||||||||
Ошибка обучения и тестирование сети |
Весовая гистограмма |
||||||||||||
|
|
||||||||||||
Коэффициенты масштабирования входных сигналов |
|||||||||||||
Показатель |
Сдвиг |
Масштаб |
|||||||||||
Среднегодовая температура воздуха |
-1,350 |
0,247 |
|||||||||||
Годовая сумма осадков |
-0,875 |
0,002 |
|||||||||||
Расход воды в год |
-0,055 |
0,457 |
|||||||||||
Смоделированная нейронная сеть обладает достаточно высокими обобщающими способностями, на что указывает хорошее совпадение исходных и предсказанных значений (верно выделенных закономерностей из входных данных) на рабочих (контрольном и тестируемом) участках. Такое качество многослойной нейронной сети обеспечило возможность корректно отображать данные, которые не использовались в обучении, и достаточно достоверно прогнозировать изменение речного стока от климатических факторов. Подтверждением результатов эффективного моделирования нейросетей для прогнозирования является анализ итоговой статистики реализации нейронных сетей на обучающей, контрольной и тестовой выборках (табл. 2). Достоверность ИНС составила – 84 %.
Таблица 2. Итоговая статистика обучения нейронной сети для прогнозирования динамики стока малой реки Везелка в зависимости от климатических факторов
Описательные статистики |
Обучающая выборка |
Контрольная выборка |
Тестовая выборка |
Математическое ожидание ошибки |
0,0259 |
-0,0305 |
0,0122 |
Стандартное отклонение ошибки |
0,1076 |
0,2214 |
0,2811 |
Математическое ожидание абсолютной ошибки |
0,0797 |
0,1535 |
0,2121 |
Корреляция |
0,976 |
0,893 |
0,828 |
Рис. 5. Динамика и прогноз стока реки Везелка
По прогнозу (рис. 5) к 2022 г. ожидаются стабильно минимальные значения стока Везелки в пределах до 0,44 м3/с. Среднее значение прогнозируемого речного стока гораздо меньше его средних значений за рассматриваемые периоды, что обусловлено снижением вклада климатических факторов при коренном антропогенном изменении параметров русла.
Заключение
Анализ влияния климатических факторов на водность малой реки показал высокую зависимость погодичных величин расхода воды от совокупного влияния температуры воздуха и атмосферных осадков, т.е. главенствующим фактором при формировании речного стока малых рек являются климатические условия. За последние тридцать лет несколько усилился отклик динамики стока малой реки от климатических факторов (с 0,86 до 0,90). Выявлены два хорошо выраженных периода динамики стока исследуемой реки: до и после преобразования русла. Во втором периоде устойчивость водной системы была разрушена, т.е. река в городских условиях при изменении параметров русла уже не в состоянии самостоятельно противостоять антропогенному давлению и требует постоянного контроля и охраны.
Прогнозная нейросетевая модель способна не только моделировать динамику водности рек, но и учитывать разнородную информацию о текущих и ожидаемых режимах функционирования водной системы. В отличие от «традиционных» методов статистики нейронные сети находят не статистически достоверное, а правдоподобное решение задачи и могут работать при определенном дефиците эмпирических данных. Эти преимущества делают нейронные сети все более и более перспективным направлением в анализе и прогнозировании динамики стока малых рек в зависимости от факторов его формирования.
Исследование выполнено при финансовой поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» в рамках научного проекта № 14.132.21.1387.
Рецензенты:
Корнилов А. Г., д.г.н., профессор, заведующий кафедрой географии и геоэкологии ФГАУ ВПО Белгородский государственный национальный исследовательский университет, г. Белгород.
Кумани М. В., д.г.н., профессор, заведующий кафедрой географии и геоэкологии ФГБОУ ВПО Курский государственный университет, г. Курск.
Библиографическая ссылка
Кузьменко Я.В., Лисецкий Ф.Н., Пичура В.И. ОЦЕНКА И ПРОГНОЗИРОВАНИЕ СТОКА МАЛЫХ РЕК В УСЛОВИЯХ АНТРОПОГЕННЫХ ВОЗДЕЙСТВИЙ И ИЗМЕНЕНИЙ КЛИМАТА // Современные проблемы науки и образования. – 2012. – № 6. ;URL: https://science-education.ru/ru/article/view?id=7640 (дата обращения: 23.02.2025).