CC BY
81
МОДЕЛЮВАННЯ І УПРАВЛІННЯ ВЕЛИКИМИ СИСТЕМАМИ
УДК 621.311
В.П. ВОЛОБОЕВ, В.П. КЛИМЕНКО
ОБ ОДНОМ ПОДХОДЕ К МОДЕЛИРОВАНИЮ ЭНЕРГОСИСТЕМ
Abstract: Correct statement of a task of modelling of distribution of streams of powers in the established operating mode of a power supply system is offered. It is shown that correctness of setting depends on the selection of the variables of system of equations (SE), which describe power system. It is proposed to use the voltage of the components of the branches of the tree of the graph of a power system as variable SE, whose correct selection is carried out with the composition of the topological matrix of the contours of graph and depends on the parameters of the components of power system, and the modification of Newton's method of the solution of SE, in which it is realized checking of correctness of task when resolving and its correct formulation.
Key words: power system, the simulation of the flow distribution of power before the steady-state regime, nodal potential method, the divergence of the iterative process of the method of the solution, bad conditionality, the correct formulation of the task of simulation, method of voltage a component of branches of a tree the graph of power system.
Анотація: Запропоновано коректну постановку задачі моделювання розподілу потоків потужностей у сталому режимі роботи енергосистеми. Показано, що коректність постановки залежить від вибору змінних системи рівнянь (СР), що описують енергосистему. Запропоновано застосовувати напруги компонент гілок дерева графа енергосистеми як змінні СР, коректний вибір яких виконується при складанні топологічної матриці контурів графа й залежить від параметрів компонентів енергосистеми, і модифікацію методу Ньютона рішення СР, у якому реалізована перевірка коректності задачі в процесі рішення і її коректне формулювання.
Ключові слова: енергосистема, моделювання розподілу потоків потужностей у сталому режимі, метод вузлових потенціалів, розбіжність ітераційного процесу, погана обумовленість, коректна постановка задачі, метод напруг компонентів гілок дерева графа енергосистеми, модифікація методу Ньютона.
Аннотация: Предложена корректная постановка задачи моделирования распределения потоков мощностей в установившемся режиме работы энергосистемы. Показано, что корректность постановки зависит от выбора переменных системы уравнений (СУ), описывающих энергосистему. Предложено применять напряжения компонент ветвей дерева графа энергосистемы в качестве переменных СУ, корректный выбор которых выполняется при составлении топологической матрицы контуров графа и зависит от параметров компонент энергосистемы, и модификацию метода Ньютона решения СУ, в котором реализована проверка корректности задачи в процессе решения и ее корректная формулировка. Ключевые слова: энергосистема, моделирование распределения потоков мощностей в установившемся режиме, метод узловых потенциалов, расходимость итерационного процесса, плохая обусловленность, корректная постановка задачи, метод напряжений компонент ветвей дерева графа энергосистемы, модификация метода Ньютона.
1. Введение
В основе анализа и планирования работы электроэнергетических систем лежит расчет распределения потоков мощностей в установившемся режиме, базирующийся на математическом моделировании электромагнитных процессов в системе. Для обеспечения оперативного управления энергосистемой ее моделирование должно выполняться в реальном масштабе времени с учетом изменения конфигурации реальной энергосистемы, условий генерации и потребления электроэнергии.
В настоящее время в России и странах СНГ известно несколько основных программных комплексов для расчёта установившихся режимов энергосистем: АНАРЭС-2000, разработанный в Институте диспетчерского управления энергетических систем (г. Новосибирск) совместно с Институтом систем энергетики им. Мелентьева СО РАН (г. Иркутск) [1, 2], комплекс программ EnergyCS, предназначенный для выполнения электротехнических расчетов при проектировании и эксплуатации электроэнергетических систем любой сложности, разработчик CSoft Development [3],
© Волобоев В.П., Клименко В.П., 2009
ISSN 1028-9763. Математичні машини і системи, 2009, № 4
комплекс «Космос» [4], разработки Института электродинамики НАН Украины (г. Киев) и целый ряд других программных комплексов [5-7].
Результаты выполненных электротехнических расчетов должны иметь достаточную степень достоверности, так как ошибки в результатах моделирования могут стать причиной аварий и даже катастроф. Достоверность результатов моделирования вытекает из корректности принятой математической модели энергосистемы и зависит от того, какие выбраны модели объектов энергосистемы и, соответственно, параметры этих моделей, а также от качества вычислительного процесса математического моделирования системы.
В данной работе рассмотрена зависимость достоверности результатов только от качества вычислительного процесса моделирования режимов энергосистемы.
2. Постановка задачи
Основные промышленные программные комплексы имеют как индивидуальные особенности, так и ряд общих черт. К общим признакам можно отнести применяемые методы составления математического описания энергосистемы и решения составленных уравнений. Как следует из литературных данных [1, 3, 4], для составления математического описания энергосистемы применяется только метод узловых потенциалов в разной модификации. В результате чего задача электротехнических расчетов сводится в конечном итоге к решению полученной системы нелинейных уравнений (СНУ) в комплексной области, для решения которой применяется метод Ньютона - Рафсона [1, 3, 4, 8]. Практически во всех комплексах программ в качестве критерия отсутствия решения моделируемой при заданных параметрах объектов энергосистемы [1, 3, 4] используется несходимость итерационного процесса по методу Ньютона.
Как следует из литературы [9], отличительной особенностью итерационного метода Ньютона - Рафсона, применяемого для решения СНУ, есть быстрая сходимость в области решения. Быстрая сходимость в окрестности решения достигается за счет применения результатов решения системы линейных алгебраических уравнений (СЛАУ), составленной определенным образом, на каждом итерационном шаге. Известно, что численное решение многих классов СЛАУ встречает значительные трудности, т.е. относится к классу некорректных задач. К числу таких задач относятся системы большой размерности с заполненными матрицами произвольной структуры, а также плохообусловленные СЛАУ. Вышеприведенные программные комплексы как раз рассчитаны на моделирование энергосистем большой размерности. Следует заметить, что для энергосистем большой размерности не всегда можно задать начальное приближение итерационного процесса в окрестности решения, особенно в случае нестандартной ситуации. Это означает, что вышеперечисленные факторы, а не отсутствие решения моделируемой энергосистемы, могут быть источником несходимости или даже расходимости итерационного процесса по методу Ньютона.
В настоящее время разработка соответствующих методов решения некорректных задач рассматривается как чисто математическая задача. Методы решения некорректных задач были предложены в работах А.Н. Тихонова [10], М.М. Лаврентьева [11], В.К. Иванова [12]. Анализ
литературы показал, что исследования методов решения некорректных задач в области моделирования энергосистем практически не проводились.
Следует отметить, что в литературе практически нет примеров системного подхода к корректной постановке технической (физической) задачи, ориентированной на компьютерные технологии, как это было в эпоху ручного счета. Как исключение из общего правила, можно отметить подход к моделированию электрических цепей, рассмотренный в статьях [13-16]. Там рассмотрен метод составления уравнений нелинейной электрической цепи, в котором уже на этапе выбора переменных уравнений, описывающих электрическую цепь, учитываются условия сходимости итерационного метода Ньютона. Кроме того, данный метод составления уравнений позволяет корректно составить СЛАУ, описывающих линейную (линеаризованную) электрическую цепь. В качестве источника информации целенаправленного выбора переменных используются параметры линейных (линеаризованных) компонент конкретной электрической цепи.
Представляет интерес рассмотрение данного метода составления уравнений электрической цепи применительно к задаче расчета установившихся режимов энергосистемы. Следует отметить, что отличительной особенностью данного метода есть то, что он применим к электрическим цепям, содержащим только двухполюсные компоненты. Поэтому, прежде всего, необходимо объекты энергосистемы представить эквивалентными схемами (моделями), состоящими из двухполюсных компонент.
3. Модели объектов энергосистемы
Согласно [2], [17], к основным объектам энергосистемы, рассматриваемой при расчете установившихся режимов, относятся производители (генераторы) и потребители (нагрузки) электроэнергии, линии передачи электроэнергии. Приведем схемы замещения отдельных объектов энергосистемы и основные соотношения для компонент, входящих в объекты.
Напомним, что участок линии передачи электроэнергии в общем случае состоит из высоковольтной линии электропередачи, трансформатора и высоковольтного выключателя. Принято считать, что в трехфазной энергосистеме значения токов и напряжений рабочих режимов симметричны. Исходя из этого, трехфазную систему заменяют однофазной. В данном случае эквивалентная схема высоковольтной линии электропередачи традиционно содержит линейные емкости Сх, С2, сопротивление Яг и индуктивность Ц (рис. 1а).
Возможны другие варианты представления линии электропередачи, но важно то, что эквивалентная схема линии электропередачи содержит линейные, независимые, напряжениеуправляемые двухполюсные компоненты и их функциональные зависимости можно представить в общем виде как
I = /(й), (1)
где I - ток, протекающий через компоненту, й - напряжение, падающее на компоненте.
Эквивалентная схема высоковольтного выключателя (ключа) приведена на рис. 1б, где в замкнутом состоянии он представлен источником напряжения Е, напряжение которого йЕ = 0 , а в разомкнутом - источником тока J (рис. 1в), ток которого ^ = 0. Иногда в качестве
эквивалентной схемы замещения ключа применяется сопротивление. В этом случае величина сопротивления зависит от состояния ключа. Если ключ находится в разомкнутом состоянии, то величина сопротивления равна сопротивлению утечки ключа, в замкнутом - остаточному сопротивлению замкнутого ключа.
1_1Л л л
с2
Е
7
а) высоковольтная линия передачи;
б) ключ К замкнут: 11Е=0; в) ключ К разомкнут:Ь=0;
Іі = Іт -'2='Ет-„
Л
й2=иЕт
г ) трансформатор
Рис. 1. Эквивалентные схемы объектов линии передачи электроэнергии Обычно для описания работы трансформатора применяют следующие соотношения:
й2 = пй!, і2 = Ї, /п,
(2)
где п - коэффициент трансформации, й!,Ї!, й2, Ї 2 - входные/выходные напряжения и токи трансформатора (рис. 1 г). Эквивалентная схема трансформатора (рис. 1 г) содержит источник напряжения Ет, зависимый от напряжения, падающего на источнике тока ,т, и источник тока ,т,
зависимый от тока источника напряжения Ет. В общем виде функциональные зависимости компонентов трансформатора можно записать как
йБ, = 1(й1г),
і.,, = Г(ЇЕГ>-
Обычно нагрузка определяется активной и реактивной мощностями РН, QН, и соответствующее ей сопротивление эквивалентной схемы имеет следующую зависимость тока Ї от напряжения й:
(Р н йа +Q „ йг) + ](Р н йг - Q н йа)
(3)
Ї = Їа + ІЇг
М 2
(4)
где Їа,Їг - активная и реактивная составляющие комплексного тока Ї, протекающего через нагрузку, йа,йг - активная и реактивная составляющие комплексного напряжения й , падающего на нагрузке, |й| - модуль напряжения й, причем й2 =и2С1 + й2г, у' - мнимая единица. Из (4)
следует, что, компонента нелинейная, независимая, напряжениеуправляемая и в общем виде может быть описана уравнением (1).
Переходя к рассмотрению эквивалентных схем источников электроэнергии, отметим, что если при заданном напряжении, падающем на источнике электроэнергии, последний обеспечивает требуемую энергосистеме активную и реактивную мощности, то он может быть представлен в схеме компонентой типа независимого источника напряжения.
Когда источник электроэнергии вырабатывает заданные активную и реактивную мощности Рг, Qг, то он представляется в схеме компонентой типа сопротивление. Зависимость между током
и напряжением описывается уравнением (4), в котором Рн, Qн заменяются на Рг, Qг, со знаком минус (отдает энергию в энергосистему).
В случае, когда источник электроэнергии вырабатывает заданную активную мощность Рг
при заданном модуле напряжения |й|, падающего на источнике, связь между током I,
вытекающим из источника в энергосистему, и напряжением й, падающем на нем, описывается следующим уравнением:
гг -г Рйа + Qйr .Рйг — Qйa
1= а + ]1г =-----а |2 Г + 3-----Г 2 ' (5)
|й[ |й[
где Р = Рн — Рг есть суммарная мощность, вырабатываемая в узле подключения генератора, а
Q = Qн — Qг - суммарная реактивная мощность. Предполагается, что источник электроэнергии
может содержать активную Рн и реактивную Qн нагрузку, которые известны, а активная йа и
реактивная йг - составляющие комплексного напряжения источника энергии определяются следующим образом:
йа =\й\cosj,
II (6)
йг = й^Пф,
где ф - фазовый угол между активной и реактивной составляющими напряжения. В этом случае
аргументами функциональной зависимости (5) выступают реактивная составляющая источника
электроэнергии Qг и фазовый угол ф. Как следует из (5) и (6), источник электроэнергии может
быть представлен в схеме сопротивлением, функциональная зависимость тока от напряжения которого задана в параметрическом виде. В общем виде функциональная зависимость тока, протекающего через компоненту, имеет следующий вид:
1= /(ф, Qг ) , (7)
а функциональная зависимость напряжения, падающего на компоненте,
й = /(ф). (8)
Следующим шагом есть построение математической модели.
4. Составление уравнений, описывающих энергосистему
Построение математической модели энергосистемы данным методом базируется на основной системе уравнений энергосистемы, куда входят уравнения, составленные на основе законов Кирхгофа, и компонентные уравнения. Для описания графа энергосистемы и соответственно уравнений на основе законов Кирхгофа применяются топологические матрицы контуров и сечений. Переменные системы уравнений выбираются из напряжений и/или токов компонент в результате анализа основной системы уравнений. При этом учитываются параметры компонентных уравнений и особенности топологических матриц, присущие конкретной энергосистеме. В конечном счете из основной системы уравнений выделяется система уравнений, соответствующая выбранным переменным, и система уравнений связи, с помощью которой вычисляются все напряжения и токи компонент. Преобразованная таким образом основная система уравнений рассматривается как математическая модель энергосистемы. Рассмотрим этапы составления математической модели энергосистемы.
Вначале составляется эквивалентная схема замещения энергосистемы, определяются компонентные уравнения и граф эквивалентной схемы. Для этого объекты энергосистемы заменяются эквивалентными схемами, состоящими из двухполюсных компонент, и устанавливаются тип, управляемость и зависимость компонент от других.
Составление топологических матриц контуров и сечений включает выбор дерева графа эквивалентной схемы и составление контуров для выбранного дерева. Дерево графа эквивалентной схемы энергосистемы выбирается таким образом, чтобы все источники напряжения включались в дерево, а все источники тока в хорды. Следует учесть, что при выборе дерева графа эквивалентной схемы в дерево включаются компоненты, которые имеют функциональные зависимости, представленные в параметрическом виде. Составление контуров выполняется присоединением ветвей хорд к дереву. Рассматривается случай, когда в качестве независимых переменных выбраны напряжения ветвей дерева. Как следует из [14, 16], сходимость метода решения при составлении уравнений схемы, а также корректное составление СЛАУ для линейной (линеаризованной) схемы будут учтены в том случае, если в контуре, образованном присоединением хорды к дереву, компоненты дерева будет иметь проводимости по величине больше, чем проводимость этой хорды.
Напряжения и токи компонентов схемы в векторах группируются в элементы, содержащие компоненты, которые входят в дерево (индекс Д), т.е. ветви, и содержащие компоненты, не входящие в дерево (индекс Х) - хорды, таким образом
IД
(9)
где й - вектор напряжений компонент, I - вектор токов компонент.
В связи с тем, что контуры образуются присоединением хорд к дереву графа схемы,
, где 1Х - единичная подматрица хорд, а
й= й Д , I = I Д
йх IХ
1
Х
топологическая матрица контуров имеет вид топологическая матрица сечений - вид \1д —Щ, где 1д - единичная подматрица ветвей, тогда
их =-Рид,
уравнения схемы в матричном виде, составленные на основе законов Кирхгофа, можно записать следующим образом:
7‘тУ (10)
‘Д— X- (11)
Индекс * обозначает транспонированную матрицу. Возможны и другие способы построения контуров, но только при рассмотренном способе в топологических матрицах можно выделить единичные подматрицы.
Элементы ветвей и хорд векторов напряжений и токов компонент (10), (11) и компонентных уравнений группируются в соответствии с типом компонент эквивалентной схемы энергосистемы. В компонентных уравнениях в качестве аргументов указываются только управляемость и
зависимость компонент, так как только они необходимы при анализе основной системы уравнений.
С учетом предложенных преобразований система уравнений (10) принимает следующий вид:
Х^ • Х Р Р* Р* Р XоЕ XоЕт Xо До Хо Дш
и^ = - 3 тЕ ЗтЕт ЗтДо ЗтДф2
Р3Е Р3Ет До Р3Дф2
иЕ
и„
и
До
ит
(12)
а (11)
1е
1 Ет
о Д
Дф2
^ЕХо РЕ3т
ЕЗ
^ЕтХо ^ЕтЗт РЕт3
РдоХо Рдо3т Рдо3
¥д х Рд з ¥д з
Д ф^о Д ф т Д ф2
Х
I,
(13)
Компонентные уравнения будут записаны в следующем виде:
иЕт = /Ет(изт)> 1 Дш = /д„( Ф2 )’
иДо = $Д,ю(ф)’ , 1 Хо = 1хп(ихп)’ (14)
1Д
где индекс Е - указывает, что данный элемент вектора относится к независимым источникам напряжения, Ет к зависимым источникам напряжения эквивалентных схем трансформаторов,
До - к компонентам, входящим в дерево, токи и напряжения которых заданы в параметрическом
виде, До - ко всем остальным компонентам, входящим в дерево, 3 - к независимым источникам
тока, Зт - к источникам тока, зависимым от токов источников напряжений трансформаторов, Хо
- ко всем остальным компонентам, входящим в хорды.
Из анализа систем уравнений (12)-(14) следует, что в случае, если известны векторы
напряжений иЕ, ид ,3 и векторы параметров ф, 0г , можно определить остальные токи и
напряжения компонент, входящих в эквивалентную схему замещения энергосистемы, выполнив следующую последовательность вычислений.
т
З
т
и,
'/Ет (-(1 и
+ + К гг ип + К
иДю))(
и Хе К * К * К* К * Г ХеЕ ХеЕт ХеДе ХеДфО
и и К * К * К * К * г иЕ иЕт ЛЩс ^ иДфО
1тДа Да 1тД фв Дфв
и Е
и Ет и Дс
и
(15)
|фО
I
Дфб = ^Дфа (ф, вг ) -1 Да = IДа (и Да ) -
1 Ха = $Ха (иХа ) -
(16)
I 1т (РЕтХа 1 Ха + КЕт1т 1 1т + КЕт1 1 Л ) ■
(17)
Из анализа системы уравнений (15) следует, что зависимые источники напряжения Ет могут образовывать контуры (циклы) через напряжения источников 1т ■ А из системы уравнений (17) следует, что зависимые источники напряжения 1т могут образовывать контуры (циклы) через источники напряжения Ет ■ Появление таких контуров означает, что напряжения зависимых источников напряжения и токов зависимых источников тока должны определяться в результате решения систем уравнений (15) и (17). Этот вопрос подробно рассмотрен в [15], там же предложен алгоритм выделения таких контуров, определения независимых переменных и составления уравнений, образованных этими контурами. Особенности эквивалентных схем объектов энергосистемы учитываются следующим образом. В случае эквивалентной схемы замещения участка линии передачи электроэнергии, приведенной на рис. 1, такие контуры (циклы) в
уравнениях (15) и (17) не образуются, т.е. элементы матриц ¥Л
*1 т Ет
0 и К
Ет1 т
Тогда, как следует из уравнений (13), выбранным переменным будет соответствовать следующая система уравнений:
а
Дф в
К^ДаХа КДаЛт КДаЛ К^ДФбХа Кдфо1т Кдфо1
I
Х а
1Г
(18)
в которой все токи необходимо представить через напряжения компонент ветвей векторов иЕ,иДв, 11 и векторов параметров ф, вг , используя системы уравнений связи (15)—(17).
Система уравнений (18) нелинейная. Она содержит как комплексные переменные иДд , так и
вещественные ф, вг . Математической моделью энергосистемы будем считать алгоритмы выбора дерева графа схемы, составления топологических матриц, группирования векторов и
т
т
0
т
топологических матриц вместе с системой уравнений (18), в которой переменными являются напряжения иДд компонент ветвей и параметры ф, вг генераторов, и систем уравнений связи (15)—(17).
Коммутация линий электропередачи или изменение условий производства электроэнергии в реальной энергосистеме учитывается в математической модели энергосистемы путем изменения параметров функциональных зависимостей или типов компонент. Эти изменения отслеживаются перевыбором дерева, пересоставлением топологических матриц контуров и сечений и перегруппировкой векторов напряжений и токов компонент в соответствии с новым набором типов компонент.
5. Модификация метода Ньютона для решения СНУ, составленных предложенным методом
Как уже упоминалось, в основных программных комплексах для решения СНУ применяется итерационный метод Ньютона. Суть метода заключается в следующем. Ограничимся рассмотрением тех эквивалентных схем энергосистемы, которые при составлении уравнений не
образуют циклы (контуры) в уравнениях (15), (17), т.е. К
0 и К
0 . Тогда СНУ (18) в
общем виде можно записать как
ї(ид ) = 0, (19)
где ид = р‘Дв ,р‘- вектор независимых переменных. Методом Ньютона принято считать следующий вычислительный итерационный процесс решения СНУ (19) [9]:
идп+1= идп - а^ид , (20)
где а - демпфирующий множитель, ЛИд - вектор приращений переменных ид на п -ом
итерационном шаге. Вектор Лид определяется из решения СЛАУ
/'(ид п) Аид п = /(ид п)’ (21)
где /(идп) - вектор невязок на п -ом итерационном шаге решения СНУ, / (идп) - матрица
частных производных СНУ (19), и имеет следующий вид:
/11 (и д п ) ■'
^(ид п )
/1т (и д п )
(22)
!ш! (иД п ) /шгп (иД п )
где т - порядок СНУ (19), /у(идп ) - частная производная г -ого уравнения по у -ой переменной.
Как следует из приведенного итерационного процесса, вопрос корректной постановки задачи может возникнуть как при решении СНУ (19) (сходимость итерационного процесса решения), так и при решении СЛАУ (21) (плохая обусловленность, система большой размерности). Это означает, что метод решения СНУ (19) и составление уравнений предлагаемым методом необходимо рассматривать как единый вычислительный процесс, так как в процессе решения
Ти Т
необходимо проверять задачу на корректность и соответствующим образом ее корректировать. Следует заметить, что в приведенном методе Ньютона рассматривается решение СНУ (19) как чисто математическая задача с полным абстрагированием от физической сущности задачи. В то же время в предложенном методе составления уравнений для корректной постановки задачи моделирования необходимо знать особенности конкретной задачи. Это означает, что метод Ньютона необходимо модифицировать таким образом, чтобы можно было применить предложенный метод составления уравнений в процессе решения.
Известно, что физической интерпретацией СЛАУ (21) электротехнической задачи есть система уравнений, описывающих линеаризованную эквивалентную схему этой задачи, т.е. схему, у которой функциональные зависимости заменяются их производными. Предлагается СЛАУ (21) получать путем составления уравнений для линеаризованной схемы. В процессе составления уравнений для линеаризованной схемы реализуется контроль задачи на корректность, и в случае необходимости выполняется перевыбор независимых переменных. Такой подход позволит рассматривать составление уравнений, описывающих схему, и их решение как единый вычислительный процесс решения задачи, что является существенным отличием предложенного подхода к моделированию энергосистемы от реализованного в программных комплексах.
Система уравнений связи (15)-(17) и система уравнений (18) применяются для составления уравнений, описывающих линеаризованную эквивалентную схему энергосистемы. При этом
следует учесть, что в системе уравнений (15)-(18) напряжения компонент и и токи компонент I
исходной схемы заменяются напряжениями компонент Аи и токами компонент Ы линеаризованной схемы. Вместо функциональных зависимостей компонент исходной схемы используются функциональные зависимости компонент линеаризованной схемы, которые имеют следующий вид:
АиДф = /Дф (фп) Аф ,
ыиЕт = /Ет(и 1п)ыит,
А 1 ф = /'«вф (фп ’ вгп )Аф + 1'Дфвв
(фп, вг )Авг ,
ЫДа = 1 Да (иДап ) Аи Да , (23)
ЫХа = 1 Ха (иХап )АиХа ' Ы- 1т 1 1т (Етп )Ы- Ет ,
где /i ( ) - численное значение производной функциональной зависимости компоненты схемы на
п -ом итерационном шаге решения СНУ. Как следует из (23), производные функциональных зависимостей компонент схемы имеют размерность проводимостей или коэффициентов усиления.
Так же, как для исходной схемы, в качестве независимых переменных системы уравнений, описывающей линеаризованную схему, используются векторы напряжений компонент линеаризованной схемы, входящих в ветви дерева Аидд , и векторы фазовых углов А ф и
реактивных мощностей линеаризованных компонент генераторов. Переменные АйДс , Ар ,
определяются из решения системы уравнений (18), которые в случае линеаризованной схемы имеют следующий вид:
АІ Д 0 Рдаха Рдг;т АІХ Х г
А і Д*о ^Д^ха Рд^т Аі, ит
= /(и Д.) ■
(24)
Учитывая (15)-(17) и (23), систему уравнений (24) нетрудно преобразовать к виду (21).
Для контроля корректности составленной СЛАУ (24) применяется следующее выражение:
АIд, - (Рдехе А■ Хе + ^ДеитА■ ит), (25)
которое, учитывая (15)—(17) и (23), нетрудно преобразовать к виду
і=і
(26)
где к - количество ветвей хорд, входящих в сечение I -ой ветви дерева, т - количество ветвей дерева в схеме, Од& - проводимость компоненты I -ой ветви дерева графа схемы, ОХ^ -
проводимость компоненты у -ой ветви хорды графа схемы, входящей в I -ое сечение. СЛАУ (25) будет считаться корректно составленной, если выполняется следующее неравенство:
°д„ г Сх,, . ОТ
В противном случае выполняются новый выбор переменных АйдС, пересоставление топологических матриц Г, Г и пересчет п -ого итерационного шага решения СНУ.
6. Тестовый пример
С целью оценки эффективности предложенного метода корректной формулировки задачи моделирования установившихся режимов энергосистемы рассмотрим тестовый пример. Тестовый пример будет промоделирован программой, в которой в первом случае составление уравнений
выполняется методом
узловых потенциалов, а во втором -
предл оже нным м етодом с контролем на
корректность в процессе решения задачи.
Количество итераций, затраченных на
получение результата,
будет служить критерием эффективности. В качестве тестового примера выбран фрагмент эквивалентной схемы энергосистемы, показанной на рис. 2.
а) фрагмент моделируемой сети;
Рис. 2. Тестовый пример
В табл. 1 приведены параметры генераторов и нагрузок схемы, в табл. 2 - параметры компонент линий передачи.
Таблица 1. Параметры генераторов и нагрузок Таблица 2. Параметры компонент линий схемы передачи
№ Имя компонента Рн Рг вн вг м иа
1 Е1 200
2 300 100 90 - 220
3 200 - 200 - -
4 К3 50 - 50 - -
5 Р4 100 50 50 220
№ Имя компонента Я10 Я2 А
1 У-1 4,3 1Е+10 2,6 -27Е+3
2 У-2 52,6 1Е+10 31,6 -337Е+3
3 У-3 15,5 1Е+10 10,8 -410Е+3
4 У-4 59,5 1Е+10 58,0 -650Е+3
5 У1-5 59,5 1Е+10 58,0 -650Е+3
В тестовом примере в качестве изменяемого параметра выбрано сопротивление Я1 линии электропередачи. Изменение сопротивления Я1 можно интерпретировать как переключение высоковольтного переключателя из закрытого состояния (сопротивление велико) в открытое (сопротивление мало). При моделировании сопротивление Я1 изменялось следующим образом:
Я1 = Я1010-п, где 1 < п < 7 .
На рис. 3 приведены результаты моделирования. Как следует из результатов
моделирования тестового примера, итерационный процесс решения системы уравнений при п = 7 , составленной методом узловых потенциалов, оказался расходящимся, в то время как система имеет решение. На решение системы уравнений, составленной предложенным методом, было затрачено меньше 15 итераций.
метод узловых потенциалов —«—предлагаемый метод
94
1/
/
/
' 7 7
1 2 3 4 5 6 7
и
Рис. 3. График зависимости количества итераций, затраченных на решение, в зависимости от метода составления уравнений, описывающих тестовый пример
6. Заключение
В связи с тем, что применение моделирования становится неотъемлемой частью анализа, планирования и оперативного управления электроэнергетическими системами, достоверность результатов моделирования является основным критерием оценки результатов моделирования. Это требует новых подходов к разработке методов и алгоритмов описания электроэнергетических систем и решения составленных описаний.
В данной работе на примере расчета установившихся режимов энергосистемы показано, что можно корректно поставить задачу моделирования на этапе составления уравнений. Корректная постановка задачи достигнута за счет выбора в качестве независимых переменных системы уравнений, описывающих энергосистему, напряжений компонент ветвей дерева графа эквивалентной схемы энергосистемы. В качестве критерия выбора переменных используются параметры функциональных зависимостей компонент эквивалентной схемы, т.е. при составлении системы уравнений, описывающих схему, учитываются конкретные особенности схемы.
Предложена модификация метода Ньютона решения СНУ, описывающих энергосистему, которая позволяет реализовать в процессе решения проверку корректности моделируемой задачи и в случае необходимости вновь корректно сформулировать задачу моделирования.
Объединение задачи составления уравнений, описывающих энергосистему, и их решение дает возможность моделировать энергосистемы с динамически изменяемыми структурами.
СПИСОК ЛИТЕРАТУРЫ
1. Программно-вычислительный комплекс АНАРЭС-2000 // http://anares.irk.ru/PDF/ar10801.pdf .
2. Модели анализа нормальных и предельно установившихся режимов комплекса АНАРЭС для ІВМ РС / О.М. Кучерова, Ю.Н. Кучеров, Щ.Н. Шепилов и др. - К.: Министерство энергетики и электрификации Украины, 1993. - 25 с.
3. Программный комплекс Епє^уСБ для проектирования электроэнергетических систем // Іійр://зд^. csoft.ru/press/pub/article_23506.html.
4. Прихно В. Л. Программный комплекс Космос оперативных расчетов режимов энергосистем на основе телеметрической информации // http://www.regimov.net/books/prihno_articles/01.html.
5. http://en-soft.ru/.
6. http://niieng.ru/.
7. http://www.digsilent.de/.
8. Конторович А.М., Дунаева Н.П. Исследование методов расчета установившихся режимов, основанных на разложении решения в ряд Тейлора. - Иркутск, 1978.
9. Ортега Д., Рейнболт В. Итерационные методы решения нелинейных систем уравнений со многими неизвестными. - М.: Мир, 1975. - 558 с.
10. Тихонов Д.Н. О решении некорректно поставленных задач и методе регуляризации // Доклады АН СССР. -1963. - № 151. - С. 501 - 504.
11. Лаврентьев М.М. О некоторых некорректных задачах математической физики. - Новосибирск: Изд-во СО АН СССР, 1962. - 92 с.
12. Иванов В.К. О некорректно поставленных задачах // Математический сборник. - 1963. - № 61 (103):2. -С. 211 - 223.
13. Волобоев В.П. Составление уравнений цепи, содержащей зависимые двухполюсники и многополюсники // Вопросы проектирования математических машин и устройств. Научный совет по кибернетике АН УСАН УССР, Киев, 1972. - С. 3 - 16.
14. Волобоев В.П. К учету сходимости численных методов при составлении уравнений цепи постоянного тока // Вопросы проектирования математических машин и устройств. Научный совет по кибернетике АН УСАН УССР. - Киев, 1972. - С. 17 - 26.
15. Волобоев В.П. О расширении класса схем, моделируемых методом напряжений ветвей дерева // Проектирование технических средств ЭВМ и систем: Сб. научн. тр. - К.: Изд-во ИК АН УССР, 1982. - С. 32 -36.
16. Волобоев В.П., Клименко В.П. Об одном подходе к моделированию сложных систем // Математичні машини і системи. - 2008. - № 4. - С. 111 - 122.
17. Сенди К. Современные методы анализа электрических систем. - М.: Энергия, 1971. - 360 с.
Стаття надійшла до редакції22.09.2009
