ГИАБ. Горный информационно-аналитический бюллетень / MIAB. Mining Informational and Analytical Bulletin, 2020;(3-1):37-45 ОРИГИНАЛЬНАЯ СТАТЬЯ / ORIGINAL PAPER
УДК 622.838:911.6:519.72 DOI: 10.25018/0236-1493-2020-31-0-37-45
РАЙОНИРОВАНИЕ ПОДРАБОТАННЫХ ТЕРРИТОРИЙ ГОРОДА БЕРЁЗОВСКИЙ
А.С. Ведерников1, П.И. Зуев1
1 Институт горного дела Уральского отделения Российской академии наук (ИГД УрО РАН)
Аннотация: районирование подработанных территорий с начала истории наблюдения за деструктивными проявлениями земной поверхности в подработанном горном массиве в виде оседаний, а в худшем случае и провалов, остаётся актуальной проблемой для населённых пунктов, находящихся вблизи либо в пределах разрабатываемых месторождений. Особую актуальность она представляет для городов «заводского типа» Урала, где исторически промышленные поселения развивались вместе и на основе горнометаллургических производств, иногда переплетаясь в территориальном плане. По прошествии многих лет обоснованное и безопасное использование этих территорий вызывает серьезные проблемы ввиду разрозненности и отсутствии информации об их состоянии. Одним из таких поселений является город Берёзовский, большая часть территории которого подработана. Закрытие и последующее затопление части шахт когда-то градообразующего золотодобывающего рудника привело к необходимости исследования территории, в ходе которого с использованием предоставленной рудником инженерно-геологической, геодезической и прочей информации была разработана методика и поэтапно произведено районирование в геоинформационной среде, установлены последствия затопления. Результаты работы являются уникальными ввиду неповторимости и разнородности параметров подземных выработок, провалов на земной поверхности, вызванных выработками, тесно сосредоточенных на исследуемой территории.
Ключевые слова: подработанная территория, геоинформационная система, дайка, провал, массив, слой, город.
Для цитирования: Ведерников А.С., Зуев П.И. Районирование подработанных территорий города Берёзовский // Горный информационно-аналитический бюллетень. — 2020. — № 3-1. — С. 37-45. DOI: 10.25018/0236-1493-2020-31-0-37-45.
Zoning of undermined urban territories in the area of the Berezovsky town
A.S. Vedernikov1, P.I. Zuev1
1 The Institute of Mining of the Ural branch of the Russian Academy of Sciences, Russia
Abstract: Zoning of undermined urban territories since the beginning of observations over disruption of undermined ground surface in the form of subsidence or sinking, in the worst case, yet remains critical for the population clusters nearby or within mineral mining areas. The problem is of the special concern for industrial towns in the Ural where historically manufacturing communities developed alongside with mining and metallurgical industries, and in the same territory. After decades and sometimes ages have passed, safe and efficient use of these areas remains problematic in view of the disembodied or missing information about their condition. One of such residential areas is the town of Berezovsky with the major portion
© А.С. Ведерников, П.И. Зуев. 2020.
of its territory being undermined. Closure and subsequent flooding of a once-urban gold mine necessitated the territorial research. During the research using the geotechnical, geodetic and other type information provided by the mine engineering service, the zoning procedure was developed and implemented in the geoinformation environment, and the consequences of the flooding were defined. The research results are unique in view of the singularity and diversity of characteristics of underground excavations and induced ground surface sinks clustered inside the study territory.
Key words: Undermined territory, geoinformation system, dyke, sink, rock mass, layer, town. For citation: Vedernikov A.S., Zuev P.I. Zoning of undermined urban territories in the area of the Berezovsky town. MIAB. Mining Inf. Anal. Bull. 2020;(3-1):37-45. [In Russ]. DOI: 10.25018/02361493-2020-31-0-37-45.
Введение
Окончание разработки месторождений и затопление рудников порождают большой комплекс работ и мероприятий. В том числе для дальнейшей эксплуатации подработанных территорий требуется определить состояние дневной поверхности подработанных территорий. Город Березовский Свердловской области практически полностью располагается в пределах разрабатываемого уже несколько столетий месторождения.
Постановка задачи
В данной работе рассматривается оценка катастрофических последствий затопления части шахт Березовского рудника для социальной и промышленной инфраструктуры города. Эта специфика взаимного расположения города и подработанной подземной добычей территории выдвигают на первый план развитие процесса сдвижения горных пород в качестве источника катастрофических явлений. Опасность его состоит в том, что на застроенной городскими объектами различного назначения подработанной территории, сохраняющей в настоящее время относительную устойчивость, при обводнении подработанных пород налегающей толщи, особенно верхней зоны выветривания, процессы сдвижения получат свое дальнейшее развитие, в том числе с развитием воронок обрушения,
провалов, трещин и плавных деформаций, превышающих предельные допустимые значения. В связи с внезапным характером проявлений деформационных процессов они представляют непосредственную опасность для населения города и сопряжены со значительным материальным ущербом [1, 2].
Вторым по значимости источником катастрофических явлений, пожалуй, является подтопление определенных районов города [3], абсолютные отметки территорий которых окажутся ниже абсолютных отметок уровня подземных вод, образовавшегося в результате ликвидации установившейся депрессионной воронки при работающем водоотливе рудника. Фактор подтопления проявится в проникновении подземных вод в подвалы, обводнении и нарушении фундаментов и в тяжелых случаях — в непосредственном затоплении и заболачивании.
Третьим фактором по значимости и вероятности проявления в качестве источника катастрофических явлений можно назвать сейсмогеодинамические процессы, представляющие собой техногенные землетрясения. Они возникают от нарушения исходного напряженно-деформированного состояния из-за горных работ и формирования вторичного техногенно измененного напряженного состояния. Накопивши-
еся упругие деформации, вызванные концентрацией напряжений и деформаций вокруг выработанных пространств, при обводнении массива могут проявиться подвижками по тектоническим нарушениям. Уровень сейсмического воздействия на здания и сооружения зависит от скорости и величины подвижки и чаще всего не вызывает серьезных разрушений. Но здания и сооружения, оказавшиеся непосредственно над активной зоной, могут получить серьезные нарушения от произошедших деформаций [4—6].
Проявления этих трех основных источников катастрофических явлений осложняются двумя взаимосвязанными специфическими факторами. Во-первых, город полностью расположен на территории месторождения. Во-вторых, месторождение разрабатывается на протяжении 265 лет с полным разнообразием технологий добычи, присущим различным периодам его существования: «Казенная добыча» золота более 150 лет велась в приповерхностной части месторождения на глубинах, не превышавших 45 — 50 м [7]. За этот период было построено более 1000 мелких шахт. Кроме того, повсеместно велась браконьерская старательская добыча мелкими горными выработками. В период гражданской войны все эти разработки подвергались затоплению.
В 20—30-е годы прошлого столетия рудник был восстановлен, построены новые шахты с современными технологиями добычи, сопровождавшимися образованием на земной поверхности воронками обрушения, провалами и формированием мульд сдвижения. В настоящее время разработка месторождения осуществляется двумя шахтами: «Южная» и «Северная», ведущими горные работы на глубинах 314 и 512 м соответственно.
Таким образом, инфраструктура города Березовского теснейшим образом переплетается с областями негативного воздействия горных разработок. В связи с этим, любые изменения в горных разработках, тем более такие радикальные, как затопление горных выработок и обводнение вмещающих пород, требуют серьезной проработки специализированными организациями.
В геологическом плане рудное поле исследуемого участка расположено в пределах Урало-Тобольского антикли-нория в опущенном блоке, ограниченном субмеридиональными разломами и крупными телами древних габбро и гиперба-зитов (с юга и севера) [8], герцинскими гранитоидами (с запада и востока). Площадь рудного поля составляет около 60 км2. Блок сложен нижнепалеозойскими слоистыми толщами (кремнисто-углистые сланцы, филлиты с прослоями туффитов, в нижней части разреза — метаморфизованные диабазы и порфи-риты), вмещающими силлы и лакколи-тообразные тела серпентинизированных перидотитов и пироксенитов. Породы прорваны многочисленными дайками плагиопорфиров, гранитпорфиров, лам-профиров различного возраста. Протяжённость даек до 20 км, мощность 2 — 40 м, простирание близмеридиональное, падение крутое. К поперечным трещинам разрыва в дайках приурочены короткие крутопадающие кварцеворудные «лестничные» жилы, менее развиты пологие широтные и диагональные.
Жилы в пределах рудного поля распространены неравномерно. Широко развиты околожильные изменения — лиственитизация гипербазитов и оса-дочно-вулканогенных пород, березити-зация жильных гранитоидов. Рудные жилы сложены крупнозернистым кварцем, карбонатом, сульфидами. Основные рудные минералы — пирит, на отдельных участках — галенит,
тетраэдрит, айкинит, халькопирит, самородное золото и другие. Золото (проба 800 — 900) распределено неравномерно в виде тонкой дисперсной вкрапленности (золотоносность ранних выделений сульфидов) и более крупных скоплений (поздние генерации сульфидов), на верхних горизонтах иногда в виде мелких самородков [8].
Отрабатываются скопления жил в дайках и во вмещающих породах. Месторождение разрабатывается подземным способом в основном с магани-зированием руды, реже — с подэтажной отбойкой и выемкой горизонтальными слоями с гидрозакладкой выемочного пространства (на нижних горизонтах). Запасы рудного золота составляют 90 т, ежегодный объем добычи — 600 кг. Таким образом, при сохранении текущих объемов производительности предприятие обеспечено сырьевой базой на 150 лет [9].
По гидрогеологическим данным, основным источником обводнения месторождения являются трещинно-жильные воды, приуроченные к тектоническим нарушениям, дайкам грано-диоритов и кварцевым жилам, которые представляют собой естественные подземные дрены и в то же время являются проводниками поверхностных вод [3]. Низкое гипсометрическое положение Березовского месторождения в замкнутой котловине обусловило создание здесь благоприятных условий для накопления подземных вод за счет инфильтрации атмосферных осадков с дренируемой водосборной площади.
Выделяются две вертикальные фильтрационные зоны:
• верхняя — зона грунтового стока;
• нижняя — зона локальных трещинных систем тектонического генезиса.
В верхней зоне проницаемость коллекторов определяется экзогенной
трещиноватостью, тектонической тре-щиноватостью и техногенной трещино-ватостью, сформированной процессами сдвижения горных пород при отработке месторождения.
На глубинах более 220 м начинается зона локальных трещинных систем тектонического генезиса. Эти системы имеют, как правило, линейный характер, сложную геометрию, разнообразные размеры и фильтрационные параметры. Очагами питания здесь являются участки выхода в зону грунтового стока. Гидравлическая связь между отдельными трещинными системами обычно отсутствует или реализуется через зону грунтового стока.
Фоновые значения водопроводимо-сти и обводненности разреза в целом имеют сравнительно низкие величины. Водопроводимость изменяется от первых единиц до 50 — 100 м2/сут, водо-притоки к скважинам в стационарном режиме не превышают 1 — 3 л/с.
Зоны разрывных тектонических нарушений и контактов интрузивных тел с толщами иного литолого-петро-графического состава характеризуются повышенной обводненностью.
Методы исследования
Целью проведенных исследований является разработка метода визуализации информации о характере и степени подработанности участков земной поверхности, позволяющего обоснованно решать вопросы о вовлечении подработанных территорий в хозяйственное использование. Решение задачи осуществляется методами геоинформатики, применение которых обеспечивает обоснованный подход к решению проблемы безопасного вовлечения в хозяйственную деятельность ранее подработанных территорий.
Актуальность данная работа получила в связи с вопросом прекраще-
ния финансирования откачки воды из отработанных шахт рудника. Администрацией города была поставлена задача — на основании существующих геодезических и гидрогеологических изысканий определить последствия заполнения водой пустот подработанных территорий города. В рамках задачи по имеющимся данным изысканий были проведены исследования и анализ подработанных территорий для структурирования параметров и прогнозирования геомеханического состояния исследуемого горного массива.
На руднике велось многолетнее наблюдение за подработанными территориями города, документировалась информация о провалах дневной поверхности. Также имелись данные о размерах учтенных подземных выработок — планы и проекции шахт, карты с выделенными участками подработки. По предварительным данным было определено, что подработанная территория занимает примерно 70 % от всей площади города. Поэтому создание геоинформационного пакета на основе этой информации стало актуальной задачей для районирования данных территорий.
Для создания пакета использовалась открытая геоинформационная система Quantum GIS (QGIS), распространяемая бесплатно и совместимая с большинством распространенных на рынке аналогичных инструментов. Функционально данная система полностью соответствует поставленной задаче [10, 11].
Как было сказано ранее, одним из источников информации явились журналы регистрации разнообразных нарушений и провалов дневной поверхности подработанной территории, которые ведутся на руднике с 1957 года. Провалы были распределены по двум категориям — временной период и зафиксированный диа-
метр воронки. За основу взята карта города Березовский, на ней отдельными слоями были вынесены гидрообъекты, растительность, жилая и производственная инфраструктура, зоны землепользования, дороги, гидросеть.
Результаты исследований
и обсуждение
В результате анализа журналов были определены параметры зафиксированных деформаций и составлена таблица для всей территории. Она содержит в себе следующие параметры: 1) дата образования провала, 2) привязка на местности, 3) привязка к вмещающей дайке, 4) размер — диаметр, 5) глубина, 6) широта, 7) долгота. На основе таблицы для статистического удобства сеть провалов была разбита на 5 диапазонов по временным отрезкам. Каждый диапазон включает в себя 5 слоев с ранжированными по диаметру провалами. Далее были созданы таблицы атрибутов, которые обеспечили визуализацию статистически удобных слоев, упомянутых ранее. Таким образом, была сформирована карта провалов, которая отображена на рис. 1. В левом верхнем углу на карту вынесены слои: 1) 5 диапазонов по временным периодам: 1957 — 1970 гг., 1971 — 1980 гг., 1981 — 1990 гг., 1991 — 2000 гг., 2011 — 2014 гг.; 2) распределение провалов по диаметру воронки провала в метрах: 0 — 12, 12 — 24, 24—36, 36 — 48, 48 — 60.
В ходе разработки пакета был создан слой даек в пределах месторождения. Информация по отработанным дайкам была получена от специалистов рудника, в её основе лежат проекционные разрезы интрузивных тел. Была произведена статистическая оценка отработки даек посредством геоинформационного обеспечения Quantum GIS. В итоге были сформированы два слоя — пространственное положение даек и степень их
Puc. 1. Карта провалов, в левом верхнем углу указаны 5 временных диапазонов, которые включают в себя по 5 слоёв с провалами, распределенными по диаметру воронки обрушения Fig. 1. Map of sinkholes, in the upper-left corner there are 5 time ranges, which include 5 layers with sinkholes distributed along the diameter of the collapse funnel
подработки. Степень подработки разграничили четырьмя категориями 0-25%, 25-50%, 50-75%, 75-100% [6]. После получения соотношения провалов и отработанных интрузивных тел, нанесённых на одну карту, была установлена их взаимосвязь на уровне техногенного генезиса.
Большинство провалов вполне ожидаемо лежат вблизи наиболее подработанных даек, что в принципе обуславливает их образование на дневной поверхности.
В ходе дальнейшего развития геоинформационного пакета был получен совмещённый план города Берё-зовский с зафиксированными на нём участками, на которых велись горные работы. Таким образом, при совмещении трёх основных тематических слоёв (зафиксированных деформа-
ций, степени подработанности даек и известных горных выработок) аналитическим методом удалось оконтурить потенциально опасный район развития геомеханических процессов техногенного происхождения (рис. 2) [12, 13].
Заключение
Исследование подтвердило возможность возникновения прогнозируемых опасных явлений при затоплении рудника. Возникновение интенсивных процессов сдвижения горных пород неизбежно вследствие обводнения породы, подверженных выветриванию.
Геоинформационный пакет включен в более масштабный проект Института горного дела УрО РАН, объединяющий объекты горнодобывающей промышленности.
Рuc. 2. Потенциально опасный район развития процессов сдвижения пород, выделен прозрачной областью на территории города
Fig. 2. A potentially dangerous area for the development of rock shifting processes is highlighted by a transparent area on the territory of the city
По результатам использования геоинформационного пакета в прогнозировании сдвижения пород на совещании с руководством города было принято решение продолжить финансирование откачки воды из отработанных пустот. Данный проект продолжает развиваться, геомеханическая картина изменяется постоянно, внесение новых данных способствует более глубокому осмыслению процессов сдвижения и формированию более детализированного представления о них, выража-
ющегося в закономерностях и предписаниях безопасности [14, 15].
Авторы выражают благодарность руководителю направления геомеханика в Институте горного дела УрО РАН — А.Д. Сашурину за советы при подготовке статьи, также благодарность выражается всем сотрудникам отдела геомеханики Института, участвовавшим в данной работе. За предоставление информационных материалов выражается благодарность ООО «Берёзовский рудник».
СПИСОК ЛИТЕРАТУРЫ
1. Turney J.E. Subsidence above inactive coal mines: information for the homeowner. Department of Natural Resources. Colorado Division of Reclamation, Mining & Safety and the Colorado Geological Survey. — Denver, USA. 1985. — 36.
2. Ефремов Е.Ю. Обоснование критерия завершения процесса воронкообра-зования. Известия Тульского Государственного Университета. Науки о земле. —
2018. №4. — С. 12 — 21.
3. Далатказин Т.Ш., Коновалова Ю.П. Прогноз последствий затопления Березовского рудника. Проблемы недропользования. — Екатеринбург. — 2017. — №3. — С. 60 — 66.
4. Сашурин А.Д. Геомеханика в горном деле: фундаментальные и прикладные исследования. Горный журнал. — Екатеринбург. — 2017. — №1. — С. 29—32.
5. Сашурин А.Д. Формирование напряженно-деформированного состояния иерархически блочного массива горных пород. Проблемы недропользования. — Екатеринбург. — 2017. — №1. — С. 38 — 44.
6. Гоигорьев Д.В. Результаты исследования массива горных пород в районе строительства Екатеринбургского метрополитена. Материалы XII Уральской молодёжной научной школы по геофизике. — Пермь: Горный институт УрО РАН,
2011. — С. 59—62.
7. Клейменов Д.А., Альбрехт В.Г., Ерохин Ю.В., Баталин А.С., Баталина А.А. Берёзовское золоторудное месторождение. История и минералы. Екатеринбург: ФГУИПП «Уральский рабочий», 2005. — 200 с. с ил.
8. Бородаевский Н.И., Бородаевская М.Б. Березовское рудное поле. — М: Гос-геолтехиздат, 1947. — 261с.
9. Уралу не хватает запасов. URL: http://news.rambler.ru/other/38549807-uralu-ne-hvataet-zapasov
10. Панжин А.А., Зуев П.И. Задачи развития горно-геологических информационных систем для учета взаимодействия объектов недропользования с природной средой. В сборнике: Уральский научный форум Всероссийская научная конференция с международным участием «Дни наук о Земле на Урале». — Екатеринбург,
2012. С. 92 — 94.
11. Цветков В.Я. Информационные модели и геоинформационные модели. Образовательные ресурсы и технологии. — 2016. № 4 (16). — С. 114—120. DOI: 10.21777/2312 — 5500—2016 — 4-114 — 120.
12. Jixiong Zhang, Qiang Zhang, Qiang Sun, Rui Gao, Deon Germain, Sami Abro. Surface subsidence control theory and application to backfill coal mining technology. Environmental Earth Sciences, 2015, Vol. 74, no 2, p. 1439.
13. Jangwon Suh, Yosoon Choi, Hyeong-Dong Park. GIS-based evaluation of mining-induced subsidence susceptibility considering 3D multiple mine drifts and estimated mined panels. Environmental Earth Sciences, 2016, Vol. 75, No. 10, p. 1.
14. Muhammad Ahsan Mahboob, Bekir Genc, Turgay Celik, Sarfraz Ali, Iqra Atif. Modeling and analysis of Lily gold mine disasters using geoinformatics. GeoJournal,
2019, р. 1.
15. Radosiaw Murdzek, Hubert Malik and Andrzej Lesniak. Ground subsidence information as a valuable layer in GIS analysis. E3S Web of Conf. 2018, Vol. 36, No. 02006, p. 8. DOI: 10.1051/e3sconf/20183602006. EES
REFERENCES
1. Turney J.E. Subsidence above inactive coal mines: information for the homeowner. Department of Natural Resources. Colorado Division of Reclamation, Mining & Safety and the Colorado Geological Survey. Denver, USA. 1985. 36.
2. Efremov E.YU. Justification of the criterion for completing the funnel formation process. Izvestiya Tul'skogo Gosudarstvennogo Universiteta. Nauki o zemle. 2018. no 4. P. 12 — 21. [In Russ]
3. Dalatkazin T.SH., Konovalova YU.P. Forecast of the consequences of flooding of the Berezovsky mine. Problems of subsurface use. Problemy nedropol'zovaniya. 2017. no 3. pp. 60 — 66. [In Russ]
4. Sashurin A.D. Geomechanics in mining: fundamental and applied research. Gornyj zhurnal. Ekaterinburg. 2017. no 1. P. 29 — 32. [In Russ]
5. Sashurin A.D. Formation of the stress-strain state of a hierarchically blocky rock mass. Problemy nedropol'zovaniya. Ekaterinburg. 2017. no 1. pp. 38 — 44. [In Russ]
6. Grigor'ev D.V. Rezul'taty issledovaniya massiva gornyh porod v rajone stroitel'stva Ekaterinburgskogo metropolitena [Results of the study of the rock mass in the construction area of the Yekaterinburg metro]. Materialy XII Ural'skoj molodyozhnoj nauchnoj shkoly po geofizike. Perm': Gornyj institut UrO RAN, 2011. P. 59 — 62. [In Russ]
7. Klejmenov D.A., Al'brekht V.G., Erohin YU.V., Batalin A.S., Batalina A.A. Beryozovskoe zolotorudnoe mestorozhdenie. Istoriya i mineraly [Berezovskoe gold Deposit. History and minerals]. Ekaterinburg: FGUIPP «Ural'skij rabochij», 2005. 200 p. s il. [In Russ]
8. Borodaevskij N.I., Borodaevskaya M.B. Berezovskoe rudnoe pole [Berezovsky ore field]. Moscow: Gosgeoltekhizdat, 1947. 261 p. [In Russ]
9. Uralu ne hvataet zapasov. URL: http://news.rambler.ru/other/38549807-uralu-ne-hvataet-zapasov [In Russ]
10. Panzhin A.A., Zuev P.I. Zadachi razvitiya gorno-geologicheskih informacionnyh sistem dlya ucheta vzaimodejstviya ob'ektov nedropol'zovaniya s prirodnoj sredoj. Vsbornike [Problems of development of mining and geological information systems for accounting of interaction of subsurface use objects with the natural environment. In the collection]: Ural'skij nauchnyj forum Vserossijskaya nauchnaya konferenciya s mezhdunarodnym uchastiem «Dni nauk o Zemle na Urale». Ekaterinburg, 2012. P. 92 — 94. [In Russ]
11. Cvetkov V.YA. Information models and geoinformation models. Obrazovatel'nye resursy i tekhnologii. 2016. no 4 (16). pp. 114 — 120. DOI: 10.21777/2312 — 5500 — 2016 — 4-114 — 120. [In Russ]
12. Jixiong Zhang, Qiang Zhang, Qiang Sun, Rui Gao, Deon Germain, Sami Abro. Surface subsidence control theory and application to backfill coal mining technology. Environmental Earth Sciences, 2015, Vol. 74, no 2, p. 1439.
13. Jangwon Suh, Yosoon Choi, Hyeong-Dong Park. GIS-based evaluation of mining-induced subsidence susceptibility considering 3D multiple mine drifts and estimated mined panels. Environmental Earth Sciences, 2016, Vol. 75, No. 10, p. 1.
14. Muhammad Ahsan Mahboob, Bekir Genc, Turgay Celik, Sarfraz Ali, Iqra Atif. Modeling and analysis of Lily gold mine disasters using geoinformatics. GeoJournal, 2019, p. 1.
15. Radostaw Murdzek, Hubert Malik and Andrzej Lesniak. Ground subsidence information as a valuable layer in GIS analysis. E3S Web of Conf. 2018, Vol. 36, No. 02006, p. 8. DOI: 10.1051/e3sconf/20183602006.
ИНФОРМАЦИЯ ОБ АВТОРАХ
Ведерников Андрей Сергеевич1 — научный сотрудник, [email protected], Зуев Павел Игоревич1 — младший научный сотрудник, [email protected], 1 Институт горного дела Уральского отделения Российской академии наук (ИГД УрО РАН), 620075 г. Екатеринбург, ГСП-219, Мамина-Сибиряка 58.
INFORMATION ABOUT THE AUTHORS
Vedernikov A.S.1, Researcher, [email protected], Zuev P.I.1, Junior Researcher, [email protected],
1 The Institute of Mining of the Ural branch of the Russian Academy of Sciences, 620075, Ekaterinburg, Russia
Получена редакцией 21.11.2019; получена после рецензии 20.01.20; принята к печати 20.03.2020. Received by the editors 21.11.2019; received after the review 20.01.20; accepted for printing 20.03.2020.