飞灰含碳量测量方法
技术领域
本发明涉及热工控制领域,具体地,涉及一种飞灰含碳量测量方法。
背景技术
在燃煤电厂的实际生产过程中,由于技术或经济上的原因,许多重要技术参数和经济参数无法用常规的传感器直接测量,严重影响燃煤电厂经济运行。飞灰含碳量是火力发电厂燃煤锅炉燃烧效率的一个重要指标,它反映了燃煤机械未完全燃烧损失的大小。目前,飞灰含碳量在线检测方法有灼烧失重法和微波检测法。
灼烧法失重法测量技术是中国电力工业标准《飞灰和炉渣可然物测定方法》及《煤的工业分析方法》中的相关方法,当含有未燃尽碳的灰样在规定的高温下经灼烧后,由于灰样中残留的碳被燃尽后使灰样的质量出现了损失,利用灰样的烧失量作为依据计算出灰样中的含碳量,
含碳量的质量(%)=[灼烧前灰样加坩埚的质量(g)—灼烧后灰样加坩埚的质量(g)]/[灼烧前灰样加坩埚的质量(g)—收灰前坩埚的质量(g)]。
图1示出了采用灼烧失重法测量飞灰含碳量的装置结构框图。如图1所示,灼烧失重法的主要工作流程是将烟道中的灰样通过振动器和取样器收集到测量单元的坩埚中,由升降机构将坩埚放置于旋转托盘,再由测量单元内部的执行机构将装有灰样的坩埚送入灼烧装置(如,电炉加热器)进行高温灼烧,电子天平实时测量收灰前、收灰后及灼烧后的重量信号,控制单元对接收到的重量信号进行计算,获得飞灰的含碳量并在控制单元的显示屏上进行显示,其中可以使用手动控制盒控制控制单元进行计算或显示。灼烧后的灰样通过***的排灰机构和真空发生器排放回烟道中。
然而灼烧失重法在现场应用中存在下列一些问题:
(1)取样问题。装置取样速度慢,所以取样时间长,灰样冷却较快,容易堵灰。
(2)微波炉加热问题。首先不能达到实验室加热温度(电力工业标准《飞灰和炉渣可然物测定方法》中提出应在810±10℃下进行灼烧,慢灰要1小时,快灰要30分钟),灼烧温度、时间不能达到要求,并且由于微波加热特性和坩埚盛灰的最低要求制约使得灰样厚度超出范围,因而只能对坩埚内的灰样外表面完全灼烧,而内部为原灰样。其次,不能按规程先烘干水分(因受烟道中温差变化因素,飞灰有时会出现结露现象),在灼烧测挥发量,所以失重变化数据中有水分变化因素。
(3)称重问题。电子天平安装在现场测量柜中进行实时称重,锅炉运行导致电子天平一直振动,无法保证正常测量精度。
(4)机械问题。此类装置一般有非常复杂的机械装置和电机、齿冷等转动部件,使用中经常出现坩埚破碎、机械故障、电机故障及坩埚工位错误等各种机械故障,由于设备故障会涉及到复杂的机械装置、电气电子装置、气路、灰路***,即使有专人维护,也难以解决错综复杂的问题。
(5)时间滞后问题。装置取样速度慢,并需要机械运转、称重、灼烧,使得检测周期时间长。
所以,灼烧失重法的在线检测装置在现场使用中具有检测滞后时间长、灰样不能烧透、称重不准确、堵灰、机械故障率高的缺陷,维护要求高、难维护。
微波检测法的原理为:因为飞灰可燃物主要成份是碳及碳的介电常数,微波测试单元可利用固定频率发射能量衡定的微波信号,飞灰中可燃物的含量越高,吸收微波能量的作用就越强。
微波检测法的***可采用无动力飞灰取样器,自动将烟道中的灰样收集到微波测试装置的测量管中,由灰位控制器自动判别收集灰位的高度。当收集到足够的灰样时,***对飞灰含碳量进行微波谐振测量。控制装置打开电磁阀接入压缩空气吹扫已分析完的灰样,根据程序设定或手动设置,飞灰可以自动经采样管道吹回烟道或者送入收灰容器,以便于化学分析化验。
由于不同物质的频率特性不同、飞灰中物质组成复杂,且煤种变化时主要是矿物质变化使得微波检测法在煤种变化时不能检测飞灰含碳量变化。
发明内容
本发明的目的是提供一种飞灰含碳量测量方法,用于实现飞灰含碳量的在线测量,提高机组运行的安全性和经济性。
为了实现上述目的,本发明提供一种飞灰含碳量测量方法,该方法包括:根据飞灰含碳量的历史数据值及对应于所述历史数据值的辅助变量值建立历史数据样本集合;在所述历史数据样本集合中选取与当前工况相近的第一样本集;根据所述第一样本集中的辅助变量值和飞灰含碳量值拟合针对辅助变量的加权系数;以及根据当前工况的辅助变量值及所述加权系数计算所述当前工况的飞灰含碳量。
优选地,所述在所述历史数据样本集合中选取与当前工况相近的第一样本集包括:计算所述当前工况的辅助变量值与所述历史数据样本集合中每一样本的辅助变量值的相关系数;以及选取大于第一阈值的相关系数所对应的样本组成所述第一样本集。
优选地,采用皮尔逊积矩相关系数法计算所述相关系数。
优选地,所述方法还包括:对所述第一样本集进行主成分分析以得到包含主要辅助变量的第二样本集;以及根据所述第二样本集来拟合所述加权系数。
优选地,使用最小二乘法拟合所述加权系数。
优选地,所述辅助变量包括:每一给煤机的给煤量、烟气含氧量、排烟温度、主蒸汽压力、主蒸流量、主蒸汽温度、总风量、省煤器入口流量值、省煤器入口压力、省煤器出口压力、发电负荷、煤质特性、一次风总压力、二次风门开度、燃尽风门开度、二次风总压与炉膛差压、空气预热器出口温度、燃烧器摆角。
优选地,所述方法还包括:更新所述历史数据样本集合。
优选地,所述更新所述历史数据样本集合包括:在所述历史样本集合中确定多个样本中心;在新出现的样本与所述多个样本中心中每一者的距离大于第二预定值的情况下,将所述新出现的样本添加至所述历史数据样本集合。
优选地,所述更新所述历史数据样本集合进一步包括:将所述新出现的样本作为一个新的样本中心以更新所述多个样本中心;以及合并更新后的多个样本中心中彼此之间的距离小于第三阈值的两个样本中心。
优选地,所述更新所述历史数据样本集合进一步包括:在所述更新后的多个样本中心中彼此之间的距离均大于第三阈值的情况下,删除更新后的历史数据样本集合中与所述多个样本中心的每一者的距离最小者。
通过上述技术方案,基于飞灰含碳量的历史数据样本集合及当前工况的辅助变量值计算当前工况的飞灰含碳量,方案简单易行无需复杂的装置,并且能够提高机组运行的安全性和经济性。
本发明的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明,但并不构成对本发明的限制。在附图中:
图1示出了采用灼烧失重法测量飞灰含碳量的装置结构框图;
图2示出了一实施例中本发明所提供的飞灰含碳量测量方法的流程图;以及
图3示出了更新历史数据样本集合的流程图。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
图2示出了一实施例中本发明所提供的飞灰含碳量测量方法的流程图。如图2所示,本发明提供一种飞灰含碳量测量方法,该方法包括:根据飞灰含碳量的历史数据值及对应于所述历史数据值的辅助变量值建立历史数据样本集合(步骤S10);在所述历史数据样本集合中选取与当前工况相近的第一样本集(步骤S20);根据所述第一样本集中的辅助变量值和飞灰含碳量值拟合针对辅助变量的加权系数(步骤S30);以及根据当前工况的辅助变量值及所述加权系数计算所述当前工况的飞灰含碳量(步骤S40)。
以下将分步骤具体介绍本发明所提供的飞灰含碳量测量方法。
步骤S10:根据飞灰含碳量的历史数据值及对应于所述历史数据值的辅助变量值建立历史数据样本集合。
其中,飞灰含碳量与现场多个变量之间存在关系,为了反映提取飞灰含碳量的多方面信息,这里选取与飞灰含碳量相关性大的辅助变量。该辅助变量可以包括:每一给煤机的给煤量(假设有五个给煤机,则应当包括给煤机A给煤量、给煤机B给煤量、给煤机C给煤量、给煤机D给煤量、给煤机E给煤量)、烟气含氧量、排烟温度、主蒸汽压力、主蒸流量、主蒸汽温度、总风量、省煤器入口流量值、省煤器入口压力、省煤器出口压力、发电负荷、煤质特性、一次风总压力、二次风门开度、燃尽风门开度、二次风总压与炉膛差压、空气预热器出口温度、燃烧器摆角等。在所建立的历史样本集合中的每一样本中包含上述辅助变量的值及对应的飞灰含碳量的数据值。
步骤S20:在所述历史数据样本集合中选取与当前工况相近的第一样本集。
在该步骤中,可以计算当前工况的辅助变量值与历史数据样本集合中每一样本的辅助变量值的相关系数,相关系数越高,说明二者相关性越高,所对应的历史样本就越与当前工况相近。这里可以采用皮尔逊积矩相关系数法计算上述的相关系数。
设当前工况的辅助变量值表示为X1=(x1,x2……xn),其中x1,x2……xn表示当前工况中不同辅助变量的值,历史样本的辅助变量值为Y1=(y1,y2……yn),其中y1,y2……yn表示某一历史样本中不同辅助变量的值,则皮尔逊积矩相关系数r的计算公式为:
式中
Sxx为变量xi对其均值的偏差平方和,Syy为变量yi对其均值的偏差的平方和,Sxy为xi、yi的偏差平方和。r就称为向量X1对Y1(或向量Y1对X1)的单项相关系数,亦可称为相关系数。
相关系数r的取值范围是-1≤r≤1,相关系数r的绝对值越接近1,向量X1和Y1之间的相关程度越高,相关系数r的绝对值越接近0,向量X1和Y1之间的相关程度越低,在所计算出来的相关系数中选取大于第一阈值的相关系数所对应的样本组成上述的第一样本集。这里,第一阈值可以根据具体情况而定,并不作特定限制。
步骤30:根据第一样本集中的辅助变量值和飞灰含碳量值拟合针对辅助变量的加权系数。优选地,这里可以针对所获得的第一样本集进行主成分分析来得到包含主要辅助变量的第二样本集,根据该第二样本集中的主要辅助变量值和飞灰含碳量值拟合针对主要辅助变量的加权系数。
在一实施例中可以采用主成分分析法对所获得的第一样本集进行主成分分析以获得第二样本集。
主成分分析法是一种降维的统计方法,它借助于一个正交变换,将其分量相关的原随机向量转化成其分量不相关的新随机向量,这在代数上表现为将原随机向量的协方差阵变换成对角形阵,在几何上表现为将原坐标系变换成新的正交坐标系,使之指向样本点散布最开的p个正交方向,然后对多维变量***进行降维处理,使之能以一个较高的精度转换成低维变量***,再通过构造适当的价值函数,进一步把低维***转化成一维***。
假设X是一个n×m的数据矩阵,其中的每一列对应于一个变量,每一行对应于第一样本集中的一个样本。矩阵X可以分解为m个向量的外积之和,即
式中,ti∈Rn被称为得分向量,pi∈Rm称为负荷向量。X的得分向量也叫做X的主元。可以写成下列矩阵形式
X=TPT (3)
其中T=(t1,t2,…,tn)称为得分矩阵,P=(p1,p2,…,pm)称为负荷矩阵。
各个得分向量之间也是正交的,即对于任何i和j,当i≠j时,满足titj=0。各个负荷向量之间也是正交的,同时每个负荷向量的长度都为1,即
当矩阵X中的向量间存在一定程度的线性相关时,数据矩阵X的变化将主要体现在最前面的几个负荷向量的方向上,数据矩阵X在最后面的几个负荷向量上的投影将会很小,它们主要是由于测量噪声引起的。这样就可以讲矩阵X进行主元分解后写成下式
式中E为误差矩阵,代表X在pk+1到pm等负荷向量方向上的变化。由于误差矩阵E主要是由于测量噪声引起的,将E忽略掉往往会起到清除测量噪声的效果,而不会引起数据中有用信息的明显损失。因而数据X可以近似地表示为
在计算出式(7)后就完成针对所获得的第一样本集所进行的进行主成分分析。
因为X的前k个主代表了X数据中的绝大多数变化,因此可以用X的前k个主元来代替那些原始输入变量进行回归分析,这样便得到下面的主元回归模型
式中,矩阵Y代表经主成分分析后的第二样本集,B=(b1 b2 … bk)T为主元回归模型参数。可利用最小二乘方法通过下式计算而得到B
则输入变量的模型参数θ为
模型参数θ即为所计算出的加权系数。使用所计算的加权系数乘以每一加权系数所对应的当前工况的辅助变量值就可以计算出当前工况的飞灰含碳量。
进一步地,为了全面反映工况情况,在出现新的样本时,需要更新历史数据样本集合。首先在历史数据样本集合中确定l(l≥1)个样本中心,计算新出现的样本与l个样本中心中每一者的距离,当每个距离均大于阈值δ时,则将新出现的样本加入到历史数据样本集合中,以该新出现的样本作为一个新的样本中心,样本中心个数为l+1。计算l+1个样本中心彼此之间的距离,如果两个样本的距离小于阈值σ,则将该两个样本中心进行合并。如果两个样本的距离均大于阈值σ,计算历史数据库中每一样本与l+1个样本中心的距离,剔除距离最小的一个样本,从而保证历史数据库中的样本数保持恒定。
在实际使用中,本发明所提供的飞灰含碳量计算方法可以通过支持C语言编程的PLC可编程控制器(如西门子的S7系列的控制器)进行实现,通过以太网或硬接线的形式读取现场数据(辅助变量值),并将计算结果显示在运行人员的监控器中。
以上结合附图详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种简单变型,这些简单变型均属于本发明的保护范围。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合,为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。
此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明的思想,其同样应当视为本发明所公开的内容。