CN113628692B - 基于数据与机理相混合的催化重整过程建模方法 - Google Patents
基于数据与机理相混合的催化重整过程建模方法 Download PDFInfo
- Publication number
- CN113628692B CN113628692B CN202110599815.6A CN202110599815A CN113628692B CN 113628692 B CN113628692 B CN 113628692B CN 202110599815 A CN202110599815 A CN 202110599815A CN 113628692 B CN113628692 B CN 113628692B
- Authority
- CN
- China
- Prior art keywords
- model
- reactor
- reaction
- catalytic reforming
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 132
- 230000007246 mechanism Effects 0.000 title claims abstract description 74
- 230000008569 process Effects 0.000 title claims abstract description 63
- 238000001833 catalytic reforming Methods 0.000 title claims abstract description 54
- 239000000203 mixture Substances 0.000 title claims description 11
- 238000013499 data model Methods 0.000 claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 26
- 238000013528 artificial neural network Methods 0.000 claims abstract description 11
- 238000012549 training Methods 0.000 claims abstract description 11
- 238000005094 computer simulation Methods 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims description 72
- 238000002407 reforming Methods 0.000 claims description 36
- 239000013598 vector Substances 0.000 claims description 33
- 239000007789 gas Substances 0.000 claims description 14
- 150000001335 aliphatic alkanes Chemical class 0.000 claims description 12
- 239000003054 catalyst Substances 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 11
- 150000004945 aromatic hydrocarbons Chemical class 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 239000001257 hydrogen Substances 0.000 claims description 9
- 229910052739 hydrogen Inorganic materials 0.000 claims description 9
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 8
- 150000001924 cycloalkanes Chemical class 0.000 claims description 8
- 238000006057 reforming reaction Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 7
- 230000002441 reversible effect Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000006356 dehydrogenation reaction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004886 process control Methods 0.000 claims description 6
- 230000000704 physical effect Effects 0.000 claims description 5
- 239000000376 reactant Substances 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 150000001336 alkenes Chemical class 0.000 claims description 3
- 238000004517 catalytic hydrocracking Methods 0.000 claims description 3
- 238000007327 hydrogenolysis reaction Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- UFWIBTONFRDIAS-UHFFFAOYSA-N Naphthalene Chemical compound C1=CC=CC2=CC=CC=C21 UFWIBTONFRDIAS-UHFFFAOYSA-N 0.000 claims 1
- 238000004590 computer program Methods 0.000 claims 1
- 150000002431 hydrogen Chemical class 0.000 claims 1
- 230000007547 defect Effects 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 description 8
- 241000196324 Embryophyta Species 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000003062 neural network model Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000007795 chemical reaction product Substances 0.000 description 2
- 238000011438 discrete method Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 241000183024 Populus tremula Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开基于数据与机理相混合的催化重整过程建模方法。建立机理模型；将机理模型进行离散，得到离散模型；采用能够求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数；构建数据模型，训练时采用神经网络学习方式；最后，训练好的数据模型以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，根据其输出与催化重整过程中反应器的输出实测值相比较，若差值大于阈值，则重新优化机理模型，反之则结束，得到所需的机理模型和数据模型相结合的催化重整过程模型。该方法弥补了机理建模与数据建模各自的不足之处，提高了催化重整过程模型的通用性与适应性。
Description
技术领域
本发明属于自动控制技术领域，涉及一种基于数据与机理相混合的催化重整过程建模方法。
背景技术
催化重整是炼油工业主要的加工工艺之一，也是解决高辛烷值汽油、芳烃原料短缺的有效方法之一。只有建立更加准确的模型，才会使催化重整装置运行的更好，也会使得产率提高，得到更好的经济收益。
重整反应的过程是重整进料与循环氢混合，并与重整反应产物进行换热，换热后进入加热炉加热到一定的温度后，进入第一个重整反应器进行重整反应，生成重整反应产物从第一个重整反应器出来后再进入加热炉加热到一定的反应温度后进入第二个重整反应器进行重整反应，依次循环直到第四个重整反应器。
针对此工艺流程在机理建模方面，国外已有很多成熟的机理建模软件，例如AspenPlus，ProⅡ，UniSim等。国外的这些建模软件，是为了针对其自己国家的实际情况所建立的，且其软件内部的机理是严格保密的，无法获取到，这就导致了其中存在的一些模型针对国内的实际情况来说不一定适用，因此建立出符合我国内实际情况的机理模型是十分必要的。机理模型的优点是可以很好的获得内部的参数，且参数易于调整，模型具有很强的适应性，但是其难点在于催化重整过程内部过于复杂且有些参数无法获得，因此机理模型想要建立的准确非常困难，此外在优化求解中还需要采用大规模的求解器，求解速度较慢。
近些年来，随着机器学习、人工智能等技术的发展，国内外很多专家学者逐渐运用此类技术来建立催化重整过程的数据模型，以此来克服传统机理建模方法的不足，其建立的数据模型可以看为一个黑箱模型，此模型不需要对机理过程的内部结构了解太多且不需要求解器，但是仅使用数据模型的话，如果数据未来的变化一旦超出一定的范围，此类模型的可靠性将会大幅度下降。
随着世界能源短缺以及越来越激烈的国际市场竞争，与发达国家同类企业相比，国内企业普遍存在着生产成本高，经济效益差等问题。因此，若可以建立较为通用且准确的模型，不仅可以节省资金，还可以提高市场竞争力。基于机理模型与数据模型各自存在的问题，本发明给出了一种基于数据与机理相混合的催化重整过程建模方法，该方法采用机理建模、神经网络学习相结合的方式，为建立出符合实际需求的催化重整过程模型提供了有利的保障。
发明内容
本发明的目的是针对现有的催化重整过程中机理建模与数据建模的各自存在的缺点，提出了一种机理模型与数据模型相混合的建模方法。
本发明依据重整动力学、热力学、物料平衡、能量平衡等基本原理，建立了较为完备的机理模型；为了使得机理模型与实际生产过程达到一致，根据实际工厂的输入数据，采用内点法的参数校正方法来对模型参数进行优化；将优化后机理模型进行离散，得到离散模型；采用能够求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数；构建数据模型，训练时采用神经网络学习方式；最后，训练好的数据模型以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，根据其输出与催化重整过程中反应器的输出实测值相比较，若差值大于阈值，则重新优化机理模型，反之则结束，得到所需的机理模型和数据模型相结合的催化重整过程模型。
该方法弥补了机理建模与数据建模各自的不足之处，提高了催化重整过程模型的通用性与适应性。
一种基于数据与机理相混合的催化重整过程建模方法，该方法主要步骤如下：
步骤1、依据重整动力学、热力学、物料平衡、能量平衡等基本原理，建立机理模型；
步骤2、为了使得机理模型与实际生产过程达到一致，根据催化重整过程中反应器的输入实测值，采用内点法的参数校正方法来对步骤(1)机理模型的模型参数进行优化；
步骤3、将步骤(2)优化后机理模型进行离散，得到离散模型；采用能够求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数；
步骤4、构建数据模型，训练时采用神经网络学习方式，以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，以催化重整过程中反应器的输出实测值为输出；其中催化重整过程中反应器的输入实测值为重整进料各集总组分含量、催化剂含量等，催化重整过程中反应器的输出实测值为反应器输出产物含量等。
步骤5、训练好的数据模型以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，根据其输出与催化重整过程中反应器的输出实测值相比较，若差值大于阈值，则返回步骤2，反之则结束，得到所需的优化后机理模型和数据模型相结合的混合模型。
所述的步骤1具体的实施步骤如下：
重整反应器是整个重整装置中的核心部分，机理建模最关键的部分就是针对重整反应器的建模，集总方法是反应器建模的必然选择。本发明将反应物料划分为33集总，其中集总的划分主要为烷烃、环烷烃、芳烃、烯烃、氢气等，根据其反应过程可以得出39个反应速率方程，如下所示：
烷烃脱氢反应(可逆)：
环烷烃脱氢反应(可逆)：
芳烃氢解反应：
烷烃加氢裂化反应：
其中为第jz个反应的反应速率；YP、YN、YA分别为烷烃、环烷烃、芳烃的摩尔流量；VC为催化剂装填量；F为进料体积流量；Kepj为第jz个反应平衡常数；/>为反应速率常数，其表达式为：
其中为反应频率因子；/>分别为反应活化能和压力指数；/>为/>压力指数下反应物中的氢分压；R'为摩尔气体常量；T为反应温度；/>为催化剂活性因子；/>为重整反应器装置因子。
本发明连续重整工艺采用常用的径向反应器，径向反应器物流通过反应器的压力降比轴向反应器小，有利于减小临氢系统的压力降。结合物料平衡和能量平衡原理，可得到重整反应器模型方程组：
其中为33集总组分的摩尔流量；R为径向反应器床层半径；H为反应器高度；液时空速LHSV＝F/Vc；/>为各集总组分的化学计量系数矩阵；/>为反应热向量；/>为恒压比热容向量；/>为重整反应的反应速率向量；
由于因此重整反应器机理模型方程组可以改写为：
其中为反应速率常数矩阵。
重整反应器机理模型见方程组(7)中的基础物性数据具体如下：
在重整模拟过程中，必要物性数据的获取是模型计算的前提条件，这些数据的准确性对模拟精度有非常大的影响。一些热力学性质数据是必须通过关联式由基础数据计算得到的。在重整温度条件下，所有反应组分均呈气态形式，可以将它们看作为理想气体状态，使用热力学一致性方程。
1)焓
Hm＝A+BT+CT2+DT3+ET4+F'T5 (8)
其中Hm为理想气体在T时的焓，T为反应温度；A、B，C，D，E，F'为导出系数。
2)熵
其中S为理想气体在T时的熵；G为导出系数。
3)恒压比热容
其中Cp为理想气体在恒压下的热容。
4)反应热
在温度T下，反应热的计算公式为：
其中代表了标准摩尔反应焓；α为化学计量系数；/>代表了标准摩尔生成焓；Cp(T)为各组分在温度T下的恒压比热容。
5)反应平衡常数计算
由于前15个反应是可逆反应，因此上述反应的温度T及标准压力为100Kpa下的平衡常数为：
其中为温度T下的标准摩尔生成吉布斯自由能，也称为标准摩尔自由焓，计算式为：
其中已知，标准摩尔生成熵/>的计算公式为：
其中为标准状态下的组分标准摩尔熵。
所述的步骤2具体的实施步骤如下：
1)为了使得机理模型与实际生产过程达到一致，根据实际工厂的输入数据，采用内点法的参数校正方法来对模型参数进行获取，设定一个待估参数向量X，X表示为反应的指前因子以及装置性能系数，其目标函数可以设定为：
其中G(X)表示催化重整过程中反应器的实测值与催化重整过程中反应器的设定值之间的关系，表示第in个约束方程，m'表示约束方程的总个数；
根据目标函数(15)构造惩罚项为允许误差为ε；
2)初始化X和惩罚函数因子νk，,设定迭代次数k＝1；
3)求解下列无约束问题：
min G(X)+νkβ(X) (16)
得到X极小值；
4)判断是否满足νkβ(Xk)＜ε，若是则结束，当前迭代次数下得到的X极小值即为最优解；反之令k＝k+1，νk＝C"νk-1，C"∈(0,1)，且返回步骤3)。
所述的步骤3具体的实施步骤如下：
方程组(7)属于既包含微分方程又包含代数方程的微分代数优化命题，基于有限元配置的全离散联立求解技术因其独特的优势非常合适求解该类问题。采用有限元配置方法对包含微分代数方程的机理模型方程组(7)进行离散，得到离散模型，即得到许多离散变量以及方程。采用求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数，例如催化剂吸附平衡常数、反应器装置因数等；
作为优选，所述有限元配置方法具有龙格-库塔离散方法的多种优势，在该方法中微分状态变量表达式如下：
其中w(z)表示当前有限元微分状态变量；表示在第iy个有限元初始处的值；/>是有限元iy的长度；/>表示在第iy个有限元上配置点qy处的一阶导数；/>为满足以下条件的K阶多项式：
其中为有限元中第ry个配置点的位置。qy＝1,...,K,ry＝1,...,K
微分方程的连续性方程表示为：
此外，对于过程控制变量和代数变量可采用拉格朗日形式的插值多项式表示。
代数变量形式如下：
过程控制变量形式如下：
其中和/>分别表示第iy个有限元置点qy处的代数变量和控制变量值，且满足z表示空间中的位置，/>表示第iy配置点所在的空间位置，/>为K阶拉格朗日多项式。
所述的步骤4中神经网络包括输入层、隐含层(中间层)和输出层，隐含层(对输入特征多层次的抽象)和输出层的神经元个数(多个线性划分)范围分别为2～4和2～6。
训练时神经网络中的各参数定义如下：
输入层单元输入向量为Pks＝(a1,a2,…,an)，目标向量Tks＝(d1,d2,…,dn)，n表示输入层单元输入向量个数；
隐含层单元输入向量Sks＝(s1,s2,…,sp)，输出向量Bks＝(b1,b2,...,bp)，p表示隐含层单元输入向量个数；
输出层单元输入向量Lks＝(l1,l2,…,lq)，输出向量Cks＝(c1,c2,…,cq)，ks＝1,2,…,m表示样本数据个数，q表示输出层单元输入向量个数；
神经网络模型的学习过程步骤如下：
a)初始化各层的连接权值和阈值，给输入层到隐含层的连接权隐含层到输出层的连接权vjt、隐含层各单元输出阈值θj、输出层各单元的输出阈值yt赋予区间(-1,1)内的随机值；is＝1,2,…,n，j＝1,2,…,p，t＝1,2,…,q；
b)选取输入样本和输出样本；
c)用输入样本、连接权、输入阈值和输出阈值计算隐含层和输出层各单元的输出:
其中sj表示隐含层单元的输入值、为输入层的输出向量、bj表示隐含层单元的输出值、lt表示输出层单元的输入值、ct表示输出层单元的输出值；
d)计算输出层各单元一般化误差然后利用隐含层到输出层的连接权vjt、隐含层的输出向量Bk＝(b1,b2,...,bp)、输出层各单元一般化误差/>计算隐含层各单元的一般化误差/>计算公式如下：
e)利用输出层各单元的一般化误差与隐含层各单元的输出值来修正隐含层到输出层的连接权vjt、输出阈值yt：
同样利用隐含层各单元的一般化误差与输入层的输入来修正输入层到隐含层的连接权wisj、输出阈值θj：
vjt(N)表示当前迭代的连接权vjt，vjt(N+1)表示下一次迭代的连接权；yt(N)表示当前迭代的输出阈值，yt(N+1)表示下一次迭代的输出阈值；
wisj(N)表示当前迭代的连接权wisj，wisj(N+1)表示下一次迭代的连接权，θj(N)表示当前迭代的阈值，θj(N+1)表示下一次迭代的阈值，N＝1,2，...，NN，其中NN表示设定的学习迭代次数；
f)选取下一个输入样本和输出样本，返回到步骤c)，直到m个训练样本训练完毕；
g)计算所有样本的累计误差E，累计误差计算方式为其中q表示输出层单元数，m表示样本数量,Et表示样本之间的误差。如果样本累计误差E小于预先设定值ε，或者当前学习迭代次数大于设定的学习迭代次数，那么学习训练结束。否则再次选取样本输入和目标输出，然后返回到步骤c)；
通过以上学习过程，获得了描述工厂实际数据输入和对机理模型的动态模拟后输出的内部参数与工厂实际的输出数据的神经网络模型。
与现有的技术相比，本发明具有如下优点：
本发明通过对离散后的机理模型动态模拟后，可以得出一些无法实际测得的反应器内部参数，如催化剂活性等，根据得到的这些动态的内部参数以及催化重整过程中反应器的实测值，运用神经网络学习的方式来建立出动态的数据模型。这样建立出来的数据与机理相混合的动态混合模型既可以解决优化求解速度较慢的问题，又可以解决数据模型无法获得内部参数的问题，两者相结合，可以既快速又准确的获得模型的输出，以此提高了催化重整过程模型的通用性与适应性。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的分析。
一种基于数据与机理相混合的催化重整过程建模方法，如图1包括以下步骤：
步骤1、依据重整动力学、热力学、物料平衡、能量平衡等基本原理，建立机理模型；
步骤2、为了使得机理模型与实际生产过程达到一致，根据催化重整过程中反应器的输入实测值，采用内点法的参数校正方法来对步骤(1)机理模型的模型参数进行优化；
步骤3、将步骤(2)优化后机理模型进行离散，得到离散模型；采用能够求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数；
步骤4、构建数据模型，训练时采用神经网络学习方式，以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，以催化重整过程中反应器的输出实测值为输出；其中催化重整过程中反应器的输入实测值为重整进料各集总组分含量、催化剂含量等，催化重整过程中反应器的输出实测值为反应器输出产物含量等。
步骤5、训练好的数据模型以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，根据其输出与催化重整过程中反应器的输出实测值相比较，若差值大于阈值，则返回步骤2，反之则结束，得到所需的优化后机理模型和数据模型相结合的混合模型。
所述的步骤1具体的实施步骤如下：
重整反应器是整个重整装置中的核心部分，机理建模最关键的部分就是针对重整反应器的建模，集总方法是反应器建模的必然选择。本发明将反应物料划分为33集总，其中集总的划分主要为烷烃、环烷烃、芳烃、烯烃、氢气等，根据其反应过程可以得出39个反应速率方程，如下所示：
烷烃脱氢反应(可逆)：
环烷烃脱氢反应(可逆)：
芳烃氢解反应：
烷烃加氢裂化反应：
其中为第jz个反应的反应速率；YP、YN、YA分别为烷烃、环烷烃、芳烃的摩尔流量；VC为催化剂装填量；F为进料体积流量；Kepj为第jz个反应平衡常数；/>为反应速率常数，其表达式为：
其中为反应频率因子；/>分别为反应活化能和压力指数；/>为/>压力指数下反应物中的氢分压；R'为摩尔气体常量；T为反应温度；/>为催化剂活性因子；/>为重整反应器装置因子。
本发明连续重整工艺采用常用的径向反应器，径向反应器物流通过反应器的压力降比轴向反应器小，有利于减小临氢系统的压力降。结合物料平衡和能量平衡原理，可得到重整反应器模型方程组：
其中为33集总组分的摩尔流量；R为径向反应器床层半径；H为反应器高度；液时空速LHSV＝F/Vc；/>为各集总组分的化学计量系数矩阵；/>为反应热向量；/>为恒压比热容向量；/>为重整反应的反应速率向量；
由于因此重整反应器机理模型方程组可以改写为：
其中为反应速率常数矩阵。
重整反应器机理模型见方程组(7)中的基础物性数据具体如下：
在重整模拟过程中，必要物性数据的获取是模型计算的前提条件，这些数据的准确性对模拟精度有非常大的影响。一些热力学性质数据是必须通过关联式由基础数据计算得到的。在重整温度条件下，所有反应组分均呈气态形式，可以将它们看作为理想气体状态，使用热力学一致性方程。
1)焓
Hm＝A+BT+CT2+DT3+ET4+F'T5 (8)
其中Hm为理想气体在T时的焓，T为反应温度；A、B，C，D，E，F'为导出系数。
2)熵
其中S为理想气体在T时的熵；G为导出系数。
3)恒压比热容
其中Cp为理想气体在恒压下的热容。
4)反应热
在温度T下，反应热的计算公式为：
其中代表了标准摩尔反应焓；α为化学计量系数；/>代表了标准摩尔生成焓；Cp(T)为各组分在温度T下的恒压比热容。
5)反应平衡常数计算
由于前15个反应是可逆反应，因此上述反应的温度T及标准压力为100Kpa下的平衡常数为：
其中为温度T下的标准摩尔生成吉布斯自由能，也称为标准摩尔自由焓，计算式为：
其中已知，标准摩尔生成熵/>的计算公式为：
其中为标准状态下的组分标准摩尔熵。
所述的步骤2具体的实施步骤如下：
2)为了使得机理模型与实际生产过程达到一致，根据实际工厂的输入数据，采用内点法的参数校正方法来对模型参数进行获取，设定一个待估参数向量X，X表示为反应的指前因子以及装置性能系数，其目标函数可以设定为：
其中G(X)表示催化重整过程中反应器的实测值与催化重整过程中反应器的设定值之间的关系，表示第in个约束方程，m'表示约束方程的总个数；
根据目标函数(15)构造惩罚项为允许误差为ε；
2)初始化X和惩罚函数因子νk，,设定迭代次数k＝1；
3)求解下列无约束问题：
min G(X)+νkβ(X) (16)
得到X极小值；
4)判断是否满足νkβ(Xk)＜ε，若是则结束，当前迭代次数下得到的X极小值即为最优解；反之令k＝k+1，νk＝C"νk-1，C"∈(0,1)，且返回步骤3)。
所述的步骤3具体的实施步骤如下：
方程组(7)属于既包含微分方程又包含代数方程的微分代数优化命题，基于有限元配置的全离散联立求解技术因其独特的优势非常合适求解该类问题。采用有限元配置方法对包含微分代数方程的机理模型方程组(7)进行离散，得到离散模型，即得到许多离散变量以及方程。采用求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数，例如催化剂吸附平衡常数、反应器装置因数等；
作为优选，所述有限元配置方法具有龙格-库塔离散方法的多种优势，在该方法中微分状态变量表达式如下：
其中w(z)表示当前有限元微分状态变量；表示在第iy个有限元初始处的值；/>是有限元iy的长度；/>表示在第iy个有限元上配置点qy处的一阶导数；/>为满足以下条件的K阶多项式：
其中为有限元中第ry个配置点的位置。qy＝1,...,K,ry＝1,...,K
微分方程的连续性方程表示为：
此外，对于过程控制变量和代数变量可采用拉格朗日形式的插值多项式表示。
代数变量形式如下：
过程控制变量形式如下：
其中和/>分别表示第iy个有限元置点qy处的代数变量和控制变量值，且满足z表示空间中的位置，/>表示第iy配置点所在的空间位置，/>为K阶拉格朗日多项式。
所述的步骤4中神经网络包括输入层、隐含层(中间层)和输出层，隐含层(对输入特征多层次的抽象)和输出层的神经元个数(多个线性划分)范围分别为2～4和2～6。
训练时神经网络中的各参数定义如下：
输入层单元输入向量为Pks＝(a1,a2,…,an)，目标向量Tks＝(d1,d2,…,dn)，n表示输入层单元输入向量个数；
隐含层单元输入向量Sks＝(s1,s2,…,sp)，输出向量Bks＝(b1,b2,...,bp)，p表示隐含层单元输入向量个数；
输出层单元输入向量Lks＝(l1,l2,…,lq)，输出向量Cks＝(c1,c2,…,cq)，ks＝1,2,…,m表示样本数据个数，q表示输出层单元输入向量个数；
神经网络模型的学习过程步骤如下：
a)初始化各层的连接权值和阈值，给输入层到隐含层的连接权隐含层到输出层的连接权vjt、隐含层各单元输出阈值θj、输出层各单元的输出阈值yt赋予区间(-1,1)内的随机值；is＝1,2,…,n，j＝1,2,…,p，t＝1,2,…,q；
b)选取输入样本和输出样本；
c)用输入样本、连接权、输入阈值和输出阈值计算隐含层和输出层各单元的输出:
其中sj表示隐含层单元的输入值、为输入层的输出向量、bj表示隐含层单元的输出值、lt表示输出层单元的输入值、ct表示输出层单元的输出值；
d)计算输出层各单元一般化误差然后利用隐含层到输出层的连接权vjt、隐含层的输出向量Bk＝(b1,b2,...,bp)、输出层各单元一般化误差/>计算隐含层各单元的一般化误差/>计算公式如下：
e)利用输出层各单元的一般化误差与隐含层各单元的输出值来修正隐含层到输出层的连接权vjt、输出阈值yt：
同样利用隐含层各单元的一般化误差与输入层的输入来修正输入层到隐含层的连接权wisj、输出阈值θj：
vjt(N)表示当前迭代的连接权vjt，vjt(N+1)表示下一次迭代的连接权；yt(N)表示当前迭代的输出阈值，yt(N+1)表示下一次迭代的输出阈值；
wisj(N)表示当前迭代的连接权wisj，wisj(N+1)表示下一次迭代的连接权，θj(N)表示当前迭代的阈值，θj(N+1)表示下一次迭代的阈值，N＝1,2，...，NN，其中NN表示设定的学习迭代次数；
f)选取下一个输入样本和输出样本，返回到步骤c)，直到m个训练样本训练完毕；
g)计算所有样本的累计误差E，累计误差计算方式为其中q表示输出层单元数，m表示样本数量,Et表示样本之间的误差。如果样本累计误差E小于预先设定值ε，或者当前学习迭代次数大于设定的学习迭代次数，那么学习训练结束。否则再次选取样本输入和目标输出，然后返回到步骤c)；
通过以上学习过程，获得了描述工厂实际数据输入和对机理模型的动态模拟后输出的内部参数与工厂实际的输出数据的神经网络模型。
以上内容是结合具体的优选实施方式对本发明所做的进一步详细说明，不能认定本发明的具体实施只限于这些说明。对于本发明所属技术领域的技术人员来说，在不脱离发明构思的前提下还可以做出一定程度的简单推演或者替换，都应当视为属于本发明的保护范围。
Claims (6)
1.基于数据与机理相混合的催化重整过程建模方法，其特征在于该方法包括以下步骤：
步骤1、建立机理模型；
步骤2、根据催化重整过程中反应器的输入实测值，采用内点法的参数校正方法来对步骤(1)机理模型的模型参数进行优化；
2-1目标函数设定为：
min G(X)
其中G(X)表示催化重整过程中反应器的实测值与催化重整过程中反应器的设定值之间的关系，表示第in个约束方程，m'表示约束方程的总个数；
根据目标函数(15)构造惩罚项为允许误差为ε；
2-2初始化X和惩罚函数因子νk，,设定迭代次数k＝1；
2-3求解下列无约束问题：
min G(X)+νkβ(X) (16)得到X极小值；
2-4判断是否满足νkβ(Xk)＜ε，若是则结束，当前迭代次数下得到的X极小值即为最优解；反之令k＝k+1，νk＝C"νk-1，C"∈(0,1)，且返回步骤2-3；
步骤3、将步骤(2)优化后机理模型进行离散，得到离散模型；采用能够求解大规模NLP问题的求解器对离散模型进行动态模拟，得到实测无法获取的反应器内部参数；
步骤4、构建数据模型，训练时采用神经网络学习方式，以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，以催化重整过程中反应器的输出实测值为输出；其中催化重整过程中反应器的输入实测值为重整进料各集总组分含量、催化剂含量，催化重整过程中反应器的输出实测值为反应器输出产物含量；
步骤5、训练好的数据模型以催化重整过程中反应器的输入实测值和动态模拟后的机理模型输出为输入，根据其输出与催化重整过程中反应器的输出实测值相比较，若差值大于阈值，则返回步骤2，反之则结束，得到所需的优化后机理模型和数据模型相结合的混合模型。
2.根据权利要求1所述的基于数据与机理相混合的催化重整过程建模方法，其特征在于步骤(1)中，将重整反应器内反应物料划分为33集总，其中集总的划分主要为烷烃、环烷烃、芳烃、烯烃、氢气，根据其反应过程可以得出39个反应速率方程，如下：
烷烃脱氢反应：
环烷烃脱氢反应：
芳烃氢解反应：
烷烃加氢裂化反应：
其中为第jz个反应的反应速率；YP、YN、YA分别为烷烃、环烷烃、芳烃的摩尔流量；VC为催化剂装填量；F为进料体积流量；Kepj为第jz个反应平衡常数；/>为反应速率常数，其表达式为：
其中为反应频率因子；/>分别为反应活化能和压力指数；/>为/>压力指数下反应物中的氢分压；R'为摩尔气体常量；T为反应温度；/>为催化剂活性因子；/>为重整反应器装置因子；
结合物料平衡和能量平衡原理，可得到重整反应器模型方程组：
其中为33集总组分的摩尔流量；R为径向反应器床层半径；H为反应器高度；液时空速LHSV＝F/Vc；/>为各集总组分的化学计量系数矩阵；/>为反应热向量；/>为恒压比热容向量；/>为重整反应的反应速率向量；
由于因此重整反应器机理模型方程组可以改写为：
其中为反应速率常数矩阵。
3.根据权利要求2所述的基于数据与机理相混合的催化重整过程建模方法，其特征在于步骤(1)中，重整反应器机理模型中的基础物性数据具体如下：
1)焓
Hm＝A+BT+CT2+DT3+ET4+F'T5 (8)
其中Hm为理想气体在T时的焓，T为反应温度；A、B，C，D，E，F'为导出系数；
2)熵
其中S为理想气体在T时的熵；G为导出系数；
3)恒压比热容
其中Cp为理想气体在恒压下的热容；
4)反应热
在温度T下，反应热的计算公式为：
其中代表了标准摩尔反应焓；α为化学计量系数；/>代表了标准摩尔生成焓；Cp(T)为各组分在温度T下的恒压比热容；
5)反应平衡常数计算
由于前15个反应是可逆反应，因此上述反应的温度T及标准压力为100Kpa下的平衡常数为：
其中为温度T下的标准摩尔生成吉布斯自由能，也称为标准摩尔自由焓，计算式为：
其中已知，标准摩尔生成熵/>的计算公式为：
其中为标准状态下的组分标准摩尔熵。
4.根据权利要求1所述的基于数据与机理相混合的催化重整过程建模方法，其特征在于步骤(3)中有限元配置方法微分状态变量表达式如下：
其中w(z)表示当前有限元微分状态变量；表示在第iy个有限元初始处的值；/>是有限元iy的长度；/>表示在第iy个有限元上配置点qy处的一阶导数；/>为满足以下条件的K阶多项式：
其中为有限元中第ry个配置点的位置；qy＝1,...,K,ry＝1,...,K
微分方程的连续性方程表示为：
此外，对于过程控制变量和代数变量可采用拉格朗日形式的插值多项式表示；
代数变量形式如下：
过程控制变量形式如下：
其中和/>分别表示第iy个有限元置点qy处的代数变量和控制变量值，且满足Z表示空间中的位置，/>表示第iy配置点所在的空间位置，/>为K阶拉格朗日多项式。
5.一种计算机可读存储介质，其上存储有计算机程序，当所述计算机程序在计算机中执行时，令计算机执行权利要求1-4中任一项所述的方法。
6.一种计算设备，包括存储器和处理器，所述存储器中存储有可执行代码，所述处理器执行所述可执行代码时，实现权利要求1-4中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110599815.6A CN113628692B (zh) | 2021-05-31 | 2021-05-31 | 基于数据与机理相混合的催化重整过程建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110599815.6A CN113628692B (zh) | 2021-05-31 | 2021-05-31 | 基于数据与机理相混合的催化重整过程建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113628692A CN113628692A (zh) | 2021-11-09 |
| CN113628692B true CN113628692B (zh) | 2024-01-30 |
Family
ID=78378073
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110599815.6A Active CN113628692B (zh) | 2021-05-31 | 2021-05-31 | 基于数据与机理相混合的催化重整过程建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113628692B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114509951B (zh) * | 2022-04-21 | 2022-08-12 | 浙江浙能航天氢能技术有限公司 | 一种基于神经网络的加氢自适应控制方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104484714A (zh) * | 2014-11-26 | 2015-04-01 | 华东理工大学 | 一种催化重整装置收率的实时预测方法 |
| CN112782979A (zh) * | 2020-12-25 | 2021-05-11 | 杭州电子科技大学 | 一种连续催化重整装置实时优化控制系统及其方法 |
| CN112857072A (zh) * | 2020-12-25 | 2021-05-28 | 杭州电子科技大学 | 一种催化重整装置中关于加热炉温度的实时控制方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190228843A1 (en) * | 2018-01-19 | 2019-07-25 | Aspen Technology, Inc. | Molecule-Based Equation Oriented Reactor Simulation System And Its Model Reduction |
-
2021
- 2021-05-31 CN CN202110599815.6A patent/CN113628692B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104484714A (zh) * | 2014-11-26 | 2015-04-01 | 华东理工大学 | 一种催化重整装置收率的实时预测方法 |
| CN112782979A (zh) * | 2020-12-25 | 2021-05-11 | 杭州电子科技大学 | 一种连续催化重整装置实时优化控制系统及其方法 |
| CN112857072A (zh) * | 2020-12-25 | 2021-05-28 | 杭州电子科技大学 | 一种催化重整装置中关于加热炉温度的实时控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113628692A (zh) | 2021-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pavão et al. | Large‐scale heat exchanger networks synthesis using simulated annealing and the novel rocket fireworks optimization | |
| Pistikopoulos et al. | Towards the integration of process design, control and scheduling: Are we getting closer? | |
| Vega et al. | Integrated design and control of chemical processes–Part I: Revision and classification | |
| CN104484714B (zh) | 一种催化重整装置收率的实时预测方法 | |
| CN104765346B (zh) | 一种炼油过程全流程建模方法 | |
| Wang et al. | A two‐layer ensemble learning framework for data‐driven soft sensor of the diesel attributes in an industrial hydrocracking process | |
| Caccavale et al. | Control and monitoring of chemical batch reactors | |
| CN111158237B (zh) | 基于神经网络的工业炉炉温多步预测控制方法 | |
| Peirelinck et al. | Domain randomization for demand response of an electric water heater | |
| CN113628692B (zh) | 基于数据与机理相混合的催化重整过程建模方法 | |
| Lin et al. | Digital-twin-based improvements to diagnosis, prognosis, strategy assessment, and discrepancy checking in a nearly autonomous management and control system | |
| CN117031950A (zh) | 深度调峰火电机组控制系统建模的方法及装置 | |
| CN103177291B (zh) | 连续搅拌反应釜的变搜索空间rna遗传算法建模方法 | |
| CN106483853A (zh) | 炼油加热炉炉膛压力的分数阶分布式动态矩阵控制方法 | |
| Sadighi et al. | Optimizing an industrial scale naphtha catalytic reforming plant using a hybrid artificial neural network and genetic algorithm technique | |
| CN104765347A (zh) | 一种渣油延迟焦化过程中收率实时预测方法 | |
| Forbes | Model structure and adjustable parameter selection for operations optimization. | |
| CN106200379A (zh) | 一种非自衡对象的分布式动态矩阵控制方法 | |
| CN108388112A (zh) | 一种批次过程二维模型预测控制方法 | |
| Li et al. | Data-driven adaptive predictive control of hydrocracking process using a covariance matrix adaption evolution strategy | |
| CN116300755A (zh) | 基于mpc的含储热的供热系统双层优化调度方法及装置 | |
| CN112857072B (zh) | 一种催化重整装置中关于加热炉温度的实时控制方法 | |
| CN115938502A (zh) | 一种基于分子级反应机理的化工产物特性预测方法及系统 | |
| CN115206448A (zh) | 一种基于ann模型的化学反应动力学计算方法 | |
| CN114187975A (zh) | 一种柴油加氢装置数据模型的优化和自动更新方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |