鍋爐過熱蒸汽出口壓力的串級專家PID控制

王恒飛 陳永良 郭春裕

(中國計量學院質量與安全工程學院，浙江 杭州 310018)

鍋爐過熱蒸汽出口壓力的串級專家PID控制

王恒飛 陳永良 郭春裕

(中國計量學院質量與安全工程學院，浙江 杭州 310018)

針對常規PID鍋爐過熱蒸汽出口壓力控制效果不佳的問題，提出并編寫了基于PCS7系統的專家PID控制器。運用Matlab軟件進行專家PID的仿真，并利用高級多功能過程控制實訓系統( SMPT 1000) 以及PCS7對控制效果進行了進一步試驗驗證，實驗中PSC7主要負責建立控制系統，實現專家控制。運行結果表明，該控制方案的調節時間短、超調量小，并且具有較好的抗干擾能力，可以實現過熱蒸汽出口壓力的穩定控制。

過熱蒸汽 專家控制 PID控制 PCS7 系統組態

0 引言

鍋爐過熱蒸汽出口壓力是蒸汽供需平衡的標志，是影響系統蒸汽平衡的重要熱力學參數，直接關系到操作的平穩與生產的安全[1]。蒸汽出口壓力控制系統具有時變、高階次、大滯后等特點，利用常規的PID很難對其實現準確的控制，于是引進先進的控制算法來解決此問題。

本文選擇專家PID 算法來解決過熱蒸汽出口壓力過渡過程時間過長這一缺點[2]。專家PID 控制是對常規PID 控制算法的一種改進，能很好地與實際工業相結合。西門子的 Simatic PCS7集散系統實現了機電一體化[3]，其中SCL模塊可以實現編寫專家PID控制算法。通過SCL編寫PID控制與專家系統結合，優化了過熱蒸汽壓力控制器的輸出，大大提高了控制系統的穩定性和控制效果。

1 過熱蒸汽出口壓力控制系統結構

過熱蒸汽出口壓力控制系統就是對出口負荷的壓力進行控制的系統。影響過熱蒸汽出口壓力的變量有燃料流量、過熱蒸汽出口流量以及過熱蒸汽出口溫度等。其中燃料流量是對過熱蒸汽出口壓力影響最顯著的變量，而且可控性好。如圖1所示是壓力串級控制系統結構圖。在串級控制系統中加入專家PID控制器PIC1104，通過主、副環的聯合控制，使整個壓力控制系統的性能提高。

圖1 過熱蒸汽出口壓力控制系統

2 PCS7系統組態

西門子PCS7將傳統的DCS 和PLC 控制系統的優點相結合，系統所有的硬件都基于統一的硬件平臺，可以根據需要選用不同的功能組件進行系統組態[4]。PCS7基于Profibus現場總線技術，把現場設備集成在過程控制系統中，現場設備通過ET200系列分布式I/O連接到Profibus總線上。系統結構圖如圖2所示。

圖2 系統總體結構

PCS7的系統組態是指使用工具軟件對計算機及軟件的各種資源進行配置，使計算機或軟件按照預先設置的指令參數，自動執行指定的任務的過程。PCS7的組態過程一般分為以下四個步驟。

① 新建項目工程

整個項目工程的新建都由工程向導(new project wizard)完成，主要進行CPU套件的選型、項目層級的設置以及創建操作員站(OS)。

② 應用服務器(AS)、OS硬件的組態

項目建立完成之后，需要進行硬件組態，實現SMPT1000(小型過程控制仿真平臺)和PCS7的數據通信。根據過熱蒸汽出口壓力控制系統的需求，將所需采集與輸出的信號線接到ET200對應的I/O接口上。其中AS組態就是將電源、CPU、CP443-1等訂貨號與實際硬件訂貨號匹配，配置CPU412-3H的MPI/DP接口。該接口支持Profibus DP通信，可實現CPU直接與遠程I/O通信。

OS組態與AS類似，需要在工程師站(ES)中組態真實的OS,并且雙方通信成功。

③ 控制算法組態

控制算法組態主要是指連續功能圖(CFC)組態。在CFC編輯器中放置控制系統所需要的功能塊，配置相關的參數。這些功能塊是具有特定功能的預置功能塊，連接相應的端子即可。本文所需要的專家PID功能塊需要使用SCL語言進行編寫，編譯后才能生成專家PID控制器。

④ 監控畫面組態

監控畫面組態就是通過WinCC人機界面組件完成的上位機界面，實現鍋爐過熱蒸汽出口壓力的監控的可視化和過程任務控制。操作員在OS界面可以以圖形化的方式控制和監視出口壓力控制系統的數據變化。

3 專家PID控制器設計

在設計專家控制器時選用了經典的控制器模型，進一步提高專家PID控制器對數據的處理能力，提高最終的控制效果。專家控制通過模擬專家的推理過程和知識，能以專家的水平解決問題。本文使用知識工程方法，應用專家系統的設計規則和實現形式構建一個實時專家智能控制系統[5-6]。專家PID控制器結構圖如圖3所示。

圖3 專家PID控制器結構圖

由圖3看出專家控制器選用的模型為[7]：

U=f(E,T,I)

(1)

式中：U=(u1,u2,…,ul)為控制器輸出集；E={e1,e2,…,em}為控制器輸入集；T={t1,t2,…,tn}為系統中的數據項集；I={i1,i2，…，ip}為推理機構的輸出集。

f為智能算子，其基本形式為：

IfEandTthen (ifIthenU)

(2)

專家PID控制器就是依據這個模型理論設計的。圖4所示的是典型的二階系統單位階躍響應誤差曲線[2]。

圖4 典型二階系統單位階躍響應誤差曲線

e(t)的采樣誤差e(k)、e(k-1)、e(k-2)分別表示前一次和前兩次采樣時刻的誤差值，則有:

Δe(k)=e(k)-e(k-1)

Δe(k-1)=e(k-1)-e(k-2)

根據誤差及其變化，并結合圖4對專家PID控制器的規則進行如下定性分析。

① 當|e(k)|>M1時，說明誤差的絕對值已經很大。不論誤差變化趨勢如何，都應考慮控制器的輸出按定值輸出，以達到迅速調整誤差、同時避免超調的目的。此時，它相當于實施開環控制?？刂破鬏敵鰹閡(k)=usp。

② 當e(k)Δe(k)>0時，說明誤差在朝誤差絕對值增大方向變化，或誤差為某一常數，未發生變化。

如果|e(k)|≥M2，說明誤差較大，可考慮由控制器實施較強的控制作用，使誤差絕對值朝減小方向變化，迅速減小誤差的絕對值?？刂破鬏敵鰹?

u(k)=u(k-1)+K1{KpΔe(k)+KIe(k)+KD[e(k)-2e(k-1)+e(k-2)]}

(3)

如果|e(k)|≥M2，說明盡管誤差朝絕對值增大方向變化，但誤差絕對值本身并不是很大，可考慮實施一般的控制作用，扭轉誤差的變化趨勢，使其朝誤差絕對值減小方向變化?？刂破鬏敵鰹?

u(k)=u(K-1)+KpΔe(k)+KIe(k)+KD[e(k)-2e(k-1)+e(k-2)]

(4)

③ 當e(k)Δe(k)<,Δe(k)Δe(k-1)>0或Δe(k)=0時，說明誤差的絕對值朝減小的方向變化，或者已經達到平衡狀態。此時，可考慮采取保持控制器輸出不變。

④ 當e(k)Δe(k)<0、Δe(k)Δe(k-1)<0時，說明誤差處于極值狀態。如果此時誤差的絕對值較大，即可考慮實施較強的控制作用；如果誤差的絕對值較小即|e(k)|

e(k)=u(k-1)+K1Kpem(k) |e(k)|>M2

(5)

e(k)=u(k-1)+K2Kpem(k) |e(k)|≤M2

(6)

⑤ 當|e(k)|≤ε(精度)時，說明誤差的絕對值很小，此時可加入積分環節，減少穩態誤差。

以上各式中，em(k)為誤差的第k個極值；u(k)為第k次控制器的輸出；K1為增益放大系數，K1>1；K2為抑制系數，0M2>0。

以上規則就是通過對誤差控制信號的瞬時響應進行模式識別。在常規PID控制器的基礎上構建專家PID控制器，可以更好地優化主控制器的輸出，充分發揮專家控制的優點，能對時變、非線性和易受到各種干擾的受控過程給出有效的控制策略。串級專家PID控制系統方塊圖如圖5所示。

圖5 控制系統方塊圖

4 仿真與驗證

4.1 對象特性辨識

在進行Matlab仿真之前，需要獲取最接近被控對象的精確數學模型。運用WinCC控件中的Visual Basic的腳本語言讀取壓力模型的輸入(燃油流量)與輸出(蒸汽壓力),再利用行列式比定階法來估計壓力模型的階次；然后利用遞推最小二乘法實現對壓力模型的系統辨識；最終確定最精確的過熱蒸汽出口壓力的數學模型[8-10]。在 Matlab 中編寫模型階次辨識語言，對多組壓力模型的輸入輸出進行計算，得到行列式比DR*(n)，從而實現對過熱蒸汽出口壓力模型階次的辨識。行列式比數據如表1所示。

表1 行列式比數據

通過數據比較可知，當n等于2時行列式比顯著增加，故壓力模型的階次可定為 2 階。

獲得壓力模型的階次之后，進一步分析論證得到最接近的壓力辨識模型的差分方程為：

y(k)=-1.011 05y(k-1)+0.025 6y(k-2)- 0.000 6u(k-1)+0.000 5u(k-2)

(7)

4.2 專家PID控制的仿真

根據式(7)的模型，先在Matlab上進行仿真，分別進行常規PID與專家PID的仿真，PID控制器是通過數字PID控制算法程序實現的[11]。整定的專家PID控制參數為：

KP=5.4KI=0.8KD=30K1=1.2K2=0.6 利用Matlab繪出各自的壓力曲線如圖6所示。

圖6 常規PID與專家PID控制的輸出比較

由圖6的輸出曲線可以看出，專家PID控制過渡時間短且基本無超調。仿真結果表明專家PID 的控制效果更加優越。

4.3 實驗驗證

由仿真結果可以明顯看出專家PID控制器的優越性，再根據前面所述的專家PID控制器的設計規則，在PCS7中運用SCL模塊編寫專家PID控制器功能塊進行實驗驗證。生成的專家控制器功能塊如圖7所示。

圖7 專家PID控制器功能塊

在PCS7中將整個鍋爐過熱蒸汽出口壓力控制系統搭建完成后，副回路控制器使用燃料流量控制器，主回路控制器分別采用常規和專家 PID控制器，準備工作完成以后運行上位機程序進行實驗。兩種PID控制的輸出結果比較如表2所示。

表2的對比結果進一步表明，專家PID 控制使控制器的輸出更加穩定，專家PID 控制在快速性、抗干擾性等方面優于常規PID 控制。

表2 常規/專家PID控制的比較

5 結束語

本文首先利用Matlab搭建的出口蒸汽壓力的數學模型，驗證了專家系統在PID控制器運用中的優越性，并通過西門子先進的PCS7系統進行進一步實驗驗證。與常規的PID控制相比較，專家PID控制提高了系統的抗干擾性能，在實際生產環節中可以起到節能減排、提高經濟效益的作用。

[1] 馬昕,張貝克.深入淺出過程控制：小鍋帶你學過控[M].北京：高等教育出版社,2013:175-177.

[2] 劉金琨.先進PID 控制及其Matlab仿真[M].北京: 電子工業出版社,2003.

[3] 西門子 A&D 集團.用于小型應用的緊湊而經濟的過程控制系統—SIMATIC PCS 7BOX[P].北京:西門子(中國)有限公司,2005.

[4] Sandro C. Control system PCS7 and M. I. S. together for the complete automation of the process in the sugar beet factory of Co. Pro. B.-Minerbio-Italy[J].Computer Aided Chemical Engineeringr,2007,24:842-845.

[5] 吳興純,趙金燕,楊秀蓮，等.智能 PID算法在爐溫度控制系統中的運用[J].機電工程, 2011(8):949-950.

[6] 郭楠,李智.專家 PID 算法在伺服系統中的應用與仿真[J].機械工程與自動化,2009(6):61-63.

[7] 韓力群.智能控制理論及應用[M].北京：機械工業出版社,2007.

[8] Ljung L,Soderstrom T.Theory and practice of recursive identification[M].MIT Press,Cambridge,Massacbusetts,London,England, 1983.

[9] 方崇智,蕭德云.過程辨識[M].北京：清華大學出版社,1988.

[10]潘立登,潘仰東.系統辨識與建模[M].北京：化學工業出版社,2004.

[11]葛寶明,林飛,李國國.先進控制理論及其應用[M].北京:機械工業出版社,2007.

Cascade Expert PID Control of Boiler Superheated Steam Outlet Pressure

To solve the problem of conventional PID, i.e., poor control effect for boiler superheated steam outlet pressure, the expert PID controller based on PCS7 system is proposed and compiled. The expert PID is simulated by adopting Matlab, and further test and verification of the control effect are conducted with advanced multifunction process control training system SMPT 1000 and PCS7, PCS7 is mainly responsible for setup the control system to implement expert control. The results of tests and operation show that the control scheme features shorter regulation time, smaller overshoot, and possesses high anti-interference capability, it is able to realize stable control for superheated steam outlet pressure.

Superheated steam Expert control PID control PCS7 System configuration

王恒飛(1990-)，男，現為中國計量學院檢測技術與自動化裝置專業在讀碩士研究生；主要從事計算機控制、電子測量及虛擬儀器技術的研究。

TH81;TP273

A

10.16086/j.cnki.issn1000-0380.201510025

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