葉凌箭，宋執(zhí)環(huán)，馬修水

(1 浙江大學(xué)寧波理工學(xué)院, 浙江 寧波 315100；2 浙江大學(xué)控制科學(xué)與工程學(xué)系, 浙江 杭州 310027)

引 言

間歇過(guò)程的優(yōu)化是一個(gè)動(dòng)態(tài)優(yōu)化問(wèn)題，數(shù)學(xué)上求解相對(duì)于穩(wěn)態(tài)過(guò)程難度更大。解析法如應(yīng)用Pontryagin 極值原理最小化Hamiltonian 函數(shù)等，對(duì)于反應(yīng)機(jī)理復(fù)雜、規(guī)模龐大的間歇過(guò)程幾乎很難實(shí)現(xiàn)。應(yīng)用數(shù)值法求解間歇過(guò)程的優(yōu)化問(wèn)題更為常見(jiàn)，如控制向量參數(shù)化法（control vector parameterization，CVP）等[1]。

對(duì)間歇生產(chǎn)過(guò)程進(jìn)行操作優(yōu)化時(shí)，如果僅在標(biāo)稱狀態(tài)下進(jìn)行一次離線優(yōu)化，優(yōu)化性能具有局限性。主要原因是工業(yè)過(guò)程在實(shí)際運(yùn)行時(shí)受到各種不確定擾動(dòng)的影響，引起了系統(tǒng)工作點(diǎn)的漂移。標(biāo)稱工作點(diǎn)的最優(yōu)輸入軌跡無(wú)法在真實(shí)工況下保持最優(yōu)性，引起了目標(biāo)函數(shù)損失。實(shí)時(shí)優(yōu)化（real-time optimization，RTO）的任務(wù)即在不確定擾動(dòng)影響下使過(guò)程系統(tǒng)重新運(yùn)行于新工況下的最優(yōu)點(diǎn)[2]。

根據(jù)RTO 策略和實(shí)現(xiàn)方法不同，Chachuat 等[3]將現(xiàn)有的RTO 方法主要分為3 大類：（1）參數(shù)估計(jì)-重優(yōu)化兩步法[4]；（2）修正適應(yīng)法（modifier adaptation）[5]；（3）直接調(diào)節(jié)輸入法。方法（1）利用過(guò)程的測(cè)量變量數(shù)據(jù)，首先基于參數(shù)化模型估計(jì)出不確定擾動(dòng)，然后重新優(yōu)化求解。方法（2）不對(duì)模型進(jìn)行更新，而是對(duì)約束條件和目標(biāo)函數(shù)進(jìn)行合理修正，使得修正后的最優(yōu)性條件（necessary conditions of optimality，NCO）和真實(shí)過(guò)程的NCO相匹配，從而找到最優(yōu)輸入。方法（3）以NCO 跟蹤法[2,6]和自優(yōu)化控制（self-optimizing control，SOC）[7]為代表，NCO 跟蹤的主要特征是測(cè)量/估計(jì)出當(dāng)前的NCO，設(shè)計(jì)控制律使得系統(tǒng)最終滿足NCO。SOC則通過(guò)跟蹤控制某些經(jīng)過(guò)離線設(shè)計(jì)的被控變量，即使不在線調(diào)整設(shè)定值，系統(tǒng)也在反饋控制器作用下自動(dòng)運(yùn)行于最優(yōu)點(diǎn)附近。以上優(yōu)化策略中，大多需要在線測(cè)量/估計(jì)目標(biāo)函數(shù)和測(cè)量變量的梯度信息，或者進(jìn)行頻繁的高負(fù)荷數(shù)學(xué)運(yùn)算。相比而言，SOC的在線實(shí)現(xiàn)最為簡(jiǎn)便，控制結(jié)構(gòu)中無(wú)須單獨(dú)的優(yōu)化層，而是通過(guò)常規(guī)控制作用實(shí)現(xiàn)RTO 效果。由于SOC 設(shè)計(jì)的被控變量（測(cè)量變量的組合）間接地獲取了梯度信息，因此無(wú)須像有限差分法對(duì)系統(tǒng)進(jìn)行激勵(lì)以獲取實(shí)驗(yàn)梯度，從而極大地加快了RTO 速度。

由于具備良好RTO 性能的同時(shí)保持了簡(jiǎn)潔的可實(shí)現(xiàn)性，SOC 自Skogestad[7]于2000年提出后就引起了相關(guān)學(xué)者的關(guān)注和研究[8-13]，大量平臺(tái)上的成功應(yīng)用不斷報(bào)道?，F(xiàn)有的SOC 主要應(yīng)用于連續(xù)過(guò)程，在間歇過(guò)程上的應(yīng)用較少。文獻(xiàn)[14-16]基于SOC 思想對(duì)間歇過(guò)程進(jìn)行了一些初步的工作，但具有如下兩個(gè)方面的局限性：（1）滿足最優(yōu)性條件的Hamiltonian 函數(shù)經(jīng)解析法推導(dǎo)，因此只適用于一些機(jī)理簡(jiǎn)單、規(guī)模微小的間歇過(guò)程；（2）要求實(shí)時(shí)測(cè)量一些關(guān)鍵的過(guò)程變量如物料濃度等，實(shí)際操作時(shí)通常很難滿足這一點(diǎn)。如果過(guò)程變量只在批次運(yùn)行結(jié)束后才被測(cè)量，只能通過(guò)批間方式實(shí)現(xiàn)RTO。目前，基于SOC 策略的間歇過(guò)程批間優(yōu)化，仍未見(jiàn) 報(bào)道。

1 連續(xù)過(guò)程自優(yōu)化控制

自優(yōu)化控制的核心思想是對(duì)系統(tǒng)的控制結(jié)構(gòu)進(jìn)行干預(yù)，通過(guò)選擇（設(shè)計(jì)）合理的被控變量實(shí)現(xiàn)RTO。擾動(dòng)產(chǎn)生時(shí)，無(wú)須進(jìn)行重新優(yōu)化求解更新設(shè)定值，而是只在常規(guī)控制器下跟蹤被控變量的恒定設(shè)定值，就能夠使系統(tǒng)自動(dòng)靠近最優(yōu)點(diǎn)運(yùn)行[7]。

在自優(yōu)化控制方法發(fā)展的十多年中，文獻(xiàn)[8-13]報(bào)道了不同方法求解最優(yōu)被控變量，如使用單獨(dú)的測(cè)量變量，或者它們的組合函數(shù)。本文采用的自優(yōu)化控制法為作者提出的NCO 近似法[13]，并將其推廣到間歇過(guò)程。

考慮如下連續(xù)過(guò)程的優(yōu)化問(wèn)題

式中，標(biāo)量J 表示需要最小化的性能指標(biāo)，u ∈ Rnu,d ∈Rnd分別是操縱變量和擾動(dòng)變量。y ∈Rny為過(guò)程輸出模型， G ∈Rng表示約束條件。假設(shè)uopt是最優(yōu)輸入，根據(jù)優(yōu)化理論，在uopt點(diǎn)應(yīng)該滿足如下KKT 條件[17]

式中，Φ 為拉格朗日函數(shù)，μ ∈ Rng為拉格朗日乘子。滿足上述KKT 條件的第一步是位于邊界的積極約束Ga為0（對(duì)應(yīng)非零拉格朗日乘子μi），Ga確定后可從式(2)約去μ 并得到[3]

式中，(·)+表示矩陣的Pseudo 偽逆矩陣。式(3)包含兩部分積極約束Ga和簡(jiǎn)約梯度?rJ，稱為連續(xù)過(guò)程優(yōu)化問(wèn)題的NCO 條件。

NCO 近似法的基本想法是利用系統(tǒng)處于最優(yōu)狀態(tài)時(shí)必須滿足NCO 條件（Ga和?rJ）為0 這一特征，將其直接作為被控變量（設(shè)定值為0）實(shí)現(xiàn)理想的自優(yōu)化控制。

假設(shè)NCO 條件（3）中的Ga部分可直接測(cè)量，則Ga直接作為被控變量進(jìn)行卡邊控制。而簡(jiǎn)約梯度?rJ 往往依賴于不可測(cè)擾動(dòng)d，因此不能在線測(cè)量。因此考慮建立?rJ 的回歸模型Z(y)，使用可測(cè)變量y估計(jì)?rJ，進(jìn)而設(shè)計(jì)被控變量c=Z(y)對(duì)其進(jìn)行間接控制，其中Z 為可測(cè)變量y 的線性或非線性函數(shù)。

以線性最小二乘擬合為例，被控變量記為c = θ0+θTy ( tf)。將過(guò)程M 個(gè)不同操作點(diǎn)的測(cè)量變量和簡(jiǎn)約梯度排列為矩陣形式

式中，下角標(biāo)i 表示第i 個(gè)操作點(diǎn)下的對(duì)應(yīng)數(shù)據(jù)，γi表示梯度值。為使自優(yōu)化控制效果在整個(gè)擾動(dòng)空間都有效，M 個(gè)操作點(diǎn)應(yīng)涵蓋所有擾動(dòng)工況。最小二乘擬合的結(jié)果為

除線性最小二乘回歸外，多項(xiàng)式回歸也可視擬合精度效果選擇性采用，也可以使用更加復(fù)雜的人工神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)回歸等。對(duì)于具有多個(gè)決策變量的系統(tǒng)，應(yīng)對(duì)每個(gè)簡(jiǎn)約梯度逐個(gè)擬合獲取多個(gè)被控變量。文獻(xiàn)[18]將自優(yōu)化控制問(wèn)題歸結(jié)為一個(gè)軟測(cè)量問(wèn)題，這里不再贅述。

2 間歇過(guò)程的批間自優(yōu)化控制

2.1 間歇過(guò)程優(yōu)化

本文研究一類參數(shù)不確定型間歇過(guò)程，將自優(yōu)化控制方法推廣到間歇過(guò)程的批間實(shí)時(shí)優(yōu)化。操作目標(biāo)表示為如下最優(yōu)化問(wèn)題

式中，tf為批次運(yùn)行時(shí)間，φ為標(biāo)量目標(biāo)函數(shù)，u(t)為連續(xù)的操縱輸入軌跡，x 為nx維狀態(tài)向量（初態(tài)x0）。J、y、d 的定義如前，其中y 測(cè)量值在批次終端獲得，即y(tf)。F、H 和G 分別為nF維動(dòng)態(tài)模型、nh維輸出方程和ng維終端約束。

最優(yōu)化問(wèn)題式(6)是一個(gè)動(dòng)態(tài)優(yōu)化問(wèn)題，可使用數(shù)值求解法如控制向量參數(shù)化法（CVP）求解得到最優(yōu)輸入軌跡uopt(t)。但由于系統(tǒng)受到不確定擾動(dòng)d的影響，因此使用標(biāo)稱模型求解得到的最優(yōu)輸入軌跡并非適用于所有工況，需根據(jù)實(shí)際操作工況進(jìn)行調(diào)整。對(duì)此，參數(shù)估計(jì)-重復(fù)優(yōu)化法[4]采用了一種直接的優(yōu)化思路，即首先根據(jù)當(dāng)前測(cè)量變量y 估計(jì)出不確定擾動(dòng)，然后更新模型并對(duì)其重新優(yōu)化求解。缺點(diǎn)是此方法的計(jì)算負(fù)荷較大，另外，由于參數(shù)估計(jì)和優(yōu)化兩者之間存在沖突，需要對(duì)系統(tǒng)進(jìn)行充 分激勵(lì)才能準(zhǔn)確估計(jì)出擾動(dòng)，因此需要折中考慮 兩者[4]。

2.2 輸入軌跡參數(shù)化

動(dòng)態(tài)優(yōu)化問(wèn)題式(1)可等價(jià)轉(zhuǎn)化為如下靜態(tài)優(yōu)化問(wèn)題[6]

操縱輸入軌跡u(t)的參數(shù)化方案對(duì)求解式(7)具有重要影響，通常做法是將批次運(yùn)行時(shí)間區(qū)間[0, tf]等間距分割成N 個(gè)時(shí)間間隔，在單個(gè)時(shí)間間隔內(nèi)u(t)為常數(shù)或多項(xiàng)式曲線，其中視為常數(shù)時(shí)即CVP法?？梢宰C明，當(dāng)N→∞時(shí)，參數(shù)化后的近似曲線可逼近任意形狀的曲線。如圖1 所示，經(jīng)CVP 法參數(shù)化后新的決策變量為u1,…,uN。

圖1 基于CVP 法的輸入軌跡參數(shù)化Fig.1 Input parameterization based on CVP (N=12)

提高CVP 法的算法精度可以通過(guò)增加N 來(lái)調(diào)節(jié)軌跡形狀，也增加了優(yōu)化問(wèn)題中的決策變量個(gè)數(shù)。因此CVP 法優(yōu)化性能的提高以增加計(jì)算復(fù)雜度為代價(jià)。這種方法的缺點(diǎn)在于，間歇過(guò)程中輸入軌跡的變化趨勢(shì)往往具有階段性，不同時(shí)間段內(nèi)的變化特性不同。對(duì)于較為平滑的部分，只需少量參數(shù)即可描述，而對(duì)較為“陡”的部分，才需較多的參數(shù)變量。對(duì)于CVP 法來(lái)說(shuō)，由于整個(gè)時(shí)間區(qū)間是等間距劃分的，因此需要大幅度增加N，以滿足最“陡”的那部分輸入軌跡曲線的需求，造成了計(jì)算資源的浪費(fèi)。

本文考慮另一種更為靈活的輸入軌跡參數(shù)化方法，下面結(jié)合圖2 進(jìn)行說(shuō)明。圖2 中均使用了非等間距的分隔方式對(duì)整個(gè)輸入軌跡曲線進(jìn)行劃分，兩者參數(shù)化后的決策變量均為u=[u1,u2,u3,t1,t2]T，圖2(a)中單個(gè)間隔內(nèi)u 均保持為常數(shù)，圖2 (b)中后兩段u 為斜坡形曲線。對(duì)于更為一般化的情況，也可以采用多項(xiàng)式曲線，n 階多項(xiàng)式由n+1 個(gè)參數(shù)確定。注意，這里3 段曲線的切換時(shí)間t1和t2也作為自由決策變量考慮，這是因?yàn)殚g歇過(guò)程在受到外界不確定擾動(dòng)的影響時(shí)，分階段軌跡曲線往往呈現(xiàn)出一種在時(shí)間尺度上的拉伸或壓縮的特性，此時(shí)使用時(shí)間參數(shù)作為自由變量更合理。

圖2 輸入軌跡參數(shù)化Fig.2 Input trajectory parameterization

圖3 使用二段2 次多項(xiàng)式曲線對(duì)圖1 所示的輸入曲線進(jìn)行了近似處理，此時(shí)輸入軌跡被參數(shù)化為7 個(gè)變量u=[a0,a1,a2,b0,b1,b2,t1]T，相比原來(lái)N=12可大大降低求解復(fù)雜度。

圖3 使用二段2 次多項(xiàng)式曲線進(jìn)行輸入軌跡參數(shù)化Fig.3 Input parameterization using two-stage second order polynomial functions

參數(shù)化法的性能好壞，很大程度上取決于參數(shù)化后的輸入軌跡是否吻合最優(yōu)輸入軌跡的形態(tài)特征，以盡量少的參數(shù)準(zhǔn)確描述。如果缺乏先驗(yàn)知識(shí)，仍然需要先借助如CVP 法求取最優(yōu)軌跡，然后再對(duì)其進(jìn)行合理的參數(shù)化。這是因?yàn)?，在系統(tǒng)受到不確定擾動(dòng)影響時(shí)，最優(yōu)輸入軌跡主要體現(xiàn)為這些參數(shù)的變化，而主要的形態(tài)特征往往不會(huì)發(fā)生較大改變，因此參數(shù)化法能大大降低優(yōu)化求解難度。本質(zhì)上來(lái)說(shuō)，參數(shù)化就是對(duì)輸入軌跡曲線施加額外的約束條件。一般來(lái)說(shuō)，如果參數(shù)化后的決策變量過(guò)少，導(dǎo)致新的軌跡曲線和目標(biāo)軌跡有較大差距，就可能引起一定程度的性能損失。此時(shí)需要設(shè)計(jì)人員進(jìn)行折中考慮，在性能損失可以承受的范圍內(nèi)，盡量減少?zèng)Q策變量數(shù)目。間歇過(guò)程的特點(diǎn)是決策變量數(shù)量龐大，各個(gè)輸入輸出通道之間相互耦合，對(duì)實(shí)現(xiàn)反饋控制是不利的。輸入軌跡參數(shù)化方案能夠大幅度減少?zèng)Q策變量數(shù)目，極大地利于反饋控制的實(shí)施，因此在本文背景下具有重要意義。

2.3 間歇過(guò)程批間自優(yōu)化

從前文分析可以看到，經(jīng)過(guò)穩(wěn)態(tài)轉(zhuǎn)化和輸入軌跡參數(shù)化后，最優(yōu)化問(wèn)題式(7)和穩(wěn)態(tài)自優(yōu)化控制問(wèn)題具有相同的形式。因此分析式(7)的最優(yōu)性條件，然后應(yīng)用NCO 近似法，可以實(shí)現(xiàn)間歇過(guò)程的批間自優(yōu)化控制。具體來(lái)說(shuō)，是建立y(tf)和最優(yōu)性條件之間的回歸模型c=Z(y(tf))，設(shè)定值恒為0，然后對(duì)其進(jìn)行反饋控制。

和連續(xù)過(guò)程自優(yōu)化的不同點(diǎn)在于：（1）連續(xù)過(guò)程的操縱變量往往是確定不變的，而間歇過(guò)程的決策變量經(jīng)過(guò)軌跡參數(shù)化得到，可以根據(jù)具體問(wèn)題進(jìn)行抉擇；（2）連續(xù)過(guò)程的控制器工作于連續(xù)或者高頻狀態(tài)，而間歇過(guò)程的控制器工作頻率為每批次進(jìn)行，以最簡(jiǎn)單的PI 控制器為例，第k+1 批次的操縱輸入uk+1的更新律為

式中，Kp和TI分別為比例系數(shù)和積分時(shí)間。

綜上，本文提出的間歇過(guò)程批間自優(yōu)化控制方法，具體實(shí)施步驟如下：

（1）使用傳統(tǒng)的數(shù)值優(yōu)化算法，在標(biāo)稱工作狀態(tài)下對(duì)間歇過(guò)程進(jìn)行優(yōu)化求解，獲得最優(yōu)輸入軌跡 ( )u t? ；

（2）觀察最優(yōu)輸入軌跡的形態(tài)特征，將其參數(shù)化為若干少量的決策變量u。然后根據(jù)參數(shù)化方案產(chǎn)生的輸入軌跡曲線u(t)，計(jì)算新的目標(biāo)函數(shù)J(u(t))。衡量J(u(t))相對(duì)于 ( ( ))J u t? 的損失程度，如果損失可以接受，則進(jìn)行下一步，否則重新調(diào)整參數(shù)化方案；

圖4 本文方法的控制方框圖Fig.4 Control diagram of proposed scheme

（3）基于間歇過(guò)程的動(dòng)態(tài)模型，以計(jì)算機(jī)仿真的方式在整個(gè)操作空間內(nèi)（包括擾動(dòng)空間和決策變量空間）進(jìn)行開(kāi)環(huán)實(shí)驗(yàn)，得到測(cè)量變量y(tf)和目標(biāo)函數(shù)J 相對(duì)于u 的梯度值?J/?u（可通過(guò)有限差分法得到）；

（4）以y(tf)為輸入，?J/?u 為輸出，建立回歸模型c=Z(y(tf))；

（5）設(shè)計(jì)控制器，被控變量為c，設(shè)定值0，控制輸入為u。

控制系統(tǒng)的在線運(yùn)行如圖4 所示，控制系統(tǒng)對(duì)第k 批次的終端輸出值yk(tf)進(jìn)行采樣，計(jì)算NCO的估計(jì)值并作為反饋控制器的被控變量ck，以批間方式跟蹤控制其設(shè)定值0。最后，使用控制輸出uk產(chǎn)生軌跡曲線uk(t)使用于下一批次。

3 仿真實(shí)例

3.1 過(guò)程描述

研究一個(gè)間歇反應(yīng)器[3]，原料A 和B 反應(yīng)生成產(chǎn)物C，同時(shí)B 自身合成生成副產(chǎn)物D，如圖5 所示。A 在反應(yīng)初始時(shí)刻一次性投放，B 在反應(yīng)過(guò)程中實(shí)時(shí)投放，流率為操縱變量F(t)

體系的模型方程為

圖5 間歇反應(yīng)器示意圖Fig.5 Schematic diagram of batch reactor

式中，cX表示物料X 的濃度，其他各符號(hào)含義參見(jiàn)表1。

表1 反應(yīng)過(guò)程參數(shù)Table 1 Parameters for reactor process

該過(guò)程的擾動(dòng)變量為主反應(yīng)常數(shù)k1，可能的波動(dòng)范圍為[0.02, 0.08]L·mol·min-1，過(guò)程運(yùn)行時(shí)k1不可測(cè)。反應(yīng)器過(guò)程的目標(biāo)函數(shù)表示為

式中，第1 項(xiàng)為產(chǎn)物C 的產(chǎn)量，第2 項(xiàng)代表了操作成本，ω 為2500 L2·min-1·mol-1。

3.2 輸入軌跡參數(shù)化

首先使用數(shù)值法求解標(biāo)稱工作點(diǎn)的最優(yōu)輸入軌跡，取CVP 法的N=15，使用Dynopt 工具箱[19]得到最優(yōu)輸入軌跡如圖6 左上角所示，此時(shí)最優(yōu)目標(biāo)函數(shù)為= 0.2412。從圖中可以觀察到，最優(yōu)輸入軌跡形狀的規(guī)律為：反應(yīng)初期反應(yīng)物B 的進(jìn)料流率F 較大，以利于反應(yīng)進(jìn)行得到更多產(chǎn)物C；反應(yīng)后期則減少B 的進(jìn)料，這是因?yàn)榉磻?yīng)生成C 需要一定時(shí)間，后續(xù)進(jìn)料來(lái)不及在反應(yīng)結(jié)束前產(chǎn)生效用，因此應(yīng)該少量進(jìn)料。

圖6 CVP 法（N=15）和輸入?yún)?shù)化法對(duì)比Fig.6 Comparison between CVP and input parameterization

根據(jù)上述規(guī)律，考慮如下一種簡(jiǎn)單的參數(shù)化方案：將整個(gè)反應(yīng)時(shí)間段分為兩個(gè)階段，第一階段[0, ts)實(shí)行恒定進(jìn)料F(t)=FB，第二階段[ts, tf]停止進(jìn)料F(t)=0。這樣整條輸入軌跡只需兩個(gè)參數(shù)[FB, ts]就完全確定。在此基礎(chǔ)上對(duì)系統(tǒng)進(jìn)行重新優(yōu)化求解，最優(yōu)目標(biāo)函數(shù)Jopt=0.2395，F(xiàn)B和ts分別為3.98×10-4L·min-1和234.6 min?？梢钥吹?，使用參數(shù)化方案將決策變量從15 減少至2 個(gè)后，只引起了目標(biāo)函數(shù)的少量損失，是非常值得的。兩者的最優(yōu)輸入軌跡和狀態(tài)變量對(duì)比如圖6 所示，可以看到cA、V、cC、cD的變化曲線幾乎是重疊的，cB略有差異，這是因?yàn)槲镔|(zhì)B 本身就是進(jìn)料。

3.3 批間自優(yōu)化控制

最優(yōu)決策變量的分布情況是假設(shè)擾動(dòng)變量已知，通過(guò)離線優(yōu)化求解得到的。而系統(tǒng)在線運(yùn)行時(shí)的真實(shí)擾動(dòng)值不可測(cè)，因此只能通過(guò)間接的方式實(shí)現(xiàn)優(yōu)化。

由于該間歇反應(yīng)器不存在積極約束部分，因此只需控制兩個(gè)梯度?J/?FB和?J/?ts為0即可保持系統(tǒng)最優(yōu)性。在擾動(dòng)變量k1∈[0.02, 0.08]的范圍內(nèi)，分別優(yōu)化求解兩個(gè)決策變量，得到它們的分布情況如圖7 所示?？梢钥吹?，F(xiàn)B的最優(yōu)值受擾動(dòng)k1變化較為敏感，而ts幾乎保持不變，變化范圍為標(biāo)稱工作點(diǎn)上下±2 min 內(nèi)波動(dòng)?？梢则?yàn)證，ts在小范圍內(nèi)變化時(shí)對(duì)目標(biāo)函數(shù)的影響不大，可以考慮取恒值ts=234.6 min 進(jìn)一步簡(jiǎn)化操作。

最終，使用NCO 近似法對(duì)該間歇反應(yīng)器進(jìn)行批間自優(yōu)化控制，只需得到可測(cè)變量y(tf)估計(jì)梯度?J/?FB的回歸模型作為被控變量。開(kāi)環(huán)狀態(tài)下，隨機(jī)變化決策變量和擾動(dòng)變量，仿真獲取100 組樣本數(shù)據(jù)，每組樣本數(shù)據(jù)包含可測(cè)變量 y(tf)和梯度?J/?FB。使用線性最小二乘擬合得到如下被控變量

圖7 最優(yōu)決策變量在整個(gè)擾動(dòng)范圍內(nèi)的分布Fig.7 Distributions of optimal decision variables in whole disturbance range

使用PI 控制器，設(shè)定c 的參考值為0 進(jìn)行批間自優(yōu)化控制。安排實(shí)驗(yàn)如下：第1～8 批次k1=0.02，第9～18 批次k1=0.08，第19～24 批次k1=0.053，實(shí)驗(yàn)結(jié)果顯示于圖8 和表2。從圖表中可以看到，方法具有良好的自優(yōu)化控制效果。擾動(dòng)發(fā)生后，由于被控變量c 在反饋控制器作用下調(diào)整Fs跟蹤設(shè)定值0，使反應(yīng)過(guò)程自動(dòng)向最優(yōu)工作點(diǎn)逼近，并且在若干批次內(nèi)收斂到各自的最終工作點(diǎn)。3 種擾動(dòng)情形中，k1= 0.02 時(shí)系統(tǒng)的目標(biāo)函數(shù)損失最大為1.18%。這部分損失主要來(lái)自兩方面：（1）從圖7中可以驗(yàn)證，k1=0.02 時(shí)ts的最優(yōu)值約為231.5 min，相對(duì)而言偏離固定值ts=234.6 min 最遠(yuǎn)，引起了一定程度損失；（2）對(duì)梯度?J/?FB擬合時(shí)引入了一定的擬合誤差，如需進(jìn)一步提高自優(yōu)化控制效果，可以考慮使用多項(xiàng)式擬合提高精度。k1= 0.08 和k1=0.053 時(shí)系統(tǒng)的目標(biāo)函數(shù)損失都很?。?.13%、0），可認(rèn)為已近似達(dá)到了最優(yōu)控制。

表2 批間自優(yōu)化控制結(jié)果和真實(shí)最優(yōu)值對(duì)比Table 2 Batch-to-batch self-optimizing control performances against true optimums

4 結(jié) 論

本文提出了一種基于NCO 近似法的間歇過(guò)程的批間自優(yōu)化控制策略。使用批次運(yùn)行終端的測(cè)量變量y(tf)，設(shè)計(jì)被控變量c=Z(y(tf))使系統(tǒng)能夠在后續(xù)批次中跟蹤設(shè)定值的同時(shí)，自動(dòng)運(yùn)行于（近似）最優(yōu)點(diǎn)。應(yīng)用本文方法時(shí)，一個(gè)重要步驟是對(duì)輸入軌跡曲線進(jìn)行準(zhǔn)確的參數(shù)化。從間歇反應(yīng)器研究的例子中可以看到，合理的參數(shù)化方案能夠在維持系統(tǒng)優(yōu)化性能的同時(shí)大大降低問(wèn)題復(fù)雜度。間歇反應(yīng)器的仿真結(jié)果表明，本文方法能有效實(shí)現(xiàn)間歇過(guò)程的批間自優(yōu)化控制。

圖8 3 種典型擾動(dòng)情形下的控制效果Fig.8 Control performances under 3 typical disturbance scenarios

[1]Srinivasan B, Palanki S, Bonvin D.Dynamic optimization of batch processes (Ⅰ): Characterization of the nominal solution [J].Comput.Chem.Eng., 2003, 27 (1): 1-26.

[2]Srinivasan B, Bonvin D, Visser E, Palanki S.Dynamic optimization of batch processes (Ⅱ): Role of measurements in handling uncertainty [J].Comput.Chem.Eng., 2003, 27 (1): 27-44.

[3]Chachuat B, Srinivasan B, Bonvin D.Adaptation strategies for real-time optimization [J].Comput.Chem.Eng., 2009, 33 (10): 1557-1567.

[4]Chen T, Liu Y, Chen J.An integrated approach to active model adaptation and on-line dynamic optimisation of batch processes [J].J.Proc.Control, 2013, 23 (10): 1350-1359.

[5]Marchetti A, Chachuat B, Bonvin D.Modifier-adaptation methodology for real-time optimization [J].Ind.Eng.Chem.Res., 2009, 48 (13): 6022-6033.

[6]Francois G, Srinivasan B, Bonvin D.Use of measurements for enforcing the necessary conditions of optimality in the presence of constraints and uncertainty [J].J.Proc.Control, 2005, 15 (6): 701-712.

[7]Skogestad S.Plantwide control: the search for the self-optimizing control structure [J].J.Proc.Control, 2000, 10 (5): 487-507.

[8]Halvorsen I J, Skogestad S, Morud J C, Alstad V.Optimal selection of controlled variables [J].Ind.Eng.Chem.Res., 2003, 42 (14): 3273-3284.

[9]Alstad V, Skogestad S.Null space method for selecting optimal measurement combinations as controlled variables [J].Ind.Eng.Chem.Res., 2007, 46 (3): 846-853.

[10]Kariwala V.Optimal measurement combination for local self-optimizing control [J].Ind.Eng.Chem.Res., 2007, 46 (11): 3629-3634.

[11]Kariwala V, Cao Y, Janardhanan S.Local self-optimizing control with average loss minimization [J].Ind.Eng.Chem.Res., 2008, 47 (4): 1150-1158.

[12]Alstad V, Skogestad S, Hori E S.Optimal measurement combinations as controlled variables [J].J.Proc.Control, 2009, 19 (1): 138-148.

[13]Ye L, Cao Y, Li Y, Song Z.Approximating necessary conditions of optimality as controlled variables [J].Ind.Eng.Chem.Res., 2013, 52 (2): 798-808.

[14]Ye L, Kariwala V, Cao Y.Dynamic optimization for batch processes with uncertainties via approximating invariant//8th IEEE Conference on Industrial Electronics and Applications [C].Melbourne, Australia, 2013: 1786-1791.

[15]Hu W, Mao J, Xiao G, Kariwala V.Selection of self-optimizing controlled variables for dynamic processes//Proceedings of 8th International Symposium on the Advanced Control of Chemical Processes [C].Singapore, 2012: 774-779.

[16]Jaschke J, Fikar M, Skogestad S.Self-optimizing invariants in dynamic optimization//50th IEEE Conference on Decision and Control and European Control Conference [C].Orlando, FL, 2011: 7753-7758.

[17]Edgar T F, Himmelblau D M, Lasdon L S.Optimization of Chemical Processes [M].2nd Ed.New York: McGraw-Hill, 2001.

[18]Ye Lingjian (葉凌箭), Ma Xiushui (馬修水).Optimal control strategy for chemical processes based on soft-sensoring technique [J].Journal of Zhejiang University: Engineering Science (浙江大學(xué)學(xué)報(bào): 工學(xué)版), 2013, (7): 1253-1257, 1266.

[19]Cizniar M, Fikar M, Latifi M A.A Matlab package for dynamic optimisation of processes//Proceedings of the 7th International Scientific-Technical Conference-PROCESS CONTROL[C].Kouty nad Desnou, Czech Republic, 2006.