УДК 543.878:541.127:519.2
В. Р. Пен, И. В. Ковалев
УПРАВЛЕНИЕ ПРОЦЕССОМ ДЕЛИГНИФИКАЦИИ СМЕШАННОГО ДРЕВЕСНОГО СЫРЬЯ
Проанализировано применение Н-фактора при варке древесины. Показана ограниченность этого критерия и невозможность его использования при варке смесей различных пород древесины. Предложена математическая модель кинетики процесса делигнификации смеси различных пород древесины и обосновано ее применение для управления варкой.
Практическая работа варочных цехов целлюлозных заводов зачастую связана с необходимостью переработки смесей различных древесных пород. При этом возникает ряд трудностей, обусловленных их различным кинетическим поведением в делигнификационных процессах. Особенно сильно совместная варка влияет на качество целлюлозы, предназначенной для химической переработки. Имея все параметры (содержание б-целлюлозы, жесткость, содержание остаточного лигнина ит.д.) в пределах стандарта, она ведет себя аномально и к последующей химической переработке оказывается непригодной. Такое поведение объясняется тем, что вся масса целлюлозы состоит из дискретных волокон, получаемых после растворения скрепляющего их лигнина. В условиях совместной варки одна часть оказывается переваренной, другая - недоваренной. Имея стандартные средние показатели, эта целлюлоза для каждого отдельного волокна дает значительную дисперсию значений контролируемых параметров. Попытка измерить их наталкивается на труднопреодолимую проблему разделения отличающихся по характеристикам целлюлозных волокон для последующего анализа. Дополнительная сложность при оценке результата делигнификации может быть обусловлена наличием синергетического эффекта, когда совместная варка дает продукт, параметры которого отличаются от тех, которые целлюлоза бы имела, будучи полученной из сваренных порознь древесных пород с последующим их смешением. Кроме того, варочный котел является очень инерционной системой, обладающей большим временем запаздывания. В сочетании с сильным экзотермическим эффектом реакции делигнификации, зависящим от скорости процесса, все эти причины приводят к тому, что температурный график варки сильно колеблется. Эти колебания температуры вызывают соответствующие колебания скорости варки, затрудняя определение момента ее окончания.
Таким образом, определение критерия, позволяющего оценивать степень завершенности процесса делигни-фикации, есть важнейшая задача управления этим процессом. Традиционный способ состоит в использовании Н-фактора, являющегося, по сути, интегральным показателем [1; 2].
На использовании интегральных показателей построен распределенный контроль объектов [3]. В частности, как параметр состояния объекта используется его среднеинтегральная инерционность Би.:
1 ^
Би = -| £(х)g(х)ёх, (1)
Ь 0
где Би - инерционность в точке х; g(x) - весовая функция; Ь - длина аппарата.
По уравнению типа (1) оценивается комплексный параметр - среднеинтегральная температура варочного реактора [1], комплексные параметры тел сложных сечений, параметры графических образов, параметры поверхностей отклика, параметры микрообъектов и т. д. [4].
Интегральные характеристики являются информационными критериями-комплексами, дающими количественную характеристику некоторых свойств и особенностей технологических процессов.
Основа интегрального параметра - интеграл вида
ь
I х(0gЦ)Ж, (2)
а
где х(г) - исследуемая функция; g(г) - весовая функция, вид которой задается в соответствии с классом явления исследуемого объекта.
Для более точного соответствия истинным значениям числовой характеристики интегральный параметр может потребовать специального формирования весовой функции процесса. Так, для оценки температурновременного режима процесса варки целлюлозы используются интегральные параметры:
9(1 )-9р
Е =[ п 10 Ж.
. (3)
. В м А---
Н =| е 9(') Ж.
к
где п - температурный коэффициент, показывающий, во сколько раз увеличивается скорость процесса варки при повышении температуры на 10 °С; 9(т) - текущее значение температуры процесса в момент времени I; 90 - базовая температура, °С; е - основание натуральных логарифмов; А, В - константы; го - начальное время нагрева (допускается t0 = 0); гк - конечное время нагрева. Интегральные параметры представляют собой комплексную оценку переменных группы С. Они обладают свойством аддитивности, т. е. при разбиении области параметров, характеризующих данный температурно-временной режим, на ряд подобластей (С1, С2, ..., Ст) общий интегральный параметр будет равен сумме частных интегральных характеристик:
Е (Со) = Е (С!) + Е (С2) +... + Е(Ст). (4)
Это свойство собственного интеграла позволяет использовать интегральные параметры для решения задач контроля, управления и оптимизации температурно-временных режимов процесса варки целлюлозы [1; 5].
Математика, механика, информатика
При варке древесины для получения необходимой точности при оценке температурных графиков технологических процессов применяют Н-фактор, который определяют с помощью уравнения Аррениуса по формуле
н = |і-
(5)
В условиях производства очень часто длительность технологических процессов изменяется. В этом случае для некоторого эталонного режима можно подсчитать Н-фактор и принять его за критерий, к которому следует стремиться. В случае отклонения технологического режима от заданного при постоянных показателях качества щепы и варочного щелока и их соотношениях можно обеспечить требуемое качество продукции, если продолжительность и температуру варки изменять таким образом, чтобы Н-фактор был близок к выбранному критерию. Для случая, когда верхний предел Н-фактора - переменная величина, вводится понятие Н-функции, которая может быть представлена в виде
н = 1 к
іс
Фактически, учитывая, что
_ Е
ке КТ (і) = к
(6)
(7)
->Ь„
Ь-
кинетическое уравнение имеет вид ЖЬ
Ж
■ = -К,Ь .
Соответственно
Ь = Ьсе іс
Ь = Ьсе
_І Кь (і)Л
(8)
(9)
(10)
(11)
Полагая, что
соответственно
ёЬ
Ж
ЖЬв
Ж
'А =-КАЬА,
= _ КЬ.
Ь = ЬА + Ьв
ёЬ
Ж
Ь = Ьс е
-| К (і)ё/
(14)
(15)
(16)
= Ьсе
Н- фактор есть интеграл по времени от константы скорости, которая в процессе варки зависит от времени.
Такой способ управления процессом за длительное время своего существования показал свою эффективность и удобство использования. Однако являясь по сути приближенным, в случае необходимости управления варкой, предназначенной для получения полуфабрикатов с точно заданными качественными показателями, приводит к погрешностям. При этом суть приближения состоит в использовании для описания кинетики делигнификации кинетического уравнения первого порядка.
Предполагая, что процесс варки идет по схеме
Причем, так как должно выполнятся условие
А)е~Н = ЬА0е-НА • ЬВ0е~Н , (17)
а из уравнения (17) следует, что
1,е-Н = ЬА0 е-Н + Ьв 0 е-Н. (18)
то требование аддитивности критерия Н не выполняется даже в случае простейшей кинетики первого порядка. И тем более оно не выполняется для процессов с более сложным кинетическим механизмом и процессов с полихрон-ной кинетикой. Поэтому при совместной варке древесины различных пород Н-фактор не может быть использован в качестве параметра, по которому оценивается степень завершенности процесса.
В качестве последнего естественно использовать степень конверсии лигнина в лигноуглеводной матрице. Для полихронных процессов она определяется из выражения
10(к, і) / (к )Жк
0ь (і) =
Ь
(19)
и формально соответствует главному требованию, предъявляемому к интегральным показателям - аддитивности. Однако невозможность явного введения в выражение (20) температуры обуславливает необходимость использования математического описания в виде системы дифференциальных уравнений. При этом управление варкой смешанного сырья ведется по модели, являющейся аддитивным объединением систем уравнений, описывающих кинетику делигнификации каждой породы по отдельности и объединенных общей концентрацией щелочи и общим уравнением ее материального баланса:
ёЬ,
Ж
= -Ь. . он • (Кь, .) +Кьк • АН0;
Таким образом, Д-фактор является просто эмпирическим параметром, реально не учитывающим кинетические особенности процесса делигнификации. Особенно важно, что, будучи аддитивным по времени
Н(іс;і,) + Н(і,;і2) = Н(іс;і2), (12)
Д-фактор не аддитивен по породному составу. То есть при совместной варке древесины различных пород
Н1 + Н2 Ф Н12. (13)
Последнее обусловлено тем, что для двух различных пород А и В
С (>, 1) <. 1 АН<2 = А^0 • ехр(-КАв ■ ();
ОН = ОН0 -аг(Ь + С0)- (20)
-а^(Ь - А) - ас,1 • (С“ - С1) --а2 • (Ь\ + С0) -„ • (4 - £2) - а^ • (С20 - С2);
Ь0 = У1°.; Ь. = У Ь. .; С0 =
1 г
= УС0.; С. = УС...
/ г
где j = 1 - лиственница, j = 2 - сосна. Ансамбли имеют индекс г = 1...3; Ь - лигнин, С - углеводы. При AНQ = 0,
AQ = 0 система (21) описывает натронную и сульфатную варку, при AQ > 0 - натронно-антрахинонную. Балансовое уравнение является суммой балансовых уравнений для сосны и лиственницы. Константы, входящие в модель, определяются при моделировании варок отдельных древесных пород.
По результатам расчетов построены кинетические кривые для совместной и раздельной варок сосны и лиственницы при температуре 445 °К (рис. 1, 2). Для сосны при совместной варке степень превращения компонентов выше, чем при раздельной. Для лиственницы - наоборот. Для натронно-антрахинонных варок характер кинетических кривых аналогичен. Изменение разности ДИ степени конверсии лигнина и углеводов совместной и раздельной варок (рис. 3) протекает в следующей динамике:
Д9 = 9 -9 . (21)
разд совм )
Для лигнина эта разница достигает наибольшего по абсолютной величине значения спустя 1 000 .1 500 с после начала варки, когда 9 = 0,5 .0,6.
В дальнейшем различия нивелируются, так как степень конверсии приближается к единице. Для углеводов кривая Д9 - Т монотонна. Характер ее изменения у разных пород практически одинаков.
о г 3 б~ 5 « Й"
я:» ю"*.е
Рис. 1. Зависимость степени конверсии от времени. Сосна: пунктир - совместная варка; сплошная - раздельная; точки - эксперимент: 1 - лигнин, 2 - углеводы
| Г 3 £ I to £"
■гно'е
Рис. 2. Зависимость степени конверсии от времени. Лиственница: пунктир - совместная варка; сплошная - раздельная; точки - эксперимент.
1 - лигнин, 2 - углеводы
Таким образом, можно констатировать, что традиционно применяемый Н-фактор не позволяет эффективно управлять процессом делигнификации смешанного древесного сырья. В то же время использование объединенной математической модели процесса с вычислением 9 обеспечивает надежное управление варкой и получение волокнистых полуфабрикатов с требуемыми параметрами.
лиственницы для натронной и натронно-антрахинонной
варок; 2 - лигнин лиственницы при натроннолй варке;
3 - лигнин лиственницы при натронно-антрахинонной варке;
4 - лигнин сосны при натронно-антрахинонной варке;
5 - лигнин сосны при натронной варке; 6 - углеводы сосны при натронной и натронно-антрахинонной варке
Библиографический список
1. Вьюков, И. Е. Математические модели и управление технологическими процессами целлюлозно-бумажной промышленности / И. Е. Вьюков, И. Ф. Зорин, В. П. Петров. М., 1975. 376 с.
2. Петров, В. П. Автоматизированные системы управления технологическими процессами / В. П. Петров, И. Ф. Зорин, С. А. Рогульская. М., 1977. 51 с.
3. Девятов, Б. Н. Теория переходных процессов в технологических аппаратах с точки зрения задач управления / Б. Н. Девятов. Новосибирск, 1964. 323 с.
4. Зорин, И. Ф. Контроль переменных при управлении процессом периодической варки целлюлозы / И. Ф. Зорин, В. П. Петров // Автоматизация, механизация, приборное и аппаратурное оформление процессов целлюлозно-бумажного производства : сб. тр. ВНИИБ. Вып. 71. Л., 1977. С. 10-18.
5. Расчет и оптимизация графика нагрева для группы варочных котлов периодического действия / В. Н. Казанский [и др.] // Автоматизация, механизация и оборудование процессов ЦБП : сб. тр. ВНИИБ. Л., 1977. С. 33-36.
V R. Pen, I. V. Kovalev
MANAGEMENT OF DELIGNIFICATION PROCESS OF THE MIXED WOOD RAW MATERIAL
Application of the H-factor at cooking wood is analysed. Limitation of this criterion and impossibility of its use are shown at cooking mixes of various breeds of wood. The mathematical kinetic model of process of delignification to a mix of various breeds of wood is offered and its application for managment of cooking is proved.