Сообщество Engee

Методика параметризации синхронных машин

Автор
avatar-gladkovgladkov
Notebook

Методика параметризации моделей синхронных машин

Проблематика

В российской электроэнергетике традиционно используется справочник Неклепаев Б.Н., Крючков И.П. "Электрическая часть станций и подстанций". Кроме всего прочего, он содержит основные параметры гидро- и турбогенераторов отечественного производства.

Современные математические модели в программных комплексах (Engee, RSDAD, PSCAD, PowerFactory, PSS/E, MATLAB/Simulink) требуют гораздо большего набора параметров, чем указано в справочнике. Хотя существуют инструменты по типу Block Parameterization Manager с различными пресет-моделями блоков, однако они обладают рядом недостатков, например, если параметра нет в справочнике, то он остается дефолтным.

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

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

Описание методики

Таблица параметров, которые нужны для параметризации математических моделей cинхронных машин с неявнополюсным (Т) и явнополюсным (Г) роторами, с указанием способа из нахождения:

Параметр Обозначение Как найти?
1 Номинальная полная мощность [1]
2 Номинальное напряжение [1]
3 Номинальная электрическая частота [1]
4 Номинальная скорость вращения [1]
5 Ток цепи возбуждения [1]
6 Напряжение цепи возбуждения [1]
7 Активное сопротивление статора [1] или для Г [2, с.40]
8 Реактивное сопротивление рассеяния статора [1] или для Г [2, с.40]
9 Синхронное реактивное сопротивление по оси d [1]
10 Синхронное реактивное сопротивление по оси q [1] или для Т [2, с.29], для Г [2, с.40]
11 Реактивное сопротивление нулевой последовательности [1] или для Г [2, с.40]
12 Переходное реактивное сопротивление по оси d [1]
13 Переходное реактивное сопротивление по оси q для Т [2, с.34]
14 Сверхпереходное реактивное сопротивление по оси d [1]
15 Сверхпереходное реактивное сопротивление по оси q для Т [2, с.34], для Г [2, с.40]
16 Переходная постоянная времени по оси d при разомкнутой обмотке статора [1] или для Г [2, с.40]
17 Сверхпереходная постоянная времени по оси d при разомкнутой обмотке статора [3, с.48] или для Г [2, с.39]
18 Переходная постоянная времени по оси q при разомкнутой обмотке статора для Т [4, с.125]
19 Переходная постоянная времени по оси d при замкнутой накоротко обмотке статора [1] или [3, с.48]
20 Сверхпереходная постоянная времени по оси d при замкнутой накоротко обмотке статора [1] или [3, с.48]
21 Переходная постоянная времени по оси q при замкнутой накоротко обмотке статора для Т [2, с.39]
22 Сверхпереходная постоянная времени по оси q при замкнутой накоротко обмотке статора для Т [2, с.39], для Г [2, с.39]

Список литературных источников

[1] Неклепаев Б.Н., Крючков И.П. Электрическая часть станций и подстанций: Справочные материалы для курсового и дипломного проектирования. – М.: Энергоатомиздат, 1989.

[2] Kimbark E.W. Power System Stability. Volume III: Synchronous Machines. – John Wiley & Sons, 1956.

[3] ГОСТ 10169-77. Генераторы электрические синхронные трехфазные с частотой вращения 3000, 1500, 1000, 750, 600, 500 и 375 об/мин. Общие технические условия. – М.: Издательство стандартов, 1977.

[4] Anderson P.M., Fouad A.A. Power System Control and Stability. 2nd Edition. – IEEE Press, 2003.

Скрипт для параметризация синхронных машин

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

Импорт json-файла с базой данных синхронных машин из справочника Неклепаева:

In [ ]:
using JSON
db = JSON.parsefile("synchronous_machines.json")["SynchronousMachines"];

Выбор синхронной машины из списка доступных:

In [ ]:
println("\nСписок синхронных машин:\n")
for (i, m) in enumerate(db)
    println("[$i] $(m["machine_model"]) ($(m["machine_type"]))")
end

print("\nВведите номер выбранной машниы: ")
idx = parse(Int, readline())
m = db[idx]

println("\nВыбрана: $(m["machine_model"])")
println("Тип: $(m["machine_type"])\n")
Список синхронных машин:

[1] Т-6-2у3 (Turbo)
[2] Т-12-2у3 (Turbo)
[3] ВГС260/99-10 (Hydro)
[4] ВГС325/89-14 (Hydro)

Введите номер выбранной машниы: stdin>  1

Выбрана: Т-6-2у3
Тип: Turbo

Подключаем и выполняем функцию для расчета параметров по описанной методике для блоков синхронных машин

In [ ]:
include("convert_parameters.jl");
calculated_parameters = convert_parameters(m)
Out[0]:
Dict{Symbol, Float64} with 25 entries:
  :V_rated      => 6.3
  :S_rated      => 7.5
  :T_p_d        => 0.77
  :X_d          => 1.651
  :X_q          => 1.651
  :T_pp_q0      => 0.0063
  :machine_type => 1.0
  :V_fd_base    => 140.0
  :w_rated      => 3000.0
  :T_pp_d       => 0.0063
  :X_l          => 0.147
  :X_p_q        => 0.9109
  :X_0          => 0.057
  :T_p_q        => 1.155
  :I_fd_base    => 249.0
  :X_pp_q       => 0.9109
  :R_a          => 0.0124
  :X_pp_d       => 0.1208
  :T_p_q0       => 2.09343
  :f_rated      => 50.0
  :T_p_d0       => 7.45
  :T_pp_q       => 0.0063
  :T_pp_d0      => 0.00890762
  :N_pole_pairs => 1.0
  :X_p_d        => 0.1708

Подключаем и выполняем функцию для вывода параметров блоков синхронных машин

In [ ]:
include("print_parameters.jl");
printed_parameters = print_parameters(calculated_parameters);
============================================================
ПАРАМЕТРЫ СИНХРОННОЙ МАШИНЫ
Тип: Турбогенератор (RoundRotor)
============================================================

1. Номинальные параметры:
   S_rated      = 7.5 МВА
   V_rated      = 6.3 кВ
   f_rated      = 50.0 Гц
   N_pole_pairs = 1.0

2. Параметры цепи возбуждения:
   I_fd_base = 249.0 А
   V_fd_base = 140.0 В

3. Основные сопротивления:
   R_a = 0.0124 о.е.
   X_l = 0.147 о.е.
   X_d = 1.651 о.е.
   X_q = 1.651 о.е.
   X_0 = 0.057 о.е.

4. Переходные и сверхпереходные сопротивления:
   X_p_d  = 0.1708 о.е.
   X_pp_d = 0.1208 о.е.
   X_p_q  = 0.9109 о.е.
   X_pp_q = 0.9109 о.е.

5. Постоянные времени по оси d:
   T_p_d0  = 7.45 с (ХХ)
   T_p_d   = 0.77 с (КЗ)
   T_pp_d0 = 0.008907615894039735 с (ХХ)
   T_pp_d  = 0.0063 с (КЗ)

6. Постоянные времени по оси q:
   T_p_q0  = 2.0934295751454606 с (ХХ)
   T_p_q   = 1.155 с (КЗ)
   T_pp_q0 = 0.0063 с (ХХ)
   T_pp_q  = 0.0063 с (КЗ)
============================================================

Выводы

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

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

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