馬永光, 毛求福
(華北電力大學(xué) 控制與計(jì)算機(jī)工程學(xué)院,河北 保定 071003)
隨著我國(guó)電力行業(yè)的迅猛發(fā)展,火力發(fā)電機(jī)組承擔(dān)著主要用電量供應(yīng)任務(wù),如何保證火電機(jī)組穩(wěn)定、安全、高效地運(yùn)行一直是電力行業(yè)研究的重點(diǎn)。
超臨界發(fā)電技術(shù)比較成熟,并且具有熱效率高、可靠性強(qiáng)、環(huán)保、調(diào)峰性能好等優(yōu)勢(shì),是我國(guó)目前火力發(fā)電技術(shù)主要的發(fā)展方向,并逐步成為電網(wǎng)的主力機(jī)組[1,2]。給水控制系統(tǒng)作為超臨界機(jī)組控制系統(tǒng)中非常重要的子系統(tǒng),其控制方案是以中間點(diǎn)溫度或者焓值作為反饋信號(hào),保證燃水比不變,進(jìn)而控制汽溫,以滿足不同負(fù)荷下對(duì)給水量的要求[3]。傳統(tǒng)的PID控制器結(jié)構(gòu)簡(jiǎn)單、穩(wěn)定性好、應(yīng)用成本低,在電廠中有著廣泛的用。但是該方法要求被控對(duì)象有精確的模型,而在實(shí)際運(yùn)行過(guò)程中,給水控制過(guò)程的復(fù)雜性和不確定性使得PID控制難以達(dá)到滿意的控制效果[4]。
給水控制系統(tǒng)的數(shù)學(xué)模型不僅具有較大的純滯后,而且其數(shù)學(xué)模型的參數(shù)會(huì)隨現(xiàn)場(chǎng)機(jī)組工況的改變而發(fā)生變化,因而對(duì)給水控制系統(tǒng)的設(shè)計(jì)不僅要滿足跟蹤性能好,還要達(dá)到其對(duì)控制對(duì)象數(shù)學(xué)模型變化具有較強(qiáng)的自適應(yīng)能力的要求[5]。對(duì)此,本文采用對(duì)模型要求低、跟蹤性能好、消除擾動(dòng)能力強(qiáng)的動(dòng)態(tài)矩陣控制算法[6]構(gòu)建超臨界機(jī)組給水串級(jí)控制系統(tǒng),將DMC控制器作為給水串級(jí)控制系統(tǒng)的外回路控制器,內(nèi)回路控制器仍采用傳統(tǒng)的PID控制器,改善了調(diào)節(jié)給水流量信號(hào)的大延遲特性,提高了給水串級(jí)控制系統(tǒng)的穩(wěn)定性、魯棒性和抗干擾性。
目前,超臨界機(jī)組多采用直流鍋爐,與汽包鍋爐有所不同,直流鍋爐的給水一次性經(jīng)過(guò)加熱、蒸發(fā)和過(guò)熱段,各段之間沒(méi)有固定的界線,其汽水流程示意圖如圖1所示。
圖1 超臨界機(jī)組直流鍋爐汽水流程示意圖
在正常運(yùn)行工況下,機(jī)組負(fù)荷一般不低于本生負(fù)荷(30%額定負(fù)荷),基本處于直流運(yùn)行階段,給水控制系統(tǒng)的主要任務(wù)是在不同負(fù)荷階段保證合適的燃水比,以控制過(guò)熱器出口溫度,實(shí)現(xiàn)過(guò)熱汽溫的粗調(diào)[7]。由于過(guò)熱汽溫對(duì)給水?dāng)_動(dòng)有很大的延遲,難以保證過(guò)熱汽溫的調(diào)節(jié)品質(zhì),因而一般采用汽水分離器出口蒸汽溫度作為中間點(diǎn)溫度來(lái)反映燃水比[8]。本文研究的350 MW超臨界機(jī)組給水系統(tǒng)采用中間點(diǎn)溫度作為燃水比的校正信號(hào),其簡(jiǎn)化原理結(jié)構(gòu)圖如圖2所示。
圖2 采用中間點(diǎn)溫度信號(hào)的給水控制方案
對(duì)圖2經(jīng)過(guò)分析簡(jiǎn)化可得基于中間點(diǎn)溫度的串級(jí)給水控制系統(tǒng),其中主回路為中間點(diǎn)溫度控制回路,副回路為給水流量控制回路,其簡(jiǎn)化原理圖如圖3所示。
圖3 基于中間點(diǎn)溫度的串級(jí)給水系統(tǒng)原理圖
動(dòng)態(tài)矩陣控制(DMC)算法是針對(duì)系統(tǒng)復(fù)雜化和不確定性逐漸發(fā)展起來(lái)的一種預(yù)測(cè)控制算法。該算法是基于對(duì)象階躍響應(yīng)的預(yù)測(cè)控制算法,適應(yīng)于漸進(jìn)穩(wěn)定的系統(tǒng),且具有計(jì)算量較少,魯棒性強(qiáng)等優(yōu)點(diǎn),在工業(yè)過(guò)程中應(yīng)用最為廣泛[9]。動(dòng)態(tài)矩陣控制算法包括預(yù)測(cè)模型、滾動(dòng)優(yōu)化和反饋校正3個(gè)部分[10],其算法基本結(jié)構(gòu)如圖4所示。
圖4 預(yù)測(cè)控制算法基本結(jié)構(gòu)圖
動(dòng)態(tài)矩陣控制算法采用被控對(duì)象的單位階躍響應(yīng)數(shù)據(jù)序列作為基本的預(yù)測(cè)模型。本文研究的超臨界機(jī)組給水系統(tǒng)是有自衡的,即當(dāng)系統(tǒng)穩(wěn)定以后,被控系統(tǒng)的輸出會(huì)基本保持不變。假設(shè)從時(shí)間t=0開(kāi)始采集的中間點(diǎn)焓值單位階躍響應(yīng)的數(shù)據(jù),經(jīng)過(guò)N個(gè)采樣周期,系統(tǒng)已經(jīng)進(jìn)入穩(wěn)態(tài),所獲得的N個(gè)數(shù)據(jù)可以記為{a1,a2,…,aN},如果在k時(shí)刻的控制增量為Δu(k),則未來(lái)N個(gè)時(shí)刻的輸出值為:
Δu(k)
(1)
Δu(k+j-1)
(2)
Δu2(k+j-1)
(3)
式中:qi、rj是權(quán)系數(shù),它們分別表示對(duì)跟蹤誤差以及控制量變化的抑制程度。
由于在實(shí)際中有模型失配、環(huán)境干擾等未知因素的存在,預(yù)測(cè)值有可能偏離實(shí)際值,此時(shí)需要利用實(shí)時(shí)信息進(jìn)行反饋校正,即可采用對(duì)輸出誤差e(k+1)加權(quán)的方式修正對(duì)未來(lái)輸出的預(yù)測(cè)。
(4)
根據(jù)基于中間點(diǎn)溫度的串級(jí)給水系統(tǒng)原理圖可知,在傳統(tǒng)的串級(jí)控制系統(tǒng)中,主回路和副回路均采用PID控制器。為改善系統(tǒng)的調(diào)節(jié)品質(zhì),將主回路控制器改為DMC控制器,副回路仍采用PID控制器,即形成DMC-PID串級(jí)控制器,其基本原理圖如圖5所示。
圖5 基于DMC-PID的串級(jí)給水系統(tǒng)原理圖
PID控制器作為一種常規(guī)控制器,其控制器參數(shù)的整定原則不再贅述。下面重點(diǎn)研究DMC控制器的參數(shù)選取原則[11]。根據(jù)動(dòng)態(tài)矩陣控制算法原理,需要選取的控制器參數(shù)主要有:采樣周期Ts、序列長(zhǎng)度N、預(yù)測(cè)長(zhǎng)度P、控制時(shí)域長(zhǎng)度M、加權(quán)系數(shù)矩陣H和R、預(yù)測(cè)誤差加權(quán)修正系數(shù)C。
(1)采樣周期Ts的選取
(2)序列長(zhǎng)度N和預(yù)測(cè)長(zhǎng)度P
對(duì)于預(yù)測(cè)控制算法來(lái)說(shuō),動(dòng)態(tài)矩陣的建立來(lái)源于被控對(duì)象的單位階躍響應(yīng)曲線,即動(dòng)態(tài)矩陣僅與階躍響應(yīng)序列有關(guān),所以無(wú)論N取多大,動(dòng)態(tài)矩陣只與前P個(gè)響應(yīng)數(shù)據(jù)有關(guān),如果決定了N的取值,取P=N即可。
(3)控制時(shí)域長(zhǎng)度M
M表示所要確定的未來(lái)施加的控制增量個(gè)數(shù),因而可選M≤P。但若M選取過(guò)小,難以保證系統(tǒng)的未來(lái)輸出能緊密跟蹤希望值;若M選取過(guò)大,會(huì)使系統(tǒng)的穩(wěn)定性和魯棒性都變差。
(4)加權(quán)系數(shù)矩陣H和R
誤差加權(quán)矩陣系數(shù)hi表示從(k+1)Ts時(shí)刻起到(k+P)Ts時(shí)刻止,各時(shí)刻預(yù)測(cè)偏差所占的比重,假定預(yù)測(cè)偏差所占比重相同,記為c1;而控制加權(quán)矩陣系數(shù)ri表示未來(lái)M個(gè)控制量各自所占比重,同樣假定預(yù)測(cè)控制量所占比重相同,記為c2。由于c1和c2共同決定DMC控制器輸出幅值的大小,因而都與被控對(duì)象的增益K有關(guān)。對(duì)于c1來(lái)說(shuō),K越大,c1應(yīng)該越??;對(duì)于c2來(lái)說(shuō),隨著K的增大,c2應(yīng)該有所增大。
(5)預(yù)測(cè)誤差加權(quán)修正系數(shù)C
預(yù)測(cè)誤差加權(quán)修正系數(shù)ci是根據(jù)當(dāng)前的預(yù)測(cè)誤差來(lái)進(jìn)行修正未來(lái)的預(yù)測(cè)輸出值,并且ci=ci-1+α(i=2,3,…,P),需要選取的參數(shù)為α,當(dāng)α增大時(shí),控制作用更加平緩。
蟻群優(yōu)化算法是一種通過(guò)人工模擬螞蟻搜索食物的過(guò)程,即通過(guò)每個(gè)個(gè)體之間的信息交流與相互協(xié)作最終找到從蟻穴到食物的最佳路徑的尋優(yōu)方法[12]。其基本流程圖如圖6所示。
圖6 蟻群優(yōu)化算法基本流程圖
鑒于該優(yōu)化算法具有較強(qiáng)的魯棒性、尋找路徑的并行性以及易于與其他算法相結(jié)合的優(yōu)越性,經(jīng)過(guò)近年的發(fā)展和不斷深入研究,已經(jīng)在連續(xù)域函數(shù)優(yōu)化方面得到較好的應(yīng)用。在控制系統(tǒng)對(duì)象模型辨識(shí)以及控制器參數(shù)優(yōu)化方面,很多學(xué)者也提出了不同的方法。本文針對(duì)DMC-PID串級(jí)給水控制系統(tǒng),運(yùn)用蟻群算法對(duì)給水系統(tǒng)模型進(jìn)行辨識(shí),并對(duì)DMC控制器參數(shù)和PID控制器參數(shù)進(jìn)行優(yōu)化。
為了獲得較好的過(guò)渡過(guò)程動(dòng)態(tài)特性,可以選取誤差絕對(duì)值時(shí)間積分性能指標(biāo)作為參數(shù)選擇的最小目標(biāo)函數(shù)。為了防止控制器輸出過(guò)大,在所選取的最小目標(biāo)函數(shù)中加入控制量的平方項(xiàng),即為:
(5)
式中:e(t)為系統(tǒng)誤差;u(t)為控制器輸出;w1、w2為權(quán)值。
本文以某電廠350 MW超臨界機(jī)組為研究對(duì)象,選用60%負(fù)荷工況下的現(xiàn)場(chǎng)運(yùn)行數(shù)據(jù),采用蟻群算法進(jìn)行系統(tǒng)辨識(shí),建立給水控制系統(tǒng)的傳遞函數(shù)模型。具體的操作步驟如下:
(1)數(shù)據(jù)預(yù)處理
由于采取的現(xiàn)場(chǎng)數(shù)據(jù)都有一定的噪聲,為了減少噪聲對(duì)辨識(shí)模型精度的影響,一般對(duì)數(shù)據(jù)進(jìn)行去零值、去粗大值以及數(shù)據(jù)平滑等處理,處理結(jié)果如圖7和圖8所示。
圖7 現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)曲線
圖8 現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)曲線
(2)辨識(shí)對(duì)象模型
圖9 給水流量辨識(shí)結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比曲線
圖10 中間點(diǎn)溫度辨識(shí)結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比曲線
辨識(shí)出的傳遞函數(shù)模型為:
在本文中主要通過(guò)MATLAB編程驗(yàn)證蟻群算法優(yōu)化的DMC-PID串級(jí)控制器與PID控制器的調(diào)節(jié)性能。下面分別對(duì)兩種控制算法的仿真試驗(yàn)過(guò)程進(jìn)行討論。
(1)給水系統(tǒng)的單位階躍擾動(dòng)試驗(yàn)
首先,運(yùn)用蟻群算法優(yōu)化出常規(guī)PID串級(jí)控制器參數(shù)為:δ1=-1.95、Ti1=9.37、Td1=2.34、δ2=0.64、Ti2=2.8;再對(duì)DMC控制器參數(shù)優(yōu)化為:N=29、c1=0.5、c2=0.03、α=2.3。給水系統(tǒng)做單位階躍擾動(dòng)試驗(yàn),仿真時(shí)間為2 000 s,仿真結(jié)果如圖11所示。
圖11 給水系統(tǒng)單位階躍響應(yīng)曲線
由仿真結(jié)果可知,在給水系統(tǒng)單位階躍響應(yīng)下,兩種控制方案的動(dòng)態(tài)性能如表1所示。
表1 單位階躍擾動(dòng)時(shí)動(dòng)態(tài)性能對(duì)比
由表1可知,在系統(tǒng)單位階躍擾動(dòng)下,采用常規(guī)PID控制方法時(shí)的超調(diào)量和調(diào)節(jié)時(shí)間明顯比采用DMC-PID控制方法時(shí)的要大得多,其他調(diào)節(jié)品質(zhì)幾乎差不多。綜合考慮,DMC-PID控制效果明顯優(yōu)于常規(guī)PID控制。
(2)給水系統(tǒng)的魯棒性試驗(yàn)
在實(shí)際過(guò)程中,由于建模誤差無(wú)法避免,機(jī)組工況時(shí)常發(fā)生改變等因素的存在,可能會(huì)造成內(nèi)部擾動(dòng),也會(huì)影響系統(tǒng)的控制性能。在保持DMC控制器與PID控制器參數(shù)不變的情況下,通過(guò)單位階躍擾動(dòng)試驗(yàn)檢驗(yàn)系統(tǒng)的魯棒性[13]。假設(shè)系統(tǒng)發(fā)生模型失配,被控對(duì)象的參數(shù)變化20%時(shí),系統(tǒng)做單位階躍擾動(dòng)試驗(yàn),仿真時(shí)間為2 000 s,仿真結(jié)果如圖12所示。
圖12 模型失配時(shí)系統(tǒng)單位階躍響應(yīng)曲線
由仿真結(jié)果可知,在給水系統(tǒng)模型失配時(shí)的單位階躍響應(yīng)下,兩種控制方案的動(dòng)態(tài)性能品質(zhì)如表2所示。
表2 模型失配時(shí)動(dòng)態(tài)性能品質(zhì)對(duì)比
由表2可知,由于給水系統(tǒng)復(fù)雜多變,系統(tǒng)模型會(huì)隨負(fù)荷發(fā)生變化,即當(dāng)模型失配時(shí),在系統(tǒng)單位階躍擾動(dòng)下,當(dāng)采用常規(guī)PID控制方法時(shí),超調(diào)量明顯增加,振蕩加劇,調(diào)節(jié)時(shí)間變長(zhǎng);而采用DMC-PID控制方法時(shí)超調(diào)量和調(diào)節(jié)時(shí)間都有所增加,相比較而言,調(diào)節(jié)品質(zhì)明顯比常規(guī)PID控制時(shí)要好,即采用DMC-PID控制的魯棒性明顯優(yōu)于常規(guī)PID控制。
(3)給水系統(tǒng)的抗干擾性試驗(yàn)
對(duì)給水系統(tǒng)進(jìn)行抗干擾性試驗(yàn),在給水系統(tǒng)單位階躍響應(yīng)試驗(yàn)中,當(dāng)仿真運(yùn)行到1 000 s時(shí),在控制器與被控對(duì)象之間加入階躍擾動(dòng)信號(hào),同時(shí)系統(tǒng)做單位階躍擾動(dòng),仿真結(jié)果如圖13所示。
圖13 給水系統(tǒng)的抗干擾階躍響應(yīng)曲線
在控制器與被控對(duì)象之間加入階躍擾動(dòng)信號(hào)時(shí),兩種控制方案的動(dòng)態(tài)性能如表3所示。
由表3可知,在控制器與被控對(duì)象之間加入階躍擾動(dòng)信號(hào)后,由于系統(tǒng)存在大慣性、大遲延的特性,因而系統(tǒng)的超調(diào)量和調(diào)節(jié)時(shí)間會(huì)增大,但同采用常規(guī)PID控制方法相比較,采用DMC-PID控制方法的擾動(dòng)輸出信號(hào)消除時(shí)間和超調(diào)量均要小一些,從總體上來(lái)說(shuō),DMC-PID控制的抗干擾能力要優(yōu)于常規(guī)PID控制。
表3 加入干擾信號(hào)時(shí)動(dòng)態(tài)性能品質(zhì)對(duì)比
通過(guò)以上仿真試驗(yàn)可知,本文針對(duì)給水控制系統(tǒng)數(shù)學(xué)模型所設(shè)計(jì)的DMC-PID控制方案是可行的,能夠滿足控制要求,并且在穩(wěn)定性、魯棒性和抗干擾性方面要優(yōu)于PID控制方法。
本文基于動(dòng)態(tài)矩陣控制理論,針對(duì)超臨界機(jī)組給水系統(tǒng)具有大慣性、大遲延特性且對(duì)象參數(shù)隨工況變化較大等因素,設(shè)計(jì)了DMC-PID串級(jí)控制系統(tǒng),并將蟻群算法應(yīng)用于控制器參數(shù)的優(yōu)化。同時(shí),運(yùn)用350 MW機(jī)組現(xiàn)場(chǎng)運(yùn)行數(shù)據(jù)建立了給水系統(tǒng)模型,并進(jìn)行了魯棒性和抗擾動(dòng)性檢測(cè)仿真試驗(yàn)。結(jié)果表明,DMC-PID串級(jí)控制方法使系統(tǒng)超調(diào)量變小、調(diào)節(jié)時(shí)間縮短、穩(wěn)定性提高,達(dá)到了很好的控制效果,且較常規(guī)PID串級(jí)控制方法有較強(qiáng)的抗干擾能力和更好的魯棒性,為進(jìn)一步應(yīng)用于實(shí)際電廠中提供了一定的理論基礎(chǔ)。