3.2 Расчет параметров режимов

Для расчета нормальных режимов используется разработанная на кафедре программа W_MAESTRO [3]. Программа составлена на алгоритмическом языке DELPHI и предназначена для расчетов нормальных установившихся режимов электрических сетей с числом узлов до 50 и числом ветвей до 70 как с приведением, так и без приведения к одной ступени напряжения. Нагрузки в узлах могут быть заданы либо постоянными мощностями, либо постоянными проводимостями, либо сочетанием постоянных мощности и проводимости. Часть узлов может быть задана с фиксированным модулем напряжения и постоянной активной мощностью. В одном из узлов (балансирующем) должен быть задан вектор напряжения по модулю фазе.

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

3.2.1.Метод и алгоритм расчета установившегося режима электрической сети

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

, (i=1, 2, ... , N), (3.1)

где N_I- число ветвей, примыкающих к узлу i, ij перетоки мощности по ветвям; N= n-1, где n - число узлов сети. При заданных мощностях_HI, _ГI во всех узлах и вектора напряжения _БУ балансирующего узла система (3.1) содержит N уравнений относительно N неизвестных напряжений в узлах сети.

Для реализации вычислений в вещественных числах от уравнений (3.1) целесообразно перейти к эквивалентной системе 2N уравнений баланса активных и реактивных мощностей в узлах сети

(3.2)

Традиционной схемой решения системы нелинейных алгебраических уравнений n-го порядка

f_I (x1,x2,...,xn)=0 , i=1,2,...,n (3.3 )

методом Ньютона предусматривается выполнение следующих шагов.

1. Задание вектора начальных приближений [] к решению системы уравнений (3.3) и точности Eps, с которой нужно получить решение.

2. Определение "невязок" уравнений при подстановке начального [] или последующих [] приближений в уравнения (3.3)

(3.4)

3. Формирование линеаризованной системы уравнений в окрестности очередного приближения []. Раскладывая каждое из уравнений (3.3) в окрестности векторав ряд Тейлора и ограничиваясь членами нулевого и первого порядков разложения

,

получим систему линейных уравнений

(3.5)

где [df/dx] - матрица Якоби, - вектор невязок итерации k;-вектор поправок к решению на К-й итерации.

4. Решение линеаризованной системы уравнений (3.5) относительно вектора и определение нового приближения

=+. (3.6)

5. Проверка условий окончания итерационного процесса и контроль сходимости. Признаком окончания итерационного процесса может служить малость всех поправок или невязокна шаге расчета. Если всеилименьше заданной величины Eps, расчет заканчивается. В противном случае расчет повторяется, начиная с п.2.

В узловых уравнениях (3.2) "невязками" являются небалансы активной и реактивной мощности ,в узлах сети, а неизвестными величинами, которые должны быть определены в результате решения этих уравнений - модули и углы векторов напряжений U. Линеаризованная система уравнений вида (3.5), соответствующая уравнениям (3.2), содержит частные производные от небалансов активнойи реактивноймощности в узлах сети по модулям и углам узловых напряжений

; ;;(i, j = 1,2,...,n) ,

поправки ,и может быть записана в виде :

(3.7)

В системе (3.7) уравнения баланса активных и реактивных мощностей записаны попарно для всех n узлов сети.

Диагональные элементы матрицы Якоби получаются как частные производные от небалансов ,в узле i по модулю и углу вектора напряжения в этом узле U_I:

(3.8)

а недиагональные элементы матрицы Якоби - это частные производные от небалансов ,в узле i по модулю и углу вектора напряжения U_I в узле j, примыкающем к узлу i:

(3.9)

При записи системы линеаризованных уравнений на шаге расчета в форме (3.7) все узлы (регулируемые {N_PU} и нерегулируемые {N_PQ}) и связи между ними представлены блоками второго порядка 2х2, алгоритм решения узловых уравнений с матрицей, состоящей из унифицированных блоков, упрощается за счет того, что отпадает необходимость учитывать различия между регулируемыми и нерегулируемыми узлами.