УДК 621.311
Н. В. Токочакова, Д. Р. Мороз, А. С. Фиков
МОДЕЛИРОВАНИЕ РЕЖИМОВ ПОТРЕБЛЕНИЯ ЭЛЕКТРИЧЕСКОЙ ЭНЕРГИИ ДЛЯ ЗАДАЧ УПРАВЛЕНИЯ ЭНЕРГЕТИЧЕСКОЙ ЭФФЕКТИВНОСТЬЮ ПРОМЫШЛЕННЫХ ПОТРЕБИТЕЛЕЙ
Предложена структурная схема управления энергетической эффективностью промышленных потребителей, основным этапом в которой является оценка ее текущего состояния. Предложен универсальный алгоритм построения расчетно-статистических моделей режимов электропотребления, позволяющий получать модели требуемого качества и адекватности.
Введение
Высокая энергоемкость национального валового продукта Республики Беларусь в условиях дефицита топливноэнергетических ресурсов (ТЭР) при непрерывном росте мировых цен на них ставит перед республикой задачу повышения эффективности использования энергетических ресурсов во всех отраслях народного хозяйства. В настоящий момент данная задача решается за счет стимулирования всех субъектов народного хозяйства к энергосбережению со стороны государства, что отражено [1]. Такое стимулирование проявляется в ежегодном установлении для каждого субъекта хозяйствования квартальных норм расхода энергоресурсов на выпуск продукции и целевого показателя по энергосбережению. Поскольку превышение установленных норм расхода энергетических ресурсов на выпуск продукции приводит к штрафным санкциям, то все потребители ТЭР заинтересованы в поиске энергосберегающих мероприятий. Кроме того, все промышленные потребители (1111) с объемом потребления ТЭР свыше 1,5 тыс. т у. т. обязаны проходить энергетический аудит, основной целью которого являются оценка текущего состояния энергетической эффективности (ЭЭФ) и разработка плана внедрения энергосберегающих мероприятий на предстоящую пятилетку. Таким образом, управление ЭЭФ ПП приобретает актуальность и особую значимость как в общегосударственном масштабе, так и для ка-
ждого ПП в рамках действующего законодательства.
Постановка задачи
Управление ЭЭФ, как и любое другое управление, включает в себя три основных этапа: первый - оценка текущего состояния ЭЭФ; второй - принятие решения о методе воздействия на ЭЭФ; третий - непосредственное воздействие на ЭЭФ (внедрение энергосберегающего мероприятия). На рис. 1
представлена структурная схема управления ЭЭФ ПП. Из рисунка видно, что управление ЭЭФ имеет циклический характер, а приостановка цикла производится в случае внедрения в производство всех научно-технических достижений, позволяющих снизить энергоемкость продукции. Непрерывное развитие научно-технического прогресса обеспечивает постоянное наличие путей повышения ЭЭФ.
Очевидно, что основным этапом управления ЭЭФ является оценка ее текущего состояния, поскольку от этого зависят все управленческие решения. Поэтому вопросу оценки текущего состояния ЭЭФ уделяется большое внимание. В соответствии с [1] основным показателем текущего состояния ЭЭФ на предприятии является достигнутая величина абсолютного или удельного расхода ТЭР на производство продукции любого назначения.
Выбор метода воздействия на ЭЭФ
Воздействие на ЭЭФ (внедрение энергосберегающего мероприятия)
Оценка текущего состояния ЭЭФ
-----------------------------
Ґ \
Процесс производства
продукции
\ У
I
Получение продукции с наименьшими затратами
а
с
с
е
егр
о
о
г
гко
с
е
н
X
е
н
I
о
н
ч
чуа
н
н
$
и
т
с
о
Рис. 1. Структурная схема управления ЭЭФ ПП: I N - информационные потоки; - физическая взаимо-
связь
Проведенные исследования для Гомельской области [2] выявили, что в подавляющем количестве случаев (92 %) расчет удельных расходов ТЭР на производство продукции и, как следствие, разработка норм расхода ТЭР производятся расчетно-аналитическим методом.
Расчетно-аналитический метод предусматривает определение удельного расхода электрической энергии (ЭЭ) расчетным путем по статьям расхода на основе электрических балансов. Удельный расход ЭЭ, определенный расчетно-аналити-
ческим методом для существующих условий работы, будет представлять собой лишь единичное (среднее) значение удельного расхода ЭЭ из множества возможных. В ходе расчетов учитываются реальные условия производства с учетом фактического состояния оборудования и режима его работы, действительных параметров производственных процессов.
Учет режимов работы каждой единицы технологического оборудования для большинства ПП предполагает большую трудоемкость расчетов и сложность сис-
темы сбора параметров режимов электропотребления. Поэтому для ПП, работающих в условиях изменяющейся производственной программы и имеющих значительный парк технологического оборудования, данный метод определения удельного расхода ЭЭ на производство продукции неприемлем [2].
Поскольку применяемые методы оценки текущего состояния ЭЭФ оказываются весьма трудоемкими и имеют значительную погрешность, то снижается качество управления ЭЭФ. Возникает задача повышения достоверности определения уровня ЭЭФ ПП в современных условиях.
Научная новизна
Режимы работы большинства современных ПП отличаются нестабильностью вследствие зависимости от поставщиков сырья и рынков сбыта продукции. Так, в рамках одного года могут быть периоды с высокой и низкой степенями загрузки производства, кроме того, может изменяться качество исход-
ного сырья и т. д. Все это приводит к появлению большого многообразия режимов потребления ЭЭ, а оценка состояния ЭЭФ по среднему значению приводит к недопустимой ошибке. Для достоверной оценки текущего состояния ЭЭФ необходим учет множества режимов потребления ЭЭ. Такое возможно при оценке текущего состояния ЭЭФ на основе расчетностатистических моделей режимов потребления ЭЭ [3].
На уровень электропотребления ПП оказывает влияние большое количество факторов (объем производства продукции, климатические условия, качество исходного сырья и конечной продукции и др.). В
[3] предлагаются мультипликативные модели режимов электропотребления, где в качестве воздействующих факторов рассмотрены объем выпускаемой продукции и общий расход ЭЭ. В качестве отклика рассмотрен удельный расход ЭЭ. Поскольку общий расход ЭЭ и объем выпускаемой продукции линейно взаимосвязаны, то авторами предложена мультипликативная модель. Однако применение такой модели для оценки ЭЭФ затруднено из-за громоздкости решения и сложности интерпретации результатов моделирования. Кроме этого, рассмотренная модель не позволяет выявлять степень влияния факто-
ров на режим электропотребления.
Режим потребления ЭЭ ПП формируется под воздействием большого количества факторов, поэтому в качестве модели режимов потребления ЭЭ предлагается регрессионная аддитивная модель электропотребления от воздействующих факторов вида
ш = £д • X, + , (1)
і=1
где Ш - электропотребление ПП, кВт-ч; Х, - фактор, влияющий на уровень потребления ЭЭ ПП; п - количество факторов, влияющих на уровень потребления ЭЭ ПП; А, - степень влияния і-го фактора на уровень потребления ЭЭ ПП; Шпост - постоянная составляющая электропотребления.
В общем виде алгоритм построения модели режимов электропотребления ПП представлен на рис. 2. Как видно из рисунка, первоочередным шагом на пути построения модели режимов электропотребления являются изучение технологического процесса производства продукции и выявление факторов, оказывающих значительное влияние на режим потребления ЭЭ ПП.
Рис. 2. Логическая схема построения модели режимов электропотребления 1111
В общем случае на режим электропотребления оказывает влияние достаточно большое количество факторов, для выявления наиболее значимых могут быть применены методы факторного анализа. Однако изучение технологического процесса в некоторых случаях позволяет выявить значимые факторы на основе аналитических зависимостей, описывающих взаимосвязь электропотребления и технологических параметров (предприятия транспорта нефти, химической промышленности и др.). Количество факторов, подлежащих включению в модель в качестве воздействующих, должно быть минимальным, а сама модель при этом должна адекватно описывать режим потребления ЭЭ.
Исследования показывают, что для описания режимов электропотребления большинства 1111, выпускающих промышленную продукцию, достаточно учета одного фактора - объемов выпускаемой продукции. Сама модель при этом обладает достаточной степенью адекватности для задач нормирования и прогнозирования электропотребления, а также оценки ЭЭФ.
Для построения модели режимов электропотребления необходима статистическая информация в достаточном объеме. Уровень дискретизации статистики должен соответствовать требованиям, предъявляемым к качеству модели, при этом он должен быть одинаков для всех факторов, входящих в модель. В случаях, когда система учета не позволяет получать статистическую информацию одинаковой дискретизации, ее уровень принимается минимальным. Значения факторов (воздействие которых на электропотребление накапливается во времени, например, объем выпуска продукции), имеющие большую дискретизацию, пересчитываются по выражению
/=]-т+1
X ,
=0-1)т+1
(2)
накоплении информации; ] - номер значения фактора, приведенного к уровню минимальной дискретизации; т - отношение дискретизации фактора х в ИБД к дискретизации фактора Х, т. е. количество накопленных значений фактора х за период накопления следующего значения фактора с минимальной дискретизацией.
Значения факторов (воздействие которых на электропотребление усредняется во времени, например, влажность сырья, температура окружающей среды и др.), имеющие большую дискретизацию, пересчитываются по выражению
і і=]-т+1
X =1 XV
т і=0-1)т+1
(3)
где Х - значение фактора при уменьшении дискретизации; х - значение фактора при
Опыт построения расчетностатистических моделей режимов потребления ЭЭ [2, 3] показывает, что в подавляющем большинстве случаев система учета расхода ЭЭ и технологических факторов позволяет накапливать ИБД с уровнем дискретизации в одни сутки.
Необходимый объем статистической информации зависит от количества факторов, включаемых в модель, и целей построения модели [4]. Наиболее полной для задач нормирования, прогнозирования расхода ЭЭ, а также оценки текущего состояния ЭЭФ считается ИБД с глубиной в один год и уровнем дискретизации в одни сутки. После накопления достаточного объема статистической информации методом наименьших квадратов производится построение регрессионной модели режимов электропотребления вида (1).
Поиск коэффициентов регрессии удобно вести в матричной форме [5, 6]. Пусть глубина ИБД равна п значений, число факторов, включаемых в модель -к, т. е. имеется п значений независимых переменных f1, /2,... , /к и соответствующие им значения зависимой переменной W, тогда уравнение (1) запишется в виде
0 = B • F + E.
(4)
где 0 - вектор значений зависимой переменной, 0 = (0,Ж2,...,Жи)т ; В - вектор искомых коэффициентов регрессии,
В = (Р1,Р2,---,вк)Т; Е - вектор случайных отклонений, Е = (є1, є 2,..., є п )т; Е - матрица значений независимых переменных размером п х (к +1):
" f11 f 12 •• f1k
F = f21 f22 • •• f2k
frn f2k • fnk
Оценка В для вектора В ведется так, чтобы разность квадратов между вектором 0 зависимых переменных и вектором
оценок 0 = ВЕ была минимальной:
0-0 = 0-B-F
3.
Решением является оценка (если ранг матрицы Е равен к +1):
B = (ft •F)-1 FT 0.
(5)
Одним из показателей адекватности модели служит коэффициент детерминации Я2, который показывает долю дисперсии, учитываемой данной моделью в общей дисперсии независимой переменной:
R 2 = 1 sM )
К®-®)’.
(б)
Однако значение R2 возрастает с ростом числа факторов к и объемом выборки n, что не говорит об улучшении адекватности модели. Поэтому для сравнения моделей вводится скорректированный (adjusted) коэффициент детерминации. Его использование более корректно для сравнения регрессий при изменении
числа переменных (регрессоров) и количества их измерений [б]:
Ri= 1 -(1 - R )2
n -1 n - k -1
(7)
При формальной оценке надежности модели различными статистическими критериями «ускользают» те реальные причинно-следственные связи между параметрами, которые определяются особенностью исследуемого процесса. Как правило, в прикладных исследованиях для выявления количественной меры влияния важна оценка именно причинных связей. Анализ структуры связей производится на основе теоретических или профессиональных предположений о наличии и направлении связи. Таким образом, количество независимых переменных, введенных в модель, может быть обоснованным с использованием аналитических зависимостей и соответствовать цели создания модели [8].
Качество построенной модели характеризуется как максимальной относительной погрешностью 5шах, так и среднеквадратичным отклонением относительной погрешности а:
5 max с
= max і оі =
100-
Wі
; (В)
і=1,2, ...,n
а =
ГК F -о)2.
1 І=1
(9)
При проведении регрессионного анализа для ограниченной по объёму совокупности параметры уравнения регрессии, коэффициенты детерминации могут быть сильно искажены действием случайных факторов. В этих условиях необходима дополнительная проверка модели на адекватность. При этом производится проверка значимости (существенности) каждого коэффициента регрессии. Значимость коэффициен-
2
S
1
тов линейной регрессии определяют с помощью 1-критерия Стьюдента [5, 6]:
tвk = |Рк
(10)
Вычисленные по формуле (9) значения сравнивают с критическими t, которые определяют по таблице Стьюдента [5, 6] с учетом принятого уровня надежности 1-а и объемом выборки.
Уровень надежности 1-а обычно принимают равным 0,95. Параметр признаётся значимым (существенным) при условии, если ^асч > ^абл . В таком случае практически невероятно, что найденное значение коэффициента регрессии вк обусловлено только случайными совпадениями.
Для проверки значимости уравнения регрессии вычисляется статистический критерий Фишера (Е-критерий) [5, 6]:
- я)
- я п - к
!(( - я )2
к -1
(11)
Рассчитанное значение Е-критерия сравнивается с табличным [5, 6] при (к-1), (п-к) степенями свободы и заданном уровне надежности 1-а. И при превышении рассчитанного значения над табличным практически невероятно, что все найденные значения коэффициентов регрессии В обусловлены только случайными совпадениями:
Е > Е1-а(-1, п - к).
(12)
В случаях, когда построенная модель не удовлетворяет требуемому качеству, процесс построения возвращается на этап определения вида модели и факторов, в нее включаемых. На этом этапе вводятся корректировки с учетом полученных ранее результатов и процесс построения модели повторяется до тех пор, пока качество не удовлетворит требованиям.
Далее производится оценка текущего состояния ЭЭФ ПП. Наиболее распространенным показателем ЭЭФ является УРЭ, поэтому для оценки текущего состояния ЭЭФ предпочтительным является переход от модели режимов электропотребления к модели УРЭ от воздействующих факторов. Кроме того, определив с помощью модели влияние каждого из факторов на УРЭ
[4], а также возможность контролируемого воздействия на сами факторы, возможно прогнозирование изменения ЭЭФ в результате внедрения энергосберегающих мероприятий.
Рассмотрим процесс построения модели режимов потребления ЭЭ на примере предприятий трубопроводного транспорта нефти:
1) изучение технологического процесса транспортировки нефти [9] по трубопроводам позволило выделить в качестве факторов, воздействующих на режимы потребления ЭЭ предприятия, следующие: Р - грузооборот нефти, тыс. т-км/сут; ёэ - эквивалентный диаметр нефтепровода, мм; V - вязкость перекачиваемой нефти, м2/с;
2) на годовом интервале производилось формирование суточных временных рядов расхода ЭЭ и воздействующих факторов:
к 1 .......■. {р}=1,2.„■
Vt }=1,2...п"’ Vэt1 =1
=1,2,...,п" =1,2,...,п".
Накопленная статистическая информация вследствие несовершенства системы учета технологических параметров содержит выбросы. Для построения модели необходима очистка временных рядов от выбросов и календарных эффектов:
к 1=1,2,.,п' , }=1,2,.. ,,п',
^11=1,2,.. ,,п', Vэt \t=1,2,...,п',
п' = п"-ЫК,
где ЫК - количество суток, исключенных из временных рядов.
Для уменьшения влияния стохастической составляющей производится сглаживание временных рядов методом скользящего среднего с периодом усреднения т:
1 ,+Т-1 ____ 1 ,+Т-1
Ж, = , Р, = -^р ,
т гт
1 І+Т-1 ____ 1 І+Т-1
^ , аэ і = -Е аэ ,■;
3) методом наименьших квадратов производится построение модели режимов электропотребления от воздействующих факторов. Результат построения модели представлен на рис. 3;
4) производится оценка качества построенной модели по выражениям (8) и (9), а проверка значимости коэффициентов регрессии с использованием Е-критерия и ,-критерия при уровне надежности а = 0,95. Результаты проверки модели приведены в табл. 1.
600000 кВт-ч/сут -
400000 300000 Ж 200000 100000 Н о
Ж = 30,0238-Е + 8812-10б-у-655059< -47432
--------1---------1---------1---------1---------1---------1---------1---------1
38500 39000 39500 40000 40500 41000 тыс. т-км/сут 42500
Р ---------►
нефти
Рис 3. Результат построения модели режимов электропотребления для предприятий транспорта
Табл. 1. Результаты проверки качества модели режимов потребления ЭЭ и значимости коэффициентов регрессии
1 = І
І =І
Наименование критерия Условие удовлетворения Значение критерия
Максимальная относительная погрешность £шах Для задач нормирования и прогнозирования показателей ЭЭФ менее 10 % 5тах = 2,50 %
Среднеквадратическое отклонение относительной погрешности ст Для задач нормирования и прогнозирования показателей ЭЭФ менее 5 % а = 1,03 %
Е-критерий Для 230 наблюдений при уровне значимости 0,95 более 3,04 Е = 1434
,-критерий Для 230 наблюдений при уровне значимости 0,975 более 1,96 І І = = б 5 4, І < = 5
Как видно из табл. 1, построенная предприятия транспорта нефти полно-
модель режимов электропотребления для стью удовлетворяет предъявляемым
требованиям. Что свидетельствует о возможности ее использования для задач оценки текущего состояния ЭЭФ.
Выводы
1. Рассмотрена актуальность проблем повышения ЭЭФ 1111 в Республике Беларусь. Предложена структурная схема управления ЭЭФ 11 , основным этапом в которой является оценка текущего стояния ЭЭФ.
2. Рассмотрены существующие методы оценки ЭЭФ ПП. Показано преимущество оценки ЭЭФ на основе расчетностатистических моделей режимов электропотребления.
3. Предложен универсальный алгоритм построения расчетно-статистических моделей режимов электропотребления, позволяющий получать модели требуемого качества и адекватности.
СПИСОК ЛИТЕРАТУРЫ
1. Об энергосбережении : Закон Респ. Беларусь от 15 июля 1998 г. № 190-3. - Минск : Комитет «Белэнергосбережение», 1998. - 15 с.
2. Токочакова, Н. В. Анализ показателей энергетической эффективности промышленных потребителей / Н. В. Токочакова, Д. Р. Мороз // Вестн. ГГТУ им. П. О. Сухого. - 2006. - № 3. -С. 66-76.
3. Оперативное прогнозирование режимов электропотребления / Р. С. Саитбаталова [и др.] // Промышленная энергетика. - 2000. - № 6. -С. 27-29.
4. Токочакова, Н. В. Математическое обеспечение задач энергетического обследования технологического процесса транспортировки нефти / Н. В. Токочакова, А. С. Фиков // Вестн. МГТУ. - 2006. - № 2. - С. 169-177.
5. Ван дер Варден, Б. Л. Математическая статистика / Б. Л. Ван дер Варден. - М. : Изд-во иностранной литературы, 1960. - 434 с.
6. Козлов, М. В. Введение в математическую статистику / М. В. Козлов, А. В. Прохоров. -М. : МГУ, 1987. - 264 с.
7. Факторный, дискриминантный и кластерный анализ / Дж. Ким [и др.] ; под ред. И. С. Енюкова. - М. : Финансы и статистика, 1989.
8. Перегудов, Ф. И. Введение в системный анализ : учеб. пособие для вузов / Ф. И. Перегудов, Ф. П. Тарасенко. - М. : Высш. шк., 1989. - 367 с.
9. Анищенко, В. А. Способ построения модели режимов электропотребления участка нефтепровода / В. А. Анищенко, Н. В. Токоча-кова, А. С. Фиков // Изв. высш. учеб. заведений и энергет. об-ний СНГ. - Энергетика. - 2006. -№ 6. - С. 44-48.
Гомельский государственный технический университет им. П. О. Сухого
Материал поступил 07.02.2007
N. V. Tokochakova, D. R. Moroz, A. S. Fikov The modeling of consumption regimes of electrical energy for the task of control industrial user’s energy effectiveness
Gomel State Polytechnic University named after P. O. Sykhoi
The structural diagram of control of industrial user’s energy effectiveness has been offered, and the estimation of its current state is the main stage. Universal algorithm of the construction calculational-statistical models of the regimes of the electric requirements, it allows to receive required quality and adequacy models.