CN112418284B - A control method and system for a SCR denitration system in a full-operating power station - Google Patents
A control method and system for a SCR denitration system in a full-operating power station Download PDFInfo
- Publication number
- CN112418284B CN112418284B CN202011276325.4A CN202011276325A CN112418284B CN 112418284 B CN112418284 B CN 112418284B CN 202011276325 A CN202011276325 A CN 202011276325A CN 112418284 B CN112418284 B CN 112418284B
- Authority
- CN
- China
- Prior art keywords
- individual
- ammonia injection
- scr denitration
- injection amount
- denitration system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/34—Chemical or biological purification of waste gases
- B01D53/74—General processes for purification of waste gases; Apparatus or devices specially adapted therefor
- B01D53/86—Catalytic processes
- B01D53/8621—Removing nitrogen compounds
- B01D53/8625—Nitrogen oxides
- B01D53/8631—Processes characterised by a specific device
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biomedical Technology (AREA)
- Environmental & Geological Engineering (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Genetics & Genomics (AREA)
- Physiology (AREA)
- Analytical Chemistry (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Treating Waste Gases (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及燃煤电站优化控制技术领域,特别涉及一种全工况电站SCR脱硝系统的控制方法及系统。The present invention relates to the technical field of coal-fired power plant optimization control, and in particular to a control method and system for a SCR denitration system of a full-operation power plant.
背景技术Background Art
在能源结构转型的过程中,新能源电力规模化接入电网对燃煤机组提出了运行灵活性要求。快速深度变负荷意味着机组运行工况大范围快速变化,锅炉工况的变化会使得燃烧产生的NOx波动加剧,这无疑加大了机组实现NOx超低排放的难度。SCR(SelectiveCatalytic Reduction)——选择性催化还原法是目前国际上技术最成熟、应用最广泛的烟气脱硝技术。SCR是在催化剂的作用下,利用还原剂NH3等来有选择性地与烟气中的NOx反应并生成无毒无污染的N2和H2O。SCR脱硝是目前主流的烟气脱硝技术,其反应是一个复杂的物理化学过程,喷氨量较多可以降低NOx排放,但会增加经济成本,并导致氨逃逸增大,影响机组安全运行。In the process of energy structure transformation, the large-scale access of new energy power to the grid has put forward operational flexibility requirements for coal-fired units. Rapid deep load change means that the operating conditions of the unit change rapidly over a large range. The change of boiler operating conditions will increase the fluctuation of NO x produced by combustion, which undoubtedly increases the difficulty of the unit to achieve ultra-low NO x emissions. SCR (Selective Catalytic Reduction) - Selective catalytic reduction is currently the most mature and widely used flue gas denitrification technology in the world. SCR uses reducing agents such as NH 3 to selectively react with NO x in flue gas under the action of catalysts to generate non-toxic and non-polluting N 2 and H 2 O. SCR denitrification is the current mainstream flue gas denitrification technology. Its reaction is a complex physical and chemical process. A large amount of ammonia injection can reduce NO x emissions, but it will increase economic costs and lead to increased ammonia slip, affecting the safe operation of the unit.
目前现场应用的SCR控制系统主要分为两种,一种是固定摩尔比控制方式,这种方式属开环控制,无法满足超低排放要求;另一种是固定出口NOx浓度控制方式,目前现场多采用串级PID控制系统,其参数是针对设计(额定)工况下的反应器特性整定得到的,当机组工况发生较大变化时,若不能对SCR系统进行有效控制,就会造成脱硝效率低或氨逃逸现象严重,难以实现最优控制,对火电厂经济环保运行极为不利。因此,采用开环控制以及串级PID的控制方式往往难以取得较好的控制效果。At present, there are two main types of SCR control systems used on site. One is the fixed molar ratio control method, which is an open-loop control method and cannot meet the ultra-low emission requirements; the other is the fixed outlet NOx concentration control method. At present, the cascade PID control system is mostly used on site, and its parameters are adjusted according to the reactor characteristics under the design (rated) working conditions. When the unit working conditions change significantly, if the SCR system cannot be effectively controlled, it will cause low denitrification efficiency or serious ammonia escape, making it difficult to achieve optimal control, which is extremely unfavorable to the economic and environmental protection operation of thermal power plants. Therefore, it is often difficult to achieve good control effects using open-loop control and cascade PID control methods.
因此,如何对脱硝系统进行优化控制,在保证达标排放的同时实现机组经济运行是燃煤电站亟待解决的问题。Therefore, how to optimize the control of the denitrification system and achieve economical operation of the unit while ensuring compliance with emission standards is an urgent problem to be solved in coal-fired power plants.
发明内容Summary of the invention
本发明的目的是提供一种全工况电站SCR脱硝系统的控制方法及系统,以实现对脱硝系统进行优化控制,在保证达标排放的同时实现机组经济运行。The purpose of the present invention is to provide a control method and system for a SCR denitration system in a power station under full-operation conditions, so as to optimize the control of the denitration system and realize economic operation of the unit while ensuring compliance with emission standards.
为实现上述目的,本发明提供了如下方案:To achieve the above object, the present invention provides the following solutions:
一种全工况电站SCR脱硝系统的控制方法,所述控制方法包括如下步骤:A control method for a SCR denitration system of a power station under full working conditions, the control method comprising the following steps:
获取SCR脱硝系统的历史运行数据;Obtain historical operating data of the SCR denitrification system;
采用聚类算法将历史运行数据划分至多个负荷区间,获得多个负荷区间的历史运行数据;A clustering algorithm is used to divide the historical operation data into multiple load intervals to obtain the historical operation data of multiple load intervals;
分别利用每个负荷区间的历史运行数据训练神经网络模型,获得每个负荷区间的SCR脱硝系统输出预测模型;SCR脱硝系统输出预测模型用于根据SCR脱硝系统当前的工作状态数据和预测时间段内的每个控制时间点的喷氨量,预测SCR脱硝系统在预测时间段内的每个预测时间点的出口NOx浓度和出口氨逃逸;所述工作状态数据包括机组负荷、SCR脱硝系统入口的NOx浓度和烟气流量;The neural network model is trained by using the historical operation data of each load interval to obtain the SCR denitration system output prediction model of each load interval; the SCR denitration system output prediction model is used to predict the outlet NOx concentration and outlet ammonia slip of the SCR denitration system at each prediction time point in the prediction time period according to the current working state data of the SCR denitration system and the ammonia injection amount at each control time point in the prediction time period; the working state data includes the unit load, the NOx concentration and the flue gas flow rate at the inlet of the SCR denitration system;
根据SCR脱硝系统的负荷指令,获取负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型;According to the load instruction of the SCR denitration system, an output prediction model of the SCR denitration system in the load interval where the unit load corresponding to the load instruction is located is obtained;
根据负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型,采用遗传算法确定使多目标优化函数最优的喷氨量,作为最优喷氨量;According to the SCR denitration system output prediction model of the load interval where the unit load corresponding to the load instruction is located, a genetic algorithm is used to determine the ammonia injection amount that optimizes the multi-objective optimization function as the optimal ammonia injection amount;
根据所述最优喷氨量对SCR脱销系统进行控制。The SCR destocking system is controlled according to the optimal ammonia injection amount.
可选的,所述神经网络模型包括输入层、输出层和隐含层,所述输入层包括4个神经元,所述输出层包括2个神经元,所述输出层包括5个神经元。Optionally, the neural network model includes an input layer, an output layer and a hidden layer, the input layer includes 4 neurons, the output layer includes 2 neurons, and the output layer includes 5 neurons.
可选的,所述根据负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型,采用遗传算法确定使多目标优化函数最优的喷氨量,作为最优喷氨量,具体包括:Optionally, the SCR denitration system output prediction model in the load interval where the unit load corresponding to the load instruction is located uses a genetic algorithm to determine the ammonia injection amount that optimizes the multi-objective optimization function as the optimal ammonia injection amount, specifically including:
以预测时间段内每个控制时间点的喷氨量为个体,初始化遗传算法的种群;The ammonia injection amount at each control time point in the prediction period is taken as an individual to initialize the population of the genetic algorithm;
判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果;Determine whether the valve opening of the SCR denitration system corresponding to the ammonia injection amount at each control time point in each individual in the population is between the lower limit value and the upper limit value of the valve opening, and obtain a first determination result;
若所述第一判断结果表示否,则将SCR脱硝系统的阀门开度小于阀门开度的下限值的个体中的喷氨量用阀门开度的下限值对应的喷氨量代替,将SCR脱硝系统的阀门开度大于阀门开度的上限值的个体中的喷氨量用阀门开度的上限值对应的喷氨量代替;If the first judgment result indicates no, the ammonia injection amount in the individual of the SCR denitration system whose valve opening is less than the lower limit of the valve opening is replaced by the ammonia injection amount corresponding to the lower limit of the valve opening, and the ammonia injection amount in the individual of the SCR denitration system whose valve opening is greater than the upper limit of the valve opening is replaced by the ammonia injection amount corresponding to the upper limit of the valve opening;
将种群中每个个体输入SCR脱硝系统输出预测模型,获得每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度和预测出口氨逃逸量;Input each individual in the population into the SCR denitrification system output prediction model to obtain the predicted outlet NOx concentration and predicted outlet ammonia slip at each prediction time point in the prediction time period corresponding to each individual;
利用公式对每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度进行校正,获得每个个体对应的预测时间段内每个预测时间点的校正后的预测出口NOx浓度;其中,为k+i·s1时刻校正后的预测出口NOx浓度,为k+i·s1时刻的预测出口NOx浓度,为k时刻的预测出口NOx浓度,为k时刻的SCR脱硝系统的实际出口NOx浓度,i表示第i个预测时间点,s1为预测步长,r为校正系数;Using the formula Correct the predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual to obtain the corrected predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual; wherein, is the predicted outlet NOx concentration after correction at time k+i·s1, is the predicted outlet NOx concentration at time k+i·s 1 , is the predicted outlet NOx concentration at time k, is the actual outlet NOx concentration of the SCR denitrification system at time k, i represents the i-th prediction time point, s 1 is the prediction step length, and r is the correction coefficient;
根据每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度、预测出口氨逃逸量和校正后的预测出口NOx浓度,分别利用第一目标函数和第二目标函数,计算每个个体的第一目标函数值和第二目标函数值;Calculate the first objective function value and the second objective function value of each individual according to the predicted outlet NOx concentration, the predicted outlet ammonia slip amount and the corrected predicted outlet NOx concentration at each predicted time point within the prediction time period corresponding to each individual using the first objective function and the second objective function respectively;
确定种群中满足第二目标函数的个体中第一目标函数最小的个体作为第L次迭代的最优个体,将第L次迭代的最优个体与第L-1次迭代的全局最优个体中第一目标函数值较大的个体设置为第L次迭代的全局最优个体;Determine the individual with the smallest first objective function among the individuals in the population that satisfy the second objective function as the optimal individual of the Lth iteration, and set the individual with the larger first objective function value between the optimal individual of the Lth iteration and the global optimal individual of the L-1th iteration as the global optimal individual of the Lth iteration;
判断迭代次数是否大于迭代次数阈值,获得第二判断结果;Determine whether the number of iterations is greater than an iteration number threshold, and obtain a second determination result;
若所述第二判断结果表示否,则令迭代次数的数值增加1,采用遗传算法中的遗传、变异和重组的方式更新种群,返回步骤“判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果”;If the second judgment result indicates no, then the value of the number of iterations is increased by 1, and the population is updated by the inheritance, mutation and recombination methods in the genetic algorithm, and the step of "determining whether the valve opening of the SCR denitrification system corresponding to the ammonia injection amount at each control time point in each individual in the population is between the lower limit and the upper limit of the valve opening to obtain the first judgment result" is returned;
若所述第二判断结果表示是,输出第L次迭代的全局最优个体作为最优喷氨量。If the second judgment result indicates yes, the global optimal individual of the Lth iteration is output as the optimal ammonia injection amount.
可选的,所述第一目标函数为:Optionally, the first objective function is:
其中,J1表示第一目标函数值,为k+i·s1时刻的预测出口氨逃逸量,M2为液氨价格,P为预测时间段的预测时间点的数量,Qgas为烟气流量,为烟气含氧量,M1为排污费价格,为k+j·s2时刻的喷氨量,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,N为机组发电量,M3为电价补贴价格,ω1为第一权重系数,ω2为第二权重系数;Where J 1 represents the first objective function value, is the predicted outlet ammonia slip at time k+i·s1, M2 is the liquid ammonia price, P is the number of predicted time points in the prediction time period, Q gas is the flue gas flow rate, is the oxygen content of flue gas, M1 is the sewage fee price, is the ammonia injection amount at time k+j·s2, j is the jth control time point, s2 is the control step size, M is the number of control time points in the prediction time period, N is the power generation of the unit, M3 is the electricity price subsidy price, ω1 is the first weight coefficient, and ω2 is the second weight coefficient;
可选的,所述第二目标函数为:Optionally, the second objective function is:
其中,J2为第二目标函数值,P为预测时间段的预测时间点的数量,r(k+i·s1)为k+i·s1时刻出口NOx浓度的期望值,||Δu(k+j·s2)||为k+j·s2时刻的喷氨量与期望喷氨量的差值,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,ω3为第三权重系数。Among them, J 2 is the second objective function value, P is the number of prediction time points in the prediction time period, r(k+i·s 1 ) is the expected value of the outlet NO x concentration at moment k+i·s 1 , ||Δu(k+j·s 2 )|| is the difference between the ammonia injection amount at moment k+j·s 2 and the expected ammonia injection amount, j is the jth control time point, s 2 is the control step size, M is the number of control time points in the prediction time period, and ω 3 is the third weight coefficient.
一种全工况电站SCR脱硝系统的控制系统,所述控制系统包括:A control system for a full-operation power station SCR denitration system, the control system comprising:
历史运行数据获取模块,用于获取SCR脱硝系统的历史运行数据;A historical operation data acquisition module is used to obtain historical operation data of the SCR denitration system;
聚类模块,用于采用聚类算法将历史运行数据划分至多个负荷区间,获得多个负荷区间的历史运行数据;A clustering module is used to divide the historical operation data into multiple load intervals by using a clustering algorithm to obtain the historical operation data of the multiple load intervals;
训练模块,用于分别利用每个负荷区间的历史运行数据训练神经网络模型,获得每个负荷区间的SCR脱硝系统输出预测模型;SCR脱硝系统输出预测模型用于根据SCR脱硝系统当前的工作状态数据和预测时间段内的每个控制时间点的喷氨量,预测SCR脱硝系统在预测时间段内的每个预测时间点的出口NOx浓度和出口氨逃逸;所述工作状态数据包括机组负荷、SCR脱硝系统入口的NOx浓度和烟气流量;A training module is used to train a neural network model using historical operation data of each load interval to obtain an output prediction model of the SCR denitration system in each load interval; the output prediction model of the SCR denitration system is used to predict the outlet NOx concentration and outlet ammonia slip of the SCR denitration system at each predicted time point in the predicted time period according to the current working state data of the SCR denitration system and the injection amount of ammonia at each control time point in the predicted time period; the working state data includes unit load, NOx concentration at the inlet of the SCR denitration system and flue gas flow;
SCR脱硝系统输出预测模型选取模块,用于根据SCR脱硝系统的负荷指令,获取负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型;An SCR denitration system output prediction model selection module is used to obtain an SCR denitration system output prediction model for a load interval where the unit load corresponding to the load instruction is located according to the load instruction of the SCR denitration system;
最优喷氨量确定模块,用于根据负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型,采用遗传算法确定使多目标优化函数最优的喷氨量,作为最优喷氨量;The optimal ammonia injection amount determination module is used to determine the ammonia injection amount that optimizes the multi-objective optimization function according to the SCR denitration system output prediction model of the load interval where the unit load corresponding to the load instruction is located, and use a genetic algorithm to determine the optimal ammonia injection amount as the optimal ammonia injection amount;
控制模块,用于根据所述最优喷氨量对SCR脱销系统进行控制。A control module is used to control the SCR destocking system according to the optimal ammonia injection amount.
可选的,所述神经网络模型包括输入层、输出层和隐含层,所述输入层包括4个神经元,所述输出层包括2个神经元,所述输出层包括5个神经元。Optionally, the neural network model includes an input layer, an output layer and a hidden layer, the input layer includes 4 neurons, the output layer includes 2 neurons, and the output layer includes 5 neurons.
可选的,所述最优喷氨量确定模块,具体包括:Optionally, the optimal ammonia injection amount determination module specifically includes:
初始化子模块,用于以预测时间段内每个控制时间点的喷氨量为个体,初始化遗传算法的种群;The initialization submodule is used to initialize the population of the genetic algorithm by taking the ammonia injection amount at each control time point in the prediction time period as an individual;
第一判断子模块,用于判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果;The first judgment submodule is used to judge whether the valve opening of the SCR denitration system corresponding to the ammonia injection amount at each control time point of each individual in the population is between the lower limit value and the upper limit value of the valve opening, and obtain a first judgment result;
个体更新子模块,用于若所述第一判断结果表示否,则将SCR脱硝系统的阀门开度小于阀门开度的下限值的个体中的喷氨量用阀门开度的下限值对应的喷氨量代替,将SCR脱硝系统的阀门开度大于阀门开度的上限值的个体中的喷氨量用阀门开度的上限值对应的喷氨量代替;an individual updating submodule, for replacing the amount of ammonia injection in the individual whose valve opening of the SCR denitration system is less than the lower limit of the valve opening with the amount of ammonia injection corresponding to the lower limit of the valve opening, and replacing the amount of ammonia injection in the individual whose valve opening of the SCR denitration system is greater than the upper limit of the valve opening with the amount of ammonia injection corresponding to the upper limit of the valve opening if the first judgment result indicates no;
预测子模块,用于将种群中每个个体输入SCR脱硝系统输出预测模型,获得每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度和预测出口氨逃逸量;The prediction submodule is used to input each individual in the population into the SCR denitration system output prediction model to obtain the predicted outlet NOx concentration and predicted outlet ammonia slip at each prediction time point in the prediction time period corresponding to each individual;
校正子模块,用于利用公式对每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度进行校正,获得每个个体对应的预测时间段内每个预测时间点的校正后的预测出口NOx浓度;其中,为k+i·s1时刻校正后的预测出口NOx浓度,为k+i·s1时刻的预测出口NOx浓度,为k时刻的预测出口NOx浓度,为k时刻的SCR脱硝系统的实际出口NOx浓度,i表示第i个预测时间点,s1为预测步长,r为校正系数;Correction submodule, used to use the formula Correct the predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual to obtain the corrected predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual; wherein, is the predicted outlet NOx concentration after correction at time k+i·s 1 , is the predicted outlet NOx concentration at time k+i·s 1 , is the predicted outlet NOx concentration at time k, is the actual outlet NOx concentration of the SCR denitrification system at time k, i represents the i-th prediction time point, s 1 is the prediction step length, and r is the correction coefficient;
目标函数值计算子模块,用于根据每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度、预测出口氨逃逸量和校正后的预测出口NOx浓度,分别利用第一目标函数和第二目标函数,计算每个个体的第一目标函数值和第二目标函数值;An objective function value calculation submodule, for calculating the first objective function value and the second objective function value of each individual according to the predicted outlet NOx concentration, the predicted outlet ammonia slip amount and the corrected predicted outlet NOx concentration at each prediction time point within the prediction time period corresponding to each individual, using the first objective function and the second objective function respectively;
最优个体确定子模块,用于确定种群中满足第二目标函数的个体中第一目标函数最小的个体作为第L次迭代的最优个体,将第L次迭代的最优个体与第L-1次迭代的全局最优个体中第一目标函数值较大的个体设置为第L次迭代的全局最优个体;The optimal individual determination submodule is used to determine the individual with the smallest first objective function among the individuals in the population that satisfy the second objective function as the optimal individual of the Lth iteration, and set the individual with the larger first objective function value between the optimal individual of the Lth iteration and the global optimal individual of the L-1th iteration as the global optimal individual of the Lth iteration;
第二判断子模块,用于判断迭代次数是否大于迭代次数阈值,获得第二判断结果;A second judgment submodule is used to judge whether the number of iterations is greater than an iteration number threshold, and obtain a second judgment result;
种群更新子模块,用于若所述第二判断结果表示否,则令迭代次数的数值增加1,采用遗传算法中的遗传、变异和重组的方式更新种群,返回步骤“判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果”;A population update submodule, for, if the second judgment result indicates no, increasing the value of the number of iterations by 1, updating the population by adopting the inheritance, mutation and recombination methods in the genetic algorithm, and returning to the step of "judging whether the valve opening of the SCR denitration system corresponding to the ammonia injection amount at each control time point of each individual in the population is between the lower limit and the upper limit of the valve opening, and obtaining the first judgment result";
最优喷氨量输出子模块,用于若所述第二判断结果表示是,输出第L次迭代的全局最优个体作为最优喷氨量。The optimal ammonia injection amount output submodule is used to output the global optimal individual of the Lth iteration as the optimal ammonia injection amount if the second judgment result indicates yes.
可选的,所述第一目标函数为:Optionally, the first objective function is:
其中,J1表示第一目标函数值,为k+i·s1时刻的预测出口氨逃逸量,M2为液氨价格,P为预测时间段的预测时间点的数量,Qgas为烟气流量,为烟气含氧量,M1为排污费价格,为k+j·s2时刻的喷氨量,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,N为机组发电量,M3为电价补贴价格,ω1为第一权重系数,ω2为第二权重系数;Where J 1 represents the first objective function value, is the predicted outlet ammonia slip at time k+i·s1, M2 is the liquid ammonia price, P is the number of predicted time points in the prediction time period, Q gas is the flue gas flow rate, is the oxygen content of flue gas, M1 is the sewage fee price, is the ammonia injection amount at time k+j·s2, j is the jth control time point, s2 is the control step size, M is the number of control time points in the prediction time period, N is the power generation of the unit, M3 is the electricity price subsidy price, ω1 is the first weight coefficient, and ω2 is the second weight coefficient;
可选的,所述第二目标函数为:Optionally, the second objective function is:
其中,J2为第二目标函数值,P为预测时间段的预测时间点的数量,r(k+i·s1)为k+i·s1时刻出口NOx浓度的期望值,||Δu(k+j·s2)||为k+j·s2时刻的喷氨量与期望喷氨量的差值,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,ω3为第三权重系数。Among them, J 2 is the second objective function value, P is the number of prediction time points in the prediction time period, r(k+i·s 1 ) is the expected value of the outlet NOx concentration at moment k+i·s 1 , ||Δu(k+j·s 2 )|| is the difference between the ammonia injection amount at moment k+j·s 2 and the expected ammonia injection amount, j is the jth control time point, s 2 is the control step size, M is the number of control time points in the prediction time period, and ω 3 is the third weight coefficient.
根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:
本发明公开了一种全工况电站SCR脱硝系统的控制方法及系统,所述控制方法:将脱硝成本加入优化目标函数,采用预测控制结构,结合神经网络和遗传算法进行模型建立和控制量寻优,实现了喷氨量的优化控制。通过构建神经网络模型预测SCR反应器的出口NOx浓度未来的变化趋势,通过采用遗传算法寻优确定SCR脱硝系统的喷氨流量,并对喷氨阀门的阀门开度进行调整,该方法能够较好地克服SCR脱硝系统大惯性、大迟延的缺点,提高喷氨量控制对机组负荷变化的响应速度,改善SCR脱硝系统的动态调节品质。The present invention discloses a control method and system for a full-operation power station SCR denitration system. The control method includes: adding denitration cost to an optimization objective function, adopting a predictive control structure, combining a neural network and a genetic algorithm to establish a model and optimize the control quantity, thereby realizing the optimization control of the ammonia injection amount. By constructing a neural network model to predict the future trend of the outlet NO x concentration of the SCR reactor, by adopting a genetic algorithm to optimize and determine the ammonia injection flow rate of the SCR denitration system, and adjusting the valve opening of the ammonia injection valve, the method can better overcome the shortcomings of the large inertia and large delay of the SCR denitration system, improve the response speed of the ammonia injection amount control to the load change of the unit, and improve the dynamic regulation quality of the SCR denitration system.
采用多目标控制方式,预测控制中既考虑了阀门开度的上下限制、也考虑了出口NOx浓度的排放浓度、氨逃逸以及脱硝系统的经济性指标,避免因执行机构饱和从而影响系统性能,能够在保证出口NOx浓度达到目标值的基础上,尽量减少喷氨量,有效降低运行费用和氨逃逸率。A multi-objective control method is adopted. The predictive control takes into account the upper and lower limits of the valve opening, the emission concentration of the outlet NOx concentration, ammonia escape, and the economic indicators of the denitrification system, so as to avoid the influence of the system performance on the saturation of the actuator. On the basis of ensuring that the outlet NOx concentration reaches the target value, the amount of ammonia injection can be minimized, and the operating cost and ammonia escape rate can be effectively reduced.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.
图1为本发明提供的一种全工况电站SCR脱硝系统的控制方法的流程图;FIG1 is a flow chart of a control method for a full-operation power station SCR denitration system provided by the present invention;
图2为本发明提供的一种全工况电站SCR脱硝系统的控制方法的原理示意图;FIG2 is a schematic diagram showing the principle of a control method for a SCR denitration system of a power station under full operating conditions provided by the present invention;
图3为本发明提供的采用遗传算法确定使多目标优化函数最优的喷氨量的原理示意图FIG. 3 is a schematic diagram of the principle of using a genetic algorithm to determine the optimal ammonia injection amount for a multi-objective optimization function provided by the present invention.
图4为本发明提供的采用遗传算法确定使多目标优化函数最优的喷氨量的流程图。FIG4 is a flow chart of the present invention for using a genetic algorithm to determine the optimal ammonia injection amount for the multi-objective optimization function.
具体实施方式DETAILED DESCRIPTION
本发明的目的是提供一种全工况电站SCR脱硝系统的控制方法及系统,以实现对脱硝系统进行优化控制,在保证达标排放的同时实现机组经济运行。The purpose of the present invention is to provide a control method and system for a SCR denitration system in a power station under full-operation conditions, so as to optimize the control of the denitration system and realize economic operation of the unit while ensuring compliance with emission standards.
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对发明作进一步详细的说明。In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.
需要说明的是,本发明实施例及附图中的术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "including" and "having" and any variations thereof in the embodiments of the present invention and the accompanying drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device including a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products or devices.
如图1-2所示,本发明提供一种全工况电站SCR脱硝系统的控制方法,所述控制方法包括如下步骤:As shown in FIG. 1-2 , the present invention provides a control method for a SCR denitration system of a power station under full operating conditions, and the control method comprises the following steps:
步骤101,获取SCR脱硝系统的历史运行数据。Step 101, obtaining historical operation data of the SCR denitration system.
历史运行数据包括所述机组负荷、SCR脱硝系统入口NOx浓度、烟气流量、喷氨量、SCR脱硝系统出口NOx浓度和氨逃逸。The historical operation data include the load of the unit, NOx concentration at the inlet of the SCR denitrification system, flue gas flow, ammonia injection amount, NOx concentration at the outlet of the SCR denitrification system and ammonia slip.
本发明采集历史运行数据时的频率为每2分钟一次。The frequency of collecting historical operation data in the present invention is once every 2 minutes.
步骤101,采用聚类算法将历史运行数据划分至多个负荷区间,获得多个负荷区间的历史运行数据。Step 101 : using a clustering algorithm to divide historical operation data into multiple load intervals, and obtaining historical operation data of multiple load intervals.
根据负荷区间将历史数据通过聚类算法进行划分,分别记为为高负荷区间、中负荷区间、低负荷区间。The historical data is divided into high load interval, medium load interval and low load interval through clustering algorithm according to the load interval.
步骤102,分别利用每个负荷区间的历史运行数据训练神经网络模型,获得每个负荷区间的SCR脱硝系统输出预测模型;SCR脱硝系统输出预测模型用于根据SCR脱硝系统当前的工作状态数据和预测时间段内的每个控制时间点的喷氨量,预测SCR脱硝系统在预测时间段内的每个预测时间点的出口NOx浓度和出口氨逃逸;所述工作状态数据包括机组负荷、SCR脱硝系统入口的NOx浓度和烟气流量。Step 102, respectively train the neural network model using the historical operation data of each load interval to obtain the SCR denitration system output prediction model for each load interval; the SCR denitration system output prediction model is used to predict the outlet NOx concentration and outlet ammonia slip of the SCR denitration system at each prediction time point within the prediction time period according to the current working status data of the SCR denitration system and the injection amount of ammonia at each control time point within the prediction time period; the working status data includes the unit load, the NOx concentration at the inlet of the SCR denitration system, and the flue gas flow rate.
根据不同负荷区间的历史数据分别训练神经网络模型,在运行控制算法的过程中,根据SCR脱硝系统实时的输入数据选择不同的神经网络模型,能更好地适应机组工况大范围变化的背景,模型预测的精确度显著提升。The neural network models are trained separately according to the historical data of different load ranges. During the operation of the control algorithm, different neural network models are selected according to the real-time input data of the SCR denitrification system. This can better adapt to the background of large-scale changes in the unit operating conditions and significantly improve the accuracy of model prediction.
所述神经网络模型为LSTM神经网络,包括输入层、输出层和隐含层,所述输入层包括4个神经元,所述输出层包括2个神经元,所述输出层包括5个神经元。The neural network model is an LSTM neural network, including an input layer, an output layer and a hidden layer. The input layer includes 4 neurons, the output layer includes 2 neurons, and the output layer includes 5 neurons.
本发明在训练之前,还对历史运行数据进行归一化预处理。训练时,设定神经网络模型训练次数为1000次,学习速率为0.05,最小误差为10-3。The present invention also performs normalization preprocessing on the historical operation data before training. During training, the neural network model training times are set to 1000 times, the learning rate is 0.05, and the minimum error is 10 -3 .
步骤103,根据SCR脱硝系统的负荷指令,获取负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型。Step 103 , according to the load instruction of the SCR denitration system, an output prediction model of the SCR denitration system in the load interval where the unit load corresponding to the load instruction is located is obtained.
步骤104,根据负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型,采用遗传算法确定使多目标优化函数最优的喷氨量,作为最优喷氨量。Step 104, according to the SCR denitration system output prediction model in the load interval where the unit load corresponding to the load instruction is located, a genetic algorithm is used to determine the ammonia injection amount that optimizes the multi-objective optimization function as the optimal ammonia injection amount.
本发明的多目标优化函数使脱硝控制系统中SCR出口NOx浓度满足要求的前提下,兼顾了出口NOx、氨逃逸和相关经济成本,在达标排放的前提下实现机组安全和经济运行。并将阀门开度的上下限制加入约束目标函数。The multi-objective optimization function of the present invention takes into account the outlet NOx , ammonia slip and related economic costs under the premise that the SCR outlet NOx concentration in the denitration control system meets the requirements, and realizes the safe and economical operation of the unit under the premise of meeting the emission standards. The upper and lower limits of the valve opening are added to the constraint objective function.
本发明在预测时域内进行在线滚动优化,通过使用遗传算法在目标函数中寻优确定控制量序列,使神经网络模型的的预测输出能够最大程度的接近出口NOx浓度以及氨逃逸的期望值。The present invention performs online rolling optimization in the prediction time domain, and uses a genetic algorithm to optimize and determine the control quantity sequence in the objective function, so that the prediction output of the neural network model can be as close to the expected values of the outlet NOx concentration and ammonia slip as possible.
本发明的遗传算法与传统的迭代算法相比,能很好的避免陷入局部极小的陷阱而出现"死循环"现象,是一种全局优化算法。Compared with the traditional iterative algorithm, the genetic algorithm of the present invention can well avoid falling into the local minimum trap and the "dead loop" phenomenon, and is a global optimization algorithm.
如图3和4所示,步骤103所述根据负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型,采用遗传算法确定使多目标优化函数最优的喷氨量,作为最优喷氨量,具体包括:As shown in FIGS. 3 and 4 , the SCR denitration system output prediction model in the load interval corresponding to the load instruction is used in step 103, and a genetic algorithm is used to determine the ammonia injection amount that optimizes the multi-objective optimization function as the optimal ammonia injection amount, which specifically includes:
步骤401,以预测时间段内每个控制时间点的喷氨量为个体,初始化遗传算法的种群;本发明设定遗传算法的种群大小为300,最大迭代次数(迭代次数阈值)为100,交配概率为0.85,变异概率为0.2。Step 401, taking the ammonia spraying amount at each control time point in the prediction time period as an individual, initialize the population of the genetic algorithm; the present invention sets the population size of the genetic algorithm to 300, the maximum number of iterations (iteration number threshold) to 100, the mating probability to 0.85, and the mutation probability to 0.2.
步骤402,判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果。Step 402, determining whether the valve opening of the SCR denitration system corresponding to the ammonia injection amount at each control time point in each individual in the population is between the lower limit and the upper limit of the valve opening, and obtaining a first determination result.
步骤403,若所述第一判断结果表示否,则将SCR脱硝系统的阀门开度小于阀门开度的下限值的个体中的喷氨量用阀门开度的下限值对应的喷氨量代替,将SCR脱硝系统的阀门开度大于阀门开度的上限值的个体中的喷氨量用阀门开度的上限值对应的喷氨量代替。本发明既考虑了阀门开度的上下限制、也考虑了出口NOx浓度的排放浓度、氨逃逸以及脱硝系统的经济性指标,避免因执行机构饱和从而影响系统性能,能够在保证出口NOx浓度达到目标值的基础上,尽量减少喷氨量,有效降低运行费用和氨逃逸率。Step 403, if the first judgment result indicates no, the amount of ammonia sprayed in the individual of the SCR denitration system whose valve opening is less than the lower limit of the valve opening is replaced by the amount of ammonia sprayed corresponding to the lower limit of the valve opening, and the amount of ammonia sprayed in the individual of the SCR denitration system whose valve opening is greater than the upper limit of the valve opening is replaced by the amount of ammonia sprayed corresponding to the upper limit of the valve opening. The present invention takes into account both the upper and lower limits of the valve opening, the emission concentration of the outlet NOx concentration, ammonia slip, and the economic index of the denitration system, avoids affecting the system performance due to the saturation of the actuator, and can minimize the amount of ammonia sprayed on the basis of ensuring that the outlet NOx concentration reaches the target value, thereby effectively reducing the operating cost and the ammonia slip rate.
步骤404,将种群中每个个体输入SCR脱硝系统输出预测模型,获得每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度和预测出口氨逃逸量。本发明对SCR脱硝系统输出预测模型输出值进行反归一化处理,得到SCR脱硝系统在每个预测时间点的出口NOx浓度以及氨逃逸的预测值。Step 404, input each individual in the population into the SCR denitration system output prediction model to obtain the predicted outlet NOx concentration and the predicted outlet ammonia slip at each prediction time point within the prediction time period corresponding to each individual. The present invention performs a denormalization process on the output value of the SCR denitration system output prediction model to obtain the predicted value of the outlet NOx concentration and ammonia slip of the SCR denitration system at each prediction time point.
步骤405,利用公式对每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度进行校正,获得每个个体对应的预测时间段内每个预测时间点的校正后的预测出口NOx浓度;其中,为k+i·s1时刻校正后的预测出口NOx浓度,为k+i·s1时刻的预测出口NOx浓度,为k时刻的预测出口NOx浓度,为k时刻的SCR脱硝系统的实际出口NOx浓度,i表示第i个预测时间点,s1为预测步长,r为校正系数。本发明的预测时间点的数量为10,控制时间点的数量为3。本发明:k时刻时,通过寻优得到的喷氨量记为u(k-1),可以得到系统实际输出ym(k-1)以及预测模型输出y(k-1)。相应的,通过u(k)可以得到ym(k)和y(k+1)。因为预测模型与实际系统之间不可避免地存在着偏差,所以将k时刻实际输出ym(k)与k时刻模型输出y(k)之间的偏差视为k时刻预测误差的估计值,并将其作为反馈校正信号补偿到k+i时刻的预测模型输出y(k+i)中,即反馈校正后的预测值yp(k+i)为:yp(k+i)=y(k+i)+r(ym(k)-y(k))。反馈校正环节考虑了上一时刻的模型预测误差,可以在一定程度上提高模型预测精度,从而提高预测控制的控制品质。Step 405, using the formula Correct the predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual to obtain the corrected predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual; wherein, is the predicted outlet NOx concentration after correction at time k+i·s1, is the predicted outlet NOx concentration at time k+i·s 1 , is the predicted outlet NOx concentration at time k, is the actual outlet NO x concentration of the SCR denitration system at time k, i represents the i-th prediction time point, s 1 is the prediction step, and r is the correction coefficient. The number of prediction time points of the present invention is 10, and the number of control time points is 3. The present invention: at time k, the amount of ammonia sprayed obtained by optimization is recorded as u(k-1), and the actual output y m (k-1) of the system and the output y(k-1) of the prediction model can be obtained. Correspondingly, y m (k) and y(k+1) can be obtained through u(k). Because there is inevitably a deviation between the prediction model and the actual system, the deviation between the actual output y m (k) at time k and the model output y(k) at time k is regarded as an estimated value of the prediction error at time k, and it is used as a feedback correction signal to compensate the prediction model output y(k+i) at time k+i, that is, the prediction value y p (k+i) after feedback correction is: y p (k+i)=y(k+i)+r(y m (k)-y(k)). The feedback correction link takes into account the model prediction error at the previous moment, which can improve the model prediction accuracy to a certain extent, thereby improving the control quality of predictive control.
步骤406,根据每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度、预测出口氨逃逸量和校正后的预测出口NOx浓度,分别利用第一目标函数和第二目标函数,计算每个个体的第一目标函数值和第二目标函数值。Step 406, according to the predicted outlet NOx concentration, the predicted outlet ammonia slip amount and the corrected predicted outlet NOx concentration at each prediction time point within the prediction time period corresponding to each individual, the first objective function and the second objective function are used to calculate the first objective function value and the second objective function value of each individual.
步骤407,确定种群中满足第二目标函数的个体中第一目标函数最小的个体作为第L次迭代的最优个体,将第L次迭代的最优个体与第L-1次迭代的全局最优个体中第一目标函数值较大的个体设置为第L次迭代的全局最优个体;Step 407, determining the individual with the smallest first objective function among the individuals in the population that satisfy the second objective function as the optimal individual of the Lth iteration, and setting the individual with the larger first objective function value between the optimal individual of the Lth iteration and the global optimal individual of the L-1th iteration as the global optimal individual of the Lth iteration;
步骤408,判断迭代次数是否大于迭代次数阈值,获得第二判断结果。Step 408, determine whether the number of iterations is greater than an iteration number threshold, and obtain a second determination result.
步骤409,若所述第二判断结果表示否,则令迭代次数的数值增加1,采用遗传算法中的遗传、变异和重组的方式更新种群,返回步骤“判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果”。Step 409, if the second judgment result indicates no, then increase the value of the number of iterations by 1, and update the population by the inheritance, mutation and recombination methods in the genetic algorithm, and return to the step of "determining whether the valve opening of the SCR denitrification system corresponding to the ammonia injection amount at each control time point in each individual in the population is between the lower limit and the upper limit of the valve opening, and obtaining the first judgment result".
步骤410,若所述第二判断结果表示是,输出第L次迭代的全局最优个体作为最优喷氨量。Step 410: If the second judgment result indicates yes, output the global optimal individual of the Lth iteration as the optimal ammonia injection amount.
其中,所述第一目标函数为:Wherein, the first objective function is:
其中,J1表示第一目标函数值,为k+i·s1时刻的预测出口氨逃逸量,M2为液氨价格,P为预测时间段的预测时间点的数量,Qgas为烟气流量,为烟气含氧量,M1为排污费价格,为k+j·s2时刻的喷氨量,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,N为机组发电量,M3为电价补贴价格,ω1为第一权重系数,ω2为第二权重系数;Where J 1 represents the first objective function value, is the predicted outlet ammonia slip at time k+i·s1, M2 is the liquid ammonia price, P is the number of predicted time points in the prediction time period, Q gas is the flue gas flow rate, is the oxygen content of flue gas, M1 is the sewage fee price, is the ammonia injection amount at time k+j·s2, j is the jth control time point, s2 is the control step size, M is the number of control time points in the prediction time period, N is the power generation of the unit, M3 is the electricity price subsidy price, ω1 is the first weight coefficient, and ω2 is the second weight coefficient;
所述第二目标函数为:The second objective function is:
其中,J2为第二目标函数值,P为预测时间段的预测时间点的数量,r(k+i·s1)为k+i·s1时刻出口NOx浓度的期望值,||Δu(k+j·s2)||为k+j·s2时刻的喷氨量与期望喷氨量的差值,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,ω3为第三权重系数。Among them, J 2 is the second objective function value, P is the number of prediction time points in the prediction time period, r(k+i·s 1 ) is the expected value of the outlet NO x concentration at moment k+i·s 1 , ||Δu(k+j·s 2 )|| is the difference between the ammonia injection amount at moment k+j·s 2 and the expected ammonia injection amount, j is the jth control time point, s 2 is the control step size, M is the number of control time points in the prediction time period, and ω 3 is the third weight coefficient.
本发明的目标函数中包含了氨逃逸量,烟气流量,烟气含氧量,氮氧化物排污量,排污费价格,氨流量,液氨价格,机组发电量,电价补贴价格。The objective function of the present invention includes ammonia escape, flue gas flow, flue gas oxygen content, nitrogen oxide emission, emission fee price, ammonia flow, liquid ammonia price, unit power generation, and electricity price subsidy price.
步骤105,根据所述最优喷氨量对SCR脱销系统进行控制。Step 105: Control the SCR destocking system according to the optimal ammonia injection amount.
本发明还提供一种全工况电站SCR脱硝系统的控制系统,所述控制系统包括:The present invention also provides a control system for a full-operation power station SCR denitration system, the control system comprising:
历史运行数据获取模块,用于获取SCR脱硝系统的历史运行数据。The historical operation data acquisition module is used to obtain the historical operation data of the SCR denitrification system.
聚类模块,用于采用聚类算法将历史运行数据划分至多个负荷区间,获得多个负荷区间的历史运行数据。The clustering module is used to divide the historical operation data into multiple load intervals by using a clustering algorithm to obtain the historical operation data of the multiple load intervals.
训练模块,用于分别利用每个负荷区间的历史运行数据训练神经网络模型,获得每个负荷区间的SCR脱硝系统输出预测模型;SCR脱硝系统输出预测模型用于根据SCR脱硝系统当前的工作状态数据和预测时间段内的每个控制时间点的喷氨量,预测SCR脱硝系统在预测时间段内的每个预测时间点的出口NOx浓度和出口氨逃逸;所述工作状态数据包括机组负荷、SCR脱硝系统入口的NOx浓度和烟气流量。The training module is used to train the neural network model using the historical operation data of each load interval to obtain the SCR denitration system output prediction model for each load interval; the SCR denitration system output prediction model is used to predict the outlet NOx concentration and outlet ammonia slip of the SCR denitration system at each prediction time point within the prediction time period according to the current working status data of the SCR denitration system and the injection amount of ammonia at each control time point within the prediction time period; the working status data includes the unit load, the NOx concentration at the inlet of the SCR denitration system and the flue gas flow rate.
所述神经网络模型包括输入层、输出层和隐含层,所述输入层包括4个神经元,所述输出层包括2个神经元,所述输出层包括5个神经元。The neural network model includes an input layer, an output layer and a hidden layer, wherein the input layer includes 4 neurons, the output layer includes 2 neurons, and the output layer includes 5 neurons.
SCR脱硝系统输出预测模型选取模块,用于根据SCR脱硝系统的负荷指令,获取负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型。The SCR denitration system output prediction model selection module is used to obtain the SCR denitration system output prediction model for the load interval where the unit load corresponding to the load instruction is located according to the load instruction of the SCR denitration system.
最优喷氨量确定模块,用于根据负荷指令对应的机组负荷所在的负荷区间的SCR脱硝系统输出预测模型,采用遗传算法确定使多目标优化函数最优的喷氨量,作为最优喷氨量。The optimal ammonia injection amount determination module is used to determine the ammonia injection amount that optimizes the multi-objective optimization function according to the SCR denitration system output prediction model in the load range where the unit load corresponding to the load instruction is located, using a genetic algorithm as the optimal ammonia injection amount.
所述最优喷氨量确定模块,具体包括:初始化子模块,用于以预测时间段内每个控制时间点的喷氨量为个体,初始化遗传算法的种群;第一判断子模块,用于判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果;个体更新子模块,用于若所述第一判断结果表示否,则将SCR脱硝系统的阀门开度小于阀门开度的下限值的个体中的喷氨量用阀门开度的下限值对应的喷氨量代替,将SCR脱硝系统的阀门开度大于阀门开度的上限值的个体中的喷氨量用阀门开度的上限值对应的喷氨量代替;预测子模块,用于将种群中每个个体输入SCR脱硝系统输出预测模型,获得每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度和预测出口氨逃逸量;校正子模块,用于利用公式对每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度进行校正,获得每个个体对应的预测时间段内每个预测时间点的校正后的预测出口NOx浓度;其中,为k+i·s1时刻校正后的预测出口NOx浓度,为k+i·s1时刻的预测出口NOx浓度,为k时刻的预测出口NOx浓度,为k时刻的SCR脱硝系统的实际出口NOx浓度,i表示第i个预测时间点,s1为预测步长,r为校正系数;目标函数值计算子模块,用于根据每个个体对应的预测时间段内每个预测时间点的预测出口NOx浓度、预测出口氨逃逸量和校正后的预测出口NOx浓度,分别利用第一目标函数和第二目标函数,计算每个个体的第一目标函数值和第二目标函数值;最优个体确定子模块,用于确定种群中满足第二目标函数的个体中第一目标函数最小的个体作为第L次迭代的最优个体,将第L次迭代的最优个体与第L-1次迭代的全局最优个体中第一目标函数值较大的个体设置为第L次迭代的全局最优个体;第二判断子模块,用于判断迭代次数是否大于迭代次数阈值,获得第二判断结果;种群更新子模块,用于若所述第二判断结果表示否,则令迭代次数的数值增加1,采用遗传算法中的遗传、变异和重组的方式更新种群,返回步骤“判断种群中每个个体中的每个控制时间点的喷氨量对应的SCR脱硝系统的阀门开度是否在阀门开度的下限值和上限值之间,获得第一判断结果”;最优喷氨量输出子模块,用于若所述第二判断结果表示是,输出第L次迭代的全局最优个体作为最优喷氨量。The optimal ammonia injection amount determination module specifically includes: an initialization submodule, which is used to initialize the population of the genetic algorithm with the ammonia injection amount at each control time point within the prediction time period as an individual; a first judgment submodule, which is used to judge whether the valve opening of the SCR denitration system corresponding to the ammonia injection amount at each control time point in each individual in the population is between the lower limit value and the upper limit value of the valve opening, and obtain a first judgment result; an individual update submodule, which is used to replace the ammonia injection amount in the individual whose valve opening of the SCR denitration system is less than the lower limit value of the valve opening with the ammonia injection amount corresponding to the lower limit value of the valve opening, and replace the ammonia injection amount in the individual whose valve opening of the SCR denitration system is greater than the upper limit value of the valve opening with the ammonia injection amount corresponding to the upper limit value of the valve opening; a prediction submodule, which is used to input each individual in the population into the SCR denitration system output prediction model, and obtain the predicted outlet NOx concentration and predicted outlet ammonia slip at each prediction time point within the prediction time period corresponding to each individual; a correction submodule, which is used to use the formula Correct the predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual to obtain the corrected predicted outlet NOx concentration at each prediction time point in the prediction time period corresponding to each individual; wherein, is the predicted outlet NOx concentration after correction at time k+i·s1, is the predicted outlet NOx concentration at time k+i·s 1 , is the predicted outlet NOx concentration at time k, is the actual outlet NOx concentration of the SCR denitrification system at time k, i represents the i-th prediction time point, s1 is the prediction step, and r is the correction coefficient; the objective function value calculation submodule is used to calculate the predicted outlet NOx concentration, predicted outlet ammonia slip and corrected predicted outlet NOx at each prediction time point in the prediction time period corresponding to each individual. x concentration, respectively using the first objective function and the second objective function to calculate the first objective function value and the second objective function value of each individual; an optimal individual determination submodule, used to determine the individual with the smallest first objective function among the individuals in the population that meet the second objective function as the optimal individual of the Lth iteration, and set the individual with the larger first objective function value between the optimal individual of the Lth iteration and the global optimal individual of the L-1th iteration as the global optimal individual of the Lth iteration; a second judgment submodule, used to judge whether the number of iterations is greater than the iteration number threshold, and obtain a second judgment result; a population update submodule, used to increase the value of the number of iterations by 1 if the second judgment result indicates no, update the population by the inheritance, mutation and recombination method in the genetic algorithm, and return to the step of "judging whether the valve opening of the SCR denitrification system corresponding to the ammonia injection amount at each control time point in each individual in the population is between the lower limit and the upper limit of the valve opening, and obtain the first judgment result"; an optimal ammonia injection amount output submodule, used to output the global optimal individual of the Lth iteration as the optimal ammonia injection amount if the second judgment result indicates yes.
控制模块,用于根据所述最优喷氨量对SCR脱销系统进行控制。A control module is used to control the SCR destocking system according to the optimal ammonia injection amount.
其中,所述第一目标函数为:Wherein, the first objective function is:
其中,J1表示第一目标函数值,为k+i·s1时刻的预测出口氨逃逸量,M2为液氨价格,P为预测时间段的预测时间点的数量,Qgas为烟气流量,为烟气含氧量,M1为排污费价格,为k+j·s2时刻的喷氨量,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,N为机组发电量,M3为电价补贴价格,ω1为第一权重系数,ω2为第二权重系数。Where J 1 represents the first objective function value, is the predicted outlet ammonia slip at time k+i·s1, M2 is the liquid ammonia price, P is the number of predicted time points in the prediction time period, Q gas is the flue gas flow rate, is the oxygen content of flue gas, M1 is the sewage fee price, is the amount of ammonia injection at time k+j·s2, j is the jth control time point, s2 is the control step, M is the number of control time points in the prediction time period, N is the power generation of the unit, M3 is the electricity price subsidy price, ω1 is the first weight coefficient, and ω2 is the second weight coefficient.
所述第二目标函数为:The second objective function is:
其中,J2为第二目标函数值,P为预测时间段的预测时间点的数量,r(k+i·s1)为k+i·s1时刻出口NOx浓度的期望值,||Δu(k+j·s2)||为k+j·s2时刻的喷氨量与期望喷氨量的差值,j为第j个控制时间点,s2为控制步长,M为预测时间段内的控制时间点的数量,ω3为第三权重系数。Among them, J 2 is the second objective function value, P is the number of prediction time points in the prediction time period, r(k+i·s 1 ) is the expected value of the outlet NO x concentration at moment k+i·s 1 , ||Δu(k+j·s 2 )|| is the difference between the ammonia injection amount at moment k+j·s 2 and the expected ammonia injection amount, j is the jth control time point, s 2 is the control step size, M is the number of control time points in the prediction time period, and ω 3 is the third weight coefficient.
根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:
1、该方法能够较好地克服SCR脱硝系统大惯性、大迟延的缺点,提高喷氨量控制对机组负荷变化的响应速度,改善SCR脱硝系统的动态调节品质;1. This method can better overcome the shortcomings of large inertia and large delay of the SCR denitrification system, improve the response speed of ammonia injection control to unit load changes, and improve the dynamic regulation quality of the SCR denitrification system;
2、采用多目标控制方式,预测控制中既考虑了阀门开度的上下限制、阀门动作速率限制、也考虑了出口NOx浓度的排放浓度、氨逃逸以及脱硝系统的经济性指标,避免因执行机构饱和从而影响系统性能,能够在保证出口NOx浓度达到目标值的基础上,尽量减少喷氨量,有效降低运行费用和氨逃逸率;2. The multi-objective control method is adopted. The predictive control takes into account the upper and lower limits of the valve opening, the valve action rate limit, the emission concentration of the outlet NOx concentration, ammonia escape and the economic indicators of the denitrification system, so as to avoid the influence of the system performance due to the saturation of the actuator. On the basis of ensuring that the outlet NOx concentration reaches the target value, the amount of ammonia injection can be reduced as much as possible, effectively reducing the operating cost and the ammonia escape rate;
3、采用本发明的预测控制技术,对模型要求较低、易于在线计算、控制效果较好。3. The predictive control technology of the present invention has low model requirements, is easy to calculate online, and has good control effect.
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.
本文中应用了具体个例对发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。This article uses specific examples to illustrate the principles and implementation methods of the invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. The described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011276325.4A CN112418284B (en) | 2020-11-16 | 2020-11-16 | A control method and system for a SCR denitration system in a full-operating power station |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011276325.4A CN112418284B (en) | 2020-11-16 | 2020-11-16 | A control method and system for a SCR denitration system in a full-operating power station |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112418284A CN112418284A (en) | 2021-02-26 |
| CN112418284B true CN112418284B (en) | 2024-08-27 |
Family
ID=74831225
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011276325.4A Active CN112418284B (en) | 2020-11-16 | 2020-11-16 | A control method and system for a SCR denitration system in a full-operating power station |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112418284B (en) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113380338B (en) * | 2021-06-16 | 2022-06-10 | 哈电发电设备国家工程研究中心有限公司 | Method for measuring, correcting and predicting NOx concentration at inlet of cyclone separator |
| CN113433911B (en) * | 2021-06-30 | 2022-05-20 | 浙江大学 | Accurate control system and method for ammonia spraying of denitration device based on accurate concentration prediction |
| CN113485106B (en) * | 2021-07-07 | 2023-03-24 | 山西大学 | Method for controlling concentration of nitrogen oxide in thermal power plant |
| CN114053865A (en) * | 2021-11-03 | 2022-02-18 | 华能铜川照金煤电有限公司 | Generalized predictive control method for SCR denitration control system of coal-fired boiler |
| CN114609986A (en) * | 2022-03-16 | 2022-06-10 | 中国中材国际工程股份有限公司 | Cement decomposing furnace denitration regulation and control optimization system and method based on predictive control |
| CN115309129A (en) * | 2022-10-11 | 2022-11-08 | 华电电力科学研究院有限公司 | SCR denitration efficiency automatic optimization regulation and control method and system |
| CN115524976B (en) * | 2022-10-27 | 2025-07-15 | 东北电力大学 | Ammonia injection adjustment method for SCR system considering boiler combustion state |
| CN116832614B (en) * | 2023-05-24 | 2024-05-24 | 华能国际电力股份有限公司上海石洞口第二电厂 | Ammonia spraying amount control method and system for SCR denitration system |
| CN116880172B (en) * | 2023-06-29 | 2024-10-18 | 华能国际电力股份有限公司上海石洞口第二电厂 | A low-load denitrification ammonia injection optimization control method and system |
| CN117707070A (en) * | 2023-12-13 | 2024-03-15 | 东方电气自动控制工程有限公司 | A PINN-based boiler or gas turbine combustion optimization method and system |
| CN119322458A (en) * | 2024-12-18 | 2025-01-17 | 中电行唐生物质能热电有限公司 | Flue gas denitration self-adaptive control method of granular denitration catalyst |
| CN119717553B (en) * | 2025-02-28 | 2025-06-03 | 华北电力大学 | Control method of denitrification system under all working conditions based on generalized predictive control and intelligent feedforward |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107243257A (en) * | 2017-05-08 | 2017-10-13 | 浙江大学 | It is adapted to the intelligence spray ammonia control system of full load |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7536232B2 (en) * | 2004-08-27 | 2009-05-19 | Alstom Technology Ltd | Model predictive control of air pollution control processes |
| CN106599586B (en) * | 2016-12-19 | 2019-05-28 | 北京国能中电节能环保技术股份有限公司 | SCR neural network based intelligently sprays ammonia optimization method and device |
| CN106681381A (en) * | 2017-01-03 | 2017-05-17 | 华北电力大学 | SCR denitration system ammonia spraying quantity optimal control system and method based on intelligent feedforward signals |
| CN107544288B (en) * | 2017-09-26 | 2021-06-04 | 天津拓科思科技有限公司 | Denitration optimization control method and system |
| CN109062053B (en) * | 2018-08-31 | 2022-11-29 | 江苏国信靖江发电有限公司 | Denitration ammonia injection control method based on multivariate correction |
| CN111581581B (en) * | 2020-04-23 | 2021-02-02 | 大唐环境产业集团股份有限公司 | Method and system for detecting NOx concentration at SCR inlet under multi-boundary condition |
| CN111522290A (en) * | 2020-04-24 | 2020-08-11 | 大唐环境产业集团股份有限公司 | Denitration control method and system based on deep learning method |
| CN111921377B (en) * | 2020-09-25 | 2021-01-26 | 大唐环境产业集团股份有限公司 | SCR denitration ammonia injection control method and system based on mechanism and data driving |
-
2020
- 2020-11-16 CN CN202011276325.4A patent/CN112418284B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107243257A (en) * | 2017-05-08 | 2017-10-13 | 浙江大学 | It is adapted to the intelligence spray ammonia control system of full load |
Non-Patent Citations (1)
| Title |
|---|
| SCR脱硝系统多目标优化控制研究;杨婷婷等;《中国电机工程学报》;第41卷(第14期);4905-4911 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112418284A (en) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112418284B (en) | A control method and system for a SCR denitration system in a full-operating power station | |
| CN107243257B (en) | It is suitble to the intelligence spray ammonia control system of full load | |
| CN104826492B (en) | Improvement method for selective catalytic reduction flue gas denitrification and ammonia injection control system | |
| CN105629736B (en) | The fired power generating unit SCR denitration Disturbance Rejection forecast Control Algorithm of data-driven | |
| CN112580250A (en) | Thermal power generating unit denitration system based on deep learning and optimization control method | |
| CN104826493B (en) | A kind of control method of selective catalytic reduction flue gas denitrification system | |
| CN110368808A (en) | A kind of the ammonia spraying amount control method and system of SCR flue gas denitrification system | |
| CN106599586B (en) | SCR neural network based intelligently sprays ammonia optimization method and device | |
| CN113433911B (en) | Accurate control system and method for ammonia spraying of denitration device based on accurate concentration prediction | |
| CN106681381A (en) | SCR denitration system ammonia spraying quantity optimal control system and method based on intelligent feedforward signals | |
| CN111968708B (en) | SCR denitration ammonia injection amount prediction method based on random forest and LSTM neural network | |
| CN109062053A (en) | A kind of denitration spray ammonia control method based on multivariate calibration | |
| CN111522290A (en) | Denitration control method and system based on deep learning method | |
| CN111897373B (en) | Model prediction-based ammonia injection flow adjusting method for SCR denitration device | |
| CN110188383A (en) | A Modeling Method for Power Station SCR Denitrification Based on Selective Integration Model | |
| CN108803309A (en) | It is a kind of that ammonia optimization method and system are intelligently sprayed based on the SCR denitration of hard measurement and model adaptation | |
| CN111460726B (en) | Optimization method for ammonia escape of coal slime fluidized bed boiler denitration system | |
| CN115524976B (en) | Ammonia injection adjustment method for SCR system considering boiler combustion state | |
| CN110299188A (en) | SCR flue gas denitrification system GRNN modeling method based on GA variables choice | |
| CN115145152A (en) | A method for synergistic optimal control of boiler combustion and denitration process | |
| CN106842955A (en) | CO after burning with exhaust gas volumn Disturbance Rejection2Trapping system forecast Control Algorithm | |
| CN111013370A (en) | Wet desulphurization slurry supply amount prediction method based on deep neural network | |
| CN116036849A (en) | A CFB boiler flue gas denitrification automatic control method and system | |
| CN109933884B (en) | A Neural Network Inverse Control Method for SCR Denitrification System of Coal-fired Units | |
| CN109833773A (en) | A kind of NO_x Reduction by Effective ammonia flow accuracy control method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |