CC BY
21
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
УДК 004.942
Оптимизация тепло-технологического процесса сушки движущейся плотной многослойной массы фосфоритовых окатышей по критерию энергоресурсоэффективности
Бобков В. И.
Постановка задачи: проблема рационального использования топливно-энергетических ресурсов наиболее важна для энергоёмких промышленных производств, к числу которых относится термическая подготовка окомкованного рудного сырья в сложной многостадийной тепло-технологической системе (ТТС) обжиговой конвейерной машине для производства фосфоритовых окатышей. Оптимизация режимов сушки движущейся плотной многослойной массы фосфоритовых окатышей в сложной ТТС, позволит повысить энергоресурсоэффективность ТТС в результате интенсификации и оптимизации тепло-массообменных процессов (ТМП) по минимуму себестоимости затрачиваемой электрической и тепловой энергии. Целью работы является повышение энергоресу-соэффективности ТТС, с использованием математической и компьютерной модели, процедуры оптимизации тепло-технологического процесса (ТТП) сушки движущейся плотной многослойной массы фосфоритовых окатышей, отличающейся использованием дискретного динамического программирования и процедур контроля эффекта перераспределения влаги по высоте многослойной засыпки окатышей, что позволяет предотвращать появление горизонтов переувлажнения, негативно влияющих на газопроницаемость многослойной массы окатышей, приводящих к необоснованному росту энергопотребления и снижению производительности ТТС в целом. Используемые методы: методы системного анализа, дискретного динамического программирования, теории оптимального управления системами с распределенными параметрами, условной многокритериальной оптимизации, вычислительной математики. Результат: представлены содержательная и математическая постановки задачи оптимизации (ТТП) сушки движущейся плотной многослойной массы фосфоритовых окатышей в сложной, многостадийной (ТТС) конвейерной обжиговой машине, отличающаяся учетом тепло-технологических особенностей обжиговой машины, интенсивности процесса внутреннего влагопереноса в окатыше и процессов переувлажнения отдельных слоев окатышей и газа-теплоносителя, что позволяет повышать энергоэффективность посредством интенсификации ТМП динамической многослойной сушки. Новизна: разработаны математическая и компьютерная модели оптимизации ТТП сушки движущейся плотной многослойной массы фосфоритовых окатышей, позволившие оптимизировать энергоресурсоэффективность сложной, динамической ТТС производства фосфоритовых окатышей. Полученные результаты применялись для расчета энергоэффективной сушки окатышей в ТТС обжиговой конвейерной машине. Обнаружено, что в оптимальном режиме многослойной сушки окатышей отсутствует зона переувлажнения, интенсифицируются ТМП, снижается расход энергии и увеличивается качество готового продукта, уменьшается доля возврата, обеспечивая ресурсосбережение.
Ключевые слова: тепло-технологическая система, тепло-технологический процесс, сушка, оптимизация, окатыши, фосфорит, энергоресурсоэффективность, обжиговая машина.
Библиографическая ссылка на статью:
Бобков В. И. Оптимизация тепло-технологического процесса сушки движущейся плотной многослойной массы фосфоритовых окатышей по критерию энергоресурсоэффективности // Системы управления, связи и безопасности. 2018. № 2. С. 56-68. URL: http://sccs.intelgr.com/archive/2018-02/04-Bobkov.pdf
Reference for citation:
Bobkov V. I. The heat-technological process optimization of a moving dense multilayer weight of phosphorite pellets drying process by the energy-resource efficiency criterion. Systems of Control, Communication and Security, 2018, no. 2, pp. 56-68. Available at: http://sccs.intelgr.com/archive/2018-02/04-Bobkov.pdf (in Russian).
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
Актуальность
Проблема рационального использования сырьевых и топливно-энергетических ресурсов особенно важна для энергоемких промышленных производств, к которым относятся конвейерные обжиговые машины для производства фосфоритовых окатышей, представляющие собой сложные тепло-технологические системы (ТТС). Режимы функционирования ТТС, а также условия проведения и завершенность тепло-технологического процесса (ТТП) сушки движущейся на конвейере плотной многослойной массы фосфоритовых окатышей, определяют качество термической обработки окатышей. Изменение физико-химического или гранулометрического состава исходного минерального фосфатного сырья, поступающего в ТТС, требует проведения дорогостоящих экспериментальных исследований по определению тепло-физических и физико-химических свойств сырья, в соответствии с которыми осуществляется перенастройка параметров технологических режимов сложной ТТС обжиговой конвейерной машины [1, 2].
Выбор рациональных технологических схем и адаптация оптимальных технологических режимов ТТС к периодически меняющимся свойствам фосфатного сырья представляют собой сложную научно-техническую проблему, для решения которой необходимы дорогостоящие экспериментальные исследования, а порой и не осуществимые, из-за особенностей протекания ТТП во времени при больших температурах. Кроме того, все существующие в настоящее время методики расчёта энергоресурсоэффективности ТТП термической подготовки фосфатного сырья являются эмпирическими [3].
В связи с этим, проблема разработки методического и программно-информационного обеспечения системного анализа и процедур принятия решений по оптимизации энергоресурсоэффективности ТТС производства фосфоритовых окатышей является актуальной научной проблемой, успешное решение которой обеспечит повышение энергоресурсоэффективности ТТС и высокое качество готовых окатышей [4, 5].
Изменение свойств сырых фосфоритовых окатышей: влажности, прочности, гранулометрического состава, является причиной нарушения технологического режима ТТП сушки. Поэтому реализация ТТП сушки в режиме оптимума энергозатрат, на основе компьютерного моделирования имеет большое значение для практики [6]. Поскольку технологические ограничения ТТС существенно сужают диапазон варьирования параметров сушки, требуется тщательный анализ взаимовлияния всех факторов, определяющих режим сушки [7].
Содержательная постановка задачи оптимизации
Для постановки и решения задачи в работе введены обозначения, представленные в таблице 1.
Технологический режим сушки фосфоритовых окатышей определяется: температурой газа-теплоносителя на входе в движущийся плотный слой Т^о(т) и скоростью его перекрёстной подачи Жё(т). Качество и завершённость ТТП сушки характеризуются: тепловыми параметрами - температура окатышей Тт(х, у, т), зависящая от координаты х по радиусу окатыша, координаты у высо-
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
ты в вертикальной многослойной укладки окатышей и времени т её пребывания в зоне сушки при горизонтальном движении многослойной массы на конвейере обжиговой машины ХЭТС, градиент температуры дТт/дх, скорость нагрева дТт/дх, температура газа на выходе из слоя окатышей Tgh(т) и параметрами состояния - влагосодержание и(у, т), интенсивность влагопереноса 1(х, у, т), относительная степень высушивания у(у, т) [8].
Таблица 1 - Обозначения
Обозначение Физический смысл обозначения
Tg - температура газа-теплоносителя, К
Tgo - температура газа-теплоносителя на входе в многослойную массу окатышей, К
Tgh - температура газа-теплоносителя на выходе из многослойной массы окатышей, К
Wg - скорость подачи газа-теплоносителя, м/с
Tm - температура материала окатышей, К
T - время, с
X - координата радиуса окатышей
y - координата вертикальной многослойной укладки окатышей (высота слоя)
I - интенсивность влагопереноса в окатыше, кг/(м2-с)
u - влагосодержание, кг/кг
¥ - относительная степень высушивания окатышей
Q - удельный расход энергии, Дж/кг
p - себестоимость энергии, руб/Дж
i - длина конвейера ТТС обжиговой машины, м
эл - индекс для обозначения электрической энергии
теп - индекс для обозначения тепловой энергии
кон - индекс для обозначения конечного значения времени т пребывания окатышей на конвейере ТТС
MAX - индекс для обозначения максимального значения величины
i=\...k - индекс для счетчика номера вакуум-камеры
Критерием завершённости ТТП служит - минимум влагосодержания среднего по высоте 0 < у < 0,35м вертикальной многослойной укладки окаты-
0,35
шей и - | ЧУ', ткон на выходе из зоны сушки т = Ткон.
0
Себестоимость готового продукта в производстве фосфоритовых окатышей формируется главным образом за счет удельного расхода электрической 0^л и тепловой 0геп энергии [9]. Оптимальным считается технологический режим, при котором необходимое качество готового продукта, определяемое завершенностью процесса сушки и - 0, достигается при минимуме стоимости израсходованной тепловой и электрической энергии, причем учитывается, что стоимость электрической энергии Рэл в несколько раз дороже тепловой Ртеп. Задача оптимизации формулируется так: найти такое значение вектора управляющих параметров Х(Т^(т), ^(т)), чтобы значения и и Р = Рэл£эл[^(т)] + РтепО-еп^о(т)] дости-
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
гали своего наименьшего значения [10]. В силу технических особенностей обжиговой машины на вектор X управляющих параметров накладываются ограничения: на скорость Wg(т) < WgMAx и температуру 7^)(т)<Тмх газа-теплоносителя на входе в слой, а также на вектор параметров состояния ХЭТП сушки £(Тт, дТт/дх, дТт/дт, Тф, и, I, у) - температуру газа-теплоносителя на выходе Tgh < ТёкМАХ, скорость нагрева дТт/дт <^ХТМХ и градиент дГт/дх темпера-
туры в окатыше, влагосодержание и < иМАХ и интенсивность влагопереноса I <1МАХ (см. рис. 1).
Рис. 1. Структурно-операторная схема ТТП сушки фосфоритовых окатышей в зоне сушки ТТС обжиговой машины: I - оператор теплообмена, II - оператор сушки.
В качестве критерия эффективности (КЭ) рассматривается функционал:
ЩТЖ), ВД] = ai U [Т^(т), Wg(x)] + Wg(x)], (1)
где коэффициенты а и а2 обеспечивают ранжирование критериев в двухкрите-риальной задаче оптимизации [11, 12].
Математическое описание тепломассообменных процессов сушки движущейся многослойной массы фосфоритовых окатышей представляется в виде системы дифференциальных уравнений в частных производных (ДУЧП), решение которой возможно лишь численными методами, поэтому отыскать функции Тёо(т) и Wg(T) аналитическими методами, используя вариационное исчисление или принцип максимума Понтрягина, не представляется возможным.
Systems of Control, Communication and Security
sccs.intelgr.com
Применение метода дискретного динамического программирования для повышения энергоресурсоэффективности ТТС
Дискретное динамическое программирование является специализированным вычислительным методом решения задач оптимизации сложных динамических систем, позволяющим представить процесс оптимизации в виде последовательности отдельных этапов (шагов). Применение декомпозиционно-топологических принципов позволяет получать все используемые в динамическом программировании функциональные рекуррентные соотношения. Особенностями метода являются хорошая алгоритмизуемость и возможность эффективного использования современных средств вычислительной техники [13].
Рассмотрим декомпозицию ТТП сушки, разбив время пребывания окатышей в зоне сушки ТТС на к равных частей, соответствующих длинам вакуум-камер А/, получив временные промежутки Дт (см. рис. 2).
У, м
Горн формирует начальную температуру газа-теплоносителя Т%
g0
0,35
TgoW
TgOk,Wgk"
Длина зоны сушки /, м Эксгаустеры обеспечивают поток газа-теплоносителя со скоростью
Рис. 2. Дискретизация времени пребывания окатышей в зоне сушки длиной / конвейерной обжиговой машины 0 < т < ткон в зависимости от прохождения отдельных вакуум-камер длиной А/, в которых формируется вектор управляющих параметров Х(Тё0г,
0
Очевидно, что каждый вектор Si последующего состояния параметров ТТП сушки определяется его предыдущим состоянием Si.1 и вектором управляющих параметров Хг{Т^, - Sг■ = X), а КЭ на каждом г-м шаге есть функция, зависящая от параметров состояния ТТП сушки на предыдущем (г-1)-м шаге, и решения X, принятого на данном этапе ^ - ^(Хг; ^ ) (см. рис. 3).
Systems of Control, Communication and Security
sccs.intelgr.com
Такая декомпозиция называется неявной - алгоритм требует координации стоимости Рэл и Ртеп промежуточных переменных Qэл и Qтеп ТТП. Задача оптимизации КЭ (1) на каждом /-м этапе неявной декомпозиции формулируется так: найти такие значения параметров Тё0/, Жё/ оптимизирующего вектора Х/, чтобы КЭ (1) был минимальным:
ВДс, Жё/) = а1 и (То, Жё/)+а2 Р(Тёо/, (2)
02
Xi
Г ^2
Г Хз
й Z
Г Xi
I S, I I S I
■ HgMg>-
Рис. 3. Декомпозиционно-топологическая схема оптимизации ТТП сушки методом дискретного динамического программирования. 80 -начальный вектор параметров состояния исходного сырья, сырых окатышей на входе, В - высушенные окатыши на выходе из зоны сушки, Qi - расход энергии, X/ - вектор управляющих параметров, 8/ - вектор параметров состояния ТТП сушки на /-м этапе декомпозиции.
Решение задачи динамического программирования начинается с последнего этапа. В этом случае условный минимум КЭ будет определяться по формуле:
К ( Sk _ ! ) = цт Fk ( ^ _ X ),
где р* - условный минимум КЭ, X* - условное оптимальное управление на
последнем шаге / = к.
Значение критерия эффективности на предпоследнем шаге (к-1)-м шаге будет равно:
К-1 (Sk_2 ) = К_1 (^^_2, Х^-) + Р (^^- ) (3)
В выражение (3) входит вектор параметров состояния , поэтому, чтобы продолжить решение, необходимо определить, как изменится вектор параметров состояния под влиянием управления на (к-1)-м шаге:
= фк-1(8/-2; Х/-1).
Условный минимум КЭ на (к-1)-м шаге определяется по формуле:
( ^к-2 )= -1 ( -2' Хк-1) + К (Фк-1 ( ^к -2 > Хк-1))}
из которой находится условный минимум К () и условный вектор оптимального управления Х*_1. Далее переходим к шагу (к-2) и т.д.
Для /-го шага ((к-/+1)-й шаг с конца) получается уравнение Беллмана:
S,
B
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
F ( 3-1 )=m? {F ( ^ X )+F+i (ф ( ^ X ))} ■ (4)
Решая последовательно уравнение Беллмана (4) на каждом /-м шаге, по-
Ю7—т'Н 7—т'Н 7—т'Н 7—т'Н 7—т'Н
—^ к-\ —^ к-2 —^ ■ ■ ■ —^ 2 —^ 1 и условное оптимальное управление X* —ч —_2 —>...—> —Далее остаётся по известному
начальному вектору параметров состояния S0 определить состояние Si и найти минимум F1 и оптимальное управление X1 и так далее по цепочке, находя оптимальные значения параметров Tg0i, Wgi на каждом i-м этапе декомпозиции, получим доминирующую последовательность {Xi(Tg0i, Wgi)} f=1, формирующую динамический режим, оптимизирующий ТТП высушивания многослойной массы окатышей по всей зоне сушки ТТС обжиговой машины.
На каждом i-м этапе неявной декомпозиции локальный КЭ Fi представляет собой функцию двух переменных F(Tg0i, Wgi) = ai й (Tgoi, Wgi) + aiP(Tgoi, Wgi), для которой требуется определить точку минимума X*(Tg0i*, Wgi*). Для решения этой задачи условной оптимизации, применялся метод деформируемого многогранника (Нелдера-Мида). Выполнение ограничений на управляющие параметры - скорость 0 < Wgi < WgMAX и температуру 0 < Tg0i < TgMAX теплоносителя на входе в слой достигается введением в функцию Fi(Tg0i, Wgi) барьерных функций, так как это ограничение является наиболее жестким, ибо оно обеспечивает математически-корректное и физически-осуществимое решение задачи [14]. Выполнение остальных условий: Tghi < TghMAX, u < uMAX, I <IMAX, ЭГт/Эх , ЭТт/Эх < , осуществляется методом скользящего допуска и штрафных функций типа «квадрата срезки», с введением вспомогательного аргумента, обеспечивающего одинаковый порядок величин входящих в систему ограничений.
Вычислительная сложность разработанного алгоритма оптимизации локального КЭ Fi в процедуре оптимизации ТТП сушки движущейся плотной многослойной массы окатышей соответствует 0(G3), где G - количество ограничений на управляющие переменные и параметры состояния.
При решении двухкритериальной задачи оптимизации, учитываются сведения о важности частных критериев й (Tg0, Wg) и Q(Tg0, Wg). В данном случае эти частные критерии строго ранжированы по важности: во-первых, необходимо обеспечить завершённость ТТП сушки й ^0, во-вторых, при этом минимизировать удельный расход энергии Q^-min. Поэтому применялся метод последовательных уступок: оптимизировался первый, наиболее важный критерий, затем назначалась величина допустимого отклонения значения этого критерия, и оптимизировался второй по важности критерий, при условии, что значение первого критерия не должно отличаться от минимального более чем на величину установленного отклонения (уступки).
Анализ результатов оптимизации ТТП
На основе разработанных автором алгоритмов динамического программирования и оптимизации, программно-реализованных в компьютерных моде-
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
лях в среде Borland C++, проводились многочисленные вычислительные эксперименты по оптимизации ТТП многослойной сушки фосфоритовых окатышей.
В оптимальном режиме сушки скорость газа-теплоносителя на протяжении всего процесса остается постоянной максимально-возможной 1,3 м/с. Температура газа-теплоносителя медленно растет, скорость ее изменения обусловлена ограничениями на скорость нагрева, градиент температур, долю переувлажнения и интенсивность влагопереноса в окатыше рис. 4.
Время нахождения в зоне сушки т, %
Рис. 3. График искомой функции Тё0(т) температуры газа-теплоносителя на входе в слой и 7g/,(i) на выходе из слоя...........
Таким образом, формируется режим, при котором температура теплоносителя умеренно растет, интенсифицируя процесс сушки в среднем по слою, и не допуская переувлажнения материала [15]. Причем газ-теплоноситель выносит максимально-возможное (зависящее от свойств теплоносителя и материала окатышей) количество влаги из слоя. Дальнейшее повышение температуры не несет значительного уменьшения среднего по слою влагосодержания, а расход тепловой энергии увеличивается [16]. Полученные результаты, применялись для расчета оптимального режима сушки фосфоритовых окатышей на обжиговой машине конвейерного типа ОК-520/536Ф.
Заключение
В результате проведенных исследований: представлена содержательная и математическая постановки задачи оптимизации ТТП сушки движущейся плотной многослойной массы фосфоритовых окатышей в конвейерной обжиговой машине, разработаны математическая и компьютерная модели оптимизации энергоресурсоэффективности ТТС, отличающиеся учетом тепло-технологических особенностей обжиговой машины, интенсивности процесса внутреннего влагопереноса в окатыше и процессов переувлажнения отдельных слоев окатышей и газа-теплоносителя, что позволяет снижать энергозатраты, по-
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
средством интенсификации ТМП многослойной сушки. Анализ результатов полученного оптимального режима ТТП сушки показал, что в зоне сушки основная часть влаги удаляется раньше, чем материал поступает в зону подогрева и обжига и остаточная влага не в состоянии разрушить окатыш даже при интенсивном, но допустимом нагреве. Проводя сушку в оптимальном режиме, стоимость израсходованной энергии можно сократить в среднем на 8%.
В данной работе для математического описания ТТП сушки и оптимизации КЭ, используются системы ДУЧП, решение которых реализуется численными методами по разработанной автором компьютерной модели. Однако такой подход не позволяет учесть присущую теплофизическим характеристикам этих процессов неопределенность. Поэтому актуальной научной задачей является исследование этих процессов методами нечеткой математики. Итерационные вычисления над нечеткими переменными конечно-разностных уравнений приводят к проблеме возрастания (накопления) неопределенности и к «размыванию» результатов, вследствие чего они могут оказаться сложно интерпретируемыми и непригодными для дальнейшего анализа. В перспективе предлагается реализовать подход к решению ДУЧП в условиях неопределенности параметров ТТП сушки с использованием нечетких численных методов, который позволит исследовать распределение температуры сферических окатышей в движущейся на конвейере ТТС плотной многослойной массе в условиях неопределенности теп-лофизических характеристик.
Предлагаемый нечеткий подход позволит сформулировать содержательную постановку задачи нечеткой оптимизации КЭ для сложных многостадийных ТТС и определить пути решения этой современной научной и важной практической задачи.
Работа выполнена при финансовой поддержке РФФИ, грант 17-01-00189_а.
Литература
1. Мешалкин В. П., Бобков В. И. Ресурсосберегающие энергоэффективные технологии обработки фосфатного сырья // ХХ Менделеевский съезд по общей и прикладной химии. - Екатеринбург, 2016. - С. 299.
2. Леонтьев Л. И. Физико-химические особенности комплексной переработки железо-содержащих руд и техногенных отходов // ХХ Менделеевский съезд по общей и прикладной химии. - Екатеринбург, 2016. -С. 92.
3. Юсфин Ю. С., Пашков Н. Ф., Антоненко Л. К., Жак Р. М., Майзель Г. М., Базилевич Т. Н. Интенсификация производства и улучшение качества окатышей. -М.: Металлургия, 1994. - 240 с.
4. Luis P., Van der Brüggen B. Exergy analysis of energy-intensive production processes: advancing towardsa sustainable chemical industry // Journal of Chemical Technology and Biotechonology. 2014. Vol. 89. № 9. P. 1288-1291.
5. Elgharbi S., Horchani-Naifer K., Férid M. Investigation of the structural and mineralogical changes of Tunisian phosphorite during calcinations // Journal of Thermal Analysis and Calorimetry. 2015. Vol. 119. № 1. P. 265-269.
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
6. Bobkov V. I., Borisov V. V., Dli M. I., Meshalkin V. P. Multicriterial optimization of the energy efficiency of the thermal preparation of raw materials // Theoretical Foundations of Chemical Engineering. 2015. Vol. 49. № 6. P. 842-846.
7. Bobkov V. I., Borisov V. V., Dli M. I., Meshalkin V. P. Modeling the calcination of phosphorite pellets in a dense bed // Theoretical Foundations of Chemical Engineering. 2015. Vol. 49. № 2. P.176-182.
8. Rudobashta S. P. Calculation of the kinetics of drying disperse materials on the basis of analytical metods // Journal of Engineering Physics and Thermophysics. 2010. Vol. 83. № 4. P. 753-759.
9. Бобков В. И. Моделирование технологических процессов при термической подготовке дисперсного фосфатного сырья // Химическая технология. 2016. № 6. С. 263-271.
10. Panchenko S. V., Shirokikh T. V. Thermal Hydraulics of moving dispersive layer of process units // Theoretical Foundations of Chemical Engineering. 2016. Vol. 50. № 2. P. 217-224.
11. Panchenko S. V., Shirokikh T. V. Thermophysical processes in burden zone of submerged arc furnaces // Theoretical Foundation of Chemical Engineering. 2014. Vol. 48. № 1. P.77-81.
12. Meshalkin V. P., Men'shikov V. V., Panchenko S. V., Panchenko D. S., Kazak A. S. Computer-aided simulation of heat- and mass-transfer processes in an orereduction electrothermal reactor // Theoretical Foundations of Chemical Engineering. 2015. Vol. 49. № 5. P. 55-60.
13. Бобков В. И., Борисов В. В., Дли М. И. Подход к исследованию теплопроводности нечеткими численными методами в условиях неопределенности теплофизических характеристик // Системы управления, связи и безопасности. 2017. № 3. С. 73-83.
14. Meshalkin V. P., Kolesnikov V. A., Desyatov A. V., Milyutina A. D., Kolesnikov A. V. Physicochemical efficiency of electroflotation of finely divided carbon nanomaterial from aqueous solutions containing surfactants // Doklady Chemistry. 2017. Vol. 476(1). P. 219-222.
15. Pancnehko S. V., Dli M. I., Borisov V. V., Panchenko D. S. Analysis of thermalphysic processes in near-electrode zone of electrothermal reactor // Non-ferrous Metals. 2016. Vol. 2. P. 57-64.
16. Panchenko S. V., Meshalkin V. P., Dli M. I., Borisov V. V. Computervisual model of thermophysical processes in electrothermal reactor // Tsvetnye Metally. 2015. Vol. 4. P. 55-60.
References
1. Meshalkin V. P., Bobkov V. I. Resource-saving energy efficient technologies of processing of phosphatic raw materials. XX Mendeleevskii s"ezdpo obshchei i prikladnoi khimii [XX Mendeleev congress on general and applied chemical]. Ekaterinburg, 2016, p. 299. (in Russian)
2. Leont'ev L. I. Physical and chemical features of the complex processing of iron-stones and technogenic wastes. XX Mendeleevskii s"ezdpo obshchei iprikladnoi
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
khimii [XX Mendeleev congress on general and applied chemical]. Ekaterinburg, 2016, p. 92. (in Russian)
3. Iusfin Iu. S., Pashkov N. F., Antonenko L. K., Zhak R. M., Maizel' G. M., Bazilevich T. N. Intensifikatsiia proizvodstva i uluchshenie kachestva okatyshei [Intensification of production and improvement of quality of pellets]. Moscow, Metallurgiia, 1994. 240p. (in Russian)
4. Luis P., Van der Bruggen B. Exergy analysis of energy-intensive production processes: advancing towardsa sustainable chemical industry. Journal of Chemical Technology andBiotechonology, 2014, vol. 89, no. 9. pp. 1288-1291.
5. Elgharbi S., Horchani-Naifer K., Ferid M. Investigation of the structural and mineralogical changes of Tunisian phosphorite during calcinations. Journal of Thermal Analysis and Calorimetry, 2015, vol. 119, no 1, pp. 265-269.
6. Bobkov V. I., Borisov V. V., Dli M. I., Meshalkin V. P. Multicriterial optimization of the energy efficiency of the thermal preparation of raw materials. Theoretical Foundations of Chemical Engineering, 2015, vol. 49, no. 6, pp. 842-846.
7. Bobkov V. I., Borisov V. V., Dli M. I., Meshalkin V. P. Modeling the calcination of phosphorite pellets in a dense bed. Theoretical Foundations of Chemical Engineering, 2015, vol. 49, no.2, pp.176-182.
8. Rudobashta S. P. Calculation of the kinetics of drying disperse materials on the basis of analytical metods. Journal of Engineering Physics and Thermophysics, 2010, vol. 83, no. 4, pp. 753-759.
9. Bobkov V. I. Modelirovanie tekhnologicheskikh protsessov pri termicheskoi podgotovke dispersnogo fosfatnogo syr'ia [Design of technological processes at thermal preparation of dispersible phosphatic raw material]. Chemical Engineering, 2016, no. 6, pp. 263-271 (in Russian).
10. Panchenko S. V., Shirokikh T. V. Thermal Hydraulics of moving dispersive layer of process units. Theoretical Foundations of Chemical Engineering, 2016, vol. 50, no. 2, pp. 217-224.
11. Panchenko S. V., Shirokikh T. V. Thermophysical processes in burden zone of submerged arc furnaces. Theoretical Foundation of Chemical Engineering, 2014, vol. 48, no. 1, pp. 77-81.
12. Meshalkin V. P., Men'shikov V. V., Panchenko S. V., Panchenko D. S., Kazak A. S. Computer-aided simulation of heat- and mass-transfer processes in an ore-reduction electrothermal reactor. Theoretical Foundations of Chemical Engineering, 2015, vol. 49, no. 5, pp. 55-60.
13. Bobkov V. I., Borisov V. V., Dli M. I. Approach to a Heat Conductivity Research by Fuzzy Numerical Methods in the Conditions of Indeterminacy Thermal Characteristics. Systems of Control, Communication and Security, 2017, no. 3, pp. 73-83. (in Russian).
14. Meshalkin V. P., Kolesnikov V. A., Desyatov A. V., Milyutina A. D., Kolesnikov A. V. Physicochemical efficiency of electroflotation of finely divided carbon nanomaterial from aqueous solutions containing surfactants. Doklady Chemistry, 2017, vol. 476(1), pp. 219-222.
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
15. Pancnehko S. V., Dli M. I., Borisov V. V., Panchenko D. S. Analysis of thermalphysic processes in near-electrode zone of electrothermal reactor. Non-ferrous Metals, 2016, no. 2, pp. 57-64.
16. Panchenko S. V., Meshalkin V. P., Dli M. I., Borisov V. V. Computervisual model of thermophysical processes in electrothermal reactor. Tsvetnye Metally, 2015, no. 4, pp. 55-60.
Статья поступила 3 апреля 2018 г.
Информация об авторе
Бобков Владимир Иванович - кандидат технических наук. Доцент кафедры высшей математики. Национальный исследовательский университет «Московский энергетический институт», филиал в г. Смоленске. Области научных интересов: системный анализ сложных теплотехнологических систем; оптимизация теплотехнологических процессов; математическое моделирование; интеллектуальная поддержка принятия решений; интеллектуальный анализ данных. E-mail: vovabobkoff@mail.ru
Адрес: 214013, Россия, г. Смоленск, Энергетический проезд, д. 1.
Системы управления,связи и безопасности №2. 2018
Systems of Control, Communication and Security sccs.intelgr.com
The heat-technological process optimization of a moving dense multilayer weight of phosphorite pellets drying process by the energy-
resource efficiency criterion
V. I. Bobkov
Statement of the problem: the problem of fuel and energy resources rational use is the most important one for energy-intensive industrial production, including the thermal preparation of pelletized raw materials in a complex multi-stage heat-technological system (HTS) by a horizontal-grate machine for the phosphorite pellets production. The drying modes optimization of a moving dense multilayer weight ofphosphorite pellets in a complex HTS will increase the energy and resource efficiency of HTS as a result of heat-mass exchange processes (HMP) augmentation and optimization at an energy and heat minimum cost. The aim of the work is to increase the energy-resource efficiency of HTS, using the mathematical and computer model, the heat-technological process (HTP) optimization of the moving dense multilayer weight of phosphorite pellets drying process, characterized by the use of discrete dynamic programming and the control of the moisture redistribution effect by the height of the pellets multilayer backfilling, which allows to prevent the appearance of waterlogging horizons, negatively affecting the multilayer pellets weight gas permeability, that leads to unnecessary energy consumption increase and reduction of HTS performance in general. The methods used: system concept, discrete dynamic programming, optimal control theory of distributed parameters system, conditional multi-criteria optimization, computational mathematics. The result: the content and mathematical statements of the optimization problem (HTP) for moving dense multilayer weight of phosphorite pellets drying process in the complex multistage (HTS) horizontal-grate machine that is characterized by the heat-technological features of this horizontal-grate machine, the internal moisture transfer process intensity in the pellet and the waterlogging of sepatare pellets layers processes and gas-coolant, which allows to increase energy efficiency through intensification of dynamic multilayer drying HMP were submitted. Novelty: mathematical and computer models for optimizing the a moving dense multi-layer weight of phosphorite pellets drying HTP have been developed, which have made it possible to optimize the energy and resource efficiency of a complex dynamic phosphate pellet production HTS. The obtained results were used to calculate the energy-efficient drying of pellets in the HTS by a horizontal-grate machine. It was also found out that there is no waterlogging zone in the optimal mode of multilayer pellets drying, that HMP are intensified, the energy consumption is reduced and the end-product quality is increased, the share of return is reduced, that provides resource saving.
Keywords: heat-technological system, heat-technological process, drying, optimization, pellets, phosphorite, energy-resource efficiency, horizontal-grate machine.
Information about Author Vladimir Ivanovich Bobkov - Ph.D. of Engineering Sciences, Associate Professor of the Dept of Mathematics. The Branch of National Research University "Moscow Power Engineering Institute" in Smolensk. Fields of research: system analysis of complex heat engineering systems; optimization of heat engineering processes; math modeling; intellectual decision-making support; data mining. E-mail: vovabob-koff@mail.ru
Address: Russia, 214013, Smolensk, Energeticheskiy proezd, 1.
