CC BY
43
© А.В. Мокеев, С.А. Пискунов, Д.Н. Ульянов, Е.И. Хромцов УДК 621.311
ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ И НАДЕЖНОСТИ РЗА ЦИФРОВЫХ ПОДСТАНЦИЙ И ЦИФРОВЫХ РЭС
А.В. Мокеев1, Пискунов2, Д.Н. Ульянов2, Е.И. Хромцов1
Неверный (Арктический) федеральный университет, г. Архангельск, Россия 2ООО «Инженерный центр «Энергосервис», г. Архангельск, Россия
Резюме: ЦЕЛЬ. Рассматриваются вопросы интеграции технологий синхронизированных векторных измерений (СВИ) и цифровой подстанции (ЦПС) для повышения эффективности и надежности систем релейной защиты и автоматики в составе цифровых подстанций и цифровых районов распределительных сетей. МЕТОДЫ. Предложены трехуровневая система релейной защиты и автоматики (РЗА) понижающих подстанций с использованием встроенной в интеллектуальный выключатель защиты (первый уровень), локального устройства РЗА (второй уровень) и централизованного устройства РЗА (третий уровень). Рассматриваемая система РЗА подстанции основана на использовании совмещенной шины процесса и шины подстанции с преимущественным применением синхронизированных векторов (синхровекторов) тока и напряжения вместо выборочных значений токов и напряжений, а также замещением значительной части MMS-сообщений. РЕЗУЛЬТАТЫ. Использование обмена синхровекторами и GOOSE-сообщениями между защитами на подстанции позволит реализовать распределенную обработку данных между защитами отдельных присоединений и реализовать защиты с абсолютной селективностью. Приведены примеры реализации защит с абсолютной селективностью для понизительных подстанций (ППС) 35-220 кВ и распределительных пунктов (РП) 620 кВ. ЗАКЛЮЧЕНИЕ. Произведено моделирование предложенных алгоритмов /*быстродействующих устройств РЗА с абсолютной селективностью с использованием среды динамического моделирования Simulink.
Ключевые слова: надежность, синхровектор, цифровая подстанция, релейная защита, цифровые фильтры.
IMPROVING THE RPA EFFICIENCY AND RELIABILITY OF DIGITAL STEP-DOWN SUBSTATIONS AND DIGITAL GRIDS
AV. Mokeev1, SA. Piskunov2, DN. Ulyanov2, EI. Khromtsov1
1 Northern (Arctic) Federal University, Arkhangelsk, Russia
Engineering Centre «Energoservice», Arkhangelsk, Russia
Abstract: THE PURPOSE. The article deals with the integration of synchronized phasor measurements (SPM) and digital substation (DS) technologies to improve the relay protection and automation systems efficiency and reliability as part of digital substations and digital districts of distribution networks. METHODS. A three-level system of relay protection and automation (RPA) of step-down substations with the use of protection built into the smart circuit breaker (first level), a local RPA device (second level) and a centralized RPA device (third level) is proposed. The considered substation RPA system is based on the use of a combined process bus and a substation bus with the predominant use of synchronized current and voltage phasors (synchrophasors) instead of currents and voltages sample values and replacing a significant part of MMS messages. RESULTS. The use of synchrophasor and GOOSE messaging between substation protections will allow distributed data processing between the protections of individual feeders and implement inherently selective protections. Examples of the inherently selective protection implementation for step-down substations 35220 kV and distribution points (DP) 6-20 kV are given. CONCLUSION. The proposed algorithms of high-speed RPA devices with inherently selective protection were simulated using the Simulink modeling.
Keywords: reliability, synchrophasor, digital substation, relay protection, digital filters.
Введение
Одна из важнейших тенденций автоматизации подстанций и сетей связана с интеграцией технологий цифровой подстанции (ЦПС) и синхронизированных векторных измерений [1]. Основные преимущества использования СВИ связаны со снижением требований к объему передаваемой информации, так как требуется передача всего шести синхронизированных векторов тока и напряжения для расчета на любом уровне управления и защиты десятков параметров режима по конкретному присоединению. На базе синхровекторов отдельных присоединений может быть произведена оценка параметров схемы замещения энергосистемы. Это открывает новые возможности для построения распределенных систем защиты и автоматики (WAPS) подстанций и электрических сетей.
Архитектура систем РЗА понижающих ПС и РП
Интеграция технологий ЦПС и СВИ позволит повысить основные характеристики РЗА, прежде всего связанные с быстродействием, чувствительностью и надежностью функционирования.
Следует выделить следующие основные предпосылки интеграции указанных технологий:
- постоянно возрастающие возможности современной микропроцессорной техники, позволяющей реализовать более эффективные алгоритмы РЗА;
- применение современных датчиков тока и напряжения, в которых отсутствуют или минимизированы погрешности вследствие насыщения, остаточной намагниченности, феррорезонанса;
- наличие горизонтальных связей между устройствами РЗА позволяет организовать обмен быстрыми дискретными и аналоговыми сообщениями и значительно упрощает реализацию защит с абсолютной селективностью;
- идентификация параметров схемы замещения и параметров эквивалентной энергосистемы на основе данных СВИ позволяет значительно повысить чувствительность токовых защит;
- применение многофункциональных устройств с поддержкой СВИ на различных уровнях управления подстанции позволит реализовать распределенные системы РЗА с эффективным резервированием даже в случае частичной или полной деградации коммуникационной сети на подстанции.
В настоящее время большинство микропроцессорных устройств РЗА среднего напряжения включают чаще всего традиционный набор защит, прежде всего токовых ступенчатых защит. Гораздо реже используются дистанционные и дифференциальные защиты. Использование СВИ позволит более широко применять защиты с абсолютной селективностью на основе токового или дистанционного принципов. Использование дополнительных признаков (критериев) распознавания повреждения [2-4] на основе синхровекторов позволит повысить быстродействие и чувствительность устройств РЗА.
Несмотря на наличие программируемой логики и встроенных часов реального времени, логические операции в РЗА выполняются по-прежнему без привязки ко времени. В этой связи для совершенствования РЗА различного функционального назначение перспективно использование темпоральной логики (temporal logic), или по-другому, временной логики, которая используется для описания последовательностей явлений для их взаимосвязи по временной шкале [5]. Имеется опыт использования темпоральной логики в устройствах регулирования напряжения в энергосистемах [6]. Использование временной логики в устройствах защиты и автоматики, при реализации оперативных блокировок позволит повысить эффективность устройств РЗА.
На основе 6 синхровекторов тока и напряжения может быть рассчитано на любом уровне управления и защиты более 50 параметров режима присоединения, а при использовании дополнительно еще 6 эквивалентных синхровекторов - около 200 параметров режима [7]. При использовании эквивалентных синхровекторов становится возможным расчет параметров по симметричным составляющим с учетом влияния высших гармоник. На основе данных СВИ появляется возможность производить оценку параметров схемы замещения электрической сети, в том числе параметров эквивалентной энергосистемы. Кроме того, появляется возможность распознавания режимов, в том числе связанных с процессами пуска и самозапуска электродвигателей.
Интеграция технологий ЦПС и СВИ открывает новые возможности для построения распределенных многоуровневых систем защиты и автоматики подстанций и электрических сетей. Создание распределенных систем РЗА позволит более широко использовать защиты с абсолютной селективностью с реализацией известных или новых принципов защиты. Актуальными вопросами для создания указанных систем являются реализация в преобразователях аналоговых сигналов (ПАС), интеллектуальной коммутационной аппаратуре и локальных устройствах РЗА функций измерения синхровекторов тока и напряжения с требуемой точностью и быстродействием. При этом отпадает необходимость вычислять синхровекторы в каждом устройстве-подписчике, что резко понижает требуемый уровень вычислительной мощности централизованных устройств РЗА.
Рассматриваемая в статье распределенная система РЗА подстанции основана на использовании совмещенной шины процесса и шины подстанции с преимущественным применением синхровекторов ^Р) вместо выборочных значений (^У) токов и напряжений, а также замещением значительной части М/^-сообщений. Это позволит снизить трафик передачи данных и обеспечит повышение темпа передачи синхровекторов. Другой важной задачей является повышение быстродействия устройств РЗА с поддержкой СВИ [8].
Использование обмена синхровекторами и GOOSE-сообщеншми между защитами на подстанции позволит реализовать распределенную обработку данных между защитами отдельных присоединений и реализовать защиты с абсолютной селективностью. Например, в устройстве релейной защиты вводов за счет использования информации о синхровекторах токов присоединений среднего напряжения и синхровекторов напряжений на шинах подстанции возможна реализация защиты шин с абсолютной селективностью. Использование наряду с дифференциальным током дополнительной информации о синхровекторах тока присоединений позволит повысить быстродействие и чувствительность защиты шин [8].
На рис.1 приведена архитектура распределенной системы РЗА для ППС (РП). Основу автоматизации распределительного устройства среднего напряжения составляют цифровые КРУ с использованием цифровых комбинированных датчиков тока и напряжения (ЦКДТН), интеллектуальных выключателей (ИВ) со встроенной резервной защитой, цифровых датчиков положения (ЦДП), температуры (ЦДТ), датчиков для дуговой защиты (ЦДО) и т.д. Все перечисленные компоненты внутри КРУ, а также устройство РЗА, соединяются между собой посредством резервируемой низкоуровневой детерминированной шины процесса.
ЦКДТН ЦКДТН ЦКДТН
цдт|цдп| Цдо
устройство РЗА трансформатора (питающией линии)
цифровое КРУ (отходящая линия)
АСУ ТП
цифровой комбинированный цифровые датчики: датчик тока и напряжения температуры, оптический, положения и т.д.
Ш
цифровое КРУ
II ш .
ШИНА ПОДСТАНЦИИ
ИЭУ ■ РЗА
ЦКДТН ЦКДТН ЦКДТН
цдфдп| цдо
Рис. 1. Структурная схема трехуровневой системы РЗА ППС (РП)
ЦКДТН является не только источником информации для ИЭУ различного функционального назначения, но и способен выполнять функции ряда ИЭУ за счёт измерения синхровекторов тока и напряжения, использование которых позволит существенно снизить требования к микроконтроллерам или сигнальным процессорам как локальных устройств РЗА среднего напряжения, так и централизованного устройства РЗА (ЦРЗА). Функции ЦРЗА целесообразно возложить на устройство ввода (рис.1) с реализацией функций защиты шин с абсолютной селективностью на основе информации
синхровекторов тока со всех присоединений распределительного устройства среднего напряжения, а также функций резервирования защит отходящих присоединений [9].
На базе решения ЦПС согласно рис.1 становится возможной реализация как трехуровневой распределенной системы защиты и автоматики (интеллектуальный выключатель, локальное устройство РЗА, централизованная система РЗА), так и двухуровневой системы защиты (возможные варианты: интеллектуальный выключатель и централизованная система РЗА, интеллектуальный выключатель и локальное устройство РЗА).
Реализация быстродействующих систем РЗА
При использовании синхровекторов тока и напряжения отдельного присоединения могут легко быть вычислены параметры режима сети, необходимые для реализации токовых и дистанционных защит, а при использовании дополнительной информации о синхровекторах с других присоединений - дифференциальная защита трансформатора, линии, шин.
Применение токовых ступенчатых защит часто ограничено вследствие недостаточной чувствительности отдельных ступеней защиты. Чувствительность таких защит может быть существенно повышена за счет оценки параметров схемы замещения подстанции и параметров эквивалентной энергосистемы [8]. Токовые защиты при этом будут близки по чувствительности к дистанционным защитам. Но это тема отдельного исследования.
Использование СВИ и быстрых коммуникаций в рамках ЦПС позволит шире использовать дистанционный и дифференциальный принципы выполнения защиты с использованием дополнительных признаков распознавания аварий и реализацией защиты с абсолютной селективностью на основе информационного обмена между защитами [8]. Массовое производство недорогих микроконтроллеров с поддержкой операций с плавающей запятой при использовании СВИ делает еще более перспективным использование дистанционного принципа выполнения защиты. При этом при использовании дистанционного принципа за счет информационного обмена также можно реализовать защиту с абсолютной селективностью. При полной деградации коммуникационной сети становится невозможной работа защиты с абсолютной селективностью. Поэтому осуществляется переход устройств защиты в автономный режим работы с относительной селективностью. Это позволяет снизить требования к надежности коммуникаций на подстанции и между подстанциями.
Рассмотрим вопросы совершенствования защит линии (трансформатора) и защиты шин. Указанные подходы могут быть распространены и на другие присоединения подстанции и между подстанциями.
Ток срабатывания дифференциальной токовой защиты в простейшем случае должен быть отстроен от максимального тока небаланса. В случае применения традиционных трансформаторов тока ток срабатывания защиты должен превышать 40% трехфазного тока короткого замыкания в конце линии в максимальном режиме энергосистемы. Из указанного порога срабатывания и значения требуемого коэффициента чувствительности следует, что отношение минимального трехфазного тока короткого замыкания (КЗ) к максимальному току КЗ должно составлять более 0,69. Это существенно ограничивает область применения простых дифференциальных защит. Для повышения чувствительности защит чаще всего применяют дифференциальные защиты с торможением [10]. Другой способ повышения чувствительности дифференциальной защиты связан с использованием адаптивного торможения [11]. При этом используется величина торможения, вычислимая на основе оценки фактических погрешностей измерения тока по концам линии и погрешности синхронизации времени. Такой способ снижает указанное выше соотношение между токами короткого замыкания до значения 0,26. Еще один способ повышения чувствительности дифференциальных защит связан с использованием нелинейного фильтра тока намагничивания [14].
В случае же применения современных датчиков тока и использовании синхрофазоров тока более предпочтительным вариантом следует признать повышение чувствительности за счет сравнения фаз [10,12,13], в данном случае фаз, получаемых на основе синхрофазоров тока. Применение современных датчиков тока и напряжения в сочетании с использованием в защитах информации о синхровекторах тока и напряжения позволяет реализовать эффективную защиту с абсолютной селективностью. Преимущества применения цифровых датчиков тока для дифференциальных защит связаны прежде всего с существенно лучшими их метрологическими характеристиками (снижение тока небаланса), отсутствием остаточной намагниченности и насыщения магнитопровода при больших кратностях тока.
Эффективность рассматриваемой защиты питающей линии с абсолютной селективностью может быть продемонстрирована на примере, приведенном на рис. 2 и рис.
3. Представлено распределительное устройство 6 кВ одного из крупных предприятий Архангельской области. Для упрощения показана только одна секция шин РП, подключенная к шинам генераторного напряжения заводской ТЭЦ. Отходящие от шин присоединения включают 3 двигателя по 3 МВт (показано только 2), трансформатор 3,2 МВА с преобладанием двигательной нагрузки, трансформатор 1 МВА с очень низкой долей двигательной нагрузки.
Рис. 2 - Защита кабельной линии
.500 А(В)
Рис. 3. Защита шин, питающей кабельной линии и отходящих присоединений
Ориентация синхровекторов для наглядности показана относительно синхровектора тока ввода, который в свою очередь для наглядности направлен в нормальном режиме к шинам. Цифрами указаны синхровекторы в следующие моменты времени: 1 - нормальный режим, 2 - по истечении 10 мс после КЗ. Моделирование системы электроснабжения производилось с помощью MATLAB/Simulink с учетом токов подпитки двигателей 6 и 0,4 кВ.
Рис. 4. Импульсные функции аналогового и цифрового фильтров
На рис. 2 приведены синхровекторы тока и напряжения при кз в начале, середине и в конце кабельной линии. Как следует из представленного рис. 2, сравнение модулей и фаз синхровекторов тока позволяет реализовать эффективную защиту с абсолютной селективностью. Дополнительно для повышения чувствительности защиты может быть использована информация о синхровекторах напряжения на шинах генераторного напряжения и шинах РП.
Рис. 5. Импульсные функции аналогового и цифрового фильтров
101 и1(?)
0.05 0.1
101 и2()
-100
- -150
-200
0.05 0.1
■ 1С1 №)
0.05 0.1 0.15 0.2 0.25 0.3
Рис. 6. Графики токов, напряжений, модулей и аргументов синхровекторов
Для эффективного использования СВИ в устройствах РЗА требуется, как указывалось выше, использовать отличные от применяемых в УСВИ цифровые фильтры. Для данной задачи был специально синтезирован аналоговый фильтр-прототип (кривая 1, рис.5) с помощью метода синтеза, основанного на расширенных спектральных представлениях сигналов и при использовании методов нелинейной оптимизации с ограничениями [15]. Для получения описания цифрового фильтра можно воспользоваться известными способами перехода. Но, как показали исследования, для фильтров с ярко выраженной несимметричной импульсной характеристики лучше использовать прямой синтез цифровых фильтров. На рис.4 приведена импульсная функция цифрового фильтра при тех же условиях синтеза, что и для фильтра 1 на рис. 4 при частоте дискретизации 80 точек на период промышленной частоты.
На рис.6 приведены графики напряжения и тока вначале иД?) , ) и в конце линии и2(?), ;2(?), а также графики соответствующим напряжениям и токам (кривые 1) модулей (кривые 2) и аргументов (кривые 3) синхровекторов при коротком замыкании в середине
98
защищаемой линии. Для наглядности огибающая тока в конце линии, определяемая на основе модуля синхровектора тока, сдвинута влево по оси времени на значение, равное групповому времени задержки фильтра.
Как показали исследования, устойчивое срабатывание защиты кабельной линии обеспечивается с большим запасом по истечении течении 10 мс после возникновения КЗ. При этом длительность импульсной функции фильтра равна 27 мс (рис. 4).
В работе [8] показано, что за счет использования синхровекторов может быть реализована защита шин с абсолютной селективностью, сочетающая дифференциальный принцип с контролем аргументов (фаз) и модулей синхровекторов присоединений (рис.3). В дифференциальной защите шин используется только разность тока ввода и суммы токов остальных присоединений. Поэтому использование дополнительной информации об амплитудах и фазах токов присоединений на основе СВИ позволяет повысить чувствительность и быстродействие защиты.
Выводы
Перспективность использования технологии СВИ для релейной защиты понизительных подстанций и РП связана с возможностью широкого внедрения быстродействующих защит с абсолютной селективностью, а также увеличения функциональных возможностей защит. При создании локальных и распределенных систем РЗ с поддержкой технологии СВИ особое внимание следует уделять моделированию защищаемого объекта и защиты, синтезу быстродействующих алгоритмов устройств РЗА.
Литература
1. Synchrophasor Monitoring for Distribution Systems: Technical Foundations and Applications. A White Paper by the NASPI Distribution Task Team., NASPI-2018-TR-001, January 2018.
2. Шнеерсон Э.М., Либах Т. Современные методы фиксации поврежденных фаз и удаленности коротких замыканий // Релейная защита и автоматика энергосистем: сб. докл. конф., М.: ВВЦ, 2004, с. 40-42.
3. Лямец Ю.Я., Нудельман Г.С., Кержаев Д.В., Романов Ю.В. Многомерная релейная защита // Современные направления развития систем релейной защиты и автоматики энергосистем: сб. докл. XIX конф., М.: Н-ИИА, 2009, с. 105-111.
4. Нагай В.И. и др. Информационные признаки аварийного режима при наличии переходного опротивления электрической дуги // Кибернетика энергетических систем. Сб. мат. XXXIX сессии, Новочеркасск, 2018, С.284-288.
5. S. Demri, V. Goranko, M. Lange. Temporal Logics in Computer Science: Finite-State Systems Cambridge. University Press,2016. 752 p.
6. Mi Zeng-Qiang, Zhao Hong-Shan, Wang Hai-Ping, Li Ying. Linear Temporal Logic Control of Power Systems // IEEE/PES Transmission & Distribution Conference & Exposition: Asia and Pacific, Dalian, 2005, pp. 1-5.
7. Мокеев А.В. и др. Расширение функциональных возможностей синхронизированных векторных измерений // Оперативное управление в электроэнергетике. 2016,1, с. 39-47.
8. Мокеев А.В., Xромцов Е.И. Многофункциональные устройства с поддержкой СВИ для цифровых подстанций и цифровых сетей // Релейная защита и автоматизация, 2019. № 4, с. 20-25.
9. Булычев А.В. и др. Релейная защита в распределительных сетях 110/35/10 кВ в условиях цифровой трансформации электроэнергетических систем // Релейная защита и автоматизация. 2019, № 1, с. 70-76.
10. Федосеев А.М. Релейная защита электроэнергетических систем. М.: Энергоатомиздат, 1992, 528 с.
11. Циглер Г. Цифровые устройства дифференциальной защиты. Принципы и область применения. М.: Энергоиздат, 2005, 273 с.
12. Кужеков С.Л. и др. Обеспечение правильного функционирования дифференциальных защит сборных шин в условиях насыщения трансформаторов тока // Изв. Вузов. Электромеханика, 2017, № 1, С.76-83.
13. Басс Э.И., Дорогунцев В.Г. Релейная защита электроэнергетических систем. М.: МЭИ, 2002, 296 с.
14. М.Г. Попов и др. Применение дополнительных тормозных сигналов для повышения чувствительности и быстродействия дифференциальной защиты воздушных линий электропередачи // Изв. вузов. Электромеханика, 2017. Т. 60, № 6. С. 91 - 96.
15. Mokeev A.V. Filter synthesis for PMU // IEEE International Conference SIBC0N-2017. Astana, 2017, pp. 1-6.
© А.В. Мокеев, С.А. Пискунов, Д.Н. Ульянов, Е.И. Хромцов Авторы публикации
Мокеев Алексей Владимирович - д-р техн. наук, профессор кафедры «Энергетики и электротехники» Северного (Арктического) федерального университета.
Пискунов Сергей Александрович — инженер группы комплексного проектирования ООО «Инженерный центр «Энергосервис».
Ульянов Дмитрий Николаевич - заместитель генерального директора ООО «Инженерный центр «Энергосервис».
Хромцов Евгений Иосифович - аспирант, Северного (Арктического) федерального университета.
References
1. Synchrophasor Monitoring for Distribution Systems: Technical Foundations and Applications. A White Paper by the NASPIDistribution Task Team., NASPI-2018-TR-001, January 2018.
2. Shneerson EM, Libakh T. Modern methods of fixing damaged phases and remoteness of short circuits. Releinaya zashchita i avtomatika energosistem: sb. dokl. konf., M.: VVTs. 2004: 40-42.
3. Lyamets YuYa, Nudel'man GS, Kerzhaev DV, Romanov YuV. Multidimensional relay protection. Sovremennye napravleniya razvitiya sistem releinoi zashchity i avtomatiki energosistem: sb. dokl. XIXkonf., M.: N-IIA. 2009: 105-111.
4. Nagai VI, et al. Information signs of emergency mode in the presence of transient resistance of an electric arc. Kibernetika energeticheskikh sistem. Sb. mat. XXXIX sessii, Novocherkassk. 2018: 284288.
5.. Demri S, Goranko V., Lange M.. Temporal Logics in Computer Science: Finite-State Systems Cambridge. University Press,2016. 752 p.
6.Mi Zeng-Qiang, Zhao Hong-Shan, Wang Hai-Ping, Li Ying. Linear Temporal Logic Control of Power Systems. IEEE/PES Transmission & Distribution Conference & Exposition: Asia and Pacific, Dalian, 2005, pp. 1-5.
7. Mokeev AV. Expansion of functionality of synchronized phasor measurements. Operativnoe upravlenie v elektroenergetike. 2016; 1: 39-47.
8. Mokeev AV, Khromtsov EI. Multifunctional devices with SVI support for digital substations and digital networks. Releinaya zashchita i avtomatizatsiya. 2019; 4: 20-25.
9. Bulychev AV, et al. Relay protection in distribution networks 110/35/10 kV in the context of digital transformation of electric power systems. Releinaya zashchita i avtomatizatsiya. 2019; 1: 70-76.
10. Fedoseev AM. Relay protection of electrical power systems. M.: Energoatomizdat, 1992. 528s.
11. Zigler G. Digital devices of differential protection. Principles and scope. M.: Energoizdat, 2005. 273s.
12. Kuzhekov SL. Ensuring the correct functioning of differential protection of busbars in conditions of saturation of current transformers. Izv. Vuzov. Elektromekhanika. 2017; 1:76-83.
13. Bass EI, Doroguntsev VG. Relay protection of electrical power systems. M.: MEI, 2002. 296s.
14. Popov MG, et. al. Application of additional brake signals to increase the sensitivity and speed of differential protection of overhead power transmission lines. Izv. vuzov. Elektromekhanika. 2017; 60; 6: 91-96.
15. Mokeev AV. Filter synthesis for PMU. IEEE International Conference SIBCON-2017. Astana, 2017, pp. 1-6.
Authors of the publication
Aleksei V. Mokeev - Northern (Arctic) Federal University, Arkhangelsk, Russia.
Sergey A. Piskunov - Engineering Centre «Energoservice», Arkhangelsk, Russia.
Dmitry N. Ulyanov - Engineering Centre «Energoservice», Arkhangelsk, Russia.
Evgeny I. Khromtsov - Northern (Arctic) Federal University, Arkhangelsk, Russia.
Поступила в редакцию 29 ноября 2020г.
