Научная статья на тему 'Результаты моделирования процессов пыления хвостохранилищ'

Результаты моделирования процессов пыления хвостохранилищ Текст научной статьи по специальности «Физика»

CC BY
194
63
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПЫЛЯЩИЕ ПОВЕРХНОСТИ / ПОИНТЕРВАЛЬНАЯ ИНТЕНСИВНОСТЬ ПЫЛЕНИЯ / ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ / КОНЦЕНТРАЦИЯ ПЫЛИ / DUSTING SURFACES / INTERVAL DUST INTENSITY / NUMERICAL MODELLING / DUST CONCENTRATION

Аннотация научной статьи по физике, автор научной работы — Маслобоев В.А., Бакланов А.А., Амосов П.В.

Представлены и апробированы подходы (зависимость Westphal D. L. et al. и схема DEAD) по оценке интенсивности пыления (пустыни, хвостохранилища и др.). Продемонстрированпереход к определению динамической скорости u*и скорости на высоте +10 мнад пылящей поверхностьюu10, необходимых для выполнения оценок интенсивностипыления. Методический подход реализован на базе двухмерной численной модели аэродинамики атмосферы района «хвостохранилище АНОФ-2 г. Апатиты». Проанализированы результаты расчетов и обработки горизонтальной скорости на высоте+10 м над пылящей поверхностью при вариации скорости ветрового потока от 5 до 23 м/с. Выполнен сравнительный анализ и указаны особенности поинтервальной (по размерам частиц) интенсивности пыления хвостохранилища АНОФ-2 с использованием зависимости Westphal D. L. et al. и схемы DEAD. Полученные значения интенсивности пыления в нижней части диапазона скорости ветрового потока близки показателю«максимальной удельной сдуваемости пыли», используемого специалистами проектных организаций при разработке документации. Получены прогнозные диапазоны уровней суммарной концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты в виде степенных функциональных зависимостей. Отмечен эффект роста концентрации мелкодисперсной пыли в районе г. Апатиты при увеличении высоты объекта на 20 м.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MODELING RESULTS OF DUST PROCESSES OF TAILING DUMPS

Approaches to evaluation of dust intensity (of deserts, tailing dumps, etc.; dependence of Westphal D. L. et al. and DEAD scheme) are described and tested. The determination of dynamic velocity and velocity at the height of +10 m above the dusting surface, which are necessary to evaluate the dust intensity, is demonstrated. The mathematical approach is realized based on a two-dimensional model of atmosphere aerodynamics in the area “tailing dumps of ANOF-2 Apatity town”. The study provides calculations and analysis of horizontal velocity at the height of +10 m above the dusting surface at the wind speed varying from 5 to 23 m/sec. Comparative analysis is given and the peculiarities of interval (based on grains sizes) dust intensity of the tailing dumps of ANOF-2 are shown using the dependence of Westphal D. L. et al. and DEAD scheme. The received values of dust intensity at the lower range limit are close to the “maximum specific dust off” value which is used by project specialists for documentation development. The forecast ranges of total dust concentrations at the height of +2 m above the dusting surface are received in the center of model of Apatity town in the form of exponentialfunctional dependencies. The effect of fine dust concentrations growth with the height increase above the test object up to 20 m, is registered.

Текст научной работы на тему «Результаты моделирования процессов пыления хвостохранилищ»

ЭКОЛОГИЯ

УДК 622.4 : 519.67

РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ ПРОЦЕССОВ ПЫЛЕНИЯ ХВОСТОХРАНИЛИЩ

В. А. Маслобоев1, А. А. Бакланов1' 2, П. В. Амосов3

Институт проблем промышленной экологии Севера КНЦ РАН 2Всемирная метеорологическая организация 3Горный институт КНЦ РАН

Аннотация

Представлены и апробированы подходы (зависимость Westphal D. L. et al. и схема DEAD) по оценке интенсивности пыления (пустыни, хвостохранилища и др.). Продемонстрирован

переход к определению динамической скорости u* и скорости на высоте +10 м над пылящей поверхностью u10, необходимых для выполнения оценок интенсивности

пыления. Методический подход реализован на базе двухмерной численной модели аэродинамики атмосферы района «хвостохранилище АНОФ-2 - г. Апатиты». Проанализированы результаты расчетов и обработки горизонтальной скорости на высоте +10 м над пылящей поверхностью при вариации скорости ветрового потока от 5 до 23 м/с. Выполнен сравнительный анализ и указаны особенности поинтервальной (по размерам частиц) интенсивности пыления хвостохранилища АНОФ-2 с использованием зависимости Westphal D. L. et al. и схемы DEAD. Полученные значения интенсивности пыления в нижней части диапазона скорости ветрового потока близки показателю «максимальной удельной сдуваемости пыли», используемого специалистами проектных организаций при разработке документации. Получены прогнозные диапазоны уровней суммарной концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты в виде степенных функциональных зависимостей. Отмечен эффект роста концентрации мелкодисперсной пыли в районе г. Апатиты при увеличении высоты объекта на 20 м.

Ключевые слова:

пылящие поверхности, поинтервальная интенсивность пыления, численное моделирование, концентрация пыли.

MODELING RESULTS OF DUST PROCESSES OF TAILING DUMPS

V. А. Masloboev1, А. А. Baklanov1' 2, P. V. Amosov3

institute of North Industrial Ecology Problems 2World Meteorological Organization 3Mining Institute of the KSC of the RAS

Abstract

Approaches to evaluation of dust intensity (of deserts, tailing dumps, etc.; dependence of Westphal D. L. et al. and DEAD scheme) are described and tested. The determination of dynamic velocity and velocity at the height of +10 m above the dusting surface, which are necessary to evaluate the dust intensity, is demonstrated. The mathematical approach is realized based on a two-dimensional model of atmosphere aerodynamics in the area "tailing dumps of ANOF-2 - Apatity town". The study provides calculations and analysis of horizontal velocity at the height of +10 m above the dusting surface at the wind speed varying from 5 to 23 m/sec. Comparative analysis is given and the peculiarities of interval (based on grains sizes) dust intensity of the tailing dumps of ANOF-2 are shown using the dependence of Westphal D. L. et al. and DEAD scheme. The received values of dust intensity at the lower range limit are close to the "maximum specific dust off" value which is used by project specialists for documentation development. The forecast ranges of total dust concentrations at the height of +2 m above the dusting surface are received in the center of model of Apatity town in the form of exponential

functional dependencies. The effect of fine dust concentrations growth with the height increase above the test object up to 20 m, is registered. Keywords:

dusting surfaces, interval dust intensity, numerical modelling, dust concentration.

Введение

Данная работа - естественное продолжение исследований специалистов КНЦ РАН по проблеме пыления хвостохранилищ второй апатито-нефелиновой обогатительной фабрики (АНОФ-2), г. Апатиты, Мурманская обл. на базе методов численного моделирования [1-5]. На предыдущем этапе авторами выполнен анализ существующих подходов по оценке интенсивности пыления (пустыни, хвостохранилища и др.). Некоторые результаты выполненного исследования представлены в докладах на российских и международных конференциях [6-8].

Рассмотрены подходы по оценке интенсивности пыления, приемлемые для решения практических задач по исследованию влияния скорости ветрового потока и высоты пылящей поверхности на загрязнение атмосферы для района г. Апатиты от воздействия близлежащего хвостохранилища АНОФ-2. По результатам анализа сделан предварительный вывод, что наиболее приемлемыми для последующих исследований загрязнения приземного слоя атмосферы вниз по ветровому потоку могут быть признаны зависимость Westphal D. L. et al. [9] и схема DEAD [10]. Именно эти зависимости дают минимальное расхождение результатов на всем рассмотренном диапазоне скорости ветрового потока. Авторы отмечают, что оба подхода базируются на функциональной зависимости потока массы от динамической скорости на высоте пыления в 4-й и 3-й степенях соответственно.

При рассмотрении турбулентного переноса пыли с хвостохранилища АНОФ-2 представляется разумным ограничиться частицами, размер которых не превышает 70 мкм.

Авторский подход к оценке динамической скорости

На рис. 1 представлена схема модели, которая принимается для обоснования величины скорости, необходимой для оценки вертикального потока массы частиц (интенсивности пыления). Основная идея - динамическая скорость на высоте пылящей поверхности u* определяется через осредненную горизонтальную компоненту скорости ветра на высоте +10 м над пылящей поверхностью u10.

Высота Нхв+10м

Поток воздуха > / / Поверхность "Л

Референтная / хвостов Нхв\

$ высота 10м/ \ Основание модели

Рис. 1. Схема модели под обоснование методики определения интенсивности пыления

На входной границе модели рассмотрим диапазон значений референтной скорости ветра на высоте 10 м над основанием модели (от 5 до 23 м/с, например, с шагом 3 м/с). Задавая

u

ref

логарифмический профиль скорости на входной границе области (шероховатость принимается

на уровне 0.05 м), выполняем расчет аэродинамики на базе стационарной стандартной (к -в)-

модели турбулентности [3-5]. Далее выполняется обработка расчетных аэродинамических параметров вдоль пунктирной линии рис. 1: операция осреднения горизонтальной компоненты

скорости (рис. 2) на высоте +10 м над пылящей поверхностью

Результаты обработки

численных экспериментов для горизонтальной скорости на высоте +10 м над пылящей поверхностью, а также переход к динамической скорости и* обобщены в табл. 1.

Рис. 2. Расчетные распределения горизонтальной скорости на высоте +10 м над поверхностью хвостохранилища ы10

Дополнительно для последующего использования при численном решении уравнения конвективно-диффузионного переноса пыли выполнено осреднение (по области моделирования) значений коэффициентов динамической турбулентной вязкости. Результаты этой операции также представлены в табл. 1. Последующий переход к значениям коэффициентов турбулентной диффузии стандартен, он ранее использован и описан авторами [1, 3].

Таблица 1

Осредненные значения горизонтальной скорости на высоте +10 м ы10, динамической скорости и* на высоте поверхности пыления и турбулентной динамической вязкости ^

ы

Скорость ветра, м/с и*, м/с ы10, м/с ^ ,Па-с

5 0.623 8.246 0.4086

8 0.996 13.19 0.6563

11 1.370 18.14 0.9053

14 1.743 23.09 1.1554

17 2.117 28.04 1.4067

20 2.490 32.98 1.6593

23 2.863 37.93 1.9136

Гранулометрический состав (в графической форме) отвальных хвостов с поверхности устоявшегося пляжа хвостохранилища АНОФ-2 приведен в диссертационной работе А. В. Стриженка «Управление экологической безопасностью намывных техногенных массивов ОАО «Апатит» в процессе их формирования» [11, с. 97, рис. 2.8]. Анализ гранулометрического состава отобранных проб пульпы и отвальных хвостов АНОФ-2 проводился на лазерном анализаторе распределения размеров частиц НопЬа LA-950. Обработка графической информации позволила авторам получить количественные показатели по долям («вес» интервала) каждого моделируемого класса размера частиц пыли (табл. 2), что необходимо для выполнения численных экспериментов по турбулентному переносу мультидисперсной примеси. Отметим, что «весовое» распределение интервалов, рассчитанное по данным работы [11], заметно отличается от того распределения, которое авторы использовали в своих прежних работах [1-3] (оно соответствует экспериментальным данным 1970-1980-х гг.). Кроме того, в табл. 2 приведены значения пороговой скорости и* для пыли класса крупности р до 70 мкм в диаметре с шагом 10 мкм, рассчитанные по методике авторов работы [10], и скорости оседания разнокалиберной пыли, которые вычислены в Стоксовом приближении [1-3] и необходимы для последующих вычислений конвективно-диффузионного переноса.

Таблица 2

Значения «веса» интервала, пороговой скорости и скорости оседания в зависимости от срединного диаметра частиц пыли

Срединный диаметр (диапазон интервала), мкм «Вес» интервала, (доли) Пороговая скорость u?, м/с Скорость оседания wg, м/с

5(0-10) 0.022 0.951 0.00195

15 (10-20) 0.083 0.420 0.0175

25 (20-30) 0.142 0.295 0.0487

35 (30-40) 0.194 0.243 0.0955

45(40-50) 0.209 0.218 0.1580

55(50-60) 0.189 0.206 0.2360

65 (60-70) 0.161 0.201 0.3290

Результаты расчетов и их анализ Интенсивность пыления

Результаты вычислений по соотношениям схемы DEAD [10] и зависимости Westphal D. L. et al. [9] поинтервальных (по размерам частиц) вертикальных потоков массы при вариации скорости ветрового потока с использованием необходимых параметров, приведенных в табл. 1 и 2, представлены в табл. 3 и 4 соответственно.

Таблица 3

Поинтервальные интенсивности пыления, вычисленные по схеме DEAD, при вариации скорости ветрового потока, кг/(м2с)

Скорость ветра, м/с Срединный диаметр, мкм

5 15 25 35 45 55 65

5 *** 0.62810-6 0.13510-5 0.19010-5 0.205 10-5 0.18610-5 0.15810-5

8 0.13010-6 0.32910-5 0.57010-5 0.77110-5 0.82310-5 0.74110-5 0.63010"5

11 0.17010-5 0.86710-5 0.14510-5 0.19510-4 0.20810-4 0.18810-4 0.15910-4

14 0.435 10-5 0.17710-4 0.293 10-5 0.39410-4 0.42110-4 0.37910-4 0.32210-4

17 0.829'10-5 0.31110-4 0.51710-5 0.695 10-4 0.743 10-4 0.66910-4 0.56910-4

20 0.13810-4 0.50010-4 0.83110-4 0.11210-3 0.12010-3 0.10810-3 0.91710-4

23 0.21010-4 0.75110-4 0.12510-3 0.16910-3 0.18010-3 0.16310"3 0.13810-3

Таблица 4

Поинтервальные интенсивности пыления, вычисленные по зависимости [9], при вариации скорости ветрового потока, кг/(м2с)

Скорость ветра, м/с Срединный диаметр, мкм

5 15 25 35 45 55 65

5 0.958 10-7 0.36210-6 0.61910-6 0.845 10-6 0.91010-6 0.823 10-6 0.70110-6

8 0.62810-6 0.23710-5 0.405 10-5 0.55410-5 0.59710-5 0.54010-5 0.460-10-5

11 0.225 10-5 0.84710-5 0.14510-4 0.19810-4 0.21310-4 0.19310-4 0.16410-4

14 0.58910-5 0.22210-4 0.38010-4 0.51910-4 0.56010-4 0.50610-4 0.43110-4

17 0.12810-4 0.483 10-4 0.82610-4 0.11310-3 0.12210-3 0.11010-3 0.93710-4

20 0.245 10-4 0.925 10-4 0.15810-3 0.21610-3 0.233 10-3 0.21110-3 0.18010-3

23 0.42910-4 0.16210-3 0.27710-3 0.378 10-3 0.40710-3 0.368 10-3 0.31410-3

Безусловно, основная особенность расчетных интенсивностей пыления, выполненных по соотношениям схемы [10], - отсутствие числового значения при скорости ветра 5 м/с для самой мелкой пыли (срединный диаметр 5 мкм). Данный факт имеет очевидное объяснение: величина пороговой скорости (табл. 2, 3-й столбец, 2-я строка) выше значения динамической скорости (табл. 1, 2-й столбец, 2-я строка), т. е. в рамках данной модели (схема DEAD) эмиссии пыли такой крупности не происходит. Кроме того, хорошо видно, что при малых скоростях ветрового потока (не более 8 м/с) для пыли практически любой крупности значения интенсивности пыления, рассчитанные по схеме [10], несколько выше аналогичных расчетных данных по зависимости [9]. При величине скорости ветрового потока 11 м/с результаты расчетов интенсивности пыления по обоим используемым подходам наиболее близки. Дальнейшее увеличение скорости ветра (встречается существенно реже) приводит к тому, что расчетные (по зависимости [9]) интенсивности пыления заметно превышают аналогичные величины, вычисленные с использованием схемы [10].

Дополнительно были проанализированы суммарные интенсивности пыления при вариации либо скорости ветрового потока (строго по решаемой задаче), либо расчетной величины динамической скорости на высоте пыления соответственно. До величины скорости ветрового потока 11 м/с доминирует кривая интенсивности пыления по схеме [10], а при более высокой скорости - кривая по зависимости [9]. Естественно, что анализируемые кривые описываются степенными функциями AVn в 4-й (зависимость [9]) и 3-й (схема [10]) степенях или скорости ветрового потока (в рамках решаемой задачи), или расчетной величины динамической скорости на высоте пыления. Коэффициенты пропорциональности A в степенных зависимостях от динамической скорости следующие: схема [10] - 3.877 и зависимость [9] - 2.896. Размерность скорости V м/с.

Естественно, авторы понимают, что для верификации численной модели необходимо проведение трудоемких экспериментальных исследований на объекте пыления, что, к сожалению, пока не может быть реализовано по целому ряду причин. Тем не менее предпринята попытка проверить объективность полученных значений интенсивности пыления. Для этого проанализированы данные отчета «Реконструкция хвостохранилища до отметки 200 м: проектная документация. Раздел 8 «Перечень мероприятий по охране окружающей среды», посвященного реконструкции хвостохранилища АНОФ-2 до высотной отметки 200 м. Специалисты ЗАО «Механобр Инжиниринг» используют показатель «максимальная удельная сдуваемость пыли». В своих оценках авторы указанного отчета используют значение

5.2910-6 кг/(м2с), что отвечает значению динамической скорости не выше 0.6 м/с. Представляется возможным сделать два вывода: 1) предлагаемый методический подход достаточно объективен и работоспособен; 2) при высоких скоростях ветрового потока прогнозируемые значения интенсивности пыления будут существенно выше тех величин, которые принимаются проектировщиками в оценках воздействия объектов пыления на окружающую среду.

Распределение концентрации пыли

Расчетные данные интенсивности пыления, представленные в табл. 3 и 4, необходимы для решения конвективно-диффузионного уравнения распространения пыли при описании граничных условий на поверхности пыления. В силу того что указанное уравнение - линейное, достаточно будет выполнить численные эксперименты по пространственно-временному распределению концентрации пыли в исследуемой области, например, для данных табл. 4 (зависимость [9]). Зная соответствующие отношения интенсивностей пыления в табл. 3 и 4, получаем значения концентрации пыли, отвечающие эмиссии пылевых частиц по схеме DEAD. В результате прогноз пространственно-временного распределения концентрации пыли будет определяться в некотором диапазоне значений.

В качестве примера выполненных расчетов приведены графики поинтервального пространственного распределения концентрации пыли на высоте +2 м над поверхностью г. Апатиты для скорости ветрового потока 5 м/с (рис. 3). Пространственные координаты модели соответствуют примерным границам г. Апатиты, а значения интенсивности пыления вычислены по зависимости [9].

113 146000

Й

^ 3,E-07 Int 1

Z с " Г Г 1 _|nt_ 2

к § -LLI

S J Ю ■ UUI Int 3

to а ^ 1С Ш Int .4

I ас Int 5

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(U -UUI

J

I о Int_ _6

-in Int_ 7

isi

3,E-08

Ось Х, м

Рис. 3. Поинтервальное пространственное распределение концентрации пыли на высоте +2 м над поверхностью для района г. Апатиты при скорости ветрового потока 5 м/с (оценка интенсивности пыления по зависимости [9])

Далее были проанализированы суммарные концентрации пыли, соответствующие выбранным значениям скорости ветрового потока и вычисленные с помощью используемых подходов и приемов. На рис. 4а и б представлены суммарные концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты при вариации либо скорости ветрового потока (строго по решаемой задаче), либо расчетной величины динамической скорости на высоте пыления (более общая форма) соответственно. Как видно на рис. 4а и б, поведение кривых

концентраций пыли следует за поведением кривых интенсивности пыления: до величины скорости ветрового потока 11 м/с (и соответствующей динамической скорости) превалирует кривая, отвечающая схеме [10], а при более высокой скорости - кривая по зависимости [9]. Пересечение кривых происходит вблизи отметки 11 м/с (и соответствующей динамической скорости).

а

с s 1/Е-04 ю ас , , , , { , , , , 1 , , , , , , , 1 , , , 1 1 , , , , , 1 14 17 ifK----23

..........

ас

S с;

JD с | 1,Е-05 —ф— CWEST

-■-CDEAD

о. 1- ■ X

о 4 зе

1.Е-06

Скорость ветра, м/с

б

1,Е-04

1,Е-05

1,Е-06

-CWEST -CDEAD

Динамическая скорость, м/с

Схема [10] - ромб; зависимость [9] - квадрат Рис. 4. Суммарные концентрации пыли в зависимости от:

а - скорости ветрового потока на референтной высоте; б - динамической скорости на высоте пыления

Достаточно очевидно, что кривые на рис. 4а и б описываются, как и функции интенсивности пыления, посредством степенных функций со своими коэффициентами аппроксимации. Например, для кривых рис. 4б наилучшим образом подходят следующие: по схеме DEAD CD = 1.21 • 10 5 • ul'77; по зависимости [9] CW = 8.96 • 10 6 • u,2-81.

Таким образом, используемый методический подход позволяет спрогнозировать потенциальную область значений концентрации пыли, ограниченную графиками, которые соответствуют различным формализациям интенсивности пыления. На рис. 4а и б такое представление выполнено в терминах суммарной концентрации пыли в зависимости либо от скорости ветрового потока, либо от величины динамической скорости на высоте пыления.

Эффект высоты хвостохранилища

Обосновав и апробировав описанный выше «инструмент» для прогноза уровня концентрации пыли в районе г. Апатиты при вариации скорости ветрового потока, авторы,

естественно, не могли не применить описанный выше подход для различных высот хвостохранилища. Напомним, что ранее авторы попытались ответить на вопрос о влиянии изменения высоты хвостохранилища АНОФ-2 со 180 до 200 м (в абсолютных высотных отметках) на загрязнение приземного слоя атмосферы г. Апатиты [3-5]. По оценкам, сделанным на базе трехмерного численного моделирования, для мелкодисперсной пыли (до 60 мкм) при скорости ветра 6 м/с увеличение высоты хвостохранилища на 20 м должно привести к росту концентрации пыли примерно на 15 %.

Поскольку в данном исследовании для высоты хвостохранилища 200 м модель была «заточена» под скорость ветрового потока 5 м/с, авторами выполнена подготовка компьютерной модели и проведена серия расчетов для высоты хвостохранилища 180 м при той же скорости ветрового потока. При этом соблюдалось условие постоянства площади пыления за счет изменения размера пруда-отстойника. Результаты расчетов поинтервальной интенсивности пыления (зависимость [9]) и концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты обобщены в табл. 5. Графическое сравнение поинтервальной и суммарной концентрации пыли в той же точке модели изображено на рис. 5.

Как видно из данных табл. 5, увеличение высоты хвостохранилища на 20 м при референтной скорости ветра 5 м/с приводит к росту интенсивности пыления примерно на 17 % и увеличению уровня концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты примерно на 20 %.

Согласно приведенным результатам, предыдущие авторские оценки и данные текущего исследования несколько отличаются, но это различие не столь существенное. В обоих случаях отмечается безусловное увеличение концентрации мелкодисперсной пыли с ростом высоты хвостохранилища. При этом именно в представляемой работе получены конкретные величины концентрации пыли, имеющие достаточную объективность.

Таблица 5

Поинтервальные интенсивности пыления (кг/(м2с), вычисленные по зависимости [9], и концентрации пыли (кг/м3) на высоте +2 м над поверхностью в центре модели г. Апатиты при

скорости ветрового потока 5 м/с

Номер (диапазон, мкм) интервала Интенсивность пыления, кг/(м с) Концентрация пыли, кг/м3

180 м 200 м 180 м 200 м

1 (0 - 10) 8.1610-8 9.5810-8 3.0510-8 3.6510-8

2 (10 - 20) 3.0810-7 3.6210-7 1.2010-7 1.4410-7

3 (20 - 30) 5.2710-7 6.1910-7 2.2110-7 2.6610-7

4 (30 - 40) 7.2010-7 8.4510-7 3.3610-7 4.0710-7

5 (40 - 50) 7.7510-7 9.1010-7 4.1510-7 5.0410-7

6 (50 - 60) 7.0110-7 8.2310-7 4.4110-7 5.4010-7

7(60 - 70) 5.9710-7 7.0110-7 4.5010-7 5.5010-7

В заключение отметим, что представленные на рис. 5 суммарные концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты при референтной скорости ветра 5 м/с для двух значений высоты хвостохранилища (180 и 200 м) в соответствии с излагаемым методическим подходом можно интерпретировать следующим образом. Для высоты хвостохранилища 180 м прогнозируемый диапазон суммарной концентрации пыли составляет 2.0110-6 - 4.5210-6 кг/м3, а для высоты 200 м - 2.4110-6 - 5.5110-6 кг/м3. Первое число в этих парах отвечает интенсивности пыления, рассчитанной по зависимости [9], второе - по схеме [10].

Рис. 5. Поинтервальное и суммарное распределение суммарной концентрации пыли на высоте +2 м в центре модели г. Апатиты при референтной скорости ветра 5 м/с (проварьирована высота хвостохранилища 180 и 200 м)

Выводы

Продемонстрирован достаточно простой и нетрудоемкий переход к определению динамической скорости u* и скорости на высоте +10 м над пылящей поверхностью u10

по данным о скорости ветра набегающего потока и высоты дамбы. Указанные величины необходимы для выполнения прогнозных расчетов вертикального потока массы на базе созданной двухмерной численной модели аэродинамики атмосферы исследуемой области.

Переход реализован на базе двух подходов по оценке интенсивности пыления (зависимость [9] и схема [10]), которые приемлемы для решения практической задачи о влиянии скорости ветрового потока и высоты пылящей поверхности на загрязнение атмосферы для района г. Апатиты от воздействия близлежащего хвостохранилища АНОФ-2. Оба подхода базируются на функциональной зависимости потока массы от динамической скорости на высоте пыления в 4-й и 3-й степенях соответственно.

Выполнен сравнительный анализ расчетных значений поинтервальной интенсивности пыления на базе выбранных подходов при вариации скорости ветрового потока в широком диапазоне. Указан экономичный путь дальнейших численных экспериментов по прогнозу пространственно-временного распределения поинтервальной концентрации пыли, основанный на линейной природе конвективно-диффузионного уравнения и выполненном анализе расчетных вертикальных потоков массы пыли с использованием зависимости [9] и схемы [10].

Выполнены расчеты уровней концентрации приземного слоя атмосферы для района г. Апатиты с использованием принятых подходов по описанию интенсивностей пыления в широком диапазоне скоростей ветрового потока и при максимальной проектной высоте хвостохранилища АНОФ-2. Получены прогнозные диапазоны уровней суммарной концентрации пыли на высоте +2 м над поверхностью в центре модели г. Апатиты в форме аналитических зависимостей степенного вида.

Эффект высоты хвостохранилища на уровень суммарной концентрации пыли продемонстрирован для ситуации изменения высоты объекта со 180 до 200 м при скорости ветрового потока 5 м/с. В подтверждение ранее выполненных оценок доказано: увеличение высоты пылящей поверхности приводит к росту концентрации мелкодисперсной пыли примерно на 20 %.

ЛИТЕРАТУРА

1. Бакланов А. А. Численное моделирование в рудничной аэрологии. Апатиты: КФ АН СССР, 1987. 200 с.

2. Baklanov A., Rigina O. Environmental modeling of dusting from the Imming and concentration sites in the Kola Peninsula, Northwest Russia // The XI World Clear Air and Environment Congress (Durban, South Africa, IUAPPA-NACA, 14-18 September 1998). Durban, 1998. Vol. 1. 4F-3. Р. 1-18. 3. Амосов П., Бакланов А, Ригина О. Численное моделирование процессов пыления хвостохранилищ. LAP LAMBERT Academic Publishing, 2014. 109 с. 4. CFD-

model developing of dust transfer at a tailings dump / P. V. Amosov [et al.] // Proceedings of the 4-th International Conference on Hazardous and Industrial Waste Management - CRETE-2014 (Chania, Crete, Greece, 2-5 September 2014,) Executive Summaries. Chania: Technical University of Crete, 2014. Р. 279-280. 5. Численное моделирование процессов пыления хвостохранилища АНОФ-2 / В. А. Маслобоев [и др.] // Вестник МГТУ. 2014. Т. 17, № 2. С. 376384. 6. Амосов П. В., Бакланов А. А. К вопросу оценки интенсивности пыления хвостохранилищ // Математические методы в технике и технологиях - мМтТ-28: сб. тр. XXVIII Междунар. науч. конф.: в 12 т. Т. 1 / под общ. ред. А. А. Большакова. Саратов: Саратов. гос. техн. ун-т, 2015; Ярославль: Ярослав. гос. техн. ун-т; Рязань: Рязанск. гос. радиотехн. ун-т, 2015. С. 3-5. 7. Амосов П. В., Бакланов А. А. Оценка вертикального потока массы пыли на хвостохранилище на базе схемы DEAD // Проблемы недропользования. 2015. № 3. С.80-85. 8. Amosov P. V, Baklanov A. A. Assessment of dusting intensity on ANOF-2 tailing by using a Westphal D. L. dependency // Proceedings of X International Symposium on Recycling Technologies and Sustainable Development (Bor, Serbia; [organizer] University of Belgrade, Technical Faculty, 4-7 November 2015) / ed. by Zoran S. Markovic. Bor: University of Belgrade, 2015 (Bor : Tercija). P. 39-43. 9. Westphal D. L., Toon O. B., Carlson T. N. A case-study of mobilization and transport of Saharan Dust // Journal of the Atmospheric Sciences. 1988. No. 45. Р. 2145-2175. 10. Marticorena B., Bergametti G. Modeling the atmospheric dust cycle.1. Design of a soil-drived dust emission scheme // Journal of Geophysical Research-Atmospheres. 1995. Vol. 100, No. D8. Р. 16415-16430. 11. Стриженок А. В. Управление экологической безопасностью намывных техногенных массивов ОАО «Апатит» в процессе их формирования: дис. ... канд. техн. наук. 2015 // Национальный минерально-сырьевой университет «Горный»: сайт. URL: http://www.spmi.ru/system/files/lib/sci/aspirant-doctorant/avtoreferaty/2015/dissertaciya_strizhenok.pdf (дата обращения: 12.09.2015).

Сведения об авторах

Маслобоев Владимир Алексеевич - доктор технических наук, профессор, директор ФГБУН Института проблем промышленной экологии Севера КНЦ РАН; e-mail: masloboev@ksc.ru

Бакланов Александр Анатольевич - доктор физико-математических наук, профессор, главный научный сотрудник ФГБУН Института проблем промышленной экологии Севера Кольского научного центра РАН, сотрудник научного отдела Всемирной метеорологической организации; e-mail: aabaklanov@yahoo.com

Амосов Павел Васильевич - кандидат технических наук, старший научный сотрудник ФГБУН Горного института Кольского научного центра РАН; e-mail: vosoma@goi.kolasc.net.ru

Information about the authors

Vladimir A. Masloboev - Dr. Sci. (Eng.), Professor, director of the Institute of North Industrial Ecology Problems of the KSC of the RAS; e-mail: masloboev@ksc.ru

Alexander A. Baklanov - Dr. Sci. (Physics and Mathematics), Professor, chief researcher of Institute of North Industrial Ecology Problems of the KSC of the RAS; researcher at the World Meteorological Organization; e-mail: aabaklanov@yahoo.com

Pavel V. Amosov - PhD (Eng.), senior researcher of the Mining Institute of the KSC of the RAS; e-mail: vosoma@goi.kolasc.net.ru

Библиографическое описание статьи

Маслобоев В. А. Результаты моделирования процессов пыления хвостохранилищ / В. А. Маслобоев, А. А. Бакланов, П. В. Амосов // Вестник Кольского научного центра РАН. - 2016. -№ 1 (24). - С. 59-68.

Bibliographic Description

Vladimir А. Masloboev, Alexander А. Baklanov, Pavel V. Amosov. Modeling Results of Dust Processes of Tailing Dumps. Herald of the Kola Science Centre of the RAS. 2016, vol. 1 (24), pp. 59-68.

i Надоели баннеры? Вы всегда можете отключить рекламу.