CN109471420B - 基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 - Google Patents
基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 Download PDFInfo
- Publication number
- CN109471420B CN109471420B CN201811109196.2A CN201811109196A CN109471420B CN 109471420 B CN109471420 B CN 109471420B CN 201811109196 A CN201811109196 A CN 201811109196A CN 109471420 B CN109471420 B CN 109471420B
- Authority
- CN
- China
- Prior art keywords
- matrix
- new
- variables
- monitoring
- air preheater
- 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
Images
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/0243—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 model based detection method, e.g. first-principles knowledge model
Abstract
本发明公开了一种基于CVA‑SFA的智能电厂大型燃煤发电机组空气预热器控制性能监测方法。针对大型火力发电机组空气预热器因参数众多、机组运行工况多变导致的控制性能监测困难的问题,运用典型变量分析(CVA)提取空气预热器控制系统变量间的相关信息,再利用慢特征分析(SFA)算法提取相关信息中的动态信息。最后,结合变量的相关性和变化快慢信息构造空气预热器控制性能在线监测模型。该方法克服了大型火电机组空气预热器因变量众多、工况变化而导致的控制性能监测困难的问题,大大提高了动态过程控制性能在线监测的准确度,有助于火电厂对空预器控制系统进行智能及时的监测,有助于保证大型火力发电机组的安全可靠运行。
Description
技术领域
本发明属于火电过程控制系统性能监测领域,特别是涉及一种面向智能电 厂大型燃煤发电机组空气预热器运行相关性信息与动态信息的在线性能监测方 法。
背景技术
控制系统在现代化的工业过程中占据非常重要的地位,生产质量、操作安 全、物能消耗等影响经济效益的指标都直接或间接的与控制系统的性能有关。 在实际生产过程中,控制系统在投入使用初期往往性能表现良好,但运行一段 时间后,由于设备的磨损、定期保养和维护不及时等原因,可能导致控制系统 的性能下降,控制性能变差会直接影响生产质量,导致经济效益亏损,若因此 引发生产故障,还会涉及到人的生命安全甚至社会企业的财产安全,带来极大 威胁。Torrres等人对2004-2005年巴西12家工厂(石化、造纸、水泥、钢铁、 采矿等),超过700个控制回路进行检验,结果显示14%回路的阀门磨损过度,15%的阀门存在迟滞问题,16%的回路存在严重的整定问题,24%的控制器输出 存在饱和现象,41%的回路因为整定问题、耦合、扰动以及执行器的问题而存在 振荡现象。
另外,实际生产中,一个生产过程可能会有数以千计的控制回路共同作用,Eastman化学公司中的两个精馏生产设备拥有多大14000个控制回路,在HVAC 生产过程中,其控制回路的数量甚至能够达到十万个。大型火力发电机组具有 较高的复杂性,具体体现在规模庞大、设备众多、参数多样化且相互影响等方 面。此外,大规模的火力发电机组,现场具有高温、高压及高噪声等特点,不 适合人工去现场考经验判断控制系统性能的优劣与变化。
控制性能评价与监测技术是过程控制领域新兴的一项重要技术,它能够利 用设备的日常运行数据,实时监测监视系统控制性能的变化,对控制系统的问 题做出早期识别和优化。对于发电机组,由于电力系统中的用电负荷是经常变 化的,为了维持有功功率平衡,保持系统频率稳定,需要发电部门相应的改变 发电机的出力以适应用电负荷的变化,即发电机组的工况不是稳定不变的。但 是现有的控制性能评价与监测方法如主成分分析、偏最小二乘法、费舍尔判别 分析,都是基于工况稳定的理想假设下进行的,因此,将其运用在大型火电机 组空气预热器控制系统性能监测上,并不能得到很好的监测效果。
发明内容
空气预热器是安装在锅炉后部,利用锅炉尾部烟气的热量加热燃料燃烧所 需空气以提高锅炉效率的热交换装置。加热后的空气(热空气或热风)温度达 200~400℃,这对改善燃料着火和燃料完全燃烧十分有利,锅炉煤粉制备系 统也需要由其供给热风作为干燥剂和供输送煤粉之用,空气预热器是现代大、 中型锅炉必不可少的设备。影响空气预热器性能的关键问题有:振动、噪声、 漏风、腐蚀和堵灰。
本发明的目的在于针对大型火电机组空气预热器因参数众多、结构复杂、 工况多变导致的控制性能监测困难的问题,运用典型变量分析与慢特征分析 融合算法提取空气预热器控制系统变量间的相关信息与变化快慢信息,克服 了大型空气预热器因变量众多、工况变化而导致的控制性能监测困难的问题。
本发明的目的是通过以下技术方案实现:基于CVA-SFA的智能电厂大型 燃煤发电机组空气预热器控制性能监测方法,该方法包括以下步骤:
(1)获取训练数据:设空气预热器的控制系统具有J个测量变量和操作变 量,每一次采样可以得到一个J×1的观测向量yk,其中下标k为采样时刻,采 样N次后得到的数据表述为一个二维观测矩阵所述 测量变量为空气预热器运行过程中可被测量的状态参数,包括空预器二次风出 口风温、引风机耗电率、增压风机电流、增压风机动叶开度、预热器漏风率、 预热器出口氧量、送风温度、再热温度、排热温度、燃料温度、密封空气加热 后压力、密封空气加热后温度等;所述操作变量包括引风机功率、引风机电流、 锅炉变总功率、汽机变总功率等;训练数据应当选取空气预热器在正常运行状 态下的采样数据。
(2)利用CVA算法提取数据的时序相关信息,该步骤通过以下子步骤实 现:
分别用所有的过去观测向量和将来观测向量构建过去观测矩阵Yp和将来 观测矩阵Yf:
其中,M=N-f-p+1,p,f为两类时滞参数,令p=f,其值可以通过 样本自相关函数来确定:
其中:autocorr(Yj,p)表示矩阵Yp第j个列向量与其时滞p的自相关系数;
(2.2)构建Hankel矩阵:计算过去矩阵和将来矩阵的协方差矩阵∑pp,∑ff以及他们的互协方差矩阵∑fp,再利用协方差与互协方差矩阵构建Hankel矩阵 H:
(2.3)奇异值分解:对Hankel矩阵进行奇异值分解可以得到Jp组典型变 量配对,用(ai TYp,bi TYf)表示第i组典型变量配对,ai T、bi T表示第i组典型变 量配对间的相关系数:
H=UDVT (6)
U和V分别为奇异向量ui,vi组成的正交矩阵,D为奇异值矩阵,U、V中 的奇异向量只成对相关,且相关性大小由D中对应的第i个奇异值γi表征。奇异 值越大(γ1>γ2>…>γJp),典型变量间的相关性越大。
Cr表示准则值,β为判断阈值,β=0.5。
由Vr计算典型变量转换矩阵C和残差变量转换矩阵L:
再利用转换矩阵可以得到典型变量空间Z和残差空间E:
Z,Ε中的列向量zk∈r×1,εk∈Jp×1分别表示在采样时刻k的典型变量 和残差变量;Z,Ε中的行向量zt,εt包含了同一变量在不同时刻的时序信息。
(3)利用慢特征分析算法(Slow Feature Analysis,SFA)提取典型变量空 间Z的变化快慢特征,该方法主要步骤如下:
(3.1)数据标准化:对典型变量空间Z按变量进行标准化处理,计算公式 如下:
zt表同一变量在不同时刻的时序向量,mean(zt)表示zt的均值,std(zt)表 示zt的标准差。
(3.2)Z经过投影后的输出信号为sj,sj表示第j个慢特征序列。考虑线 性条件下, 表示系数向量,这等价于寻找一个从标准化输入 信号Z中提取慢特征信号s=[s1 T,s2 T,…,sr T]T的转换矩阵 即s=WZ。慢特征信号sj要满足的目标函数及约 束条件为:
目标函数:
约束条件为:
(3.3)白化:利用奇异值分解,对输入数据的协方差矩阵<ZZT>进行白化 处理可以去除数据中的相关性,使提取出的慢特征值携带不同的信息:
其中:Λ-1/2BT为白化矩阵,Ο为对应的白化后的输入信号。
W=PΛ-1/2BT (17)
(4)划分慢特征s:最慢的特征对应最小的特征值,将特征值由小到大排 列,并依据特征值大小将前l个特征划分为s中变化较慢的特征,用sd表示;将 后(r-l)个特征划分为s中变化较快的特征,用se表示。划分依据l的确定方法 为,首先利用慢特征值s的变化快慢表示过程变量的变化快慢:
R=W-1
将提取出的慢特征值中缓慢度比输入数据缓慢度还要大的特征划分为快特 征,一共有Me个这样的快特征:
这里card{·}表示集合{·}中元素个数。根据式(19)确定的Me值,对应将 矩阵Ω也划分成两部分:
(5)计算动态监测指标:从典型变量空间的第一个样本点开始,每个样本 点可以得到一组动态监测指标(Sd 2,Se 2)。
(6)确定基于动态监测指标的控制限:利用核密度估计的方法,先估计出 动态监测指标Sd 2的概率密度函数p(x),对于给定显著性水平α,Sd 2的控制限 Sd 2 UCL的计算方式为:
以同样的方法可以计算出Se 2的控制限Se 2 UCL。
(7)对残差变量建立统计指标Qk:
其中:εi,k表示残差空间E中第i个变量在时刻k的取值,即E中第i行第k列 的元素。
同样用核密度估计的方法计算Qk统计指标的控制限QUCL:
(8)在线监测控制性能:基于步骤(2)到(4)建立的CVA-SFA模型、 步骤(5)到步骤(7)所得的三个监测统计量在线监测空气预热器控制系统的 性能状态。该步骤由以下子步骤来实现:
(8.1)获取新在线数据以及新数据预处理:采集到新的一段观测数据后,其中,下表new表示新观测数据,首先按照步 骤(2)将Ynew拓展成过去矩阵,并根据步骤(2)中获得的均值和标准差对过 去矩阵进行标准化处理得到Ypnew。
(8.2)提取出新观测数据的典型变量和残差变量:标准化处理后,利用步 骤(2)确定的转换矩阵Vr和L计算出新观测数据的典型变量空间Znew和残差空 间Enew。
(8.3)提取新观测数据的典型变量空间Znew中的慢特征:首先按照步骤(3.1) 中确定的均值和方差对Znew进行标准化处理,之后利用步骤(3.4)中确定的慢 特征转换矩阵W,提取出标准化Znew的慢特征snew,并按照之前的划分参数将 snew划分成sdnew和sfnew。
(8.4)计算新监测统计指标:根据建立的模型以及步骤(5)(7)中确定的 计算方法,计算两组典型变量空间下的监测统计指标Sd 2 new,Se 2 new和一组残差 空间监测指标Qknew:
(8.5)在线判断空气预热器控制性能状态:实时比较三个监测指标与其各 自的统计控制限,若三个监测指标都位于统计控制限之内,表明控制系统正常 工作;若有一个或以上监测指标超出正常控制限,表明控制系统有异常状况发 生。
本发明的有益效果在于:本发明针对大型火电机组空气预热器因参数众多、 结构复杂、工况多变导致的控制性能监测困难的问题,提出了基于CVA-SFA的 智能电厂大型燃煤发电机组空气预热器控制性能监测方法,该方法运用典型变 量分析提取空气预热器控制系统变量间的相关信息,再利用慢特征分析算法, 提取相关信息中的动态信息。最后,结合变量的相关性和变化快慢信息构造空 气预热器控制性能在线监测模型。该方法克服了大型空气预热器因变量众多、 工况变化而导致的控制性能监测困难的问题,大大提高了动态过程控制性能在 线监测的准确度,有助于火电厂对空气预热器系统进行有效及时的监测,有助 于保证大型火力发电机组的安全可靠运行,同时达到提高其生产效益的生产需求。
附图说明
图1是本发明面向智能电厂大型燃煤发电机组空气预热器在线性能监测方 法的流程图,(a)为离线建模过程流程图,(b)为在线监测过程流程图;
图2是本发明CVA-SFA方法用于统计过程监测的结果图,(a)为正常情况 下监测结果图,(b)为异常情况下监测结果图。
具体实施方式
下面结合附图及具体实例,对本发明作进一步详细说明。
本发明以浙能集团下属嘉华电厂#5号机组空气预热器为例,该机组的功率 为60万千瓦,为大型火电机组,包括63个过程变量,包括空预器二次风出口 风温、引风机耗电率、增压风机电流、增压风机动叶开度、预热器漏风率、预 热器出口氧量、送风温度、再热温度、排热温度、燃料温度、密封空气加热后 压力、密封空气加热后温度、引风机功率、引风机电流、锅炉变总功率、汽机 变总功率等。
应该理解,本发明不止局限于上述实例的火电发电过程,凡是熟悉本领域 的技术人员在不违背本发明的前提下还可以做出等同变型或替换,这些等同的 变型或替换均包含在本申请权利要求所限定的范围内。
如图1所示,本发明是基于CVA-SFA的智能电厂大型燃煤发电机组空气预 热器控制性能监测方法,包括以下步骤:
(1)获取训练数据:设空气预热器的控制系统具有J个测量变量和操作变 量,每一次采样可以得到一个J×1的观测向量yk,其中下标k为时间指标,采 样N次后得到的数据表述为一个二维观测矩阵本实 例中,采样周期为10分钟,共4654个样本,63个过程变量,所测变量为机组 运行过程中空气预热器的流量、功率、温度、压力、开度等状态变量。
(2)利用CVA算法提取数据的时序相关信息,该步骤通过以下子步骤实 现:
分别用所有的过去观测向量和将来观测向量构建过去观测矩阵Yp和将来 观测矩阵Yf:
其中,M=N-f-p+1,p,f为两类时滞参数,令p=f,其值可以通过 样本自相关函数来确定:
其中:autocorr(Yj,p)表示矩阵Yp第j个列向量与其时滞p的自相关系数;
(2.2)构建Hankel矩阵:计算过去矩阵和将来矩阵的协方差矩阵∑pp,∑ff以及他们的互协方差矩阵∑fp,再利用协方差与互协方差矩阵构建Hankel矩阵 H:
(2.3)奇异值分解:对Hankel矩阵进行奇异值分解可以得到Jp组典型变 量配对,用(ai TYp,bi TYf)表示第i组典型变量配对,ai T、bi T表示第i组典型变 量配对间的相关系数:
H=UDVT (6)
U和V分别为奇异向量ui,vi组成的正交矩阵,D为奇异值矩阵,U、V中 的奇异向量只成对相关,且相关性大小由D中对应的第i个奇异值γi表征。奇异 值越大(γ1>γ2>…>γJp),典型变量间的相关性越大。
Cr表示准则值,β为判断阈值,β=0.5。
由Vr计算典型变量转换矩阵C和残差变量转换矩阵L:
再利用转换矩阵可以得到典型变量空间Z和残差空间E:
Z,Ε中的列向量zk∈r×1,εk∈Jp×1分别表示在采样时刻k的典型变量 和残差变量;Z,Ε中的行向量zt,εt包含了同一变量在不同时刻的时序信息。
(3)利用慢特征分析算法(Slow Feature Analysis,SFA)提取典型变量空 间Z的变化快慢特征,该方法主要步骤如下:
(3.1)数据标准化:对典型变量空间Z按变量进行标准化处理,计算公式 如下:
zt表同一变量在不同时刻的时序向量,mean(zt)表示zt的均值,std(zt)表 示zt的标准差。
(3.2)Z经过投影后的输出信号为sj,sj表示第j个慢特征序列。考虑线 性条件下, 表示系数向量,这等价于寻找一个从标准化输入 信号Z中提取慢特征信号s=[s1 T,s2 T,…,sr T]T的转换矩阵 即s=WZ。慢特征信号sj要满足的目标函数及约 束条件为:
目标函数:
约束条件为:
(3.3)白化:利用奇异值分解,对输入数据的协方差矩阵<ZZT>进行白化 处理可以去除数据中的相关性,使提取出的慢特征值携带不同的信息:
其中:Λ-1/2BT为白化矩阵,Ο为对应的白化后的输入信号。
W=PΛ-1/2BT (17)
(4)划分慢特征s:最慢的特征对应最小的特征值,将特征值由小到大排 列,并依据特征值大小将前l个特征划分为s中变化较慢的特征,用sd表示;将 后(r-l)个特征划分为s中变化较快的特征,用se表示。划分依据l的确定方法 为,首先利用慢特征值s的变化快慢表示过程变量的变化快慢:
R=W-1
将提取出的慢特征值中缓慢度比输入数据缓慢度还要大的特征划分为快特 征,一共有Me个这样的快特征:
这里card{·}表示集合{·}中元素个数。根据式(19)确定的Me值,对应将 矩阵Ω也划分成两部分:
(5)计算动态监测指标:从典型变量空间的第一个样本点开始,每个样本 点可以得到一组动态监测指标(Sd 2,Se 2)。
(6)确定基于动态监测指标的控制限:利用核密度估计的方法,先估计出 动态监测指标Sd 2的概率密度函数p(x),对于给定显著性水平α,Sd 2的控制限 Sd 2 UCL的计算方式为:
以同样的方法可以计算出Se 2的控制限Se 2 UCL。
(7)对残差变量建立统计指标Qk:
其中:εi,k表示残差空间E中第i个变量在时刻k的取值,即E中第i行第k列 的元素。
同样用核密度估计的方法计算Qk统计指标的控制限QUCL:
(8)在线监测控制性能:基于步骤(2)到(4)建立的CVA-SFA模型、 步骤(5)到步骤(7)所得的三个监测统计量可以在线监测空气预热器控制系 统的性能状态。该步骤由以下子步骤来实现:
(8.1)获取新在线数据以及新数据预处理:采集到新的一段观测数据后,其中,下表new表示新观测数据,首先按照步 骤(2)将Ynew拓展成过去矩阵,并根据步骤(2)中获得的均值和标准差对过 去矩阵进行标准化处理得到Ypnew。本实例中,新数据共有两份,数据一为正常 工况下采集的数据,采样周期为10分钟,共3100个样本,63个过程变量,数 据二为发生异常工况下记录的数据,采样周期为10分钟,共1862个样本,63 个过程变量,所测变量为机组运行过程中空预器的流量、振动、温度、压力、 阀门开度等;
(8.2)提取出新观测数据的典型变量和残差变量:标准化处理后,利用步 骤(2)确定的转换矩阵Vr和L计算出新观测数据的典型变量空间Znew和残差空 间Enew。
(8.3)提取新观测数据的典型变量空间Znew中的慢特征:首先按照步骤(3.1) 中确定的均值和方差对Znew进行标准化处理,之后利用步骤(3.4)中确定的慢 特征转换矩阵W,提取出标准化Znew的慢特征snew,并按照之前的划分参数将 snew划分成sdnew和sfnew。
(8.4)计算新监测统计指标:根据建立的模型以及步骤(5)(7)中确定的 计算方法,计算两组典型变量空间下的监测统计指标Sd 2 new,Se 2 new和一组残差 空间监测指标Qknew:
(8.5)在线判断空气预热器控制性能状态:实时比较三个监测指标与其各 自的统计控制限,若三个监测指标都位于统计控制限之内,表明控制系统正常 工作;若有一个或以上监测指标超出正常控制限,表明控制系统有异常状况发 生。图2(a)中,三组统计量及对应控制线中,仅个别点的统计量超过了控制 线,在置信水平α=0.05的条件下,可以认为新工况数据是正常的,即控制系统 表现正常;图2(b)中,统计量Sd 2,Se 2均在第944-961、1080–1104、1222–1264 三段明显超出阈值线,统计量Qk在第326个采样点左右第一次明显超限后一直 维持超限状态,并同样在第944–961、1080–1104、1222–1264三段明显超限, 据此可以判断控制系统发生异常,这时就可以采用适当的故障诊断方法,比如 贡献图方法分析隔离出可能的故障变量。
本发明运用典型变量分析提取空气预热器控制系统变量间的相关信息,再 利用慢特征分析算法,提取相关信息中的动态特征,由该方法提取的特征可以 反应控制器的调节作用。最后,结合变量的相关性和变化快慢信息构造空气预 热器控制性能在线监测模型,该方法克服了大型空气预热器因变量众多、工况 变化而导致的控制性能监测困难的问题,大大提高了动态过程控制性能在线监 测的准确度,有助于火电厂对空气预热器系统进行有效及时的监测,有助于保 证大型火力发电机组的安全可靠运行,同时达到提高其生产效益的生产需求。
Claims (1)
1.基于CVA-SFA的智能电厂大型燃煤发电机组空气预热器控制性能监测方法,其特征在于,该方法包括以下步骤:
(1)获取训练数据:设空气预热器的控制系统具有J个测量变量和操作变量,每一次采样可以得到一个J×1的观测向量yk,其中下标k为采样时刻,采样N次后得到的数据表述为一个二维观测矩阵所述测量变量为空气预热器运行过程中可被测量的状态参数,包括空预器二次风出口风温、引风机耗电率、增压风机电流、增压风机动叶开度、预热器漏风率、预热器出口氧量、送风温度、再热温度、排热温度、燃料温度、密封空气加热后压力、密封空气加热后温度;所述操作变量包括引风机功率、引风机电流、锅炉变总功率、汽机变总功率;训练数据应当选取空气预热器在正常运行状态下的采样数据;
(2)利用CVA算法提取数据的时序相关信息,该步骤通过以下子步骤实现:
分别用所有的过去观测向量和将来观测向量构建过去观测矩阵Yp和将来观测矩阵Yf:
其中,M=N-f-p+1,p,f为两类时滞参数,令p=f,其值可以通过样本自相关函数来确定:
其中:autocorr(Yj,p)表示矩阵Yp第j个列向量与其时滞p的自相关系数;
(2.2)构建Hankel矩阵:计算过去观测矩阵和将来观测矩阵的协方差矩阵∑pp,∑ff以及他们的互协方差矩阵∑fp,再利用协方差与互协方差矩阵构建Hankel矩阵H:
(2.3)奇异值分解:对Hankel矩阵进行奇异值分解可以得到Jp组典型变量配对,用(ai TYp,bi TYf)表示第i组典型变量配对,ai T、bi T表示第i组典型变量配对间的相关系数:
H=UDVT (6)
U和V分别为奇异向量ui,vi组成的正交矩阵,D为奇异值矩阵,U、V中的奇异向量只成对相关,且相关性大小由D中对应的第i个奇异值γi表征;奇异值越大,典型变量间的相关性越大;
Cr表示准则值,β为判断阈值,β=0.5;
由Vr计算典型变量转换矩阵C和残差变量转换矩阵L:
再利用转换矩阵可以得到典型变量空间Z和残差空间E:
Z,Ε中的列向量zk∈r×1,εk∈Jp×1分别表示在采样时刻k的典型变量和残差变量;Z,Ε中的行向量zt,εt包含了同一变量在不同时刻的时序信息;
(3)利用慢特征分析算法(Slow Feature Analysis,SFA)提取典型变量空间Z的变化快慢特征,该方法主要步骤如下:
(3.1)数据标准化:对典型变量空间Z按变量进行标准化处理,计算公式如下:
zt表同一变量在不同时刻的时序向量,mean(zt)表示zt的均值,std(zt)表示zt的标准差;
(3.2)Z经过投影后的输出信号为sj,sj表示第j个慢特征序列;考虑线性条件下,表示系数向量,这等价于寻找一个从标准化输入信号Z中提取慢特征信号s=[s1 T,s2 T,…,sr T]T的转换矩阵即s=WZ;慢特征信号sj要满足的目标函数及约束条件为:
目标函数:
约束条件为:
(3.3)白化:利用奇异值分解,对输入数据的协方差矩阵〈ZZT〉进行白化处理可以去除数据中的相关性,使提取出的慢特征值携带不同的信息:
其中:Λ-1/2BT为白化矩阵,Ο为对应的白化后的输入信号;
W=PΛ-1/2BT (17)
(4)划分慢特征s:最慢的特征对应最小的特征值,将特征值由小到大排列,并依据特征值大小将前l个特征划分为s中变化较慢的特征,用sd表示;将后(r-l)个特征划分为s中变化较快的特征,用se表示;划分依据l的确定方法为,首先利用慢特征值s的变化快慢表示过程变量的变化快慢:
将提取出的慢特征值中缓慢度比输入数据缓慢度还要大的特征划分为快特征,一共有Me个这样的快特征:
这里card{·}表示集合{·}中元素个数;根据式(19)确定的Me值,对应将矩阵Ω也划分成两部分:
(5)计算动态监测指标:从典型变量空间的第一个样本点开始,每个样本点可以得到一组动态监测指标(Sd 2,Se 2);
(6)确定基于动态监测指标的控制限:利用核密度估计的方法,先估计出动态监测指标Sd 2的概率密度函数p(x),对于给定显著性水平α,Sd 2的控制限Sd 2 UCL的计算方式为:
以同样的方法可以计算出Se 2的控制限Se 2 UCL;
(7)对残差变量建立统计指标Qk:
其中:εi,k表示残差空间E中第i个变量在时刻k的取值,即E中第i行第k列的元素;
同样用核密度估计的方法计算Qk统计指标的控制限QUCL:
(8)在线监测控制性能:基于步骤(2)到(4)建立的CVA-SFA模型、步骤(5)到步骤(7)所得的三个监测统计量在线监测空气预热器控制系统的性能状态;该步骤由以下子步骤来实现:
(8.1)获取新在线数据以及新数据预处理:采集到新的一段观测数据后,其中,下标new表示新观测数据,首先按照步骤(2)将Ynew拓展成过去观测矩阵,并根据步骤(2)中获得的均值和标准差对过去观测矩阵进行标准化处理得到Ypnew;
(8.2)提取出新观测数据的典型变量和残差变量:标准化处理后,利用步骤(2)确定的转换矩阵Vr和L计算出新观测数据的典型变量空间Znew和残差空间Enew;
(8.3)提取新观测数据的典型变量空间Znew中的慢特征:首先按照步骤(3.1)中确定的均值和方差对Znew进行标准化处理,之后利用步骤(3.4)中确定的慢特征转换矩阵W,提取出标准化Znew的慢特征snew,并按照之前的划分参数将snew划分成sdnew和sfnew;
(8.5)在线判断空气预热器控制性能状态:实时比较三个监测指标与其各自的统计控制限,若三个监测指标都位于统计控制限之内,表明控制系统正常工作;若有一个或以上监测指标超出正常控制限,表明控制系统有异常状况发生。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811109196.2A CN109471420B (zh) | 2018-09-21 | 2018-09-21 | 基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811109196.2A CN109471420B (zh) | 2018-09-21 | 2018-09-21 | 基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109471420A CN109471420A (zh) | 2019-03-15 |
CN109471420B true CN109471420B (zh) | 2020-08-14 |
Family
ID=65663241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811109196.2A Active CN109471420B (zh) | 2018-09-21 | 2018-09-21 | 基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109471420B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110222393B (zh) * | 2019-05-28 | 2020-12-18 | 浙江大学 | 基于细粒度风电发电状态划分的风机叶片结冰异常监测方法 |
CN111949003B (zh) * | 2020-07-17 | 2021-09-03 | 浙江浙能技术研究院有限公司 | 一种基于SFA与Hellinger距离的闭环控制回路性能评价方法 |
CN112255978A (zh) * | 2020-10-11 | 2021-01-22 | 国网湖北省电力有限公司电力科学研究院 | 一种燃煤锅炉智慧运维系统及方法 |
CN113376551B (zh) * | 2021-06-23 | 2022-11-08 | 北京市腾河智慧能源科技有限公司 | 串联故障电弧的检测方法及系统、设备、存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902981A (zh) * | 2012-09-13 | 2013-01-30 | 中国科学院自动化研究所 | 基于慢特征分析的暴力视频检测方法 |
CN103970092A (zh) * | 2014-04-13 | 2014-08-06 | 北京工业大学 | 一种基于自适应fcm的多阶段发酵过程故障监测方法 |
CN104598681A (zh) * | 2015-01-14 | 2015-05-06 | 清华大学 | 基于缓慢特征分析的过程监控方法和系统 |
CN105320764A (zh) * | 2015-10-26 | 2016-02-10 | 天津大学 | 一种基于增量慢特征的3d模型检索方法及其检索装置 |
CN105955241A (zh) * | 2016-06-03 | 2016-09-21 | 北京科技大学 | 一种基于联合数据驱动生产过程的质量故障定位方法 |
CN106482826A (zh) * | 2016-09-28 | 2017-03-08 | 浙江大学 | 一种用于油气锅炉燃烧振动的多信号融合测量诊断方法 |
CN107122611A (zh) * | 2017-04-28 | 2017-09-01 | 中国石油大学(华东) | 青霉素发酵过程质量相关故障检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3367061B1 (en) * | 2017-02-28 | 2020-07-15 | Honda Research Institute Europe GmbH | Navigation system based on slow feature gradients |
-
2018
- 2018-09-21 CN CN201811109196.2A patent/CN109471420B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902981A (zh) * | 2012-09-13 | 2013-01-30 | 中国科学院自动化研究所 | 基于慢特征分析的暴力视频检测方法 |
CN103970092A (zh) * | 2014-04-13 | 2014-08-06 | 北京工业大学 | 一种基于自适应fcm的多阶段发酵过程故障监测方法 |
CN104598681A (zh) * | 2015-01-14 | 2015-05-06 | 清华大学 | 基于缓慢特征分析的过程监控方法和系统 |
CN105320764A (zh) * | 2015-10-26 | 2016-02-10 | 天津大学 | 一种基于增量慢特征的3d模型检索方法及其检索装置 |
CN105955241A (zh) * | 2016-06-03 | 2016-09-21 | 北京科技大学 | 一种基于联合数据驱动生产过程的质量故障定位方法 |
CN106482826A (zh) * | 2016-09-28 | 2017-03-08 | 浙江大学 | 一种用于油气锅炉燃烧振动的多信号融合测量诊断方法 |
CN107122611A (zh) * | 2017-04-28 | 2017-09-01 | 中国石油大学(华东) | 青霉素发酵过程质量相关故障检测方法 |
Non-Patent Citations (2)
Title |
---|
Control Performance Monitoring with Temporal Features and Dissimilarity Analysis for Nonstationary Dynamic Processes;Chunhui Zhao;《IFAC Papers Online》;20180727;第357-362页 * |
基于动态单类随机森林的非线性过程监控方法;曹玉苹等;《化工学报》;20170430(第04期);第1459-1465页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109471420A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109471420B (zh) | 基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 | |
CN109491358B (zh) | 一种面向百万千瓦超超临界机组锅炉动态信息的控制性能监测方法 | |
CN109538311B (zh) | 面向高端发电装备中汽轮机的控制性能实时监测方法 | |
CN108446529B (zh) | 基于广义互熵—dpca算法的有机朗肯循环系统故障检测方法 | |
Jiang et al. | Data reconciliation and gross error detection for operational data in power plants | |
CN111159844B (zh) | 一种电站燃气轮机排气温度的异常检测方法 | |
CN113011010A (zh) | 基于结构机理和运行数据的锅炉故障诊断方法及诊断系统 | |
CN109238760B (zh) | 基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法 | |
Zeqiu et al. | Modeling and optimization of the steam turbine network of an ethylene plant | |
CN104748807A (zh) | 一种基于流量校正的电站主蒸汽流量在线计算方法 | |
CN115730191A (zh) | 一种基于注意力机制的磨煤机故障预警方法 | |
Blanco et al. | New investigation on diagnosing steam production systems from multivariate time series applied to thermal power plants | |
CN110942258B (zh) | 一种性能驱动的工业过程异常监测方法 | |
CN103207567B (zh) | 一种低误报率的改进主元分析过程监测方法及其监测系统 | |
Lindner et al. | Data-driven fault detection with process topology for fault identification | |
CN115510904A (zh) | 基于时序预测的锅炉受热面积灰监测方法 | |
Wang et al. | Deep-learning modeling and control optimization framework for intelligent thermal power plants: A practice on superheated steam temperature | |
CN106896795A (zh) | 一种火电厂变负荷协调控制品质评价方法 | |
CN109270917B (zh) | 一种面向智能电厂汽轮机轴承的闭环控制系统故障退化状态预测方法 | |
CN111949003B (zh) | 一种基于SFA与Hellinger距离的闭环控制回路性能评价方法 | |
CN103697958A (zh) | 燃煤机组汽包出口饱和蒸汽质量流量的实时计量方法 | |
CN113298133B (zh) | 超临界机组锅炉爆管故障诊断方法 | |
Wang et al. | Power plant data filtering based on Gaussian naive Bayesian classification and prediction error method | |
CN105279553A (zh) | 一种高加给水系统故障程度识别方法 | |
CN115130564A (zh) | 一种基于并行gmm-lcksvd的高压加热器在线监测方法 |
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 |