本发明属于火电脱销技术领域，涉及一种scr(选择性催化还原)脱硝智能喷氨优化方法及系统，尤其涉及一种基于矩阵式测量仪表和模型自适应pid的scr脱硝智能喷氨优化方法及系统。
背景技术：
目前，为了实现氮氧化物超低排放，大部分燃煤火电厂都安装了scr脱硝装置，采用cems系统(烟气在线监测系统)采集入口、出口氮氧化物浓度，然后再进行pid(比例-积分-微分)反馈控制。
pid反馈控制如图1所示，氨流量计算值为(入口氮氧化物浓度测量值-出口氮氧化物浓度测量值)×烟气量×氨氮摩尔比，其中氨氮摩尔比基本为固定值，入口氮氧化物浓度，出口氮氧化物浓度，烟气量通过仪表测量得到。
上述方案存在下述问题：
1.入口、出口氮氧化物浓度采用cems系统测量，cems系统采样管线比较长，造成测量纯迟延较大，测量值有2-3分钟延迟；
2.现有cems系统采用单点取样测量，因此导致测量数据不能代表整个截面平均浓度；
3.采用pid控制，pid参数在初始设定好后就不再改变，所以在机组负荷运行工况变化时，脱硝系统调整不及时，容易出现超标排放；
4.为了保证达标排放，电厂运行人员将pid参数的设定值设定得非常低，导致喷氨量过大，控制系统超调量较大，系统响应速率较低。不仅浪费还原剂，还增加了后续设备堵塞的风险。
因此，开发一种scr智能优化喷氨系统，实现精准喷氨对于scr脱硝装置安全经济运行具有重要意义。
技术实现要素：
针对现有技术中存在的技术问题，本发明的目的在于提供一种基于模型自适应pid的scr脱硝智能喷氨优化方法及系统。
一种基于模型自适应pid的scr脱硝智能喷氨优化方法，其步骤包括：
1)基于燃煤火电机组燃烧系统的运行数据，建立入口氮氧化物浓度模型；然后利用该入口氮氧化物浓度模型和当前采集的运行数据，得到入口氮氧化物浓度的预测值；
2)根据入口氮氧化物浓度预测值计算当前运行工况下喷氨量作为前馈信号，将该喷氨量前馈信号发送给scr脱硝系统；
3)测量燃煤火电机组燃烧系统的出口氮氧化物浓度的实测值，将该实测值与出口氮氧化物浓度设定值做偏差后输入自适应pid控制器，自适应pid控制器基于该偏差得到喷氨量反馈控制信号并发送给scr脱硝系统；
4)scr脱硝系统根据该喷氨量前馈信号和喷氨量反馈控制信号控制scr脱硝反应器的喷氨量。
进一步的，通过主元分析算法对所述运行数据进行划分，确定出与入口氮氧化物关联紧密的参数；然后将确定出的与入口氮氧化物关联紧密的参数作为训练数据，采用深度学习算法或者最小二乘支持向量机算法，建立入口氮氧化物浓度模型。
进一步的，所述运行数据包括机组负荷、一次风量、二次风量、总煤量、磨煤机运行台数。
进一步的，得到所述喷氨量前馈信号的方法为：将所述预测值乘以烟气量，然后再乘以氨氮摩尔比计算出氨气量需求量，作为scr脱硝系统的所述喷氨量前馈信号。
进一步的，所述自适应pid控制器动态辨识scr喷氨系统模型，根据模型变化自适应调整pid参数；自适应pid控制器参数可以采用遗传算法或者pso算法优化得到。进一步的，采用多点矩阵式测量仪表测量所述出口氮氧化物浓度的实测值。
一种基于模型自适应pid的scr脱硝智能喷氨优化系统，其特征在于，包括喷氨量前馈信号产生单元、自适应pid控制器和scr脱硝系统；其中，
所述喷氨量前馈信号产生单元，连接煤火电机组燃烧系统运行数据传感器，基于燃煤火电机组燃烧系统运行数据传感器采集的运行数据，建立入口氮氧化物浓度模型；并利用该入口氮氧化物浓度模型和当前采集的运行数据，得到入口氮氧化物浓度的预测值，然后根据入口氮氧化物浓度预测值计算当前运行工况下的喷氨量前馈信号；
所述自适应pid控制器，连接出口氮氧化物浓度测量仪表，用于根据燃煤火电机组燃烧系统的出口氮氧化物浓度的实测值与出口氮氧化物浓度设定值的偏差得到喷氨量反馈控制信号；
所述scr脱硝系统，分别连接所述喷氨量前馈信号产生单元和所述自适应pid控制器以接收所述喷氨量前馈信号和喷氨量反馈控制信号，用于根据该喷氨量前馈信号和喷氨量反馈控制信号控制scr脱硝反应器的喷氨量。
本发明首先采集燃煤火电机组燃烧系统一次风、二次风、负荷等数据，然后基于采集数据建立入口氮氧化物浓度模型，然后根据当前采集数据和入口氮氧化物浓度模型得到入口氮氧化物浓度预测值，然后根据入口氮氧化物浓度预测值计算当前运行工况下喷氨量前馈，将前馈加入到氮氧化物浓度反馈控制回路，实现精准喷氨。
入口氮氧化物浓度预测方法：
根据锅炉运行数据，初步选择机组负荷、一次风量、二次风量、总煤量、磨煤机运行台数参数及对应的入口氮氧化物浓度，首先通过主元分析算法对所选参数进行分析，确定出与入口氮氧化物关联紧密的参数，即确定出主要参数和次要参数。然后，将确定出的与入口氮氧化物关联紧密的参数作为训练数据，采用深度学习算法或者最小二乘支持向量机算法，建立入口氮氧化物浓度模型，从而利用该入口氮氧化物浓度模型实现对入口氮氧化物浓度的预测。其中，与入口氮氧化物关联紧密的参数包括总风量、总煤量、风煤比、一次风与二次风比值。
然后，基于运行数据建立scr喷氨系统模型，检测到运行工况变化时，自适应调整pid参数，从而实现优化控制。
scr喷氨系统模型是根据入口氮氧化物浓度和出口氮氧化物浓度建立的脱硝模型，基于机理与数据相结合建立，数据为入口氮氧化物浓度和出口氮氧化物浓度参数。基于当前喷氨系统模型参数，采用遗传算法或者pso优化算法，可得到当前工况下的最佳pid控制参数。在实际运行中，系统不断在线辨识scr喷氨模型，当监测到模型变化比较大的时候，重新计算pid参数，从而实现pid参数自适应调整，始终保持较好的控制参数。
本发明也可通过现场测试获得典型工况的系统特性参数。然后，根据不同典型工况，确定几组pid控制器参数。在实际运行中，自动检测运行工况，然后确定采用当前工况对应的较优pid参数，即根据运行工况不同进行切换，实现分段pid。
本发明入口、出口氮氧化物浓度采用矩阵式cems测量仪表进行测量，可以测得浓度的平均值。
综上，本发明的创新之处在于：
(1)针对氮氧化物测量延迟大的问题，本发明提出了采用软测量的方案。软测量是通过建立锅炉运行参数，比如一次风、二次风、负荷等数据与氮氧化物浓度的模型，通过计算的方式得到氮氧化物浓度的值，从而可以避免实际测量仪表的延迟，从而实现提前控制。
(2)针对现有仪表单点取样测量，测量结果不具有代表性的问题，提出了采用矩阵式测量仪表，矩阵式测量是通过多点取样然后再混合，测量结果能反应浓度平均值。
(3)针对控制算法调节特性差的问题，提出了采用智能前馈加反馈相结合的控制方案，智能前馈主要是通过预测的氮氧化物浓度计算得到。反馈方面主要是改进控制器。
通过本发明可实现：
1.脱硝出口氮氧化物浓度偏差可以控制在5～10mg/nm3之间，可以保证在变工况情况下稳定达标排放。
2.脱硝控制系统的系统响应速率高，投资较小，改造工期短。
附图说明
图1为现有系统控制流程图；
图2为本发明的控制流程图；
图3为本发明系统控制平台架构图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本方案主要是提出一种入口氮氧化物浓度软测量与工况自适应pid相结合的脱硝优化控制方法流程。
本发明首先采集燃煤火电机组燃烧系统一次风、二次风、负荷等数据及与这些参数对应的入口氮氧化物浓度，建立入口氮氧化物浓度模型。
如图2所示，基于锅炉一次风量、二次风量、给煤量、制粉系统运行台数等参数，根据上述入口氮氧化物浓度模型，得到入口氮氧化物浓度预测值，再乘以烟气量，然后再乘以氨氮摩尔比计算出氨气量需求量，作为scr脱硝系统前馈信号，前馈信号可以保证负荷变化时，可以快速调节喷氨量。
为了准确测量氮氧化物浓度，采用多点矩阵式测量仪表。多点矩阵式测量仪表是通过抽取烟道中3-5个截面的烟气，然后进行混合，混合后进入cems测量得到出口氮氧化物浓度实测值。这种方法测量的数据准确性较高。出口氮氧化物浓度设定与实测值做偏差后进入自适应pid控制器，作为喷氨量反馈控制信号。自适应pid控制器参数根据scr喷氨系统模型自适应调整。
其中，“烟气量”是指实际测量的锅炉烟气量，单位是nm3/h。“氨氮摩尔比”是指脱硝中氨气和氮氧化物浓度的比值，是脱硝中常用概念，一般为0.7-0.9。“氨气量需求量”的计算公式是：氨气需求量＝入口氮氧化物×烟气量×氨氮摩尔比。
为了保证该算法可以在现场应用，并对控制算法进行保密，本方案的实施基于脱硝优化控制平台实现。该平台的核心是一个高性能的优化控制器，通过数据采集卡件从dcs系统获得计算所需参数，经过计算后，再返回原dcs系统，实现闭环控制。控制器与原dcs系统通信采用modbus、rs485等通信方式，可以与国内主流dcs系统进行双向通信。脱硝优化控制器与dcs现场控制器可以实现无扰切换。控制平台如图3所示，优化控制器中主要包括系统通信模块和核心算法计算模块，系统通信模块主要负责实现与dcs进行数据输入，输出。核心算法模块主要是实现优化控制算法。系统运行时，首先通过系统通信模块采集运行数据，然后输入到核心算法模块，核心算法模块再输出到dcs系统，从而实现闭环控制。
表1为采用本发明方法和现有方法在实际运行中的对比，从这些对比中可以看出，采用本发明方法，可以实现scr脱硝装置的精准喷氨，提高系统运行效率，降低成本。
表1为实际运行效果对比
动态工况运行效果对比分析
1改造前运行效果
【1】机组升负荷工况
变工况下反应器出口nox含量的动态特性是考察脱硝系统自动控制效果最重要的部分之一。此部分分析升负荷过程中反应器出口nox含量的动态特性。
负荷由365mw附近升至395mw附近过程中，脱硝系统甲、乙侧反应器出口nox含量的统计分析如表2。
表2为365mw升至395mw负荷段脱硝性能数据
【2】机组降负荷工况
以下是降负荷过程中反应器出口nox含量的动态特性的分析。
负荷由500mw附近降至330mw附近过程中，脱硝系统甲、乙侧反应器出口nox含量的统计分析如表3。
表3为500mw降至330mw负荷段脱硝性能数据
2改造后运行效果
【1】机组升负荷工况
以下为机组负荷从400mw上升至500mw时的运行过程，一般在此升负荷区间，会启动一台磨煤机。磨的启停期间对反应器入口nox含量的短时影响很大，反应器出口nox含量超标也多半是受此影响。
负荷由400mw升至500mw过程中，脱硝系统甲、乙侧反应器出口nox含量的统计分析如表4。
表4为负荷由400mw升至500mw运行区间脱硝性能数据
负荷由550mw升至600mw运行区间，脱硝系统甲、乙侧反应器出口nox含量的统计分析如表5。
表5为负荷由550mw升至600mw运行区间脱硝性能数据
【2】机组降负荷工况
此部分将用同样的方式对降负荷过程中脱硝反应器的性能进行了分析。
负荷由600mw降至550mw运行区间，脱硝系统甲、乙侧反应器出口nox含量的统计分析如表6。
表6为负荷由600mw降至550mw运行区间脱硝性能数据
负荷由550mw降至500mw运行区间，脱硝系统甲、乙侧反应器出口nox含量统计分析如表7。
表7为负荷由550mw降至500mw运行区间脱硝性能数据
负荷由500mw降至400mw运行区间，脱硝系统甲、乙侧反应器出口nox含量统计分析如表8。
表8为负荷由500mw降至400mw运行区间脱硝性能数据
以上所述仅为本发明的较佳实施例，并不用以限制本发明，本领域的技术人员应可理解，凡在本发明的精神和原则之内所作的任何修改，等同替换或改进等，均应包含在本发明的保护范围之内，保护范围以权利要求书所界定者为准。