2.2 Прогнозирование индекса техногенной нагрузки в Республике Татарстан

Прогнозирование индекса техногенной нагрузки в сфере экологического состояния на 2013-2015 год на территории Республики Татарстан позволит оценить уровень негативного антропогенного воздействия на окружающую среду. Программа прогностического исследования уровня техногенной нагрузки на экологическом состоянии Республики Татарстан (Приложение А).

Анализ степени экологического состояния в целом проводится с использованием метода экстраполяции показателей индексов техногенной нагрузки.

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

Итак, с помощью Функции ПРЕДСКАЗ прогнозируется индекс техногенной нагрузки для Западного экономического района за 2013 и 2015 год. Для этого к таблице данных добавляется два столбца: «Период прогноза» и «Прогноз индекса техногенной нагрузки». Получаются следующие значения (Рисунок 1).

Рисунок 1- Прогнозное значение индекса техногенной нагрузки для Западного экономического района на 2014 и 2016 год

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

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

y=(1-α)×y_i +α×(y_i-1^сгл +T_i-1) (1)

где T_i=(1-β)×(y_i^сгл –y_i-1^сгл)+β×T_i-1;

α=0,4;

β=0,5;

y^сгл₁=y_1.

T₁=(y₂-y₁-y₄-y₃)/2 (2)

В результате применения метода Холта получаются сглаженные значения (Рисунок 2).

Рисунок 2- Сглаженные значения индекса техногенной нагрузки для Западного экономического района

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

1.Алгебраические многочлены

Функция регрессии в виде алгебраического многочлена имеет вид:

f(Х)=b₀+b₁X+b₂X²+…+b_mX^m, (3)

где b₀,b₁,b₂…,b_m - определяются на основе исходных данных.

Такая функция называется полиномиальной.

Для вычисления коэффициентов необходимо предварительно вычислить значения квадратов периодов. Далее коэффициенты многочленов с помощью функции ЛИНЕЙН.

Расчеты многочленов второй степени для Западного экономического района (Приложение Б, Рисунок Б.1).

Далее находится скорректированный коэффициент детерминации. Он вычисляется по следующей формуле:

1-(7-1)×(1-G4)/(7-2), (4)

где G4-коэффициент детерминации;

7-количество точек данных.

Рассчитанный скорректированный коэффициент детерминации многочлена второй степени равен 0,99561726.

Исходя из полученных данных строится график (по сглаживанию Холта и многочлену второй степени) (Рисунок 3).

Рисунок 3- График алгебраического многочлена для Западного экономического района

2.Логарифмическая функция

Логарифмическая функция применяется, если прогнозируемая переменная сначала имеет быстрый рост, который затем замедляется. Эта функция имеет вид:

Y=b₀+b₁ln(X) (5)

К исходным данным добавляется столбец, в котором вычисляются значение функции Ln(X). Коэффициенты b₀ и b₁ вычисляются с помощью функции ЛИНЕЙН.

Расчеты логарифмической функции для Западного экономического района (Приложение Б, Рисунок Б.2).

По аналогии с предыдущим примером вычисляется скорректированный коэффициент детерминации. Его значение равно 0,943973.

Исходя из полученных данных так же строится график(по сглаживанию Холта и логарифмической функции)(Рисунок 4):

Рисунок 4- График логарифмической функции для Западного экономического района

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

3.Гиперболическая функция

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

Y=b₀+b₁/X (6)

Добавляется столбец «Обратная функция», в котором записываются значения 1/X. Коэффициенты b₁ и b₀ вычисляются с помощью функции ЛИНЕЙН.

Расчеты гиперболической функции для Западного экономического района (Приложение Б, Рисунок Б.3).

По аналогии с предыдущим вычисляется скорректированный коэффициент детерминации. Его значение равно 0,748458.

Исходя из полученных данных строится график (по сглаживанию Холта и гиперболической функции) (Рисунок 5).

Рисунок 5- График гиперболической функции для Западного экономического района

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

4.Экспоненциальная функция

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

Y=b₀ (b₁)^X (7)

Коэффициенты b₀ и b₁ вычисляются с помощью функции ЛГРФПРИБЛ.

Расчеты экспоненциальной функции для Западного экономического района (Приложение Б, Рисунок Б.4).

По аналогии с предыдущим примером вычисляется скорректированный коэффициент детерминации. Его значение равно 0,975686.

Строится график (по сглаживанию Холта и экспоненциальной функции) (Рисунок 6).

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

Рисунок 6- График экспоненциальной функции для Западного экономического района

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

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

1.Выделение сезонной составляющей в аддитивной модели.

К исходным данным добавляются столбцы s_i и S_i, Корректирующий коэффициент m, сезонные коэффициенты, Функция прогнозирования (Рисунок 7).

Рисунок 7 - Расчеты сезонной составляющей в аддитивной модели для Западного экономического района

Значение si находится как индексы техногенной нагрузки минус многочлен второй степени; Si – находится по формуле:

Si=СРЗНАЧ(Е2;СМЕЩ(Е2;3;0);СМЕЩ(Е2;6;0);СМЕЩ(Е2;9;0)) (8)

Значение корректирующего коэффициента:

m=СУММ(F2:F4)/3 (9)

Сезонные коэффициенты:

F2-$G$2 (10)

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

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

Рисунок 8 - График сезонной составляющей в аддитивной модели для Западного экономического района

2.Выделение сезонной составляющей в мультипликативной модели.

Для мультипликативной модели проводятся похожие этапы- расчеты, как и для аддитивной.

Значение si находится отношением индексов техногенной нагрузки и многочлена второй степени; m – находится по формуле:

m=3/СУММ(Q2:Q4) (11)

Сезонный коэффициент:

Q2×$R$2 (12)

Функция прогнозирования находится как произведение сезонного коэффициента на многочлен второй степени (Рисунок 9).

Рисунок 9- Расчеты сезонной составляющей мультипликативной модели для Западного экономического района

Рисунок 10- График сезонной составляющей аддитивной модели для Западного экономического района

Следующий этап-построение функции прогнозирования и вычисление прогнозных значений.

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

1.Коэффициент детерминации.

Для его вычисления к исходным данным расчета добавляется столбец «Остатки», е², столбец «Квадрат отклонения», «Коэффициент детерминации» (Рисунок 11). В ячейку M3 записывается формула:

=1-СУММ(К2:К9)/СУММ(L2:L9) (13)

Рисунок 11 - Расчеты коэффициента детерминации аддитивной модели для Западного экономического района

2.Среднее абсолютное отклонение (показывает средний абсолютный уровень остатков). В ячейку М5 записывается формула:

=СРЗНАЧ(ABS(J2:J9)) (14)

3.Средняя абсолютная ошибка в процентах (показывает средний относительный уровень остатков). Ее величина не должна превышать 7,5%. В ячейку М7 записывается формула:

=СРЗНАЧ(ABS(J2:J9/C2:C9))×100 (15)

4.Тест на независимость остатков (может принимать значения от 0 до 4). В ячейку М9 для вычисления записывается формула:

=СУММКВРАЗН(J3:J9;J2:J8)/CУММКВ(J2:J9) (16)

Аналогичные расчеты проводятся и для мультипликативной модели (Рисунок 12).

Рисунок 12 - Расчеты коэффициента детерминации мультипликативной модели для Западного экономического района

Сравнивая все полученные показатели, выбираем ту модель, которая с большей вероятностью отражает исходные данные (модель, у которой Коэффициент детерминации больше, Среднее абсолютное отклонение меньше, Средняя абсолютная ошибка в процентах меньше 7,5% и Тест на независимость остатков дает результат либо меньше 0,5 либо 3,5).Указанным условиям удовлетворяет мультипликативная модель (Таблица 6).

Таблица 6 - Сравнение аддитивной и мультипликативной моделей

Аддитивная модель	Мультипликативная модель
Коэффициент детерминации	Коэффициент детерминации
0,997468836	0,997186405
Среднее абсолютное отклонение	Среднее абсолютное отклонение
0,000119967	0,000127571
Средняя абсолютная ошибка в процентах	Средняя абсолютная ошибка в процентах
0,741462132	0,775209095
Тест на независимость остатков	Тест на независимость остатков
1,02956029	1,082012895

Таким образом, исходя из расчетов сезонной составляющей в мультипликативной модели для Западного экономического района, индексы техногенной нагрузки для 2013-2015 года будут следующими: 2013-0,017591781, 2014-0,018611942, 2015-0,1996571(Рисунок 9).

Аналогично таким же образом прогнозируются индексы техногенной нагрузки на 2013-2015 год для Южного, Северного, Восточного, Нефтяного, Столичного и Камского экономического района.

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

Экономический район	2013	2014	2015
Западный	0,017	0,018	0,019
Южный	0,038	0,039	0,040
Северный	0,050	0,052	0,054
Восточный	0,069	0,072	0,074
Нефтяной	0,168	0,174	0,180
Столичный	0,273	0,266	0,259
Камский	0,500	0,493	0,485

Таблица 7 – Спрогнозированный индекс техногенной нагрузки экономических районов Республики Татарстан на 2013-2015 год

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

Индекс техногенной нагрузки в Нефтяном экономическом районе к 2015 году по сравнению с предыдущими годами увеличится, а Столичный и Камский экономические районы к 2015 уменьшат свой показатель, но, все же, он останется самым высоким показателем техногенной нагрузки по Республике Татарстан. В остальных экономических районах индекс техногенной нагрузки имеет наименьшее значение и пока не представляет угрозы. Динамика показателя индекса в этих районах не сильно растет и остается на допустимом уровне. А для других, что касается отмеченных выше районов с наибольшим индексом, необходимо предусмотреть меры, разработать определенные мероприятия для снижения уровня техногенной нагрузки.