陽極焙燒爐爐溫智能控制算法研究

李 麗

（中國煤炭科工集團太原研究院有限公司）

在工業焙燒爐控制過程中，由于系統本身存在非線性、大時滯、大時變及溫度、壓力、氣氛之間的相互耦合等，且該系統為典型的分布參數系統［1］，因此難以建立簡化的數學模型。若采用傳統控制方法，難以克服系統的滯后、時變和干擾，無法保證控制系統的魯棒性，要實現對生產過程的最佳控制是很難的。

焙燒過程分為4個階段：低溫預熱階段、中溫焙燒階段、高溫燒結階段和冷卻階段［2～4］。 應針對焙燒不同階段，制定相應的控制策略。 陽極焙燒爐分為4個區域： 預熱區 （1P～3P爐室）、 加熱區（4P～6P爐室）、自然冷卻區（7P～10P爐室）和強制冷卻區（11P爐室）。 由于加熱區控制效果的好壞對焙燒制品的質量和能耗的影響最大，因此筆者重點研究加熱區溫度控制，提出采用零極點配置PID預測控制算法來解決加熱區爐溫控制的難題。

1 加熱區控制特性分析

加熱區是整個焙燒過程中最重要的焙燒階段，包括4P～6P共3個爐室，每個爐室配備一個機電一體化燃燒架，燃燒架利用重油燃料燃燒的熱量對陽極進行強制加熱，生陽極在此區被加熱到工藝要求的溫度。

加熱區數學模型輸入和輸出之間的關系如圖1所示。 不難看出， 這是一個雙輸入雙輸出系統。 兩個輸出的作用是相互耦合的，要找到兩者之間的關系實現解耦控制較為困難，因此對模型進行簡化。

圖1 加熱區焙燒過程的輸入輸出關系

影響加熱區爐溫的擾動因素很多， 主要有：排煙總管抽力的變化和排煙閥門定位的準確性f1，爐體變形漏風量f2，重油熱值變化及油壓、油溫波動f3，入口煙氣溫度的變化量f4，陽極制品的溫度和陽極制品蓄熱回收程度f5，爐體散熱量。

在加熱區影響爐溫的主要因素是重油流量，排煙閥開度對爐溫的影響是通過負壓的改變來體現的，而且在加熱區負壓相對于重油流量而言對爐溫的影響要小。 調節爐室負壓的作用，第一是保證燃料的充分燃燒，第二是作為預熱區爐室的控制量來保證預熱區爐室升溫。 因此，當調節重油流量能夠滿足升溫要求時，將負壓作為可測擾動量f7來處理，對其進行前饋補償，系統簡化為圖2所示的模型， 此時系統變為一個單入單出系統。 圖2中的Wo2（s）為爐室負壓數學模型。 當只調節重油流量難以滿足升溫要求時，將負壓作為控制量來處理，系統簡化為圖3所示的模型，此時系統變為一個雙入單出系統。

圖2 負壓作前饋補償的加熱區焙燒過程的輸入輸出關系

圖3 負壓作控制量的加熱區焙燒過程的輸入輸出關系

由圖3可知，在加熱區，擾動f1和f2對爐溫的影響主要由負壓調節器構成的控制回路來克服，擾動f3～f6對爐溫的影響主要由爐溫調節器構成的控制回路來克服。

2 加熱區陽極焙燒爐模型的建立

2.1 模型的類型和結構

焙燒爐的數學模型實質上為一個分布參數的慢時變系統，為簡化模型，在加熱區可以將工作點附近的模型線性化， 用一個以爐溫為輸出、重油流量為輸入的集中參數離散差分方程模型來表示。 考慮到負壓變化對爐溫有很大的影響，故在建模時， 可將負壓作為可測干擾量來處理，再考慮到現場隨機噪聲干擾的影響，這樣就可以用帶可測負壓干擾、有控制項的自回歸滑動平均模型（CARMA）來描述被控對象的特性，即：

其中，u（k）、q（k）為k時刻系統的輸入和輸出；v（k）、ξ（k）為k時刻可測干擾和隨機噪聲；d1、d2為控制量和可測干擾通道的時延；z-1為單位后移算子；A、B、C為傳遞系數。

2.2 模型參數的辨識

根據對被控對象機理模型的研究，被控對象可近似認為是一個二階時變滯后系統。 因此，可以確定傳遞函數Gu（z-1）的模型階次為2，一般有模型階次na≥nb，na≥nc。 因此，不妨設na=nb=nc=n=2，若nb=0或nc=0，則辨識出模型中的相應系數為0。

系統的時延參數d1和d2可以用以下方法確定。 確定系統輸入變量u（k）的時延τ1和擾動變量v（k）的時延τ2的大致范圍，則d1=τ1/Ts和d2=τ2/Ts的范圍可確定，Ts為時間。 根據對現場數據的分析，確定0＜d1≤3，0＜d2≤2。d1和d2分別取不同的值，利用時變遺忘因子的遞推最小二乘法對模型參數進行辨識。 使損失函數J取得最小值的時延d1和d2，即最佳時延。

采集控制系統歷史數據庫中的數據，利用遞推最小二乘法，用仿真軟件編制程序進行離線建模。 應用損失函數檢驗法和殘差白性檢驗法驗證模型的階次并確定系統時延。 經過大量的仿真分析，當n=2，d1=4，d2=2時損失函數J取得最小值。

因此可求得當Ts為30 s時脈沖通道傳遞函數Gu（z-1）和擾動通道傳遞函數Gv（z-1）分別為：

3 零極點配置PID預測控制

從前文可知，陽極焙燒爐的運行是一個具有大滯后的被控工業生產過程， 且其周圍環境復雜，條件惡劣，受到揮發分燃燒、漏風及爐體散熱等多種因素影響，因此需要設計一個能夠應付各種情況，并能獲得較好控制效果的控制器，以保證生產的平穩進行。 基于這一考慮，筆者采用了零極點配置PID預測控制器， 并對該控制器進行了改進，實現了負壓和爐溫的協調控制。

3.1 前饋補償功能的零極點配置PID預測控制器設計

建立陽極焙燒爐在加熱區的數學模型時，已把負壓波動當作可測干擾來處理。 根據復合控制理論，在系統中設計一前饋控制器，組成前饋與反饋復合控制器， 可有效抑制可測干擾的影響。有前饋補償的零極點配置PID預測控制系統的結 構如圖4所示。

圖4 有前饋補償的零極點配置PID預測控制結構框圖

圖4中，ysp（k）為爐溫設定值；uv（k）為前饋控制器的輸出值； Gv（z-1）為擾動通道的傳遞函數；Gξ（z-1）為隨機噪聲通道模型傳遞函數；GM（z-1）為前饋控制器的傳遞函數；Gc（z-1）為控制器的傳遞函數；uo（k）=u（k）-uv（k），為控制器Gc（z-1）和前饋控制器GM（z-1）輸出量的疊加值；uo′（k）為經過輸出限幅單元的輸出。

從圖4中可以看出， 擾動作用通過兩條通道作用于系統，一條是擾動通道，另一條是前饋補償通道，前饋控制器的設計應該滿足，當有擾動時，上述兩條通道對輸出爐溫波動的影響相互抵消，所以有：

當滿足式（4）時，擾動作用對輸出的影響被前饋控制器完全補償了。 前饋控制器的輸出為：

若要實現對擾動作用的動態靜態補償，則要求B（z-1）的零點位于單位圓內，即要求被控過程為一最小相位系統。 此外，還要求d2≥d1，否則要求前饋控制器給出前導控制作用，這在物理上是無法實現的。 根據陽極焙燒爐的數學模型可知，系統控制通道和擾動通道均為最小相位系統，滿足前饋控制器穩定的要求。 因此時d1=4，d2=2，無法獲取k+2時刻的擾動，故在式（5）中令d1=d2，得到可實現的前饋控制器的輸出為：

因此所設計的前饋控制器，不能實現對爐室負壓波動所引起爐溫變化的完全補償。

3.2 前饋補償和協調控制功能的零極點配置PID預測控制

由于負壓對于其他爐室來說是一個可測不可控的擾動量，若將負壓看作其他爐室控制的一個可調節的輸入量，則必然存在一個協調控制的問題。

筆者通過對陽極焙燒工藝的詳細分析，設計了可以協調控制各爐室爐溫和負壓的控制器。 下面以4P、5P爐室的智能控制算法為例來介紹。

設計的具有前饋補償和協調控制功能的智能控制器的基本思想是，當uo（k）∈［umin，umax］時，將負壓v（k）當作可測不可控的擾動量來處理，由前饋控制器的輸出uv（k）來對其進行補償；當uo（k）超出了其工藝允許范圍時，令控制量uo（k）等于下邊界值umin或上邊界值umax， 并將負壓v（k）當作可以調節的控制量來處理，此時負壓v（k）的調節算法如下：

a. 如uo（k）

b. 將此時刻的爐溫設定值ysp（k）和umin代入焙燒爐的數學模型y（k）=Gu（z-1）u（k）+Gv（z-1）v（k）中，求得實際輸出的爐溫控制值vsp（k）。

c. 對vsp（k）進行輸出限幅。

圖5 4P爐室負壓協調控制器結構框圖

圖5中傳遞函數Wc（z-1）采用增量式數字PID設4P爐室通過調節負壓達到設定值v′sp（k），排煙閥門開度控制信號為u4P（k）。 同理，加熱區5P和6P爐室如需要協調控制負壓時排煙閥門的開度控制信號為u5P（k）和u6P（k）。 執行負壓協調算法時，K1閉合K2斷開；不執行負壓協調算法時，K2閉合K1斷開。

e. 因為工藝1P爐室的升溫要求較加熱區低，特別是在初級階段，瀝青處于軟化狀態，并不進行縮聚反應，適當減低升溫速率，延長升溫時間，對炭塊的質量影響不大［3］。因此，為保證加熱區爐室的升溫可以適當地降低預熱區的負壓控制要求，將各爐室對排煙閥門的開度要求按照一定的權重系數相疊加，以疊加值作為控制排煙閥門開度的實際控制量，可用下式來表示：

當uo（k）>umax時，uo（k）=umax，采用上述相同的算法，在此不再敘述。

調節負壓對預熱區升溫過程影響很大，故λ1、λ2較大，由于負壓波動對6P爐室爐溫影響較小，故權重系數λ5要設置得小一點。

由前面的討論可知，陽極焙燒爐為一個參數慢時變系統，需要在線采集輸入輸出數據對過程的模型進行在線修正。 筆者采用時變遺忘因子的遞推最小二乘法進行模型的在線辨識與修正，周期為30 s。

4 仿真分析

對4P爐室的升溫工藝曲線進行仿真，升溫起點溫度設定為830 ℃， 終點溫度設定為1 030 ℃，升溫時間為25 h，保溫時間為3 h，控制周期為0.5 min。 實際系統中干擾量爐室負壓應為一隨機變化量， 為易于反映前饋控制器的補償效果，取負壓v（k）為跳變的方波信號，范圍為-90～-60 Pa，跳變周期為50 min，如圖6所示。圖7為加入前饋控制器前后的溫升曲線。

圖6 爐室負壓階躍擾動信號

圖7 4P爐室加入前饋控制器前后的溫升曲線

圖8 4P爐室智能算法與普通PID算法控制效果比較

從圖8可以看出，加入前饋控制器之后，大幅降低了輸出波動，爐溫能較精確地跟蹤設定的溫升曲線，加入負壓階躍擾動后，未加前饋控制器時爐溫與設定溫度偏差在［-25 ℃，25 ℃］范圍內波動，加入前饋控制器后爐溫與設定溫度偏差在［-5 ℃，5 ℃］的范圍內波動，而且對模型參數的變化也有很好的抑制能力。

5 結束語

詳細分析了加熱區控制系統的特點，基于工作點附近的CARMA模型，采用帶時變遺忘因子的遞推最小二乘法辨識得到加熱區陽極焙燒爐的數學模型，并以此模型為基礎，在加熱區爐室設計了具有前饋補償和協調控制功能的零極點配置PID預測智能控制器。 通過仿真驗證證明了筆者提出的智能控制算法可以較好地解決爐室負壓波動對爐溫造成的干擾問題，從而實現加熱區爐溫的精確控制。