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

Алгоритмы - шпоры (Final)

.doc
Скачиваний:
81
Добавлен:
31.03.2015
Размер:
2.78 Mб
Скачать

1’ Области рационального использования средств ВТ. Характеристики каждой из областей.

1’ Эксплуатационные расчёты:

а) прогнозирование графика нагрузки.

б) оценка состояния (формирование модели режима на основе телеизмерений,

расчёт режима)

в) расчёт токов КЗ

г) расчет переходных процессов

д) оптимизационные расчеты

2’ Проектные расчёты:

а) расчёты режимов

б) расчет токов КЗ

в) расчет устойчивости

г) выбор трассы линии, выбор конфигурации сети, определение осуществимости.

3’ Планирование развития ЭС:

Решается сложная оптимизационная задача, но с учетом ограничений:

- по устойчивости

- по надежности

- по качеству

4’ Научные исследования.

- от расчёта устойчивости до синтеза структуры автоматического регулирования

5’ Применение компьютеров в контурах управления (противоаварийная автоматика, …), т. е. выполняется задача информационно-вычислительного обеспечения диспетчера.

Несмотря на высокое развитие, факторы экономного использования памяти не потеряли своей актуальности.

2’ Техническая постановка задачи расчету установившегося режима. Самостоятельное и подчиненное значение расчетов УР. Результаты расчета УР

Получение параметров УР имеет самостоятельное значение, т.е. позволяет ответить на вопрос, осуществим ли режим (возможна ли передача мощности от источника к потребителю, меньше ли допустимых значений токи в элементах сети, допустимы ли значения напряжений в узловых точках сети).

Помимо этого, результаты расчёта УР могут являться данными для других расчётов (подчиненное значение):

- проверка статической устойчивости

- проверка динамической устойчивости

- расчет токов КЗ

- решение задач по эквивалентированию.

В результате расчета УР находим:

- напряжения в узлах ЭС

- потокораспределение в сети

- суммарные потери активной и реактивной мощности

- некоторые специаьлные программные комплексы позволяют также:

а) учесть изменение частоты

б) осуществить проверку апериодической статической устойчивости

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

3’ Классификация элементов ЭЭС (активные, пассивные, узловые, линейные). Модели всех элементов ЭЭС. Классификация узлов.

Для выполнения расчетов УР ЭС необходима информация о схеме сети и её параметрах, о потребителях (нагрузках подстанций) и источниках электроэнергии (электростанциях), т. К. реальными элементами энергосистем являются:

1) генераторы электростанций (СГ), синхронные компенсаторы (СК), статические тиристорные компенсаторы (СТК) и другие источники.

2) нагрузки подстанций, представленные в виде мощности нагрузки Sн, управляемого шунтирующего реактора (УШР) или в виде проводимости шунта нагрузки Yш. 3) электрическая сеть – элементы, передающие и преобразующие электроэнергию, т. е. ЛЭП, трансформаторы и автотрансформаторы, а также нерегулируемые статические устройства продольной и поперечной компенсации, схемы замещения которых содержат активные, индуктивные и емкостные сопротивления, образующие продольные или поперечные (связанные с нейтралью) ветви. Все реальные элементы энергосистемы можно разделить на активные и пассивные. Активные – элементы, генерирующие и потребляющие электроэнергию, которые в расчетах УР представляются в виде источников тока (задающих токов). В ряде случаев ветвь намагничивания Т или АТ вводится постоянным отбором мощности, задаваемым аналогично нагрузке, в виде источника тока. Пассивные – элементы, параметры схем замещения которых в расчетах УР принимаются постоянными, не зависящими от напряжений и токов в этих элементах

Все активные элементы – узловые, т.е. поперечные ветви.

Пассивные эл-ты – узловые (поперечные) или продольные ветви.

Для расчёта УР все реальные эл-ты должны быть смоделированы. Расчетная схема эл. системы составляется на 1 фазу с нейтралью.

4’ Модель ветви, модель узла. Генераторные, нагрузочные узлы. Базисный уз.

Для расчёта УР все реальные эл-ты должны быть смоделированы. Расчетная схема эл. системы составляется на 1 фазу с нейтралью. Обычно все элементы, классификация которых приведена выше, ипользуются в современных ПВК для расчета режимов и укладываются в 2 модели элементов ЭС.

1) Модель ветви. Т.к. в состав продольных элементов входят не только линии электропередач, но и трансформаторы или АТ, то модель ветви представляется как модель ЛЭП (П-образная сх. замещ-я) + модель Т или АТ (Г-образная сх. замещ-я). Если схема замещения электрической системы не приводится к одной ступени напряжения, то в схемы зам-я тр-рных эл-тов вводятся идеальные к-ты трансф-ции.

2) Модель узла. Каждый узел может характеризоваться в расчете УР следующей парой параметров:

- активной и реактивной генерируемой мощностью Pг, Qг

- комплексным напряжением (модулем Ui, фазой δ)

- активной и реактивной мощностью нагрузки Pн, Qн

- активной и реактивной проводимостью на землю (Yш=g+-jb)

Генераторный узел м. задаваться 2 путями:

- опорный ген. узел – задаются значения Pг, Uг и ограничения по Qг. неопорный ген. узел – задаются Pг, Qг.Нагрузочный узел в расчетах режимов м. задаваться:

- статическими характеристиками активной и реактивной мощности по напряжению

- при расчете на перспективу – постоянной мощностью

- в виде шунтов на землю с постоянной проводимостью

Базисный узел – узел, в котором фиксируется модуль напряжения Uбаз и фазовый угол δ=0.

5’ Обобщённая статическая характеристика по напряжению (СХН). Возможность представления нагрузки различными способами на основе СХН.

а) представление нагрузки постоянным сопротивлением:

б) представление нагрузки постоянным током:

в) представление нагрузки постоянной мощностью:

6’ Представление исходных данных об ЭЭС для расчета УР в ней.

Таблица по формату f8.0, включает в себя массивы узлов и ветвей.

Массив узлов:

- код массива (0201)

- номер узла

- Uном – номинальное напряжение в узле

- Pг – активная генерируемая мощность

- Qг – реактивная генерируемая мощность

- Pн – активная мощность нагрузки

- Qн – реактивная мощность нагрузки

- Umod – модуль напряжения в узле

- Qmin – ограничение по генерируемой мощности

- Qmax - ограничение по генерируемой мощности

Массив ветвей:

- код массива (0301)

- Nn - номер узла начала ветви

- Nk - номер узла конца ветви

- R – активное сопротивление ветви

- x – реактивное сопротивление ветви

- Yc – проводимость ветви на землю

- Tk – коэффициент трансформации

- Tf – коэффициент трансформации

7 Ядро информационно-вычислительной системы. База данных. Модели данных. СУБД и ее функции.

СУБД – система

управления базами данных.

Основные модели баз данных:

  1. Иерархическая (дерево)

  2. Сетевая (самопроизвольная связь)

  3. Реляционная (таблица). Достоинство: простота модификации и логичная структура. Недостаток: Медленнодействующая.

СУБД – совокупность языковых и программных средств, обеспечивающих доступ, управление, создание баз данных.

Функции СУБД:

  1. Формирование БД (начальное заполнение)

  2. Модификация структуры БД и изменение записей.

  3. Сортировка и упорядочение данных.

  4. Генерация отчетов

  5. Обработка стандартных запросов

  6. Интерфейс с другими программными продуктами

8 Основные уравнения, описывающие установившийся режим ЭЭС. Узловое уравнение состояния ЭЭС (в форме баланса токов, мощности).

(k=m-n)

Из (1) и (2) =>

- Узловое уравнение в форме баланса токов.

- Узловое уравнение в форме баланса мощности.

Y – матрица узловых проводимостей, размерность (n*n)

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

напряжений. (столбец неизвестных (nx1) )

S - столбец заданных мощностей  

Uд^-1 - квадратная матрица (nxn), в которой присутствуют элементы только на главной диагонали.

Система (2) – нелинейна.

Допущения:

  1. Вся система приведена к одному номинальному напряжению

  2. Все мощности в узлах заданы (все генераторные узлы – неопорные)

Примечение: ^ - сопряженное значение

· - комплексное значение

9’/11Базисный и балансирующий узлы. Требования к выбору балансирующего узла.

Балансирующий узел => не задана мощность => после решения

системы (n-1) уравнений можно будет найти Sбал куда будут включены все небалансы.

Базисный узел – столбец.

Балансирующий узел – строка.

Требования к выбору балансирующего узла:

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

  2. В балансирующем узле должны быть резервы активной и реактивной мощности.

  3. Напряжение в этом узле не должно зависеть от режима.

  4. Этот узел должен быть наиболее мощной станцией в ЭС, либо это д.б. часть другой ЭС, высшей ступени иерархии (сэквивалентированная).

Про базисный:Если - решение системы уравнений,то

- так же решение этой системы.

Доказательство:

Следовательно, чтобы избежать неоднозначности системы решений, необходимо задать базисный узел, где δбаз=0.

10 Базисный узел. Необходимость задания базисного узла в расчетах УР.

1. Если - решение системы уравнений, то и

- так же решение этой системы.

Доказательство:

Следовательно, чтобы избежать неоднозначности системы решений, необходимо задать базисный узел, где δбаз=0.

2.

Если мы задали мощности => мы задали потери, а потери, в свою очередь, зависят от напряжения, которое является неизвестным.

=> Во всех узлах задать мощность нельзя.

В базисном узле необходимо задать модуль напряжения и его фазу.

=> Вычеркиваем столбец базисного узла => (n-1) неизвестное напряжение

=>Необходимо вычеркнуть уравнение балансирующего узла, в котором задана мощность.

12Узловое уравнение состояния ЭЭС в форме баланса токов или мощности при условии совмещения базисного и балансирующего узлов.

Узловое уравнение в форме баланса токов.

Узловое уравнение в форме баланса мощности.

Y – квадратная матрица узловых проводимостей, размерность

(n-1)*(n-1)

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

напряжений. (столбец неизвестных (n-1)x1

^

S - столбец комплексно-сопряженных значений мощности в узлах схемы зам. (n-1)*1

^ ^

Uд^-1= diag(1/Ui) - диагональная матрица комплексно сопряженных значений в узлах схемы зам. (nxn), в которой присутствуют элементы только на главной диагонали.

Yjбаз – столбец проводимостей ветвей, связанных с базисным узлом. (n-1)*1

Uбаз – напряжение в базисном узле.

13-14’ Формирование матрицы узловых проводимостей. Алгоритм и программная реализация

1)Если среди элементов, ограничивающих ветвь i-j нет «земли», то она формирует 4 эл-та матрицы Y.

Пример:

Yij = Y31 ; Y13 = Y31 = - Y1-3

Yii = Y33 → +Y1-3

Yjj = Y11 → +Y1-3

nn

nk

Z=R+jX

Yc

3

1

Z1-3

Yc1-3

2

0

Z2-0

1

2

Z1-2

Yc1-2

0

1

Z1-0

4

3

Z3-4

2

3

Z2-3

2)Если среди номеров узлов, ограничивающих ветвь, есть ветвь земли, то проводимость этой ветви входит только в 1 диаг. элемент

Yii = Y22 → +Y2-0

Y=

Y1-2+Y1-3 +Y1-0 +YC1-2/2 +

YC1-23/2

- Y1-2

- Y1-3

0

- Y1-2

Y1-2+Y2-3 +Y2-0+YC1-2/2

- Y2-3

0

- Y1-3

- Y2-3

Y1-3+Y2-3 +Y3-4+YC1-3/2

- Y3-4

0

0

- Y3-4

Y3-4

Программная реализация (↓ - переход на новую строку)

Subroutine matry (nn,nk,z,yc,m,n,y,l)

Dimension nn(1),nk(m),yc(n)

Complex z(m), y(l,l)

Do i=1,n ↓ do j=1,n ↓ y(i,j)=(0,0) ↓ end do ↓ end do ↓ do i=1,n ↓ y(i,j)=(0,0) ↓ do k=1,m ↓ i=nn(k) ↓ j=nk(k) ↓ c=(1,0)/z(k) ↓ if(i*j.ne.0)then ↓ y(i,j)=-c ↓ y(j,i)=-c ↓ y(i,i)=y(i,i)+c+0.5*cmplx(0,yc(k)) ↓ else ↓ j=i+j ↓ end if ↓ y(j,j)=y(j,j)+c+0.5*cmplx(0,yc(k)) ↓ end do ↓ end do

15-16’ Учет трансформаторных связей при формировании матрицы узловых проводимостей

17’ Прямые и итерационные методы решения систем линейных алгебраических уравнений (СЛАУ). Необходимость использования этих методов для решения задач расчета установившихся режимов.

1) прямые (точные) методы – дают решение за конечное число арифметических итераций, которое определяется вычислительной схемой метода а также порядком и структурой матрицы коэффициентов СЛАУ. (Гаусса)

2) Итерационные – это методы, в которых за счет многократного выполнения однообразных вычислений (итераций) получается решение с наперед заданной точностью. Чем выше точность, тем больше количество итераций. (м. Зейделя, min ф-ции ошибки, релаксационные методы)

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

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

В результате получаем СЛАУ, которую необходимо решать на каждой итерации.

18’ Метод Гаусса с обратным ходом. Вычислительная схема пр и обр хода. +/- метода

1 этап: Прямой ход: n шагов – матрица коэффициентов приводится к треугольному виду:

2 этап: Обратный ход – находятся Xn-1,….X1

1 этап: , а для для К-го шага:

2 этап:

«+» : 1) прямой (точный) метод, получаем решение за конечное число шагов. 2) если диагональн элементы >> недиагональных, то хорошая точность расчета

«-» : 1)вычисления по формулам (2) и (3) дают большую погрешность, если диагональн элементы → 0 (малы)

2)вычисления по формулам (4) и (5) могут дать большие погрешности, если вычитаются близкие величины (при detA → 0). Неточность задания исходных данных может привести к большой погрешности, ПРИМЕР: сис-ма уравнений 1000x1+x2=1002; 2001x1+2x2=2005; решая ее получим x1= 1, x2= 2. Изменив исх данные с ε=0,1% система будет: 1001x1+x2=1002; 2001x1+2x2=2005; решив новую систему получим: x1= -1, x2= 2003.

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

19’ Алгоритмическая и программная реализация метода Гаусса с обратным ходом.

Алгоритм:

Прямой ход:

Обратный ход:

Программная реализация:

Subroutine Gauss (A,B,n,l)

Complex A(l,l), B(n)

do k = 1,n-1

b(k) = b(k) / a(k,k)

do j = k+1, n

a(k,j) = a(k,j) / a(k,k)

b(j) = b(j) - a(j,k) ∙ b(k)

do i = k+1, n

a(i,j) = a(i,j) – a(i,k) ∙ a(k,j)

end do

end do end do

! обратный ход

b(n) = b(n) / a(n,n)

do i = n-1, 1, -1

do j = i+1, n

b(i) = b(i) – a(i,j) ∙b

end do

end do

end

20’ Принципы учета слабой заполненности сетевых матриц при использовании метода Гаусса.

1) Порядок исключения ненулевых элементов в методе Гаусса должен быть таким, чтобы в процессе исключения появились минимальное кол-во новых ненулевых элементов;

Итог: 1) Надо либо менять порядок следования исходных узлов

2) Либо не изменяя порядка исключения переименовывая узлы (см билет 23)

Порядок исключения: в первую очередь должны исключаться узлы менее всего связанные с другими узлами (имеющие наименьшее число инцидентных связей с другими элементами)

2) Хранить в памяти машины нужно только ненулевые элементы и только с ними нужно производить арифметические операции (см билет 24)

21’ Порядок исключения неизвестных в методе Гаусса с обратным ходом.

1) Порядок исключения ненулевых элементов в методе Гаусса должен быть таким, чтобы в процессе исключения появились минимальное кол-во новых ненулевых элементов;

Итог: 1) Надо либо менять порядок следования исходных узлов

2) Либо не изменяя порядка исключения переименовывая узлы (см билет 23)

Порядок исключения: в первую очередь должны исключаться узлы менее всего связанные с другими узлами (имеющие наименьшее число инцидентных связей с другими элементами)

2) Хранить в памяти машины нужно только ненулевые элементы и только с ними нужно производить арифметические операции (см билет 24)

22’ Коэффициент заполненности матриц. Хранение ненулевых элементов матриц.

- кол-во элементов; n – кол-во узлов; m – ветвей

В больших схемах

Вывод: чем больше кол-во узлов схемы тем меньше коэффициент заполненности

Хранить необходимо только ненулевые элементы матрицы Y, причем только верхний Δ, так как матрица Y – симметрична

(см билет 24)

Требования к схемам хранения матриц.

1) Компактность

2) Простота формирования

3) Простота использования

а) простота выборки элементов

б) гибкость изменения хранимой информации

(т.к. поле 1 шага появл. новые ненул элементы)

I. Задаем только ненул. элементы:

VALUE

Y12

Y43

Y33

Y23

- значения

ILINE

1

4

3

2

- номер строки

ICOL

2

3

3

3

- номер столбца

Экономия от n2 до 4n чисел

23’ Алгоритмы упорядочения, их классификация.

2) компромисс между программно не очень сложной перенумерацией и экономией памяти машины

3) процесс упорядочивания проходит до начала исключения узлов

4) осуществляется на каждом шаге прямого хода

5) нумерация, предполагающая миним. требуемый объем памяти, мин. суммарное кол-во ненулевых элементов

При программной перенумерации происходит только имитация процесса исключения неизвестных

24’ Хранение слабозаполненных матриц. Схемы упаковки матриц. Требования к схемам хранения матриц.

1) Компактность

2) Простота формирования

3) Простота использования

а) простота выборки элементов

б) гибкость изменения хранимой информации

(т.к. поле 1 шага появл. новые ненул элементы)

I. Задаем только ненул. элементы:

VALUE

Y12

Y43

Y33

Y23

- значения

ILINE

1

4

3

2

- номер строки

ICOL

2

3

3

3

- номер столбца

Экономия от n2 до 4n чисел

II. Организуем отдельно массив диагональных элементов, для которых не храним номер строки и столбца:

DIAG

Y11

Y22

Y33

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

VALUE

Y12

Y43

Y33

ILINE

1

4

2

ICOL

2

3

3

III.

VALE

Y12

Y15

Y21

NZERO – количество ненул. элементов, связанных с данным узлом, искл. ILINE

ICOL

2

5

1

DIAG

Y11

Y22

Y33

NZERO

2

2

3

IV.

DIAG

Y11

Y22

Y33

Чтобы не перестраивать весь массив, появл. элементы дописываются в конец, но перекомпилируется массив INEXT

ILINE

1

1

2

ICOL

2

5

1

INEXT

2

2

3

VALUE

Y12

Y15

Y21

25’ Алгоритм формирования матрицы У в компактной форме.

У11

У12

У13

У15

×

У22

У23

×

×

У33

У34

У35

×

У44

×

×

У55

diag

У11

У22

У33

У44

У55

nzero

3

1

2

0

0

Только те эл-ты, которые есть в верхнем Δ

value

У12

У13

У15

У23

У34

У35

icol

2

3

5

3

4

5

nadr

1

4

5

7

7

Создается доп массив NADR – номер ячейки массива VALUE, с которого начинаются проводимости, связанные с этим узлом

7 – потому что 6 эл-тов всего в VALUE, показывает что больше ни с чем не связано.

26’ Алгоритмическая и программная реализация формирования матрицы Y в компактной форме.См.25 Программа:

Subroutine ysz (nn, nk, z, yc, y, diag, nzero, value, icol, nadr, n, m)

Complex z(m), yc(m), diag(n), value (1)

Dimension nzero(n), icol(1), nadr(n)

do i=1,m

diag(i)=(0.0)

nzero(i)=0

end do; do i=1,m

i1=nn(i)

i2=nk(i)

y(i)=(1,0.)/z(i)

if(i1*i2.ne.0) then j=i2

if (i2.gt.i1) j=i1

nzero(j)= nzero(j)+1

diag(i1)= diag(i1)+y(i)+0.5*yc(i)

else i2=i1+i2

endif

diag(i2)= diag(i2)+y(i)+0.5*yc(i)

end do

nadr(1)=1 do i=2,n

nadr(i)= nadr(i-1)+ nzero(i-1)

end do do i=1,m

i1=nn(i)

i2=nk(i)

if(i1*i2.ne.0) then j=i2

if (i2.gt.i1) j=i1

k=nadr(j)

nadr(j)= nadr(j)+1

value (k)=-y(i)

icol(k)= i1+i2=j

endif

end do end return

27’ Итерационные методы расчета УР. Алгоритм расчета УР методом Зейделя.

1)Поузловые методы – м., где исх переменные ищутся из каждого след. ур-я последовательно.

2)Одноврем. реш. ур. – все ур. системы реш-ся одновременно на каждой итерации.

Бал.токов: метод Гаусса, обращения, факторизации матрицы

Бал. Мощностей: метод Ньютона

Эффективность:

1)кол-во вычислений на кажд. Итерации;

2)суммарное кол-во итераций;

3)суммарное время вычислений.

Нет абсолютно эффективных и абсолютно сходящихся методов

28’ Метод Зейделя применительно к решению нелинейного узлового ур-я в форме баланса токов.

YU=Uд-1S-YjбазUбаз примечание: над Uд-1S надо ^

Y11U1+

…..

+Y1nUn=

S1/U1-

Y1базUбаз

….

…..

Yn1U1+

…..

+YnnUn=

Sn/Un-

YnбазUбаз

Правая линейная часть, левая нелинейная часть.

U1(k)=1/Y11*( S1/U1(k-1)- Y1базUбаз-Y12U1 (k-1)- Y13U3 (k-1)-….. -Y1nUn (k-1))

U2(k)=1/Y22*( S2/U2(k-1)- Y2базUбаз-Y21U2 (k)- Y23U3 (k-1)-….. –Y2nUn (k-1))

………..

Un(k)=1/Ynn*( Sn/Un(k-1)- YnбазUбаз-Yn1U1 (k)- Yn2U2 (k)-….. –Yn-1nUn-1 (k))

Итерационная формула Зейделя:

Ui(k)=1/Yii*( Si/Ui(k-1)- YiбазUбаз(i-1) YijUj (k)-Σ YijUj (k-1))

Sбал=UбалΣYбалjUj

Должно быть задано макс. кол-во итерацмй,номер балансирующего узла,точность

max ׀Ui (k)- Ui (k-1)׀≤εk≈0.5÷0.05

max ׀Pнб׀≤εр ≈ 0.5 МВт max ׀Qнб׀≤εQ ≈ 1 Мвар

Примечание: везде над S/U надо поставить ^