CC BY
25
УДК 519.711.3, 004.942, 669.711
О.В.ТИТОВ, аспирант, [email protected]
Санкт-Петербургский государственный горный институт (технический университет)
O.V.TITOV, post-graduate student, [email protected] Saint Petersburg State Mining Institute (Technical University)
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОАНИЕ ОБЖИГА ИЗВЕСТКОВО-НЕФЕЛИНОВОЙ И ЦЕМЕНТНОЙ ШИХТЫ В ТРУБЧАТЫХ
ВРАЩАЮЩИХСЯ ПЕЧАХ
Статья посвящена разработке математических моделей процессов спекания извест-ково-нефелиновой и обжига цементной шихты в трубчатых вращающихся печах с учетом цепных завес. Цепные завесы устанавливаются во вращающихся печах для интенсификации процессов теплообмена между обжигаемым материалом и газовым потоком, сушки материала и для лучшего пылеулавливания. Рассмотрено влияние цепных завес на процессы теплообмена и пылеулавливания. Приведены уравнения, позволяющие описать эти процессы в математической модели трубчатой вращающейся печи. Полученная математическая модель процессов спекания и обжига реализована в программном комплексе ReactOp.
Ключевые слова: математическая модель; шихта; трубчатая печь; цепная завеса.
MATHEMATICAL MODELLING OF ROASTNG OF CALC-NEPHELINE AND CEMENT MIXTURE IN TUBULAR ROTARY
KILNS
The article is devoted to developing a mathematical models of the sintering processes of calc-nepheline and firing the cement mixture in a tubular rotary kilns based chain curtains. Chain curtains are installed in a rotary kilns for the intensification of heat exchange between the material and burn the gas flow, drying material, dust and carbon capture. The influence of chain curtains on the processes of heat transfer and dust separation are investigated. The developed equations allow to describe these processes in a mathematical model of a tubular rotary kiln. The resulting mathematical model of sintering and firing is implemented in software package ReactOp.
Key words: mathematical model, furnace feed, tubular rotary kiln, chain curtains.
Отличительной характеристикой таких процессов как спекание глиноземной и обжиг цементной шихты в трубчатых печах являются высокие температуры обжига, значительное время пребывания шихты в аппарате и высокие скорости движения газовой фазы в печи [3, 4]. В процессах спекания и обжига происходит изменение химического состава реагирующих фаз и значительное увеличение объема газовой фазы. Протекающие в печах спекания и обжига процессы являются, в основном эндотерми-
320
ческими и требуют подвода большого количества тепла. Основными эндотермическими процессами, происходящими в печи при вводе в нее шихты в процессах спекания и обжига при производстве глинозема, и цементного клинкера являются процессы декарбонизации, при которых происходит разложение карбонатов кальция и магния, а также процессы сушки и дегидратации кристаллогидратов. Для определения оптимальных условий проведения таких процессов необходима математическая модель ос-
новных физико-химических превращений. Наличие такой модели позволит описать количественно весь ход физико-химических превращений по длине аппарата, определить оптимальные условия для отдельных стадий процесса и разработать оптимальную технологическую схему и систему автоматизации для поддержания выбранного оптимального режима.
Основными теплообменными элементами трубчатых печей являются цепные завесы [1, 2].
Цепные завесы обычно устанавливают в загрузочной части печи, где температура газов не превышает 550 К (рис.1). В начале загрузочной части печи, где поступающая в печь шихта текучая, она налипает на поверхность теплообменных устройств и непосредственно соприкасается с горячими газами. По мере сушки шихта постепенно теряет способность налипать на элементы теплообменных устройств, которые в этом случае выполняют функции регенераторов, воспринимая тепло от газов и передавая его материалу. В районе расположения тепло-обменных устройств продвижение материала вдоль корпуса печи замедляется, вследствие чего повышается заполнение печи материалом, а также увеличивается поверхность теплообмена. В то же время улучшается перемешивание материала, в результате чего выравнивается его температура по сечению печи, увеличиваются разности и перепад температур между газами и обрабатываемым материалом, материалом и поверхностью теплообменников, что также способствует интенсификации теплообмена. Цепная завеса положительно влияет не только на теплообмен материала с горячими газами, но и на образование пыли и ее улавливание.
Величина конвективного теплообмена в зоне цепных завес описывается уравнением:
бкон - а кон Кр Д
(1)
где 2кон - конвективный теплообмен кДж/(м с-К); акон -коэффициент конвективного теплообмена, кДж/(м2-с-К); Кр - коэффициент плотности навески, %; Д£ - разность температур, К.
а
Рис. 1. Цепные завесы вращающихся печей:
а - конструкция части цепной завесы (1 - рядовое звено;
2 - концевое звено); б - узел крепления цепей к корпусу печи (1 - кронштейн; 2, 3 - болт с гайкой; 4 - корпус печи;
5 - концевое звено цепи; 6 - футеровка)
В соответствии с уравнением (1), конвективный теплообмен ^кон) должен увеличиваться при росте Кр, но в связи с тем, что при этом уменьшаются скорость и поток газа через цепи, то снижается коэффициент акон, что приводит к уменьшению Qкон. На участке пластичного материала при высоких Кр и тМ / тц может произойти полное перекрытие цепного сечения и тогда Qкон « 0.
Скорость газового потока оказывает значительное влияние на пылеунос из печи, который пропорционален скорости газа в 34 степени. Это приводит к значительному увеличению пылеуноса при росте скорости газа в пространстве над материалом.
В горячей части цепной завесы увеличение слоя материала в подцепном пространстве приводит к увеличению регенеративного теплообмена, но наряду с этим увеличивается и пылевынос за счет возрастания скорости газа. Поэтому в горячей части цепной зоны нельзя допускать чрезмерного увеличения слоя материала, необходимо перераспределение скорости в цепях и подцепном пространстве за счет снижения плотности навески. Для расчета основных параметров завесы со свободно висящими цепями рассчитывается максимально возможная длина цепи [1]:
К = ^п / п ,
(2)
Рис.2 Изменение содержания пыли по длине печи при спекании глиноземной шихты
Рис.3. Температурные профили газовой и твердой фазы при спекании глиноземной шихты с учетом цепной завесы
где DП - диаметр печи, м; п - число винтовых линий навески.
Общая поверхность цепей на участке длиной Lуч
Sобщ = naLуЧLц^ПП sin(90 - р)), (3)
где а - количество цепей, навешиваемых на один крепеж; Lуч, - длина участка, м, S1ПМ -
поверхность 1 м цепи, м2; h - шаг навески, м; р - угол наклона винтовой линии, град.
Поверхностный коэффициент плотности навески Кр рассчитывается как отношение площади поверхности цепей Sобщ к
площади поверхности корпуса на этом участке (м2/м2):
KF = So6rn, Zуч),
(4)
где Sобщ - общая поверхность цепей на рассматриваемом участке, м2; Оп - диаметр печи, м; Lуч - длина участка навески цепей, м.
Объемный коэффициент плотности навески Ку, рассчитывается как отношение площади поверхности цепей S общ к объему участка, на котором они навешены, м2/м3,
Kv = Яобщ/(^п2 Ly4)/4.
(5)
С учетом принятых допущений и закономерностей, описывающих процессы теплообмена и пылеуноса в присутствии цепных завес, математическая модель, описывающая изменение концентраций компонентов твердой и газовой фаз, будет иметь вид
322
dci 1 M
—L = w ;
dl us~~i dT 1 M
^ = —Y w. Q -B (T - T ) ± dl us ji jj Л s oJ
±Kf(Ts -Tg)±Kсер(Гs -Tg);
щ=i
dl u„
Kf T - Tg )± Kcep(Ts - Tg ); (6)
dC gaspili = GsZ + Z'
dl
100
K.F = Kcep , еСЛИ X ^ lcep K.F , еСЛИ X > ^ ,
cep iF '
сер :
где I - длина печи, м; - концентрация /-го компонента, кмоль/м3; и& - скорость материала, м/с; wij - скорость изменения концентрации /-го компонента в 7-ой реакции, кмоль/(м -с); Тё и Т - соответственно температуры газовой и твердой фазы, К; 2' - количество пыли, циркулирующее в печи, %; GS -
расход твердого материала в слое, кг/кг.
Для программной реализации разработанной математической модели был использовании программный комплекс ReаctOp [5] и система реакций, происходящих при получении компонентов клинкера и спекании шихты для производства глинозема [6].
Результаты моделирования процессов спекания глиноземной шихты и цементной шихты с учетом всех приведенных соотношений по теплообмену и пылеуносу представлены на рис.2, 3.
Выводы
1. Созданная математическая модель с учетом цепной завесы, позволяет прогнозировать состав конечного продукта спекания в зависимости от минералогического состава исходной шихты и условий работы печи.
2. Разработанная математическая модель позволяет отследить изменение составов и температурные профили твердой и газовой фаз по длине печи, выявить наиболее энергоемкие зоны печи, их смещение в зависимости от состава исходной шихты, прогнозировать точное количество тепловой энергии, необходимой для получения качественного спека из шихт различного состава.
3. На основе разработанной математической модели процесса спекания может быть определен оптимальный температурный профиль, обеспечивающий заданный конечный состав шихты.
ЛИТЕРАТУРА
1. Борисов И.Н. Некоторые особенности поведения шламов на цепях при высушивании / И.Н.Борисов, Л.С.Дурнева // Цемент и его применение. 2002. № 5. С.26-28.
2. Борисов И.Н. Пути повышения эффективности работы цепных завес вращающихся печей // Цемент и его применение. 2003. № 3. С. 17-20.
3. Сизяков В.М. Научно-технический прогресс в производстве глинозема / В.М.Сизяков, Н.С.Шморгуненко // Цветные металлы. М., 1981. № 8. С.43-45.
4. Сизяков В.М. Повышение качества белитового шлама при комплексной переработке нефелинов / В.М.Сизяков, П.В.Яшунин, А.И.Алексеев // Цветная металлургия. 1980. № 13. С.24-26.
5. Шариков Ю.В. Моделирование процессов и объектов в металлургии / Ю.В.Шариков, И.Н.Белоглазов, А.Ю. Фирсов. СПГГИ. СПб. 2006.
6. Шариков Ю.В. Кинетика химических превращений при термообработке шихты в трубчатых вращающихся печах при производстве глинозема и цементного клинкера // Цветные металлы. 2011. № 4.
REFERENCES
1. Borisov I.N. Some features of shludge behavior on chains at drying / I.N.Borisov, L.S.Durneva // Cement and its applying. 2002. № 5. P. 26-28.
2. Borisov I.N. Direction of effectivitis increasing of operabilitis chain curtains of rotary kilns / Cement and its applying. 2002. № 5. P. 17-20.
3. Sizyakov V.M. Scintific and technical prpogress in alumina production / Sizyakov V.M,. Shmorgunenko N.S. // Non-ferrous metals. Moscow, 1981. № 8. P.43-45.
4. Sizyakov V.M. Quality increasing of belite mixture at complex nepheline treatment / V.M.Sizyakov, P.V.Yakushin, A.I.Alekseev // Non-ferrous metallurgy. 1980. № 13. P.24-26.
5. Sharikov Y.V. Modelling of the processes and objects in metallurgy/ RPC SPMI. Saint-Petersburg, 2006.
6. Sharikov Y.V. Kinetics of the chemical transformation at thermotreatment of the furnace feed in tubular rotary kilns at production of alumina and cement clinkers / Y.V.Sharikov, O.V.Titov // Non-ferrous metals. 2011. № 4.
