Винокур В. М., Кавалеров Б. В., Шигапов А. А.,
г. Пермь

УДК 621.311.001.57

АЛГОРИТМИЧЕСКАЯ БАЗА ДЛЯ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ
СИСТЕМ ЭЛЕКТРОСНАБЖЕНИЯ ПРОИЗВОЛЬНОЙ КОНФИГУРАЦИИ

(презентация)

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

Задачи математического моделирования

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

В зависимости от потребностей исследователей, требований заказчика алгоритмическая база позволяет создавать компьютерные программы различного назначения. Перечислим лишь неко-торые из них:

Алгоритмическая база предоставляет разработчику весь необходимый арсенал методов и средств, реализованный в соответствии с форматом заказчика и основанный на строго обоснованных критериях практической эксплуатации в условиях промышленного производства. Алгоритмическая база «МЭЭС» формирует гибкую систему, которая может автоматически настраиваться на конкретные потребности заказчика и модифицировать свою структуру в соответствии с текущими ситуациями. Сегодня мы имеем несколько форматов, которые предложили заказчики, данные форматы легли в основу поля возможностей, формируемого при создании инструментальных средств и расширяемого по мере развития алгоритмической базы. Одним из главных требований являются адаптируемость и возможность к дальнейшему развитию, к расширению класса используемых методов и их свободному включению в инструментарий.

Разработка алгоритмической базы потребовала серьезной научно-исследовательской работы. В результате на сегодняшний день мы имеем необходимый арсенал методов и алгоритмов для компьютерного моделирования произвольных электросистем. Охарактеризуем объект исследования.

Основные принципы построения алгоритмической базы

Среди особенностей, присущих электроэнергетическим объектам, в первую очередь следует отметить значительную сложность большинства из них, которая проявляется в многоэлементности и иерархичности структуры, обилии степеней свободы, разнообразии параметров, характеризующих состояние [1]. Наличие у электроэнергетических объектов указанных признаков позволяет причислить их к объектам кибернетического типа [2]. Кроме того, необходимо отметить, что стандартные технические решения, обычно применяемые для электростанций больших и средних мощностей, оказываются неприемлемыми для малых (мини-энергосистем). В то же время каких-либо специальных норм технологического проектирования для них не существует, что требует от проектировщика творческого подхода к решению возникающих технических проблем [3]. В связи с этим задачи, которые необходимо решать при эксплуатации существующих, а также при создании новых образцов, оказываются весьма сложными. Решение таких задач на интуитивном уровне недопустимо.

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

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

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

Энергосистема является сложной системой. Основными отличительными признаками сложной системы являются следующие признаки [4].

  1. Наличие большого количества взаимно связанных и взаимодействующих между собой элементов.
  2. Сложность функций, выполняемых системой и направленных на достижение заданной цели функционирования.
  3. Возможность разбиения системы на подсистемы, цели функционирования которых подчинены общей цели функционирования всей системы.
  4. Наличие управления.

Процесс функционирования сложной системы можно представить как совокупность действий и элементов, подчиненных единой цели.

Математическая модель сложной системы состоит из математических моделей элементов и математической модели взаимодействия между элементами [4].

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

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

Была решена задача структурной декомпозиции. Выбор уровня декомпозиции осуществлялся, исходя из следующих соображений. Целесообразно, как уже отмечалось, получить по возможности наиболее универсальную математическую модель, ориентированную на широкий класс электроэнергетических объектов. Возможность получения такой модели обусловлена тем, что электроэнергетические объекты состоят в основном из ограниченного набора типовых элементов. К их числу относятся вращающиеся электрические машины, трансформаторы, линии электропередачи, реакторы, батареи конденсаторов и т.п. Таким образом, наиболее целесообразным и естественным представляется выбрать в качестве структурных элементов мини-энергосистем вышеперечисленные наиболее распространенные типовые элементы электрических систем и сетей.

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

Обобщенные модели структурных элементов

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

  (1)

где I, pI – вектор статорных токов и вектор производных статорных токов элемента; U – вектор напряжений, приложенных между внешними зажимами элемента; A, B – матрицы, размерность которых зависит от системы координат в которых моделируется структурный элемент (d,q или a,b,c), а также от того, полные это уравнения или упрощенные; H – вектор, определяющий воздействие на элемент со стороны средств регулирования электрических параметров; знак «–» в уравнении (1) означает, что структурный элемент является источником электрической мощности (генератором), знак «+» – потребителем; компоненты вектора H для пассивных элементов электрической системы – нулевые, для синхронных машин при отсутствии РВ (регулирования возбуждения) – постоянные, при наличии РВ – меняются на каждом шаге расчета.

Уравнения (1) решаются на каждом шаге расчета методом численного интегрирования относительно внешних токов структурного элемента. При этом следует отметить, что для пассивных элементов выражение (1) содержит всю систему дифференциальных уравнений элемента. Для электрических машин уравнения (1) следует интегрировать совместно с уравнениями роторных контуров и уравнениями, описывающими механическое состояние элемента. Таким образом, для электрических машин вектора I и pI не совпадают по размерности: вектор Iсодержит не только внешние, но и внутренние (роторные) токи, расчет которых производится на каждом шаге численного интегрирования полной системы дифференциальных уравнений электрической машины. Знаки при слагаемых уравнения (1) соответствуют нормальному направлению токов в электрических машинах (от генератора к узлу).

Взаимодействие структурных элементов

После приведения математических моделей всех элементов к единой обобщенной форме записи (1) строится математическая модель взаимодействия структурных элементов.

Электрическая система изображается в виде совокупности соединенных между собой двухполюсников (ветвей), то есть элементов с двумя зажимами. Наиболее общим понятием, соответствующим представлению системы в виде соединения ветвей, является понятие графа [5]. Представление системы электроснабжения в виде схемы является полным, то есть дает достаточную информацию для решения задачи моделирования, если оно содержит как соединения входящих в нее элементов, так и характеристики этих элементов:

На основании анализа возможных методов моделирования взаимодействия элементов [7] выбор был сделан в пользу методов, основанных на решении уравнений связи.

Для построения алгоритма решения уравнений связи и автоматизации построения расчетных схем электрических систем используется матрично-топологическое направление теории цепей.

Проиллюстрируем получение уравнений связи электрически взаимодействующих элементов. Принцип формирования систем уравнений связи, как отмечается в [2], базируется на условии инвариантности граничных переменных, согласно которому в общем для элементов узле в произвольный момент времени каждая скалярная граничная переменная, определяемая в одном элементе, тождественно равна себе самой, определенной в любом другом из этих элементов, а алгебраическая сумма векторных граничных переменных при тех же условиях равна нулю. При этом скалярной является переменная, не имеющая направления (в нашем случае – потенциал), а векторной – граничная переменная, характеризующаяся направлением (в нашем случае – электрический ток). Уравнения связи могут быть составлены как баланс активной и реактивной мощностей в узловой точке [6]:

, , (2)

что равноценно уравнению 1-го закона Кирхгофа для узловой точки.

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

(3)

Аналогично составляются уравнения связи и для более сложных структурных схем.

На основании системы уравнений (3) формируется векторно-матричное уравнение для расчета узловых напряжений на каждом шаге расчета:

, (4)

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

У каждого генератора система d, q координат (Парка-Горева) связана с его ротором и поэтому при их согласовании необходимо использовать следующие матрицы преобразований. Матрица преобразования:

  , (5)

где δ ij – угол между осями Парка-Горева двух различных генераторов.

Из структуры матрицы инцидентности и структуры матрицы видно, что их элементы определяются конфигурацией системы и выбором систем координат, к которым выполняется приведение. Вектор НY  и матрица Y определяются составом элементов системы.

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

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

Основные вычислительные трудности при реализации алгоритма (1) связаны с обращением матрицы , размер которой равен 2n, где n – число узлов моделируемой системы электроснабжения.

Использование матричных операций позволяет создать алгоритмы цифрового моделирования сложных переходов в любой их последовательности для самых различных эквивалентных схем мини-энергосистем.

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

Созданная алгоритмическая база «МЭЭС» обеспечивает формирование различных вариантов компьютерных программ для моделирования мини-энергосистем произвольной конфигурации и состава структурных элементов. Реализуется моделирование всех основных режимов работы, в том числе, статических, динамических, несимметричных.

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

Литература

  1. Коротков Б.А., Попков Е.Н. Алгоритмы имитационного моделирования переходных процессов в электрических системах. – Л.: Изд-во Ленинградского университета, 1987. – 280 с.
  2. Веников В.А. Кибернетические модели электрических систем. – М.: Энергоатомиздат, 1982. – 328 c .
  3. Бут Д.А. Синтез автономных электроэнергетических систем // Электричество. – 1994. – № 1. – С. 1–17.
  4. Бусленко Н.П. Моделирование сложных систем. – М.: Наука, 1978. – 440 с.
  5. Оре О. Теория графов. – М.: Наука, 1980. – 336 с.
  6. Веретенников Л.П. Исследование процессов в судовых электроэнергетических системах. Теория и методы. – Л.: Судостроение, 1975. – 376 с.
  7. Кавалеров Б.В. Выбор способа моделирования взаимодействия структурных элементов системы электроснабжения // Наука, техника, инновации: Материалы докладов региональной научной конференции 05 – 08 декабря 2002 г., ч.2 – Новосибирск, 2002. – С. 137–138.