Оценка технического состояния трансформаторных вводов на основе нечетких алгоритмов

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

q=Ка,

где К - градуировочный коэффициент измерительного устройства, включенного в испытательную схему.

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

- индикатор ЧР, измерение с помощью которого производится во время испытания изоляции приложенным повышенным напряжением, а также в случаях, когда в условиях эксплуатации ввод оборудован стационарным устройством присоединения;

- сигнализатор ЧР для автоматического непрерывного контроля в условиях эксплуатации;

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

С целью ранней диагностики применяют методы измерений, выявляющие ЧР с интенсивностью меньшей, чем интенсивность критических разрядов (порядка 10-9 Кл). Для сигнализации предаварийного состояния достаточно выявить разряды критической интенсивности. При этом необходим частый или даже непрерывный контроль. Имеются данные об успешном выявлении дефектов изоляции ввода 750 кВ путем измерения ЧР с интенсивностью порядка 10-6 Кл. При испытаниях герметичных вводов 110 кВ были измерены разряды интенсивностью (3-6)10-6 Кл.

Однако необходимо отметить следующие недостатки метода непосредственного измерения ЧР.

Во-первых, сигналы ЧР очень слабые, а при измерениях ЧР в условиях эксплуатации источники помех, как правило, не могут быть устранены, тем самым сильно искажают результаты измерений. Ниже 15 кГц возможно влияние высших гармоник промышленной и комбинированных частот, выше 2 МГц возможно снижение уровня сигнала от ЧР. В верхней области частот часто имеются помехи от мощных местных радиопередатчиков. Основным источником неустранимых помех при измерениях в эксплуатационных условиях являются коронные разряды на проводах, арматуре и оборудовании (основной уровень помех). В распределительных устройствах иногда наблюдается также высокий дополнительный уровень помех, который обычно является следствием ЧР, внешних по отношению к контролируемому вводу. К ним относятся разряды между шинами и головками проходных трансформаторов тока при отсутствии между ними перемычки, разряды между элементами токопроводов блочных трансформаторов, разряды на заостренных краях арматуры или на концах ножей отключенных разъединителей и т.п.

Во-вторых, если за период времени, равный продолжительности реакции устройства, на входе измерительного элемента 1 (рис.3.6) будут действовать несколько импульсов, то их энергия суммируется и на выходе появится один эквивалентный импульс. При этом не только будет утеряна такая характеристика, как количество разрядов (или их средняя частота следования), но и будут искажены данные о заряде импульсов, ибо на выходе усилителя амплитуда импульсов будет иметь случайное значение, зависящее не только от заряда, но и от интервала между импульсами.

10. Оценка технического состояния трансформаторных вводов на основе теории нечетких множеств

Традиционным методам диагностики трансформаторных вводов присущи многие недостатки, так как они не учитывают существующую неполноту и нечеткость информации о состоянии ввода.

Анализ опытов диагностики вводов обнаруживает, что для большинства случаев существует такое решающее правило "если параметр X не выше нормы Xн, то ..., а если параметр X выше нормы Xн , то ... ". Норма Xн обычно определяется из результата статистической обработки аварийных состояний по параметру X, однако чем больше мощность и класс напряжения силового трансформатора, тем скуднее такая статистика в связи с трудностью ее реализации. Это привело к тому, что норма Xн - нечеткая величина и нет резких границ или жестких граней, отделяющих одно состояние ввода от другого.

Другое непростое обстоятельство при диагностике вводов по традиционным методам заключается в том, что количество контролируемых параметров значительно, поэтому количество комбинаций этих параметров (при разных условиях) довольно большое. А опыт диагностики в виде решающих правил "если параметр X1..., параметр X2..., ..., параметр Xп..., то ..." не охватывает полный набор комбинаций параметров. Отсюда возникает вопрос: если комбинация параметров находится вне существующих решающих правил, то какое решение необходимо принимать в таком случае ?

Для преодоления этих трудностей целесообразно применять теорию нечетких множеств Заде [9]. Эта теория была предложена для количественного анализа таких гуманистических систем как лингвистика, экономика, политика. Однако в последнее время отмечено применение теории нечетких множеств в технике для решения задач проектирования и управления в медицинской диагностике.

В теории канторовских множеств произвольное подмножество А универсального множества U однозначно определяется своим индикатором:

Заде расширил класс подмножества U, введя понятие нечеткого, "расплывчатого" множества. Нечеткому множеству соответствует обобщенный индикатор (числовые функция со значениями из всего отрезка [0,1]), получивший название функции принадлежности µ(x). Тогда справедливо следующее определение: нечеткое множе-ство содержит элементы с функциями принадлежности, принимаю-щими любые значения на интервале [0,1], хотя бы одна из которых отлична от единицы.

Ниже приведем основные определения и операции с нечеткими множествами, которые необходимые нам при решении задачи диагностики трансформаторных вводов [8].

Объединением нечетких множеств А и В в пространстве X называется нечеткое множество А и В с функцией принадлежности, определяемой соотношением вида

,

или, в дизъюнктивной форме,

,

где символ ”V “ означает максимум.

Пересечением нечетких множеств А и В в пространстве V называется нечеткое множество А ^ В с функцией принадлежности, определяемой соотношением вида

,

или, в конъюнктивной форме,

,

где символ “ ^ ” означает минимум.

Согласно общего алгоритма [10] методика диагностики вводов на основе теории нечетких множества реализуется в несколько этапов.

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

На втором этапе необходимо задавать вид функции принадлежности нечетких термов при разных контролируемых параметрах по мнениям опытного персонала по диагностике вводов.

На третьем этапе необходимо зафиксировать значения измеренных параметров состояния ввода

На четвертом этапе определяем функции принадлежности нечетких термов при фиксированных значениях измеренных параметров .

На пятом этапе вычисляем значения многопараметрических функций принадлежности при векторе состояния для всех диагнозов d пользуясь логическими уравнениями из базы знаний. При этом логические операции И (^) и ИЛИ (V) над функциями принадлежности заменяются на операции min и max.

На шестом этапе определяем решение d, для которого

Этому решению и будет соответствовать искомый диагноз ввода с вектором параметров состояния

Аппаратная реализация этой методики реализуется с помощью элементов вычислительной техники. При этом вычислительная часть легко реализуется путем выполнения операций min и mах. Для изложения базы знаний будем обозначать: tgδ1 - тангенс угла диэлектрических потерь основной изоляции (%); tgδ3 - тангенс угла диэлектрических потерь наружных слоев изоляции (%); tgδм↓70C - тангенс угла диэлектрических потерь масла, измеренный при подъеме температуры в точке 70°С (%); tgδм↑70C -тангенс угла диэлектрических потерь масла, измеренный при снижении температуры в точке 70°С, после предварительной выдержки масла при температуре 90-100°С в течение 15-20 минут (%); Ai- значение концентрации растворенного в масле 1 -го газа или суммы концентраций газов (% об. 10-4= мкл / л), где i - Н2, С2Н2, ∑СХНУ = СН4 + С2Н2 + С2Н4 + С2H6; Vi - абсолютная скорость нарастания концентрации i-го газа или суммы газов (% об. 10-4 / мес),

где Аiп- значение концентрации i -го газа при предыдущем измерении (все обозначения нижним индексом с буквой п означают предыдущее измерение), Т - промежуток времени с момента предыдущего измерения (мес.); V∑- относительная скорость нарастания концентрации суммы всех газов (% / мес),

где Аiн- начальное (безопасное) значение концентрации i -го газа (% об. 10 ) определяется из табл. 7, N - граничное значение контролируемого параметра (Ai, tgδ, V, V∑ и т.д.) определяется из табл.3.5-3.7.

На основании [11] приведем правила принятия решений, при этом мы постоянно будем обращаться к таблицам 3.2-3.6.

Правило 1: Вводы эксплуатируются с обычно принятой периодичностью контроля при выполнении следующих условий:

тангенс угла диэлектрических потерь основной изоляции tgδ1 удовлетворяют условиям табл. 3.5;

-значение tgδ1 возросло не более чем в 1.2 раза;

-тангенс угла диэлектрических потерь наружных слоев изоляции tgδ3 не превышает значений, приведенных в табл. 3.5;

-значение tgδ3 возросло не более чем в 1.5 раза по сравнению с предыдущими измерениями;

-значение tgδ3 возросло не более чем в 2 раза по сравнению с данными при вводе в эксплуатацию;

-тангенс угла диэлектрических потерь масла на спаде температуры (tgδм↓70C) не превышает значений, приведенных в табл. 3.5 (если такие измерения проводились);

-концентрация содержания любого из газов Н2, С2Н2, СХНУ не превышает значений, приведенных в табл.3.5, а соотношение АCO/АCO2<0.2;

-относительная скорость увеличения концентрации суммы горючих газов Н2 + ∑СХНУ в месяц (V∑).

V∑< 10 (%/мес.) при периодичности контроля через 1 год,

V∑< 4 (%/мес.) при периодичности контроля через 3 года.

Эти условия в наглядном виде приведены в таблице 3.2.

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

-имеется tgδ3 > 1.5 N, по табл. 3.5;

-значение tgδ1 увеличивалась более чем в 1.2 раза по сравнению с предыдущим измерением, а tgδ3 превышает значения, приведенные в табл. 3.5;

- концентрация растворенных газов превышает значения по табл. 3.8, а тангенс наружных слоев изоляции tgδ3 превышает значения,приведенные в табл. 3.5;

- концентрация газов Н2, С2Н2 и ∑СХНУ превышают или равны значениям, приведенным в табл. 3.5 одновременно;

-абсолютная скорость нарастания любого из газов Н2, С2Н2, и ∑СХНУ превышает значения, приведенные в табл. 3.6;

- тангенс угла диэлектрических потерь масла, измеренный при подъеме и на спаде температуры превышают значения, приведенные в табл. 3.7 и 3.5 соответственно (если проводился отбор пробы масла);

- концентрация любого из газов Н2, С2Н2 и ∑СХНУ превышает в 2 и более раза значения, приведенные в табл. 3.5;

- величина tgδ3 превышает значения, приведенные в табл. 3.5 и значения tgδ3п при предыдущем измерении, при этом относительная скорость нарастания суммы горючих газов в месяц Vг превышает значений по табл. 3.2.

Эти условия в наглядном виде приведены в таблице 3.3. Отмечаем особый случай, когда tgδ1<0. При этом если tg1 превышает значения, приведенные в табл.4 то считают, что ввод подлежит отбраковке, а если tgδ3 меньше нормы по табл. 3.5 то необходимо убедиться в истинности результатов измерений. Для этого необходимо предпринять меры по исключению влияния погодных условий и загрязнения внешней изоляции.

Правило 3: Вводы с предполагаемым наличием дефекта, у которых значение контролируемых параметров отклоняются от нормы, требуют уточнения диагноза эксплуатационного состояния. Для этого необходимо увеличить объем испытаний, провести дополнительные измерения или повторный контроль при сокращенной периодичности. Правило формируется в виде "если ... то ...".

-если значение tgδ1 лежит в пределах величин, приведенных в табл. 3.5, а tgδ3 превышает указанное в табл. 3.5, но не более чем в 1.5 раза, а концентрации газов не превышают значений табл. 3.8, то необходимо произвести отбор пробы масла для измерения tgδм↑70C и tgδм↓70C измеренные при подъеме и спаде температуры, если эти тангенсы не превышают значений в табл. 3.7 и 3.5 соответственно, то допустима проведение повторного контроля через год;

-  если значение tgδ1 лежит в пределах величин, приведенных в табл. 3.5, а tgδ3 превышает указанное в табл. 3.5, но не более чем в 1.5 раза, а концентрация газов превышает значения по табл. 3.8, но не более граничных значений по табл. 3.5, а относительная скорость нарастания суммы горючих газов в месяц V1 < N (табл. 3.2) то допустимо проведение повторного контроля через 0.5 года;

-  если концентрации растворенных газов достигают граничных значений только за счет содержания Н2 и СH4, но не более удвоенного значения по табл. 3.5 и при этом АCO/АCO2 < 0.15, то допустимо проведение повторного контроля через 0.5 года;

-если установлено, что относительная скорость нарастания суммы горючих газов в месяц V∑ не превышает значений по таб. 3.2, а концентрация каждого из газов Н2, С2Н2 и ∑СХНУ не превышает значений по табл. 3.5, или концентрация газов превышает граничное значение, но не более удвоенного граничного значения по табл. 3.5 только по отдельным газам, то необходимо произвести повторный контроль не позднее, чем через время

где N - граничные значения концентраций газов по табл. 3.5; Аi ,АiП - концентрации газов, полученные при последнем и предыдущем измерении соответственно; Т - период между последним и предыдущим измерением; Тсл - минимальное значение из рассчитанных по газам Н2, ∑СХHУ по выражению в правой части выше приведенного неравенства.

Эти условия в наглядном виде приведены в таблице 3.4.

В соответствии с базой знаний (табл.3.2-3.4) будем определять следующие технические состояния:

Таблица 3.9

Исходя из базы знаний, целесообразно ввести следующие входные параметры с соответствующими возможными диапазонами изменения [11]:

Таблица 3.10

В случае, когда при измерении x1 <0, то необходимо проверить тщательно результаты других измерений и повторно производить измерение параметра x1. Если это подтверждается, то ввод подлежит отбраковке.

Задача диагностики состоит в том, чтобы каждому сочетанию значений факторов поставить в соответствие одно из решений dm.

Параметры x1-x18, определенные выше, будем рассматривать как лингвистические переменные. Кроме того, введем еще одну лингвистическую переменную: d - опасность повреждения ввода, которая измеряется уровнями d1 - d5.

Для оценки значений лингвистических переменных x1, x2, x5, x6, x17, x11, x12, x13, x14, x15, x16, x17 будем использовать два терма: Н - низкий, В – высокий. Для оценки значений лингвистических переменных x3, x4, x8, x9, x10, x18 будем использовать три терма: Н - низкий, С - средний, В – высокий. Каждый из этих термов задает нечеткое ограничение на множество, заданное с помощью соответствующей функции принадлежности.

Предполагаем, что функции принадлежности параметров x1,, x2, x5, x6, x17, x11, x12, x13, x14, x15, x16, x17 имеют одинаковый вид для каждого терма Н или В.(рис. 3.7.); функции принадлежности параметров x3, x4, x8, x9, x10, x18 имеют одинаковый вид для каждого терма Н, С или В (рис. .8.).

Из таблиц 3.2 – 3.4 формулируем следующие нечеткие высказывания:

1)ЕСЛИ (x1= Н) и (x2 = Н) и (x1- Н) и [(x1= Н) или (x4 = С)] и

(x5 = Н) и (x7 = Н) и (x8 = Н) и (x9 = Н) и (x10 = И) и (x11=Н) и [(x18= Н) или (x18 = С)],

то d = d1

2)ЕСЛИ [x3 = В),

или [(x3 = С) и (x2 = В)],

или {( x3 = С) и [(x1= В) или (x13 = В) или (x14=B)]}э

или [(x8 = С) и (x9 = С) и (x10 = С)],

или (x15 = В),

или (x16 = В),

или (x17 = В),

или [(x6 = В) и (x7 = В)],

или (x8 = В),

или (x9 = В),

или (x11 = В),

или {(x3 = С) и [(x4 = С) или (x4 = В)] и (x11 =В)},

то d = d2

3)ЕСЛИ [(x3=С) и (x1= Н) и (x12 = Н) и (x13 = Н) и (x14 = Н) и(x6 = Н) и (x7 = Н)],

то d = d3

4)ЕСЛИ [(x3 = С) и (x1 = Н) и (x12 = В) и (x13 = В) и (x14 = В) и(x8 = Н) и (x9 = Н) и (x10 = Н) и (x11 = Н)] или [(x18 = Н) и (x8=С) и (x10 = С) и (x13 = Н)],

то d = d4

5)ЕСЛИ [(x11 = В) и (x8 = Н) и (x9 = Н) и (x10 - Н)], или {( x11 =В) и [(x8 = С) или (x9 = С) или (x10=С)]},

то d = d5

Пользуясь функциями принадлежности, запишем эти логические высказывания в виде логических уравнений. При этом заменяем слово "и" операцией “^” (для краткости будем использовать знак "•", слово "или" операцией “V”.

Согласно общего алгоритма [10], решению задачи диагностики соответствует тот диагноз, который имеет максимальное значение функции принадлежности:

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

Отметим, что из выше приведенных правил ЕСЛИ...ТО... можно получить однопарамстрические правила:

ЕСЛИ (x3=В), то d = d2;

ЕСЛИ (x8 = В), то d = d2;

ЕСЛИ (x9 = В), то d = d2;

ЕСЛИ (x10 = В), то d = d2;

ЕСЛИ (x15=В), то d = d2;

ЕСЛИ (x16 = В), то d = d2;

ЕСЛИ (x17 = В), то d = d2;

(при этом необходимо учитывать, что параметр xз измеряется по мостовой схеме, а параметры x8, x9, x10, x15, x16, x17 измеряются по методу хроматографического анализа);

двухпараметрические правила:

ЕСЛИ [(x3=С) и (x2 = В)], то d = d2;

ЕСЛИ [(x6 =В) и (x7= В)], то d = d2;

трехпараметрические правила:

ЕСЛИ [(x8 = С) и (x9 = С) и (x10 = С)], то d = d2;

ЕСЛИ {(xз = С) и [(x4=С) или (x4 = В)] и (x11 = В)}, то d = d2;

четырехпараметрические правила:

ЕСЛИ {(x3 =С) и [(x12=В) или (x13 = В) или (x14 = В)]},

то d = d2

ЕСЛИ [(x18 = Н) и (x8 = С) и (x10 = С) и (x13 = Н)]

то d = d4;

ЕСЛИ [(x11 = В) и (x8 = Н) и (x9 = Н) и (x10 = Н)],

то d = d5;

ЕСЛИ {( x11=В) и [(x8=С) или (x9 = С) или (x10 = С)]},

то d = d5;

семипараметрическое правило:

ЕСЛИ [(x3 = С) и (x1 = Н) и (x12 = Н) и (x13 = Н) и (x14 = Н) и

(x6 = Н) и (x7 = Н)], то d = d3;

девятипараметрическое правило: ЕСЛИ [(x3=С) и (x1= Н) и (x12 =В) и (x13= В) и (x14 = В) и (x8 = Н) и (x9 = Н) и (x0 = Н) и (x11 = Н)],

то d = d4

одиннадцатипараметрическое правило:

ЕСЛИ (x1 = Н) и (x2 = Н) и (x3 = Н) и [(x4 = Н) или (x4 = С)] и (x5 = Н) и (x7= Н) и (x8 = Н) и (x9 = Н) и (x10 = Н) и (x11 = Н) и [(x18 = Н) или

(x18 = С)],

то d = d1;

Отсюда следует целесообразность контроля в первую очередь по однопараметрическим правилам сначала параметра x3. Если x3 = В (при этом будем условно говорить, что параметр x3 больше принадлежит терму В, то сразу принимаем решение d = d2, иначе проверяем один из параметров x8 , x9 , x10, x15, x16, x17. Если один из этих параметров больше принадлежит терму В, то также сразу принимаем решение d2, иначе проверяем по двухпараметрическим правилам путем дополнительного рассмотрения параметра x2 (учитываем, что параметр x2 измеряется одновременно с параметром x3 при контроле на подстанции по методу измерения tg). Если условия по двухпараметрическим правилам не выполняются, то переходим к трехпараметрическим правилам и т.д.

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

Традиционная диагностика по правилам 1, 2, 3 является частным случаем предложенной выше методики с применением теории нечетких множеств.

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

Ниже приведем один практический пример, при котором правила традиционной диагностики не позволяют принять решение, а основанные на теории нечетких множеств - позволяют.

Данные измерений для диагностики ввода 110 кВ с маслом типа ГК приведены в таблицах 3.11 и 3.12.

Подставляя данные из таблиц 3.11 и 3.12 в формулы для x с учетом таблиц 3.5 – 3.8, получим значения параметров x. Проверка по одно-, двух-, трех-, четырех-, семи-, девяти- и одиннадцатипараметрическому правилу не дает ответа на диагноз. Это значит, что традиционная диагностика неприемлема в этой ситуации.

Подставляя в формулы найденные значения x, получим значения функций принадлежности. Все эти величины приведены в таблице 3.13.

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

Отсюда следует, что среди пяти диагнозов максимальное значение функции принадлежности имеет диагноз d2. Следовательно, ввод подлежит немедленной отбраковке.

Таблица 3.11 Данные измерений  (%)

Таблица 3.12 Данные измерений по методу ХАРГ (% об. )

Список литературы

1. Рыжкин В.Я. Тепловые электрические станции. - М.: Энергия, 1992.-276с.

2. Тепловые и атомные станции. Книга 3. Справочник. Под ред. В.А. Григорьева, В.М. Зорина. 2-е издание, переработанное –М.: «Энергия», 1989, - 600 с.

3. Рыжкова Л.Д., Козулин В.С. Электрооборудование станций и подстанций. Учебник для техникумов. – 2-е изд., перераб. - М.: «Энергия», 1980, - 600 с.

4. Электрическая часть электрических станций и подстанций. Под ред. А.А. Васильева. Учебник для вузов - М., «Энергия», 1980. - 608 с.

5. Электрическая часть электрических станций и подстанций. Справочные материалы для курсового и дипломного проектирования. Под ред. Б.Н. Неклепаева. Изд.2 -е , перер. М., «Энергия», 1972. - 336 с.

6. Околович М. Н. Проектирование электрических станций: Учебник для вузов. - М.: Энергоатомиздат, 1982. - 400 с.

7. Бажанов С.А., Воскресенский В.Ф. Монтаж и эксплуатация маслонаполненных вводов.-М.:Энергоатомиздат, 1981.-104с.

8. Штовба С.Д. Введение в теорию нечетких множеств и нечеткую логику.

9. Заде Л.А. Понятие лингвистической переменной и ее применение к принятию приблеженных решений. М., 1976.

10. Ронштейн А.П. Медицинская диагнстика на нечеткой логике. – Винница: Континент – ПРИМ, 1996. – 132с.

11. Галузевий керівний документ. Маслонаповнені вводи напругою 110 – 750 кВ. Типова інструкція з експлуатації. – Київ, 2004

12. Методические указания по диагностике состояния изоляции высоковольтных вводов 110 – 750 кВ. М.,1990.

13. Методичні вказівки з техніко – економічного обгрунтування інвестиційних проектів електричних станцій. Уклад.: Є.Г. Скловська, К.Г. Тодорович. – К.: Політехніка, 2002. – 24с.

14. Охрана труда в єлектроустановках: Учебник для вузов / Под ред. Б.А. Князевского. – 3-е изд., перераб. И доп. – М.: Єнергоатомиздат, 1983. – 336 с.

Размещено на

Страницы: 1, 2, 3