Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Современные методы анализа электрических систем

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
19.49 Mб
Скачать

Впрактических случаях

т= 3; ttiik=4;

т= 2; тш = 3,3.

Если <т»/сг<0,2, то ограничения, связанные с пропускной способностью линии, практически несущественны.

Таким образом, по экономической оценке определяет­ ся tnik, после чего из выражения (15-43) получаем про­ пускную способность

 

= tnik^ik — E/ft,

(15-51)

где Fib находится по (15-42)

и (15-31):

 

 

= S Си,

Сы, а (Кг - Ti -

X t), (15-52)

а

— из выражения (15-42):

 

 

 

< 4 = 2

(Сы, /^)*

 

15-3,д. Соединительная линия (связь) между дву­ мя системами. Предполагая, что обе системы имеют оди­ наковый состав электростанций, для определения сред­ неквадратичного отклонения резерва и дефицита мощно­ сти можно исходить из следующего:

Ft =

aF; о1=

<з]/‘а ;

(15-53)

Е2 =

(1—a)F;

aa = e j/l —а,

 

где

а = - * 1___

аК г + К г

Ограничение, определяемое уравнением (15-37), можно записать следующим образом:

/ а 6,+ / 1—<**»=£

(15-54)

Это уравнение есть уравнение прямой на плоскости, различные положения которой показаны на рис. 15-3 для коэффициента надежности f =т = 3.

Ограничение, обусловленное пропускной способно­ стью связи, в этом случае рассмотреть просто, так как

коэффициент распределения тока С равен 1 или 0. Оче­ видно, что

:

Если базисным узлом является узел 2, то дефицит мощности системы 1 покрывает система 2. Ограничение представляется в виде прямой, параллельной оси коор-

Рис.

15-3. Изменение мощности

при

объединении двух

систем.

Плотность

вероятности

дефицита

мощности-

г2=*53,99% г3= 4,43%; г4=

 

 

=0,13%; a=Ki/K

 

динат 2. Аналогично в случае базисного узла 1 ограни­ чивающая прямая параллельна оси координат 1. Рас­ стояния обеих ограничивающих линий от начала коор­ динат равны:

_ F n

+ Kn _Fx + Kn .

'“21 ---

--

F12

Ч~ ^12

^2 Н~ Кп

 

т12

 

 

 

°1 2

Согласно § 15-3, г примем:

т^ — т^ — А.

Тогда

^ . = 40,—ф .

Поэтому

"Т1=в(4 — 3]/сГ) у а,;

^ • = ( 4 — 3 У Т ^ а ) У Т = а .

Из величин Кч\ и /Си следует выбрать большую. Для а, меньших 0,5 (за исключением очень малых значений а), большее значение имеет величина Кп- При дефиците

332

мощности в соответствии с графиком на рис. 15-3 для малых значений а передача мощности из малой системы

в■большую сильно ограничена.

Вобласти а=0,5, ..., 1, большие значения дает /С12. Результаты расчетов приведены на рис. 15-4, по которо­ му можно определить экономию резервной мощности и

Рис. 15-4. Зависимость экономически целесо­

образной пропускной способности (/Сл) меж­ системной связи от вели­ чины а (Д/С — экономия резерва мощности).

требуемую пропускную способность. Экономия резерв­ ной мощности выражается следующей зависимостью:

^ = ( ^ + К ' П ^ а - l)f.

(15-55)

Подставляя значение J — 3, получаем:

^ З О / а - К П ^ - ! ) .

15-3,д. Методы анализа для простых систем. Иссле­ дования, изложенные в § 15-3,г, можно продолжить для объединения трех систем. Можно определить, что суммарная пропускная способность линий всегда больше, чем мощность присоединения. Среднеквадратичное от­ клонение дефицита мощности отдельных присоединений равно 0,8—1,3 от общей величины среднеквадратичного отклонения. Для мощностей отдельных линий, состав­ ляющих контур, оно равно 0,5—0,9 от среднеквадратич­ ного отклонения общего дефицита.

Необходимо, чтобы отношение мощностей отдельных

систем энергообъединения не

Превосходило 1 :3

или

1:4, так как иначе не будет

удовлетворительным

ис­

пользование межеистемной связи.

 

S33

Плановая передача мощности и восполнение дефици­ та мощности обычно осуществляются совместно. Резуль­ тирующий поток мощности по передаче образуется пу­ тем обычного сложения Пусть в г-й системе /С, есть располагаемая мощность, Г, — пиковая нагрузка, тогда можно определить величину общего резерва мощности

К — Т.

Плановая передача мощности из г-й системы:

PH = J - K i - T i .

(15-56)

Тем самым с помощью методов, изложенных в гл. 6, определятся плановые потоки мощности. Однако эту ве­ личину следует рассматривать как идеальную, ее необ­ ходимо уточнить путем учета мощности, требуемой для покрытия дефицитов, возникающих в отдельных систе­ мах. В этом случае математическое ожидание дефицита мощности уже не надо определять; оно содержится в плановой составляющей. По величине среднеквадра­ тичных отклонений дефицитов мощностей можно опре­ делить приращения мощностей по передачам, которые могут достигать 0,5—0,9 от общей величины среднеквад­ ратичного отклонения. Такие приращения необходимо учитывать.

Например, в энергообъединении общей мощностью 20 Гвт среднеквадратичное отклонение общего дефицита равно 2%, что по отношению к плановой передаче мощ­ ности может дать приращение в 200—360 Мет.

Такие приближенные оценки могут быть использова­ ны при проектировании объединений электрических си­ стем.

15-4. ИССЛЕДОВАНИЕ ЭЛЕКТРИЧЕСКИХ СЕТЕЙ МЕТОДОМ МОНТЕ-КАРЛО

Исследуя «надежность» сети, необходимо изучать ее поведение при случайно происходящих процессах. Для этой цели удобно пользоваться методом статистических испытаний, известным под названием метода Монте-Кар­ ло (Л. 32, 154—156]. Стохастическую модель сети можно составить из следующих частей:

1) нз нагрузок, представленных нормально распре­ деленными случайными величинами с математическими ожиданиями, равными максимальным нагрузкам узлов;

334

2) из функции распределения аварийного снижения мощности генераторов каждого узла сети, эти функции можно задать с учетом стохастических закономерностей повреждений отдельных генераторов;

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

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

а) дефицит мощности в системе; б) ограничение пропускной способности соединитель­

ной линии:

6-1) в силу достижения максимально передаваемой мощности при отсутствии повреждений в сети;

6-2) при повреждении 1, 2, 3, 4-х элементов сети; в) снижение напряжений на нагрузках, которое мо­

жет возникнуть:

в-1) вследствие дефицита реактивной мощности; в-2) вследствие неблагоприятных условий передачи

электроэнергии; г) нарушение динамической устойчивости.

Наряду с испытаниями, моделирующими поврежде­ ния в системе, требуется также определение дефицита мощности или энергии. При исследовании важно пра­ вильно выбрать последовательность расчетов. Если, на­ пример, имеет место случай «а», то становятся излиш­ ними исследования случаев «б»—«г». В случае «6-2* целесообразно отметить элемент сети, приводящий к ограничению передачи мощности. Случаи «в» и «г» требуют дальнейших исследований, так как даже при использовании современных электронных вычислитель­ ных машин они представляют еще трудную задачу.

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

Точность метода Монте-Карло характеризуется раз­

ностью между математическим ожиданием | случайной величины | (с примерно нормальным распределением)

335

исследуемой системы и частостью а, определяемой путем статистических испытаний [Л. 32]. Эта разность удовлет­ воряет оценке

 

(15-57)

где а — среднеквадратичное отклонение

т — коэффи­

циент надежности [см. зависимость (15-4)], по которому

можно определить, с какой вероятностью

может

быть

удовлетворено

приведенное

неравенство

(например,

<2(1,95) =0,915,

или Q(3) =0,9986); N — число необходи­

мых испытаний.

 

 

 

 

Если вероятность повреждений q мала и при каждом

испытании одинакова, то из (15-6) a ^ Y q ,

а число тре­

буемых испытаний:

 

 

 

 

 

 

(15-58)

Пусть, например, <7=0,04;

е/<7=0,1, т = 2, тогда

N —

= 1 000.

Вероятности повреждений основных элементов сети (данные Electricite de France, EdF):

ЛЭП высокого напряжения: <7=0,001; трансформатор высокого напряжения: <7=0,003 (пе

риод времени 1/10 года). Вероятность повреждения машин:

вгод ввода в эксплуатацию <7 = 0 ,1 ;

вследующем году <7 = 0,08;

востальные годы <7=0,05.

В других энергетических системах вероятность по­ вреждения машин также колеблется в пределах 0,05—0,1.

Вероятность повреждения прототипной машины в 1—2 года ввода в эксплуатацию <7=0,3.

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

Посредством сечения сети (см. гл. 2) система разби­ вается на две части и исследуется пропускная способ­ ность линий между ними. По распределению, заданному

336

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

Необходимая передаваемая мощность между двумя частями сети, разделенными сечением, определяется ме­ тодом Монте-Карло. Путем проведения 5000 статистиче­ ских испытаний, обеспечивающих достаточную точность, на ЭЦВМ типа ШМ-7094 для системы, состоящей из 20 узлов, за 5 мин определяется распределение вероятно­ стей для требуемой пропускной способности, которое сравнивается с данными, полученными в результате расчетов потокораспределения мощностей.

Выбирая примерно 200 различных сечений, можно достаточно подробно исследовать систему. Время ма­ шинного расчета всего исследования составляет при этом 5x200= 1000 м ин^ 16 ч.

П Р И Л О Ж Е Н И Е

НЕКОТОРЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ МАТРИЦ

ПЛ. ОСНОВНЫЕ СООТНОШЕНИЯ О б о з н а ч е н и е м а т р и ц ы : А; в случае т строк и п столб-

ЦОВ Агп п* В случае квадратной матрицы m = n

обозначение А п. Для

квадратной матрицы число п называют порядком матрицы.

Элементы матрицы:

 

 

. aln

 

 

 

#п

#12 •

 

 

А =

#21

#22 •

#2П

 

( П Ы )

 

_#ml

#m2 •

• #mn _

 

 

(Элементы можно обозначать также прописной

буквой Ац.)

Т р а н с п о н и р о в а н н а я

с о п р я ж е н н а я

м а т р и ц а по­

лучается заменой

строк

столбцами, а отдельных

элементов — на

комплексно-сопряженные с ними. Обозначение: А*. Если А*=[Ьц], то bji через элементы матрицы А выражаются следующим образом:

dij—bjib

Матрица с действительными

элементами с и м м е т р и ч н а ,

если

 

(П1-2)

А=А*.

Матрицу с комплексными элементами называют эрмитовой, если

А =А *.

(П1-3)

Для с т о л б ц о в о й м а т р и ц ы

( в е к т о р а - с т о л б ц а )

приме­

няется обозначение

 

 

 

# 1

 

для с т р о ч н о й м а т р и ц ы ( в е к т о р а - с т р о к и ) — обозначение

а* = [л,.........а„].

338

С к а л я р н о е п р о и з в е д е н и е д в у х в е к т о р о в

 

 

п

 

 

 

 

 

a*b =

2 a%bi.

 

 

(П1-4)

 

 

i=1

 

 

 

 

Д и а д н о е ( в н е ш н е е )

п р о и з в е д е н и е д в у х в е к т о р о в

(коротко: диада)

axbx

ахЬг .

 

 

 

 

 

.

alf)n

 

ab* =

й2Ьх

и2Ь2 •

.

агЬп

(Ш-5)

 

 

 

 

 

_

arfil

дпЬ2 .

апЬп _

 

С у м м а д в у х м а т р и ц :

 

 

 

 

 

А + В = С

определяется для матриц, имеющих одинаковое число строк и столб­ цов,

Сгj = di з Ь гj.

(П1-6)

П р о и з в е д е н и е д в у х м а т р и ц :

АВ = С

имеет смысл, если число столбцов матрицы А равно числу строк ма­ трицы В, и определяется согласно правилу

п

CU = S aihPh1'

(П1-7)

k-\

 

В общем случае произведение не обладает свойством перестано­

вочности (коммутативности).

матрица, у кото­

Д и а г о н а л ь н о й м а т р и ц е й называется

рой не равны нулю только элементы, находящиеся на главной диаго­ нали Диагональная матрица обозначается так:

<С^1» &2>

ИЛИ

 

 

diag (ах, а2, ... , ап)

 

(см. разд, 6-1,е).

 

 

Е д и н и ч н а я м а т р и ц а :

 

 

1= <1, 1, ...,

1>.

(П1-8)

Обозначение единичной матрицы п-го порядка Iw.

единичной

Е д и н и ч н ы е в е к т о р ы — это

отдельные столбцы

матрицы. Обозначение первого столбца е, еп; обозначение последнего столбца f, fn.

“ 1 ~

-

0

 

0

 

.

 

е =

; f==

;

(П1-9)

 

 

0

 

_ 0 _

_

1

_

339

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

Н-ь •••, iv

Минорная матрица, у которой опущена t-я строка и у-й столбец, обо­ значается:

At j.

Гиперматрицами называются матрицы, элементами которых в свою очередь являются матрицы:

 

ГАП

А121

 

 

 

LА21

A22J

 

 

Здесь Аи не следует путать с минорными матрицами.

 

Обозначение о п р е д е л и т е л я

к в а д р а т н о й

м а т р и ц ы :

|А] =

det А =

I й"

'

й'п

(П1-10)

 

 

| йп1

. . .

апп

 

Определитель квадратной матрицы:

 

 

 

 

а

 

 

 

 

IА | =

£ sig n

од, а

. д ,

 

1=1

где суммирование надо вести по всевозможным перестановкахм чисел

ви о% ...» Оп из

ряда 1,

2,

, . n; sign<r=l, если перестановка о ь

О2 , . . Оп четная, и — 1, если перестановка нечетная.

Аналитическая

ф о р м а

т е о р е м ы Л а п л а с а :

 

k

х+а*

 

2

|а,|= Е (-1)“ '

< П М | >

 

 

 

Эта формула

указывает,

что из верхних

k строк матрицы Ап следует

взять ^ ^ ^

различных

миноров

£-го порядка и умножить

на соответствующие дополнительные миноры (пk)-ro порядка. Все полученные произведения надо сложить с соответствующими зна­ ками.

Теорема

Коши — Вине (см. § 4-1).

О п р е д е л и т е л ь г и п е р м а т р и ц ы :

[А* В,

I А, I • |Ва — А аА Г^! I, если | А! I Ф

 

} (П1-12)

IА2 в 2

 

 

Ы*о. J

. | В2 1• | А! — BiB2

Ag I* если

| В2 1

 

Алгебраическое дополнение матрицы А к элементу ац\

 

At i = ( - 1)<+МА « | -

(П1-13)

П р и с о е д и н е н н а я м а т р и ц а

к матрице

А определяется так:

 

adjA = [ А Н];

 

 

A (?dj А) = (adj А) А = / А 11.

(П1-14)

340

Соседние файлы в папке книги