Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
18
Вінницький національний технічний університет
УДК 519.876.5:621.311
МАТЕМАТИЧНІ МОДЕЛІ ДЕТЕРМІНІЗАЦІЇ ПРОЦЕСІВ В СИСТЕМАХ ЕЛЕКТРОПОСТАЧАННЯ
01.05.02 “Математичне моделювання та обчислювальні методи”
АВТОРЕФЕРАТ
дисертації на здобуття наукового ступеня
кандидата технічних наук
Вінниця
Дисертацією є рукопис.
Робота виконана у Вінницькому національному технічному університеті Міністерства освіти і науки України.
Науковий керівник: доктор технічних наук, професор
Мокін Борис Іванович,
Вінницький національний технічний університет,
завідувач кафедри електромеханічних систем
автоматизації
Офіційні опоненти: доктор технічних наук, професор
Квєтний Роман Наумович,
Вінницький національний технічний університет,
завідувач кафедри автоматики та інформаційно-
вимірювальної техніки
кандидат технічних наук, доцент
Нагул Володимир Іванович
Південно-Західна електроенергетична система,
відособлена структурна одиниця Державного
підприємства Національна енергетична компанія
“Укренерго”, м. Вінниця
начальник виробничо-технічної служби
Провідна установа: Державний науково-дослідний інститут інформаційної
інфраструктури Державного комітету звязку та
інформатизації і НАН України, м. Львів,
відділ інформаційних технологій
Захист відбудеться 16.06.2004 р. о 14 годині на засіданні спеціалізованої вченої ради Д 05.052.01 у Вінницькому національному технічному університеті за адресою: 21021, м. Вінниця, Хмельницьке шосе, 95.
З дисертацією можна ознайомитись у бібліотеці Вінницького національного технічного університету за адресою: 21021, м. Вінниця, Хмельницьке шосе, 95.
Автореферат розісланий 7.05.2004 р.
Вчений секретар
cпеціалізованої вченої ради Захарченко С. М.
ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ
Успішне вирішення задачі поточного диспетчерського контролю електроспоживання є необхідною умовою проведення діючої політики енергозбереження, раціонального керування електроспоживанням, контролю за виконанням графіків обмежень електроспоживання, тощо. Розвязання цієї задачі дозволить виявити вузли мережі з позаобліковим електроспоживанням, що забезпечить успішну боротьбу з комерційними втратами електроенергії.
Традиційний шлях розвязання вищеперерахованих задач потребує оснащення телеметричними джерелами інформації переважної більшості трансформаторних підстанцій електромережі. Це означає, що на кожній із цих підстанцій кожне приєднання 110-35-10-6 кВ повинно бути обладнано сенсорами телесигналізації, напруги і потужності, кожна підстанція повинна мати апаратуру прийому, обробки та передачі телеметричної інформації, а також надійний канал зв'язку з оперативно-диспетчерською службою підприємства електромереж.
Однак сьогодні у вітчизняних мережах телевимірюваннями охоплена незначна частина підстанцій, тому реалізація апаратного шляху досягнення повного диспетчерського контролю за електроспоживанням вимагатиме багатьох років роботи і значних додаткових капіталовкладень.
В існуючих оперативно-інформаційних комплексах (ОІК), як альтернативу відсутнім телевимірам, часто використовують введені оператором параметри електроспоживання, що базуються на попередніх статистичних даних. Але, зважаючи на те, що принципово неможливо підібрати вручну параметри режиму так, щоб вони в кожний момент часу були сумісні за законами Кірхгофа з параметрами, що безперервно змінюються, отриманими в результаті телевимірів, приходимо до висновку, що застосування ручних введень параметрів електроспоживання не може бути повноцінною альтернативою відсутнім телевимірам.
Очевидно, що найбільш доцільною альтернативою відсутнім телевимірам є відтворення параметрів електроспоживання в усіх вузлах мережі, у яких відсутні телевиміри, у темпі процесів, що протікають у цих вузлах, за допомогою математичних моделей. Тому, розробка математичних моделей та методів, які дозволяють розвязувати подібні задачі є актуальною.
В межах цього дослідження розвязувались такі задачі:
Одержані наукові результати впроваджені в програмному комплексі “Аналітична Система Відтворення Електроспоживання” (АСВЕС), який з 1997р. успішно функціонує в оперативно-диспетчерській службі Вінницьких центральних високовольтних електромереж (ВЦВЕМ). Крім того, проводиться поетапне впровадження АСВЕС в оперативно-диспетчерській службі Вінницьких східних високовольтних електромереж (ВСВЕМ). Впровадження підтверджується відповідним актом.
Можливості АСВЕС дозволяють відтворювати і відображати параметри поточного електроспоживання в усіх вузлах без телевимірів, а також забезпечувати розрахунок і контроль балансів активної потужності й енергії підприємств електромереж, РЕМ і окремих значних споживачів із достатньою точністю навіть за умови дуже обмеженої кількості реальних вимірів.
Слід особливо відзначити, що АСВЕС дає можливість безперервно вести поточні розрахунки технічних втрат в мережі, що дозволяє визначати місця позаоблікового електроспоживання.
У вступі обгрунтовано актуальність теми, визначено мету та задачі досліджень, відзначено наукову новизну та практичну цінність роботи.
В першому розділі розглядаються математичні методи вивчення та аналізу процесів в складних системах, дається огляд відомих математичних моделей для розв'язання задач відтворення параметрів електроспоживання, обгрунтовуються задачі дослідження.
Задача детермінізації процесів в системах є в певному розумінні оберненою задачею по відношенню до задачі ідентифікації. Для такої постановки задачі доцільно користуватися більш універсальним, ніж в традиційній теорії систем, означенням поняття “система”, а саме: “Система це довільна множина, на якій визначена деяка структура”. При цьому структурою С на множині S називають відношення, в якому знаходяться між собою елементи цієї множини. Саме на цій основі і сформована аксіоматика теорії детермінізації процесів в складних системах, яка опублікована в роботі [4].
Вона виходить із постулату:
Для будь-якої технічної системи S існує непустий клас систем , якому вона належить, такий, що:
Введемо для класу декомпозицію DK цього класу, як пару виду DK = (U, P), де U універсум елементів, P= ||pij||, i=1, ….w; j=1, …..n простір можливих властивостей всіх елементів універсуму U, а ||pij|| матриця властивостей елементів універсуму U, w кількість елементів універсуму U, n максимальна кількість властивостей елемента із U.
Виділимо в класі систему S, для якої необхідно розвязати задачу детермінізації, тобто відтворити невідомі параметри. На основі декомпозиції DK цього класу для системи S можна ввести простір моделей системи S на декомпозиції DK класу , як декартів добуток виду = CP, де: C = (E) простір структур можливих моделей системи S.
Носій цього простору Eρ(U) є множиною елементів моделі системи S, де ρ(U) множина всіх підмножин універсуму U; (E) множина бінарних відношень на E таких, що їх транзитивне замикання зв'язне, тобто
(E) = {c | (c E)a,bE ((ab)((a,b)(b,a) ))}, де транзитивне замикання бінарного відношення c(E), яке є можливою структурою моделі системи S.
Кожен елемент eE моделі системи S, крім того, що він має власні параметри та зв'язки з іншими елементами моделі системи, може мати також і канали зв'язку з оточуючим середовищем, які не входять до структури моделі системи S. Через ці канали оточуюче середовище впливає на систему за допомогою вхідних параметрів, а сама система впливає на оточуюче середовище за допомогою вихідних параметрів. Таким чином кожний елемент моделі можна записати як трійку e=(xe ,ye ,pe), де xe = (xe1 , xe2 ,….. xek) вектор входу елементу e, k максимальна кількість вхідних параметрів елементу із E, ye = (ye1 , ye2 ,….. yeg) вектор виходу елементу e, g максимальна кількість вихідних параметрів елементу із E; ,pe e-тий рядок матриці ||pij||.
Законом Z, що діє на просторі моделей ,називається відображення виду Z:XY, де X={x | x=(x,….,xe ,….,xh)} простір можливих вхідних параметрів, Y={y | y=(y,….,ye ,…..,yh)} простір можливих вихідних параметрів, h кількість елементів множини E.
Вектори x,y, матриця ||pij|| та відображення Z залежно від природи причинно-наслідкових зв'язків можуть бути представлені в детермінованій, стохастичній або нечіткій формах.
Якщо відображення Z представлене в стохастичній формі, то, незалежно від форми представлення вхідного вектору x та матриці властивостей ||pij||, вихідний вектор буде мати стохастичний характер. У випадку, якщо вектор x або матриця ||pij|| представлені в нечіткій формі, то на виході відображення Z будемо мати нечіткий вектор y. Акцентуємо увагу на тому, що навіть для фіксованого детермінованого числового аргументу x, його образ при стохастичному відображенні Z буде випадковою величиною, а при нечіткому відображенні нечіткою величиною.
Нечітке відображення Z, незалежно від форми представлення вхідного вектору x та матриці властивостей ||pij||, призведе до нечіткого характеру вихідного вектору. В окремому випадку, якщо вектор x або матриця ||pij|| представлені в стохастичній формі, то на виході відображення Z будемо мати нечіткий випадковий вектор y.
У випадку детермінованих апріорі відомих законів , незалежно від форми представлення векторів x,y та матриці ||pij|| (детермінованій, стохастичній, нечіткій), задача відтворення невідомих координат векторів x, y та невідомих елементів матриці ||pij|| називається детермінізацією системи S. Очевидно, що в просторі детермінованих параметрів детермінований закон Z є бієктивним відображенням або бієктивним звуженням деякого довільного відображення, в просторі стохастичних параметрів характеризуватиме тренд, а в просторі нечітких параметрів буде функцією нечітких аргументів, яка задає “розмитий тренд”.
Таким чином, задачу відтворення параметрів електроспоживання можна сформулювати як окремий випадок задачі детермінізації процесів в складних системах.
В результаті огляду та аналізу відомих математичних моделей для розв'язання задач відтворення параметрів електроспоживання були виявлені деякі невирішені проблеми:
В результаті цих міркувань були обгрунтовані мета та задачі дослідження.
В другому розділі досліджено загальні умови, за якими можливе відтворення стану складних технічних систем, більш детально проаналізовані умови відтворення стану систем з лінійними та лінеаризованими функціональними зв'язками, доведена теорема, яка дозволяє визначити можливість відтворення стану лінійних систем без попередніх обчислень.
Згідно означень, наведених в першому розділі, кожна модель системи S є пара виду =(c,||pij||), для якої має місце співвідношення:
y=Z (x,)=Z (x,(c,||pij||)), (1)
де yY вектор вихідних параметрів, xX вектор вхідних параметрів.
Нехай для деякої сукупності відомих координат векторів x та y, частково відомих елементів матриці властивостей ||pij|| при заданій структурі C існує алгоритм відтворення всіх невідомих координат векторів x та y, а також невідомих елементів матриці властивостей ||pij||. Така сукупність відомих параметрів називається детермінізуючою. Будемо вважати, що кожен елемент eE моделі системи S може бути чи не бути зв'язаним з об'єктом, який називається джерелом інформації. У першому випадку вхід, вихід, власні параметри та зв'язки такого елементу з іншими елементами стають визначеними.
Для кожного елементу eE введемо адитивну функцію інформаційних витрат ω(e,c) з від'ємним інгредієнтом. Якщо елемент е, не зв'язаний з джерелом інформації, то ω(e,c)=0.
Цілком зрозуміло, що може існувати велика кількість варіантів розташування джерел інформації на елементах множини , які дають детермінізуючу комбінацію параметрів для заданої моделі системи S. Нехай існує таких варіантів. Кожен з цих варіантів характеризується множиною елементів E, зв'язаних з джерелом інформації, та функцією інформаційних витрат: Ω() = ∑ω(e,c). Очевидно, що серед цих варіантів можливо вибрати оптимальний , виходячи з критерію: . Якщо функція інформаційних витрат обмежена значенням Ωlim, то може так статись, що {1,2,3,…,}(Ω()>Ωlim), тобто кожний варіант розташування джерел інформації в , який дає детермінізуючу комбінацію параметрів, не відповідає обмеженню. Це буде означати, що однозначне відтворення всіх невідомих координат векторів x та y, а також невідомих елементів матриці властивостей ||pij|| неможливе. У цьому випадку можлива квазідетермінізація, тобто відтворення невідомих параметрів з прийнятною похибкою.
Рівняння (1) можна переписати у вигляді: (yв, yн) = Z((xв, хн),(c, ||pij||)), де yв відомі компоненти вектора вихідних параметрів, yн невідомі компоненти вектора вихідних параметрів, xв відомі компоненти вектора вхідних параметрів, хн невідомі компоненти вектора вхідних параметрів. Якщо у випадку відсутності детермінізуючої сукупності відомих параметрів yв, xв, кожен з векторів yн, хн невідомих параметрів можна розкласти на дві частини на параметри, які можна однозначно відтворити (позначимо їх yн+, хн+) та на параметри, які неможливо однозначно відтворити (позначимо їх yн-, хн-), то матимемо
(yв, yн+, yн-) = Z((xв, хн+, хн-),(c, ||pij||)). (2)
Таким чином задача квазідетермінізації полягає в тому, щоб якимось чином визначити саме параметри yн-, хн- з прийнятною в межах задачі похибкою.
Із загальних міркувань випливають три можливі види квазідетермінізації:
Вона полягає в тому, що дослідник, враховуючи свій власний досвід, інтуїцію та думку експертів, задає деякі евристичні припущення у вигляді детермінованих функціональних співвідношень між параметрами yн-, хн-, які так доповнюють реально діючі закони Z, що деякі доступні при обмеженні Ωlim комбінації параметрів стають детермінізуючими. Ці припущення в процесі дослідження реальної системи можуть надалі уточнюватись, а можлива похибка відтворення зменшуватись.
Цей тип квазідетермінізації характеризується тим, що параметри yн-, хн- (тобто параметри, які неможливо однозначно відтворити при заданому обмеженні Ωlim) представляються у вигляді статистичних оцінок випадкових величин і вираз (2) є звичайним детермінованим відображенням випадкових та детермінованих величин.
Вона полягає в тому, що координати векторів yн- та хн- представляються у вигляді нечітких множин (наприклад, регулярних нечітких чисел). В результаті, замість звичайної системи рівнянь (2) може бути отримана система нечітких рівнянь.
Розглянемо умови однозначної детермінізації для випадку повністю відомих елементів матриці властивостей ||pij||, відомої структури С та лінійних законів Z.
Назвемо елементи системи S з повністю відомими параметрами фіксованими елементами, а елементи з повністю невідомими параметрами вільними елементами.
Умови однозначної детермінізації визначаються сформульованою і доведеною в другому розділі теоремою:
Нехай існує система S, яка має фіксовані та вільні елементи. Нехай також кількість фіксованих та вільних елементів співпадає та більша за 1. В цьому випадку, якщо на графі цієї системи фіксовані та вільні вершини розташовані таким чином, що їх неможливо поділити на пари ф-в, ланцюги між якими попарно не перетинаються, то система рівнянь, що побудована для цього графа несумісна.
Ця теорема дозволяє визначити умови однозначної детермінізації системи S без попередніх обчислень. Для цього цілком достатньо провести топологічний аналіз графу системи S.
В третьому розділі досліджується особливість проблеми неповноти вхідних даних при детермінізації режимів електричних мереж, яка повязана з відсутністю вимірів кутів зсуву фаз між напругами різних вузлів мережі. Для зменшення впливу цієї невизначеності на процес відтворення режиму електроспоживання в просторі комплексних електротехнічних параметрів застосовується лінійна неевклідова метрика.
Розглянемо тепер деякі суттєві особливості забезпечення телевимірами електричних мереж. Як відомо, пристрої для виміру комплексних значень потужності в лініях (тобто її активної та реактивної складової) відносно прості за принципом дії та недорогі, оскільки кут зсуву фаз визначається між струмом та напругою, що виміряються на одній і тій же лінії. Що стосується пристроїв для виміру комплексних значень напруги на шинах підстанцій, то в цьому випадку складність полягає в тому, що кут зсуву фаз необхідно визначати між напругами різних підстанцій. Це потребує таких великих капіталовкладень, що робить практично недоцільними такі виміри. Отже ми маємо можливість одержувати від телеметричних пристроїв тільки модулі напруг у вузлах мережі. А це в свою чергу означає, що при застосуванні теорії детермінізації процесів в складних системах для відтворення режимів електричних мереж ми будемо мати справу не з однозначною детермінізацією, а з квазідетермінізацією. Тому дуже важлива розробка математичних моделей та методів, які дозволяють знизити рівень невизначеності вимірів напруг на шинах підстанцій.
В цьому розділі запропоновано новий спосіб детермінізації цієї невизначеності. Його суть полягає в тому, що в просторі комплексних електротехнічних параметрів додатково вводиться лінійна неевклідова метрика у вигляді так званого лінійного неевклідового модулю вектора напруги U# =β|Ua|+|Up|, де Ua активна складова вектору напруги, Up реактивна складова вектору напруги, коефіцієнт приведення.
Використання лінійного неевклідового модуля вектора напруги регламентується сформульованою і доведеною в цьому розділі теоремою:
Для мереж 110-35 кВ при індуктивному характері навантаження лінійна неевклідова втрата напруги в лінії дорівнює 0 тоді та лише тоді, коли виконується умова R / X = (1 βtg φ)/( β+ tg φ), де Rпитомий активний опір лінії, Xпитомий індуктивний опір лінії, φ кут зсуву фаз між напругою і струмом, β коефіцієнт приведення.
Проаналізуємо межі використання цієї теореми для задачі відтворення параметрів електроспоживання.
Очевидно, що R, X та φ є незалежними від дослідника величинами, а впливати він може тільки на коефіцієнт приведення (β ). Розрахунки показали, що для повітряних ліній для діапазону перерізів від 50 до 240 ммта значень Cosφ від 0.85 до 0.99, β може приймати значення від 0.887 до 2.141.
Таким чином, вибираючи для кожної лінії певного перерізу та з деяким значенням коефіцієнту потужності відповідне значення β , ми можемо забезпечити виконання умови ΔU#=0.
Якщо умова теореми не виконується, то мають місце співвідношення:
β > (X - R tgφ)/(X tgφ + R)ΔU#>0 та β < (X - R tgφ)/(X tgφ + R)ΔU#<0.
В цьому розділі виконано порівняльний аналіз відносних звичайної та неевклідової втрат напруги за однакових режимних умов. Виявилося, що навіть у випадку, якщо дослідник неточно спрогнозує Cosφ в лінії та відповідно неточно визначить коефіцієнт приведення β , то ΔU# і в цьому випадку буде у 6 разів менше за ΔU.
В четвертому розділі запропоновано розвиток теорії нечітких множин за рахунок введення нетрадиційного класу операцій та принципів узагальнення, будується аксіоматика цього класу та формулюються принципи побудови моделей нечіткої квазідетермінізації.
Необхідність введення нового класу операцій та принципів узагальнення в теорію нечітких множин пояснюється тим, що використання максимінного принципу узагальнення (ММПУ) вимагає повної інформації про функції належності нечітко визначених параметрів задачі, а це для певного кола задач є практично нездійсненною процедурою. Однією з таких задач є нечітка квазідетермінізація процесів в системах електропостачання. Для такої системи навіть найдосвідченіший експерт може визначити для невідомих нечітких параметрів системи лише їх носії або окремі множини -рівня. Побудова на цій основі повних функцій належності невідомих нечітких параметрів є справою ризикованою та ненадійною.
Новий клас (ми будемо надалі називати ці операції та принципи узагальнення -слабкими) є менш вимогливим до повноти даних про функції належності нечітких множин і дозволяє формувати модель нечіткої квазідетермінізації навіть якщо відомі лише окремі множини -рівня.
Сформулюємо більш загальний ніж ММПУ та менш вимогливий до повноти даних про функції належності принцип -слабкого узагальнення для чітких відображень нечітких множин.
Нечітка множина f() в Y є -слабким образом нечіткої множини в X при чіткому відображенні f:XY за означенням тоді та лише тоді, коли yY(f(A)(y) xf-1(y)(A(x))), де (0,1], f-1(y) прообраз елемента уY при чіткому відображенні f:XY.
Аналізуючи запропонований принцип -слабкого узагальнення, легко побачити, що він визначає не один образ нечіткої множини, а ціле сімейство можливих образів, серед яких утримуються і образи, одержані за допомогою ММПУ. Цей принцип дозволяє сформувати образ нечіткого параметра в умовах неповністю заданої функції належності цього параметра.
Якщо у виразі (2) вважати апріорі відомими структуру С та матрицю властивостей ||pij||, його можна записати у вигляді: (yв, yн+, yн-) = F((xв, хн+, хн-), де F функціональна залежність між векторами вхідних та вихідних даних, яка враховує структуру системи та матрицю властивостей її елементів.
Зрозуміло, що в залежності від характеру функції F, отриманий вираз може бути представлений як система рівнянь (лінійних, квадратних, кубічних, поліноміальних, тригонометричних, логарифмічних, тощо) і конкретні алгоритми розвязання такої системи будуть значно відрізнятися один від одного. Тому сформулюємо основні принципи побудови моделей нечіткої квазідетермінізації.
Очевидно, що відомі параметри (yв, xв) та невідомі параметри (yн+, хн+,), які можна однозначно відтворити, немає ніякого сенсу представляти у вигляді нечітких множин. Це викликало б значне ускладнення системи рівнянь. Нечіткими множинами слід представляти лише невідомі параметри (yн-, хн-), які неможливо однозначно відтворити.
Для того, щоб була можливість застосувати принципи -слабкого узагальнення, необхідно отримати від експертів деякий спектр множин -рівнів (,,,….,n) для всіх координат одного з векторів невідомих параметрів (наприклад, вихідного yн-). Далі для кожного з -рівнів складається система нечітких рівнянь. При цьому, якщо задані нечіткі множини будуть визначені на числовій прямій R, то це будуть інтервальні рівняння, які розвязуються набагато легше.
В результаті буде отримано аналогічний спектр множин -рівнів (,,,….,n) для всіх координат вхідного (хн-) вектора невідомих параметрів.
Дуже важливим є те, що досить часто експерт не може визначитись з формами функцій належності нечітких координат одного вектора, але може значно легше це зробити за отриманим спектром множин -рівнів іншого вектора. Саме в цих випадках принципи -слабкого узагальнення практично не мають альтернативи.
В пятому розділі синтезується алгоритм розташування джерел інформації у вузлах електричної мережі з урахуванням умов відтворення параметрів режиму, розробляються математичні моделі евристичної та нечіткої квазідетермінізації та синтезуються відповідні алгоритми, аналізується ефективність програмної реалізації задачі відтворення режиму електричної мережі.
В результаті розгляду алгоритму розташування джерел інформації у вузлах електричної мережі за умови, що фінансових обмежень на кількість джерел інформації не існує, встановлено, що джерела інформації доцільно у першу чергу встановлювати у тих вузлах мережі (вершинах графу), які мають найбільший степінь вершини (тобто кількість ребер, що відходять від цієї вершини), а також на всіх лініях, що відходять від цих вузлів. Оскільки, в реальних електричних мережах найбільший степінь вершини мають вузли живлення (тобто шини живильних підстанцій), тому на першому кроку алгоритму, телеметричні пристрої встановлюються в усіх вузлах живлення мережі. Далі проводиться визначення фіксованих та вільних вузлів і за теоремою, яка доведена в другому розділі, перевіряються умови однозначної детермінізації. Аналогічні кроки повторюються до тих пір, доки кількість невідомих параметрів не зрівняється з кількістю рівнянь за методом вузлових потенціалів, які потрібні для відтворення параметрів режиму.
Зрозуміло, що у випадку існування фінансових обмежень на кількість джерел інформації, цілком можливо, що кількість невідомих параметрів залишиться більшою за кількість рівнянь за методом вузлових потенціалів, і ми будемо мати справу з квазідетермінізацією. Крім того, внаслідок суттєвої невизначеності вхідних даних, яка розглядалась в третьому розділі, нам доведеться будувати модель саме квазідетермінізації навіть у випадку відсутності фінансових обмежень.
Спочатку зупинимося на моделі евристичної квазідетермінізації. Вона базується на двох евристичних припущеннях.
Перше припущення, яке повинно усунути невизначеність, що виникає внаслідок фінансових обмежень на кількість джерел інформації, полягає в тому, що для всієї підмножини вільних вузлів мережі, параметри яких не можуть бути однозначно відтворені, в кожний момент часу співвідношення між фактичними значеннями повних навантажень цих вузлів повинні дорівнювати співвідношенням між типовими значеннями повних навантажень цих вузлів. Типові значення повних навантажень вузлів визначаються з типових графіків навантажень всіх вузлів мережі, які формуються на основі вимірів навантажень для характерних діб кожного сезону.
Інакше кажучи: для довільної пари вільних вузлів в будь-який момент часу t має місце співвідношення Sфі(t)/Smі(t)= Sфj(t)/Smj(t) , де Sфі(t), Sфj(t) фактичні значення повного навантаження вузлів для моменту часу t; Smі(t), Smj(t) типові значення повного навантаження вузлів для моменту часу t.
Друге припущення, яке повинно усунути невизначеність, що виникає внаслідок відсутністю вимірів куту зсуву фаз напругами різних підстанцій, полягає в тому, що для довільного вузла мережі лінійний неевклідовий модуль вектора напруги, розрахований при середньостатистичних навантаженнях в мережі, буде дорівнювати лінійному неевклідовому модулю напруги цього же вузла, розрахованому при фактичних навантаженнях. Звичайно при цьому ми будемо вважати, що Cosφ в мережі майже не змінюються. Завдяки цьому припущенню ми маємо можливість за відомими U# та U знайти =Ua - jUp, де Ua = (βU#+√((1+ β)U- ( U#)))/(1+ β), Up = (U# - β√((1+ β)U-( U#)))/(1+ β).
Базуючись на висунутих вище евристичних припущеннях, в цьому розділі синтезовано алгоритм евристичної квазідетермінізації для задачі відтворення режиму електричної мережі.
В пятому розділі також синтезовано модель та алгоритм нечіткої квазідетермінізації, що базуються на загальних принципах, викладених в четвертому розділі. При цьому слід мати на увазі, що для усунення невизначеності, повязаної з відсутністю вимірів куту зсуву фаз напругами різних підстанцій, доцільно як і раніше застосовувати друге евристичне припущення, тому ми будемо мати модель комбінованої евристично-нечіткої квазідетермінізації.
Моделі та алгоритми, які наведені в цьому розділі, лягли в основу аналітичної системи відтворення електроспоживання (АСВЕС), що була введена в дослідну експлуатацію в листопаді 1997 року в оперативно-диспетчерській службі Вінницьких центральних високовольтних електромереж (ВЦВЕМ). АСВЕС реалізована для операційного середовища Windows 98, 2000 в системі програмування Visual C++.
Можливості системи відтворення режиму електроспоживання дозволяють: відтворювати та відображати параметри поточного електроспоживання вузлів, не охоплених телевимірами; на основі виміряних та відтворених параметрів електроспоживання вузлів мережі, забезпечувати розрахунок інтегральних показників електроспоживання ВЦВЕМ, а саме балансів активної потужності суміжних з ВЦВЕМ підприємств електромереж, власних районів електромереж і окремих значних споживачів; здійснювати контроль за виконанням обмежень на електроспоживання; тощо. Дуже важливим є те, що АСВЕС дає можливість безперервно вести поточні розрахунки технічних втрат в мережі, що дозволяє визначати місця позаоблікового електроспоживання.
На протязі всього часу дослідної експлуатації АСВЕС, як фахівцями ВЦВЕМ, так і автором проводилась постійна робота по перевірці похибки відтворення режиму мережі, тобто адекватності алгоритму евристичної квазідетермінізації. Зазначимо, що в інформаційно-вимірювальну систему (ІВС), яка постачає вхідними даними АСВЕС, на протязі дослідної експлуатації постійно додавалися нові телеметричні пристрої, тому максимальна похибка відтворення режиму мережі відповідно знижувалась.
Графік залежності максимальної похибки інтегральних показників електроспоживання від відносної кількості трансформаторних підстанцій, обладнаних телеметричними пристроями, зображено на рисунку.
Виявлено, що максимальною складовою похибки відтворення режиму мережі є похибка, яка утворюється внаслідок похибок ІВС. Ця складова “відповідає” приблизно за 60% загальної похибки. Наступною за величиною складовою є похибка, яка утворюється внаслідок порушення евристичного припущення щодо пропорційності навантажень у вузлах їх типовим навантаженням. Ця складова “відповідає” приблизно за 30% загальної похибки. Нарешті, останньою складовою є похибка, яка утворюється внаслідок порушення евристичного припущення щодо стабільності значень Cosφ в лініях. Ця складова “відповідає” приблизно за 10% загальної похибки.
Встановлено, що для суттєвого зниження загальної похибки відтворення режиму мережі необхідні такі заходи:
АНОТАЦІЯ
Кацив С.Ш. Математичні моделі детермінізації процесів в системах електропостачання. Рукопис.
Дисертація на здобуття наукового ступеня кандидата технічних наук за спеціальністю 01.05.02 “Математичне моделювання та обчислювальні методи”. Вінницький національний технічний університет, Вінниця, 2004.
Дисертація присвячена новим підходам в математичному моделюванні задач відтворення параметрів режимів електричних мереж в умовах недостатньої кількості телеметричних вимірів. Визначені умови розвязку таких задач. Завдяки використанню в просторі комплексних електротехнічних параметрів лінійної неевклідової метрики зменшена невизначеність вхідних даних таких задач. В теорії нечітких множин введено нетрадиційний клас слабких операцій та принципів узагальнення. Розроблені моделі та алгоритми евристичної та нечіткої квазідетермінізації.
Ключові слова: детермінізація, квазідетермінізація, лінійна неевклідова метрика, нечіткі множини, слабкі операції, принципи узагальнення, процеси, математичні моделі, системи електропостачання.
Кацыв С.Ш. Математические модели детерминизации процессов в системах электроснабжения. - Рукопись.
Диссертация на соискание ученой степени кандидата технических наук по специальности 01.05.02 - “Математическое моделирование и вычислительные методы”. - Винницкий национальный технический университет, Винница, 2004.
Диссертация посвящена новым подходам в математическом моделировании задач воссоздания параметров режимов электрических сетей в условиях недостаточного количества телеметрических измерений.
Традиционный путь решения этой задачи требует оснащения телеметрическими источниками информации подавляющего большинства трансформаторных подстанций электросети, что требует значительных капиталовложений.
Наиболее целесообразная альтернатива отсутствующим телеизмерениям это воссоздание параметров электропотребления во всех узлах сети, в которых отсутствуют измерения.
Серьезное внимание в работе уделено развитию существующего математического аппарата теории нечетких множеств. Для этого предлагается новый нетрадиционный класс операций и принципов обобщения так называемые α-слабые операции. Этот класс менее требователен к полноте данных о функциях принадлежности нечетких множеств и позволяет оперировать ними даже если известны лишь отдельные множества α-уровня. В работе разработана строгая непротиворечивая аксиоматика основных операций этого класса.
α-слабые операции и принципы обобщения определяют не один образ нечеткого множества, а целое семейство возможных образов, среди которых содержатся и образы, полученные с помощью традиционных операций. Это позволяет формировать нечеткие математические модели в условиях не полностью заданных функций принадлежности нечетких множеств.
Научные результаты, полученные в этой работе, имеют практическую ценность. В первую очередь, это выражается в возможности решать задачу воссоздания состояния для технических систем различной природы (например, систем электроснабжения, газоснабжения, теплоснабжения, водоснабжения и т.д.). При этом, теорема, определяющая условия воссоздания состояния систем, позволяет без предварительных вычислений задать места установки источников информации на графе рассматриваемой системы.
Математическая модель нечеткой квазидетерминизации, основанная на α-слабых операциях и принципах обобщения дает возможность резко снизить требования к экспертной информации и повысить эффективность алгоритмом, т.к. дает возможность перейти от нечетких уравнений к интервальным уравнениям.
Полученные научные результаты внедрены в программном комплексе “Аналитическая Система Воссоздания Электропотребления” (АСВЭП), который с 1997 г. успешно функционирует в оперативно-диспетчерской службе Винницких центральных высоковольтных электросетей (ВЦВЭС). В настоящее время проводится поэтапное внедрение АСВЭП в оперативно-диспетчерской службе Винницких восточных высоковольтных электросетей (ВВВЭС). Планируется также внедрение программного комплекса и в оперативно-диспетчерской службе ОАО “АК Винницаоблэнерго”.
Возможности программного комплекса АСВЭП позволяют: воссоздавать и отображать параметры текущего электропотребления во всех узлах сети, в которых отсутствуют телеизмерения; обеспечивать расчет и контроль текущих интегральных показателей электропотребления (балансов активной мощности и энергии предприятий электросетей и отдельных крупных потребителей); формировать, накапливать и просматривать архивы параметров режима отдельных подстанций и интегральных показателей электропотребления.
Следует особо отметить, что АСВЭП дает возможность непрерывно вести текущие расчеты технических потерь в сети, что позволяет определять места неучтенного электропотребления.
Ключевые слова: детерминизация, квазидетерминизация, линейная неэвклидова метрика, нечеткие множества, слабые операции, принципы обобщения, процессы, математические модели, системы электроснабжения.
Katsyv S. Mathematical models of process determinization in electric supply systems. Manuscript.
Thesis for a Ph.D. science degree by speciality 01.05.02 Mathematical modeling and computing methods. Vinnytsia National Technical University, Vinnytsia, 2004.
The thesis dedicated to a new approach in mathematical modeling of the reproducing of regime parameters of electrical nets in the conditions of deficient telemetric measuring. The conditions of the resolution of some tasks are define. Thanks to use in the complex electrical parameters space the linear unevklid metric, the indetermination of data-in is lower. In fuzzy sets theory the untraditional class of the weak operations and principles of generalization is introduced. The algorithms of the heuristic and fuzzy kvazydeterminization are developed.
Key words: determinization, kvazydeterminization, linear unevklid metric, fuzzy sets, weak operations, principles of generalization, process, mathematical models, electric supply systems.