CC BY
111
УДК 621.311 DOI: 10.30724/1998-9903-2020-22-3-60-67
МОДЕЛИРОВАНИЕ ГАЗОТУРБИННОЙ УСТАНОВКИ С ПРОГНОСТИЧЕСКИМИ РЕГУЛЯТОРАМИ НАПРЯЖЕНИЯ И СКОРОСТИ
Ю. Н. Булатов1, А. В. Крюков2' 3, Нгуен Ван Хуан 3
1Братский государственный университет, г. Братск, Россия 2 Иркутский государственный университет путей сообщения, г. Иркутск, Россия 3 Иркутский национальный исследовательский технический университет,
г. Иркутск, Россия
ORCID: https://orcid.org/0000-0002-3716-5357, [email protected]
Резюме: В статье представлены результаты компьютерных исследований, цель которых состояла в разработке динамических моделей двухвальных газотурбинных установок (ГТУ) малой мощности, оснащенных автоматическими регуляторами возбуждения и скорости с прогностическими звеньями. При использовании ГТУ в изолированных системах электроснабжения отклонения частоты и напряжения могут быть значительными, что неблагоприятно влияет на показатели качества электроэнергии. Настройка регуляторов вводимой в эксплуатацию ГТУ требует сложных расчетов. Применение алгоритмов, позволяющих с минимальными затратами усовершенствовать классические регуляторы, позволяет эффективно решать задачи их настройки. К их числу можно отнести прогностические алгоритмы, позволяющие настроить систему автоматического управления с помощью одного параметра - времени прогноза. Это особенно актуально в условиях необходимости ускоренного ввода в эксплуатацию ГТУ. Прекдставленные в статье модели реализованы с использованием пакета Simulink системы MATLAB. Исследовано влияние прогностических алгоритмов на качество управления напряжением и скоростью вращения ротора генератора ГТУ в режимах подключения и сброса нагрузки. Результаты компьютерного моделирования дали возможность сформулировать следующие выводы: применение прогностических регуляторов при подключении к ГТУ дополнительной нагрузки позволяет убрать колебательность, уменьшить провалы напряжения, снизить время переходного процесса на 2,5 с; при резком сбросе нагрузки за счет применения прогностических алгоритмов удается полностью убрать колебательность, уменьшить перерегулирование скорости вращения ротора и перенапряжения на зажимах генератора, а также заметно снизить время переходного процесса для скорости по сравнению с классическими регуляторами; использование прогностическиих алгоритмов позволяет получить приемлемые показатели качества переходных процессов без применения сложных процедур настройки регуляторов.
Ключевые слова: распределенная генерация, газотурбинные установки, автоматические регуляторы возбуждения и скорости, прогностические алгоритмы.
Для цитирования: Булатов Ю.Н., Крюков А.В., Нгуен Ван Хуан. Моделирование газотурбинной установки с прогностическими регуляторами напряжения и скорости // Известия высших учебных заведений. ПРОБЛЕМЫ ЭНЕРГЕТИКИ. 2020. Т. 22. № 3. С. 60-67. doi:10.30724/1998-9903-2020-22-3-60-67.
MODELING A GAS-TURBINE UNIT WITH PROGNOSTIC REGULATORS OF VOLTAGE AND SPEED
YuN. Bulatov1, AV. Kryukov 2'3, Nguyen Van Huan 3
1 Bratsk State University, Bratsk, Russia 2 Irkutsk State Transport University, Irkutsk, Russia 3 Irkutsk National Research Technical University, Irkutsk, Russia
ORCID: https://orcid.org/0000-0002-3 716-5357, [email protected]
Abstract: Object. The article presents the results of computer research, the purpose of which was to develop dynamic models of low-power twin-shaft gas turbine units (GTU) equipped with automatic excitation and speed regulators with predictive links. Methods. Setting up the regulators of the commissioned gas turbine requires complex calculations. The use of algorithms that make it possible to improve the classical regulators with minimal costs allows us to effectively solve the problems of their settings. These include prognostic algorithms that allow you to configure the automatic control system using one parameter - forecast time. The models presented in the article are implemented using the Simulink package of the MATLAB system. Results. The influence of prognostic algorithms on the quality of controlling the voltage and rotor speed of a gas turbine generator in the modes of connecting and dumping loads is studied. Studies have shown that increasing the coefficient of the amplifier of the auto-predictive speed regulator can significantly reduce overshoot and the transition process, and also has a positive effect on voltage regulation. Conclusions. The results made it possible to formulate the following conclusions: the use of prognostic regulators when connecting an additional load to the gas turbine allows you to remove the oscillation, reduce voltage dips, reduce the transient time by 2.5 s; with a sharp load shedding due to the use of prognostic algorithms, it is possible to completely remove the oscillation, reduce the overshoot of the rotor speed and overvoltage at the generator terminals, and also significantly reduce the transient time for speed compared to classical regulators; the use ofprognostic algorithms allows us to obtain acceptable quality indicators of transients without the use of complex regulator tuning procedures.
Key words: distributed generation, gas turbine units, automatic regulators of excitation and speed, prognostic algorithms.
For citation: Bulatov YuN, Kryukov AV, Nguyen Van Huan. Modeling a gas-turbine unit with prognostic regulators of voltage and speed. Power engineering: research, equipment, technology. 2020;22(3):60-67. doi: 10.30724/1998-9903-2020-22-3-60-67.
Использование установок распределённой генерации (РГ), находящихся в непосредственной близости от потребителей электроэнергии, позволяет уменьшить потери мощности и повысить надежность электроснабжения. Это особенно актуально при модернизации изолированных систем электроснабжения (СЭС).
Широкое практическое применение для электроснабжения потребителей находят многовальные газотурбинные установки (ГТУ), имеющие лучшие характеристики, чем у одновальных [1-3]:
- при изменениях нагрузки в широких пределахах экономичность многовальных ГТУ
значительно выше;
- момент механической инерции многовальных ГТУ существенно меньше;
- электромеханическая постоянная времени многовальных ГТУ в 2.. .4 раза ниже.
В работе рассматривается двухвальная ГТУ, структурная схема которой показана на рис. 1. Путем изменения подачи топлива в камеру сгорания с помощью автоматического регулятора скорости (АРС) поддерживается заданная частота вращения силовой турбины и ротора генератора; при этом мощность силовой турбины изменяется с некоторым запаздыванием, обусловленным моментом инерции вращающихся масс.
Рис. 1. Структурная схема двухвальной ГТУ: К - компрессор; Т - турбина высокого давления; КС - камера сгорания; СТ - силовая турбина; Р - редуктор; ДЧ - датчик частоты вращения; ОВ - обмотка возбуждения; ТН - трансформатор напряжения; АРВ - автоматический регулятор возбуждения; АРС - автоматический регулятор
Введение
Синхронный генератор
газ
скорости 61
При использовании ГТУ в изолированных СЭС отклонения частоты и напряжения от номинального значения могут быть значительными, что неблагоприятно влияет на показатели качества электроэнергии. Настройка регуляторов вновь вводимой в эксплуатацию ГТУ традиционными методами требует многовариантных и сложных расчетов. Применение алгоритмов, позволяющих с минимальными затратами усовершенствовать классические регуляторы позволяет более эффективно решать задачи по их настройке. К их числу можно отнести прогностические алгоритмы [4, 5], позволяющие настроить систему автоматического управления с помощью одного параметра - времени прогноза. Это особенно актуально в условиях необходимости ускоренного ввода в эксплуатацию ГТУ.
В статье приведены результаты компьютерных исследований влияния прогностических алгоритмов на качество управления напряжением и скоростью вращения ротора генератора ГТУ в режимах подключения и сброса нагрузки потребителей. Моделирование ГТУ с автоматическими регуляторами проводилось в системе MATLAB.
Применение прогностического АРВ и автопрогностического АРС позволяет значительно уменьшить перерегулирование и время переходного процесса напряжения и частоты вращения ротора синхронного генератора по сравнению с классическими регуляторами без использования сложных процедур настройки.
Описание модели ГТУ и используемых прогностических регуляторов
При исследовании динамических режимов работы ГТУ возникает задача выбора модели турбины; при этом могут использоваться модели ГТУ различной степени детализации [6-9]. В основном это модели одновальных ГТУ. Для применения целого ряда моделей ГТУ требуются экспериментальные характеристики [10-12]. Сложность использования этих моделей состоит в том, что необходим большой объем данных о конкретной ГТУ, которых нет в открытом доступе. Проблемным также являются трудоемкость построения полной модели ГТУ и сложность подбора ее параметров для обеспечения стабильной работы во всех изучаемых режимах [10].
Динамические свойства двухвальной ГТУ достаточно точно можно учесть с помощью модели [1], схема которой представлена на рис. 2. Модель реализована с использованием пакета Simulink системы MATLAB. Моделировалась ГТУ номинальной мощностью 2,5 МВт на базе двухвальной газовой турбины с редуктором со стороны генератора. Номинальная частота вращения силовой турбины составляет 5500 об/мин, а ротора синхронного генератора - 1000 об/мин. В представленной модели учитываются расход топлива и изменение кинетической энергии газогенератора (турбины и компрессора).
Камера сгорания (блок Combustion chamber на рис. 2), обуславливающая некоторое запаздывание при изменении подачи и скорости сгорания топлива, моделировалась апериодическим звеном первого порядка с передаточной функцией
1
T„S +1
где Tcc - постоянная времени камеры сгорания, которая принималась равной 0,01 с.
Газогенератор моделировался интегрирующим звеном, охваченным жесткой отрицательной обратнойй связью [1], постоянная времени которого принималась равной 1,99 с исходя из мощности и частоты вращения турбины с моментом инерции 454 кг-м2 , приведенного к частоте вращения ротора генератора. Механическая мощность турбины определялась с учетом доли мощности турбины газогенератора по отношению к номинальной мощности силовой турбины с помощью коэффициента Кг (при моделировании принимался равным 0,55). Передаточный коэффициент редуктора ^ составлял 5,56. Инерционность редуктора учитывалась в модели генератора.
Turbine
Рис. 2. Схема модели двухвальной газовой турбины в Simulink 62
Синхронный генератор моделировался стандартным блоком пакета SymPowerSystems - Synchronous Machine pu Fundamental. При моделировании использовались следующие параметры: индуктивное сопротивление по продольной оси Xd = 2,34 о.е.; ЭДС E = 1,25 о.е.; напряжение U = 1 о.е.; постоянная механической инерции Т = 1,77 с и др.
Эквивалентная постоянная механической инерции ГТУ (турбины и генератора) составляла
Tje = 3,76 с.
Схема используемой модели ГТУ, работающей на выделенную нагрузку, представлена на рис. 3. В модели предусматривалось внесение возмущения в виде отключения или подключения активно-индуктивной нагрузки мощностью 0,645 + j0,17 МВА с помощью выключателя (блок Breaker). Подробное описание используемых в модели регуляторов скорости и напряжения приведено в [5,13,14].
На рис. 4 и 5 представлены схемы используемых Simulink-моделей прогностического АРВ и автопрогностического АРС [5,13]. В блоке Amplifier (рис.4) учитывалось
K
последовательное соединение электронного усилителя с передаточной функцией -a— и
Tas +1
линейного прогнозирующего звена с передаточной функцией Tps +1. При моделировании
для АРВ принимались следующие числовые значения параметров: Ka = 1; Ta = 0,001 с, Tp = 0,525 с.
Automatic prognostic speed regulator
delta
Gas turbine
[—-0
ARS 0_
Ug
If Uare
Uare Vf
Prognostic voltage Excitation System
Ug delta
Discrete Ts=1e-05 s.
43
43
Synchronous Machine pu Fundamental
A Vabc
B a b
C c
A B-VWT- C
Bff
C
0.645+j0.17 MVA
Рис.3. Схема модели ГТУ в MATLAB
Р
Рис. 4. Схема Simulink-модели прогностического АРВ: Ug - напряжение генератора; SetUg - заданное значение напряжения генератора; If- ток возбуждения генератора; <т - скорость вращения ротора генератора; Set &т - заданное значение скорости; k0u, k1u, kUf, k0<, k1ra - коэффициенты настройки АРВ
Tp
Рис. 5. Схема Simulink-модели автопрогностического АРС:
1.8+j0.4 MVA
Kp, К,, Kd - коэффициенты настройки АРС
Постоянная времени прогнозирующего звена АРС определлась частотой собственных колебаний ГТУ по выражению [5]
j, _ 1 _ 2п __2п__(1)
LP
fp \dPg
.ж E • Ug .cosb.ж
dbg Tje y Xd Tje
где fp - частота собственных колебаний ротора генератора при неизменности по модулю и
углу вектора напряжения на его шинах - парциальная частота, являющаяся параметром генератора, комплексно характеризующим его инерционность и жесткость связи с СЭС [15]; œg - номинальная частота вращения ротора генератора, рад/с; Тje - эквивалентная
постоянная механической инерции агрегата, с; X d - индуктивное сопротивление
генератора по продольной оси, о.е.; Ug - номинальное напряжение генератора, о.е.;
5 - угол между напряжением и ЭДС Eq генератора, зависящий от нагрузочного режима, эл.
град.
Частота собственных колебаний моделируемой ГТУ для указанных выше параметров в зависимости от угла нагрузки 5 представлена в модели автопрогностического АРС (рис.5) следующей функцией:
fp(5) = 0,614-Vcos 5 .
Результаты моделирования
Исследования проводились для выявления влияния прогностических алгоритмов на качество управления напряжением и скоростью вращения ротора генератора ГТУ в следующих режимах:
- подключение дополнительной нагрузки (28 % от номинальной мощности установки);
- резкий сброс нагрузки (30 % от номинальной мощности установки).
При моделировании коэффициенты настройки классических и прогностических регуляторов задавались, исходя из практических соображений, следующими: для АРВ kou = 5, kiu = 1, kf 0, koffl = 1, = 1; для АРС Kp = 1, K = 1, Kd = 0,1; при этом коэффициент усилителя Ka для классического и автопрогностического АРС принимал значения 1 и 50.
Результаты моделирования указанных режимов работы ГТУ с классическими и прогностическими АРВ и АРС представлены на рис. 6-8.
Полученные результаты показали, что увеличение коэффициента усилителя Ka АРС позволяет значительно уменьшить перерегулирование и время переходного процесса для скорости вращения ротора, а также оказывает положительное влияние на регулирование напряжения (рис. 6 и 7).
Напряжение, о.е.
1.02 ]
0.96 0.94 А. 0.9 0.33 0.36
2
1
Время, с
Ю 11 13 14 13 16 17 ]] Г_ 13 14 is 16 ]7
а) б)
Рис. 6. Осциллограммы скорости вращения ротора (а) и напряжения (б) генератора при подключении дополнительной нагрузки (для АРС коэффициент усилителя Ka=1): 1 - АРВ и АРС без прогнозирующих звеньев; 2 - прогностический АРВ и автопрогностический АРС
а) б)
Рис. 7. Осциллограммы скорости вращения ротора (а) и напряжения (б) генератора при подключении дополнительной нагрузки (для АРС коэффициент усилителя Ka=50): 1 - АРВ и АРС без прогнозирующих звеньев; 2 - прогностический АРВ и автопрогностический АРС
Напряжение, о.е.
V к
V /
10 11 12 13 14 15 16 17 10 11 12 В 14 1S 16 17
а) б)
Рис.8. Осциллограммы скорости вращения ротора (а) и напряжения (б) генератора при сбросе нагрузки на 30% (для АРС коэффициент усилителя Ka=1): 1 - АРВ и АРС без прогнозирующих звеньев; 2 - прогностический АРВ и автопрогностический АРС
Применение прогностического АРВ и автопрогностического АРС дает возможность:
- при подключении к ГТУ дополнительной нагрузки убрать колебательность (рис. 6 и 7), уменьшить провалы напряжения и скорости вращения ротора генератора (рис. 6, б); снизить время переходного процесса для скорости вращения ротора генератора на 2,5 с и получить время переходного процесса для напряжения и скорости практически одинаковыми (рис. 6, а и б);
- при резком сбросе нагрузки ГТУ полностью убрать колебательность (рис. 8), уменьшить перерегулирование скорости вращения ротора и перенапряжение на зажимах генератора, а также снизить на 2,5 с время переходного процесса для скорости по сравнению с классическими регуляторами.
Заключение
Разработаны компьютерные модели двухвальной ГТУ малой мощности с прогностическим АРВ и автопрогностиическим АРС. Исследовано влияние прогностических алгоритмов на качество управления напряжением и скоростью вращения ротора генератора ГТУ в режимах подключения и сброса нагрузки потребителей. Результаты компьютерного моделирования позволяют сделать следующие выводы:
1. Применение прогностического АРВ и автопрогностического АРС при подключении к ГТУ дополнительной нагрузки позволяет убрать колебательность, уменьшить провал напряжения и скорости вращения ротора генератора, снизить время переходного процесса на 2,5 с;
2. При резком сбросе нагрузки ГТУ прогностические алгоритмы дают возможность полностью убрать колебательность, уменьшить перерегулирование скорости вращения ротора и перенапряжения на зажимах генератора, а также снизить время переходного процесса по сравнению с классическими регуляторами.
3. Использование прогностическиих алгоритмов позволяет получить достаточно высокие показатели качества переходного процесса напряжения и скорости вращения ротора генератора ГТУ без применения сложных расчетов и процедур настройки регуляторов.
Литература
1. Меркурьев Г.В., Шаргин Ю.М. Устойчивость энергосистем. В 2-х т. СПб.: НОУ «Центр подготовки кадров энергетики», 2008. Т. 2. 376 с.
2. Асаинов Д.Н., Гусев Ю.П. Исследование динамической устойчивости и электродинамической стойкости ГТУ // Вестник МЭИ. 2010. №2. C. 55-61.
3. Илюшин П.В., Филиппов С.П., Новиков Н.Л. Требования к маневренности газотурбинных и газопоршневых генерирующих установок // В сборнике: Методические вопросы исследования надежности больших систем энергетики в 2-х книгах. 2019. С. 343-352.
4. Пикина Г.А. Принцип управления по прогнозу и возможность настройки систем регулирования одним параметром // Новое в российской электроэнергетике. 2014. № 3. С. 5-13.
5. Bulatov Yu.N., Kryukov A.V., Nguyen Van Huan. Automatic prognostic regulators of distributed generators // International Multi-Conference on Industrial Engineering and Modern Technologies (Vladivostok, 03-04 October 2018). Vladivostok: IEEE, 2018. V. 463. pp. 1-4. https://doi.org/10.1109/FarEastCon.2018.8602718
6. Centeno P., Egido I., Domingo C., et al. Review of gas turbine models for power system stability studies. - 9th Spanish Portuguese Congress on Electrical Engineering, 2005.
7. Rowen W.I. Simplified Mathematical Representations of singleshaft gas turbines in mechanical drive services. - Turbo Mach. Int., 1992. V. 33, N. 5. pp. 26-32.
8. Опарин Д.А., Кавалеров Б.В. О моделировании газотурбинных установок при управлении электростанциями малой и средней мощности // Вестник ПНИПУ. 2014. №12. С.5-13.
9. Mohammad Reza Bank Tavakoli, Behrooz Vahidi, Wolfgang Gawlik. An Educational Guide to Extract the Parameters of Heavy Duty Gas Turbines Model in Dynamic Studies Based on Operational Data // IEEE Trans. Power Syst. 2009. V. 24. N. 3. pp. 1366-1374.
10. Бахмисов О.В., Кузнецов О.Н. Методика моделирования газотурбинных и парогазовых установок большой мощности при исследовании процессов в ЭЭС // Электричество. 2016. № 5. C.27-34.
11. Meegahapola L. and Flynn D. Characterization of Gas Turbine Lean Blowout During Frequency Excursions in Power Networks. IEEE Trans. Power Syst. 2014. N. 99. pp. 1-11.
12. Shalan H., Hassan M., Bahgat A. Comparative Study on Modeling of Gas Turbines in Combined Cycle Power Plants. Proc. of the 14th International Middle East Power Systems Conference (MEPC0N'10). 2010. pp. 970-976.
13. Bulatov Yu.N., Kryukov A.V. Optimization of automatic regulator settings of the distributed generation plants on the basis of genetic algorithm // 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM). IEEE Conference Publications. 2016, pp. 1-6.
14. Булатов Ю.Н., Крюков А.В., Нгуен Ван Хуан. Методика настройки прогностических регуляторов установок распределенной генерации // Известия высших учебных заведений. ПРОБЛЕМЫ ЭНЕРГЕТИКИ. 2016. № 11-12. С. 84-95.
15. Бушуев В.В, Лизалек Н.Н., Новиков Н.Л. Динамические свойства энергосистем. М.: Энергоатомиздат, 1995. 320 с.
Авторы публикации
Булатов Юрий Николаевич - канд. техн. наук, доцент, заведующий кафедрой электроэнергетики и электротехники, Братский государственный университет.
Крюков Андрей Васильевич - д-р техн. наук, профессор кафедры электроэнергетики транспорта, Иркутский государственный университет путей сообщения; профессор кафедры электроснабжения и электротехники, Иркутский национальный исследовательский технический университет.
Нгуен Ван Хуан - магистр техники и технологии, аспирант, Иркутский национальный исследовательский технический университет.
References
1. Merkuryev GV, Shargin YuM. Ustojchivost' energosistem. In 2 volumes. St. Petersburg: NOU «Centr podgotovki kadrov energetiki», 2008. V. 2. 376 p.
2. Asainov DN, Gusev YuP. Research of dynamic stability and electrodynamic resistence of the
gas turbine plants. Bulletin of the Moscow Power Engineering Institute. 2010; 2: 55-61.
3. Ilyushin PV, Filippov SP, Novikov NL. Trebovaniya k manevrennosti gazoturbinnyh i gazoporshnevyh generiruyushchih ustanovok. Metodicheskie voprosy issledovaniya nadezhnosti bol'shih sistem energetiki. 2019;1:343-352.
4. Pikina GA. Princip upravleniya po prognozu i vozmozhnost' nastrojki sistem regulirovaniya odnim parametrom. New in Russian electric power industry. 2014;3:5-13.
5. Bulatov YuN, Kryukov AV, Nguyen Van Huan. Automatic prognostic regulators of distributed generators. International Multi-Conference on Industrial Engineering and Modern Technologies (October 03-04, 2018). Vladivostok: IEEE, 2018;463:1-4. https://doi.org/10.1109/FarEastCon.2018.8602718
6. Centeno P, Egido I, Domingo C, .et al. Review of gas turbine models for power system stability studies. 9th Spanish Portuguese Congress on Electrical Engineering, 2005.
7. Rowen WI. Simplified Mathematical Representations of singleshaft gas turbines in mechanical drive services. Turbo Mach. Int., 1992;33(5):26-32.
8. Oparin DA, Kavalerov BV. O modelirovanii gazoturbinnyh ustanovok pri upravlenii elektrostanciyami maloj i srednej moshchnosti. Bulletin of Perm National Research Polytechnic University. 2014; 12: 5-13.
9. Mohammad Reza Bank Tavakoli, Behrooz Vahidi, Wolfgang Gawlik. An Educational Guide to Extract the Parameters of Heavy Duty Gas Turbines Model in Dynamic Studies Based on Operational Data. IEEE Trans. Power Syst. 2009;24 (3):1366-1374.
10. Bakhmisov OV, Kuznetsov ON. Methods of Simulation of Heavy Duty Gas Turbine Units and Combined Cycle Power Plants for Power Systems Studies. Electricity. 2016;5: 27-34.
11. Meegahapola L. and Flynn D. Characterization of Gas Turbine Lean Blowout During Frequency Excursions in Power Networks. IEEE Trans. Power Syst. 2014; 99: 1-11.
12. Shalan H., Hassan M., Bahgat A. Comparative Study on Modeling of Gas Turbines in Combined Cycle Power Plants. Proc. of the 14th International Middle East Power Systems Conference (MEPCON'IO). 2010, pp. 970-976.
13. Bulatov Yu.N, Kryukov AV. Optimization of automatic regulator settings of the distributed generation plants on the basis of genetic algorithm. 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM). IEEE Conference Publications. 2016. pp. 1-6. doi: 10.1109 / ICIEAM.2016.7911456.
14. Bulatov YuN, Kryukov AV, Huan NV. Methodology of determining forecasting controllers of distributed generation plants. Proceedings of the higher educational institutions. ENERGY SECTOR PROBLEMS. 2016;11-12:84-95.
15. Bushuev VV, Lizalek NN, Novikov NL. Dinamicheskie svojstva energosistem. Moscow: Energoatomizdat. 1995. 320 p.
Authors of the publication
Yuri N. Bulatov - Bratsk State University, Bratsk, Russia. Email: [email protected].
Andrey V. Kryukov - Irkutsk State Transport University, Irkutsk, Russia. Irkutsk National Research Technical University, Irkutsk, Russia.
Van Huan Nguyen - Irkutsk National Research Technical University, Irkutsk, Russia.
Поступила в редакцию 23 апреля 2020 г.
