УДК 621.31+519.2
МОДЕЛИРОВАНИЕ СИСТЕМЫ КОНТРОЛЯ ЗАТРАТ ДЛЯ ПОДДЕРЖКИ ПРИНЯТИЯ РЕШЕНИЙ ПРИ УПРАВЛЕНИИ ЭЛЕКТРОСЕТЕВЫМ
ПРЕДПРИЯТИЕМ
А.М. КУМАРИТОВ, А.Э. ДЗГОЕВ, О.Р. БАБОЧИЕВ, И.М. ХУЗМИЕВ
Северо-Кавказский горно-металлургический институт (Государственный технологический университет) СКГМИ (ГТУ), г. Владикавказ
С использованием метода планирования эксперимента получена полезная адекватная модель, связывающая финансовые затраты по передаче данных о потреблении электроэнергии с количеством опросов разных категорий измерительных приборов учета электроэнергии. Определены категории измерительных приборов учета расхода электроэнергии и выявлена категория, которая имеет большую зависимость затрат от количества опросов.
Ключевые слова: математическая статистика, планирование эксперимента, оптимизация затрат на передачу данных электропотребления, настройка схемы опроса измерительных приборов.
Введение
Автоматизированная система контроля и учета электроэнергии (АСКУЭ) широко применяется в электроустановках промышленных и электросетевых предприятий, а также в жилищно-коммунальном хозяйстве - там, где необходим точный и своевременный учет перетоков и потребления электроэнергии. Применение АСКУЭ увеличивает эффективность процесса коммерческого и технического учета электроэнергии, повышая финансовые показатели эксплуатирующего субъекта.
Степень финансовой эффективности применения внедренных систем автоматизации напрямую зависит от затрат на их использование. Так, при оценке затрат на эксплуатацию АСКУЭ, в случаях использования дорогостоящего оборудования, каналов связи и неоптимальных тарифов на передачу данных о потреблении электроэнергии, констатируется превышение эксплуатационных затрат над финансовыми издержками процесса учета в условиях отсутствия АСКУЭ, что, соответственно, делает использование системы экономически нецелесообразным. Поэтому одно из направлений развития АСКУЭ - это минимизация стоимости внедрения и эксплуатации системы без ухудшения функциональных характеристик.
Одна из основных затратных составляющих применения АСКУЭ - это использование платных каналов связи сторонних провайдеров телекоммуникационных услуг. Особенно это актуально для горных территорий, где в ряде случаев отсутствует инфраструктура дешёвых и быстрых сетей связи, и единственным приемлемым вариантом остается сотовая связь.
Расходы на связь в АСКУЭ зависят от многих факторов. Основным фактором является количественно-временная характеристика использования каналов связи, определяющая величину объемов передаваемых данных и длительность интервалов времени передачи данных. Задача АСКУЭ - сбор информации с измерительных приборов учета. От того, насколько эффективно будет настроена схема опроса объектов учета, зависит как объем собранной информации (собираемость данных), с одной стороны, так и количественно-временная характеристика использования каналов
© А.М. Кумаритов, А.Э. Дзгоев, О.Р. Бабочиев, И.М. Хузмиев Проблемы энергетики, 2015, № 9-10
связи, с другой стороны. Для улучшения собираемости данных с измерительных приборов эксплуатирующие субъекты вынуждены наращивать технические и материальные ресурсы (число обслуживаемых серверов, пропускную способность каналов связи, посредством которых осуществляется передача данных и т.д.), что связано с существенными финансовыми затратами.
Деятельность подобного рода систем осуществляется в соответствии с требованиями к их функционированию, предъявляемыми законодательством РФ и нормативными документами в электроэнергетической отрасли: в пункте 6 «Организация сбора и использования коммерческой информации от АСКУЭ для расчетов на оптовом рынке» положения «Об организации коммерческого учета электроэнергии и мощности на оптовом рынке и приложениях»; в пункте 9 «Общие технические требования к системе учета электроэнергии» типовой инструкции по учету электроэнергии при ее производстве, передаче и распределении РД34.09.101-94; в пункте «Организация системы измерений и сбора информации о фактическом производстве электрической энергии и мощности и об их потреблении на оптовом рынке» 3 статьи 33 Федерального закона «Об электроэнергетике» от 30.12.2012 N 291-ФЗ.
Другой важный фактор - это ценовая политика провайдеров (операторов), гибкость тарифной системы и т.п.
Указанные два фактора могут быть взаимоувязаны путем проведения системного анализа и принятия по его результату эффективных управленческих решений с конечной целью оптимизации технико-экономических показателей использования АСКУЭ.
Таким образом, важной и актуальной прикладной задачей для электросетевых организаций и иных субъектов, использующих АСКУЭ, является нахождение баланса между увеличением собираемости данных о потреблении электроэнергии и допустимыми финансовыми затратами с целью максимизации экономической эффективности применения систем автоматизации учета.
Методика исследования
Данная задача решена в результате использования метода планирования эксперимента [1, 2]. Для этой цели был использован план Бокса (В3), обладающий Б-оптимальностью и являющийся наиболее мощным критерием оптимальности [3].
Выбор плана Бокса, в отличие от критериев Л-оптимальных, .-оптимальных, С-оптимальных, ^-оптимальных, связан с тем, что он относится к разряду О-оптимальных критериев, имеющих максимальное значение определителя матрицы М или минимальное значение М-1 [2].
Используемый программный продукт
Все необходимые формулы для расчётов приведены в терминах программного продукта «Ма{Неай».
Для расчёта коэффициентов уравнения регрессии был использован МНК (Метод Наименьших Квадратов) с учётом линеаризации независимых переменных в квадратичной форме. Коэффициенты регрессии рассчитаны по формуле
В = (ХТХ) ~'ХТ1, (1)
где В - коэффициенты регрессионной модели; X - матрица независимых переменных; У - матрица-столбец зависимой переменной; (хТX] - левая часть системы нормальных уравнений.
Основные результаты исследования и их обсуждение
Для составления матрицы планирования эксперимента решено было разделить общее количество измерительных приборов учёта электропотребления на три категории: в безразмерном масштабе - Х1, Х2, Х3; в размерном масштабе - А, В, С.
Безразмерный масштаб
Х1 (Категория Х1) - критические измерительные приборы учета электропотребления, опрашиваемые каждый час или каждые 30 минут.
К ним относятся точки учета, преимущественно расположенные на крупных объектах, подстанциях, информация о потреблении электроэнергии с которых имеет чрезвычайно важное значение как для измерения и оперативного контроля за потреблением электроэнергии, так и для надежного функционирования оптового рынка электроэнергии.
Х2 (Категория Х2) - измерительные приборы учета электропотребления, опрашиваемые один раз в сутки.
Точки учета, опрос которых осуществляется раз в сутки, можно отнести к так называемому мелкомоторному сектору потребителей электроэнергии. В данную категорию входят трансформаторные подстанции (ТП), малые и средние промышленные предприятия и т.д.
Х3 (Категория Х3) - измерительные приборы учета, опрашиваемые один раз в трое суток.
К точкам учета, опрашиваемым один раз в трое суток, относятся приборы учета, принадлежащие потребителям частного сектора (многоквартирные дома, частные хозяйства и т.д. ).
Размерный масштаб
А - интервал опроса измерительных приборов учета электропотребления от 0 до 60 мин, с частотой 1 час или 30 минут.
В - интервал измерительных приборов учета электропотребления от 0 до 24 часов, с частотой 1 раз в сутки.
С - интервал опроса измерительных приборов учета от 0 до 3 суток. Х1, Х2, Х3 - независимые переменные в безразмерном масштабе,
А, В, С - независимые переменные в размерном масштабе,
AS, BS, CS - средние значения,
У - финансовые затраты (руб.) на установление связи точек учета с сервером АСКУЭ посредством применяемых каналов связи.
Используя план Бокса, в матрице планирования проведено 14 опытов.
Так как N = 2к + 2 • к , то N = 23 + 2 • 3 = 14 опытов, где к - число независимых переменных Х.
Матрица планирования и результаты эксперимента в размерном и безразмерном масштабах представлены в табл. 1.
Далее представлены расчеты соответствующих значений опытов (с 1 -го по 14-ый), результаты которых внесены в табл. 1.
Лщп:= 1 4^:= 60 АБ := ^ + Атах = 30,5 ДА := Л^х -^ = 29,5
А - АЧ А — АЧ
._ ^^тт ^ _ _1 ._ тах _1
тт'= ДА " тах'= ДА "
:= 1 В_ := 24 ВБ := Втп + Втх = 12,5 ДВ := - В8 = 11,5
Х2 •= Вт1п ~ В$ _ 1 у ^ = Втах ~ В$ _ 1
т1п' АВ тах' АВ
0™:=! Стах:=72 С5:= Стп +Стах =36,5 АС:= С^ -СЯ = 35,5 С — С — с$
уп ._ Ст1п Со _ 1 уп ._ Стах _ 1
т1п'= АС " тах'= АС "
Таблица 1
Матрица планирования и результаты эксперимента в размерном и безразмерном масштабах
Размерный Безразмерный Затраты
масштаб масштаб Функциональная
№ А, В, С, Х1 Х2 Х3 У, роль матрицы
опыта мин. час. сут. руб.
1 1 1 1 -1 -1 -1 163.3
2 60 1 1 +1 -1 -1 176.3
3 1 24 1 -1 +1 -1 191.4
4 60 24 1 +1 +1 -1 352.0 линейная часть
5 1 1 72 -1 -1 +1 135.2 матрицы
6 60 1 72 +1 -1 +1 195.0
7 1 24 72 -1 +1 +1 278.2
8 60 24 72 +1 +1 +1 266.2
9 1 1 36,5 -1 0 0 135.2
10 60 1 36,5 +1 0 0 127.4
11 30,5 12,5 36,5 0 -1 0 113.9 нелинейная
12 30,5 12,5 36,5 0 +1 0 117.5 часть
13 30,5 12,5 1 0 0 -1 112.0 матрицы
14 30,5 12,5 72 0 0 +1 115.2
Безразмерный масштаб
Проведен расчет коэффициентов регрессионной модели в безразмерном масштабе.
N := 14 - число опытов,
к := 10 - число коэффициентов регрессии, так как решено было использовать нелинейный полином второй степени с тремя независимыми переменными XI, Х2, Х3.
X0 : =
(1 1 (-1 ^ (-11
1 1 -1
1 -1 1
1 1 1
1 -1 -1
1 1 -1
1 -1 1
XI : = X2 :=
1 1 1
1 -1 0
1 1 0
1 0 -1
1 0 1
1 0 0
V1 V V 0 V V 0 V
X3 : =
(-1 1 -1 -1 -1 1 1 1 1 0 0 0 0 -1
V 1 V
Y : =
(163.3^1
176.3
191.4 352.0 135.2 195.0 278.2 266.2 135.2
127.4 113.9
117.5 112.0
V 115.2,
где X0 - фиктивная переменная (для формирования свободного члена коэффициента регрессии B0);
X1, X2, X3 (независимые переменные в безразмерном масштабе) и Y (зависимая переменная) указаны выше в тексте данной статьи.
Для формирования матрицы независимых переменных (Х) в программном продукте MathCad была использована функция augment.
Г 2 2 2
X:= а^шеп(Х0,Х1 ,Х2,Х3,Х1 ,Х2 ,Х3 ,(Х1-Х2),(Х1-Х3),(Х2-Х3)_. Далее проведен расчет коэффициентов регрессионной модели по формуле (1).
B := (XT-х) XT-Y =
0
0 70.45
1 21.36
2 42.16
3 -0.52
4 60.85
5 45.25
6 43.15
7 9.475
8 -15.725
9 1.3
Коэффициенты регрессионной модели
Ранги независимых переменных (по модулю)
Определены ранги независимых переменных по силе их влияния на зависимую переменную уравнения регрессии нелинейного полинома второй степени с тремя независимыми переменными.
Ранги независимых переменных (X.1, X2, X3) рассчитаны для того, чтобы понять силу их влияния на зависимую переменную (Г). Результат представлен в табл. 2. В2 =42,16 В =21,36 В3 =-0,52
Таблица 2
Ранги независимых переменных
Ранг Сила влияния Безразмерный масштаб
I |42,16| X2
II |21,36| X1
III |-0,52| X3
По силе влияния на затраты по передаче данных независимые переменные ранжируются следующим образом: Х2, XI, Х3.
Далее была применена формула для расчетных значений зависимой переменной:
УЯ := X • В. (2)
Для определения тесноты связи между зависимыми переменными (Г) и расчетными значениями зависимой переменной (УЯ) рассчитан коэффициент корреляции [4]:
г := согг(УЯ,У) = 0,914 - коэффициент корреляции между расчетными и экспериментальными данными
Полная модель в безразмерном масштабе независимых переменных Разработана полная регрессионная модель в безразмерном масштабе переменных:
У := В0 + В1 • X1 + В2 • X2 + В3 • X3 + В4 • XI2 + В5 • X22 +
В6 • X32 + В7 • X!- X2 + В8 • X1- X3 + В9 • X2 • X3. (3)
Оценка значимости коэффициента корреляции была произведена с помощью критерия Стъюдента.
JN - 2
-— = 7,805 - расчетное значение критерия Стьюдента,
1 - r
1; := д? (0,975, N - 2) = 2,179 - табличное значение критерия Стьюдента. В связи с тем, что й > ^ коэффициент корреляции признан статистически значимым с 95% надежностью. На этом же основании уравнение регрессии может быть признано адекватным экспериментальным данным.
Оценка адекватности модели с помощью коэффициента детерминации Я2 (ЯЯ) Одной из эффективных оценок адекватности регрессионной модели, мерой качества уравнения регрессии, характеристикой прогностической силы анализируемой
регрессионной модели является коэффициент детерминации (R2). Величина R2 показывает, какая часть (доля) вариации зависимой переменной обусловлена вариацией
объясняющей переменной. Чем ближе R2 к единице, тем лучше регрессионная модель аппроксимирует эмпирические данные [5].
Zy
YSR := —— = 177,057 - среднее значение зависимой переменной,
N
Z (YR - YSR)2
RR := —--— = 0,835 - коэффициент детерминации (определенности),
(Y - YSR)2
•v/RR = 0,914 - коэффициент корреляции.
В связи с тем, что л/RR = 0,914 равняется коэффициенту корреляции между векторами Y и YR (r := corr(YR,Y) = 0,914), принимается гипотеза об адекватности регрессионной модели экспериментальным данным как по t-критерию Стьюдента, так и по коэффициенту детерминации.
Размерный масштаб
Далее представлены расчеты разработки регрессионной модели в размерном масштабе независимых переменных.
Таблица 3
Матрица планирования и результаты эксперимента в размерном и безразмерном масштабах
Размерный Безразмерный Затраты
масштаб масштаб Функциональная
№ A, B, C, X1 X2 X3 Y, роль матрицы
опыта мин. час. сут. руб.
1 1 1 1 -1 -1 -1 163.3
2 60 1 1 +1 -1 -1 176.3
3 1 24 1 -1 +1 -1 191.4
4 60 24 1 +1 +1 -1 352.0 линейная часть
5 1 1 72 -1 -1 +1 135.2 матрицы
6 60 1 72 +1 -1 +1 195.0
7 1 24 72 -1 +1 +1 278.2
8 60 24 72 +1 +1 +1 266.2
9 1 1 36,5 -1 0 0 135.2
10 60 1 36,5 +1 0 0 127.4
11 30,5 12,5 36,5 0 -1 0 113.9 нелинейная часть
12 30,5 12,5 36,5 0 +1 0 117.5 матрицы
13 30,5 12,5 1 0 0 -1 112.0
14 30,5 12,5 72 0 0 +1 115.2
X0 : =
A : =
AW
1 " ( 1 " ( 1 " ' 163.3"
60 1 1 176.3
1 24 1 191.4
60 24 1 352.0
1 1 72 135.2
60 1 72 195.0
1 60 B : = 24 24 C : = aw 72 72 Y : = 278.2 266.2
1 12.5 36.5 135.2
60 12.5 36.5 127.4
30.5 1 36.5 113.9
30.5 24 36.5 117.5
30.5 30.5; 12.5 v 12.5 , 1 V 72 , 112.0 v 115.2,
где A, B, C независимые переменные в размерном масштабе.
Для формирования матрицы независимых переменных (Х) в программном продукте MathCad была использована функция augment.
Г 2 2 2 -> -> ->"
£:= augment_X0,A,B,C,A ,B ,C ,(A-B),(AC),(B C)
Далее по формуле (1) рассчитаны коэффициенты регрессионного уравнения (регрессионной модели)
B1 := (XT-х) XT-Y =
0
0 162.582
1 -3.342
2 -5.856
3 -2.096
4 0.07
5 0.342
6 0.034
7 0.028
8 -0.015
9 3.184-10-3
Коэффициенты регрессионной модели в размерном масштабе независимых переменных
Далее по формуле (2) проведен расчет значений зависимой переменной: Ж := X • 51.
Полная модель в размерном масштабе независимых переменных Полная модель в размерном масштабе независимых переменных представлена формулой
Y ;= B0 + B1 • A + B2 • B + B3 • С + B4 • A2 + B5 • B2 +
+B6 • С2 + B7 • A • B + B8 • A • С + B9 • B • С. (4)
r := corr(YR,Y) = 0,914 - коэффициент корреляции между расчетными (ГЛ) и экспериментальными данными (Г).
Выводы
По результатам планирования эксперимента и создания на его основе регрессионной модели была установлена зависимость затрат на опрос точек учета от количества опросов и категорий измерительных приборов учета электроэнергии.
В ходе ранжирования трех категорий точек учета (X1, X2, X3) были получены важные для дальнейшего исследования результаты.
Из них следует, что наиболее влияющими на затраты являются точки учета, принадлежащие к мелкомоторному сектору B (Категория Х2), на втором месте расположились точки учета сектора А (Категория Х1), являющиеся критичными точками учета , и замыкает тройку - точки учета частного сектора C (Категория Х3). Из чего можно сделать вывод, что в первую очередь необходимо оптимизировать опрос точек учета мелкомоторного сектора (B).
Полученная зависимость затрат от типа точки опроса объясняется следующими причинами.
По группе «С» затраты наименьшие в связи с наиболее оптимальной схемой опроса, реализованной в АСКУЭ сектора населения. В процессе сбора данных задействованы механизмы сжатия и пакетирования информации. Данные передаются единовременно по большим совокупностям индивидуальных приборов учета, при этом применяется тарификация с оплатой исключительно объемов переданных данных без учета времени на передачу и количества попыток установления связи удаленной точки с сервером АСКУЭ. Также на снижении затрат сказывается низкая периодичность сбора данных. В результате, несмотря на наибольшее число точек учета, задействованных в данной группе, затраты на эксплуатацию каналов связи из расчета на одну точку учета минимальны относительно остальных групп.
По группе «А» основным фактором затрат является большая периодичность опроса - один час и чаще. При этом для повышения надежности сбора используется более дорогостоящий способ соединения по стандартному GSM каналу с повременной тарификацией вызова и с зависимостью затрат от количества попыток соединения. Относительные затраты из расчета на одну точку учета по данной группе максимальны, однако ввиду малого количества точек учета величина суммарных затрат группы «А » занимает промежуточную позицию.
По группе «В» суммарные затраты максимальны, основной причиной чего является значительное число точек учета, оборудованных индивидуальными терминалами связи с персональным абонентским номером. При этом часть терминалов, в связи с недостаточным качеством предоставления услуги доступа в Интернет, функционирует в стандартном режиме передачи данных. Необходимость ежесуточного обновления данных и неоптимальная конфигурация расписания опроса также способствуют увеличению затрат.
Summary
Using the experiment design method an useful adequate model connecting the financial expenditure on data transference in terms of the electric energy consumption with
the counting amount of various measuring devices was developed. The categories of the measuring devices registrating the energy expenditure were determined and the category having the greatest dependence of expenditure on the counting amount was revealed.
Keywords: mathematical statistics , experiment design , optimization of data transfer costs of electricity, setting up schemes measuring devices.
Литература
1. Налимов В.В., Чернова Н.А. Статистические методы планирования экстремальных экспериментов. / М.: Издательство «Наука». Главная редакция физико-математической литературы. 1965.
2. Новик Ф.С., Арсов Я.Б. Оптимизация процессов технологии металлов методами планирования экспериментов / М.: Машиностроение; София: Техника, 1980. 304 с.
3. Винарский М.С., Лурье М.В. Планирование эксперимента в технологических исследованиях./ «Техшка», 1975. 168 с.
4. Кумаритов А.М., Дзгоев А.Э., Бетрозов М.С., Бочкор Е.Н. Анализ малой выборки экспериментальных данных при управлении энергоснабжением и энергосбережением региона // Известия высших учебных заведений. Проблемы энергетики. № 9-10. 2013. С. 81-88.
5. Кремер Н.Ш., Путко Б.А. Эконометрика: Учебник для вузов / Под ред. проф. Н.Ш. Кремера. М.: ЮНИТИ-ДАНА, 2002. 311 с.
Поступила в редакцию 11 февраля 2016 г.
Кумаритов Алан Мелитонович - д-р техн. наук, профессор, заведующий кафедрой «Информационные системы в экономике» Северо-Кавказского горно-металлургического института (Государственного технологического университета) СКГМИ (ГТУ), г. Владикавказ. E-mail: [email protected].
Дзгоев Алан Эдуардович - канд. техн. наук, доцент кафедры «Информационные системы в экономике» Северо-Кавказского горно-металлургического института (Государственного технологического университета) СКГМИ (ГТУ), г. Владикавказ. E-mail: [email protected].
Бабочиев Олег Рамазанович - аспирант кафедры «Информационные системы в экономике» СевероКавказского горно-металлургического института (Государственного технологического университета) СКГМИ (ГТУ), г. Владикавказ. E-mail: [email protected].
Хузмиев Игорь Маратович - канд. техн. наук, заместитель начальника службы реализации услуг и учета электроэнергии Северо-Осетинского филиала, ОАО "МРСК Северного Кавказа". Тел.: 8 (918) 820-23-24. E-mail: [email protected].