close

Вход

Забыли?

вход по аккаунту

?

Теоретико-графовый метод анализа информативности кинетических экспериментов при определении параметров.

код для вставкиСкачать
МАТЕМАТИКА и ФИЗИКА
1126
УДК 544.4
ТЕОРЕТИКО-ГРАФОВЫЙ МЕТОД АНАЛИЗА ИНФОРМАТИВНОСТИ
КИНЕТИЧЕСКИХ ЭКСПЕРИМЕНТОВ ПРИ ОПРЕДЕЛЕНИИ ПАРАМЕТРОВ
© С. И. Спивак1*, А. С. Исмагилова2, Р. А. Гибаева2
1Башкирский государственный университет
Россия, Республика Башкортостан, 450076 г. Уфа, ул. Заки Валиди, 32.
2Башкирский
государственный университет, Нефтекамский филиал
Россия, Республика Башкортостан, 452681 г. Нефтекамск, ул. Социалистическая, 65.
Тел.: +7 (347) 229 96 35.
*Email: semen.spivak@mail.ru
Предметом исследования настоящей работы являются обратные задачи идентификации механизмов сложных химических реакций. Нахождение базиса нелинейных параметрических функций кинетических констант сводится к анализу графа механизма сложной химической реакции. Рассмотрен механизм синтеза винилхлорида на катализаторе «сулема–
активированный уголь».
Ключевые слова: Обратная задача, кинетические параметры модели, информативность, граф механизма реакции.
Введение
Обратная задача состоит в определении констант скоростей элементарных стадий, входящих в механизм
сложной химической реакции, на основе экспериментальных данных о концентрациях участвующих в реакции
веществ [1]. Недостаточная информативность реально доступных массивов кинетических измерений приводит
к неоднозначности решения обратной задачи. Важно знать, какие из индивидуальных констант однозначно
определяются по заданной структуре эксперимента, какие определяются только в виде функциональных комбинаций, сколько независимых комбинаций, каков их явный вид. Если число таких комбинаций меньше общего
числа констант, то в технологических целях удобно иметь дело с моделями, содержащими меньшее число параметров.
В работах [2–4] проведено детальное математическое исследование проблемы информативности. Построена методология анализа информативности кинетических измерений при решении обратных задач, позволяющая выделять число и вид нелинейных параметрических функций (НПФ) констант скоростей реакций, допускающих однозначную оценку по кинетическому эксперименту разных типов. В качестве математического аппарата используется теория неявных функций и теория непрерывных групповых преобразований. Трудность
состоит в том, что использование этой теории связано со сложными аналитическими вычислениями, применение которых сильно затрудняется из-за большой размерности в реальных задачах построения кинетических
моделей.
В [5] предложен метод, позволяющий упростить исследование на информативность кинетических моделей
сложных реакций. Этот метод основан на теоретико-графовом анализе независимых маршрутов графа Вольперта [6]. Анализ информативности сложной схемы реакций разделяется на анализ информативности для тех схем,
которые отвечают за протекание реакции по каждому из независимых маршрутов.
В настоящей работе решается задача нахождения базиса НПФ кинетических констант при помощи графа
механизма реакции, введенного М. И.Темкиным [7].
Теоретико-графовый метод нахождения НПФ механизма синтеза винилхлорида
на катализаторе «сулема-активированный уголь»
В графе Темкина вершины представляют промежуточные вещества, ребра – стадии [7]. Напомним, что
вещества, концентрации которых не могут быть определены в ходе протекания реакции, называют промежуточными. Интерпретация механизма реакции на графе Темкина удобна в случаях, когда реакции линейны относительно промежуточных веществ.
Для механизма синтеза винилхлорида на катализаторе «сулема-активированный уголь» [1]
1 ) HgCl
2
 HCl  С 2 H
k1
k1
2
 HgCl
k 10
k2
2 ) HgCl
3 ) HgCl
2
2
 HCl  СCl  HgCl
k 20
 HCl  С 2 H
2
2
C2H
2
1 ) Y1  X 1  Y 2
 HCl
k 10
2 ) Y1  X
 2 HCl
k3
2
 HCl  C 2 H 3 Cl  HgCl
2
 HCl
k2
2
 Y3
k 20
k3
3 ) Y 2  X 2  X 3  Y1
Вестник Башкирского университета. 2014. Т. 19. №4
ISSN 1998-4812
k4
4 ) HgCl
2
 2 HCl  С 2 H 2  C 2 H 3 Cl  HgCl
2
1127
k4
4 ) Y 3  X 1  X 3  Y1
 HCl
граф Темкина изображен на рис. 1.
Рис. 1. Граф Темкина.
Приведем некоторые понятия из теории графов [8]. Совокупность ребер, продолжающих друг друга, называют цепью. Дерево – цепь, проходящая через все вершины графа к фиксированной вершине (корню).
Построим всевозможные деревья графа исследуемого механизма (рис. 2).
Рис. 2. Деревья графа механизма синтеза винилхлорида.
В графе механизма реакции вес ребра равен скорости расходования (образования) промежуточного вещества в той или иной стадии. Выражения для весов ребер получаются, если скорости прямой и обратной реакций
разделить на концентрации участвующих в реакции промежуточных веществ. Веса ребер, ведущих в фиксированную вершину дерева, для параллельных реакций складываются, а последовательные перемножаются. Величина дерева равна произведению весов всех составляющих его ребер.
Введем обозначения. Пусть [ x 1 , x 2 , x 3 ]  [ C 2 H 2 , C C l , C 2 H 3 C l ] – концентрации измеряемых веществ,
[ y 1 , y 2 , y 3 ]  [ HgCl
2
 HCl , HgCl
2
C2H
 HCl , HgCl
2
2
 2 HCl ] – концентрации промежуточных веществ.


Промежуточные вещества связаны между собой уравнением материального баланса. Через w i , w i обозначим
скорость i-й реакции в прямом и обратном направлении соответственно, 1  i  4 , v i – вес дуги, U – вес графа,
U 
U
j
,
где
U
j
–
вес
дерева,
j-го
т.е.
дерева
с
корнем
в
вершине
j
j
( j  H g C l 2  H C l , H g C l 2  C 2 H 2  H C l , H g C l 2  2 H C l ). Запишем скорости реакций в соответствии с зако-
ном действующих масс:

w 1  k 1 x1 y 1 ,

w1  k 1 0 y 2 ,


w2  k 2 x2 y 1 ,
w 2  k 20 y 3 ,
w3  k3 x2 y 2 ,



w 4  k 4 x1 y 3 , где k i – константа скорости реакции. Веса ребер графа: v 1  k 1 x 1 , v 1  k 1 0 , v 2  k 2 x 2 ,

v 2  k 20 ,

v3  k 3 x2 ,

v 4  k 4 x 1 . Веса деревьев и графа: U 1  v 1 v 2 , U

2


 v1 v 2 , U
3


 v1 v 4 , U
4

 v2 v3 ,
МАТЕМАТИКА и ФИЗИКА
1128








U  v 1 v 2  v 1 v 2  v 1 v 4  v 2 v 3 . Используя правило Мэзона для нахождения концентраций промежуточных
веществ, можно выписать стационарное кинетическое уравнение. Скорость протекания реакции:
2
R 
R
k

2



2
v1 v1 v 2 v 2 v 3 v 4








v1 v 2  v1 v 2  v1 v 4  v 2 v 3
2
'3
' 3
k 1 k 2 k 3 k 4 k 1 0 k 2 0 x1 x 2

k 1 k 2 0 x1  k 2 k 1 0 x 2  k 1 k 4 x1  k 2 k 3 x 2
'
'
'2
' 2

( x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 0 k 2 x 2 k 3 k 4 )( k 2 x 2 x 1  k 1 0 k 2 0  k 4 x 1  k 3 x 2 k 2 0 )
2
2
2
2

2
2
2
( k 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1  k 1 0 k 2 0  k 4 k 1  k 3 x 2 k 2 0 )
2
2
2

2
2

( k 1 x 1 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )
2
2
2
2


2
2
2

( k 1 x 1 k 1 0 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )
2
2
2
2

2
2
2
( k 1 x 1 k 1 0 k 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 k 2 0 )
2
2
2

2
2
2
( k 1 x 1 k 1 0 k 2 x 2 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0  k 4 x 1 k 1  k 3 x 2 )
2
2
2
2


2

2
2

( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  x 2 k 2 0 )
2
2
2
2

2
2
( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  k 4 x 1 k 1  k 3 x 2 k 2 0 )  ( k 1 x 1 k 1 0 k 2 x 2 k 2 0 k 3 k 4 )( k 2 x 2 k 1 x 1  k 1 0 k 2 0  x 1 k 1  k 3 x 2 k 2 0 )
2
2
2

2
2
2
где   k 2 x 2 (1   2 ) k 1 x1 (1   1 )  k 10 k 20  k 4 x1 (1   1 ) k 2 x 2 (1   2 )  k 3 x 2 (1   2 ) k 20 .
Существует ненулевая матрица
1


 1 1

0



0


k4
Q  
 (1   1 ) k 1

0

0

1



k1

0


зависящая от k и  , такая, что
R
k

2



1


1 2

k3

(1   2 ) k 2 


0


0

0


0


1


k2

0
,
Q  0 .
Система, отвечающая последнему соотношению, будет иметь вид:
.
ISSN 1998-4812
Вестник Башкирского университета. 2014. Т. 19. №4
1
p
1   1 k1
1
p
1   2 k 2


k4
p

(1   1 ) k 1  k 4
k3
p
(1   2 ) k 2  k 3
1 p
k1  1

1
1129
 0,
p
k 2  2
 0.
Базис НПФ состоит из двух компонент и может быть представлен в виде:
1 
k4
k1
 k 1 (1   1 ) ,  2 
k3
k2
 k 2 (1   2 )
.
Таким образом, в механизме синтеза винилхлорида на катализаторе «сулема-активированный уголь» доступная экспериментальная информация позволяет определить две независимые параметрические функции.
Описанный в настоящей работе подход позволяет существенно упростить исследование на информативность кинетических моделей сложных реакций.
ЛИТЕРАТУРА
1.
2.
3.
4.
5.
6.
7.
8.
Яблонский Г. С., Быков В. И., Горбань А. Н. Кинетические модели каталитических реакций. Новосибирск: Наука, 1983. 255 с.
Спивак С. И., Горский В. Г. Неединственность решения задачи восстановления кинетических констант // Доклады Академии наук.
1981. Т. 257. №2. С. 412–415.
Горский В. Г., Спивак С. И. Исследование идентифицируемости параметров – один из важнейших этапо построения математических
моделей в химии // Журнал структурной химии. 1988. Т. 29. №6. С. 119.
Кудашев В. Р., Спивак С. И. Информативность кинетических измерений при определении параметров математических моделей нестационарной химической кинетики // Теоретические основы химической технологии. 1992. Т. 26. №6. С. 872–879.
Спивак С. И., Исмагилова А. С. Информативность кинетических измерений и обратные задачи химической кинетики // Доклады
Академии наук. 2013. Т. 451. №3. С. 296–298.
Вольперт А. И., Худяев С. И. Анализ в классах разрывных функций и уравнения математической физики. М: Наука, 1975. 394 с.
Механизм и кинетика сложных каталитических реакций // Лекции, прочитанные на первом симпозиуме Международного конгресса
по катализу. М.: Наука, 1970. С. 57–76.
Оре О. Графы и их применение. М.: КомКнига. 2006. C. 132.
Поступила в редакцию 29.09.2014 г.
МАТЕМАТИКА и ФИЗИКА
1130
GRAPH-THEORETICAL METHOD OF INFORMATIVITY ANALYSIS OF
KINETIC EXPERIMENTS DETERMINING PARAMETERS
© S. I. Spivak1*, A. S. Ismagilova2, R. A. Gibaeva2
1Bashkirsky State University
32 Zaki Validi St., 450076 Ufa, Republic of Bashkortostan, Russia.
2Bashkirsky
State University, Neftekamsk branch
65 Socialisticheskaya St., 452681 Neftekamsk, Republic of Bashkortostan, Russia.
Phone: +7 (347) 229 96 35.
*Email: semen.spivak@mail.ru
This work deals with inverse problems of chemical kinetics. Systems of non-linear ordinary differential, differential algebraic,
or algebraic equations, depending on the conditions of the experiment, namely, non-stationary, stationary, or equilibrium cases, constitute the mathematical object of the study. The inverse problem consists of determining the rate constants of elementary steps based
on experimental data on the concentrations of compounds involved in the reaction. The insufficient informativity of the actually
available arrays of kinetic measurement data results in the non-uniqueness of inverse problem solution. It is important to know which
individual constants are unambiguously determined in an experiment with predefined structure, which ones are only determined as
functional combinations, how many independent combinations are there, and what their explicit form is. It is clear that problems of
this kind are highly laborious in computational respect. The difficulty is that the use of this theory involves complex analytical calculations, the use of which is hampered considerably due to high dimensionality in real tasks of kinetic model building. The main result
involves building a methodology for analyzing the informativity of kinetic measurements in the solution of inverse problems, which
would allow one to determine the number and form of independent combinations of reaction rate constants. The developed algorithms became the basis of software for the solution of inverse problems of chemical kinetics.
Keywords: inverse problem, the kinetic parameters of the model, informativity, graph of the reaction mechanism.
Published in Russian. Do not hesitate to contact us at bulletin_bsu@mail.ru if you need translation of the article.
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
Yablonskii G. S., Bykov V. I., Gorban' A. N. Kineticheskie modeli kataliticheskikh reaktsii [Kinetic Models of Catalytic Reactions].
Novosibirsk: Nauka, 1983.
Spivak S. I., Gorskii V. G. Doklady Akademii nauk. 1981. Vol. 257. No. 2. Pp. 412–415.
Gorskii V. G., Spivak S. I. Zhurnal strukturnoi khimii. 1988. Vol. 29. No. 6. Pp. 119.
Kudashev V. R., Spivak S. I. Teoreticheskie osnovy khimicheskoi tekhnologii. 1992. Vol. 26. No. 6. Pp. 872–879.
Spivak S. I., Ismagilova A. S. Doklady Akademii nauk. 2013. Vol. 451. No. 3. Pp. 296–298.
Vol'pert A. I., Khudyaev S. I. Analiz v klassakh razryvnykh funktsii i uravneniya matematicheskoi fiziki [Analysis in Classes of Discontinuous Functions and Equations of Mathematical Physics]. Moscow: Nauka, 1975.
Mekhanizm i kinetika slozhnykh kataliticheskikh reaktsii Lektsii, prochitannye na pervom simpoziume Mezhdunarodnogo kongressa
po katalizu. Moscow: Nauka, 1970. Pp. 57–76.
Ore O. Grafy i ikh primenenie [Graphs and Their Application]. Moscow: KomKniga. 2006. Pp. 132.
Received 29.09.2014.
Документ
Категория
Без категории
Просмотров
9
Размер файла
394 Кб
Теги
анализа, кинетическая, метод, информативность, эксперимент, определение, графовых, параметры, теоретико
1/--страниц
Пожаловаться на содержимое документа