專利名稱：一種內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)及方法
技術(shù)領域：
本發(fā)明涉及精餾節(jié)能控制系統(tǒng)與方法設計問題，尤其提出了一種內(nèi)部熱耦合精餾 塔的自適應非線性控制系統(tǒng)及方法。
背景技術(shù)：
精餾過程是化工過程中的一種核心過程，精餾塔是其中的一個關(guān)鍵單元。長久以 來，精餾塔因為高耗能，低能效的問題成為國際精餾領域研究的焦點。目前針對精餾過程的 能耗問題主要有兩方面的解決方案一種設計新型結(jié)構(gòu)，利用熱量耦合實現(xiàn)能量重復利用 達到節(jié)能目的，一種設計高效精餾過程控制策略，提高產(chǎn)品生產(chǎn)質(zhì)量減少廢料從而達到節(jié) 能目的。盡管有較多實驗研究證明內(nèi)部熱耦合精餾塔能夠顯著提高能源利用率，但是由于 內(nèi)部熱耦合精餾塔的精餾段與提餾段之間存在極強的耦合性且該塔具有十分復雜的強非 線性，該塔的控制策略設計顯得尤為困難。傳統(tǒng)的PID，內(nèi)膜控制方案等已經(jīng)不能滿足要求，在內(nèi)部熱耦合精餾過程控制當 中，這些方案已經(jīng)很難使精餾過程穩(wěn)定。而基于線性辨識模型的預測控制方案只能工作在 穩(wěn)態(tài)工作點附近，稍微增大干擾幅度，或者設定值階躍變化系統(tǒng)控制質(zhì)量則出現(xiàn)明顯下降。 實現(xiàn)內(nèi)部熱耦合精餾塔高效節(jié)能過程的有效的非線性控制方案具有十分重要的意義。
發(fā)明內(nèi)容
為了克服現(xiàn)有的內(nèi)部熱耦合精餾塔的控制方法的抑制干擾能力差、控制效果差的 不足，本發(fā)明提供一種抑制干擾能力良好、控制效果好的內(nèi)部熱耦合精餾塔的自適應非線 性控制系統(tǒng)及方法。本發(fā)明解決其技術(shù)問題所采用的技術(shù)方案是一種內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)，包括與內(nèi)部熱耦合精餾塔直接 連接的現(xiàn)場智能儀表和DCS系統(tǒng)，所述DCS系統(tǒng)包括存儲裝置、控制站和上位機，所述現(xiàn)場 智能儀表與存儲裝置、控制站和上位機相連，所述的上位機包括用以計算輸出內(nèi)部熱耦合 精餾塔控制變量值的自適應非線性控制器，所述自適應非線性控制器包括組分推斷模塊，用以從現(xiàn)場智能儀表獲取溫度，壓強數(shù)據(jù)，計算內(nèi)部熱耦合精餾塔 各塊塔板的組分濃度，并將組分濃度計算結(jié)果存儲在歷史數(shù)據(jù)庫當中，采用式(1) (2)得 到 其中，k為當前采樣時刻，下腳標i為塔板編號，1為塔頂編號，f為進料板編號，η 為塔底編號，Xi (k)為k采樣時刻液相輕組分濃度，PJk)為k采樣時刻精餾段壓強、PS提餾段壓強，Ti(k)為k采樣時刻塔內(nèi)各塊塔板的溫度，α為相對揮發(fā)度，a、b、c為安東尼常數(shù)；模型參數(shù)自適應校正模塊，用以實現(xiàn)模型參數(shù)的在線更新，在線擬合模型函數(shù)如 下⑶、⑷ 其中，友為第i塊塔板處液相組分濃度預估濃度，Xfflin,r>Xfflax,r>kr, Xmin,s、Xmax,s、ks為 擬合參數(shù)，sr, Ss分別為內(nèi)部熱耦合精餾塔精餾段、提餾段液相組分濃度分布的位置；自適應非線性控制律求解模塊，用以根據(jù)當前組分濃度數(shù)據(jù)，參考軌跡和當前時 刻操作變量值求取當前的控制變量的理想改變值，采用式(5)-(14)得到 其中，k為當前采樣時刻，t為采樣周期，Xi (k)、Yi (k)分別為k采樣時刻第i塊塔 板輕組分液相輕組分濃度和汽相輕組分濃度，QiGO為第i塊塔板之間的熱耦合量，UA為傳 熱速率，Xi^ (K)為k采樣時刻第i+f-Ι塊塔板液相輕組分濃度，q(k)為k采樣時刻進料熱狀況，PJk)為當k采樣時刻精餾段壓強，F(xiàn)為進料流率，Zf為進料組分濃度，V1GO^fGO分 別為k采樣時刻塔頂和進料板的汽相流率，L1 (k)、Lf_i (k)、Ln(k)分別為k采樣時刻塔頂、第 f-Ι塊塔板和塔底的液相流率，H為持液量，λ為汽化潛熱，Y1*, Xn*分別為塔頂?shù)钠噍p組 分濃度Y1、塔底的液相輕組分濃度Xn的設定值，Xj (k)、Xn(k)、Xf^1 (k)分別為k采樣時刻第 j塊塔板、塔底、第f-Ι塊塔板的液相輕組分濃度，Y1 (k)、Yf (k)分別為k采樣時刻塔頂和進 料板的汽相輕組分濃度，K1, K2, K3, K4為控制律參數(shù)，S；, S；分別為精餾段提餾段拐點參考 軌跡，SJk) ,Sr(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔精餾段液相組分濃度分布的位 置，Ss(k)、Ss(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔提餾段液相組分濃度分布的位置， Aq(k), APJk)，分別為當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強 的當前理想改變值。作為優(yōu)選的一種方案所述的上位機還包括人機界面模塊，用以設定采樣周期t， 控制律參數(shù)K1, K2, K3, K4和塔頂?shù)钠噍p組分濃度的設定值ΥΛ塔底的液相輕組分濃度的 設定值χΛ并顯示控制器的輸出曲線和被控變量即塔頂塔底液相輕組分濃度的記錄曲線。一種內(nèi)部熱耦合精餾塔的自適應非線性控制方法，所述的控制方法包括以下步 驟1)確定采樣周期t，并將t值，相對揮發(fā)度α，提餾段壓強Ps，安東尼常數(shù)a、b、C、 保存在歷史數(shù)據(jù)庫當中；2)設定控制律參數(shù)K1, K2, K3, K4和塔頂?shù)钠噍p組分濃度的設定值ΥΛ塔底的液 相輕組分濃度的設定值Xn*;3)從現(xiàn)場智能儀表獲取k采樣時刻精餾段壓強&提餾段壓強Ps，以及各塔板溫度 Ti,計算液相輕組分濃度值，采用式(1) (2)得到式 其中，k為當前采樣時刻，下腳標i為塔板編號，1為塔頂編號，f為進料板編號，η 為塔底編號，Xi (k)為k采樣時刻液相輕組分濃度，PJk)為k采樣時刻精餾段壓強、PS提餾 段壓強,TiGO為k采樣時刻塔內(nèi)各塊塔板的溫度，α為相對揮發(fā)度，a、b、c為安東尼常數(shù)；4)用歷史數(shù)據(jù)庫中組分推斷模塊計算出的組分濃度數(shù)據(jù)，在線校正模塊實現(xiàn)模型 參數(shù)，擬合模型函數(shù)如下(3)、(4)X1 = Xmmr1=1,2,……fl(3) 其中，友為第i塊塔板處液相組分濃度預估濃度，Xfflinjr, Xfflaxjr, kr, Xfflinjs, Xfflaxjs, ks為 擬合參數(shù)，sr, Ss分別為內(nèi)部熱耦合精餾塔精餾段，提餾段液相組分濃度分布的位置；5)根據(jù)當前組分濃度數(shù)據(jù)，模型函數(shù)和當前時刻操作變量值求取當前的控制變量 的理想改變值，采用式(5)-(14)得到 其中，k為當前采樣時刻，t為采樣周期，Xi (k)、Ti (k)分別為k采樣時刻第i塊塔 板輕組分液相輕組分濃度和汽相輕組分濃度，QiGO為第i塊塔板之間的熱耦合量，UA為傳 熱速率，Xi^GO為k采樣時刻第i+f-l塊塔板液相輕組分濃度，q(k)為k采樣時刻進料熱 狀況，PJk)為當k采樣時刻精餾段壓強，F(xiàn)為進料流率，Zf為進料組分濃度，V1GO^fGO分 別為k采樣時刻塔頂和進料板的汽相流率，L1 (k)、Lf_i (k)、Ln(k)分別為k采樣時刻塔頂、第 f-Ι塊塔板和塔底的液相流率，，H為持液量，λ為汽化潛熱，Y廣X1；分別為塔頂?shù)钠噍p組 分濃度Y1、塔底的液相輕組分濃度Xn的設定值，Xj (k)、Xn(k)、Xf^1 (k)分別為k采樣時刻第 j塊塔板、塔底、第f-Ι塊塔板的液相輕組分濃度，Y1 (k)、Yf (k)分別為k采樣時刻塔頂和進 料板的汽相輕組分濃度，K1, K2, K3, K4為控制律參數(shù)，S；, S；分別為精餾段提餾段拐點參考 軌跡，SJk) ,Sr(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔精餾段液相組分濃度分布的位 置，Ss(k)、Ss(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔提餾段液相組分濃度分布的位置， Aq(k), APJk)，分別為當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強 的當前理想改變值；6)將當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強的當前理想改變值Aq(k)，APr(k)輸送給DCS系統(tǒng)中的控制站，調(diào)整內(nèi)部熱耦合精餾塔的進料熱狀 況值和精餾段壓強值。進一步，所述歷史數(shù)據(jù)庫為DCS系統(tǒng)中的存儲裝置，控制站讀取歷史數(shù)據(jù)庫，顯示 內(nèi)部熱耦合精餾塔工作過程狀態(tài)。本發(fā)明的有益效果主要表現(xiàn)在1.所述的非線性控制方案建立在高精度非線性 模型基礎上，能夠及時抑制干擾作用；2.控制方案較好地處理了耦合問題，能夠快速準確 地跟蹤設定值變化。
圖1是本發(fā)明所提出的內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)的結(jié)構(gòu)圖。圖2是上位機控制器實現(xiàn)方法的原理圖。
具體實施例方式下面根據(jù)附圖具體說明本發(fā)明。實施例1參照圖1，圖2，一種內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)，控制系統(tǒng)包括 與內(nèi)部熱耦合精餾塔直接連接的現(xiàn)場智能儀表2和DCS系統(tǒng)，所述DCS系統(tǒng)包括存儲裝置 4、控制站5和上位機6，其中現(xiàn)場智能儀表2、存儲裝置4、控制站5和上位機6通過現(xiàn)場總 線和數(shù)據(jù)接口 3依次相連；所述的智能儀表用以實現(xiàn)檢測模塊7的功能即檢測得到內(nèi)部熱 耦合精餾塔的溫度、壓強數(shù)據(jù)，所述的數(shù)據(jù)接口用以實現(xiàn)I/O模塊8的功能即數(shù)據(jù)的輸入輸 出ο所述的上位機包括用以計算輸出內(nèi)部熱耦合精餾塔控制變量值的自適應非線性 控制器，所述自適應非線性控制器包括組分推斷模塊9、模型參數(shù)自適應校正模塊10和自 適應非線性控制律求解模塊11 ；組分推斷模塊9，用以從現(xiàn)場智能儀表獲取溫度，壓強數(shù)據(jù)，計算內(nèi)部熱耦合精餾 塔各塊塔板的組分濃度，并將組分濃度計算結(jié)果存儲在歷史數(shù)據(jù)庫當中，采用式(1) (2)得 到 其中，k為當前采樣時刻，下腳標i為塔板編號，1為塔頂編號，f為進料板編號， η為塔底編號，Xi (k)為k采樣時刻液相輕組分濃度，PJk)為k采樣時刻精餾段壓強、Ps提 餾段壓強，TiGO為k采樣時刻塔內(nèi)各塊塔板的溫度，α為相對揮發(fā)度，a、b、c為安東尼 (Antonie)常數(shù)。模型參數(shù)自適應校正模塊10，用以實現(xiàn)模型參數(shù)的在線更新，在線擬合模型函數(shù) 如下
(3)
其中友為第i塊塔板處液相組分濃度預估濃度，Xfflinjr, Xfflaxjr, kr, Xfflinjs, Xfflaxjs, ks為 擬合參數(shù)，Sr, Ss分別為內(nèi)部熱耦合精餾塔精餾段，提餾段液相組分濃度分布的位置；自適應非線性控制律求解模塊11，用以根據(jù)當前組分濃度數(shù)據(jù)，參考軌跡和當前 時刻操作變量值求取當前的控制變量的理想改變值，采用式(5)-(14)得到 其中，k為當前采樣時刻，t為采樣周期，Xi (k)、Yi (k)分別為k采樣時刻第i塊塔 板輕組分液相輕組分濃度和汽相輕組分濃度，QiGO為第i塊塔板之間的熱耦合量，UA為傳 熱速率，Xi^GO為k采樣時刻第i+f-l塊塔板液相輕組分濃度，q(k)為k采樣時刻進料熱 狀況，PJk)為當k采樣時刻精餾段壓強，F(xiàn)為進料流率，Zf為進料組分濃度，V1GO^fGO分 別為k采樣時刻塔頂和進料板的汽相流率，L1 (k)、Lf_i (k)、Ln(k)分別為k采樣時刻塔頂、第 f-Ι塊塔板和塔底的液相流率，，H為持液量，λ為汽化潛熱，Y廣X1；分別為塔頂?shù)钠噍p組 分濃度Y1、塔底的液相輕組分濃度Xn的設定值，X^K)、Xn(k)、Xf_i (k)分別為k采樣時刻第j 塊塔板、塔底、第f-ι塊塔板的液相輕組分濃度，Y1GO、Yf(k)分別為k采樣時刻塔頂和進料板的汽相輕組分濃度，K1, K2, K3, K4為控制律參數(shù)，K1, K3 e [2,200]，K2, K4 e [20，2000]，根 據(jù)具體操作對象特性具體調(diào)節(jié)，S；, S；分別為精餾段提餾段拐點參考軌跡，Sr (k)、Sr (i)分 別為k和i采樣時刻內(nèi)部熱耦合精餾塔精餾段液相組分濃度分布的位置，Ss (k)、Ss(i)分別 為k和i采樣時刻內(nèi)部熱耦合精餾塔提餾段液相組分濃度分布的位置，Aq(k)，APJk)，分 別為當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強的當前理想改變值。人機界面模塊12，用于設定采樣周期T，控制律參數(shù)K1, K2, K3, K4和塔頂?shù)钠噍p 組分濃度的設定值Y久塔底的液相輕組分濃度的設定值ΧΛ并顯示控制器的輸出曲線和被 控變量即塔頂塔底液相輕組分濃度的記錄曲線。實施例2參照圖1和圖2，一種內(nèi)部熱耦合精餾塔的自適應非線性控制方法，所述的控制方 法包括以下步驟1)確定采樣周期t，并將t值，相對揮發(fā)度α，提餾段壓強Ps，安東尼(Antonie)常 數(shù)a、b、C、保存在歷史數(shù)據(jù)庫當中；2)設定控制律參數(shù)K1, K2, K3, K4和塔頂?shù)钠噍p組分濃度的設定值ΥΛ塔底的液 相輕組分濃度的設定值Xn*;3)從現(xiàn)場智能儀表獲取k采樣時刻精餾段壓強&提餾段壓強Ps，以及各塔板溫度 Ti,計算液相輕組分濃度值，采用式(1) (2)得到 其中，k為當前采樣時刻，下腳標i為塔板編號，1為塔頂編號，f為進料板編號，η 為塔底編號，Xi (k)為k采樣時刻液相輕組分濃度，PJk)為k采樣時刻精餾段壓強、PS提餾 段壓強,TiGO為k采樣時刻塔內(nèi)各塊塔板的溫度，α為相對揮發(fā)度，a、b、c為安東尼常數(shù)；4)用歷史數(shù)據(jù)庫中組分濃度實現(xiàn)模型參數(shù)的在線更新，在線擬合模型函數(shù)如下 其中，友為第i塊塔板處液相組分濃度預估濃度，Xfflinjr, Xfflaxjr, kr, Xfflinjs, Xfflaxjs, ks為 擬合參數(shù)，sr, Ss分別為內(nèi)部熱耦合精餾塔精餾段，提餾段液相組分濃度分布的位置；5)根據(jù)當前組分濃度數(shù)據(jù)，模型函數(shù)和當前時刻操作變量值求取當前的控制變量 的理想改變值，采用式(5)-(14)得到 其中，k為當前采樣時刻，t為采樣周期，Xi (k)、Yi (k)分別為k采樣時刻第i塊塔 板輕組分液相輕組分濃度和汽相輕組分濃度，QiGO為第i塊塔板之間的熱耦合量，UA為傳 熱速率，Xi^GO為k采樣時刻第i+f-l塊塔板液相輕組分濃度，q(k)為k采樣時刻進料熱 狀況，PJk)為當k采樣時刻精餾段壓強，F(xiàn)為進料流率，Zf為進料組分濃度，V1GO^fGO分 別為k采樣時刻塔頂和進料板的汽相流率，L1 (k)、Lf_i (k)、Ln(k)分別為k采樣時刻塔頂、第 f-Ι塊塔板和塔底的液相流率，H為持液量，λ為汽化潛熱，Y1*, Xn*分別為塔頂?shù)钠噍p組 分濃度Y1、塔底的液相輕組分濃度Xn的設定值，Xj (k)、Xn(K)、Xf^1 (k)分別為k采樣時刻第 j塊塔板、塔底、第f-Ι塊塔板的液相輕組分濃度，Y1 (k)、Yf (k)分別為k采樣時刻塔頂和進 料板的汽相輕組分濃度，K1, K2, K3, K4為控制律參數(shù)，K1, K3 e [2，200]，K2, K4 e [20，2000]， 根據(jù)具體操作對象特性具體調(diào)節(jié)，S；, S/分別為精餾段提餾段拐點參考軌跡，Sr(k),Sr(i) 分別為k和i采樣時刻內(nèi)部熱耦合精餾塔精餾段液相組分濃度分布的位置，Ss (k)、Ss(i)分 別為k和i采樣時刻內(nèi)部熱耦合精餾塔提餾段液相組分濃度分布的位置，Aq(k)，APr(k)， 分別為當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強的當前理想改變 值；6)將當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強的當前理 想改變值Aq(k)，APr(k)輸送給DCS系統(tǒng)中的控制站，調(diào)整內(nèi)部熱耦合精餾塔的進料熱狀 況值和精餾段壓強值。所述的歷史數(shù)據(jù)庫為DCS系統(tǒng)中的存儲裝置4，所述的DCS系統(tǒng)包括數(shù)據(jù)接口 3、 存儲裝置4和控制站5，其中控制站5可以讀取歷史數(shù)據(jù)庫，顯示內(nèi)部熱耦合精餾塔工作過程狀態(tài)。
上述實施例用來解釋說明本發(fā)明，而不是對本發(fā)明進行限制，在本發(fā)明的精神和 權(quán)利要求的保護范圍內(nèi)，對本發(fā)明作出的任何修改和改變，都落入本發(fā)明的保護范圍。
權(quán)利要求
一種內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)，包括與內(nèi)部熱耦合精餾塔直接連接的現(xiàn)場智能儀表和DCS系統(tǒng)，所述DCS系統(tǒng)包括存儲裝置、控制站和上位機，所述現(xiàn)場智能儀表與存儲裝置、控制站和上位機相連，其特征在于所述的上位機包括用以計算輸出內(nèi)部熱耦合精餾塔控制變量值的自適應非線性控制器，所述自適應非線性控制器包括組分推斷模塊，用以從現(xiàn)場智能儀表獲取溫度，壓強數(shù)據(jù)，計算內(nèi)部熱耦合精餾塔各塊塔板的組分濃度，并將組分濃度計算結(jié)果存儲在歷史數(shù)據(jù)庫當中，采用式(1)(2)得到 <mrow><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>P</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>×</mo><mi>α</mi><mo>×</mo><msup> <mn>10</mn> <mrow><mfrac> <mrow><msub> <mi>T</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mi>c</mi> </mrow> <mi>b</mi></mfrac><mo>-</mo><mi>a</mi> </mrow></msup><mo>-</mo><mn>1</mn> </mrow> <mrow><mi>α</mi><mo>-</mo><mn>1</mn> </mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>f</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>P</mi> <mi>s</mi></msub><mo>×</mo><mi>α</mi><mo>×</mo><msup> <mn>10</mn> <mrow><mfrac> <mrow><msub> <mi>T</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mi>c</mi> </mrow> <mi>b</mi></mfrac><mo>-</mo><mi>a</mi> </mrow></msup><mo>-</mo><mn>1</mn> </mrow> <mrow><mi>α</mi><mo>-</mo><mn>1</mn> </mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mi>f</mi><mo>,</mo><mi>f</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow> </mrow>其中，k為當前采樣時刻，下腳標i為塔板編號，1為塔頂編號，f為進料板編號，n為塔底編號，Xi(k)為k采樣時刻液相輕組分濃度，Pr(k)為k采樣時刻精餾段壓強、Ps提餾段壓強，Ti(k)為k采樣時刻塔內(nèi)各塊塔板的溫度，α為相對揮發(fā)度，a、b、c為安東尼常數(shù)；模型參數(shù)自適應校正模塊，用以實現(xiàn)模型參數(shù)的在線更新，在線擬合模型函數(shù)如下(3)、(4) <mrow><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mi>i</mi></msub><mo>=</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>-</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub> </mrow> <mrow><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <mrow><mo>-</mo><msub> <mi>k</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub><mi>S</mi><mi>r</mi> </msub> <mo>)</mo></mrow> </mrow></msup> </mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>f</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mover><mi>X</mi><mo>^</mo> </mover> <mi>i</mi></msub><mo>=</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>+</mo><mfrac> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>-</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub> </mrow> <mrow><mn>1</mn><mo>+</mo><msup> <mi>e</mi> <mrow><mo>-</mo><msub> <mi>k</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub><mi>S</mi><mi>s</mi> </msub> <mo>)</mo></mrow> </mrow></msup> </mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mi>f</mi><mo>,</mo><mi>f</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo></mrow> </mrow>其中，為第i塊塔板處液相組分濃度預估濃度，Xmin，r、Xmax，r、kr，Xmin，s、Xmax，s、ks為擬合參數(shù)，Sr、Ss分別為內(nèi)部熱耦合精餾塔精餾段、提餾段液相組分濃度分布的位置；自適應非線性控制律求解模塊，用以根據(jù)當前組分濃度數(shù)據(jù)，參考軌跡和當前時刻操作變量值求取當前的控制變量的理想改變值，采用式(5)-(14)得到 <mrow><msub> <mi>Y</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>αX</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mrow> <mo>(</mo> <mi>α</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mn>1</mn> </mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>Q</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><mi>UA</mi><mo>×</mo><mi>b</mi><mrow> <mo>(</mo> <mfrac><mn>1</mn><mrow> <mi>a</mi> <mo>-</mo> <mi>ln</mi> <mo>{</mo> <mrow><mo>(</mo><msub> <mi>P</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mi>ΔPr</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>×</mo><mi>t</mi><mo>)</mo> </mrow> <mo>/</mo> <mo>[</mo> <msub><mi>X</mi><mi>i</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>+</mo> <mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub> <mi>X</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mo>/</mo> <mi>α</mi> <mo>]</mo> <mo>}</mo></mrow> </mfrac></mrow> </mrow>i＝1，2，……，f-1(6) <mrow><mo>-</mo><mfrac> <mn>1</mn> <mrow><mi>a</mi><mo>-</mo><mi>ln</mi><mo>{</mo><msub> <mi>p</mi> <mi>S</mi></msub><mo>/</mo><mo>[</mo><msub> <mi>X</mi> <mrow><mi>i</mi><mo>+</mo><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub><mi>X</mi><mrow> <mi>i</mi> <mo>+</mo> <mi>f</mi> <mo>-</mo> <mn>1</mn></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>/</mo><mi>α</mi><mo>]</mo><mo>}</mo> </mrow></mfrac><mo>)</mo> </mrow>V1(k)＝F(1-q(k)-Δq(k)×t)(7)Ln(k)＝F(q(k)+Δq(k)×t) (8) <mrow><msub> <mi>L</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mfrac> <mrow><msub> <mi>Q</mi> <mi>i</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mi>λ</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo></mrow> </mrow>Vf(k)＝V1(k)+Lf-1(k) (10) <mrow><msubsup> <mi>S</mi> <mi>r</mi> <mo>*</mo></msubsup><mo>=</mo><mn>1</mn><mo>+</mo><mfrac> <mn>1</mn> <msub><mi>k</mi><mi>r</mi> </msub></mfrac><mi>ln</mi><mrow> <mo>(</mo> <mfrac><mrow> <msub><mi>X</mi><mrow> <mi>max</mi> <mo>,</mo> <mi>r</mi></mrow> </msub> <mo>-</mo> <msubsup><mi>Y</mi><mn>1</mn><mo>*</mo> </msubsup> <mo>/</mo> <mo>[</mo> <mi>α</mi> <mo>-</mo> <mrow><mo>(</mo><mi>α</mi><mo>-</mo><mn>1</mn><mo>)</mo> </mrow> <msubsup><mi>Y</mi><mn>1</mn><mo>*</mo> </msubsup> <mo>]</mo></mrow><mrow> <msubsup><mi>Y</mi><mn>1</mn><mo>*</mo> </msubsup> <mo>/</mo> <mo>[</mo> <mi>α</mi> <mo>-</mo> <mrow><mo>(</mo><mi>α</mi><mo>-</mo><mn>1</mn><mo>)</mo> </mrow> <msubsup><mi>Y</mi><mn>1</mn><mo>*</mo> </msubsup> <mo>]</mo> <mo>-</mo> <msub><mi>X</mi><mrow> <mi>min</mi> <mo>,</mo> <mi>r</mi></mrow> </msub></mrow> </mfrac> <mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>f</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo></mrow> </mrow> <mrow><msubsup> <mi>S</mi> <mi>s</mi> <mo>*</mo></msubsup><mo>=</mo><mi>n</mi><mo>+</mo><mfrac> <mn>1</mn> <msub><mi>k</mi><mi>s</mi> </msub></mfrac><mi>ln</mi><mrow> <mo>(</mo> <mfrac><mrow> <msub><mi>X</mi><mrow> <mi>max</mi> <mo>,</mo> <mi>s</mi></mrow> </msub> <mo>-</mo> <msubsup><mi>X</mi><mi>n</mi><mo>*</mo> </msubsup></mrow><mrow> <msubsup><mi>X</mi><mi>n</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>X</mi><mrow> <mi>min</mi> <mo>,</mo> <mi>s</mi></mrow> </msub></mrow> </mfrac> <mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mi>f</mi><mo>,</mo><mi>f</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo></mrow> </mrow> <mrow><mfrac> <mrow><mo>-</mo><msub> <mi>V</mi> <mn>1</mn></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>Y</mi> <mn>1</mn></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>V</mi> <mi>f</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>Y</mi> <mi>f</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>L</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>X</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mi>H</mi><munderover> <mi>Σ</mi> <mrow><mi>j</mi><mo>=</mo><mn>1</mn> </mrow> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mfrac> <mrow><msub> <mi>k</mi> <mi>r</mi></msub><mrow> <mo>(</mo> <msub><mi>X</mi><mrow> <mi>min</mi> <mo>,</mo> <mi>r</mi></mrow> </msub> <mo>-</mo> <msub><mi>X</mi><mi>j</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mrow> <mo>(</mo> <msub><mi>X</mi><mi>j</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>X</mi><mrow> <mi>max</mi> <mo>,</mo> <mi>r</mi></mrow> </msub> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>r</mi> </mrow></msub><mo>-</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>r</mi> </mrow></msub> </mrow></mfrac> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo></mrow> </mrow> <mrow><mo>=</mo><msub> <mi>K</mi> <mn>1</mn></msub><mrow> <mo>(</mo> <msubsup><mi>S</mi><mi>r</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>S</mi><mi>r</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>+</mo><msub> <mi>K</mi> <mn>2</mn></msub><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>k</mi></munderover><mrow> <mo>(</mo> <msubsup><mi>S</mi><mi>r</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>S</mi><mi>r</mi> </msub> <mrow><mo>(</mo><mi>i</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mi>t</mi> </mrow> <mrow><mfrac> <mrow><msub> <mrow><mo>-</mo><mi>V</mi> </mrow> <mi>f</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>Y</mi> <mi>f</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>L</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>X</mi> <mrow><mi>f</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>FZ</mi> <mi>f</mi></msub><mo>-</mo><msub> <mi>L</mi> <mi>n</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>X</mi> <mi>n</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><mi>H</mi><munderover> <mi>Σ</mi> <mrow><mi>j</mi><mo>=</mo><mi>f</mi> </mrow> <mi>n</mi></munderover><mfrac> <mrow><msub> <mi>k</mi> <mi>s</mi></msub><mrow> <mo>(</mo> <msub><mi>X</mi><mrow> <mi>min</mi> <mo>,</mo> <mi>s</mi></mrow> </msub> <mo>-</mo> <msub><mi>X</mi><mi>j</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mrow> <mo>(</mo> <msub><mi>X</mi><mi>j</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>X</mi><mrow> <mi>max</mi> <mo>,</mo> <mi>s</mi></mrow> </msub> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>X</mi> <mrow><mi>max</mi><mo>,</mo><mi>s</mi> </mrow></msub><mo>-</mo><msub> <mi>X</mi> <mrow><mi>min</mi><mo>,</mo><mi>s</mi> </mrow></msub> </mrow></mfrac> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo></mrow> </mrow> <mrow><mo>=</mo><msub> <mi>K</mi> <mn>3</mn></msub><mrow> <mo>(</mo> <msubsup><mi>S</mi><mi>s</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>S</mi><mi>s</mi> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>+</mo><msub> <mi>K</mi> <mn>4</mn></msub><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>k</mi></munderover><mrow> <mo>(</mo> <msubsup><mi>S</mi><mi>s</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>S</mi><mi>s</mi> </msub> <mrow><mo>(</mo><mi>i</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mi>t</mi> </mrow>其中，k為當前采樣時刻，t為采樣周期，Xi(k)、Yi(k)分別為k采樣時刻第i塊塔板輕組分液相輕組分濃度和汽相輕組分濃度，Qi(k)為第i塊塔板之間的熱耦合量，UA為傳熱速率，Xi+f-1(k)為k采樣時刻第i+f-1塊塔板液相輕組分濃度，q(k)為k采樣時刻進料熱狀況，Pr(k)為當k采樣時刻精餾段壓強，F(xiàn)為進料流率，Zf為進料組分濃度，V1(k)、Vf(k)分別為k采樣時刻塔頂和進料板的汽相流率，L1(k)、Lf-1(k)、Ln(k)分別為k采樣時刻塔頂、第f-1塊塔板和塔底的液相流率，H為持液量，λ為汽化潛熱，Y1*、Xn*分別為塔頂?shù)钠噍p組分濃度Y1、塔底的液相輕組分濃度Xn的設定值，Xj(k)、Xn(k)、Xf-1(k)分別為k采樣時刻第j塊塔板、塔底、第f-1塊塔板的液相輕組分濃度，Y1(k)、Yf(k)分別為k采樣時刻塔頂和進料板的汽相輕組分濃度，K1，K2，K3，K4為控制律參數(shù)，Sr*，Ss*分別為精餾段提餾段拐點參考軌跡，Sr(k)、Sr(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔精餾段液相組分濃度分布的位置，Ss(k)、Ss(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔提餾段液相組分濃度分布的位置，Δq(k)，ΔPr(k)，分別為當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強的當前理想改變值。FDA0000022884190000015.tif
2.如權(quán)利要求1所述的內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)，其特征在于所 述的上位機還包括人機界面模塊，用以設定采樣周期t，控制律參數(shù)K1, K2, K3, K4和塔頂?shù)?汽相輕組分濃度的設定值ΥΛ塔底的液相輕組分濃度的設定值X:，并顯示控制器的輸出曲 線和被控變量即塔頂塔底液相輕組分濃度的記錄曲線。
3.一種用如權(quán)利要求1所述的內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)實現(xiàn)的自 適應非線性控制方法，其特征在于所述的控制方法包括以下步驟1)確定采樣周期t，并將t值，相對揮發(fā)度α，提餾段壓強Ps，安東尼常數(shù)a、b、C、保存 在歷史數(shù)據(jù)庫當中；2)設定控制律參數(shù)K1,K2, K3, K4和塔頂?shù)钠噍p組分濃度的設定值ΥΛ塔底的液相輕 組分濃度的設定值Xn*;3)從現(xiàn)場智能儀表獲取k采樣時刻精餾段壓強&提餾段壓強Ps，以及各塔板溫度Ti, 計算液相輕組分濃度值，采用式(1) (2)得到其中，k為當前采樣時刻，下腳標i為塔板編號，1為塔頂編號，f為進料板編號，η為塔 底編號，Xi (k)為k采樣時刻液相輕組分濃度，PJk)為k采樣時刻精餾段壓強、Ps提餾段壓 強，YiGO為k采樣時刻塔內(nèi)各塊塔板的溫度，α為相對揮發(fā)度，a、b、c為安東尼常數(shù)；4)用歷史數(shù)據(jù)庫中組分推斷模塊計算出的組分濃度數(shù)據(jù)，在線校正模塊實現(xiàn)模型參 數(shù)，擬合模型函數(shù)如下(3)、(4) 其中，友為第i塊塔板處液相組分濃度預估濃度，Xfflinjr, Xfflaxjr, kr, Xfflinjs, Xfflaxjs, ks為擬合 參數(shù)，Sr, Ss分別為內(nèi)部熱耦合精餾塔精餾段，提餾段液相組分濃度分布的位置；5)根據(jù)當前組分濃度數(shù)據(jù)，模型函數(shù)和當前時刻操作變量值求取當前的控制變量的理 想改變值，采用式(5)-(14)得到 其中，k為當前采樣時刻，t為采樣周期，Xi (k)、Yi (k)分別為k采樣時刻第i塊塔板輕 組分液相輕組分濃度和汽相輕組分濃度，QiGO為第i塊塔板之間的熱耦合量，UA為傳熱速 率，Xi+f-!(k)為k采樣時刻第i+f-Ι塊塔板液相輕組分濃度，q(k)為k采樣時刻進料熱狀 況，Pr (k)為當k采樣時刻精餾段壓強，F(xiàn)為進料流率，Zf為進料組分濃度，Vjk)、Vf(k)分 別為k采樣時刻塔頂和進料板的汽相流率，L1 (k)、Lf_i (k)、Ln(k)分別為k采樣時刻塔頂、第 f-Ι塊塔板和塔底的液相流率，，H為持液量，λ為汽化潛熱，Y廣X1；分別為塔頂?shù)钠噍p組 分濃度Y1、塔底的液相輕組分濃度Xn的設定值，Xj (k)、Xn(K)、Xf^1 (k)分別為k采樣時刻第 j塊塔板、塔底、第f-Ι塊塔板的液相輕組分濃度，Y1 (k)、Yf (k)分別為k采樣時刻塔頂和進 料板的汽相輕組分濃度，K1, k2，K3, K4為控制律參數(shù)，S；, S；分別為精餾段提餾段拐點參考 軌跡，SJk) ,Sr(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔精餾段液相組分濃度分布的位 置，Ss(k)、Ss(i)分別為k和i采樣時刻內(nèi)部熱耦合精餾塔提餾段液相組分濃度分布的位置， Aq(k), APJk)，分別為當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強 的當前理想改變值；6)將當前時刻內(nèi)部熱耦合精餾塔控制變量即進料熱狀況和精餾段壓強的當前理想改 變值Aq(k)，APr(k)輸送給DCS系統(tǒng)中的控制站，調(diào)整內(nèi)部熱耦合精餾塔的進料熱狀況值 和精餾段壓強值。
4.如權(quán)利要求3所述的自適應非線性控制方法，其特征在于所述歷史數(shù)據(jù)庫為DCS 系統(tǒng)中的存儲裝置，控制站讀取歷史數(shù)據(jù)庫，顯示內(nèi)部熱耦合精餾塔工作過程狀態(tài)。
全文摘要
一種內(nèi)部熱耦合精餾塔的自適應非線性控制系統(tǒng)，包括與內(nèi)部熱耦合精餾塔直接連接的現(xiàn)場智能儀表和DCS系統(tǒng)，所述DCS系統(tǒng)包括存儲裝置、控制器和上位機，所述現(xiàn)場智能儀表與存儲裝置、控制站和上位機相連，所述的上位機包括用以計算輸出內(nèi)部熱耦合精餾塔控制變量值的自適應非線性控制器，所述自適應非線性控制器包括組分推斷模塊，模型參數(shù)自適應校正模塊，自適應非線性控制律求解模塊。同時提供了一種內(nèi)部熱耦合精餾塔的自適應非線性控制方法。本發(fā)明能夠很好地處理內(nèi)部熱耦合節(jié)能精餾的強非線性特征，具有高效的在線運算速度，并具有非常好的伺服跟蹤控制效果，干擾抑制效果良好。
文檔編號G05B19/418GK101887261SQ20101021426
公開日2010年11月17日 申請日期2010年6月30日 優(yōu)先權(quán)日2010年6月30日
發(fā)明者劉興高, 周葉翔 申請人:浙江大學