Идентификация линейной обобщенной нагрузки в трехфазной трехпроводной сети в задаче распределения фактических вкладов в точке общего присоединения Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.31 ^к 10.20998/2074-272Х.2016.6.11

Ю.Л. Саенко, Д.Н. Калюжный

ИДЕНТИФИКАЦИЯ ЛИНЕЙНОЙ ОБОБЩЕННОЙ НАГРУЗКИ В ТРЕХФАЗНОЙ ТРЕХПРОВОДНОЙ СЕТИ В ЗАДАЧЕ РАСПРЕДЕЛЕНИЯ ФАКТИЧЕСКИХ ВКЛАДОВ В ТОЧКЕ ОБЩЕГО ПРИСОЕДИНЕНИЯ

На основi матричного методу анал1зу електричних кт визначена структура е^валентног схеми замщення узагаль-неного навантаження. Запропоновано споаб наближеного визначення параметрiв е^валентног схеми замщення узагальненого навантаження. Виконано анажз помилок наближеного визначення параметрiв е^валентног схеми замщення узагальненого навантаження стосовно задачi розподту фактичных внесшв у точщ загального приеднання.

Бiбл. 7, табл. 1, рис. 4.

Ключовi слова: еквiвалентна схема замщення лшшного узагальненого навантаження, фактичний внесок, точка загального приеднання.

На основе матричного метода анализа электрических цепей определена структура эквивалентной схемы замещения обобщенной нагрузки. Предложен способ приближенного определения параметров эквивалентной схемы замещения обобщенной нагрузки. Выполнен анализ ошибок приближенного определения параметров эквивалентной схемы замещения обобщенной нагрузки применительно к задаче распределения фактических вкладов в точке общего присоединения. Библ. 7, табл. 1, рис. 4.

Ключевые слова: эквивалентная схема замещения линейной обобщенной нагрузки, фактический вклад, точка общего присоединения.

Введение. Несимметрия и отклонение напряжений являются одними из основных характеристик электрической энергии, которые в большинстве случаев не соответствуют ее качеству [1]. Последствия потребления и передачи электрической энергии пониженного качества общеизвестны [2] и, в конечном счете, сводятся к финансовым убыткам [3]. Определение фактических вкладов в понижение качества электрической энергии в точке общего присоединения (ТОП) является решением задачи о распределении компенсации экономических потерь потерпевшей стороне между всеми участниками рынка электрической энергии.

Постановка задачи. Существующие методы распределения фактических вкладов [4], независимо от выбранного математического пространства своей реализации, имеют одну общую проблему. Она заключается в адекватном представлении и задании параметров схемы замещения обобщенной нагрузки. Так, в пространстве симметричных составляющих обобщенную нагрузку представляют в виде совокупности источника тока или ЭДС со своим сопротивлением, что соответствует схемам замещения Нортона или Тевенена. В пространстве фазных составляющих можно ограничиться пассивными элементами, где однозначного мнения по поводу структуры схемы замещения применительно к задаче распределения фактических вкладов нет.

При несимметричной нагрузке использование метода симметричных составляющих не дает никаких преимуществ по сравнению с пространством фазных составляющих. Исходя из этого, в статье рассматривается вопрос определения структуры и параметров схемы замещения линейной искажающей обобщенной нагрузки в трехфазных трехпроводных сетях в рамках пространства фазных составляющих.

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

фазных составляющих для задачи распределения фактических вкладов в токе общего присоединения.

Результаты исследования. Рассмотрим систему электроснабжения (СЭС) с выделенной обобщенной нагрузкой относительно ТОП, в которой производятся измерения показателей качества электрической энергии (рис. 1,а). Будем считать, что все элементы электроэнергетической системы и приемники электрической энергии обобщенной нагрузки могут вызывать только линейные (синусоидальные) искажения напряжения. Тогда схему замещения такой СЭС можно представить в виде соединения активных и пассивных элементов (рис. 1,6). Дополнительно примем допущение о пренебрежении взаимными сопротивлениями и проводимостями отдельных элементов СЭС.

Обобщенная нагрузка

Рис. 1. Система электроснабжения: а - структурная схема относительно ТОП; 6 - условная схема замещения электроэнергетической системы и обобщенной нагрузки

© Ю.Л. Саенко, Д.Н. Калюжный

Y * аа Yab х Фa J a

Yba Ybb Фь J b

Согласно методу узловых потенциалов общая схема замещения электроэнергетической системы и обобщенной нагрузки описывается следующим матричным уравнением:

(1)

где а - порядковый номер последнего узла со стороны схемы замещения электроэнергетической системы; Ь - количество независимых узлов в общей схеме замещения электроэнергетической системы и обобщенной нагрузки.

Решение уравнения (1) относительно подматрицы Ф а с учетом отсутствия активных элементов в схеме замещения обобщенной нагрузки имеет следующий вид:

(Уаа + Уед ) Ф а = У а , (2)

где Уед - матрица узловых проводимостей эквивалентной обобщенной нагрузки, которая характеризует структуру ее схемы замещения.

Матрица Уед определяется следующим образом:

Yeq - Yab х Ybb х Yba •

YBC &

Juab • yab + 0 • YBC + UAC • yac - -A ; [UBA • YAB + UBC • YBC + 0 • YAC - —B •

(4)

( Yab , Ybc и Yac ), является неопределенной и совместной. Ее общее решение имеет следующий вид:

■Ía -Uac • Yac)/u

[ Y-BC -~{¡-C + UAC • YAC УUBC-

Yab -i

-AB ■

(5)

Y

Для определения неизвестных проводимостей .АВ, ВС и ас можно прибегнуть к переопределению системы (5) путем использования данных двух независимых, например, следующих друг за другом измерений. Применяя метод наименьших квадратов к такой системе уравнений, неизвестные проводимости будут определяться следующим образом:

Y - U

UT и г

U T1,

(6)

(3)

Особенностью матрицы Уед в рассматриваемом

нами случае является то, что в ней все элементы, за исключением элементов блока размерностью [3x3], нулевые. Из этого следует, что эквивалентная схема замещения обобщенной нагрузки должна иметь три независимых узла. Этому условию соответствует структура схемы замещения, изображенная на рис. 2.

Рис. 2. Структура эквивалентной схемы замещения обобщенной нагрузки

Рассмотрим задачу определения параметров схемы замещения обобщенной нагрузки. Исходными данными для этого являются измеряемые в ТОП фазные токи (I а , 1_в и 1_с ) и линейные напряжения

( иав , ивс и иас ).

Параметры схем замещения и измеряемые параметры режима связаны между собой следующим образом:

где и - матрица измеренных напряжений; I - матри-

^

ца-столбец измеренных токов; У - матрица-столбец искомых проводимостей.

Использование изложенного подхода определения неизвестных проводимостей требует соблюдения следующих условий. Измеренные параметры режима работы должны отображать два различных состояния СЭС, вызванные изменениями только со стороны электроэнергетической системы. С практической точки зрения такая ситуация в СЭС маловероятна, поэтому многие исследователи, например [5, 6], предлагают проводить активный эксперимент на основе принудительного коммутирования батареи конденсаторов в ТОП. Учитывая методику измерения ПКЭ, где минимальный интервал измерения параметров режима работы начинается с 0,2 с, этот подход оказывается неприемлемым. Кроме этого, результаты математического моделирования показывают высокую чувствительность решения (6) к незначительным изменениям состояния обобщенной нагрузки, а также величине возмущающего воздействия со стороны электроэнергетической системы. Как следствие, адекватное определение параметров схемы замещения обобщенной нагрузки согласно выражению (6) в условиях реальной эксплуатации является технически не реализуемым.

Рассмотрим возможность приближенного определения параметров схемы замещения обобщенной нагрузки. С этой целью обратимся к ее схеме замещения, представленной в виде соединения элементов в звезду (рис. 3).

Система (4), состоящая из двух алгебраических линейных уравнений и содержащая три неизвестных

Рис. 3. Эквивалентная схема замещения обобщенной нагрузки с соединением элементов в звезду

в

Рассматриваемая схема замещения с четырьмя узлами описывается следующей системой трех уравнений:

Ua ■ Ya = La ; Ub ■ Yb = Lb ;

(7)

U,

■ Yc = Lc •

Как известно, точное измерение фазных напряжений в трехфазных трехпроводных сетях, из-за физического отсутствия нейтральной точки, в общем случае выполнить невозможно. Поэтому для приближенного определения фазных напряжений воспользуемся измеряемыми значениями линейных напряжений:

UA =

UAB1 ■

, - у30

+ kAB2 ■е

j30

)/ V3;

UB = Uba + UA; UC = Ucb + UB ,

(8)

где Ц.АВ1 и Ц.АВ2 - линейные напряжения прямой и

обратной последовательностей соответственно.

Ошибка определения фазных напряжений для схемы замещения рис. 3 будет определяться разно-

стью

точного

(u A + U В + UC )/3

(рис. 4):

N

и приближенного

< = <n -(u

A + UB + kC

)/3-

(9)

_ ив

Рис. 4. Точные и приближенные значения фазных напряжений для эквивалентной схемы замещения обобщенной нагрузки с соединением элементов в звезду (рис. 3)

Таким образом, приближенные значения прово-димостей УА, УВ и УС , вычисленные согласно (7) с учетом (8), связаны со своими точными значениями следующим равенством:

' (10)

Yph = Yph - Yph■<k

ph '

где символ «ph» соответствует фазе А, В или С.

Для решения задачи распределения ФВ основополагающим является разделение параметров схемы замещения на неискажающие и искажающие части [7]. При этом наиболее чувствительными к данному разделению, с точки зрения ошибок определения фазных проводимостей, оказываются неискажающие (nondisturbing) и искажающие (disturbing) части,

используемые для распределения ФВ по несимметрии напряжений:

nondis ph

:Ь + Yb + Yc У3

Ydis = Y - Ynondis. —ph ~ —ph —ph •

(11)

Исходя из этого, проведем анализ ошибок определения неискажающих и искажающих частей параметров эквивалентной схемы замещения обобщенной нагрузки для двух характерных случаев, когда на шинах сбалансированной и несбалансированной нагрузок имеет место несимметричное напряжение.

При сбалансированной нагрузке с несимметричным напряжением ошибки определения У

nondis ph

Y% составляют:

SY

nondis = Yph

(

ph

= Yph Uph-

\

1 1 1

—+—+—

v kA UB kA;

Yph

(

3

1 1 1

kA + йВВ+йс

(12)

• (13)

значений потенциалов точки N

Упростим множитель

(k А2 У k

положении, что

(

а1

( \

1 1 1

U А + U В + UC

<< 1:

в пред-

1 1 1

-+-+-

V k А1 + Ша2 UB kA;

uа2 _+ ub2 +_ uc2

(14)

(12 (B112 (C 1}) ()2

3 ■и А2

Тогда с учетом (14) приближенное значение

„v nondis

sYph составит:

SY

nondis ph

U

а2

(15)

-ph V •

( A1i

Принимая во внимание, что множитель

k

llkА1J

'-АН илм на два порядка уменьшает результат произведения У^' лр в (15), величину погрешности

SY

nondis ph

можно считать стремящейся к нулю. Так,

например, при \л<т равном 10 % от модуля фазного

значения прямой последовательности (

U А1

) и

U А2

равном 4 % от

k А1J

ошибка определения YpTd's

сбалансированной нагрузки будет иметь порядок - 0,004 • , что в относительных единицах соответствует -0,4 %.

Ошибка ЗУУ^ с учетом (14) будет приближенно определяться следующим образом:

3

sïdL - Yph

( \ 1 + ua2

(16)

Uph Ai J

Учитывая, что UA2/(UAiF величина второго

1 Uph

порядка малости относительно записать:

лр

5Lph ~Lph ■

Uph

Как видно из (17), значение ошибки зУр* будет определяться величиной лр. Так, например, при

\лр\ = 0,1

U

ph

порядок ошибки sYph составит

ph

idis

0,1 ^ри при истинном значении = 0.

При несбалансированной нагрузке с несимметричным напряжением имеем:

(

„v nondis

sLph

лр

U

■ ej(a~pA J + YB ■

лр

U

eJ(-Рв ) + Y

лр

U

eJ'('

а-рс )

(18)

sLpS = Lph ■л\Uph -1 ■

f

1a лр

—=-

V U A

■ eJ(a-\A ).

+ h

л

лр ■ e^B ) + yc лр ■ ^'(арJ

иB UC /

(19)

где « = ащ{лр); РрН = ш%Црк).

Допустим, что значения фазных напряжений (ЦА, Ц в и ЦС ) приближенно можно считать симметричными. Тогда с учетом возможности произвольного поворота системы векторов имеем:

Y.

лр

U

eJ0 + Y

i-B

а

лр

лр

U-

e J120

+Y

лр

U

QJ 240

U

■ ■ kî ■ eJ0 + if ■ ej120 + YCs ■ eJ240

(20)

ph\

частей проводимостей (lAs , YdBS и yC ) таковы, что

eJ 240 в сумме дадут вектор

модуль ошибки зУ"0"Р'* для случая несбалансированной нагрузки может быть оценен следующим не-

dis

Y

dis

Y

dis

то

равенством:

„v nondis

sLph

\лр

U

ph

Y dis —ph max

(21)

где Ур* тах - максимальная фазная искажающая проводимость, которая для крайнего случая (несимметричная нагрузка, соединенная по схеме звезда (рис. 3,) задана одним плечом) равна:

Y

dis

=2 Y.

можем

(17)

-ph max з — ph'

Таким образом, для \л<т = 0,1

Uph

(22)

имеем

SY

nondis ph

< 0,067

Y

ph

dis

Рассмотрим ошибку зУр* с учетом того, что второе слагаемое в (19) представляет собой зУ"°"'* (21):

SY

dis

лр

■1

— si nondis ) ph -S1ph

(23)

-рк ~ и \^рк —рк Г Црк

С учетом (22) модуль максимального значения ошибки зУр'* может быть оценен таким неравенст-

вом:

SY

dis

ph max

\лср

Uph

Y

Количественно < 0,17

при \л(р

ph ' = 0,1

(24)

U

ph

Si dis —ph max

Lph

Из проведенного анализа ошибок si

nondis ph

dis ph

sY ph следует, что неискажающие части проводимо-

Если предположить, что аргументы искажающих

dis vdis vdis

A , LB и Le

j 0 /120

их умножение на единичные векторы e , e и

стей сбалансированной и несбалансированной обобщенной нагрузки могут быть определены с достаточно высокой точностью. Определение искажающих частей характеризуется значительными ошибками, соизмеримыми с точными значениями искомых величин.

Выводы. Структура эквивалентной схемы замещения линейной обобщенной нагрузки в трехфазных трехпроводных сетях определяется тремя независимыми узлами, что соответствует схеме соединения элементов в треугольник. Точное определение параметров этой схемы замещения на основе методики измерения показателей качества электрической энергии является технически не реализуемым. Предложен способ приближенного определения параметров эквивалентной схемы замещения линейной обобщенной нагрузки, основанный на приближенном определении фазных напряжений. Ошибки вычисления неиска-жающих частей параметров эквивалентной схемы замещения линейной обобщенной нагрузки характеризуются незначительными величинами, которыми можно пренебречь. Ошибки вычисления искажающих частей могут достигать величин, соизмеримых с точными значениями.

СПИСОК ЛИТЕРАТУРЫ 1. Kartashev I., Tulsky V., Shamonov R. Assessing of impact of disturbing load on power quality // 17th International

x

<

и

+

+

<

Conference on Electricity Distribution. Barcelona, 12-15 May, 2003. pp. 1-5.

2. Жежеленко И.В., Шидловский А.К., Пивняк Г.Г., Саен-ко Ю.Л., Нойбергер Н.А. Электромагнитная совместимость потребителей. - М.: Машиностроение, 2012. - 349 с.

3. Чэпмэн Д. Цена низкого качества электроэнергии // Энергосбережение. - 2004. - №1. - С. 66-69.

4. Sayenko Yu., Kalyuzhniy D. Analytical methods for determination of the factual contributions impact of the objects connected to power system on the distortion of symmetry and sinusoidal waveform of voltages. Przeglad Elektrotechniczny. -2015. - vol.11. - pp. 81-85. doi: 10.15199/48.2015.11.23.

5. Гамазин С.И., Петрович В.А. Определение фактического вклада потребителя в искажении параметров качества электрической энергии // Промышленная энергетика. - 2003. -№1. - С. 32-38.

6. Aiqiang Pan, Jian Zhou. Power quality analysis and harmonic tracing in city grid based on big monitoring data // 23rd International Conference on Electricity Distribution. Lyon, 1518 June 2015. pp. 1-4.

7. Sayenko Y., Sukhonos M., Kalyuzhniy D., Bolgov V. Mathematical model for real-time assessment of contributions of disturbing sources to power quality level at a Point of Common Coupling // 10th Electric Power Quality and Supply Reliability (PQ) Conference. Tallinn, 29-31 August, 2016, pp. 29-35. doi: 10.1109/pq.2016.7724085.

REFERENCES

1. Kartashev I., Tulsky V., Shamonov R. Assessing of impact of disturbing load on power quality. 17th International Conference on Electricity Distribution. Barcelona, 12-15 May, 2003. pp. 1-5.

2. Zhezhelenko I.V., Shidlovskiy A.K., Pivnyak G.G., Sayenko Yu.L., Noyberger N.A. Elektromagnitnaya sovmestimost potre-biteley [Electromagnetic compatibility of customers]. Moscow, Mashinostroenie Publ., 2012. 349 p. (Rus).

3. Chepmen D. Price of low power quality. Energosberezhenie - Energy Saving, 2004, no.1, pp. 66-69. (Rus).

4. Sayenko Yu., Kalyuzhniy D. Analytical methods for determination of the factual contributions impact of the objects connected to power system on the distortion of symmetry and sinusoidal waveform of voltages. Przeglad Elektrotechniczny, 2015, vol.11, pp. 81-85. doi: 10.15199/48.2015.11.23.

5. Gamazin S.I., Petrovich V.A. Determining the actual contribution of the consumer in the distortion of the quality parameters of electric energy. Promyishlennaya energetika, 2003, no.1, pp. 38-42. (Rus).

6. Aiqiang Pan, Jian Zhou. Power quality analysis and harmonic tracing in city grid based on big monitoring data. 23rd International Conference on Electricity Distribution. Lyon, 1518 June 2015. pp. 1-4.

7. Sayenko Y., Sukhonos M., Kalyuzhniy D., Bolgov V. Mathematical model for real-time assessment of contributions of disturbing sources to power quality level at a Point of Common Coupling. 10th Electric Power Quality and Supply Reliability (PQ) Conference, Tallinn, 29-31 August, 2016, pp. 29-35. doi: 10.1109/pq.2016.7724085.

Поступила (received) 29.09.2016

Саенко Юрий Леонидович1, д.т.н., проф., Калюжный Дмитрий Николаевич2, к.т.н., доц.,

1 Приазовский государственный технический университет, 87500, Мариуполь, ул. Университетская, 7,

тел/phone +38 0629 446551, e-mail: [email protected]

2 Харьковский национальный университет городского хозяйства им. А.Н. Бекетова, 61002, Харьков, ул. Революции, 12,

тел/phone +38 050 5606835, e-mail: [email protected]

Yu.L. Sayenko1, D.N. Kalyuzhniy2

1 Pryazovskyi State Technical University,

7, Universytets'ka Str., Mariupol, 87500, Ukraine.

2 O.M. Beketov National University of Urban Economy in Kharkiv,

12, Revolution Str., Kharkiv, 61002, Ukraine. Identification of generalized linear load in three-phase three-wire network in the problem of the distribution of actual contributions at the point of common coupling. Purpose. To perform structural and parametric identification of generalized load equivalent circuit of three-phase three-wire load in the network in the space of phase components. Methodology. Underlying structural identification methods are matrix analysis of electrical circuits. Parametric identification is based on the basic laws of electrical engineering. Results. The structure of a generalized load equivalent circuit is composed in three independent nodes. An approximate method for determining its parameters is proposed. The estimation error determination undistorted and distorted parts of the parameters of generalized load equivalent circuit. Originality. Approximate determination of equivalent circuit parameters are based on the results of a single measurement of voltages and phase currents. Practical value. The proposed replacement structure and a method for determining its parameters of the circuit can be used in the problem of the distribution of actual contributions at the point of common coupling. References 7, tables 1, figures 4. Key words: equivalent circuit of linear generalized load, actual contribution, point of common coupling.