Программа определения места повреждения на воздушной линии электропередачи, страница 8

Матрица Y поврежденной ЛЭП будем считать функцией двух переменных: расстояния l до места повреждения и переходного сопротивления r короткого за­мыкания. Принцип ОМП состоит в подборе таких значений l и r, при которых расчётные значения токов I по концам линии отличаются от зарегистрированных аппаратурой I_изм на возможно меньшую величину, и.е. при которых длина вектора невязок

₂

минимальна.

На значения l и r наложим ограничения в виде допустимой вилки значений. Таким образом, задача ОМП сводится к поиску минимума функции двух пере­менных при наличии линейных ограничений типа неравенств. Для минимизации используется процедура BCPol математической библиотеки Фортрана, а блок-схема расчёта целевой функции представлена рис. 2.1.

		Рис. 2.1 Алгоритм определения места повреждения

 

Предварительно, до запуска процедуры минимизации, должна быть сфор­мирована матрица Y неповреждённой линии.

2.2 Формирование матрицы Y собственных и взаимных проводимостей        поврежденной воздушной линии электропередачи

          Поврежденную воздушную линию электропередачи можно представить в виде трех многополюсников объединенных в один многополюсник типа «Сборка».

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

 

Рис.2.2 Схема замещения поврежденной ЛЭП

Так как линия принимается короткой, следовательно, ее обобщенные пара­метры могут быть вычислены путем умножения погонных параметров на длину. До повреждения линия должна представлять собою обратимый двухполюсник. Соответствующие полюсам блоки матрицы Y должны удовлетворять следующим условиям:

Матрица Y собственных и взаимных проводимостей поврежденной ЛЭП формируется из матрицы Y неповрежденной ЛЭП, при этом матрица Y неповре­жденной ЛЭП должна находится в Y файле.

Алгоритм формирования матрицы Y поврежденной ЛЭП заключается в следующем.

1.  По матрице Y неповрежденной ЛЭП формируем матрицы Y многопо­люсников, соответствующих начальному и конечному участкам. Прини­маем, что собственные  проводимости прямо, а взаимные – обратно про­порциональны длине участков, то есть  .

2.    Пересчитываем матрицы Y начального и конечного участков, а также дополнительного двухполюсника вредителя с учетом разъемов, через  которые они подключатся к месту повреждения. Многополюсник-вреди­тель позволяет учесть, например, конечную проводимость дуги в месте повреждения. Число фаз многополюсника вредителя принимается рав­ным 2, так как учитываем только одно сопротивление короткого замыка­ния, а следовательно рассматриваем только двухфазные и однофазные короткие замыкания.

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

4.  Исключаем висячий и земляной полюса.

Опишем программную реализацию этого алгоритма.

Формирование матрицы проводимостей поврежденного многополюсника сосредоточено в модуле «DamageLine»

module DamageLine

use UnitName – Модуль описания дескрипторов рабочих файлов

use Global – Модуль описания глобальных переменных

use MatrADim – Модуль описания разъемов

implicit none

public:: MPDamLine

integer*2::  &

  Con1UN(2), &  – Пользовательские  номера разъёмов подключения на­чального участка к месту повреждения

  Con2UN(2), &  – Пользовательские  номера разъёмов подключения конеч­ного участка к месту повреждения

  ConFUN(2), &  – Пользовательские  номера разъёмов подключения полю­сов вредителя к месту повреждения

  kfm  – Количество фаз многополюсника

private:: Con1UN,Con2UN,ConFUN, DamLine, AddYMP

          Contains