CN106769646B - 一种基于反演思想的精馏过程双变量扰动原因诊断方法 - Google Patents

一种基于反演思想的精馏过程双变量扰动原因诊断方法 Download PDF

Info

Publication number
CN106769646B
CN106769646B CN201610993039.7A CN201610993039A CN106769646B CN 106769646 B CN106769646 B CN 106769646B CN 201610993039 A CN201610993039 A CN 201610993039A CN 106769646 B CN106769646 B CN 106769646B
Authority
CN
China
Prior art keywords
disturbance
inversion
cause
bivariate
rectification process
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
Application number
CN201610993039.7A
Other languages
English (en)
Other versions
CN106769646A (zh
Inventor
田文德
梁慧婷
崔哲
李乐宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao University of Science and Technology
Original Assignee
Qingdao University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qingdao University of Science and Technology filed Critical Qingdao University of Science and Technology
Priority to CN201610993039.7A priority Critical patent/CN106769646B/zh
Publication of CN106769646A publication Critical patent/CN106769646A/zh
Application granted granted Critical
Publication of CN106769646B publication Critical patent/CN106769646B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于反演思想的精馏过程双变量扰动原因诊断方法。本发明通过双变量扰动响应曲线之间分布的规律，提取特征量，通过特征量建立与扰动原因之间的数学模型，运算实现双变量扰动的定量辨识，建立的反演模型可以深入分析扰动原因，提高了扰动诊断的准确性，并且在只有扰动数据的情况下可以快速、准确地得出扰动量的大小。

Description

一种基于反演思想的精馏过程双变量扰动原因诊断方法
技术领域
本发明涉及一种精馏过程故障诊断技术领域，特别涉及一种基于反演思想的精馏过程双变量扰动原因诊断方法。
背景技术
精馏是石油化工生产过程中应用最广泛的操作之一，对精馏设备和操作所作的微小的扰动会产生巨大的经济损失。为确保精馏装置安全、稳定运行，及时发现异常状态中参数的裂化趋势，需要及时识别异常原因，从源头上预防和控制精馏事故的发生。同时精馏过程中存在大量扰动，使得扰动原因难以确定，扰动量信息诊断精度难以提高。
对于双变量的故障诊断，普遍使用的是多变量统计的方法。如主元分析(PCA)、偏最小二乘(PLS)算法等，以及这些算法的改进算法，例如动态PCA、核PCA及独立成分分析(ICA)与动态时间规整(DTW)、累积和控制图(CUSUM)或支持向量机(SVM)等的结合。以上方法只是在故障发生之后诊断出故障的位置或是实时侦测故障是否发生，不能深入分析扰动原因。本文将反演的思想应用于精馏故障诊断领域，在只有扰动数据或曲线的情况下，使用此方法能快速、准确地实现扰动量的定量分析。
发明内容
针对上述现有技术中存在的问题，本发明提供一种基于反演思想的精馏过程扰动原因诊断的方法，可深入分析扰动原因，快速、准确地实现扰动量的定量辨识。
为实现上述目的，本发明提供如下技术方案：
一种基于反演思想的精馏过程双变量扰动原因诊断的方法，包括扰动原因反演模型的建立，其具体实施步骤如下：
(1)基于精馏过程机理，动态模拟出双变量同时扰动时的正常及异常样本；
(2)绘制出塔顶组分(苯浓度)随进料量、进料组成及回流比扰动下的响应曲线；
(3)根据扰动曲线的分布规律提取特征量，建立扰动特征量与扰动原因的反演模型。
进一步，步骤(3)中建立扰动特征量与扰动原因的反演模型具体如下：
根据模拟出的扰动响应曲线的分布规律，找出两个能特征表示曲线的量(曲线在某点的切线斜率和扰动曲线和正常曲线之间的欧式距离)，根据遗传算法搜索最优函数形式，将反演问题转化为函数优化问题，得到反演模型。
本发明的有益效果如下：
(1)本发明提供了一种基于反演思想的精馏过程双变量扰动原因诊断方法，实现了由扰动数据直接反向推理出扰动量的目的，更深入的分析了扰动原因。
附图说明
图1为进料量和进料组成扰动时苯浓度的响应曲线图；
图2为进料量和回流比扰动时苯浓度的响应曲线图；
图3为进料组成和回流比扰动时苯浓度的响应曲线图；
具体实施方式
为了使本领域的人员更好地理解本发明的技术方案，下面对本发明的技术方案进行清楚、完整的描述，基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的其它类同实施例，都应当属于本申请保护的范围。
实施例一：
一种基于反演思想的精馏过程双变量扰动原因诊断方法，包括扰动原因反演模型的建立，其步骤具体如下：
(1)基于精馏过程动态机理，基于MATLAB平台动态模拟出双变量同时扰动时的正常及异常样本；
(2)绘制出塔顶轻组分(苯浓度)随进料量、进料组成及回流比扰动下的响应曲线；
(3)根据扰动曲线的分布规律提取特征量，建立扰动特征量与扰动原因的反演模型。
步骤(2)中塔顶轻组分(苯浓度)随进料量、进料组成及回流比扰动下的响应曲线具体如下：
设置试验时间为40h，采样间隔为5s，异常类型设定为进料流量与进料组成、进料流量与回流比、进料组成与回流比双变量同时扰动，扰动幅度设置为±1～10％以内。
步骤(3)中建立扰动特征量与扰动原因的反演模型具体如下：
如图1、2、3所示，可以看出，六种扰动类型对苯的影响大体趋势相同。随着扰动量的增加扰动曲线排布在标准曲线的上方，最后的浓度越来越大并趋近于1，随着扰动量的减少扰动曲线出现在标准曲线的下方，且呈现一定的规律，根据模拟出的扰动响应曲线的分布规律，找出两个能特征表示曲线的量(曲线在1h处的切线斜率和扰动曲线和正常曲线之间的欧式距离)，根据遗传算法搜索最优函数形式，将反演问题转化为函数优化问题，得到反演模型。具体反演模型如式(1)～(12):
Figure BDA0001150642890000041
Figure BDA0001150642890000042
Figure BDA0001150642890000043
ΔzF2＝3.598-0.00291*d2+0.000175*d2 2-1.082E-5*d2 3+2.748E-7*d2 4
-2.56E-9*d2 5-46.6*k2+150.035*k2 2 (4)
ΔF3＝0.79+0.938*d3-1.0225*d3 2+0.268*d3 3+206.55*ln(x3)+10298.64*(ln(x3))2 (5)
Figure BDA0001150642890000044
Figure BDA0001150642890000045
Figure BDA0001150642890000046
Figure BDA0001150642890000047
Figure BDA0001150642890000048
Figure BDA0001150642890000049
Figure BDA00011506428900000410
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims (2)

1.一种基于反演思想的精馏过程双变量扰动原因诊断方法，其特征在于，包括如下步骤：
1)基于精馏机理动态模拟扰动正常、扰动样本，绘制出扰动响应曲线；
2)依据扰动响应曲线的规律性，选择各扰动曲线与标准曲线之间的欧氏距离和扰动曲线在某点的切线斜率作为特征量定量表示各扰动；
3)基于遗传算法自搜索特征量与扰动量之间的函数关系式，建立扰动原因反演模型，实现定量分析扰动原因的目的。
2.根据权利要求1所述的一种基于反演思想的精馏过程双变量扰动原因诊断方法，其特征在于，步骤3)中所述扰动原因反演模型基于以下方程实现扰动原因的定量分析：
Figure FDA0002268608770000011
Figure FDA0002268608770000012
Figure FDA0002268608770000013
Figure FDA0002268608770000014
Figure FDA0002268608770000015
Figure FDA0002268608770000016
Figure FDA0002268608770000017
Figure FDA0002268608770000018
Figure FDA0002268608770000019
Figure FDA0002268608770000022
Figure FDA0002268608770000023
其中，d为欧式距离，k为切线斜率，ΔF，ΔzF，ΔR分别为进料量、进料组成和回流比扰动量。
CN201610993039.7A 2016-11-11 2016-11-11 一种基于反演思想的精馏过程双变量扰动原因诊断方法 Active CN106769646B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610993039.7A CN106769646B (zh) 2016-11-11 2016-11-11 一种基于反演思想的精馏过程双变量扰动原因诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610993039.7A CN106769646B (zh) 2016-11-11 2016-11-11 一种基于反演思想的精馏过程双变量扰动原因诊断方法

Publications (2)

Publication Number Publication Date
CN106769646A CN106769646A (zh) 2017-05-31
CN106769646B true CN106769646B (zh) 2020-03-03

Family

ID=58973400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610993039.7A Active CN106769646B (zh) 2016-11-11 2016-11-11 一种基于反演思想的精馏过程双变量扰动原因诊断方法

Country Status (1)

Country Link
CN (1) CN106769646B (zh)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102339040B (zh) * 2010-07-15 2013-06-19 清华大学 精馏塔优化控制方法
CN102522759B (zh) * 2011-12-13 2013-11-06 南方电网科学研究院有限责任公司 一种基于扰动后电压轨迹的电力系统暂态失稳实时判别方法
CN102760208B (zh) * 2012-07-03 2015-07-29 清华大学 基于模拟疫苗的动态人工免疫故障诊断方法
CN102903268A (zh) * 2012-08-23 2013-01-30 上海国际化建工程咨询公司 合成气化工及石油化工生产过程动态模拟仿真培训系统
CN103048927B (zh) * 2012-12-28 2015-03-25 浙江大学 用于精馏系统的模型预测控制方法
CN103488085A (zh) * 2013-09-22 2014-01-01 上海交通大学 甲醇四塔精馏系统多目标优化控制方法
CN103678903A (zh) * 2013-12-06 2014-03-26 杭州电子科技大学 一种动态精馏塔建模方法
CN104731057B (zh) * 2015-01-28 2017-10-20 河北工业大学 一种基于有效相对增益矩阵法的萃取精馏工艺控制方法

Also Published As

Publication number Publication date
CN106769646A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106982196B (zh) 一种异常访问检测方法及设备
CN106770015B (zh) 一种基于主成分分析相似判别的油品性质检测方法
WO2019050966A3 (en) AUTOMATED SAMPLE WORKFLOW PORTALING AND DATA ANALYSIS
CN112129741B (zh) 绝缘油老化分析方法、装置、计算机设备和存储介质
CN104990893B (zh) 一种基于相似判别的汽油辛烷值检测方法
CN106649755B (zh) 一种多维度实时变电设备数据的阈值自适应设置异常检测方法
CN105608334B (zh) 一种评价烟用辅材参数影响程度的方法
CN104316660A (zh) 快速检测豆瓣酱风味、等级及产品溯源的方法
CN103488135A (zh) 一种用于半导体生产加工过程监控的统计过程控制方法
CN108844941B (zh) 一种基于拉曼光谱和pca-hca的不同品位磷矿的鉴别和分类方法
CN107403196A (zh) 基于谱聚类分析的即时学习建模方法
US20190228122A1 (en) Method of fast identifying the distribution rule of wind speed
MX2021002598A (es) Metodo y dispositivo electronico para monitorear la fabricacion de un producto de metal, programa informatico relacionado e instalacion.
US6973390B2 (en) Method and system for analyzing wafer yield against uses of a semiconductor tool
CN106769646B (zh) 一种基于反演思想的精馏过程双变量扰动原因诊断方法
Jeong Semiconductor wafer defect classification using support vector machine with weighted dynamic time warping kernel function
CN108760268B (zh) 一种基于信息熵的立磨运行数据阶跃故障诊断方法
CN110298767B (zh) 一种火力发电厂时间序列变量异常监测方法及系统
CN109472038A (zh) 一种钢铁冶炼工艺参数曲线形态匹配的方法及系统
MY193844A (en) Plant operation support system and plant operation support method
Luo et al. A prediction method based on improved ridge regression
CN104078381B (zh) 一种量测机台监测图规格界限设定的方法
Yang et al. A simple approach for monitoring business service time variation
CN109783698A (zh) 基于Merkle-tree的工业生产数据实体识别方法
CN105447511A (zh) 一种基于Adaboost Haar-Like特征的SVM目标检测方法

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