一种基于动态运行数据分析的燃烧优化控制方法

技术领域

本发明涉及相关性分析和燃烧优化控制技术领域，具体涉及一种基于动态运行数据分析的燃烧优化控制方法。

背景技术

自从“智能电厂”成为发电行业追逐的热点，传统火电行业智能发电技术研发逐渐从面上的先进信息技术的应用扩展到底层核心难点技术的研发层次。发电厂DCS系统所存储的运行数据则是智能发电技术研发的重要基础和条件。可以在一定规模发电集团的技术研究院建设“智能发电技术中心”,集中进行智能化技术的研发；各个发电厂和发电机组向该中心传送实时数据和信息,并从该中心获得智能控制、智能运维、智能管理等方面的技术支持。但到目前为止，对机组运行历史数据的智能化运用的理论和技术均不充分。近年来研究人员还提出了将整个锅炉工况划分为不同的分区，通过数据驱动的混合策略来解决多目标燃烧优化问题；提出了基于模糊关联规则的运行优化目标值确定方法，将数据挖掘技术引入火电厂优化过程；基于深度Q网络(DQN)和长短时记忆(LSTM)模块的集成燃烧优化系统ThermalNet；用计算流体动力学(CFD)模拟生成一些数据作为人工神经网络建模的训练样本，加上运行历史运行数据，建立了用于预测锅炉运行和排放特性的人工神经网络(ANN)模型；采用动态数据挖掘模型用于调节门流量特性计算等一些方法用以提高机组运行的效率和智能化程度。由于机组运行始终处于动态变化过程中，认为同一时段内系统参数的均值必须在机组处于相对稳定的工况下才能真实反映机组运行状态，因此需要界定机组的稳定工况，计算并推选样本数据集。

在中国公开了(申请号为：201310574516.2，授权公告号:CN 103574581 B)一种火电机组NOx燃烧优化方法及系统，采集锅炉运行参数；计算所述锅炉运行参数与锅炉效率和NOx排放量之间的相关系数，选择符合相关系数条件的所述锅炉运行参数；以及根据所选择的符合所述相关系数条件的所述锅炉运行参数计算所述锅炉效率和所述NOx排放量，根据符合所述锅炉效率和所述NOx排放量条件的所述锅炉运行参数调整锅炉的运行参数。

针对上述专利一种火电机组NOx燃烧优化方法及系统，发现以下不足：由于机组运行参数变化越剧烈，以变化参数计算的机组经济指标偏差就越大。因此，在分析机组的运行经济性时，大部分文献均认为，只有稳定工况下的数据才是有效和具有参考价值的。机组历史数据库中的数据记录了机组的所有运行工况，因此需要对历史数据进行有效工况的检测与辨识。

发明内容

本发明针对上述的问题提供一种结合稳态分量和波动数据进行相关性分析的基于动态运行数据分析的燃烧优化控制方法。

本发明是这样实现的：一种基于动态运行数据分析的燃烧优化控制方法，其特征在于：所述该方法包括以下步骤：

步骤1：获取机组过去一段时间内DCS存储的相关历史运行数据；

步骤2：设全部运行数据集为：A(m,k)，其中m为采样序号，m＝1,2,3......M，M为总数据采样数；k代表运行参数，k＝1,2,3......K；将不同负荷、气温区间的数据包归在不同的数据子集中，运行数据表示为：X(i,j,k,l)i代表划分的负荷区间数，i＝1,2,3......I个；j代表划分的气温区间数，j＝1,2,3......J个；k代表运行参数，k＝1,2,3......K；l代表该子集内数据包的个数，也即工况点个数，l＝1,2,3......L(i，j)个，各个子集内工况点个数是存在差异的，M＝∑L(i,j)；由式(1)计算不同负荷、不同气温区间的运行数据子集内的所有运行参数的平均值：

步骤3：对计算得到的

步骤4：获得运行数据的稳态分量后，从全部运行数据集A(m,k)中减去对应的稳态分量

步骤5：计算出各参数与总煤量和/或烟气NOx浓度之间的相关系数；

步骤6：根据得到的相关系数，分析正相关系数大于0.1和负相关系数小于-0.1的运行参数对燃烧经济性和污染物生成的影响，给出优化建议。

优选的，所述步骤1中，获取机组过去一段时间内DCS存储的相关历史运行数据，可以将机组整个负荷、气温区间各划分为20-30个区间。

优选的，所述步骤2中选取的运行参数涵盖锅炉运行主要监测和控制参数，包括实际负荷、主蒸汽参数中流量、压力、温度、给水流量、总燃料量、各层(角)燃烧器给煤量、总风量、送风量、引风量、全部二次风门开度、烟气含氧量、排烟温度、NOx浓度、送风机入口风温。

优选的，所述步骤4中，根据不同运行参数的数值相差较大，进一步计算相对波动分量ΔA(m,k)如下：

其中，(i,j)为m采样时刻机组所处的负荷和气温区间序号。

优选的，所述步骤5中，采用线性相关系数法计算出各参数与总煤量和/或烟气NOx浓度之间的相关系数：

其中：相关系数r(k,k')是针对全部k个采样参数，当r(k,k')＝r(k',k)且r(k,k)＝1；分析每一个采样参数与其他所有采样参数之间的相关系数的变化情况。

优选的，所述步骤6中，选择DCS燃料量的大小直接作为锅炉效率和机组经济性评价参考指标，加上NOx排放值，联合作为优化工况筛选依据；

其中L'(i',j')是L(i,j)的子集，k

根据获得的L'(i',j')后，再使用公式(1)计算各个区段的均值，作为优化控制参数的基础。

本发明的有益效果：

1、通过将一段时间内全部机组动态运行数据根据负荷、气温划分二维区间的基础上，计算不同负荷、气温区间的运行数据的平均值，并计算出所有运行参数动态分量之间的相关系数，得出了波动性数据，通过波动性数据揭示了运行参数之间相关变化规律的特性，为操作员控制锅炉运行数据提供参考，有助于提高锅炉+运行的稳定性，减少锅炉报警次数。

2、通过选用机组含动态工况的运行数据，对工况按负荷、气温划分得更加细密，使用优化规律的锅炉运行结果显示锅炉效率明显提高、NOx排放量明显降低，使得到的稳态分量和波动分量更加准确，计算出的相关系数更加精准，得到的优化方法更适用。

附图说明

图1为本发明的操作流程图。

图2为主蒸汽流量随时间变化的示意图。

图3为参数实际运行值、平均运行值以及相对脉动值中总煤量和总风量比较图。

图4为参数实际运行值、平均运行值以及相对脉动值中烟气NOx浓度和烟气氧量比较图。

图5为参数实际运行值、平均运行值以及相对脉动值中主蒸汽温度和B后墙燃尽风出口电动调节门1位置比较图。

具体实施方式

以下结合附图对本发明做进一步概况。

一种基于动态运行数据分析的燃烧优化控制方法，

步骤1：获取机组过去一年的DCS存储的相关历史运行数据；

选取机组过去一年时间段内DCS存储的历史数据，以1-3min为周期,将机组整个负荷、气温区间各划分为20-30个区间，对于600MW级机组，基本上10MW一个负荷段，一个负荷区段内的负荷变化不大，可以认为负荷基本不变；

步骤2：设全部运行数据集为：A(m,k)，其中m为采样序号，m＝1,2,3......M，M为总数据采样数；k代表运行参数，k＝1,2,3......K；将不同负荷、气温区间的数据包归在不同的数据子集中，运行数据表示为：X(i,j,k,l)i代表划分的负荷区间数，i＝1,2,3......I个；j代表划分的气温区间数，j＝1,2,3......J个；k代表运行参数，k＝1,2,3......K；l代表该子集内数据包的个数，也即工况点个数，l＝1,2,3......L(i，j)个，各个子集内工况点个数是存在差异的，M＝∑L(i,j)；选取的运行参数涵盖锅炉运行主要监测和控制参数，包括实际负荷、主蒸汽参数中流量、压力、温度、给水流量、总燃料量、各层(角)燃烧器给煤量、总风量、送风量、引风量、全部二次风门开度、烟气含氧量、排烟温度、NOx浓度、送风机入口风温；

由式(1)计算不同负荷、不同气温区间的运行数据子集内的所有运行参数的平均值：

步骤3：对计算得到的

步骤4：获得运行数据的稳态分量后，从全部运行数据集A(m,k)中减去对应的稳态分量

根据不同运行参数的数值相差较大，进一步计算相对波动分量ΔA(m,k)如下：

其中，(i,j)为m采样时刻机组所处的负荷和气温区间序号；

步骤5：计算出各参数与总煤量和/或烟气NOx浓度之间的相关系数；

采用线性相关系数法计算出各参数与总煤量和/或烟气NOx浓度之间的相关系数：

其中：相关系数r(k,k')是针对全部k个采样参数，当r(k,k')＝r(k',k)且r(k,k)＝1；分析每一个采样参数与其他所有采样参数之间的相关系数的变化情况；

步骤6：根据得到的相关系数，分析正相关系数大于0.1和负相关系数小于-0.1的运行参数对燃烧经济性和污染物生成的影响；给出优化建议；DCS燃料量的大小直接作为锅炉效率和机组经济性评价参考指标；加上NOx排放值，联合作为优化工况筛选依据；

其中L'(i',j')是L(i,j)的子集，k

根据获得的L'(i',j')后，在使用公式(1)计算各个区段的均值，作为优化控制参数的基础。

采用某电厂350MW燃煤汽轮发电机组上进行了数据运行分析和优化运行试验，选用机组锅炉为超临界参数变压运行直流炉，东方锅炉(集团)股份有限公司设计制造，锅炉型号:DG1150/25.4-n2型。型式为单炉膛、前后墙对冲燃烧方式，尾部双烟道结构，采用挡板调节再热器温、一次再热、平衡通风、固态排渣、全钢构架、全悬吊结构、门型露天布置。锅炉炉膛布置膜式水冷壁，炉膛出口处布置屏式过热器，水平烟道内装设高、低温两级过热器，尾部竖井交错布置两级省煤器和两级空气预热器，水平烟道和转向室均为膜式壁顶棚管包墙管。锅炉宽度15101.2mm，锅炉深度13678.8mm。采用中速磨煤机、冷一次风机正压直吹、负压炉膛、平衡通风制粉燃烧系统。配5台中速磨煤机，其中4台运行，1台备用。煤粉细度Rg＝21％。以1分钟为采样周期，运行数据采集数量为24000个，大约为16天的运行数据。

选取主蒸汽流量作为锅炉和机组出力的标记参数，图2是主蒸汽流量随时间变化的数据图。图3、图4和图5中给出了典型参数实际运行值、平均运行值以及相对脉动值比较，包括总煤量(a)、总风量(b)、烟气NOx浓度(c)、烟气氧量(d)、主蒸汽温度(e)、B后墙燃尽风出口电动调节门1位置(f)。脉动值均在各图的下方显示。总煤量、总风量、B后墙燃尽风出口电动调节门1位置随负荷变化很明显；烟气氧量随负荷有一定变化；烟气NOx浓度、主蒸汽温度随负荷变化不明显。它们的脉动值均接近0均值的随机数。

表1为左侧全部分析参数与总煤量；中侧为#2炉A侧反应器入口烟气NOx浓度；右侧为总煤量+#2炉A侧反应器入口烟气NOx浓度的相关系数。

对表1的数据进行分析，发现运行中的总风量、烟气氧量、送风机A、B动叶调节门控制指令、引风量、#2炉B侧反应器入口烟气NOx浓度、#2A、#2B侧氨气混合器前氨气流量、烟气氧量、D组燃烧器A侧、B侧、C组燃烧器B侧、A侧二次风母管调节门位置等参量与煤量或烟气NOx浓度的正相关系数较大，对这一类参数进行适当的降低将有助于降低总煤耗量，降低污染物排放，以及提高机组运行的经济性等。

在表1中正相关系数大于0.1的阈值，代表两个参数之间呈现明显的同向变化；给煤机A、B、E给煤率、#1、#2SCR反应器加氨稀释空气流量，意味着燃料消耗量将会响应的减少；如表1中，#2炉A侧、#2炉B侧反应器入口烟气NOx浓度、A侧、B侧高温过热器集箱蒸汽压力，则表示参数值与总煤量呈负相关，为了降低总燃料量消耗，运行中这些参数应适当增加。

实际运行中，不会单纯追求运行经济性，或者氮氧化物排放，通常是综合优化经济性和污染物生成。根据动态数据分析得到的运行参数优化控制规律，减去各个参数的运行平均值得到的结果的正负，即优化调整的方向；与表1中右侧综合考虑煤量和NOx综合优化得到的正相关系数大于0.1的主要控制参数的调整方向一致，证明动态运行数据的动态分量的相关性分析得到的结果与燃烧优化有内在的联系。

工作原理：将一段时间内全部机组动态运行数据根据负荷、气温划分二维区间的基础上，计算不同负荷、气温区间的运行数据的平均值，然后从实际运行数据中减去其平均值，获得运行数据的动态分量，并计算出所有运行参数动态分量之间的相关系数。给出了总煤量、NOx生成量、以及总煤量和NOx生成量综合参数与所有运行参数动态分量之间的相关系数计算值，分析了正相关系数大于0.1和负相关系数小于-0.1的运行参数对燃烧经济性和污染物生成的影响，以及运行优化调整的建议。

以上所述仅为本发明的实施方式而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理内所作的任何修改、等同替换、改进等，均应包括在本发明的权利要求范围之内。