Прогнозирование затрат электрической энергии в зависимости от объемов изготовленной продукции Текст научной статьи по специальности «Экономика и бизнес»

ЭНЕРГОСБЕРЕЖЕНИЕ

УДК 519.216.3:621.317.38

Левченко Сергій Андрійович, кандидат технічних наук, доцент

Запорізька державна інженерна академія. Кафедра електротехніки та енергетичного менеджменту, м. Запоріжжя, Україна. Пр Леніна 226, м. Запоріжжя, 69006. Тел.: +38 (063) 6045676.

E-mail: [email protected]

ПРОГНОЗУВАННЯ ВИТРАТ ЕЛЕКТРИЧНОЇ ЕНЕРГІЇ В ЗАЛЕЖНОСТІ ВІД ОБ’ЄМІВ

ВИГОТОВЛЕНОЇ ПРОДУКЦІЇ

В статті розглянуто прогнозування електроспоживання за допомогою кореляційно-регресійного аналізу та теорії ймовірностей на перспективу на прикладі киснево-компресорного цеху ПАТ “Запоріжсталь-АГ". Ключові слова: рівняння регресії, довірчі інтервали, нормальний розподіл.

Левченко Сергей Андреевич, кандидат технических наук, доцент

Запорожская государственная инженерная академия. Кафедра электротехники и энергетического менеджмента, г. Запорожье, Украина. Пр. Ленина 226, г Запорожье, 69006. Тел.: +38 (066) 6045676. E-mail:[email protected]

ПРОГНОЗИРОВАНИЕ ЗАТРАТ ЭЛЕКТРИЧЕСКОЙ ЭНЕРГИИ В ЗАВИСИМОСТИ ОТ ОБЪЕМОВ

ИЗГОТОВЛЕННОЙ ПРОДУКЦИИ

В статье рассмотрено прогнозирование электропотребления с помощью корреляционно-регрессионного анализа и теории вероятностей на перспективу на примере кислородно-компрессорного цеха ПАО “Запорожсталь-АГ".

Ключевые слова: уравнение регрессии, доверительные интервалы, нормальное распределение.

Sergiy Levchenko, Ph. D. in Engineering Science

Zaporizhzhya State Engineering Academy. Department of Electrotechnics and energy management, Zaporizhzhya, Ukraine. Leninaprosp. 226, Zaporizhye, 69006. Tel. +38 (063) 6045676. E-mail:[email protected]

PREDICTION POWER CONSUMPTION DEPENDING ON VOLUME OF MANUFACTURED

PRODUCTION

This article considers forecasting power consumption with using correlation regression analysis and probability theory on perspective on the example of the oxygen compressor plant of PSC "Zaporizhstal-AG".

Keywords: regression equation, confidential intervals, normal distribution.

Вступ

При нестабільності економічних умов практична потреба в об'єктивному співставленні обсягів виробництва та питомих витрат електроенергії на підприємствах велика, оскільки обумовлена необхідністю, по-перше, прогнозувати електроспоживання при значних змінах обсягів виробництва для формування заявки в енергосистему і, по-друге, оцінювати раціональність використання електроенергії.

В умовах функціонування оптового ринку електроенергії (ОРЕ), для великих і середніх споживачів електричної енергії гостро стоїть питання точності планування електроспоживання. Складання достовірної заявки на споживання електроенергії повинно дозволити знизити витрати на придбану електричну енергію.

Практика показує, що для прогнозування споживання електроенергії не існує загального, єдиного методу: кожне виробництво містить індивідуальні технологічні цикли, які, сумуючись, утворюють унікальний часовий процес. Однак у всіх виробничих циклах споживання енергії можна знайти спільні риси, тим самим утворюючи методичну базу для виконання точного прогнозу.

При наявності великої статистичної бази з електроспоживання та об’ємів виготовленої продукції найбільш доцільним методом прогнозування буде кореляційно-регресійний аналіз з можливістю визначити ймовірність такого прогнозу.

Аналіз літературних даних і постановка проблеми

Питаннями прогнозування електроспоживання займалися такі вчені, як Праховник А. В., Пащук В. Я., Копцев Л. А., Філіпов С. П., Кулік М. М. та інші.

12

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

ЭНЕРГОСБЕРЕЖЕНИЕ

Перед тим, як подати заявки з енергоспоживання на ОРЕ на перспективний період, підприємства мають дати попередню оцінку величині електроспоживання. Проте, найчастіше прогнозовані величини не відповідають дійсним фактичним витратам.

Для того, щоб застрахувати підприємство від можливих збитків, пов’язаних з понад договірними витратами електроенергії або перевищенням розміру замовлення, необхідно підвищити точність побудови прогнозу енергоспоживання.

Мета і задачі дослідження

При достатньому об’ємі статистичних даних, проблему підвищення точності прогнозування графіків електричного навантаження доцільно вирішити шляхом виявлення залежності між об’ємами виготовленої продукції та відповідними витратами електроенергії за допомогою кореляційно-регресійного аналізу.

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

Прогнозування електроспоживання за допомогою методів кореляційно-регресійного аналізу та теорії ймовірностей

В якості експериментальних даних, необхідних для проведення дослідження, використано статистичні дані киснево-компресорного цеху ПАТ “Запоріжсталь-АГ” з об’ємів виготовленого кисню та іншої продукції — стисненого повітря, азоту і аргону, та відповідних витрат електричної енергії за період з 2009 до вересня 2013 року за місяцями.

За допомогою засобів MS Excel в якості регресійного рівняння для математичного опису залежності між об’ ємами виготовленої продукції та відповідними витратами електричної енергії було підібрано степеневу функцію, яка найбільш точно відповідає статистичним даним:

У x = a • (1)

де x - об’єм виготовленої продукції;

у - витрати електричної енергії на виготовлення x об’єму продукції; a - коефіцієнт рівняння регресії, який показує витрати електроенергії на холостому ході; b - коефіцієнт рівнянні регресії, який показує середню зміну витрат електроенергії із зміною об’ єму виготовленої продукції на одну одиницю; є - випадкова помилка, що підлягає мінімізації.

Для оцінки параметрів регресійного рівняння використано підхід, заснований на лінеаризації моделі, який полягає в тому, що за допомогою належних перетворень вихідних змінних досліджувану нелінійну залежність представляють у вигляді лінійного співвідношення між перетвореними змінними, а саме — про логарифмувавши обидві частини рівняння та замінив змінні:

у * = a* + b * • x * +s*,

(2)

де у* = lgух ; x* = lgx; а* = lga; b* = b; є* = /g-є.

Далі за методом найменших квадратів даються оцінки параметрів даної лінійної функції:

*

a

b *

£ (lg Уі) •£ (lg x.)2 -£ (lg x. • !g у) •£ (lg x.)

n •£ (lg xi)2 - (£ (lg xi))2 £ (lg xi • !g Уг )- a * •£ (lg xi)

£ (lgxi)2

(3)

Параметри степеневої функції знаходяться відповідно, як 10a

та 10b .

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

13

ЭНЕРГОСБЕРЕЖЕНИЕ

Для оцінки суттєвості рівняння регресії в цілому розраховується F-критерій Фішера, необхідний для перевірки нульової гіпотези Но', факторна і залишкова дисперсії, обумовлені рівнянням регресії, не відрізняються одна від одної.

Рфакш _ £ (yx - У)2 • (п - 2)

Р.зал £ (У - Уx )2

(4)

Якщо розрахований критерій більше табличного, гіпотеза відхиляється.

Для оцінки тісноти зв’язку між енерговитратами та об’ємами продукції розраховується лінійний показник - індекс кореляції, значення якого знаходяться в межах 0 < R < 1,

R

1

II

£( у - Ух)2 £( у - y)2

(5)

Для оцінки долі дисперсії, поясненої рівнянням регресії, розраховується коефіцієнт детермінації, який свідчить про правильність підібраної моделі регресії і дорівнює квадрату індексу кореляції R2.

Якість вибраної моделі характеризує середня помилка апроксимації, значення якої мають знаходитися в межах 5-7%,

1

А _ -■£

п

(У - Ух)

У

•100 .

(6)

Для оцінки суттєвості коефіцієнта регресії його величина порівнюється з його стандартною помилкою, тобто визначається фактичне значення t-критерію Стьюдента,

tb _-b-, яке потім порівнюється з табличним значенням при певному рівні значущості а і

mb

числі ступенів свободи (n - 2).

Величина стандартної помилки спільно з t-розподілом Стьюдента при n - 2 ступенях свободи застосовується для перевірки суттєвості коефіцієнта регресії та для розрахунку його довірчих інтервалів,

ть _

£( у - y x)2

S2

---=———--------_ ---------- . (7)

І(п-2)•£(x-x)2 ]£(x-x)2

Якщо фактичне значення t-критерію перевищую табличне, тоді гіпотезу про незначимість коефіцієнту регресії можна відхилити.

Довірчий інтервал для коефіцієнту регресії визначається, як b ± t • mb.

Стандартна помилка параметру а визначається, як,

та _

£ (у - yx)2 £•

п - 2 п •£ (x - x)2 ^

S2 •■

£-

п

£(x - x)

2

(8)

Процедура оцінки значимості даного параметру така саме, як і для параметру b; обчислюється t-критерій, tа _-^, його величина порівнюється з табличним значенням при

т„

п - 2 ступінях свободи.

Стандартна помилка прогнозованого за лінією регресії значення визначається за формулою, де xk — значення змінної x,

т>. _ S •

1 + (xk - x )2

п £(x-x)2,

(9)

2

2

14

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

ЭНЕРГОСБЕРЕЖЕНИЕ

Для прогнозного значення У Хк довірчі інтервали визначаються наступним чином:

У.

хк

Л • тУх

< У < У + t

— У — у хк a

т ;

(10)

Далі за теорією ймовірностей визначається достовірність даного прогнозу, для чого в якості випадкової величини береться значення питомого електроспоживання.

Для знаходження ймовірності потрапляння величини питомого енергоспоживання в довірчі інтервали, розраховані за кореляційно-регресійним методом, спочатку за допомогою критерію х2 (“хі квадрат”) Пірсона перевіряється гіпотеза про нормальний розподіл даної величини.

Перевірка виконується за наступним алгоритмом.

1. Методом добутків, розділивши генеральну сукупність значень питомого електроспоживання на інтервали, обчислюються вибіркова середня хВ та вибіркове середнє

квадратичне відхилення ов:

х в = M * • h + C ; (11)

Db = М 2* - (М1*)2 J. h2 , (12)

де h - крок (різниця між двома сусідніми варіантами); C - хибний нуль (варіанта, яка розміщена приблизно в * V ni • иг * М1 = умовний момент першого порядку; М 2 n середині варіаційного ряду); V ni • иг2 й = умовний момент n

другого порядку; хі - С иг = умовна варіанта; h

n - емпірична частота потрапляння питомого енергоспоживання в інтервал. Для зменшення помилки, викликаної групуванням даних за інтервалами, при розрахунку дисперсії і відповідно середнього квадратичного відхилення роблять поправку Шепарда:

°В =^Db - (112) • h2 . (13)

2. Обчислюються теоретичні частоти потрапляння питомого електроспоживання в інтервал:

n • h ( )

Фг) ,

n =■

де n - об’єм вибірки (сума всіх частот); h - крок (різниця між двома сусідніми варіантами);

Хг - хв .

ОВ ’

(р(и) =

л/2"

ж

3. Емпіричні та теоретичні частоти порівнюються між собою, для чого: а) знаходять спостережуване значення критерію Пірсона:

(14)

(15)

(16)

иг =

1

и

е

X

2

спост

v (ni - n')2 ^ n'

(17)

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

15

ЭНЕРГОСБЕРЕЖЕНИЕ

б) за таблицею критичних точок розподілу х, за заданим рівнем значущості а і числу ступенів свободи k=S-3, де S - число груп вибірки, знаходять критичну точку %кр2 (a; k)

області прийняття нульової гіпотези.

Якщо - немає підстав відхилити гіпотезу про нормальний розподіл генеральної сукупності, тобто, емпіричні та теоретичні частоти відрізняються несуттєво {випадково).

Після підтвердження гіпотези, знаючи густину та функцію розподілу, можна визначити ймовірність потрапляння величини питомого електроспоживання в заданий довірчий інтервал {а; в):

P{a< X < Р)

Р ( x-a )2

Г е - [(2ст)]2 dx

a

Ф(Р-а л ф(а-а)

Ф(----) -Ф(---)-.

а

а

(18)

де а - вибіркова середня ^з;

ф(Р—а) і ф(а—а) - функції Лапласа, значення яких знаходяться за таблицею. а а

Отже, для киснево-компресорного цеху було виявлено наступні залежності витрат електричної енергії на виготовлення об’ємів кисню та іншої продукції відповідно:

у = 0,345x1,047 та y = 0,186x0,972.

Розрахунок необхідних коефіцієнтів та перевірка рівняння регресії в цілому показали правильність підбору функції та її параметрів.

Довірчі інтервали для даних рівнянь відповідно виглядають, як:

у - 2,005 • my x < у * < у + 2,005 • mу x;

Xk ’ У X У У xk

де стандартні помилки прогнозованих значень:

my = 0,024 •.

yx

+ (Xk - 4,55)2 56 0,342

m - = 0,011 •

yx

+ (Xk - 5,125)2 56 0,167

Графік залежності прогнозного енергоспоживання за даними Бізнес плану на 2013 -2016 роки від необхідних для виготовлення об’ємів кисню з урахуванням планових витрат та довірчих інтервалів за місяцями зображено на рис. 1, а для іншої продукції - на рис. 2.

Розраховані ймовірності потрапляння прогнозного питомого енергоспоживання в отримані довірчі інтервали: 0,87 та 0,92 відповідно підтверджують високу якість прогнозу.

Аналіз обох графіків показав, що планові витрати електричної енергії наближаються до верхнього довірчого інтервалу прогнозного значення. Проте, ймовірність їх потрапляння в даний інтервал буде значно меншою, що свідчить про доцільність використання прогнозних значень.

Висновки

Завдяки застосованому кореляційно-регресійному аналізу та теорії ймовірностей прогнозовані значення енергоспоживання максимально наближено до фактичних, що, в свою чергу, призведе до мінімізації грошових витрат, пов’ язаних з неточністю прогнозування.

Список використаних джерел

1. Єлісєєва, І. І. Економетріка: навч. посібник / С. В. Куришева, Т. В. Костєєва та ін.; під заг. ред. 1.1. Єлісєєвої. - 2-е вид. - М.: Фінанси і статистика, 2007. - 576 с.

2. Гмурман В. Є. Теорія ймовірностей і математична статистика: навч. посібник для вузів / В.Є. Гмурман - М., Вища школа, 2003. - 479 с.

3. Уварова Л. Г. Прогнозування графіків електричних навантажень підприємства / Л. Г. Уварова, С. А. Левченко//Вісник національного технічного університету «ХПІ». - 2011. - N 53. - С. 122-128.

16

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

ЭНЕРГОСБЕРЕЖЕНИЕ

“І-----------1-----------1-----------1------------r

10 20 30 40 50

2013-2016 роки за місяцями

Рис. 1. Графік електричного навантаження в залежності від об’ємів виготовленого

кисню на 2013-2016 роки

20 ЗО

20 13-2016 роки за місяцями

Рис. 2. Графік електричного навантаження в залежності від об’ємів виготовленої іншої продукції на 2013-2016 роки

Referenses

1. Eliseeva, 1.1. Econometrics [Ekonometrika]. 2007, 576 p.

2. Gmurman, V. E. Probability theory and mathematical statistics [Teoriya ymovirnostey i matematychna statystyka]. 2003, 479 p.

3. Uvarova, L. G. & Levchenko S. A. (2011), Prediction charts electrical loads enterprise [Prognozuvannya graflklv elektrichnih navantazhen pldpriyemstva], Proceedings of the National Technical University "HPT',.no. 53,

P. 122-128.

Поступила в редакцию 23.09 2014 г.

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

17