CC BY
72
Оригинальная статья / Original article УДК: 621.311
DOI: 10.21285/1814-3520-2016-7-141-150
ОПРЕДЕЛЕНИЕ МЕСТ КОРОТКИХ ЗАМЫКАНИЙ НА ДЛИННЫХ ЛИНИЯХ ЭЛЕКТРОПЕРЕДАЧИ С УЧЕТОМ РАСПРЕДЕЛЕННОГО ХАРАКТЕРА ПАРАМЕТРОВ ЛИНИЙ
© С.Н. Маточкин1, С.Г. Тигунцев2, А.Б. Осак3
12
' Иркутский национальный исследовательский технический университет, 664074, Россия, г. Иркутск, ул. Лермонтова, 83. 3Институт систем энергетики им. Л.А. Мелентьева СО РАН, 664033' Россия, г. Иркутск, ул. Лермонтова, 130.
Резюме. Цель. Одним из путей повышения точности при определении места короткого замыкания (ОМКЗ) на длинных (магистральных и сверхдальних) воздушных линиях электропередачи высокого и сверхвысокого напряжения (220-750 кВ) является учет распределенного характера параметров линии, поэтому целью исследования стало создание эффективной методики определения места короткого замыкания. Материал и методы. Разработана методика определения места короткого замыкания на длинных трехфазных воздушных линиях электропередачи. В основе методики, учитывающей распределенный характер параметров системы, использовано численное интегрирование для телеграфных уравнений, моделирующих возникающие электромагнитные процессы. Место короткого замыкания определяется путем построения графиков изменения напряжения вдоль линии с двух ее концов. Результаты. Проверка методики проведена на трехфазной модели линии с использованием программы расчета режимов в фазных координатах. Результаты представлены на графиках и сведены в таблицу, на их основе сделан вывод о высокой точности метода. Заключение. Методика позволяет определять место всех видов коротких замыканий для любых напряжений, при этом она не чувствительна к переходному сопротивлению. Точность ОМКЗ зависит только от точности измерения токов и напряжений, а также от точности задания параметров схемы замещения линии.
Ключевые слова: воздушные линии электропередачи, электрические сети, определение мест коротких замыканий на воздушных линиях, дифференциальные уравнения для линий электропередачи, короткое замыкание, распределенные параметры.
Формат цитирования: Маточкин С.Н., Тигунцев С.Г., Осак А.Б. Определение мест коротких замыканий на длинных линиях электропередачи с учетом распределенного характера параметров линий // Вестник Иркутского государственного технического университета. 2016. № 7. С. 141-150. DOI: 10.21285/1814-3520-2016-7-141-150
FAULT LOCATION ON LONG POWER TRANSMISSION LINES WITH ALLOWANCE FOR THE DISTRIBUTED
NATURE OF LINE PARAMETERS
S.N. Matochkin, S.G. Tiguntsev, A.B. Osak
Irkutsk National Research Technical University, 83, Lermontov St., Irkutsk, 664074, Russia. Melentiev Energy Systems Institute SB RAS, 130, Lermontov St., Irkutsk, 664033, Russia.
Abstract. Purpose. One of the methods improving fault location accuracy on long (trunk and ultra-long) overhead transmission lines of high and extra high voltage (220-750 kV) is the account of the distributed nature of line parameters. Therefore, the purpose of the study is the creation of an effective procedure for fault location determination. Material and methods. A procedure is designed for fault location on long three-phase overhead power transmission lines. Considering distributed nature of system parameters it is based on the use of numerical integration for telegraph equations describing electromagnetic processes that arise in the line. The fault location is determined by plotting a voltage curve along the transmission line from its two ends. Results. The procedure has been tested on a three-phase model of the transmission line using a program for calculating modes in phase coordinates. The results are shown in the graphs and summarized in
1
Маточкин Сергей Николаевич, аспирант, e-mail: [email protected] Matochkin Sergey, Postgraduate, e-mail: [email protected]
2Тигунцев Степан Георгиевич, кандидат технических наук, доцент кафедры электрических станций, сетей и систем, e-mail: [email protected]
Tiguntsev Stepan, Candidate of Engineering, Associate Professor of the Department of Electrical Stations,
Networks and Systems, e-mail: [email protected]
3Осак Алексей Борисович, научный сотрудник, e-mail: [email protected]
Osak Aleksei, Researcher, e-mail: [email protected]
tables. The conclusion on the high accuracy of the method is made on their basis. Conclusion. The procedure enables determination of the location of all types of short circuits, for any voltages. Moreover, it is not sensitive to transition resistance. The accuracy of fault location determination (FLD) depends on the accuracy of current and voltage measurements, as well as on the accuracy of transmission line equivalent circuit parameter setting.
Keywords: overhead power transmission lines, electrical networks, fault location on overhead lines, differential equations for power transmission lines, short circuit, distributed parameters
For citation: Matochkin S.N., Tiguntsev S.G., Osak A.B. Fault location on long power transmission lines with allowance for the distributed nature of line parameters. Proceedings of Irkutsk State Technical University. 2016, no. 7, pp. 141-150 (in Russian). DOI: 10.21285/1814-3520-2016-7-141-150
Введение
Одним из путей повышения точности при определении места короткого замыкания (ОМКЗ) на длинных (магистральных и сверхдальних) воздушных линях электропередачи высокого и сверхвысокого напряжения (220-750 кВ) является учет распределенного характера параметров линии.
В общем случае трехфазные линии характеризуются наличием взаимоиндукции между проводами и емкости между проводами и землей (рис. 1).
U'a U'b U'c
U''a U''b U''c
= _
E''a
Рис. 1. Схема замещения трехфазной воздушной линии электропередачи Fig. 1. Equivalent circuit of a three-phase overhead electricity transmission line
Линия в трехфазном виде с учетом различия продольных и поперечных параметров наиболее полно описывается дифференциальными уравнениями:
du di m = n di
■—t = ri + L — + I M —mm dx kk k dt m = l km dt
di m = n du m = n d(u - U )
■— = g u + I g (u - u ) + C + I C —k-m
dx k k km k m k dt Л km dt
(1)
m = 1
где k = a, Ь, c - обозначения проводов фаз; Lk, gk, Ck - собственные параметры проводов фаз на единицу длины с учетом влияния земли; Mkm и С^ - взаимная индуктивность и емкость между ^м и m-м проводами на единицу длины с учетом влияния земли.
Наличие частных производных в уравнениях обусловлено тем, что напряжение u(x, ^ и ток i(x, ^ зависят от двух аргументов. Если в линии существует режим гармонических колебаний, то можно применять символический метод анализа. Так как комплексные амплитуды напряжения и тока не зависят от времени, то вид телеграфных уравнений существенно упрощается:
dU dx
dU
dx dU
a = 7 I + 7 I + 7 I
7aIa + 7abIb + 7ac1c
= 7 b Ib + 7 bc Ie + Zab I a
c = 7CIC + 7bcIb + 7Ja
dx
- ^ = YUa + YabUab + YJJac dx
- ^ = YU + + Y.U
dx
- ^ = YU + + yjut dx
(2)
h h ah ab bc bc
c c ac ac
bc bc
Решение системы дифференциальных уравнений (2) для трехфазной линии электропередачи в соответствии с теорией волновых процессов - задача громоздкая и для практики малоприменимая. Аналитического решения этих уравнений для целей ОМКЗ авторам не известно.
В тоже время телеграфные уравнения, полученные для однофазной линии, можно применить по отдельности к прямой, обратной и нулевой последовательностям схемы замещения трехфазной линии электропередачи, но при этом вносится методическая погрешность в ОМКЗ, обусловленная усреднением параметров линии. С учетом этого длинные трехфазные линии обычно принимают симметричными и описывают дифференциальными (телеграфными) уравнениями для одной фазы в симметричных составляющих. Для таких уравнений имеется аналитическое решение. Известны методы ОМКЗ на длинных линиях, представленных в схемах замещения симметричных составляющих, которые используют аналитические решения дифференциальных уравнений [1-3].
Создание эффективной методики определения места короткого замыкания
Уравнения (2) можно решить методом численного интегрирования [4]. С учетом этого разработан алгоритм расчетного ОМКЗ, в котором по измеренным при коротком замыкании фазным токам и напряжениям по концам линии, по известным погонным продольным и поперечным параметрам трехфазной линии определяется место КЗ путем построения графиков пошагового изменения напряжения от концов линии к месту короткого замыкания. При этом фазные напряжения в конце каждого шага (участка линии, длина которого может быть любой, вплоть до пролета между опорами) определяются как разница между напряжением в начале участка и падением напряжения на участке, который представлен трехфазной П-образной схемой замещения, а фазные токи в конце участка определяются как разница токов в начале участка и токов в поперечных проводимостях линии.
В предлагаемом методе предварительно формируется модель линии как значения продольных и поперечных параметров N участков П-образной схемы замещения линии в трехфазном виде:
7 = ij
7 7AA ij 7 7AB ij 7 7AC ij Y iaa ij Y aAB IJ Y aAC IJ
7 7BA ij 7 7BB ij 7 7BC ij ; Y = IJ Y aAB IJ Y ABB ij Y ABC IJ
7 7ca ij 7 7CB ij 7 7CC ij Y CA ij Y ACB ij Y YCC ij
где 2ААу,2ВВу,2ССу - значения собственных продольных сопротивлений фаз участка Н линии,
Ом; £^у,2АСц,2ВАу,2ВСу,2САц,2СВц - значения взаимных продольных сопротивлений фаз участка Н линии, Ом; - значения собственных поперечных емкостных прово-
димостей фаз половины участка Н линии, Сим; - значения
взаимных поперечных емкостных проводимостей фаз половины участка Н линии, Сим (рис. 2).
Рис. 2. Участок П-образной схемы замещения воздушной линии в трехфазном виде Fig. 2. A section of a U-shaped three-phase equivalent circuit of an overhead line
Значения собственных и взаимных сопротивлений емкостных проводимостей фаз на «землю» и взаимных емкостных проводимостей между проводами фаз определяются по общеизвестным выражениям [5].
Далее, после получения измеренных при коротком замыкании значений фазных напряжений на шинах и токов с двух концов линии (' и "), задаются поочередно точки j в конце каждого участка вдоль линии, рассчитываются и сохраняются для двух концов линии значения комплексных фазных напряжений в каждой j-й точке по выражениям:
! ! . ! ! . |
Ù =(1-Y *Z )"1*(1 + Y*Z MJ.-Z *i ; j ij г ij ij 1 ij ij
. !! !! !! , !! !! !! !! . || U = (1 - Y *Z ) *(1 + Y *Z VU.-Z *I ,
j У Г У Г i У D
ûAl . !! û;
где Ù. = Ubi ÛCI , U. = UBi u;
i i
- значения комплексных фазных напряжений в каждой /-й точке ли-
нии с одного и другого конца линии, для / = 1 значения напряжений, измеренных на шинах од- значения комплексных фазных напряжений в каж-
ного и другого конца линии, В; U. =
j
Bj
Cj
дой у-й точке линии с одного и другого конца линии, В; I.. =
LAIJ
LBij
Cij
- значения комплексных фаз-
ных токов на участке /-/ с одного и другого конца линии, для / = 1; значения комплексных фаз-
ных токов, измеренных с одного и другого конца линии, А; Ъ =
у
7 7 7
7AA ij 7AB ij 7AC ij
7 7 7
7BA ij 7BB ij 7BC ij
7 7 7
7CA ij 7CB ij 7CC ij
H
7 = ij
7 7 7
7AAIJ 7ABIJ 7ACIJ
7 7 7
7BAIJ 7BBij 7BCij
7 7 7
7CAij 7CBij 7CCij
- значения продольных собственных и взаимных сопротивлений
участков /-/ схемы замещения линии с одного и другого конца линии, Ом.
Рассчитываются значения фазных токов в поперечных емкостных проводимостях в /-й и /-й точках участка линии по выражениям:
.! ! . ! .! I .1 .11 !!.!!.!! !! . !!
I = У *и ; I = У *и ; I = У *и ; I = У .
10 1] 1 ,)0 1] ] 10 1] 1 ,]0 1] ]
Если в /-м узле включена отпайка, то рассчитываются значения фазных токов в поперечных емкостных проводимостях и в отпайке в /-й точке линии по выражениям:
I =(y'+Y )*Ù ; I = (Y + Y )*U .
J0 ij J отп j J0 ij J отп j
Рассчитываются и сохраняются значения фазных токов в продольных сопротивлениях в каждом (//+1)-ом участке линии по выражениям:
.! . ! .! .! .!! .!!.!! . !!
I =1 -I -I ; I =1 -I -I ,
ij + 1 ij Ю jo ij + 1 ij iO jo
где У. =
¡J
y' yAA ij y' yAB ij y' yAC ij tt y" yaa ij y'' yAB ij y'' yAC ij
y' YAB ij y' 1BB ij y' ybc ij , y = ij Y yAB ij yBB ij YBC ij
y' YCA ij y' yCB ij y' yCC ij y" yCA ij YCB ij YCC ij
- значения поперечных собственных и
взаимных емкостных проводимостей половины участка /-/ схемы замещения линии с одного и
другого конца линии, Сим; I =
Ai
LBi
Ci
I =
i0
Ai
LBi
Ci
- значения рассчитанных фазных токов в попе-
речных емкостных проводимостях в начале каждого //-го участка линии с одного и другого кон-
ца линии, А; I =
Aj
Bj
, I = jo
Aj
Bj
- значения рассчитанных фазных токов в поперечных емкост-
ных проводимостях в конце каждого //-го участка линии с одного и другого конца линии, А;
41
41
Y
отп J
Y A отп J 0 0
0 Y B отп J 0
0 0 Y Y отп
- фазные значения проводимостей отпайки, включающие в
себя проводимости линии и трансформатора от отпайки до нагрузки и нагрузки отпайки, Сим.
Далее из сохраненных значений комплексных фазных напряжений и и и выделяются модули, по которым строятся графики с двумя осями зависимости модулей напряжений от номера участка (от расстояния). Точка пересечения графиков соответствует точке короткого замыкания.
Обсуждение результатов
Таким образом, представленный метод позволяет определять место всех видов коротких замыканий для линий любых напряжений, при этом он не чувствителен к переходному сопротивлению КЗ.
Линия может быть задана как однородная, т.е. параметры участков вдоль всей длины не изменяются, так и разнородная, т.е. учитываются особенности каждого участка линии, что позволяет с высокой точностью учитывать транспозицию линии и отпайки на линии.
Линия может быть как одноцепная, так и многоцепная, причем цепи могут быть разных напряжений.
В предлагаемом методе не требуется выполнять синхронизацию измеренных величин токов и напряжений по концам линии.
Проверка метода выполнена на условной трехфазной модели линии напряжением 500 кВ, длиной 300 км, представленной 300-ми участками П-образной схемы замещения (см. рис. 2). Расчеты режимов с КЗ проводились в программе расчета электрических режимов в фазных координатах, созданной на кафедре электрических станций, сетей и систем Иркутского национального исследовательского технического университета [6].
Математическим описанием установившихся режимов (УР) электрической системы являются уравнения состояния линейной электрической цепи на основе законов Ома и Кирхгофа. Количество уравнений состояния таких систем велико, поэтому для составления системы используется формализованный подход представления конфигурации схемы замещения с помощью элементов теории графов и алгебры матриц. Параметры элементов сети при формализованном подходе отображаются диагональной матрицей сопротивления ветвей:
Zв = diag(Zl), i = 1,...,Ш.
Взаимная индуктивная, а также магнитная связь между элементами учитывается включением в матрицу недиагональных элементов.
Далее, используя формализованный подход, для расчета УР в фазных координатах формируют первую М и вторую N матрицы инциденций.
Матрица М характеризует соединения ветвей в узлах, число строк которой равно числу вершин графа (узлов), а число столбцов - числу ребер (ветвей):
М = (т), i = 1..п; j = 1...т;
ту = 1, если узел / является начальной вершиной ветви у; ту = -1, если узел / является конечной вершиной ветви у; ту = 0, если узел / не является вершиной ветви у.
Матрица N соединяет ветви в независимые контуры; число строк в ней равно числу независимых контуров к, а число столбцов - числу ветвей т:
N = Ц), ! = 1...к; ) = 1...Ш,
п/ = 1, если ветвь входит в контур и их направления совпадают; п/ = -1, если ветвь входит в контур и их направления противоположны; п/ = 0, если ветвь не входит в контур.
Каждая строка матрицы N показывает, какие ветви входят в контур, каждый столбец - в состав каких контуров входит ветвь.
Далее определяют матрицу контурных сопротивлений:
матрицу контурных проводимостей:
ZK = N*ZB*NT ;
Y = Z
-1.
матрицу контурных токов:
матрицу токов в ветвях:
IK = Z-1 *E ;
к к к'
Ib = n*Ik ;
матрицу падений напряжения в ветвях:
dUB = Zb*IB - EB ;
матрицу падений напряжения в ветвях дерева графа относительно балансирующего узла:
dUa = Zb3 *Ib3 - Eb
матрицу напряжений в узлах схемы:
Uy = (МТ )-1*dU
a-
В расчете на модели задавались точки короткого замыкания на разных расстояниях от концов линии с разными переходными сопротивлениями 1п. Результаты некоторых из них представлены на графиках (рис. 3-6).
Рис. 3. Zn = 0, короткое замыкание на 270 км Fig. 3. Zn = 0, short circuit at the point of 270 km
Рис. 4. Zn = 0, короткое замыкание на 75 км Fig. 4. Zn = 0, short circuit at the point of 75 km
Рис. 5. Zn = 10, короткое замыкание на 225 км Fig. 5. Zn = 10, short circuit at the point of 225 km
Рис. 6. Zn = 30, короткое замыкание на 30 км Fig. 6. Zn = 30, short circuit at the point of 30 km
Полные результаты для различных смоделированных коротких замыканий сведены в
табл.1.
шя Энергетика
M Power Engineering
Таблица 1
Результаты определения места короткого замыкания для различных смоделированных точек КЗ на модели линии
Table 1
Results of fault location determination (FLD) for various simulated points of short circuits
on the line model
Zn = 0 Ом / Zn = 0 Ohm
Место КЗ / Fault location Результаты ОМКЗ / FLD Results Погрешность / Error
Относительное расстояние / Relative distance км/km км/km км/km %
0,1 30 30 0 0
0,25 75 75 0 0
0,5 150 150 0 0
0,75 225 225 0 0
0,9 270 270 0 0
Zn = 10 Ом / Zn = 10 Ohm
0,1 30 30 0 0
0,25 75 75 0 0
0,5 150 150 0 0
0,75 225 225 0 0
0,9 270 270 0 0
Zn = 30 Ом / Zn = 30 Ohm
0,1 30 30 0 0
0,25 75 75 0 0
0,5 150 150 0 0
0,75 225 225 0 0
0,9 270 270 0 0
Заключение
Представленная методика позволяет определять место всех видов коротких замыканий для любых напряжений, при этом она не чувствительна к переходному сопротивлению. Точность ОМКЗ зависит только от точности измерения токов и напряжений, а также от точности задания параметров схемы замещения линии.
Библиографический список
1. Висящев А.Н. Приборы и методы определения места повреждения на линиях электропередачи: в 2 ч. Иркутск: Изд-во ИрГТУ, 2001.Ч. 1. 188 с.
2. Висящев А.Н., Устинов А.А. Определение места повреждения на линиях электропередачи с учетом волновых процессов // Повышение эффективности производства и использования энергии в условиях Сибири: материалы Всерос. науч.-практ. конф. Иркутск: Изд-во ИрГТУ, 2006. C. 311-322.
3. Шалыт Г.М. Определение мест повреждений в электрических сетях. М.: Энергоиздат, 1982. 312 с.
4. Вершинин Е.А., Трофимов С.Н. Два метода исследования проводных структур: материалы VI Всерос. конф. молодых ученых по математическому моделированию и информационным технологиям (Кемерово, 29-31 октября 2005) [Электронный ресурс]. URL: http://www.ict.nsc.ru/ws/YM2005/9472/ (28.03.2016).
5. Ульянов С.А. Электромагнитные переходные процессы в электрических системах. М.: Энергия, 1970г. 519 с.
6. Осак А.Б. Практические методы расчета токов короткого замыкания в ЭЭС при сложных видах замыканий с использованием фазных координат // Технические науки - от теории к практике. 2016. № 2 (50). С. 181-192.
References
1. Visyashchev A.N. Pribory i metody opredeleniya mesta povrezhdeniya na liniyakh elektroperedachi [Devices and methods for fault location determination on power transmission lines]. Irkutsk, IrGTU Publ., 2001, 188 p. (in Russian).
2. Visyashchev A.N., Ustinov A.A. Opredelenite mesta povrezhdeniya na liniyakh elektroperedachi s uchetom volnovikh processov [Wave process-wise determination of fault locations on power transmission lines]. Materialy Vserossiiskoi nauchno-prakticheskoi konferentsii "Povyshenie effektivnosti proizvodstva i ispol'zovaniya energii v usloviyakh Sibiri" [Materials of the All-Russian Scientific and Practical Conference "Improving Efficiency of Energy Production and Use in Siberia"]. Irkutsk, IrGTU Publ., 2006, pp. 311-322 (in Russian).
3. Shalyt G.M. Opredelenie mest povrezhdenii v elektricheskikh setyakh [Identification of fault locations in electric networks]. Moscow, Energoizdat Publ., 1982, -312 p.
4. Vershinin E.A., Trofimov S.N. Dva metoda issledovaniya provodnykh struktur [Two research methods of wire structures]. Kemerovo, KGU Publ, 2005. Availabe at: http://www.ict.nsc.ru/ws/YM2005/9472/ (accessed 28 Merch 2016).
5. Ul'yanov S.A. Elektromagnitnye perekhodnye protsessy v elektricheskikh sistemakh [Electromagnetic transition processes in electrical power systems]. Moscow, Energiya Publ., 1970, 519 p. (in Russian).
6. Osak A.B. Prakticheskie metody rascheta tokov korotkogo zamykaniya v EES pri slozhnykh vidakh zamykanii s ispol'zovaniem faznykh koordinat [Practical methods of short circuit current calculation in the electrical power system under complex types of short circuits with the use of phase coordinates]. Tekhnicheskie nauki - ot teorii kpraktike [Engineering sciences - from theory to practice]. 2016, vol. 50, no. 2, pp. 181-192 (in Russian).
Критерии авторства
Маточкин С.Н., Тигунцев С.Г. и Осак А.Б. имеют равные авторские права. Ответственность за плагиат несет Ма-точкин С.Н.
Authorship criteria
Matochkin S.N., Tiguntsev S.G. and Osak A.B. have equal copyrights. Matochkin S.N. bears the responsibility for avoiding plagiarism.
Конфликт интересов
Авторы заявляют об отсутствии конфликта интересов.
Conflict of interests
The authors declare that there is no conflict of interests regarding the publication of this article.
Статья поступила 19.04.2016 г. The article was received 19 April 2016
