Методика расчета оценки ветровой эрозии и пыления

Утверждена приказом
Председателя Комитета
государственного энергетического
надзора Министерства индустрии и
новых технологий
Республики Казахстан
от «24» ноября 2010 года
№112-П
Методика расчета оценки ветровой эрозии и пыления золоотвалов
тепловых электростанций
Содержание
1
2
3
4
5
5.1
Введение
Область применения
Нормативные ссылки
Термины и определения
Обозначения и сокращения
Основные положения
Механизм и основные показатели ветровой эрозии золошлакоотвала
1
2
2
3
4
4
ТЭС
4
11
13
18
5.2 Характеристика золоотвала как площадного пылящего объекта
5.3 Расчетная оценка ветровой эрозии золовых шлаков ТЭС
Приложение 1 Значения поправочных коэффициентов К3 и К4
Приложение 2 Приведение экспериментально определенной удельной
сдуваемости к скорости ветра на высоте флюгера
Приложение 3 Удельная сдуваемость золовых отложений на отвалах
действующих ГРЭС
Приложение 4 Алгоритм и пример расчетного определения ветровой
эрозии золошлакоотвала
Библиография
19
21
22
29
Введение
Одним из существенных воздействий тепловых электростанций (ТЭС) на
окружающую среду является вынос в атмосферу пылевых частиц с поверхности
золошлакоотвалов (ЗШО) ТЭС в результате их ветровой эрозии и последующее
их осаждение на почве. Пыление возникает из-за несовершенства проектных
решений и технологии складирования золошлаков, нарушения правил
эксплуатации указанных объектов, являющихся источником неорганизованных
выбросов загрязнителей в атмосферу.
2
Оценка пылевых выбросов с поверхности ЗШО ТЭС рекомендуется
экспериментальными, аналоговыми, расчетными методами. Наиболее полная и
достоверная оценка пылящих свойств золоотвалов рекомендуется полевыми
исследованиями при правильной их организации, с использованием
специальной аппаратуры. Универсальный способ оценки пылящих свойств
ЗШО - расчетный метод, основанный на определенной физической модели
процесса ветровой эрозии.
Настоящая Методика расчета оценки ветровой эрозии и пыления
золоотвалов ТЭС (далее - Методика) составлена на основе теоретических,
лабораторных и полевых исследований, обобщения данных выполненных
наблюдений, исследований, проектных разработок в отечественной и
зарубежной энергетике, смежных отраслях. [1]; [2]; [3]; [5].
1. Область применения
Настоящая Методика рекомендуется:
- ТЭС и энергетическим объединеним для оценки воздействия ЗШО ТЭС
на состояние окружающей природной среды и разработки мероприятий по
ограничению ее загрязнения, для определения валового выброса пылевых
частиц в атмосферу с указанных поверхностных источников с целью
формирования отчетности и расчетов платежей за выбросы;
- проектным организациям при разработке прогнозных оценок
воздействия новых ЗШО на состояние воздушного бассейна и оценки
эффективности проектных технологических и природоохранных решений;
- контролирующим органам для оценки загрязняющих свойств объектов
ТЭС и эффективности природоохранной деятельности ТЭС, уточнения размера
платы за выбросы в атмосферу, экологической экспертизы создаваемых и
реконструируемых объектов. [3], [5].
2. Нормативные ссылки
1. Концепция экологической безопасности Республики Казахстан на
2004÷2015 годы от 3 декабря 2003 года № 1241.
2. Экологический кодекс Республики Казахстан (с изменениями и
дополнениями по состоянию на 17.07.2009 г.
3. Санитарно-эпидемиологические правила и нормы РК «Санитарноэпидемиологические требования к проектированию производственных
объектов», № 3792 от 17.08.2005.
4. СНиП РК 2.04-03-2002 «Строительная теплотехника» Астана, 2002.
3. Термины и определения
Адгезия – (от лат. adhaesio – прилипание) в физике – сцепление
поверхностей разнородных твёрдых и (или) жидких тел. Адгезия обусловлена
межмолекулярным взаимодействием в поверхностном слое и характеризуется
3
удельной работой, необходимой для разделения поверхностей.
Адгезия
существенно влияет на природу трения соприкасающихся поверхностей;
Батиметрическая съемка – процесс сбора данных о глубинах в зоне
исследования. Батиметрические данные используются для общего анализа
рельефа дна и картосоставления, решения прикладных задач (анализ
заполнения ЗШО, его заносимости, контроль дноуглубительных работ и др.);
Ветровая эрозия – разрушение почв и искусственных пород ветром,
включающее в себя отрыв и вынос обломков материала и сопровождающееся
их отложением. Ветровая эрозия разделяется на два типа – повседневная
ветровая эрозия и пыльные бури;
Вторичное пыление золоотвала – отрыв от слоя осаждаемых
эродированных частиц;
Гранулометрический (механический) состав складируемого материала –
относительное содержание в искусственной смеси частиц различных размеров
независимо от их химического или минералогического состава;
Диффузия – процесс переноса материи или энергии из области с высокой
концентрацией в область с низкой концентрацией;
Захоронение отходов – складирование отходов в местах, специально
установленных для их безопасного хранения в течение неограниченного срока;
Золошлаки – продукты сжигания твердого топлива;
Золошлакоотвал – место для складирования золы и шлака, образующихся
при сжигании твердого топлива на ТЭС;
Мониторинг окружающей среды – система регулярных наблюдений и
контроля, для оценки состояния окружающей среды, анализа происходящих в
ней процессов и своевременного выявления тенденций ее изменений;
Первичное пыление золоотвала – отрыв пылевых частиц от намытого
слоя;
Природопользование – осуществление пользования природными
ресурсами и (или) эмиссии в окружающую среду;
Экологическая экспертиза – установление соответствия хозяйственной и
иной деятельности нормативам качества окружающей среды и экологическим
требованиям, в целях предупреждения неблагоприятных воздействий этой
деятельности на окружающую среду;
Эродируемость золошлакоотвала – способность верхнего слоя ЗШО
разрушаться под воздействием ветра и атмосферных осадков, сопровождаемая
перемещением и переотложением продуктов разрушения [1], [2], [6], [7], [8].
4. Обозначения и сокращения
отн.ед.
мкм
ЗШО
РД
СНиП
ТЭС
относительные единицы
микрометр = 10-6м
Золошлакоотвал
Руководящий документ
Строительные нормы и правила
Тепловые электрические станции
4
5. Основные положения
5.1. Механизм и основные показатели ветровой эрозии
золошлакоотвала ТЭС
5.1.1. Гидрозолошлакоотвалы ТЭС - это хранилища отходов
энергетического производства – золы и шлака, для которых по мере заполнения
рекомендуется проведение мероприятий по консервации или рекультивации и
дальнейшая передача в хозяйственное пользование. В сухие периоды года
большинство ЗШО (кроме случаев складирования цементирующихся зол) источники интенсивного пыления в результате ветрового воздействия.
Ветровая эрозия золошлаковых полей – это разрушение сложившейся в
процессе складирования и статического хранения структуры золы уноса ТЭС
под воздействием ветровоздушного потока. Она определяется характером и
интенсивностью воздействия ветрового потока, свойствами складируемого
материала, технологией складирования, а также конструкцией объекта.
5.1.2. Рекомендуется выделение трех составляющих процесса ветровой
эрозии:
- отрыв и взлет частицы с поверхности;
- перемещение ее в запыленном потоке над поверхностью золошлакового
поля;
- рассеивание эродированных частиц золы за пределами золоотвала после
схода пылевого облака с дамбы.
5.1.3. Воздействие ветровоздушного потока на каждую отдельную
частицу на поверхности слоя связано с несколькими одновременно
действующими механизмами, наиболее важные из которых:
- лобовое аэродинамическое давление, побуждающее к сдвигу по
направлению ветра вдоль поверхности;
- перепад статического давления, возникающий при обтекании частицы и
создающий подъемное усилие;
- турбулентная диффузия в ветровом потоке, создающая переменные,
пульсирующие по величине и направлению усилия на частицу и ослабляющие
гравитационные и адгезионные связи частицы со слоем.
5.1.4. Сдвиговая или динамическая скорость потока U* – важнейший
параметр, характеризующий условия ветровой эрозии поверхности
пылевидного материала и определяющий величину создаваемого потоком
усилия отрыва частицы от поверхности (либо сдвига ее по поверхности).
Непосредственно этим параметром характеризуется степень интенсивности
турбулентного пульсационного движения в пограничном слое. Значение U*
определяется из логарифмического уравнения Кармана-Прандтля (1),
учитывающего высотный градиент скорости незапыленного ветровоздушного
потока, связанный с тормозящим действием подстилающей поверхности [8]:
Uz 1
Z
U  k  ln( d )  8,5
*
П
(1)
5
где: UZ – скорость ветра на высоте флюгера Z, м/с;
Z
– высота установки флюгера, м;
k
– постоянная Кармана (~0,4);
dП – диаметр частицы, м;
U* – динамическая скорость потока
Значения U* для золовых частиц с агрегатной плотностью ρа = 2,6 г/см3
показаны на рисунке 1.
Рисунок 1. Динамическая скорость ветровоздушного потока при различной
скорости ветра на уровне флюгера (Z = 10м).
Минимальное значение U*t, (м/с), соответствующая началу пыления,
определяется по графику (рисунок 2) или по формуле (2):
U  A  (  g  d П ) 0,5
*t
(2)
где: δ = ρП / ρВ
ρВ - плотность воздуха, кг/м3;
ρП - агрегатная плотность пылевых частиц, кг/м3;
g
- ускорение силы тяжести, м/с
dП - диаметр частицы, м;
А - эмпирический коэффициент, принимаемый равным 0,08 ÷ 0,12
Для
каждого
отдельного
случая
складирования
золошлаков
рекомендуется минимальная скорость ветра UКР на уровне флюгера,
соответствующая пороговой динамической скорости U*t, при которой
начинается пыление поверхности золоотвала. Определение этой величины
рекомендуется расчетным путем, исходя из зависимости (1) по формуле (3),
которой соответствует график на рисунке 1, либо наблюдениями по данному
объекту [5]; [7]; [8].
6
Рисунок 2. Значение пороговой скорости U*t при А=0,1
 z
U KP  [2,5  ln 
 dП

  8,5]  U *t

(3)
где: Z
– высота установки флюгера, м;
dП – диаметр частицы, м;
U*t – пороговая (минимальная) динамическая скорость
5.1.5. Поведение золошлаковых частиц в слое при продольном ветровом
воздействии характеризуется физическими свойствами материала и
индивидуальными свойствами каждой частицы (размер, плотность,
конфигурация, парусность), условиями ее физико-химического воздействия в
слое со смежными частицами (адгезионные свойства, химическое воздействие),
а также структурой слоя частиц (влажность, гранулометрический состав) и
условиями его формирования (насыпной, намывной).
5.1.6. Влияние гранулометрического состава отложений на пылящие
свойства слоя экстремально. Наибольшая сдуваемость характерна для слоя с
преобладанием частиц размером 0,05 ÷ 0,2 мм. Малопылящий характер
свойственно поверхности золошлаковых отложений, образованная частицами
dП < 0,05 мм, что обусловлено взаимодействием между частицами при их более
плотной упаковке. С увеличением доли частиц размером более 0,25 ÷ 0,3 мм (за
счет фракционирования при намыве слоя, либо по мере эрозионного истощения
поверхности сухого участка) интенсивность сдува золошлаков резко снижается.
5.1.7. Определяющим фактором химического взаимодействия золовых
частиц с образованием различной степени устойчивости сцементированных
конгломератов признается содержание оксида кальция CaO.
При содержании CaO < 10 % в золоотвале при намыве и хранении нет
структурной перестройки зернового состава поступающих золошлаков. При
содержании CaO = 10 ÷ 25 % существует частичная структурная перестройка,
приводящая к образованию на золошлаковом поле локальных участков слабо
сцементированных отложений. Золошлаки с содержанием CaO > 25 %
самоцементируются в монолит, стойкий к атмосферным воздействиям.
5.1.8. Интенсивность ветровой эрозии поверхности зависима от способа
формирования слоя – для намытого слоя золошлакоотвала она на порядок
ниже, чем для насыпного и неуплотненного слоя. Эти сведения свидетельство
7
существенного различия условий пыления сухих свеженамытых участков и зон
вторичного пыления, поверхность которых образована осажденными
эродированными частицами.
5.1.9. Капиллярное увлажнение поверхности слоя, применение орошения
в качестве способа пылеподавления, осадки, талые воды, изменяя влажность
материала, существенно факторы, оказывающие влияние на пыление
поверхности золоотвала. Пыление характерно для участков, где влажность
поверхности слоя составляет 1 ÷ 3 % (воздушно-сухой материал), и полностью
отсутствует при влажности более 6 ÷ 7 %.
5.1.10. Интегральная характеристика эродируемости частиц представлена
удельной сдуваемостью материала mO, определяемой экспериментальным
путем продувками в аэродинамической трубе проб золы, отобранных на
золоотвале, с моделированием условий намыва золы. Для приближенной
оценки порядка величины mO (г/м2с) рекомендуется эмпирическая зависимость:
mO  c  U* (U*2  U*2t )
(4)
где: с – универсальная постоянная (с = 100)
U* – динамическая скорость потока
U*t – пороговая (минимальная) динамическая скорость
5.1.11. Известны три формы движения эродированных пылевых частиц
под действием воздушного потока:
- сальтация (скачкообразное перемещение);
- во взвешенном состоянии;
- безотрывное перемещение частиц по поверхности (перекатывание).
Зола, сдуваемая с поверхности золоотвала и выносимая за его пределы, –
полидисперсная с пылевыми частицами субмикронных размеров до 500 мкм
(витающие – размером до 40 мкм и гравитирующие, которые участвуют в
сальтирующем движении, – размером до 500 мкм). Крупные частицы размером
свыше 500 мкм перемещаются по поверхности под действием лобового усилия
воздушного потока безотрывно. Эти частицы практически не выносятся за
пределы золошлакового поля. Витающие частицы увлекаются ветровым
потоком и рассеиваются на значительном удалении от золошлакоотвала.
При сальтации частица, получая начальный подъемный пульс от
воздушного потока, перемещается вниз под воздействием силы тяжести и силы
трения о воздух. Критерием разделения процессов сальтации и переноса взвеси
служит отношение динамической скорости U*t к скорости гравитационного
оседания частиц Vt . Граничный размер пылевой частицы dГР (минимальный
размер сальтирующей частицы), разделяющий витающие и сальтирующие
частицы, соответствует значению комплекса
U *t
 1,0
Vt
(5)
8
Рекомендуется определение значения граничного размера пылевой
частицы dГР по графику на рисунке 3.
Рисунок 3. Значение граничного размера частиц dГР при различной агрегатной
плотности золы.
5.1.12. Каждому значению скорости ветра на высоте флюгера Uz
соответствует свой предельный (максимальный) размер эродируемой частицы
dmax , определяемый из соотношения U*t / Ut по графику на Рисунке 4.
Рисунок 4. Значение предельного размера эродируемых золовых частиц при
различной скорости ветра на уровне флюгера для различной плотности
пылевых частиц ρП: 1 - ρП = 3,0 г/см3; 2- ρП = 2,65 г/см3; 3- ρП = 2,0 г/см3.
5.1.13. Определение высоты подъема эродированной частицы h (без учета
воздействия восходящих воздушных потоков и турбулентных вихревых
образований) рекомендуется по формуле:
2
U ЭФ
1  1
h
  [1    ln(1  a)  1] , м
g а  а
a  0,0383 
U *  U ЭФ
  dП
где: Uэф - скорость ветра на уровне оси пылевого облака, равная 0,8 Uz
g
- ускорение силы тяжести, м/с
dП - диаметр частицы, м;
(6)
9
5.1.14. Определение пролета (м) сальтирующей
золошлаковым полем рекомендуется зависимостью:
LП 
h
 U ЭФ ,
Vt
частицы
над
(7)
где: h
- высота подъема эродированной частицы
Vt - скорость гравитационного оседания частиц
Uэф - скорость ветра на уровне оси пылевого облака, равной 0,8 Uz
5.1.15. На длине перелета одиночной пылевой частицы от наветренной
границы пылящего участка Lпр, происходит нарастание массы эродированных
частиц, переносимых ветровоздушным потоком, где отсутствует осаждение
оторванных от поверхности частиц.
На последующем участке процесс выдувания частиц из слоя
характеризуется переходом эродированных сальтирующих частиц из потока на
поверхность, а так же участием во вторичном пылении. Условия отрыва от слоя
осаждаемых эродированных частиц (вторичного пыления) существенно
отличны от характера отрыва пылевых частиц от намытого слоя (первичного
пыления). Интенсивность сдува частиц при вторичном пылении существенно
выше. Поэтому, весь эродированный материал первичного сдува, включая
витающие и сальтирующие частицы (со всей пылящей неэкранированной
прудом поверхности золоотвала) выносится к границе золошлакового поля (к
ограждающей дамбе отвала). Участок непосредственно перед ограждающей
дамбой, сама дамба, ее низовой откос, дренажные и нагорные канавы - это зоны
обеспыливания ветрового потока. В зоне обеспыливания ветрового потока
действуют гравитационные силы и турбулентная диффузия.
5.1.16. Перенос золовых частиц, поступающих в атмосферу с открытой
поверхности складируемого материала, и переход их на подстилающую
поверхность в прилежащей зоне - это механизм, отличный от процессов
рассеивания в атмосфере твердых частиц с дымовыми газами ТЭС и
твердофазных выпадений из дымового факела.
Обеспыливание ветровоздушного потока (мг/м3) - это эмпирическая
экспоненциальная зависимость (8):
 X  O  c  aX ,
(8)
где: μ0 – запыленность на границе золоотвала, мг/м3
μх – запыленность на расстоянии Х от дамбы, мг/м3
α – коэффициент затухания, 1/м; средний коэффициент затухания по
данным разных съемок составляет αср = 6,2·103
Угол бокового раскрытия факела пыли после схода с дамбы по мере его
распространения равен 10° [5]; [7]; [8].
5.1.17. При определении максимального текущего выноса золовых частиц
с золоотвала и максимальных приземных концентраций рекомендуется
10
выполнение расчета по максимальной скорости ветра на уровне флюгера с
повторяемостью не менее 5%.
5.2. Характеристика золоотвала как площадного пылящего объекта
5.2.1. Для золошлакоотвала не характерна сплошная пылящая
поверхность. Выделим три основных возможных зон пыления:
- сухие пляжи на золошлаковых полях;
- дамбы, сложенные из золошлаков (неприсыпанных);
- отложения пыли в аэродинамической тени дамбы (вторичное
пыление).
Поверхность золошлакового поля в значительной мере характеризуется
способом намыва. Сбросом золошлаковой пульпы из распределительного
пульпопровода по периметру отвала обеспечивается фракционирование
материала – осаждение наиболее крупных фракций у внутреннего откоса
дамбы. Наличие отстойного пруда намывного золоотвала - это условие полного
осаждения наиболее мелких частиц золы. Неоднородность фракционного
состава поверхностного слоя золошлакового поля - это следствие
неравномерного распределения золошлакового материала по фронту намыва:
наиболее крупные фракции осаждаются у первых по ходу пульпы выпусков.
Рекомендуются способы направленного формирования неоднородности
фракционного состава осаждаемого материала – в частности, раздельный
способ намыва, при котором наиболее крупные и тяжелые фракции (в основном
шламовые частицы) отводятся из нижней точки распределительного
пульпопровода и сбрасываются в периферийной части золоотвала, оставшаяся
часть пульпы отводится к центральной части отвала.
В зоне надводных отложений золошлакоотвала увлажнение намытого
материала рекомендуется проводить за счет капиллярного подъема воды,
периодического изменения уровня пруда-осветлителя, осадков, применения
специальных дождевальных систем поверхностного увлажнения.
Нарушение поступления влаги к поверхности золошлакового поля
активно способствует образованию сухих пляжей.
Сухой пляж золоотвала - это сложное сочетание участков интенсивного,
среднего и малого пыления. Интенсивное пыление характерно на участках с
толстым слоем золы песчаной (несвязанной) фракции. Степень однородности
намываемых золошлаков по дисперсному составу - существенный фактор в
дифференциации золошлакового материала при намыве и формировании
пылящих сухих пляжей. Золошлаки высокой степени однородности
(характерной для экибастузских углей) слабо дифференцируются при намыве в
золоотвал - причина опасности пыления по большей части поверхности
золошлакового поля. При высоком содержании шлаковых фракций (угли
подмосковного бассейна, сжигание углей в топках с жидким шлакоудалением)
на пляжах намыва осаждаются крупные фракции (более 0,5 мм в диаметре),
устойчивые к ветровой эрозии. Тонкодисперсная зола складируется в зоне
отстойного пруда и капиллярного смачивания [5].
11
5.2.2. Эродируемость поверхности слоя пылевых частиц зависима от
условий формирования слоя и эродирующих свойств частиц. Наиболее важные
параметры поверхности, подвергаемой ветровому воздействию:
- гранулометрический состав складируемого материала;
- влажность поверхностного слоя пылевых частиц – технологическая (за
счет увлажнения путем орошения, остаточная влажность при колебании уровня
воды в золоотвалах, увлажнение поверхностного слоя капиллярной влагой из
подслойного водяного объема) и вследствие атмосферных осадков;
- плотность поверхностного слоя, определяемая способом намыва золы
(подводный или надводный, рассредоточенный или сосредоточенный),
использованием механических средств уплотнения поверхностного слоя;
- снеговой покров;
- растительный покров на отработанных ЗШО [5]; [7].
5.2.3. Эродируемость объекта рекомендуется определять с учетом
следующих дополнительных факторов (помимо рассмотренных в пунктах 5.2.1;
5.2.2 свойств материала и состояния поверхностного слоя):
- защищенности объекта от ветрового воздействия рельефом
прилегающей местности;
- конструкции объекта – наличия и высоты ограждающих систем (дамбы,
защитные стенки, лесопосадки по периметру золоотвала);
- эксплуатационных характеристик – уровня поверхности слоя
относительно дамбы, окружающей территории, площади сухих пляжей в
гидрозолоотвалах.
Пруд-осветлитель - это накопитель наиболее мелких фракций золы. При
снижении уровня воды относительно расчетного осушенные отложения возможные источники интенсивного пыления. Пруд практически полностью
поглощает гравитирующие эродированные и, частично, витающие частицы,
поступающие с ветровым потоком на акваторию.
При обтекании ветровоздушным потоком дамбы происходит частичное
его обеспыливание – выпадение значительной части эродированных
золошлаковых частиц, преимущественно сальтирующих.
Рельеф примыкающей местности и высотное профилирование
золошлакоотвала - существенные факторы влияния на характер и
интенсивность ветровых потоков вдоль поверхности золошлакового поля. В
ходе эксплуатации золоотвала по мере заполнения его и наращивания
ограждающих сооружений высотная характеристика объекта по отношению к
примыкающей местности существенно трансформируется [5]; [7]; [8].
5.3. Расчетная оценка ветровой эрозии золовых шлаков ТЭС
5.3.1. При оценке пылевых выбросов золоотвалов ТЭС рекомендуется
принятие следующих исходных положений:
- учитывая вынос всех эродируемых витающих частиц ветровоздушным
потоком за пределы золошлакоотвала, осаждение их на поверхности прудаосветлителя и дамбы маловероятно;
12
- сальтирующие частицы полностью экранируются отстойным прудом,
при оценке пылящих свойств объектов средний размер частицы dП
рекомендуется принимать по наиболее характерному (средневзвешенному)
размеру эродируемых частиц в диапазоне от dmax до 0;
- при оценке загрязнения воздушного бассейна пылящим
золошлакоотвалом рекомендуется расчет как для линейного источника,
располагающегося на границе отвала с подветренной стороны, при этом
высотный размер данного источника равно высоте пылящего облака;
- однородность гранулометрического состава материала на пылящей
поверхности - одно из условий расчета;
- оценка пыления рекомендуется для нормальных, характерных для
данного региона метеорологических условий, не учитывая возможность
аномальных атмосферных явлений (смерч, ураганные ветры), которые могут
привести к нерасчетному катастрофическому сдуванию складируемого
мелкозернистого и пылевидного материала;
- для существующих ЗШО рекомендуются следующие параметры:
площадь и положение пылящих участков золошлакового поля;
гранулометрический состав на пылящих участках; период пыления (скорость
ветра, соответствующая началу пыления золоотвала, продолжительность
периода, в течение которого возможно пыление); высота пылевого облака над
поверхностью золоотвала и на гребне дамбы;
Для проектируемых ЗШО рекомендуются:
- площадь пруда-осветлителя – по проектными материалами (до 70 ÷ 80 %
площади золошлакоотвала);
- площадь пылящего участка – до 20 % поверхности сухих пляжей;
- сосредоточение пылящих участков у дамбы с подветренной стороны;
- средняя влажность пылящих частиц WП = 2,5 %;
- средняя скорость ветровоздушного потока, осуществляющего перенос
эродируемых частиц над поверхностью золошлакоотвала (на оси пылевого
облака) - 0,8Uср ;
- высота пылевого облака на гребне дамбы - 2h, т.е. удвоенная высота
подъема эродированных частиц над золошлаковым полем.
При отсутствии отдельных характеристик ветровой эрозии для
существующих ЗШО рекомендуется принятие значений соответствующих
характеристик для проектируемых объектов.
5.3.2. Учет конструктивных, планировочных и природно-климатических
факторов рекомендуется введением ряда поправочных коэффициентов K 1÷ K4
к расчетной величине сдува:
K1 - обеспыливание пылевого потока за счет осаждения золовых частиц
при обтекании дамбы. Значение K1 рекомендуется в зависимости от
превышения гребня дамбы относительно уровня поверхности золошлакового
поля по графику на рисунке 5;
K2 - состояние поверхностного слоя (коркообразование, агрегатирование
золовых частиц в слое в результате химического взаимодействия) в
зависимости от содержания окиси кальция в золе:
K2 = 1,0 при CaO < 10 %;
13
K2 = 1 - 1,6129∙10-3(СаО)2 при CaO = 10 ÷ 25 %; K2 = 0 при CaO > 25 %;
K3 - защищенность объекта от ветрового воздействия (влияние высотных
элементов рельефа, специальных ветрозащитных сооружений, лесопосадок) и
закрепления поверхности зольного пляжа. Выбор значения коэффициента K3
рекомендуется по таблице 1 Приложения 1 к настоящей Методике. При
одновременном воздействии нескольких факторов защищенности величина
коэффициент K3 - это произведение соответствующих коэффициентов.
K4 - применение оперативных методов пылеподавления (орошение
пылящей поверхности водой и др.). Выбор значения коэффициента K4
рекомендуется по таблице 2 Приложения 1 к настоящей Методике.
Рисунок 5. Зависимость коэффициента переноса золовых частиц при обтекании
дамбы K1 от скорости ветра Uz при различных превышениях гребня дамбы над
зольным пляжем h.
5.3.3. Гранулометрический состав поверхности пылящих участков отвала
рекомендуется определять экспериментальным путем для данного конкретного
складируемого материала при характерной влажности, либо по справочным
данным для пойменной золошлаковой зоны отвала [9].
5.3.4. При оценке ветровой эрозии (т/год) и среднегодового текущего
выноса частиц (г/с) с поверхности золоотвала рекомендуется определение
отдельных характеристик эродируемого материала и ветрового режима
следующим образом:
- средняя скорость ветра на уровне флюгера UСР в зоне размещения
золоотвала принимается как средневзвешенное значение в диапазоне от
скорости ветра (начало сдува золовых частиц UКР) до максимальной скорости
ветра UМАХ с учетом повторяемости градаций скоростей;
- средний размер эродируемых частиц dСР принимается как
средневзвешенная величина в диапазоне от dМАХ до 0:
d CP 
 (d  a)i
d
i
где: dì – средний размер частиц i-той фракции
(9)
14
a – доля частиц i-той фракции;
i – число фракций эродируемых частиц
Максимальный размер эродируемых частиц dMAX рекомендуется
определять по средней скорости ветра в пылеопасный период U’CP,
первоначальное значение которой рекомендуется равным 5 ÷ 7 м/с, далее
рекомендуется последовательное определение: d MAX  d CP  U *t  U KP  U CP
При расхождении значений U’CP и UCP более чем на 0,2 м/с
рекомендуется повторный расчет, после корректировки значения U’CP .
5.3.5. Расчет продолжительности периода возможной ветровой эрозии
поверхности золоотвала τП рекомендуется, исходя из двух основных временных
характеристик:
- относительной продолжительности пылеопасного ветрового режима
τ’Uкр, когда скорость ветра на флюгере Uz больше скорости начала пыления UКP
(определяемой по средневзвешенному размеру эродируемых частиц);
- относительной продолжительности периода возможного пыления τ’ПОВ
по состоянию поверхности золоотвала, исключающего из рассматриваемого
периода τ продолжительность периодов укрытия золовых пляжей устойчивым
снеговым покровом, увлажнения осадками и талыми водами:
'
'
'
'
 ПОВ
 1  ( CH
  OC
  T'   ШТ
)
(10)
где: τ’СН - относительная продолжительность устойчивого снегового покрова;
τ’ОС - относительная продолжительность осадков в виде дождя и мокрого
снега;
τ’Т - относительная продолжительность увлажнения поверхности золоотвала
талыми водами;
τ’ШТ - относительная продолжительность штиля
В целом определение продолжительности периода (ч) возможной
ветровой эрозии (пыления) поверхности золоотвала рекомендуется наложением
двух рассмотренных временных характеристик:
'
'
 ПЫЛ     ПОВ
  Ukp
или
'
'
'
'
 ПЫЛ    [1  ( СН
  OC
  T'   ШТ
)]   Ukp
,
(11)
где: τ - продолжительность рассматриваемого периода, ч
5.3.6. Текущий вынос золовых частиц с поверхности золоотвала (г/с).
При скорости ветра выше критической UКР (ниже которой m0 =0).
М ВЫН  М ВЫН ВИТ  М ВЫН САЛЬТ ,
(12)
где: МВЫН ВИТ - масса выносимых за границы золошлакоотвала витающих
частиц, г/с
МВЫН САЛЬТ - масса выносимых за границы золошлакоотвала сальтирующих
15
частиц, г/с
М ВЫН  тО  (nВИТ  S  nСАЛЬТ  SЭФ  К1 )  К 2  К3  К 4
(13)
где: mO - удельная сдуваемость материала пылящей поверхности при данном
значении скорости потока на высоте флюгера, г/м2с, определяемая
экспериментально продувкой проб золы с пылящих участков золоотвала в
аэродинамической трубе с моделированием условий намыва золы и с
приведением к скорости потока на высоте флюгера в соответствии с
приложением 2 к настоящей Методике, либо по формуле (4);
nвит ,nсальт - доля витающих и сальтирующих частиц в сдуваемой золе;
Sэф - эффективная площадь пылящей поверхности (м2), на которой
завершается нарастание в потоке массы сальтирующих частиц:
SЭФ  S  SЭК
(14)
где: S - полная площадь пылящей поверхности золоотвала, м2;
Sэк - часть пылящего участка золоотвала, экранируемая отстойным прудом
(находящаяся с наветренной стороны относительно пруда), м2
В приложении 3 к настоящей Методике на рисунке 6 приводится
удельная сдуваемость золовых отложений на отвалах ГРЭС использующих
уголь различных месторождений.
5.3.7. Годовой вынос золовых частиц или эродируемость объекта (т/год)
рекомендуется определением суммы выноса золовых частиц по каждому
направлению ветра
М ЭР   М ЭРi
(15)
где: MЭРi - годовой вынос золовых частиц по каждому направлению ветра:
М ЭРi  mOCP  (nВИТ  S  пСАЛЬТ  SЭФ  К1 )  К 2  К3  К 4   ПЫЛ  3,6  10 6 , т/год
(16)
где: тСР0 - удельная сдуваемость в г/м с, соответствующая средней скорости
ветра UСР;
К1;К2;К3;К4 - поправочные коэффициенты (пункты 5.3.2)
τПЫЛ - продолжительность периода пыления, ч.
16
Приложение 1 к «Методике расчета
оценки ветровой эрозии и пыления
золоотвалов тепловых
электростанций»
Значения поправочных коэффициентов К3 и К4
Таблица 1 - Значения поправочного коэффициента К3, учитывающего
защищенность золошлакоотвала от ветрового воздействия (влияние высотных
элементов рельефа, специальных ветрозащитных сооружений, лесопосадок) и
закрепления поверхности зольного пляжа
Факторы защищенности отвала от пыления
1. Закрытие отвала высотными элементами рельефа:
- с одной стороны
- с двух сторон
- с трех сторон
2. Сооружение сплошных барьеров по периметру дамб отвала
(решетчатые ограды, прокладка пульпопроводов по гребню дамбы,
лесополосы вдоль границ отвала)
3. Относительное увеличение высоты дамбы, расположенной с
подветренной стороны по отношению к преобладающим ветрам
4. Закрепление поверхности зольного пляжа вяжущими веществами
(коркообразование)
5. Закрепление поверхности зольного пляжа шлаком
6. Закрепление поверхности зольного пляжа защитным слоем из
связного грунта (суглинок, глина)
К3
0,6
0,3
0,15
0,7
0,85
0,1
0,05
0,02
Примечание: При одновременном воздействии нескольких факторов
защищенности коэффициент K3 определяется перемножением соответствующих
коэффициентов.
Таблица 2 - Значения поправочного коэффициента К4, учитывающего
защищенность золошлакоотвала при применении оперативных методов
пылеподавления
Метод оперативного пылеподавления
1. Поднятие уровня воды в пруде-осветлителе выше уровня золового
пляжа
2. Периодическое орошение сухих пляжей стационарными
дождевальными установками или поливочными машинами
3. Смачивание сухих пляжей осветленной водой, подаваемой по
резервному пульпопроводу разводящей сети
4. Учащение переключения пульповыпусков в теплое время
К4
0÷0,2
0,1÷0,5
0,2÷0,3
0,7
17
Приложение 2 к «Методике расчета
оценки ветровой эрозии и пыления
золоотвалов тепловых
электростанций»
Приведение экспериментально определенной удельной сдуваемости
к скорости ветра на высоте флюгера
Для оценки эродируемости частиц складируемого материала
используется устанавливаемый экспериментальным путем интегральный
параметр - удельная сдуваемость материала m0, которая является базовой
характеристикой для расчета ветровой эрозии поверхности золошлакоотвала и
определяется продувкой проб материала в аэродинамической трубе.
В аэродинамической трубе (так же, как в натурных условиях) над
обдуваемой поверхностью золошлакового материала формируется высотный
градиент скорости. В закрытом канале, используемом для экспериментального
определения удельной сдуваемости золошлаковых частиц, динамическая
скорость потока определяется через величину сдвигового напряжения на стенке
канала:
U*  (
 W 0,5
)
В
 W   B  u 2  0,045 Re0, 25
(17)
(18)
где: τW - величина сдвигового напряжения на стенке канала, определяемая по
формуле (18);
ρВ - плотность воздуха, кг/м3;
u - скорость на оси канала, м/с;
Re - число Рейнольдса, определяемое по формуле (19);
R - гидравлический радиус канала, определяемый по формуле (20), м;
a; b - стороны прямоугольного канала, м;
ν - кинематическая вязкость воздуха, м2/с
uR
(19)
v
a b
(20)
R
0,5  (a  b)
Располагая расчетной зависимостью U* = f(UZ) (1) для натурного
золошлакоотвала и экспериментальным графиком функции m0 = f(U*),
устанавливается значение удельной сдуваемости для материала, складируемого
на золошлакоотвале, при любом значении скорости ветра, на любой высоте.
Re 
18
Приложение 3 к «Методике расчета
оценки ветровой эрозии и пыления
золоотвалов тепловых
электростанций»
Удельная сдуваемость золовых отложений на отвалах Рефтинской и
Южноуральской ГРЭС
Скорость ветра U*t , м/с
Рисунок 6. Зависимость удельной сдуваемости намытой золы mO от скорости
ветра на высоте флюгера U*t (10 м)
19
Приложение 4 к «Методике расчета
оценки ветровой эрозии и пыления
золоотвалов тепловых
электростанций»
Алгоритм и пример расчетного определения ветровой эрозии
золошлакоотвала
Таблица 1 - Пример расчетного определения ветровой эрозии золошлакоотвала.
Наименование параметра
Расчетная формула, источник
Величина параметра в
примере расчета
1. Основные исходные данные
Для
Для проектируемых Проектируемый золоотвал
действующих
отвалов
ТЭС
отвалов
1.1. Содержание CaO в
фондовые
справочные
Содержание CaO < 10 %
золошлаках
материалы ТЭС
материалы
> 10 мм – 3,4;
10 ÷ 5 мм – 4,4;
1.2. Гранулометри- ческий
5 ÷ 2 мм – 1,65;
состав (в %) золошлаков
2 ÷ 1 мм – 3,4;
на поверхности
данные
справочные
1 ÷ 0,5 мм – 4,35;
золошлаковой зоны сухого лабораторного материалы по виду
0,5 ÷ 0,25 мм – 18,5;
пляжа, агрегатная
анализа
сжигаемого топлива 0,25 ÷ 0,125 мм – 25,8;
плотность
0,125 ÷ 0,063 мм – 21,05;
ρП, кг/м3
0,063 ÷ 0,04 мм – 8,9;
< 0,04мм – 8,55.
ρП =2400 кг/м3
1.3. Характеристика
данные метеонаблюдений
ветрового режима:
0 ÷ 1 м/с – 38,05;
2 ÷ 3 м/с – 33,2;
1.3.1. Повторяемость в
4 ÷ 5 м/с – 17,5;
течение года скорости
данные метеонаблюдений
6 ÷ 7 м/с – 7,76;
ветра различных градаций
8 ÷ 9 м/с – 2,6;
РUi , %
10 ÷ 11 м/c – 0,54;
> 11м/с – 0,35
1.3.2. Повторяемость
С – 20,8; СВ – 12;
различных направлений
В – 7,2;
данные метеонаблюдений или по
ветра в течение года (по
ЮВ – 5,8;
данным [4]
восьмирумбовой розе
Ю – 12, ЮЗ – 10,2;
ветров) Pi, %
З – 12; СЗ – 20;
1.4. Относительная
продолжительность
периодов ограничения
пыления по
данные метеонаблюдений
56,0
состоянию поверхности в
% годового фонда
времени, в том числе:
20
Продолжение таблицы 1
Наименование параметра
Расчетная формула, источник
Величина параметра в
примере расчета
1.4.1.Устойчивый
снеговой покров τСН
1.4.2. Увлажнение талыми
водами τТ
данные метеонаблюдений или по
данным [4]
42,5
данные метеонаблюдений
2,7
1.4.3. Осадки τОС
1.5. Площадь золового
поля SОТВ , м2
1.6. Площадь прудаосветлителя SПР , м2
1.7. Площадь пылящих
участков, S, м2
2.1 .Предельный размер
эродируемых частиц
dmax , мкм
2.2. Средневзвешенный
размер dПСР (мкм)
эродируемых частиц в
пылящем слое (при dП <
dmax)
2.3. Граничный размер
эродируемых частиц,
разделяющий
сальтирующие и
витающие частицы dГР ,
мкм
данные метеонаблюдений или по
10,8
данным [4]
Батиметричеспроект
6525000
кая съемка
БатиметричесSПР =(0,6÷ 0,9)·
2680000
кая съемка
SОТВ
БатиметричесS=0,2·( SОТВ – SПР )
769000
кая съемка
2. Характеристика эродируемых частиц
график на рисунке 4 или по формуле:
U CP
 A  (  g  d max ) 0,5
Z
[2,5  ln(
)  8,5]
при U‫׳‬CP = 6,2 м/с; dmax =
d max
160
Предварительно принимается
значение средней скорости ветра
пылеопасного периода U‫׳‬CP
 (d П а) i
d ПСР 
 аi
i - количество градаций размера
0,098
частиц;
а - весовая доля соответствующей
градации
график на рисунке 3
dГР =31,0
График на рисунке 2 или по формуле:
2.4. Пороговая
динамическая скорость
ветрового потока для
средневзвешенного
размера эродируемых
частиц в слое U*t , м/с
2.5. Доля витающих и
сальтирующих частиц в
общей массе
эродируемого золового
материала (dП < dmax )
U *t  0,1  (  g  d ПСР ) 0,5
где: g- ускорение силы тяжести, м/c2;

П
В
U*t = 0,135
ρП - плотность частиц, кг/м3,
ρВ - плотность воздуха, кг/м3
 а d max dгр
п ВИТ 
1  а  d max
nСАЛ 
а
dгр 0
1  а  d max
nСАЛ = 0,87
nВИТ = 0,13
21
Продолжение таблицы 1
Наименование параметра
Расчетная формула, источник
Величина параметра в
примере расчета
3. Характеристика ветрового режима
3.1. Скорость ветра на График на рисунке 1 или по формуле:
уровне флюгера UКР,
10
U KP  [2,5  ln(
)  8,5]
соответствующая U*t
z
0
,
098
U KP  [2,5  ln( СР )  8,5]  U *t
(скорости начала пыления
dП
 0,166
поверхностного слоя), м/с
 (U  P) i ;
U CP 
 Pi
3.2. Среднегодовая
скорость ветра на уровне где i - количество градаций скорости
UCP =6,2 = U’CP , пересчет
в диапазоне от UКР до Umax ,
флюгера в пылеопасном
не требуется
ветровом режиме (U > U- средняя скорость ветра в пределах
градации
UКР), UCP , м/с
Р- повторяемость различных
направлений ветра, %
3.3. Относительная
продолжительность
ветрового режима (U >
Данные метеонаблюдений
τ′Uп = 0,08
UКР), способствующего
выдуванию золовых
частиц τ′Uп
6,2
График на Рисунке 1 или по формуле: U * 


10
3.4. Динамическая
U CP
[2,5  ln 
  8,5]
U


3
*
скорость потока U* ,
 0,098  10 
 z 
[2,5  ln  СР   8,5]
соответствующая UCP , м/с
 0,166
 dП 
4. Характеристика золоотвала
4.1. Эффективная площадь
пылящей поверхности SЭФ
, м2
4.2. Продолжительность периода
возможного пыления
Р
 ПЫЛ
τРПЫЛ в течение
рассматриваемого периода
τ по каждому
направлению, ч
4.3. Поправочные
коэффициенты к величине
выноса золовых частиц:
К1
SСЭФ = 619000;
SСВЭФ = 425000;
SСЗЭФ = 653000;
SЮЭФ = 469000;
SЭФ= S – SЭК
SЮВЭФ = 260000;
SЮЗЭФ = 569000;
SЗЭФ = 769000;
SВЭФ = 635000
При τ = 8760 часов:
τСПЫЛ = 30,7; τСВПЫЛ =
17,73;
'
'
   [1  ( CH
  OC
  T'   ШТ )]   U'  P
τСЗПЫЛ = 29,6;
τЮПЫЛ = 17,73; τЮВПЫЛ =
8,57; τЮЗПЫЛ = 15,07;
τЗПЫЛ = 17,73; τВПЫЛ =
10,64.
По графику на Рисунке 5
К1 =0,775 (превышение
График на рисунке 5
гребня дамбы над
пляжем составляет 1 м)
22
Продолжение таблицы 1
Наименование параметра
К2
Расчетная формула, источник
СаО < 10 %
СаО = 10 ÷ 25 %
К2 = 1,0;
К2 = 1 – 103(0,127
СаО)2;
К2 = 0
СаО > 25 %
К3
Таблица 1
К4
Таблица 2
5. Параметры ветровой эрозии
5.1. Удельная сдуваемость Данные экспериментальной продувки
(приложение 2)
золовых частиц, mO,
Приближенная оценка по формуле:
соответствующая
2
mo  100  U *2  (U *2  U *2t )
величине UCP , г/м с
5.2. Эродируемость
золоотвала по каждому из
направлений ветра, МЭРi ,
т/год
M i  mO  (nВИТ  S  nСАЛ  S ЭФ  К1 ) 

К 24  ПЫЛ
3600
Величина параметра в
примере расчета
К2 = 1,0
К3 = 0,15 (отвал укрыт
от ветра холмами с трех
сторон)
К4 = 0,7 (увлажнение
пляжа переключением
пульпопроводов)
mO=100·0,1662·(0,1662–
0,1352)= 0,0257
МСЭР = 154,29;
МСВЭР = 66,58;
МСЗЭР = 155,35;
МЮЭР = 71,69;
МЮВЭР = 22,92;
МЮЗЭР = 70,8;
МЗЭР = 106,52;
МВЭР = 54,59.
5.3.Эродируемость
Σ(МЭР)i
МЭР =702,74
объекта МЭР , т/год
6. Расчет текущего пылевого выноса и рассеивания золовых частиц в атмосфере
6.1. Направление ветра
по заданию
юго-восточный
6.2.Скорость ветра Uz , м/с
по заданию
6,2
6.3. Предельный размер
при U’CP = 6,2 м/с dmax =
эродируемых частиц dmax,
график рисунок 4
160
мкм
 (d П  а) i
6.4.Средневзвешенный
d ПСР 
СР
размер d П (мкм)
 ai
dСРП = 0,098
эродируемых частиц в
i - количество градаций размера
пылящем слое
частиц;
(при dП < dmax)
a - весовая доля градации
6.5. Граничный размер
эродируемых частиц, dГР ,
мкм
6.6. Доля витающих и
сальтирующих частиц в
общей массе
эродируемого золового
материала (dП < dmax)
по графику рисунок 3
п ВИТ 
а
пСАЛ 
dГР = 31,0
d max  dгр
1  а  d max
а
dгр 0
1  а  d max
nВИТ = 0,13
nСАЛ = 0,87
23
Продолжение таблицы 1
Величина параметра в
примере расчета
График на Рисунке 1 или по формуле:
6,2
U* 


10
6.7. Динамическая
[2,5  ln 
  8,5]
U

3
z
скорость потока U*,
U* 
 0,098  10 
 z 
соответствующая UZ , м/с
 0,166
[2,5  ln  СР   8,5]
d
 П 
Наименование параметра
Расчетная формула, источник
График на рисунке 2 или по формуле:
6.8. Пороговая
динамическая скорость
ветра для
средневзвешенного
размера эродируемых
частиц в слое U*t м/с
U *t  0,1  (  g  d ПСР ) 0,5
где: g- ускорение силы тяжести, м/c2;

П
В
U*t =0,135
ρП - плотность частиц, кг/м3,
ρВ - плотность воздуха, кг/м3
6.9. Эффективная площадь
пылящей поверхности SЭФ
SЭФ =SЮВЭФ
SЮВЭФ =260000
, м2
6.10. Удельная
Данные экспериментальной продувки
сдуваемость золовых
(приложение 2)
mO =100·0,1662·(0,1662–
частиц mO ,
Приближенная оценка по формуле:
0,1352)= 0,0257
соответствующая
2
2
2
mo  100  U *  (U *  U *t )
величине UZ , г/м2с
Т
6.11.Текущий вынос
М ВЫН

Т
золовых частиц М ВЫН ,
МТВЫН = 742,8
тек
 тО  (пВИТ  S  nСАЛ  S ЭФ  К1 )  К 24
г/с
для действующих ЗШО – по
6.12.Ширина пылевого
наблюдениям; для проектируемых –
LП = 2460
облака LП , м
длина дамбы, с которой сходит
пылевое облако.
6.13. Эффективная
скорость ветра,
непосредственно
UЭФ = 0,8·UСР
UЭФ = 0,8·6,2 = 4,96
воздействующего на
эродированную частицу,
UЭФ , м/с
4,96 2
h(
)
9,8  0,173
1


2
 U ЭФ

1
(1  0,173 )  ln( 1  0,173)  1
  [(1  )  ln( 1  a)  1]
h  


6.14. Высота подъема
g  a 
a

h = 1,196
эродированных частиц над
золошлаковым полем h , м
а
U
а  0,0383  U *  ЭФ
4,96
 dП
0,0383  0,166 

1860  0,098 10 3
 0,173
6.15. Высота пылевого
h0=2 h
h0 = 2·1,196 = 2,392
облака на дамбе h0
24
Продолжение таблицы 1
Наименование параметра
6.16. Начальная
концентрация пылевых
частиц на сходе с дамбы
μ0, мг/м3
6.17. Приземная
концентрация пылевых
частиц μx (мг/м3) на
удалении x (м) от дамбы
Расчетная формула, источник
0 
М ВЫН
LП  h0  U ЭФ
Величина параметра в
примере расчета
742,8 10 3
0 
2460  2,392  4,96
 0  25,45
 х   0  е  ах ;
где a - коэффициент
рассеивания, равный 6,2·10-3
при x = 200 м μ = 7,36;
при x = 500 м μ = 1,15;
при x = 830 м μ = 0,148;
при x = 1000 м μ = 0,0516
Рисунок 7. Расчетная схема золоотвала, рассматриваемого в примере.
Примечание: Роза ветров пылеопасного периода принимается для каждого
объекта индивидуально по данным [4] или результатам местных наблюдений.
25
Библиография
1. Концепция экологической безопасности Республики Казахстан на
2004÷2015 годы от 3 декабря 2003 года №1241
2. Экологический кодекс Республики Казахстан (с изменениями и
дополнениями по состоянию на 17.07.2009 г.)
3. Санитарно-эпидемиологические правила и нормы РК «Санитарноэпидемиологические требования к проектированию производственных
объектов», № 3792 от 17.08.05
4. СНиП РК 2.04-03-2002 «Строительная теплотехника» Астана, 2002
5. РД 153-34.0-02.106-98 «Методика расчетной оценки ветровой эрозии и
пыления золошлакоотвала ТЭС», УралОГРЭС, 1999
6. Клименко А.В. Зорин В.М. «Теплоэнергетика и теплотехника»,
Справочная серия в 4 книгах М, Издательство МЭИ, 2000
7. Иванов Б. Н. «Мир физической гидродинамики: От проблем
турбулентности до физики космоса». Изд.2, М:., URSS, 2010
8. Биркгоф Г. «Гидродинамика». М.: Издательство иностранной
литературы, 1963
9. «Состав и свойства золы и шлака ТЭС», Справочное пособие под
редакцией В.А.Мелентьева, Л.,1985