摘要：針對(duì)目前火電廠熱控系統(tǒng)中大量dcs參數(shù)運(yùn)行調(diào)試現(xiàn)狀，提出了采用matalb進(jìn)行火電廠熱控系統(tǒng)dcs系統(tǒng)參數(shù)優(yōu)化的實(shí)現(xiàn)方案，介紹了matalb提供的非線性優(yōu)化工具模塊的使用技術(shù)和控制系統(tǒng)的參數(shù)優(yōu)化實(shí)例。
1引言
目前火電廠熱控系統(tǒng)已廣泛采用分散控制系統(tǒng)dcb。dcs控制系統(tǒng)基本是由一些簡(jiǎn)單的標(biāo)準(zhǔn)結(jié)構(gòu)組態(tài)而成，如將基本的和各種變形的pid調(diào)節(jié)器、微分器、一階慣性環(huán)節(jié)等’組態(tài)構(gòu)成的鍋爐三沖量給水控制系統(tǒng)、鍋爐過(guò)熱汽溫調(diào)節(jié)系統(tǒng)，以及機(jī)組協(xié)調(diào)控制系統(tǒng)等。dc$所構(gòu)成的控制系統(tǒng)的拓?fù)浣Y(jié)構(gòu)特點(diǎn)是各基本單元簡(jiǎn)單而標(biāo)準(zhǔn)化，復(fù)雜功能的實(shí)現(xiàn)通過(guò)用標(biāo)準(zhǔn)基本單元的復(fù)雜連接而完成。這使得dcs環(huán)境下的控制系統(tǒng)具有可任意組態(tài)的特點(diǎn)，并使得dcs系統(tǒng)在實(shí)際使用到具體工程時(shí)控制系統(tǒng)結(jié)構(gòu)和參數(shù)才*確定下來(lái)。由于各個(gè)具體工程系統(tǒng)的設(shè)計(jì)與安裝及運(yùn)行往往都不*一致，影響了許多成功的控制器參數(shù)優(yōu)化技術(shù)和軟件包在不同dcs實(shí)際工程環(huán)境的推廣應(yīng)用。并且由于實(shí)際工程調(diào)試時(shí)間緊，控制系統(tǒng)運(yùn)行責(zé)任重大，工程經(jīng)初步調(diào)試和168小時(shí)投運(yùn)后，控制系統(tǒng)幾乎就不再做參數(shù)優(yōu)化試驗(yàn)，控制系統(tǒng)的優(yōu)化性能常常沒(méi)有充分發(fā)揮[1]。
—般說(shuō)來(lái)，控制系統(tǒng)的優(yōu)化問(wèn)題大致可分為兩類：一類是調(diào)節(jié)器的型式已經(jīng)確定（如dcs中的pi、pid調(diào)節(jié)規(guī)律）時(shí)，通過(guò)調(diào)整控制器參數(shù)，使某目標(biāo)函數(shù)達(dá)到*：另一類為控制對(duì)象已知，尋找*控制作用使某目標(biāo)函數(shù)達(dá)到*，包括尋找控制器的結(jié)構(gòu)、型式及參數(shù)等。火電廠熱控系統(tǒng)采用dcs時(shí)的結(jié)構(gòu)和調(diào)節(jié)器的類型是確定的，故控制系統(tǒng)的優(yōu)十七就是整定調(diào)節(jié)器參數(shù)。調(diào)節(jié)器參數(shù)的整定除可以采用傳統(tǒng)的頻域分析法外，也可用參數(shù)優(yōu)化的方法予以解決[2]。
參數(shù)優(yōu)化有兩種途徑：一種是間接尋優(yōu)，在得到目標(biāo)函數(shù)解析式的基礎(chǔ)上，根據(jù)該目標(biāo)函數(shù)取極值的充分必要條件，計(jì)算出參數(shù)的*解；另一種為直接尋優(yōu)，即直接在參數(shù)空間中，按照一定的規(guī)律進(jìn)行探索尋優(yōu)，尋得的目標(biāo)函數(shù)即為*的參數(shù)點(diǎn)?；痣姀S熱控系統(tǒng)由于存在大量的非線性和純遲延，解析分析計(jì)算較困難，本文將介紹利用matalb軟件包中非線性優(yōu)化工具ncd（nonlinem·，controldesign），進(jìn)行火電廠熱控系統(tǒng)dcs中pid控制系統(tǒng)參數(shù)優(yōu)化的方法。在介紹matalb提供的非線性優(yōu)化工具模塊的使用技術(shù)的基礎(chǔ)上，結(jié)合系統(tǒng)非線性、純遲延和變參數(shù)特點(diǎn)，給出將該技術(shù)應(yīng)用于火電廠風(fēng)量控制雙執(zhí)行機(jī)構(gòu)系統(tǒng)的參數(shù)優(yōu)化實(shí)例。
2采用matalb實(shí)現(xiàn)火電廠熱控非線性系統(tǒng)優(yōu)化
matalb已成為通用的控制系統(tǒng)分析工具[3兒它不但有用于動(dòng)態(tài)系統(tǒng)仿真的simulink工具箱，還開發(fā)了一個(gè)于非線性控制和優(yōu)化設(shè)計(jì)的工具箱ncd（nonlinearcontroldesign）。借助于ncd工具箱，不但可以利用simulink進(jìn)行系統(tǒng)的動(dòng)態(tài)仿真與分，析，而且可以利用simulink進(jìn)行系統(tǒng)參數(shù)的優(yōu)化設(shè)計(jì)。非線性系統(tǒng)優(yōu)化設(shè)計(jì)可以參考demo中的典型示例進(jìn)行。在創(chuàng)建火電廠熱控應(yīng)用系統(tǒng)時(shí)，可將典型示例的模塊拷貝到該應(yīng)用系統(tǒng)中使用。
某電廠300mw機(jī)組采用西門子t／mpdcs系統(tǒng)。機(jī)組投運(yùn)agc時(shí)需要對(duì)各級(jí)子系統(tǒng)進(jìn)行優(yōu)化，以便能在安全約束下較快地響應(yīng)電網(wǎng)負(fù)荷需求、dcs中的送（引）風(fēng)系統(tǒng)為雙執(zhí)行機(jī)構(gòu)，其簡(jiǎn)化仿真系統(tǒng)如圖1。送（引）風(fēng)a、b兩側(cè)回路由偏置系數(shù)c2協(xié)調(diào)，實(shí)現(xiàn)多輸出控制（mocs）。圖中fd-a和m—b分別表示a側(cè)和b側(cè)的調(diào)節(jié)通道，仿真中用有區(qū)別的小慣性環(huán)節(jié)實(shí)現(xiàn)。風(fēng)量對(duì)象的模型圖中用plant表示（plant的實(shí)際特性如圖2）。圖1上部方框（u’，data-inputl）為在電廠采集的實(shí)際運(yùn)行的風(fēng)量指令數(shù)據(jù)，可經(jīng)開關(guān)叨換送入plant模型得到仿真輸出風(fēng)量信號(hào)，將該仿真輸出風(fēng)量信號(hào)與實(shí)際的輸出風(fēng)量信號(hào)進(jìn)行比較可以驗(yàn)證實(shí)驗(yàn)?zāi)Ｐ?，并提供?yōu)化結(jié)果可用與否的參考。fd644ini為matlab優(yōu)化的初始化文件，而ncd-outpu就是matlab的非線性仿真優(yōu)化工具。
利用ncd工具箱?；痣姀S熱控系統(tǒng)的參數(shù)優(yōu)十t任務(wù)可很方便的進(jìn)行。ncd的使用方法是：在matalbcommand窗口下鍵人ncdblock，即可彈出nonlinearcontroldesignblockset模塊庫(kù)。這里有一個(gè)用于優(yōu)化的模塊ncdoutput．把它連接到要優(yōu)化的用simulink建立的應(yīng)用系統(tǒng)框圖模型的輸出端，啟動(dòng)仿真，則系統(tǒng)的參數(shù)優(yōu)化設(shè)計(jì)自動(dòng)進(jìn)行。
為了介紹優(yōu)化設(shè)計(jì)的方法和步驟，下面結(jié)合風(fēng)量系統(tǒng)實(shí)例來(lái)說(shuō)明。
3ncd優(yōu)化技術(shù)
設(shè)圖1火電廠風(fēng)量控制對(duì)象plant的傳遞函數(shù)可表達(dá)式為：
此處設(shè)系統(tǒng)包含飽和環(huán)節(jié)和速度限制環(huán)節(jié)兩個(gè)非線性環(huán)節(jié)，以及一個(gè)純遲延環(huán)節(jié)。再設(shè)系統(tǒng)參數(shù)具有不確定因素：a2在40∽50之間變化，a1在1．5∽4．5范圍內(nèi)變化。運(yùn)行中可變參數(shù)平衡點(diǎn)為a2=43，a1=3。
控制系統(tǒng)采用pid調(diào)節(jié)時(shí)的優(yōu)化設(shè)計(jì)任務(wù)為：設(shè)計(jì)一組pid控制器參數(shù)，使該單位反饋閉環(huán)系統(tǒng)滿足：zui大超調(diào)量不大于20％；上升時(shí)間不大于los；調(diào)整時(shí)間不大于30s。并且還希望閉環(huán)系統(tǒng)對(duì)不確定因素a1、a2具有魯棒性。
優(yōu)化設(shè)計(jì)的方法和步驟主要如下：
（1）用simulink建立控制對(duì)象模型（包括速度執(zhí)行機(jī)構(gòu)的限幅環(huán)節(jié)、飽和環(huán)節(jié)和對(duì)象傳遞函數(shù)）如圖2；
（2）在matlabcommand窗口中輸人命令ncdblock；將ncdoutput模塊與step階躍輸入、pid等一起構(gòu)成控制系統(tǒng)，建立如圖l所示的用于參數(shù)優(yōu)化設(shè)計(jì)的simulink模型方框圖。其中的控制對(duì)象為圖2的子系統(tǒng)，pid為simulink提供的標(biāo)準(zhǔn)模塊。實(shí)際應(yīng)用中可根據(jù)dcs算法定制。
（3）雙占ncdoutput模塊，彈出ncdblockset約束窗口，如圖3所示
（4）選擇options菜單：
①通過(guò)stepresponse命令定義階躍響應(yīng)性能限制：調(diào)整時(shí)間30s；上升時(shí)間10s；穩(wěn)態(tài)誤差百分?jǐn)?shù)5；超調(diào)量百分?jǐn)?shù)20：振蕩負(fù)幅值百分?jǐn)?shù)1。階躍響應(yīng)性能限制也可以直接用鼠標(biāo)在ncdblockset約束窗口設(shè)置。
②選擇timerange命令，設(shè)置x軸0—100s。
③選擇輸出坐標(biāo)y—axis命令，設(shè)置y軸幅值0．131—1．321。
（5）選擇optimization（優(yōu)化）菜單
①選擇pammeters命令，定義待優(yōu)化調(diào)整變量及有關(guān)參數(shù)如下：待調(diào)整優(yōu)化變量kp，ki，kd及其上下限，變量允差（0．001）和約束允差（0．001）。
②選擇uncertainty命令，定義不確定變量及有關(guān)參數(shù)：不確定變量為ala2及其上下限[al／240]、[a1*250]
（6）利用初始化模塊fd644ini對(duì)系統(tǒng)模型參數(shù)初始化：kp=0．6；ki=0．05；kd=2；a2=43；al=3
（7）選擇optimization菜單start命令，開始對(duì)調(diào)整變量的優(yōu)化。優(yōu)化過(guò)程實(shí)質(zhì)上是利用matlab優(yōu)化工具箱函數(shù)constr．m對(duì)調(diào)整變量kp，ki，kd進(jìn)行優(yōu)化計(jì)算，直到系統(tǒng)階躍響應(yīng)指標(biāo)的約束條件得到滿足為止的過(guò)程。優(yōu)化時(shí)ncdblockset約束窗口不斷顯示階躍響應(yīng)曲線，matlabcommand窗口也不斷顯示有關(guān)信息，一旦性能約束條件滿足
，優(yōu)化過(guò)程停止。
（8）優(yōu)化結(jié)束后，在matlabcommand窗口，鍵人優(yōu)化整定的變量名，即得到優(yōu)化的參數(shù)值：
kp=2．3449；ki=0．1115；kd=10．3382
ncdblockset約束窗口如圖3，初始階躍響應(yīng)為曲線1，控制器的優(yōu)化的系統(tǒng)階躍響應(yīng)為曲線2。將該組優(yōu)化的pid參數(shù)與dcs系統(tǒng)中現(xiàn)有的參數(shù)比較，則可以根據(jù)實(shí)際工況修正原參數(shù)，從而得到優(yōu)化的系統(tǒng)。
若風(fēng)機(jī)a、b側(cè)的pid參數(shù)需要分別優(yōu)化時(shí)，則可將圖l中：pid的參數(shù)分別用：kpl、ki1、kdl、kp2、ki2、kd2表示，在參數(shù)初始化時(shí)作相應(yīng)考慮，即可得到優(yōu)化結(jié)果。若將圖二的復(fù)雜模型變換為電廠熱控中常用的對(duì)象模型，即將plant的模型僅用多個(gè)慣性環(huán)節(jié)串聯(lián)表示時(shí)，優(yōu)化的方法和步驟仍基本相同，只是免去了參數(shù)的不確定性和非線性。
4結(jié)束語(yǔ)
本文在分析現(xiàn)有dcs中pid控制系統(tǒng)的參數(shù)優(yōu)化技術(shù)的基礎(chǔ)上，提出了采用matlab軟件包進(jìn)行火電廠dcs中pid參數(shù)非線性優(yōu)化的方法，給出了采用matlab的非線性優(yōu)化模塊進(jìn)行控制系統(tǒng)優(yōu)化的實(shí)例，結(jié)合火電廠送風(fēng)雙執(zhí)行機(jī)構(gòu)系統(tǒng)進(jìn)行了討論。應(yīng)用該優(yōu)化技術(shù)時(shí)，應(yīng)注意仿真尋優(yōu)對(duì)初值的依賴，ncd優(yōu)化給出的是初值附近的局部*解而不是全局*解。在某電廠300mw機(jī)組ccs系統(tǒng)細(xì)調(diào)實(shí)現(xiàn)機(jī)組agc的工作中，采用mayubncd的輔助優(yōu)化方法取得了成效。
參考文獻(xiàn)
1李平康，dcs中的pid控制器參數(shù)自整定技術(shù)初探，華北電力技術(shù)，1999．10
2韓樸、朱希彥，自動(dòng)控制系統(tǒng)數(shù)字仿真，中國(guó)電力出版社，1996．1，p：213-254
3龔劍、朱亮，matlab5．x入門與提高，清華大學(xué)出版社，2000．3
4李遵基，熱工自動(dòng)控制系統(tǒng)，中國(guó)電力出版社，1997．10，p·126-129