王麗,李蓓,張帆*,陳金偉
(1.北京化工大學(xué)機電工程學(xué)院,北京 100029;2.上海交通大學(xué)動力機械及工程教育部重點實驗室,上海 200240)
隨著社會經(jīng)濟發(fā)展,能源消費的快速增長給能源供應(yīng)和環(huán)境保護帶來了巨大壓力。我國在中短期內(nèi)傳統(tǒng)能源仍占據(jù)主導(dǎo)地位,到2050年一次能源需求中煤炭占比將大幅下降,而天然氣在能源需求中的占比將有所提升[1]??紤]到我國能源結(jié)構(gòu)及控制碳排放[2-3],對傳統(tǒng)能源的高效低成本利用對于促進能源結(jié)構(gòu)向清潔化、多元化過渡具有重要意義。
燃料電池具有發(fā)電效率較高、對環(huán)境污染小等特點,使其成為一種具有較好前景的能源利用方式[4]。其中固體氧化物燃料電池(Solid Oxide Fuel Cell,SOFC)工作溫度較高,可以產(chǎn)生高品位排氣,并與燃氣輪機(Gas Turbine,GT)組成混合發(fā)電系統(tǒng),進一步提高發(fā)電效率,實現(xiàn)能源梯級利用。SOFCGT 系統(tǒng)具有變量多、耦合強、非線性特性顯著等特點,在偏離設(shè)計工況點運行的過程中,極易出現(xiàn)系統(tǒng)性能突然降低以及設(shè)備運行安全隱患等嚴(yán)重問題,因此需要設(shè)計一套有效的控制系統(tǒng),以保證SOFC-GT 系統(tǒng)在變工況條件下運行的安全性和運行效率。在運行過程中,影響電池性能最重要的因素是溫度。
許多學(xué)者對SOFC-GT 混合系統(tǒng)的安全高效運行進行了廣泛的研究。Pezzini 等[5]針對SOFC-GT混合動力系統(tǒng)設(shè)計了基于狀態(tài)空間模型的多變量控制器,調(diào)節(jié)控制渦輪速度和陰極氣流。Chen 等[6]利用比例-積分-微分(PID)控制器設(shè)計了5 種不同的溫度控制系統(tǒng),證明了單側(cè)陰極溫度控制策略的有效性。Emami 等[7]應(yīng)用一組滿足SOFC-GT 混合系統(tǒng)的時間延遲和增益不確定性的魯棒性PI 控制器控制SOFC-GT 混合系統(tǒng)性能。Chen 等[8]設(shè)計了一種具有陽極及陰極入口溫度控制回路,采用2 個PI控制器控制SOFC 的溫度性能。Wang 等[9]提出了多回路PI 控制器與協(xié)調(diào)保護控制相結(jié)合的控制方法,消除了SOFC-GT 的不穩(wěn)定性。Bao 等[10]提出了一種單神經(jīng)元自適應(yīng)PI 反饋控制器,在系統(tǒng)遇到不確定性和擾動時自動調(diào)整控制參數(shù)。Emami等[11]采用分散式PI 控制SOFC-GT 系統(tǒng),控制渦輪速度和燃料電池質(zhì)量流量。
在傳統(tǒng)控制方法研究的基礎(chǔ)上,現(xiàn)階段很多學(xué)者引入模糊邏輯控制、滑??刂?、模型預(yù)測控制(MPC)等先進控制理論對SOFC 混合系統(tǒng)進行控制。Vigneysh 等[12]介紹了SOFC 與光伏系統(tǒng)耦合的混合系統(tǒng),利用模糊邏輯控制器來控制系統(tǒng)電壓變化。Wu 等[13]提出了一種基于高階滑模觀測器的復(fù)合非線性控制器用于估計系統(tǒng)內(nèi)部溫度,并應(yīng)用反饋線性化控制器和前饋控制器的復(fù)合控制器控制溫度梯度。Tsai等[14]研究使用在線前饋人工神經(jīng)網(wǎng)絡(luò)來控制渦輪速度和燃料電池陰極氣流。Xiong等[15]提出了基于前饋與模型預(yù)測與GPR 流量預(yù)測模型的SOFC 流量控制原理。Ferrari 等[16]和Rossi等[17]設(shè)計了MPC 控制器,測試正常運行條件下的線性化控制器設(shè)計對性能嚴(yán)重下降影響系統(tǒng)的控制效果。Tsai等[18]提出了一種多模型自適應(yīng)預(yù)測控制算法,控制渦輪轉(zhuǎn)速和陰極氣流,在不同操作條件下 實 現(xiàn) 更 優(yōu) 控 制。Ferrai 等[19]利 用MPC 實 現(xiàn) 了SOFC-GT 系統(tǒng)溫度梯度的有效控制。胡小夫等[20]研究了SOFC 與微型燃氣輪機(MGT)聯(lián)合發(fā)電系統(tǒng)的發(fā)展趨勢。
上述研究中,大部分學(xué)者針對SOFC-GT系統(tǒng)的多個運行參數(shù)進行控制,對于溫度控制多選擇PI 控制SOFC 單側(cè)溫度,實際過程中SOFC 的陽極通道和陰極通道存在一定的溫差,在負載跟蹤過程中會產(chǎn)生很大的熱應(yīng)力。因此,采用單側(cè)溫度控制是不準(zhǔn)確的。而如果設(shè)計多回路分散式PI 控制器,則難以應(yīng)對系統(tǒng)的多變量耦合及大范圍變工況運行時的強非線性等問題。現(xiàn)階段的研究中,已證明了MPC能夠?qū)崿F(xiàn)SOFC-GT 系統(tǒng)的變工況調(diào)節(jié),滿足SOFC溫度調(diào)控等方面的需求。
本文以SOFC-GT混合動力系統(tǒng)為研究對象,系統(tǒng)引入引射器,實現(xiàn)陽極與陰極高溫排氣回流[21]。應(yīng)用增益調(diào)度MPC 控制策略,滿足運行過程中負荷跟蹤及SOFC 溫度小范圍調(diào)整的要求。MPC 能夠?qū)崿F(xiàn)對控制SOFC-GT混合動力系統(tǒng)的多變量控制,能夠顯式處理約束,適用于非線性系統(tǒng)或受控變量高度耦合系統(tǒng)。但MPC 是一種基于模型的控制策略,對于預(yù)測模型有較高的要求,當(dāng)模型失配時,預(yù)測效果下降,控制性能也隨之下降。因此,采用了MPC 與增益調(diào)度結(jié)合的策略,根據(jù)實際運行工況切換MPC 控制器。為了驗證所設(shè)計控制器的性能,將其與常規(guī)線性MPC 控制器和分散式SIMC-PI 控制器進行比較,結(jié)果表明增益調(diào)度MPC 控制以較短的調(diào)節(jié)時間和較小的超調(diào)實現(xiàn)更準(zhǔn)確的溫度控制。
圖1 為SOFC-GT 混合動力系統(tǒng)循環(huán)裝置結(jié)構(gòu)。由圖1可見,在陽極循環(huán)系統(tǒng)中,首先高壓燃料進入燃料換熱器中預(yù)熱,產(chǎn)生的高溫高壓燃料作為一次流送入陽極引射器中以回流部分低壓陽極排氣。回流的低壓陽極排氣溫度較高且混雜水蒸氣,可以對一次流燃料進行預(yù)熱,防止發(fā)生積碳現(xiàn)象。一次流燃料與低壓陽極排氣混合進入重整換熱器中反應(yīng)生成富氫燃料,送入SOFC陽極參與電化學(xué)反應(yīng)。
圖1 SOFC-GT 混合動力系統(tǒng)循環(huán)裝置[21]Fig.1 Circulation system of the SOFC-GT hybrid power system[21]
陰極循環(huán)系統(tǒng)中,空氣首先通過壓氣機和引風(fēng)機加壓,高壓氣體進入空氣換熱器中加熱形成高溫高壓氣體后進入陰極引射器中。同時,部分低壓陰極排氣回流到陰極引射器中,兩者混合送入SOFC陰極參加電化學(xué)反應(yīng)。回流的低壓陰極排氣溫度較高,可對一次流空氣進行預(yù)熱。
電化學(xué)反應(yīng)之后的陽極排氣和陰極排氣中,仍有部分燃料未參與反應(yīng),該部分燃料會被送到后燃室中與陰極的剩余排氣進行燃燒反應(yīng)。燃燒后生成的高溫氣體進入重整換熱器為加熱燃料提供熱量。透平排氣作為熱流依次進入空氣和燃料換熱器中以加熱空氣和燃料,實現(xiàn)了能量的多級利用,也提高了SOFC-GT混合動力系統(tǒng)的工作效率。
建立不同工況下雙輸入雙輸出的SOFC-GT 系統(tǒng)預(yù)測模型,輸入為后燃室燃料流量和燃氣輪機的轉(zhuǎn)速,輸出為SOFC 的陽極入口溫度和陰極入口溫度。當(dāng)燃氣輪機轉(zhuǎn)速或后燃室燃料流量發(fā)生階躍變化時,針對#1—#5 工況(分別對應(yīng)60%,70%,80%,90%和100%額定工況)下對SOFC 進行動態(tài)仿真分析,獲得對應(yīng)的SOFC 陰極和陽極入口的溫度響應(yīng)。針對各工況下的仿真數(shù)據(jù)進行縮放處理,將燃氣輪機轉(zhuǎn)速變化數(shù)據(jù)縮小至原數(shù)據(jù)的1/1 000,并將后燃室燃料量數(shù)據(jù)放大100 倍。其次,將每時刻的數(shù)據(jù)與初始時刻數(shù)據(jù)的差值通過傳遞函數(shù)進行辨識。通過辨識得到相應(yīng)工況下的傳遞函數(shù)矩陣。SOFC-GT混合動力系統(tǒng)邊界條件見表1。
針對上述5種工況的輸入輸出對設(shè)置不同的零極點數(shù)目,辨識出一階/二階模型,再對二者進行比較,選擇與原系統(tǒng)模型擬合精度更高的一種作為辨識模型。#1—#5工況下傳遞函數(shù)矩陣G(s)分別為
辨識出的傳遞函數(shù)模型與機理模型的擬合程度均在90%以上,基于傳遞函數(shù)設(shè)計的控制策略能夠較為精準(zhǔn)地實現(xiàn)對SOFC-GT 混合動力系統(tǒng)的控制。
為了使SOFC-GT 混合動力系統(tǒng)運行在高效合理狀態(tài),設(shè)計無偏模型預(yù)測控制算法對系統(tǒng)進行控制。本研究采用基于積分作用的無偏MPC 算法,實現(xiàn)SOFC-GT混合動力系統(tǒng)的功率跟蹤和溫度控制。
為不失一般性,考慮到受擾動影響的系統(tǒng),其標(biāo)稱離散狀態(tài)空間模型可以表示為
式中:x∈Rnx表示維數(shù)為nx的狀態(tài)量;u∈Rn u,表示維數(shù)為nu的輸入量;y∈Rny表示維數(shù)為ny的輸出量;A,B,C,D是狀態(tài)空間模型的系統(tǒng)矩陣,其維度根據(jù)系統(tǒng)的狀態(tài)量、輸入量和輸出量決定。
考慮系統(tǒng)受到不可測擾動d的影響,可表示為
式中:Bd,Dd是不可測擾動d(k)進入系統(tǒng)的擾動矩陣。
對于系統(tǒng)式(6)—(7)控制要求為將系統(tǒng)的輸出量y驅(qū)動至設(shè)定值yr,為此可構(gòu)造時的目標(biāo)函數(shù)
為克服系統(tǒng)受控過程中受到的不確定影響,實現(xiàn)無偏控制,需要對MPC 算法式(11)—(15)進行改進,通過偏差變量的形式在算法中引入積分信號消除不可測擾動在系統(tǒng)穩(wěn)態(tài)時產(chǎn)生的影響,
含有積分信號的擴增模型式(25)—(26)將作為預(yù)測模型,替代原預(yù)測模型,作為系統(tǒng)未來狀態(tài)量和輸出量的預(yù)測基礎(chǔ)。
基于擴增模型,可以獲得在k時刻的預(yù)測值為
式中:y(k+i|k),i= 1,…,N表示k時刻系統(tǒng)的預(yù)測輸出;x(k|k)為k時刻狀態(tài)值。相應(yīng)的標(biāo)函數(shù)改寫為
本研究中,觀測器增益設(shè)計為穩(wěn)態(tài)Kalman 濾波增益。利用系統(tǒng)歷史實測輸出值和輸入值,可以通過穩(wěn)態(tài)Kalman 濾波器對系統(tǒng)當(dāng)前時刻的狀態(tài)量進行估計,即
式中:L為穩(wěn)態(tài)Kalman濾波增益。
因此,在得到狀態(tài)量估計值x?(k)后,可以得到擴增狀態(tài)量x?(k)的估計值x??(k)
最終,在MPC 滾動優(yōu)化過程中,對于k時刻,僅采用第1 個取值Δu作為控制輸入的增量,最終得到k時刻后燃室燃料流量和燃氣輪機轉(zhuǎn)速的輸入值。
傳統(tǒng)MPC 控制器設(shè)計是基于系統(tǒng)某一穩(wěn)態(tài)工作點進行線性模型設(shè)計,在大工況范圍內(nèi)變化過程中不能保證系統(tǒng)在整個工況范圍內(nèi)均具有良好的控制效果。增益調(diào)度的設(shè)計方法充分利用傳統(tǒng)MPC 控制器在穩(wěn)態(tài)工作點附近具有良好控制效果的優(yōu)點,實現(xiàn)SOFC-GT混合動力系統(tǒng)在大工況范圍內(nèi)跟蹤功率設(shè)定。
結(jié)合增益調(diào)度的無偏預(yù)測控制算法重點在于采用多個不同工況的局部線性子模型來逼近全局非線性動態(tài),并為每個子模型設(shè)計滿足其穩(wěn)定性和魯棒性的MPC 控制器,通過切換控制器的方法實現(xiàn)在不同工況辨識線性模型,減小預(yù)測模型失配帶來的影響??刂七^程中為了提高運算效率,非當(dāng)前工況范圍內(nèi)的MPC 控制器不計算最優(yōu)控制序列,但為了實現(xiàn)控制器之間的無擾動切換,防止控制器切換時操縱變量突然變化,其他控制器繼續(xù)執(zhí)行狀態(tài)估計。增益調(diào)度模型預(yù)測控制結(jié)構(gòu)如圖2所示。
圖2 基于增益調(diào)度無偏模型預(yù)測控制策略控制結(jié)構(gòu)Fig.2 Control system structure unbiased model predictive control strategy based on gain scheduling
在SOFC-GT混合動力系統(tǒng)控制過程中,分別針對#1—#5 工況設(shè)計MPC 控制器。SOFC-GT 混合動力系統(tǒng)運行過程中總功率會實時變化,總功率的變化反映出系統(tǒng)工況的變動,根據(jù)產(chǎn)生的調(diào)度信號切換到最合適當(dāng)前工況的控制器,避免運行過程中出現(xiàn)模型失配現(xiàn)象。調(diào)度過程中,當(dāng)功率≤213.2 kW時選擇#1 控制器(針對#1 工況設(shè)計);功率>213.2 kW 且≤246.0 kW 時選擇#2 控制器(針對#2 工況設(shè)計);功率>246.0 kW 且≤278.8 kW 時選擇#3控制器(針對#3 工況設(shè)計);功率>278.8 kW 且≤311.0 kW時選擇#4 控制器(針對#4 工況設(shè)計);功率>311.0 kW 時選擇#5控制器(針對#5工況設(shè)計)。轉(zhuǎn)換值選擇相鄰工況下SOFC-GT混合系統(tǒng)總功率的中點值。
針對SOFC-GT混合動力系統(tǒng)的動態(tài)特性,采用所設(shè)計的增益調(diào)度無偏模型預(yù)測控制進行變工況運行時的溫度調(diào)節(jié)。以SOFC-GT 混合動力系統(tǒng)的機理模型為被控對象,以不同工況辨識得到的線性模型作為預(yù)測模型設(shè)計各局部子控制器,被控變量分別為燃料電池陽極和陰極的入口溫度,控制變量為后燃燒室燃料量和燃氣輪機轉(zhuǎn)速指令,調(diào)度變量為混合系統(tǒng)輸出總功率,該值可以反映系統(tǒng)運行的不同工況。局部控制器的參數(shù)選擇見表2。
表2 MPC控制器參數(shù)設(shè)置Table 2 Parameter settings of an MPC controller
為驗證所設(shè)計控制器的性能,將本研究提出的算法與額定工況下設(shè)計的單模型預(yù)測控制器和基于SIMC 準(zhǔn)則整定的分散式PI 控制器進行了對比分析。在系統(tǒng)工況持續(xù)變化情況下,選擇#5 工況下設(shè)計的單模型預(yù)測控制器進行系統(tǒng)控制。工況變化期間,模型預(yù)測控制器保持不變。單模型預(yù)測控制器的參數(shù)設(shè)置見表2。采用分散式PI控制回路控制陽極通道和陰極通道的進口溫度,通過SIMC 準(zhǔn)則進行PI控制器參數(shù)整定,具體參數(shù)見表3。
表3 PI控制器參數(shù)設(shè)置Table 3 Parameter settings of a PI controller
為實現(xiàn)系統(tǒng)負載跟蹤響應(yīng),保證SOFC-GT混合動力系統(tǒng)持續(xù)高效運行,需要對SOFC-GT混合動力系統(tǒng)輸出功率進行快速跟蹤控制。本仿真采用大范圍連續(xù)變工況運行,SOFC-GT 混合動力系統(tǒng)從#5工況連續(xù)降負荷運行至#1 工況,設(shè)置目標(biāo)功率為斜坡信號,在第1 000 s開始以0.132的斜率下降,直到第2 000 s下降至#1工況后保持功率穩(wěn)定。
圖3 為3 種控制策略中SOFC-GT,SOFC 和GT的輸出功率動態(tài)響應(yīng)。3 種控制策略下,單模型預(yù)測控制和增益調(diào)度模型預(yù)測控制的功率響應(yīng)近似,PI 控制中功率響應(yīng)變化幅度較大,在3 500 s 時,SOFC-GT 功率的PI 控制過程中出現(xiàn)大幅度振蕩現(xiàn)象,且出現(xiàn)較大的超調(diào)現(xiàn)象,穩(wěn)定時間很長。增益調(diào)度模型控制策略通過SOFC-GT 的實際輸出功率來判斷是否切換子模型,與功率的變化吻合,證明了增益調(diào)度模型預(yù)測控制的切換效果較好。由系統(tǒng)中SOFC和GT的功率響應(yīng)可以看出,GT的功率響應(yīng)時間短,SOFC 的功率響應(yīng)時間較長,這是因為燃氣輪機轉(zhuǎn)軸的機械慣性較小,而燃料電池內(nèi)部的熱慣性大,導(dǎo)致其響應(yīng)時間較長,因此整個混合系統(tǒng)的功率響應(yīng)時間較長。
圖3 3種控制策略中SOFC-GT,SOFC和GT的輸出功率動態(tài)響應(yīng)Fig.3 Dynamic response of output power of SOFC-GT,SOFC and GT in three control strategies
分析3 種控制策略的性能時,隨著功率設(shè)定值的降低,SOFC 內(nèi)部電化學(xué)反應(yīng)產(chǎn)生的熱能開始減少,SOFC 進口處溫度偏離了期望值,3 種控制系統(tǒng)開始產(chǎn)生動作。圖4表示后燃室燃料流量動態(tài)響應(yīng)變化,由圖4可見系統(tǒng)功率設(shè)定變化時,在PI控制系統(tǒng)中,后燃室燃料流量出現(xiàn)了一定程度的振蕩,恢復(fù)穩(wěn)定時間很長;在單模性預(yù)測控制系統(tǒng)中,后燃室燃料流量從初始值迅速階躍到0.002 92 kg/s,出現(xiàn)小幅度震蕩,最后穩(wěn)定;在增益調(diào)度模型預(yù)測控制系統(tǒng)中,后燃室燃料流量從初始值迅速階躍到0.006 76 kg/s,最后穩(wěn)定在0.001 10 kg/s。
圖4 后燃室燃料流量動態(tài)響應(yīng)Fig.4 Dynamic response of afterburner fuel flow
圖5 為在SOFC 的溫度控制中燃機轉(zhuǎn)速的動態(tài)響應(yīng)。當(dāng)系統(tǒng)功率設(shè)定發(fā)生變化的時,在增益調(diào)度模型預(yù)測控制系統(tǒng)當(dāng)中,轉(zhuǎn)速從初始值70 000 r/min降至53 820 r/min,繼而提升到58 271 r/min;單模型預(yù)測控制系統(tǒng)中,轉(zhuǎn)速從初始值降至52 630 r/min,然后穩(wěn)定在58 271 r/min,其變化幅度要比增益調(diào)度模型預(yù)測控制中大;在PI 控制過程中,轉(zhuǎn)速從初始值降至48 600 r/min,然后穩(wěn)定在57 000 r/min。因為在系統(tǒng)辨識過程中進行了縮放處理,還原縮放后偏差被放大,且PI控制器控制過程較慢,SIMC-PI在仿真時間設(shè)置為20 000 s 時尚未達到穩(wěn)定狀態(tài),穩(wěn)定時間過長,在仿真時間足夠長的情況下,3 種方法最終穩(wěn)定在同一點。
SOFC陽極進口溫度控制效果如圖6所示。3種控制策略中,PI 控制出現(xiàn)了大幅度振蕩,控制溫度時間較長,在12 000 s 后才能達到穩(wěn)定狀態(tài),這是因為PI 控制方法對于SOFC 進口溫度的控制是通過2個控制回路建立的,2 個回路會產(chǎn)生相互影響。由于陰極進口溫度的比例常數(shù)較大,導(dǎo)致控制過程中出現(xiàn)一定程度的振蕩。單模型預(yù)測控制在1 000 s時出現(xiàn)了模型不匹配問題,導(dǎo)致控制出現(xiàn)一定幅度的波動。由于工況變化范圍太大,單模型預(yù)測控制在10 000 s 左右時才達到穩(wěn)定狀態(tài)。增益調(diào)度模型預(yù)測控制波動幅度最小,在3 000 s 左右時候就能夠?qū)⒛P涂刂圃诜€(wěn)定狀態(tài),且控制過程中,陽極進口溫度的變化速率保持在3 K/min 以內(nèi)。
圖6 SOFC陽極進口溫度控制效果比較Fig.6 Control effect of the SOFC anode inlet temperature
SOFC 陰極進口溫度控制效果如圖7 所示。PI控制系統(tǒng)對后燃室中燃料流量的調(diào)節(jié)速度較慢,在12 000 s 左右將系統(tǒng)控制在穩(wěn)定狀態(tài),控制時間較長,并且在調(diào)節(jié)過程中出現(xiàn)了較大的超調(diào)現(xiàn)象。這是由于在陰極進口溫度控制回路中,PI 控制系統(tǒng)的比例常數(shù)較大導(dǎo)致的。增益調(diào)度模型預(yù)測控制和單模型預(yù)測控制在溫度偏離期望值時可迅速反應(yīng),并將溫度調(diào)節(jié)到期望值附近。但在單模型預(yù)測控制過程中,因為模型失配問題,陰極入口溫度有一定幅度的波動,導(dǎo)致系統(tǒng)不穩(wěn)定。增益調(diào)度模型預(yù)測控制方法可以在2 000 s 左右將陰極進口溫度控制在穩(wěn)定狀態(tài),且控制過程中溫度的變化速度保持在3 K/min以內(nèi)。
圖7 SOFC陰極進口溫度控制效果比較Fig.7 Control effect of the SOFC cathode inlet temperature
溫度是影響燃料電池工作性能的重要因素,為使SOFC 能夠安全穩(wěn)定的運行,必須把SOFC 的陽極進口溫度和陰極進口溫度控制在一定范圍內(nèi),而SOFC-GT 系統(tǒng)結(jié)構(gòu)復(fù)雜,具有強耦合、強非線性的特點,傳統(tǒng)PI 控制很難實現(xiàn)復(fù)雜系統(tǒng)中的控制目標(biāo)。本文針對SOFC-GT混合動力系統(tǒng)中SOFC進口溫度的控制,設(shè)計了一種增益調(diào)度模型預(yù)測控制方法。在SOFC-GT 混合動力系統(tǒng)大范圍變工況運行時,增益調(diào)度模型預(yù)測控制可以很好地適應(yīng)工況變化導(dǎo)致的模型失配問題,減小非線性帶來的控制性能下降問題,與雙回路PI控制以及常規(guī)線性MPC控制相比,能夠更好地跟蹤功率設(shè)定,控制SOFC 進口溫度恒定。