閻 群,劉 波,路 輝,崔家瑞,黃若愚,王義軒,李 擎,曹 斌

(1.北京科技大學(xué) 自動化學(xué)院,北京 100083;2.貴陽鋁鎂設(shè)計研究院有限公司,貴州 貴陽 550003)

氧化鋁濃度是鋁電解生產(chǎn)過程中至關(guān)重要的參數(shù)[1]。鋁電解生產(chǎn)時電解質(zhì)中氧化鋁濃度一般需要控制在1.5%～3%[2],過高容易在槽底生成氧化鋁沉淀,降低鋁電解的電流效率并且影響電解槽壽命;而過低則容易導(dǎo)致鋁電解過程陽極效應(yīng)頻發(fā),造成大量不必要的能源損耗和環(huán)境污染。所以,實現(xiàn)鋁電解生產(chǎn)過程中氧化鋁濃度精確控制對于鋁電解槽的物料平衡、穩(wěn)定高效運行具有十分重要的意義。

目前,國內(nèi)外針對氧化鋁濃度控制方法的研究已有許多,主要包括自適應(yīng)控制[3]、槽電阻斜率控制[4]、模糊控制[5-6],以及多種方法的綜合控制[7-8]等。文獻(xiàn)[3]和文獻(xiàn)[4]是利用槽電阻與氧化鋁濃度的關(guān)系,通過不同下料速率的切換來實現(xiàn)氧化鋁濃度的控制。這種控制思路能夠有效估計氧化鋁濃度,使電解槽處于基本穩(wěn)定的工作狀態(tài),但存在控制動作相對粗獷、氧化鋁濃度波動大等問題。文獻(xiàn)[5]提出的模糊控制方法有效利用歷史經(jīng)驗和專家知識,制定并不斷修改模糊控制規(guī)則,實現(xiàn)較高的氧化鋁濃度控制精度,但是存在模糊規(guī)則設(shè)置非常依賴專家知識和歷史經(jīng)驗的問題。文獻(xiàn)[7-8]結(jié)合多種控制方法的優(yōu)點,進(jìn)一步提高了氧化鋁濃度控制的精度,但存在控制器結(jié)構(gòu)復(fù)雜、計算所需時間長、控制成本高,難以在現(xiàn)有槽控機(jī)硬件系統(tǒng)上實現(xiàn)等問題。

為了建立能夠完整描述鋁電解過程非線性特性、便于與槽控機(jī)硬件系統(tǒng)結(jié)合,實現(xiàn)不依賴專家經(jīng)驗的氧化鋁濃度的精確控制。本團(tuán)隊在氧化鋁濃度預(yù)測[9]、氧化鋁濃度軟測量[10-11]、分布式陽極電流監(jiān)測[12-13]、數(shù)字化鋁電解槽系統(tǒng)框架設(shè)計[14]等方面做了大量研究工作。在前期工作基礎(chǔ)上,本文提出了一種基于非線性模型預(yù)測控制的氧化鋁濃度精確控制方法,其控制系統(tǒng)結(jié)構(gòu)簡圖如圖1所示。

圖1 模型預(yù)測控制系統(tǒng)結(jié)構(gòu)簡圖

1 基于機(jī)理的內(nèi)模輸入輸出變量選擇

鋁電解過程中表征鋁電解槽運行狀態(tài)的變量參數(shù)較多,如:下料量(下料間隔)、工作電壓、分布式電流、電解溫度、鋁水平、電解質(zhì)水平、分子比等。在鋁電解穩(wěn)定運行過程中,影響氧化鋁濃度變化的主要有下料過程和極距調(diào)整過程。傳統(tǒng)鋁電解槽點式下料控制方式主要采用定時下料或者欠、過量切換下料[2],即固定下料量,通過改變下料間隔來實現(xiàn)下料控制。傳統(tǒng)極距控制是基于槽電阻的區(qū)間控制[2],即根據(jù)工作電壓與目標(biāo)電壓區(qū)間的比較來判斷是否需要進(jìn)行極距調(diào)整。

針對氧化鋁濃度精確控制問題,通過綜合考慮鋁電解下料、極距調(diào)整等過程中下料量、工作電壓、陽極分布式電流、電解溫度、鋁水平、電解質(zhì)水平、分子比等關(guān)鍵性能指標(biāo)變量的可控性、可測性,以及硬件裝置的可實現(xiàn)性來選取建模變量。相關(guān)參數(shù)的特征、硬件在線檢測可實現(xiàn)性等如表1所示。

表1 鋁電解槽過程關(guān)鍵參數(shù)特征及選取情況

根據(jù)選取規(guī)則,最終選取了下料量和工作電壓作為NMPC內(nèi)模的輸入變量,氧化鋁濃度為輸出變量。

2 基于Hammerstein結(jié)構(gòu)的氧化鋁濃度子空間建模

針對氧化鋁濃度變化過程存在的非線性特性,采用Hammerstein模型結(jié)構(gòu)作為鋁電解過程氧化鋁濃度控制系統(tǒng)的內(nèi)模,如式(1)所示:

(1)

式中:t=1,…,N-1,u=[u1u2]T為內(nèi)模的輸入變量(下料量和工作電壓);

yt——模型輸出變量(氧化鋁濃度);

xt——系統(tǒng)的狀態(tài);

f(ut)=〔f1(u1(t))f2(u2(t))〕T;

vt,wt——零均值高斯白噪聲序列向量。

本文采用文獻(xiàn)[15]提出的基于最小二乘支持向量機(jī)的非線性Hammerstein系統(tǒng)子空間辨識算法。該方法在子空間辨識(N4SID)算法的基礎(chǔ)上進(jìn)行了非線性擴(kuò)展,使其能夠?qū)崿F(xiàn)對于Hammerstein模型結(jié)構(gòu)的辨識。其辨識流程如圖2所示。

圖2 建模流程

3 基于NMPC的氧化鋁濃度控制系統(tǒng)設(shè)計

針對基于Hammerstein結(jié)構(gòu)內(nèi)模的特殊性,在求解控制律的時候通過逆運算移除控制系統(tǒng)的非線性項[16],需要考慮逆運算過程中的求根問題,設(shè)計圖3所示的控制器結(jié)構(gòu)。

圖3 非線性模型預(yù)測控制系統(tǒng)結(jié)構(gòu)

圖中,P為被控對象(鋁電解槽),F和Gm串聯(lián)構(gòu)成Hammerstein結(jié)構(gòu)的內(nèi)模,Gc和F-1串聯(lián)構(gòu)成控制器,f(u)為中間變量。F為體現(xiàn)內(nèi)模非線性特性的多項式函數(shù),F-1為F的逆運算。

控制系統(tǒng)性能指標(biāo)函數(shù)定義為:

(2)

s.t.1)umin≤u≤umax

2)f(u)min≤f(u)≤f(u)max

s.t.umin≤u≤umax

3)rmin≤r≤rmax

式中:Np、Nc——分別為預(yù)測步長和控制步長;

△f(u)t+k-1——t+k-1時刻中間變量f(u)的增量(其中u=[u1u2]T);

rt+k和yt+k|t——分別為t+k時刻系統(tǒng)的設(shè)定輸出值和預(yù)測輸出值;

Q和R——分別為輸出項和控制項的權(quán)重。

針對某400 kA鋁電解槽,控制量u為下料量(u1)和工作電壓(u2),r為氧化鋁濃度設(shè)定值,需要滿足物理條件約束(u1min≤u1≤u1max,u2min≤u2≤u2max,rmin≤r≤rmax)。由于非線性項為一元高次函數(shù),存在無解和多解的情況,因此在性能指標(biāo)中引入對f(u)的約束(f(u)min≤f(u)≤f(u)max),保證求逆過程獲得具有實際物理意義的最優(yōu)實根。

控制系統(tǒng)的控制策略為:在每一個采樣時刻,將系統(tǒng)當(dāng)前時刻的狀態(tài),控制系統(tǒng)設(shè)定值軌跡和內(nèi)模代入性能指標(biāo)函數(shù)中,根據(jù)兩步法求解最佳控制律,將控制律的第一個分量作用于控制系統(tǒng)。在下一個采樣時刻,重復(fù)該過程。具體求解步驟如下:

Step1:根據(jù)式(1)的狀態(tài)空間方程進(jìn)行遞推,得到y(tǒng)(k+1|k)、y(k+2|k)…y(k+n|k)的遞推表達(dá)式(3),其中y(k+n|k)表示k時刻預(yù)測k+n時刻系統(tǒng)的輸出;

y(k+1|k)=Cx(k+1|k)+Df(u)(k+1|k)=CAx(k)+CBf(u)(k|k)+Df(u)(k+1|k)

y(k+2|k)=Cx(k+2|k)+Df(u)(k+2|k)=CA2x(k)+CABf(u)(k|k)+CBf(u)(k+1|k)+Df(u)(k+1|k)

y(k+n|k)=Cx(k+n-1|k)+Df(u)(k+n|k)=CAnx(k)+CAn-1Bf(u)(k|k)+

CAn-2Bf(u)(k+1|k)+…+CBf(u)(k+n-1|k)+Df(u)(k+n|k)

(3)

整理得:

Yk=Mx(k)+NUk

其中:

Yk=〔y(k+1|k)y(k+2|k)…y(k+n|k)〕T

M=〔CACA2CA3…CAn〕T

Uk=〔f(u)(k|k)f(u)(k+1|k)…f(u)(k+n|k)〕T

Step2:將遞推式(3)代入式(2),可得性能指標(biāo)函數(shù)的標(biāo)準(zhǔn)二次型:

(4)

式中:H=2(NTQN+R)

fT=2(Mx(k)-r)TQN

Step3:采用兩步法求解最佳控制律。首先求解step2的二次型表達(dá)式(4),得到最優(yōu)f(u);然后根據(jù)f(u)的反函數(shù)求解最佳控制律u。

4 仿真實驗與結(jié)果分析

本文對某鋁廠400 kA鋁電解槽進(jìn)行了數(shù)據(jù)采集,基于實際采集數(shù)據(jù),建立了氧化鋁濃度內(nèi)模,并應(yīng)用非線性模型預(yù)測控制實現(xiàn)對氧化鋁濃度的控制,其中控制部分考慮了無干擾情況和存在下料量不準(zhǔn)確、電壓波動等干擾的情況。

4.1 數(shù)據(jù)采集

實驗數(shù)據(jù)包括電解槽生產(chǎn)日報數(shù)據(jù)和現(xiàn)場采集數(shù)據(jù)。工作電壓、下料量通過日報數(shù)據(jù)獲取,下料間隔通過現(xiàn)場記錄獲取,氧化鋁濃度數(shù)據(jù)通過現(xiàn)場采集電解質(zhì)樣品進(jìn)行實驗室化驗獲取?？紤]到人工采集電解質(zhì)樣品的可操作性,將人工采樣間隔設(shè)置為20 s,采樣起始時刻設(shè)置為下料時刻、下料時刻+5 s、下料時刻+10 s、下料時刻+15 s。再通過插值實現(xiàn)下料周期內(nèi)電解質(zhì)樣品數(shù)據(jù)的擴(kuò)充。電解質(zhì)現(xiàn)場采集過程如圖4所示。

圖4 電解質(zhì)樣品現(xiàn)場采集圖

4.2 建模及模型誤差分析

對現(xiàn)場采集的數(shù)據(jù)進(jìn)行預(yù)處理后,令經(jīng)驗參數(shù)b=0.8,系統(tǒng)階次k=2,多項式階次m=2,建立的氧化鋁濃度模型參數(shù)為:

C=〔0.5421 -0.0337〕,D=〔0.2049 0.0426〕,

建立的內(nèi)模的均方根誤差(RMSE)為0.0930,且預(yù)測輸出的變化趨勢與實際值趨于一致,能夠較完整描述氧化鋁濃度的變化過程。

4.3 控制效果分析

某400 kA鋁電解槽單次下料時間間隔為90 s左右,由于存在下料干擾等因素,考慮30 s的冗余,設(shè)置控制系統(tǒng)的預(yù)測時長為60 s,預(yù)測時域Np=12。同時,在綜合考慮系統(tǒng)動態(tài)響應(yīng)速度、穩(wěn)態(tài)誤差和計算復(fù)雜度的基礎(chǔ)上,設(shè)置控制時域Nc=1、輸出項權(quán)重Q=1和控制項權(quán)重R=0.05。

4.3.1 無干擾情況

氧化鋁濃度初始值設(shè)定為1.5%,在100 s時刻和200 s時刻分別將氧化鋁濃度設(shè)定值改為2%和2.5%,其控制量和輸出響應(yīng)分別如圖5a～圖5d所示。

由圖5可以看出,實際的氧化鋁濃度能夠很快跟隨設(shè)定值的變化,動態(tài)響應(yīng)較快,而且穩(wěn)態(tài)誤差較小。

圖5 無干擾情況下的模型預(yù)測控制效果

4.3.2 下料量不準(zhǔn)確情況

考慮鋁電解槽實際生產(chǎn)中單次下料的下料量并不一定非常準(zhǔn)確,存在下料器堵塞以及人工與槽控機(jī)沖突導(dǎo)致短時內(nèi)多次下料的情況。在100 s時刻將下料量分別增加和減少20%,模擬執(zhí)行器故障和下料沖突的情況。以增加20%下料量為例,控制效果如圖6所示。

圖6 t=100s時刻下料量增加20%的模型預(yù)測控制效果

從圖6可以看出,在100 s時刻,控制系統(tǒng)會自動減小控制下料量,氧化鋁濃度能夠很快回到設(shè)定值,實現(xiàn)快速調(diào)整。此外,由圖6a計算可得下料量減少的平均比例為19.5%,與引入的干擾基本一致,說明該控制系統(tǒng)能夠有效抑制下料不準(zhǔn)確的情況,有較強(qiáng)的魯棒性。

4.3.3 電壓波動情況

考慮到實際生產(chǎn)中,鋁液、電解質(zhì)波動等會引起工作電壓波動,在采樣時刻130～140 s和170～180 s分別引入工作電壓波動,控制效果如圖7所示。

圖7 考慮極距調(diào)整和電壓波動的模型預(yù)測控制效果

由圖7可知,相對于穩(wěn)態(tài)時刻,在干擾引入時刻及持續(xù)時間內(nèi),氧化鋁濃度波動較小,干擾消失后,系統(tǒng)能夠快速恢復(fù)穩(wěn)定,證明了該控制系統(tǒng)具有較強(qiáng)的抗干擾能力。

5 結(jié) 語

針對現(xiàn)有氧化鋁濃度控制的不足,提出了一種基于非線性模型預(yù)測控制的氧化鋁濃度控制方法。仿真實驗表明:本文提出的控制方法能在不依賴于專家經(jīng)驗的基礎(chǔ)上通過下料和極距調(diào)整實現(xiàn)氧化鋁濃度的精確控制,為實際鋁電解槽的氧化鋁濃度控制提供一定的指導(dǎo)作用。目前只進(jìn)行針對鋁電解槽氧化鋁濃度非線性模型預(yù)測控制的仿真研究,后續(xù)將進(jìn)一步進(jìn)行鋁電解槽氧化鋁濃度控制的實際工業(yè)實驗。