ОЦЕНКА СТРУКТУРНЫХ И ТЕХНОЛОГИЧЕСКИХ ХАРАКТЕРИСТИК УГЛЕЙ ПО ДАННЫМ ДЕРИВАТОГРАФИИ И ПИРОЛИЗА
Шевкопляс В.Н., Лящук С.Н., Бутузова Л.Ф. Химия твердого топлива 2004 год №5 стр. 68-75.
Основная задача, которая решается при исследовании термических превращений твердых горючих ископаемых (ТГИ) - это поиск рациональных и эффективных путей их использования. Эта задача требует выявления оптимальных условий переработки и экспресс-методов оценки характера термических превращений ТГИ. На сегодня прикладная углехимия располагает множеством методик, позволяющих прогнозировать свойства ТГИ и их поведение в различных термических процессах. Прошли проверку временем ряд подразделений углей в классификации, которые используются на практике для оценки их технологических параметров. Основными из них являются генетические, промышленные и промышленно-генетические классификации, охватывающие угли всего ряда метаморфизма от бурых и каменных до антрацитов.
В последние годы в литературе опубликован ряд работ, использующих математико-статистический анализ для нахождения зависимостей между технологическими показателями углей, их техническим, элементным составом и данными физико-химических методов исследования. Предложены уравнения, которые позволили выявить наличие линейной связи между выходом летучих веществ и содержанием органической серы для углей низкой стадии метаморфизма [1]; уравнения для прогнозирования выхода H2S из высокосернистых углей при гидрогенизации [2]; для моделирования структурных характеристик органической массы угля (ОМУ) по показателям энтальпии образования и атомизации индивидуальных соединений [3]; для оценки вклада отдельных атомов в исследуемые свойства углей, а также для построения классификатора углей по их технологическим показателям [4].
Однако многие из используемых в химической технологии топлив методов основаны на сложившихся стереотипах в понимании структуры ТГИ, которые не выходят за рамки традиционных представлений и не отвечают современным требованиям. Возникла необходимость существенного пересмотра некоторых из них, тем более что для этого имеются объективные предпосылки, а именно: накоплен большой массив экспериментальных данных, полученных с использованием традиционных и современных методов анализа. Статистическая обработка имеющихся данных позволит создать научный базис для оценки качества и свойств углей, расширить границы наших представлений о составе и структуре углеродных материалов и, следовательно, сделает возможным получить более объективную информацию об их технологических характеристиках.
Одним из вышеупомянутых методов является «метод определения выхода продуктов полукоксования» (ГОСТ 3168-93, ИСО 647-74), который предусматривает термическое разложение топлив при температуре 520 оС без доступа воздуха с образованием пирогенетической воды, смолы, газа и полукокса [5]. Выход продуктов полукоксования – это показатель, отражающий химическую природу, строение и свойства органической массы углей. Он является одним из классификационных параметров, широко использующихся в промышленности для оценки технологических характеристик ТГИ с целью более эффективного их использования в различных процессах конверсии [6-8]. Основным недостатком данного метода является фиксированная температура проведения пиролиза - 520оС, которая, по данным дериватографии, не соответствует температуре собственно разложения наименее устойчивой части органической массы (Тмак. на кривой ДТГ). Для углей Б и Д эта температура завышена, а для углей марки Т и антрацитов - занижена. Следовательно, метод предназначен, главным образом, для оценки технологических характеристик углей средней стадии метаморфизма (угли марок Г, Ж, К и ОС) и не может быть распространен на угли всего ряда метаморфизма.
Целью данной работы было установление зависимостей «состав-свойство» для полного ряда ТГИ, основываясь на результатах определения выхода продуктов пиролиза и данных дериватографии. Решение поставленной задачи требовало выбора оптимальной температуры проведения пиролиза, которая была бы в равной степени приемлемой для всех углей и позволяла связать их стадию метаморфизма с технологическими характеристиками. На первом этапе исследований был выбран метод дериватографии, который позволяет определить искомую температуру. По виду дериватограмм (кривые ДТА, ДТГ и ТГ) возможно проследить за ходом термодеструкции ТГИ в интервале температур от 20 до 950 оС и охарактеризовать процесс при каждой конкретной температуре.
В качестве объектов исследования были использованы бурые угли Днепровского бассейна, каменные угли и антрациты Донецкого бассейна, характеристика которых представлена в таблице 1.
Таблица 1
Характеристика объектов исследования |
||||||||||
№ п/п |
Шахта, марка угля, |
Технический анализ, % |
Элементный анализ, (daf), % |
|
||||||
Wa |
Ad |
Sdt |
Vdaf |
C |
H |
O+N |
|
|
||
1. |
ш.Новомиргородская, Б |
42,0 |
11,9 |
2,9 |
63,2 |
65,4 |
6,1 |
28,5 |
|
|
2. |
р. Константиновский, Б |
47,4 |
11,2 |
4,1 |
61,0 |
66,8 |
5,8 |
27,4 |
|
|
3. |
Кураховская, Д, l2 |
3,0 |
5,3 |
1,0 |
37,0 |
76,2 |
5,1 |
14,9 |
|
|
4. |
Кураховская, Д, l4 |
1,9 |
8,6 |
5,6 |
43,0 |
79,0 |
5,2 |
13,0 |
|
|
5. |
Лидиевская, Г, l3 |
1,6 |
0,4 |
2,3 |
38,0 |
82,2 |
5,4 |
10,9 |
|
|
6. |
Лидиевская, Г, l1 |
1,2 |
3,8 |
1,0 |
35,0 |
82,7 |
5,3 |
11,0 |
|
|
7. |
им. Засядько, Ж k8 |
- |
6,3 |
4,1 |
30,5 |
85,4 |
5,2 |
5,3 |
|
|
7. |
им. Засядько, Ж, l1 |
- |
7,2 |
1,1 |
32,7 |
86,1 |
5,4 |
7,4 |
|
|
9. |
Комсомолец, Ж, m2 |
- |
10,4 |
1,2 |
29,5 |
87,4 |
5,3 |
6,1 |
|
|
10. |
им. Батова, К, k8 |
- |
6,0 |
3,0 |
22,4 |
88,6 |
4,8 |
3,6 |
|
|
11. |
Ясиновская-Глубокая, ОС, l6 |
- |
9,0 |
1,7 |
21,0 |
88,9 |
4,5 |
4,8 |
|
|
12. |
им. 60-летия Сов. Украины, Т, h8 |
- |
5,7 |
1,2 |
12,9 |
90,3 |
4,3 |
4,2 |
|
|
13. |
№ 20, А, h8 |
- |
7,0 |
1,7 |
4,0 |
93,6 |
2,3 |
2,4 |
|
|
14. |
Прогресс, А, h8 |
- |
7,3 |
1,5 |
4,2 |
93,7 |
1,9 |
2,9 |
|
|
Пиролиз образцов углей (фракция 0,2 мм) проводили в дериватографе Q-1500 системы Паулик-Паулик-Эрдеи в платиновом тигле с крышкой в атмосфере инертного газа аргона. Скорость нагрева составляла 10 оС/мин., навеска угля - 500 мг.
По анализу кривых ДТГ и ТГ определяли температуры максимумов и температурные интервалы основного термического разложения углей (таблица 2) [9]. Приведенные в таблице значения Тн. и Тк. характеризуют температуру начала и окончания эндоэффекта на кривой ДТГ, а Тмак. - температуру максимума этого эффекта. Из полученных данных следует, что температурный интервал основного термического разложения смещается в область более высоких температур с увеличением степени метаморфизма угля. Следовательно, благодаря различию структуры углей метаморфического ряда, данная температурная область является для каждого угля специфической, хотя и отображает общие закономерности протекания термических преобразований. Температурный максимум Тмак., соответствующий наибольшей скорости потери массы для всех изученных образцов, является той температурной точкой, которая отделяет температурную область преимущественной деструкции (нисходящая ветвь ДТГ) от области, где совместно с реакциями вторичной деструкции уже начинают протекать реакции конденсации (восходящая ветвь ДТГ). Таким образом, Тмак. является той характеристической температурной точкой, которая дает объективную информацию о структуре угля.
Таблица 2
Характеристика температурного интервала основного термического разложения угля |
|||||||
№ п/п |
Марка угля, пласт |
Сdaf, % |
Данные ДТГ, оС |
|
|||
Тн. |
Тмак. |
Тк. |
|
|
|||
1. |
Б |
65,4 |
200 |
375 |
450 |
|
|
2. |
Б |
66,8 |
195 |
370 |
445 |
|
|
3. |
Д, l2 |
76,2 |
305 |
405 |
480 |
|
|
4. |
Д, l4 |
79,0 |
320 |
405 |
480 |
|
|
5. |
Г, l3 |
82,2 |
295 |
425 |
550 |
|
|
6. |
Г, l1 |
82,7 |
335 |
425 |
510 |
|
|
7. |
Ж, k8 |
85,4 |
365 |
445 |
545 |
|
|
8. |
Ж, l1 |
86,1 |
360 |
440 |
550 |
|
|
9. |
Ж, m2 |
87,2 |
340 |
455 |
515 |
|
|
10. |
К, k8 |
88,6 |
400 |
475 |
565 |
|
|
11. |
ОС, l6 |
88,9 |
410 |
475 |
525 |
|
|
12. |
Т, h8 |
90,3 |
430 |
550 |
640 |
|
|
13. |
A, h8 |
93,6 |
520 |
625 |
750 |
|
|
14. |
А, h8 |
93,7 |
525 |
625 |
745 |
|
|
Для проверки этой гипотезы был проведен пиролиз углей ряда метаморфизма при температуре Тмак., подсчитан материальный баланс процесса и найдены уравнения регрессии, связывающие параметр Сdaf, температуру Тмак. с технологическим параметром угля - выходом смолы пиролиза (СМ).
Пиролиз углей осуществляли следующим образом. Навеску испытуемого образца массой 5-10 г (фракция 0-0,5 мм) загружали в реактор с неподвижным слоем (объем 20 см3), вакуумировали и нагревали со скоростью 25 оС/мин. до температуры максимума основного термического разложения (Тмак.), которая была ранее определена для каждого угля дериватографическим методом, и выдерживали при данной температуре в течение 10 мин. Это время было достаточным, чтобы при данных условиях эксперимента образующиеся первичные продукты пиролиза угля полностью вышли из зоны реакции [10]. Пиролиз проводили в закрытой системе без давления. По окончании эксперимента весовым методом определяли выход пирогенетической воды, смолы пиролиза, газа, твердого остатка и рассчитывали материальный баланс процесса. Расхождения между результатами двух параллельных определений не превышали по абсолютной величине для смолы пиролиза - ±0,5 %; пирогенетической воды - ±0,4 % и полукокса - ±0,7 % [5].
Количественная характеристика продуктов пиролиза для углей ряда метаморфизма приведена в таблице 3. Из таблицы следует, что в случае проведения пиролиза при Тмак, прослеживается взаимосвязь между Тмак. и выходом продуктов пиролиза для всех марок углей, включая угли высокой стадии метаморфизма (угли марки Т и А). Кроме того, наблюдается снижение выхода смолы пиролиза и общей степени конверсии в ряду метаморфизма от бурых углей до антрацитов.
Таблица 3
Материальный баланс пиролиза углей ряда метаморфизма, проведенного при температуре максимума основного термического разложения углей (Тмак, оС) |
|||||||
№ п/п |
Марка угля, пласт |
Продукты пиролиза, % |
Конверсия, % |
|
|||
Вода |
Смола |
Твердые |
Газы |
|
|||
2. |
Б |
4,6 |
6,9 |
80,6 |
7,9 |
19,4 |
|
1. |
Б |
5,0 |
5,3 |
81,2 |
8,5 |
18,7 |
|
3. |
Д, l2 |
7,2 |
3,4 |
83,9 |
5,5 |
16,1 |
|
4. |
Д, l4 |
3,4 |
3,0 |
86,8 |
6,8 |
13,8 |
|
5. |
Г, l3 |
3,6 |
3,1 |
89,2 |
4,1 |
10,7 |
|
6. |
Г, l1 |
2,9 |
2,8 |
89,1 |
5,2 |
10,9 |
|
7. |
Ж, k8 |
2,8 |
2,8 |
90,5 |
3,9 |
9,5 |
|
8. |
Ж, l1 |
2,0 |
2,3 |
93,7 |
2,0 |
6,3 |
|
9. |
Ж, m2 |
5,4 |
1,9 |
91,4 |
1,3 |
8,6 |
|
10. |
К, k8 |
5,5 |
1,2 |
87,6 |
5,7 |
12,4 |
|
11. |
ОС, l6 |
1,5 |
0,9 |
96,1 |
1,5 |
3,9 |
|
12. |
Т, h8 |
1,9 |
1,1 |
95,9 |
1,1 |
4,1 |
|
13. |
A, h8 |
1,3 |
0,1 |
97,6 |
1,0 |
2,4 |
|
14. |
А, h8 |
1,0 |
0,1 |
98,4 |
0,5 |
1,6 |
|
На рис. 1-2 показана зависимость выхода смолы пиролиза (СМ) от содержания углерода в угле Сdaf и от Тмак.. На рис. 3 приведена зависимость между содержанием углерода в угле Сdaf и величиной Тмак.
Рис. 1 Зависимость выхода смолы пиролиза (СМ) из углей ряда метаморфизма от содержания Сdaf, %
Для построения математической модели, адекватно описывающей взаимозависимость экспериментальных значений Cdaf, Tmax и CM [набор функций Cdaf = f (Tmax), Cdaf = f (CM), Tmax = f (Cdaf), Tmax = f (CM), CM = f (Tmax), а также CM = f (Cdaf)] нами использованы уравнения двухпараметровой регрессии, где аргументы представлены в виде квадратичных функций соответствующих параметров:
Y = Ao + A1X1 + A2X12
Известно, что подобными функциями, как правило, удается успешно описать широкий круг нелинейных физико-химических зависимостей [11-12] и построить достаточно надежные прогнозирующие модели.
В таблице 4 приведены результаты обработки полученных экспериментальных данных по этим уравнениям, где r – коэффициент корреляции, а So – стандартное отклонение. Рассчитанные средние величины ?ср. относительных расхождений между экспериментальными и вычисленными значениями функций Еср (%) = (Yэксп. - Yвыч)?100/ Yвыч для уравнений 1-6 составляли 1,8, 2,1, 19,9, 23,8, 4,1 и 3,7 % соответственно.
Зависимости Cdaf = f (Tmax), CM = f (Tmax) и, соответственно Tmax = f (Cdaf), Tmax = f (CM) - нелинейные (рис. 2-3), поэтому они могут быть описаны только квадратичными уравнениями 1, 3, 5, 6 (таблица 4). Полученные при этом коэффициенты корреляции во всех случаях достаточно высоки, но величины Еср существенно отличаются - от очень низких (уравнение 1) до достаточно больших (уравнение 4).
Таблица 4
Результаты регрессионного анализа изученных зависимостей. |
|||||||||
№ уравнения |
Параметры уравнения |
Коэффициенты уравнения |
r |
So |
|
||||
Y |
X1 |
X2 |
Ao |
A1 |
A2 |
|
|||
1 |
Cdaf |
Tmax |
Tmax2 |
-140,4 ± 25,5 |
0,83 ± 0,10 |
-0,0007 ± 0,0001 |
0,972 |
2,247</TD< tr> |
|
2 |
Cdaf |
CM |
CM2 |
94,5 ± 1,0 |
-4,51 ± 0,32 |
*) |
0,970 |
2,233 |
|
3 |
CM |
Tmax |
Tmax2 |
45,2 ± 7,6 |
-0,16 ± 0,03 |
0,00010 ± 0,00002 |
0,945 |
0,674 |
|
4 |
CM |
Cdaf |
Cdaf 2 |
19,9 ± 1,3 |
-0,21 ± 0,01 |
*) |
0,970 |
0,480 |
|
5 |
Tmax |
Cdaf |
Cdaf 2 |
3119 ± 550 |
-76,4 ± 13,9 |
0,53 ± 0,09 |
0,963 |
24,05 |
|
6 |
Tmax |
CM |
CM2 |
611 ± 17 |
-91,4 ± 11,8 |
8,47 ± 1,71 |
0,953 |
27,08 |
|
Примечание: *) - незначимый параметр (dА2 > A2)
Таким образом, рассматриваемые уравнения, позволяют определить с достаточно высокой степенью точности значения Tmax только лишь на основе данных элементного анализа углей (содержание углерода Cdaf). С другой стороны, исходя из значений Tmax, легко оценить как содержание углерода в конкретном угле, так и выход смолы пиролиза (СМ), не прибегая к соответствующему анализу. Однако, в последнем случае, значительные расхождения между экспериментальными и вычисленными значениями СМ не позволяют сделать это с высокой точностью.
Что касается зависимостей Cdaf = f (CM) и, соответственно, CM = f (Cdaf), то для изученного набора углей они хорошо описываются линейными однопараметровыми уравнениями 2 и 4 с высокими значениями коэффициента корреляции (таблица 4).
Таким образом, в настоящей работе предложен новый подход, позволяющий усовершенствовать метод определения выхода продуктов пиролиза из бурых, каменных углей и антрацитов. Метод предусматривает нагревание угля при температуре максимума основного термического разложения (Тмак.), которая предварительно определяется методом дериватографии. Показано, что температура Тмак. на кривой ДТГ является той температурной точкой, которую можно применять для оценки структурных и технологических характеристик углей. Данная температурная точка позволяет связывать воедино угли всего ряда метаморфизма и может быть использована при создании новой промышленной классификации углей [6].
1. Бутузова Л.Ф., Турчанина О.Н., Збыковский Е.И., Бутузов Г.Н. Статистические характеристики состава и свойств низкометаморфизованных углей Донбасса разных генетических типов // Углехимический журнал. – 2003. - №3-4. – С. 7-12
2. Осипов А.М., Лящук С.Н., Бойко З.В., Грищук С.В. Влияние состава сернистых углей Украины на выделение сероводорода при гидрогенизации // Химия твердого топлива.-2003.-№ 5.-С.68-75.
3. Иванов И.А., Гюльмалиев А.М, Гладун Т.Г., Гагарин С.Г., Головин Г.С. // Применение расширенного базиса структурных параметров для расчета энтальпии и атомизации соединений, моделирующих строение углей // Химия твердого топлива.-2003.-№ 1.-С.30-41.
4. Иванов И.А., Гагарин С.Г., Гюльмалиев А.М. Регрессивные зависимости технологических показателей углей от их элементного состава // Химия твердого топлива.-2003.-№ 5.-С.9-16.
5. ГОСТ 3168-93 (ИСО 647-74). Твердое топливо. Методы определения выхода продуктов полукоксования Взамен ГОСТ3168-75 СЭВ Введен 01.01.87.-М.: Изд-во стандартов, 1987.-12 с.
6. Авгушевич И.В., Броновец Т.М., Еремин И.В., Медведев А.В., Чурбаков В.Ф. Аналитическая химия и технический анализ угля.-М.: Недра, 1987.-336 с.
7. Камнева А.И., Платонов В.В. Теоретические основы химической технологи горючих ископаемых.-М.: Химия, 1990.-С. 155-156.
8. Тайц Е.М., Андреева И.А. Методы анализа и испытание углей.-М.: Недра, 1983.-301 с.
Скляр М.Г., Тютюнников Ю.Б. Химия твердых горючих ископаемых.-Киев: Выща школа, 1985.-247 с.
9. Аронов С.Г., Скляр М.Г., Тютюнников Ю.Б. Комплексная химико-технологическая переработка углей.-Киев: Техніка, 1968.-262 с.
10. Батунер Л.М., Позин М.Е. Математические методы в химической технике.-Л.: Химия, 1971.-824 с.
11. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. М.: Статистика, 1973. 561 с.
Патент 68874 А Україна, МПК С10В53/04 Спосіб визначення виходу продуктів піролізу з вугілля / Шевкопляс В.М., Галушко Л.Я. (UA).-№ 20031110274; Заявлено 14.11.03; Опубл. 16.08.04. –Бюл. № 8.-3 с.