Russian Federation
Razrabotana matematicheskaya model' sistemy aspiracii cementnoy mel'nicy, kotoraya mozhet byt' ispol'zovana dlya ee rascheta i pri sozdanii avtomatizirovannyh sistem upravleniya.
cementnaya mel'nica, aspiraciya, matematicheskaya model', upravlenie.
Введение. Технологическая аспирация цементных мельниц является важным условием увеличения их производительности, снижения удельных энергозатрат и повышения качества цемента [1]. В связи с модернизацией цементных заводов и переводом мельниц в замкнутый цикл повышаются требования к эффективности работы их аспирационных систем. При расчете и проектировании систем аспирации необходимо учитывать аэродинамическое сопротивление мельницы и других элементов аспирационного тракта, температуру и запыленность отсасываемого воздуха, подсосы наружного воздуха, характеристику аспирационного вентилятора и другие особенности процесса аспирации. Влияние этих факторов на режим аспирации исследовалось в ряде работ [2–6], однако методы комплексного аэродинамического расчета систем аспирации остаются недостаточно разработанными. В связи с этим подбор аспирационных вентиляторов и режимов их работы, как правило, ведется без должного обоснования, что приводит к перерасходу электроэнергии и снижению эффективности аспирации. Для разработки методов расчета и управления системой аспирации необходимо создание ее математической модели, связывающей между собой конструктивно–режимные параметры системы с управляющими факторами [7].
Основная часть. Аэродинамика основного элемента аспирационного тракта - мельницы (рис. 1) осложняется рядом особенностей, которые носят недетерминированный характер.
|
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
6 |
|
6 |
|
6 |
Рис. 1. Схема аспирационного тракта (1 – мельница, 2 – аспирационная шахта, 3 – группа циклонов,
4 – рукавный фильтр, 5 – аспирационный вентилятор, 6 –места подсосов наружного воздуха
Так, щели межкамерных перегородок, которые вносят существенный вклад в сопротивление мельницы, случайным образом перекрываются мелющими телами и забиваются измельчаемым материалом, а параметры аспирационного воздуха (температура, запыленность, влажность), которые также влияют на сопротивление мельницы, сложным и практически непредсказуемым образом меняются вдоль мельницы.
Поэтому зависимость сопротивления мельницы DP0 от массового расхода проходящего через нее аспирационного воздуха G0 целесообразно определять используя понятие эквивалентного отверстия [8], т.е. круглого отверстия в тонкой стенке, имеющего при тех же условиях сопротивление, равное сопротивлению мельницы:
. (1)
Здесь – плотность воздуха, просасываемого через мельницу, T0 – его средняя абсолютная температура, b = PM/R, P = 101325 Па, M = 0,029 кг/моль, R = 8,314 Дж/(моль∙K), b = 353,4 кг∙K/м3, Z0 – средняя по мельнице концентрация взвешенной пыли, A – площадь эквивалентного отверстия, m =0,65 – коэффициент расхода через него [9]. Величина A может быть найдена в результате аппроксимации полученной опытным путем зависимости DP0 от G0 функцией:
(2)
По найденному методом наименьших квадратов (МНК) значению параметра kМ площадь эквивалентного отверстия мельницы определяется соотношением:
(3)
Для мельниц размерами 3´14 м
A = 0,15…0,19 м2 [9]
Соотношение (1) для аэродинамического сопротивления мельницы может быть представлено в виде:
, (4)
где
(5)
В производственных условиях аспирация цементных мельниц сопровождается значительными подсосами наружного воздуха в местах соединения элементов аспирационного тракта. Так, перед аспирационной шахтой расход подсасываемого воздуха может составлять , перед циклонами – , перед фильтром и перед вентилятором . Суммарный подсос может превышать расход просасываемого через мельницу воздуха в несколько раз [10]. Количество наружного воздуха, подсасываемого через неплотности i–го сопряжения аспирационного тракта зависит от площади этих неплотностей Fi, разрежения в месте их расположения DPi и плотности наружного воздуха :
, (6)
где T – температура наружного воздуха, – коэффициент местного сопротивления щелевидного отверстия при наличии проходящего потока [11].
Подсосы приводят не только к увеличению расхода аспирационного воздуха, но и к снижению его температуры и запыленности. Влияние изменений этих параметров необходимо учитывать при определении сопротивлений элементов аспирационного тракта. Температура и объемный расход пылевоздушной смеси в аспирационной шахте определяется соотношениями:
,
(7)
Запыленный воздушный поток, поступающий в аспирационную шахту, дополнительно захватывает с собой часть выгружаемого из мельницы цемента, поэтому его начальная запыленность Z1 значительно превосходит концентрацию пыли на выходе из мельницы и может достигать 330…770 г/м3 [12]. Концентрация пыли на выходе аспирационной шахты определяется массовым расходом пылеуноса Gу, который зависит от расхода аспирационного воздуха G0 и достигает для мельниц открытого цикла измельчения от 2 до 5% производительности мельницы Gм [13]. Полагая эту зависимость в интервале 1 < G0 < 9 кг/с линейной, получим
(8)
где – относительная доля пылеуноса, параметры kу и B определяются в результате обработки опытных данных.
Эффективность пылеосаждения в аспирационной шахте, происходящего в результате интенсивной коагуляции и седиментации укрупненных частиц, достигает 90%. В этом случае входная Z1 и выходная Z2 концентрации цементной пыли в аспирационной шахте связаны соотношением Z1 = 9Z2, а для средней концентрации пыли получим: , где
(9)
Аэродинамическое сопротивление аспирационной шахты DPa обусловлено поворотом и расширением аспирационного потока в разгрузочном устройстве мельницы, весом выносимой пыли, разностью температур аспирационного и наружного воздуха, а также сужением потока аспирационного воздуха при его входе в группу циклонов
(10)
Здесь Fp – площадь свободного поперечного сечения разгрузочного устройства, g = 9,81 м/c2 – ускорение силы тяжести, H, Fa – высота и площадь поперечного сечения аспирационной шахты; и – коэффициенты местного сопротивления переходов из разгрузочного устройства в аспирационную шахту и из нее в группу циклонов.
В качестве первой ступени очистки аспирационного воздуха чаще всего используются группы циклонов ЦН–15 или ЦН–15у диаметром мм. Коэффициент сопротивления группы таких циклонов с учетом запыленности аспирационного воздуха равен , а эффективность очистки достигает значения , если скорость пылевоздушной смеси, отнесенная к площади поперечного сечения цилиндрической части циклонов (условная скорость) близка к оптимальной [11]:
, (11)
где – сумарная площадь поперечных сечений циклонов, n = 2..4 – их число.
При заданном числе циклонов их диаметр определяются соотношением:
(12)
Сопротивление группы циклонов определяется формулой:
(13)
Здесь – средняя плотность пылевоздушной смеси в циклонах, , Tц – плотность воздуха и его температура, Zцс – средняя концентрация пыли в циклонах.
Соотношения для определения Tц и Zсц выводятся из уравнений теплового и материального баланса:
(14)
(15)
С учетом условия (10) и соотношений (12, 13) формула для сопротивления группы циклонов принимает вид:
(16)
Для второй ступени очистки аспирационного воздуха используются рукавные фильтры серий ФРКИ – 360 с производительностью до 32400 м3/ч.
Сопротивление рукавного фильтра может составлять Па, а эффективность достигает значений .
Температура и запыленность аспирационного воздуха на входе в рукавный фильтр, которые снижаются из–за подсоса наружного воздуха в количестве G3, могут быть найдены по формулам:
(17)
(18)
Заключительным элементом аспирационного тракта является вентилятор. Расход, температура и запыленность воздуха, выбрасываемого вентилятором в атмосферу с учетом подсасываемого им наружного воздуха G4, можно рассчитать по формулам:
, (19)
, (20)
, (21)
(22)
где – полный расход подсасываемого воздуха.
Вентилятор может использоваться в качестве управляющего элемента системы аспирации, поскольку при изменении частоты его вращения меняется объем подсосов и температура аспирационного воздуха, что наряду с частотой вращения влияет на режим работы вентилятора. Поэтому математическая модель системы управления аспирацией цементной мельницы должна включать в себя аналитическое выражение характеристики вентилятора с учетом ее зависимости от частоты вращения и температуры отсасываемого воздуха [14]. В качестве такого выражения предлагается следующая аппроксимация характеристик вентиляторов:
,(23)
где ns, Ts – частота вращения и абсолютная температура перемещаемого воздуха при стендовых испытаниях вентилятора, n и T – рабочие значения этих параметров, значения коэффициентов A,B и C определяются в результате статистической обработки опытных характеристик вентиляторов [15].
Так, аппроксимации характеристик дымососов Д–12 и Д–15.5, полученных при ns = 970 об/мин и Ts = 473 К (ts = 200 °C), имеют вид:
(24)
(25)
Потребляемую мощность вентилятора можно найти по формуле
кВт, (26)
где – коэффициент полезного действия (КПД) вентилятора, зависящий от его производительности. Эти зависимости для дымососов Д–12 и Д–15,5 можно аппроксимировать в следующем виде:
(27)
(28)
Систему аспирации цементной мельницы можно представить в виде сети, состоящей из пяти вложенных друг в друга и условно замкнутых через внешнюю среду контуров (рис. 2).
В соответствии со вторым законом Кирхгофа система уравнений, описывающая аспирационную сеть, имеет вид:
(29)
где , i = 1, 2, 3, 4 – сопротивления неплотностей i-го сопряжения элементов аспирационного тракта:
, (30)
где Gi – массовый расход наружного воздуха, проходящего через неплотности i-го сопряжения элементов аспирационного тракта,
(31)
Условие сохранения массового расхода аспирационного воздуха учитывалось при выводе соотношений для сопротивлений отдельных элементов системы аспирации.
Вычитая из первого уравнения системы (29) второе, а со второго – третье и т.д. приведем ее к следующему виду:
(32)
|
DP0 |
|
DPa |
|
DPц |
|
DPф |
|
Pv |
|
G1 |
|
G2 |
|
G3 |
|
G4 |
|
G0 |
|
Gв |
|
G0+G1 |
|
DP2 |
|
DP3 |
|
DP4 |
|
DP1 |
|
G0+G1+G2 |
|
G0+G1+G2+G3 |
|
G0 |
|
G0 |
Рис. 2. Схема распределения воздушных ( ) и пылевоздушных ( ) потоков в системе аспирации цементной мельницы
Для численного решения системы уравнений (32) приведем ее к виду:
(33)
Уравнения (31) связывают вычисляемые массовые расходы воздуха, проходящего через мельницу G0 и через неплотности аспирационного тракта (подсосы G1–G4), с конструктивно–технологическими параметрами мельницы, аспирационного тракта и температурой внешнего воздуха. Эти уравнения вместе с соотношениями для сопротивлений элементов аспирационного тракта (4, 10, 15) и характеристикой вентилятора (22) составляют математическую модель системы аспирации. Путем численной реализации этой модели исследовалось влияние управляющих воздействий (изменения частоты вращения вентилятора, режима регенерации рукавного фильтра и других) с учетом возмущающих воздействий (изменений площади неплотностей и температуры воздуха в мельнице и снаружи) на характеристику аспирационного тракта и основные показатели аспирационного процесса: скорость воздушного потока внутри мельницы Ua, относительный уровень подсосов Kп = 100Gп/G0, % , энергозатраты на работу вентилятора, температура запыленного воздуха, поступающего в циклоны и фильтр, запыленность аспирационных выбросов. Рассмотрим численную реализацию математической модели системы аспирации на примере мельницы открытого цикла измельчения размерами 3´14 м. Производительность мельницы GМ = 43 т/ч = 12 кг/c, степень ее заполнения j = 0.27, площадь эквивалентного отверстия A = 0,17 м2, площадь свободного сечения мельницы Fсв = 4,66 м2, температура и запыленность аспирационного воздуха на выходе из мельницы T0 = 110 °C = 383 K, Z0 = 0,4 кг/м3, ky = 0,01, B = -0,01; концентрация пыли на входе в аспирационную шахту Z1=0,7 кг/м3 , а ее высота H = 14 м, число циклонов n = 2, их эффективность ; сопротивление рукавного фильтра DPф = 1400 Па, его эффективность hф = 0.998, в качестве аспирационных вентиляторов рассматривались дымососы Д–12 и Д–15.5.
Наименее энергозатратная аспирация цементных мельниц возможна при полном отсутствии подсосов. Однако, характеристика герметизированного аспирационного тракта (рис. 3, 1) показывает , что поддерживать такой режим аспирации (Ua = 0,7 м/c, G0= 3 кг/c, Q0 = 11,7 тыс. м3/ч, PV = 3100 Па, Nв = 14,4 кВт при hв = 0,7) можно лишь с помощью вентиляторов высокого давления (разрежения).
В реальных условиях полное устранение подсосов практически невозможно, поскольку их влияние проявляется уже при площади всех щелей и неплотностей F » 0,025 м2. При наличии подсосов характеристика аспирационного тракта становится более пологой (рис. 3, 2), что позволяет использовать в качестве аспирационных вентиляторы среднего давления.
|
4000 |
|
3500 |
|
3000 |
|
2500 |
|
5 |
|
10 |
|
15 |
|
Р, Па |
|
Gв, кг/c |
|
4 |
|
1 |
|
2 |
|
3 |
Рис. 3. Характеристики аспирационного тракта (1–без подсосов, 2 – с подсосами, F = 0,076) и вентиляторов
(3 – Д-12, n = 970 об/мин, t = 75°C, 4 – Д-15,5, n = 970 об/мин, t = 77 ° C)
Так, в теплый период года ( t=25°C) при общей площади неплотностей сопряжений аспирационного тракта F = 0.076 м2 (F1=0,05, F2=0,01, F3=0,003 и F4 = 0,014 м2 ) с помощью дымососа Д-12 можно поддерживать режим аспирации (рис. 3) с параметрами: Ua= 0,7 м/с, Gв = 5,17 кг/c, Qв = 18,3 тыс. м3/ч, Nв = 25 кВт, hв =0,64, коэффициент пылеуноса p=0,02, относительный объем подсосов Kп = 71%, температура воздуха, поступающего в циклоны и рукавный фильтр tц = 83 °С, tф = 82 °С, запыленность аспирационных выбросов Zв = 19 мг/м3. При понижении температуры наружного воздуха до 5 °С (холодный период года) уровень подсосов сохраняется, а интенсивность вентиляции мельницы и потребляемая вентилятором мощность повышаются (Uа =0,76 м/c, Nв = 27,2 кВт). Снижение скорости вращения вентилятора на 15 об/мин. позволяет восстановить значения этих параметров. Рассмотренный режим работы системы аспирации для реальных производственных условий по–видимому является оптимальным. Широко распространенное применение в аналогичных условиях более мощного дымососа Д-15,5, сопровождаемое значительным увеличением объема подсасываемого воздуха и потребляемой вентилятором мощности (рис. 3), прежде всего обусловлено неудовлетворительным состоянием аспирационного тракта.
Выводы. 1. Необходимым условием снижение энергозатрат на аспирацию цементных мельниц является герметизация сопряжений элементов аспирационного тракта.
2. Для реализации энергосберегающего режима работы системы аспирации необходимо регулирование характеристики аспирационного вентилятора.
3. Разработанная математическая модель может быть использована для обоснованного подбора аспирационных вентиляторов, а также при разработке автоматизированных систем управления аспирацией цементных мельниц.
1. Bogdanov V.S., Nesmeyanov N.P., Pirockiy V.Z., Morozov A.I. Mehanicheskoe oborudovanie predpriyatiy stroitel'nyh materialov. Belgorod: Izd-vo BelGTASM, 1998. 180 s.
2. Pirockiy V.Z., Bogdanov V.S., Sevost'yanov V.S. Aspiraciya cementnyh mel'nic. M.: VNIIESM, 1984. Vyp. 1. S. 1-51.
3. Kryhtin G.S., Zharko V.I. Aerodinamika cementnyh mel'nic. Nauchnye soobscheniya NIICementa. Vyp. 23. M. Stroyizdat, 1968. S.42-46.
4. Minko V.A., Logachev I.N. Raschet aspi-racii i sistem CPU. Belgord: Izd-vo BelG-TASM, 1994. 53 s.
5. Shaptala V.G. Matematicheskoe modelirovanie v prikladnyh zadachah mehaniki dvuhfaznyh potokov. Belgorod: Izd-vo BelGTASM. 1996. 102 S.
6. Logachev I.N., Logachev K.I. Aerodina-micheskie osnovy aspiracii. S.-Pb.: Himiya, 2005. 659s.
7. Rubanov V.G., Filatov A.G. Modeliro-vanie sistem. Belgorod: Izd-vo BGTU, 2006. 349 s.
8. Rudnichnaya ventilyaciya. Spravochnik. Pod red. K.Z. Ushakova. M.: Nedra, 1988. 440 s.
9. Tihovidov B.D., Balera N.D. Vybor kriteriya aspiracionnyh kachestv trubnoy sharovoy mel'nicy. Cement. 1983. № 11. S. 14-16.
10. Satarin V.I., Perli S.V. Dvizhenie i obespylivanie gazov v cementnom proizvod-stve. M. Gosstroyizdat. 1960. 265s.
11. Idel'chik I.E. Spravochnik po gidravlicheskim soprotivleniyam. M., Mashinostroenie, 1975. 559 s.
12. Lapshin A.B., Koziko G.G., Musatyan S.A. Obespylivanie aspiracionnogo vozduha ce-mentnyh mel'nic. Cement. 1989. № 5. S. 20-21.
13. Aspiraciya cementnyh mel'nic [Elektronnyy resurs] rezhim dostupa: http://ecologenergy.com (data obrascheniya 29.01.2017)
14. Rubanov V.G. Teoriya avtomaticheskogo upravleniya (matematicheskie modeli, analiz i sintez lineynyh sistem). Belgorod: Izd-vo BGTU. 2009. 199 s.
15. Aerodinamicheskie harakteristiki dymososov VD i D. [Elektronnyy resurs] rezhim dostupa: http://ventilator.spb.ru (data obrascheniya 29.01.2017).
