- •Лекции по надежности электроэнергетических систем.
- •Лекция №1. Надежность и ее свойства. Категории потребителей по надежности. Причины и характер отказов объектов. Средства обеспечения надежности. Измерение надежности.
- •Лекция №2. Классификация отказов оборудования. Показатели надежности невосстанавливаемого элемента. Классификация отказов оборудования
- •Лекция №3. Показатели надежности восстанавливаемого элемента. Комплексные показатели надежности.
- •Лекция №4. Расчет показателей надежности одноэлементной схемы.
- •Лекция №8. Вероятности послеаварийных состояний сложных схем электрических соединений. Метод расчета надежности схем на основе формулы полной вероятности.
Лекция №4. Расчет показателей надежности одноэлементной схемы.
Функционирование элементов систем электроснабжения наиболее адекватно можно описать случайным процессом.
Процесс называется марковским, если для каждого момента времени вероятность любого состояния элемента в будущем зависит только от состояния его в настоящий момент времени и не зависит от того, каким образом элемент пришел в это состояние.
Процесс отказов и восстановлений обладает важным свойством: время пребывания элементов в состоянии безотказной работы или восстановления не зависит от того, сколько времени он уже провел в этом состоянии, и распределено по экспоненциальному закону.
Если процесс, протекающий в системе со счетным множеством состояний и непрерывным временем, является марковским, то его можно описать с помощью обыкновенных дифференциальных уравнений, в которых неизвестными функциями являются вероятности состояний.
Рассмотрим элемент, который может находиться в двух состояниях: 0 – исправной работы и 1 – отказа (элемент восстанавливается). Определим вероятности состояний p0(t) и p1(t) элемента в произвольный момент времени t при различных начальных условиях.
Эту задачу решим при условии, что поток отказов простейший с интенсивностью , закон распределения времени между отказами (частота отказов):
.
Время восстановления описывается также экспоненциальным законом распределения с параметром :
Сначала необходимо составить направленный граф состояний элементов (рис. 1), затем составить систему дифференциальных уравнений для вероятностей состояний.
Рис. 10. Расчетный граф состояний одноэлементной схемы.
Правило составления дифференциальных уравнений:
В левой части
каждого уравнения стоит производная
,
а в правой части – столько членов,
сколько ребер связано непосредственно
с данным состоянием; если ребро
оканчивается в данном состоянии, то
член имеет знак плюс, если начинается
из данного состояния – знак минус.
Каждый член равен произведению
интенсивности потока событий, переводящего
элемент или систему по данному ребру,
и вероятности того состояния, из которого
начинается ребро.
.
Система дифференциальных уравнений может быть использована для определения вероятностей безотказной работы электрических систем, функции и коэффициента готовности, вероятности нахождения в ремонте (восстановлении) среднего времени пребывания системы в любом состоянии на относительно коротких интервалах времени, когда необходим учет начальных условий (состояний элементов).
Решив эту систему уравнений при начальных условиях p0(0)=1, p1(0)=0, т.е. предположив, что в момент t=0 элемент находится в рабочем состоянии, получим:
Если в начальный момент времени элемент был в состоянии отказа (восстановления) p0(0)=0, p1(0)=1, то
Для установившегося
значения вероятности
,
в практических расчетах
коэффициент
готовности выражается формулой:
,
а коэффициент вынужденного простоя
.
Степень затухания экспоненты:
.
Это выражение получено в предположении, что время безотказной работы практически всех элементов электрических систем много больше времени аварийного восстановления. Формулы для практических расчетов можно преобразовать к виду:
при p0(0)=1, p1(0)=0:
при p0(0)=0, p1(0)=1:
Коэффициент
готовности или коэффициент вынужденного
простоя можно интерпретировать как
среднюю вероятность нахождения систем
соответственно в рабочем состоянии и
в состоянии отказа. Чем меньше среднее
время восстановления элемента, тем
больше коэффициент затухания
,
тем быстрее процесс стремиться к
установившемуся значению вероятности,
т.е. к установившемся значениям kг
и kп.
Обычно
при определении показателей надежности
электроснабжения, связанных с
технико-экономической оценкой последствий
от недоотпуска электроэнергии для
достаточно длительных интервалов
времени (сезон, год), без большой
погрешности вероятности состояний
системы можно считать по установившимся
средним вероятностям при
.
Такого рода режимы с точки зрения
надежности называются предельными
режимами элемента или системы.
Рис. 11. Зависимости коэффициента готовности и коэффициента вынужденного простоя от времени при различных начальных условиях.
Лекция №5. Расчет показателей надежности системы, состоящей из последовательных элементов. Расчет показателей надежности системы, состоящей из параллельно соединенных элементов. Расчет показателей надежности системы с произвольным соединением элементов. Расчет показателей надежности с учетом преднамеренных отключений элементов.
Расчет показателей надежности системы, состоящей из последовательных элементов.
Система, состоящая из n последовательных восстанавливаемых элементов, отказывает в тех случаях, когда отказывает любой из элементов. Данная система имеет два состояния: 0 – все элементы в рабочем состоянии, 1 – один из элементов в отказовом состоянии. Система дифференциальных уравнений выглядит следующим образом:
Решив эту систему уравнений при начальных условиях p0(0)=1, p1(0)=0, т.е. предположив, что в момент t=0 элемент находится в рабочем состоянии, получим:
Если в начальный момент времени элемент был в состоянии отказа (восстановления) p0(0)=0, p1(0)=1, то
Для установившегося значения вероятности , в практических расчетах коэффициент готовности выражается формулой:
,
а коэффициент вынужденного простоя
.
Вероятность отказового состояния за время T:
Если элементы
неоднородные, т.е.
,
то
Расчет показателей надежности системы, состоящей из параллельно соединенных элементов.
Параллельное соединение в смысле надежности элементов означает, что при отказе одного из элементов система продолжает выполнять свои функции, т.е. предполагается, что автоматическое резервирование каждого элемента с пропускной способностью, достаточной для полного обеспечения мощности потребителя.
В общем случае, когда система при таком резервировании состоит из n независимых элементов, число возможных состояний системы будет 2n, если считать, что каждый элемент может быть в двух состояниях – рабочем и отказовом – и отказ системы наступает когда все элементы окажутся в отказовом состоянии.
Рассмотрим параллельное соединение двух элементов (две цепи линии электропередачи, двухтрансформаторные подстанции и т.д.). Такая система может находиться в четырех состояниях: 1 – оба элемента в рабочем состоянии; 2 – первый элемент в отказовом состоянии, второй – в рабочем; 3 – второй элемент в отказовом состоянии, первый в рабочем; 4 – оба элемента в отказовом состоянии. Соответствующие вероятности этих состояний будут p1(t), p2(t), p3(t), p4(t).
Граф переходов для этой системы приведен на рис 12.
Рис. 12. Расчетный граф состояний системы.
Система дифференциальных уравнений, описывающих вероятности состояний во времени следующая:
Решение для стационарного состояния выглядит следующим образом:
Коэффициент готовности определяется следующим образом:
Коэффициент вынужденного простоя системы:
При рассмотрении одноэлементной системы было показано, что коэффициент затухания экспоненты обратно пропорционален времени восстановления элемента:
При условии, что
:
,
где
имеет смысл эквивалентного времени восстановления двух параллельно соединенных элементов.
Рассматривая систему с двумя параллельно соединенными элементами как один эквивалентный элемент, можем записать:
,
откуда
Параметр потока
отказов системы из двух резервирующих
друг друга элементов равен произведению
параметра потока отказов первого на
среднюю вероятность отказового состояния
второго плюс произведение параметра
потока отказов второго на среднюю
вероятность отказа первого элемента.
Иначе говоря, слагаемые
имеют смысл удельного числа отказов
системы при отказе соответственно
первого и второго элементов.
Для определения параметра потока отказов в случае n параллельно соединенных элементов необходимо рассмотреть столько слагаемых, сколько элементов в системе:
.
Расчет показателей надежности системы с произвольным соединением элементов.
Основным методом при исследовании надежности сложных ЭЭС используется метод структурного анализа схем электрических соединений. При расчете показателей надежности с помощью структурных схем анализируются не все возможные состояния схемы, а только состояния отказа того минимального набора элементов, которые приводят к отказу системы относительно рассматриваемого узла (минимальные сечения). Это дает возможность использовать известные приемы расчета, в частности применить формулу для суммы вероятностей совместимых событий – событий отказа сечений.
Вероятность отказа схемы относительно расчетного узла нагрузки:
|
(1.24) |
где
– количество элементов в сечении,
– вероятность отказа i-го сечения:
|
(1.25) |
здесь
– вероятность отказа j-го элемента
i-го сечения;
– число элементов в i-ом сечении;
|
(1.26) |
В расчет необходимо вводить сечения с числом элементов не более 2-3 в зависимости от необходимой точности расчетов. Здесь k1, k2, k3 – количество одно-, двух и трехэлементных сечений схемы относительно узла нагрузки соответственно.
Определение минимальных сечений по матрице минимальных путей.
Минимальный путь – минимальный набор элементов, которые обеспечивают нормальное функционирование схемы (передачу электроэнергии) от источника питания до узла нагрузки.
Сечения формируются по матрице минимальных путей исходя из следующих соображений: если элемент входит во все пути, то он является одноэлементным сечением, если хотя бы один из двух элементов i и j поочередно входит во все пути и i и j не являются одноэлементными сечениями, то они образуют двухэлементное сечение. Аналогичные рассуждения можно провести относительно трех, четырех и т.д. элементов.
Поэтому для формирования минимальных сечений целесообразно пути представить в виде матрицы путей П, столбцы которой соответствуют элементам схемы, строки – путям. Элементам этой матрицы присваивается значение 1, если элемент входит в соответствующий путь, 0 – если элемент не входит в путь.
По матрице путей можно определить все сечения в схеме относительно рассматриваемого узла нагрузки, анализируя сочетания столбцов матрицы П. Сечения рассматриваются в порядке возрастания числа элементов в них (одноэлементные, двухэлементные, трехэлементные и т.д.).
Для выявления
одноэлементных сечений по матрице путей
достаточно выявить ее столбец, состоящий
из одних единиц (пересечение всех путей
в одном элементе). Двухэлементные сечения
получаются в результате логического
сложения двух любых столбцов матрицы
путей (пересечение всех путей в двух
элементах). Если в результате получается
столбец, состоящий из одних единиц, то
эти два элемента образуют двухэлементное
сечение.
Трехэлементное сечение получается в
результате логического сложения трех
столбцов (пересечение всех путей в трех
элементах) и т.д.
Правила логического сложения:
Расчет показателей надежности с учетом преднамеренных отключений элементов.
Вероятность наложения аварийного отключения одного элемента на преднамеренное отключение другого (а не наоборот) определяется с учетом начальных условий. Следовательно, вероятность его аварии за время преднамеренного отключения рассматриваемого элемента меньше средней вероятности его аварийного отключения:
|
(1.27) |
Если
продолжительность преднамеренного
отключения принять равной
,
то
|
(1.28) |
где
-
коэффициент, зависящий (рис.13) от
соотношения времени аварийного
восстановления резервирующего
эквивалентного элемента и времени
преднамеренного отключения рассматриваемого
элемента.
Рис.13.
Зависимость понижающего коэффициента
k
от соотношения времени преднамеренного
отключения элемента
и эквивалентного времени восстановления
оставшейся части системы
.
В том случае, если
расчеты показателей надежности
выполняются на длительных интервалах
времени
,
то в практических расчетах можно не
учитывать начальные условия, и фактор
уменьшения вероятности совпадения
преднамеренного отключения одного
элемента и аварийного отключения другого
учитывается введением коэффициента
|
(1.29) |
Рис.14. Фрагмент ЭЭС.
На рис.14 изображен фрагмент системы электроснабжения, в которой имеются две параллельные ветви, состоящие из nв элементов каждая. Один из элементов цепи выводится в плановый ремонт на время Tпр. Вероятность отказа системы до ремонта, во время ремонта и после него описываются следующими выражениями:
|
|
|
|
|
|
|
|
|
|
Примерный вид зависимость вероятности отказа системы во времени – на рис. 15.
Рис. 15. Зависимость вероятности отказа системы во времени при отключении одного из элементов для проведения планового ремонта.
Система с последовательным соединением элементов:
В нерезервированных системах электроснабжения для уменьшения вероятности отключенного состояния и уменьшения количества перерывов электроснабжения стремятся по возможности совместить преднамеренные отключения элементов для ремонтов. Коэффициент вынужденного простоя такой системы с учетом преднамеренных отключений определяется по формуле:
где
- наибольшая вероятность преднамеренного
отключения одного из n
элементов системы.
Результирующий параметр потока отказов и эквивалентное время восстановления равны:
где
- наибольшая частота отключения одного
из n элементов схемы.
Система с параллельным соединением элементов:
Коэффициент вынужденного простоя:
Параметр потока отказов системы:
.
Лекция №6. Отражение логики функционирования электрических схем в расчетах показателей надежности. Метод расчета функциональной надежности ЭЭС. Комплексные показатели надежности. Пример расчета недоотпуска электроэнергии.
Отражение логики функционирования электрических схем в расчетах показателей надежности.
Отказы большей части элементов электрических систем: воздушных и кабельных линий, генераторов, трансформаторов, сборных шин, выключателей и разъединителей (в статическом состоянии) – обычно сопровождаются коротким замыканием. Отказы такого вида не локализуются в самих отказавших элементах, а приводят к необходимости работы релейной защиты и автоматики, воздействующей на коммутационные аппараты (выключатели, в низковольтных сетях – автоматы и предохранители), которые локализуют отказ элемента и ограничивают зону и продолжительность его воздействия на другие элементы.
Коммутационные аппараты могут сами отказывать при необходимости отключения (включения), т.е. при заявке на срабатывание; кроме того, не каждый элемент в электрической системе оборудован с двух сторон автоматическим коммутационным аппаратом, способным локализовать его отказ, поэтому зона его воздействия на систему в зависимости от схемы коммутации может быть весьма значительна, даже при отказе одного независимого элемента. При этом вынужденные отключения других элементов, обусловленные отказом рассматриваемого, классифицируются как зависимые простои. Продолжительность такого рода вынужденных отключений может быть меньше, чем продолжительность аварийного восстановления.
В расчетах структурной надежности предусматривается определение показателей относительно узлов нагрузки системы для состояний их полного отказа, т.е. на основе сечений схемы. Ранее рассматривались только основные сечения схемы, по которым электроэнергия передается непосредственно в узел нагрузки. Взаимосвязь между отказами элементов может порождать ситуации, когда отказы элементов, по которым не передается и не может передаваться электроэнергия в узел нагрузки, также приводят к состоянию полного отказа схемы относительно узла. Минимальные совокупности элементов, по которым непосредственно электроэнергия не передается или передается по части из них, но отказ которых приводит к потере питания узла, называются дополнительными сечениями схемы. Относительно узла формирование условий отказовых состояний в оценках структурной надежности сводится к определению основных и дополнительных сечений схемы.
Для каждого элемента основного сечения определяются все элементы системы, от которых на него может перейти отказ, т.е. вычисляются условные параметры потоков отказов и средние условные вероятности, последние – по продолжительности воздействия зависимых отказов. Затем последовательно вместо элемента основного сечения вводятся элементы, от которых может перейти отказ, тем самым получаются дополнительные сечения. Причем эти элементы вводятся своими условными параметрами потока отказов и соответствующей продолжительностью воздействия. Если основное сечение – одноэлементное, то параметры потоков и вероятности отказов основного и всех дополнительных сечений суммируются. Число дополнительных сечений, как правило, гораздо больше, чем основных, зависимость между отказами элементов как бы осуществляет “размножение” сечений.
Элементы дополнительных сечений образуются из различных списков. По продолжительности воздействия отказов и коэффициентам связи эти списки обозначаются как В, П, А. Так, если элемент не отделен от узла никакими коммутационными аппаратами, то его отказ приводит к отказу узла на продолжительность восстановления tв, а коэффициент перехода отказа равен единице (список В); если отделен, например, неавтоматическим коммутационным аппаратом, то его отказ приводит к отказу узла на продолжительность оперативных переключений tп, коэффициент перехода также равен единице (список П); если отделен автоматическим коммутационным аппаратом; то коэффициент перехода равен вероятности отказа qа выключателя, а продолжительность воздействия tп, а (список А).
При образовании дополнительных сечений с использованием указанных списков необходимо постоянно осуществлять контроль, не является ли полученное сечение основным и не включается ли какой-либо элемент в дополнительное более одного раза.
Метод расчета функциональной надежности ЭЭС.
При отключении отдельных элементов система находится на более низком уровне функционирования по сравнению с нормальным. При этом параметры режимов элементов и систем могут выходить за пределы, допустимые по нормам, что обуславливает отказ, отключение элементов или части системы, т.е. переход системы на еще более низкий уровень функционирования. Чтобы избежать этого, вводятся ограничения по режимам. К ним относятся: форсировка регулирующих и компенсирующих устройств, включение ненагруженных или малонагруженных резервных элементов, коммутационные изменения в схеме и т.д. Если эти мероприятия не обеспечивают введение режима элементов или всей системы в пределы допустимых норм, то отключается часть потребителей.
Рассмотрим определение значения недоотпуска электроэнергии для резервированной системы (рис.16), состоящей из двух элементов, один из которых отключен (j-й элемент), а в другом возникает перегрузка (i-й элемент).
Рис. 16. Схема резервированной сети.
Вероятность перегрузки γIi характеризует относительную долю времени, когда ветвь i будет перегружена. Если считать, что при любой перегрузке сверх расчетной в системе на соответствующее значение будет снижена нагрузка узлов, т.е. ветвь будет полностью разгружена, то недоотпущенная электроэнергия в этом состоянии равна:
|
|
где
– вероятность состояния;
;
Uномi – номинальное напряжение в
i-й ветви.
Однако в состоянии
перегрузки ветвь может быть разгружена
не полностью, а лишь на величину
,
превышающую допустимое значение тока
для этой ветви. Тогда недоотпуск
электроэнергии потребителям будет
минимальным, т.е. эта ситуация соответствует
нижней границе недоотпуска:
|
|
где
при нормальном законе распределения
|
|
На практике недоотпуск электроэнергии может находиться в диапазоне между верхним и нижним расчетными значениями. Весь анализ будет проводиться только для одного значения недоотпуска электроэнергии, под которым подразумевается его нижняя оценка.
При анализе ремонтных состояний ЭЭС необходимо оценивать функциональную надежность объектов ЭЭС как при "стационарном" состоянии, т.е. при отсутствии отказов элементов, так и при наложении отказов других элементов на плановые ремонты.
Комплексные показатели надежности.
Из комплексных показателей надежности практическое применение имеют два показателя: математическое ожидание недоотпущенной электроэнергии за период времени и экономический ущерб от перерывов электроснабжения.
Математическое ожидание недоотпущенной электроэнергии потребителям узла нагрузки вследствие полного прекращения его питания от системы:
|
(1.35)
|
или
|
(1.36)
|
где
– математическое ожидание активной
мощности потребителей узла;
– математическое ожидание недоотпуска
электроэнергии потребителям; Т – период
времени, для которого производится
расчет; Рнб, Тнб – максимум
суммарной мощности узла и время
использования максимума мощности
соответственно;
–
вероятность отказа относительно
расчетного узла.
Нарушение нормального функционирования объекта ЭЭС вследствие отказов элементов вызывает экономические санкции к собственнику объекта и как к покупателю электроэнергии и как к продавцу ее.
Как покупатель электроэнергии – вследствие отказа на своем объекте вызывает штрафные санкции со стороны поставщика, обусловленные снижением экономичности работы его оборудования (пережогом топлива, повышением уровня запасом воды на ГЭС, запиранием мощности и т.д.), составляющую Шпр.
Как продавец электроэнергии – вследствие недополучения прибыли от продажи электроэнергии (обозначим эту составляющую ПР) и штрафных санкций со стороны покупателя. Эта составляющая называется ущербом у потребителя У. В данном случае величина ущерба принимается равной Уд, согласно условиям договора между поставщиком и потребителем.
В общем случае недополученная прибыль объекта ЭЭС равна:
|
|
В зависимости от условий договора объекта ЭЭС с продавцом и покупателем электроэнергии первая и третья составляющие могут принимать любые значения, в том числе и нулевые, но ПР в принципе не может принимать нулевые значения, т.к. электроэнергия относится к физическим субстанциям.
В ряде случаев при оценке необходимой степени обеспечения надежности электроснабжения потребителей пользуются некоторыми средними экономическими показателями. Это связано с тем, что к электрическим сетям присоединяются различные промышленные предприятия, и точная оценка технико-экономического ущерба при перерыве их электроснабжения представляет значительные трудности. Например, при оценке ущерба потребителям можно использовать значение ущерба y0, отнесенные к 1 кВт*ч недоотпущенной электроэнергии.
При этом в случае прекращения электроснабжения потребителей ущерб потребителям равен:
У=у0Энпр. |
|
Удельный экономический ущерб y0 может многократно превышать стоимость недоотпущенной потребителю электроэнергии.
Пример расчета недоотпуска электроэнергии для резервированной схемы.
Исходные данные:
Нагрузка потребителя электроэнергии изменяется по равномерному закону в диапазоне от 20 до 100 МВт. Допустимая на один элемент – 45 МВт.
Показатели надежности элементов:
Показатель |
Значение |
|
1 |
|
2 |
|
400 |
|
500 |
|
3 |
|
2 |
|
600 |
|
700 |
Требуется определить недоотпуск электроэнергии потребителю за год.
Расчет показателей структурной надежности системы:
Вероятность отказа элементов:
Понижающие коэффициенты для учета преднамеренных отключений элементов:
Вероятность полного прекращения питания потребителя:
Расчет показателей функциональной надежности системы:
Все элементы находятся в работоспособном состоянии:
,
,
Первый элемент отказал или находится в плановом ремонте:
,
,
Второй элемент отказал или находится в плановом ремонте:
, ,
Суммарный недоотпуск электроэнергии потребителям:
Лекция №7. Режимные характеристики сложных ЭЭС. Метод пересчета параметров режимов и обобщенных параметров при изменении состояния системы. Основные методы структурного анализа систем. Пример расчета структурной и функциональной надежности системы.
Режимные характеристики сложных ЭЭС.
Для расчетов
токораспределения и напряжений узлов
используются обобщенные параметры
схем, матрица узловых сопротивлений
,
коэффициентов распределения C. Из
теории электрических сетей известно,
что токи ветвей и напряжения узлов
относительно напряжения базового узла
определяется следующим образом:
|
|
где
– матрица задающих токов узлов
(положительное направление принимается
к узлу).
Все обобщенные параметры целесообразно выразить через матрицу узловых сопротивлений Z, которая получается в результате обращения матрицы узловых проводимостей Yу:
Z=Yу-1=(MуMt)-1, |
|
где M, у – соответственно первая матрица инциденций и матица проводимостей ветвей (обычно диагональная); t – знак транспонирования.
Порядок обращаемой матрицы равен числу независимых узлов схемы n. Матрица коэффициентов распределения:
С=ГZ.
|
|
Таким образом, исходной информацией, характеризующей схему, являются: полностью заполненная матрица узловых сопротивлений Z размером nXn (в общем случае она получается обращением матрицы Yу порядка n); слабозаполненная матрица Гt размером nXm, где m – число ветвей схемы (она легко получается из матрицы Мt заменой элементов +1 и –1 на значения соответствующих проводимостей ветвей у с сохранением знаков, отражающих выбранные направления); диагональная матрица у (если не учитываются взаимные индуктивности между ветвями).
Метод пересчета параметров режимов и обобщенных параметров при изменении состояния системы.
Для оценки недоотпуска электроэнергии при любых моделях схем и режимов электропотребления необходим анализ значительного числа состояний. При этом оказывается нерациональным расчет параметров режимов с получением матрицы Z(r) для каждого состояния обращением матрицы Y(r). Основные причины этого:
конфигурация систем в ремонтных состояниях незначительно отличаются одна от другой, так как отключаются и включаются одновременно обычно небольшое количество элементов;
обращение матрицы Yу(r) является очень трудоемкой операцией, к тому же для каждого состояния необходима информация о всех элементах полностью заполненной матрицы Z(r), в противном случае не удается выявить элементы (узлы и ветви) с недопустимыми параметрами режимов.
С ростом размера анализируемой системы нерациональность такого подхода трансформируется, по существу, в невозможность практической реализации с помощью современных средств вычислительной техники.
Необходим метод пересчета матрицы Z, при котором выполнялись бы следующие требования:
использование единой информации о режимах и схеме для пересчета параметров режимов при большом количестве состояний системы. В качестве такой информации целесообразно использовать информацию для нормального состояния системы;
возможность ускоренного пересчета параметров режима и обобщенных параметров схемы с выявлением тех элементов, режимы которых выходят за допустимые пределы;
ускоренный пересчет любой части матриц обобщенных параметров без пересчета всей матрицы. Это имеет существенное значение при анализе схем большого размера.
При изменении состояния системы нагрузки отключаемых или включаемых ветвей, если узлы сохраняются в схеме, перераспределяются между оставшимися в работе, т.е. происходит приращение нагрузок ветвей.
Если, например,
отключаются r ветвей, то нагрузки
оставшихся ветвей
можно представить в виде:
|
|
где
– приращение нагрузок ветвей,
обусловленное изменением схемы и
перераспределением нагрузок коммутируемых
ветвей.
Значение этого
приращения для каждой ветви зависит от
нагрузки
отключаемых ветвей в доремонтном
состоянии, схемы системы и обобщенных
параметров, поэтому:
|
|
где gI – матрица распределения нагрузок отключаемых или включаемых ветвей, определяемая по изменению обобщенных параметров схемы.
Аналогичные рассуждения можно провести относительно напряжений узлов V. При отключении r ветвей
|
|
где gV – матрица распределения напряжений узлов, инцидентных отключаемым или включаемым ветвям, определяемая по изменению обобщенных параметров схемы;
– матрица разности
напряжений узлов, инцидентных коммутируемым
ветвям в доремонтном состоянии схемы.
Матрицы gV и gI (обобщенные параметры изменения схем) зависят только от исходных обобщенных параметров схем и изменений конфигурации.
В общем случае ремонтное состояние системы может реализовываться не только отключением r ветвей, но и включением s ветвей, так как при ремонтах производятся переключения в схеме.
Основным обобщенным параметром схемы системы, по которому определяются режимы нагрузок ветвей и напряжений узлов, является матрица узловых сопротивлений. Матрицы gV и gI получаются в результате пересчета элементов матрицы Z.
Анализ числовых
значений матриц
и
позволяет выявить зону влияния нагрузок
отключаемых ветвей на оставшиеся в
работе элементы схемы и наметить
мероприятия для введения режимов в
допустимые пределы в каждом состоянии.
Основные методы структурного анализа систем.
При оценке надежности ЭЭС невозможно обойтись без специальных методов и алгоритмов анализа структуры системы, на основе которых выявляются основные и дополнительные сечения схемы, базовые пути питания потребителей электроэнергии, анализируется схема выдачи мощности электрическими станциями и др. В основном все алгоритмы анализа системы построены на графах.
Рассмотрим два основных метода - поиск в ширину и поиск в глубину, на основании которых можно решить любой круг задач, связанных с анализом графов.
Поиск в ширину:
Пусть задан граф G=(V,E) и фиксированная начальная вершина s. Алгоритм поиска в ширину перечисляет все достижимые из s (если идти по ребрам) вершины в порядке возрастания расстояния от s. Расстоянием считается длина (число ребер) кратчайшего пути. В процессе поиска из графа выделяется часть, называемая “деревом поиска в ширину” с корнем s. Она содержит все достижимые из s вершины (и только их). Для каждой из них путь из корня в дереве поиска будет одним из кратчайших путей (из начальной вершины) в графе. Алгоритм применим как к ориентированным, так и к неориентированным графам.
Название объясняется тем, что в процессе поиска мы идем вширь, а не вглубь (сначала просматриваем все соседние вершины, затем соседей соседей и т.д.).
Построенный таким образом подграф Gπ графа G представляет собой дерево, в котором для каждой вершины ν имеется единственный простой путь из s в ν. Этот путь будет кратчайшим путем из s в ν в графе G.
Поиск в глубину:
Стратегия поиска в глубину такова: идти “вглубь”, пока это возможно (есть не пройденные исходящие ребра), и возвращаться и искать другой путь, когда таких ребер нет. Так делается, пока не обнаружены все вершины, достижимые из исходной. (Если после этого остаются необнаруженные вершины, можно выбрать одну из них и повторять процесс, и делать так до тех пор, пока мы не обнаружим все вершины графа).
Как и при поиске в ширину, впервые обнаружив вершину ν, смежную с u, мы отмечаем это событие. Получается дерево – или несколько деревьев, если поиск повторяется из нескольких вершин. Говоря дальше о поиске в глубину, всегда предполагается, что поиск повторяется.
Как алгоритм поиска в ширину, так и алгоритм поиска в глубину служат для формирования кратчайших путей из одной вершины графа в другую. Набор кратчайших путей, например, от узла нагрузки до источника(ов) питания необходим для определения минимальных сечений относительно рассматриваемой точки. Простейший способ определения основных сечений схемы – вынужденное пересечение нескольких путей в одном, двух и т.д. элементах. Так, например, если два независимых пути пересекаются в одном элементе, то этот элемент является одноэлементным сечением; если четыре независимых пути пересекаются в двух элементах, то эти элементы формируют двухэлементное сечение относительно узла нагрузки и т.д.
