CC BY
221
УДК 621.311
Горячев В.Я., Михайлов С.А., Шатова Ю.А.
Пензенский государственный университет
МАТЕМАТИЧЕСКИЕ МЕТОДЫ ОПРЕДЕЛЕНИЯ ЦЕНТРА РАСПРЕДЕЛЕННЫХ ПО ПОВЕРХНОСТИ НАГРУЗОК
Аннотация. Обзор существующих аналитических способов определения местоположения центров электрических нагрузок. Предлагается метод определения координат центра электрических нагрузок на базе кривой нормального закона распределения плотности вероятности. Приводятся результаты расчетов в программе, написанной в среде MathLab.
Ключевые слова. Программирование, математическая модель, графический интерфейс.
Эффективность работы систем электроснабжения в значительной степени зависит от взаимного расположения подстанций и потребителей электрической энергии. Это связано, первую очередь, с потерями на передачу электрической энергии и со стоимостью работ по электрификации определенного региона. С другой стороны, эффективность работы электрических сетей зависит от наличия информационных систем и систем управления объектами. Разработка последних определяется правильным проектированием систем электроснабжения.
При проектировании современных систем электроснабжения всегда возникали задачи определения числа, расположения источников питания, распределения приемников электроэнергии по источникам питания и др. Решать перечисленные задачи становится все сложнее. Это объясняется тем, что проектировщикам приходится оперировать с большим количеством исходных данных, объем которых постоянно увеличивается. В первую очередь это относится к возросшему числу приемников электрической энергии. Большой объем данных и постоянный его рост привели к широкому внедрению вычислительной техники в проектную практику, что потребовало разработки иных подходов к проектированию. Первое представление о характере распределения нагрузок по территории объекта получают с помощью картограммы нагрузок [1] . В основе метода лежит принцип определения «центра масс». Для решения задач определения места оптимального расположения подстанции в методе нахождения «центра масс» нагрузок используется представление нагрузки в форме равновысоких цилиндров, объем которых пропорционален мощности соответствующей нагрузки. Возможности метода ограничены, так как позволяет определить место расположения лишь о д-ной подстанции для питания всех потребителей данного региона [3] . Практика проектирования систем электроснабжения показывает, что данная методика не дает однозначного ответа на вопрос оптимального расположения нескольких подстанций, если имеется некоторое количество потребителей, распределенных по площади, но частично «сгруппированных». Такая ситуация имеет место в реальной действительности, так как в пределах одной сети электроснабжения имеется ряд центров, вокруг которых и расположены потребители.
Для определения мест расположения нескольких подстанций в одной системе электроснабжения иногда используется метод эквипотенциальных поверхностей. При этом каждой нагрузке, расположенной в некоторой точке, присваивается потенциал пропорциональный мощности нагрузки. Центром нагрузки в этом случае считается точка с наивысшим потенциалом. Однако действенность такого метода весьма сомнительна, так как превалирующей всегда является нагрузка с максимальной мощностью.
В основу предлагаемого метода определения центров нагрузок для оптимального расположения подстанций предполагается в качестве базисной функции использовать уравнение поверхности :
f ( X, у)
1
2na2
(x-a)1 +(у-b)2 2a2
(2)
где: a - радиус рассеяния значений функции, a и b
координатах ху .
Определим объем фигуры, ограниченной поверхностью сначала проинтегрируем исходную функцию по x пределах
g(y)
■=ю х=ю (x-af+(y-bf
f f(x, y)dx = -----2 f e a dx
_•* 2na _■*
1
>/2na
e
(y-b)2 2a2
(3)
Проинтегрируем полученную функцию по у .
- координаты вершины поверхности в
f(x,у) и плоскостью ху , для чего от -ю до +Ю .
у=ю
v = f g(у)Ф/
у=-ю
1
•J2na
f
(у-b)2
e а dy = 1
(4)
Таким образом, объем тела, ограниченного поверхностью f(x, у) и плоскостью ху равен единице.
При a = 0 приведенная функция в точке х = a и у = b имеет значение, равное бесконечности, однако объем, ограниченный функцией и плоскостью ху равен единице. При а^ю данная функция в точке х = a и у = b имеет значение, стремящееся к нулю, однако объем, ограниченный функцией и плоскостью ху остается равным единице. При увеличении а происходит распределение объема по поверхности плоскости ху . При этом объем под поверхностью остается неизменным.
Коэффициент а является радиусом окружности, образованной точками перегиба поверхности. Это доказывается и тем, что равенство второй производной от базовой функции нулю, дает уравнение окружности с радиусом а. Поэтому эту величину будем называть радиусом рассеяния.
а б
Рисунок 1
На рисунке 1 представлена криволинейная поверхность при & = 10 (рис. 1а) и & = 20 (рис. 1б). Ордината верхней точки левой поверхности равна примерно шести единицам, а ордината правой поверхности в пять раз меньше. Это объясняется, естественно, тем, что объем фигуры, образованной поверхностью и плоскостью xy не зависит от радиуса рассеяния.
Таким образом, максимальное значение функции при постоянном объеме фигуры зависит от значения & .
при x = a и y = b , f(a,b) = /тз
1
2&2
Или & =
V2^-fm
График зависимости & = p(fmx) приведен на рисунке 2
1
Полученная зависимость в дальнейшем
позволит вычислять максимальное значение функции по значению & и наоборот.
Это свойство функции предлагается использовать для нахождения мест расположения подстанций для питания нагрузок, распределенных по поверхности из условия минимальных затрат на строительство.
Каждую нагрузку на поверхности будем представлять фигурами, представляющими собой тела вращения ограниченные поверхностями, образованными произведением базовой функции на мощность нагрузки. Для потребителя мощностью P и координатами a и b при данном значении & поверхность описывается уравнением
f (x, y) = P
1
2я&2
((x-a)2 +(y~br 2&2
Удельной плотностью мощности будем называть величину, равную сти на единицу поверхности.
Объем тела, ограниченного данной поверхностью и плоскостью
e
среднему xy будет
значению мощночисленно равен
мощности. Ранее указывалось на то, что при увеличении & степень рассеяния мощности по поверхности увеличивается. Другими словами удельная мощность в каждой точке поверхности уменьшается.
В том случае, когда количество потребителей велико, объем вычислений, необходимых для нахождения координат подстанции, увеличивается пропорционально этому количеству. С целью со-
кращения расчетов отдельные группы потребителей можно заменить одним эквивалентным потребителем, мощность которого равна сумме мощностей приемников. Для этого необходимо определить радиус рассеяния, при котором поверхность трех нагрузок будет практически повторять поверхность одновершинной фигуры. Проведен анализ зависимости ординат вершины поверхности, равной по объему сумме трех нагрузок, от значения радиуса рассеяния а и высоты «одновершинной» поверхности, объем под которой равен сумме мощностей нагрузки от значения радиуса рассеяния а .
На рисунке 3 кривая 1 представляет собой зависимость высоты поверхности полученной по предлагаемой выше методике от радиуса рассеяния. Расстояние между местами расположения нагрузок составляет 20 метров. Кривая 2 представляет зависимость высоты эквивалентной поверхности, объем под которой равен сумме мощностей нагрузки. Поведение кривых, представленных на рисунке 3, показывает на то, что при радиусе рассеяния равном примерно 30 метрам координаты вершин практически равны и отличаются друг от друга не более чем на 10 процентов. Из сказанного следует вывод о том, что при радиусе рассеяния, равном примерно полутора максимальным расстояниям между потребителями, совокупность потребителей можно заменить одним эквивалентным, мощностью, равной сумме мощностей группируемых потребителей. Последнее положение приводит к мысли о том, что при работе с реальными потребителями, распределенными по поверхности определенного региона, для эффективного использования метода необходимо провести анализ нагрузок на предмет их плотности распределения с учетом их мощности.
О 100 200 300 400 500 600 700 600 900 1 000
X
Рисунок 4
Определение центра масс по известной до настоящего времени теории заключается практически в вычислении координат центра круга, площадь которого равна сумме площадей кругов, представляющих нагрузки с их координатами [4].
Рассмотрим эффективность нового метода на примере девяти нагрузок. На рисунке 4 нагрузки представлены «холмами», распределенными по поверхности при коэффициенте распределения равном а = 6 . При коэффициенте распределения а = 70 распределение удельной мощности по поверхности меняется.
На рисунке 5 явно выражены 3 вершины эквивалентных нагрузок, координаты которых определяют центры расположения подстанций. В том случае, когда возникает необходимость найти место расположения одной подстанции для питания всех потребителей, рекомендуется выбрать другой коэффициент распределения, зависящий от степени разброса потребителей.
Очевидно то, что представленный способ определения мест расположения подстанций в значительной степени упрощает процедуру вычислений при использовании вычислительной техники при произвольном количестве потребителей. Такой метод может быть использован не только для определения оптимальных мест расположения питающих подстанций, но и для решения других задач, связанных с выбором оптимального расположения объектов на поверхности. Для окончательного вывода об эффективности использования метода необходимо провести анализ экономической эффективности количества подстанций.
Предлагаемый подход к выбору оптимальных с точки зрения затрат на строительство мест расположения подстанций является универсальным. Он может быть использован и для определения мест расположения компенсационных установок для сокращения мощности на передачу электрической энергии. Анализ режимов работы реальных разветвленных электрических сетей показывает то, что большое количество электрической энергии теряется на «передачу» реактивной мощности [5]. Потребителями реактивной мощности являются предприятия с большим количеством станков, в которых используются асинхронные двигатели. К потребителям реактивной мощности относятся и трансформаторы подстанций, работающие в номинальном или недогруженном режиме. Имеется ряд и других потре-
бителей реактивной мощности. Для сокращения потерь, связанных с передачей реактивной мощности, необходимо использовать компенсационные установки. В настоящее время предприятия не заинтересованы в улучшении коэффициента мощности своих потребителей, из-за чего эффективность электроэнергетических сетей снижается. Поэтому единственным решением этой проблемы является установка компенсационных устройств на подстанциях различного уровня. Для того чтобы работа компенсационных установок была эффективной, необходимо на первом этапе определить место установки компенсаторов. И для этого может быть использована предложенная методика. Вторым этапом должна быть разработка автоматических систем регулирования мощности, генерируемой установками.
Однако использование локальных автоматизированных компенсаторов не может преодолеть дефицит генерируемой мощности, так как кроме локальных потребителей реактивной мощности существуют, как было сказано выше, распределенные по сети потребители. Для обеспечения эффективности работы всех компенсирующих устройств сети необходимо разрабатывать систему дистанционного автоматического управления, работа которой основана на анализе распределения реактивных мощностей по сети электроснабжения.
Представленный метод определения центров нагрузок эффективен в случае использования средств вычислительной техники. Метод распределенных или удельных нагрузок является теоретической основой построения информационных и управляющих систем электрических цепей. Наилучший результат он обеспечивает при проектировании новых сетей электроснабжения частично сгруппированных потребителей. Для использования этого метода авторами разработана программа, с помощью которой по исходным данным можно быстро получить координаты подстанций для дальнейшего расчета их мощностей проектирования распределительных сетей.
ЛИТЕРАТУРА
1. Фёдоров А.А. Основы электроснабжения промышленных предприятий /Фёдоров А.А., Каменева В.В. - М.: Энергоатомиздат, 1984 - 466 с.
2. Каменева В.В. Область рассеяния центра электрических нагрузок. - М.:МЭИ, 1971
3. Каменева В.В. К вопросу определения местоположения главных понизительных или распределительных подстанций промышленных предприятий./Каменева В.В., Киреева Э.А., - Электричество,
1972 - 73с.
4. Идельчик В.И. Электрические системы и сети. -М.: Энергоатомиздат, 1989 - 592 с.
5. Железко Ю.С. Компенсация реактивной мощности и повышение качества электроэнергии. -М.: Энергоатомиздат, 1985 г. -224 с.
