CN109144035A - 一种基于支持向量的化工过程监测方法 - Google Patents
一种基于支持向量的化工过程监测方法 Download PDFInfo
- Publication number
- CN109144035A CN109144035A CN201811135141.9A CN201811135141A CN109144035A CN 109144035 A CN109144035 A CN 109144035A CN 201811135141 A CN201811135141 A CN 201811135141A CN 109144035 A CN109144035 A CN 109144035A
- Authority
- CN
- China
- Prior art keywords
- variable
- indicate
- data
- indicates
- monitoring
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0224—Process history based detection method, e.g. whereby history implies the availability of large amounts of data
- G05B23/024—Quantitative history assessment, e.g. mathematical relationships between available data; Functions therefor; Principal component analysis [PCA]; Partial least square [PLS]; Statistical classifiers, e.g. Bayesian networks, linear regression or correlation analysis; Neural networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于支持向量的化工过程监测方法,包括如下步骤:步骤1:收集化工过程中传感器数据,进行整理并执行偏最小二乘算法,计算得到相应的参数;步骤2:使用步骤1中得到的支持向量对化工过程进行监测,判断是否出现故障。本发明提出了一种新型的基于支持向量的化工过程监测方法。该方法通过采集化工过程中的数据建立过程的预测模型,然后基于支持向量创建一种新的过程监测指标。改进了传统过程监测指标不能很好地监测非线性过程的缺点,改善了监测的灵敏性。
Description
技术领域
本发明属于自动化技术领域,涉及一种基于支持向量的化工过程监测方法。
背景技术
作为现代工业过程的重要组成部分,过程监控发挥着重要作用。已经出现了大量的过程监测方法来检测故障。作为最常用的多变量过程监测方法之一,偏最小二乘已被广泛用于过程监测。传统的偏最小二乘方法是假设过程数据符合高斯分布,并且这些变量之间的关系是线性的,这并不符合实际情况。然而偏最小二乘仍然是一种有效且简单的数据提取方法,因此,应当构建新的统计监测指标来充分发挥偏最小二乘的优势。
发明内容
本发明的目的是针对传统的过程监测技术不能很好的监测非高斯过程数据的不足之处,提出了一种新型的基于支持向量的化工过程监测方法。该方法将原始数据投影到更高维的特征空间来解决非线性问题,引入拉格朗日方法来减小计算复杂度。基于支持向量创建新的过程监测指标,加强了过程监测系统的灵敏度,提高了故障诊断率。其具体技术方案如下:
步骤1:收集化工过程中传感器数据,进行整理并执行偏最小二乘算法,计算得到相应的参数。具体步骤是:
1.1收集化工过程中产生的数据,数据分为两类,过程数据X和质量数据Y,包含N个样本。
X=[x1,x2...xm],x1,x2...xm∈RN×1
Y=[y1,y2...yp],x1,x2...xp∈RN×1
其中,x1,x2...xm等分别表示化工过程中的反应物浓度,压力,温度…阀门开度等变量,y1,y2...yp等分别表示产物浓度,产物纯度…产物温度等与质量相关的变量。m表示过程变量的个数,p表示质量变量的个数。
1.2由于过程中采集到的数据可能存在单位不一致,数值相差过大等情况,因此对步骤1.1中采集的数据进行标准化。
其中,表示过程数据的均值,表示质量数据的均值。表示标准化后的过程数据,表示标准化后的质量数据。
1.3计算标准化过程变量的潜在变量t和标准化质量变量的潜在变量u,首先初始化u为Y的任意一列,
其中,w表示标准化过程变量与潜在变量u对应的权重向量,u表示标准化质量变量对应的潜在变量,‘表示转置符号。
1.4重复步骤1.3,直到潜在变量t收敛,然后缩减标准化过程变量和标准化质量变量
1.5根据步骤1.3和1.4,计算需要的参数。
B=W(P'W)-1C'
其中,W表示权重矩阵,P和C表示负载向量,B表示回归系数矩阵。T表示标准化过程变量的支持向量。U表示u的聚合矩阵。
1.6经过步骤1.2-步骤1.5,标准化过程变量和标准化质量变量分别可以被分解为:
其中,E表示标准化过程变量的残差矩阵,F表示标准化质量变量的残差矩阵。
1.7经过以上步骤,可以得到标准化质量变量的预测模型。
其中,表示质量变量Y的预测值。
步骤2:使用步骤一中得到的支持向量对化工过程进行监测,判断是否出现故障。具体步骤为:
2.1将每个样本映射到特征空间中的更高维的超球面,超球面的中心a和半径R,则超球体的最小体积可以得到:
其中,yi表示第i个样本映射到高维特征空间后的值,V表示超球面体积与误差之间的折衷参数。ζi表示松弛变量。
2.2上面的问题可以使用另一种更容易计算的方法,即引入拉格朗日乘子的方式。
其中,xi,xj分别表示第i,j个样本的过程变量;K(xi,xj)表示第i,j个样本对应的核函数。αi,αj均表示拉格朗日乘子。
2.3计算超球面的半径。
其中,xs表示任意的一个支持向量。
2.4计算新的数据样本到中心的距离,构建新的监测统计量。
其中,xnew表示新的样本数据。D表示xnew到超球面中心的距离。
有益效果
本发明提出了一种新型的基于支持向量的化工过程监测方法。该方法通过采集化工过程中的数据建立过程的预测模型,然后基于支持向量创建一种新的过程监测指标。改进了传统过程监测指标不能很好地监测非线性过程的缺点,改善了监测的灵敏性。
具体实施方式
下面对本发明作进一步说明。
以田纳西-伊斯曼过程为例:
田纳西-伊斯曼过程由反应器,冷凝器,分离器,压缩机和汽提塔五部分组成,包含12个操作变量和41个测量变量。有21个故障变量。
步骤1:收集化工过程中传感器数据,进行整理并执行偏最小二乘算法,计算得到相应的参数。具体步骤是:
1.1收集化工过程中产生的数据,数据分为两类,过程数据X和质量数据Y,包含N个样本。
X=[x1,x2...xm],x1,x2...xm∈RN×1
Y=[y1,y2...yp],x1,x2...xp∈RN×1
其中,x1,x2...xm等分别表示化工过程中的反应物浓度,压力,温度…阀门开度等变量,y1,y2...yp等分别表示产物浓度,产物纯度…产物温度等与质量相关的变量。
1.2由于过程中采集到的数据可能存在单位不一致,数值相差过大等情况,因此对步骤1.1中采集的数据进行标准化。
其中,表示过程数据的均值,表示质量数据的均值。表示标准化后的过程数据,表示标准化后的质量数据。
1.3计算标准化过程变量的潜在变量t和标准化质量变量的潜在变量u,首先初始化u为Y的任意一列,
其中,w表示标准化过程变量与潜在变量u对应的权重向量,u表示标准化质量变量与潜在变量t对应的权重向量,‘表示转置符号。
1.4重复步骤1.3,直到潜在变量t收敛,然后缩减标准化过程变量和标准化质量变量
1.5根据步骤1.3和1.4,计算需要的参数。
B=W(P'W)-1C'
其中,W表示权重矩阵,P和C表示负载向量,B表示回归系数矩阵。T表示标准化过程变量的支持向量。
1.6经过步骤1.2-步骤1.5,标准化过程变量和标准化质量变量分别可以被分解为:
其中,E表示标准化过程变量的残差矩阵,F表示标准化质量变量的残差矩阵。
1.7经过以上步骤,可以得到标准化质量变量的预测模型。
其中,表示质量变量Y的预测值。
步骤二:使用步骤一中得到的支持向量对化工过程进行监测,判断是否出现故障。具体步骤为:
2.1将每个样本映射到特征空间中的更高维的超球面,超球面的中心a和半径R,则超球体的最小体积可以得到:
其中,yi表示第i个样本映射到高维特征空间后的值,C表示超球面体积与误差之间的折衷参数。ζi表示松弛变量。
2.2上面的问题可以使用另一种更容易计算的方法,即引入拉格朗日乘子的方式。
其中,xi,xj分别表示第i,j个样本的过程变量;K(xi,xj)表示第i,j个样本对应的核函数。αi,αj均表示拉格朗日乘子。
2.3计算超球面的半径。
其中,xs表示任意的一个支持向量。
2.4计算新的数据样本到中心的距离,构建新的监测统计量。
其中,xnew表示新的样本数据。D表示xnew到超球面中心的距离。
Claims (3)
1.一种基于支持向量的化工过程监测方法,包括如下步骤:
步骤1:收集化工过程中传感器数据,进行整理并执行偏最小二乘算法,计算得到相应的参数;
步骤2:使用步骤1中得到的支持向量对化工过程进行监测,判断是否出现故障。
2.如权利要求1所述基于支持向量的化工过程监测方法,其特征在于:
所述步骤1具体如下:
1.1收集化工过程中产生的数据,数据分为两类,过程数据X和质量数据Y,包含N个样本。
X=[x1,x2...xm],x1,x2...xm∈RN×1
Y=[y1,y2...yp],x1,x2...xp∈RN×1
其中,x1,x2...xm等分别表示化工过程中的反应物浓度,压力,温度…阀门开度等变量,y1,y2...yp等分别表示产物浓度,产物纯度…产物温度等与质量相关的变量。m表示过程变量的个数,p表示质量变量的个数;
1.2对步骤1.1中采集的数据进行标准化:
其中,表示过程数据的均值,表示质量数据的均值。表示标准化后的过程数据,表示标准化后的质量数据;
1.3计算标准化过程变量的潜在变量t和标准化质量变量的潜在变量u,首先初始化u为Y的任意一列,
其中,w表示标准化过程变量与潜在变量u对应的权重向量,u表示标准化质量变量对应的潜在变量,‘表示转置符号;
1.4重复步骤1.3,直到潜在变量t收敛,然后缩减标准化过程变量和标准化质量变量
1.5根据步骤1.3和1.4,计算需要的参数:
B=W(P'W)-1C'
其中,W表示权重矩阵,P和C表示负载向量,B表示回归系数矩阵。T表示标准化过程变量的支持向量。U表示u的聚合矩阵;
1.6经过步骤1.2-步骤1.5,标准化过程变量和标准化质量变量分别可以被分解为:
其中,E表示标准化过程变量的残差矩阵,F表示标准化质量变量的残差矩阵;
1.7经过以上步骤,可以得到标准化质量变量的预测模型:
其中,表示质量变量Y的预测值。
3.如权利要求2所述基于支持向量的化工过程监测方法,其特征在于:
所述步骤2具体为:
2.1将每个样本映射到特征空间中的更高维的超球面,超球面的中心a和半径R,则超球体的最小体积可以得到:
其中,yi表示第i个样本映射到高维特征空间后的值,V表示超球面体积与误差之间的折衷参数,ζi表示松弛变量;
2.2引入拉格朗日乘子的方式:
其中,xi,xj分别表示第i,j个样本的过程变量;K(xi,xj)表示第i,j个样本对应的核函数。αi,αj均表示拉格朗日乘子;
2.3计算超球面的半径:
其中,xs表示任意的一个支持向量;
2.4计算新的数据样本到中心的距离,构建新的监测统计量:
其中,xnew表示新的样本数据,D表示xnew到超球面中心的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811135141.9A CN109144035A (zh) | 2018-09-27 | 2018-09-27 | 一种基于支持向量的化工过程监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811135141.9A CN109144035A (zh) | 2018-09-27 | 2018-09-27 | 一种基于支持向量的化工过程监测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109144035A true CN109144035A (zh) | 2019-01-04 |
Family
ID=64813157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811135141.9A Pending CN109144035A (zh) | 2018-09-27 | 2018-09-27 | 一种基于支持向量的化工过程监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109144035A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111679651A (zh) * | 2020-06-08 | 2020-09-18 | 中国人民解放军火箭军工程大学 | 用于故障引起的变结构变参数系统的辨识方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101446831A (zh) * | 2008-12-30 | 2009-06-03 | 东北大学 | 一种分散的过程监测方法 |
CN101458522A (zh) * | 2009-01-08 | 2009-06-17 | 浙江大学 | 基于主元分析和支持向量数据描述的多工况过程监控方法 |
CN104504288A (zh) * | 2015-01-12 | 2015-04-08 | 江南大学 | 基于多向支持向量聚类的非线性多阶段间歇过程软测量方法 |
CN108062566A (zh) * | 2017-12-13 | 2018-05-22 | 北京工业大学 | 一种基于多核潜在特征提取的智能集成软测量方法 |
-
2018
- 2018-09-27 CN CN201811135141.9A patent/CN109144035A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101446831A (zh) * | 2008-12-30 | 2009-06-03 | 东北大学 | 一种分散的过程监测方法 |
CN101458522A (zh) * | 2009-01-08 | 2009-06-17 | 浙江大学 | 基于主元分析和支持向量数据描述的多工况过程监控方法 |
CN104504288A (zh) * | 2015-01-12 | 2015-04-08 | 江南大学 | 基于多向支持向量聚类的非线性多阶段间歇过程软测量方法 |
CN108062566A (zh) * | 2017-12-13 | 2018-05-22 | 北京工业大学 | 一种基于多核潜在特征提取的智能集成软测量方法 |
Non-Patent Citations (4)
Title |
---|
张铃君: "基于改进偏最小二乘的过程监测方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)》 * |
沈非凡: "基于轨迹分析的工业过程建模与故障检测", 《中国博士学位论文全文数据库 信息科技辑》 * |
王建平: "基于高阶偏最小二乘的间歇过程监控方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)》 * |
衷路生 等: "基于得分矩阵和支持向量数据描述(SM-SVDD)的过程监测", 《科学技术与工程》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111679651A (zh) * | 2020-06-08 | 2020-09-18 | 中国人民解放军火箭军工程大学 | 用于故障引起的变结构变参数系统的辨识方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107632592B (zh) | 基于高效递推核主元分析的非线性时变过程故障监测方法 | |
CN104699077B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 | |
CN108388234B (zh) | 一种基于相关性划分多变量块pca模型的故障监测方法 | |
CN109407652B (zh) | 基于主辅pca模型的多变量工业过程故障检测方法 | |
CN101169623A (zh) | 基于核主元分析贡献图的非线性过程故障辨识方法 | |
CN105404280A (zh) | 基于自回归动态隐变量模型的工业过程故障检测方法 | |
CN108375965B (zh) | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 | |
CN107403196B (zh) | 基于谱聚类分析的即时学习建模对脱丁烷塔底丁烷浓度进行预测的方法 | |
CN103389701B (zh) | 基于分布式数据模型的厂级过程故障检测与诊断方法 | |
CN104267668A (zh) | 基于贝叶斯方法的航天阀门零件加工过程故障诊断方法 | |
CN104536439B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障诊断方法 | |
CN106094749B (zh) | 基于改进核熵成分分析的非线性故障检测方法及应用 | |
CN108508865A (zh) | 一种基于分散式osc-pls回归模型的故障检测方法 | |
CN108345284A (zh) | 一种基于两变量块的质量相关故障检测方法 | |
CN115096627B (zh) | 一种液压成形智能装备制造过程故障诊断与运维方法及系统 | |
CN111397902A (zh) | 一种基于特征对齐卷积神经网络的滚动轴承故障诊断方法 | |
CN114757269A (zh) | 一种基于局部子空间-邻域保持嵌入的复杂过程精细化故障检测方法 | |
CN109144035A (zh) | 一种基于支持向量的化工过程监测方法 | |
CN108762242A (zh) | 一种基于多块典型相关分析模型的分布式故障检测方法 | |
CN111639304A (zh) | 基于Xgboost回归模型的CSTR故障定位方法 | |
CN113222046B (zh) | 基于过滤策略的特征对齐自编码器故障分类方法 | |
CN114611067A (zh) | 一种基于典型变量相异性分析的化工过程缓变故障检测方法 | |
CN112149054B (zh) | 基于时序扩展的正交邻域保持嵌入模型的构建与应用 | |
CN112748331A (zh) | 一种基于ds证据融合的断路器机械故障识别方法及装置 | |
CN116224887A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190104 |