Обширные площади территории Якутии находятся в различных климатических зонах. В республике хорошо развитая речная сеть, которая играет огромную роль в жизнедеятельности и жизнеобеспечении региона как одна из основных транспортных магистралей при перевозке пассажиров и грузов в короткий навигационный период, длительность которого зависит от географической и природно-климатической специфики. Ежегодные колебания уровней и расходов воды в реке представляют собой случайные величины, поэтому для построения прогностических моделей часто используются методы математической статистики. Статистические модели позволяют учесть некоторые общие закономерности межгодовых колебаний климата, многолетних колебаний среднегодовых расходов воды, его цикличности, синхронности и асинхронности в различных районах и построить прогнозные модели на ближайшие годы. Исследование синхронности колебаний среднегодовых расходов воды в реках имеет важное значение также при решении многих научных и практических гидрологических задач, связанных с комплексным использованием водных ресурсов при планировании и проектировании крупных энергетических и водохозяйственных систем. Синхронность колебаний может быть полезна для установления закономерности изменчивости годового стока вод в арктической зоне, из-за недостаточного количества гидрологических постов. Представленные результаты могут стать основой для экстраполяции выявленных закономерностей расходов воды северных рек на более значительные временные интервалы, а возможно и на другие регионы со схожими природно-климатическими условиями.
Целью данной работы является решение задачи, связанной с анализом, адаптацией существующих и разработкой новых подходов при исследованиях закономерностей многолетних колебаний среднегодовых расходов воды в реках, полученных по многолетним данным наблюдений и статистики Гидрометцентра, используемых для планирования и проектирования крупных энергетических и водохозяйственных систем на территории Республики Саха (Якутия).
В создание методов и моделей, основанных на вероятностной теории описания речного стока, большой вклад был сделан российскими исследователями (Н.А. Картвелишвили, М.Ф. Менкель, В.А. Лобанов, С.Н. Крицкий, Д.А. Бураков и др.). Результаты исследований многолетних изменений основных параметров водного баланса рек и водосборов учитывались и использовались при построении моделей в работах [1–3]. В исследованиях [4–6] рассматривались, анализировались и классифицировались экстремальные гидрологические явления на территории России и других стран, был определен вклад природных и антропогенных влияний в их формирование. На основании полученных прогнозных моделей были построены сценарии развития опасных гидрологических ситуаций, что давало возможность устранения и смягчения негативных последствий, а также использования полученных результатов при проектировании и управлении хозяйственными и промышленными объектами.
Материалы и методы исследования
В данной работе на основе метода парных корреляций и исследования трендов среднегодовых расходов воды выполнен анализ синхронности и пространственной однородности колебаний расходов воды рек Яны, Колымы и Вилюя. Выбор рек производился с учетом следующих критериев: наличие достаточно продолжительного ряда гидрологических наблюдений (более 75 лет); «непересечение» водосборов рек; отсутствие в бассейнах исследуемых водотоков существенных антропогенных воздействий на сток. Исследуемая территория состоит из двух регионов: регион I – реки бассейна Яны; регион II охватывает реки бассейна Колымы и Вилюя.
Река Колыма образуется от слияния рек Аян-Юрях и Кулу, берущих начало на Охотско-Колымском нагорье. Длина 2129 км (от истока реки Кеньеличи, правой составляющей реки Кулу, – 2513 км), из них около 1,4 тыс. км на территории Магаданской области, остальное – на территории Якутии. Площадь бассейна 643 тыс. км2. Впадает в Колымский залив Восточно-Сибирского моря тремя главными протоками: Каменная Колыма – правая, судоходная, Походская Колыма и Чукочья. Длина дельты 110 км, площадь 3000 км2.
Река Яна образуется слиянием рек Сартанг и Дулгалах, которые стекают со склонов Верхоянского хребта. Длина – 872 км, площадь бассейна – 238 000 км2. При впадении в Янский залив моря Лаптевых река образует дельту площадью 10200 км2. При разработке методов и моделей прогнозирования поведения временных рядов за некоторый период времени выполняются следующие шаги [5, 6]:
– Необходимо провести предварительный анализ данных;
– Разработка новых или адаптация известных моделей, т.е. выбор кривых, описывающих явление, и численное оценивание параметров модели;
– Проверка адекватности и оценка точности используемых моделей;
– Оценка точечного и интервального прогноза.
Одними из методов вероятностно-статистического аппарата, которые подходят для прогнозирования изменений среднегодовых расходов воды в реках является метод гармонического анализа временных рядов [7, 8]. Построим прогностические модели среднегодовых расходов воды на реках Яна, Колыма, Вилюй до 2011 г. и на их основе составим числовой прогноз изменения до 2015 г. В данной работе необходимо выполнить прогнозирование временного ряда Ht+1, Ht+2,..., по наблюдениям H1, H2,..., Ht где t – время. Предполагается, что присущие данному процессу особенности и закономерности сохраняются в течение исследуемого периода, что является корректным при моделировании краткосрочных либо оперативных прогнозов. Для временных рядов используется часто статистический подход, который состоит в том, что временной ряд, описывающий развитие процесса, можно представить в виде трех составляющих компонент: функцию тренда (тенденция развития), циклическую компоненту, остаточную компоненту (белый шум).
H(t) = f(t) + u(t) + e(t),
где f(t) – функция тренда (тенденция развития), u(t) – циклическая компонента, e(t) – остаточная компонента.
Предварительный анализ данных состоит в выполнении сглаживания данных наблюдения и определении наличия тренда временного ряда. При определении вида функции тренда и вычислении коэффициентов применяется метод наименьших квадратов. Этап удаления тренда необходимо сопровождать исследованием временного ряда на стационарность. Оценка зависимостей между последовательными значениями одного и того же ряда производится используя коэффициент автокорреляции. Считая значения рассматриваемого временного ряда случайными, предполагаем значения между их уровнями независимыми. Величина коэффициента автокорреляции может быть использована при определении периода колебаний. Построение статистической модели было выполнено, используя методику обработки данных, в основе которой применялась нелинейная многопараметрическая регрессия [7].
Начальным этапом моделирования данных является построение регрессионной модели. Для получения прогнозных значений среднегодовых расходов воды необходимо провести композицию составляющих модели и обратное преобразование. Адекватность полученной модели определяется сравнением данных наблюдения с результатами прогноза.
Прогноз среднегодовых расходов воды на реке Яна в районе Верхоянска (данные с 1935 г.). Судоходство на участке Батагай – Верхоянск нерегулярно и возможно лишь в полноводный период. Для получения достаточно точных и надежных прогнозов необходимо изучить состояние процесса и выделить основные компоненты, которые входят в состав исследуемого временного ряда (тренд и циклические колебания – отклонения уровней отдельных периодов времени от тенденции динамики). Для аппроксимации временного тренда наиболее простым и чаще всего используемым подходом является его полиномиальное представление, когда соответствующая составляющая временного ряда – функция является полиномом, коэффициенты которого могут быть найдены методом наименьших квадратов.
Результаты исследования и их обсуждение
Стандартным методом выделения циклических составляющих временного ряда является разложение его в ряд Фурье. Используя методы предварительной обработки данных и анализ Фурье, мы получили математическую гармоническую модель, описывающую исходные данные, на основе которой делаем прогнозирование.
На рис. 1 представлены среднегодовые расходы воды, исходные данные и гармоническая составляющая на реке Яна в районе Верхоянска.
Прогноз среднегодовых расходов воды на реке Колыма в районе Среднеколымска (данные с 1934 г.). Колыма судоходна от устья реки Бахапча (регулярное судоходство – от Сеймчана); навигация 4–5 мес. Аналогичным образом производится выявление наличия тренда, сглаживание временных рядов для рек Колыма в районе Среднеколымска и Вилюй в районе с. Сунтар. Подборка коэффициентов и выбор моделей тренда осуществляется на основании метода наименьших квадратов. После удаления тренда необходимо провести анализ временного ряда на стационарность. К ряду остатков применяем гармонический анализ. Суммарный прогноз строится с учетом тренда и экстраполированных значений ряда Фурье.
а)
б)
Рис. 1. а) исходные данные (синим) и тренд (красным) среднегодовых расходов воды на реке Яна в районе Верхоянска; б) исходные данные (синим) и гармоническая составляющая (красным) на реке Колыма в районе Среднеколымска
Таблица 1
Таблица относительной погрешности прогностических моделей
Участки рек |
2012 |
Ош. |
2013 |
Ош. |
2014 |
Ош. |
2015 |
Ош. |
||||
Факт |
Прог. |
% |
Факт. |
Прог. |
% |
Факт. |
Прог. |
Факт |
Прог. |
% |
||
Яна в районе Верхоянска |
209 |
204 |
2 |
166 |
190 |
14 |
125 |
161 |
28 |
188 |
171 |
9 |
Колыма в районе Среднеколымска |
2460 |
2671 |
8 |
3410 |
2748 |
19 |
3450 |
2839 |
17 |
2040 |
2639 |
29 |
Вилюй в районе Сунтар |
629 |
684 |
8 |
565 |
581 |
3 |
859 |
659 |
23 |
1070 |
816 |
24 |
Рис. 2. Тренды временных рядов среднегодовых расходов воды рек Яна, Вилюй и Колыма в районе Верхоянска, Среднеколымска и Сунтар соответственно
Среднегодовые расходы воды на реке Вилюй в районе Сунтар (данные с 1932 г). Река Вилюй – это самый длинный приток реки Лены. Регулярное судоходство осуществляется до посёлка Сунтар. Циклическая составляющая ряда остатков строится на основании результатов гармонического анализа. Суммарный прогноз строится с учетом тренда и экстраполированных значений ряда Фурье. Результаты расчетов и оценка погрешностей представлены в табл. 1.
Адекватность построенных моделей прогнозирования оценивается путем сравнения фактических и предсказанных значений. Полученные значения относительной погрешности говорят о среднем уровне точности построенных моделей.
Степень синхронности или асинхронности колебаний среднегодовых расходов воды рек была определена с помощью метода парной корреляции, результаты представлены в табл. 2.
Таблица 2
Коэффициенты корреляции среднегодовых расходов воды рек Яна, Колыма и Вилюй в районе Верхоянска, Среднеколымска и Сунтар соответственно
№ п/п |
Участки рек |
1 |
2 |
3 |
1 |
Яна в районе Верхоянска |
1 |
0,34 |
0,17 |
2 |
Колыма в районе Среднеколымска |
0,34 |
1 |
0,04 |
3 |
Вилюй в районе Сунтар |
0,17 |
0,04 |
1 |
Анализ коэффициентов корреляции среднегодовых расходов воды рек Яна, Колыма и Вилюй в районе Верхоянска, Среднеколымска и Сунтар соответственно показывает, что корреляция между ними очень слабая, синхронность данных рек практически отсутствует. При установлении синхронности, асинхронности и цикличности колебаний расходов воды рек используются различные показатели и методы. В данной работе кроме метода парной корреляции использовался также графический – сравнение временных трендов среднегодовых расходов воды, который показал асинхронность рек Вилюй и Колыма, коэффициент парной корреляции этих рек равен 0,04. На графике можно отметить, что наиболее маловодные годы реки Колыма совпадают многоводному периоду на реке Вилюй и наоборот (рис. 2).
Заключение
Одним из методов вероятностно-статистического аппарата, который подходит для прогнозирования изменений многолетних колебаний среднегодовых расходов воды, является метод гармонического анализа временных рядов. Предложенный подход позволяет получать прогнозные оценки многолетних колебаний среднегодовых расходов воды на определенный промежуток времени и получать модели изучаемого процесса (тренда, гармонической и случайной составляющих). Степень синхронности или асинхронности колебаний среднегодовых расходов воды рек определялась с помощью метода парной корреляции. Установлено, что в многолетних колебаниях среднегодовых расходов воды исследованных рек синхронность выражена слабо. Но в некоторых случаях коэффициент парной корреляции может оказаться недостаточной характеристикой синхронности колебаний водности, так как коэффициент корреляции определяется на основании всего ряда среднегодовых расходов воды и является усредненной характеристикой. Поэтому в данной работе кроме метода парной корреляции использовался также графический – сравнение временных трендов среднегодовых расходов воды, который показал асинхронность рек Вилюй и Колыма. Предложенный метод может быть полезен для прогнозирования закономерностей колебания уровня воды на реках Севера.
Библиографическая ссылка
Капитонова Т.А., Слепцов О.И., Кусатов К.И., Стручкова Г.П. АНАЛИЗ МНОГОЛЕТНИХ КОЛЕБАНИЙ СРЕДНЕГОДОВЫХ РАСХОДОВ ВОДЫ РЕК ЯКУТИИ // Успехи современного естествознания. – 2018. – № 11-1. – С. 109-114;URL: https://natural-sciences.ru/ru/article/view?id=36913 (дата обращения: 16.08.2024).