CC BY
42
Вестник РУДН. Серия: Экология и безопасность жизнедеятельности
RUDN Journal of Ecology and Life Safety
2018 Том 26 № 3 354-366
http://journals.rudn.ru/ecology
DOI 10.22363/2313-2310-2018-26-3-354-366 УДК 57.049
Прогнозирование экологических последствий загрязненности воздуха свалочным газом в жилой застройке
Санкт-Петербургский университет технологий управления и экономики Российская Федерация, 190103, Санкт-Петербург, Лермонтовский пр., д. 44, лит. А
В статье выполнено прогнозирование экологических последствий загрязненности компонентами свалочного газа воздуха в жилой застройке, расположенной на территории вблизи полигонов захоронения отходов. Предложена методика расчета рассеивания свалочного газа с учетом влияния застройки, основанная на численном решении уравнений диффузии, теплообмена, Рейнольдса, показательного распределения, дополненных соответствующими моделями турбулентности, начальными и граничными условиями. Методика позволяет учитывать пространственную неоднородность полигонов и свалок как источников загрязняющего вещества, а также процесс проникновения загрязненного воздуха в помещения. Выполнена программная реализация и верификация методики на примере полигона «Новоселки» (г. Санкт-Петербург) и прилегающих к нему складских зданий и жилого комплекса для наиболее неблагоприятного ветрового режима территории с учетом осреднения по времени. Проведен эколого-токсикологический анализ свалочного газа с точки зрения основных его компонентов, а также анализ канцерогенного и неканцерогенного рисков при хроническом воздействии на абстрактного человека, проживающего в помещении, в которое проникает загрязненный свалочным газом воздух. Результаты работы могут быть использованы для определения величины предельно допустимого выброса компонентов свалочного газа, размеров санитарно-за-щитных зон полигонов и свалок, выявления наиболее загрязненных мест на улице и в помещении, оценки риска для здоровья жильцов в результате хронического действия ядовитых веществ, синдрома «больного здания», вызванного проникновением загрязненного воздуха в помещения через неплотности ограждающих конструкций, при решении прикладных задач экологической безопасности строительства, жилищно-коммунального хозяйства и общественного здравоохранения.
Ключевые слова: свалочный газ, полигоны захоронения отходов, синдром «больных зданий», канцероген, хронический риск
Введение. Свалочный газ представляет собой газовую смесь, образующуюся в результате разложения на полигонах и свалках твердых бытовых и промышленных отходов (ТБО и ПО), представляющую опасность для окружающей среды. В составе свалочного газа помимо метана и углекислого газа могут присутствовать аммиак, бензол, сероводород, бенз(а)пирен, фенол, ксилол, диоксины, меркаптаны, эфиры, алкинбензолы и другие экотоксиканты. В условиях уплотнительной
А.В. Сауц
© Сауц А.В., 2018
This work is licensed under a Creative Commons Attribution 4.0 International License
городской застроики расположение здании, в том числе жилых и торговых комплексов, нередко оказывается в зоне негативного влияния полигонов и свалок. На практике соблюдение границ санитарно-защитных зон часто не обеспечивает соблюдение требований к качеству воздуха в застройке. В результате в зоне влияния полигонов создается неблагоприятная экологическая обстановка. Стандартная методика [1], действующая в настоящее время на территории РФ, а также ее прототип ОНД-86, реализованные в программном комплексе УПРЗА «Эколог», не учитывают процесс проникновения загрязненного воздуха в помещения, не предназначены для расчета рассеивания загрязняющих веществ от пространственно-неоднородных источников и прогнозирования экологических последствий загрязненности воздуха.
Цель исследования — прогнозирование экологических последствий загрязненности воздуха свалочным газом в жилой застройке.
Методы и результаты исследования. Расчет диффузии пассивной примеси концентрацией с, мг/м3, производится по уравнению
д-С+Щсу-k Vc) dt ( kxyz Vc) dz
—+V(cV -kxyz Vc)--^ = I, (1)
где I — время, с; V = V (и; V; м>) — вектор скорости движения воздуха, м/с; I — эмис
кхуг=ку; К
сия вещества, мг/(м2с); к ? = diag{kx; ку; к7} — кинематическая турбулентная вязкость,
м2/с; ж, — скорость оседания или подъема примеси в воздухе, м/с.
При расчете поля скоростей движения наружного воздуха приняты допущения, что среда является изотермической и несжимаемой. Для расчетов применяются следующие подходы.
1. При большом удалении полигона от здания используется система уравнений Навье — Стокса для несжимаемой вязкой жидкости:
P^V+P(V? -V)V = v[n(vV+VVT)-p ] dt , (2)
vV = о
где р — плотность, кг/м3; п — динамическая вязкость, Пас; р — давление, Па.
Для определения турбулентной структуры потока система (2) дополняется метеорологической моделью турбулентности Дж. Смагоринского:
кх _ ку _Ах АУкбаз/ А 2 к _ кбаз ^2/ А 2 ,
где Д — масштаб ячейки, м; кбаз — базовый коэффициент турбулентной вязкости, м2/с, определяемый по формулам
_ |кф +еА2\Defl VI -Ш /Л, Ш <1
кбаз Ъ, т > 1 . (4)
Масштаб А прямоугольных ячеек размерами Ах, Ау, Аг определяется как: — при умеренной анизотропии сетки
А = ^АхАуАг;
(5)
при сильной анизотропии сетки
А = тах{Ах; Ау; Аг},
(6)
где кф = 1—15 м2/с; е = 0,1—0,4 — константа Смагоринского; Бв/ — диссипативная или деформационная функция:
Бв/ =
/^..л2 ^.л2 /дмл2
ди дх
+
ду \дУу
+
дг
1
+2
ду ди
— + —
дх ду
2
+
( дм ду ^2 ( ди дм ^2
- + -ду дг
+
+
дг дх
(7)
Ш — критерий Ричардсона, определяемый по формуле
ш=£ д-
дг
+
[у ф(г) -У а
/[Бв/ ]2,
(8)
где Уф(г) = —йТ^/йг — градиент фоновой температуры, °/м; Т — температура воздуха в расчетной точке, К; Тф — фоновая температура, К; уа = 0,01 °/м — адиабатический градиент; Ш « 0 при Т« Тф.
2. При проведении аэродинамических расчетов непосредственно вблизи застройки модель Смагоринского может давать ложные результаты. В связи с этим вблизи препятствий формируются дополнительные расчетные подобласти (рис. 1), в которых для расчетов используется уравнение Рейнольдса, дополненное моделью турбулентности «к-е»:
Ре ду+ро^-Уу=у
УК = 0
(п + п,)(+УУТ)-р
дк — / — — т\2
Р0 ^ + Р0УУк-У[(П + П, /^к)Ук] = 0,5Л, (( + УУ ) -р0£ ,
де £ / Т \2 £2
Р0 ^ + Р0УУе = У[(П + П, /)Уе] = 0,5Се1 ( + УУ ) -р0Се2 к
П, =Р0 Сц к 2/ е
(9)
где к — кинетическая энергия турбулентности; е — скорость диссипации турбулентности; Сп = 0,09; Се1 = 1,44; Се1 = 1,92; ок = 1; ое = 1,3.
Кинематическая турбулентная вязкость определяется по выражению
кхуг = П,/Р0.
(10)
Рис. 1. Построение дополнительной расчетной подобласти вблизи зданий [Figure 1. The construction of additional settlement subareas in close proximity to buildings]
В помещениях зданий расчет вентиляционных течений производится с помощью системы уравнений Навье — Стокса и теплопроводности в приближении Буссинеска — Обербека:
Ро dV + Po(V?-V)V = v
n(VK + VV )-p
+ Ро g[Pr (T - To) -1]
дТ dt VV = 0
pCp -V(XVT) + pcpVVT = Q
(11)
где р0 — плотность воздуха, кг/м3, при некой равновесной температуре Т0, К; вт — коэффициент термического расширения воздуха, К-1; — ускорение свободного падения, м/с2; ср — изобарная теплоемкость, Дж/(кг К); 1 — коэффициент теплопроводности Вт/(м-К); Q — источник или сток тепла, Вт/м2.
Для численного решения приведенных дифференциальных уравнений задаются соответствующие начальные и граничные условия.
Влияние полигона или свалки как источника примеси определяется с помощью уравнения [2]
I = 11574,07^/(ТтеПл^оЛ), (12)
где т — эмиссия свалочного газа, т/год (рассчитывается по методике АКХ имени К.Д. Памфилова); Z — высота массива отходов, м; Ттепл — период года, дней, во время которого Тф > 273,15 К; VПоЛ — объем полигона или свалки, м3.
Эмиссию веществ, поступающих в помещение за счет неплотности ограждающей конструкции (оконный блок), можно найти как [2]
1 = у^ритДрАрн^ (13)
где у — коэффициент пересчета разовой концентрации в среднесуточную [3]; сприт — разовая концентрация загрязняющего вещества в приточном воздухе, мг/м3; Др — перепад давлений воздуха между наружной и внутренней поверхностями ограждающих конструкций, Па; рн — плотность наружного воздуха, кг/м3; Ru — сопротивление воздухопроницанию ограждающих конструкций, кг/(м2 с).
Эмиссию загрязняющих веществ, поступающих в помещение через открытое окно или вентиляционный клапан, можно найти как
С = УСприт. (14)
На входной границе наружной расчетной области профиль скорости ветра задается уравнением Кармана:
и = и*1п[(г - Нср + г/О/г)]/?, (15)
где и* — скорость ветра, м/с; Нср — средняя высота препятствий на пути ветра, м; г0 — шероховатость поверхности, м; ? = 0,4 — коэффициент Кармана; С — коэффициент сопротивления, м.
Скорость поступления приточного воздуха принимается равномерной по всей площади оконного блока и определяется как
V = Ар/(рнШи). (16)
Постановка задачи принята стационарной, ее решение выполняется с использованием метода установления. Программная реализация модели выполнена в среде Сот8о1Ми1Ирйу81с8 3.5а для наиболее неблагоприятного метеорежима для полигона «Новеселки» (г. Санкт-Петербург), на расстоянии 500 м от которого расположен автомобильный завод «Ниссан Мэнуфэкчуринг Рус», а на расстоянии 5 км — жилой комплекс «Дом на Нижне-Каменской улице» (рис. 2).
Рис. 2. Расположение полигона «Новоселки», завода «Ниссан Мэнуфэкчуринг Рус» и жилого комплекса «Дом на Нижне-Каменской улице» [Figure 2. The location of the "Novoselki" landfill, plant "Nissan Manufacturing Rus" and residential complex "Apartment house on the Nizhne-Kamenskaya street"]
Полученные результаты и их обсуждение. На рис. 3 приведены результаты расчетов концентраций аммиака и бензола в зоне влияния полигона «Новоселки».
Более детальный результат расчетов концентраций аммиака и бензола вблизи жилого комплекса приведен на рис. 4.
Расчеты загрязнения воздуха аммиаком и бензолом также произведены для помещений здания на первом этаже жилого дома с наветренной и подветренной сторон (рис. 5). Результаты расчетов приведены на рис. 6 и 7.
Результаты расчетов концентраций сравнивались с замерами, произведенными с газоанализатором УГ-2. В табл. 1 сравнение разовых приземных концентраций в наружном воздухе приведено в долях ПДК максимальной разовой (ПДКмр), в табл. 2 — среднесуточных концентраций в центре помещения с наветренной стороны в долях ПДК среднесуточной (ПДКсс).
с Ш
Рис. 3. Поле разовых приземных концентраций аммиака (а) и бензола (б), мг/м3, в районе полигона «Новоселки» (11.03.2014) [Figure 3. Field of single surface concentrations ammonia (a) and benzene (б), mg per m3, near the landfill "Novoselk" (11.03.2014)]
Рис. 4. Поле разовых приземных концентраций аммиака (а) и бензола (б), мг/м3, вблизи жилого комплекса (11.03.2014) [Figure 4. Field a single ground level concentrations ammonia (a) and benzene (б), mg per m3, near the residential compound (11.03.2014)]
Рис. 5. Места расположения исследуемых помещений в многоквартирном доме: 1 — с наветренной стороны; 2 — с подветренной стороны [Figure 5. Locations of the studied premises in an apartment building: 1 — on the windward side; 2 — on the leeward side]
Рис. 6. Поле среднесуточных концентраций аммиака (а) и бензола (б), мг/м3, в помещении с наветренной стороны (11.03.2014) [Figure 6. Field of average daily concentrations of ammonia (a) and benzene (б), mg per m3, in the room from the windward side (11.03.2014)]
Рис. 7. Поле среднесуточных концентраций аммиака (а) и бензола (б), мг/м3, в помещении с подветренной стороны (11.03.2014) [Figure 7. Field of average daily concentrations of ammonia (a) and benzene (б), mg per m3 in the room from the leeward side (11.03.2014)]
Таблица 1
Значения разовых приземных концентраций компонентов свалочного газа в наружном воздухе, доли ПДКмр, полученные различными способами
Компонент свалочного газа ПДКМр [4] Вблизи складов автомобильного завода Вблизи жилого комплекса
Экспериментальное Расчетное Экспериментальное Расчётное
Аммиак 0,2 61,5 ± 15,4 75,0 24,5 ± 6,1 34,5
Бензол 0,3 39,7 ± 9,9 47,3 16,74 ± 4,2 1,7
Table 1
The value of a single near-surface concentrations of components of landfill gas in the outer air, the share of the maximum single MPC, obtained in different ways
Component of landfill gas Maximum single MPC [4] Near the warehouses of the automobile plant Near the residential complex
Experimental Calculated Experimental Calculated
Ammonia 0.2 61.5 ± 15.4 75.0 24.5 ± 6.1 34.5
Benzene 0.3 39.7 ± 9.9 47.3 16.74 ± 4.2 1.7
Таблица 2
Значения среднесуточных концентраций компонентов свалочного газа в воздухе в центре помещения с наветренной стороны, доли ПДКсс, полученные различными способами
Компонент свалочного газа ПДКсс [4] Значения, доли ПДКсс
Экспериментальное Расчетное
Аммиак 0,04 32,5 ± 8,1 32,5
Бензол 0,1 10,0 ± 2,5 12,0
Table 2
The average daily concentrations of landfill gas components in the air in the center of the room from the windward side, the proportion of MPC average daily, obtained in various ways
Component of landfill gas MPC average daily [4] Values of the share of the average daily MPC
Experimental Calculated
Ammonia 0.04 32.5 ± 8.1 32.5
Benzene 0.1 10.0 ± 2.5 12.0
Теперь проанализируем последствия загрязненности воздуха свалочным газом в застройке. При длительном воздействии на человеческий организм метан способен вызывать асфиксию, изменение кровяного давления и головную боль [5]. Углекислый газ у человека и животных также способен вызывать асфиксию и кроме того возбуждающе действовать на дыхательные центры в головном мозге, что приводит к увеличению частоты и глубины дыхания [6]. Аммиак у людей и животных в больших концентрациях способен вызывать раздражение дыхательных путей, пульмонологические заболевания, тошноту, рвоту, повышенную нервную возбудимость, бессонницу и т.д. Сероводород при длительном вдыхании вызывает головную боль, слезотечение, раздражение в носу, при длительном воздействии — усталость, тошноту, головную боль [5]. Перечисленные симптомы, наблюдаемые у жильцов исследуемого дома, являются наиболее типичными при возникновении синдрома «больного здания». При этом без проведения соответствующих расчетов или замеров будет достаточно трудно установить истинную причину его появления и развития, так как некоторые компоненты свалочного газа не имеют запаха или имеют высокий обонятельный порог чувствительности (кроме аммиака и сероводорода).
Оценить неканцерогенный риск R при хроническом воздействии экотокси-кантов можно с помощью беспороговой модели воздействия [4]:
bC
r = 1 - о,84 K зПДКсс, (17)
где С — среднесуточная концентрация загрязняющего вещества, мг/м3; b — коэффициент изоэффективности, величина которого зависит от класса опасности вещества (b = 2,4 для 1 класса, b = 1,31 для 2 класса, b = 1 для 3 класса, b = 0,86 для 4 класса); K — коэффициент запаса, величина которого зависит от класса опасности вещества (K = 7,5 для 1 класса, K3 = 6 для 2 класса, Кз = 4,5 для 3 класса, K3 = 3 для 4 класса).
Некоторые компоненты биогаза, например бензол, обладают канцерогенным действием. Оценить индивидуальный канцерогенный риск CR которых можно, дополнив построенную математическую модель следующим уравнением [5]
CR = 1 - в-^1^0, (18)
где SF — фактор канцерогенного потенциала, кг день/мг; ЛАББ — среднесуточная доза в течение жизни, мг/(кг день), определяемая как
ЛАБО = , (19)
365ВЖ-АТ
где c — концентрация вещества в загрязненной среде, мг/м3; CR — скорость поступления воздействующей среды (воздуха), м3/день; ED — продолжительность воздействия, лет; EF — частота воздействия, дней/год; BW — масса тела человека, кг; AT— период усреднения экспозиции (для канцерогенов AT = 70 лет).
Расчет рисков по имеющимся полям концентраций компонентов биогаза выполнен в программной среде Surfer 11. В качестве примера на рис. 8 и 9 приведе-
ны результаты расчетов для исследуемого помещения с наветренной стороны. При этом для расчета канцерогенного риска, вызванного парами бензола при ингаляционном воздействии, рассматривается среднестатистический абстрактный взрослый житель массой 70 кг, вдыхающий 20 м3 воздуха в день на протяжении 337 дней в году (предполагается, что отпуск жилец проводит вне дома), проживающий в квартире в течение 70 лет.
Рис. 8. Распределение неканцерагенного риска при хроническом воздействии аммиака (а) и бензола (б) в помещении с наветренной стороны (центральный разрез вдоль помещения) [Figure 8. Distribution of non-carcinogenic risk in chronic exposure to ammonia (a) and benzene (б) in the room from the windward side (central incision along the room)]
Рис. 9. Распределение канцерогенного риска, вызванного парами бензола при ингаляционном воздействии в помещении с наветренной стороны (центральный разрез вдоль помещения) [Figure 9. Distribution of carcinogenic risk caused by benzene vapors under inhalation exposure in the room from the windward side (central incision along the room)]
В помещении с подветренной стороны диапазон канцерогенного риска допустим для профессиональных групп, но не допустим для жильцов (10~4 < СЯ < < 10~3) [7]. Диапазон неканцерогенного риска при хроническом воздействии аммиака и бензола не приемлем ни для профессиональных групп, ни для жильцов (Я > 10~3) [7].
Заключение. В данной работе выполнено прогнозирование экологических последствий загрязненности воздуха свалочным газом в жилой застройке. Использование предложенной методики расчета позволяет выявить наиболее загрязнен-
ные места, выполнить оценку экологических последствий загрязненности воздуха, рассчитать величину предельно допустимого выброса (ПДВ), размеров санитарно-защитных зон, что представляет практическую значимость в решении задач экологической безопасности строительства, ЖКХ, общественного здравоохранения.
Список литературы
[1] Приказ Министерства природных ресурсов и экологии РФ от 06 июня 2017 г. № 273 «Об утверждении методов расчетов рассеивания выбросов вредных (загрязняющих) веществ в атмосферном воздухе».
[2] Сауц А.В. Обеспечение экологически безопасного воздушного режима зданий, расположенных вблизи полигонов твердых бытовых и промышленных отходов: дис. ... канд. техн. наук. СПб., 2014. 140 с.
[3] Теверовский Е.И., Терновский И.А. Допустимые выбросы радиоактивных и вредных химических веществ в приземный слой атмосферы. М.: Энергоатомиздат, 1980. 240 с.
[4] ГН 2.1.6.3492-17. Предельно допустимые концентрации (ПДК) загрязняющих веществ в атмосферном воздухе городских и сельских поселений.
[5] Балахчина Т.К. Оценка воздействия свалочного газа с полигонов твердых бытовых отходов на человека // Физиология. Медицина. Экология человека. 2012. № 2. С. 41—57.
[6] Лазарев Н.В., Левина Э.Н. и др. Вредные вещества в промышленности: справочник для химиков, инженеров и врачей. 7-е изд., пер. и доп.: в 3 т. Т. III. Неорганические и элемен-торганические вещества. Ленинград: Химия, 1977. 608 с.
[7] Плуготаренко Н.К., Чижова В.И., Петров В.В. и др. Применение методик оценки риска для анализа качества окружающей среды // Технологии техносферной безопасности. 2014. № 2 (54). URL: http://www.agps-2006.narod.ru/ttb/2014-3/30-03-14.ttb.pdf (дата обращения: 25.03.2018).
[8] Р 2.1.10.1920-04. Руководство по оценке риска для здоровья населения при воздействии химических веществ, загрязняющих окружающую среду.
История статьи:
Дата поступления в редакцию: 14.06.2018 Дата принятия к печати: 30.07.2018
Для цитирования:
Сауц А.В. Прогнозирование экологических последствий загрязненности воздуха свалочным газом в жилой застройке // Вестник Российского университета дружбы народов. Серия: Экология и безопасность жизнедеятельности. 2018. Т. 26. № 3. С. 354—366. DOI 10.22363/2313-2310-2018-26-3-354-366
Сведения об авторе:
Сауц Артур Валерьевич — кандидат технических наук, доцент кафедры маркетинга и социальных коммуникаций ЧОУ ВО «Санкт-Петербургский университет технологий управления и экономики». Контактная информация: e-mail: mining_expert@mail.com; eLIBRARY SPIN-код: 5191-0602
Forecasting the environmental effects of air pollution with landfill gas in residential areas
A.V. Sauts
Saint-Petersburg University of Management Technologies and Economics 44 Lermontovsky prospect, litera A, Saint Petersburg, 190103, Russian Federation
In the article the prediction of ecological consequences of contamination of components of the landfill gas air to a residential building, located on the territory close to the landfill. A method for calculating the dispersion of landfill gas taking into account the influence of the building, based on the numerical solution of differential equations, heat transfer, Reynolds, exponential distribution, supplemented by corresponding models of turbulence, initial and boundary conditions. The technique allows to take into account the spatial heterogeneity of landfills and rolling pins as the source of the pollutant, as well as the process of penetration of polluted air into the interference. Performed software implementation and verification of the method on the example of the range "Novoselki" (Saint Petersburg) and the adjacent warehouse buildings and residential complex for the most unfavourable wind conditions taking into account the averaging over time. The ecological and toxicological analysis of landfill gas taking into account its main components, and also the analysis of carcinogenic and non-carcinogenic risk at chronic influence on the abstract person living in the room in which the air polluted by landfill gas penetrates is carried out. The results can be used to determine the amount of maximum permissible emissions of landfill gas components, the scale moat sanitary protection zones of landfills and dumps, detecting the identification of the most polluted places on the street and indoors, the health risk assessment of residents as a result of chronic action of toxic substances, the syndrome of "sick building" caused by penetration of contaminated air into the room through openings in the walling, the solution of applied problems of ecological safety of construction, housing and communal services and public health.
Keywords: landfill gas, landfill sites, syndrome of "sick buildings", carcinogen, chronic risk
References
[1] Order of the Ministry of Natural Resources and Environment of the Russian Federation of June 06, 2017 No. 273 On approval of calculation methodsfor dispersion of emissions of harmful (polluting) substances in the atmospheric air.
[2] Sauts AV Ensuring environmentally safe air regime of buildings located near landfills for solid household and industrial waste (dissertation of cand. tech. sciences). Saint Petersburg; 2014.
[3] Teverovsky EI, Ternovsky IA. Permissible emissions of radioactive and harmful chemicals into the surface layer of the atmosphere. Moscow: Energoatomizdat Publ.; 1980.
[4] Hygienic standard GN 2.1.6.3492-17. Maximum Permissible Concentrations (MPCs) of Pollutants in the Air of Urban and Rural Settlements.
[5] Balakhchina TK. Assessment of the impact of landfill gas from solid waste landfills on humans. Physiology. The medicine. Human ecology. 2012;(2): 41—57.
[6] Lazarev N, Levina E et al. Harmful substances in industry. Handbook for chemists, engineers and doctors. 7th ed., trans. and add: in 3 vol. Vol. III. Inorganic and elemental substances. Leningrad: Chemistry Publ.; 1977.
[7] Plugotarenko NK, Chizhova VI, Petrov VV et al. Application of risk assessment techniques for environmental quality analysis. Technosphere Safety Technologies. 2014;2(54). http://www.agps-2006.narod.ru/ttb/2014-3/30-03-14.ttb.pdf (access date: 25.03.2018).
[8] Guideline R 2.1.10.1920-04. Guidelines for assessing the risk to public health from exposed chemicals that pollute the environment.
Article history:
Received: 14.06.2018 Revised: 30.07.2018
For citation:
Sauts AV. Forecasting the environmental effects of air pollution with landfill gas in residential areas. RUDN Journal ofEcology and Life Safety. 2018;26(3): 354-366. DOI 10.22363/2313-23102018-26-3-354-366
Bio Note:
Arthur Valerevich Sauts — Candidate of Technical Sciences, Associate Professor of Marketing and Social Communications Department, Private Educational Institution of Higher Education "Saint Petersburg University of Management Technologies and Economics". Contact information: e-mail: mining_expert@mail.com; eLIBRARY SPIN-code: 5191-0602
