Модель оптимизации параметров размола шаровой барабанной мельницы ШК-32 Текст научной статьи по специальности «Физика»

 Энергетика

S5S Power Engineering

Оригинальная статья / Original article УДК: 536.24

DOI: 10.21285/1814-3520-2016-10-128-134

МОДЕЛЬ ОПТИМИЗАЦИИ ПАРАМЕТРОВ РАЗМОЛА ШАРОВОЙ БАРАБАННОЙ МЕЛЬНИЦЫ ШК-32

© С.Д. Иванов1, А.Н. Кудряшов2, В.В. Ощепков3

Иркутский национальный исследовательский технический университет, 664074, Россия, г. Иркутск, ул. Лермонтова, 83.

РЕЗЮМЕ. ЦЕЛЬ. Проведение исследования влияния диаметра угольных частиц на расход электроэнергии на пылеприготовление шаровой барабанной мельницы (ШБМ) ШК-32 с целью оптимизации ее работы. МЕТОДЫ. Использованы численные алгоритмы многопараметрической оптимизации транспорта пылевоздушной смеси в ШБМ. РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ. Получена целевая функция для минимизации энергозатрат на пылеприготовление в ШБМ ШК-32. Данная функция показывает возможность снижения расхода электроэнергии на 15%. ВЫВОДЫ. Важным фактором снижения энергозатрат на размол и пневмотранспорт является выбор оптимальной степени заполнения мельницы шарами.

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

Формат цитирования: Иванов С.Д., Кудряшов А.Н., Ощепков В.В. Модель оптимизации параметров размола шаровой барабанной мельницы ШК-32 // Вестник Иркутского государственного технического университета. 2016. Т. 20. № 10. С. 128-134. DOI: 10.21285/1814-3520-2016-10-128-134

OPTIMIZATION MODEL OF BC-32 BALL PULVERIZER MILL GRINDING PARAMETERS S.D. Ivanov, A. N. Kudryashov, V.V. Oshchepkov

Irkutsk National Research Technical University, 83, Lermontov St., Irkutsk, 664074, Russia.

ABSTRACT. The PURPOSE of this paper is to study the effect of the coal particle diameter on power consumption for pulverization by the ball pulverizer mill BC-32 (ball conical) in order to optimize its operation. METHODS. Numerical algorithms of multi-parameter optimization of coal-air mixture transport in the ball pulverizer mill are used. RESULTS AND THEIR DISCUSSION. A target function is obtained to minimize energy consumption for pulverization in the BC-32 ball pulverizer mill. This function shows the possibility to reduce energy consumption by 15%. CONCLUSIONS. The choice of an optimal degree of filling the mill with balls is an important factor in the reduction of energy consumption for grinding and pneumatic transport.

Keywords: optimization, ball pulverizer mill, energy consumption for coal pulverization, pneumatic transport power, grinding power, balls filling degree, average diameter of coal dust particles

For citation: Ivanov S.D., Kudryashov A.N., Oshchepkov V.V. Optimization model of BC-32 ball pulverizer mill grinding parameters. Proceedings of Irkutsk State Technical University. 2016. vol. 20, no.10. pp. 128-134. (In Russian) DOI: 10.21285/1814-3520-2016-10-128-134

Введение

Технологические процессы измель- Так, на ТЭС расход электрической энергии чения угля с помощью шаровых барабан- на пылеприготовление составляет около ных мельниц (ШБМ) довольно энергоемки. 25% от общего расхода электроэнергии на

1

Иванов Сергей Дмитриевич, кандидат физико-математических наук, доцент кафедры теплоэнергетики Института энергетики, e-mail: ivsd55@yandex.ru

Ivanov Sergey, Candidate of Physical and Mathematical sciences, Associate Professor of the Department of Heat-Power Engineering of the Institute of Power Engineering, e-mail: ivsd55@yandex.ru

2Кудряшов Александр Николаевич, кандидат технических наук, доцент, заведующий кафедрой теплоэнергетики Института энергетики, e-mail: kan@istu.irk.ru

Kudryashov Alexander, Candidate of technical sciences, Associate Professor, Head of the Department of Heat-Power Engineering of the Institute of Power Engineering, e-mail: kan@istu.irk.ru

3Ощепков Василий Владимирович, аспирант, e-mail: voui@mail.ru Oshchepkov Vasily, Postgraduate, e-mail: voui@mail.ru

собственные нужды, или около 2% от общей выработки электрической энергии. Поэтому оптимизация рабочих параметров ШБМ с целью снижения этих расходов актуальна.

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

вынос пыли из барабана уменьшается, происходит переизмельчение угля и понижение производительности мельницы. С усилением скорости потока воздуха производительность мельницы растет, но при этом повышаются затраты на вентиляцию вследствие увеличения среднего размера уносимых частиц угля [1].

Материал, цель и методы исследования

В данной статье рассматривается метод оптимизации рабочих параметров мельницы ШК-32 с учетом среднего диаметра частиц угольной пыли.

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

N (й)

Э = Э + Э = _1 +

Э Эр ^ ЭШБМ ^

B

N

+ 'ШБМ (d) = Ф(d),

B

где Мр - мощность, затрачиваемая на размол угля в ШБМ; МШБМ - мощность, затрачиваемая мельничным вентилятором на продувку пыли; й - средний размер угольной частицы после размола; В - производительность мельницы по сухой пыли после ШБМ.

По экспериментальным данным мощности Мб, потребляемой на вращение барабана [2], можно найти мощность Мр по уравнению из [3]

Np = N6 - QMex - N60,

(1)

где = 0,7^ - тепло, выделяющееся в результате работы мелющих органов;

1,86 _

N0 =-ВЬпб8б - мощность на вращение пустого барабана (без мелющих эле-

ментов); г]б - КПД привода; О - эквивалентный диаметр мельницы; ^ - длина мельницы; п - скорость вращения барабана; Бб - толщина стенки барабана.

Мощность, затрачиваемую на пневмотранспорт угольной пыли через ШБМ, также можно вычислить по экспериментально измеренным данным по формуле из [3]:

+ 1 АРшБМ (й)

N„

( d ) = B-

ß

РЛм

см'1МВ

где л = Б/Ов - массовая расходная концентрация пыли; - расход воздуха после сушки топлива в ШБМ [2]; АрШБМ - перепад давления в ШБМ [2];

рсм=арч+ (1 -а)рв - плотность смеси воздуха и пыли; а - объемная доля частиц в воздухе; рч - плотность частиц угольной

пыли; р - плотность воздуха.

Средний диаметр частиц угольной пыли был вычислен по уравнению Розина - Раммлера:

R = е~

(2)

где ^ - остаток на каком-либо сите; Ь и п - опытные коэффициенты, характеризующие соответственно тонкость помола и равномерность зернового состава (показатель полидисперсности), постоянные для данного топлива и данного метода размола величины (определяются рассевом навесок пыли на двух ситах, чаще используют сита

Ьх

с размером ячеек 90 и 200 мкм); x - размер частиц.

Экспериментальные значения ^до и R2oo после сепаратора представлены в [2], по ним определили соответствующие значения параметров Ь и п.

Продифференцировав уравнение (2) по размерам частиц x, получили функцию плотности распределения частиц пыли по размерам:

F(b,n,x) = - dR/ dx = nbx

n—1 - bxn

Средний диаметр частиц

xmax

d = J xF(b, n, x)dx

где xmax - максимальный диаметр частиц (в рассматриваемом случае 0,0005 м).

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

Рис. 1. Экспериментальная зависимость полного удельного расхода электроэнергии на пылеприготовление от среднего диаметра частиц угля Fig. 1. Experimental dependence of the total specific energy consumption for coal pulverization

on the average coal particle diameter

Результаты исследования и их обсуждение

Как видно из рис. 1, имеется минимум энергозатрат при определенном диаметре частиц. Далее рассмотрим пути его теоретического снижения.

Теоретическая мощность вращения барабана ЫТб определяется по формуле из [3]:

NT = 11 (0,127D3Ln6yiy6 + 1,86DLn6S6), (3)

где Y - насыпной удельный вес шаров; уб - степень заполнения ШБМ шарами.

На рис. 2 представлено сравнение экспериментальной мощности ЫЭб, потребляемой мельницей, и теоретической Ытб,

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

Теоретическая мощность размола Ыб с учетом выражений (1) и (3) принимает вид

NT = NT — ОТ — N

N р N б Омех N 60 ■

(4)

0

1400

1000

600

200

: Nß, кВт • A * *- • • A • -a- •-m 1 • Ï * ♦ A V3 N6

..............A «.......... Ä * A VT

A A A A A

d: M KM

50

60

70

80

90

100

110

Рис. 2. Сравнение экспериментальной (•, сплошная прямая) и теоретической (▲, пунктирная прямая) мощностей вращения барабана мельницы Fig. 2. Comparison of experimental (•, full line) and theoretical (▲, dotted line) powers of mill cylinder rotation

Чтобы получить целевую функцию, необходимо найти связь мощности ИТр со

средним диаметром частиц угольной пыли.

Полезная мощность размола идет на образование новых поверхностей пыли в результате удара (трением пренебрегаем) [4]:

Я = ^ ± ф,

пов ъ г-г /jr '

2Ерч г=0

(5)

где а - предел прочности на сжатие частицы угля; Е - модуль упругости материа-

ла (модуль Юнга); п = 1п \ — | - число

1п 2 ^ ё )

этапов измельчения от начального диаметра б0 до конечного ё ; ф = 0,98 - доля энергии, приходящаяся на измельчение ударом.

Приняв число п = 5-30, можно аппроксимировать сумму в формуле (5).

Данная аппроксимация представлена на рис. 3 и имеет следующий функциональный вид:

ф( d ) =-1,463 + 2,294 Г ¿0

V d

>0,7

Рис. 3. Аппроксимация ступенчатой функции в (6) непрерывной зависимостью Fig. 3. Step function approximation in (6) by a continuous dependence

Максимальный коэффициент полезного действия процесса размола, учитывая уравнения (4) и (5), можно записать так:

max

=

N

NТр°

0,0012,

(6)

Т ,опт

где Ы - теоретическая оптимальная мощность размола при оптимальной степени загрузки = 0,2 [3].

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

Np(d) =

Ba2

In 2 I d

2 Eprf

(7)

i =0

Теоретическую мощность, затрачиваемую на пневмотранспорт угольной пыли через ШБМ, можно найти из теоретического перепада давления на ШБМ крТптм:

N Т

(d)

= B

ß+1 лртшш (d)

ß

РамЛм

(8)

см' 1МВ

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

АрШШМ (d) =APmp(d) + Дринж (d)■

(9)

В [5] получены критериальные уравнения для этих составляющих перепада давлений в ШБМ. В несколько модернизированном виде перепад давлений, связанный с транспортом газовзвеси, равен:

. 1,214 • 101

APmV (d) = ---Х

Re

1,5

ß f D J pewe (d)

1 + ßV L

(10)

40е (й)

где ^= —- скорость воздуха;

Re„ =

nD 2p 4Ge (d)

nDVe

- число Рейнольдса; -

динамическая вязкость воздуха.

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

ЛриНж(d) = 2,78• 107ß5 X

ч20 2

/■ — \ 0,379

X

Vd J

D Y p w2 (11)

DJ (1+ß>pw2,

где й* = (ув/-меи) - характерный диаметр

частиц при Яе = ¡К = 1! Wви - скорость витания частиц; ув - кинематическая

вязкость воздуха.

Средняя скорость витания частиц определяется следующим образом [5]:

-ч

Weu = J WeU1(x)F(b,n,x)dx 0

xmax

+ J Weu 2 (x)F(b n x)dx,

+

x1

где w (х) = — gx2—--скорость вита-

18 Ре^е

ния при Reeü < 0,2; Reeu =

w„, x

- число

Рейнольдса по скорости витания;

Weu2( Х) =

0,1

Re, x

К

- скорость витания при

Ree

4000;

Яеви1 = (0,568 + 0,133Яе2вМ+0,00045 ЯеАвМ)-число Рейнольдса по скорости витания WвИ 2; Яеви0 = WвИo х/ ^в - число Рейнольдса

V

e

X

по скорости витания;

КиО = U55^/gXp7pУР ; Х - диаметр

частиц; х, = (О,2• 18р^2 /(Рчg)/■33 - граничное значение диаметра частиц при Реви = 0,2, м; хтах - максимальный диаметр частиц (в рассматриваемом случае 0,0005 м).

Подставляя (10) и (11) в (9), а затем в (8) и учитывая (7), получаем функционал для минимизации энергозатрат на размол и пневмотранспорт газовзвеси в ШБМ:

Ф^) =

а

МПм

2 ЕРЛ.

1,214-101'

Re15

Ф( d )■

M I D

|1 + M V L

2,78 • 107 M'5( 1 + m)

( я -A0'379

d_ d*

w„

На рис. 4 представлены экспериментальная зависимость полного удельного расхода электроэнергии на пылеприго-товление Э и теоретического функционала Ф от среднего диаметра частиц угля.

3

20

140

100

60

20

: Э, кВт-ч/т 3

* •

_ » • • • \ Ф

• с/, мкм

40

60

80

100

120

Рис. 4. Зависимость энергозатрат от среднего диаметра частиц угля Fig. 4. Energy consumption dependence on average coal particle diameter

Заключение

Из рис. 4 видно, что по экспериментальным точкам минимальный расход электроэнергии равен 84,3 кВтч/т, при этом значение среднего диаметра равно 75,1 мкм, а по теоретическому расчету расход электроэнергии - 71,4 кВтч/т, средний диаметр - 72,3 мкм. Разница в энергозатратах между теоретическим и экспериментальным расчетами составляет порядка 15%.

В эксперименте [2] мощность ЫЭ,

потребляемая на вращение барабана, изменялась в диапазоне от 720 до 950 кВт в зависимости от степени заполнения ШБМ шарами = 0,191; 0,231; 0,263; 0,3;

0,318, при этом щбопт ~ 0,2 для данного типа брони. Теоретическая мощность, рассчитанная по уравнению (5), равна 660 кВт. Имеются два фактора, которые влияют на увеличение энергозатрат при размоле топлива:

- степень заполнения ШБМ шарами;

- мощность на вращение барабана, измеренная в эксперименте, больше мощности, рассчитанной по уравнению (3) во всех измерениях, {щы)> Щ

(см. рис. 2).

Сформулируем основные результаты:

- получена функция для минимизации энергозатрат на размол угля и пневмотранспорта в ШБМ ШК-32;

- предварительные оценки показывают, что теоретически рассчитанные энергозатраты ниже экспериментальных на 15%.

Библиографический список

1. Левит Г.Т. Пылеприготовление на тепловых электростанциях. М.: Энергоатомиздат, 1990. 384 с.

2. Елизаров В.В., Сеннов В.С. Отчет по тепловым испытаниям котла ТП-81 и пылесистемы 5А Иркутской ТЭЦ-9. Иркутск: РЭУ «Иркутскэнерго», 1970. 234 с.

3. Соколов Н.В., Кисельгоф М.Л. Расчет и проектирование пылеприготовительных установок котельных агрегатов. Л.: ОНТИ ЦКТИ, 1971. 312 с.

4. Guangbin Dr.Yu., Кузнецова М.М., Мараховская М.Б., Алексина А.А. Определение энергозатрат на

измельчение твердых материалов // Журнал технической физики. 2015. Т. 85. № 5. С. 145-147. 5. Иванов С.Д., Кудряшов А.Н., Ощепков В.В. Математическая модель аэродинамического сопротивления шаровой барабанной мельницы при транспорте угольной газовзвеси // Горение топлива: теория, эксперимент, приложения: материалы IX Всероссийской конференции с международным участием (Новосибирск, 16-18 ноября 2015 г.). Новосибирск, 2015. С. 106-110.

References

1. Levit G.T. Pyleprigotovlenie na teplovyh ehlek-trostanciyah. [Pulverization in thermal power plants]. Moscow, Energoatomizdat Publ., 1990, 384 p. (In Russian)

2. Elizarov V.V., Sennov V.P. Otchet po teplovym ispytaniyam kotla TP-81 i pylesistemy 5A Irkutskoj TEC 9 [Report on the thermal tests of TP-81 boiler and 5A coal pulverization system of Irkutsk CHP 9]. Irkutsk, REU "Irkutskenergo" Publ., 1975, 234 p. (In Russian)

3. Sokolov N. V., Kisel'gof M. L. Raschet i proektiro-vanie pyleprigotovitel'nyh ustanovok kotel'nyh agrega-tov [Calculation and design of boiler unit pulverization plants]. Leningrad, ONTI CKTI Publ., 1971, 312 p. (In Russian)

4. Guangbin Dr. Yu., Kuznecova M.M., Marahovskaya M.B., Aleksina A.A. Opredelenie energozatrat na iz-

Критерии авторства

Иванов С.Д., Кудряшов А.Н., Ощепков В.В. исследовали влияние диаметра угольных частиц на расход электроэнергии на пылеприготовление шаровой барабанной мельницы ШК-32 с целью оптимизации ее работы, провели обобщение и написали рукопись. Ощепков В.В. несет ответственность за плагиат.

Конфликт интересов

Авторы заявляют об отсутствии конфликта интересов.

Статья поступила 16.08.2016 г.

mel'chenie tverdyh materialov [Determination of energy consumption for solid material grinding]. Zhurnal tekhnicheskoj fiziki [Technical Physics]. 2015, vol. 85, no. 5, pp. 145-147. (In Russian) 5. Ivanov S.D., Kudryashov A.N., Oshchepkov V.V. Matematicheskaya model' aehrodinamicheskogo sopro-tivleniya sharovoj barabannoj mel'nicy pri transporte ugol'noj gazovzvesi [Mathematical model of ball cylindrical mill aerodynamic resistance under the coal gas suspension transport]. Materialy IX Vserossijskoj kon-ferencii s mezhdunarodnym uchastiem "Gorenie topliva: teoriya, eksperiment, prilozheniya" [Proceedings of the IX All-Russian Conference with International Participation "Fuel combustion: theory, experiment, application"]. Novosibirsk, 2015, pp. 106-110. (In Russian)

Authorship criteria

Ivanov S.D., Kudryashov A.N., Oshchepkov V.V. studied the influence of the coal particle diameter on power consumption for the dust preparation by the ball pulverizer mill BC-32 in order to optimize its operation. They summarized the material and wrote the manuscript. Oshchepkov V.V. bears the responsibility for plagiarism.

Conflict of interests

The authors declare that there is no conflict of interests regarding the publication of this article.

The article was received 16 August 2016