基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法技术

本发明专利技术公开了一种基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法。本发明专利技术针对磨煤机过程变量众多，控制器调节作用复杂的问题，运用典型变量分析提取过程的时序相关关系的信息，同时，结合慢特征分析算法，提取过程动静态特征，综合闭环系统调节作用与运行状态相关的时序相关关系和动静态特征，在不同的子空间内建立了动静态在线监测指标对磨煤机进行过程监测。该方法的优势在于能够充分反应闭环系统的调节作用，有利于对复杂过程特性的了解，增强了对磨煤机在线过程监测的可靠性和可信度，有助于火电厂对磨煤机运行状态做出准确判断，及时发现故障，从而保证了磨煤机的安全可靠运行。

【技术实现步骤摘要】
基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法
本专利技术属于火电过程过程监测领域，特别是涉及一种考虑闭环系统的调节作用，综合分析过程变量的时序相关关系以及过程的动静特征，对磨煤机的动静协同的在线过程监测方法。
技术介绍
我国电力工业快速发展，并且以主要火电为主，而煤是火力发电厂的主要燃料。据资料显示，我国每千瓦时平均供电耗煤比发达国家约高50g标准煤，这说明机组实际运行状态和设计状态有着极大的区别。2003年欧洲和北美发生的大规模停电事故提醒我们，现阶段的人类生活离不开电能，电能已成为人类生活的一部分，停电给国民经济造成的损失远远超过电力系统本身的损失，因此发电设备安全可靠的而运行显得尤为重要。随着我国电力工业的迅猛发展，火电发电技术的不断发展，发电机组的容量不断增加，复杂程度不断提高，需要处理判断的运行信息逐渐增多，因此火电厂一些关键设备的可靠性、可用性与安全性的问题日益突出。随着火力发电高参数、大容量机组的不断发展和应用，其相应的辅机设备也逐步趋向大型化，且在火电厂生产过程中的作用也愈发重要。根据近几年机组参与可靠性统计的数据显示，电力设备中的五大辅机，如磨煤机、给水泵、送风机、引风机、高压加热器，如果发生故障，将影响整个机组安全、经济运行，甚至还可能给现场人员带来安全危害，它的重要性甚至不亚于锅炉、汽轮机、发电机三大主机。为了发现和预防设备故障，保证设备安全、稳定的连续运行、减少事故的发生，保障运行安全、经济型。需要对设备运行状态进行监测分析。火电厂磨煤机作为锅炉燃烧制粉系统的核心设备，它是将煤块破碎并磨成煤粉的机械设备，其工作状况对整个电厂系统运行的安全和经济性有重要的影响。如磨煤机的碾磨压力、通风量、磨煤机中载荷等，它们的大小和异常都将影响生产效率。磨煤机通常必须在恶劣环境下工作，其工作过程是一个复杂的过程，且磨煤机发生故障的机率较高，而且一旦发生故障将会带来人生安全和经济效益的双重损失。例如2003年6月19日，甘肃火电工程公司平凉电厂磨煤机内部吊装分离器倒锥体(935.2kg)重心偏移，致1人死亡；2016年2月25日，大唐长山热电厂在1号炉C磨煤机发生一起热烟风灼烫事故，造成3人死亡。因此，因此确保火电厂磨煤机的安全可靠运行至关重要。针对磨煤机系统复杂、过程变量多样化的问题，前人对此已从不同的角度做了相应的研究与探讨，提出了相应的在线过程监测方法。总体来说，主要包括基于解析模型的方法和基于数据驱动的方法。基于模型的方法的不足主要有：基于模型的状态监测方法主要应用于线性时不变系统，然而大多数系统和设备都是非线性的、时变的.这种方法中对于传感器输出的信号变化是由于系统设备参数、结构的变化所致，还是由于传感器自身状态、参数的变化所致，很难辨别出来。基于数据驱动的流程工业过程监测研究与应用于上个世纪90年代逐渐兴起，由于工业自动化程度的加深和信息化水平的不断提高使得大量的过程数据得以存储并利用。因此，基于数据驱动的过程监测方法得到越来越多研究者的关注。然而，现有的数据分析方法，如主成分分析、偏最小二乘法、费舍尔判别分析，都仅仅考虑了过程的静态信息，且忽视了过程中闭环系统的调节作用对过程状态的影响。因此，将其运用在磨煤机上，并不能得到很好的监测效果。本专利技术的内容深入考虑了磨煤机实际过程运行中闭环系统的调节作用对于监测性能的影响，综合分析过程变量的时序相关关系以及过程的动静态特征，提出了一种新的面向磨煤机的动静特征协同分析的在线监测方法。
技术实现思路
本专利技术的目的在于针对现有的针对磨煤机的过程监测方法的不足，提出一种基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法，该方法充分考虑到闭环系统调节作用对过程特性的影响，主要体现在相关关系的变化以及过程的动态特征，首先运用典型变量分析提取过程的时序相关信息，同时，结合慢特征分析算法，提取过程动态信息，在不同的子空间内建立了动静态在线监测指标对磨煤机进行过程检测，实现对磨煤机运行状态的准确监测。本专利技术的目的通过以下技术方案实现：一种基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法，磨煤机是将煤块破碎并磨成粉末的机械，是火电厂中重要的辅机之一，磨煤机常见的故障有进出口漏风、出口温度高、磨煤机着火等，该方法包括以下步骤:(1)获取待分析数据：设一个磨煤机生产过程具有J个测量变量和操作变量，每一次采样可以得到一个1×J的向量，采样N次后得到的数据表述为一个二维矩阵X(N×J)，所述测量变量为运行过程中可被测量的过程变量，包括温度、转速、压力、阀门开度等。(2)基于典型变量分析的磨煤机运行数据时序相关信息提取，该步骤由以下子步骤来实现：(2.1)提取磨煤机运行数据的典型变量，分析时序相关关系：将每个时间点t处的测量值x(t)(1×J)分别用l个过去测量值x(t-1),x(t-2),…,x(t-l)，以及h个将来测量值x(t),x(t+1),…,x(t+h)的展开。其中，下标p代表过去，f代表将来。xp(t)表示过去测量值的集合，xf(t)表示将来测量值的集合；其中，l和h根据与当前时刻x(t)的相关程度确定，具体如下：与当前时刻t间隔l+1个采样间隔的过去时刻的测量值xp(t+l+1)和当前时刻x(t)的相关程度表示为：其中α是一个阈值，0<α<0.5。autocorr(Xi,p)是第i个过程变量加p个时滞后的自相关系数。和当前时刻的相关程度，当Al+1其小于α时，可认为xp(t+l+1)和x(t)不存在相关关系，不参与构建矩阵xp(t)。选择能满足上式的最小的值确定为l，同时，确定出h＝l，从而构建出xp(t)和xf(t)。当某测量值与当前时刻t的时间间隔大于所确定的l和h时，变量间的相关关系可以忽略。(2.2)将不同时刻展开的向量组成过去和将来的矩阵：Xp＝[xp(t+1),xp(t+2),…,xp(t+m)](4)Xf＝[xf(t+1),xf(t+2),…,xf(t+m)](5)其中下标p代表过去，f代表将来，m＝N-l-h+1。(2.3)求取不同时刻t的典型变量z(t)以及残差变量e(t)：其中，∑p,f表示过去矩阵Xp和将来矩阵Xf的协方差矩阵，∑p,p表示过去矩阵Xp的协方差矩阵，上标1/2表示对矩阵内元素进行开方操作，∑f,f表示将来矩阵Xf的协方差矩阵，上标1/2表示对矩阵内元素进行开方操作。表示式(6)等号左侧进行SVD分解，得到正交矩阵U和V，以及对角矩阵Λ，其对角线上的系数是主成分相关系数γ1≥…≥γr，下表r表示对角矩阵Λ对角线上元素个数。其中J代表求取相关变量的转换矩阵，L代表求取残差标量的转换矩阵，I表示单位矩阵，Uk包含了矩阵U的k列，k的选择根据下式决定：其中β是个阈值，0.5≤β≤1。Z＝JXP(10)E＝LXP(11)其中矩阵Z和矩阵E分别表示主成分空间和残差空间。Z和E的每列分别是典型变量z(t)和残差变量e(t)。(3)基于慢特征分析的对磨煤机运行的协同监测，该步骤由以下子步骤来实现：(3.1)分别对主成分空间Z和残差空间E进行SFA建模sc＝WcZ(12)se＝WeE(13)其中sc是提取出的主成分空间Z的慢特征，Wc本文档来自技高网...

【技术保护点】
1.一种基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法，其特征在于，该方法包括以下步骤：(1)获取待分析数据：设一个磨煤机生产过程具有J个测量变量和操作变量，每一次采样可以得到一个1×J的向量，采样N次后得到的数据表述为一个二维矩阵X(N×J)，所述测量变量为运行过程中可被测量的过程变量，包括温度、转速、压力、阀门开度等。(2)基于典型变量分析的磨煤机运行数据时序相关信息提取，该步骤由以下子步骤来实现：(2.1)提取磨煤机运行数据的典型变量，分析时序相关关系：将每个时间点t处的测量值x(t)(1×J)分别用l个过去测量值x(t‑1)，x(t‑2)，…，x(t‑l)，以及h个将来测量值x(t)，x(t+1)，…，x(t+h)的展开。

【技术特征摘要】
1.一种基于典型相关分析与慢特征分析的智能电厂燃煤发电机组磨煤机的在线监测方法，其特征在于，该方法包括以下步骤：(1)获取待分析数据：设一个磨煤机生产过程具有J个测量变量和操作变量，每一次采样可以得到一个1×J的向量，采样N次后得到的数据表述为一个二维矩阵X(N×J)，所述测量变量为运行过程中可被测量的过程变量，包括温度、转速、压力、阀门开度等。(2)基于典型变量分析的磨煤机运行数据时序相关信息提取，该步骤由以下子步骤来实现：(2.1)提取磨煤机运行数据的典型变量，分析时序相关关系：将每个时间点t处的测量值x(t)(1×J)分别用l个过去测量值x(t-1)，x(t-2)，…，x(t-l)，以及h个将来测量值x(t)，x(t+1)，…，x(t+h)的展开。其中，下标p代表过去，f代表将来。xp(t)表示过去测量值的集合，xf(t)表示将来测量值的集合；其中，l和h根据与当前时刻x(t)的相关程度确定，具体如下：与当前时刻t间隔l+1个采样间隔的过去时刻的测量值xp(t+l+1)和当前时刻x(t)的相关程度表示为：其中α是一个阈值，0＜α＜0.5。autocorr(Xi，p)是第i个过程变量加p个时滞后的自相关系数。和当前时刻的相关程度，当Al+1其小于α时，可认为xp(t+l+1)和x(t)不存在相关关系，不参与构建矩阵xp(t)。选择能满足上式的最小的值确定为l，同时，确定出h＝l，从而构建出xp(t)和xf(t)。当某测量值与当前时刻t的时间间隔大于所确定的l和h时，变量间的相关关系可以忽略。(2.2)将不同时刻展开的向量组成过去和将来的矩阵：Xp＝[xp(t+1)，xp(t+2)，…，xp(t+m)](4)Xf＝[xf(t+1)，xf(t+2)，…，xf(t+m)](5)其中下标p代表过去，f代表将来，m＝N-l-h+1。(2.3)求取不同时刻t的典型变量z(t)以及残差变量e(t)：其中，∑p，f表示过去矩阵Xp和将来矩阵Xf的协方差矩阵，∑p，p表示过去矩阵Xp的协方差矩阵，上标1/2表示对矩阵内元素进行开方操作，∑f，f表示将来矩阵Xf的协方差矩阵，上标1/2表示对矩阵内元素进行开方操作。表示式(6)等号左侧进行SVD分解，得到正交矩阵U和V，以及对角矩阵Λ，其对角线上的系数是主成分相关系数γ1≥…≥γr，下表r表示对角矩阵Λ对角线上元素个数。其中J代表求取相关变量的转换矩阵，L代表求取残差标量的转换矩阵，I表示单位矩阵，Uk包含了矩阵U的k列，k的选择根据下式决定：其中β是个阈值，0.5≤β≤1。Z＝JXP(10)E＝LXP(11)其中矩阵Z和矩阵E分别表示主成分空间和残差空间。Z和E的每列分别是典型变量z(t)和残差变量e(t)。(3)基于慢特征分析的对磨煤机运行的协同监测，该步骤由以下子步骤来实现：(3.1)分别对主成分空间Z和残差空间E进行SFA建模sc＝WcZ(12)se＝WeE(13)其中sc是提取出的主成分空间Z的慢特征，Wc是主成分空间Z的转换矩阵；se是提取出的残差空间E的慢特征，We是残差空间E的转换矩阵。(3.2)选取主成分空间中慢特征sc的主慢特征。变化缓慢的特征能够表示过程变化总体趋势，而变化较快的那些特征可以看作是噪声；将全部慢特征按照变化快慢进行排序，选择Rm(Rm＜Rc)个变化慢的慢特征作为主慢特征sc，d，剩余的慢特征可看作噪声，其中Rc为全部慢特征数；主慢特征数Rm的选择从重构的角度考虑，具体如下：过程变量xj可以通过慢特征s重构：其中，是的第j列，是将的某几个元素用0替代，元素个数为Rm，Rm的确定方法如下：降噪重构的过程变量的缓慢程度能够代表重构的过程变量保留了多少重要的信息；应含有尽可能少的变化快的噪声...

【专利技术属性】
技术研发人员：赵春晖，田畅，范海东，陈积明，孙优贤，李清毅，沙万里，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33