Совершенствование средств управления

advertisement
На правах рукописи
ДУДИКОВ Юрий Сергеевич
СОВЕРШЕНСТВОВАНИЕ СРЕДСТВ УПРАВЛЕНИЯ
ЭКСПЛУАТАЦИЕЙ СИСТЕМ ЭЛЕКТРОСНАБЖЕНИЯ НА
ОСНОВЕ ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ
Специальность 05.09.03 – «Электротехнические комплексы и системы»
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
Самара – 2009
Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Самарском государственном техническом университете на кафедре «Автоматизированные электроэнергетические
системы».
Научный руководитель:
–
доктор технических наук, профессор
Гольдштейн Валерий Геннадьевич,
Самарский государственный
технический университет
Официальные оппоненты:
–
доктор технических наук, профессор
Таджибаев Алексей Ибрагимович,
Петербургский энергетический институт
повышения квалификации,
г. Санкт-Петербург
–
кандидат технических наук, доцент
Сенько Владислав Владимирович,
Тольяттинский государственный
университет, г. Тольятти
Ведущее предприятие:
–
Южно-Российский государственный технический университет (Новочеркасский
политехнический институт)
Защита состоится «17» февраля 2009 г. в 13 час. 00 мин. на заседании диссертационного совета Д 212.217.04 ГОУВПО Самарский государственный технический университет (СамГТУ) по адресу: г. Самара, Первомайская ул., д. 18,
1-й корпус, ауд. 4.
Отзывы по данной работе в двух экземплярах, заверенные печатью, просим направлять по адресу: Россия, 443100, г. Самара, Молодогвардейская ул.
244, Главный корпус, ГОУВПО Самарский государственный технический университет, ученому секретарю диссертационного совета Д 212.217.04,
тел.: (846) 278-44-96, факс (846) 278-44-00, e-mail: aees@samgtu.ru.
С диссертацией можно ознакомиться в библиотеке СамГТУ.
Автореферат разослан «15» января 2009 г.
Ученый секретарь диссертационного
совета Д 212.217.04, кандидат
технических наук, доцент
2
Е.А. Кротков
ОБЩАЯ ХАРАКТЕРИСТИКА ДИССЕРТАЦИОННОЙ РАБОТЫ
Актуальность темы. В процессе современного развития систем электроснабжения (СЭС) возникает необходимость в пересмотре традиционных подходов к решению ряда задач таких, как проектирование электрических схем
энергообъектов, анализ их надежности, организация диспетчерского управления (отслеживание текущего состояния схемы) и эксплуатации, получение
объективной информации об измеренных (ИП) и расчетных параметрах (РП)
режимов, расчет и анализ топологически близких режимов, подготовка персонала и проведение тренировок посредством тренажеров и аниматоров, диагностика оборудования. Для эффективного решения этих задач необходима гибкая
схемная графика, при применении которой возникает специфичная ситуация,
когда диспетчерские (ДС) и коммутационные (КС) схемы содержат до нескольких тысяч элементов, что значительно превышает топологические размеры традиционных расчетных схем.
Даже для современных вычислительных устройств решение режимных задач при этом усложняется. Хотя эти схемы содержат наиболее полную информацию, для анализа режимов все-таки необходимы расчетные схемы замещения, а его результаты должны возвращаться на ДС (КС). Кроме этого текущие
изменения в реальной сети на расчетной схеме своевременно не отражаются.
Можно констатировать актуальность возникающих при этом проблем:
- создание модели, которая наиболее полно отражает реальное состояние сети,
- гибкий учет всех происходящих в ней изменений,
- реализация ее преобразования в расчетную схему замещения и выполнение
на ней расчета режима,
- переход и отражение результатов анализа режима на исходную схему для их
оперативного использования в процессах управления, в тренажерах и др.
Проведенный анализ показал возможность коррекции традиционных
подходов к проблемам анализа, оценки и контроля состояния СЭС на основе совместного применения положений теории «распознавания образов»
(ТРО) и виртуальной имитационной модели (ВИМ) на основе табличных операторов, которые являются в этом направлении одними из наиболее эффективных
методологических инструментов.
Важность решения задач совершенствования методов и средств управления
состоянием СЭС неоднократно отмечалась на международных и отечественных
семинарах, конференциях, форумах, посвященных проблемам их эксплуатации, а
также в правительственных и отраслевых решениях. Сказанное выше подтверждает актуальность проблемы и темы диссертации.
Представленная диссертационная работа поддерживается грантом Ученого Совета СамГТУ № 5 от 2008 г.
Цель работы заключается в повышении эффективности и оперативности
решения эксплуатационных задач управления СЭС в части анализа режимов,
оценки, контроля состояния и реализации процессов тренажирования на основе совершенствования табличных методов.
3
Для реализации этой цели в настоящей работе решаются следующие задачи.
Научные задачи.
1. Научное обоснование развития возможностей тренажеров и имитаторов для
СЭС на основе табличных методов (ТМ).
2. Разработка методов моделирования режимов и подрежимов СЭС большого
топологического объема на основе применения табличных операторов, адекватных системам уравнений режимов СЭС.
3. Оценка состояния СЭС по данным телеизмерений (ТИ) и систем учета электроэнергии с учетом различий и наличия дефектов в составе ТИ с помощью ТМ.
Практические задачи.
1. Практическая реализация методов анализа установившихся (УР), аварийных
режимов и подрежимов КЗ в тренажерах и имитаторах СЭС.
2. Создание программного модуля для выполнения многовариантных расчетов
режимов в тренажерах и оценки состояния СЭС.
3. Построение алгоритма и разработка программного модуля для определения
параметров состояния СЭС по данным ТИ и систем учета электроэнергии с
различной полнотой и достоверностью информации с помощью таблично
топологических методов.
Основные методы научных исследований. Научные исследования в диссертационной работе основаны на ТРО, топологического анализа, множеств,
многомерных пространств и направленных графов, теории ИМ и графоаналитического моделирования режимов для изменяющихся квазистационарных состояний СЭС и др. Оценка корректности моделирования проведена путем сравнения
с результатами, полученными по данным эксплуатации СЭС, с помощью измерений, экспериментов и комплексных тестов работы моделей.
Основные положения, выносимые на защиту.
1. Научное обоснование и реализация взаимных переходов от реальной схемы
к виртуальным имитационным моделям, построенным на основе их графических представлений.
2. Совершенствование систем управления и тренажирования режимами и подрежимами с помощью табличных модификаций методов анализа УР и аварийных режимов СЭС на базе метода наложения.
3. Инвариантная оценка состояния СЭС с поиском и устранением дефектов
телеинформации.
Научная новизна работы определяется следующими результатами выполненных исследований.
1. Развитие систем тренажирования по технологическим режимам СЭС с помощью виртуальных имитационных моделей и теории распознавания образов в рамках единого информационно-топологического пространства.
2. Совершенствование методов моделирования режимов и подрежимов для
оперативного управления эксплуатацией СЭС непосредственно по табличным структурам, являющимися их виртуальными отображениями.
4
3. Методика оперативного восстановления параметров режимов СЭС при
наличии полной, избыточной и недостаточной информации по данным ТИ с
достоверизацией «портрета» этих режимов на основе отбраковки ошибочных данных в исходной и расчетной информации.
Практическая ценность работы определяется результатами следующих
разработок.
1. Практическая реализация в режимных и коммутационных тренажерах и проверка работоспособности методов анализа УР, аварийных режимов и подрежимов КЗ СЭС различного объема.
2. Реализация программного комплекса для оценки параметров режимов СЭС
по данным ТИ и систем учета электроэнергии с различной степенью полноты исходной информации.
3. Применение предложенных методов и программных средств в производственных, учебных и научно-исследовательских организациях.
Реализация в промышленности, проектной практике, учебном процессе и внедрение результатов.
Результаты представленной диссертационной работы реализованы в ряде
проектов для предприятий СЭС. В частности, за 2006-2008 г.г. они внедрены в
ПО филиала ОАО «МРСК Волги», Новокуйбышевском нефтеперерабатывающем заводе и др.
Материалы, изложенные в диссертационной работе и публикациях [1-12], выполненных автором лично и в соавторстве, используются в учебном процессе Петербургского энергетического института повышения квалификации, Самарского,
Ульяновского, и государственных технических университетов Тольяттинского
государственного университета, Самарского университета путей сообщения и др.
Обоснованность и достоверность научных положений, выводов и рекомендаций, а также практических выводов и программных продуктов базируется на корректном применении математических методов топологии, теории
графов и подтверждается адекватным поведением моделей по сравнению с
процессами на реальных энергообъектах, а также результатами измерений и
внедрения в составе ИАК "ПЕГАС" и программном комплексе «МОДУС».
Апробация работы. Основные положения диссертации докладывались и
обсуждались на 5 международных и 5 всероссийских научно-технических конференциях, а также научно-технических семинарах кафедры «АЭЭС» СамГТУ.
Публикации. По теме диссертации опубликовано 12 печатных работ (2
работы в журнале по списку ВАК).
Структура и объем диссертации. Работа состоит из введения, четырех
глав, заключения, списка литературы и приложений. Общий объем работы составляет 144 страницы. Библиография включает 105 наименований.
КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается научная и практическая актуальность, формулируются идея, цель, научные и практические задачи диссертационной работы.
5
Показана научная новизна и практическая значимость работы, приводятся основные положения, выносимые на защиту, сведения об апробации и внедрении
результатов работы.
В первой главе освещено современное состояние вопросов имитационного
графоаналитического динамического представления топологических структур
СЭС на основе табличного моделирования. Приведен подробный анализ современного состояния теории и методов анализа режимов СЭС. Выработан единый
подход к задачам оценивания параметров режимов СЭС на основе ИМ и ТРО. Дана также оценка основных результатов, которые изложены в работах отечественных и зарубежных авторов, применяются в практических приложениях и предлагаются к практическому применению по актуальным направлениям математического моделирования схем СЭС и их электрических режимов.
Большой вклад в решение вопросов построения ВИМ, построения эффективных алгоритмов формирования перехода от графических изображений схем
СЭС к виртуальным моделям в виде матриц и таблиц внесли такие известные
ученые и исследователи, как: Бартоломей П.И., Браун Дж., Гамм А.З., Гераскин О.Т., Голуб И.И., Гольдштейн В.Г., Демстер А., Егоров А.О., Жуков С.А.,
Жуков Л.А., Идельчик В.И., Качанова Н.А., Кохлас Дж., Крылов В.А., Крумм
Л.А., Мельников Н.А., Паздерин А.В., Стратан И.А., Строев В.А., Счастливый
Г.Г., Цукерник Л.В., Уоссермен Ф., Ушаков И.А., и др.
Основным в работе современного диспетчера СЭС и функциями тренажера является: наблюдение за состоянием объекта; диагностика и анализ различных ситуаций; планирование и выполнение действий по нормализации ситуации и принятие решений. Причем вся информация об операциях накапливается
в соответствующих журналах и базах данных.
Выполнение всех этих функций, равно как и анализ режимов опирается на
создание ВИМ. В работе выработан общий подход к решаемым задачам на
основе ТМ и ИМ. ВИМ режимов и топологии СЭС использует концепцию единого информационно-топологического пространства (ЕИТП), которая реализуется с
помощью ТМ и широко распространенных в настоящее время в эксплуатации СЭС
графических редакторов. Эта концепция включает в себя методологию, информационное и алгоритмическое обеспечение решения широкого круга технологических
задач качественного и количественного анализа режимов и оценки состояния СЭС,
которые выполняются с помощью информационно – аналитических и графических
программных комплексов.
Здесь для схем, моделирующих СЭС, основным инструментом служат графические представления, построенные с помощью стандартных графических
отображений элементов СЭС. В процессе их задания непосредственным образом
строится ВИМ СЭС в виде таблицы, содержащей диспетчерские и цифровые обозначения (топологию), формальные данные (марки оборудования, длины ЛЭП и
др.) элементов и расчетные параметры их схем замещения. Эти параметры автоматически формируются на основе сведений из базы данных по оборудованию.
6
На основе проведенного анализа известных работ и задач, возникающих
при управлении СЭС, сформулированы основные положения виртуального ИМ
СЭС и их режимов с помощью теоретико-множественного подхода, являющегося теоретической основой построения ЕИТП.
В качестве математического описания ЕИТП, не претендуя на исчерпывающую полноту, можно предложить имитационную модель СЭС в виде глобального информационно-топологического множества
М →{Мг, Мту, Мтк, Пд, Хо, Пс, Пр, Пи,…},
(1)
представленного здесь следующими подмножествами:
Мг – виртуальная графическая схема замещения, состоящая из множества графических изображений элементов СЭС, соединенных между собой в соответствие с реальным состоянием СЭС;
Мту, Мтк – табличные структуры в виде топологических списков, построенные,
соответственно, по узловому и контурному принципам;
Пд, Хо – паспортные данные и характеристики оборудования;
Пс – параметры схем замещения элементов СЭС, определенные в соответствие
с Пд и Хо;
Пр – параметры режимов СЭС, найденные с помощью расчетных процедур;
Пи, Пу – ИП режимов СЭС, переданные в информационную модель, соответственно, с помощью средств ТИ, телесигнализации, а также средств учета электроэнергии, в частности, АСКУЭ.
Основным назначением множественного описания (1) является реализация
ряда основных информационных и расчетных процессов ИМ, в полном объеме,
которые представлены в диссертации.
ТМ описания электрической схемы заключается в представлении информации о схеме сети в виде списковых табличных операторов. Основными из
них являются таблицы: узловых соединений (Ту); узловых характеристик (Тух);
контурных характеристик I и II рода (Ткх1, Ткх2) и др. Эти таблицы строятся
автоматически редактором «МОДУС”» при формировании графического изображения схемы, являясь составной частью ИМ, они дают возможность получить
решений уравнений состояния СЭС без их составления и подготовки.
Во второй главе реализуются табличные модификации основных методов
анализа и расчета УР и токов КЗ и, в частности, с помощью метода сопряженных
градиентов (МСГ). При использовании тренажеров в промышленных масштабах
объем данных, которые необходимо моделировать и анализировать, очень большой. Так, например, при моделировании городской сети в КС число ТП и РП вместе с коммутационными аппаратами составляет несколько тысяч.
Управление нормальными и аварийными режимами в этом случае ограничивается тем, оперативно-диспетчерский персонал (ОДП) не сможет использовать
большинство собираемой информации из-за ее разнообразия, огромного количества и отсутствия ЕИТП, что делает бессмысленным само внедрение АСУ-систем.
Фактически они выполняют только функции измерения и хранения данных. В работе показано, что результаты расчета УР для СЭС, имеющих несколько ступеней
7
номинальных напряжений, могут быть получены с использованием итерационных
процедур метода Гаусса-Зейделя, Ньютона-Рафсона и градиентного метода с помощью ТМ без формирования соответствующих уравнений.
Уравнения режимов в таких сетях составлены с учетом того, что: задано положительное направление от узла начала ветви (УНВ) p к узлу конца ветви (УКВ)
q; продольное сопротивление ветви Z вычислено для ступени номинального
напряжения сети, где находится УНВ p; идеальный трансформатор принимается
без потерь и находится всегда в УКВ; известен коэффициент трансформации K pq .
Сущность табличной модификации расчета УР по методу Гаусса-Зейделя
состоит в использовании узловых уравнений, которые для произвольного узла
р могут быть реализованы с помощью табличного оператора Тух:


 Y   K qp U
 Y 
Т ух  K pq U
q pq
q qp
q
q

,
 
U
(2)
p


2
Т ух  Y pq   K qp Yqp 
q
q

В выражении (2) в числителе и знаменателе присутствуют по две суммы.
В каждой паре сумм первая составляется по всем ветвям, для которых узел р
является УНВ; вторая – по всем ветвям, для которых узел р является УКВ.
Самым целесообразным способом построения вычислительного процесса в
этом случае представляется способ суммирования, при последовательном просмотре всех строк Тух. При этом для каждого пакета узловых характеристик p производится набор сумм, участвующих в выражении (2) и вычисляется новое значение напряжения узла p. По окончании просмотра Тух оказываются определены
напряжения во всех узлах схемы.
При построении итерационного процесса решения УР по методу Ньютона-Рафсона система нелинейных уравнений W (X) = 0, где переменными являются модули U и фазы δ напряжений узлов, а функция W соответствует узловым небалансам активной ∆P и реактивной ∆Q мощности на каждой i-ой итерации последовательно заменяется линейной системой, реализуемой с помощью табличного оператора Ту
δ 
 ΔP 
 W

Ту 
(3)
 X i   W X i 1  , где X    , W   
U
 X

 
ΔQ 
Ее решения дают значения неизвестных Xi, более близкие к решению названной нелинейной системы, чем предыдущие приближения Xi-1.
Элементы матрицы Якоби – это частные производные небалансов активной и реактивной мощности по модулям и фазам напряжений узлов, , которые
в расчетной процедуре по методу Ньютона-Рафсона реализуются в табличной
форме с помощью оператора Ту:
ΔP
ΔP
ΔQ
ΔQ
 W H N  
Ту 
,N=
,L=
,K=
.

 , где H =

δ
U
δ
U
 X  L K  

8

При построении итерационного процесса по градиентному методу в системе
нелинейных УУР, W (X) = 0, где вектор Х и W и составлены из вещественных U1,
I1 и мнимых U2, I2 значений узловых напряжений и токов, а функция W соответствует вещественным и мнимым значениям узловых небалансов по току.
Вводится неотрицательный функционал =(W^WТ), равный квадрату эвклидовой длины вектора W. Здесь обозначены символами «^» - скалярное произведение
векторов и «Т» - транспонирование.
На каждой i+1-ой итерации новые значения Xi+1 находятся в виде
Xi+1 = Xi - Ту gradX φ.t
(4)
с шагом t по вектору-антиградиенту, вычисленному в Xi, то есть по направлению наискорейшего спуска к абсолютному минимуму , соответствующему
положению электрического равновесия по первому закону Кирхгофа для токов
в узлах. Это дает значения неизвестных Xi+1, более близкие к решению линейной системы, чем исходное приближение.
~
 ~Iнб ~ 

 ~ ~
J ~ 

grad ~I   2Ту 
 Iнб   2Ту (Y  Iнб )  2D(
 Iнб ) ;
(5)
нб
U




 U

~
~
Iнб - вектор узловых токов небаланса; J - вектор узловых задающих токов.
В работе приведены сравнительный анализ рассмотренных методов, а
также принципы построения и основные позиции алгоритмов по выражениям
(2) – (5) с помощью табличных операторов Тух и Ту В них по определению вычисления проводятся только по реально существующим связям (ветвям), то
есть арифметические операции с нулевыми элементами в отличие от процедур,
построенных на основе матричных методов, полностью исключены.
Для расчета режима КЗ в узловой форме в работе предлагается модификация
МСГ непосредственно по Ту. При решении УР в комплексной форме производится их эквивалентное преобразование к вещественному виду и умножение
основной матрицы А системы на транспонированную матрицу, чтобы она стала
симметричной и положительно определенной. Предлагаемая модификация
МСГ в эквивалентной вещественной форме позволяет сократить количество
вычислительных операций в 2-4 раза по сравнению с комплексной формой.
Рассмотренные подходы к расчетам режимов КЗ и УР СЭС с использованием табличного способа описания ее топологической структуры реализованы
в программном комплексе «ПЕГАС», разработанном в СамГТУ. В частности
для фрагмента схемы электроснабжения Новокуйбышевского нефтеперерабатывающего завода выполнялись тестовые расчеты УР и режимов КЗ предложенными модификациями.
В рамках положений ИМ разработан испытательный программный комплекс (ИПК) для анализа режимов КЗ в схемах большого топологического объема на основе табличного подхода. Практическая проверка работоспособности и
оценка возможностей модификации МСГ, проведена по результатам серийных
расчетов режимов КЗ в сетях, в которых число узлов у > 1000.
9
Для этого с помощью датчика случайных чисел ИПК строит имитационные топологические модели СЭС в виде таблиц соединений и параметров ветвей. Кроме
того, моделируются и формируются следующие основные эксплуатационные признаки, характеристические соотношения и положения построения сетевых структур.
1. Различные уровни напряжения (≤
3) сети с произвольными процентными соотношениями состава по ступеням и связям между ними, количеству ветвей с источниками.
2. Произвольное соотношение между числом ветвей в и узлов у (от
разомкнутой до сложно замкнутой
схемы - 1,0 < в/у < 2,5).
3. Количество трансформаторных
и поперечных ветвей, в частности,
генераторных с характерными диаРис.1 Зависимости отношения числа итераций к числу узлов для схем сетей количе- пазонами сопротивлений для конством узлов у до 1000 и при в/у 1,5 < в/у < кретных классов напряжения.
2,5
4. Возможность выполнения расчета режима КЗ в вещественной и комплексной
форме. В последнем случае – при произвольном соотношении активного и реактивного сопротивлений ветвей.
После окончания процесса расчета конкретного режима, когда реккуррентно
вычисляемая невязка узловых небалансов по току становилась меньше наперед
заданной погрешности, производилась дополнительная независимая от решения
по МСГ проверка полученного решения по условиям оценки положения электрического равновесия. Иначе говоря, по рассчитанным напряжениям в узлах определялись токи в ветвях и в последовательном процессе просмотра строк Ту формировались балансы токов в узлах.
На рис. 1 представлены результаты расчетов режимов КЗ в схемах, имеющих
до 1000 узлов при различных показателях замкнутости схемы. Для объективной
демонстрации численной устойчивости алгоритма верхняя граница по количеству
узлов была выбрана значительно больше, чем это требуется в настоящее время.
По результатам проведенных расчетов можно сделать следующие выводы.
1. Ни в одном из вычислительных экспериментов в диапазоне до 1000 узлов не
было зафиксировано случаев численной неустойчивости.
2. Отмечена характерная тенденция относительного уменьшения числа шагов в
алгоритме МСГ при увеличении числа узлов.
3. Слабая зависимость числа шагов вычислительного процесса от диапазона возможных значений сопротивлений ветвей схемы, коэффициента замкнутости
схемы, числа генераторных и трансформаторных ветвей.
В третьей главе с помощью методов табличного анализа решаются задачи,
часто возникающие в процессах тренажирования и управления СЭС. Они связаны
10
с расчетами подрежимов, отличающихся от режима, принятого за основной некоторыми изменениями параметров схемы и режима.
Характерными случаями этих изменений являются отключение существующих и подключение новых ветвей схемы, перемещение точки КЗ по ветви, сети и др. При этом нередко в расчетах возникают численная неустойчивость или полное расхождение результатов с предполагаемыми. Решение задач, связанных с этими изменениями возможно двумя способами: первый из
них – изменение топологии сети. В результате перестраивается ВИМ, и по ней
любым методом расчета режимов находится новое решение. В качестве второго способа, прежде всего, для режимов КЗ в сложнозамкнутых сетях для анализа
подрежимов предлагается использовать метод наложения.
В частности, хорошо известно, что с помощью этого метода действительный режим КЗ можно рассматривать как результат наложения собственно аварийного режима на предшествующий доаварийный.
При решении систем уравнений узловых напряжений и контурных токов,
составленных в комплексной форме, с основными матрицами узловых проводимостей Yу и контурных сопротивлений Zк для них применяют общепринятое LU разложение на произведение двух треугольных матриц. С учетом слабого заполнения оно реализуется в так называемой упорядоченной форме (в виде табличных списков) для минимизации появления новых нулевых элементов.
Факт изменения основных матриц систем уравнений узловых напряжений
и контурных токов в подрежимах очевиден. Поэтому в работе предлагается
методика расчета подрежимов без изменений матриц Yу и Zк основного режима,
но с их косвенным учетом с помощью метода наложения. При этом становится
возможным самую трудоемкую операцию L-U разложение для анализируемой
расчетной схемы выполнить только один раз. В работе даны расчетные процедуры учета в правых частях основной системы уравнений физических условий и
изменений, связанных с моделированием подрежимов.
В частности, по результатам решения систем уравнений с помощью L-U
разложения, можно найти комплексный вектор токов ветвей I и исходного режима, а далее - токораспределение I для подрежима, отличающегося, наприр
мер, от исходного режима совместным отключением и подключением ветвей
j, k, l, f. Для этого на исходный режим (рис. 2) накладывается искусственный
 , E ,
(рис. 3), в котором в отключаемых ветвях включены источники э.д.с. E
рj
рk
создающие в этих ветвях совместным действием токи I рj , I рk , равные по величине токам исходного режима I иj , I иk , но противоположные по направлению. В результате наложения искусственного режима на исходный, токи в ветвях j, k будут отсутствовать. В подключаемых ветвях вводятся источники тока J l ,
J суммарное действие которых должно обеспечивать на каждом из них напряf
11
 , равное тому, которое было между узлами подклю , U
жение подключения U
иf
иl
чения каждой новой ветви в исходном режиме.
Рис. 2 Исходный режим
Рис. 3 Искусственный режим
При произвольном числе отключаемых nо и подключаемых nп ветвей система
уравнений для определения значений э.д.с. и токов соответственно источников
напряжения и тока может быть записана в матричном виде в комплексной форме
  K J  I  I ; 
Ypо E
pо
m рп
pо
u 
(6)





K н Epо  Z рп J рп  U рп  U u .

где I - ток, создаваемый источником напряжения в искусственном режиме при
ро
отключении линии, равный по величине исходному току I и , но противополож - напряжение на зажимах источника тока, вводиный ему по направлению; U
рп
 ; Yро и
мого в схему при подключении линии, равное исходному напряжению U
и
Zрп - подматрицы узловых проводимостей и контурных сопротивлений при отключении и подключении линий соответственно; J рп и E ро - источники тока и
э.д.с при подключении и отключении линии, соответственно.
Число уравнений в системе (6) равно (nо+ nп). Размерности векторов E ро , I ро
 , a также подматриц Yро и Zрп определяются соответственно числами nо
и J рп , U
рп
и nп. Подматрица Кт, имеющая размерность (nо, nп), определяет коэффициенты
связей по току ветвей с источниками тока и ветвей с источниками напряжения,
или, иначе говоря, между подключаемыми и отключаемыми ветвями. Аналогично
подматрица Кн, имеющая размерность (nп, nо), определяет коэффициенты связей
по напряжению ветвей с источниками напряжения и ветвей с источниками тока,
или между отключаемыми и подключаемыми ветвями.
Определение элементов подматриц Yро и Zрп производится по результатам
расчетов парциальных режимов (ПР). ПР определяются по тому же треугольному разложению основной матрицы сети. При этом для наиболее простого случая,
когда отключается одна ветвь, рассматривается искусственный режим, в котором
все источники сети отключены, а в рассматриваемой сети включен единичный
12
источник эдс, направленный против тока предшествующего режима. Этот режим
и получил название парциального.
Теперь после решения для него основной системы уравнений, по соотношению
предшествующего и парциального тока данной ветви определяется коэффициент
пропорциональности, позволяющий получить значение эдс источника из граничного условия по равенству тока отключаемой ветви току предшествующего режима, а
также все токи и напряжения от этого одного искусственного источника. Накладывая этот режим на исходный, получаем ноль – отсутствие тока в отключенной ветви
и, соответственно, токи и напряжения режима в новом состоянии сети.
Аналогично источники тока определяются в ПР и для подключения ветвей,
только в этом случае подключаются единичные источники тока, и определение
граничных условий производится по напряжениям.
С использованием предлагаемых принципов в работе решаются и другие задачи: изменение параметров ветвей, в том числе и коэффициентов трансформации, подключение и отключение источников в сети, и изменение их режимного
состояния, перемещение точки короткого замыкания по узлам и ветвям сети и т.д.
В результате решения системы (6) определяются значения э.д.с. и токов источников напряжения и тока, являющиеся коэффициентами пропорциональности изменения величин токов парциальных режимов и токораспределение I c
подрежима совместного подключения и отключения ветвей:
 I  J I  ...
I c  I u  E
(7)
pj j
пl l
Предлагаемая методика анализа подрежимов КЗ реализована для ряда схем
структурных предприятий ОАО «МРСК Волги».
В четвертой главе констатируется, что в настоящее время в СЭС измерительные датчики и преобразователи, устройства сбора телеинформации, средства
связи и пр. имеют большой физический износ и морально устарели. Кроме того,
на подстанциях, вторичные цепи нередко работают с перегрузкой. Все это приводит к тому, что ОДП довольно часто не имеет всей необходимой информации
о текущем режиме СЭС, или она неточна или недостоверна. В результате на
ОДП значительно увеличивается аналитическая нагрузка, т.к. необходимо не
только анализировать сам режим и, но также оценивать и учитывать возможную
неточность имеющейся у него информации. При этом ОДП может принимать
неверные решения, приводящие даже к аварийным ситуациям.
Решение проблемы восстановления параметров режима по данным ТИ с
оценкой достоверности измеренных и расчетных данных основывается на использовании таблично-топологического узлового метода (ТТУМ), реализованного в ЕИТП.
Из существующих методов оценки режимов СЭС широко распространены
матричные методы. Основным моментом в этом случае является построение
матриц инциденций и операции с ними. При реализации матричных алгоритмов
необходимо применять специальные меры по формированию, упаковке и обработке матриц инциденций, что делает их громоздкими. Кроме этого, эти методы
13
в построении портрета режима используют весь набор приходящих данных, в
том числе с ошибками, которые интегрально влияют на конечный результат.
Предлагаемый метод позволяет определить искомые параметры режима в
комплексной форме по измеренным модулям напряжения и по составляющим
активной (P) и реактивной (Q) мощности на основе одной из заданных форм
исходной информации и уравнений четырехполюсника.
1. Модули напряжений КУ (U) и потоки P и Q мощности по ветвям дерева
схемы сети, связанных с данным КУ (форма PQU).
2. Потоки P и Q и модули токов (I) по ветвям дерева схемы, связанных с данным КУ (форма PQI).
3. Модули U в КУ, потоки Р и модули I в ветвях дерева схемы, связанных с
данным КУ (форма PUI).
4. Модули U в КУ, потоки Q и модули I в ветвях дерева схемы, связанных с
данным КУ (форма QUI).
Форма PQU
1. КУ на входе четырехполюсника



S

1ˆ 2 


(8)
U 2  Т уд  DU 1  B  


U
1


2. КУ на выходе четырехполюсника

 
 1
2
(9)
U 2  Т уд  ( A  U 2  B S12ˆ )

 U1

Форма PQI
 для этого случая совпадает с (8).
1. КУ на входе четырехполюсника. Выражение U
2
2. КУ на выходе четырехполюсника.
 S


U

1 12ˆ
 Т 
(10)
U
2
уд 
2

A
S

B
I


ˆ
2
 12

В выражениях (8-10) Туд – это оператор табличной реализации расчета уз определовых напряжений по ветвям дерева схемы. Для форм PUI и QUI U
2
ляется по выражениям, записанным для двух предыдущих вариантов задания
исходной информации, но только с учетом правильного выбора направления
потоков соответствующих мощностей. При этом можно выделить ветви, для
которых условные направления потоков Р и Q совпадают. При произвольном
направлении потоков Р и Q между двумя соседними узлами эти мощности
необходимо определять, выполнив четыре расчета при всех возможных сочетаниях знаков. Расчет с правильным выбором знаков даст совпадение расчетного модуля напряжения узла конца с заданным.
Названные уравнения записаны относительно входа и выхода четырехполюсника, что соответствует условному началу (НВ) и концу (КВ) ветви. При
14
этом для расчета режима текущей ветви дерева возможны только два случая: а)
модуль измеренного (ИН) и комплекс рассчитанного напряжения (РН) находятся в узле условного НВ (по ходу дерева); б) модуль ИН находится в конце, а
комплекс РН - в узле условного НВ. В обоих случаях необходимо определить
комплекс напряжения U 2 в узле, являющемся условным КВ.
В диссертации представлены эффективные алгоритмы выбора состава дерева и расстановка расчетных КУ по одному из деревьев электрической схемы
с произвольным выбором одного из них в качестве узла корня дерева.
Полученные расчетные выражения являются универсальными, так как при
решении данной задачи с помощью уравнений режима четырехполюсника могут быть представлены схемы замещения любых элементов системы. В работе
предложены расчетные выражения для различных схем замещения ветвей.
Предлагаемый метод практически не зависит от полноты исходной информации. В условиях избыточности всей схемы вся информация по узлу, а
также о потоках мощности в ветвях используется для контроля результатов и
самих измерений. В ненаблюдаемых участках схемы используются псевдоизмерения, которые формируются на основе данных контрольных замеров. Их
точность определяется весовыми коэффициентами, которые назначаются по
результатам экспертных оценок с использованием статистических методов
оценки погрешностей измерений замеров. При задании псевдоизмерений их
точность устанавливают в несколько раз меньше точности реального телеизмерения подобного вида в зависимости от качества исходной информации.
В процессе решения задач восстановления режимов СЭС по данным ТИ и
учета электроэнергии встает кардинальная проблема определения соответствия
портрета – образа расчетного режима измеренным данным. При использовании
последних необходимо учитывать их несовершенство, вызванное наличием систематических и спорадических (случайных) ошибок, отклонений и сбоев.
По результатам анализа статистических данных о нарушениях и сбоях
в работе ТИ структурных подразделений ОАО МРСК Волги в работе пр иведены классификация ошибок поступающих ТИ и способы их устранения.
С позиций ТРО основной задачей анализа СЭС является приближение так
называемого виртуального портрета режима к реальному режиму сети. Другими
словами необходимо установить сходства и различия между ИП и РП сети. Переход от реального режима к его виртуальному портрету осуществляется на основе ИМ, реализуемого с помощью ТМ и установления критериев соответствия
между реальной сетью и ее аналогом.
На основе ТРО в работе сформированы формализованные критерии к оценке
правдоподобия образа. Далее в таблице 1 приведены критерии правдоподобия
образов – портретов режимов СЭС. Эти критерии позволяют проверить достоверность получаемых оценок всей системы в целом и по каждому узлу, параметру
в отдельности.
15
Таблица 1
Критерии правдоподобия портретов режимов СЭС
№
п/п
1
2
3
4
5
6
7
8
9
10
16
Формульное представление
Дифференциальный
| ξi - ξj | ≤ εп, где
по
параметру а)  | n | ;  | n | ;
i
j
i
j
(напряжению в узле,
току, потоку и поте- б) ξi=Re (ni ) ; ξi=Re ( n j ) ;
рям мощности в
в) ξj=Im ( n j ) ; ξi=Im ( n j ) ;
ветви)
Интегральный - двоичный (по факту
Σу[|Re (n k ) | + |Im (n k ) |]≤ εд
приращения напряжений всех узлов)
Вид критерия
Дифференциальный
узловой по мощностям (токам)
φк =| ∆ Sk |2 ≤ εSк
Интегральный узловой по мощностям
(токам)
Дифференциальный
по электроэнергии,
поступившей и отпущенной в узле
  k  ε Sk
Sk   Skq
q
q
|ΣkWk ПОСТ – ΣkWk ОТП| ≤ εW
Интегральный
по
электроэнергии, поΣW ПОСТ i - ΣW ОТП j ≤ εΣW
ступившей и отпущенной в сеть
Дифференциальный
по потерям электроΣk∆Wi - Σk∆Wj ≤ ε∆W
энергии в отдельной
ветви
Интегральный
по
потерям
электроΣ∆Wi - Σ∆Wj ≤ εΣ∆W
энергии во всей сети
Интегральный
по
суммарной эвклидовой норме квадратов Σку (U ку изм - Uку расч)2 ≤ εдU2
разностей ИН и РН
для всех КУ.
Дифференциальный
по модулю разности |U ку изм - Uку расч| ≤ ε д∆U
ИН и РН в КУ.
Примечание
εп – допустимая погрешность
(ДП) по параметру п в текущих
образах (портретах) режима
ПРi и ПРj
εд – допустимая мера отличия
от нуля суммы по всем узлам k
текущих приращений узловых
параметров.
εSк –ДП узловых квадратичных
сумм потоков мощности (токов) по всем q ветвям, сходящимся в узле k.
εSк – ДП – эвклидова длина
векторов небалансов мощности (токов) по всем узлам.
εW – ДП между поступившей Wк
ПОСТ и отпущенной Wк ОТП электроэнергией при составлении
баланса по узлу k.
εΣW – ДП между WПОСТ и WОТП
электроэнергией в текущих
образах (портретах) режимов
ПРi и ПРj при составлении баланса по всем узлам сети.
ε∆W – ДП по потерям электроэнергии в текущих образах
(портретах) режимов ПРi и ПРj
εΣ∆W – ДП по потерям электроэнергии всей сети в текущих
образах режимов ПРi и ПРj
ε дU2 - ДП по эвклидовой норме разностей ИН и РН
ε д∆U – ДП разности ИН и РН в
КУ
Продолжение таблицы 1
11
12
Интегральный
по
разнице (п.9) измеренных (ИМ) и расчетных (РМ) мощностей для всех КУ.
Дифференциальный
по модулю разности
ИМ и РМ в КУ.
Σ pq (S pq изм - S pq расч)2 ≤
ε д∆S2
ε д∆S – ДП по сумме квадратов
разности ИМ и РМ по ветвям.
|S pq изм –S pq расч|≤ ε д∆S
ε д∆S - ДП разности ИМ и РМ
по ветвям в КУ.
Наиболее частыми ситуациями являются нарушения и сбои в информации
по потоку мощности одной ветви, напряжения по одному или в целом по всем
присоединениям КУ. Значения ИМ Р и Q по ветвям, сходящимся в конкретном
КУ, имеют оценку достоверности по критериям 3, 4 и 11, 12.
Для оценки достоверности и качества модулей напряжений в КУ предлагается метод репит-контроля. В его интегральном варианте рассматривается режим всей сети и используется ТТУМ. В дифференциальном варианте комплекс
напряжения в КУ в КВ может быть однозначно определен по известному комплексу напряжения в НВ и измеренному комплексу мощности в конце.
Рассматривая эти оценки в целом и по каждому КУ в отдельности сделаны
выводы о качестве конкретного пакета ТИ.
1. Не выполнение условий проверки по всем узлам, начиная с опорного КУ –
сбойный узел. Повтор расчета от других узлов.
2. Отсутствие отклонений в параметрах в опорном КУ и наличие в некоторой
их части дерева схемы – коррекция отбракованных измерений путем сравнения
расчетов от расчета различных опорных узлов.
3. Выполнение условий проверки по всей совокупности КУ – отсутствие дефектов ТИ в КУ.
Решение задачи оценки параметров режима СЭС по данным ТИ проводилось для ряда схем Самарского производственного отделения, филиала ОАО
«МРСК Волги».
В заключении отражены основные результаты исследований в соответствии с поставленными задачами, решение которых обеспечило достижение
цели диссертационной работы.
1. Научно обоснована концепция взаимных переходов между реальной схемой
электроснабжения и ВИМ, на основе таблично-топологических методов и
ЕИТП, которые в 3-4 раза повышают оперативность управления топологически близкими режимами, и возрастает эффективность подготовки ОДП СЭС.
2. Для режимных и коммутационных тренажеров научно обоснованы модифицированные методы решения задач расчета УР и аварийных режимов, позволяющие оперативно решать задачи качественного и количественного анализа
и управления СЭС.
3. Реализованы алгоритмы анализа УР и аварийных режимов СЭС в составе ИАК
Пегас и графического редактора МОДУС, по точности выполнения расчетов
17
сопоставимые с известными матричными аналогами, но сочетающие достоинства табличных представлений и возможности анализа непосредственно по реальному графическому отображению СЭС.
4. Расширены возможности тренажирования и диспетчерского управления аварийными режимами и подрежимами СЭС различного топологического объема
за счет применения принципа наложения, и проведена их опытная проверка.
5. Разработан ИПК по расчету аварийных режимов СЭС различного топологического объема с помощью табличного имитационного моделирования и модификации МСГ, обеспечивающего численную устойчивость и уменьшение числа итераций при увеличении числа узлов схемы.
6. Предложена методика оценки параметров режимов СЭС на основе ТТУМ,
позволяющая повысить оперативность выполнения расчетов в 2-3 раза и до 9598% увеличить достоверность отбраковки ошибок в исходной и рассчитанной
информации по данным ТИ за счет предложенных критериев правдоподобия.
7. Разработан программный пакет по анализу топологии и оценки параметров режима СЭС по данным ТИ и систем учета электроэнергии наряду с эффективными алгоритмами детекции и отбраковки ошибок в исходной информации.
В приложении приведены акты внедрения работы в предприятиях, проектных
организациях и учебном процессе.
Основное содержание работы отражено в следующих публикациях
В изданиях по списку ВАК:
1. Дудиков Ю.С. Методика расчета подрежимов коротких замыканий с постоянной матрицей узловых проводимостей в электрических сетях [Текст] / Ю.С.
Дудиков // Известия вузов. Электромеханика. - 2007. - №6. - С. 83-87.
2. Дудиков Ю.С. Проблема восстановления параметров режима электрической
сети по данным ТИ и ТС, связанная с детекцией ошибок в исходной информации / В.Г. Гольдштейн, Ю.С. Дудиков [Текст] // Известия вузов. Электромеханика. Спец. Выпуск «Электроснабжение». - 2007. – С. 80-81.
В других научных изданиях:
3. Дудиков Ю.С. Теоретико-множественный поход к имитационному моделированию электрических сетей [Текст] / Ю.С. Дудиков, А.В. Романов // «Тинчуринские чтения»: Материалы II-ой молодежной межд. научн. конф. – Казань:
КГЭУ. - 2008. - С.146-147.
4. Дудиков Ю.С. К вопросу о реализации метода сопряженных градиентов в
задачах анализа и оптимизации режимов электрических сетей [Текст] / Ю.С.
Дудиков // «Математическое моделирование и краевые задачи»: Труды IV-ой
Всероссийской научной конференции с международным участием. Часть 2 –
Самара: СамГТУ. - 2007. - С. 62-65.
5. Дудиков Ю.С. Применение метода наложения в расчетах подрежимов коротких замыканий [Текст] / В.Г. Гольдштейн, Ю.С. Дудиков // «Энергетика: состо-
18
яние, проблемы, перспективы»: Труды Всерос. научн.-техн. конф. – Оренбург:
ОГУ. – 2007. - С. 29-34.
6. Дудиков Ю.С. Учет оперативных изменений в электрической сети при решении
задач анализа режимов с неизменной матрицой узловых проводимостей / В.Г.
Гольдштейн, Ю.С. Дудиков [Текст] // «Энергетика и энергоэффективные технологии»: Труды II Межд. научн.-техн. конф. – Липецк: ЛГТУ. - 2007. - С. 211-215.
7. Дудиков Ю.С. Табличный способ построения моделей электрических сетей /
В.Г. Гольдштейн, Ю.С. Дудиков, Е.В. Подшивалова [Текст] // «Энергетика:
экология, надежность, безопасность»: Материалы XIII Всерос. научн. конф. Томск: ТПУ. - 2007. - С. 53-57.
8. Дудиков Ю.С. Табличный способ построения моделей электрических сетей /
В.Г. Гольдштейн, Ю.С. Дудиков, Е.В. Подшивалова [Текст] // «Наука. Технологии. Инновации»: Материалы Всерос. научн. конф. молодых ученых. - Новосибирск. Часть 3: НГТУ. - 2007. - С. 210-212.
9. Дудиков Ю.С. Применение табличного метода для дорасчета режима электрической сети по данным оперативно-информационных комплексов [Текст] / Ю.С.
Дудиков // Вестник СамГТУ. Технические науки. – 2007. - №2(20). - С.141-147.
10. Дудиков Ю.С. Применение метода сопряженных градиентов при табличной
реализации режимов электрических сетей [Текст] / В.Г. Гольдштейн, Ю.С. Дудиков // «Оптимизация режимов работы электротехнических комплексов»:
Межвузовский сборник научных трудов. – Красноярск: СБУ. - 2008. - С. 180 189.
11. Дудиков Ю.С. Метод наложения в задачах расчета режимов коротких замыканий электротехнических комплексов [Текст] / В.Г. Гольдштейн, Ю.С. Дудиков,
Е.В. Подшивалова // «Радиоэлектроника, электротехника и энергетика»: Тезисы
докладов 14-й межд. научн.-техн. конф. студентов и аспирантов. 28–29 февраля
2008 г. – М.: МЭИ, 2008. 3 том. – С. 259 - 260.
12. Дудиков Ю.С. Восстановление режима электрической сети в условиях избыточной информации по данным телеизмерений [Текст] / В.Г. Гольдштейн, Ю.С.
Дудиков // «Перенапряжения и надежность эксплуатации электрооборудования»:
Международная научно-техническая конференция, Выпуск 6. – Спб.: ПЭИПК,
2008. – С. 135-140.
Личный вклад автора. Основные научные результаты, включенные в диссертацию, опубликованы в 12 работах. Статьи [1, 3, 9] написаны лично. В
остальных - автору принадлежат: общая постановка научных проблем, путей и
методов построения основных решений [2, 4, 6]; в работе [5] - построение математических моделей и алгоритмов, реализация математических моделей; в работах [7, 8, 11] – формализация процессов построения моделей; в работах [9, 10,
12] – постановка задач, выполнение расчетов и обобщение и обсуждение результатов расчётов и исследований. Кроме того, во всех публикациях выполнены
19
редактирование работ при представлении их в печать и их коррекция по замечаниям рецензентов.
20
Автореферат отпечатан с разрешения диссертационного совета Д212.217.04 ГОУ ВПО Самарский
государственный технический университет
(протокол № 13 от 25 декабря 2008 г.)
Заказ № 960 Тираж 100 экз.
Отпечатано на ризографе.
ГОУВПО Самарский государственный технический университет
Отдел типографии и оперативной печати
443100 г. Самара ул. Молодогвардейская, 244
21
Download