УДК 621.311
УЧЕТ КЛИМАТИЧЕСКИХ УСЛОВИЙ РЕСПУБЛИКИ БЕЛАРУСЬ ПРИ ПРОГНОЗИРОВАНИИ И НОРМИРОВАНИИ ПОТРЕБЛЕНИЯ ЭЛЕКТРИЧЕСКОЙ ЭНЕРГИИ ПРЕДПРИЯТИЯМИ ТРУБОПРОВОДНОГО ТРАНСПОРТА НЕФТИ
А. С. ФИКОВ, С. В. БОРСУК
Учреждение образования «Гомельский государственный технический университет имени П. О. Сухого»,
Республика Беларусь
Введение
В Концепции Программы [1] социально-экономического развития Республики Беларусь на 2006-2010 гг. в числе важнейших проблем отмечается высокая энергоемкость выпускаемой продукции, что приводит к перерасходу ТЭР. Одним из инструментов государственной политики в области энергосбережения является нормирование расхода ТЭР на выпуск продукции. В целях регламентирования единых методических подходов к нормированию и прогнозированию расходов топливно-энергетических ресурсов на производство единицы продукции введено в действие положение о нормировании расхода топлива, тепловой и электрической энергии в народном хозяйстве Республики Беларусь
[2]. Согласно этому документу, основной задачей нормирования расхода ТЭР является обеспечение применения при планировании производства продукции технически обоснованных норм расхода топлива, тепловой и электрической энергии (ЭЭ).
В настоящее время для предприятий трубопроводного транспорта нефти Республики Беларусь прогнозирование и нормирование потребления ЭЭ производится с использованием регрессионной модели зависимости электропотребления от технологических факторов [3]: грузооборота нефти, ее вязкости или температуры, эквивалентного диаметра нефтепровода:
W = р рр +р, •у +рй < +Po, катн/сут (1)
где Рр - коэффициент регрессии, кВт-ч/тыс. т-км; Р - грузооборот нефти, тыс. т-км/сут;
2 2 - коэффициент регрессии, кВт-ч-с/(м -сут); V - вязкость нефти, м /с; Ра -
коэффициент регрессии, кВт-ч/(м-сут); ёэ - эквивалентный диаметр нефтепровода, м; Р0
- свободный член уравнения регрессии, кВт-ч/сут.
Количество независимых переменных, включенных в модель, может варьироваться в зависимости от условий функционирования нефтепровода. Так, например, эквивалентный диаметр нефтепровода используется в модели только для прогнозирования и нормирования расхода ЭЭ в условиях проведения реконструкции линейной части нефтепровода.
Постановка задачи
Для упрощения процедур нормирования и прогнозирования авторами [3] предложено сезонное отклонение электропотребления от среднегодового значения учитывать поправочным коэффициентом, что позволяет определять квартальный расход ЭЭ, используя минимальное количество независимых переменных. Данный коэффициент для /-го квартала находится по следующему выражению:
где - фактическое потребление ЭЭ за і-й квартал, кВт-ч; - расчетное
электропотребление за і-й квартал, без учета сезонного изменения электропотребления (определяется по однофакторной модели зависимости электропотребления от грузооборота нефти), кВт-ч.
Для нормирования и прогнозирования ЭЭ в условиях стабильной динамики вязкости с учетом поправочного коэффициента модель имеет вид:
Следует отметить, что при таком подходе определения поправочного коэффициента в прогнозное значение электропотребления закладывается влияние сезонного изменения вязкости нефти предыдущего года. Такой подход к нормированию и прогнозированию расхода ЭЭ дает хорошие результаты в условиях стабильной динамики вязкости нефти
[3], поскольку в прогнозном периоде сохраняется величина сезонного изменения электропотребления, вызванного изменением вязкости нефти. Вязкость нефти, транспортируемой по территории Республики Беларусь, зависит от химического состава (месторождения) и ее температуры. Прогнозирование вязкости нефти в функции химического состава затруднительно, поскольку точно неизвестно из каких месторождений и в каких пропорциях будет поступать нефть на территорию республики в прогнозном периоде. В свою очередь температура нефти на территории Республики Беларусь формируется под воздействием температуры окружающей среды, которая, даже для таких значительных промежутков времени как квартал, не является постоянной величиной и на момент разработки норм не известна. Таким образом, нормирование и прогнозирование расхода ЭЭ ведется в условиях неопределенной исходной информации.
При разработке норм потребления ТЭР необходимо учитывать нормативные условия производства продукции и оказания услуг, регламентируемые техническими нормативными правовыми актами. К таким условиям, например, относятся: освещенность рабочего места, нормы расхода горячей воды на одного работника, температуры теплоносителя и окружающей среды. Применительно к технологическому процессу транспортировки нефти в качестве нормативных условий может быть принята температура окружающей среды (воздуха), регламентируемая [4]. Ставится задача разработать способ расчета нормативных поправочных коэффициентов сезонного отклонения электропотребления от среднего за год значения (НИК), исходя из среднеквартальных температур окружающей среды.
Предлагаемые способы решения
Для решения поставленной задачи в регрессионной модели (1) вязкость нефти целесообразно заменить на температуру воздуха. Модель электропотребления участка нефтепровода для определения НИК в условиях отсутствия реконструкции нефтепровода примет вид:
где РР - коэффициент регрессии при грузообороте, кВт-ч/тыс.т-км; Р - грузооборот нефти, тыс. т-км/сут; Р{ - коэффициент регрессии при температуре воздуха, кВт-ч/сут-°С; 1ос- температура воздуха, °С; Р0 - свободный член уравнения регрессии, кВт-ч/сут.
Алгоритм построения модели аналогичен рассмотренному в [3], при периоде усреднения статистической информации 90 суток.
Ж = kv-(р рР + Р0), кВт-ч/сут.
(3)
Ж = Р р- Р + Р (-ґ0с + р 0, кВт-ч/сут,
(4)
Замена в модели (1) вязкости нефти на температуру окружающей без изменения вида модели становится возможной вследствие линейности характера взаимозависимости параметров toс и V (рис. 2).
о
>
-15
-5
5
15
25
Рис. 2. Зависимость среднесуточной вязкости нефти участка нефтепровода «Мозырь - Адамова Застава» от среднесуточной температуры воздуха в г. Мозыре
Величины НИК предложено определять по выражению
90-Шос
Рр-Рср + Рх ^ср.кв, +Р
90-ЖХос х°ргод
Р р-Рср + Рх 'Хср.год +Р.
о. е.,
(5)
где W °с сркв/ - среднесуточное электропотребление при температуре окружающей среды,
равной среднеквартальной, кВт-ч/сут; Wt°'c Хсргод/ - среднесуточное электропотребление при температуре окружающей среды, равной средней за год по данным [4], кВт-ч/сут; Рср -среднесуточный грузооборот за год, тыс. т-км/сут; Хсркв. - среднеквартальная температура наружного воздуха /-го квартала, °С; Хсрго воздуха, °С.
Среднеквартальная температура Хср кв. количеству суток в месяце:
- средняя за год температура наружного определяется как средневзвешенная по
-т
у= ср.мес у з
°С,
(6)
т.
з =1 3
где
.срмес з - среднемесячная температура наружного воздуха за у-й месяц в /-м квартале, принимается по данным [4], °С; ту. - количество суток ву-м месяце, сут.
Ирогнозирование электропотребления в зависимости от временного цикла можно классифицировать на оперативное, краткосрочное, внутримесячное, текущее и долгосрочное, годовое, перспективное [3]. При этом нормирование расхода ТЭР относится к задаче текущего и долгосрочного прогнозирования. При прогнозировании расхода ЭЭ участка нефтепровода за осенне-зимний или весенне-летний период в основу построения модели (4) ложится статистика за соответствующий период предыдущего года. При этом значения НИК (нормативное отклонение среднего электропотребления за квартал от среднего за осенне-зимний или весенне-летний период) будут отличаться от полученных по годовой статистической информации. Определим аналитическую взаимосвязь НИК, полученных по годовой и сезонной статистике. Также покажем, что
0
0
3
сумма НПК (^с), определенных по сезонной статистике, равна 2, а сумма НПК (kV),
определенных по годовой статистике, равна 4.
Предполагая поквартальное равенство количества суток перекачки нефти по участку нефтепровода, среднесуточное электропотребление за год при неизменном среднесуточном грузообороте можно записать в виде:
Ж
ср.год
1 4
= 4ЕЖср.г, карч/сут
4 і=1
(7)
где Жср, - среднесуточное электропотребление за /-й квартал, кВт-ч/сут.
Аналогично, среднесуточные значения электропотребления по сезонам года примут вид:
Жср2,3 = 0,5(Жср.2 + Жср.3 ) , кВт^ч/сут;
ЖСр1,4 = 0,5(ЖСр.1 + ЖСр.4 ), кВт^ч/сут,
(8)
(9)
где Жср 1, Жср 2, Жср 3, Жср 4 - среднесуточные значения электропотребления соответственно
за первый, второй, третий и четвертый кварталы, кВт-ч/сут.
Учитывая, что среднесуточные за квартал значения электропотребления можно выразить через среднесуточное за год электропотребление и соответствующий НПК
( К. = К Ср.год ), формулы (8) и (9) можно переписать:
Жср2 3 = 0,5 • Жср.год •(кV2 + кс3), к&тч/суг;
• (^1 + kv 4 ), кВт^ч/сут,
Ж 14 = 0,5 • Ж
ср ’ ср.год
(10)
(11)
где kv1, kv2, kv3, кс4- НПК соответственно первого, второго, третьего и четвертого
кварталов, определенные по годовой статистической информации.
Величины НПК, характеризующие отклонение среднеквартального электропотребления от среднего за осенне-зимний или весенне-летний период, находятся по выражениям:
кс = кс1 -
Ж
"ср.1
Ж
1,4 ’ 42
Ж
к° = р2
ср
Жс
2,3 5 43
Ж
к с = р.3
ср
Жс
2,3 5 44
Ж
к с = р.4
ср
Жс
(12)
ср
Подставляя выражения (12) в (10) и (11) получим две системы уравнений, отражающие взаимосвязь между НПК, полученными по сезонной и годовой статистической информации:
кс = 2 •кс1
к с 1 + к V 4
к с = . 2 • кс 4
кс1 + к V 4
2 • к
кс = -
с2
кс 2 + кс 3
кс 2 + кс 3
(13)
Сложив левые и правые части уравнения первой системы, получим:
кс + кс = 2 •кс 1 +. ^1 + к с4 = к + к
1 ^ ^С4
2 • к
4
кс1 + кс 4
=2.
(14)
Справедливость равенства ксс2 + ксс3 = 2 доказывается аналогично.
Величины НПК, характеризующие отклонение среднеквартального электропотребления от среднего за год, находятся по выражению:
*„= . (15)
Сумма коэффициентов к ^ с учетом (7) и (15) запишется в виде:
Е к-=Е Ж,., ЖР„Д=4Е Ж,., / Е Ж*=4. (16)
Использование систем (13) позволяет получить значения коэффициентов кс непосредственно через значения коэффициентов к ^, минуя построение регрессионных моделей электропотребления по сезонной статистической информации. Тождества (14) и (16) позволяют производить проверку расчета коэффициентов к v и к °.
Исследовательская часть
Произведем расчет НПК для участков нефтепровода РУП «Гомельтранснефть Дружба»: «Унеча - Мозырь», «Мозырь - Адамова Застава», «Мозырь - Броды». На основе суточной статистической информации за 2006 г. сформирована информационная база данных, включающая: расход ЭЭ (Ж), грузооборот нефти (Р), температуру воздуха (7) в г. Мозыре за 2006 г. (табл. 1).
Таблица 1
Фрагмент информационной базы данных по режимам электропотребления участка нефтепровода «Мозырь - Адамова Застава»
Дата W, кВт^ч/сут Р, тыс. т^км/сут и О С ©
02.01.2006 928830 62061 -0,09
03.01.2006 987400 61854 -1,45
04.01.2006 987020 61828 -4,41
05.01.2006 992470 63409 -5,46
06.01.2006 1000420 63488 -6,30
С использованием метода наименьших квадратов и скользящего среднего с периодом усреднения 90 сут построены модели электропотребления участков нефтепровода вида
(4). Найденные коэффициенты уравнения регрессии сведены в табл. 2.
Таблица 2
Коэффициенты регрессии моделей электропотребления участков нефтепровода
РУП «Гомельтранснефть Дружба»
Наименование участка Коэффициенты рег рессии
Рр» кВт^ч/ тыс. т^км Р.» кВт^ч/сут^°С Ро, кВт^ч/сут
Унеча - Мозырь 18,30 -1509 -231767
Мозырь - Адамова Застава 31,74 -3498 -1016632
Мозырь - Броды 84,19 -159 -86197
Важным этапом построения модели является проверка качества уравнения регрессии. При увеличении грузооборота нефти затраты ЭЭ на транспорт нефти возрастают, что
выражается положительным коэффициентом регрессии при грузообороте. Знак коэффициента регрессии при температуре должен быть отрицательный, поскольку температура окружающей среды связана с вязкостью нефти обратно пропорциональной зависимостью. Рост вязкости нефти вызывает увеличение затрат ЭЭ на ее транспортировку. Знаки при полученных коэффициентах регрессии не противоречат физике протекания технологического процесса.
Проверка статистической значимости коэффициентов регрессии с т объясняющими переменными произведена на основе Х-статистики [5]. При уровне значимости а =0,05
наблюдаемое значение Х -статистики сравнивалось с критической точкой 7^2 п-т-1
распределения Стьюдента. Для участка «Мозырь - Адамова Застава» параметр Хр
составляет 50,29, что превышает критическое значение Ха/2п-т-1 (п = 252), равное 1,648.
Коэффициент регрессии является статистически значимым. Результаты проверки значимости коэффициентов регрессии сведены в табл. 3.
Таблица 3
Проверка статистической значимости коэффициентов уравнения регрессии
Наименование участка Расчетная .-статистика Табличное значение .-критерия
.р .о
Унеча - Мозырь 38,82 -24,68 -12,20 1,648
Мозырь - Адамова Застава 50,29 -44,30 -26,39
Мозырь - Броды 30,71 -7,83 -16,28
Общая проверка значимости уравнения регрессии произведена с использованием критерия Фишера [3]. Рассчитанное значение сравнивалось с табличным [5] при (к -1), (п - к) степенями свободы и уровне надежности 0,95. При превышении рассчитанным значением табличного, практически невероятно, что все найденные значения коэффициентов регрессии обусловлены только случайными совпадениями. Для участка «Мозырь - Адамова Застава» параметр F составляет 1477, что превышает критическое значение Р1-а(к-1п-к) (п = 252), равное 3,0. Следовательно, можно сказать об общей
статистической значимости коэффициентов уравнения регрессии. Результаты проверки общей значимости коэффициентов регрессии сведены в табл. 4.
Таблица 4
Проверка общей статистической значимости коэффициентов уравнения регрессии
Наименование участка Расчетная ^-статистика Табличное значение ^-критерия
Унеча - Мозырь 1065 3,0
Мозырь - Адамова Застава 1477
Мозырь - Броды 689
Качество построенной модели характеризуется как максимальной относительной
5П
а [3]:
так и среднеквадратичным отклонением относительной погрешности
5тах = тах\ 5, =
р
100-*-
%,
(17)
I =1,2
п
а =
п - к -
7-1(8,-8)
1 ,=1
%. (18)
Проверка качества моделей электропотребления участков нефтепровода приведена в табл. 5. Результаты проверки показывают высокую точность моделирования электропотребления по наиболее энергоемким участкам нефтепровода «Унеча - Мозырь» и «Мозырь - Адамова Застава». Для участка «Мозырь - Броды» точность моделирования несколько ниже, однако с учетом малой энергоемкости участка (5 % от суммарного технологического потребления ЭЭ) может считаться удовлетворительной.
Таблица 5
Проверка качества построения уравнения регрессии
Наименование участка 5тах, % С| %
Унеча - Мозырь 3,03 1,39
Мозырь - Адамова Застава 1,58 0,81
Мозырь - Броды 6,97 2,53
Полученные величины НПК при среднем грузообороте нефти по участкам нефтепровода приведены в табл. 6.
Таблица 6
Значения нормативных поправочных коэффициентов участков нефтепровода
РУП «Гомельтранснефть Дружба»
Наименование участка ^1 к 2 К 3 к,4 X к«
Унеча - Мозырь 1,034 0,981 0,969 1,015 4,0
Мозырь - Адамова Застава 1,042 0,976 0,961 1,019 4,0
Мозырь - Броды 1,023 0,986 0,978 1,010 4,0
Заключение
1. Предложен новый способ расчета нормативных поправочных коэффициентов (НПК), характеризующих нормативное отклонение среднесуточного электропотребления участка нефтепровода за квартал от среднесуточного за год, основанный на построении двухфакторной аддитивной регрессионной модели зависимости электропотребления участка нефтепровода от грузооборота нефти и температуры окружающей среды, отличающийся учетом температуры окружающей среды за многолетний период наблюдений, позволяющий более объективно прогнозировать и нормировать квартальный удельный расход ЭЭ на транспортировку нефти в условиях неопределенной исходной информации.
2. Предложены выражения, позволяющие получить значения НПК, характеризующие отклонение среднесуточного электропотребления участка нефтепровода за квартал от среднесуточного за осенне-зимний или весенне-летний период (к,), непосредственно через значения НПК, характеризующие отклонение среднесуточного электропотребления участка нефтепровода за квартал от среднесуточного за год (к,), без построения регрессионных моделей электропотребления за осенне-зимний или весенне-летний период.
Литература
1. О Республиканской программе энергосбережения на 2006-2010 годы : постановление Совета Министров Респ. Беларусь, 2 февр. 2006 г., № 137 // Нац. реестр правовых актов Респ. Беларусь. - 2006. - № 24. - 5/17219.
2. Об утверждении Положения о нормировании расхода топлива, тепловой и электрической энергии в народном хозяйстве Республики Беларусь : постановление Комитета по энергоэффективности при Совете Министров Республики Беларусь от 19 ноября 2002 г. № 9 // Нац. реестр правовых актов Респ. Беларусь. - 2006. - № 24. -5/17219.
3. Анищенко, В. А. Оценка и нормирование показателей энергоэффективности предприятий трубопроводного транспорта нефти / В. А. Анищенко, Н. В. Токочакова. -Гомель : ГГТУ им. П. О. Сухого, 2007. - 233 с.
4. Строительная климатология: СНБ 2.04.05 - 2000. - Введ. 1.07.2001. - Минск: Минстройархитектуры, 2001. - 37 с.
5. Бородич, С. А. Вводный курс эконометрики : учеб. пособие / С. А. Бородич. - Минск : БГУ, 2000.- 354 с.
Получено 11.05.2007 г.