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

доктора химических наук
Пен, Владимир Робертович
город
Санкт-Петербург
год
2010
специальность ВАК РФ
05.21.03
Автореферат по технологии, машинам и оборудованию лесозаготовок, лесного хозяйства, деревопереработки и химической переработки биомассы дерева на тему «Кинетика делигнификации древесины»

Автореферат диссертации по теме "Кинетика делигнификации древесины"

Пен Владимир Робертович

КИНЕТИКА ДЕЛИГНИФИКАЦИИ ДРЕВЕСИНЫ

05.21.03 - Технология и оборудование химической переработки биомассы дерева; химия древесины

Автореферат диссертации на соискание ученой степени доктора химических наук

2 8 ОКТ 7010

Санкт-Петербург 2010

004611923

Работа выполнена в Сибирском государственном технологическом университете

Официальные оппоненты:

Доктор химических наук, профессор Дейнеко Иван Павлович

Доктор химических наук, профессор Боголицын Константин Григорьевич

Доктор технических наук, профессор Вураско Алеся Валерьевна

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

Защита состоится 26 октября 2010 г. в н часов на заседании диссертационного совета Д 212.220.01 в Санкт-Петербургской государственной лесотехнической академии им. С.М.Кирова по адресу: 194021, Санкт-Петербург, Институтский пер., д.5

С диссертацией можно ознакомиться в библиотеке Санкт-Петербургской государственной лесотехнической академии им. С.М.Кирова

Отзывы на автореферат в двух экземплярах с подписями, заверенными печатью, просим направлять по адресу: 194021, г.Санкт-Петербург, Институтский пер., д. 5, Санкт-Петербургская государственная лесотехническая академия им.С.М.Кирова, Ученый Совет

Автореферат разослан Я сентября 2010 г.

Ученый секретарь диссертационного совета, доктор технических наук, профессор 1Л\J7jyyVl~ Калинин

Актуальность проблемы. Теория промышленных способов делигнификации древесины составляет обширный и важный раздел химии растительного сырья. Многолетними усилиями большого числа исследовательских коллективов и отдельных ученых воссоздана относительно полная картина химических превращений компонентов древесины при щелочных, сульфитных и окислительных методах производства технической целлюлозы. В значительно меньшей степени разработана химическая кинетика процессов делигнификации. Одной из причин этого является недостаточная информативность методов классической химической кинетики при исследовании такого «неклассического» объекта, как древесина с её сложной макромолекулярной структурой и анизотропией свойств. Разработанные в диссертации новые методы и предложенные идеи, такие, как математическая теория полихронной кинетики делигнификации, применение статистического моделирования методом Монте-Карло для оценки влияния топологии макромолекулы лигнина на кинетику процесса деструкции, метод получения аппроксимационных кинетических уравнений и метод оценки кинетических параметров, в сочетании с мощными возможностями современной вычислительной техники, позволили существенно продвинуться вперед в этой области.

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

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

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

На основе представлений о полихронности кинетики делигнифика-ции в диссертации проведен анализ четырех способов варки целлюлозы -натронного, сульфатного, натронно-антрахинонного, окислительного. Реакции компонентов древесины с гидроксидом натрия протекают, как независимые стадии при всех щелочных способах делигнификации. Кинетические параметры натронной варки используются при моделировании других щелочных способов получения целлюлозы. Сульфатная варка в настоящее время является основным промышленным способом делигнификации, на её долю приходится более 90 % мирового производства целлюлозы. Вместе с тем, значительная техногенная нагрузка на окружающую среду, обусловленная применением серосодержащих соединений, стимулирует поиски более экологически безопасных способов делигнификации. Натронная варка с добавками каталитических количеств антрахинона уже успешно реализована на ряде предприятий и считается одним из наиболее перспективных бессерных способов производства целлюлозы. По этой же причине в качестве альтернативы сульфатной варке рассматриваются окислительные способы делигнификации. В диссертации приведены результаты кинетического исследования одного из новых перспективных процессов -окисления лигнина надуксусной кислотой, генерируемой в системе «пе-роксид водорода - уксусная кислота - вода - катализатор».

Исследования проводились по Государственной Hill «Комплексное использование и воспроизводство древесного сырья» (1992-1997 г.г.) и Федеральной целевой НТП «Исследования и разработки по приоритетным направлениям развития науки и техники гражданского назначения», подпрограмма «Комплексное использование древесного сырья» (1998-2001 г.г.).

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

Задачи исследования:

- с использованием метода Монте-Карло и теории графов создать стохастическую имитационную модель деструкции и проанализировать влияние топологической структуры и химической неоднородности макромолекул лигнина на кинетику процесса делигнификации;

- разработать математическую теорию полихронной кинетики процессов делигнификации древесины;

- выполнить экспериментальное исследование кинетики щелочных и окислительных методов делигнификации;

- получить математические модели кинетики изученных процессов делигнификации древесины;

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

- разработать метод кинетического анализа процесса варки различных древесных пород для оценки степени неоднородности получаемой целлюлозы;

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

Научная новизна. Автором впервые получены следующие научные результаты.

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

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

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

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

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

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

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

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

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

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

Книги «Математический аппарат исследования структуры и деструкции полимеров», «Кинетика делигнификации древесины» и «Теоретические основы делигнификации» используются в качестве учебных пособий студентами СибГТУ и родственных вузов.

Основные положения, выносимые на защиту. Автор защищает:

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

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

- математическую теорию полиэкспоненциальной кинетики процесса делигнификации;

- метод получения математических моделей химической кинетики для систем с распределенными кинетическими параметрами;

- метод оценки параметров математичеких моделей, описывающих кинетику процессов в системах с полиэкспоненциальной кинетикой;

- кинетическую теорию и математическую модель окислительной делигнификации;

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

Апробация работы. Результаты исследований докладывались и обсуждались: на Всесоюзной конференции «Теоретические проблемы новых методов делигнификации» (Братск, 1985); на Всесоюзной конференции «Современные проблемы химической технологии» (Красноярск, 1986); на 7-й Всесоюзной конференции по химии и использованию лигнина (Рига, 1987); на 6-й Всесоюзной конференции по физике и химии целлюлозы (Минск, 1990); на международной конференции «Ьапх-98» (Красноярск, 1998); на Всероссийской конференции с международным участием «Дос-

тижения науки и техники - развитию сибирских регионов» (Красноярск, 1999); на региональных конференциях «Проблемы химико-лесного комплекса» (Красноярск, 1998-2009); на Международной конференции «Научные исследования высшей школы по приоритетным направлениям науки и техники» (Савонна-Барселона-Тенерифе-Савонна, 2009); на Всероссийской заочной электронной научной конференции «Новые материалы и химические технологии» (Вебсайт Российской академии естествознания www.rae.ru. 2009); на III Общероссийской конференции «Актуальные вопросы науки и образования» (Москва, 2010).

Публикации. По теме диссертации опубликовано 42 работы. В их числе 3 монографии, 32 статьи, 7 тезисов докладов.

Объем и структура диссертации. Диссертация состоит из введения, девяти глав основного текста, выводов, библиографического списка (440 наименований) и приложения, включает 64 таблицы, 67 рисунков, изложена на 390 страницах.

Содержание работы

1 Анализ представлений в области строения и кинетики деструкции полимеров древесины

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

а) усложнение математического описания с включением большого числа переменных факторов и подгоночных параметров;

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

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

о

где п(т) - число разрушенных связей к моменту г, щ - начальное общее число связей; ]{к) - функция распределения связей по константам скоро-

(1.1)

ста кг, 0(пй,к, т) - решение кинетического уравнения для отдельного типа связи в макромолекуле.

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

2 Развитие теории полихронной кинетики

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

а -

тх

(2.1)

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

КмАХ

с(г)= | С(к,1)/(к)с1к> (2.2)

Кмт

где есть решение соответствующего кинетического уравнения, а /(к) - некая функция распределения, учитывающая кинетическую неоднородность вступающего в реакцию вещества.

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

где *,.(г) — наблюдаемая эффективная константа скорости.

Кинетическая кривая, соответственно, описывается уравнением

с,(0= (2-4)

о

В случае реакции первого порядка с1с1 , , .

= (2.5)

и для полихрониых процессов справедливо следующее уравнение:

г

о

Фактически выражения (2.2) и (2.6) для реакций первого порядка — это две формы записи одного и того же явления, так как для (2.5) имеет место отношение

{*/(*) ехр(-А*)г/Л

т • (2.7)

Дня (2.3) справедливо более общее выражение

л

Кт . (2.8)

о

При этом предполагается, что У(к) от времени не зависит, т.е. скорость релаксации достаточна для восстановления _/(£) в ходе реакции, но недостаточна для выравнивания значений констант и приведения системы к виду, когда /(к) = ¿(¿).

Соответственно, среднее к определяется, как 00

к = _ (2.9)

-00

Перенумеруем локальные состояния реагентов в ансамбле индексом j = 1,2,N. Зависимость концентрации подансамблей в состоянии } от времени определяется уравнением

¿СЛ О

¿^=-I "Vе, (О+Е - *А<о=-I • (2-Ю)

Здесь \Vjj- - вероятность перехода подансамбля из состояния _/ в состояние _/' Ц - константа скорости химической реакции в подансамбле, находящемся в состоянии}. Два первых слагаемых описывают релаксацию,

тогда как последнее слагаемое описывает химическое превращение в продукты реакции. Суммируя уравнения (2.9), получаем:

Ткгм

— —^СДО; (2.11)

с(?)=£с,(о. (2Л2)

}=1

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

С(Г) = Х }соехр(^)/«^. (2.13)

|'=1 к,

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

С(0 = Е | С0 ехрН*)/(*)Л = ¿С0 скР(Ч,1)т(Кы - К,) (2-14) '=1 к, ы

(где К, - ;-е значение к на отрезке [0,&„<„]), или

С,(0 = С0,-ехр(-А;;г);

лг (2.15)

С(0=1С,(0-/=1

Если г, - доля начальной концентрации г-го ансамбля в общем количестве ансамблей, то

г.^тмш-к.у, (216)

С/о -

Знание функции плотности распределения /(к) позволяет из (2.2) найти вид уравнения для С((). Так как (2.2) по отношению к распределению /{к) является интегральным уравнением Фредгольма первого рода, то определение Дк) есть некорректная по Адамару задача, требующая специальных методов решения.

В общем случае необходимо рассматривать систему, для которой скорость релаксации может быть настолько незначительной, что исходное распределение не успевает восстанавливаться. Очевидно, что распределение подансамблей оказывается зависящим от времени и необходимо вводить плотность распределения Лк, (), для которой

С(/)= ¡С(к,0/(Шк. (2.17)

Киш

Для таких систем среднее' значение к оказывается зависящим от времени:

k(t)= ¡kf(k,t)dk_

(2.18)

В реальной системе подсистемы образуют непрерывное распределение, поэтому такая простейшая модель в общем случае не верна. При движении по поверхности потенциальной энергии в таких системах реализуется множество траекторий, каждая их которых характеризуется собственной константой скорости к. Можно ввести в рассмотрение пространство констант скоростей, движение в котором по множеству траекторий реализует случайный процесс к = (к,: /б5) случайных величин к,. Исходя из такого представления, для изменения константы скорости уравнения (2.5) можно записать:

^ = -ак + т ■£((), (2.19)

где ант- постоянные; ф) - «белый шум», то есть случайная функция, обладающая корреляционной функцией

*£(/) = £(*). (2.20)

Уравнению (2.19) соответствует второе уравнение Колмогорова

(2.21)

df д , .г 2 б/.

— = —(-akf -т —)

dt дк аг'

где а и т - константы; - двухвременная плотность распределе-

ния по константам скорости.

В диссертации получено решение этого уравнения;

Яа,;гА)=- 1

ехр[- 2о(г-01

•ехр

{-[¿г-£,ехр(-а(г-/))]2

2—[1 2а

-ехр(-2 я(г-О)]

(2.22)

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

но, что для этого периода можно постулировать неизменность /(г,кп во времени. Тогда д//3г = 0 и, соответственно,

ак дк

При этом

со

/(Т,к- 1,к.) =.№; \/(к)с1к = 1; Лк)>0,

-со

При интегрировании (2.23) получается уравнение

т1 — + <зЛ/" = О ¿к -1

(2.23)

решение которого имеет вид

/(*) = -

1

2 а

-ехр

(2.24)

(2.25)

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

(2.26)

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

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

с,(0=/с,(*,0Я*)<Л;

о

(2.27)

Су(0 = %{к,№к)с1к.

О !

Здесь Лк) — плотность распределения, зависящая от вектора к, компонентами которого являются п переменных к, определяемых числом п рассматриваемых в кинетической модели делигнификации веществ; С?,{&,/)

— решение соответствующего кинетического уравнения для /-го компонента древесины.

В соответствии с (2.2) уравнение кинетической кривой имеет вид:

Кшх 1

C« = Z J -exp

' ——Jlx

2a,

(k-kj

'2a,

dk,

(2.28)

Если используется традиционный подход, основанный на системе обыкновенных дифференциальных уравнений, то

(2.29)

>

dt

константа скорости Щ) зависит от времени и вычисляется по уравнению

и

*(0 =

2 а,

-exp

2a,

dt

2a,

-exp

{k-kj

2 a,

(2.30)

G(k„t)dk,

Применение выражений (2.13...2.16) дает возможность моделировать кинетику делигнификации в виде системы обыкновенных дифференциальных уравнений. Модель такого типа имеет традиционный для формальной кинетики вид и позволяет естественным образом использовать уравнение Аррениуса и строить зависимость эффективных констант скоростей от температуры.

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

Лк).

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

нетических ансамблей, так и константы скоростей соответствующих им реакций.

Суммарная площадь S всех N прямоугольных распределений равна

я = 1 ДА, •/(*,-), ! = U,iV, (2.31)

i

поэтому минимизируемый функционал F имеет вид

м

F = -/(¿,;)-l)2 ~> min 7 (2.32)

i=l i

где j = 1, ... , М соответствует числу уровней температуры, при которых проведены кинетические эксперименты.

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

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

3 Статистическое имитационное моделирование деструкции и делигнификации

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

В диссертации для анализа деструкции макромолекул с разными ва-

риантами их топологической структуры использовано моделирование методом статистических испытаний (метод Монте-Карло). Моделирование базируется на следующих исходных предпосылках:

- макромолекулы лигнина состоят из элементарных звеньев - фе-нилпропановых единиц (ФПЕ), соединенных разнообразными по химической природе и прочности связями;

- распределение связей разного типа в пределах макромолекулы случайное;

- наиболее вероятная топологическая структура макромолекул лигнина in situ - сетка с относительно редкими узлами ветвления и редкими внутримолекулярными «сшивками»;

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

- jr-меры, степень полимеризации которых оказывается меньше некоторого критического значения хкр, переходят в раствор (процесс делиг-нификации).

Моделирование процесса деструкции макромолекул разной топологической структуры включает следующие операции:

- генерацию графа, имитирующего заданную топологическую структуру макромолекулы, узлами которого являются фенилпропановые единицы, а ребрами - связи между ФПЕ;

- нумерацию ребер графа;

- многократное «разыгрывание» (случайный выбор) номеров удаляемых ребер (т. е. номеров разрываемых связей);

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

- построение дифференциальных кривых распределения фрагментов деструкции по молекулярным массам (кривых ММР).

Для выполнения имитационного моделирования разработана оригинальная компьютерная программа.

В ходе исследования были рассмотрены такие структуры: цепная (вариант А); дендровидная с разветвлениями от каждого узла (вариант Б); дендровидная с редкими разветвлениями (вариант В); дендровидно-сшитая (варианты Г и Д).

Цепная структура макромолекул (вариант А) не характерна для лигнина. Она представляет интерес по той причине, что её поведение при статистической (случайной) деструкции достаточно полно описано математически. Это предоставило возможность проверить применимость и информативность метода Монте-Карло при решении поставленной задачи. Выполненное исследование подтвердило хорошее соответствие между результатами статистического и аналитического моделирования процесса.

Дендровидная структура с разветвлением каждого узла (вариант Б) является предельным случаем графа-дерева без циклов (рисунок 3.1). Такая топологическая структура также не характерна для лигнина. Её анализ в сопоставлении с цепной структурой позволил наглядно проследить, как отражается на кривых ММР продуктов деструкции присутствие разветвляющих узлов. Наиболее характерной особенностью ММР продуктов разрушения макромолекулы с дендровидной структурой оказалась бимо-дальность распределения, особенно заметная при умеренной степени деструкции. Очевидная причина этого в том, что крупные фрагменты макромолекулы могут образоваться только при удалении ребер графа, расположенных не очень далеко от корня. Число таких ребер относительно невелико, следовательно, и вероятность их случайного удаления мала. Напротив, число ребер на «периферии» графа значительно больше. Соответственно значительно выше и вероятность отрыва низкомолекулярных фрагментов в начальном периоде деструкции. Постепенно, по мере деградации крупных фрагментов, кривая ММР становится унимодальной.

Рисунок 3.1- Фрагмент дендровидной топологической структуры макромолекулы с разветвлениями от каждого узла (вариант Б)

Рисунок 3.2 - Фрагмент дендровидной топологической структуры макромолекул с редкими разветвлениями (вариант В, сплошные линии), с редкой внутримолекулярной сшивкой (вариант Г, сплошные и пунктирные линии) и более частой внутримолекулярной сшивкой (вариант Д, сплошные, пунктирные и двойные пунктирные линии)

Дендровидная структура с редкими разветвлениями (рисунок 3.2, вариант В) является промежуточной по отношению к двум рассмотренным

выше моделям, она содержит, наряду с разветвляющими ФПЕ, цепные участки разной протяженности. Кривые ММР в этом случае также занимают промежуточное положение: при небольшой степени деструкции они бимодальны, но «пик», соответствующий высокомолекулярной фракции, быстро исчезает, и уже при степени деструкции (доле разорванных связей) у = 0,15 кривая ММР становится унимодальной с узким распределением по

степени полимеризации.

Дендровидно-сшитая структура наиболее близка к реальной топологической структуре нативного лигнина. Модели для анализа получены добавлением внутримолекулярных сшивок к дендровидным структурам с редкими разветвлениями. Места сшивок разыгрывались и были, таким образом, случайными. Уже небольшое число внутримолекулярных сшивок -4 на 100 ФПЕ- заметно отражается на процессе деструкции (рисунок 3.3, вариант Г). Если в структурах без циклов разрыв каждой связи сопровождается деполимеризацией, то в рассматриваемых структурах к деполимеризации приводит только разрыв связей, не входящих в цикл. В результате этого при одинаковой доле разорванных связей (в сравнении со структурами без циклов) наличие циклов обусловливает уменьшение доли низкомолекулярных фракций и увеличение доли и степени полимеризации высокомолекулярных фракций. Дальнейшее увеличение доли связей, включенных во внутримолекулярную сшивку, до 10 на 100 ФПЕ приводит к углублению указанных различий (рисунок 3.3, вариант Д).

0.02

200

Рисунок 3.3 - Кривые ММР продуктов деструкции полимера (х-меров) с дендровидно-сшитой топологической структурой (слева вариант Г, справа вариант Д): 1 - степень деструкции у= 0,004; 2 - у~ 0,019; 3 - у~ 0,077; 4-у =0,308.

Для описания динамики, разрыва связей между ФПЕ (деструкции) можно принять первый кинетический порядок с константой скорости к': 1-у - ехр(-й'г),

г = 1Ц1 -у).

(3.2)

Здесь г- характеристическое время, т. е. продолжительность достижения заданной степени деструкции.

-1п(1-й)

г, мин

Рисунок 3.4 - Полулогарифмические анаморфозы кинетических кривых (вариант топологической структуры Д, однотипные связи между ФПЕ): 1 - динамика деструкции; 2 - динамика делигнификации.

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

Результаты имитационного моделирования позволяют перейти от оценок степени деструкции лигнина к оценкам степени делигнификации. Предполагается, что в раствор переходят х-меры, степень полимеризации которых меньше зс,р. Массовая доля растворившихся фракций в (то есть степень делигнификации) при фиксированной степени деструкции может быть определена как отношение площади под соответствующей кривой ММР левее абсциссы х^ к общей площади под этой кривой. На рисунке 3.4 (линия 2) показана динамика делигнификации, вычисленная для лигнина с топологической структурой, соответствующей варианту Д. Стандартная математическая обработка показала, что этот процесс близок к первому кинетическому порядку, т. е. может быть описан уравнением

1-0 = ехр(-£"-г), (3-3)

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

При анализе влияния топологической структуры макромолекул лигнина на характер деструкции не обсуждалась природа связей между ФПЕ. Тип связей определяет только скорость деструкции, но не вид кривых ММР. Если макромолекулы содержат связи разной прочности и они распределены случайно, то при одинаковой степени деструкции кривые ММР

также одинаковы, независимо от типа связей между ФПЕ. Скорость же процессов деструкции и делигнификации, естественно, зависит от характера связей.

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

4 9 Г '1

ио - Е = Е К ехр(-£,. г)|, (3.4)

ы /=1

где «о/ ~ число связей г'-го типа в исходной макромолекуле, причем

ч

5>о,- = ло;

;=1

и, - число разорванных связей г-го типа; V - константа скорости разрыва связей /-го типа.

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

Г,в

Рисунок 3.5 - Динамика деструкции и делигнификации (вариант топологической структуры Д, разнотипные связи между ФПЕ): I - динамика деструкции; 2 - динамика делигнификации.

О 100 200 300 400 500

г, мин Результаты имитацион-

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

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

Во-вторых, распределение кинетических параметров при делигни-

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

4 Кинетика натронной варки

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

В самом общем виде натронный процесс может быть представлен схемой реакций

I + ЛаСЮ-> I,; (4.1)

С + ЫаОН Ср, (4.2)

где Ьр и Ср~ переходящие в раствор продукты деструкции лигнина и углеводов. Для определения кинетических параметров была выполнена серия варок при значительном избытке активной щелочи в сравнении с необходимым теоретическим ее расходом, что позволило принять в первом приближении нулевой частный порядок по ИаОН. В этом случае система уравнений скоростей реакций (4.1) и (4.2) записывается в виде двух независимых дифференциальных уравнений:

-кс[С].

dr LL (4.3)

d[C]_

dr

В случае полихронного процесса, очевидно,

кшх

L(r)= J L(kL,T)f(kL)dk!;

о (4.4)

кш\

С(г)= j C(kc,z)f(kc)dkc, о

где L(ku г) и С(кс, т) - решения (4.3); Цг) и С(г) - экспериментально наблюдаемые кинетические кривые; fikC) и _f[kc) - функции плотности распределения констант скоростей для лигнина и углеводов. Для решения интегральных уравнений Фредгольма первого рода (4.4) относительно функции распределения констант скоростей был использован стандартный алгоритм метода регуляризации Тихонова. Результат одного из вычислений

представлен на рисунке 4.1. Остальные распределения выглядят аналогично.

/№> 165, 160 155 150 145 140 1Л5 по

6 Kxio3

Рисунок 4.1- Плотность функции распределения лигнина сосны при 393К: f-плотность функции распределения; к- значение константы скорости

Математическое моделирование процесса с переменной концентрацией щелочи выполнялось на основе предложенной выше теории. Для определения числа ансамблей использовался дифференциальный метод. Оказалось, что компоненты древесины при натронной варке ведут себя так, как если бы они состояли из трех видов лигнина Ц, ¿2 и и трех видов углеводов Сь С-1 и Сз, различающихся реакционной способностью. Этот результат хорошо согласуется с известным фактом существования трех основных видов связей между фенилпропановыми звеньями макромолекулы лигнина, определяющих сложность ее структуры. Поэтому можно предположить, что в силу нерегулярности строения полимера, найденные методом Тихонова распределения являются суммой трех распределений, каждое из которых определяется дисперсией кинетических свойств одного из трех типов связей в макромолекуле лигнина. Исходя из этого минимизируемый функционал имеет вид:

(4.5)

где у = 3 , г-число точек, /¡'- значение плотности распределения, найденное по методу Тихонова, /°{тк,сгк)- уравнение для плотности распределения. Модель кинетики делигнификации, сформулированная в виде системы интегральных уравнений, позволяет с достаточной точностью описывать кинетические кривые процесса делигнификации. Однако моделям этого типа присущ очевидный недостаток, связанный с тем, что явное введение температуры в такую модель затруднено. Поэтому возникает необходимость переформулировать задачу так, чтобы получить модель, позволяющую учитывать влияние температуры явным образом. Таким свой-

ством обладают модели формальной кинетики, записанные в виде систем обыкновенных дифференциальных уравнений.

Для процесса делигнификации система уравнений, аппроксимирующая натронную варку сосны и лиственницы, запишется в виде:

^-кШон]

ах

^ = ОН]

аЖ-кс\С\он] (45)

ат

[ОЯ]= «4 + «, {[С„, ]+ [С„]+ [сп}+ а2 {[./,„]+ [1п}+ [/-„]}

где а,), а,,а2- экспериментально определяемые параметры.

Для определения констант скоростей кц и ка систему уравнений (4.5) решаличисленно. Варьируя значения констант, минимизировали сумму квадратов отклонений расчетных значений [I] и [С] от их экспериментальных значений.

Численные значения найденных величин приведены в таблицах 4.1 и 4.2 (для древесины сосны). Об удовлетворительной точности аппроксимации кинетических кривых можно судить по рисунку 4.2. На это же указывают результаты статистического анализа — малые остаточные дисперсии Лег и высокие значения критерия эффективности Аналогичные результаты получены для древесины лиственницы.

Таблица 4.1 - Константы скорости, остаточные дисперсии и критерии эффективности: натронная варка, лигнин сосны

Т,К Константы скорости, дм г с 2 •У ост **

Ша

393 -12,5 -15,9 -18,8 1,91 148

413 -11,1 -13,4 -17,9 2,14 98

433 -9,88 -П,4 -13,5 1,44 188

453 -8,85 -9,94 -11,4 1,85 120

473 -7,40 -8,48 -9,95 0,96 81

Таблица 4.2 - Константы скорости, остаточные дисперсии и критерии эффективности: натронная варка, углеводы сосны

т,к Константы скорости, дм"г 'с 1 с2 л ост ^эф

1п£С1 1п кс2 \пксз

393 -8,01 -12,0 -19,3 4,14 60

413 -11,6 -13,9 -16,2 1,44 122

433 -10,4 -12,0 -13,2 1,68 144

453 -9,27 -10,6 -11,9 0,90 102

473 -7,86 -8,73 -9,09 0,73 170

Рисунок 4.2 - Динамика изменения концентрации лигнина (а) и углеводов (б) при натронной варке сосны: [£], [С] - концентрации лигнина и углеводов, г/дм\ т- время, с. Температура, К: 1 - 393; 2 - 433; 3 - 473; точки - результаты экспериментов, линии - аппроксимация

5 Кинетика сульфатной варки

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

Наиболее лабильны по отношению к действию щелочи эфирные связи Са~О С4 в фенольных арилпропановых звеньях макромолекул лиг-

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

Структуры с Ср —О—(24 связями между фенилпропановыми звеньями и с фенольными гидроксилами превращаются в условиях щелочной делиг-нификации в относительно устойчивые структуры типа хинонметидов.

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

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

В соответствии со сформулированными предположениями формально-кинетическая схема реакций, протекающих при сульфатной варке, представлена в виде:

Ьг+тон Ьы[ (5.1)

и

Ч (5.2)

Ьы + Ьзш \

к-2

кз

¿М* (5.3)

Ьш + ЫаОЕ (5.4)

Ср—> С0!с+££Г, (5.5)

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

Рисунок 5.1 - Динамика изменения концентрации компонентов сульфатной варки сосны (а) и лиственницы (б) при температуре 433 К: 1 - углеводы; 2 - лигнин; 5 - активная щелочь; точки - результаты экспериментов, линии -аппроксимация

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

Кинетика сульфатной варки древесины сосны и лиственницы изучена в диапазоне температур от 393 до 473 Л". Расчет кинетических параметров выполнялся по тому же алгоритму, что и при анализе результатов натронной варки. О качестве аппроксимации можно судить по рисунку 5.1. Адекватность математической модели реальному процессу и ее высокая эффективность подтверждена статистическим анализом.

Сравнение кинетических параметров натронной и сульфатной варки выявило ряд особенностей. При обоих способах делипшфикации поведение компонентов древесины во многом идентично. Сохранилось три группы кинетических ансамблей лигнина и и три группы ансамблей углеводов С\, Сг и С3. Массовые доли всех групп остались практически неизменными. Почти не изменились также величины энергии активации процессов для ¿3, С], С2 и С3. Наиболее интересным результатом оказалось снижение величины энергии активации только группы кинетических ансамблей ¿2 примерно на 20 кДж/моль в сравнении с натронной варкой. Это наблюдение косвенно подтверждает справедливость предположений, сформулированных в начале раздела.

Концентрации,

г*НГ3,г

6 Кинетика натрошю-антрахиноиной варки

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

В диссертации предложен следующий кинетический механизм де-лигнификации:

кц

1/+ДЬОН -> 1р}\

к^

и+КЬОН < Ь^ : к-1 к4

■"ХМ + КаОН -1р2;

Ьум + ЛНО 1а:

к-5

и -> 1р2 + АО \

С + АО —> Ст + АНО;

Ср + АО —Сж + АНО :

км

ЛНО + Ц -> М;

ки

и + КаОН -> Ц3 ,

где Ьа - аддукт лигнина и антрагидрохинона; М — продукты реакции «отравления» катализатора.

Анализ системы дифференциальных уравнений скоростей стадий позволил вывести два следствия. При малых начальных концентрациях [57Г] и антрахинона лимитирующими стадиями являются реакции образования аддуктов и скорость делигнификации должна быть линейной функцией этих концентраций. При высоких концентрациях и [ЛО] лимитирующей стадией становится образование хинонметидной структуры, протекающее без участия этих соединений, поэтому сульфатная и натрон-но-антрахинонная варки должны протекать в этом случае с одинаковыми скоростями. Выполнение обоих следствий подтверждено экспериментами.

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

с1т

¿г

4Ср2] йх

V ;=1 ;=1 7

Н<?]«ю3. ли3 2»0

(6.1)

140

ив Рисунок 6.1 - Динамика измена нения концентраций компота нентов субстрата натронно-антрахинонной варки древесины сосны: 1 — углеводы; 2 — !0 лигнин; 3 - активная щелочь; 40 4 _ антрахинон о

2го Математическую обра-

240 ботку результатов эксперименте тального исследования наш тронно-антрахинонной варки древесины сосны и лиственни-1~° цы проводили по изложенной ранее методике. Статистиче-40 ский анализ подтвердил адек-о ватность математического описания результатам экспериментов. О точности аппроксимации можно судить по рисунку 6.1, на котором приведена часть результатов.

7 Полиэкспоненциальная кинетика делигнификации и принцип линейности свободных энергий

Деструкция лигнина при щелочной делигнификации происходит путем разрыва главным образом эфирных связей между фенилпропановыми структурными единицами. Если считать основной причиной полихронно-сти кинетики при щелочной делигнификации дисперсию энергии связей между ФПЕ из-за наличия разнотипных связей и влияния их окружения в макромолекуле лигнина, то следует ожидать, что для групп кинетических ансамблей лигнина в реакциях щелочного гидролиза связей будет выполняться принцип линейности свободных энергий. Анализ следствий, вытекающих из этой гипотезы, позволил получить экспериментальное подтверждение обсуждаемого явления. Одним из таких следствий является обнаруженное в ходе экспериментов существование кинетического компенсационного эффекта. В пользу гипотезы свидетельствует также установленная экспериментально линейность зависимостей, вытекающих из логарифмической формы уравнения Эйринга (рисунок 7.1) и изокинетиче-ских соотношений (рисунок 7.2): ЬН*=А+Т-

где АН* и АЗ* - соответственно энтальпия и энтропия активации; Г - температура; А - эмпирический коэффициент.

-1а {¡с/Т\

26 24

•» Л

20 П 16 и

<

150

250

А*?, ДжлюлъК

300 350

V*—1

2,0

2,4 2,6

-150

■100

-50

Г'

■ 103,1С1

О 50

ДЧх, Дж/люлъ-К

Рисунок 7.1 (слева) - Диаграммы Эйринга для лигнина (сплошные линии -сосна, пунктирные - лиственница; 1 - 2 - ¿2; 3 - ¿3).

Рисунок 7.2 (справа) - Изокинетическая диаграмма для лигнина сосны (светлые точки) и лиственницы.

8 Кинетика окислительной делигнификацин

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

В диссертации предложена и обоснована следующая схема образования надуксусной кислоты:

к,

Н202 + Cat ~—Г CatO+Н20; k-i

АсОН + CatO 1—" АсООН + Cat. к-2

к?

CatO -► Cat + 02;

h

1Ш2+АсОН 1—* АсООН+Н20: k-t

ks

Н202 -- Н2С) + 02

Здесь знак CatO использован для обозначения пероксокомплекса.

В соответствии с этой схемой кинетика образования надуксусной кислоты представлена следующей системой дифференциальных уравнений скоростей:

= [Я2 02 }[Cat] + к_х [<CatO][H2 О) -

ат

~ к4[Н202][Ас0Щ + к_4[АсООН][Н2О] - к5[Н202];

d[ АсООН] = к^Ас0щ^Сащ _ к_2 [AcOOH][Cat] + dr

+ к4 [Н202 ][АсОН] - к_4 [АсООН][Н2О];

= к, [Н2 02 }[Cat] - к_} [СаЮ)[Н2 О] -

dr

- к2 [CatO][AcOH] + к_2[Cat][AcOOH] - к3[СаЮ].

50 100 150 200 Продолжит етьность, мин

Рисунок 8.1-Динамика изменения текущих концентраций пероксида водорода

(а) и надуксусной кислоты

(б)

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

[АсООН], «о/ц-дм* 0,6

50 100 150 :оо Продолжтигшосп., ш.н

250

д й- — .....

it. ы

к

• ЗЬЗ А ■ 353 А

* 342 К о 333 К *321 А *313 К

г 11Г*, С

Рисунок 8.2 - Кинетические кривые делигнификации сосновой древесины; точки результаты эксперимента; линии - расчетные значения

Больший интерес для процесса делигнификации представляет динамика коцентрации надуксусной кислоты в растворе. Как видно из рисунка 8.1 (б), кривые [АсООН] ~ г экстремальны. Это обусловлено обратимостью реакций образования надуксусной кислоты и одновременным протеканием процессов разложения пероксида водорода и пероксокомплексов с выделением молекулярного кислорода, которые необратимы.

Эксперименты по изучению кинетики делигнификации включали серию изотермических варок сосновой и осиновой древесины при различных температурах в диапазоне от 313 до 363 К. Часть результатов показана точками на рисунке 8.2.

Отметим наиболее интересные особенности кинетических кривых растворения лигнина.

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

Во-вторых, по мере повышения температуры варки снижается максимально достижимая степень конверсии лигнина. Это особенно заметно при температурах выше 340 К. Причиной этого является, очевидно, инактивация (кислотная конденсация) лигнина.

В-третьих, при варке осиновой древесины явление инактивации лигнина не обнаружено. Аналогичная особенность осинового лигнина проявляется и при других способах делигнификации.

Модель кинетики делигнификации представлена приведенной выше системой дифференциальных уравнений скоростей реакций компонентов варочного раствора, дополненной уравнениями скоростей делигнификации и инактивации лигнина:

щр]

(¿Т

= к6[Ц[АсООН]"-,

ат

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

9 Моделирование совместной щелочной варки древесины разных пород

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

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

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

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

(9-1)

= -\Срт-\}{ОН-\, (9-2)

[ОН-]=[ОН-]0 -5:{в1(Л1[1^)]}-хи(40Е[ср&-)]1

1 I ' ) 1 I / ) (9.3)

где г = 1, 2 или 3; индекс) обозначает принадлежность компонента сосне или лиственнице; / - время.

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

9.1 Прочностные свойства целлюлозы

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

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

лголозе (в диапазоне от 4 до 7 %) могут быть аппроксимированы линейными уравнениями регрессии общего вида

Г^йш+МГ, (9.4)

где У, - свойства целлюлозы; Х- массовая доля лигнина в целлюлозе, %; Ь0 и Ь\ - коэффициенты регрессии.

С целью прогнозирования свойств целлюлозы при варке смешанного древесного сырья нужно:

- используя математическую модель вида (9.1) - (9.3), вычислить выходы целлюлозы из древесины в любой момент варки отдельно для каждой породы, как отношение

[С]+Щ [СЬ+Мо'

а также массовые доли лигнина в целлюлозе — как отношение

щ .

[С] + [Ь]'

- вычислить прочностные свойства целлюлозы, получившейся из древесины каждой породы, по уравнениям вида (9.4);

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

9.2 Критерий управления варочным процессом

Интегральные характеристики являются информационными критериями-комплексами, дающими количественную характеристику некоторых свойств и особенностей технологических процессов. Основа интегрального параметра — интеграл вида:

ь

а

где х(0 - исследуемая функция; - весовая функция, вид которой задается в соответствии с классом явления исследуемого объекта. Интегральные параметры представляют собой комплексную оценку переменных. Они обладают свойством аддитивности, т. е. при разбиении области параметров, характеризующих данный температурно-временной режим, на ряд подобластей (С,(С2,...,Ст), общий интегральный параметр будет равен сумме частных интегральных характеристик: )= ЛС,) + Л С3)+... + ПСт).

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

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

Я = ]к0е ^Л.

(9.5)

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

Фактически, с учетом того, что Е

к0еКт^=к, (9.6)

Я-фактор есть интеграл по времени от константы скорости, которая зависит от температуры и может изменяться в ходе варки.

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

ОД. Л

откуда

- = (9.7)

[¿] = [^]ехр

V 'о

(9.8)

или, что то же самое, .

[¿] = Г1оКн. (9.9)

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

Яа0;О + Я(^2) = Я(/0;Ц (910)

Я-фактор не аддитивен по породному составу древесного сырья, т.е. при совместной варке древесины разных пород А и В имеет место неравенство

НА + НВФНАВ. (9.11)

Это обусловлено тем, что для двух различных пород А а В уравнения скоростей (9.7) должны быть записаны в виде

Л

О^—ки 1 (9Л2)

Поскольку суммарный липшн

[¿] = [£Л + [1в], (9-13)

общая скорость процесса должна также описываться уравнением

Ш

-)т* (9Л4)

Это означает, что для аддитивности //-фактора должно выполняться условие

е~н =е~Нл -е'"в. (9.15)

В то же время из (9.13) следует

[Ьй}е-Н ^Ь^е-" +[Ьва}е-н. (9.16)

Соответственно, требование аддитивности /У-фактора по отношению к породному составу не выполняется даже в случае простейшей кинетики первого порядка. Тем более оно не выполняется для процессов с более сложным кинетическим механизмом и процессов с полихронной кинетикой.

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

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

Выводы

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

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

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

4. Решена проблема построения последовательной математической теории кинетики процессов деструкции основных высокомолекулярных компонентов древесины и ее делигнификации на основе принципа полихронности, теории Крамерса и теории случайных процессов.

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

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

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

8. Установлено, что при всех вариантах щелочных варок поведение кинетических ансамблей лигнина и углеводов древесины подчиняется принципу линейности свободных энергий.

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

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

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

По теме диссертации опубликовано 42 работы. Основное содержание

диссертации отражено в следующих публикациях. Монографии

1. Пен, В.Р. Кинетика делигнификации древесины [Текст] / В.Р. Пен, Р.З. Пен - Красноярск: СибГТУ, 1998. - 200 с.

2. Пен, В.Р. Математический аппарат исследования структуры и деструкции полимеров [Текст] / В.Р. Пен. - Красноярск: СибГТУ, 1999. - 137 с.

3. Пен, В.Р. Теоретические основы делигнификации [Текст] / В.Р. Пен, Р.З. Пен. - Красноярск: ООО Красноярский писатель, 2007. - 348 с.

Прочие публикации

4. Пен, В.Р. Полихронная кинетика деструкции полисахаридов в процессе щелочной делигнификации [Текст] / В.Р. Пен // Тез. докл. 6-й Всесоюз. конф. по физике и химии целлюлозы. - Минск, 1990. - С. 180.

5. Пен, В.Р. Кинетика натронно-антрахинонной делигнификации древесины [Текст] / В.Р. Пен, Р.З. Пен, И.Л. Шапиро // Журнал прикладной химии. - 1997. -Т. 70. - Вып. 2. - С. 330-333.

6. Реп, V.R. Isokinetic correlation in alkaline delignification of larch wood [Текст] / V.R. Pen, M.Yu. Katrukhina // Larix-98: World Resources for Breeding. Resistance and Utilisation. IUFRO Interdivis. - Krasnoyarsk, 1998. - P. 73.

7. Pen, V.R. Kinetics of the soda and soda-anthraquinone pulping of larch wood in a wide temperature range [Текст] / V.R. Pen, I.L. Shapiro // Larix-98:

World Resources for Breeding, Resistance and Utilisation. IUFRO Interdivis. - Krasnoyarsk, 1998. - P. 73.

8. Особенности кинетики окислительной делигнификации лиственной древесины [Текст] / В.Р. Пен [и др.] Н Тез. докл. конф. «Проблемы химико-лесного комплекса». - Красноярск, 1998. - С. 79.

9. Ступенчатая окислительно-щелочная делигнификация древесины без применения соединений серы и хлора [Текст] / В.Р. Пен [и др.] // Тез. докл. конф. «Проблемы химико-лесного комплекса». - Красноярск, 1998. - С. 124.

10. Пен, В.Р. Стохастическое моделирование кинетики деструкции полимеров [Текст] / В.Р. Пен, С.И. Левченко // Тез. докл. конф. «Проблемы химико-лесного комплекса». - Красноярск, 1998. - С. 176.

11. Пен, В.Р. Компьютерное имитационное моделирование деструкции полимеров методом Монте-Карло [Текст] / В.Р. Пен, С.И. Левченко. // Тез. докл. конф. «Проблемы химико-лесного комплекса». - Красноярск, 1998. -С.248.

12. Пен, В.Р. Кинетика делигнификации древесины. 9. Стохастическая модель деструкции лигнина как полидисперсного полимера [Текст] /

B.Р. Пен // Химия растительного сырья. - 1998. - № 3. - С. 119-122.

13. Пен, В.Р. Имитационное моделирование деструкции макромолекул лигнина с разной топологической структурой [Текст] / В.Р. Пен, Р.З. Пен Н Журнал прикладной химии. - 1999. - Т. 72. - Вып. 7. - С. 11941198.

14. Оптимальный режим получения целлюлозы без применения серы [Текст] / В.Р. Пен [и др.] // Тез. докл. Всероссийской конф. с международным участием «Достижения науки и техники - развитию сибирских регионов». - Красноярск, 1999. - Часть 2. - С. 271-272.

15. Пен, В.Р. Пероксобораты как катализаторы при делигнификации древесины пероксидом водорода [Текст] / В.Р. Пен, Е.Ю. Беляев, И.В. Ганзен // Тез. докл. Всероссийской конф. с международным участием «Достижения науки и техники - развитию сибирских регионов». - Красноярск, 1999. - Часть 2. - С. 272-273.

16. Кинетика делигнификации лиственной древесины перуксусной кислотой [Текст] / В.Р. Пен и др. // Вестник Сибирского гос. технол. унив-та. - 1999. -№ 1.-С. 76-83.

17. Пен, В.Р. Стохастическая модель делигнификации [Текст] / В.Р. Пен // Вестник Сибирского гос. технол. унив-та. -1999. - № 1. - С. 84-97.

18. Пен, В.Р. Кинетика делигнификации древесины. 10. Некоторые особенности кинетики сульфатной и натронно-антрахинонной варки [Текст] 1 В.Р. Пен [и др.] // Химия растительного сырья. - 2001. - № 1. -

C. 37-41.

19. Пен, В.Р. Оценка плотности функции распределения для полиэкспоненциальной кинетики деструкции полимеров древесины [Текст] / В.Р. Пен // Химия растительного сырья. - 2004. - № 2. - С. 97-100.

20. Пен, В.Р. Об одном способе получения математических моделей кинетики деструкции полимеров древесины [Текст] / В.Р. Пен // Химия растительного сырья. - 2004. - № 2. - С. 101-103.

21. Пен, В.Р Способ получения математических моделей кинетики деструкции полимеров на основе принципа полихронности. [Текст] / В.Р. Пен, С.И. Левченко // Известия вузов. Химия и химическая технология.

- 2005. - Т.48. - Вып. 8. - С. 99-105.

22. Пен, В.Р. Влияние среды на кинетику разложения хилоловых эфиров [Текст] / В.Р. Пен, С.И. Левченко // Известия вузов. Химия и химическая технология. - 2006. - Т.49. - Вып. 7. - С. 73-76.

23. Пен, В.Р. Использование балансовых уравнений для учета стехиометрии в моделировании кинетики реакций полимеров нерегулярного строения [Текст] / В.Р. Пен // Вестник КрасГАУ. - 2006. - С. 147-150.

24. Пен, В.Р. О критерии управления процессом варки древесины [Текст] / В.Р. Пен, О.Г. Морозова // Вестник КрасГАУ. - 2006. - С. 150155.

25. Пен, В.Р. Управление процессом делигнификации смешанного древесного сырья [Текст] / В.Р. Пен, И.В. Ковалев // Вестник СибГАУ. -2007.-№1,-С. 57-61.

26. Пен, В.Р. Кинетика деструкции полимеров с распределенными кинетическими параметрами как стохастический процесс [Текст] / В.Р. Пен, И.В. Ковалев, С.И. Левченко // Вестник СибГАУ. - 2008. - №4. - С. 2830.

27. Пен, В.Р. Низкотемпературная окислительная делигнификация древесины. 15. Варка с рециклом щелока [Текст] / В.Р. Пен, Н.В. Каретникова, И.Л. Шапиро, A.B. Бывшев // Химия растительного сырья. - 2008.

- № 2. -С. 11-14.

28. Пен, В.Р. Учет стехиометрии в кинетике делигнификации древесины [Текст] / В.Р. Пен, С.И. Шишканова, И.Л. Шапиро // Химия растительного сырья. - 2008. - № 2. - С. 137-138.

29. Пен, В.Р. Критерий управления щелочной варкой смесей древесных пород [Текст] i В.Р. Пен, С.И. Левченко, С.Л. Карпенко // Химия растительного сырья. - 2008. - № 2. - С. 139-140.

30. Пен, В.Р. Распределение Больцмана и полиэкспоненциальная кинетика натронной варки древесины [Текст] / В.Р. Пен, С.И. Левченко, И.Л. Шапиро 1! Химия растительного сырья - 2008. - № 2. - С. 141-142.

31. Пен, В.Р. О критерии управления процессом щелочной варки древесины [Текст] / В.Р. Пен, С.И. Шишканова, С.Л. Карпенко // Химия растительного сырья. - 2008. 2. - С. 143-145.

32. Пен, В.Р. Полидисперсность целлюлозы из сосны и лиственницы

[Текст] / В.Р. Пен, И.Л. Шапиро // Международный журнал прикладных и фундаментальных исследований. -2009. - № 4. - С. 81-82.

33. Пен, В.Р. Теплотворная способность пероксидного щелока [Текст] / В.Р. Пен, Н.В. Каретникова // Успехи современного естествознания. -2009,- №8.-С. 104-105.

34. Пен, В.Р. Кинетика разложения пероксида водорода в присутствии вольфрамата, молибдата, серной кислоты и их смесей [Текст] / В.Р. Пен, Н.В. Каретникова, И.Л. Шапиро, И.В. Мирошниченко // Успехи современного естествознания. - 2010. - № 9. - С. 212-213.

35. Пен, В.Р. Синергизм в каталитических системах делигнификации древесины перуксусной кислотой [Текст] / В.Р. Пен [и др.] // Успехи современного естествознания. - 2010. - № 9. - С. 213-214.

Подписано в печать 07.09.2010. Формат 60><84 1/16. Уч-изд. 2.5 . Тираж 110 экз. Изд.№ 2/6(2010)

Заказ № 862._

Редакционно-издательский центр СибГТУ 660049, Красноярск, пр. Мира, 82.