《基于代理模型模式匹配的锅炉实时操作优化》由会员分享,可在线阅读,更多相关《基于代理模型模式匹配的锅炉实时操作优化(15页珍藏版)》请在金锄头文库上搜索。
1、 基于代理模型模式匹配的锅炉实时操作优化 钟崴,林雪茹,林小杰,周懿(1.浙江大学 浙江省清洁能源与碳中和重点实验室,浙江 杭州 310027;2.浙江大学 能源工程学院,浙江 杭州 310027;3.浙江大学 嘉兴研究院,浙江 嘉兴 314024)发展风力、光伏可再生能源发电是推动我国能源改革的必要手段1-4.我国火电机组占主导地位5,因此研究如何优化火电机组的运行操作具有必要性与迫切性.燃煤锅炉是主要的耗能设备,如何优化其运行达到高效、低碳是核心问题.于晨龙等6-9对锅炉汽机的最佳运行组合及负荷优化问题展开了研究.在实际的执行环节,电厂锅炉的调控模式较粗犷,依赖人的经验展开优化方案的执行,
2、执行效果不佳或难以实施.系统复杂、工况多变,动态调整耗时大,易造成生产波动,降低经济效益.有必要将操作经验数字化,并基于模型进行决策,以辅助人员开展调控工作.目前,学术界展开以聚类法、回归模型法、在线监测法为代表的锅炉操作数据模型化研究.Song等10-12通过聚类的方法实现燃烧效率的最大化,但聚类方法的计算效率较低.吴家标等13-16提出回归模型法,解决了聚类法的计算效率难题,但难以适应工况变化和实际的操作模式.Nikula等17-18将在线监测法应用于锅炉性能优化,可以监测工况是否属于正常状态或优化状态,但未考虑参数的耦合性能,给出改进方法.随着深度学习的发展,Zou等19-21将深度强化
3、学习应用于锅炉的燃烧过程,可以考虑工况的时变和参数的耦合,但深度强化学习模型的可解释性和在连续过程的应用方面研究尚未成熟.现有的锅炉操作优化研究存在可执行性、工况适应性和计算效率不足的问题,需要进一步探索实时优化22场景下的方法.本文提出适用于多工况在线应用的代理模型模式匹配(pattern matching of agent model, PMAM)建模框架,以提升燃煤锅炉的效率,优化锅炉运行.主要贡献如下.1) 提出新的模式匹配方法,解决锅炉操作优化的可执行性问题.2) 提出基于模式匹配方法的优化操作库构建方法,解决锅炉操作优化的工况适应性问题.3) 提出基于神经网络算法的在线应用代理模型
4、,解决锅炉操作优化的计算效率问题.1 锅炉实时操作优化建模方法燃煤锅炉转化煤为蒸汽驱动发电,需要控制煤量、风量参数,保证安全并实现需要的蒸汽量.为了提高锅炉效率,需要根据输入的工况数据进行优化操作建议.采用反平衡方法23计算锅炉效率,输出给煤机、引风机、送风机的变频数据,以达到效率最高的目标.如图1所示为PMAM锅炉实时操作优化建模方法.该方法基于改进的模式匹配算法及神经网络算法,分为预建模、在线应用2个部分.预建模部分基于考虑滞后及时间窗的模式匹配优化模型,引入注意力机制参数、状态参数区间频率法、调控最小3层方案优化机制,得到各个工况的优化操作,构建工况参数与优化操作的最优代理模型.在线应用
5、部分基于工况参数与优化操作参数的最优代理模型,实现在线优化计算.图1 基于代理模型模式匹配的锅炉实时操作优化计算框架Fig.1 Boiler real-time operation optimization calculation framework based on pattern matching of agent model1.1 数据预处理1.1.1 滞后性计算 锅炉操作中的滞后性主要体现在给煤机的频率与锅炉主蒸汽流量的关系.基于相关性分析,开展滞后性计算.算法1:滞后性计算输入:样本个数为n,样本维度为1的样本Y1、Y2,最大滞后时长为tdelayMax输出:滞后时长tdelay1)
6、 将样本Y1、Y2归一化,得到Y1n 、Y2n.2) 构造数据集合dcorr ,第1列、第2列分别为Y1n、Y2n.3)foriinrange(-tdelayMax,tdelayMax):4) 删掉带有空值的行,得到修改后的数据集合dcorr.5) 计算dcorr的相关性系数,得到维度的相关性系数矩阵.6) 取dcorr 第1列第3行开始的数据,表征Y1与经过时间平移处理的Y2之间的相关性.相关性系数绝对值最大的时间移动值为tdelay ,即Y1滞后于Y2的时间为tdelay.设Y1、Y2均为时间长度为n的一维时间序列数据,研究Y1滞后于Y2的时间tdelay,则计算不同时间平移后Y2与Y1之
7、间的相关性,相关性越大,则认为平移的时间长度越优.数据归一化方法采用最大最小化方法,具体公式如下:数据集合dcorr的构造方法如图2所示.将Y2数据集进行时间轴平)移,得到图2的框选部分,时间长度为(n-tdelayMax的数据集即为dcorr.图2 滞后性计算中的数据重构示意图Fig.2 Schematic diagram of data reconstruction in hysteresis calculation1.1.2 基于主成分分析及k-均值改进的模糊C均值聚类模型 设置数据集X=XM,XO,XR,XCR.其中XM为工况表征参数,包括主蒸汽温度、压力与煤质数据等;XO为目标表征参
8、数;XR为关键的状态表征参数;XC为可调控参数;X的总时间维度为n0.根据滞后时间,对数据集X进行时间轴对齐重构,得到新的数据集合为X=XM,XO,XR,XCR,X的总时间维度为n.由于不同的参数数据范围不同,对数据集进行最大最小化的归一化处理.为了提高后续模式匹配优化的计算效率,处理后的数据使用主成分分析模型(principal component analysis, PCA)24降低维度.PCA降低数据维度的算法如下.算法2:PCA降低数据维度输入:样本个数为n,样本特征数为q的样本X,降维后数据特征占比为RlowV输出:低维样本DlowX1) 样本的每一列去均值化,得到DmeanMove
9、d.2) 求协方差矩阵DcovMat.3) 求协方差矩阵的特征值与特征向量.4) 将特征值从大到小排序,并依次取出,直到取出的特征值加和大于RlowV,取得对应的特征向量DredE.5) 计算降维后的数据DlowX=DmeanMovedDredE.6) 返回低维样本DlowX、协方差矩阵DcovMat.针对PCA降维后的数据,采用k-均值(kmeans)改进的模糊C均值聚类模型(fuzzy C-means, FCM)25进行工况聚类.记u=uijcn为隶属度矩阵,c为聚类数,聚类中心为v=v1,v2,,其中1ic,2cn,q为数据特征数.FCM目标函数如下.式中:m为模糊指数,控制每个距离平方
10、误差的占比,m1.采用k-means方法改进FCM方法中的初始模糊矩阵生成步骤,可以降低聚类结果的随机性,加快算法的收敛.改进FCM方法的计算步骤如下.1)初始化聚类参数:c、m、更新迭代阈值、迭代次数l、最大迭代次数lmax.其中,更新迭代阈值为2次迭代的模糊矩阵偏差绝对值容忍度.2)利用k-means聚类法生成初始的模糊矩阵U.算法3:k-means聚类生成初始的模糊矩阵输入:样本X、类别数目c输出:模糊矩阵U1) 初始化,令l= 0,随机选择c个样本点作为初始的聚类中心ml=m1l,mcl.2) 对样本进行聚类,计算每个样本到类中心的距离.将每个样本分配到与其最近的中心的所属类,生成聚类
11、结果cl.3) 计算聚类结果cl中每个类的样本均值,作为新的聚类中心4) 若符合停止条件(划分结果不再改变或得到迭代次数上限),则终止,输出c*=cl;否则,l=l+1,返回步骤2).5) 根据c*,计算U.6) 返回U.3)计算聚类中心:4)更新模糊矩阵:5)若Ul+1-Ul或llmax,则终止计算;否则l=l+1,跳转至步骤3).1.2 考虑时间窗口的模式匹配优化模型经过模型聚类后,得到分类后的数据集:为了避免测量误差和异常操作的影响,考虑时间维度的因素,将锅炉操作参数模式匹配优化视为带时间序列属性的问题,选取一段时间窗口的数据进行模式匹配优化计算.对于类别i第z个时间点,设时间窗长度为t
12、deta.步骤1)5)为工况注意力机制参数层,考虑工况参数相对于目标参数的重要性,筛选相似工况.步骤6)、7)为状态参数频率区间层,利用状态参数与锅炉效率的高效区域进行优化.步骤8)为调控最小层,考虑实际的调控习惯,使得建议的操作参数与当前操作参数之间的偏离较小.1)获取匹配分数下限值cthresholdLow、匹配分数上限值cthresholdUp、关键参数的优先级顺序corderR.2)根据工况表征参数XMi,对目标表征参数XOi计 算相关性系数,取绝对值得到注意力参数cak,i.3)筛选工况表征参数中变动偏差符合匹配分数上、下限的数据,得到更新后的XMinew.式中:Q为XMi的总时间维
13、度,G为XMi的总变量维度.4)计算z-tdeta+1,z时间段的工况表征参数与XMinew相似度得分,Q为XMinew的总时间维度,q为其中的一个时间点,1qQ:5)按相似度得分排序更新XMinew,得到对应时间的关键状态表征参数XRinew、目标表征参数XOinew 、操作参数XCinew.6)对于rlen(corderR),令r= 1: 取出第r个关键状态参数的数据条XRinew,:,r,计算最大值cmax与最小值cmin.根据区间单元长度cfreq,将XRinew,:,r划分为多个区间,统计每个区间的频率.将连续的区间合并为长度为clend的长区间,ceil为向上取整函数,round
14、为四舍五入取整函数.统计连续长区间的频率,按高到低排序,即为筛选后的关键状态参数的区间.根据目标参数的优化方向,取在该长区间的50%数据,得到筛选后的数据,令XRinew等于筛选后的数据;r=r+1.7)根据输出的关键状态表征参数XRinew,更新得到对应时间的工况表征参数XMinew、目标表征参数XOinew 、操作参数XCinew.8)计算z时间点的操作参数XCi,z与XCinew相似度得分,并按从高到低排序,输出得分最高的XCi,z为工况i第z个时间点的优化操作X*Ci,z.1.3 最优代理模型利用3层方案模式匹配机制,分别对历史生产数据进行优化,生成优化操作库.为了解决直接搜索优化操作
15、库数据量过大、搜索效率过低的问题,构建基于神经网络算法的最优代理模型.该代理模型以工况参数为输入,以优化后的操作、效率和关键参数为输出.式中:f为最优代理函数,为类别i第z个时间点的工况数据,X*Ci,z、X*Oi,z、X*Ri,z分别为类别i第z个时间点的优化可操作参数、目标参数、关键状态参数.在线应用时,采用孤立森林26(isolation forest)法进行偏离检测.若当前工况非偏离点,则调用式(13)计算优化的操作参数;否则更新预建模的聚类中心与最优代理模型,输出当前的操作参数为优化方案.2 案例选取以浙江省某热电公司为案例进行研究.该公司的机组规模为4台蒸汽质量流量为220 t/h的高温超高压循环流化床锅炉.选取1#锅炉数据展开计算实验,并以2#锅炉数据开展验证.选择的工况参数、目标参数、关键状态参数、操作参数如表1所示.表1 锅炉操作优化计算的相关参数表Tab.1 Table of parameters related to boiler operation optimization calculation选择历史数据时间段为2022-02-16 22:412022-04-29 21:29,时间间隔为1 min,共103 586条数据构成原