Опытный размол бурого угля в стендовой среднеходной мельнице Текст научной статьи по специальности «Технологии материалов»

УДК 621.926

ОПЫТНЫЙ РАЗМОЛ БУРОГО УГЛЯ В СТЕНДОВОЙ СРЕДНЕХОДНОЙ МЕЛЬНИЦЕ

ШУВАЛОВ С.И., д-р техн. наук, МИХЕЕВ Г.Г., канд. техн. наук

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

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

EXPERIMENTAL BROWN COAL GRINDING IN BOWL MILL

S.I. SHUVALOV, Ph.D., G.G. MIKHEEV, Ph.D.

The work represents the results of experiment-calculated research of damp brown coal grinding in bowl mill, and the results of kinetic research of coal grinding on mill grindability table as well. The authors suggest the mathematical model for dispersed makeup and substance mass flow compiling in ventilated roller with separator. This model describes the processes of big particles grinding, agglomeration of smaller particles, particle classification in upper grindability place and in the main separator. There was carried out the evaluation of damp coal influence, grindability bowl types and separator division effectiveness over mill productivity.

Key words: coal-preparation system, grindability apparatus, grinding intensity, mathematical model, carbonic particle agglomeration process.

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

Для оценки влияния влажности угля на условия размола в валковой среднеходной мельнице и выбора наиболее пригодной конструкции отдельных узлов пы-лесистемы на стендовой установке ИГЭУ [1] был проведен опытный размол бурого угля березовского месторождения при изменении его влажности от воздушно-сухого состояния К = 12 %) до предельной влажности wp = 31 %.

В первой серии опытов проводилось исследование кинетики измельчения угля на размольном столе. До начала опытов с мельницы был снят сепаратор пыли и демонтирован сопловой аппарат. Исходный материал подавался на размольный стол по центральной течке, проходил под валками и собирался в бункере провала. В опытах использовался уголь (табл. 1) после дробилки и отсев на сите 6,3 мм.

Таблица 1. Дисперсные составы опытных партий углей

5, мкм Полные остатки R,%

уголь после дробилки отсев на сите 6,3 мм

7000 11,6 0

5000 13,9 5

2500 37,8 38

1600 51,7 56

710 70 72,1

400 77,5 80

200 85,5 85

90 91,3 93

0 100 100

ОАО «Тяжмаш» оснащает среднеходные мельницы коническими валками, ведущие в этой отрасли фирмы «Бабкок» и «Леше» устанавливают тороидальные валки. Для оценки влияния формы валков опыты на стендовой мельнице проводились с обоими типами валков.

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

Анализ зависимости этого параметра от влажности угля, поступающего в мельницу (рис. 1), показывает, что при влажности угля wp = 12-25 % более высокую интенсивность измельчения обеспечивают тороидальные валки, при более высокой влажности их преимущество быстро убывает, а при влажности выше 30 % предпочтительны конические валки.

Рис. 1. Влияние влажности угля на кратность измельчения: 1 -конические валки, дробленка; 2 - тороидальные валки, дроблен' ка; 3 - конические валки, отсев; 4 - тороидальные валки, отсев

Присутствие в исходном материале даже небольшого количества относительно крупных частиц заметно ухудшает интенсивность размола. Отсев из дробленого угля 10 % частиц круп-

і

нее 6,3 мм (рис. 1, кривые 3 и 4) существенно повышает кратность измельчения во всем диапазоне изменения влажности угля.

Оценка интенсивности измельчения по кратности ки является достаточно грубой. Она очень чувствительна к погрешностям дисперсного анализа и не учитывает изменение структуры дисперсного состава. Вместе с тем было замечено, что после прохождения через мельницу угля с влажностью 31 % содержание в продуктах размола частиц мельче 90 мкм стало ниже, чем в исходном материале. Этот результат свидетельствует о том, что наряду с разрушением частиц исходного материала на рабочем столе происходит объединение мелких частиц в более крупные агломераты. Кроме того, при влажности угля выше 25 % на размольном столе с тороидальными валками начинают появляться пластинчатые наслоения, укрупняющиеся со временем до 20-25 мм. Эти наслоения обладают малой механической прочностью и при ситовом анализе разрушаются, однако при сушке угля в промышленной мельнице такие образования могут упрочниться и вызвать серьезные затруднения в работе мельницы. При размоле влажного угля коническими валками таких наслоений не наблюдалось.

Для более детального исследования процесса формирования дисперсного состава на рабочем столе мельницы и влияния особенностей размола влажных углей на производительность пылесистемы была разработана математическая модель среднеходной мельницы и проведено численное моделирование процесса.

В теплоэнергетике принято дисперсный состав пыли описывать с помощью кривой полных остатков Щб), определяющей процентное содержание в поли-дисперсной смеси частиц крупнее размера б. Для полных остатков выполняются граничные условия

К (тах ) = 0, К (б™ ) = 100%.

Здесь бтах, бт|п - максимальный и минимальный размер частиц в смеси. Обычно принимают бтт=0.

Для описания процессов измельчения и классификации более удобно пользоваться плотностью распределения частиц по размерам 1(б), которая связана с полными остатками соотношениями

К (б)-100 Т' . < (б) = -■

8

При этом из условия нормировки следует

тах | і

= 1.

Наиболее часто кинетику измельчения описывают с помощью селективной и распределительной функций [2]:

(1)

Селективная функция Э(б) определяет интенсивность разрушения частиц размером б в единицу времени или в течение одного цикла нагружения. Распределительная функция В(б,£) соответствует доле продуктов разрушения частиц начальным размером б, оказавшихся после разрушения мельче размера £. Для этих функций должны выполняться следующие соотношения:

Б (б)> 0, Б (0) = 0;

В (б, ^)< В (б,^) при ^ > ^2; (2)

В (б, б) = 0; В (б,0) = 1.

Формула (1) предполагает, что при нагружении более крупные частицы разрушаются с образованием более мелких осколков, при этом полные остатки на всех размерах уменьшаются. Для учета возможности укрупнения мелких частиц аналогично (1) введем селективную функцию агломерации А(б), определяющую интенсивность образования новых агломератов с участием частиц размером б, и распределительную функцию агломерации Н(б,5) [3], значения которой равны весовой доле частиц классов интервала размеров [б, б+с1б], перешедших во вновь образованные агломераты размером крупнее £. Очевидно, что должны выполняться следующие соотношения:

А (б)> 0,

Н(6,^)> Н(б,^2) ПРИ ^1 < ^2, (3)

Н(б, б) = 1, Н(б, бтах ) = 0.

Пусть в начальный момент времени 1 = 0 дисперсный состав размалываемого материала соответствует плотности распределения 10(б). За период времени С из класса частиц размера [б, б+Сб] вследствие разрушения и агломерации уйдет материал в количестве

^ (б) [Б (б) + А (б)] Ьб^ . В этот же класс за счет

осколков от разрушения более крупных частиц попадет материал в количестве

тах

і Б (5)

8В (5,8). ----—----.

58

ґ0 (5)<15<. Также в этот класс

перейдут агломераты из более мелких частиц в

количестве -

| а(5)^Н(8))о (5)^5л.

Знак «-»

вызван тем, что производная

дН(5,8)

58

отрица-

тельна. Окончательно уравнение измельчения и агломерации примет вид

= -[Б (8) + а (8)] .0 (8)+ )ах Б (5)дВ8 (5)0

і а (5)

дН (58))

----------іп

58

(5) <15.

(4)

Распределительная функция измельчения В(8,5) определяет долю осколков размером менее 5, образовавшихся после разрушения частицы размером 8. Результаты наших исследований [4-6] показали, что при достаточном энергоподводе происходит хрупкое разрушение частиц с равномерным распределением продуктов разрушения по размерам. Это соответствует линейному распределению, которое может быть представлено зависимостью

в о8,5)=§.

(5)

Такое распределение наблюдается при измельчении различных материалов в шаровых

8

8

8

барабанных, вибрационных и струйных мельницах, когда одиночные частицы разрушаются стесненным или свободным ударом. Вместе с тем механизм раздавливания частиц в слое на столе валковой мельницы вносит определенную специфику в формирование дисперсного состава измельчаемого материала. Если в мельницах другого типа подведенная при ударе к частице энергия оказывается большей, чем требуется для хрупкого разрушения, то ее излишек расходуется в основном на увеличение кинетической энергии образующихся осколков. В валковой мельнице этот путь диссипации избытка энергии крайне ограничен. Образующиеся после первичного разрушения частицы осколки не могут быстро уйти из-под валков и при избытке энергии повторно подвергаются нагружению.

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

Селективная функция измельчения S(5) определяет интенсивность разрушения частиц размером 5 в единицу времени или за цикл нагружения. При численном интегрировании значение S(5i) означает долю материала, перешедшего из интервала размеров [5i, 5i+i] в более мелкие интервалы.

Наиболее часто селективную функцию аппроксимируют степенной зависимостью

S (5) = а5р, (6)

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

В шаровых барабанных и вибрационных мельницах интенсивность измельчения частиц определенного размера зависит от конструкции и режима работы мельницы и от прочностных свойств материала. В валковых мельницах в качестве определяющего фактора дополнительно должен быть введен некоторый параметр, характеризующий долевое содержание частиц других классов крупности. Согласно [7], в (6) введем коэффициент

kf (8) =

(7)

где бт - медианный размер полидисперсной смеси на размольном столе; у - параметр, определяемый из эксперимента.

Тогда выражение для селективной функции будет иметь вид

(8)

Процесс агломерации угольных частиц изучен в значительно меньшей степени, чем их разрушение. По аналогии с (7) будем считать, что интенсивность агломерации частиц размером б определяется выражением

A (б) = цбю, (9)

где ц и ю - коэффициенты, зависящие от свойств материала и условий измельчения.

Поскольку процесс образования агломератов связан с взаимодействием поверхностей частиц, то предположим, что частицы размером б переходят в более крупные классы пропорционально суммарной поверхности частиц этих классов. В этом случае распределительная функция Н(б, А), определяющая долю частиц начальным

размером б, переходящих в агломераты крупнее А, должна иметь вид

H(8, ^)=-

(10)

Уравнение (4) требует численного решения. Разделим весь диапазон изменения размеров на n классов таким образом, чтобы 81 =8m¡n, Sn+i=Smax. Представим непрерывный процесс размола в виде последовательных циклов с периодом At. Уравнение (4) примет вид

f (8) = f0 (8,)[1 -(S(8,) + A(8i))At] +

i-1 n (11)

+At^ S (8k ) bkif0 (8k ) + At ^ A (8k ) hkif0 (8k ) ■ k=1 k=i+1

Здесь коэффициенты bk¡ определяют долю образующихся осколков размером [8¡, 8¡+1] при измельчении частиц размером [8k, 8k+1]; hk¡ определяют долю частиц, перешедших при агломерации из диапазона [8k, 8k+1] в диапазон [8¡, 8¡+1].

При описании дисперсных составов в виде набора n узких фракций будем считать, что характерным размером ¡-й фракции является среднее арифметическое границ интервала фракции:

• 8,+1

2

Выражение (8) принимает вид

S (8) = а8в

(12)

(13)

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

ьк1^б|-б:+1 . (14)

8k+1 8n+1

Коэффициенты hk¡ равны f0 (8k )

hki=Ti

(15)

Уравнение (12) принимает вид

8

*1 (8,) = *0 (8)

і -1 і

(

і -

а8в( 8

8т

V

ц8ю

А!

+А!

1 -1 8«іЧУ 8|-8і+1 * (88 )

8=1

8т) 88+1 8п+1

(16)

+А!ц

8?^ (88^ °, 8=|+1 5 *

0 (8к ) 8к

При отсутствии агломерации (16) переходит в уравнение

*1 (8,) = *о (8,)

( (8 У'

1 - а 8, V V 8т ], 1

+А!

1 -1 ТТЯ ( 88 У 8 - 8;

£а8*1 і8

8=1

Ї+1

8т ) 88+1 - 8п+1

*о (88)

(17)

Параметрами идентификации модели измельчения и агломерации являются коэффициенты а, р, у, ц и ю, а также количество эквивалентных нагружений при прохождении материала под валком Пг. При отсутствии агломерации, в соответствии с (17), в модели чистого измельчения остаются параметры а, р, у и Пг.

При анализе результатов опытного размола было отмечено, что наибольшую зависимость от влажности имеют параметры а и ц, отражающие интенсивность процессов измельчения и агломерации. В соответствии с этим на первой стадии обработки экспериментальных данных для всего массива независимо от влажности угля и типов размольных валков было найдено, что наблюдаемая в опытах структура дисперсного состава размолотых частиц достигается при следующих значениях коэффициентов: Пе = 7; р = 0,65; у = -0,25; ю = -0,85. Оставшиеся параметры а и ц определялись путем минимизации суммы квадратов расхождений значений полных остатков, полученных ситовым анализом и рассчитанных по (16) (рис. 2, 3).

Рис. 2. Зависимость коэффициента пропорциональности селективной функции измельчения от влажности угля: 1 -конические валки; 2 - тороидальные валки

При использовании конических валков с увеличением влажности параметр а сначала несколько возрастает с 0,00154 при \^р = 12 % до 0,0019 при вдР = 18,8 %, а затем снижается практически до первоначального значения а = 0,00149 при wp = 25 % и а = 0,00147 при wp = 31 %. Тороидальные валки во всем

диапазоне изменения влажности обеспечивают несколько большую скорость измельчения: а = 0,00176 при wp = 12 %, а = 0,00244 при ад/5 = 25 % и а = 0,00152 при чР = 31 %.

Рис. 3. Зависимость коэффициента пропорциональности селективной функции агломерации от влажности угля: 1 - конические валки; 2 - тороидальные валки

Однако конические валки менее способствуют процессу агломерации. При размоле воздушно-сухой пробы с wp = 12 % коэффициент Ц, пропорциональный скорости агломерации, однозначно равен нулю. При увеличении влажности угля до 25 % значение параметра несколько увеличивается до 0,56-0,64 и только при более высокой влажности wp = 31 % коэффициент ц быстро возрастает до 7,35.

При размоле угля тороидальными валками агломерация продуктов разрушения более интенсивна. Значение коэффициента пропорциональности ц более монотонно возрастает с 0,76 при wp = 12 % до 3,34 при wp = 25 % и до 13,15 при wp = 31 %.

Вторая серия опытов включала в себя измельчение доведенного до воздушно-сухого состояния с wp = 12 % бурого угля в вентилируемой мельнице с сепаратором пыли конструкции ИГЭУ (табл. 2).

Таблица 2. Результаты опытов по измельчению бурого угля в вентилируемой мельнице с сепаратором, wp =12 %

Параметры Опытные данные

Производительность, кг/ч 144 144 136 90 120 116

Вентиляция, м 3/ч 780 580 450 275 360 360

Высота зонь , мм 100 100 100 200 200 200

Лопатки сепаратора, град 45 45 45 30 45 30

Дисперсные составы, К,%

710 2 1 0,2 0 0 0

400 12 7,7 1,5 0 0 0

Размер ячеек сит, 200 40 32 10 0,5 2,3 3,6

мкм 90 69 56 28 19 28 34

0 100 100 100 100 100 100

Опыты показали, что получение более мелкой готовой пыли может быть достигнуто за счет снижения вентиляции мельницы или уменьшения границы разделения сепаратора (за

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

Для оценки производительности мельницы в зависимости от влажности угля, эффективности разделения сепаратора и тонкости готовой пыли по математической модели вентилируемой среднеходной мельницы с сепаратором пыли [1], в которой процесс разрушения и агломерации описывался уравнением (17), были проведены вариантные расчеты. При этом принималось, что максимальная производительность мельницы ограничивается одинаковым для всех случаев максимально допустимым потоком сырого угля и возврата на размольный стол мельницы. Для большей наглядности результаты оценки производительности пылесистемы представлены в относительном виде, где в качестве нормирующей величины принята производительность мельницы с коническими валками и сепаратором типа ТКЗ-ВТИ при размоле угля с влажностью wр = 12 % до тонкости Р90 = 30 %.

Анализ зависимостей производительности мельницы от тонкости готовой пыли при размоле воздушно-сухого угля с влажностью wр = 12 % (рис. 4) показывает, что при угрублении готовой пыли с 1^90 = 30 % до Р90 = 35 % максимальная производительность мельницы с коническими валками и сепаратором ТКЗ-ВТИ увеличивается на 8 %; дальнейшее угрубление пыли до Р90 = 40 % позволяет поднять производительность мельницы на 18 %; уменьшение тонкости пыли до Р90 = 15 % снижает производительность мельницы на 26 %.

В 3 4

\хг

\2_ ^90,%

10

15

20

25

30

35

40

45

Рис. 4. Влияние тонкости готовой пыли на производительность мельницы, \^р = 12 %: 1 - конические валки и сепаратор ТКЗ-ВТИ;

2 - конические валки и сепаратор ИГЭУ; 3 - тороидальные валки и сепаратор ТКЗ-ВТИ; 4 - тороидальные валки и сепаратор ИГЭУ

Установка на эту мельницу более эффективного статического сепаратора ИГЭУ позволяет увеличить производительность мельницы при одинаковой тонкости готовой пыли на 2 % при Р90 = 30 % и на 7 % при Р90 = 15 %.

Применение тороидальных валков с сепаратором ТКЗ-ВТИ при Р90 = 30 % увеличивает производительность мельницы на 9 %, а при Р90 = 40 % - на 11 %. Вместе с тем по мере утонения пыли преимущество тороидальных валков снижается и при Р90 < 15 % более эффективными становятся конические валки.

Применение сепаратора ИГЭУ усиливает положительный эффект от использования тороидальных валков. При Р90 = 30 % производительность мельницы становится выше на 13,8 % по сравнению с базовым вариантом. При Р90 = 35 % и Р90 = 40 % прирост производительности составляет 16 % и 19,4 %, соответственно. При Р90 = 15 % установка сепаратора ИГЭУ и тороидальных валков увеличивает производительность мельницы на 10,8 %.

Анализ зависимостей производительности мельницы с коническими валками и сепаратором ТКЗ-ВТИ от влажности угля (рис. 5) показывает, что с

увеличением влажности от 12 % до 18,6 % производительность мельницы возрастает на 13,517,7 %. Дальнейшее повышение влажности до 25 % приводит к снижению производительности по сравнению с размолом сухого угля на 7-13 %. При влажности 31 % производительность мельницы падает на треть при Р90 = 40 % и более чем в три раза при Р90 = 15 %.

Рис. 5. Зависимости производительности мельницы с коническими валками и сепаратором ТКЗ-ВТИ от тонкости готовой пыли: 1 - ^ = 12 %; 2 - ^ = 18,6 %; 3 - wp = 25 %; 4 - wp = 31 %

Анализ зависимостей производительности мельницы от тонкости пыли при размоле угля с влажностью 31 % (рис. 6) показывает, что при тонкости готовой пыли 1^90 = 30 % производительность мельницы с тороидальными валками и сепаратором ТКЗ-ВТИ снижается более чем в два раза, по сравнению с размолом сухого угля.

В этих условиях более перспективным является использование мельниц с коническими валками и статическим сепаратором конструкции ИГЭУ.

При 1^90 = 30 % производительность такой мельницы на 20 % выше, чем у мельницы с коническими валками и сепаратором ТКЗ-ВТИ и на 33 % выше, чем у мельницы с тороидальными валками и сепаратором ТКЗ-ВТИ.

Рис. 6. Зависимости производительности мельницы от тонкости готовой пыли, шр = 31%: 1 - конические валки и сепаратор ТКЗ-ВТИ; 2 - конические валки и сепаратор ИГЭУ; 3 - тороидальные валки и сепаратор ТКЗ-ВТИ

Следует обратить внимание еще на одну особенность размола влажных углей. При влажности угля 31 % для получения пыли с Р90 = 30 % статический сепаратор ТКЗ-ВТИ должен быть настроен на остаточно низкую границу 10-11 мкм. Более эффективный сепаратор ИГЭУ должен разделять мельничный продукт по более высокой границе 22-23 мкм.

1,5

0,5

Результаты расчетно-экспериментального исследования измельчения бурого угля в валковой сред-неходной мельнице показали, что в достаточно широком диапазоне изменения влажности для получения относительно крупной пыли с Р90 > 15 % более предпочтительны мельницы с тороидальными валками. Однако существует некоторое пороговое значение влажности угля, выше которого могут возникнуть существенные затруднения при его размоле. Мельницы с тороидальными валками при более высокой интенсивности измельчения не исключают образования крупных агломератов, которые могут нарушить стабильность работы пылесистемы. При влажности, когда процессы агломерации становятся более существенными, предпочтительнее использовать конические валки.

Вопрос о выборе типа валков для получения более мелкой пыли с 1^90 < 10 % является открытым. Следует ожидать, что в этом случае также могут оказаться более эффективными конические валки.

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

Представленные результаты могут быть использованы при проектировании промышленных среднеход-ных мельниц для размола материалов, склонных к агломерации, и получения тонкодисперсных порошков.

Список литературы

1. Шувалов С.И., Михеев Г.Г., Кишкин Ю.Н.

Выбор конструкции и типоразмера сепаратора для валковой среднеходной мельницы // Вестник ИГЭУ. -2003. - Вып. 3. - С. 22-26.

2. Гарднер Р.П., Аустин Л.Г. Исследование измельчения в мельнице периодического действия / В кн.: Труды Европейского совещания по измельчению. -М.: Стройиздат, 1966. - С. 219-248.

3. Шувалов С.И., Смирнов А.С., Краснов Е.В. Измельчение и агломерации частиц в шаровых барабанных мельницах // Вестник ИГЭУ. - 2003. - Вып. 1. -

С. 65-67.

4. Шувалов С.И. Закономерности преобразования дисперсного состава материала в процессах его измельчения в шаровых барабанных мельницах: Дис... канд. техн. наук. - Иваново, 1983.

5. Михеев Г.Г. Интенсификация процессов измельчения сыпучих материалов в среднеходных валковых мельницах: Дис. канд. техн. наук. - Иваново, 1984.

6. Шувалов С.И. Структурная и режимная оптимизация процессов фракционирования порошков: Дис. д-ра техн. наук. - Иваново, 1995.

7. Лебедев Д.Е. Распределение энергии по фракциям материала при его измельчении и его влияние на прогнозирование фракционного состава: Дис. канд. техн. наук. - Иваново, 2001.

Шувалов Сергей Ильич,

ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор кафедры тепловых электрических станций, телефон (4932) 26-99-13, e-mail: [email protected]

Михеев Геннадий Григорьевич,

ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры тепловых электрических станций, телефон (4932) 26-99-13, e-mail: [email protected]