CN115155310B - 一种scr脱硝系统喷氨精准优化方法 - Google Patents
一种scr脱硝系统喷氨精准优化方法 Download PDFInfo
- Publication number
- CN115155310B CN115155310B CN202210828002.4A CN202210828002A CN115155310B CN 115155310 B CN115155310 B CN 115155310B CN 202210828002 A CN202210828002 A CN 202210828002A CN 115155310 B CN115155310 B CN 115155310B
- Authority
- CN
- China
- Prior art keywords
- nox
- concentration
- time
- nox concentration
- scr
- 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
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 52
- 229910021529 ammonia Inorganic materials 0.000 title claims abstract description 47
- 238000005457 optimization Methods 0.000 title claims abstract description 21
- 238000005507 spraying Methods 0.000 title claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 73
- 239000003245 coal Substances 0.000 claims abstract description 72
- 230000004044 response Effects 0.000 claims abstract description 41
- 238000002485 combustion reaction Methods 0.000 claims abstract description 36
- 238000002347 injection Methods 0.000 claims abstract description 30
- 239000007924 injection Substances 0.000 claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 23
- 230000008859 change Effects 0.000 claims description 43
- 230000008569 process Effects 0.000 claims description 23
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 18
- 239000003546 flue gas Substances 0.000 claims description 18
- 239000000779 smoke Substances 0.000 claims description 13
- 230000001105 regulatory effect Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 7
- 239000007789 gas Substances 0.000 claims description 6
- 230000001684 chronic effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000001052 transient effect Effects 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 238000010998 test method Methods 0.000 claims description 3
- 238000000738 capillary electrophoresis-mass spectrometry Methods 0.000 claims 2
- 238000010531 catalytic reduction reaction Methods 0.000 abstract description 2
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 403
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000001595 flow curve Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Classifications
-
- 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
-
- 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
- B01D2257/00—Components to be removed
- B01D2257/40—Nitrogen compounds
- B01D2257/404—Nitrogen oxides other than dinitrogen oxide
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D2258/00—Sources of waste gases
- B01D2258/02—Other waste gases
- B01D2258/0283—Flue gases
Landscapes
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Analytical Chemistry (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Treating Waste Gases (AREA)
- Regulation And Control Of Combustion (AREA)
Abstract
本发明公开了一种SCR脱硝系统喷氨精准优化方法。本发明通过试验测定NOx浓度测量延迟时间,试验测定风、煤调节时NOx响应时间,全负荷工况正交试验,获得影响NOx浓度的关键特征数据和全可调范围的运行数据。修正DCS系统NOx浓度测量时间,构建表征锅炉运行动态特性的数据变量,构建表征NOx生成浓度基准的数据变量;结合锅炉运行历史数据构建能实时预测锅炉动态工况下燃烧NOx浓度的数据结构,建立燃烧NOx浓度软测量模型,较NOx浓度测量系统提前约一个延迟时间软测量NOx浓度。通过NOx浓度偏置法对DCS控制逻辑改造,实现SCR精准喷氨。本发明解决了脱硝系统控制延迟的问题,实现了喷氨系统实时精准控制,改善了NOx浓度瞬时超标,同时延长喷氨阀门寿命。
Description
技术领域
本发明属于燃煤电站烟气脱硝领域,具体涉及一种SCR脱硝系统喷氨精准优化方法。
背景技术
燃煤电厂广泛采用选择性催化还原技术(SCR技术)来降低烟气中氮氧化物(NOx)的排放浓度。SCR脱硝系统通过喷入氨气来还原烟气中的NOx,因此喷氨量要与烟气中的NOx浓度及目标控制NOx浓度相匹配。由于脱硝反应存在一定的反应时间且SCR入口NOx浓度存在一定幅度的波动性,目前电厂主要采用“PID+前馈+反馈”的控制方式来控制喷氨量,以SCR入口处烟气NOx浓度为前馈量、SCR出口处烟气NOx浓度为反馈量。针对NOx浓度的测量,CEMS(烟气连续监测系统)基本上都是采用抽气采样分析的方法,从SCR入口烟道到CEMS分析仪通常距离较长,这会导致测量的NOx浓度会存在一段抽气时间的延迟,延迟时间达几十秒甚至数分钟,从而导致SCR测量控制过程存在大延迟的问题。在系统不能及时检测SCR入口NOx浓度变化时,常规的控制方式易出现对SCR入口NOx变化的反应不及时而造成SCR出口NOx浓度瞬间大幅波动的情况,喷氨流量曲线和NOx排放曲线出现大量的“毛刺”或“噪声”现象。电厂往往为了保证NOx排放不超标,只能将NOx排放控制目标值设置到远低于标准目标值,过高的脱硝效率和瞬间的过量喷氨给锅炉运行的安全性带来较大的不利影响。
SCR入口为高温高尘烟气环境,NOx等烟气组分很难像温度一样直接测量,所以NOx浓度测量延迟问题目前尚难以通过直接测量手段解决,通过软测量方法实时预测NOx浓度是一种较为认可的手段。锅炉燃烧NOx(理论上即是SCR入口NOx,本文将模型软测量值或实际值定义为燃烧NOx、CEMS实测值定义为SCR入口NOx)的生成过程受多个变量共同影响,这些变量间同时具有强耦合性,使得锅炉动态运行时燃烧NOx生成量的软测量尤为困难。目前的研究大多是基于稳态工况下建立燃烧NOx排放量的预测模型,很少涉及变工况的情况,而变工况时通常是燃烧NOx波动较大的时候,且变工况恰是锅炉运行常态。目前对CEMS测量延迟问题虽有文献提及,但有效的解决措施难以见到,关键难点在于如何实现锅炉动态工况下软测量的燃烧NOx曲线具有可靠的准确性、跟随性、实时性和稳定性。
发明内容
本发明的目的是针对现有技术的不足,提出一种SCR脱硝系统喷氨精准优化方法,通过更科学的数据结构建立精准实时预测动态工况下燃烧NOx浓度的软测量模型,解决CEMS测量SCR入口NOx浓度存在的时间延迟问题,实现了喷氨系统实时精准控制,对改善NOx浓度瞬时超标、延长喷氨阀门寿命及提高整体系统运行安全性、稳定性都具有重要作用。该方法考虑CEMS测量延迟时间、锅炉动态特性表征、锅炉静态及缓变状态下NOx浓度基准等,建立动态工况条件下的燃烧NOx浓度实时软测量模型,较CEMS系统提前约一个延迟时间获得SCR入口NOx浓度。燃烧NOx浓度软测量值与CEMS实测值的差值作为偏置量,加入到SCR喷氨控制系统的前馈逻辑中,实现DCS原有逻辑的最小改动以及优化系统的最简单投撤与可靠运行。
为实现上述目的,本发明提供的技术方案如下:
一种SCR脱硝系统喷氨精准优化方法,该方法包括以下步骤:
步骤1:试验测定SCR入口NOx浓度测量延迟时间;试验测定风量、煤量等调节时NOx浓度响应时间;开展影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,获得全可调范围的运行数据;
步骤2:将DCS系统的SCR入口NOx浓度测量值进行延迟时间的修正;构建表征锅炉运行动态特性的数据变量;构建表征燃烧NOx浓度基准的数据变量;结合影响NOx浓度生成的锅炉运行实时参数变量,构建实时预测锅炉动态工况下燃烧NOx浓度的软测量模型数据结构;
步骤3:根据软测量模型数据结构,采集锅炉运行历史数据,利用预测算法和优化算法建立能够在动态工况下预测燃烧NOx浓度的软测量模型;
步骤4:搭建与DCS系统实时数据交互的外置服务器,根据锅炉实时运行数据进行实时的软测量燃烧NOx浓度,即相比较SCR入口NOx浓度提前约一个延迟时间;计算燃烧NOx浓度与SCR入口NOx浓度的差值,即NOx浓度偏置值,并送入DCS系统,从而克服DCS系统不能进行复杂计算的问题;
步骤5:在DCS系统的喷氨控制逻辑中,将前馈变量SCR入口NOx浓度加上NOx浓度偏置值,实现SCR精准喷氨。
进一步,所述的步骤1中试验测定SCR入口NOx浓度延迟时间,具体实现如下:
首先调整锅炉工况使SCR入口NOx浓度发生变化的过程中,使用便携式烟气分析仪在SCR入口处的同一个NOx测点位置连续测量;其次修正便携式烟气分析仪测量系统自身测量延迟时间,得到时间段T内烟气中实际NOx浓度曲线,获取DCS系统中同一时间段T内SCR入口的NOx浓度曲线;最后利用相关性系数分析法计算两条NOx浓度曲线的相关性系数,采用迭代法将SCR入口NOx逐步前移并计算两条NOx浓度曲线的相关性系数。
进一步,所述的步骤1中试验测定SCR入口NOx浓度延迟时间,具体实现如下:
(1)假设时间t,t=0,1,2,3,······k,k为预估最大可能延迟时间,单位s;
(2)将SCR入口的NOx浓度曲线前移t,计算前移后SCR入口的NOx浓度曲线与实际NOx浓度曲线的相关性系数;
(3)绘制相关性系数-时间t曲线;
(4)相关性系数-时间t曲线的峰值点对应的时间t即为CEMS测量系统的测量的延迟时间tc。
进一步,所述的便携式烟气分析仪测量系统自身测量延迟时间,指包括采样枪和管路在内的便携式烟气分析仪测量系统的测量延迟时间;将该测量系统与NO标气瓶连接,打开标气瓶阀门直到烟气分析仪上示数随时间变化基本稳定,从打开阀门至烟气分析仪示数接近稳定的时间即为便携式烟气分析仪测量系统自身测量延迟时间。
进一步,所述的步骤1中试验测定风量、煤量调节时NOx浓度响应时间,具体实现如下:
锅炉的风量、煤量变化对应炉膛氧量、配风方式、负荷的变化,与NOx的生成有较强的正相关性,从工况改变指令下达、到给煤量变化、到炉膛燃烧时煤量风量变化、再到燃烧生成的NOx从炉膛到烟道到SCR入口测点位置的变化具有一定的延迟性,该延迟时间即为NOx浓度响应时间;同样通过相关性系数分析法,将风量调整过程中的实际NOx浓度与风量绘制相关性系数-时间t曲线;将煤量调整过程中的实际NOx浓度与煤量绘制相关性系数-时间t曲线;取相关性系数-时间t曲线上的峰值点对应的时间为风量、煤量调节时NOx浓度响应时间。
进一步,所述的步骤1中开展影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,具体实现如下:
锅炉在常规运行方式中,所有影响NOx生成的因素未必都能达到最大幅度的调整,尤其是组合式调整;因此采用正交试验法将负荷、燃烧器各风门、磨煤机分离器挡板、风煤比、燃尽风率的参数进行最大可调整范围的正交组合,从而获得最大可调范围的运行数据,为优化点寻优提供更全面的数据支撑。
进一步,所述的步骤2中DCS系统NOx浓度测量时间修正,具体实现如下:
由于CEMS测量系统存在延迟,故当前时刻DCS系统中NOx浓度测量值实际为延迟时间tc前烟道中的NOx浓度,因此将从DCS系统中采集的SCR入口NOx浓度前移tc时间,实现SCR入口NOx浓度与运行工况精准对应。
进一步,所述的步骤2中构建表征锅炉运行动态特性的数据变量,构建表征NOx生成浓度基准的数据变量,具体实现如下:
所述的步骤2中构建表征锅炉运行动态特性的数据变量,根据测定的风量、煤量调节NOx浓度的响应时间,将锅炉负荷、燃烧器各风门开度、各磨煤机分离器挡板、各磨煤机给煤量、各磨煤机一次风量、风煤比、燃尽风率的过去时刻的数据作为软测量模型输入端参数来表征变工况的幅度及速度情况,即表征影响NOx浓度变化的关键参数的动态变化过程,实现动态工况下的燃烧NOx软测量;
所述步骤2中构建表征NOx生成浓度基准的数据变量,以当前时刻的SCR入口NOx浓度作为输入变量,即1个SCR入口NOx浓度延迟时间tc前时刻的实际NOx浓度作为输入变量,从而修正包括煤种、设备或天气的非瞬时变化的工况参数对NOx浓度的固有的或动态的慢性影响;在此基准上软测量NOx浓度能够紧紧跟随实测的NOx浓度;
所述的过去时刻的数据,以NOx浓度响应时间为基准,取一个或多个n倍的NOx浓度响应时间前为过去时刻,同时为简化工作量,将锅炉负荷、磨煤机分离器挡板、磨煤机给煤量的变量的响应时间统一为煤量变化的NOx浓度响应时间;磨煤机一次风量、燃烧器风门、风煤比、燃尽风率的变量的响应时间统一为风量变化的NOx浓度响应时间。
进一步,所述的步骤2中构建实时预测锅炉动态工况下燃烧NOx生成浓度的软测量模型数据结构,具体实现如下:
将前述DCS系统NOx浓度测量时间修正后的实际NOx浓度作为软测量模型输出变量;将表征锅炉运行动态特性的数据变量、表征NOx生成浓度基准的数据变量、影响NOx生成的锅炉运行实时参数变量共同作为软测量模型的输入变量。
进一步,所述的步骤5在DCS系统的喷氨控制逻辑上,将原有的前馈变量SCR入口NOx浓度添加NOx浓度偏置值,实现DCS系统原有喷氨控制逻辑上的最小改动和优化系统的最简单投撤。
本发明的有益效果如下:
(1)延迟时间低。由于CEMS系统测量的SCR入口NOx浓度会存在一段较长时间的延迟,使得喷氨量不能实时跟踪工况变化。本发明通过试验测量得到CEMS测量系统的延迟时间,并将DCS系统SCR入口NOx浓度测量时间根据延迟时间进行修正,较CEMS系统提前约一个延迟时间获得实际NOx浓度。实际应用表明DCS前馈的NOx浓度信号的延迟时间由原来的约1分钟减少到10秒以内,投运后喷氨流量曲线和NOx排放浓度曲线较投运前显著平稳。
(2)调节范围广。目前对喷氨优化的研究大多是基于稳态工况下建立NOx浓度预测模型,很少涉及变工况的情况,而变工况时通常是NOx波动较大的时候,且变工况恰是锅炉运行常态。本发明基于试验结果建立了能体现实际锅炉动态工况的数据结构,再基于影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,结合AGC动态运行,获取相对全面的锅炉运行历史数据,最终使得模型在动态工况条件下依然能够正常调节,调节范围广。
(3)预测精度高。本发明基于前一时刻的SCR入口NOx浓度进行预测,能够修正如煤种或设备等非瞬时变化的工况参数的固有的或动态的慢性影响,在此基准上预测的NOx浓度能够紧紧跟随实际的NOx浓度,使得软测量NOx浓度无论在趋势上还是数值上都能紧紧跟随实际NOx浓度的变化。
(4)改造量小。常规思路一般将软测量值代替测量值去进行DCS逻辑改造,本发明通过取NOx浓度软测量值与CEMS实测值的差值作为偏置量,加入到SCR喷氨控制系统的前馈逻辑中,一方面防止了软测量过程中SCR入口NOx浓度出现瞬间变化,另一方面避免了对喷氨逻辑大幅修改而增加控制系统复杂度,只在原有的喷氨控制逻辑的前馈信号中加一个偏置计算,实现DCS原有逻辑的最小改动和优化系统的最简单投撤。
附图说明
图1为本发明提供的SCR脱硝系统喷氨精准优化系统与方法的流程图。
图2为本发明相关性系数分析法判断延迟时间的流程图。
图3为本发明提供的SCR脱硝系统喷氨精准优化系统与方法应用逻辑示意图。
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示,一种SCR脱硝系统喷氨精准优化方法,该方法包括如下步骤:
步骤1:试验测定SCR入口NOx浓度测量延迟时间;试验测定因风量、煤量等变量调节时NOx浓度响应时间;设计影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,获得全可调范围的运行数据。
上述试验结果为建立更合理有效的燃烧NOx浓度软测量模型的数据结构提供符合实际设备特性的数据及其关系支持。
步骤2:将DCS系统的SCR入口NOx浓度测量值进行延迟时间修正;构建表征锅炉运行动态特性的数据变量;构建表征燃烧NOx浓度基准的数据变量。
进一步结合影响NOx生成的锅炉运行实时参数变量,构建实时预测锅炉动态工况下燃烧NOx浓度的软测量模型数据结构。
步骤3:根据软测量模型数据结构,采集锅炉运行历史数据(包括正交试验期间的稳态、动态工况数据和其它部分时间段的动态工况数据),利用预测算法和优化算法建立能够在动态工况下预测燃烧NOx浓度的软测量模型。
步骤4:搭建与DCS系统实时数据交互的外置服务器,根据锅炉实时运行数据实时软测量燃烧NOx浓度(较SCR入口NOx浓度提前约一个延迟时间),计算燃烧NOx浓度与SCR入口NOx浓度差值,即NOx浓度偏置值,并送入DCS系统。
步骤5:在DCS系统的喷氨控制逻辑中,将前馈变量SCR入口NOx浓度加上NOx浓度偏置值,实现SCR精准喷氨。
进一步,所述的关键调控设备是指影响NOx浓度的关键设备。
进一步,所述步骤1中试验测定SCR入口NOx浓度延迟时间,在锅炉工况调整使SCR入口NOx浓度发生变化过程中,使用便携式烟气分析仪在SCR入口处的同一个NOx测点位置连续测量。修正便携式烟气分析仪测量系统自身测量延迟时间,得到时间段T内烟气中实际NOx浓度曲线,获取DCS系统中同一时间段T内SCR入口的NOx浓度曲线。利用相关性系数分析法(本发明优选斯皮尔曼相关性系数,亦可用其它相关性系数)计算两条NOx浓度曲线的相关性系数,采用迭代法将SCR入口NOx逐步前移并计算两条NOx浓度曲线的相关性系数;计算过程如图2所示,具体为:
(1)假设时间t,t=0,1,2,3,······k,k为预估最大可能延迟时间,单位s;
(2)将SCR入口的NOx浓度曲线前移t,计算前移后SCR入口的NOx浓度曲线与实际NOx浓度曲线的相关性系数;
(3)绘制相关性系数-时间t曲线;
(4)相关性系数-时间t曲线的峰值点对应的时间t即为CEMS测量系统的测量延迟时间tc。
所述斯皮尔曼相关性系数计算公式如下:
由于实际应用中,变量间的连结是无关紧要的,故可简化计算为:
其中,ρs为斯皮尔曼相关性系数,n为数据的数量,di为两个数据次序的差值:
di=rg(Xi)-rg(Yi) (3)
其中,rg(Xi)和rg(Yi)分别表示两列数据按照升序(降序)排列后,第i个数据在原始排列中的位置。本发明中两列数据是SCR入口NOx浓度曲线和实际NOx浓度曲线。
进一步,所述便携式烟气分析仪测量系统自身测量延迟时间,指包括采样枪和管路在内的便携式烟气分析仪测量系统的测量延迟时间。将该测量系统与NO标气瓶连接,打开标气瓶阀门直到烟气分析仪上示数随时间变化基本稳定,从打开阀门至烟气分析仪示数接近稳定的时间即为便携式烟气分析仪测量系统自身测量延迟时间。
进一步,所述步骤1中试验测定风量、煤量等变量调节时NOx浓度响应时间,锅炉的风量、煤量变化对应炉膛氧量、配风方式、负荷等的变化,与NOx的生成有较强的正相关性,从工况改变指令下达、到给煤量变化、到炉膛燃烧时煤量风量变化、再到燃烧生成的NOx从炉膛到烟道到SCR入口测点位置的变化具有一定的延迟性,该延迟时间即为NOx浓度响应时间。同样通过上述相关性系数分析法,将风量调整过程中的实际NOx浓度(SCR入口NOx浓度进行延迟时间修正后)与风量绘制相关性系数-时间t曲线;绘制煤量调整过程中的实际NOx浓度与煤量的相关性系数-时间t曲线;取相关性系数-时间t曲线上的峰值点对应的时间为风量、煤量调节时NOx浓度响应时间。
绘制风量调整过程中的实际NOx浓度与风量的相关性系数-时间t曲线具体实现如下:
(1)假设调整风量响应时间t1,t1=0,1,2,3,······k1,k1为预估最大可能响应时间,单位s;
(2)将实际NOx浓度曲线前移t1,计算前移后的实际NOx浓度曲线与风量曲线的相关性系数;
(3)绘制风量调整过程中的实际NOx浓度与风量的相关性系数-时间t曲线;
(4)相关性系数-时间t曲线的峰值点对应的时间t即为风量调节时NOx浓度响应时间。
绘制煤量调整过程中的实际NOx浓度与煤量的相关性系数-时间t曲线具体实现如下:
(1)假设调整煤量响应时间t2,t2=0,1,2,3,······k2,k2为预估最大可能响应时间,单位s;
(2)将实际NOx浓度曲线前移t2,计算前移后的实际NOx浓度曲线与煤量曲线的相关性系数;
(3)绘制煤量调整过程中的实际NOx浓度与风量的相关性系数-时间t曲线;
(4)相关性系数-时间t曲线的峰值点对应的时间t即为煤量调节时NOx浓度响应时间。
进一步,所述步骤1中设计影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,锅炉在常规运行方式中,所有影响NOx生成的因素未必都能达到最大幅度的调整,尤其是组合式调整,因此采用正交试验法将负荷、燃烧器各风门、磨煤机分离器挡板、风煤比、燃尽风率等参数进行最大可调整范围的正交组合,从而获得最大可调范围的运行数据,为优化点寻优提供更全面的数据支撑。
进一步,所述步骤2中DCS系统NOx浓度测量时间修正,由于上述CEMS测量系统存在延迟,故当前时刻DCS系统中NOx浓度测量值实际为延迟时间tc前烟道中的NOx浓度,因此将从DCS系统中采集的SCR入口NOx浓度前移tc时间,实现SCR入口NOx浓度与运行工况精准对应。
进一步,所述步骤2中构建表征锅炉运行动态特性的数据变量,根据测定的风量、煤量等调节NOx浓度的响应时间,将锅炉负荷、燃烧器各风门开度、各磨煤机分离器挡板、各磨煤机给煤量、各磨煤机一次风量、风煤比、燃尽风率等变量的过去时刻的数据作为软测量模型输入端参数来表征变工况的幅度及速度情况,即表征影响NOx浓度变化的关键参数的动态变化过程,实现动态工况下的燃烧NOx软测量。
所述的过去时刻的数据,以NOx浓度响应时间为基准,取一个或多个n倍的NOx浓度响应时间前为过去时刻,n取值为1~5。为简化工作量,可以将锅炉负荷、磨煤机分离器挡板、磨煤机给煤量等变量的响应时间统一为煤量变化的NOx浓度响应时间,磨煤机一次风量、燃烧器风门、风煤比、燃尽风率等变量的响应时间统一为风量变化的NOx浓度响应时间,亦可精细化测试每个变量变化的NOx浓度响应时间。
进一步,所述步骤2中构建表征NOx生成浓度基准的数据变量,以当前时刻的SCR入口NOx浓度(即1个SCR入口NOx浓度延迟时间tc前时刻的实际NOx浓度)作为输入变量,能够修正如煤种、设备或天气等非瞬时变化的工况参数对NOx浓度的固有的或动态的慢性影响。在此基准上软测量NOx浓度能够紧紧跟随实测的NOx浓度。
所述步骤2中构建实时预测锅炉动态工况下燃烧NOx生成浓度的软测量模型数据结构,将前述DCS系统NOx浓度测量时间修正后的实际NOx浓度作为软测量模型输出变量。将表征锅炉运行动态特性的数据变量、表征NOx生成浓度基准的数据变量、影响NOx生成的锅炉运行实时参数变量共同作为软测量模型的输入变量。
进一步,所述步骤3中锅炉运行历史数据,包括影响NOx浓度的关键调控设备及负荷全可调范围的正交试验数据及AGC状态下全负荷工况运行数据。
所述步骤3中利用预测算法及优化算法建立能够在动态工况下预测燃烧NOx浓度的软测量模型,采用神经网络、支持向量机、XGBoost等预测算法建立燃烧NOx浓度软测量模型,建模过程中采用遗传算法等优化算法对预测算法的关键参数进行寻优,建模总体过程与常规相关文献基本一致。本发明实际应用的为XGBoost+遗传算法。
所述步骤4搭建与DCS数据交互的外置服务器,外置服务器与DCS系统通过RS485进行通讯,服务器在DCS系统中实时读取软测量模型所需的锅炉运行工况参数。测量燃烧NOx浓度(较SCR入口NOx浓度提前约一个延迟时间),并计算燃烧NOx浓度与SCR入口NOx浓度差值,即NOx浓度偏置值。返回DCS系统,从而克服DCS系统不能进行复杂计算的问题。
如图3所示,所述步骤5在DCS系统的喷氨控制逻辑上,将原有的前馈变量SCR入口NOx浓度加上NOx浓度偏置值。采用偏置值的方法主要是为了一方面防止软测量过程中SCR入口NOx浓度出现瞬间变化,另一方面避免对喷氨逻辑大幅修改而增加控制系统复杂度,只在原有的喷氨控制逻辑的前馈信号中加一个偏置计算,实现DCS原有逻辑的最小改动和优化系统的最简单投撤。
以上所述仅为本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种SCR脱硝系统喷氨精准优化方法,其特征在于包括以下步骤:
步骤1:试验测定SCR入口NOx浓度测量延迟时间;试验测定风量、煤量调节时NOx浓度响应时间;开展影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,获得全可调范围的运行数据;
步骤2:将DCS系统的SCR入口NOx浓度测量值进行延迟时间的修正;构建表征锅炉运行动态特性的数据变量;构建表征燃烧NOx浓度基准的数据变量;结合影响NOx浓度生成的锅炉运行实时参数变量,构建实时预测锅炉动态工况下燃烧NOx浓度的软测量模型数据结构;
步骤3:根据软测量模型数据结构,采集锅炉运行历史数据,利用预测算法和优化算法建立能够在动态工况下预测燃烧NOx浓度的软测量模型;
步骤4:搭建与DCS系统实时数据交互的外置服务器,根据锅炉实时运行数据进行实时的软测量燃烧NOx浓度,即相比较SCR入口NOx浓度提前一个延迟时间;计算燃烧NOx浓度与SCR入口NOx浓度的差值,即NOx浓度偏置值,并送入DCS系统,从而克服DCS系统不能进行复杂计算的问题;
步骤5:在DCS系统的喷氨控制逻辑中,将前馈变量SCR入口NOx浓度加上NOx浓度偏置值,实现SCR精准喷氨;
步骤1中试验测定SCR入口NOx浓度延迟时间,具体实现如下:
首先调整锅炉工况使SCR入口NOx浓度发生变化的过程中,使用便携式烟气分析仪在SCR入口处的同一个NOx测点位置连续测量;其次修正便携式烟气分析仪测量系统自身测量延迟时间,得到时间段T内烟气中实际NOx浓度曲线,获取DCS系统中同一时间段T内SCR入口的NOx浓度曲线;最后利用相关性系数分析法计算两条NOx浓度曲线的相关性系数,采用迭代法将SCR入口NOx逐步前移并计算两条NOx浓度曲线的相关性系数;
步骤1中试验测定SCR入口NOx浓度延迟时间,具体实现如下:
(1)假设时间t,t=0,1,2,3,······k,k为预估最大可能延迟时间,单位s;
(2)将SCR入口的NOx浓度曲线前移t,计算前移后SCR入口的NOx浓度曲线与实际NOx浓度曲线的相关性系数;
(3)绘制相关性系数-时间t曲线;
(4)相关性系数-时间t曲线的峰值点对应的时间t即为CEMS测量系统的测量的延迟时间tc。
2.根据权利要求1所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于所述便携式烟气分析仪测量系统自身测量延迟时间,指包括采样枪和管路在内的便携式烟气分析仪测量系统的测量延迟时间;将该测量系统与NO标气瓶连接,打开标气瓶阀门直到烟气分析仪上示数随时间变化基本稳定,从打开阀门至烟气分析仪示数接近稳定的时间即为便携式烟气分析仪测量系统自身测量延迟时间。
3.根据权利要求1所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于所述步骤1中试验测定风量、煤量调节时NOx浓度响应时间,具体实现如下:
锅炉的风量、煤量变化对应炉膛氧量、配风方式、负荷的变化,与NOx的生成有较强的正相关性,从工况改变指令下达、到给煤量变化、到炉膛燃烧时煤量风量变化、再到燃烧生成的NOx从炉膛到烟道到SCR入口测点位置的变化具有一定的延迟性,该延迟时间即为NOx浓度响应时间;同样通过相关性系数分析法,将风量调整过程中的实际NOx浓度与风量绘制相关性系数-时间t曲线;将煤量调整过程中的实际NOx浓度与煤量绘制相关性系数-时间t曲线;取相关性系数-时间t曲线上的峰值点对应的时间为风量、煤量调节时NOx浓度响应时间。
4.根据权利要求3所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于步骤1中开展影响NOx浓度的关键调控设备及负荷全可调范围的正交试验,具体实现如下:
锅炉在常规运行方式中,所有影响NOx生成的因素未必都能达到最大幅度的调整;因此采用正交试验法将负荷、燃烧器各风门、磨煤机分离器挡板、风煤比、燃尽风率的参数进行最大可调整范围的正交组合,从而获得最大可调范围的运行数据,为优化点寻优提供更全面的数据支撑。
5.根据权利要求4所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于步骤2中DCS系统NOx浓度测量时间修正,具体实现如下:
由于CEMS测量系统存在延迟,故当前时刻DCS系统中NOx浓度测量值实际为延迟时间tc前烟道中的NOx浓度,因此将从DCS系统中采集的SCR入口NOx浓度前移tc时间,实现SCR入口NOx浓度与运行工况精准对应。
6.根据权利要求5所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于步骤2中构建表征锅炉运行动态特性的数据变量,构建表征NOx生成浓度基准的数据变量,具体实现如下:
所述的步骤2中构建表征锅炉运行动态特性的数据变量,根据测定的风量、煤量调节NOx浓度的响应时间,将锅炉负荷、燃烧器各风门开度、各磨煤机分离器挡板、各磨煤机给煤量、各磨煤机一次风量、风煤比、燃尽风率的过去时刻的数据作为软测量模型输入端参数来表征变工况的幅度及速度情况,即表征影响NOx浓度变化的关键参数的动态变化过程,实现动态工况下的燃烧NOx软测量;
所述步骤2中构建表征NOx生成浓度基准的数据变量,以当前时刻的SCR入口NOx浓度作为输入变量,即1个SCR入口NOx浓度延迟时间tc前时刻的实际NOx浓度作为输入变量,从而修正包括煤种、设备或天气的非瞬时变化的工况参数对NOx浓度的固有的或动态的慢性影响;在此基准上软测量NOx浓度能够紧紧跟随实测的NOx浓度;
所述的过去时刻的数据,以NOx浓度响应时间为基准,取一个或多个n倍的NOx浓度响应时间前为过去时刻,同时为简化工作量,将锅炉负荷、磨煤机分离器挡板、磨煤机给煤量的变量的响应时间统一为煤量变化的NOx浓度响应时间;磨煤机一次风量、燃烧器风门、风煤比、燃尽风率的变量的响应时间统一为风量变化的NOx浓度响应时间。
7.根据权利要求6所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于所述步骤2中构建实时预测锅炉动态工况下燃烧NOx生成浓度的软测量模型数据结构,具体实现如下:
将前述DCS系统NOx浓度测量时间修正后的实际NOx浓度作为软测量模型输出变量;将表征锅炉运行动态特性的数据变量、表征NOx生成浓度基准的数据变量、影响NOx生成的锅炉运行实时参数变量共同作为软测量模型的输入变量。
8.根据权利要求7所述的一种SCR脱硝系统喷氨精准优化方法,其特征在于步骤5在DCS系统的喷氨控制逻辑上,将原有的前馈变量SCR入口NOx浓度添加NOx浓度偏置值,实现DCS系统原有喷氨控制逻辑上的最小改动和优化系统的最简单投撤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210828002.4A CN115155310B (zh) | 2022-07-13 | 2022-07-13 | 一种scr脱硝系统喷氨精准优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210828002.4A CN115155310B (zh) | 2022-07-13 | 2022-07-13 | 一种scr脱硝系统喷氨精准优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115155310A CN115155310A (zh) | 2022-10-11 |
CN115155310B true CN115155310B (zh) | 2024-04-26 |
Family
ID=83494153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210828002.4A Active CN115155310B (zh) | 2022-07-13 | 2022-07-13 | 一种scr脱硝系统喷氨精准优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115155310B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08155267A (ja) * | 1994-12-01 | 1996-06-18 | Tokyo Electric Power Co Inc:The | 脱硝制御装置 |
CN102000482A (zh) * | 2010-11-15 | 2011-04-06 | 无锡科立泰科技有限公司 | 氮氧化物去除控制系统和方法 |
CN103034208A (zh) * | 2012-12-17 | 2013-04-10 | 江苏方天电力技术有限公司 | 燃煤发电机组脱硝电价监控方法 |
CN106681381A (zh) * | 2017-01-03 | 2017-05-17 | 华北电力大学 | 一种基于智能前馈信号的scr脱硝系统喷氨量优化控制系统及方法 |
CN110501901A (zh) * | 2019-07-08 | 2019-11-26 | 华能南京金陵发电有限公司 | 基于全负荷精准脱硝要求的scr喷氨调节优化的方法 |
CN112580250A (zh) * | 2020-11-12 | 2021-03-30 | 山东纳鑫电力科技有限公司 | 基于深度学习的火电机组脱硝系统及优化控制方法 |
CN113433911A (zh) * | 2021-06-30 | 2021-09-24 | 浙江大学 | 基于浓度精准预测的脱硝装置喷氨精准控制系统与方法 |
-
2022
- 2022-07-13 CN CN202210828002.4A patent/CN115155310B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08155267A (ja) * | 1994-12-01 | 1996-06-18 | Tokyo Electric Power Co Inc:The | 脱硝制御装置 |
CN102000482A (zh) * | 2010-11-15 | 2011-04-06 | 无锡科立泰科技有限公司 | 氮氧化物去除控制系统和方法 |
CN103034208A (zh) * | 2012-12-17 | 2013-04-10 | 江苏方天电力技术有限公司 | 燃煤发电机组脱硝电价监控方法 |
CN106681381A (zh) * | 2017-01-03 | 2017-05-17 | 华北电力大学 | 一种基于智能前馈信号的scr脱硝系统喷氨量优化控制系统及方法 |
CN110501901A (zh) * | 2019-07-08 | 2019-11-26 | 华能南京金陵发电有限公司 | 基于全负荷精准脱硝要求的scr喷氨调节优化的方法 |
CN112580250A (zh) * | 2020-11-12 | 2021-03-30 | 山东纳鑫电力科技有限公司 | 基于深度学习的火电机组脱硝系统及优化控制方法 |
CN113433911A (zh) * | 2021-06-30 | 2021-09-24 | 浙江大学 | 基于浓度精准预测的脱硝装置喷氨精准控制系统与方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115155310A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104826492B (zh) | 一种选择性催化还原法烟气脱硝喷氨控制系统的改进方法 | |
CN109343349B (zh) | 一种基于喷氨量补偿器的scr烟气脱硝优化控制系统和方法 | |
CN105797576B (zh) | 一种燃煤机组脱硝喷氨控制方法 | |
CN107544288B (zh) | 一种脱硝优化控制方法和系统 | |
CN104826493B (zh) | 一种选择性催化还原法烟气脱硝系统的控制方法 | |
CN108664006A (zh) | 一种基于分区控制和先进控制的scr脱硝智能喷氨优化方法及系统 | |
CN112580250A (zh) | 基于深度学习的火电机组脱硝系统及优化控制方法 | |
CN113433911B (zh) | 基于浓度精准预测的脱硝装置喷氨精准控制系统与方法 | |
CN106681381A (zh) | 一种基于智能前馈信号的scr脱硝系统喷氨量优化控制系统及方法 | |
CN108628177A (zh) | 一种基于模型自适应pid的scr脱硝智能喷氨优化方法及系统 | |
CN108837698A (zh) | 基于先进测量仪表和先进控制算法的scr脱硝喷氨优化方法及系统 | |
CN113418207B (zh) | 一种电站炉膛燃烧监测诊断装置及方法 | |
CN111897373B (zh) | 一种基于模型预测的scr脱硝装置喷氨流量调节方法 | |
CN112664975B (zh) | 一种适用于煤粉锅炉的风量控制方法 | |
CN108803309A (zh) | 一种基于软测量和模型自适应的scr脱硝智能喷氨优化方法及系统 | |
CN110501901B (zh) | 基于全负荷精准脱硝要求的scr喷氨调节优化的方法 | |
CN107831651A (zh) | 煤粉电站锅炉scr与sncr脱硝系统协同控制方法 | |
CN102460018B (zh) | 尤其矿物燃料锅炉燃烧室内燃烧过程的控制方法和燃烧系统 | |
CN115685743A (zh) | 一种智能控制燃煤锅炉及其智能预测调控烟气排放方法 | |
CN103593578A (zh) | 焦炉加热燃烧过程烟道吸力反馈设定方法 | |
CN112221347A (zh) | 一种scr脱硝系统精准喷氨控制方法 | |
CN115155310B (zh) | 一种scr脱硝系统喷氨精准优化方法 | |
CN111644054A (zh) | 一种发电厂脱硝喷氨控制系统 | |
CN115656461A (zh) | 一种基于煤质软测量的煤电机组实时碳排放量监测方法 | |
CN109260948B (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 |