Кинетика и катализ, 2020, T. 61, № 6, стр. 783-788
Метод решения обратной задачи химической кинетики многостадийных реакций
Н. И. Кольцов *
ФГБОУ ВО Чувашский государственный университет им. И.Н. Ульянова, химико-фармацевтический факультет
428015 Чебоксары, Московский просп., 15, Россия
* E-mail: koltsovni@mail.ru
Поступила в редакцию 16.12.2019
После доработки 17.02.2020
Принята к публикации 21.02.2020
Аннотация
Изложен метод решения обратной задачи химической кинетики по данным нестационарных экспериментов для многостадийных реакций. Метод основан на применении сплайнов, учитывающих релаксационные особенности реакций с монотонной кинетикой. Приведены примеры использования метода для определения интервалов физичных значений констант скоростей стадий двух- и трехстадийной реакций. Оценено влияние структуры сплайна на погрешность определения констант скоростей стадий. Показана эффективность применения метода для решения обратной задачи химической кинетики для многостадийных реакций с учетом ошибки определения экспериментальных данных.
ВВЕДЕНИЕ
Решение обратной задачи химической кинетики, связанной с определением кинетических параметров реакций по экспериментальным данным, является актуальной проблемой. Основными особенностями этой задачи являются погрешность экспериментальных данных и неоднозначность решений. Существующие к ее решению подходы, основанные, как правило, на алгоритмах оптимизации и регуляризации, позволяют определить интервалы изменений значений некоторых констант и комплексов констант скоростей стадий реакций [1–23]. В работе [24] изложен метод решения обратной задачи, позволяющий определять значения отдельных констант скоростей стадий и возможные интервалы их изменений для некоторых классов каталитических реакций без использования этих алгоритмов и учета различной информативности участков релаксационных зависимостей. Представлялось интересным разработать подход к решению обратной задачи для более широких классов многостадийных химических реакций с учетом различной информативности различных участков релаксационных кривых.
В данной статье приведен такой метод, который с помощью сплайнов [25] с разным числом точек на различных участках экспериментальных кинетических кривых учитывает релаксационные особенности протекания реакций в нестационарных условиях. Этот метод позволяет по нестационарным данным с учетом погрешности эксперимента найти интервалы физичных значений констант скоростей элементарных стадий для многостадийных реакций.
ТЕОРЕТИЧЕСКАЯ ЧАСТЬ
Рассмотрим произвольную химическую реакцию, протекающую через стадии вида
где Ak – реагенты, k = 1, …, K – номер реагента, a±ij ≥ 0 − стехиометрические коэффициенты реагента Ak в стадии i = 1, …, s. Динамика такой реакции в открытом безградиентном изотермическом реакторе в рамках закона действующих масс описывается системой K обыкновенных дифференциальных уравнений (ОДУ)(1)
$\begin{gathered} A_{k}^{'} = \Sigma ({{a}_{{ - ik}}} - {{a}_{{ + ik}}})({{r}_{i}} - {{r}_{{ - i}}}){\text{ }} + {{q}^{0}}A_{k}^{0}--q{{A}_{k}}, \\ k = 1, \ldots ,K, \\ \end{gathered} $Если для реакции выполняются линейные стехиометрические законы сохранения (ЛСЗС) вида
где αjk и Cj – константы, то размерность системы (1) можно уменьшить на число этих ЛСЗС. Как показано в работе [26], точное число ЛСЗС равно где Rs − ранг матрицы A = (aik − a−ik), Rk − ранг по комплексам констант скоростей стадий. Выразим c помощью ЛСЗС число зависимых реагентов Ns через остальные (независимые) и исключим их из системы (1)−(2). Если Ns = 0, то система (1)−(2) не меняется. С учетом этого далее будем считать, что уравнения (1)−(2) включают только независимые реагенты и их число равно К.Пусть наблюдаемые концентрации всех реагентов характеризуются монотонным изменением во времени. Составим таблицу их значений Akn от времени tn с шагом h, где n = 1, 2, 3, … − номер измерения. Выберем из этой таблицы набор из N точек (таких наборов может быть несколько, и можно выбрать любой). Каждый такой набор соответствует некоторому кусочно-линейному приближению (линейному сплайну), заменяющему кривые концентраций вписанными непрерывными ломанными линиями, состоящими из N − 1 отрезков прямых (фрагментов). Выберем опорные (используемые в дальнейших расчетах) значения концентраций реагентов для каждого фрагмента как средние арифметические между их границами
(5)
$A_{{kn}}^{*} = {{({{A}_{{kn}}} + {{A}_{{kn}}}_{{ - 1}})} \mathord{\left/ {\vphantom {{({{A}_{{kn}}} + {{A}_{{kn}}}_{{ - 1}})} 2}} \right. \kern-0em} 2},\,\,\,\,n = 1,2{\text{ }}, \ldots ,N,$(6)
$A_{{kn}}^{'} \approx ({{A}_{{kn}}} - {{{{A}_{{kn}}}_{{ - 1}})} \mathord{\left/ {\vphantom {{{{A}_{{kn}}}_{{ - 1}})} {{{H}_{n}}}}} \right. \kern-0em} {{{H}_{n}}}},\,\,\,n = 1,2, \ldots ,N,$(7)
$\begin{gathered} \Sigma ({{a}_{{ - ik}}} - {{a}_{{ik}}})(r_{{in}}^{*} - r_{{ - in}}^{*})--qA_{{kn}}^{*} \approx \\ \approx {{q}^{0}}A_{{kn}}^{0} - ({{A}_{{kn}}} - {{{{A}_{{kn}}}_{{ - 1}})} \mathord{\left/ {\vphantom {{{{A}_{{kn}}}_{{ - 1}})} {{{H}_{n}}}}} \right. \kern-0em} {{{H}_{n}}}}, \\ \end{gathered} $При K × N = 2s система (8) имеет единственное решение
(9)
${{k}_{{ \pm i}}} = {{{{D}_{{ \pm i}}}} \mathord{\left/ {\vphantom {{{{D}_{{ \pm i}}}} D}} \right. \kern-0em} D},\,\,\,\,i = 1, \ldots ,s,$(10)
${{\Delta }_{{ + i}}}\Delta > 0,\,\,\,\,{{\Delta }_{{ - i}}}\Delta \geqslant 0,\,\,\,\,i = 1, \ldots ,s.$При K × N < 2s существует бесконечно много решений, и константы скоростей стадий определяются неоднозначно. Условия физичности аналогичны (10). Иначе, т.е. если условия (8)−(10) не выполняются, система (7) не имеет физичных решений при данном выборе точек для расчетов. В этих случаях необходимо выбрать другой набор экспериментальных точек. Полный перебор всех возможных комбинаций экспериментальных точек позволяет установить интервалы возможных изменений констант скоростей стадий. Если ни один из наборов экспериментальных точек не дает физичных значений констант скоростей стадий, то обратная задача не разрешима и необходимы новые экспериментальные данные или альтернативная стадийная схема реакции.
Пример 1. Пусть реакция А = С + D протекает по двухстадийной схеме
Для нее уравнения (2) запишутся как
(1.1)
$\begin{gathered} A{\kern 1pt} {\kern 1pt} ' = - {{k}_{ + }}_{1}A + {{k}_{{ - 1}}}B + {{q}^{0}}{{A}^{0}} - qA, \\ B{\kern 1pt} ' = {{k}_{ + }}_{1}A - {{k}_{{ - 1}}}B - {{k}_{ + }}_{2}B{\text{ }} + {\text{ }}{{k}_{{ - 2}}}CD + {{q}^{0}}{{B}^{0}} - qB. \\ \end{gathered} $(1.2)
$\begin{gathered} C{\kern 1pt} ' = {{k}_{ + }}_{2}B - {{k}_{{ - 2}}}CD + {{q}^{0}}{{C}^{0}} - qC, \\ D{\kern 1pt} ' = {{k}_{ + }}_{2}B - {{k}_{{ - 2}}}CD + {{q}^{0}}{{D}^{0}} - qD. \\ \end{gathered} $Система (1.1)−(1.2) характеризуется двумя ЛСЗС A + B + C = 1 и A + B + D = 1, т.е. в реакции (II) участвуют два зависимых и два независимых реагента. Выберем в качестве независимых реагентов, например, А и С. Исключим с помощью этих ЛСЗС зависимые реагенты, тогда система (1.1)−(1.2) примет вид
(1.4)
$C{\kern 1pt} ' = {{k}_{ + }}_{2}(1 - A - С) - {{k}_{{ - 2}}}{{C}^{2}} + {{q}^{0}}{{C}^{0}} - qC.$Зададим, например, $k_{{ + 1}}^{*}$ = 1, $k_{{ - 1}}^{*}$ = 1, $k_{{ + 2}}^{*}$ = 1, $k_{{ - 2}}^{*}$ = 1 (“истинные” значения констант скоростей стадий) и н. у. реакции q = q0 = 1, A0 = 1, B0 = 0, C0 = 0, D0 = 0. Численно проинтегрируем систему (1.3)−(1.4) с шагом h на интервале t ∈ [0,tmax]. Получим N точек по каждому реагенту и примем их за исходные данные. Выделим границы фрагментов быстрой, средней и медленной релаксаций. Выберем на двух информативных фрагментах две опорные точки (A1, С1) и (A2, С2), удовлетворяющие вместе с н. у. условиям физичности (10). Подставим эти значения в (1.3)−(1.4) и получим систему линейных по константам уравнений
(1.5)
$A_{1}^{'} = - {{k}_{ + }}_{1}{{A}_{1}} + {{k}_{{ - 1}}}(1 - {{A}_{1}} - {{С}_{1}}) + {{q}^{0}}{{A}^{0}} - q{{A}_{1}},$(1.6)
$C_{1}^{'} = {{k}_{ + }}_{2}(1 - {{A}_{1}} - {{С}_{1}}) - {{k}_{{ - 2}}}C_{1}^{2} + {{q}^{0}}{{C}^{0}} - q{{C}_{1}},$(1.7)
$A_{2}^{'} = - {{k}_{ + }}_{1}{{A}_{2}} + {{k}_{{ - 1}}}(1 - {{A}_{2}} - {{С}_{2}}){\text{ }} + {{q}^{0}}{{A}^{0}} - q{{A}_{2}},$(1.8)
$C_{2}^{'} = {{k}_{ + }}_{2}(1 - {{A}_{2}} - {{С}_{2}}) - {{k}_{{ - 2}}}C_{2}^{2} + {{q}^{0}}{{C}^{0}} - q{{C}_{2}}.$Результаты решения этой системы без учета шума (ошибки эксперимента) для разного числа измерений и границ информативных фрагментов приведены в табл. 1.
Таблица 1.
N | k+1 | k−1 | k+2 | k−2 | Ф1 | Ф2 | R, % | E, % |
---|---|---|---|---|---|---|---|---|
3 | 0.4677 | −0.4879 | 2.1659 | 10.5642 | 1−2 | 2−3 | 1.51 | 3.19 |
5 | 0.7438 | 0.4267 | 1.7183 | 11.1214 | 1−2 | 2−5 | 0.94 | 2.34 |
7 | 0.8975 | 1.0533 | 1.4612 | 10.3054 | 1−2 | 2−7 | 1.15 | 3.27 |
13 | 0.8975 | 1.0533 | 1.4612 | 10.3054 | 1−3 | 3−13 | 0.85 | 3.17 |
31 | 0.8975 | 1.0533 | 1.4612 | 10.3054 | 1−6 | 6−31 | 0.55 | 3.11 |
Из табл. 1 следует, что с ростом числа экспериментальных точек N найденные физичные значения констант стабилизируются (метод устойчив). Ошибки метода не превышают погрешностей измерений <5%. С учетом неоднозначности выбора фрагментов сплайнов решением обратной задачи (без учета шума) являются интервалы значений констант, определенных по пяти и более точкам k+1 ∈ [0.7438, 0.8975], k−1 ∈ [0.4267, 1.0533], k+2 ∈ ∈ [1.4612, 1.7183], k−2 ∈ [10.3054, 11.1214]. Эти интервалы (кроме последнего) близки к “истинным” значениям констант скоростей стадий.
Для оценки результативности метода варьировали “истинные” значения констант для 7−13 точек. Например, для “истинных” значений $k_{{ + 1}}^{*}$ = 1, $k_{{ - 1}}^{*}$ = 0.1, $k_{{ + 2}}^{*}$ = 2, $k_{{ - 2}}^{*}$ = 0.2 метод дает k+1 = 1.0041, k−1 = 0.6996, k+2 = 3.0096, k−2 = 5.5848 (R = 1.21, E = 4.56). Для значений $k_{{ + 1}}^{*}$ = 2, $k_{{ - 1}}^{*}$ = 0.5, $k_{{ + 2}}^{*}$ = 3, $k_{{ - 2}}^{*}$ = 1 получим k+1 = 1.7180, k−1 = 0.5997, k+2 = = 5.6507, k−2 = 7.8756 (R = 1.31, E = 4.90). Даже для “жестких” значений $k_{{ + 1}}^{*}$ = 1, $k_{{ - 1}}^{*}$ = 0.1, $k_{{ + 2}}^{*}$ = 10, $k_{{ - 2}}^{*}$ = = 100 получим k+1 = 0.9796, k−1 = 0.4975, k+2 = = 6.0535, k−2 = 59.8843 (R = 1.22, E = 4.60). Как видно, метод достаточно результативен. Для оценки помехоустойчивости метода определили влияние уровня шума в исходных данных на погрешности вычислений (табл. 2).
Таблица 2.
Sn, % | k+1 | k−1 | k+2 | k−2 | R, % | E, % |
---|---|---|---|---|---|---|
0 | 0.8975 | 1.0533 | 1.4612 | 10.3054 | 0.85 | 3.17 |
1 | 0.9164 | 1.1515 | 1.4391 | 9.5460 | 0.85 | 3.19 |
2 | 0.9362 | 1.2551 | 1.4176 | 8.8192 | 0.87 | 3.26 |
3 | 0.9569 | 1.3644 | 1.3968 | 8.1232 | 0.90 | 3.38 |
4 | 0.9787 | 1.4798 | 1.3766 | 7.4565 | 0.94 | 3.54 |
5 | 1.0016 | 1.6020 | 1.3570 | 6.8175 | 1.00 | 3.74 |
10 | 1.1373 | 2.3377 | 1.2670 | 3.9923 | 1.37 | 5.12 |
15 | 1.3238 | 3.3736 | 1.1882 | 1.6837 | 1.82 | 6.84 |
20 | 1.5998 | 4.9376 | 1.1185 | –0.2162 | 2.32 | 8.71 |
Из табл. 2 следует, что с ростом уровня шума от 1 до 15% найденные значения констант остаются физичными и решениями обратной задачи (с учетом шума) являются интервалы k+1 ∈ [0.8975, 1.3238], k−1 ∈ [1.0533, 3.3736], k+2 ∈ [1.1882, 1.4612], k−2 ∈ [1.6837, 10.3054], которые ближе к “истинным” значениям констант скоростей стадий, чем интервалы (без учета шума).
Пример 2. Пусть реакция А = C протекает через три стадии
(III)
$1)\,\,{\mathbf{A}} = {\mathbf{B}},\,\,\,{\text{ }}2)\,\,{\mathbf{A}} = 2{\mathbf{C}},\,\,\,\,{\text{ }}3)\,\,{\mathbf{B}}{\text{ }} + {\text{ }}{\mathbf{C}} = {\mathbf{A}}.$Для нее уравнения (1) запишутся как
(2.1)
$\begin{gathered} A{\kern 1pt} ' = - {{k}_{ + }}_{1}A + {{k}_{{ - 1}}}B - {{k}_{ + }}_{2}A + {{k}_{{ - 2}}}{{C}^{2}} + \\ + \,\,{{k}_{ + }}_{3}BC - {{k}_{{ - 3}}}A + {{q}^{0}}{{A}^{0}} - qA, \\ \end{gathered} $(2.2)
$B{\kern 1pt} ' = {{k}_{ + }}_{1}A - {{k}_{{ - 1}}}B - {{k}_{ + }}_{3}BC + {{k}_{{ - 3}}}A + {{q}^{0}}{{B}^{0}} - qB.$(2.3)
$\begin{gathered} C{\kern 1pt} ' = {\text{ }}2{{k}_{ + }}_{2}A - 2{{k}_{{ - 2}}}{{C}^{2}} - {{k}_{ + }}_{3}BC + \\ + \,\,{{k}_{{ - 3}}}A + {{q}^{0}}{{C}^{0}} - qC.~~~ \\ \end{gathered} $В этой системе уравнений нет ЛСЗС и все реагенты независимы. Зададим “истинные” значения констант $k_{{ + 1}}^{*}$ = 1, $k_{{ - 1}}^{*}$ = 1, $k_{{ + 2}}^{*}$ = 1, $k_{{ - 2}}^{*}$ = 1, $k_{{ + 3}}^{*}$ = 1, $k_{{ - 3}}^{*}$ = 1 и н. у. q = q0 = 1, A0 = 1, B0 = 0, C0 = 0. Выберем удовлетворяющие условиям (10) две опорные точки (A1,С1) и (A2,С2) на информативных фрагментах и подставим их значения в уравнения (2.1)–(2.3):
(2.4)
$\begin{gathered} A_{1}^{'} = - {{k}_{ + }}_{1}{{A}_{1}} + {{k}_{{ - 1}}}{{B}_{1}} - {{k}_{ + }}_{2}{{A}_{1}} + {{k}_{{ - 2}}}C_{1}^{2} + \\ + \,\,{{k}_{ + }}_{3}{{B}_{1}}{{C}_{1}} - {{k}_{{ - 3}}}{{A}_{1}} + {{q}^{0}}{{A}^{0}} - q{{A}_{1}},~~~ \\ \end{gathered} $(2.5)
$B_{1}^{'} = {{k}_{ + }}_{1}{{A}_{1}} - {{k}_{{ - 1}}}{{B}_{1}} - {{k}_{ + }}_{3}{{B}_{1}}{{C}_{1}} + {{k}_{{ - 3}}}{{A}_{1}} + {{q}^{0}}{{B}^{0}} - q{{B}_{1}},$(2.6)
$\begin{gathered} C_{1}^{'} = {\text{ }}2{{k}_{ + }}_{2}{{A}_{1}} - 2{{k}_{{ - 2}}}C_{1}^{2} - {{k}_{ + }}_{3}{{B}_{1}}{{C}_{1}} + \\ + \,\,{{k}_{{ - 3}}}{{A}_{1}} + {{q}^{0}}{{C}^{0}} - q{{C}_{1}}, \\ \end{gathered} $(2.7)
$\begin{gathered} A_{2}^{'} = - {{k}_{ + }}_{1}{{A}_{2}} + {{k}_{{ - 1}}}{{B}_{2}} - {{k}_{ + }}_{2}{{A}_{2}} + {{k}_{{ - 2}}}C_{2}^{2} + \\ + \,\,{{k}_{ + }}_{3}{{B}_{2}}{{C}_{2}} - {{k}_{{ - 3}}}{{A}_{2}} + {{q}^{0}}{{A}^{0}} - q{{A}_{2}}, \\ \end{gathered} $(2.8)
$\begin{gathered} B_{2}^{'} = {{k}_{ + }}_{1}{{A}_{2}} - {{k}_{{ - 1}}}{{B}_{2}} - {{k}_{ + }}_{3}{{B}_{2}}{{C}_{2}} + \\ + \,\,{{k}_{{ - 3}}}{{A}_{2}} + {{q}^{0}}{{B}^{0}} - q{{B}_{2}}, \\ \end{gathered} $(2.9)
$\begin{gathered} C_{2}^{'} = {\text{ }}2{{k}_{ + }}_{2}{{A}_{2}} - 2{{k}_{{ - 2}}}C_{2}^{2} - {{k}_{ + }}_{3}{{B}_{2}}{{C}_{2}} + \\ + \,\,{{k}_{{ - 3}}}{{A}_{2}} + {{q}^{0}}{{C}^{0}} - q{{C}_{2}}. \\ \end{gathered} $Результаты решения этой системы без учета шума для разного числа измерений и границ информативных фрагментов приведены в табл. 3.
Таблица 3.
N | k+1 | k−1 | k+2 | k−2 | k+3 | k−3 | Ф1 | Ф 2 | R, % | E, % |
---|---|---|---|---|---|---|---|---|---|---|
3 | 0.5822 | 0.4564 | 0.6157 | 0.4356 | 0.2625 | 0.2625 | 1−2 | 2−3 | 2.11 | 4.47 |
5 | 0.8851 | 1.0229 | 0.8792 | 1.0981 | 1.1255 | 0.8971 | 1−2 | 2−5 | 1.42 | 3.56 |
9 | 0.8851 | 1.0229 | 0.8792 | 1.0981 | 1.1255 | 0.8971 | 1−3 | 3−9 | 1.09 | 3.48 |
17 | 0.9683 | 1.2922 | 0.9433 | 1.4686 | 1.6052 | 0.9507 | 1−4 | 4−17 | 1.10 | 4.69 |
33 | 0.9683 | 1.2922 | 0.9433 | 1.4686 | 1.6052 | 0.9507 | 1−7 | 7−33 | 0.79 | 4.63 |
Как видно из этой таблицы, найденные (без учета шума) константы близки к “истинным” значениям и устойчиво находятся по пяти и более точкам в интервалах k+1 ∈ [0.8851, 0.9683], k−1 ∈ ∈[1.0229, 1.2922], k+2 ∈ [0.8972, 0.9433], k−2 ∈ ∈ [1.0981, 1.4686], k+3 ∈ [1.12555, 1.6052], k−3 ∈ [0.8971, 0.9507]. Для оценки помехоустойчивости метода определили влияние уровня шума в исходных данных на погрешности вычислений (табл. 4).
Таблица 4.
Sn, % | k+1 | k−1 | k+2 | k−2 | k+3 | k−3 | R, % | E, % |
---|---|---|---|---|---|---|---|---|
0 | 0.9683 | 1.2922 | 0.9433 | 1.4686 | 1.6052 | 0.9507 | 0.79 | 4.63 |
1 | 0.8361 | 1.0961 | 0.8892 | 1.3553 | 1.9725 | 1.9725 | 0.84 | 4.91 |
2 | 0.7043 | 0.8998 | 0.8353 | 1.2439 | 2.3336 | 1.1961 | 0.89 | 5.20 |
3 | 0.5727 | 0.7035 | 0.7817 | 1.1342 | 2.6886 | 1.3179 | 0.94 | 5.51 |
4 | 0.4415 | 0.5070 | 0.7283 | 1.0263 | 3.0377 | 1.4391 | 1.00 | 5.83 |
5 | 0.3107 | 0.3105 | 0.6752 | 0.9201 | 3.3809 | 1.5598 | 1.06 | 6.17 |
10 | –0.3388 | –0.6729 | 0.4133 | 0.4141 | 5.0130 | 2.1546 | 1.36 | 7.96 |
Из табл. 4 видно, что с ростом уровня шума от 1 до 5% найденные константы остаются физичными и решениями обратной задачи (с учетом шума) являются интервалы k+1 ∈ [0.3107, 0.9683], k−1 ∈ [0.3105, 1.2922], k+2 ∈ [0.6752, 0.9433], k−2 ∈ ∈ [0.9201, 1.4686], k+3 ∈ [1.6052, 3.3809], k−3 ∈ [0.9507, 1.5598], которые близки к “истинным” значениям констант скоростей стадий.
Пример 3. В работе [27] было экспериментально установлено, что реакция гидроалюминирования олефинов для алюминийорганического соединения ${\text{HAlBu}}_{2}^{i}$ (диизобутилалюминийгидрид) в присутствии катализатора Cp2ZrCl2 протекает по схеме
(IV)
$1)\,\,{\mathbf{A}} = {\mathbf{2B}},\,\,\,\,{\text{ }}2)\,\,{\mathbf{B}} + {\mathbf{C}} = {\text{ }}{\mathbf{D}},$(3.1)
$\begin{gathered} A{\kern 1pt} ' = - {{r}_{{ + 1}}} + {{r}_{{ - 1}}},\,\,\,\,B{\kern 1pt} ' = {\text{ }}2{{r}_{{ + 1}}} - 2{{r}_{{ - 1}}} - {{r}_{{ + 2}}} + {{r}_{ - }}_{2}, \\ C{\kern 1pt} ' = - {{r}_{{ + 2}}} + {{r}_{ - }}_{2},\,\,\,D{\kern 1pt} ' = {{r}_{{ + 2}}} - {{r}_{ - }}_{2}, \\ \end{gathered} $где r+1 = k+1A, r−1 = k−1B2, r+2 = k+2BC, r−2 = k−2D. Выпишем ЛСЗС 2A + B + D = 2 и C + D = 1, выберем независимые реагенты А и D, выразим через них B и C и исключим их из выражения (3.1):
(3.2)
$\begin{gathered} A{\kern 1pt} ' = - {{k}_{ + }}_{1}A + {{k}_{{ - 1}}}{{(2 - 2A - D)}^{2}},~ \\ D{\kern 1pt} ' = {{k}_{{ + 2}}}(2 - 2A - D)(1 - D) - {{k}_{ - }}_{2}D. \\ \end{gathered} $Проинтегрируем систему (3.2) при истинных значениях констант с шагом h = 30 на интервале t ∈ [0,200] с н. у. A0 = 1, B0 = D0 = 0, C0 = 1. Получим N = 7 “экспериментальных” точек, разделим их на два фрагмента, выберем опорные точки в середине этих фрагментов и вычислим концентрации (A1,D1), (A2,D2) и скорости $A_{1}^{'},D_{1}^{'}$ и $A_{2}^{'},D_{2}^{'}$ в этих точках. Подставим эти значения в выражения (3.2) и получим линейную систему из четырех уравнений для определения четырех констант скоростей стадий:
(3.3)
$\begin{gathered} A_{1}^{'} = - {{k}_{ + }}_{1}{{A}_{1}} + {{k}_{{ - 1}}}{{(2 - 2{{A}_{1}} - {{D}_{1}})}^{2}},~ \\ D_{1}^{'} = {{k}_{ + }}_{2}(2 - 2{{A}_{1}} - {{D}_{1}})(1 - {{D}_{1}}) - {{k}_{ - }}_{2}{{D}_{1}}, \\ A_{2}^{'} = - {{k}_{ + }}_{1}{{A}_{2}} + {{k}_{ - }}_{1}{{(2 - 2{{A}_{2}} - {{D}_{2}})}^{2}},~ \\ D_{2}^{'} = {{k}_{ + }}_{2}(2 - 2{{A}_{2}} - {{D}_{2}})(1 - {{D}_{2}}) - {{k}_{ - }}_{2}{{D}_{2}}. \\ \end{gathered} $Результаты решения этой системы приведены в табл. 5.
Таблица 5.
N | Sn, % | k+1 | k−1 | k+2 | k−2 | Ф1 | Ф2 | R, % | E, % |
---|---|---|---|---|---|---|---|---|---|
7 | 0 | 0.0221 | 0.0126 | 0.0578 | 0.0018 | 1−2 | 2−7 | 2.38 | 0.86 |
7 | 1 | 0.0216 | 0.0125 | 0.0590 | 0.0019 | 1−2 | 2−7 | 2.41 | 0.86 |
7 | 5 | 0.0200 | 0.0120 | 0.0644 | 0.0024 | 1−2 | 2−7 | 2.54 | 0.86 |
7 | 10 | 0.0183 | 0.0120 | 0.0727 | 0.0032 | 1−2 | 2−7 | 2.68 | 0.89 |
7 | 20 | 0.0157 | 0.0114 | 0.0986 | 0.0054 | 1−2 | 2−7 | 2.91 | 1.19 |
Как видно из табл. 5, рассчитанные значения констант c учетом уровня шума 20% находятся в интервалах k+1 ∈ [0.0157, 0.0221], k−1 ∈ [0.0114, 0.0126], k+2 ∈ [0.0578, 0.0986], k−2 ∈ [0.0018, 0.0054], которые описывают кинетику данной реакции со стандартной ошибкой, не превышающей 3%.
ЗАКЛЮЧЕНИЕ
Описан метод решения обратной задачи химической кинетики, который позволяет определять интервалы значений констант скоростей элементарных стадий многостадийных реакций по нестационарным концентрациям реагентов с учетом погрешности измерений. Особенность метода состоит в том, что разные участки релаксационных кривых описываются линейными сплайнами с различным числом экспериментальных точек. Этот подход позволяет учитывать различную информативность исходных данных и находить интервалы физичных значений констант скоростей элементарных стадий исследуемой реакции. Эффективность метода показана на примерах двух- и трехстадийной реакций с учетом уровня шума в исходных экспериментальных данных.
Список литературы
Писаренко В.Н., Погорелов А.Г. Планирование кинетических исследований. М.: Наука, 1969. 176 с.
Слинько М.Г., Спивак С.И., Тимошенко В.И. // Кинетика и катализ. 1972. Т. 13. № 3. С. 1570.
Клибанов М.В., Спивак С.И., Тимошенко В.И., Слинько М.Г. // ДАН СССР. 1973. Т. 208. № 5. С. 1387.
Яблонский Г.С., Спивак С.И. Математические модели химической кинетики. М.: Знание, 1977. 200 с.
Методы исследования быстрых реакций / Под ред. Хеммис Г. М.: Мир, 1977. 716 с.
Дейч А.М. Методы идентификации динамических систем. М.: Энергия, 1979. 240 с.
Клибанов М.В. Обратные задачи математической физики. Куйбышев: КГУ, 1980. 42 с.
Лукашенок В.Н. // Теорет. основы хим. технологии. 1980. Т. 14. № 1. С. 86.
Спивак С.И., Горский В.Г. // ДАН СССР. 1981. Т. 257. № 2. С. 412.
Яблонский Г.С., Быков В.И., Горбань А.Н. Кинетические модели каталитических реакций. Новосибирск: Наука, 1983. 255 с.
Горский В.Г. Планирование кинетических экспериментов. М.: Наука, 1984. 241 с.
Полак Л.С., Гольденберг М.Я., Левицкий А.А. Вычислительные методы в химической кинетике. М.: Наука, 1984. 280 с.
Тихонов А.Н. Методы решения некорректных задач. М.: Наука, 1986. 288 с.
Погорелов А.Г. Обратные задачи нестационарной химической кинетики: системный подход. М.: Наука, 1988. 392 с.
Федотов А.М. Некорректные задачи со случайными ошибками в данных. Новосибирск: Наука, 1990. 279 с.
Морозов В.А., Гребенников А.И. Методы решения некорректно поставленных задач: алгоритмический аспект. М.: Изд-во МГУ, 1992. 319 с.
Денисов А.М. Введение в теорию обратных задач. М.: Изд-во МГУ, 1994. 208с.
Кольцов Н.И., Алексеев Б.В., Федотов В.Х. Всесоюз. конф. Актуальные проблемы нефтехимии. Тез. докл. Уфа, 1985. С. 105.
Кольцов Н.И., Алексеев Б.В. // ДАН СССР. 1988. Т. 298. № 2. С. 401.
Кольцов Н.И., Алексеев Б.В. // Журн. физ. химии. 1989. Т. 63. № 4. С. 1125.
Кольцов Н.И., Алексеев Б.В. // Журн. физ. химии. 1989. Т. 63. № 6. С. 1623.
Кольцов Н.И. Математическое моделирование каталитических реакций. Чебоксары: Изд-во Чуваш. ун-та, 2007. 294 с.
Исмагилова А.С., Спивак С.И. Обратные задачи химической кинетики, Saarbrucken: Lap Lambert Academic Publishing, 2013. 118 с.
Федотов В.Х., Кольцов Н.И. // Хим. физика. 2016. Т. 35. № 10. С. 9.
Роджерс Д., Адамс Дж. Математические основы машинной графики. М.: Мир, 2001. 604 с.
Алексеев Б.В., Кольцов Н.И., Федотов В.Х. // Журн. физ. химии. 1992. Т. 66. № 12. С. 3219.
Вильданова Р.Ф. Дис. … канд. хим. наук. Уфа: Баш. ГУ, 2007.
Дополнительные материалы отсутствуют.
Инструменты
Кинетика и катализ