适用于甲醇制烯烃反应器温度压力的监测方法技术

本发明专利技术公开了适用于甲醇制烯烃反应器温度压力的监测方法，初步筛选DMTO反应器的变量；计算MIC值，并选择辅助变量；建立SVR模型；计算变量预测值与真实值之间的残差，对残差进行标准化处理，得到离线建模数据；建立PCA模型，将残差数据划分为主元空间和残差空间，在主元空间中提取变量间信息；通过PCA模型计算监测统计量及控制限；基于SVR模型对实时数据计算模型预测值，计算其与真实值的残差；对残差矩阵进行标准化，得到矩阵；对DMTO反应器进行实时监测。本发明专利技术可以提升DMTO装置反应器的实时检测精度，并可以分析出异常变化情况，进而提升了DMTO装置的稳定性和质量，具有一定的适用性。适用性。适用性。

【技术实现步骤摘要】
适用于甲醇制烯烃反应器温度压力的监测方法


[0001]本专利技术属于煤化工监测
，具体涉及一种适用于甲醇制烯烃反应器温度压力的监测方法。

技术介绍

[0002]在石油化工行业中，乙烯和丙烯的生产能力一定程度上代表了石化工业水平，其生产路线主要是石脑油蒸汽裂解，现有研究共同研发了新型技术
‑
甲醇制烯烃技术，即DMTO，开创了由煤制烯烃的新技术路径，改变了烯烃生产完全依赖于石油的状况，为甲醇合成低碳烯烃开辟了一条新途径。DMTO的成功应用对促进煤化工与石油化工协调发展，实现煤炭资源清洁高效利用，保障能源安全，缓解石油资源供应紧张局面，具有重大意义。
[0003]二甲醚的转化率会随之降低，导致乙烯和丙烯的选择性降低。另一方面，温度的升高会加剧副反应速率，增加催化剂的积碳速率。因此，反应器的温度和压力通常希望被控制在较小范围内来确保过程的长周期稳定运行，同时也可提高烯烃选择性，保证DMTO装置的产品质量。在工业生产过程中，反应器的温度、压力等关键变量也是操作人员十分关注的对象。煤化工企业对生产过程的监测主要是通过DCS监测单个变量的变化来实现，其原理是观测单个变量数值是否超过限定值，当变量数值超出限定值范围时DCS系统则会报警。然而，通常情况下温度或压力的异常波动不会触及集散控制系统的报警限，基于单变量的故障识别方法也忽略了变量之间的相关性，反应器温度、压力的异常波动前兆很有可能隐藏在其他变量数据中。因此，操作人员很难通过DCS实时或提前获知反应器温度、压力的异常变化。
[0004]随着工厂中传感器的广泛应用和数据传输技术的发展，大量的包含过程信息的历史数据被储存下来，使数据驱动过程监测方法取得了快速的发展。相比于针对单个变量数值的过程监测手段，数据驱动的过程监测方法可以从多变量的角度挖掘过程运行的内在本质特征，及时识别出生产过程中的过程偏离，为操作人员提供过程评价和如何应对过程偏离提供参考。目前，基于数据驱动的过程监测算法开发的智能监测系统已在炼油和石化等行业成功应用。但是目前DMTO装置的温度压力的监测存在的检测准确度和及时性不足，以至于DMTO装置的运行和产品质量有待进一步提升。

技术实现思路

[0005]本专利技术的目的是提供一种适用于甲醇制烯烃反应器温度压力的监测方法，解决了目前DMTO装置的温度压力的监测存在的检测准确度和及时性不足，以至于DMTO装置的运行和产品质量有待进一步提升的问题。
[0006]本专利技术所采用的技术方案是：
[0007]适用于甲醇制烯烃反应器温度压力的监测方法，S1：初步筛选DMTO反应器的变量；S2，计算与反应器温度压力相关变量的MIC值，并选择辅助变量；S3，基于辅助变量建立SVR模型；S4，基于SVR模型计算每个变量预测值与真实值之间的残差，并对残差进行标准化处理，得到离线建模数据；S5，利用残差建立PCA模型，将残差数据划分为主元空间和残差空
间，在主元空间中提取用于数据降维的变量间信息；S6，通过PCA模型计算监测统计量及控制限；S7，基于SVR模型对采集的实时数据计算模型预测值，计算其与真实值的残差；S8，使用离线建模数据对残差矩阵进行标准化，得到标准化后的矩阵；S9，对DMTO反应器进行实时监测。
[0008]本专利技术的特点还在于：
[0009]S9具体为：
[0010]通过离线建模时得到的载荷矩阵将标准化后的残差映射到离线建模时建立的PCA模型的主元空间中，并计算残差在主元空间上的监测统计量，通过与离线建模时的监测控制限进行对比，当连续三个点的监测统计量大于控制限时，则判定系统出现故障，系统发出报警并记录故障数据；反之则系统处于正常状态，继续对下一个时刻的数据进行监测。
[0011]S1具体为：
[0012]根据具体工艺流程初步筛选DMTO反应器温度压力相关变量作为监测变量。
[0013]S2具体为：
[0014]计算与反应器温度压力相关变量的MIC值，根据MIC值选择合适的辅助变量，然后采集生产过程一个完整周期的正常工况数据组成离线数据集X∈R
n
×
m
,n为数据集样本个数，m为监测变量个数，MIC计算方法下所示：
[0015]首先计算两个变量间的互信息：
[0016][0017]其中，I(X；Y)为变量X、Y的互信息，p
(x)
和p
(y)
是各自变量的边际概率密度函数，p
(x,y)
为变量X与变量Y之间的联合概率密度函数，记录其中最大的互信息值为I(D,x,y)；为了比较不同维度之间的值，需要对互信息值进行归一化处理，归一化后的值区间在0
‑
1之间，归一化及最大互信息值的计算公式如下：
[0018][0019][0020]其中，a和b是在x轴和y轴上划分格子的个数，B
(n)
为网格划分的上限个数，一般取值B
(n)
＝n
0.6
，n为数据规模。
[0021]S3具体为：
[0022]SVR支持向量回归的原理如下：
[0023]M＝{(x1,y1),(x2,y2),......,(x
m
,y
m
)}
ꢀꢀ
(4)；
[0024]对于给定训练集样本M，希望得到一个回归模型使得f(x)与y尽可能接近，w,b是待确定的回归参数，如果能够忍受f(x)与y之间最大偏差为ε，支持向量回归的求解问题就可以变为如式(5)所示形式：
[0025][0026]其中C是正则化参数，l
e
是不敏感损失函数如式(6)所示：
[0027][0028]引入松弛变量ξ
i
和后，式(5)可重写如式(7)所示：
[0029][0030]通过拉格朗日乘子μ
i
的引入，可以将原问题转变为对应的对偶问题来求解：
[0031][0032]令L对w,b,ξ
i
,的偏导数为0，即可对w,b进行确定，得到最终的支持向量回归解形式如式(9)所示：
[0033][0034]引入核函数，从而能够完美解决数据的非线性问题，又避免了计算的复杂度过高；核函数的表达式如公式(10)所示。
[0035]k(x
i
,x
j
)＝Φ(x
i
)
T
Φ(x
j
)
ꢀꢀ
(10)；
[0036]其中Φ()表示非线性变换，用于将输入空间映射到高维特征空间中。
[0037]S5中，PCA具体计算步骤如下所示：
[0038]记标准化后得到的数据集为X，利用协方差矩阵计算各维度之间的相关性系数矩阵Cov(X)，协方差的公式为：
[0039][0040]然后对相关性系数矩阵Cov(YX)进行奇异值分解，可以得到其特征值和特征向量，将特征向量排成一个矩阵本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.适用于甲醇制烯烃反应器温度压力的监测方法，其特征在于，S1：初步筛选DMTO反应器的变量；S2，计算与反应器温度压力相关变量的MIC值，并选择辅助变量；S3，基于辅助变量建立SVR模型；S4，基于SVR模型计算每个变量预测值与真实值之间的残差，并对残差进行标准化处理，得到离线建模数据；S5，利用残差建立PCA模型，将残差数据划分为主元空间和残差空间，在主元空间中提取用于数据降维的变量间信息；S6，通过PCA模型计算监测统计量及控制限；S7，基于SVR模型对采集的实时数据计算模型预测值，计算其与真实值的残差；S8，使用离线建模数据对残差矩阵进行标准化，得到标准化后的矩阵；S9，对DMTO反应器进行实时监测。2.根据权利要求1所述的适用于甲醇制烯烃反应器温度压力的监测方法，其特征在于，步骤9具体为：通过离线建模时得到的载荷矩阵将标准化后的残差映射到离线建模时建立的PCA模型的主元空间中，并计算残差在主元空间上的监测统计量，通过与离线建模时的监测控制限进行对比，当连续三个点的监测统计量大于控制限时，则判定系统出现故障，系统发出报警并记录故障数据；反之则系统处于正常状态，继续对下一个时刻的数据进行监测。3.根据权利要求1所述的适用于甲醇制烯烃反应器温度压力的监测方法，其特征在于，步骤1具体为：根据具体工艺流程初步筛选DMTO反应器温度压力相关变量作为监测变量。4.根据权利要求1所述的适用于甲醇制烯烃反应器温度压力的监测方法，其特征在于，步骤2具体为：计算与反应器温度压力相关变量的MIC值，根据MIC值选择辅助变量，然后采集生产过程一个完整周期的正常工况数据组成离线数据集X∈R
n
×
m
,n为数据集样本个数，m为监测变量个数，MIC计算方法下所示：首先计算两个变量间的互信息：其中，I(X；Y)为变量X、Y的互信息，p
(x)
.和p
(y)
是各自变量的边际概率密度函数，p
(x,y)
为变量X与变量Y之间的联合概率密度函数，记录其中最大的互信息值为I(D,x,y)；为了比较不同维度之间的值，需要对互信息值进行归一化处理，归一化后的值区间在0
‑
1之间，归一化及最大互信息值的计算公式如下：(2)；其中，a和b是在x轴和y轴上划分格子的个数，B
(n)
为网格划分的上限个数，一般取值B
(n)
＝n
0.6
，n为数据规模。5.根据权利要求1所述的适用于甲醇制烯烃反应器温度压力的监测方法，其特征在于，步骤3具体为：SVR支持向量回归的原理如下：M＝{(x1,y1),(x2,y2),......,(x
m
,y
m
)}
ꢀꢀ
(4)；
对于给定训练集样本M，希望得到一个回归模型使得f(x)与y尽可能接近，w,b是待确定的回归参数，如果能够忍受f(x)与y之间...

【专利技术属性】
技术研发人员：赵泽盟，史元腾，李超，锁良晨，吴晓茹，
申请(专利权)人：中煤能源研究院有限责任公司，
类型：发明
国别省市：