一种分子水平反应动力学模型构建方法及装置、存储介质与流程

1.本公开涉及分子炼油技术领域，尤其涉及一种分子水平反应动力学模型构建方法及装置、存储介质及设备。


背景技术：

2.目前，炼油过程装置建模方法主要采用集总模型方法，依据将物料中各组分的宏观物性、结构特征等动力学性质进行相似性归类，再进行各集总组分间反应网络的构建和反应参数的计算。然而，每一个集总组分实际由大量的纯分子组成，实际上在建模过程中将其视为具有均匀物性的、虚拟的单一组分。因此，集总模型在新物料和新催化剂的扩展适应性上存在明显缺陷。


技术实现要素：

3.为了解决上述技术问题或者至少部分地解决上述技术问题，本公开的实施例提供了一种分子水平反应动力学模型构建方法及装置、存储介质及设备。
4.第一方面，本公开的实施例提供了一种分子水平反应动力学模型构建方法，所述方法包括：
5.获取原料分子组成矩阵，其中，所述原料分子组成矩阵包括每种原料分子的结构导向集总表示和含量；
6.基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，其中，所述反应规则包括每种原料分子在其对应的反应路径中结构导向集总表示的变化；
7.将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量，并根据产物分子的含量计算每种产物的收率；
8.计算每种产物的收率与实际反应产物的收率之间的相对差值，并根据所述相对差值对预设的反应速率常数进行调整，将所述相对差值之和满足预设条件时的反应速率常数作为分子水平反应动力学模型的参数，使得分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵。
9.在一种可能的实施方式中，所述分子水平反应动力学模型为常减压蒸馏模型、渣油加氢模型、催化裂化模型、催化裂解模型、延迟焦化模型、加氢裂化模型、催化重整模型、烷基化模型、汽油加氢模型、柴油加氢模型、蜡油加氢模型、汽柴油加氢模型、气体分馏模型、芳烃抽提模型和制氢模型中的其中一种。
10.在一种可能的实施方式中，所述基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，包括：
11.按照预设的反应规则对每种原料分子的结构导向集总表示进行遍历，得到每种原料分子对应的反应路径以及由多条反应路径组成的反应网络。
12.在一种可能的实施方式中，所述将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数初值和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量，包括：
13.将当前反应的反应时长平均分为多个微元反应段；
14.将当前反应的实际生产单元入口的反应温度t、压力p、焓值h和比热容c作为第一个微元反应段入口的初始温度、初始压力、初始焓值和比热容；
15.对于每一个微元反应段，利用龙格库塔算法求解每条反应规则对应的反应速率常数下的常微分方程组，得到该微元反应段出口的产物浓度及浓度随时间t的变化；
16.根据和微元反应段入口的初始焓值h，计算焓值h随时间t的变化，以得到该微元反应段出口的焓值；
17.根据和比热容c，计算温度t随时间t的变化，以得到该微元反应段出口的温度；
18.根据理想气体常数r和温度t随时间的变化计算压力p随时间t的变化，以得到该微元反应段出口的压力；
19.将该微元反应段出口的温度、压力、焓值、各个产物浓度作为下一个微元反应段的进口参数，计算下一个微元反应段出口的参数，直到最后一个微元反应段，将最后一个微元反应段出口的产物分子的含量作为该生产单元的产物分子的含量。
20.在一种可能的实施方式中，所述反应时长为原料分子在生产单元中停留的时长。
21.在一种可能的实施方式中，通过以下表达式，计算每种产物的收率与实际反应产物的收率之间的相对差值，并将相对差值进行求和，包括：
[0022][0023]
err为每种产物的计算收率与实际收率之间的相对差值求和，为第i种产物的计算收率，为第i种产物的实际收率，为第i种产物的计算收率与实际收率之间的相对差值，n为产物的种类总数。
[0024]
在一种可能的实施方式中，通过以下表达式，计算所述反应速率常数：
[0025][0026]
其中，k为反应规则的反应速率常数，ka、kb、kc分别为与催化剂、反应温度、反应压力相关的反应动力学参数，e为反应活化能，t为反应温度，p为反应压力，pk为反应压力对反应速率常数影响的常数。
[0027]
在一种可能的实施方式中，根据所述相对差值对预设的反应速率常数进行调整，将所述相对差值之和满足预设条件时的反应速率常数作为分子水平反应动力学模型的参数，包括：
[0028]
当所述相对差值大于预设阈值，调整动力学参数，得到调整后动力学参数，基于所述调整后动力学参数计算反应速率常数，并根据反应速率常数和每种原料分子的含量重新求解所述常微分方程，得到每种产物分子的含量；
[0029]
其中，预设条件指根据实际应用场景所能接受的最大误差。
[0030]
在一种可能的实施方式中，所述分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵，包括：
[0031]
基于龙格库塔算法，根据所述相对差值之和满足预设条件时的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量。
[0032]
在一种可能的实施方式中，所述获取原料分子组成矩阵，包括：
[0033]
基于结构导向集总分子表征方法对每种原料分子进行向量表征，得到每种原料分子的结构导向集总表示；
[0034]
将每种原料分子的结构导向集总表示与含量作为完整向量；
[0035]
将所有原料分子的完整向量组合为原料分子组成矩阵。
[0036]
第二方面，本公开的实施例提供了一种分子水平反应动力学模型构建装置，包括：
[0037]
获取模块，其用于获取原料分子组成矩阵，其中，所述原料分子组成矩阵包括每种原料分子的结构导向集总表示和含量；
[0038]
生成模块，其用于基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，其中，所述反应规则包括每种原料分子在其对应的反应路径中结构导向集总表示的变化；
[0039]
求解模块，其用于将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量；
[0040]
调整模块，其用于计算每种产物的收率与实际反应产物的收率之间的相对差值，并根据所述相对差值对预设的反应速率常数进行调整，将所述相对差值之和满足预设条件时的反应速率常数作为分子水平反应动力学模型的参数，使得分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵。
[0041]
第三方面，本公开的实施例提供了一种电子设备，包括处理器、通信接口、存储器和通信总线，其中，处理器、通信接口和存储器通过通信总线完成相互间的通信；
[0042]
存储器，用于存放计算机程序；
[0043]
处理器，用于执行存储器上所存放的程序时，实现上述的分子水平反应动力学模型构建方法。
[0044]
第四方面，本公开的实施例提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的分子水平反应动力学模型构建方法。
[0045]
本公开实施例提供的上述技术方案与现有技术相比至少具有如下优点的部分或全部：
[0046]
本公开实施例所述的分子水平反应动力学模型构建方法，基于原料分子水平表征，在不依赖原料组分的条件下，针对不同的过程反应机理，实现对反应进程中分子组成的转化跟踪，对炼化过程各装置反应过程具有清晰解释，实现产物性质预测。
附图说明
[0047]
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。
[0048]
为了更清楚地说明本公开实施例或现有技术中的技术方案，下面将对实施例或相
关技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0049]
图1示意性示出了根据本公开实施例的分子水平反应动力学模型构建方法流程示意图；
[0050]
图2示意性示出了根据本公开实施例的分子水平反应动力学模型构建装置的结构框图；
[0051]
图3示意性示出了根据本公开实施例的电子设备的结构框图。
具体实施方式
[0052]
为使本公开实施例的目的、技术方案和优点更加清楚，下面将结合本公开实施例中的附图，对本公开实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本公开的一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本公开保护的范围。
[0053]
参见图1，本公开的实施例提供了一种分子水平反应动力学模型构建方法，所述方法包括：
[0054]
s1，获取原料分子组成矩阵，其中，所述原料分子组成矩阵包括每种原料分子的结构导向集总表示和含量；
[0055]
在一些实施例中，原料分子组成矩阵通过以下步骤获取：包括原料分子组成结构导向集总表示矩阵和原料分子组成含量矩阵，其中，结构导向集总表示如下表1所示。
[0056]
表1
[0057]
a6a4a2n6n5n4n3n2n1rbrmeihaansrsannnrnnorokoniv
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
[0058]
用上述24种特征结构根据分子导向集总表示法中对应的规则表示原料分子组成，每个分子转化为24个片段单元的数量构成的一维结构向量，分子向量记为a=[a
11
, a
12
,
ꢀ……a1n
]。a中各个分量表示该分子中对应表1中特征结构的数量，原料用一个维数为n
×
24分子矩阵和一个维数为1
×
n分子含量向量表示，分子矩阵记为b=，分子含量向量记为c=[c1,c2,
ꢀ……cn
]，其中n为原料中的分子种类，分子矩阵b每一行为一个分子向量，分子含量向量c的每个分量对应分子矩阵b每一行分子向量的含量，表示原料中对应分子的摩尔分数；
[0059]
a6：在所有芳香族分子中都出现的一个六碳芳香环，可以单独存在。
[0060]
a4：一种附在a6（或另一个a4环）上的四碳芳香环，这是结构增量，用以构建聚合的多环结构，不能单独存在。
[0061]
a2：二碳芳香族结构增量，a2用于附加到多环芳烃的“海湾区域”形成新的多环芳烃。
[0062]
n6和n5：六碳和五碳环烷烃。
[0063]
n4、n3、n2、n1：额外的包含四个、三个、两个和一个碳的脂肪环结构增量，它们必须附着在其他脂肪环或芳香环结构中,不能单独存在。
[0064]
r：连接在环结构上的所有烷基结构包含的碳原子个数，或是没有环结构存在时脂
肪族分子中的碳原子个数。
[0065]
ih：引入氢元素相关的结构增量，来描述分子饱和度。如果没有环结构，ih=1表示链烷烃，ih=0表示单烯烃，ih=-1表示双烯烃；如果有环存在，ih=-1表示环烯烃。
[0066]
br：表示侧链烷基、直链烷基或烯烃上分支节点的个数，它不能区分甲基、乙基、丙基支链，所以统一假定只存在甲基支链，实际炼油过程中，考虑支链种类对反应的影响意义不大，而以上假设能够表现分支个数的影响，忽略分支支链种类，如甲基、乙基、丙基支链的影响，能够满足实际需要。
[0067]
me：确定烷基结构中直接连接在芳香环或脂肪环中的碳原子的甲基基团的数量。特别地，当r=1或me=r-1时，其他结构增量能够确定环结构上的甲基数目，根据约定me不在用于表示甲基的数量。
[0068]
aa：任意两个非结构增量环（a6、n6或n5）间的联苯桥连结构。
[0069]
ns、nn和no：位于脂肪环或脂肪链中且与两个碳原子相连的硫、氮、氧原子。ns、nn和no分别指用s原子、-nh-基团和o原子替代一个-ch
2-。
[0070]
rs、rn和ro：在一个碳原子和氢原子之间插入一个s原子、含n的-nh-基团或o原子，分别构成硫醇、胺或醇类基团。
[0071]
an：在芳香环中用氮基团替代碳，如吡啶和喹啉。an基团用=n-替代了=ch-。
[0072]
ko：替代了-ch
2-或-ch3形成酮或醛基。
[0073]
ni、v：在卟啉类分子中出现。
[0074]
s2，基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，其中，所述反应规则包括每种原料分子在其对应的反应路径中结构导向集总表示的变化；
[0075]
s3，将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量，并根据产物分子的含量计算每种产物的收率；
[0076]
在一些实施例中，将当前反应划分为多个微元反应段中，根据反应停留时间划分微元反应段，在每个微元反应段里面，根据所述常微分方程计算反应网络中每个分子物质的量浓度随时间的变化，其中，常微分方程是根据反应网络自动生成的，描述了反应网络里面所有分子的浓度随时间变化的方程，另外，反应网络是根据原料分子和反应规则自动生成的。
[0077]
计算每种产物的收率与实际反应产物的收率之间的相对差值，并根据所述相对差值对预设的反应速率常数进行调整，将所述相对差值之和满足预设条件时的反应速率常数作为分子水平反应动力学模型的参数，使得分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵，其中，利用内德-米德（nelder
–
mead）非线性优化算法对预设的反应速率常数进行调整。
[0078]
在一些实施例中，在根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络之后，所述方法还包括：
[0079]
将每个反应路径的每条反应路径的产物分子的结构导向集总表示与预设分子集进行对比，仅保留在预设分子集中存在的产物分子及其对应的反应路径，作为有效产物分子和有效反应路径，以计算每种有效产物分子的含量，其中，所述预设分子集为基于结构导
向集总方法，根据炼油过程各装置实际发生的反应及分子产物预定义得到的，用于对模型基于反应规则生成的产物进行初步判断。
[0080]
在一种可能的实施方式中，所述分子水平反应动力学模型为蒸馏模型、渣油加氢模型、催化裂化模型、催化裂解模型、延迟焦化模型、加氢裂化模型、催化重整模型、烷基化模型、汽油加氢模型、柴油加氢模型、蜡油加氢模型、汽柴油加氢模型、气分模型、芳烃抽提模型和制氢模型中的其中一种。
[0081]
在一些实施例中，对于每一种模型，根据模型内实际发生的化学反应制定相对应的反应规则，得到每一个模型中对应原料分子组成中每种单分子对应的反应路径。
[0082]
在一些实施例中，原料输入时，首先，依据原料实沸点（tbp）蒸馏曲线数据进行原油馏分的切割计算，即利用馏分的沸程差别，通过蒸馏工艺从原油中分离出石脑油、柴油，以及蜡油和渣油等，分别对应不同模型反应规则。每种原料分子按照对应模型反应规则集合中的反应规则进行反应，得到每种分子对应的反应路径。每种分子按照反应规则进行第一次反应生成中间产物后，则继续作为反应物判断是否满足另一条反应规则继续按照反应规则进行后续反应，直至该中间产物不符合反应规则集合中的任一反应规则，则该中间产物为最终产物，所经历反应的汇总即该分子的反应路径，生成对应反应网络。
[0083]
在一些实施例中，所述基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，包括：
[0084]
按照预设的反应规则对每种原料分子的结构导向集总表示进行遍历，得到每种原料分子对应的反应路径以及由多条反应路径组成的反应网络。
[0085]
在一些实施例中，所述将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量，包括：
[0086]
将当前反应的反应时长平均分为n个微元反应时长段，每个微元反应段的反应时长为；
[0087]
将当前反应的实际生产单元入口的反应温度t、压力p、焓值h和比热容c作为第一个微元反应段入口的初始温度、初始压力、初始焓值和比热容；
[0088]
对于每一个微元反应段，利用龙格库塔算法求解每条反应规则对应的反应速率常数ki下的常微分方程组，得到该微元反应段出口的反应物浓度及浓度c随时间的变化，其中，常微分方程组包括：零级反应，；一级反应，；二级反应，；
[0089]
根据和微元反应段入口焓值，计算焓值随时间的变化；
[0090]
根据和比热容，计算温度随时间的变化，以得到该微元反应段出口的温度；
[0091]
根据理想气体常数r和温度随时间的变化计算压力随时间的变化，以得到该微元反应段出口的压力；
[0092]
将该微元反应段出口的温度、压力、焓值、各个产物浓度作为下一个微元反应段的进口参数，计算下一个微元反应段出口的参数，直到最后一个微元反应段，将最后一个微元反应段出口的产物含量作为该生产单元的产物的含量。
[0093]
在一些实施例中，通过以下表达式，计算每种产物的收率与实际反应产物的收率之间的相对差值，并将相对差值进行求和，包括：
[0094][0095]
其中，err为每种产物的计算收率与实际收率之间的相对差值求和，为第i种产物的计算收率（质量分数），为第i种产物的实际收率（质量分数），为第i种产物的计算收率与实际收率之间的相对差值，n为产物的种类总数。
[0096]
在一些实施例中，通过以下表达式，计算所述反应速率常数：
[0097][0098]
其中，k为反应规则的反应速率常数，ka、kb、kc分别为与催化剂、反应温度、反应压力相关的反应动力学参数，e为反应活化能，t为反应温度，p为反应压力，pk为反应压力对反应速率常数影响的常数。
[0099]
在一些实施例中，所述根据所述相对差值对预设的反应速率常数进行调整，包括：
[0100]
响应于所述相对差值大于预设阈值，调整动力学参数，并根据调整后动力学参数计算的反应速率常数和每种原料分子的含量重新求解所述常微分方程，得到每种产物分子的含量。
[0101]
在一些实施例中，所述分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵，包括：
[0102]
基于拉道radau iia算法，根据所述相对差值之和满足预设条件时的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量。
[0103]
在一些实施例中，所述原料分子组成通过以下步骤得到：
[0104]
从原油分子组成中借助切割方法得到，或者，
[0105]
对原油进行蒸馏切割后得到的一次或二次加工产物的分子组成，其中，加工产物的分子组成通过气相色谱-质谱、全二维气相色谱法、四级杆气相色谱-质谱仪检测法、气相色谱或场电离-飞行时间质谱检测法、近红外光谱法、核磁共振波谱法、拉曼光谱法、傅立叶变换离子回旋共振质谱法、静电场轨道阱质谱法和离子淌度质谱法中的一种或多种确定。
[0106]
在一些实施例中，所述获取原料分子组成矩阵，包括：
[0107]
基于结构导向集总分子表征方法对每种原料分子进行向量表征，得到每种原料分子的结构导向集总表示；
[0108]
将每种原料分子的结构导向集总表示与含量作为完整向量；
[0109]
将所有原料分子的完整向量组合为原料分子组成矩阵。
[0110]
本公开的分子水平反应动力学模型构建方法基于结构导向集总方法，通过制定不同装置分子尺度反应规则，基于微元反应段划分建立分子水平炼化装置模拟模型，实现炼油加工过程复杂反应网络自动生成，实现物料平衡、能量平衡，并根据分子集的预设验证产物分子的有效性，大幅减少计算量的同时，也使反应更加贴合实际，根据实际装置原料分子组分及含量，实现对反应进程中分子组成的转化跟踪，对炼化过程各装置反应过程具有清晰解释，实现产物性质预测。此外，该模型可以为分子水平炼油过程中任意一个装置，具备
一定的通用性，分子水平的计算增加了模型的计算精度，其计算结果能够很好地预测产品收率，对炼化过程生产起到很好指导作用，提高炼化过程各装置的生产效率。
[0111]
参见图2，本公开的实施例提供了一种分子水平反应动力学模型构建装置，包括：
[0112]
获取模块11，其用于获取原料分子组成矩阵，其中，所述原料分子组成矩阵包括每种原料分子的结构导向集总表示和含量；
[0113]
生成模块12，其用于基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，其中，所述反应规则包括每种原料分子在其对应的反应路径中结构导向集总表示的变化；
[0114]
求解模块13，其用于将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量；
[0115]
调整模块14，其用于计算每种产物的收率与实际反应产物的收率之间的相对差值，并根据所述相对差值对预设的反应速率常数进行调整，将所述相对差值之和满足预设条件时的反应速率常数作为分子水平反应动力学模型的参数，使得分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵。
[0116]
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程，在此不再赘述。
[0117]
对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本发明方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
[0118]
上述实施例中，获取模块11、生成模块12、求解模块13和调整模块14中的任意多个可以合并在一个模块中实现，或者其中的任意一个模块可以被拆分成多个模块。或者，这些模块中的一个或多个模块的至少部分功能可以与其他模块的至少部分功能相结合，并在一个模块中实现。获取模块11、生成模块12、求解模块13和调整模块14中的至少一个可以至少被部分地实现为硬件电路，例如现场可编程门阵列（fpga）、可编程逻辑阵列（pla）、片上系统、基板上的系统、封装上的系统、专用集成电路（asic），或可以通过对电路进行集成或封装的任何其他的合理方式等硬件或固件来实现，或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者，获取模块11、生成模块12、求解模块13和调整模块14中的至少一个可以至少被部分地实现为计算机程序模块，当该计算机程序模块被运行时，可以执行相应的功能。
[0119]
参见图3，本公开的实施例提供的电子设备，包括处理器1110、通信接口1120、存储器1130和通信总线1140，其中，处理器1110，通信接口1120，存储器1130通过通信总线1140完成相互间的通信；
[0120]
存储器1130，用于存放计算机程序；
[0121]
处理器1110，用于执行存储器1130上所存放的程序时，实现如下所示分子水平反应动力学模型构建方法：
[0122]
获取原料分子组成矩阵，其中，所述原料分子组成矩阵包括每种原料分子的结构导向集总表示和含量；
[0123]
基于预设的反应规则，根据每种原料分子的结构导向集总表示生成多条反应路径以及由多条反应路径组成的反应网络，其中，所述反应规则包括每种原料分子在其对应的反应路径中结构导向集总表示的变化；
[0124]
将当前反应划分为多个微元反应段，并对每一个微元反应段根据反应网络构建常微分方程，根据预设的反应速率常数和每种原料分子的含量求解所述常微分方程，得到每种产物分子的含量；
[0125]
计算每种产物的收率与实际反应产物的收率之间的相对差值，并根据所述相对差值对预设的反应速率常数进行调整，将所述相对差值之和满足预设条件时的反应速率常数作为分子水平反应动力学模型的参数，使得分子水平反应动力学模型根据反应原料的分子组成矩阵预测反应产物的分子组成矩阵。
[0126]
上述的通信总线1140可以是外设部件互连标准(peripheralcomponent interconnect，简称pci)总线或扩展工业标准结构(extendedindustry standard architecture，简称eisa)总线等。该通信总线1140可以分为地址总线、数据总线、控制总线等。为便于表示，图中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。
[0127]
通信接口1120用于上述电子设备与其他设备之间的通信。
[0128]
存储器1130可以包括随机存取存储器(random accessmemory，简称ram)，也可以包括非易失性存储器(non-volatilememory)，例如至少一个磁盘存储器。可选的，存储器1130还可以是至少一个位于远离前述处理器1110的存储装置。
[0129]
上述的处理器1110可以是通用处理器，包括中央处理器(centralprocessing unit，简称cpu)、网络处理器(networkprocessor，简称np)等；还可以是数字信号处理器(digitalsignal processing，简称dsp)、专用集成电路(applicationspecific integrated circuit，简称asic)、现场可编程门阵列(field-programmable gate array，简称fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
[0130]
本公开的实施例还提供了一种计算机可读存储介质。上述计算机可读存储介质上存储有计算机程序，上述计算机程序被处理器执行时实现如上所述的分子水平反应动力学模型构建方法。
[0131]
该计算机可读存储介质可以是上述实施例中描述的设备/装置中所包含的；也可以是单独存在，而未装配入该设备/装置中。上述计算机可读存储介质承载有一个或者多个程序，当上述一个或者多个程序被执行时，实现根据本公开实施例的分子水平反应动力学模型构建方法。
[0132]
根据本公开的实施例，计算机可读存储介质可以是非易失性的计算机可读存储介质，例如可以包括但不限于：便携式计算机磁盘、硬盘、随机访问存储器（ram）、只读存储器（rom）、可擦式可编程只读存储器（eprom或闪存）、便携式紧凑磁盘只读存储器（cd-rom）、光存储器件、磁存储器件、或者上述的任意合适的组合。在本公开中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
[0133]
需要说明的是，在本文中，诸如“第一”和“第二”等之类的关系术语仅仅用来将一
个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0134]
以上所述仅是本公开的具体实施方式，使本领域技术人员能够理解或实现本公开。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本公开的精神或范围的情况下，在其它实施例中实现。因此，本公开将不会被限制于本文所示的这些实施例，而是要符合与本文所申请的原理和新颖特点相一致的最宽的范围。