Електротехнiчнi комплекси та системи. Силова електронка
УДК 621.314 doi: 10.20998/2074-272X.2019.5.04
Г.Г. Жемеров, Д.С. Крылов, А.В. Машура
ЭНЕРГОЭФФЕКТИВНОСТЬ СИСТЕМЫ ЭЛЕКТРОСНАБЖЕНИЯ МЕТРОПОЛИТЕНА С РЕКУПЕРАЦИЕЙ ЭЛЕКТРОЭНЕРГИИ ПРИ ТОРМОЖЕНИИ
Мета. Метою статтг е оцшка ККД системи електропостачання метропол1тену, в якш використовуеться чоти-рьохквадрантний прив1д пост1йного струму з рекуперацею енерги в мережу живлення в режим1 гальмування. Методика. Для проведення досшджень використовувалася теор1я електричних кщ математичне моделювання в пакет1 Matlab. Результаты. Отримана теоретична залежшсть ККД СЕ з двонаправленим потоком енерги в1д коефщента резистивного короткого замикання на клемах навантаження. Теоретичний результат перев1рений моделюванням Наукова новизна. Розроблена екв1валентна схема системи електропостачання метрополтену з чотирьохквадрант-ним приводом постшного струму i можлив1стю рекуперацй енерги в мережу живлення в режим1 гальмування, визна-чеш и параметри, заданий графЫ руху електропогзда. Практичне значення. Використання отриманих залежностей i результат1в моделювання дозволить визначити напрямок перспективного розвитку системи електропостачання метрополтену, оптим1зувати иенергоефективтсть. Бiбл. 8, табл. 1, рис. 5.
Ключовi слова: система електропостачання, енерпя, коефщент повернення енергй, коефщент корисно" дй, рекупера-щя енергй".
Цель. Целью статьи является оценка КПД системы электроснабжения метрополитена, в которой используется четырёхквадрантный привод постоянного тока с рекуперацией энергии в питающую сеть в режиме торможения. Методика. Для проведения исследований использовалась теория электрических цепей, математическое моделирование в пакете Matlab. Результаты. Получена теоретическая зависимость КПД СЭ с двунаправленным потоком энергии от коэффициента резистивного короткого замыкания на клеммах нагрузки. Теоретический результат проверен моделированием. Научная новизна. Разработана эквивалентная схема системы электроснабжения метрополитена с четырёхквадрантным приводом постоянного тока и возможностью рекуперации энергии в питающую сеть в режиме торможения, определены её параметры, задан график движения электропоезда. Практическое значение. Использование полученных зависимостей и результатов моделирования позволит определить направление перспективного развития системы электроснабжения метрополитена, оптимизировать её энергоэффективность. Библ. 8, табл. 1, рис. 5.
Ключевые слова: система электроснабжения, энергия, коэффициент возврата энергии, коэффициент полезного действия, рекуперация энергии.
Введение. Функционирование транспортной системы современного крупного города невозможно без использования метрополитена, обеспечивающего значительную часть пассажирских перевозок. Его высокая надёжность сочетается с не самой высокой энергоэффективностью, что обусловлено применением коллекторного электропривода постоянного тока последовательного возбуждения без возможности возврата энергии в питающую сеть. Повышению коэффициента полезного действия (КПД) системы электроснабжения метрополитена посвящено много научных трудов [1-3]. Одним из решений проблемы энергосбережения является использование четырёхквад-рантного электропривода постоянного тока, что даёт возможность организовать двунаправленный поток электрической энергии между источником и нагрузкой. Это позволит выводить накопленную в движущемся составе энергию обратно в промышленную сеть трёхфазного переменного тока, что, в свою очередь, должно повысить КПД всей системы электроснабжения метрополитена. Однако, как показано в [3], эффект энергосбережения от использования рекуперативного торможения не всегда очевиден. Он зависит от конфигурации используемой системы электроснабжения и режимов работы электропривода. При определённых условиях возможен эффект снижения суммарного КПД системы вследствие возникновения дополнительных потерь при возврате энергии в сеть.
Целью работы является оценка КПД системы электроснабжения метрополитена, в которой исполь-
зуется четырёхквадрантный привод постоянного тока с рекуперацией энергии в питающую сеть в режиме торможения.
Структура тяговой подстанции. Традиционная тяговая подстанция для преобразования переменного напряжения в постоянное использует неуправляемые диодные выпрямители, что не позволяет осуществить возврат энергии в питающую сеть. Для реализации возможного повышения эффективности схемы вместо диодных мостов необходимо использовать четырех-квадрантный тиристорный выпрямитель, представленный на рис. 1.
Сеть 6(10) kV представлена трехфазной симметричной системой синусоидальных напряжений uSA, uSB, uSC. Параметры сети учтены активным сопротивлением R0. Параметры линии, соединяющей тяговую подстанцию и преобразовательный трехфазный трансформатор 6(10) kV / 0.71 kV, определяются активным сопротивлением Rj. Сетевые обмотки трансформатора T1 подключаются к сети 6(10) kV, а вентильные - к шестипульсному четы-рехквадрантному мостовому выпрямителю VS1 -VS12. Потери в тиристорном мосте представлены эквивалентным источником противо-ЭДС уровня 1 V в прямом и обратном направлениях и перенесены в цепь постоянного тока. Параметры линии от трансформатора T1 до выпрямительных мостов соответствуют активным сопротивлением R2. Нагрузка представлена двигателем постоянного тока (ДПТ)
© Г.Г. Жемеров, Д.С. Крылов, А.В. Машура
с независимым возбуждением. Параметры линии от управляемого выпрямителя (УВ) до двигателя постоянного тока учтены сопротивлением контактного рельса Якя. Индуктивности в линии питания,
присутствующие там по факту, не участвуют в формировании потерь при передаче энергии, поэтому они перенесены в нагрузку и в эквивалентной схеме они не показаны.
Рис. 1. Эквивалентная схема СЭ МП с четырехквадрашным УВ
Для адекватной оценки энергоэффективности системы электроснабжения метрополитена (СЭ МП) необходимо знать график движения поездов, который, согласно [1-3] содержит следующие интервалы: интервал разгона от нулевой до номинальной скорости (время tac) в среднем составляет 20-30 s; время торможения от номинальной до нулевой скорости (tbr) в среднем равно 40-50 s; интервал стоянки поезда (время tst) обычно равен 25 s; интервал движения с номинальной скоростью (tmov) составляет 110-130 s. Учитывая, что время движения подвижного состава между двумя станциями в среднем составляет три минуты [1], в соответствии с [3] примем следующие значения интервалов движения: tac = 25 s, tmov = 115 s, tbr =45 s, tst = 25 s.
График изменения тока, напряжения и мощности нагрузки для указанных интервалов движения состава в рассматриваемой СЭ МП может иметь вид, представленный на рис. 2.
Uи
h h 0,41н
0.11, н
-Il н
Pi h
0,4Pl m 0,1 Pl h
-Pl h
V
s
<-*r+-
Lc
'tumi, Л
tbr
П7
.1
s
-к
s
Рис. 2. График движения состава между станциями
В интервале разгона поезда а система автоматического регулирования линейно изменяет скорость состава от нулевой до номинальной, что, при незави-
симом возбуждении МПТ, соответствует линейному росту напряжения нагрузки от нулевого до номинального значения. Ток нагрузки при этом ограничивается на номинальном уровне. Развиваемая на этом этапе мощность также линейно нарастает до номинального значения. В интервале движения tmov с номинальной скоростью, к нагрузке будет приложено номинальное напряжение, а поезд преодолевает лобовое сопротивление воздуха и силу трения, развивая до 30-50 % от номинальной мощности. На этом же уровне будет находиться и потребляемый нагрузкой ток. В режиме торможения система автоматического регулирования обеспечивает плавное линейное снижение скорости до нулевого уровня за время tbr. Напряжение нагрузки также будет линейно снижаться до нулевого уровня. При этом для обеспечения возврата в питающую сеть запасённой поездом кинетической энергии, необходимо обеспечить смену полярности тока нагрузки с ограничением его на уровне, не превышающем номинального значения. С началом рекуперации ток нагрузки переходит на реверсивную вентильную группу и поддерживается отрицательным вплоть до момента достижения равенства возвращаемой и потребляемой составом энергии на этапе торможения. После этого изменяющаяся по линейному закону мощность нагрузки снова станет положительной, а ток нагрузки переходит на положительную вентильную группу в результате переключения мостов выпрямителя. В интервале стоянки tst электропривод поезда энергии не потребляет. На всех этапах движения состав потребляет энергию собственных нужд, которая расходуется на обогрев, освещение и вентиляцию вагонов, её величина может доходить до 10 % от номинальной. Это учтено на графиках, представленных на рис. 2.
Развиваемая на этапах движения мощность зависит от физических параметров поезда, от его скорости и массы. Масса подвижного состава определяется числом вагонов и количеством пассажиров в каждом вагоне. В соответствии с [4], поезд состоит из пяти вагонов массой 33 t каждый. В вагон помещается от 200 до 300 пассажиров со средним весом 60-70 kg. Таким образом, можно считать, что масса поезда составляет 200-250 t. Номинальная скорость движения равна 25 m/s или 90 km/h. Согласно [4], действующие на линиях метрополитена поезда оснащены электро-
приводом номинальной мощностью до 2 МШ. В эквивалентной схеме для дальнейших расчётов и моделирования был выбран двигатель типа ЫРВООКБ с номинальной мощностью 2.013 МШ и номинальным током 3053 А. Моменты сопротивления и инерции поезда приведены к его ротору.
Для расчёта энергии потерь при двунаправленном потоке необходимо задаться параметрами схемы СЭ МП, приведенной на рис. 1. Характеристики питающей сети определяются параметрами трёхфазного трансформатора питающей подстанции типа ТМН 4000/35/6 [1], для которого сопротивление фазы Я0 = 0.1 О [1]. Параметры линии 1 (см. рис.1) определяются расстоянием между тяговой подстанцией и преобразовательным трансформатором, которое, в среднем, составляет от 1 до 3 кш [1]. Используемый в линии 1 алюминиевый трехжильный кабель имеет величину сопротивления фазы Ль равное 0.3 О/кш, а его сечение выбирается по току, который может потреблять рассматриваемый электропривод, и равно 95 шш2 [1]. Преобразовательный трансформатор серии ТСЗП-2500/10У3 6(10)/0.71 к¥ имеет номинальную мощность 2.509 МШ и потери короткого замыкания 20 кШ. Суммарное эквивалентное сопротивление его фазы Яту для него будет равно 2 шО. Параметры линии 2 определяются расстоянием между преобразовательным трансформатором т1 и выпрямителем, которое принято равным 50 ш. При этом сечение медного кабеля будет равно 1000 шш2, величина сопротивление фазы Я2 - 0.9 шО. Стальной контактный рельс ЯКя имеет стандартное сечение 6600 шш2 и сопротивление, равное 9 шО/кш. Его длина может изменяться в диапазоне от 1 до 3 кш, в зависимости от расположения поезда на прогоне между станциями. Активное сопротивление выбранной ранее машины постоянного тока (МПТ) RJ равно 8 шО.
КПД СЭ МП с двунаправленным потоком энергии. Оценим КПД системы электроснабжения, представленной на рис. 1. Согласно [3], максимально возможный КПД СЭ с двунаправленным потоком энергии определяется по формуле:
rmax— (2 rmax-^ ) '
- кЕ
1 - кЕ
rmaxo
(1)
0 < кЕ — < 1.
(2)
s—
где Р^ и Р^ - мощность источника в прямом и обратном потоках энергии соответственно.
Величину максимально возможного значения КПД СЭ МП в прямом потоке энергии г/тах^ определяется по выражению [3]:
-1 11 - 1
Vmax— — , + J . "
4 к
(3)
где к5с - коэффициент короткого замыкания, определяемый отношением мощности короткого замыкания на клеммах нагрузки к полезной активной мощности нагрузки:
kSC —
(4)
usf
где Ряс - мощность резистивного короткого замыкания СЭ МП при отключенной нагрузке; - среднее значение полезной активной мощности нагрузки в интервале повторяемости графика движения поезда согласно рис. 2.
Значения Ри^ и кЕ зависят от графика движения поезда, интервалов задания, скоростей разгона и торможения поезда. Мощность резистивного короткого замыкания Р8с зависит от конфигурации СЭ МП и может быть определена из соотношения:
3U 2
sm
PSC
2R
(5)
■Z
где ^тах^ и т]тах^ - максимально возможное значение КПД трехфазной СЭ в прямом и обратном потоках энергии соответственно; кЕ - коэффициент возврата энергии из нагрузки в источник, определяемый по выражению из [3]:
SC
где Usm - амплитуда синусоидального фазного напряжения источника электропитания; RE - эквивалентное активное сопротивление СЭ МП, приведенной на рис. 1.
Величина активного эквивалентного сопротивления системы электроснабжения, согласно рис. 1, включает в себя следующие составляющие:
RZ — R0 + R1 + RTV + R2 + RRF + RKR + RJ , (6) где Ro - сопротивление фазы источника переменного напряжения 6(10) kV, приведенное к вторичной обмотке преобразовательного трансформатора (ПТ);
Rj - сопротивление фазы участка линии 1, приведенное ко вторичной обмотке ПТ; RTV - суммарное сопротивление фазы ПТ; R2 - сопротивление фазы участка линии 2 от трансформатора до выпрямителя; RRF - сопротивление управляемого выпрямителя; RKR - сопротивление контактного рельса; RJ - сопротивление якорной цепи ДПТ.
Величина максимально возможного значения КПД СЭ МП в обратном потоке энергии r¡maxможет быть определена по следующему выражению [3]:
(7)
1 + кЕ kSC
Найдём величину максимально возможного КПД СЭ с двунаправленным потоком энергии и определим возможный диапазон его изменений, используя приведенные выше выражения.
Для определения коэффициента возврата энергии из нагрузки в источник кЕ, согласно (2), необходимо знать PS^_ и PS—. Их значения можно определить из графика движения поезда, приведенного на рис. 2. Вычислив площади под кривой графика изменения мощности для прямого и обратного потоков энергии, получим значения PS— = 50.3 MW, PS^ = 108.7 MW и, в соответствии с (2), кЕ = 0.5.
Среднее значение полезной активной мощности нагрузки найдём интегрированием графика мгновенной мощности в интервале повторяемости движения поезда. Получили значение Puf = 1.44 MW.
Для нахождения мощности резистивного короткого замыкания Р5С определим составляющие эквивалентного активного сопротивления системы электропитания ЯЕ и возможный диапазон их изменений.
Согласно перечисленным выше данным, приведенное сопротивление фазы источника Я0 может быть вычислено по выражению:
Яо = кЯо, (8)
где к = 1/кй2 - коэффициент приведения параметров элементов первичной обмотки преобразовательного трансформатора ко вторичной, равный 0.014. Значение сопротивления Яо равно 1.4 шО.
Аналогично вычисляется приведенное сопротивление фазы линии 1 Я1:
Я1 = кЯ1. (9)
Остальные составляющие выражения (6) и возможный диапазон их изменения были определены выше, значения параметров сведены в табл. 1, в соответствии с которой сопротивление Я1 лежит в диапазоне от 4.2 шО до 12.6 шО. Активное эквивалентное сопротивление СЭ ЯЕ, приведенной на рис. 1, будет иметь значения в диапазоне от 27 шО до 44 шО.
Таблица 1
Сопротивления схемы СЭ МП и диапазон их изменений
2 _i
1 _ kEksckadd ^
Параметр Значение
Ro , mQ 1.4
1000 m 4.2
R1, mQ 2000 m 8.4
3000 m 12.6
Rjy, mQ 3
R2, mQ 1
Rrf, mQ 1
1000 m 9
RKR, mQ 2000 m 13.5
3000 m 18
Rj, mQ 8
1 +
0.5 + -J0.25 _ k
sc
_1
-_ кг
_1
add ^
Vreal —
1 _ кг
-•(10)
График зависимости реального КПД СЭ с двунаправленным потоком энергии от коэффициента короткого замыкания на клеммах нагрузки ksc приведен на рис. 3 пунктирной линией.
Мощность короткого замыкания РС рассчитанная по (5), в зависимости от параметров схемы, имеет значение от 34 до 56 МЖ Коэффициент короткого замыкания, рассчитанный по (4), в зависимости от активного эквивалентного сопротивления, лежит в диапазоне от 25 до 40.
В реальной системе электроснабжения могут присутствовать дополнительные потери электроэнергии, которые можно учесть в теоретических расчётах введением коэффициента дополнительных потерь ка¿¿. В этом случае, КПД СЭ МП может быть рассчитано из [3] по выражению:
Рис. 3. КПД СЭ МП
Моделирование СЭ МП с двунаправленным потоком энергии. Для экспериментальной проверки теоретических результатов была разработана МаЛаЬ-модель СЭ МП с четырехквадрантным управляемым выпрямителем, имитирующая схему, приведенную на рис. 1. Ма^аЬ-модель изображена на рис. 4. Она состоит из следующих блоков:
• силовая схема - блоки 1, 3, 4, 5, 7, 8, 10;
• система управления тиристорными УВ - блок 6;
• регуляторы момента, тока и скорости - блоки 9, 13;
• датчик тока и напряжения - блок 2;
• расчетчик - блок 11;
• многолучевой осциллограф - блок 12.
Назначение блоков силовой схемы: 1 - промышленная сеть; 3 - кабели, соединяющей тяговую подстанцию и преобразовательный трехфазный трансформатор 6(10) к¥ / 0.71 кУ, который обозначен блоком 4;
5 - кабели, идущие от трансформатора 4 до выпрямительных мостов 7; 8 - стальной контактный рельс, соединяющий УВ с двигателем постоянного тока 10.
Параметры элементов силовой цепи в модели были заданы в строгом соответствии с данными СЭ МП, приведенными ранее. Данные модели двигателя постоянного тока соответствую таковым для машины типа ЫРВООКБ. Механическая часть электропривода была приведена к ротору машины постоянного тока, а кинетическая энергия, запасаемая составом при движении - к энергии эквивалентного маховика. Задаваемая в блоке 9 нагрузка МПТ учитывает как потери собственных нужд, так и потери на трение и лобовое сопротивление воздуха движущемуся составу.
Система управления выпрямителем построена по вертикальному принципу и имеет арккосинусоидаль-ную характеристику фазосмещающего устройства. Логика переключения мостов отслеживает сигнал задания с выхода регулятора и мгновенную величину тока нагрузки, принимая решение на перевод импульсов в зависимости от их суперпозиции.
^ГУ-ч
v_puy—W [kpd")-» u
► * —►
V'abc A labt
В a b С с
Three-Ftiase Source
ппп Regulator
Hl
13
Thyryslor bridges 7
Рис. 4. MatLab -модель СЭ МП
Система авторегулирования выполнена по замкнутому принципу с использованием двухконтурного подчинённого регулятора тока-скорости, настроенного на технический оптимум. Это обеспечило качественное поддержание заданной скорости в соответствии с графиком движения поезда.
В результате моделирования получены: осциллограммы изменения тока, напряжения и скорости потока энергии для интервалов движения состава в рассматриваемой СЭ МП, приведенные на рис. 5. Построен график зависимости реального КПД СЭ с двунаправленным потоком энергии от коэффициента короткого замыкания на клеммах нагрузки, который, для наглядности, был приведен на рис. 3, совместно с теоретическим графиком.
1000 то
m да
200
о -100 бооо
401» 2000
0
■2000 -дан
4
2
1 О -!
V
s
MW
2(1 40 60 SO ШО 120 140 МП ¡80 Рис. 5. Осциллограммы работы схемы
200
Выводы.
1. Разработана эквивалентная схема системы электроснабжения метрополитена с четырёхквадрантным
приводом постоянного тока и возможностью рекуперации энергии в питающую сеть в режиме торможения, определены её параметры, задан график движения электропоезда.
2. С использованием известных формул для разработанной системы электроснабжения получена теоретическая зависимость КПД СЭ с двунаправленным потоком энергии от коэффициента резистивного короткого замыкания на клеммах нагрузки.
3. Построена MatLab -модель СЭ МП с четырех-квадрантным электроприводом постоянного тока на базе шестипульсного мостового управляемого выпрямителя с возможностью реализации двунаправленного потока электрической энергии между источником и нагрузкой.
4. С использованием MatLab-модели снята реальная зависимость КПД СЭ с двунаправленным потоком энергии от коэффициента короткого замыкания на клеммах нагрузки, повторяющая динамику теоретической кривой. Некоторое расхождение между теоретической кривой и данными полученной модели может быть объяснено учетом в модели дополнительных потерь электроэнергии, таких, как трение, сопротивления воздуха, коммутационные потери и т. д.
СПИСОК ЛИТЕРАТУРЫ
1. Жемеров Г.Г., Ильина Н.А, Тугай Д.В., Холод О.И. Системы электроснабжения метрополитена с современными полупроводниковыми преобразователями и накопителями энергии // Електротехнжа i електромехашка. - 2013. - №1. -С. 41-49. doi: 10.20998/2074-272X.2013.1.10.
2. Жемеров Г.Г., Тугай Д.В., Холод О.И. Энергоэффективность систем электроснабжения подвижного состава метрополитена // Техшчна електродинамжа. - 2014. - №1. - С. 67-74.
3. Жемеров Г.Г., 1лына Н.О., Машура А.В., Тугай Д.В. Енергоефектившсть рекуперативного гальмування в системах електропостачання з регульованими електроприводами // Техшчна електродинамжа. - 2018. - №6. - С. 73-76. doi: doi: 10.15407/techned2018.06.073.
4. Электровагоны типов 81-718/719, 81-718.0/719.0 и 81718.2/719.2. Режим доступа: http://vagon.metro.ru/passenger/81-718.html.
5. Костин Н.А., Никитенко А.В. Автономность рекуперативного торможения - основа надежной энергоэффектив-
s
s
ной рекуперации на электроподвижном составе постоянного тока // З^зничний транспорт Украши. - 2014. - №3. - С. 15-23.
6. Щербак Я.В., Нерубацький В.П. Ан^з застосування рекуперативного гальмування на затзницях Украши // Заль зничний транспорт Украши. - 2011. - №2. - С. 30-34.
7. Саблин О.И. Анализ качества рекуперируемой электроэнергии в системе электрического транспорта// Вюник НТУ «ХП1». - 2013. - №38. - С. 187-190.
8. Посохов 1.М. Дослщження ринку вагонобудування Украши та конкурентоспроможносп промислових тдприемств залiзничного транспорту на свгговому ринку та ринку краш СНД // Вiсник НТУ «ХП1». - 2015. - №60. - С. 115-118.
REFERENCES
1. Zhemerov G.G., Ilyina N.A., Tugay D.V., Kholod O.I. Subway power systems with modern semiconductor converters and energy storage devices. Electrical engineering & electrome-chanics, 2013, no.1, pp. 41-49. (Rus) doi: 10.20998/2074-272X.2013.1.10.
2. Zhemerov G.G., Tugay D.V., Kholod O.I. Energy efficiency of power supple systems for a subway. Technical electrodynamics, 2014, no.1, pp. 67-74. (Rus).
3. Zhemerov G.G., Ilina N.A., Mashura A.V., Tugay D.V. Efficiency of regenerative braking in energy supply systems with electric regulated drives. Technical electrodynamics, 2018, no.8, pp. 73-76. (Ukr). doi: 10.15407/techned2018.06.073.
4. Electric carriage types 81-718/719, 81-718.0/719.0 and 81718.2/719.2. Available at: http://vagon.metro.ru/passenger/81-718.html (Accessed 12 May 2016). (Rus).
5. Kostin N.A., Nikitenko A.V. Autonomy of recuperative braking - the basis of reliable energy-efficient recovery on the electromotive force constant current. Railway Transport of Ukraine, 2014, no.3, pp. 15-23. (Rus).
6. Shcherbak Ya.V., Nerubatsky V.P. Analysis of the use of recuperative inhibition on the railways of Ukraine. Railway Transport of Ukraine, 2011, no.2, pp. 30-34. (Ukr).
7. Sablin O.I. Analysis of the quality of the recovered electricity in the electric transport system. Bulletin of NTU «KhPI», 2013, no.38, pp. 187-190. (Rus).
8. Posokhov I.M. Research of the car-building market of Ukraine and competitiveness of the industrial enterprises of railway transport in the world market and CIS markets. Bulletin of NTU «KhPI», 2015, no.60, pp. 115-118. (Ukr).
Поступила (received) 05.06.2019
Жемеров Георгий Георгиевич1, д.т.н., проф., Крылов Денис Сергеевич1, к.т.н., доц., Машура Артём Вячеславович1, аспирант, 1 Национальный технический университет «Харьковский политехнический институт», 61002, Харьков, ул. Кирпичева, 2, тел/phone +380 577076312, e-mail: [email protected]
G.G. Zhemerov1, D.S. Krylov1, A. V. Mashura1 1 National Technical University «Kharkiv Polytechnic Institute», 2, Kyrpychova Str., Kharkiv, 61002, Ukraine. Energy efficiency of the subway electrical supply system with electrical energy recovery at braking. Purpose. The purpose of the paper is to assess the efficiency of the subway power supply system, which uses a four-quadrant DC drive with energy recovery in the supply network in the braking mode. Methodology. We have applied the theory of electrical circuits and mathematical simulation in Matlab package. Results. The theoretical dependence of the efficiency of the electrical supply system with a bidirectional flow of energy on the coefficient of resistive short circuit at the load terminals has been obtained. The theoretical result is verified by modeling. Originality. The equivalent circuit of the subway power supply system with a four-quadrant DC drive and the possibility of energy recovery to the supply network in braking mode is developed, its parameters are determined, and the schedule of the electric train movement was set. Practical value. The use of the obtained dependencies and simulation results will allow to determine the direction of the future development of the subway power supply system and optimize its energy efficiency. References 8, tables 1, figures 5. Key words: power supply system, energy, energy return coefficient, efficiency, energy recovery.