CN109538311B - 面向高端发电装备中汽轮机的控制性能实时监测方法 - Google Patents
面向高端发电装备中汽轮机的控制性能实时监测方法 Download PDFInfo
- Publication number
- CN109538311B CN109538311B CN201811110373.9A CN201811110373A CN109538311B CN 109538311 B CN109538311 B CN 109538311B CN 201811110373 A CN201811110373 A CN 201811110373A CN 109538311 B CN109538311 B CN 109538311B
- Authority
- CN
- China
- Prior art keywords
- matrix
- new
- slow
- monitoring
- steam turbine
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000010248 power generation Methods 0.000 title claims abstract description 14
- 230000008859 change Effects 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 95
- 239000013598 vector Substances 0.000 claims description 39
- 230000009466 transformation Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 10
- 230000002087 whitening effect Effects 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 239000003245 coal Substances 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 229910000897 Babbitt (metal) Inorganic materials 0.000 claims description 3
- 238000005311 autocorrelation function Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 238000004886 process control Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 description 14
- 230000008901 benefit Effects 0.000 description 5
- 239000000284 extract Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004821 distillation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009434 installation 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
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/003—Arrangements for testing or measuring
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明公开了一种面向高端发电装备中汽轮机的控制性能实时监测方法。汽轮机性能指标与火电厂机组的负荷、运行参数之间存在着复杂的关系,本发明针对高端发电装备中汽轮机因参数众多、工况多变导致的控制性能监测困难的问题,运用典型变量分析提取汽轮机控制***变量间的相关信息,再利用慢特征分析算法,提取相关信息中的动态信息。最后,结合变量的相关性和变化快慢信息构造汽轮机控制性能在线监测模型。该方法克服了大型汽轮机因变量众多、工况变化而导致的控制性能监测困难的问题,大大提高了动态过程控制性能在线监测的准确度,有助于火电厂对汽轮机控制***进行有效及时的监测,对保证高端发电装备的安全可靠运行意义重大。
Description
技术领域
本发明属于火电过程控制***性能监测领域,特别是涉及高端发电装备中汽轮机运行相关性信息与动态信息的在线性能监测方法。
背景技术
控制***在现代化的工业过程中占据非常重要的地位,生产质量、操作安全、物能消耗等影响经济效益的指标都直接或间接的与控制***的性能有关。在实际生产过程中,控制***在投入使用初期往往性能表现良好,但运行一段时间后,由于设备的磨损、定期保养和维护不及时等原因,可能导致控制***的性能下降,控制性能变差会直接影响生产质量,导致经济效益亏损,若因此引发生产故障,还会涉及到人的生命安全甚至社会企业的财产安全,带来极大威胁。Torrres等人对2004-2005年巴西12家工厂(石化、造纸、水泥、钢铁、采矿等),超过700个控制回路进行检验,结果显示14%回路的阀门磨损过度,15%的阀门存在迟滞问题,16%的回路存在严重的整定问题,24%的控制器输出存在饱和现象,41%的回路因为整定问题、耦合、扰动以及执行器的问题而存在振荡现象。
另外,实际生产中,一个生产过程可能会有数以千计的控制回路共同作用,Eastman化学公司中的两个精馏生产设备拥有多大14000个控制回路,在HVAC生产过程中,其控制回路的数量甚至能够达到十万个。高端发电装备具有较高的复杂性,具体体现在规模庞大、设备众多、参数多样化且相互影响等方面。此外,大规模的发电机组,现场具有高温、高压及高噪声等特点。
控制性能评价与监测技术是过程控制领域新兴的一项重要技术,它能够利用设备的日常运行数据,实时监测监视***控制性能的变化,对控制***的问题做出早期识别和优化。对于发电机组,由于电力***中的用电负荷是经常变化的,为了维持有功功率平衡,保持***频率稳定,需要发电部门相应的改变发电机的出力以适应用电负荷的变化,即发电机组的工况不是稳定不变的。但是现有的控制性能评价与监测方法如主成分分析、偏最小二乘法、费舍尔判别分析,都是基于工况稳定的理想假设下进行的,因此,将其运用在大型火力发电机组汽轮机控制***性能监测上,并不能得到很好的监测效果。
发明内容
汽轮机是以连续流动的蒸汽作为工质,将蒸汽的热能转换为机械能的一种旋转式热力发动机。一般由喷管(或静叶片)、动叶片、叶轮、机轴轴承和汽缸等组成。汽轮机具有转速高、运行平稳可靠、单机功率可制成很大(已超过1000MW)和便于与发电机直接联接等优点,但整机构造复杂,对设计、制造、安装、运行和检修技术均有很严格的要求,且必需配置相应蒸汽参数和容量的锅炉。除已广泛用作现代中、大型火力发电厂、核电站和大型舰船的主要发动机外,还可用于大型化工或钢铁等联合企业。
本发明的目的在于针对大型火力发电机组汽轮机因参数众多、结构复杂、工况多变导致的控制性能监测困难的问题,运用典型变量分析与慢特征分析融合算法提取汽轮机控制***变量间的相关信息与变化快慢信息,克服了大型汽轮机因变量众多、工况变化而导致的控制性能监测困难的问题。
本发明的目的是通过以下技术方案实现:高端发电装备中汽轮机的控制性能实时监测方法,该方法包括以下步骤:
(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)
D=diag(γ1,γ2,…,γJp)
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和残差空间Ε中的慢特征sZ,sE。以提取典型变量空间Z中慢特征sZ为例,该方法主要步骤如下:
(3.1)数据标准化:对典型变量空间Z按变量进行标准化处理,计算公式如下:
zt表同一变量在不同时刻的时序向量,mean(zt)表示zt的均值,std(zt)表示zt的标准差。
(3.2)Z经过投影后的输出信号为sZj,sZj表示sZ第j个慢特征序列。考虑线性条件下,表示系数向量,这等价于寻找一个从标准化输入信号Z中提取慢特征信号sZ=[sZ1 T,sZ2 T,...,sZr T]T的转换矩阵即sZ=WZZ。慢特征信号sZj要满足的目标函数及约束条件为:
目标函数:
约束条件为:
(3.3)白化:利用奇异值分解,对输入数据的协方差矩阵<ZZT>进行白化处理可以去除数据中的相关性,使提取出的慢特征值携带不同的信息:
其中:ΛZ -1/2BT为白化矩阵,ΟZ为对应的白化后的输入信号。
WZ=PΛZ -1/2BT (17)
所述残差空间Ε中的慢特征sE的提取方法与上述典型变量空间Z中慢特征sZ的提取方法相同。
(4)划分慢特征sZ:最慢的特征对应最小的特征值,将特征值由小到大排列,并依据特征值大小将前l个特征划分为sZ中变化较慢的特征,用sZ,d表示;将后(r-l)个特征划分为sZ中变化较快的特征,用sZ,e表示。划分依据l的确定方法为,首先利用慢特征值sZ的变化快慢表示过程变量的变化快慢:
将提取出的慢特征值中缓慢度比输入数据缓慢度还要大的特征划分为快特征,一共有Me个这样的快特征:
这里card{·}表示集合{·}中元素个数。根据式(19)确定的Me值,对应将矩阵ΩZ也划分成两部分:
(5)计算动态监测指标:从典型变量空间的第一个样本点开始,每个样本点可以得到一组动态监测指标(SZ,d 2,SZ,e 2)。
(7)按照步骤(3)到步骤(6)所述方法,提取残差空间Ε的慢特征sE并将sE划分成两部分sE,d,sE,e,建立监测指标SE,d 2,SE,e 2并计算控制限 与对典型变量空间Z的处理方式相同,不再赘述。
(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)中确定的慢特征转换矩阵WZ,提取出标准化Znew的慢特征sZnew,并按照之前的划分参数将sZnew划分成sZ,d new和sZ,e new,同样根据WE可以得到Enew,进一步得到sE,d new和sE,e new。
(8.5)在线判断汽轮机控制性能状态:实时比较四个监测指标与其各自的统计控制限,若四个监测指标都位于统计控制限之内,表明控制***正常工作;若有一个或以上监测指标超出正常控制限,表明控制***有异常状况发生。
本发明的有益效果在于:本发明针对大型火力发电机组汽轮机因参数众多、结构复杂、工况多变导致的控制性能监测困难的问题,提出了一种面向高端发电装备中汽轮机的控制性能实时监测方法,该方法运用典型变量分析提取汽轮机控制***变量间的相关信息,再利用慢特征分析算法,提取相关信息和残差信息中的动态信息。最后,结合变量的相关性和变化快慢信息构造汽轮机控制性能在线监测模型。该方法克服了大型汽轮机因变量众多、工况变化而导致的控制性能监测困难的问题,大大提高了动态过程控制性能在线监测的准确度,有助于火电厂对汽轮机控制***进行有效及时的监测,有助于保证大型发电机组的安全可靠运行,同时达到提高其生产效益的生产需求。
附图说明
图1是本发明面向高端发电装备中汽轮机的控制性能实时监测方法的流程图,(a)为离线建模过程流程图,(b)为在线监测过程流程图;
图2是本发明CVA-SFA方法用于统计过程监测的结果图,(a)为正常情况下监测结果图,(b)为异常情况下监测结果图。
具体实施方式
下面结合附图及具体实例,对本发明作进一步详细说明。
本发明以浙能集团下属嘉华电厂#5号机组汽轮机为例,该机组的功率为60万千瓦,为大型火力发电机组,包括60个过程变量,包括轴承金属温度、轴承轴向振动、发电机有功功率、汽机转速、励磁三相温度、凝结水流量(三分钟平均)、给水压力、给水流量、给煤机给煤量等,以及一些阀门开度。
应该理解,本发明不止局限于上述实例的火电发电过程,凡是熟悉本领域的技术人员在不违背本发明的前提下还可以做出等同变型或替换,这些等同的变型或替换均包含在本申请权利要求所限定的范围内。
如图1所示,本发明是一种面向高端发电装备中汽轮机的控制性能实时监测方法,包括以下步骤:
(1)获取训练数据:设汽轮机的控制***具有J个测量变量和操作变量,每一次采样可以得到一个J×1的观测向量yk,其中下标k为时间指标,采样N次后得到的数据表述为一个二维观测矩阵本实例中,采样周期为10分钟,共4566个样本,60个过程变量,所测变量为汽轮机运行过程中的流量、振动、温度、压力、阀门开度等;
(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)
D=diag(γ1,γ2,…,γJp)
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和残差空间Ε中的慢特征sZ,sE。以提取典型变量空间Z中慢特征sZ为例,该方法主要步骤如下:
(3.1)数据标准化:对典型变量空间Z按变量进行标准化处理,计算公式如下:
zt表同一变量在不同时刻的时序向量,mean(zt)表示zt的均值,std(zt)表示zt的标准差。
(3.2)Z经过投影后的输出信号为sZj,sZj表示sZ第j个慢特征序列。考虑线性条件下,表示系数向量,这等价于寻找一个从标准化输入信号Z中提取慢特征信号sZ=[sZ1 T,sZ2 T,…,sZr T]T的转换矩阵即sZ=WZZ。慢特征信号sZj要满足的目标函数及约束条件为:
目标函数:
约束条件为:
(3.3)白化:利用奇异值分解,对输入数据的协方差矩阵<ZZT>进行白化处理可以去除数据中的相关性,使提取出的慢特征值携带不同的信息:
其中:ΛZ -1/2BT为白化矩阵,ΟZ为对应的白化后的输入信号。
WZ=PΛZ -1/2BT (17)
所述残差空间Ε中的慢特征sE的提取方法与上述典型变量空间Z中慢特征sZ的提取方法相同。
(4)划分慢特征sZ:最慢的特征对应最小的特征值,将特征值由小到大排列,并依据特征值大小将前l个特征划分为sZ中变化较慢的特征,用sZ,d表示;将后(r-l)个特征划分为sZ中变化较快的特征,用sZ,e表示。划分依据l的确定方法为,首先利用慢特征值sZ的变化快慢表示过程变量的变化快慢:
将提取出的慢特征值中缓慢度比输入数据缓慢度还要大的特征划分为快特征,一共有Me个这样的快特征:
这里card{·}表示集合{·}中元素个数。根据式(19)确定的Me值,对应将矩阵ΩZ也划分成两部分:
(5)计算动态监测指标:从典型变量空间的第一个样本点开始,每个样本点可以得到一组动态监测指标(SZ,d 2,SZ,e 2)。
(7)按照步骤(3)到步骤(6)所述方法,提取残差空间Ε的慢特征sE并将sE划分成两部分sE,d,sE,e,建立监测指标SE,d 2,SE,e 2并计算控制限 与对典型变量空间Z的处理方式相同,不再赘述。
(8)在线监测控制性能:基于步骤(2)到(4)建立的CVA-SFA模型、步骤(5)到步骤(7)所得的四个监测统计量可以在线监测汽轮机控制***的性能状态。该步骤由以下子步骤来实现:
(8.1)获取新在线数据以及新数据预处理:采集到新的一段观测数据后,其中,下表new表示新观测数据,首先按照步骤(2)将Ynew拓展成过去矩阵,并根据步骤(2)中获得的均值和标准差对过去矩阵进行标准化处理得到Ypnew。本实例中,新数据共有两份,数据一为正常工况下采集的数据,采样周期为10分钟,共2271个样本,60个过程变量,数据二为发生异常工况下记录的数据,采样周期为10分钟,共2563个样本,60个过程变量,所测变量为汽轮机运行过程中的流量、振动、温度、转速、电流等;
(8.2)提取出新观测数据的典型变量和残差变量:标准化处理后,利用步骤(2)确定的转换矩阵Vr和L计算出新观测数据的典型变量空间Znew和残差空间Enew。
(8.3)提取新观测数据的典型变量空间Znew中的慢特征:首先按照步骤(3.1)中确定的均值和方差对Znew进行标准化处理,之后利用步骤(3.4)中确定的慢特征转换矩阵WZ,提取出标准化Znew的慢特征sZnew,并按照之前的划分参数将sZnew划分成sZ,d new和sZ,e new,同样根据WE可以得到Enew,进一步得到sE,d new和sE,e new。
(8.5)在线判断汽轮机控制性能状态:实时比较四个监测指标与其各自的统计控制限,若四个监测指标都位于统计控制限之内,表明控制***正常工作;若有一个或以上监测指标超出正常控制限,表明控制***有异常状况发生。图2(a)中,四组统计量及对应控制线中,仅个别点的统计量超过了控制线,在置信水平α=0.05的条件下,可以认为新工况数据是正常的,即控制***表现正常;图2(b)中,四组统计量SZ,d 2,SZ,e 2,SE,d 2,SE,e 2均在第1330–1430、1808–1825、1943–1972、2093–2128四段明显超出阈值线,统计量Qk在第326个采样点左右第一次明显超限后一直维持超限状态,据此可以判断在这些时段内控制***有异常发生,这时就可以采用适当的故障诊断方法,比如贡献图方法分析隔离出可能的故障变量。
本发明运用典型变量分析提取汽轮机控制***变量间的相关信息,再利用慢特征分析算法,提取相关信息中的动态特征,由该方法提取的特征可以反应控制器的调节作用。最后,结合变量的相关性和变化快慢信息构造汽轮机控制性能在线监测模型,该方法克服了大型汽轮机因变量众多、工况变化而导致的控制性能监测困难的问题,大大提高了动态过程控制性能在线监测的准确度,有助于火电厂对汽轮机***进行有效及时的监测,有助于保证大型火力发电机组的安全可靠运行,同时达到提高其生产效益的生产需求。
Claims (1)
1.一种面向高端发电装备中汽轮机的控制性能实时监测方法,其特征在于,该方法包括以下步骤:
(1)获取训练数据:设汽轮机的控制***具有J个测量变量和操作变量,每一次采样得到一个J×1的观测向量yk,其中下标k为时间指标,采样N次后得到的数据表述为一个二维观测矩阵所述测量变量为汽轮机运行过程中可被测量的状态参数,包括轴承金属温度、轴承轴向振动、发电机有功功率、汽机转速和励磁三相温度;所述操作变量包括凝结水流量、给水压力、给水流量和给煤机给煤量;训练数据应当选取汽轮机在正常运行状态下的采样数据;
(2)利用CVA算法提取数据的时序相关信息,该步骤通过以下子步骤实现:
分别用所有的过去观测向量和将来观测向量构建过去观测矩阵Yp和将来观测矩阵Yf:
其中,M=N-f-p+1,p,f为两类时滞参数,令p=f,其值通过样本自相关函数Ap+1来确定:
其中: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)利用慢特征分析算法分别提取典型变量空间Z和残差空间Ε中的慢特征sZ,sE,典型变量空间Z中慢特征sZ的提取方法如下:
(3.1)数据标准化:对典型变量空间Z按变量进行标准化处理,计算公式如下:
zt表同一变量在不同时刻的时序向量,mean(zt)表示zt的均值,std(zt)表示zt的标准差;
(3.2)Z经过投影后的输出信号为sZj,sZj表示sZ第j个慢特征序列;考虑线性条件下,表示系数向量,这等价于寻找一个从标准化输入信号Z中提取慢特征信号sZ=psZ1 T,sZ2 T,…,sZr T]T的转换矩阵即sZ=WZZ;慢特征信号sZj要满足的目标函数及约束条件为:
目标函数:
约束条件为:
(3.3)白化:利用奇异值分解,对输入数据的协方差矩阵<ZZT>进行白化处理去除数据中的相关性,使提取出的慢特征值携带不同的信息:
其中:ΛZ -1/2BT为白化矩阵,ΟZ为对应的白化后的输入信号;
WZ=PΛZ -1/2BT (17)
所述残差空间Ε中的慢特征sE的提取方法与上述典型变量空间Z中慢特征sZ的提取方法相同;
(4)划分慢特征sZ:最慢的特征对应最小的特征值,将特征值由小到大排列,并依据特征值大小将前l个特征划分为sZ中变化较慢的特征,用sZ,d表示;将后(r-l)个特征划分为sZ中变化较快的特征,用sZ,e表示;划分依据l的确定方法为,首先利用慢特征值sZ的变化快慢表示过程变量的变化快慢:
将提取出的慢特征值中缓慢度比输入数据缓慢度还要大的特征划分为快特征,一共有Me个这样的快特征:
这里card{·}表示集合{·}中元素个数;根据式(19)确定的Me值,对应将矩阵ΩZ也划分成两部分:
(5)计算动态监测指标:从典型变量空间的第一个样本点开始,每个样本点得到一组动态监测指标SZ,d 2,SZ,e 2;
(8)在线监测控制性能:基于步骤(2)到(4)建立的公式、步骤(5)到步骤(7)所得的四个监测统计量在线监测汽轮机控制***的性能状态,该步骤由以下子步骤来实现:
(8.1)获取新在线数据以及新数据预处理:采集到新的一段观测数据后,其中,下标new表示新观测数据,首先按照步骤(2)将Ynew拓展成过去矩阵,并根据步骤(2)中获得的均值和标准差对过去矩阵进行标准化处理得到Ypnew;
(8.2)提取出新观测数据的典型变量和残差变量:标准化处理后,利用步骤(2)确定的转换矩阵Vr和L计算出新观测数据的典型变量空间Znew和残差空间Enew;
(8.3)提取新观测数据的典型变量空间Znew中的慢特征:首先按照步骤(3.1)中确定的均值和方差对Znew进行标准化处理,之后利用步骤(3.4)中确定的慢特征转换矩阵WZ,提取出标准化Znew的慢特征sZnew,并按照之前的划分参数将sZnew划分成sZ,d new和sZ,e new,同样根据WE得到Enew,进一步得到sE,d new和sE,e new;
(8.5)在线判断汽轮机控制性能状态:实时比较四个监测指标与其各自的统计控制限,若四个监测指标都位于统计控制限之内,表明控制***正常工作;若有一个或以上监测指标超出正常控制限,表明控制***有异常状况发生。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811110373.9A CN109538311B (zh) | 2018-09-21 | 2018-09-21 | 面向高端发电装备中汽轮机的控制性能实时监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811110373.9A CN109538311B (zh) | 2018-09-21 | 2018-09-21 | 面向高端发电装备中汽轮机的控制性能实时监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109538311A CN109538311A (zh) | 2019-03-29 |
CN109538311B true CN109538311B (zh) | 2020-08-04 |
Family
ID=65843543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811110373.9A Active CN109538311B (zh) | 2018-09-21 | 2018-09-21 | 面向高端发电装备中汽轮机的控制性能实时监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109538311B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111188761A (zh) * | 2019-12-31 | 2020-05-22 | 杭州哲达科技股份有限公司 | 一种基于Fourier-CVA模型面向机泵设备的监测方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111199246B (zh) * | 2019-12-24 | 2024-04-05 | 泉州装备制造研究所 | 一种工况分类方法 |
CN111474920B (zh) * | 2020-04-29 | 2021-04-06 | 电子科技大学 | 一种基于ica与acva的故障检测方法 |
CN113065583B (zh) * | 2021-03-13 | 2023-11-14 | 宁波大学科学技术学院 | 一种基于在线非线性判别特征分析的精馏过程异常监测方法 |
CN114396317B (zh) * | 2021-12-01 | 2022-12-16 | 上海发电设备成套设计研究院有限责任公司 | 核电汽轮机多目标多维度在线联合监控方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1494020A1 (en) * | 2003-06-30 | 2005-01-05 | Siemens Westinghouse Power Corporation | Method and apparatus for measuring on line failure of turbine thermal barrier coatings |
US20160033580A1 (en) * | 2012-05-29 | 2016-02-04 | Board Of Regents Of The University Of Nebraska | Detecting Faults in Turbine Generators |
CN106097151A (zh) * | 2016-06-24 | 2016-11-09 | 清华大学 | 一种基于数据协调的降低电厂数据不确定度的方法 |
CN106680012A (zh) * | 2017-01-25 | 2017-05-17 | 浙江大学 | 一种面向大型燃煤发电机组非平稳过程的故障检测方法和诊断方法 |
CN107766457A (zh) * | 2017-09-27 | 2018-03-06 | 华能国际电力股份有限公司上海石洞口第二电厂 | 一种火电厂运营管理***及其任务执行方法 |
US20180171774A1 (en) * | 2016-12-21 | 2018-06-21 | Schlumberger Technology Corporation | Drillstring sticking management framework |
CN108334674A (zh) * | 2018-01-17 | 2018-07-27 | 浙江大学 | 一种基于参数关联性智能分析的汽轮机高压缸运行状态监测方法 |
-
2018
- 2018-09-21 CN CN201811110373.9A patent/CN109538311B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1494020A1 (en) * | 2003-06-30 | 2005-01-05 | Siemens Westinghouse Power Corporation | Method and apparatus for measuring on line failure of turbine thermal barrier coatings |
US20160033580A1 (en) * | 2012-05-29 | 2016-02-04 | Board Of Regents Of The University Of Nebraska | Detecting Faults in Turbine Generators |
CN106097151A (zh) * | 2016-06-24 | 2016-11-09 | 清华大学 | 一种基于数据协调的降低电厂数据不确定度的方法 |
US20180171774A1 (en) * | 2016-12-21 | 2018-06-21 | Schlumberger Technology Corporation | Drillstring sticking management framework |
CN106680012A (zh) * | 2017-01-25 | 2017-05-17 | 浙江大学 | 一种面向大型燃煤发电机组非平稳过程的故障检测方法和诊断方法 |
CN107766457A (zh) * | 2017-09-27 | 2018-03-06 | 华能国际电力股份有限公司上海石洞口第二电厂 | 一种火电厂运营管理***及其任务执行方法 |
CN108334674A (zh) * | 2018-01-17 | 2018-07-27 | 浙江大学 | 一种基于参数关联性智能分析的汽轮机高压缸运行状态监测方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111188761A (zh) * | 2019-12-31 | 2020-05-22 | 杭州哲达科技股份有限公司 | 一种基于Fourier-CVA模型面向机泵设备的监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109538311A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109538311B (zh) | 面向高端发电装备中汽轮机的控制性能实时监测方法 | |
CN109491358B (zh) | 一种面向百万千瓦超超临界机组锅炉动态信息的控制性能监测方法 | |
CN109471420B (zh) | 基于cva-sfa的智能电厂大型燃煤发电机组空气预热器控制性能监测方法 | |
CN110262450B (zh) | 面向汽轮机的多种故障特性协同分析的故障预测方法 | |
CN112179655B (zh) | 一种基于阈值分级的汽轮发电机故障预警方法 | |
CN108446529B (zh) | 基于广义互熵—dpca算法的有机朗肯循环***故障检测方法 | |
US20040215353A1 (en) | Process and device for evaluating the performance of a process control system | |
CN111159844B (zh) | 一种电站燃气轮机排气温度的异常检测方法 | |
CN112836941B (zh) | 一种煤电机组汽轮机高加***在线健康状况评估方法 | |
CHEN et al. | Fault diagnosis in chemical process based on self-organizing map integrated with fisher discriminant analysis | |
CN115730191A (zh) | 一种基于注意力机制的磨煤机故障预警方法 | |
CN109184821B (zh) | 一种大型发电机组汽轮机的闭环信息分析的在线监测方法 | |
CN102136038A (zh) | 基于状态空间的多特征设备状态评价方法及应用 | |
CN114116832A (zh) | 一种基于数据驱动的配电网异常识别方法 | |
CN109270917B (zh) | 一种面向智能电厂汽轮机轴承的闭环控制***故障退化状态预测方法 | |
Xinyi et al. | Fault diagnosis in chemical process based on self-organizing map integrated with fisher discriminant analysis | |
CN111949003B (zh) | 一种基于SFA与Hellinger距离的闭环控制回路性能评价方法 | |
CN105279553A (zh) | 一种高加给水***故障程度识别方法 | |
CN115130564A (zh) | 一种基于并行gmm-lcksvd的高压加热器在线监测方法 | |
Wang et al. | Power plant data filtering based on Gaussian naive Bayesian classification and prediction error method | |
CN111695581B (zh) | 一种基于自回归极限学习机的非线性动态过程监测方法 | |
Li et al. | Physics-guided deep learning method for tool condition monitoring in smart machining system | |
CN113298133A (zh) | 超临界机组锅炉爆管故障诊断方法 | |
Ma et al. | Process monitoring of the pneumatic control valve using canonical variate analysis | |
CN114819743B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |