Математические методы анализа и прогнозирования в электроэнергетике. Учебное пособие
.pdfВ многомерном ряду представлена динамика нескольких показателей, характеризующих один процесс (явление). Например, показатели деятельности предприятия, наблюдаемые помесячно, можно представить в виде многомерного ряда (табл. 1.2). В качестве показателей можно взять показатели регионального электропотребления:
Y – объемы потребления электроэнергии в регионе; Х1 – объем валового внутреннего продукта;
Х2 – объем покупки электроэнергии по свободным и регулируемым долгосрочным договорам;
Х3 – объемы покупки электроэнергии на рынке на сутки вперед;
Х4 – объемы покупки электроэнергии на балансирующем рынке;
Х5 – температура окружающей среды.
Таблица 1.2
Многомерный временной ряд
Момент времени |
Номер момента |
|
Значения показателей |
|
||
(дата) |
времени, t |
Y |
X1 |
X2 |
… |
X5 |
Январь 2020 г. |
1 |
y1 |
x11 |
x12 |
… |
x15 |
Февраль 2020 г. |
2 |
y2 |
x21 |
x22 |
… |
x25 |
Март 2020 г. |
3 |
y3 |
x31 |
x32 |
… |
x35 |
… |
… |
… |
… |
… |
… |
… |
|
|
|
|
|
|
|
… |
n |
yn |
xn1 |
xn2 |
… |
xn5 |
В данном пособии изложена методика изучения абсолютных моментных одномерных временных рядов с равноотстоящими уровнями (далее – временных рядов).
1.3. Моделинестационарныхвременныхрядов
Значения уровней ряда формируются под влиянием различных факторов: общей тенденции, сезонных и циклических закономерностей изучаемого процесса, а также случайных воздействий.
11
Тренд представляет собой долговременную систематическую линейную или нелинейную компоненту, закономерно изменяющуюся во времени (тенденцию).
Сезонная составляющая – это периодически повторяющаяся компонента, обусловленная временем года.
Циклическая составляющая – это периодически повторяющаяся компонента, обусловленная технологическим процессом или экономическими закономерностями.
Эти три вида регулярных (систематических) компонент часто имеются в рядах одновременно. Например, потребление электроэнергии может возрастать из года в год (тренд), но при этом оно может содержать и сезонную составляющую (на спрос и объемы электропотребления могут оказывать влияние сезонные факторы). На объем потребления электроэнергии также могут влиять технологические циклы (циклическая составляющая).
Необязательно, чтобы в формировании значений уровней ряда участвовали все четыре вида факторов (табл. 1.3). Однако случайная составляющая присутствует всегда и является обязательной составляющей временного ряда, т.к. случайность непременно присутствует в производственном процессе.
Случайная компонента является действием большого числа относительно слабых, второстепенных факторов. Влияние каждого из второстепенных факторов на значение уровней ряда незначительно, но ощущается их суммарное воздействие.
Таблица 1.3
Факторы, влияющие на значение уровней временного ряда
Компонента |
Классификация |
Описание |
|
Систематическая |
Общая устойчивая долговременная |
Тренд, yt |
||
|
|
тенденция |
|
|
|
Сезонная, st |
Систематическая |
Достаточно регулярные периодиче- |
|
|
ские изменения, происходящие в каж- |
|
|
дом 12-месячном периоде из года в год |
|
|
|
12
|
|
Окончание табл. 1.3 |
|
|
|
Компонента |
Классификация |
Описание |
|
|
|
Циклическая, ct |
Систематическая |
Повторяющиеся подъемы и спады, |
|
|
проходящие 4 фазы: пик, рецессия, |
|
|
депрессия, подъем |
|
|
|
Случайная, et |
Нерегулярная |
Непредвиденные воздействия внеш- |
|
|
них факторов |
|
|
|
Перечисленные в табл. 1.3 компоненты уровней временного ряда носят теоретический характер и оцениваются статистическими методами. Разделение временного ряда на составляющие компоненты называется фильтрацией.
В зависимости от связи компонент временного ряда между собой могут быть построены следующие модели.
Модель временного ряда, которая представляет каждое значение уровня ряда как сумму указанных выше компонент, называется аддитивной. Согласно этой модели, любое значение временного ряда представляется в виде суммы
|
|
yt = yt +st +ct +et, |
(1.1) |
гдеyt – значение уровня временного ряда;
e
yt, st, ct, t – соответственно значения трендовой, сезонной, циклической и случайной компонент в момент времени t.
Значения et также называются остатками (случайными возмущениями) модели.
Аддитивная модель применяется в тех случаях, когда анализируемый временной ряд имеет приблизительно одинаковую амплитуду колебаний на протяжении всей длины ряда (рис. 1.3).
Наиболее общей является мультипликативная модель временного ряда, используемая при анализе ежемесячных, ежеквартальных и ежегодных данных по электропотреблению и потому чаще всего применяемая в статистических исследованиях.
13
|
1800 |
|
|
|
|
|
|
|
|
|
1600 |
|
|
|
|
|
|
|
|
t |
1400 |
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
1200 |
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
1 |
5 |
9 |
13 |
17 |
21 |
25 |
29 |
33 |
|
|
|
|
t |
|
|
|
|
|
Рис. 1.3. График временного ряда с постоянной амплитудой |
|||||||||
|
|
|
колебаний уровней |
|
|
|
|||
В мультипликативной модели значения его уровней представляется в виде произведения
|
|
yt = yt st ct et. |
(1.2) |
Мультипликативная модель применяется в тех случаях, когда анализируемый временной ряд имеет возрастающую (уменьшающуюся) амплитуду колебаний на протяжении всей длины ряда (рис. 1.4).
yt |
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
1500 |
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
500 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
1 |
5 |
9 |
13 |
17 |
21 |
25 |
29 |
33 |
|
|
|
|
t |
|
|
|
|
Рис. 1.4. График временного ряда с непостоянной амплитудой колебаний уровней
14
Для проверки адекватности построенной модели вместо мультипликативной удобнее использовать смешанную модель с аддитивными остатками:
|
|
yt = yt st ct +et. |
(1.3) |
15
2. ПРЕДВАРИТЕЛЬНЫЙАНАЛИЗ ВРЕМЕННОГОРЯДА
2.1. Графическийанализвременногоряда
Анализ временного ряда начинается с построения и анализа его графика. График временного ряда – ломаная линия. На горизонтальной оси откладываются моменты времени t, на вертикальной оси – значения соответствующих уровней yt временного ряда.
Визуальный анализ графика временного ряда позволяет определить наличие и характер тренда, наличие сезонных и циклических колебаний, предполагаемых точек смены тенденции, выбросов (аномальных значений). Предварительные выводы, сделанные по виду графика, впоследствии проверяются статистическими расчетами.
Выявление сезонных (циклических) колебаний удобно проводить методом «наложения». Когда на горизонтальной оси откладываются периоды времени предполагаемого периода колебаний (например, месяцы или кварталы года), а по вертикальной оси – значения соответствующих уровней за несколько лет. При наличии сезонной (циклической) компоненты графики будут подобными (в части наличия пиков или спадов).
Пример.
Имеются месячные значения показателя «Электропотребление региона» (МВт ч) за 2018–2023 гг. (табл. 2.1). Провести графический анализ ряда.
|
|
|
|
|
|
|
Таблица 2.1 |
|
|
|
|
Исходные данные |
|
|
|||
|
|
|
|
|
|
|
|
|
|
Период |
|
Электропо- |
|
Период |
|
Электропо- |
|
|
|
|
требление, |
|
|
|
требление, |
|
Год |
Месяц |
t |
Год |
Месяц |
t |
|||
МВт ч, (yt) |
МВт ч, (yt) |
|||||||
|
|
|
|
|
|
|||
2018 |
январь |
1 |
407,1 |
2021 |
январь |
37 |
727,3 |
|
2018 |
февраль |
2 |
451,9 |
2021 |
февраль |
38 |
868,7 |
|
|
|
|
|
|
|
|
|
|
16
Окончание табл. 2.1
|
Период |
|
Электропо- |
|
Год |
Месяц |
t |
требление, |
|
МВт ч, (yt) |
||||
|
|
|
||
2018 |
март |
3 |
495,3 |
|
2018 |
апрель |
4 |
516,9 |
|
|
|
|
|
|
2018 |
май |
5 |
523,8 |
|
|
|
|
|
|
2018 |
июнь |
6 |
567,7 |
|
|
|
|
|
|
2018 |
июль |
7 |
479,7 |
|
|
|
|
|
|
2018 |
август |
8 |
463,7 |
|
2018 |
сентябрь |
9 |
439,7 |
|
2018 |
октябрь |
10 |
441,2 |
|
|
|
|
|
|
2018 |
ноябрь |
11 |
470,3 |
|
|
|
|
|
|
2018 |
декабрь |
12 |
543,7 |
|
|
|
|
|
|
2019 |
январь |
13 |
468,1 |
|
|
|
|
|
|
2019 |
февраль |
14 |
486,6 |
|
2019 |
март |
15 |
514,4 |
|
|
|
|
|
|
2019 |
апрель |
16 |
601,0 |
|
|
|
|
|
|
2019 |
май |
17 |
503,4 |
|
|
|
|
|
|
2019 |
июнь |
18 |
548,9 |
|
|
|
|
|
|
2019 |
июль |
19 |
601,5 |
|
|
|
|
|
|
2019 |
август |
20 |
568,6 |
|
2019 |
сентябрь |
21 |
571,3 |
|
|
|
|
|
|
2019 |
октябрь |
22 |
632,3 |
|
|
|
|
|
|
2019 |
ноябрь |
23 |
589,7 |
|
|
|
|
|
|
2019 |
декабрь |
24 |
796,3 |
|
|
|
|
|
|
2020 |
январь |
25 |
530,2 |
|
|
|
|
|
|
2020 |
февраль |
26 |
620,0 |
|
2020 |
март |
27 |
612,3 |
|
|
|
|
|
|
2020 |
апрель |
28 |
697,8 |
|
|
|
|
|
|
2020 |
май |
29 |
647,9 |
|
|
|
|
|
|
2020 |
июнь |
30 |
663,2 |
|
|
|
|
|
|
2020 |
июль |
31 |
716,0 |
|
|
|
|
|
|
2020 |
август |
32 |
736,0 |
|
2020 |
сентябрь |
33 |
731,7 |
|
|
|
|
|
|
2020 |
октябрь |
34 |
674,1 |
|
|
|
|
|
|
2020 |
ноябрь |
35 |
712,9 |
|
|
|
|
|
|
2020 |
декабрь |
36 |
1059,9 |
|
|
|
|
|
|
Период |
|
Электропо- |
|
Год |
Месяц |
t |
требление, |
|
МВт ч, (yt) |
||||
|
|
|
||
2021 |
март |
39 |
965,2 |
|
2021 |
апрель |
40 |
1041,2 |
|
|
|
|
|
|
2021 |
май |
41 |
954,3 |
|
|
|
|
|
|
2021 |
июнь |
42 |
1058,5 |
|
|
|
|
|
|
2021 |
июль |
43 |
1150,0 |
|
|
|
|
|
|
2021 |
август |
44 |
1072,0 |
|
2021 |
сентябрь |
45 |
1053,0 |
|
2021 |
октябрь |
46 |
1539,9 |
|
|
|
|
|
|
2021 |
ноябрь |
47 |
1507,3 |
|
|
|
|
|
|
2021 |
декабрь |
48 |
2013,0 |
|
|
|
|
|
|
2022 |
январь |
49 |
1345,3 |
|
|
|
|
|
|
2022 |
февраль |
50 |
1587,0 |
|
2022 |
март |
51 |
1655,2 |
|
|
|
|
|
|
2022 |
апрель |
52 |
1723,0 |
|
|
|
|
|
|
2022 |
май |
53 |
1701,5 |
|
|
|
|
|
|
2022 |
июнь |
54 |
1923,9 |
|
|
|
|
|
|
2022 |
июль |
55 |
1852,6 |
|
|
|
|
|
|
2022 |
август |
56 |
1987,6 |
|
2022 |
сентябрь |
57 |
2003,3 |
|
|
|
|
|
|
2022 |
октябрь |
58 |
1971,8 |
|
|
|
|
|
|
2022 |
ноябрь |
59 |
2027,0 |
|
|
|
|
|
|
2022 |
декабрь |
60 |
2492,7 |
|
|
|
|
|
|
2023 |
январь |
61 |
1793,8 |
|
|
|
|
|
|
2023 |
февраль |
62 |
2146,9 |
|
2023 |
март |
63 |
2422,3 |
|
|
|
|
|
|
2023 |
апрель |
64 |
2380,8 |
|
|
|
|
|
|
2023 |
май |
65 |
2315,4 |
|
|
|
|
|
|
2023 |
июнь |
66 |
2632,8 |
|
|
|
|
|
|
2023 |
июль |
67 |
2584,1 |
|
|
|
|
|
|
2023 |
август |
68 |
2606,7 |
|
2023 |
сентябрь |
69 |
2547,5 |
|
|
|
|
|
|
2023 |
октябрь |
70 |
2485,6 |
|
|
|
|
|
|
2023 |
ноябрь |
71 |
2611,8 |
|
|
|
|
|
|
2023 |
декабрь |
72 |
3406,5 |
|
|
|
|
|
17
Для проведения расчетов исходный ряд запишем в виде табл. 2.2.
|
|
|
|
|
Таблица 2.2 |
|
|
|
Временной ряд |
|
|
|
|
|
|
|
|
|
|
|
t |
yt |
t |
yt |
t |
|
yt |
1 |
407,1 |
25 |
530,2 |
49 |
|
1345,3 |
|
|
|
|
|
|
|
2 |
451,9 |
26 |
620,0 |
50 |
|
1587,0 |
|
|
|
|
|
|
|
3 |
495,3 |
27 |
612,3 |
51 |
|
1655,2 |
|
|
|
|
|
|
|
4 |
516,9 |
28 |
697,8 |
52 |
|
1723,0 |
|
|
|
|
|
|
|
5 |
523,8 |
29 |
647,9 |
53 |
|
1701,5 |
|
|
|
|
|
|
|
6 |
567,7 |
30 |
663,2 |
54 |
|
1923,9 |
|
|
|
|
|
|
|
7 |
479,7 |
31 |
716,0 |
55 |
|
1852,6 |
|
|
|
|
|
|
|
8 |
463,7 |
32 |
736,0 |
56 |
|
1987,6 |
|
|
|
|
|
|
|
9 |
439,7 |
33 |
731,7 |
57 |
|
2003,3 |
|
|
|
|
|
|
|
10 |
441,2 |
34 |
674,1 |
58 |
|
1971,8 |
|
|
|
|
|
|
|
11 |
470,3 |
35 |
712,9 |
59 |
|
2027,0 |
12 |
543,7 |
36 |
1059,9 |
60 |
|
2492,7 |
|
|
|
|
|
|
|
13 |
468,1 |
37 |
727,3 |
61 |
|
1793,8 |
|
|
|
|
|
|
|
14 |
486,6 |
38 |
868,7 |
62 |
|
2146,9 |
15 |
514,4 |
39 |
965,2 |
63 |
|
2422,3 |
|
|
|
|
|
|
|
16 |
601,0 |
40 |
1041,2 |
64 |
|
2380,8 |
|
|
|
|
|
|
|
17 |
503,4 |
41 |
954,3 |
65 |
|
2315,4 |
|
|
|
|
|
|
|
18 |
548,9 |
42 |
1058,5 |
66 |
|
2632,8 |
|
|
|
|
|
|
|
19 |
601,5 |
43 |
1150,0 |
67 |
|
2584,1 |
|
|
|
|
|
|
|
20 |
568,6 |
44 |
1072,0 |
68 |
|
2606,7 |
|
|
|
|
|
|
|
21 |
571,3 |
45 |
1053,0 |
69 |
|
2547,5 |
|
|
|
|
|
|
|
22 |
632,3 |
46 |
1539,9 |
70 |
|
2485,6 |
|
|
|
|
|
|
|
23 |
589,7 |
47 |
1507,3 |
71 |
|
2611,8 |
|
|
|
|
|
|
|
24 |
796,3 |
48 |
2013,0 |
72 |
|
3406,5 |
|
|
|
|
|
|
|
График временного ряда изучаемого показателя «Электропотребление региона» изображен на рис. 2.1.
На рис. 2.1 видны периодические колебания графика. Наиболее выраженные колебания повторяются через каждые 12 моментов времени, что в нашем случае соответствует 12 месяцам. Между этими 12 месяцами также есть незначительные
18
колебания, повторяющиеся с определенной периодичностью. Таким образом, можно сделать предположение, что данный временной ряд не случаен и эти колебания – годовые, полугодовые и квартальные изменения заявленной мощности. Также можно заметить, что, начиная с момента времени t = 32, график возрастает быстрее, а амплитуда колебаний увеличивается, что является косвенными признаками смены тенденции.
|
4000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Электропотребление, МВт ч |
3500 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2500 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1500 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
500 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
6 |
11 |
16 |
21 |
26 |
31 |
36 |
41 |
46 |
51 |
56 |
61 |
66 |
71 |
|
|
|
|
|
|
t – номер месяца |
|
|
|
|
|
||||
Рис. 2.1. График временного ряда электропотребления
На рис. 2.2 построен график электропотребления методом наложения. На графике прослеживается наличие сезонности, причем для последних трех лет сезонность проявляется сильнее.
ч |
4000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
МВт |
3500 |
|
|
|
|
|
|
|
|
|
|
|
3000 |
|
|
|
|
|
|
|
|
|
|
|
|
Электропотребление, |
|
|
|
|
|
|
|
|
|
|
|
|
2500 |
|
|
|
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|
|
|
|
1500 |
|
|
|
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
|
|
|
|
|
|
500 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
|
|
|
|
|
|
t – номер месяца |
|
|
|
|
|||
2018 |
2019 |
2020 |
2021 |
2022 |
2023 |
Рис. 2.2. График временного ряда электропотребления
19
2.2. Выявлениеаномальныхзначений
Уровни временного ряда могут содержать аномальные значения («выбросы»), которые могут возникать из-за ошибок двух типов.
Ошибки 1-го рода (устранимые) возникают при сборе, записи и передаче данных.
Ошибки 2-го рода (неустранимые) вызываются объективными причинами (аварии, катастрофы, кризисы, природные явления и др.).
Выбросы искажают результаты моделирования. Поэтому выявление и тестирование выбросов – обязательный этап предварительного анализа временного ряда.
При наличии аномальных значений их необходимо заменить реальными близлежащими или расчетными значениями. Чаще всего их заменяют средней арифметической двух соседних уровней ряда либо соответствующим значением по кривой, аппроксимирующей данный временной ряд. В случае ошибок 2-го рода помимо «сглаживания» выбросов производится анализ как самих значений, так и вызвавших их причин.
В процессе визуального анализа выявляются точки графика, «подозрительные» на выброс. Для тестирования подозрительного значения yt можно использовать метод Ирвина.
Рассчитывается наблюдаемое значение критерия (статистика):
u = |
|
|
yt −yt−1 |
|
, |
(2.1) |
|
|
|||||
|
|
|||||
|
|
|
|
|||
t |
|
sy |
|
|
|
|
|
|
|
|
|
|
|
гдесреднеквадратическое отклонение sy рассчитывается с помощью формул:
|
|
|
|
n |
|
|
n |
|
|
|
|
|
|
∑(yt −y)2 |
|
|
∑yt |
|
|
s |
y |
= |
|
t=1 |
; y = |
t=1 |
. |
(2.2) |
|
n −1 |
|
||||||||
|
|
|
|
|
n |
|
|||
|
|
|
|
|
|
|
|||
20
