多變量FOPID-MAC在機(jī)爐協(xié)調(diào)控制中的應(yīng)用

Application of Multivariable Fractional Order PID-MAC in Boiler-turbine Coordinated Control

郭　偉1,2,3　王漢杰3　夏友亮3　郁振波3

(江蘇省大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心1，江蘇 南京　210044；

江蘇省氣象能源利用與控制工程技術(shù)研究中心2，江蘇 南京　210044；南京信息工程大學(xué)信息與控制學(xué)院3，江蘇 南京　210044)

多變量FOPID-MAC在機(jī)爐協(xié)調(diào)控制中的應(yīng)用

Application of Multivariable Fractional Order PID-MAC in Boiler-turbine Coordinated Control

郭偉1,2,3王漢杰3夏友亮3郁振波3

(江蘇省大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心1，江蘇 南京210044；

江蘇省氣象能源利用與控制工程技術(shù)研究中心2，江蘇 南京210044；南京信息工程大學(xué)信息與控制學(xué)院3，江蘇 南京210044)

摘要：為了克服火電單元機(jī)組機(jī)爐協(xié)調(diào)控制中的多變量、強(qiáng)耦合、強(qiáng)干擾、大時(shí)滯等問(wèn)題，提出了一種在狀態(tài)空間方程形式下的多變量模型算法控制(MAC)與分?jǐn)?shù)階PID(FOPID) 控制相結(jié)合的新型模型預(yù)測(cè)控算法MFOPID-MAC，并利用李雅普諾夫第二方法證明該算法推導(dǎo)出的控制律可以保證系統(tǒng)閉環(huán)穩(wěn)定。將改進(jìn)的多變量FOPID-MAC算法應(yīng)用于火電單元機(jī)組的機(jī)爐協(xié)調(diào)控制系統(tǒng)，并進(jìn)行階躍響應(yīng)、抗干擾仿真實(shí)驗(yàn)與絕對(duì)誤差積分指標(biāo)分析。仿真結(jié)果充分顯示多變量FOPID-MAC控制方法具有良好的設(shè)定值跟蹤特性與靜態(tài)特性、較強(qiáng)的魯棒性與滿(mǎn)意的抗干擾性。

關(guān)鍵詞：模型算法控制(MAC)分?jǐn)?shù)階PID控制多變量模型性能分析Matlab仿真

Abstract：In order to overcome the problems in boiler-turbine coordinated control for fossil fired power unit, e.g., multiple variables, strong coupling and disturbance, and large time delay, etc., the new type of multi-variable model predictive control algorithm MFOPID-MAC which combines the multi-variable model algorithm control (MAC) and fractional-order PID (FOPID) is proposed. By adopting Lyapunov the second method, it is verified that the control law derived from the algorithm can ensure the closed loop stability of the system. Step response and anti-disturbance simulation experiments of the improved MFOPID-MAC are carried out in boiler-turbine coordinated control system of fossil fired power unit, and the integral absolute error index is analyzed. The results of simulation show that the multi-variable FOPID-MAC method has satisfactory set-point tracking and static characteristics and good robustness as well as disturbance rejection capability.

Keywords：Model algorithmic control(MAC)FOPID controlMultivariable modelPerformance analysisMatlab simulation

0引言

火電單元機(jī)組協(xié)調(diào)控制是典型的多輸入多輸出系統(tǒng)，伴隨著電網(wǎng)和單元機(jī)組容量的不斷擴(kuò)大，用戶(hù)對(duì)用電質(zhì)量要求的不斷提高，研究多變量控制策略在熱工過(guò)程中的應(yīng)用愈加重要[1-4]。目前,國(guó)內(nèi)火電機(jī)組的協(xié)調(diào)控制方法有傳統(tǒng)PID控制、動(dòng)態(tài)矩陣控制策略、模糊解耦合策略、采用對(duì)象特性的在線(xiàn)預(yù)估補(bǔ)償方法等[5-7]。

本文將狀態(tài)空間方程形式下的多變量模型算法控制與分?jǐn)?shù)階PID控制算法相結(jié)合，提出了一種新型的多變量控制策略——多變量PID模型算法控制(multi-variable fractional-order PID model algorithmic control,MFOPID-MAC),并以某300 MW燃煤直流鍋爐再熱機(jī)組為模型，在Matlab中進(jìn)行仿真實(shí)驗(yàn)。研究表明，MFOPID-MAC控制器在不同負(fù)荷條件下都能保證被控量快速平穩(wěn)地追蹤設(shè)定值，靜態(tài)性能優(yōu)異，且對(duì)外部干擾具有良好的抑制效果。

1基于狀態(tài)空間形式的多變量MAC算法

模型算法控制(model algorithmic control,MAC)是預(yù)測(cè)控制中具有代表性的典型控制算法，具有預(yù)測(cè)模型、反饋校正和滾動(dòng)優(yōu)化的基本特征[9-10]。設(shè)被控對(duì)象有m個(gè)控制輸入、n個(gè)輸出，建模時(shí)域?yàn)镹，控制時(shí)域?yàn)镸，預(yù)測(cè)時(shí)域?yàn)镻，則多變量MAC預(yù)測(cè)模型取如下的狀態(tài)空間方程形式：

(1)

式中：Xm(k)=[xm1(k),xm2(k),…,xmn(k)]T；U(k)=[u1(k),u2(k),…,un(k)]T；Ym(k)=[ym1(k),ym2(k),…,ymn(k)]T；Am、Bm、Cm分別為預(yù)測(cè)模型狀態(tài)方程的系統(tǒng)矩陣、輸入矩陣和輸出矩陣。

遞推推導(dǎo)求得預(yù)測(cè)模型輸出為：

Ym(k+P)=CmXm(k+P)=CmAmPXm(k)+

(2)

在n個(gè)輸出的對(duì)象中，若取參考軌跡為一階指數(shù)形式，則(k+i)時(shí)刻參考軌跡為：

(3)

為了減少由外部擾動(dòng)和模型誤差引起的多步預(yù)測(cè)誤差，在多變量系統(tǒng)中，采用如下反饋校正：

E(k+1)=hE(k)=Yp(k)-Ym(k)

(4)

式中:E(k+1)=[e1(k+1),e2(k+1),…,en(k+1)]T，ei(k+1)=ei(k)=ypi(k)-ymi(k),i=1,2,…,n;Yp(k)=[yp1(k),yp2(k),…,ypn(k)]T;Ym(k)=[ym1(k),ym2(k),…,ymn(k)]T;h為系統(tǒng)的誤差校正矩陣。

那么，未來(lái)P時(shí)刻預(yù)測(cè)模型被修正為:

(5)

多變量MAC目標(biāo)函數(shù)二次型性能指標(biāo)形式如下：

J(k)=[Yp(k+P)-Yr(k+P)]T×Q×

(6)

(7)

2多變量FOPID-MAC控制算法與魯棒性分析

2.1　分?jǐn)?shù)階PID控制理論

1999年,Podlubny I設(shè)計(jì)了分?jǐn)?shù)階PIλDμ控制器[11]。它比傳統(tǒng)PID控制器多了2個(gè)可調(diào)參數(shù)，使得控制器設(shè)計(jì)更加多變靈活，性能更優(yōu)越。

設(shè)λ為積分階數(shù)(λ>0) ，μ為微分階數(shù)(μ>0),Kp、Ki、Kd為控制器參數(shù)，分?jǐn)?shù)階PID控制器的表達(dá)式為:

u(t)=Kpe(t)+KiDt-λe(t)+KdDtμe(t)

(8)

通過(guò)拉氏變換得到控制器的傳遞函數(shù)為[12]：

(9)

將式(8)寫(xiě)成離散形式:

u(k)=Kpe(t)+

(10)

易得：

u(k-1)=Kpe(k-1)+

(11)

將式(10)與式(11)相減，得到分?jǐn)?shù)階PID的增量式表達(dá)式：

Δu(k)=u(k)-u(k-1)=

(12)

(13)

2.2　多變量FOPID-MAC控制算法

現(xiàn)將FOPID控制和多變量MAC控制策略結(jié)合起來(lái)，使FOPID-MAC控制器既擁有分?jǐn)?shù)階比例、積分和微分的廣義結(jié)構(gòu)特點(diǎn)，又擁有MAC的預(yù)測(cè)功能。將多變量FOPID-MAC的目標(biāo)函數(shù)取為：

(14)

U(k+i)=[u1(k+i),…,um(k+i)]T

i=0,1,2…,P-1

而(k+i)時(shí)刻誤差可表示為:

E(k+i)=Yp(k+i)-Yr(k+i)=

Ym(k+i)+Yp(k)-Ym(k)-

Yr(k+i)=GiU(k)+Di(k)

(15)

定義:

(16)

(17)

(18)

(19)

2.3　多變量FOPID-MAC控制算法的穩(wěn)定性分析

(20)

將式(20)代入式(1)，可得閉環(huán)系統(tǒng)狀態(tài)如下：

Xm(k+1)=AmXm(k)-BmKXm(k)

(21)

定理1設(shè)P、Q為對(duì)稱(chēng)正定矩陣，如果滿(mǎn)足如下黎卡提方程，則多變量FOPID-MAC控制律可以保證系統(tǒng)漸近穩(wěn)定。

(Am-BmK)TP(Am-BmK)-P=-Q

(22)

證明將閉環(huán)系統(tǒng)的李雅普諾夫函數(shù)定義為如下二次函數(shù)的形式:

(23)

其中矩陣P是對(duì)稱(chēng)正定矩陣。那么李雅普諾夫函數(shù)V[Xm(k)]的增量為:

由假設(shè)(22)可知:

因此，李雅普諾夫函數(shù)在多變量FOPID-MAC控制律中單調(diào)遞減，因此閉環(huán)系統(tǒng)漸近穩(wěn)定。

3仿真研究

3.1　火電單元機(jī)組的數(shù)學(xué)模型

火電單元機(jī)組由爐、機(jī)、電、控等若干個(gè)部分構(gòu)成，具有多變量、強(qiáng)耦合、強(qiáng)干擾、大滯后的特點(diǎn)[14]。

如圖1所示，從負(fù)荷(功率)控制看，火電單元機(jī)組(包括鍋爐、汽輪機(jī)各個(gè)子系統(tǒng)在內(nèi))可以看作是一個(gè)具有雙輸入雙輸出的多變量被控系統(tǒng)[15]。

圖1中:N為輸出功率，P為機(jī)前壓力，它們是被控量；uT為主汽門(mén)調(diào)節(jié)閥開(kāi)度，uB為鍋爐負(fù)荷(燃燒率)指令，它們是控制量，GNuT(s)、GPuT(s)、GNuB(s)、GPuB(s)分別是四個(gè)回路的傳遞函數(shù)。其中除了GPuB(s)為負(fù)(uT增加，P減小)外，其余傳遞函數(shù)均為正。

圖1　火電單元機(jī)組協(xié)調(diào)控制受控對(duì)象簡(jiǎn)化框圖

被控對(duì)象的近似數(shù)學(xué)模型為:

(24)

本節(jié)的控制系統(tǒng)以某300 MW燃煤直流鍋爐再熱機(jī)組為模型，將其在100%和70%負(fù)荷下工況近似為如下的線(xiàn)性化模型[14]。

① 100%負(fù)荷

② 70%負(fù)荷

3.2　仿真實(shí)驗(yàn)與結(jié)果分析

為了證明多變量FOPID-MAC控制方法控制系統(tǒng)的動(dòng)態(tài)特性、靜態(tài)特性以及對(duì)參數(shù)攝動(dòng)的魯棒性，在確保控制器結(jié)構(gòu)與參數(shù)都一致的情況下，分別對(duì)100%負(fù)荷與70%負(fù)荷下工況進(jìn)行階躍響應(yīng)的追蹤實(shí)驗(yàn)和抗干擾實(shí)驗(yàn)；然后再和FOPID與MAC兩種控制策略進(jìn)行對(duì)比；最后還比較了在階躍實(shí)驗(yàn)中，三種控制策略的絕對(duì)誤差積分(integral absolute error,IAE)性能指標(biāo)。

3.2.1階躍實(shí)驗(yàn)

在階躍輸入下進(jìn)行仿真，圖2為功率設(shè)定值N0作15 MW(5%)階躍變化時(shí)，利用多變量FOPID-MAC、MAC、FOPID三種控制策略的仿真結(jié)果。

圖2　N0作15 MW(5%)階躍變化時(shí)仿真曲線(xiàn)

從圖2可以看出，多變量FOPID-MAC策略擁有優(yōu)越的模型適應(yīng)性和對(duì)參數(shù)攝動(dòng)的魯棒性，無(wú)論是在100%負(fù)荷或是在70%負(fù)荷下，大約經(jīng)125 s后，功率皆可趨于穩(wěn)定。與此同時(shí)，主汽壓P的動(dòng)態(tài)偏差最大值小于0.14 MPa。

在此情況下，將多變量FOPID-MAC與FOPID和MAC兩種方法相比較。比較結(jié)果顯示FOPID方法動(dòng)態(tài)響應(yīng)速度最快，但也帶來(lái)了明顯的超調(diào)。MAC方法靜態(tài)性能良好且無(wú)明顯超調(diào)，但響應(yīng)速度較慢。FOPID-MAC方法雖然在70%負(fù)荷條件下略有超調(diào)，但總體來(lái)說(shuō)，它結(jié)合了其他兩種方法的優(yōu)點(diǎn)，具有良好的動(dòng)態(tài)與靜態(tài)性能。

圖3為當(dāng)主汽壓設(shè)定值P0發(fā)生1 MPa的階躍變化時(shí)系統(tǒng)的響應(yīng)曲線(xiàn)。多變量FOPID-MAC控制策略下，主汽壓P經(jīng)過(guò)120 s就趨于穩(wěn)定。但是功率N的最大動(dòng)態(tài)偏差較大，為2.5 MW左右。此時(shí)，也將多變量FOPID-MAC與FOPID與MAC兩種方法進(jìn)行比較。

由圖3可得，多變量FOPID-MAC方法不會(huì)出現(xiàn)如FOPID方法那樣明顯的最大動(dòng)態(tài)偏差與超調(diào)量,且其響應(yīng)速度較MAC方法快，靜態(tài)性能也很好，穩(wěn)態(tài)誤差趨近于0。

圖3　P0作1 MPa階躍變化時(shí)仿真曲線(xiàn)

通過(guò)上述仿真曲線(xiàn)的對(duì)比分析可知，多變量FOPID-MAC控制策略，其被控量能快速跟蹤設(shè)定值，在70%負(fù)荷條件下，功率有略微超調(diào)卻沒(méi)有振蕩。當(dāng)其中一個(gè)被控量的設(shè)定值階躍時(shí)，另一個(gè)被控量的變化不劇烈，這表明多變量FOPID-MAC控制方法有效地削弱了多變量之間的耦合作用，保證了系統(tǒng)控制品質(zhì)。在100%負(fù)荷與70%負(fù)荷下，系統(tǒng)的動(dòng)態(tài)指標(biāo)差距不大，這說(shuō)明利用多變量FOPID-MAC方法所設(shè)計(jì)的機(jī)爐協(xié)調(diào)控制系統(tǒng)，對(duì)受控對(duì)象模型參數(shù)的變動(dòng)有很好的適應(yīng)性和魯棒性。

3.2.2抗干擾實(shí)驗(yàn)

現(xiàn)將來(lái)自于電網(wǎng)的干擾表示為作用于功率N上的附加方波信號(hào)ΔN，幅值為5 MW，在500 s時(shí)出現(xiàn)。在ΔN擾動(dòng)下，系統(tǒng)過(guò)渡過(guò)程曲線(xiàn)如圖4(100%負(fù)荷)、圖5(70%負(fù)荷)所示。

圖4　功率受擾時(shí)系統(tǒng)過(guò)渡過(guò)程曲線(xiàn)(100%負(fù)荷)

從仿真結(jié)果可以看出，多變量FOPID-MAC方法能快速對(duì)擾動(dòng)作出響應(yīng)，并迅速恢復(fù)穩(wěn)態(tài)值。與MAC方法相比，多變量FOPID-MAC由于引入了廣義的分?jǐn)?shù)階PID參數(shù)，對(duì)于耦合量P雖然帶來(lái)了小幅振蕩，但更重要的是加快了對(duì)擾動(dòng)的響應(yīng)速度。在不同負(fù)荷下，MFOPID-MAC抗擾性能相差較小，這表明控制系統(tǒng)對(duì)干擾和對(duì)像模型不確定性具有良好的適應(yīng)性。與FOPID方法相比，多變量FOPID-MAC避免了大幅的振蕩與超調(diào)，且能快速恢復(fù)到新的穩(wěn)定值。

由上述比較可得，多變量FOPID-MAC控制策略下的抗干擾性能優(yōu)于未經(jīng)優(yōu)化的FOPID與MAC控制策略，而且能迅速恢復(fù)穩(wěn)態(tài)值；在不同負(fù)荷條件下，多變量FOPID-MAC的抗干擾性能相差不大，這顯示出該控制策略對(duì)外擾和對(duì)象模型的不確定性有良好的適應(yīng)性和魯棒性。

圖5　功率受擾時(shí)系統(tǒng)過(guò)渡過(guò)程曲線(xiàn)(70%負(fù)荷)

3.2.3IAE指標(biāo)分析

在輸出功率N0作15 MW階躍響應(yīng)與機(jī)前主汽壓P0做1 MPa的響應(yīng)時(shí)，該控制系統(tǒng)在100%負(fù)荷時(shí)每一條通道的絕對(duì)誤差積分(IAE)指標(biāo)的計(jì)算結(jié)果如表1所示。

表1　標(biāo)稱(chēng)情況下各回路IAE指標(biāo)

由表1可得，MAC控制方法在耦合通道中IAE指標(biāo)較小，而在其他情況下，多變量FOPID-MAC控制策略下的IAE指標(biāo)要優(yōu)于其他控制策略。

4結(jié)束語(yǔ)

本文在多變量MAC與分?jǐn)?shù)階PID控制算法的基礎(chǔ)上，提出了新型的多變量FOPID-MAC模型算法控制策略，并以某300 MW燃煤直流鍋爐再熱機(jī)組為模型，在Matlab中進(jìn)行階躍仿真實(shí)驗(yàn)、抗干擾仿真實(shí)驗(yàn)以及絕對(duì)誤差積分性能分析，并與MAC和分?jǐn)?shù)階PID控制

方法進(jìn)行對(duì)比。由比較結(jié)果可以看出，經(jīng)改進(jìn)的算法解決了火電機(jī)組耦合性強(qiáng)、時(shí)延大等對(duì)象特性，且靜態(tài)誤差小，具有較強(qiáng)的魯棒性和滿(mǎn)意的抗干擾性能。因此，多變量FOPID-MAC模型算法控制策略是解決火電機(jī)組機(jī)爐協(xié)調(diào)控制難題的一種有效方法。

參考文獻(xiàn)

[1] 錢(qián)朝明，王會(huì)，吳永存，等.1 000 MW火電機(jī)組協(xié)調(diào)控制策略?xún)?yōu)化[J].中國(guó)電力，2013(3):68-72.

[2] 樓冠男，譚文，鄭勤玲.汽包鍋爐單元機(jī)組協(xié)調(diào)系統(tǒng)的線(xiàn)性自抗擾控制[J].中國(guó)電機(jī)工程學(xué)報(bào)，2011，31(23)：94-100.

[3] 陳敏,徐翃霞.1 000 MW超超臨界機(jī)組協(xié)調(diào)優(yōu)化控制策略及應(yīng)用[J].自動(dòng)化儀表,2014(6):50-53.

[4] 姚峻，祝建飛，金峰.1 000 MW機(jī)組節(jié)能型協(xié)調(diào)控制系統(tǒng)的設(shè)計(jì)與應(yīng)用[J].中國(guó)電力，2010(6):79-84.

[5] 周磊,費(fèi)樹(shù)岷,黃家才.一種新型的預(yù)測(cè)函數(shù)PID控制算法的研究[J].自動(dòng)化儀表,2013(9):63-66.

[6] 李玉紅，王東風(fēng)，劉紅軍，等.機(jī)爐協(xié)調(diào)系統(tǒng)的動(dòng)態(tài)矩陣控制仿真研究[J].電力科學(xué)與工程，2004(4)：11-14.

[7] 劉晴,張正剛,丁維明.神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)PID控制在氣化爐中的應(yīng)用[J].自動(dòng)化儀表,2014(5):60-62.

[8] 郭偉，姚少杰.基于時(shí)域的PID動(dòng)態(tài)矩陣控制算法改進(jìn)[J].儀器儀表學(xué)報(bào)，2007，28(12)：2174-2178.

[9] 席裕庚.預(yù)測(cè)控制[M].北京：國(guó)防工業(yè)出版社，1993：10-18.

[10]Rouhani R，Mehra R K.Model algorithmic control(MAC)；basic theoretical properties[J].Automatica，1982，18(4)：401-414.

[11]Podlubny I.Fractional-order systems and PIλDμ controllers[J].IEEE Transactions on Automatic Control，1999，44(1):208-21.

[12]郭偉，倪家健，李濤，等.基于時(shí)域的分?jǐn)?shù)階PID預(yù)測(cè)函數(shù)勵(lì)磁控制器[J].儀器儀表學(xué)報(bào)，2011，32(11):2461-2467.

[13]Tan W，Marquez H J，Chen T.Multivariable robust controller design for a boiler system[J].Control Systems Technology，IEEE Transactions on，2002，10(5)：735-742.

[14]吳國(guó)垣.熱力過(guò)程TITO系統(tǒng)的PID控制研究[D].北京：清華大學(xué)，2004.

[15]胡克定，汪自勤，錢(qián)鐘韓.單元機(jī)組協(xié)調(diào)控制一種頻域設(shè)計(jì)方法[J].東南大學(xué)學(xué)報(bào)，1989，19(1)：69-77.

中圖分類(lèi)號(hào)：TH89;TK323

文獻(xiàn)標(biāo)志碼：A

DOI:10.16086/j.cnki.issn1000-0380.201508018

江蘇省高校大學(xué)生創(chuàng)新創(chuàng)業(yè)項(xiàng)目(編號(hào)：201410300034Z)。

修改稿收到日期：2014-12-11。

第一作者郭偉(1960-)，男，1987年畢業(yè)于山東工業(yè)大學(xué)工業(yè)自動(dòng)化專(zhuān)業(yè)，獲碩士學(xué)位，教授；主要從事預(yù)測(cè)控制和電力系統(tǒng)控制等方面的教學(xué)與科研工作。