CN100461043C - 工业聚丙烯生产熔融指数检测故障诊断***及方法 - Google Patents
工业聚丙烯生产熔融指数检测故障诊断***及方法 Download PDFInfo
- Publication number
- CN100461043C CN100461043C CNB2006101554147A CN200610155414A CN100461043C CN 100461043 C CN100461043 C CN 100461043C CN B2006101554147 A CNB2006101554147 A CN B2006101554147A CN 200610155414 A CN200610155414 A CN 200610155414A CN 100461043 C CN100461043 C CN 100461043C
- Authority
- CN
- China
- Prior art keywords
- module
- data
- training
- calculate
- sigma
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 84
- 239000004743 Polypropylene Substances 0.000 title claims abstract description 33
- -1 polypropylene Polymers 0.000 title claims abstract description 33
- 229920001155 polypropylene Polymers 0.000 title claims abstract description 33
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 32
- 238000003745 diagnosis Methods 0.000 title claims abstract description 30
- 238000001514 detection method Methods 0.000 title claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 52
- 238000004458 analytical method Methods 0.000 claims abstract description 50
- 238000012706 support-vector machine Methods 0.000 claims abstract description 45
- 238000012549 training Methods 0.000 claims description 76
- 239000011159 matrix material Substances 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 48
- 230000004927 fusion Effects 0.000 claims description 15
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 claims description 15
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 10
- 239000000284 extract Substances 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 8
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 7
- 230000005856 abnormality Effects 0.000 claims description 7
- 238000013459 approach Methods 0.000 claims description 7
- 239000003054 catalyst Substances 0.000 claims description 7
- 229910052739 hydrogen Inorganic materials 0.000 claims description 7
- 239000001257 hydrogen Substances 0.000 claims description 7
- 239000007788 liquid Substances 0.000 claims description 7
- 238000006467 substitution reaction Methods 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 2
- 230000006870 function Effects 0.000 description 40
- 238000006116 polymerization reaction Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000002844 melting Methods 0.000 description 3
- 230000008018 melting Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000002405 diagnostic procedure Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 229920003023 plastic Polymers 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- HGAZMNJKRQFZKS-UHFFFAOYSA-N chloroethene;ethenyl acetate Chemical compound ClC=C.CC(=O)OC=C HGAZMNJKRQFZKS-UHFFFAOYSA-N 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 229920005672 polyolefin resin Polymers 0.000 description 1
- 239000004800 polyvinyl chloride Substances 0.000 description 1
- 229920000915 polyvinyl chloride Polymers 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920003002 synthetic resin Polymers 0.000 description 1
- 239000000057 synthetic resin Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种工业聚丙烯生产熔融指数检测故障诊断***,包括与工业过程对象连接的现场智能仪表、DCS***以及上位机,所述的DCS***由数据接口、控制站、数据库构成;智能仪表、DCS***、上位机依次相连,所述的上位机包括标准化处理模块、小波分解模块、主元分析功能模块、小波重构功能模块、支持向量机分类器功能模块以及故障判断模块。以及提出了一种故障诊断方法。本发明提供一种同时考虑工业过程数据的复共线性、非线性特性和多尺度特性、能够得到良好的诊断效果的的工业聚丙烯生产熔融指数检测故障诊断***及方法。
Description
(一)技术领域
本发明涉及工业过程故障诊断领域,特别地,涉及一种工业聚丙烯生产熔融指数检测故障诊断***及方法。
(二)背景技术
聚丙烯是以丙烯单体为主聚合而成的一种合成树脂,是塑料工业中的重要产品。在目前我国的聚烯烃树脂中,成为仅次于聚乙烯和聚氯乙烯的第三大塑料。在聚丙烯生产过程中,熔融指数(MI)是反映产品质量的一个重要指标,是生产质量控制和牌号切换的重要依据。但MI只能离线检测,一般离线分析至少需要近2小时,耗资而且耗时,特别是离线分析的2小时期间将无法及时了解聚丙烯生产过程的状态。因此,选取与熔融指数密切相关的易测变量作为二次变量,从中分析熔融指数,检测生产过程是否正常,对丙烯聚合生产过程进行监控就显得异常重要。
聚丙烯生产过程从本质上来说是具有多尺度特性的,比如微观、宏观聚合动力学的空间多尺度特性,过程变量总是在不同采样率或尺度下得到的,相应的控制或操作也是在不同的时间和空间尺度下发生的。只有全面考虑到过程的复共线性和非线性特性,同时考虑丙烯聚合生产过程的多尺度特性,才能对聚丙烯生产进行准确、可靠的故障诊断。但是目前的故障诊断只考虑了聚丙烯生产过程的复共线性和非线性特性,而没有考虑到过程的多尺度特性,往往难以得到较好的故障诊断效果。
(三)发明内容
为了克服已有的聚丙烯生产熔融指数检测故障诊断***的没有考虑过程的多尺度特性、难以得到较好的诊断效果的不足,本发明提供一种同时考虑工业过程数据的复共线性、非线性特性和多尺度特性、能够得到良好的诊断效果的工业聚丙烯生产熔融指数检测故障诊断***及方法。
本发明解决其技术问题所采用的技术方案是:
一种工业聚丙烯生产熔融指数检测故障诊断***,包括与聚丙烯生产过程对象连接的现场智能仪表、DCS***以及上位机,所述的DCS***由数据接口、控制站、数据库构成;智能仪表、DCS***、上位机依次相连,所述的上位机包括:
标准化处理模块,用于对数据库中采集***正常时关键变量数据进行标准化处理,各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
1)计算均值:
2)计算方差:
3)标准化:
其中,TX为训练样本,N为训练样本数,TX为训练样本的均值;
小波分解功能模块,用于采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息,依照分解层数,采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj。细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子;
④从而得到原始信息的分解信息,采用下式进行计算:
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k采用Mallat算法计算;
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
主元分析模块,用于进行主元分析提取主成分,依照主元分析方差提取率,采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑X;
②对∑X进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,
对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F;
小波重构功能模块,用于进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;
支持向量机分类器功能模块,用于核函数采用径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),依照置信概率,将训练过程化为如下二次规划求解问题:
得到分类函数,即如下函数的符号函数:
其中,αi(i=1,…,N)是拉各朗日乘子,xi(i=1,…,N)是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的的方向,b为决定超平面位置的参数,δ为核参数;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态;
信号采集模块,用于设定每次采样的时间间隙,采集现场智能仪表的信号;
待诊断数据确定模块,用于将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;
故障诊断模块,用于对待检测数据VX用训练时得到的TX和进行标准化处理,并将标准化处理后的数据作为小波分解模块的输入,用训练时相同的参数对输入数据进行小波分解,得到的系数作为中主元分析模块的输入;用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到小波重构模块;将对应的数据相加即得到原待测数据的主成分,并将所得到的主成分输入到支持向量机分类器模块;将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态。
作为优选的一种方案:所述的上位机还包括:判别模型更新模块,用于定期将过程状态正常的点添加到训练集VX中,输出到标准化处理模块、小波分解模块、主元分析模块、小波重构功能模块,并更新支持向量机分类器的分类模型。
作为优选的另一种方案:所述的上位机还包括:结果显示模块,用于将故障诊断结果传给DCS***,并在DCS的控制站显示过程状态,并通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示。
作为优选的再一种方案:所述的关键变量包括主催化剂流率f4、辅催化剂流率f5、三股丙烯进料流率(f1、f2、f3)、釜内流体温度T、釜内流体压强P、釜内液位l和釜内氢气体积浓度α。
一种工业聚丙烯生产熔融指数检测故障诊断方法,所述的故障诊断方法包括以下步骤:
(1)、确定故障诊断所用的关键变量,从DCS数据库的历史数据库中分别采集***正常和故障时所述变量的数据作为训练样本TX;
(2)、在小波分解模块、主元分析模块和支持向量机分类器模块中,分别设置小波分解层数、主元分析方差提取率、支持向量机核参数和置信概率参数,设定DCS中的采样周期;
(3)、训练样本TX在上位机中,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
3.1)计算均值:
3.2)计算方差:
3.3)标准化:
其中,N为训练样本数,N为训练样本数,TX为训练样本的均值;
(4)、用于采用Mallat塔式分解算法将采集的原始信号分解为一系列近似信息和细节信息,采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj。细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子;
④从而得到原始信息的分解信息,采用下式进行计算:
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k采用Mallat算法计算;
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
(5)、进行主元分析提取主成分,采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑X;
②对∑X进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F;
(6)、进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;
(7)、核函数采用径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),将训练过程化为如下二次规划求解问题:
得到分类函数,即如下函数的符号函数:
其中,αi(i=1,…,N)是拉各朗日乘子,xi(i=1,…,N)是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的的方向,b为决定超平面位置的参数,δ为核参数;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态;(8)、将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;对待检测数据VX用训练时得到的TX和进行标准化处理,并将标准化处理后的数据作为小波分解模块的输入,用训练时相同的参数对输入数据进行小波分解,得到的系数作为中主元分析模块的输入;用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到小波重构模块;将对应的数据相加即得到原待测数据的主成分,并将所得到的主成分输入到支持向量机分类器模块;将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态。
作为优选的一种方案:所述的故障诊断方法还包括:(9)、定期将过程状态正常的点添加到训练集VX中,重复(3)~(7)的训练过程,以便及时更新支持向量机分类器的分类模型。
作为优选的另一种方案:在所述的(8)中,计算判别函数值,并在上位机的人机界面上显示过程的状态,上位机将故障诊断结果传给DCS***,并在DCS的控制站显示过程状态,同时通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示。
作为优选的再一种方案:所述的关键变量包括主催化剂流率f4、辅催化剂流率f5、三股丙烯进料流率(f1、f2、f3)、釜内流体温度T、釜内流体压强P、釜内液位l和釜内氢气体积浓度α。
本发明同时考虑聚丙烯生产过程数据的复共线性、非线性特性和多尺度特性,将主元分析、支持向量机和小波分析相结合,对工业过程进行故障诊断。其中,主元分析用于处理工业过程生产数据的复相关性,支持向量机分类器用于解决非线性分类问题,小波分析用于获取过程在不同尺度下的信息。
本发明的有益效果主要表现在:同时兼顾聚丙烯生产过程数据中的复相关性、多尺度特性和非线性特性等特性,充分利用了主元分析,多尺度***理论和支持向量机各自的特点,将主元分析的解相关性能力,小波分析对信息不同尺度下的强分解与重构能力以及支持向量机的多变量非线性映射能力很好地结合了起来,发挥了各自的优势,使得故障诊断更加可靠有效,能更好的指导生产,提高生产效益。
(四)附图说明
图1是本发明所提出的故障诊断***的硬件结构图。
图2是本发明所提出的故障诊断***功能模块图。
图3是小波分析的Mallat算法分解图。
图4是聚丙烯生产流程简图。
图5是MS-PAC-SVM检测效果图。
图6是本发明上位机的原理框图。
(五)具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
参照图1、图2、图3、图4、图5以及图6,工业聚丙烯生产熔融指数检测故障诊断***,包括与聚丙烯生产过程对象1连接的现场智能仪表2、DCS***以及上位机6,所述的DCS***由数据接口3、控制站4、数据库5构成;智能仪表2、DCS***、上位机6通过现场总线依次相连,所述的上位机6包括:
标准化处理模块7,用于对数据进行标准化处理,各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
1)计算均值:
2)计算方差:
3)标准化:
其中,TX为训练样本,N为训练样本数,TX为训练样本的均值;
小波分解功能模块8,用于采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息,根据小波分解层数,采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj。细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子;
④从而得到原始信息的分解信息,采用下式进行计算:
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k采用Mallat算法计算;
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
主元分析模块9,用于进行主元分析提取主成分,依照设置的主元分析方差提取率,采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑X;
②对∑X进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F;
小波重构功能模块10,用于进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;
支持向量机分类器功能模块11,用于核函数采用径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),依照设置的置信概率,将训练过程化为如下二次规划求解问题:
得到分类函数,即如下函数的符号函数:
其中,αi(i=1,…,N)是拉各朗日乘子,xi(i=1,…,N)是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的的方向,b为决定超平面位置的参数,δ为核参数;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态;
信号采集模块12,用于设定每次采样的时间间隙,采集现场智能仪表的信号;
待诊断数据确定模块13,用于将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;
故障诊断模块14,用于对待检测数据VX用训练时得到的TX和进行标准化处理,并将标准化处理后的数据作为小波分解模块的输入,用训练时相同的参数对输入数据进行小波分解,得到的系数作为中主元分析模块的输入;用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到小波重构模块;将对应的数据相加即得到原待测数据的主成分,并将所得到的主成分输入到支持向量机分类器模块;将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态。
所述的上位机还包括:判别模型更新模块15,用于定期将过程状态正常的点添加到训练集VX中,输出到标准化处理模块7、小波分解模块8、主元分析模块9、小波重构功能模块10,并更新支持向量机分类器模块11的分类模型。
所述的上位机还包括:结果显示模块16,用于将故障诊断结果传给DCS,并在DCS的控制站显示过程状态,并通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示。
本实施例的工业聚丙烯生产熔融指数检测故障诊断***的硬件结构图如附图1所示,所述的故障诊断***核心由包括标准化模块7、小波分解模块8、主元分析模块9、小波重构模块10、支持向量机分类器模块11等五大功能模块和人机界面的上位机6构成,此外还包括:现场智能仪表2,DCS***和现场总线。所述的DCS***由数据接口3、控制站4、数据库5构成;丙烯聚合生产过程1、智能仪表2、DCS***、上位机6通过现场总线依次相连,实现信息流的上传和下达。故障诊断***在上位机6上运行,可以方便地与底层***进行信息交换,及时应对***故障。
本实施例述的故障诊断***的功能模块图如附图2所示,主要包括标准化处理模块7、小波分解模块8、主元分析模块9、小波重构模块10、支持向量机分类器模块11等五大功能模块。
所述的故障诊断方法按照如下步骤进行实施:
1、从DCS数据库5的历史数据库中采集***正常时以下九个变量的数据作为训练样本TX:主催化剂流率f4、辅催化剂流率f5、三股丙烯进料流率(f1、f2、f3)、釜内流体温度T、釜内流体压强P、釜内液位l和釜内氢气体积浓度α;
2、在上位机6的小波分解模块8、主元分析模块9和支持向量机分类器模块11中,分别设置小波分解层数、主元分析方差提取率、支持向量机核参数和置信概率等参数,并设定DCS中的采样周期;
3、训练样本TX在上位机6中,依次经过标准化处理7、小波分解8、主元分析9、小波重构10、支持向量机11等功能模块,采用以下步骤来完成诊断***的训练:
1)上位机6的标准化处理功能模块7,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X。采用以下步骤来实现:
①计算均值:
②计算方差:
③标准化:
其中N为训练样本数,,N为训练样本数,TX为训练样本的均值;
上位机6的标准化处理功能模块7所进行的标准化处理能消除各变量因为量纲不同造成的影响。
2)在上位机6的小波分解功能模块8中,采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息。所述的上位机6中小波分解模块8的小波分析采用db3小波,分解层数取为3—7。采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj。细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子。
④从而得到原始信息的分解信息,采用下式进行计算:
其中第一项代表近似信息,第二项为细节信息。逼近因子aJ,k与细节因子dj,k采用Mallat算法计算。其中所采用的Mallat算法的塔式分解图,如图2所示。
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
实际聚丙烯生产过程从本质上来说是具有多尺度特性的,各个尺度下的信息含量和所体现的***特性是不同的,如果直接对***进行建模,将忽略这种不同,从而导致会导致结果上的偏差。用小波分解提取各个尺度的信息,能更充分的挖掘过程数据所携带的信息,提升结果的精确度。
3)上位机6的主元分析模块9,进行主元分析,提取主成分。所述的主元分析总方差提取率大于80%,计算过程采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑X;
②对∑X进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F。
很显然,分析***在一个低维空间要比在一个高维空间容易得多。主元分析在力求数据信息丢失最少的原则下,对高维的变量空间降维,以得到丙烯聚合生产过程变量体系的少数几个线性组合,并且这几个线性组合所构成的综合变量将尽可能多地保留原丙烯聚合生产过程变量变异方面的信息。
4)上位机6的小波重构功能模块10,进行小波重构。
根据小波理论,将丙烯聚合生产过程各尺度下得到的主元相加,即得到丙烯聚合生产过程总的主元。
5)训练上位机6中的支持向量机分类器功能模块11的分类模型。
所述的上位机6中的支持向量机分类器功能模块9的核函数,采用径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),将训练过程化为如下二次规划求解问题:
从而得到分类函数,即如下函数的符号函数:
其中,αi(i=1,…,N)是拉各朗日乘子,xi(i=1,…,N)是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的的方向,b为决定超平面位置的参数,δ为核参数;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态。
支持向量机基于统计学习理论,采用结构风险最小化准则,很好地解决了工业聚丙烯生产熔融指数检测过程中存在的小样本、局部极小点、高维数等难题,用于工业聚丙烯生产熔融指数检测的分类问题能提高分类精度。
4、***开始投运:
1)用定时器,设置好每次采样的时间间隔;
2)现场智能仪表2检测过程数据并传送到DCS数据库5的实时数据库中;
3)上位机6在每个定时周期从DCS数据库5的实时数据库中,得到最新的变量数据,作为待诊断数据VX;
5)上位机6的小波分解模块8,用训练时同样的参数对输入数据进行小波分解,得到的系数作为上位机6中主元分析模块9的输入;
6)上位机6的主元分析模块9,用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到上位机6的小波重构模块10;
7)上位机6的小波重构模块10,将对应的数据相加即得到原待测数据的主成分,并将所得到的组成分输入到上位机6的支持向量机分类器模块11;
8)上位机6的支持向量机分类器模块11,将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态,并在上位机6的人机界面上显示过程的状态
9)上位机6将故障诊断结果传给DCS,并在DCS的控制站4显示过程状态,同时通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示,使得现场操作工可以及时应对。
5、分类器模型更新
在***投运过程中,定期将工业聚丙烯生产过程状态正常的点添加到训练集TX中,重复步骤3的训练过程,以便及时更新上位机6的支持向量机分类器11的分类模型,使分类器模型具有较好的分类效果。
下面详细说明本发明的一个具体实施例。
以聚丙烯生产HYPOL工艺实际工业生产为例。图三给出了典型的Hypol连续搅拌釜(CSTR)法生产聚丙烯的工艺流程图,前2釜是CSTR反应器、后2釜是流化床反应器(FBR)。选取主催化剂流率、辅催化剂流率、三股丙烯进料流率、釜内流体温度、釜内流体压强、釜内液位、釜内氢气体积浓度九个易测操作变量作为模型的输入量,从生产过程的DCS***中获取这九个参数的数据作为训练样本,其中五十个正常的样本作为训练集,另二十二个样本点作为测试集数据验证诊断效果。小波分解层数为5,PCA提取主成分个数为7,支持向量机核参数为5,置信概率0.98,采样周期为2小时。图5为MS-PAC-SVM检测效果图,图中只画出了前两个主成分的分布。表1列出了图5相对应的测试集中实际故障点和本***检测出的故障点,可以看出仅15号故障点漏报,误报率为0。显然,本***具有较高的诊断准确性。
实际故障点 | 1,2,12,15,16 |
检测故障点 | 1,2,12,16 |
表1。
实施例2
参照图1、图2、图3以及图4,一种工业聚丙烯生产熔融指数检测故障诊断方法,所述的故障诊断方法包括以下步骤:
(1)、确定故障诊断所用的关键变量,从DCS数据库的历史数据库中分别采集***正常和故障时所述变量的数据作为训练样本TX;
(2)、在小波分解模块8、主元分析模块9和支持向量机分类器模块11中,分别设置小波分解层数、主元分析方差提取率、支持向量机核参数和置信概率等参数,设定DCS中的采样周期;
(3)、训练样本TX在上位机6中,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
3.1)计算均值:
3.2)计算方差:
3.3)标准化:
其中,N为训练样本数。
(4)、用于采用Mallat塔式分解算法将采集的原始信号分解为一系列近似信息和细节信息,采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj。细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子;
④从而得到原始信息的分解信息,采用下式进行计算:
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k采用Mallat算法计算;
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
(5)、进行主元分析提取主成分,采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑X;
②对∑X进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F;
(6)、进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;
(7)、核函数采用径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),将训练过程化为如下二次规划求解问题:
得到分类函数,即如下函数的符号函数:
其中,αi(i=1,…,N)是拉各朗日乘子,xi(i=1,…,N)是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的的方向,b为决定超平面位置的参数,δ为核参数;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态;(8)、将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;对待检测数据VX用训练时得到的TX和进行标准化处理,并将标准化处理后的数据作为小波分解模块的输入,用训练时相同的参数对输入数据进行小波分解,得到的系数作为中主元分析模块的输入;用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到小波重构模块;将对应的数据相加即得到原待测数据的主成分,并将所得到的主成分输入到支持向量机分类器模块;将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态。
所述的故障诊断方法还包括:(9)、定期将过程状态正常的点添加到训练集VX中,重复(3)~(7)的训练过程,以便及时更新支持向量机分类器的分类模型。
在所述的(8)中,计算判别函数值,并在上位机的人机界面上显示过程的状态,上位机将故障诊断结果传给DCS***,并在DCS的控制站显示过程状态,同时通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示。
所述的关键变量包括主催化剂流率f4、辅催化剂流率f5、三股丙烯进料流率(f1、f2、f3)、釜内流体温度T、釜内流体压强P、釜内液位l和釜内氢气体积浓度α。
Claims (7)
1.一种工业聚丙烯生产熔融指数检测故障诊断***,包括与聚丙烯生产过程对象连接的现场智能仪表、DCS***以及上位机,所述的DCS***由数据接口、控制站、数据库构成;智能仪表、DCS***、上位机依次相连,其特征在于:所述的上位机包括:
标准化处理模块,用于对数据库中采集***正常时关键变量数据进行标准化处理,各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
1)计算均值:
2)计算方差:
3)标准化:
其中,TX为训练样本,N为训练样本数,TX为训练样本的均值;
小波分解功能模块,用于采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息,根据小波分解层数,采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj,细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子;
④从而得到原始信息的分解信息,采用下式进行计算:
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k采用Mallat算法计算;
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
主元分析模块,用于进行主元分析提取主成分,依照设置的主元分析方差提取率,采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑x;
②对∑x进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F;
小波重构功能模块,用于进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;
支持向量机分类器功能模块,用于将核函数设为径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),依照设置的置信概率,将训练过程化为如下二次规划求解问题:
得到分类函数,即如下函数的符号函数:
其中,αi,i=1,…,N是拉各朗日乘子,xi,i=1,…,N是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的方向,b为决定超平面位置的参数,σ为基宽度;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态;
信号采集模块,用于设定每次采样的时间间隙,采集现场智能仪表的信号;
待诊断数据确定模块,用于将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;
故障诊断模块,用于对待诊断数据VX用训练时得到的TX和进行标准化处理,并将标准化处理后的数据作为小波分解模块的输入,用训练时相同的参数对输入数据进行小波分解,得到的系数作为主元分析模块的输入;用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到小波重构模块;将对应的数据相加即得到原待测数据的主成分,并将所得到的组成分输入到支持向量机分类器模块;将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态。2、如权利要求1所述的工业聚丙烯生产熔融指数检测故障诊断***,其特征在于:所述的上位机还包括:
判别模型更新模块,用于定期将过程状态正常的点添加到训练集中,输出到标准化处理模块、小波分解模块、主元分析模块、小波重构功能模块,并更新支持向量机分类器的分类模型。
3.如权利要求1或2所述的工业聚丙烯生产熔融指数检测故障诊断***,其特征在于:所述的上位机还包括:
结果显示模块,用于将故障诊断结果传给DCS***,并在DCS的控制站显示过程状态,并通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示。
4.如权利要求3所述的工业聚丙烯生产熔融指数检测故障诊断***,其特征在于:所述关键变量包括主催化剂流率f4、辅催化剂流率f5、三股丙烯进料流率f1、f2、f3、釜内流体温度T、釜内流体压强P、釜内液位l和釜内氢气体积浓度。
5.一种用如权利要求1所述的工业聚丙烯生产熔融指数检测故障诊断***实现的故障诊断方法,其特征在于:所述的故障诊断方法包括以下步骤:
1)、确定故障诊断所用的关键变量,从DCS数据库的历史数据库中分别采集***正常和故障时所述变量的数据作为训练样本TX;
2)、在小波分解模块、主元分析模块和支持向量机分类器模块中,分别设置小波分解层数、主元分析方差提取率、支持向量机核参数和置信概率参数,并设定DCS中的采样周期;
3)、训练样本TX在上位机中,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
3.1)计算均值:
3.2)计算方差:
3.3)标准化:
其中,TX为训练样本,N为训练样本数,TX为训练样本的均值;
4)、用于采用Mallat塔式分解算法将采集的原始信号分解为一系列近似信息和细节信息,采用以下步骤来实现:
①将原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度;
③计算细节空间Wj。细节空间Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,采用下式来计算:
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)} (5)
其中j是尺度因子,k是平移因子;
④从而得到原始信息的分解信息,采用下式进行计算:
具中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k采用Mallat算法计算;
所述的近似信息AJf(t)与细节信息Djf(t)(j=1,2,...,J),定义如下:
5)、进行主元分析提取主成分,采用协方差奇异值分解的方法,采用以下步骤来实现:
①计算X的协方差阵,记为∑x;
②对∑x进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;
③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;
④选取特征向量矩阵U的前k列,作为变换矩阵T;
⑤计算主元,通过式子F=T×X计算得到主元F;
6)、进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;
7)、核函数采用径向基函数K(xi,x)=exp(-‖x-xi‖/σ2),将训练过程化为如下二次规划求解问题:
得到分类函数,即如下函数的符号函数:
其中,αi,i=1,…,N是拉各朗日乘子,xi,i=1,…,N是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的方向,b为决定超平面位置的参数,σ为基宽度;
定义当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态;
6.如权利要求5所述的故障诊断方法,其特征在于:所述的故障诊断方法还包括:9)、定期将过程状态正常的点添加到训练集中,重复(3)~(7)的训练过程,以便及时更新支持向量机分类器的分类模型。
7.如权利要求5或6所述的故障诊断方法,其特征在于:在所述的8)中,计算判别函数值,并在上位机的人机界面上显示过程的状态,上位机将故障诊断结果传给DCS***,并在DCS的控制站显示过程状态,同时通过DCS***和现场总线将过程状态信息传递到现场操作站进行显示。
8.如权利要求7所述的故障诊断方法,其特征在于:所述的关键变量包括主催化剂流率f4、辅催化剂流率f5、三股丙烯进料流率f1、f2、f3、釜内流体温度T、釜内流体压强P、釜内液位l和釜内氢气体积浓度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101554147A CN100461043C (zh) | 2006-12-22 | 2006-12-22 | 工业聚丙烯生产熔融指数检测故障诊断***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101554147A CN100461043C (zh) | 2006-12-22 | 2006-12-22 | 工业聚丙烯生产熔融指数检测故障诊断***及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101046687A CN101046687A (zh) | 2007-10-03 |
CN100461043C true CN100461043C (zh) | 2009-02-11 |
Family
ID=38771351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101554147A Expired - Fee Related CN100461043C (zh) | 2006-12-22 | 2006-12-22 | 工业聚丙烯生产熔融指数检测故障诊断***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100461043C (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470426B (zh) * | 2007-12-27 | 2011-02-16 | 北京北方微电子基地设备工艺研究中心有限责任公司 | 一种故障检测的方法和*** |
JP6085862B2 (ja) * | 2012-10-10 | 2017-03-01 | 陳山 鵬 | マルチ・コンディション・モニターを用いた状態監視方法および状態監視装置システム |
CN103675006B (zh) * | 2013-09-22 | 2015-09-09 | 浙江大学 | 最小二乘的工业熔融指数软测量仪表和方法 |
CN103674778B (zh) * | 2013-09-22 | 2015-08-26 | 浙江大学 | Rbf粒子群优化的工业熔融指数软测量仪表和方法 |
CN105867345B (zh) * | 2016-03-24 | 2018-04-27 | 浙江科技学院 | 一种多变量化工过程的故障源和故障传播路径定位方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1241721A (zh) * | 1999-03-22 | 2000-01-19 | 西安交通大学 | 基于非线性频谱分析的故障检测与诊断方法 |
US20010048367A1 (en) * | 2000-06-02 | 2001-12-06 | Richard Wimmer | System and method for diagnosing fault conditions associated with powering an electrical load |
CN1415685A (zh) * | 2001-10-30 | 2003-05-07 | 中国石油化工股份有限公司 | 控制聚酰胺热熔胶熔融指数的方法 |
US20030196492A1 (en) * | 2002-04-17 | 2003-10-23 | Remboski Donald J. | Fault detection system having audio analysis and method of using the same |
CN1451689A (zh) * | 2003-05-09 | 2003-10-29 | 东华大学 | 超高熔融指数聚丙烯及其制备方法 |
-
2006
- 2006-12-22 CN CNB2006101554147A patent/CN100461043C/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1241721A (zh) * | 1999-03-22 | 2000-01-19 | 西安交通大学 | 基于非线性频谱分析的故障检测与诊断方法 |
US20010048367A1 (en) * | 2000-06-02 | 2001-12-06 | Richard Wimmer | System and method for diagnosing fault conditions associated with powering an electrical load |
CN1415685A (zh) * | 2001-10-30 | 2003-05-07 | 中国石油化工股份有限公司 | 控制聚酰胺热熔胶熔融指数的方法 |
US20030196492A1 (en) * | 2002-04-17 | 2003-10-23 | Remboski Donald J. | Fault detection system having audio analysis and method of using the same |
CN1451689A (zh) * | 2003-05-09 | 2003-10-29 | 东华大学 | 超高熔融指数聚丙烯及其制备方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101046687A (zh) | 2007-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100461044C (zh) | 一种丙烯聚合生产中熔融指数检测的故障诊断***及方法 | |
CN100480926C (zh) | 一种基于小波分析的工业生产过程故障诊断***及方法 | |
CN100517141C (zh) | 丙烯聚合生产数据检测及故障诊断***及方法 | |
CN100470417C (zh) | 工业生产过程小样本条件下的故障诊断***及方法 | |
CN109459993B (zh) | 一种流程工业过程在线自适应故障监测与诊断方法 | |
CN103914064B (zh) | 基于多分类器和d-s证据融合的工业过程故障诊断方法 | |
CN106872657B (zh) | 一种多变量水质参数时序数据异常事件检测方法 | |
US8868985B2 (en) | Supervised fault learning using rule-generated samples for machine condition monitoring | |
Mounce et al. | Pattern matching and associative artificial neural networks for water distribution system time series data analysis | |
CN101305327A (zh) | 与多元统计分析一起用于过程中的故障检测和隔离及异常情况预防的统计特征 | |
CN105955219A (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN100461043C (zh) | 工业聚丙烯生产熔融指数检测故障诊断***及方法 | |
CN107065834B (zh) | 湿法冶金过程中浓密机的故障诊断方法 | |
CN103901880A (zh) | 基于多分类器和d-s证据融合的工业过程故障检测方法 | |
Dorgo et al. | Decision trees for informative process alarm definition and alarm-based fault classification | |
CN111580506A (zh) | 基于信息融合的工业过程故障诊断方法 | |
CN109298633A (zh) | 基于自适应分块非负矩阵分解的化工生产过程故障监测方法 | |
CN104914850A (zh) | 基于切换线性动态***模型的工业过程故障诊断方法 | |
CN116108371B (zh) | 基于级联异常生成网络的云服务异常诊断方法与*** | |
CN106647650A (zh) | 基于变量加权pca模型的分散式工业过程监测方法 | |
CN105334823A (zh) | 基于有监督的线性动态***模型的工业过程故障检测方法 | |
CN102601881A (zh) | 一种橡胶硬度的在线质量监控和预测模型更新的方法 | |
CN110009126A (zh) | 基于pls模型与pca贡献度融合的在线报警分析方法 | |
CN111122811A (zh) | 一种oica和rnn融合模型的污水处理过程故障监测方法 | |
CN116380445A (zh) | 基于振动波形的设备状态诊断方法及相关装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090211 Termination date: 20111222 |