CN113009890B - 一种基于bmo-lssvm的火电厂scr烟气脱硝控制系统的优化控制方法 - Google Patents
一种基于bmo-lssvm的火电厂scr烟气脱硝控制系统的优化控制方法 Download PDFInfo
- Publication number
- CN113009890B CN113009890B CN202110109418.6A CN202110109418A CN113009890B CN 113009890 B CN113009890 B CN 113009890B CN 202110109418 A CN202110109418 A CN 202110109418A CN 113009890 B CN113009890 B CN 113009890B
- Authority
- CN
- China
- Prior art keywords
- lssvm
- bmo
- flue gas
- scr
- barnacle
- 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
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 title claims abstract description 68
- 239000003546 flue gas Substances 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 56
- 108700018263 Brassica oleracea SCR Proteins 0.000 title claims abstract description 13
- 238000005457 optimization Methods 0.000 title description 35
- 238000012549 training Methods 0.000 claims abstract description 32
- 238000012360 testing method Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000010606 normalization Methods 0.000 claims abstract description 7
- HKSZLNNOFSGOKW-HMWZOHBLSA-N staurosporine Chemical compound C12=C3N4C5=CC=CC=C5C3=C3CNC(=O)C3=C2C2=CC=CC=C2N1[C@@H]1C[C@H](NC)[C@H](OC)[C@@]4(C)O1 HKSZLNNOFSGOKW-HMWZOHBLSA-N 0.000 claims abstract description 4
- 241000238586 Cirripedia Species 0.000 claims description 73
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 claims description 52
- 229910021529 ammonia Inorganic materials 0.000 claims description 25
- 230000013011 mating Effects 0.000 claims description 15
- 239000000243 solution Substances 0.000 claims description 14
- 210000004392 genitalia Anatomy 0.000 claims description 12
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 8
- 230000004720 fertilization Effects 0.000 claims description 8
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 229910000069 nitrogen hydride Inorganic materials 0.000 claims description 4
- 241000894007 species Species 0.000 claims description 2
- 238000001035 drying Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 33
- 238000012843 least square support vector machine Methods 0.000 description 28
- BIGPRXCJEDHCLP-UHFFFAOYSA-N ammonium bisulfate Chemical compound [NH4+].OS([O-])(=O)=O BIGPRXCJEDHCLP-UHFFFAOYSA-N 0.000 description 6
- 239000003054 catalyst Substances 0.000 description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N nitrogen Substances N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- XYXNTHIYBIDHGM-UHFFFAOYSA-N ammonium thiosulfate Chemical compound [NH4+].[NH4+].[O-]S([O-])(=O)=S XYXNTHIYBIDHGM-UHFFFAOYSA-N 0.000 description 2
- 238000010531 catalytic reduction reaction Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- VHUUQVKOLVNVRT-UHFFFAOYSA-N Ammonium hydroxide Chemical compound [NH4+].[OH-] VHUUQVKOLVNVRT-UHFFFAOYSA-N 0.000 description 1
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 description 1
- 230000002378 acidificating effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 1
- 235000011114 ammonium hydroxide Nutrition 0.000 description 1
- BFNBIHQBYMNNAN-UHFFFAOYSA-N ammonium sulfate Chemical compound N.N.OS(O)(=O)=O BFNBIHQBYMNNAN-UHFFFAOYSA-N 0.000 description 1
- 229910052921 ammonium sulfate Inorganic materials 0.000 description 1
- 235000011130 ammonium sulphate Nutrition 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000004202 carbamide Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 230000009849 deactivation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000010196 hermaphroditism Effects 0.000 description 1
- 230000002779 inactivation Effects 0.000 description 1
- 230000007758 mating behavior Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000000276 sedentary effect Effects 0.000 description 1
- 238000007086 side reaction Methods 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000009182 swimming Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- 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
-
- 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/8696—Controlling the catalytic process
-
- 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/90—Injecting reactants
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32339—Object oriented modeling, design, analysis, implementation, simulation language
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02C—CAPTURE, STORAGE, SEQUESTRATION OR DISPOSAL OF GREENHOUSE GASES [GHG]
- Y02C20/00—Capture or disposal of greenhouse gases
- Y02C20/10—Capture or disposal of greenhouse gases of nitrous oxide (N2O)
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Biomedical Technology (AREA)
- Analytical Chemistry (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Health & Medical Sciences (AREA)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Treating Waste Gases (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提供了一种基于BMO‑LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,包括如下步骤:Step1:动态采集包含SCR烟气脱硝系统模型的输入指标输出指标;Step2:对样本数据集进行去燥处理和归一化处理,分成训练集和测试集:Step3:建立LSSVM模型,以训练集LSSVM模型进行训练;Step4:使用BMO对LSSVM模型的核函数参数和惩罚系数进行寻优,并将获得的最优核函数参数σ以及惩罚系数C赋给LSSVM模型,得到BMO‑LSSVM模型;Step5:用训练集P01对BMO‑LSSVM模型进行训练;Stsp6:利用测试集检验,若超出设定的阈值,则重复Step4‑6;Step7:根据预测结果,及时调整SCR反应器出口的NOx浓度。本申请在保证脱硝效率的前提下,以解决现火电厂SCR脱硝出口NOx浓度波动过大,且短时刻时有超出最低排放标准的现象的问题。
Description
技术领域
本发明涉及燃煤电站烟气脱硝领域,具体涉及一种基于BMO算法优化LSSVM的SCR烟气脱硝控制系统优化方法。
背景技术
最新修订的《火电厂大气污染物排放标准》对NOx的排放制定了非常严格的要求,要求 NOx排放在50mg/m3以内。只依靠低氮燃烧不能达到排放指标要求,因此燃烧后控制NOx,即在烟道尾部加装脱硝装置是目前大型燃煤电厂的普遍选择。
在脱硝技术中,SCR(选择性催化还原)因其脱硝效率高、烟气处理量大、技术成熟成为燃煤电站的主要选择。SCR是指在催化剂和氧气的存在下,在320℃~427℃温度范围下,还原剂(无水氨、氨水或尿素)有选择性地与烟气中的NOx反应生成无害的氮和水,从而去除烟气中的NOx。
目前火电厂普遍存在脱硝出口NOx浓度波动过大,且短时刻时有超出最低排放标准,即 50mg/m3的情况发生。在短时刻出现超出排放标准的情况下,电厂工作人员常手动调节喷氨阀门,以控制反应器出口NOx浓度不超标,但此时即会造成喷氨量过多,导致氨逃逸问题加剧,同时增加副反应产生的可能,生成过多的具有强腐蚀性以及强粘性的硫酸氢铵以及次硫酸铵,加剧空预器堵塞问题,严重影响整个脱硝机组的安全稳定运行,并且增加电厂经济运行成本。
发明内容
本申请提供一种基于BMO-LSSVM的火电厂SCR脱硝控制系统的优化控制方法,在保证脱硝效率的前提下,以解决现火电厂SCR脱硝出口NOx浓度波动过大,且短时刻时有超出最低排放标准的现象的问题。
为解决上述技术问题,本发明采用的技术方案为:
一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,
包括如下步骤:
Step1:动态采集包含SCR烟气脱硝系统模型的输入指标X和输出指标Y的样本数据集XY(X,Y),动态收集的间隔时间为t;
Step2:对样本数据集XY(X,Y)进行去噪处理,形成去噪数据集P,对去噪数据集P进行归一化处理形成归一化数据集P0,将归一化数据集P0分成训练集P01和测试集记为P02;
Step3:建立LSSVM模型,以训练集P01对LSSVM模型进行训练,得到基于LSSVM的 SCR烟气脱硝控制系统;
Step4:使用BMO对LSSVM模型的核函数参数σ和惩罚系数C进行寻优,并将获得的最优核函数参数σ以及惩罚系数C赋给LSSVM模型,得到BMO-LSSVM模型;
Step5:用训练集P01对BMO-LSSVM模型进行训练,得到基于BMO-LSSVM的SCR烟气脱硝控制模型;
Stsp6:利用测试集P02对基于BMO-LSSVM的SCR烟气脱硝控制模型的预测效果进行检验,若SCR反应器出口NOx浓度的预测值与实际值之间的误差在设定的阈值内,将基于BMO-LSSVM的SCR烟气脱硝控制模型用于优化SCR烟气脱硝控制系统,若超出设定的阈值,则重复Step4-6;
Step7:根据预测的SCR反应器出口NOx浓度y获取最佳喷氨控制量,及时调整SCR反应器出口的NOx浓度。
优选的,其中,step1中输入指标X包括机组负荷x1、SCR入口NOx浓度x2、入口烟气流量x3、入口烟气NH3/NOx摩尔比x4、入口烟气温度x5、加氨流量x6、出口NOx设定值x7、氨逃逸量x8以及前一时刻出口NOx浓度x9,输出指标Y为SCR反应器出口 NOx浓度y。
优选的,Step2中,对数据进行归一化处理采用如下公式:
P0=(P-Pmin)/(Pmax-Pmin) (1)
式中,P0为去噪数据集P归一化后得到的数据集,P为去噪数据集,Pmax、Pmin分别为去燥数据集P中的最大值和最小值。
进一步的,Step2中对数据去噪操作包括:
去除数据集中的机组负荷低于满负荷后25%的数据,去除数据集中存在缺失值的无效数据,以及机组负荷不在负荷范围内的错误数据。
其中,Step3中LSSVM的模型采用下式:
其中,K(x,xi)是核函数,用来代替高维空间上的内积运算,xi是核函数中心,x是训练样本的输入值,αi表示拉格朗日乘子,αi>0,i=1,2,...,N,b是偏置常数。
优选的,K(x,xi)选用RBF核函数:
其中,σ为核函数参数。
其中,Step4包括如下步骤:
Step4-1:初始化藤壶种群D,设置藤壶的数量并设置最大迭代次数;假设候选解为藤壶,其中种群向量表达为:
其中,D是初始群种,N为控制变量个数,n为种群大小即藤壶数量。每个控制变量都有相应的上下界约束:
其中,ubi,lbi分别表示第i个控制变量的上界和下界;
Step4-2:对初始种群D进行评估,基于适应度函数计算每个藤壶的适应度值,然后对各个解进行排序,当前最优解在D的最顶部,设T为当前最优解,适应度函数为:
Step4-3:设置生殖器长度pl;
Step4-4:进行交配父代的选择:
barnacle_d=rand perm(n) (7)
barnacle_m=rand perm(n) (8)
式中,barnacle_d和barnacle_m为交配的父代和母代,n为种群的大小,rand perm(n)表示随机选取整个藤壶种群中的一个藤壶;
Step4-5:从藤壶的父代产生新的子代;
如果选择z≤pl,选择的藤壶是在父亲的藤壶生殖器长度的范围内,发生了利用过程,根据子代生殖公式产生新的子代:
如果选择z>pl,即当要交配的藤壶的选择超出了预期的pl值时,依据远程受精公式发生远程受精:
Step4-6:调整每个控制变量的边界,计算每个藤壶的适应度值,并进行排序,更新T,若达到最大迭代次数,则输出全局最优解,即得到最优的核函数参数σ和惩罚系数C;否则重复Step4-4—Step4-6。
优选的,Step4-1中藤壶数量为60,最大迭代次数为100。
其中,Step4-5中子代生殖公式为:
其中,Step4-5中远程受精公式为:
式中,rand()为[0,1]内的随机数。
本发明的有益效果:
本申请提供一种基于BMO-LSSVM的火电厂SCR脱硝控制系统的优化控制方法,在保证脱硝效率的前提下,以解决现火电厂SCR脱硝出口NOx浓度波动过大,且短时刻时有超出最低排放标准的现象的问题。
(1)本发明利用藤壶交配优化算法(BMO)优化最小二乘支持向量机(LSSVM)建立SCR烟气脱硝控制系统的优化控制模型,该方法一定程度上克服了传统PID控制的不及时,造成短时刻时有超出最低排放标准的情况的发生,该方法利用核函数将数据映射到高维空间进行回归,考虑了数据间的非线性关系;
(2)本发明通过BMO算法优化LSSVM,选取最优核函数参数以及惩罚系数,使得模型优化效果达到最优,准确的预测反应器出口NOx浓度,及时准确的获取最佳喷氨控制量,从而优化出口NOx浓度,算法可消除原LSSVM模型中核函数参数以及惩罚系数随机选择或人为设置对模型优化效果的影响,稳定性好、精度高。
(3)本发明提供一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,预测SCR出口NOx浓度,及时准确的获取最佳喷氨控制量,从而优化出口NOx浓度,减少电厂人员在短时刻出现超出最低排放标准时,选择手动调节喷氨调节阀,加大喷氨量,从而增加具有强酸性和强粘性的硫酸氢铵(NH4HSO4)以及次硫酸铵((NH4)2SO2) 的生成的情况的发生,减少催化剂失活的问题以及空预器堵塞的情况,促进整个烟气脱硝系统的安全稳定运行。
附图说明
图1为本发明提供的基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化方法流程示意图;
图2为本发明提供的BMO算法优化LSSVM的核函数参数σ和惩罚系数C的流程示意图;
图3为本发明提供的BMO算法在藤壶数目为10时的藤壶交配过程的选择过程示意图;
图4为本发明提供的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化方法选取某电厂800组实际历史运行数据,优化前后的反应器出口NOx浓度曲线对比图;
图5为本发明提供的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化方法在变负荷情况下,优化前后的反应器出口NOx浓度曲线对比图;
图6为本发明提供的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化方法在变负荷情况下,优化前后的脱硝效率对比曲线。
具体实施方式
下面结合附图以及具体实施方法对本发明一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法作进一步详细说明。
如图1所示,一种基于BMO-LSSVM的火电厂SCR脱硝控制系统优化方法,包括如下步骤:
Step1:选取样本数据信息:根据SCR烟气脱硝工程的实际情况并经理论分析,选取机组负荷x1、SCR入口NOx浓度x2、入口烟气流量x3、入口烟气NH3/NOx摩尔比x4、入口烟气温度x5、加氨流量x6、出口NOx设定值x7、氨逃逸量x8以及前一时刻出口NOx浓度x9共九个因素作为SCR烟气脱硝系统模型的输入量,SCR反应器出口NOx浓度y作为输出量; Step2:选取并处理数据集:从某电厂实际历史数据中选取样本数据,采样间隔为1min,经分析处理,使用python去除数据集中的机组负荷低于满负荷后25%的数据,去除数据集中存在缺失值的无效数据,以及机组负荷不在负荷范围内的错误数据,选取具有代表性的数据构成数据集P,该组数据涉及低负荷(满负荷的25%开始)、中负荷(满负荷的50%开始)、高负荷 (满负荷的75%开始)三种负荷状态,包含多个稳定状态和多个动态状态,稳定状态即某段时间内稳定于某个负荷的状态,动态状态即负荷不断变化的状态。将数据集P分为训练集和测试集,并将数据集进行归一化处理得到数据集P0,将归一化后的训练集记为P01,归一化后的测试集记为P02;
Step3:建立LSSVM优化控制模型,以处理好的训练集P01对LSSVM模型进行训练,训练好的LSSVM模型可用于优化SCR烟气脱硝控制系统;
Step4:使用藤壶交配优化算法(BMO)对最小二乘支持向量机(LSSVM)的核函数参数σ和惩罚系数C进行寻优,并将获得的最优核函数参数以及惩罚系数赋给LSSVM,得到BMO-LSSVM优化控制模型;
Step5:用处理好的训练集P01重新训练经优化得到的BMO-LSSVM模型;
Stsp6:利用处理好的测试集P02对训练好的BMO-LSSVM模型的预测效果进行检验,计算预测结果与实际测量值之间的误差,若在设定的误差范围内,则可将训练好的模型用于优化SCR烟气脱硝控制系统,若超出设定的误差范围,则重复Step4-6;
Step7:根据预测的出口NOx浓度及时准确的获取最佳喷氨控制量,对SCR出口NOx浓度进行优化。
具体的,所述Step2中,对数据进行归一化处理采用如下公式:
P0=(P-Pmin)/(Pmax-Pmin) (1)
式中,P0为归一化后得到的数据集,P为原始数据集,Pmax、Pmin分别为原始数据集中的最大值和最小值,最终将样本点数据映射到[0,1]范围之间。
具体的,所述Step3中,包括如下具体步骤,首先构建决策函数如下:
在n维空间中,将最优化问题表示为:
其中,ei是松弛变量,C是惩罚系数。利用上式构建Lagrange函数,如下所示:
其中,αi>0(i=1,2,...,N)是拉格朗日乘子。
根据KKT最优条件:
消去ω和ei,可得到:
其中,Q=[1,...,1]T,Ki=K(x,xi)为矩阵K中的元素。将得到LSSVM的函数估计:
其中,K(x,xi)是核函数,用来代替高维空间上的内积运算,xi是核函数中心,x是训练样本的输入值,本文选用结构简单且泛化性能好的RBF为核函数,即:
其中,σ为核函数参数。
藤壶交配优化(BMO)算法模拟自然界中藤壶的交配行为来解决优化问题。藤壶是自侏罗纪时代就存在的微生物,出生时即会游泳,当成年时,它们会附着在水中的物体上,并长出壳。大多数藤壶雌雄同体,这意味着它们有雄性和雌性繁殖。为了应对潮汐变化和久坐不动的生活方式,它们的生殖器可以达到身长的7~8倍,生殖器所能接触到的所有邻居和潜在的配偶竞争者就构成了其交配种群,因此生殖器长度变化在确定交配种群大小和局部交配竞争中起到重要作用。
如图2所示,所述Step4中,BMO算法优化LSSVM的核函数参数σ和惩罚系数C的具体步骤如下:
Step4-1:初始化藤壶种群Di,设置藤壶的数量并设置最大迭代次数;
本实施例中藤壶数量设为60,最大迭代次数设为100。
假设候选解为藤壶,其中种群向量可以表达为:
其中,N为控制变量个数,即问题维度,本方法需优化LSSVM模型的核函数参数σ和惩罚系数C,故控制变量个数为2,n为种群大小即藤壶数量。每个控制变量都有相应的上下界约束:
其中,ubi,lbi分别表示第i个控制变量的上界和下界。
Step4-2:对初始种群D进行评估,计算每个藤壶的适应度值,然后对各个解进行排序,当前最优解在D的最顶部,设T为当前最优解。适应度函数为:
Step4-3:设置生殖器长度pl,本实施例中将藤壶的最大生殖器长度设为其体长的七倍,即pl=7,那么在某一次迭代中,藤壶1只能与藤壶2~7中的一个进行交配,如果藤壶1选择了藤壶8,那么将不会执行正常的交配过程,子代的产生是通过远程受精过程进行的,为方便展示,如图3所示,展示了10个藤壶交配过程的选择;
Step4-4:根据下式进行选择:
barnacle_d=rand perm(n) (7)
barnacle_m=rand perm(n) (8)
式中,barnacle_d和barnacle_m为交配的父代,n为种群的大小,rand perm(n)表示随机选取整个藤壶种群中的一个藤壶。
Step4-5:从藤壶的父代产生新的子代。
如果选择z≤pl,即选择的藤壶是在父亲的藤壶生殖器长度的范围内,就发生了利用过程,根据下式产生新的子代变量:
式中,p为[0,1]内均匀分布的随机数,q=1-p,和分别为等式(7)和(8) 选择的父亲藤壶和母亲藤壶的变量,p和q代表了父亲和母亲的特征在下一代中所占的百分比。因此,子代根据0到1之间的随机数概率继承父亲和母亲的行为。
如果选择z>pl,即当要交配的藤壶的选择超出了预期的pl值时,就会发生远程受精:
式中,rand()为[0,1]内的随机数。
等式(10)为藤壶子代进化的简单方法,新的子代根据母亲藤壶生成以进行探索,这是由于新的子代是由母亲藤壶产生的,因为它接受了水中由其他藤壶释放的精子。
Step4-6:调整每个变量的边界,计算每个藤壶的适应度值,并进行排序,更新T,若达到最大迭代次数,则输出全局最优解,即得到最优的核函数参数σ和惩罚系数C,否则重复Step4-4—Step4-6。
如图4所示,为本实施例选取某电厂具有代表性的800组实际历史运行数据进行实验,使用BMO-LSSVM优化后的SCR反应器出口NOx浓度曲线,并与原传统PID控制下反应器出口NOx浓度曲线进行对比。
PID控制下,短时刻SCR出口NOx浓度出现较多超出50mg/m3的情况,最高达58.94mg/m3,相比较之下,采用BMO-LSSVM模型进行优化后,SCR反应器出口NOx浓度短时刻未出现超出50mg/m3的情况,最高为41.6mg/m3。
如图5所示,为本实施例在变负荷下的,采用BMO-LSSVM优化后的SCR出口NOx浓度,并与原传统PID控制下进行对比。
火电厂SCR烟气脱硝系统的整个控制过程具有惯性大、非线性、大滞后性等特征,反应机理十分复杂,使用传统的PID方法进行控制的效果存在一定的问题。
在变工况下,当机组负荷为300MW时,SCR出口NOx浓度为47.1mg/m3,即有超出最低排放标准的趋势,在机组负荷为350MW时,SCR反应器出口浓度已超标,为54.8mg/m3,而环保部门是对出口净烟气NOx浓度进行考核,整个尾气处理的系统存在大滞后的特点,那么在环保部门对出口净烟气NOx浓度进行监测时,将很有可能出现无法达到排放标准的情况,那么在实际情况下,工作人员常在将出现超出排放标准的时候,加大喷氨量,则增加了氨逃逸的可能,并会增加NH4HSO4以及(NH4)2SO2的生成,从而影响整个脱硝机组的的安全问题。在机组负荷增加时,BMO-LSSVM模型依旧可以较好的进行控制,在高负荷350MW 时,SCR反应器出口NOx浓度为36.1mg/m3,可见,采用BMO-LSSVM模型进行控制,可以更好的达到脱硝的目的。
如图6所示,为本实施例在变负荷下的,采用BMO-LSSVM优化后的脱硝效率,并与原传统PID控制下进行对比。在变工况下,传统PID控制在低负荷100WM时,脱硝效率为88.89%,在高负荷350MW时,脱硝效率仅为67.77%,而采用BMO-LSSVM优化控制模型进行优化后,低负荷100MW时,脱硝效率可提升至96.74%,在高负荷350MW时,脱硝效率可提升至88.91%。
本发明需一个或多个处理器、存储器以及程序,其中存储器需存储上述所需数据,且一个或多个程序存储在其中,并由处理器执行,一个或多个程序用于执行基于BM0算法优化 LSSVM的SCR烟气脱硝控制系统的优化控制方法的执行。
本发明提供一种基于藤壶交配优化算法(BMO)优化的最小二乘支持向量机(LSSVM) 的SCR烟气脱硝控制系统的优化控制方法,选取某电厂实际历史运行数据为样本数据集,采样间隔为1min,经分析处理,去除无效、错误数据,选取具有代表性的数据构成数据集P,该组数据涉及低、中、高三种负荷状态,包含多个稳定状态和多个动态状态。并将数据集P 分为训练集和测试集,并将数据集进行归一化处理得到数据集P0,将归一化后的训练集记为 P01,归一化后的测试集记为P02。根据SCR烟气脱硝工程的实际情况并经理论分析,选取机组负荷x1、SCR入口NOx浓度x2、入口烟气流量x3、入口烟气NH3/NOx摩尔比x4、入口烟气温度x5、加氨流量x6、出口NOx设定值x7、氨逃逸量x8以及前一时刻出口NOx浓度x9共九个因素作为SCR烟气脱硝系统模型的输入量,SCR反应器出口NOx浓度y作为输出量。之后建立LSSVM模型,并使用训练集P01进行训练;采用BMO算法对LSSVM的的核函数参数σ和惩罚系数C进行寻优,并将获得的最优核函数参数σ以及惩罚系数C赋给LSSVM,得到 BMO-LSSVM优化控制模型;并使用训练集P01对BMO-LSSVM模型进行训练,利用处理好的测试集P02对训练好的BMO-LSSVM模型的预测效果进行检验,计算预测结果与实际测量值之间的误差,若在设定的误差范围内,则可将训练好的模型用于优化SCR烟气脱硝控制系统,若超出设定的误差范围,则重新进行核函数参数以及惩罚系数的寻优,以及后续操作,即重复上述Step4-6;之后根据预测的出口NOx浓度及时准确的获取最佳喷氨控制量,对SCR 出口NOx浓度进行优化。
本发明表明可以通过BMO算法优化LSSVM对火电厂SCR烟气脱硝控制系统的优化控制,采用电厂实际历史运行数据进行训练并测试,所选数据集涉及低、中、高三种负荷状态,包含多个稳定状态和多个动态状态,可以为目前的SCR脱硝过程中,反应器出口NOx浓度波动过大,短时刻时有超出最低排放指标的状况发生的问题提供一个解决方法。本发明利用BMO算法优化LSSVM建立SCR烟气脱硝控制系统的优化控制模型,该方法一定程度上克服了传统PID控制的不及时,造成短时刻时有超出最低排放标准的情况的发生,该方法利用核函数将数据映射到高维空间进行回归,考虑了数据间的非线性关系;选取最优核函数参数以及惩罚系数,使得模型优化效果达到最优,准确的预测反应器出口NOx浓度,及时准确的获取最佳喷氨控制量,从而优化出口NOx浓度,算法可消除原LSSVM模型中核函数参数以及惩罚系数随机选择或人为设置对模型优化效果的影响,稳定性好、精度高;
本发明提供一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,预测SCR出口NOx浓度,及时准确的获取最佳喷氨控制量,从而优化出口NOx浓度,减少电厂人员在短时刻出现超出最低排放标准时,选择手动调节喷氨调节阀,加大喷氨量,从而增加具有强酸性和强粘性的硫酸氢铵(NH4HSO4)以及次硫酸铵((NH4)2SO2)的生成的情况的发生,减少催化剂失活的问题以及空预器堵塞的情况,促进整个火电厂的烟气脱硝系统的安全稳定运行。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下被实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
在本发明的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
类似地,应当理解,为了精简本公开并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多特征。更确切地说,如权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
如在此所使用的那样,除非另行规定,使用序数词“第一”、“第二”、“第三”等等来描述普通对象仅仅表示涉及类似对象的不同实例,并且并不意图暗示这样被描述的对象必须具有时间上、空间上、排序方面或者以任意其它方式的给定顺序。
尽管根据有限数量的实施例描述了本发明,但是受益于上面的描述,本技术领域内的技术人员明白,在由此描述的本发明的范围内,可以设想其它实施例。此外,应当注意,本说明书中使用的语言主要是为了可读性和教导的目的而选择的,而不是为了解释或者限定本发明的主题而选择的。因此,在不偏离所附权利要求书的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。对于本发明的范围,对本发明所做的公开是说明性的,而非限制性的,本发明的范围由所附权利要求书限定。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,包括如下步骤:
Step1:动态采集包含SCR烟气脱硝系统模型的输入指标X和输出指标Y的样本数据集XY(X,Y),动态收集的间隔时间为t;
Step2:对样本数据集XY(X,Y)进行去噪处理,形成去噪数据集P,对去噪数据集P进行归一化处理形成归一化数据集P0,将归一化数据集P0分成训练集P01和测试集记为P02;
Step3:建立LSSVM模型,以训练集P01对LSSVM模型进行训练,得到基于LSSVM的SCR烟气脱硝控制系统;
Step4:使用BMO对LSSVM模型的核函数参数σ和惩罚系数C进行寻优,并将获得的最优核函数参数σ以及惩罚系数C赋给LSSVM模型,得到BMO-LSSVM模型;
Step5:用训练集P01对BMO-LSSVM模型进行训练,得到基于BMO-LSSVM的SCR烟气脱硝控制模型;
Stsp6:利用测试集P02对基于BMO-LSSVM的SCR烟气脱硝控制模型的预测效果进行检验,若SCR反应器出口NOx浓度的预测值与实际值之间的误差在设定的阈值内,将基于BMO-LSSVM的SCR烟气脱硝控制模型用于优化SCR烟气脱硝控制系统,若超出设定的阈值,则重复Step4-6;
Step7:根据预测的SCR反应器出口NOx浓度y获取最佳喷氨控制量,及时调整SCR反应器出口的NOx浓度。
2.根据权利要求1所述的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,Step1 中输入指标X包括机组负荷x1、SCR入口NOx浓度x2、入口烟气流量x3、入口烟气NH3/NOx摩尔比x4、入口烟气温度x5、加氨流量x6、出口NOx设定值x7、氨逃逸量x8以及前一时刻出口NOx浓度x9,输出指标Y为SCR反应器出口NOx浓度y。
3.根据权利要求1所述的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,Step2中,对数据进行归一化处理采用如下公式:
P0=(P-Pmin)/(Pmax-Pmin) (1)
式中,P0为去噪数据集P归一化后得到的数据集,P为去噪数据集,Pmax、Pmin分别为去燥数据集P中的最大值和最小值。
4.根据权利要求1所述的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,Step2中对数据去噪操作包括:
去除数据集中的机组负荷低于满负荷后25%的数据,去除数据集中存在缺失值的无效数据,以及机组负荷不在负荷范围内的错误数据。
7.根据权利要求1所述的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,Step4包括如下步骤:
Step4-1:初始化藤壶种群D,设置藤壶的数量并设置最大迭代次数;假设候选解为藤壶,其中种群向量表达为:
其中,D是初始群种,N为控制变量个数,n为种群大小即藤壶数量, 每个控制变量都有相应的上下界约束:
其中,ubi,lbi分别表示第i个控制变量的上界和下界;
Step4-2:对初始种群D进行评估,基于适应度函数计算每个藤壶的适应度值,然后对各个解进行排序,当前最优解在D的最顶部,设T为当前最优解,适应度函数为:
Step4-3:设置生殖器长度pl;
Step4-4:进行交配父代的选择:
barnacle_d=rand perm(n) (7)
barnacle_m=rand perm(n) (8)
式中,barnacle_d和barnacle_m为交配的父代和母代,n为种群的大小,rand perm(n)表示随机选取整个藤壶种群中的一个藤壶;
Step4-5:从藤壶的父代产生新的子代;
如果选择z≤pl,选择的藤壶是在父亲的藤壶生殖器长度的范围内,发生了利用过程,根据子代生殖公式产生新的子代:
如果选择z>pl,即当要交配的藤壶的选择超出了预期的pl值时,依据远程受精公式发生远程受精:
Step4-6:调整每个控制变量的边界,计算每个藤壶的适应度值,并进行排序,更新T,若达到最大迭代次数,则输出全局最优解,即得到最优的核函数参数σ和惩罚系数C;否则重复Step4-4—Step4-6。
8.根据权利要求7所述的一种基于BMO-LSSVM的火电厂SCR烟气脱硝控制系统的优化控制方法,其特征在于,Step4-1中藤壶数量为60,最大迭代次数为100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110109418.6A CN113009890B (zh) | 2021-01-25 | 2021-01-25 | 一种基于bmo-lssvm的火电厂scr烟气脱硝控制系统的优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110109418.6A CN113009890B (zh) | 2021-01-25 | 2021-01-25 | 一种基于bmo-lssvm的火电厂scr烟气脱硝控制系统的优化控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113009890A CN113009890A (zh) | 2021-06-22 |
CN113009890B true CN113009890B (zh) | 2021-12-31 |
Family
ID=76384960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110109418.6A Active CN113009890B (zh) | 2021-01-25 | 2021-01-25 | 一种基于bmo-lssvm的火电厂scr烟气脱硝控制系统的优化控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113009890B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113554225B (zh) * | 2021-07-22 | 2023-09-26 | 浙江理工大学 | 一种基于藤壶交配算法的流水车间人力调度优化方法 |
CN113611375B (zh) * | 2021-08-09 | 2023-10-10 | 成都佳华物链云科技有限公司 | 火电厂系统中的数据确定方法、装置、设备以及存储介质 |
CN114534458B (zh) * | 2022-03-09 | 2023-04-07 | 山东爱亿普环保科技股份有限公司 | 一种喷射脱硝还原剂的控制方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109670625A (zh) * | 2018-08-28 | 2019-04-23 | 江苏方天电力技术有限公司 | 基于无迹卡尔曼滤波最小二乘支持向量机的NOx排放浓度预测方法 |
CN110094251A (zh) * | 2019-05-05 | 2019-08-06 | 东南大学 | 基于分时段多模型建模的scr催化剂性能劣化分析方法 |
CN112085277A (zh) * | 2020-09-10 | 2020-12-15 | 南京工业大学 | 一种基于机器学习的scr脱硝系统预测模型优化方法 |
CN112085196A (zh) * | 2020-09-10 | 2020-12-15 | 南京工业大学 | 一种基于sca算法优化bp神经网络的scr脱硝系统的喷氨调节方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102012224184A1 (de) * | 2012-12-21 | 2014-06-26 | Sms Siemag Ag | Verfahren zur Vorhersage, Steuerung und/oder Regelung von Stahlwerksprozessen |
-
2021
- 2021-01-25 CN CN202110109418.6A patent/CN113009890B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109670625A (zh) * | 2018-08-28 | 2019-04-23 | 江苏方天电力技术有限公司 | 基于无迹卡尔曼滤波最小二乘支持向量机的NOx排放浓度预测方法 |
CN110094251A (zh) * | 2019-05-05 | 2019-08-06 | 东南大学 | 基于分时段多模型建模的scr催化剂性能劣化分析方法 |
CN112085277A (zh) * | 2020-09-10 | 2020-12-15 | 南京工业大学 | 一种基于机器学习的scr脱硝系统预测模型优化方法 |
CN112085196A (zh) * | 2020-09-10 | 2020-12-15 | 南京工业大学 | 一种基于sca算法优化bp神经网络的scr脱硝系统的喷氨调节方法 |
Non-Patent Citations (2)
Title |
---|
"Gray Level Image Contrast Enhancement Using Barnacles Mating Optimizer";Ahmed Shameem 等;《IEEE Access》;20200914;第8卷;第169196-169214页 * |
"基于SLRMD与多分类SVM的变压器故障诊断方法";刘宇芳 等;《机械设计与制造工程》;20191231;第48卷(第12期);第71-75页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113009890A (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113009890B (zh) | 一种基于bmo-lssvm的火电厂scr烟气脱硝控制系统的优化控制方法 | |
CN114191953B (zh) | 一种基于卷积神经网络和XGBoost的烟气脱硫脱硝控制方法 | |
CN104636600B (zh) | 基于极限学习机的高含硫天然气净化工艺建模、优化方法 | |
CN111460726B (zh) | 一种煤泥流化床锅炉脱硝系统氨逃逸的优化方法 | |
CN110299188A (zh) | 基于ga变量选择的scr烟气脱硝系统grnn建模方法 | |
CN107168065A (zh) | 一种用于选择性催化还原脱硝装置的控制方法及系统 | |
CN110264012B (zh) | 基于经验模态分解的可再生能源功率组合预测方法及系统 | |
CN112085196A (zh) | 一种基于sca算法优化bp神经网络的scr脱硝系统的喷氨调节方法 | |
CN112488315A (zh) | 一种基于深度强化学习和遗传算法的分批调度优化方法 | |
CN111368988B (zh) | 一种利用稀疏性的深度学习训练硬件加速器 | |
CN106680428A (zh) | 一种用于脱硝控制系统的软测量方法 | |
CN114791571B (zh) | 一种基于改进cso-lstm网络的锂离子电池寿命预测方法及装置 | |
CN116688754A (zh) | 船舶烟气脱硫自动控制系统及其方法 | |
CN118039022B (zh) | 二氧化硫浓度预测方法、装置及计算机设备 | |
CN115755624A (zh) | 一种基于进化算法的燃煤锅炉多目标优化方法 | |
CN116484933A (zh) | 一种基于改进布谷鸟搜索算法的bp神经网络优化方法 | |
Altuğ et al. | Predicting tailpipe NOx emission using supervised learning algorithms | |
CN108509692A (zh) | 一种基于MiUKFNN算法的高含硫天然气脱硫工艺建模方法 | |
CN116603390B (zh) | 船舶尾气洗涤系统及其方法 | |
CN116205135A (zh) | 基于数据分解和神经网络的so2排放预测方法及系统 | |
CN114860398B (zh) | 一种智能的云平台的任务调度方法、装置及设备 | |
CN101478159B (zh) | 一种暂态稳定约束电力系统潮流优化方法 | |
CN106777468A (zh) | 高含硫天然气脱硫工艺强跟踪演化建模方法 | |
CN116596049A (zh) | 基于强化学习的ai处理器架构优化方法、装置及存储介质 | |
CN116312869A (zh) | 催化裂化再生烟气氮氧化物预测方法、装置及系统 |
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 |