Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Математические методы анализа и прогнозирования в электроэнергетике. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
691 Кб
Скачать

В многомерном ряду представлена динамика нескольких показателей, характеризующих один процесс (явление). Например, показатели деятельности предприятия, наблюдаемые помесячно, можно представить в виде многомерного ряда (табл.  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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]