朱竹軍，白建云，劉林仙

(1.山西平朔煤矸石發(fā)電有限公司，山西 朔州 036800；2.山西大學(xué)自動(dòng)化系，山西 太原 030013)

0 引言

國(guó)內(nèi)外大多學(xué)者對(duì)大型火力發(fā)電機(jī)組選擇性非催化還原(selective non-catalytic reduction，SNCR)脫硝系統(tǒng)的相關(guān)研究，主要集中在NOX控制機(jī)理和工藝[1-5]。AfsinGungor研究了循環(huán)流化床中的二氧化硫和氮氧化物排放量之間的關(guān)系[1]，Javed.M等研究了SNCR技術(shù)對(duì)燃燒產(chǎn)生的氮氧化物的控制[2]。還有少部分國(guó)內(nèi)學(xué)者對(duì)火力發(fā)電機(jī)組SNCR脫硝系統(tǒng)控制的研究集中在控制方法的設(shè)計(jì)和仿真[6-9]，而國(guó)外公開文獻(xiàn)中未曾發(fā)現(xiàn)此方面的研究。國(guó)內(nèi)顧志國(guó)等設(shè)計(jì)了一種帶死區(qū)的無模型自適應(yīng)控制策略，但該控制策略只進(jìn)行了仿真，并沒有應(yīng)用于生產(chǎn)現(xiàn)場(chǎng)[6]。大多數(shù)的仿真在設(shè)計(jì)時(shí)并沒有考慮到實(shí)際運(yùn)行中的各種復(fù)雜多變的因素，以及現(xiàn)有的分散控制系統(tǒng)是否支持高級(jí)的算法等問題，從而導(dǎo)致設(shè)計(jì)的控制方法并不適用于現(xiàn)場(chǎng)實(shí)施。SNCR脫硝系統(tǒng)NOX脫除過程的大遲延、大慣性、多因素耦合等特性使常規(guī)控制方法無法滿足其控制要求，而一些先進(jìn)控制算法因不能建立比較準(zhǔn)確的數(shù)學(xué)模型等問題，很難在現(xiàn)場(chǎng)實(shí)施，加大了實(shí)現(xiàn)SNCR脫硝系統(tǒng)自動(dòng)控制的難度。

針對(duì)SNCR脫硝系統(tǒng)被控對(duì)象的特點(diǎn)，提出了一種基于專家模糊的SNCR脫硝系統(tǒng)自動(dòng)控制策略。該控制策略與DCS相結(jié)合，成功應(yīng)用于某電廠300 MW循環(huán)流化床鍋爐機(jī)組的SNCR脫硝系統(tǒng)，實(shí)現(xiàn)了機(jī)組SNCR脫硝系統(tǒng)的自動(dòng)、穩(wěn)定運(yùn)行。

1 SNCR脫硝工藝簡(jiǎn)介

某300 MW循環(huán)流化床鍋爐采用SNCR脫硝技術(shù)實(shí)現(xiàn)NOX的脫除，具體工藝如圖1所示。由圖1可知，為保證出口煙氣分析系統(tǒng)(continuous emission monitoring system，CEMS)測(cè)量的NOX值達(dá)到環(huán)保標(biāo)準(zhǔn)，需要在旋風(fēng)分離器附近噴入適量的氨水來降低NOX含量。氨水噴入量由變頻泵A通過調(diào)節(jié)頻率來控制。為保證噴頭有較好的霧化、穿透效果，通過變頻泵B使混合器出口壓力達(dá)到設(shè)定值，使之有足夠的噴射壓力，從而保證噴射效果。

圖1 SNCR脫硝工藝簡(jiǎn)圖

2 自動(dòng)控制策略設(shè)計(jì)

2.1 基于專家經(jīng)驗(yàn)的設(shè)定值形成結(jié)構(gòu)

在還原劑量設(shè)定值形成過程中，整個(gè)控制回路沒有使用常規(guī)的PID控制，而是借鑒了PID與專家控制相結(jié)合的動(dòng)態(tài)疊加方法。在PID控制過程中，如果被控量有下降的趨勢(shì)，控制量仍在累加；而采用基于專家經(jīng)驗(yàn)的動(dòng)態(tài)疊加方法，在控制量有下降趨勢(shì)時(shí)，會(huì)進(jìn)行判定，不再繼續(xù)進(jìn)行控制量的累加。其他情況與此類似。

SNCR脫硝系統(tǒng)控制策略如圖2所示。

圖2 SNCR脫硝系統(tǒng)控制策略圖

設(shè)定值形成基本結(jié)構(gòu)如圖3所示

圖3 設(shè)定值形成基本結(jié)構(gòu)圖

y(k)=y(k-1)+x

(1)

當(dāng)滿足條件時(shí)，發(fā)送一個(gè)計(jì)算機(jī)掃描周期的脈沖信號(hào)，將計(jì)算得到的還原劑量疊加到初始值；然后經(jīng)過一定的時(shí)間(即一個(gè)判定周期)，再次判斷是否滿足還原劑量變化的條件。

2.2 數(shù)據(jù)的預(yù)處理

由于現(xiàn)場(chǎng)測(cè)量得到的數(shù)據(jù)往往帶有隨機(jī)擾動(dòng)，可以通過濾波來消除干擾對(duì)控制的影響?？紤]到控制方法中需要用到的測(cè)量信號(hào)，選取了帶有隨機(jī)擾動(dòng)的被控量NOX、負(fù)荷、風(fēng)煤比、氧量、石灰石的測(cè)量值及其微分進(jìn)行濾波處理，以消除干擾對(duì)控制的影響。實(shí)際組態(tài)中，可以通過調(diào)節(jié)微分時(shí)間或加入慣性環(huán)節(jié)實(shí)現(xiàn)濾波，同時(shí)需要設(shè)置合適的參數(shù)，將隨機(jī)擾動(dòng)導(dǎo)致的數(shù)值變化與是否確實(shí)發(fā)生變化加以區(qū)分。在數(shù)據(jù)處理時(shí)，對(duì)參數(shù)進(jìn)行調(diào)整，將NOX、負(fù)荷、風(fēng)煤比、石灰石和氧量隨機(jī)干擾導(dǎo)致的變化率限制在±0.1之內(nèi)：±0.1之內(nèi)時(shí)，認(rèn)為是隨機(jī)干擾導(dǎo)致的波動(dòng)；超過±0.1時(shí)，認(rèn)為NOX、負(fù)荷等確實(shí)發(fā)生了變化。

2.3 基于專家經(jīng)驗(yàn)的區(qū)間控制

設(shè)計(jì)一個(gè)區(qū)域范圍/死區(qū)范圍，忽略被控量NOX在區(qū)域范圍內(nèi)的變化，增強(qiáng)了魯棒性；當(dāng)被控量NOX測(cè)量值在所設(shè)定區(qū)域范圍內(nèi)時(shí)，則還原劑設(shè)定值不進(jìn)行調(diào)整。當(dāng)被控量NOX測(cè)量值不在所設(shè)定區(qū)域范圍內(nèi)時(shí)，還原劑設(shè)定值則按偏離區(qū)域的程度采取相應(yīng)的專家算法來進(jìn)行調(diào)整。

①當(dāng)被控量NOX超過所設(shè)區(qū)域范圍時(shí)，還原劑調(diào)節(jié)量按照一定的專家經(jīng)驗(yàn)算法確定，具體計(jì)算方法如下：

ΔY=K1X1+K2X2+K3X3+K4X4

(2)

式中:ΔY為控制器的疊加值;X1為負(fù)荷；X2為氧量;X3為偏差；X4為偏差變化率；K1、K2、K3、K4為權(quán)重。

根據(jù)負(fù)荷、氧量、NOX偏差值、NOX變化率的變化量，決定增減還原劑的調(diào)節(jié)量。

②當(dāng)被控量NOX超過區(qū)域范圍時(shí)，將會(huì)判斷是否應(yīng)該進(jìn)行控制量的增減。

允許增加還原劑量的條件如下。

①NOX偏差值(NOX測(cè)量值減去設(shè)定值)大于所設(shè)死區(qū)值；

②NOX變化率大于-0.1；

③處于自動(dòng)調(diào)節(jié)狀態(tài)。

允許減少還原劑量的條件如下。

①NOX偏差值(NOX測(cè)量值減去設(shè)定值)小于所設(shè)死區(qū)值；

②NOX變化率小于0.1；

③處于自動(dòng)調(diào)節(jié)狀態(tài)；

④風(fēng)煤比變化率小于0.1；

⑤負(fù)荷變化率小于0.1。

當(dāng)滿足這些增減條件時(shí)，將會(huì)發(fā)送一個(gè)掃描周期的脈沖來增減由式(1)計(jì)算得到的還原劑量。給其噴入的還原劑一定的反應(yīng)時(shí)間，經(jīng)過一段時(shí)間反應(yīng)后，再次判定是否滿足條件，以決定是否增減還原劑量。

2.4 基于專家經(jīng)驗(yàn)的穩(wěn)態(tài)優(yōu)化

穩(wěn)態(tài)優(yōu)化基本原理如圖4所示。

圖4 穩(wěn)態(tài)優(yōu)化基本原理圖

圖4中：PV為測(cè)量值；SV為設(shè)定值。

為了減少還原劑消耗量，在到達(dá)NOX設(shè)定區(qū)間的情況下，進(jìn)一步對(duì)還原劑消耗量進(jìn)行優(yōu)化。如果被控量NOX在所設(shè)定區(qū)域范圍內(nèi)且NOX設(shè)定值大于NOX測(cè)量值，對(duì)還原劑設(shè)定值進(jìn)行微調(diào)，進(jìn)一步減少還原劑消耗量，使被控量NOX測(cè)量值在所設(shè)定區(qū)域范圍內(nèi)。最終使被控量NOX設(shè)定值小于NOX測(cè)量值，停止調(diào)節(jié)。

具體調(diào)節(jié)規(guī)則如下。

①將負(fù)荷經(jīng)函數(shù)F(x)折算成所需要的還原劑量。

②當(dāng)NOX實(shí)測(cè)值在所設(shè)定區(qū)域范圍內(nèi)，并且NOX沒有上升趨勢(shì)(即NOX變化率小于0.1)時(shí)，經(jīng)過判斷，減少由式(1)計(jì)算得到的還原劑量。經(jīng)過一定的反應(yīng)時(shí)間后，再判斷是否應(yīng)該進(jìn)一步優(yōu)化。

當(dāng)達(dá)到一定條件時(shí)，優(yōu)化作用消失。具體條件如下。

①負(fù)荷、風(fēng)煤比變化率大于0.1時(shí)，優(yōu)化作用消失，待穩(wěn)定后進(jìn)行條件判定。

②還原劑量小于各工況下所設(shè)的最低值時(shí)，優(yōu)化作用消失。

2.5 基于解耦規(guī)則的前饋控制

2.5.1 前饋控制的經(jīng)驗(yàn)算法

通過研究分析，發(fā)現(xiàn)負(fù)荷、風(fēng)煤比、氧量和石灰石對(duì)NOX的生成影響較大，所以采用負(fù)荷、風(fēng)煤比、石灰石作為前饋調(diào)節(jié)；而調(diào)節(jié)量的多少，滿足什么條件可進(jìn)行調(diào)節(jié)具有一定規(guī)則。具體介紹如下。

①負(fù)荷前饋調(diào)節(jié)規(guī)則。

將負(fù)荷經(jīng)函數(shù)F(x1)折算成所需要的還原劑基準(zhǔn)量，同時(shí)將氧量經(jīng)函數(shù)F(x2)折算成修正量，對(duì)還原劑基準(zhǔn)量進(jìn)行修正，即:

ΔF=F(x1)F(x2)

(3)

式中：ΔY為控制器的疊加值；x1為負(fù)荷；x2為氧量。

②風(fēng)煤比前饋調(diào)節(jié)規(guī)則。

將負(fù)荷經(jīng)函數(shù)F(x3)折算成所需要的還原劑基準(zhǔn)量，將風(fēng)煤比變化率經(jīng)函數(shù)F(x4)折算成還原劑基準(zhǔn)量的修正值,即：

ΔY=F(x3)F(x4)

(4)

式中：ΔY為控制器的疊加值；x3為負(fù)荷；x4為風(fēng)煤比變化率。

③石灰石前饋調(diào)節(jié)規(guī)則。

將負(fù)荷經(jīng)函數(shù)F(x5)折算成所需要的還原劑基準(zhǔn)量，同時(shí)將石灰石變化率經(jīng)函數(shù)F(x6)折算成還原劑基準(zhǔn)量的修正值，即:

ΔY=F(x5)F(x6)

(5)

式中：ΔY為控制器的疊加值；x5為負(fù)荷；x6為石灰石變化率。

2.5.2 前饋控制的解耦規(guī)則

當(dāng)負(fù)荷測(cè)量值變化時(shí)，風(fēng)煤比也隨之發(fā)生變化。當(dāng)負(fù)荷增大時(shí)，如果風(fēng)煤比下降幅度不是很大，則此時(shí)負(fù)荷對(duì)NOX的影響比風(fēng)煤比大；當(dāng)負(fù)荷減小時(shí)，如果風(fēng)煤比上升幅度較大，則此時(shí)風(fēng)煤比的變化對(duì)NOX的影響比負(fù)荷大。針對(duì)彼此之間的耦合關(guān)系，在控制策略中對(duì)其進(jìn)行了解耦。

(1)負(fù)荷作為前饋。

①允許增加還原劑量的條件有：

負(fù)荷變化率大于等于0.1；

NOX測(cè)量值大于等于NOX設(shè)定值；

NOX變化率大于等于-0.1；

風(fēng)煤比變化率大于等于-0.1。

②允許減少還原劑量的條件有:

負(fù)荷變化率小于等于-0.1；

NOX實(shí)測(cè)值小于等于NOX設(shè)定值；

NOX變化率小于等于0.1；

風(fēng)煤比變化率小于等于0.1。

根據(jù)歷史數(shù)據(jù)分析，負(fù)荷的變化相對(duì)較慢，選擇每分鐘進(jìn)行一次判斷。到達(dá)判定時(shí)間同時(shí)滿足以上條件時(shí)，發(fā)送一個(gè)脈沖，按照式(2)得到疊加量，進(jìn)行還原劑量調(diào)整。

(2)風(fēng)煤比作為前饋。

風(fēng)煤比作為前饋時(shí)的控制方式與負(fù)荷作為前饋時(shí)的控制方式相似，在此不再詳述。

根據(jù)歷史數(shù)據(jù)分析，風(fēng)煤比的變化相對(duì)較快，所以選擇每30 s進(jìn)行一次判斷。到達(dá)判定時(shí)間同時(shí)滿足以上條件時(shí)，發(fā)送一個(gè)脈沖，按照式(3)得到的疊加量進(jìn)行還原劑量調(diào)整。

2.6 基于模糊的還原劑低限/高限保護(hù)

在逐步降低還原劑量的同時(shí)，給其設(shè)置一個(gè)最低流量限制。其目的是要有一定的還原劑余量，使NOX的變化不太敏感。在逐步增加還原劑量的同時(shí)，給其設(shè)置一個(gè)最高流量限制，防止噴入過量的還原劑。其目的是要保證經(jīng)濟(jì)效益以及降低氨逃逸。選取負(fù)荷和氧量信號(hào)作為模糊輸入，根據(jù)歷史數(shù)據(jù)以及運(yùn)行人員的經(jīng)驗(yàn)，分別建立了模糊規(guī)則。根據(jù)輸入的數(shù)值、模糊規(guī)則庫分別計(jì)算還原劑的低限值和還原劑的高限值。

2.7 基于專家經(jīng)驗(yàn)的還原劑低限保護(hù)自適應(yīng)控制

本文提出了基于專家經(jīng)驗(yàn)的還原劑低限保護(hù)自適應(yīng)控制方法。在建立模糊規(guī)則時(shí)，最低限值選取的是正常運(yùn)行工況下的一般值。為了提高經(jīng)濟(jì)性，當(dāng)NOX生成較少時(shí)，需要采用自適應(yīng)的方法進(jìn)一步減少低限值，從而減少還原劑量。

當(dāng)還原劑量為當(dāng)前工況下對(duì)應(yīng)還原劑的低限值，且NOX設(shè)定值大于NOX測(cè)量值5 mg/Nm3時(shí)，經(jīng)過一段時(shí)間的判斷，如果這種情況一直沒有變化，則在還原劑低限值的基礎(chǔ)上進(jìn)一步減少還原劑量。然后再進(jìn)行判斷，直到NOX測(cè)量值大于NOX設(shè)定值2 mg/Nm3時(shí)，逐步恢復(fù)到原來的還原劑低限值。

2.8 基于模糊的NOX快速保護(hù)

當(dāng)燃燒工況大幅度變化導(dǎo)致NOX急劇變化時(shí)，常規(guī)控制已不能滿足控制要求。當(dāng)被控量NOX超過設(shè)定的上限值時(shí)，將程序自動(dòng)切換至基于模糊的被控量NOX快速保護(hù)。NOX快速保護(hù)算法如圖5所示。

圖5 NOX快速保護(hù)算法示意圖

由圖5可知，取NOX測(cè)量值和NOX變化率作為模糊塊的輸入。根據(jù)歷史數(shù)據(jù)以及運(yùn)行人員的經(jīng)驗(yàn)建立模糊規(guī)則，根據(jù)輸入值的大小計(jì)算出NOX所需的還原劑基準(zhǔn)值。另取負(fù)荷和氧量信號(hào)作為模糊輸入建立模糊規(guī)則，得到還原劑基準(zhǔn)值的修正值。最后得到NOX的還原劑保護(hù)量。

3 仿真驗(yàn)證

DCS仿真曲線如圖6所示。

圖6 仿真曲線

試驗(yàn)機(jī)組使用的控制系統(tǒng)為新華DCS系統(tǒng)。為了將該控制策略應(yīng)用于試驗(yàn)機(jī)組中，在仿真驗(yàn)證該控制策略的可行性時(shí)，直接在新華DCS系統(tǒng)上搭建控制對(duì)象模型，并采用控制策略程序，這樣仿真時(shí)能夠更貼近實(shí)際發(fā)現(xiàn)問題，且更方便將自動(dòng)控制程序移植到生產(chǎn)現(xiàn)場(chǎng)的DCS中。

由圖6可知，該控制策略在仿真時(shí)滿足控制要求，NOX測(cè)量值能較好地跟隨NOX設(shè)定值的變化；當(dāng)NOX因各種原因突變時(shí)，氨水量會(huì)迅速增加，使NOX快速下降。該仿真驗(yàn)證得到的可行性結(jié)果為該控制策略在試驗(yàn)機(jī)組上的實(shí)施提供了基礎(chǔ)。

4 實(shí)際運(yùn)行效果與經(jīng)濟(jì)效果分析

4.1 實(shí)際運(yùn)行效果

通過組態(tài)調(diào)試，將該控制程序成功應(yīng)用于300 MW循環(huán)流化床機(jī)組SNCR脫硝系統(tǒng)。脫硝系統(tǒng)投入自動(dòng)運(yùn)行以后，系統(tǒng)的氨水量較好地跟隨NOX的變化。當(dāng)設(shè)定值為35 mg/Nm3時(shí)，NOX排放量在30～40 mg/Nm3之間波動(dòng)，NOX排放量與設(shè)定值之間的每小時(shí)平均誤差在±2 mg/Nm3之內(nèi)。該控制策略實(shí)現(xiàn)了SNCR脫硝系統(tǒng)的自動(dòng)控制，解決了NOX突變氨水量調(diào)節(jié)跟不上以及氨水過噴的問題。

4.2 經(jīng)濟(jì)指標(biāo)

統(tǒng)計(jì)自動(dòng)化系統(tǒng)投入前后各一個(gè)星期的氨水消耗量和發(fā)電量數(shù)據(jù)，計(jì)算每300 MW需要消耗多少氨水。相關(guān)經(jīng)濟(jì)指標(biāo)數(shù)據(jù)如表1和表2所示。

表1 系統(tǒng)投入前經(jīng)濟(jì)指標(biāo)

經(jīng)計(jì)算，在手動(dòng)調(diào)節(jié)下，每300 MW·h消耗0.91 t氨水;在自動(dòng)調(diào)節(jié)下，每300 MW·h消耗0.75 t氨水，氨水消耗量每300 MW·h比原來降低17.6%。

表2 系統(tǒng)投入后經(jīng)濟(jì)指標(biāo)

5 結(jié)束語

SNCR脫硝系統(tǒng)具有大遲延、大慣性等特性。其被控對(duì)象的復(fù)雜性使得精確模型較難獲取。這些特性導(dǎo)致常規(guī)的PID控制，以及一些先進(jìn)控制應(yīng)用受限。結(jié)合運(yùn)行人員的手動(dòng)調(diào)節(jié)規(guī)律，本文提出了一種基于專家模糊的SNCR脫硝系統(tǒng)自動(dòng)控制策略。將此控制策略組態(tài)應(yīng)用于300 MW CFB機(jī)組的新華DCS系統(tǒng)。經(jīng)過調(diào)試，該策略實(shí)現(xiàn)了SNCR脫硝系統(tǒng)的穩(wěn)定及經(jīng)濟(jì)運(yùn)行，運(yùn)行人員工作量大大降低，氨水消耗量每300 MW·h，較原先降低17.6%，電廠經(jīng)濟(jì)效益顯著提高。