基于平衡理論的模擬移動(dòng)床工藝參數(shù)魯棒尋優(yōu)

魏朋，陳珺，王志國，劉飛

（江南大學(xué)輕工過程先進(jìn)控制教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 無錫 214122）

引 言

模擬移動(dòng)床(simulated moving bed,SMB)是連續(xù)色譜分離技術(shù)的主要代表，其根據(jù)進(jìn)料液中各組分在固定相和流動(dòng)相中吸附和分配系數(shù)的不同實(shí)現(xiàn)彼此的分離[1]。模擬移動(dòng)床工藝由Broughton 等[2]于20世紀(jì)60年代首先提出，因其結(jié)合了固定床和真實(shí)移動(dòng)床的優(yōu)點(diǎn)，具有傳質(zhì)驅(qū)動(dòng)力大、分離能力強(qiáng)、可連續(xù)化生產(chǎn)和避免床層磨損等優(yōu)勢[3-5]。同時(shí)，與其他熱分離方法如蒸餾法、結(jié)晶法等相比，模擬移動(dòng)床消耗的能量更少，且具有經(jīng)濟(jì)環(huán)保、操作簡單的優(yōu)點(diǎn)[6]。因此，模擬移動(dòng)床在石油化工、制糖工業(yè)、生物制藥、精細(xì)化工等領(lǐng)域都有應(yīng)用[7-9]。

模擬移動(dòng)床的分離效果除了和固相吸附劑、洗脫液極性、色譜柱的填裝和尺寸有關(guān)外[10-11]，還和工藝參數(shù)，如切換時(shí)間和四個(gè)分離區(qū)域的流量等有關(guān)[12]。關(guān)于模擬移動(dòng)床工藝參數(shù)的尋優(yōu)一直是分離領(lǐng)域的研究重點(diǎn)之一。Shen等[13]在對二甲苯的分離中，提出了一種動(dòng)態(tài)順序二次規(guī)劃優(yōu)化算法以尋找最佳的工藝參數(shù)，使整個(gè)工藝的生產(chǎn)率得以提高。Li 等[14]在從工業(yè)糖漿中純化木寡糖時(shí)，采用NSGA-Ⅱ解決了具有不同目標(biāo)函數(shù)和約束的各種優(yōu)化問題。Matos 等[15]在雙萘酚對映體的分離過程中，提出了并行粒子群優(yōu)化算法，實(shí)現(xiàn)了生產(chǎn)率的最大化和洗脫液消耗的最小化。胡蓉等[16]采用多目標(biāo)教學(xué)優(yōu)化算法對C8 芳烴吸附分離過程進(jìn)行優(yōu)化，實(shí)現(xiàn)了對二甲苯收率最大化和解吸劑用量最小化。Degerman等[17]在通過離子交換色譜法純化免疫球蛋白G 時(shí)，在研究的工藝參數(shù)上進(jìn)行了敏感性和不確定性分析，分別用于確定哪些模型參數(shù)對產(chǎn)品性能有較大影響和產(chǎn)品因質(zhì)量問題失敗的概率，然后進(jìn)行魯棒優(yōu)化，從而確定更大的設(shè)計(jì)空間。Borg 等[18]分析了實(shí)驗(yàn)條件的不確定性對模型參數(shù)的影響，接著使用逆方法將模型校準(zhǔn)到實(shí)驗(yàn)數(shù)據(jù)來估計(jì)模型參數(shù)的分布，最后用兩個(gè)案例證明了該方法能用于分析不同的工藝參數(shù)，還能以一種更魯棒的方式設(shè)計(jì)出最佳過程。Nestola等[19]分析了當(dāng)顆?？紫堵屎蜕V柱體積都受到不確定因素影響時(shí)，在嚴(yán)格優(yōu)化理論的框架下提出了一種穩(wěn)健設(shè)計(jì)的一般程序，該程序使得工藝參數(shù)的最佳值能應(yīng)對模型參數(shù)在其不確定區(qū)間內(nèi)的所有變化。文獻(xiàn)[13-16]雖然可以根據(jù)不同的性能指標(biāo)以多目標(biāo)優(yōu)化算法得到一組最優(yōu)的工藝參數(shù)，但都是在模型參數(shù)固定不變的情況下得到的。實(shí)際中，由于工業(yè)現(xiàn)場情況存在眾多不確定因素，如進(jìn)料濃度變化、等溫線參數(shù)變化、循環(huán)泵的穩(wěn)定性和色譜柱的填充不一致等[20-21]，其導(dǎo)致的模型參數(shù)的變化會(huì)使所得初始工藝參數(shù)不滿足分離的性能指標(biāo)。因此，尋找一組魯棒的模擬移動(dòng)床工藝參數(shù)值，對確保生產(chǎn)過程平穩(wěn)運(yùn)行很有必要。文獻(xiàn)[17-19]在模擬移動(dòng)床的魯棒設(shè)計(jì)過程中分析了模型參數(shù)的不確定性對分離性能的影響，但未給出一種簡便的魯棒設(shè)計(jì)方法。

本文提出一種基于色譜平衡理論的模擬移動(dòng)床魯棒工藝參數(shù)尋優(yōu)方法。通過分析模型參數(shù)如何受到干擾的不確定性影響，從而將工作點(diǎn)移至故障更少、可靠性更高的參數(shù)區(qū)域。以果糖和葡萄糖的分離過程為例，首先對模擬移動(dòng)床系統(tǒng)建模，再進(jìn)一步分析在吸附劑失活和色譜柱堵塞情況下分離性能的變化，然后通過平衡理論尋優(yōu)獲得一組魯棒的工藝參數(shù)，最后驗(yàn)證在新工藝參數(shù)下能否有效克服擾動(dòng)對其分離性能的影響。

1 模擬移動(dòng)床吸附分離過程

圖1 所示為一個(gè)8 柱4 區(qū)的模擬移動(dòng)床色譜裝置結(jié)構(gòu)示意圖，每個(gè)色譜柱設(shè)有樣品進(jìn)出口，并通過多位閥和管道相互串聯(lián)。整個(gè)裝置分為四個(gè)區(qū)域，每個(gè)區(qū)域的流速不同且在混合物的分離中具有特定的功能。假設(shè)A 為強(qiáng)吸附組分，B 為弱吸附組分，各區(qū)劃分和功能如下[5]。

Ⅰ區(qū)是洗脫區(qū)，位于洗脫液進(jìn)口和提取液出口之間，其功能主要是從固定相中解吸組分A，實(shí)現(xiàn)固相吸附劑的再生，在Ⅰ區(qū)應(yīng)將固定相清洗干凈，使其不能將組分A帶入Ⅳ區(qū)。

Ⅱ區(qū)是精餾區(qū)，位于提取液出口和進(jìn)料口之間，其主要作用在于使固相吸附劑中的組分B 被組分A 不斷置換出來，使組分A 反復(fù)吸附、解吸而濃縮，從而使提取液中只包含組分A，不包含組分B。

Ⅲ區(qū)是吸附區(qū)，位于進(jìn)料口和提余液出口之間，其功能主要是吸附組分A，防止組分A 進(jìn)入提余液中，以保證組分A的高回收率。

Ⅳ區(qū)是過渡區(qū)，位于提余液出口和洗脫液進(jìn)口之間，其主要功能是使物料在Ⅲ區(qū)沒有被吸附的組分B被吸附。

模擬移動(dòng)床的具體分離過程如下：如圖1(a)所示，在初始時(shí)刻t=t0，待分離物料和洗脫液分別從柱4和柱5、柱8和柱1之間進(jìn)入，提取液和提余液分別從柱2 和柱3、柱6 和柱7 之間流出。經(jīng)過一個(gè)切換周期ts，模擬移動(dòng)床沿流動(dòng)相流動(dòng)方向切換多位閥以改變物料進(jìn)出口的位置，如圖1(b)所示。此時(shí)，待分離物料和洗脫液分別從柱5 和柱6、柱1 和柱2 之間進(jìn)入，提取液和提余液分別從柱3 和柱4、柱7 和柱8之間流出。多位閥經(jīng)過8次切換即1個(gè)循環(huán)，物料的進(jìn)出口位置又回到初始時(shí)刻的位置。模擬移動(dòng)床經(jīng)過多次切換之后，兩組分的色譜波峰會(huì)出現(xiàn)明顯分離。

圖1 傳統(tǒng)模擬移動(dòng)床操作示意圖Fig.1 Schematic diagram of conventional SMB operation

2 模擬移動(dòng)床的數(shù)學(xué)模型

模擬移動(dòng)床的數(shù)學(xué)模型由兩部分組成：單個(gè)色譜柱模型和節(jié)點(diǎn)平衡模型。本文選取色譜柱的平衡擴(kuò)散模型，該模型主要有以下假設(shè)[22]：

(1)吸附過程中溫度不變；

(2)組分在固定相與流動(dòng)相間的平衡是瞬間完成的；

(3)軸向擴(kuò)散和傳質(zhì)阻力的影響集中為表觀軸向擴(kuò)散系數(shù)；

(4)每個(gè)色譜柱填裝均勻，空隙率一致；

(5)沿色譜柱的流速恒定。

2.1 色譜柱中物料守恒方程

流動(dòng)相的質(zhì)量守恒方程為：

式中，ci和qi分別表示組分i(i=A 或B)在流動(dòng)相與固定相中的濃度；v表示色譜柱中流動(dòng)相的流動(dòng)速度；ε表示床層的空隙率；Da表示表觀軸向擴(kuò)散系數(shù)；z、t分別表示空間和時(shí)間坐標(biāo)。

2.2 相間吸附平衡方程

吸附平衡方程描述溶質(zhì)分子的吸附和解吸，本文選取線性吸附等溫線，即：

式中，Hi表示組分i的Henry常數(shù)。

2.3 初始和邊界條件

色譜柱分離過程一般都是從空柱開始的，因此，初始條件為：

2.4 節(jié)點(diǎn)平衡模型

洗脫液節(jié)點(diǎn)處：

2.5 性能指標(biāo)

評估分離過程的性能指標(biāo)有：純度、回收率、生產(chǎn)率、洗脫液消耗等。

純度反映產(chǎn)品質(zhì)量，其定義為在一個(gè)切換周期內(nèi)提取液和提余液中目標(biāo)組分的平均濃度與總濃度的比值，則提取液中組分A的純度為:

3 基于平衡理論的參數(shù)尋優(yōu)

上述模型由一組偏微分方程和代數(shù)方程相耦合而成。偏微分方程描述系統(tǒng)的動(dòng)態(tài)過程，一般直接對偏微分方程求解難度較大，因此常先對偏微分方程進(jìn)行離散化處理，將其轉(zhuǎn)化成常微分方程組，然后求解常微分方程組的數(shù)值解。已有學(xué)者采取不同的數(shù)值方法求解偏微分方程，例如有限元法[11]、有限體積法[23]、正交配置法[24]、時(shí)空守恒元解元法[25]和有限元正交配置法[26]等。

有限元正交配置法是有限元法和正交配置法的結(jié)合，首先將色譜柱沿軸向劃分為I個(gè)有限元，然后在每個(gè)有限元上應(yīng)用正交配置法，每個(gè)有限元中的內(nèi)配置點(diǎn)個(gè)數(shù)為N[27]，如圖2所示。為了避免重復(fù)計(jì)算，假定濃度函數(shù)及其一階導(dǎo)數(shù)在相鄰兩個(gè)有限元的邊界處是連續(xù)的。色譜柱的邊界條件分別對應(yīng)第一個(gè)有限元的第一個(gè)配置點(diǎn)和最后一個(gè)有限元的最后一個(gè)配置點(diǎn)。

圖2 有限元上的正交配置Fig.2 Orthogonal collocation on finite element

根據(jù)式(1)，對第i個(gè)有限元的第n個(gè)內(nèi)部配置點(diǎn)，可將方程離散為式（14）：

式中，i=1,2,…,I，n=2,3,…,N；h為每個(gè)有限元的長度，A和B為通過正交配置法獲得的系數(shù)矩陣，F(xiàn)為相比，定義為：F=(1 -ε)/ε。

有限元邊界處的連續(xù)條件為：

本文采用5 個(gè)有限元，每個(gè)有限元的內(nèi)配置點(diǎn)個(gè)數(shù)為3，可將模型離散為42個(gè)微分代數(shù)方程，然后進(jìn)行求解。

平衡理論是目前應(yīng)用廣泛的模擬移動(dòng)床工藝參數(shù)尋優(yōu)方法，主要用于確定各區(qū)的最佳流量以滿足混合物的有效分離，它由Storti 等[28]于1993 年提出。該理論定義了流量比mj的概念，即某區(qū)液相凈流量與固相凈流量的比值，化簡后如下：

值得注意的是，平衡理論只考慮了Ⅱ區(qū)和Ⅲ區(qū)流量比的影響，Ⅰ區(qū)和Ⅳ區(qū)的流量比是事先確定的，從而將m2-m3平面劃分為不同的分離區(qū)域，如圖3 所示，其中由w、a和b三點(diǎn)確定的三角形區(qū)域?yàn)閮山M分的完全分離區(qū)域。

根據(jù)圖3 中不同的分離區(qū)域，可大概確定各區(qū)的流量，為模擬移動(dòng)床模型的求解提供工藝參數(shù)的理論指導(dǎo)，再結(jié)合前述的模型求解方法即有限元正交配置法，則可計(jì)算出相應(yīng)的性能指標(biāo)。

圖3 線性吸附等溫線描述的系統(tǒng)在m2-m3平面上的不同分離區(qū)域Fig.3 Different separation regions of the system described by linear adsorption isotherms on the m2-m3 plane

在實(shí)際中，系統(tǒng)常受到非受控不確定因素的干擾。因此，工藝參數(shù)的設(shè)計(jì)尤為重要，需使生產(chǎn)過程在干擾方面是魯棒的。由于干擾的性質(zhì)，系統(tǒng)的魯棒性不可能在所有情況下都能得到保證。所以，在一組樣本數(shù)據(jù)中，需指定某一置信水平，使生產(chǎn)能以該置信水平滿足性能指標(biāo)的要求。模擬移動(dòng)床的尋優(yōu)求解流程如圖4所示。

圖4 模擬移動(dòng)床的尋優(yōu)求解流程圖Fig.4 Flow chart of optimizing solution for simulated moving bed

4 仿真實(shí)驗(yàn)

本文選擇果葡糖漿的模擬移動(dòng)床分離過程驗(yàn)證所提方法的有效性。果糖是一種高甜度的糖，被廣泛添加到果汁飲料當(dāng)中，既能夠增加甜度，還能夠保持果味的原來風(fēng)味。此外，果糖也可以應(yīng)用在醫(yī)學(xué)領(lǐng)域，比如被制成注射液，還可以作為重要醫(yī)藥的合成起始物[29]。雖然在果葡糖漿的分離過程中洗脫液的消耗是很重要的一個(gè)性能指標(biāo)，也有學(xué)者表明SSMB 和ISMB 可以減少洗脫液的消耗[12,30]，但本文以果葡糖漿的分離過程為例不會(huì)影響驗(yàn)證所提方法的有效性。文獻(xiàn)[6]介紹了從果葡糖漿中分離果糖，根據(jù)工藝要求，當(dāng)提取液中果糖的純度達(dá)到95%以上時(shí)，得到的產(chǎn)品才能滿足質(zhì)量要求；當(dāng)進(jìn)料中果糖的回收率達(dá)到95%以上時(shí)，才能滿足生產(chǎn)的經(jīng)濟(jì)效益。由于文獻(xiàn)[6]未能給出軸向擴(kuò)散系數(shù)的值，故采用一典型的數(shù)值，具體的模型參數(shù)和初始工藝參數(shù)如表1 所示，其中組分A 為果糖，組分B為葡萄糖。

表1 模擬移動(dòng)床模型參數(shù)及工藝參數(shù)Table 1 Model parameters and process parameters of simulated moving bed

4.1 模擬移動(dòng)床運(yùn)行過程分析

為了更直觀地介紹模擬移動(dòng)床的分離過程和驗(yàn)證所采用模型的準(zhǔn)確性，通過實(shí)驗(yàn)分析了模擬移動(dòng)床在一個(gè)循環(huán)內(nèi)的軸向濃度變化，如圖5 所示?？梢钥闯雒拷?jīng)過一次切換，組分A、B 的濃度都在增大，且組分A 向提取口移動(dòng)，組分B 向提余口移動(dòng)。隨著切換次數(shù)的增加，組分向各自出口的移動(dòng)速度變慢，最后基本上不再移動(dòng)且濃度不再變化，此時(shí)系統(tǒng)到達(dá)了循環(huán)穩(wěn)態(tài)。

圖5 模擬移動(dòng)床的運(yùn)行過程(各子圖均為每個(gè)運(yùn)行周期末尾時(shí)刻的濃度分布)Fig.5 Running process of the simulate moving bed

圖6 為模擬移動(dòng)床經(jīng)過72 次切換后到達(dá)的循環(huán)穩(wěn)態(tài)。需要注意的是，此時(shí)到達(dá)的循環(huán)穩(wěn)態(tài)是周期性的，其周期等于模擬移動(dòng)床的切換周期。系統(tǒng)在到達(dá)循環(huán)穩(wěn)態(tài)后，就可以分別在提取口和提余口得到高純度的分離產(chǎn)品。根據(jù)本文方法，在表1 中的初始工藝參數(shù)下，實(shí)驗(yàn)得到的提取口組分A 的平均純度為99.31%，提余口組分B 的平均純度為99.83%，然后持續(xù)進(jìn)原料與洗脫液，系統(tǒng)可實(shí)現(xiàn)連續(xù)化生產(chǎn)。實(shí)驗(yàn)結(jié)果表明，模擬移動(dòng)床可有效分離果葡糖漿混合物，且能夠分別在兩出口得到高純度的果糖和葡萄糖產(chǎn)品。

圖6 模擬移動(dòng)床的穩(wěn)態(tài)濃度分布(穩(wěn)態(tài)下每個(gè)切換周期末尾時(shí)刻的濃度分布)Fig.6 Steady-state concentration distribution of the simulate moving bed

4.2 擾動(dòng)作用下的魯棒設(shè)計(jì)

在模擬移動(dòng)床的運(yùn)行過程中會(huì)存在一些難以避免的擾動(dòng)，這些擾動(dòng)會(huì)對模型參數(shù)造成不確定性的影響，如吸附劑失活、色譜柱堵塞、循環(huán)泵誤差和溫度變化等擾動(dòng)會(huì)分別對色譜柱的空隙率ε、有效柱長L、各分離區(qū)的體積流量Qj和Henry系數(shù)Hi造成不確定性的影響?，F(xiàn)假設(shè)吸附劑失活和色譜柱堵塞為主要擾動(dòng)，并假設(shè)受其影響的ε和L是相互獨(dú)立且服從正態(tài)分布的隨機(jī)變量，這使得最大變化的可能性較小。現(xiàn)取300組數(shù)據(jù)為樣本，其中ε和L的均值分別為0.38 和53.6，標(biāo)準(zhǔn)差分別為均值3%的波動(dòng)量即0.0114 和1.608，計(jì)算可得ε和L落在各自區(qū)間的置信度分別為98.97%和99.57%。

在給定的初始工藝參數(shù)下(即QD，QE，QF，QR保持不變)，隨機(jī)擾動(dòng)對模擬移動(dòng)床的分離性能產(chǎn)生的影響如圖7 所示。從圖7 可知，根據(jù)工藝要求，許多擾動(dòng)使工藝無法滿足純度或回收率的要求，或兩者都有。就產(chǎn)品的純度而言，落在虛線以下即純度小于95%的點(diǎn)在300 組樣本數(shù)據(jù)中占57 個(gè)，生產(chǎn)失敗的概率高達(dá)19%，這在實(shí)際生產(chǎn)中是絕對不允許的。因此，需要尋找魯棒的工藝參數(shù)，使分離在擾動(dòng)作用下不會(huì)冒生產(chǎn)失敗的風(fēng)險(xiǎn)。

圖7 隨機(jī)擾動(dòng)作用下生成的300組正態(tài)分布值(L,ε)運(yùn)行在初始工藝參數(shù)下得到的產(chǎn)品純度和回收率(圖(a)和(b)中的藍(lán)點(diǎn)代表各自的性能達(dá)標(biāo)，紅點(diǎn)代表未達(dá)標(biāo)；圖(c)中的紅點(diǎn)代表模型參數(shù)的原始值，藍(lán)點(diǎn)代表正態(tài)分布值)Fig.7 The product purity and recovery rate obtained by running 300 sets of normal distribution values(L,ε)generated under the random disturbance action under the initial process parameters

圖8 所示為模擬移動(dòng)床工作在初始工藝參數(shù)時(shí)，300 組樣本數(shù)據(jù)(L,ε)在m2-m3平面所對應(yīng)的點(diǎn)的分布，從圖8 可知有75 個(gè)點(diǎn)落在了完全分離區(qū)域的外面，概率為34.67%。在分析這些樣本數(shù)據(jù)生成的性能指標(biāo)時(shí),發(fā)現(xiàn)導(dǎo)致提取液中產(chǎn)品純度不合格的點(diǎn)(L,ε)所對應(yīng)的m2均小于HB，即使是m2≥HB時(shí)，m2靠近HB的點(diǎn)所對應(yīng)的提取液純度也相對較低。

圖8 初始工藝參數(shù)下樣本點(diǎn)(L,ε)對應(yīng)的流量比在m2-m3平面上的分布(紅點(diǎn)代表落在了完全分離區(qū)域的外面，藍(lán)點(diǎn)代表落在了完全分離區(qū)域的里面)Fig.8 Distribution of the flow ratio corresponding to the sample point(L,ε)on the m2-m3 plane at the initial process parameters

文獻(xiàn)[31]結(jié)果表明，在保證Ⅰ區(qū)和Ⅳ區(qū)流量比不變的情況下，若提取液流量QE減小，則提取液中強(qiáng)吸附組分的純度會(huì)得到提高，但QE減小過多，會(huì)導(dǎo)致強(qiáng)吸附組分損失在提余液中，影響工藝的生產(chǎn)率和回收率。增加進(jìn)料QF,雖然可以提高工藝的生產(chǎn)率和減少洗脫液的消耗，但若進(jìn)料量的增加超出了色譜柱的吸附容量，則會(huì)有更多的強(qiáng)吸附組分隨提余液排出，使強(qiáng)吸附組分的收率降低。若進(jìn)一步增加進(jìn)料，則含有雜質(zhì)的再生洗脫液會(huì)穿過過渡區(qū)，降低提取液中強(qiáng)吸附組分的純度。

根據(jù)文獻(xiàn)[31]中結(jié)果表明的QE和QF對分離性能的影響，再結(jié)合模擬移動(dòng)床的節(jié)點(diǎn)模型，需要適當(dāng)減小QE和QF，使模擬移動(dòng)床工作在新的工藝參數(shù)時(shí)，有盡可能多的樣本點(diǎn)(L,ε)對應(yīng)的流量比落在完全分離區(qū)域內(nèi)，且離區(qū)域邊緣有一定的距離，才能保證產(chǎn)品的純度和回收率。根據(jù)式(7)～式(10)和式(19)計(jì)算，本次工藝參數(shù)的魯棒值選取為：QD= 0.0414 ml/s，QE= 0.0306 ml/s，QF= 0.0148 ml/s，QR= 0.0256 ml/s。

圖9所示為模擬移動(dòng)床工作在新的魯棒工藝參數(shù)時(shí)得到的流量比分布，從圖中可以看出300 組樣本數(shù)據(jù)中有291個(gè)點(diǎn)即97%的點(diǎn)落在了完全分離區(qū)域內(nèi)，符合預(yù)期。

圖9 魯棒工藝參數(shù)下樣本點(diǎn)(L,ε)對應(yīng)的流量比在m2 - m3平面上的分布(紅點(diǎn)代表落在了完全分離區(qū)域的外面，藍(lán)點(diǎn)代表落在了完全分離區(qū)域的里面)Fig.9 Distribution of the flow ratio corresponding to the sample point(L,ε)on the m2 - m3 plane at the robust process parameters

模擬移動(dòng)床運(yùn)行在魯棒工藝參數(shù)下得到的產(chǎn)品純度和回收率分布如圖10 所示，從圖可得，產(chǎn)品純度低于95%的點(diǎn)只有3 個(gè)，即批次失敗的概率僅有1%，而且總體純度水平較之前有了很大的提高；此外，產(chǎn)品回收率達(dá)到95%的置信水平為93.67%，滿足了原料的回收要求?？傊?，選取的魯棒工藝參數(shù)能夠很好地克服擾動(dòng)的存在，使整個(gè)工藝過程能夠穩(wěn)定生產(chǎn)。與此同時(shí)，整個(gè)過程的生產(chǎn)率有所降低，也增加了洗脫液的消耗量，但這點(diǎn)損耗與流程工業(yè)的平穩(wěn)生產(chǎn)相比，是完全可以接受的。

圖10 隨機(jī)擾動(dòng)作用下生成的300組正態(tài)分布值(L,ε)運(yùn)行在魯棒工藝參數(shù)下得到的產(chǎn)品純度和回收率(圖(a)和(b)中的藍(lán)點(diǎn)代表各自的性能達(dá)標(biāo)，紅點(diǎn)代表未達(dá)標(biāo)；圖(c)中的紅點(diǎn)代表模型參數(shù)的原始值，藍(lán)點(diǎn)代表正態(tài)分布值)Fig.10 The product purity and recovery rate obtained by running 300 sets of normal distribution values(L,ε)generated under the random disturbance action under the robust process parameters

5 結(jié) 論

本文提出一種基于色譜平衡理論的模擬移動(dòng)床魯棒工藝參數(shù)尋優(yōu)方法。在采用有限元正交配置法求解模型的基礎(chǔ)上，分析了擾動(dòng)對工藝分離性能的影響，再以平衡理論為依據(jù)，獲得了一組魯棒工藝參數(shù)。通過以糖醇分離行業(yè)常見的果糖和葡萄糖的分離過程為例，假定吸附劑失活和色譜柱堵塞為主要影響模型參數(shù)的擾動(dòng)，實(shí)驗(yàn)結(jié)果表明所提方法得到的新工藝參數(shù)具有良好的魯棒性，同時(shí)工藝能以99%的置信水平滿足產(chǎn)品的純度要求，93.67%的置信水平滿足原料的回收要求。此外，該方法相比其他多目標(biāo)優(yōu)化算法，具有算法簡便、尋優(yōu)過程花費(fèi)時(shí)間少的優(yōu)勢。本文提出的研究方法除了可用于研究吸附劑失活和色譜柱堵塞對分離性能的影響外，還可用于研究循環(huán)泵誤差、溫度變化等不確定性因素對分離性能的影響，以設(shè)計(jì)更加魯棒的工藝過程。

符 號(hào) 說 明

ci——組分i在流動(dòng)相中的濃度，g/ml

ci,k——分別表示組分i在洗脫液、提取液、進(jìn)料液和提余液中的濃度(k=D,E,F,R)，g/ml

cˉi——組分i在出口處的平均濃度，g/ml

cini,j,couti,j——分別表示組分i在j區(qū)的入口濃度和出口濃度，g/ml

Da——表觀軸向擴(kuò)散系數(shù)

F——相比

Hi——組分i的Henry常數(shù)

h——每個(gè)有限元的長度，cm

I——有限元的個(gè)數(shù)

L——色譜柱的柱長，cm

mj——j區(qū)的流量比

N——限元的內(nèi)配置點(diǎn)個(gè)數(shù)

Qj——洗脫液、提取液、進(jìn)料液、提余液的體積流量(j=D，E，F(xiàn)，R)，ml/s

Qj——分離區(qū)域Ⅰ、Ⅱ、Ⅲ、Ⅳ的體積流量(j=Ⅰ,Ⅱ,Ⅲ,Ⅳ)，ml/s

qi——組分i在固定相中的濃度，g/ml

t——時(shí)間坐標(biāo)

ts——模擬移動(dòng)床的切換周期，s

v——譜柱中流動(dòng)相的流動(dòng)速度，cm/s

z——空間坐標(biāo)

ε——床層空隙率