шш^ш
УДК 620.9 DOI: 10.30724/1998-9903-2019-21-6-51-60
МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ПРОЦЕССА ТЕПЛООБМЕНА ЧАСТИЦ ГОРЮЧЕГО СЛАНЦА И ЗОЛЬНЫМ ТЕПЛОНОСИТЕЛЕМ В УСЛОВИЯХ ПОЛУКОКСОВАНИЯ ВО ВРАЩАЮЩЕМСЯ БАРАБАННОМ РЕАКТОРЕ УТТ
А.А. Селиванов1, П.А.Батраков 2, А.Н. Мракин1, И.А. Вдовенко1, О.В. Афанасьева 3
1 Саратовский государственный технический университет имени Гагарина Ю.А., г. Саратов, Россия. e-mail: [email protected] 2Омский государственный технически университет, г. Омск, Россия e-mail: [email protected] 3Казанский государственный энергетический университет, г. Казань, Россия
e-mail: [email protected]
Резюме: В статье рассматриваются вопросы моделирования процесса теплообмена в барабанном реакторе установки полукоксования горючего сланца с твердым теплоносителем между частицами золы и сланца в условиях постоянного выделения парогазовых продуктов. Выявлена доминирующая роль конвективной составляющей процесса и получены кривые прогрева частиц сланца с учетом эндотермического эффекта.
Установлено, что наиболее мелкая фракция топлива прогревается в течение 30-60 секунд и в дальнейшем служит источником теплоты для остальных фракций. В интервале 400-600 секунд завершается прогрев самой крупной фракции. Дальнейшее пребывание сланцезольной смеси в реакторе вызвано отставанием скорости процесса термического разложения от теплообмена, а снижение температуры на 20-250С наличием эндотермического эффекта разложения органического вещества.
Ключевые слова: горючий сланец, полукоксование, теплообмен, твердый теплоноситель, барабанный реактор.
Благодарности: Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта №17-08-00295 «А».
Для цитирования: Селиванов А.А., Батраков П.А., Мракин А.Н., Вдовенко И.А., Афанасьева О.В. Математическое описание процесса теплообмена частиц горючего сланца и зольным теплоносителем в условиях полукоксования во вращающемся барабанном реакторе УТТ // Известия высших учебных заведений. ПРОБЛЕМЫ ЭНЕРГЕТИКИ. 2019. Т. 21. № 6. С. 51-60. doi:10.30724/1998-9903-2019-21-6-51-60.
THE MATHEMATICAL DESCRIPTION OF PROCESS OF HEAT EXCHANGE PARTICLES OF OIL SHALE AND ASH THE COOLANT IN THE SEMI-COKING CONDITIONS IN A ROTATING DRUM REACTOR INSTALLATION WITH SOLID
HEAT CARRIER
AA Selivanov1, PA Batrakov2, AN Markin1, IA Vdovenko1, OV Afanaseva3
1Gagarin Yu.A. Saratov State Technical University, Saratov, Russia 2Omsk State Technical University, Omsk, Russia 3Kazan State Power Engineering University, Kazan, Russia
Abstract: The article deals with the modeling of heat transfer process in a drum reactor of a semi-coking unit of oil shale with a solid coolant between ash and shale particles in the conditions of constant separation of vapor-gas products. The dominant role of the convective component of the process is revealed and the curves of heating of shale particles are obtained taking into account the endothermic effect.
It is established that the smallest fraction of fuel warms up within 30-60 seconds and further serves as a source of heat for other fractions. In the interval of 400-600 seconds, the heating of the largest fraction is completed. The further stay of the shale mixture in the reactor is caused by the lagging rate of the thermal decomposition process from heat exchange, and the temperature decrease by 20-250C by the presence of the endothermic effect of the decomposition of organic matter.
Key words: oil shale, semi-coking, heat exchange, solid coolant, drum reactor.
For citation: Selivanov AA, Batrakov PA, Markin AN, Vdovenko IA, Afanaseva OV. The mathematical description of process of heat exchange particles of oil shale and ash the coolant in the semi-coking conditions in a rotating drum reactor installation with solid heat carrier. Power engineering: research, equipment, technology. 2019;21(6):51-60. doi:10.30724/1998-9903-2019-21-6-51-60.
Введение
Проектирование и строительство эффективных тепломассобменных устройств, работающих с оптимальными показателями, требует продолжительных и дорогостоящих исследований с первичной реализацией процессов на стендовых и полупромышленных аналогах. Важной составляющей является математическое описание, учитывающее основные закономерности теплового взаимодействия фаз.
Термическая переработка полидисперсных сыпучих топлив при постоянном перемешивании с инертным теплоносителем имеет достаточно широкое распространение в химической промышленности и энергетике. В особенности может быть выделена активно развивающаяся отрасль переработки горючего сланца на основе использования установок типа УТТ в Эстонии [1, 2].
При этом многочисленные вопросы остаются при детальном изучении механизма протекания процессов передачи теплоты в условиях взаимодействия сланца и зольного теплоносителя. Экспериментальные и теоретические исследования по протеканию тепломассообменных процессов при термической переработке мелкодисперсного топлива приведены в [1-5]. Обращает на себя внимание различие в условиях проведения физических экспериментов, конструктивные особенности аппаратов, а также особенности состава органической и минеральной части горючего сланца в зависимости от места добычи. Кроме того имеется ряд результатов, полученный аналитическим способом при принятых допущениях. Использование этих данным для инженерных и практических расчетов требует определенных уточнений и дополнений. Сложность в описании процесса теплового взаимодействия сланца и зольного теплоносителя также вносит полидисперсность и различие составов выделяющихся парогазовых продуктов с течением времени. Очень важна уже имеющаяся информация по взаимодействию топлива с твердым и газовым теплоносителем.
Таким образом, можно полагать об отсутствии на сегодняшний день математической модели процесса теплообмена сланца и зольного теплоносителя в барабанном реакторе УТТ с учетом особенностей состава самого топлива и продуктов термического взаимодействия, неоднородности в гранулометрическом составе и форме частиц, конструктивных особенностей аппарата с учетов вращения и др.
Теоретические основы и результаты
При описании закономерностей теплового взаимодействия частиц сланца и теплоносителя особый интерес представляют работы [1, 2, 6]. Неравномерность прогрева частиц по времени происходит по причине реагирования двух неоднородных компонентов (гранулометрический состав дисперсных материалов). На (рис.1 )показана схема взаимодействия частиц с учетом непрерывно выделяющихся продуктов полукоксования.
Исходя из (рис.1) перенос теплоты от частиц зольного теплоносителя (ЗТ) топливным частицам осуществляется через парогазовую прослойку (ПГС) с постоянно изменяющимися свойствами и составом. В свою очередь контактное взаимодействие частиц носит кратковременный характер ввиду вращения барабанного реактора.
Известно, прогрев твердых частиц определяется одновременным действием двух процессов: теплопередача от теплоносителя к поверхности и распространение ее вглубь частиц и зависит от размера, формы частиц, пористости, физических свойств теплоносителя и скорости движения газовой среды [7].
Рис. 1. Схема взаимодействия частиц сланца и золы в барабанном реакторе (БР)
Критерий Био позволит определить наличие внутреннего температурного градиента, который является определяющим для нахождения скорость проникновения теплоты от поверхности к центру частицы(1)
а • I
Ы =-, (1)
2 • 1
где I - определяющий размер, м; а -коэффициент теплоотдачи, 1 - коэффициент теплопроводности, Вт/(м-°С); размер
Как известно, сланец имеет плиточную структуру строения, что, в свою очередь, влияет на минимизацию расстояния от поверхности частиц до центра и значительно сокращает время прогрева.
При определении площади поверхности частиц неправильной формы используют эквивалентный диаметр йэ, м. Для учета формы частиц используется коэффициент формы Ф, определяющий отклонение формы реальных частиц от шара. В нашем случае важно учитывать тепловоспринимающую поверхность сланцевых и зольных частиц. При соотношениях сторон параллелепипеда фактор формы составит (2)
п • сС^ С^
ф. =
э = 0,8727-^-,
V 4а
(2)
где - отношение площадей поверхности сферы к площади поверхности реальной
С 2
частицы /-той фракции одинакового объема; э - эквивалентный диаметр частиц, м ;
а - длина ребра частиц правильной формы, м.
В табл. 1, 2 приведены значения с и Ф частиц геометрически правильной формы и осредненные значения ф для частиц неправильной формы.
а
Таблица 1
Значения ^ и Ф частиц различных форм
№ Форма частицы Эквивалентные диаметры Фактор формы, Ф,
dv,i и
1 2 3 4 5
1 - ( а х а х 2 а) 1,5632• а 1,7841 • а 0,8762
2 - ( а х а х 0, 5 а) 0,9847• а 1,1284• а 0,8727
3 - - 0,5759• а 0,8740 • а 0,6589
^ — а, высотой И = 0,1а
4 Пластина шестиугольная с ребром а, высотой И = 0,1 а 0,7917 • а 1,3583 • а 0,5828
Средние значения фактора формы для топливных частиц Таблица 2
№ Вид топлива Ф2 Ф1
1 сланец - размер фракции (^ = 2,5 4- 11,2 мм) 0,426 0,653
2 сланец - размер фракции^ = 5 4 25 мм) 0,46^0,68 0,755
3 сланец - размер фракции^ = 34 4 62,5 мм) 0,758 0,871
4 каменный уголь (^ = 6 4 11,25 мм) 0,536 0,732
5 металлургический кокс (^ = 6 4 11,25 мм) 0,403 0,635
В работе Симонова В.Ф. «Определение кинетических констант выделения газовых продуктов при термодеструкции сернистого горючего сланца» были рассмотрены фракционный состав сланца для полукоксования, где 84 % от общей массы частицы имеют размер < 5 мм. Данные фракции сланца как уже было сказано выше имеют плиточную форму: И = 0,54а, Ь = 1,43а, Ь = а. При этом важнейшим показателем является
И
минимальное расстояние от центра части к поверхности. = 0,94мм. В
2
ориентировочном интервале значений а (15^75 Вт/м2 °С) и известному значению , Хсл = 0,266, Вт/м°Свыполняется соотношение В- - 0,1_ 0,2. Данный факт подтверждают
авторы [1]. Таким образом, для мелкой и средней фракции сланца, занимаемой подавляющее большинство при переработке в схеме УТТ, можно предполагать отсутствие внутреннего температурного градиента.
Вторым основным процессом является перенос теплоты от поверхности частиц ЗТ через ПГС с постоянно изменяющимся составом к поверхности сланцевых частиц.
Данные о работе пилотных и промышленных реакторов типа УТТ, об особенностях движения сыпучего материала в БР, а также о структурно механических характеристиках Кашпирского сланца и сыпучих продуктах полукоксования приведены в работах [8,9] -табл.3, 4
Таблица 3
Технические характеристики установок УТТ [2]
Установка Реактор Частота вращения, об/мин. Время пребывания Скорость движения твердого материала, м/с Производитель ность по сланцу, т/ч
Отношение диаметров барабана (внутр./ внешн.) Длина, м твердого материала, мин. ПГС, с.
УТТ-200 1,64/2 4,5 1,83 20,6 18 0,35 6,34
УТТ-500 2,52/3 6,0 1,18 21,4 15 0,57 20,2
УТТ-3000 4,38/5 14,0 0,92 20,7 21,2 1,1 139,0
Таблица 4
Структурно механические характеристики Кашпирского сланца и его продуктов
Материал Угол естественного откоса, град. Плотность, кг/м3 Пористость, % по объему
динамический статический насыпная истинная кажущаяся
Сланец 45,4 42,1 819 2180 1830 15,3
КЗО (сланцевый полукокс) 42,7 39,5 823 2350 1760 23,7
Зольный теплоноситель 42,2 38,7 788 2480 1690 32,2
Зола уноса 41,8 37,9 762 2410 1650 33,4
На основе характеристик табл. 3, 4,данных о степени заполнения и конструктивных параметрах БР [9], а также результатов исследования процесса выделения ПГС [2] при известных расходных характеристиках по сланцу и золе получим соотношения для расчета объема каждого элементарного участка (рис. 2.) - (3) и времени пребывания сланцезольной смеси в БР - (4):
/ Л
V =
V Рсл рз У
(1 - Е)
1 п • Б,
2 4
2
б 7 3
• I, м
(3)
т=2V
(4)
где т?- - временной интервал расчета, с; V - свободный внутренний объем барабанного реактора, м3; £сл , ^ - расход сланца и ЗТ, кг/с; е - порозность движущегося слоя, согласно; р , р - кажущаяся плотность сланца и ЗТ, кг/м3; - внутренний диаметр БР,
П 3
м; 2 V - объем выделившихся компонентов ПГС, м /с. 1 1
Рис. 2. Распределение материальных потоков в процессе полукоксования в БР
Время пребывания объема в БР зависит от особенностей процесса выделения ПГС в сегменте толщиной I, м. (рис.2).
Лучистый теплообмен в установках подобного типа не превышает 7%, следовательно, влияние его по сравнению с конвективным незначительно [10].
Вторым важным обстоятельством является анализ контактной составляющей теплообмена. Типовые модели теплового взаимодействия твердых мелкодисперсных материалов имеют существенные различия с картиной взаимодействия сланца и теплоносителя в реакторе в условиях выделения ПГС [11]. В частности они предполагают продолжительный тепловой контакт частиц, который отсутствует в барабанном реакторе ввиду постоянного перемешивания и поступательного движения вперед, а также активной фильтрации продуктов полукоксования. Это дает право полагать о подавляющем влиянии конвективного переноса теплоты, а условия теплового взаимодействия могут быть представлены на основе результатов работ [1, 2] (Рис. 1).
П
1
Ми = 2 +1,3 • Рт'15 + 0,66 • Кс2т15 • Р^0'31
Таким образом, при движении смеси частиц сланца и ЗТ будет происходить фильтрация ПГС через разрыхленный слой частиц со значением порозности слоя превышающим рекомендуемые для движущегося плотного слоя. В работе [9] приведены данные по взаимодействию твердого теплоносителя (в форме шаров) и газа пропускаемого через засыпку. Значения критерия Рейнольдса 12< Яе < 2100. Уравнение, характеризующее значение критерия Нуссельта, имеет вид (5)
' Рг ^^
--(5)
V Ргст У
Похожие выражения приведены в работах А.П. Баскакова, В.Н. Тимофеева [10]. Скорость фильтрации газа через засыпку будет определяться с учетом постоянно растущего количества продуктов полукоксования.
V
" =—, (6)
ПГС / • 8
где / = • I - поверхностная площадь сегмента участка твердой фазы при движении
материала в реакторе (рис.2), м2; - суммарный объем ПГС на расчетном участке, м3/с. Число Рейнольдса вычисляется по формуле (7)
4 • w
Яе =-—, (7)
/о
где /0 = /(1 - 8) - удельная поверхность частиц в 1 м3 слоя, м2/м3; V - кинематический
коэффициент вязкости ПГС, м2/с.
Объемный расход ПГС будет складываться из значений удельных объемных расходов газов, выделяющихся на каждом участке за время в соответствии с температурой в заданном промежутке с учетом массового расхода (8), м3/с:
ТГ П лет уд
^лет = 2 8г- • иу , (8)
лет уд
где 8 - выход ПГС на /-том участке, кг/с; и - удельный объем ПГС на /-м участке (9),
м3/кг:
А уд Л
уд п и =2
V
. М
V з У
т
V т0 У
(9)
где из - удельный объем /-го компонента, м /кг.
Тепловой поток как со стороны зольных частиц к ПГС, так и от ПГС к сланцевым частицам будет определяться по классическому уравнению Ньютона-Рихмана (10).Частицы сланца (ьтой фракции), в заданный временной интервал Лт, способны воспринимать от парогазовой среды количество теплоты Л(?г, кДж, определяемое выражением
Ла = а, • Г, (С - / ) Лт •10-3, (10)
где, - соответственно коэффициенты теплоотдачи, Вт/(м2-°С); Б; -тепловоспринимающие поверхности, м2; . - температуры поверхностей частиц,°С,
каждой /-й фракции. С учетом данных по фракционному составу и геометрической форме, выражение для определения тепловоспринимающей поверхности примет вид (11)
6
р, = / 'т, =-;—тт •т,. (11)
' г г р, • •ф, г
1
Обсуждение
Расчеты по разработанному математическому описанию проводились с учетом результатов [8, 12, 15], где получены данные по константам выделения продуктов полукоксования для различных уровней прогрева материала с соответствующим выходом ПГС.
Для проверки правильности выбранных соотношений для описания особенностей теплового взаимодействия сланца и ЗТ с учетом ПГС необходимо сопоставление результатов численного и физического экспериментов. В качестве исходных данных были приняты результаты [13, 14].
Таблица 5
Сравнение результатов математической модели с экспериментальными результатами [10]
Показатели № опыта (численный эксперимент/ физический эксперимент)
1 2 3 4
Время пребывания смеси, с 1198/1188 1036/1144 512/799 769/954
Газы дезоксидации (СО2, Н2О, Н2Б), % масс 33,27/32,27 32,93/32,00 31,98/28,7 32,1/40,4
Легкие газы (Н2, СО, СН4, С^), % масс 9,69/8,1 9,65/11,7 9,4/16,9 9,51/35,3
Смола, % масс. 57,03/59,6 57,41/56,3 58,6/54,4 56,3/24,3
Среднеопытная погрешность, % 1,71 1,36 4,99 22,03
Общая погрешность, % 7,52
Согласно таблице 5, суммарная относительная погрешность по всем опытам находится в пределах 10%, а отдельно по опытам 25%, что подтверждает адекватность выбранных экспериментальных данных и эмпирических соотношений.
Данные рисунка 3 определяют интенсивность прогрева различных фракций сланца и снижение температуры теплоносителя, полученные в результате численного эксперимента на математической модели процесса полукоксования сланца и ЗТ (теплообмен осложненный выделением ПГС) [6,9,12].
0 200 400 600 800 1000 1200
Время полукоксования в реакторе, сек
Рис. 3 Кривые прогрева сланца полифракционного состава: 1 - ЗТd=1 мм; 2 - фракция 0-3 мм; 3 - фракция 3-7 мм; 4 - фракция 7-11 мм; 5 - фракция 11-15 мм
Как видно из рис. 3 имеет место перегрев (выше принятой температуры полукоксования 5000 С) мелких фракций, которые в процессе движения материала по барабанному реактору играют роль дополнительного источника теплоты для более крупных фракций. В интервале 400-600 секунд завершается прогрев самой крупной фракции. Кроме того впоследствии происходит понижение температуры сланцезольной смеси за счет учтенного в модели эндотермического эффекта разложения органической массы сланца с соответствующими кинетическими константами реакций [8].
57
Выводы:
1. Показано практическое отсутствие внутреннего температурного градиента у большинства частиц сланца при прогреве.
2. Выявлена доминирующая роль конвективного переноса теплоты от ЗТ к частицам сланца через ПГС с изменяющимся составом и свойствами.
3. В результате проведение численных экспериментов получены данные прогрева горючего сланца полифракционного состава. Показано, что наиболее мелкая фракция топлива прогревается в течение 30-60 секунд и в дальнейшем служит источником теплоты для остальных фракций.
4. Показано влияние эндотермического эффекта разложения органического вещества ГС, что приводит к снижению температуры потока на 20-25 °С.
Литература
1.Блохин, А.И.,. Зарецкий М.И, Стельмах Г.П., и др. Энерготехнологическая переработка топлив твердым теплоносителем. М.: Светлый стан, 2005. 336 с.
2.Волков Э.П. Энерготехнологическое использование сланца. М.: Издательский дом МЭИ, 2015. 480 с.
3.Иорудас, К.А., Стельмах Г.П. Тепло и массообмен в реакторе установок с твердым теплоносителем: Тезисы доклада 2 Всесоюзной научно-технической конференции. "Гидроупругость и долговечность конструкций энергетического оборудования" // Изд. Лит. ССР. Июль 1990. С. 213-214.
4.Хасхачих В.В. Экспериментальное исследование тепло-массообменных процессов переработки горючего сланца в установке с твердым теплоносителем (УТТ) // В сборнике: Состояние и перспективы развития электро- и теплотехнологии 18 Бенардосовские чтения: материалы международной научно-технической конференции. 2015. С. 220-223.
5.Печенегов, Ю.Я., Мракин А.Н. Влияние стефановского потока на теплообмен в системе "газ-твердая частица" при термохимической конверсии твердого топлива // Инженерно-физический журнал. 2017. Т. 90. № 5. С. 1152-1154.
6.Al-Ayed, OS. Distillation curves under the influence of temperature and particle size of Ellajjun oil shale // International Green Energy Conference. Canada: N. p., 2005. pp. 1-4.
7.AdnanAl-Harahsheha, Omar Al-Ayed, Moh'dAl-Harahsheh, Rajab Abu-El-Halawah. Heating rate effect on fractional yield and composition of oil retorted from El-lajjun oil shale // Journal of Analytical and Applied Pyrolysis, V.89. Issue 2. November 2010. pp. 239-243.
8.Симонов В.Ф., Селиванов А.А., Мракин А.Н., и др. Определение кинетических констант выделения газовых продуктов при термодеструкции сернистого горючего сланца // Вестник Саратовского государственного технического университета. 2015. Т. 4. №1 (81). С. 219-224.
9.Selivanov A.A., Batrakov P.A., Mrakin A.N., et al. Mathematical modeling of the carbon burn-out process from the shale semi-coke surface of aerofountain technological furnace in the solid heat-carrier oil shale semi-coking scheme // Procedia Engineering. 2016. 152. pp. 134-139.
10.A.P.Baskakov, B.Leckner. Radiative heat transfer in circulating fluidized bed furnaces / Powder Technology. 1997. V. 90. Issue 3. pp 213-218
11.A. I. Nafsun, F. Herz, E. Specht, H. et al. Thermal bed mixing in rotary drums for different operational parameters // Chemical Engineering Science. 2017. V.16. pp. 346-353.
12.Simonov V.F., Mrakin A.N., Selivanov A.A., et al., Development of an algorithm for calculating the parameters of a drum reactor for low-temperature carbonization of oil shale with ash-based heat-transfer agent // Chemical and Petroleum Engineering. 2017. V. 52. pp. 675-681.
13.Казаков, Е.И., Мамай В.И., Мишанин В.А., и др. Определение основных технологических параметров термической переработки сланца Перелюб-Благодатовского месторождения на стендовой установке с твердым теплоносителем (УТТ) // Межвузовский научный сборник СПИ - Саратов: СПИ, 1982. С. 101-105.
14.Volkov E., Stelmakh G. The stages of research on creating commercial units for processing oil shale fines. Development of the "Galoter" process in 1944-1999 // Oil Shale. 1999. V.16. N.2. pp. 161-185.
15.Симонов В.Ф. О влиянии скорости нагрева горючих сланцев на динамику выделения летучих продуктов и их качественный состав // Вестник СГТУ. 2011. Т. 4. № 1 (59). С. 181-185.
Авторы публикации
Селиванов Алексей Александрович - канд. техн. наук, доцент кафедры «Промышленная теплотехника», Саратовский государственный технический университет имени Гагарина Ю.А.
Батраков Пётр Андреевич - канд. техн. наук, доцент кафедры «Теплоэнергетика», Омский государственный технический университет.
Афанасьева Ольга Валерьевна - канд. техн. наук, доцент кафедры ЭМС, Казанский государственный энергетический университет.
Мракин Антон Николаевич - канд. техн. наук, доцент кафедры «Промышленная теплотехника», Саратовский государственный технический университет имени Гагарина Ю.А.
Вдовенко Иван Анатольевич - канд. техн. наук, доцент кафедры «Промышленная теплотехника», Саратовский государственный технический университет имени Гагарина Ю.А.
Referenses
1.Blokhin, AI, Zaretsky MI, Stelmach GP, et al. Energy technology processing of fuels with solid heat carrier. M.: Svetly Stan, 2005. 336 p.
2.Volkov EP. Energy technology use of slate. M : Publishing House MEI, 2015. 480 p.
3.Iorudas KA, Stelmakh GP. Heat and mass transfer in the reactor of plants with solid coolant: Abstracts of the 2 scientific and technical conf. "Hydroelasticity and durability ofpower equipment designs". Ed. Lit. SSR. July 1990.pp. 213-214.
4.Khaskhikhikh VV. An experimental study of heat and mass transfer processes for processing oil shale in a plant with a solid coolant (UTT). In the collection: Current State and Prospects for the Development of Electrical and Thermotechnology 18 Benardos Readings: Materials of the International Scientific and Technical Conference. 2015. pp. 220-223.
5.Pechenegov YuYa, Mrakin AN. The effect of the Stefan stream on heat transfer in the gas-solid particle system during thermochemical conversion of solid fuel. Engineering Physics Journal. 2017;90(5): 1152-1154.
6.Al-Ayed. Distillation curves under the influence of temperature and particle size of Ellajjun oil shale. International Green Energy Conference. Canada: 2005. pp. 1-4.
7.AdnanAl-Harahsheha, Omar Al-Ayed, Moh'dAl-Harahsheh, Rajab Abu-El-Halawah. Heating rate effect on fractional yield and composition of oil retorted from El-lajjun oil shale. Journal of Analytical and AppliedPyrolysis2010. V.89. Issue 2. pp.239-243.
8.Simonov VF, Selivanov AA, Mrakin AN, et al. Determination of the kinetic constants for the evolution of gaseous products during thermal decomposition of sour fuel. Bulletin of the Saratov State Technical University. 2015;4(1):(81):219-224.
9.Selivanov AA, Batrakov PA, Mrakin AN. Mathematical modeling of the carbon burn-out process from the shale semi-coke surface of aerofountain technological furnace in the solid heat-carrier oil shale semi-coking scheme . Procedia Engineering. 2016;152:134-139.
10.Baskakov AP, Leckner B. Radiative heat transfer in circulating fluidized bed furnaces. Powder Technology. 1997;90(3):213-218.
11. Nafsun A, Herz F, Specht E, et al. Thermal bed mixing in rotary drums for different operational parameters. Chemical Engineering Science. 2017;160 (16) :346-353.
12.Simonov VF, Mrakin AN, Selivanov AA, et al. Development of an algorithm for calculating the parameters of a drum reactor for low-temperature carbonization of oil shale with ash-based heat-transfer agent . Chemical and Petroleum Engineering. 2017;52:675-681.
13.Kazako EI, Mamai VI, mishanin VA, et al. Determination of main technological parameters of thermal processing of oil shale Perelyub-Blagodatovka field for bench set with solid heat carrier (UTT) interuniversity scientific collection of SLEEP - Saratov: SPI, 1982. pp. 101-105.
14.Volkov E, Stelmakh G. The stages of research on creating commercial units for processing oil shale fines. Development of the "Galoter" process in 1944-1999. Oil Shale. 1999;16(2):161-185. 4. А.
15.Simonov VF, Morev AA. About the influence of the heating rate of oil shale on the dynamics of the evolution of volatile products and their qualitative composition . Bulletin of SSTU. 2011;4(1): (59):181-185.
Authors of the publication Alexey A. Selivanov - Gagarin Yu.A. Saratov State Technical University, Saratov, Russia. Petr A Batrakov - Omsk State Technical University, Omsk, Russia. Olga V. Afanaseva - Kazan State Power Engineering University, Kazan, Russia. Anton N. Mrakin - Gagarin Yu.A. Saratov State Technical University, Saratov, Russia. Ivan A Vdovenko - Gagarin Yu.A. Saratov State Technical University, Saratov, Russia. Поступила в редакцию 15 июля 2019г.