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

Размещено на http://www.allbest.ru/

УДК 519.2:330.007

ФГОУ ВПО МГУП

АНАЛИЗ МНОГОЛЕТНИХ КОЛЕБАНИЙ ГОДОВОГО СТОКА р. ДОН

Н.В. Муращенкова

О.С. Фалеева

г. Москва

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

Совместный анализ разностно-интегральных кривых наблюденного и условно-естественного стока позволил выявить относительную роль климатического и антропогенного фактора в многолетних изменениях стока [4].

Естественный (с 1936 года условно-естественный) сток воды р. Дон с 1981-1900 гг. был близок к среднему многолетнему, с 1901-1914 гг. наблюдалась фаза пониженного, а с 1915-1933 гг. - повышенного. С 1934-1940 гг. наблюдалась фаза резкого понижения стока, сменившаяся резким повышением в 1941-1943 гг., затем периоды средней водности чередовались с периодами понижения водности в 1948-1955 и 1971-1977 гг. После резкого повышения стока с 1978-1981 гг. сток близок к среднему многолетнему значению.

Средняя многолетняя величина условно-естественного стока равна 27,3 км3 . Причем, условно-естественный сток воды в наиболее выраженный маловодный период 1933-1939 гг. составил 17,5 км3, а в наиболее выраженный многоводный период 1940-1942 гг. - 45,1 км3 .

На фоне выделенных выше периодов чередования маловодных и многоводных лет, обусловленных климатическими причинами, для наблюденного стока весь период с 1949- 2000 гг. в целом маловодный, в то время как восстановленный сток в этот период незначительно колебался около естественной нормы, за исключением маловодного периода 1972-1976 гг. и многоводного периода 1977-1981 гг.

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

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

До 1951 года среднегодовые потери стока не превышали 1,8 км3, в период заполнения Цимлянского водохранилища 1952-1955 гг. они составили в среднем около 6,3 км3 .

В период 1956-1965 гг. потери составили в среднем 3,2 км3 . В период 1966-1980 гг. - 4,7 км3 , в 1981-1990 гг. потери еще увеличились и составили в среднем 6,7 км3 , при этом максимальные годовые потери 8,6 км3 были в 1990 году. Затем потери снизились и составили в 1991-1999 году 3,9 км3, минимальные годовые потери 1,7 км3 были в 1996 году.

Средний условно-естественный сток в периоды 1881-1951 гг. и 1952-2000 гг. различается незначительно и составляет, соответственно 27,8 и 26,5 км3. Наблюденный сток в 1952-2000 году составил 21,6 км3 , то есть на 21% ниже нормы естественного стока.

В последние 19 лет условно-естественный сток был близок к естественной норме и составил 26,7 км3 , а наблюденный сток составил 21,1 км3 - на 21% меньше. В многоводный период 1977-1981 гг. средний условно-естественный сток составил 32,9 км3, а наблюденный сток составил 28,5 км3 - на 13% меньше.

В таблице приведены основные статистические характеристики рассматриваемых временных рядов стока по периодам.

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

Основные статистические параметры годового стока р. Дон

Примечание: (числитель - наблюденный сток, знаменатель - условно-естественный, Wср - среднемноголетний объем, W - среднеквадратичное отклонение, СV - коэффициент вариации, R(1) - коэффициент корреляции стока смежных лет).

Модель авторегрессии (АР) можно представить в виде

Хt = Ф1Хt-1 + Ф2Хt-2 + …+ ФpХt-p + at,

где Хt - совокупность случайных величин в определенные моменты времени t, предсказанная по совокупности значений Хt-p в предыдущие моменты времени (t-p); Фp - параметр модели (предиктор величины Хt по величинам Хt-p); at - случайная величина, характеризующая белый шум; p - порядок модели авторегрессии АР(p).

Модель скользящего среднего (СС) представляет собой

Хt = at - 1аt-1 + 2аt-2 + …+ qаt-q,

где q - параметр модели (предиктор величины Хt по величинам аt-q); q - порядок модели скользящего среднего СС(q).

Смешанная модель авторегрессии - скользящего среднего (АРСС(p,q))

Хt = Ф1Хt-1 + Ф2Хt-2 + …+ ФpХt-p + at - 1аt-1 + 2аt-2 + …+ qаt-q,

Однако рассмотренные выше модели применимы только к стационарному ряду, а исследуемые гидрологические ряды являются нестационарными.

Для приведения нестационарного ряда к стационарному вводится оператор разности = 1 - В. Тогда Хt = Хt - Хt-1, и модель (3) можно представить в виде: климатический антропогенный сток водность

Ф(В) d Хt = (В) at,

где Ф(В), (В) - стационарные операторы авторегрессии и скользящего среднего порядка p и q, соответственно; d - порядок разности модели.

Разность порядка d ряда {Хt} вычисляется по формуле

Wt d Хt = (1 - B)d Хt = ( Хt = ,

и является стационарным обратимым процессом АРСС(p,q).

Так как процесс {Хt} можно получить d-кратным суммированием (интегрированием) процесса {Wt}, то процсс, задаваемой моделью (4) называют процессом АРПСС(p,d,q), где П означает проинтегрированной. Процесс АРПСС(p,d,q) может использоваться для моделирования как для стационарных, так и для нестационарных временных рядов.

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

Пусть исходный динамический ряд имеет вид

.

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

Выражение

представляет собой экспоненциальную среднюю k-го порядка для ряда yt, где - параметр сглаживания. В расчетах для определения экспоненциальной средней используется рекуррентная формула [5]

.

Использование соотношения (8) предполагает задание начальных условий S0[1], S0[2], …, S0[k]. Для этого можно воспользоваться формулой Брауна-Мейера, связывающей коэффициенты прогнозирующего полинома с экспоненциальным средним соответствующих порядков

,

где р = 1,2,…, n+1; - оценка коэффициентов; = 1-. Для определения вида порядка модели необходимо рассмотреть выборочные автокорреляционную функцию (АКФ) и частную автокорреляционную функцию (ЧАКФ). При этом используются следующие практические критерии [3]:

если АКФ экспоненциально затухает, а ЧАКФ имеет выброс лишь при сдвиге по времени, равном 1, то процесс наилучшим образом описывается моделью авторегрессии первого порядка (АР (1));

если АКФ имеет форму затухающей синусоидальной волны или экспоненциально затухает, а ЧАКФ имеет выброс только при сдвиге по времени, равном 1 и 2, то процесс описывается моделью авторегрессии второго порядка (АР(2));

если АКФ имеет выброс при сдвиге по времени, равном 1, а ЧАКФ экспоненциально затухает, то подходящей является модель скользящего среднего первого порядка (СС (1));

если АКФ имеет выбросы при сдвиге по времени, равном 1 и 2, а ЧАКФ имеет форму синусоидальной волны или экспоненциально затухает, то используется модель скользящего среднего второго порядка (СС(2));

если АКФ и ЧАКФ экспоненциально затухают, начиная со сдвига по времени, равного 1, то используется модель авторегрессии и скользящего среднего первого порядка (АРСС (1,1)).

Так как для АКФ и ЧАКФ временного ряда годового стока р. Дон не имеют выбросов и четко выраженных тенденций к затуханию (рис. 2), то, следовательно, ряд является нестационарным, и модели авторегрессии и скользящего среднего не подходят для его описания, так как они применимы только к стационарным рядам [1]. Для приведения ряда к стационарному необходимо трансформировать его, взяв первую разность. АКФ трансформированного временного ряда имеет выброс при сдвиге по времени, равном 1, остальные значения статистически незначимы, а ЧАКФ экспоненциально затухает (рис. 3). Следовательно, рассматриваемый процесс условно-естественного стока р. Дон может быть идентифицирован моделью АРПСС (0,1,1) с параметром 1 = 0,62851. Анализ остатков между расчетными и наблюденными величинами свидетельствует в пользу правомерности использования данной модели, так как остатки похожи на белый шум. Для более детального описания процесса стока применялись модели с сезонной составляющей (мультипликативные) с различными «сезонными» периодами. Наилучшие результаты дает модель с лагом 6, с параметрами Фp = 0,14605, 1 = 0,69984, Фp = -0,8804, 1= - 0,7796. Прогнозируемый годовой сток р. Дон по данной модели принимает значения от 26,5 до 28,6 км3/год.

При использовании модели экспоненциального сглаживания были получены следующие результаты. Наиболее адекватной оказалась аддитивная модель с сезонной составляющей и с затухающим трендом с различными параметрами. Наиболее адекватной представляется модель с параметрами: = 0,5, = 0,3, = 0,1, сезонный лаг - 11. Прогнозируемые значения годового стока р. Дон при применении данной модели колеблются от 23,5 до 31,0 км3/год.

Рис. 1. Интегральная кривая функция годового стока р. Дон, стока р. Дон Рис. 2. Автокорреляционная ряда условно-естественного стока р. Дон за период 1881-2000

Рис. 3. АКФ(а) и ЧАКФ (б) преобразованного по первой разности ряда многолетних колебаний стока р. Дон за период 1881-2000 гг.

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

Библиографический список

1. Боровиков В. П. , Ивченко Г. И. Прогнозироание в системе STATISICA в среде WINDOWS. М.: Финансы и статистика. 1999. 382 с.

2. Исмайылов Г. Х., Перминов А. В. Возможные подходы к оценке будущих значений гидрологических характеристик //В сб. Проблемы научного обеспечения развития эколого-экономического потенциала России. М.: МГУП. 2004. С. 44-49.

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

4. Отчет ГОИН Комплексное исследование гидролого-экологического состояния устьевых областей рек южных морей России. М., 2002.

5. Рабочая книга по прогнозированию. М.: Мысль, 1982. 430 с.

Размещено на Allbest.ru