автореферат диссертации по химической технологии, 05.17.08, диссертация на тему:Совершенствование техники и технологии процесса термической переработки древесных отходов

На правах рукописи

Грачев Андрей Николаевич

СОВЕРШЕНСТВОВАНИЕ ТЕХНИКИ И ТЕХНОЛОГИИ ПРОЦЕССА ТЕРМИЧЕСКОЙ ПЕРЕРАБОТКИ ДРЕВЕСНЫХ ОТХОДОВ

05.17.08- Процессы и аппараты химических технологий

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Иваново 2005

Работа выполнена в Казанском государственном технологическом университете на кафедре переработки древесных материалов.

Научный руководитель

Официальные оппоненты

Ведущая организация

-доктор технических наук, профессор, заслуженный изобретатель РФ Сафин Рушан Гареевич. -доктор технических наук, профессор Овчинников Лев Николаевич; -доктор химических наук Почивалов Константин Васильевич. -Государственный научно-

исследовательский институт

химических продуктов (г. Казань).

Защита диссертации состоится 14 марта 2005 года в_часов на заседании

диссертационного совета Д 212.063.05 при Ивановском государственном химико-технологическом университете по адресу: 153000, г. Иваново, пр. Ф. Энгельса, 7.

С диссертацией можно ознакомиться в библиотеке Ивановского государственного химико-технологического университета.

Автореферат разослан_февраля 2005 г.

Ученый секретарь диссертационного совета

А/

Зуева Г.А.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Одним из видов возобновляемых энергетических ресурсов являются отходы деревообрабатывающих предприятий. Однако в настоящее время уровень использования древесной биомассы для получения энергии на отечественных деревообрабатывающих предприятиях не превышает 9% от общего энергетического потенциала древесных отходов. Основными проблемами, лимитирующими процесс сжигания древесных отходов, являются: во-первых, неустойчивость и низкая эффективность процесса горения древесных отходов; во-вторых, отходы деревообрабатывающих и мебельных предприятий образуют при сжигании высокотоксичные соединения.

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

Работа выполнялась в соответствии с планом НИОКР АН РТ (договор подряда № 07-7.5-229/2004).

Цель работы:

1) разработка эффективной схемы процесса энергетической переработки древесных отходов;

2) разработка методики расчета экологически безопасной установки для термической переработки древесных отходов с эффективной схемой утилизации тепла;

3) математическое моделирование технологического процесса энергетической переработки древесных отходов;

4) промышленная реализация результатов работы.

Научная новизна.

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

2) Впервые разработана методика расчета процесса термической переработки древесных отходов с учетом процессов предварительной сушки за счет тепла низкотемпературных газов, выбрасываемых в традиционных схемах в атмосферу, и очистки дымовых газов от токсичных соединений.

3) Разработан графический способ выбора рациональных условий проведения процесса в зависимости от свойств древесных отходов в виде диаграммы.

4) Разработана конструкция экспериментального стенда для исследования процесса горения древесных отходов, а также методика и соответствующее программное обеспечение для проведения и обработки экспериментальных исследований по кинетике горения.

Практическая ценность. На основе результатов моделирования, дающих возможность оценить влияние свойств древесных отходов на режимные параметры процесса термической переработки древесных отходов и выявить кинетические характеристики процесса горения, разработана методика расчета конструктивных и режимных параметров процесса сжигания древесных отходов. Полученная методика расчета была использована при разработке установки для термической переработки древесных отходов. Использование данной установки позволило: повысить эффективность сжигания древесных отходов на 15-30 %, снизить до допустимых норм выделение токсичных веществ в окружающую среду, а также обеспечить комплексное использование сырьевых ресурсов на предприятии. Представленные в графическом виде анализ результатов моделирования и разработанная диаграмма позволяют выявить пути усовершенствования аналогичных технологических процессов. Разработанные на базе экспериментального стенда алгоритмы и программное обеспечение могут быть использованы для всестороннего автоматизированного контроля над процессом горения в подобных устройствах.

Реализация работы. Результаты теоретических и экспериментальных исследований реализованы при создании методики расчета и проектировании промышленной установки для термической переработки древесных отходов, а также в учебном процессе на базе экспериментального стенда для изучения процесса сжигания древесных частиц. Внедрение установки для термической переработки древесных отходов в технологический цикл ЗАО «Ласкрафт» (г. Казань) позволило получить суммарный годовой экономический эффект в размере 325 тыс. руб. Использование экспериментального стенда для исследования процесса сжигания древесных частиц в учебном процессе позволяет студентам в ходе практических занятий по дисциплине «Технологические процессы и оборудование деревообрабатывающих производств» (ТПОДП) осуществлять всестороннее изучение процессов горения и газоочистки. Результаты исследований были также использованы при разработке: углевыжигательной печи, установки для сжигания газовых выбросов, вытяжного устройства для аккумуляторов и устройства для нанесения защитных покрытий на поверхности труб.

Автор защищает:

1) предложенную технологическую схему установки для термической переработки древесных отходов;

2) методику расчета установки для термической переработки древесных отходов;

3) результаты математического моделирования и экспериментального исследования;

4) конструкцию экспериментального стенда для исследования процесса горения древесных отходов, а также методику и соответствующее программное обеспечение для проведения и обработки результатов экспериментальных исследований по кинетике горения;

5) конструкцию промышленной установки для термической переработки древесных отходов.

Апробация работы» Основные научные положения и результаты работы докладывались на: IV республиканской научной конференции «Актуальные экологические проблемы Республики Татарстан», Казань, 2000 г.; Всероссийских научно-практических конференциях «Химико-лесной комплекс- проблемы и решения», Красноярск, 2002-2004 г.; Международной научно-технической конференции «Лес -2004». Брянск, 2004 г.; Всероссийской научно-технической конференции «Актуальные проблемы развития лесного комплекса», Вологда, 2004 г.; Международных научных конференциях «Математические методы в технике и технологиях», 2002, 2004г.; Международной научной конференции «Энерго-ресурсосберегающие технологии и оборудование, экологически безопасные производства», Иваново, 2004 г.; IV международном симпозиуме «Строение, свойства и качество древесины -2004», Санкт-Петербург, 2004 г.

Публикации. По теме диссертации опубликовано 15 печатных работ. Из них 1 статья и 4 патента на изобретения.

Объем и структура работы. Диссертационная работа состоит из введения, четырех глав, заключения, библиографического списка и приложений. Работа содержит 143 страницы текста основного содержания, 65 рисунков и 2 таблицы. Библиографический список включает 180 работ российских и зарубежных авторов.

На всех этапах работы в качестве научного консультанта принимал участие к.т.н., доцент Башкиров В.Н.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

решения которых были использованы численные методы. Данная модель включает в себя: математическую модель процесса сжигания древесного топлива, математическую модель процесса предварительной сушки древесного топлива, математическую модель процесса очистки газовых выбросов.

Математическая модель процесса сжигания древесного топлива. Продолжительность сгорания древесной частицы тч, можно представить в виде суммы продолжительностей процессов прогрева Т1Я) и сушки Тс, термического разложения с горением летучих X, и выгорания угля т,. Тогда выражение для определения Тч можно записать в следующем виде

Стадию прогрева древесной частицы можно описать дифференциальным уравнением нестационарной теплопроводности с соответствующими граничными условиями.

Изменение влагосодержания частицы в период постоянной скорости сушки можно определить с помощью уравнения

(1ш_МР"-Р>рч (2)

с!т р0у„

В процессе термического разложения массу летучих, выделившихся к данному моменту времени можно определить как произведение доли летучих веществ выделившихся при разложении древесины к моменту

времени т, на массу частицы. При этом долю летучих веществ ул(т), выделившихся при разложении древесной частицы к моменту времени можно определить на основании выражения

Ул(т) = У

лО

1 -ехр

-|к0 ех|

4

—1

(3)

Изменение размеров частицы в процессе выгорания углерода может быть записано в виде соотношения

аз

= -о.

2М.

(4)

с!т р.

Расчетные выражения скорости выгорания углерода в кислородной зоне необходимо выбирать в зависимости от значения критерия Семенова ве.

Если 8е<0,4, то выгорание частицы будет протекать по схеме негорящего пограничного слоя. Скорость выгорания коксовой частицы в этом случае можно определить по выражению

(5)

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

В этом случае выражение для определения скорости выгорания углерода запишется в виде

в. =

ЛТ

N.

N »

-— (Р, +—)

1 + К ' 2

(6)

1 + К,

При значениях критерия Семенова более 2 выгорание окиси углерода происходит в пределах пограничной пленки и скорость выгорания коксовой частицы можно определить по уравнению

(7)

Для восстановительной зоны выражение скорости выгорания углерода запишется в виде:

(8)

в. =-а„

Величина N в выражениях (5)-(8) представляет собой отношение константы скорости химической реакции к диффузионному коэффициенту массоотдачи, то есть

N =-¡- = N11, а„

(9)

и имеет смысл диффузионно-химического критерия подобия.

Зависимость константы скорости химической реакции от температуры определяется на основании закона Аррениуса. Температура в слое в зависимости от размера частицы может быть определена из уравнения теплового баланса,, которое запишется в виде выражения:

лб

КТяру-(5>слМ0г+53у02) 6'У°2'Р'|,РМ.

Оог-

(10)

Парциальные давления участвующих в реакции компонентов, в зависимости от текущего размера частицы, можно определить по выражениям: для кислородной зоны

Р2(5) = ^

а„„

\5оу

е у„ (б'а -51) Р.(8) = —• " -Р;

аС1 о

и для восстановительной зоны

у() (451 -25')

У о 5!, +25;,

(И) (12)

Причем значение размера частицы на границе кислородной зоны 6, , можно получить с помощью решения уравнения (12) относительно 8, которое при Р)=0 запишется в виде

Скорость продвижения частицы и ее положение в слое можно определить с помощью выражения

— = —(1-е)—. с1т 5„ (Зт

(16)

Совместное решение уравнений (4)-(16) позволяет рассчитать продолжительность выгорания коксовой частицы.

Математическая модель процесса предварительной сушки древесного топлива. Влагосодержание древесных отходов на выходе из сушильного

бункера можно определить при соблюдении теплового баланса

-со.

р, - с, +Р:

со,,-со, 03

+ С, I —-- > Г --

100 ~ 100

(Сг1ш(1_з.))+СЛь. 100 ""100

(17)

где значение температуры мокрого термометра можно определить по выражению

Значение равновесного влагосодержания древесной частицы в зависимости от параметров сушильного агента может быть определено по выражению

шр=13,094.(г,;Г,3.е*

(19)

При этом относительная влажность газов на выходе из сушильного бункера определяется из соотношения:

(20)

где давление насыщенного водяного пара в зависимости от температуры определяется уравнением Антуана. С другой стороны влагосодержание газа на выходе из сушильного бункера запишется в виде

X =-

V V + ДШ

у н,о I 1и> т " "

(21)

Совместное решение уравнений (19)-(21) относительно позволяет

определить значение температуры и относительную влажность топочных газов на выходе из сушильного бункера. Продолжительность высокотемпературной сушки древесной частицы можно определить по выражениям:

для первого периода сушки

Хф ' 72(1с-100) I, а + 400Х)' (22)

для второго периода сушки

_0,9Ьрг.«о„-а),у 1 Ь V .а>„-ю, " 72ас-юо) {а+ШГ*-^- (23)

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

М* = Вм(Чуг + д\у)(у., -?.) = ь(х„ -Х„). (24)

Плотность орошения вычисляется из соотношения

р И ' (25)

Коэффициент массоотдачи по газовой фазе в абсорберах с регулярной насадкой находят из соотношения

Р^ОЛбТ^е/^г/'-^у4'. (26)

Коэффициент массоотдачи в жидкой фазе может быть определен по выражению:

Необходимая высота слоя насадки рассчитывается по соотношению:

= 0,785«1Ч|/, ' (28)

Гидравлическое сопротивление орошаемой насадки рассчитывается по уравнению:

[ 9................

(29)

Др, = X ] уехр(2,303Уи)

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

В третьей главе изложены результаты математического моделирования и экспериментальных исследований процесса термической переработки древесных отходов, а также проведена проверка полученной математической модели на адекватность реальному процессу. В качестве модельного материала использовалась древесина сосны. Для проверки формализованной модели процесса термической переработки древесных отходов на адекватность реальному процессу был создан экспериментальный стенд, представленный на

рис. 1. Результаты экспериментальных исследований и результаты, полученные математическим моделированием в идентичных условиях.

Рис. 1. Схема экспериментального стенда для исследования процесса сжигания древесных

частиц

1- камера сгорания; 2-система очистки; 3- блок отбора проб; 4- колосниковая решетка; 5-дозатор; 6- нагревательный элемент, 7-ротаметр; 8-вентилятор, 9- компьютер, 10- аналого-цифровой преобразователь; 11-блок управления; 12,13 термопары; 14- цифровая видеокамера; 15-барботер; 16- насадочный абсорбер; 17-сепаратор; 18-приемная ванна.

представлены в виде графических зависимостей (рис. 2-7), на которых сплошными линиями изображены данные, полученные расчетным путем, точками - экспериментальные значения. Как видно из представленной на рис. 2. зависимости начальная влажность древесной частицы оказывает существенное влияние на продолжительность процесса сжигания в целом. Причем степень влияния распределяется следующим образом: при влажности до 30% время сжигания древесины меняется незначительно по линейной

Рис. 2 Зависимость времени сжигания Рис 3 Зависимость относительной древесной частицы от начальной температуры процесса горения от начальной

влажности влажности древесной частицы

зависимости вида т » аШ+Ь, а при влажности свыше 30% с возрастанием влажности увеличение времени сгорания более значимое и характеризуется

зависимостью вида т +с1. Снижение температурного уровня процесса с

увеличением влажности древесных частиц представлено на рис. 3. в виде

Рис 4. Зависимость продолжительности Рис 5. Зависимость относительной температуры процесса сжигания древесной частицы от процесса сжигания от коэффициента избытка

коэффициента избытка воздуха а воздуха

графической зависимости относительной температуры процесса горения от начальной влажности древесной частицы. Следует отметить, что значение влажности древесных частиц оказывает существенное влияние на продолжительность стадии прогрева и термического разложения древесины, практически не оказывая влияния на длительность стадии выгорания коксового остатка. Влияние коэффициента избытка воздуха на продолжительность и температурный уровень процесса представлено в виде графических зависимостей на рис. 4. и рис. 5. Максимальный температурный уровень с минимальной продолжительностью процесса сжигания достигается при значении коэффициента избытка воздуха равном 1,15.

Изменение продолжительности стадии выгорании коксового остатка в зависимости от размера частицы представлено на рис. 6. Из приведенной на рис. 7. зависимости скорости выгорания коксового остатка от текущего размера частицы в противоточном слое видно, что скорость выгорания коксового остатка по мере нагревания топочных газов в слое интенсивно увеличивается до тех пор, пока концентрация кислорода в топочных газах достаточно велика. Дальнейшее уменьшение концентрации кислорода в топочных газах приводит к снижению скорости выгорания коксовых частиц.

(1.6658 0,0086

Рис. 6. Зависимость продолжительности

выгорания коксового остатка в противоточном слое от размера частицы

Рис 7 Зависимость скорости выгорания коксового остатка от текущего размера частицы в противоточном слое

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

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

1- транспортер; 2-сушильный бункер, 3-шнековый питатель, 4- газогенераторная топка, 5-водогрейный котел; 6- система очистки

термической переработке древесных отходов, а также учебно-методическая программная среда, предназначенная для проведения практических занятий на базе экспериментального стенда (рис. 1.) с использованием компьютера по курсу «ТПОДП». Полученные данные были использованы при проектировании и изготовлении промышленной установки для термической переработки древесных отходов, которая представлена на рис. 8. Данная установка внедрена на ЗАО «Ласкрафт» с целью утилизации отходов столярного цеха и получения тепловой энергии. Анализ результатов испытаний установки, показал высокую эффективность процесса сжигания древесных отходов по данной технологической схеме. Применение данной установки на ЗАО «Ласкрафт» позволяет: снизить расход топлива при сохранении номинальной мощности котла, повысить стабильность процесса при сжигании топлива с высоким влагосодержанием, повысить энергоэффективность процесса сжигания в целом

за счет более полного использования тепла отходящих газов, а также снизить выбросы токсичных компонентов дымовых газов в воздушный бассейн за счет применения многоступенчатой системы очистки газовых выбросов

Результаты проведенных исследований были использованы также при разработке: углевыжигательной печи, установки для сжигания газовых

выбросов (рис 9.), вытяжного устройства для аккумуляторов, устройства для нанесения защитных покрытий на поверхности труб (рис. 10.)

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

1. Проведен анализ современного состояния технологического процесса термической переработки древесных отходов

2. Выявлена актуальность совершенствования технологической схемы процесса термической переработки древесных отходов.

3 На основе физической картины процесса и принятых допущений создана математическая модель процесса термической переработки древесных отходов с учетом предварительной сушки древесного топлива и абсорбционной очистки отходящих топочных газов

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

5. Получены зависимости, позволяющие оценить влияние режимных параметров и конструктивных характеристик на эффективность процесса термической переработки древесных отходов. Разработана методика графического расчета параметров сушильного бункера и других элементов оборудования энергетической переработки древесных отходов.

6. В ходе моделирования на основе обобщенного экспериментального материала получены эмпирические выражения по определению

Рис 9 Схема установки для сжигания ™с '^Устройство для нанесения защитных газовых выбросов (Пат № 2184909) покрытий на поверхности труб (Пат №

2170623)

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

7. Результаты теоретических и экспериментальных исследований процесса были использованы при проектировании и разработке ряда опытно-промышленных и промышленных установок, одна из которых внедрена на ЗАО «Ласкрафт» с экономическим эффектом в размере 325 тыс. руб. Промышленные испытания и технико-экономические исследования показали высокую эффективность применения предложенных технических решений в условиях производства.

8. Разработанные на базе экспериментального стенда алгоритмы и программное обеспечение могут быть использованы для всестороннего автоматизированного контроля над процессом горения в подобных устройствах.

9. Использование экспериментального стенда для исследования процесса сжигания древесных частиц в учебном процессе позволяет студентам в ходе практических занятий по дисциплине «ТПОДП» осуществлять всестороннее изучение процессов горения и газоочистки.

Основные обозначения со- влажность, %;т- время, сек; В,,-расход топлива, кг/сек;5- размер частицы, м;р-плотность, кг/м3; Р-площадь, м2; р, Р -парциальное, общее давление, Па; у-мольная доля, кМоль/кМоль; у,У-объем, удельный объем, м3, м3/кг;<р-отн. влажность газовой среды; Я- унив. газовая постоянная, Дж/(Кмоль К); с- массовая теплоемкость, Дж/(кг К); е -порозность; М-молярная масса, кг/кМоль; Н -высота, м; Э-коэффициент диффузии м/с; Т, (температура, К(°С); у-удельный вес, кг/м3; Др-сопротивление, Па; а- удельная поверхность, м2/кг; х- влагосодержание газа, кг/кг; )3, а„-коэффициент массоотдачи м/с; Е-энергия активации, кДж/моль; к-константа скорости химической реакции, Усек; -скорость выгорания углерода, кМоль/(м2 с); Скудельная теплота, кДж/кг; !- энтальпия, кДж/кг (кДж/м3); V-доля содержания компонента;«-коэффициент избытка воздуха;9-скорость течения, м/с; г-теплота парообразования, кДж/кг; Ь-толщина частицы, мм; Т-коэффициент, зависящий от характеристик насадки; II-плотность орошения, м/с; МаГ,-расход извлекаемого компонента, кг/с; Х-ср. отн. концентрация компонента в жидкой фазе; У-ср. отн. концентрация компонента в газовой фазе;(3-диаметр, м; <3Х1Ш-уд. теплота сгорания, кДж/кг; ДW- удельное количество удаляемой влаги, кг/кг;

Индексы: ч- частица; пр- прогрев; с-сушка; л- термическое разложение, летучие; у- выгорание коксового остатка, углерод; п- пар; в- вода, воздух; н-насыщенный, начальный; к- конечный; 1- кислород ; 2- углекислый газ; 3-окись углерода; 4-реакция горения окиси углерода; 0- при нормальных условиях, начальный ; к.з.- кислородная зона; в.з.- восстановительная зона; сл-слой; пот- потерь; г-газ; пол - полный; пов -поверхность ; ка - котел; сух-

сухой; м- мокрый, меньший; б- бункер, большее; пер- переходный; ср-среднее; ж- жидкость; нсд- насадка; э-эквивалентный; у- газовая фаза; х- жидкая фаза.

Основные положения диссертации опубликованы в следующих работах:

1. Грачев А.Н., Башкиров В.Н., Сафин Р.Г. Использование методов приближения при моделировании процесса термической переработки древесных отходов // Химия и химическая технология. -2004. -Т. 47. -№10. -С. 137-140.

2. Патент РФ № 2232348. Установка для термической переработки твердых отходов /Сафин Р.Г., Башкиров В.Н., Грачев А.Н. и др. -Бюл. № 18 . 2004.

3. Патент РФ № 2161349. Вытяжное устройство для аккумуляторов. /Сафин Р.Г., Лашков В.А., Грачев А.Н. и др. -Бюл. № 36, 2000.

4. Патент РФ № 2170623. Устройство для нанесения защитных покрытий на поверхности фуб. /Сафин P.I ., Лашков В.А., Грачев А.И. и др. -Бюл. № 25, 2001.

5. Патент РФ № 2184909. Установка для сжигания газовых выбросов /Сафин Р.Г., Лашков В.А.. Грачев А.Н. и др. -Бюл. № 34.2002.

6. Грачев А.Н, Сафин Р.Г. и др. Комплексное решение экологических задач на автотранспортных предприятиях республики Татарстан.// Актуальные экологические проблемы Республики Татарстан: Мат. IV респ. науч. конф.. -Казань: Новое знание, 2000. С. 131-132.

7. Грачев А.Н., Сафин Р.Г., Башкиров В.Н. Термическое обезвреживание твердых промышленных отходов // ММТТ-15. -Тамбов: Тамбовский гос. техн. ун-т, ТА 2002. -С.85-86.

8. Грачев А.Н., Сафин Р.Г., Башкиров В.Н. Использование отходов деревообрабатывающих предприятий в качестве дополнительного источника тепловой энергии // Всероссийская научно-практическая конференция «Химико -лесной комплекс - проблемы и решения». -Красноярск: 2002. -Т. II. -С. 286288.

9. Башкиров В.Н., Сафин Р.Г., Грачев А.Н. Расчет количества дымовых газов и определение их состава при сжигании отходов деревообрабатывающих предприятий // Тезисы докладов всероссийской научно-практической конференции «Лесной и химический комплекс - проблемы и решения». -Красноярск: 2003. -С. 189-193.

10. Валеев И.А., Сафин P.P., Грачев А.Н., и др. Ресурсосберегающая технология переработки древесных отходов // «Лес -2004».- Брянск: 2004. -С. 121-123.

11. Грачев А.Н, Башкиров В.Н., Сафин Р.Г. Расчетно-экспериментальное исследование времени сгорания влажных древесных отходов // МТТ-17. -Кострома: изд. КГТУ, 2004. -Т.9. -С. 108-110.

12. Грачев А.Н., Башкиров В.Н. Сафин Р.Г. и др. Повышение энергетической эффективности процесса термической переработки отходов в деревообрабатывающих производствах // Международная научная конференция

Of./7-PS. 2/

«Энерго-ресурсосберегающие технологии и оборудование, экологически безопасные производства». -Иваново: 2004. -Т.2. -С. 47.

13. Сафин P.I., Валеев И.А., Грачев А.Н., и др. Математическое моделирование процесса пиролиза // «Строение, свойства и качество древесины -2004». -Санкт-Петербург: СПбЛТА, 2004. -T.I. -C.342-343.

14. Ьашкиров В.Н., Сафин P.P., 1 рачев А.Н., и др. Перспективы развития лесного комплекса // Материалы Всероссийской научно-технической конференции «Актуальные проблемы развития лесного комплекса». -Вологда: 2004. -С.33-35.

15. Сафин P.P., Валеев И.А., Грачев А.Н. и др. Установка для пирогенетической переработки древесных отходов // Всероссийская научно-практическая конференция «Лесной и химический комплексы- проблемы и решения». -Красноярск: 2004. -С. 65-70.

Соискатель

А.Н. Грачев

Заказ Тираж 100 экз.

Офсетная лаборатория Казанского государственного технологического университета 420015,Казанц.К Маркса, 68

590

ВВЕДЕНИЕ.

Глава 1. ТЕОРЕТИЧЕСКИЕ ОСНОВЫ И АППАРАТУРНОЕ ОФОРМЛЕНИЕ ПРОЦЕССА ТЕРМИЧЕСКОЙ

ПЕРЕРАБОТКИ ДРЕВЕСНЫХ ОТХОДОВ.

1.1. Классификация и характеристика древесных отходов как топлива.

1.2. Современное состояние технологии и аппаратурного оформления процесса сжигания древесных отходов.

1.3. Анализ систем очистки дымовых газов.

1.4. Анализ теоретических основ сжигания древесных отходов.

Выводы.

Постановка задачи.

Глава 2. МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ПРОЦЕССА ТЕРМИЧЕСКОЙ

ПЕРЕРАБОТКИ ДРЕВЕСНЫХ ОТХОДОВ.

2.1. Физическая картина процесса сжигания древесных отходов.

2.2. Формализация процесса.

2.3. Методика расчета установки для термической переработки древесных отходов.

2.3.1. Математическая модель процесса сжигания древесного топлива.

2.3.2. Математическая модель процесса предварительной сушки древесного топлива.

2.3.3. Математическая модель процесса очистки газовых выбросов

2.4. Алгоритм расчета процесса термической переработки древесных отходов.

Выводы.

Глава 3. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И

ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ПРОЦЕССА w ТЕРМИЧЕСКОЙ ПЕРЕРАБОТКИ ДРЕВЕСНЫХ ОТХОДОВ.

3.1. Описание экспериментального стенда для исследования процесса сжигания древесных частиц.

3.2. Методика проведения экспериментов и обработка экспериментальных данных.

3.3. Анализ результатов математического моделирования и tv экспериментальных данных.

Выводы.

Глава 4. ПРОМЫШЛЕННАЯ РЕАЛИЗАЦИЯ РЕЗУЛЬТАТОВ

ИССЛЕДОВАНИЯ ПРОЦЕССОВ, ПРОТЕКАЮЩИХ ПРИ ТЕРМИЧЕСКОЙ ПЕРЕРАБОТКЕ ДРЕВЕСНЫХ

ОТХОДОВ.

4.1. Исследование структурно-механических свойств древесных отходов столярного цеха.

4.2. Разработка установок для сжигания древесных отходов.

4.2.1. Опытная установка для сжигания древесных отходов.

4.2.2. Промышленная установка для сжигания древесных отходов.

4.3. Реализация результатов исследований в учебном процессе.

4.4. Разработка новых технологических процессов и оборудования на основе проведенных исследований.

Выводы.

Актуальность темы. Главными приоритетами программы развития отечественной энергетики являются повышение эффективности использования энергетических ресурсов и снижение негативного воздействия на окружающую среду /98, 99/, что достигается разработкой новых технологий и использованием возобновляемых ресурсов.

Одним из видов возобновляемых энергетических ресурсов являются отходы деревообрабатывающих предприятий. Однако в настоящее время уровень использования древесной биомассы для получения энергии на отечественных деревообрабатывающих предприятиях не превышает 9% от общего энергетического потенциала древесных отходов /104/. При этом комплексное использование отходов деревообработки в энергетическом балансе страны дает возможность значительно экономить традиционные виды топлива, а также снизить выделение парниковых газов в атмосферу.

Основными проблемами, лимитирующими процесс сжигания древесных отходов, являются: во-первых, неустойчивость и низкая эффективность процесса горения древесных отходов; во- вторых, отходы деревообрабатывающих и мебельных предприятий образуют при сжигании высокотоксичные соединения.

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

Работа выполнялась в соответствии с планом НИОКР АН РТ (договор подряда № 07-7.5-229/2004).

Цель работы:

1) разработка эффективной схемы процесса энергетической переработки древесных отходов;

2) разработка методики расчета экологически безопасной установки для термической переработки древесных отходов с эффективной схемой утилизации тепла;

3) математическое моделирование технологического процесса энергетической переработки древесных отходов;

4) промышленная реализация результатов работы.

Научная новизна.

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

2) Впервые разработана методика расчета процесса термической переработки древесных отходов с учетом процессов предварительной сушки за счет тепла низкотемпературных газов, выбрасываемых в традиционных схемах в атмосферу, и очистки дымовых газов от токсичных соединений.

3) Разработан графический способ выбора рациональных условий проведения процесса в зависимости от свойств древесных отходов в виде диаграммы.

4) Разработана конструкция экспериментального стенда для исследования процесса горения древесных отходов, а также методика и соответствующее программное обеспечение для проведения и обработки экспериментальных исследований по кинетике горения.

Практическая ценность. На основе результатов моделирования, дающих возможность оценить влияние свойств древесных отходов на режимные параметры процесса термической переработки древесных отходов и выявить кинетические характеристики процесса горения, разработана методика расчета конструктивных и режимных параметров процесса сжигания древесных отходов. Полученная методика расчета была использована при разработке установки для термической переработки древесных отходов. Использование данной установки позволило: повысить эффективность сжигания древесных отходов на 15-30 %, снизить до допустимых норм выделение токсичных веществ в окружающую среду, а также обеспечить комплексное использование сырьевых ресурсов на предприятии. Представленные в графическом виде анализ результатов моделирования и разработанная диаграмма позволяют выявить пути усовершенствования аналогичных технологических процессов. Разработанные на базе экспериментального стенда алгоритмы и программное обеспечение могут быть использованы для всестороннего автоматизированного контроля над процессом горения в подобных устройствах.

Реализация работы. Результаты теоретических и экспериментальных исследований реализованы при создании методики расчета и проектировании промышленной установки для термической переработки древесных отходов, а также в учебном процессе на базе экспериментального стенда для изучения процесса сжигания древесных частиц. Внедрение установки для термической переработки древесных отходов в технологический цикл ЗАО «Ласкрафт»(г. Казань) позволило получить суммарный годовой экономический эффект в размере 325 тыс. руб. Использование экспериментального стенда для исследования процесса сжигания древесных частиц в учебном процессе позволяет студентам в ходе практических занятий по дисциплине «Технологические процессы и оборудование деревообрабатывающих производств» (ТПОДП) осуществлять всестороннее изучение процессов горения и газоочистки. Результаты исследований были также использованы при разработке: углевыжи-гательной печи, установки для сжигания газовых выбросов, вытяжного устройства для аккумуляторов и устройства для нанесения защитных покрытий на поверхности труб. Автор защищает:

1) предложенную технологическую схему установки для сжигания древесных отходов;

2) методику расчета установки для сжигания древесных отходов;

3) результаты математического моделирования и экспериментального исследования;

4) конструкцию экспериментального стенда для исследования процесса горения древесных отходов, а также методику и соответствующее программное обеспечение для проведения и обработки экспериментальных исследований по кинетике горения;

5) конструкцию промышленной установки для сжигания древесных отходов.

Апробация работы. Основные научные положения и результаты работы докладывались на: IV республиканской научной конференции «Актуальные экологические проблемы Республики Татарстан», Казань, 2000 г.; Всероссийских научно-практических конференциях «Химико -лесной комплекс-проблемы и решения», Красноярск, 2002-2004 г.; Международной научно-технической конференции «Лес -2004», Брянск, 2004 г.; Всероссийской научно-технической конференции «Актуальные проблемы развития лесного комплекса». Вологда, 2004 г.; Международных научных конференциях «Математические методы в технике и технологиях». 2004г., 2002г.; Международной научной конференции «Энерго-ресурсосберегающие технологии и оборудование, экологически безопасные производства», Иваново, 2004 г.; IV международном симпозиуме «Строение, свойства и качество древесины -2004», Санкт-Петербург, 2004 г.

Публикации. По теме диссертации опубликовано 15 печатных работ. Из них 1 статья и 4 патента на изобретения.

Объем и структура работы. Диссертационная работа состоит из введения, четырех глав, заключения, библиографического списка и приложений. Работа содержит 143 страницы текста основного содержания, 65 рисунков и 3 таблицы. Библиографический список включает 180 работ российских и зарубежных авторов. В первой главе проведен анализ современного состояния теории и техники по изучаемому вопросу. Во второй главе представлены математическое описание процесса термической переработки древесных отхо

Выводы

На основе результатов теоретических и экспериментальных исследований разработана инженерная методика расчета процессов, протекающих при сжигании древесных отходов, а также учебно-методическая программная среда, предназначенная для проведения практических занятий на базе экспериментального стенда с использованием компьютера по курсу «Технологические процессы и оборудование деревообрабатывающих производств».

Полученные данные были использованы при проектировании и изготовлении ряда опытно-промышленных и промышленных установок, одна из которых внедрена на ЗАО «Ласкрафт». Промышленные испытания и технико-экономические исследования показали эффективность применения предложенных технических решений в условиях производства.

ЗАКЛЮЧЕНИЕ

Настоящая работа содержит научно-обоснованные решения, направленные на совершенствование процесса термической переработки древесных отходов с целью получения тепловой энергии.

Аналитический обзор теории и практики по данному вопросу показал, что энергетическое использование древесных отходов лимитируется повышенной влажностью и нестабильностью состава. Термическая переработка на предприятиях деревообрабатывающей промышленности преимущественно осуществляется в топочных устройствах слоевого типа. Предварительное снижение влажности древесных отходов и применение абсорбционный системы очистки в топочных устройствах позволит решить ряд актуальных проблем, возникающих при утилизации древесных отходов. Однако отсутствие обобщенных методик расчета процесса термической переработки сдерживает его промышленное использование. Поэтому целесообразным представляется: исследование взаимосвязанных процессов термической переработки, предварительной сушки топлива и абсорбционной очистки отходящих дымовых газов; разработка обобщенной математической модели всего комплекса явлений процесса на основе известных положений теории тепло- и массообмена; создание усовершенствованных схем энергетической переработки древесных отходов; выдача рекомендаций по практическому применению результатов исследования.

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

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

Результаты теоретических и экспериментальных исследований реализованы при создании методики расчета и проектировании промышленной установки для сжигания древесных отходов, а также в учебном процессе в виде экспериментального стенда для изучения процесса сжигания древесных частиц. Использование установки для сжигания древесных отходов в технологическом цикле ЗАО «Ласкрафт» позволило получить суммарный годовой экономический эффект в размере 325 тыс. руб. Использование экспериментального стенда для исследования процесса сжигания древесных частиц в учебном процессе позволяет студентам в ходе практических занятий по дисциплине «Технологические процессы и оборудование деревообрабатывающих производств» осуществлять всестороннее изучение процессов горения и газоочистки. Результаты исследований были также использованы при разработке: углевыжигательной печи, установки для сжигания газовых выбросов, вытяжного устройства для аккумуляторов и устройства для нанесения защитных покрытий на поверхности труб, новизна решений которых подтверждена четырьмя патентами РФ.

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ

Зольность, %; коэффициенты в уравнении Антуана; коэффициент теплообмена излучением и конвекцией, Дж/(м2 с К); коэффициент избытка воздуха; коэффициент температуропроводности м2/с; удельная поверхность, м2/кг расход топлива, кг/сек; коэффициент массоотдачи м/с; определяющий размер частицы, м; толщина частицы, мм; коэффициент диффузии м/с; диаметр, м; порозность; энергия активации, кДж/моль; площадь поверхности, м2; относительная влажность газовой среды; текущая масса, кг; У скорость выгорания углерода, кМоль/(м с); высота, м; удельная энтальпия, кДж/кг, кДж/м3; поток массы кг/(м2 с) коэффициент теплопередачи, массопередачи, Вт/(м2 К), л кмоль/(м с (ед.дв.с.)); константа скорости химической реакции, 1/сек; масса извлекаемого компонента, кг/с; диффузионно-химический критерий подобия; количество вещества, моль; мощность, Вт; коэффициент полезного действия, %; парциальное давление, Па; сопротивление, Па; плотность, кг/м3; общее давление, Па; удельная теплота, кДж/кг; удельная мощность внутреннего источника, Вт/м3; теплонапряжение колосниковой решетки, кВт/(м2 ч) теплота парообразования, кДж/кг; универсальная газовая постоянная Дж/(Кмоль К) объем, удельный объем, м3, м3/кг; доля объемного содержания компонента; скорость течения фазы, м/с; влажность, %; удельное количество удаляемой влаги, кг/кг; мольная доля; плотность орошения, м/с; параметр зависящий от формы частицы; молярная масса, г/моль; теплоемкость удельная Дж/(кг К) текущее время, продолжительность, сек; текущая температура, °С(К); средняя относительная концентрация компонента в жидкой фазе, кМоль/ кМоль; средняя относительная концентрация компонента в газовой фазе, кМоль/ кМоль; текущая координата, влагосодержание газа, кг/кг;

X - коэффициент теплопроводности Вт/(м К)

- коэффициент, зависящий от характеристик насадки;

Индексы: ч- частица; пр- прогрев; с-сушка; л- термическое разложение, летучие; у- выгорание коксового остатка, углерод; п- пар; в- вода, воздух; н- насыщенный; н- начальный; к- конечный; хим- химический; 1- кислород ; 2-углекислый газ; 3- окись углерода; 0- при нормальных условиях, начальный ; к.з.- кислородная зона; в.з- восстановительная зона; сл- слой; пот- потерь; г-газ; пол - полный; пов -поверхность ; ка — котел; сух- сухой; м- мокрый, меньший; б- бункер, большее; пер- переходный; ср-среднее; ж- жидкость; нсд- насадка; э-эквивалентный; у- газовая фаза; х- жидкая фаза.

1. Агроскин А. А. Физика угля. М.: Недра, 1965.-352 с.

2. Ананьин П.И., Петри В.П. Высокотемпературная сушка древесины. М.: Гослесбумиздат, 1963. - 305 с.

3. Ахназарова C.JL, Кафаров В.В. Оптимизация эксперимента в химии и химической технологии. М.: Высш. школа, 1978.- 319с.

4. Аэров М.Э., Тодес О.М., Наринский Д.А. Аппараты со стационарным зернистым слоем. Л.: Химия, 1979. - 176 с.

5. Аэродинамический расчет котельных установок (нормативный метод) Под ред. С.И. Мочана. Изд. 3-е. Л., Энергия, 1977. 256 с.

6. Бабий В.И., Куваев Ю.Ф. Горение угольной пыли и расчет пыле-угольного факела. М.: Энергоатомиздат, 1986. - 208 с.

7. Баженов В.А., Карасев Е.И., Мерсов Е.Д. Технология и оборудование производства древесных плит и пластиков. М.: Экология, 1992. - 416 с.

8. Бахмачевский Б.И. Теплотехника. -М.: Энергия, 1963. -725с.

9. Белосельский Б.С., Соляков В.К. Энергетическое топливо. М.: Энергия, 1980.-168 с.

10. Боровиков A.M., Уголев Б.И. Справочник по древесине. М.: Лесная пром-сть, 1989. - 296 с.

11. Бретшнайдер С. Свойства газов и жидкостей. Инженерные методы расчёта. М.; Л.: Химия, 1966. 535 с.

12. Быстров А.Ф. Основы для эффективного использования древесныхотходов деревообрабатывающего предприятия // Деревообрабатывающая промышленность. -1999. -№5, -С. 27-28.

13. Ведерников H.A. Клеточная стенка древесины. -Рига: Зинатне. 1972. -510 с.

14. Вержбицкий В.М. Численные методы (математический анализ и обыкновенные дифференциальные уравнения). -М.:Высш. шк., -2001. -382 с.

15. Виленкин С.Я. Статистическая обработка результатов исследования случайных функций. -М.: Энергия, 1979.-320 с.

16. Вильяме Ф.А. Теория горения. М.: Наука, 1976.

17. Внуков А.К. Защита атмосферы от выбросов энергообъектов: Справ. -М.: Энергоатомиздат, 1992.-176 с.

18. Волков Е.А. Численные методы. М.: Наука, 1987. - 248 с.

19. Вопросы аэродинамики и теплопередачи в котельно-топочных процессах. Под ред. Г. Ф. Кнорре. Л., Машгиз, 1958.

20. Вулис Л.А Тепловой режим горения. -М.: Госэнергоиздат, 1958, -340с

21. Гаврилов Л.Ф., Малкин Б.М. Загрязнение и очистка поверхностей нагрева котельных установок. М.: Энергия, 1980. - 328 с.

22. Гидравлический расчет котельных агрегатов (нормативный метод). -М.: Энергия, -1978. -256с.

23. Головков С.И. Энергетическое использование древесных отходов. -М.: Лесная пром-ть. -1987.-224с.

24. Голубев Л.Г. Древесиноведение. -Казань: Изд-во КГТУ, -2000. 64с.

25. ГОСТ 11022-75 Угли бурые, каменные, антрацит и сланцы горючие. Метод определения зольности. Взамен ГОСТ 11022-64; Введ. с 01.01.76. -М.: Издательство стандартов, 1975. - 6с.

26. ГОСТ 147-74 (СТ СЭВ 1463-78) Топливо твердое. Метод определения высшей теплоты сгорания и вычисление низшей теплоты сгорания. -Введ. с 01.01.75. М.: Издательство стандартов, 1985. - 20 с.

27. ГОСТ 2408.1-88 (СТ СЭВ 6040-87) Топливо твердое. Методы определения углерода и водорода. Взамен ГОСТ 6389-81, ГОСТ 2408.175; Введ. с 01.07.89. - М.: Издательство стандартов, 1989. - 22 с.

28. ГОСТ 6382-80 (СТ СЭВ 2033-79) Угли бурые, каменные, антрацит и сланцы горючие. Метод определения выхода летучих веществ. Взамен ГОСТ 6382-75; Введ. с 01.01.83. - М.: Издательство стандартов, 1987.- 6 с.

29. ГОСТ 8606-72 (СТ СЭВ 1462-78, СТ СЭВ 2230-80) Топливо твердое. Метод определения серы. Взамен ГОСТ 8606-68; Введ. с 01.07.73. -М.: Издательство стандартов, 1983. - 9 с.

30. Грачев А.Н, Башкиров В.Н., Сафин Р.Г. К расчету состава дымовых газов при сжигании древесных отходов. // Аннотации сообщений научной сессии КГТУ.- Казань: 2003. -С.116.

31. Грачев А.Н, Башкиров В.Н., Сафин Р.Г. Очистка дымовых газов при сжигании древесных отходов // Аннотации сообщений научной сессии КГТУ. -Казань: -2003. С.117.

32. Грачев А.Н., Башкиров В.Н., Сафин Р.Г. Использование методов приближения при моделировании процесса термической переработки древесных отходов // Химия и химическая технология. -2004. -Т. 47.10, -С. 137-140.

33. Грачев А.Н., Сафин Р.Г., Башкиров В.Н. Термическое обезвреживание твердых промышленных отходов// Математические методы в технике и технологиях.- ММТТ-15. -Тамбов: Тамбовский гос. техн. ун-т, -Т.4.-2002.-С.85-86.

34. Грачев А.Н.,Сафин Р.Г., Башкиров В.Н. Использование отходов деревообрабатывающих предприятий в качестве дополнительного источника тепловой энергии // Всероссийская научно-практическая конференция «Химико -лесной комплекс- проблемы и решения».

35. Красноярск: -2002. -Т. И. -С. 286-288.

36. Григорьев JI.H. МТС Емеля. //Деревообрабатывающая промышленность. -2000. -№4, -С. 32-33.

37. Джирис, Р, Винтербак Ю., Энгберг Ю., Охман М. Изменения свойств древесины и их влияние на качество гранул. КТНВ, Studsvikbibliote-ket. 2002. ТВ 2-12.

38. Доброхотов В.И. К проблеме воздействия энергетики на окружающую среду // Теплоэнергетика. 1995. - № 2. - С. 2-5.

39. Дрейпер П., Смит Г. Прикладной регрессионный анализ. М.: Статистика, 1973, 392 с. 96. W. Francis. Fuels and Fuel Technology. OxfordLondon, 1965.-120 c.

40. Дыбок В.В. Получение синтетических моторных топлив при утилизации древесных отходов. //Лесная промышленность. 1999. №1, с 1820.

41. Дытнерский Ю.И. Процессы и аппараты химической технологии. Изд. 3-е. В 2-х кн.: Часть 2. Массообменные процессы и аппараты. М.: Химия, 2002.-368 с.

42. Житков А.В. Утилизация древесной коры. М.: Лесн. пром-ть, 1985.-136с.

43. Защита окружающей среды от техногенных воздействий: Учебное пособие. Под общ. ред. Невской Г.Ф.М: Изд-во МГОУ, 1993, 218 с.

44. Зельдович Я.Б., Барнеблатт Г.И., Либрович В.Б. и др. Математическая теория горения и взрыва. -М.: Наука, -1980. -560 с.52