黃冬，趙民帥，羅雄麟

?

乙烯精餾塔控制中降液管時(shí)滯效應(yīng)影響分析

黃冬，趙民帥，羅雄麟

（中國(guó)石油大學(xué)（北京）自動(dòng)化系，北京 102249）

在精餾塔動(dòng)態(tài)建模中忽略了降液管容積對(duì)液相流動(dòng)及傳質(zhì)的滯后，導(dǎo)致模型與實(shí)際精餾塔存在明顯差異。針對(duì)某實(shí)際乙烯精餾塔，通過機(jī)理分析建模，建立理論降液管模型，并在原精餾塔模型基礎(chǔ)上構(gòu)建了考慮降液管時(shí)滯效應(yīng)的乙烯精餾塔動(dòng)態(tài)模型。通過仿真，將該模型與原模型的動(dòng)態(tài)特性進(jìn)行比較，其差異性說明考慮降液管能夠更加準(zhǔn)確地把握精餾過程的動(dòng)態(tài)特性。另外，對(duì)兩種動(dòng)態(tài)模型分別設(shè)計(jì)控制器，整定得到的控制器參數(shù)差別很大，說明忽略降液管的影響造成在此基礎(chǔ)上設(shè)計(jì)的控制器可能不適用于實(shí)際裝置。因此，在動(dòng)態(tài)建模中考慮降液管能夠更加準(zhǔn)確地對(duì)精餾塔進(jìn)行分析、控制和優(yōu)化，具有一定的現(xiàn)實(shí)意義。

過程系統(tǒng)；動(dòng)態(tài)仿真；計(jì)算機(jī)模擬；蒸餾；降液管

引 言

精餾塔是化工過程中重要的傳熱傳質(zhì)單元[1-2]，精餾塔的模型化也受到廣泛關(guān)注[3-5]。模型化主要是以一定的形式相對(duì)準(zhǔn)確地反映精餾塔的內(nèi)部特性和規(guī)律，為精餾塔的設(shè)計(jì)、分析、控制和優(yōu)化提供重要手段[6-10]。精餾塔的機(jī)理模型分為平衡級(jí)模型和非平衡級(jí)模型。平衡級(jí)模型應(yīng)用比較廣泛，而且理論成熟[11]；非平衡級(jí)模型從20世紀(jì)80年代開始興起，但主要局限于學(xué)術(shù)研究方面[12]。精餾塔模型也可分為穩(wěn)態(tài)模型和動(dòng)態(tài)模型，穩(wěn)態(tài)模型主要應(yīng)用于精餾塔的穩(wěn)態(tài)計(jì)算[13]，動(dòng)態(tài)模型可應(yīng)用于動(dòng)態(tài)優(yōu)化和控制設(shè)計(jì)[14]。

精餾塔動(dòng)態(tài)模型的研究中，一般只考慮理論板上液相的動(dòng)態(tài)特性，而忽略了降液管對(duì)液相變化的影響[15-16]。對(duì)于降液管的研究，Zhang等[17]針對(duì)降液管液相混合程度進(jìn)行了計(jì)算，認(rèn)為降液管內(nèi)液體處于不完全混合狀態(tài)。Weiler等[18]研究了降液管液相不完全混合對(duì)精餾塔板效率的影響。黃雪雷等[19]對(duì)板式塔弓形降液管液相流場(chǎng)進(jìn)行了數(shù)值模擬。曹振恒等[20]將傳統(tǒng)篩板塔的弓形降液管改為均分布于塔板篩孔間的多降液管結(jié)構(gòu)，提高了塔板的傳質(zhì)效率。這些研究針對(duì)單個(gè)降液管結(jié)構(gòu)等方面的計(jì)算，未考慮降液管的存在對(duì)于精餾塔動(dòng)態(tài)特性的影響。

本工作通過分析降液管的結(jié)構(gòu)建立降液管的動(dòng)態(tài)機(jī)理模型，并引入到乙烯精餾塔模型中，研究降液管的時(shí)滯效應(yīng)對(duì)乙烯精餾塔溫度和質(zhì)量的動(dòng)態(tài)性能影響，同時(shí)討論了有無降液管時(shí)塔頂和塔底溫度控制器參數(shù)的差異。

1 問題描述

乙烯精餾塔是乙烯生產(chǎn)裝置中的重要操作單元，由于其負(fù)荷大、產(chǎn)品質(zhì)量高等特點(diǎn)，對(duì)操作的正確性和準(zhǔn)確性要求很高[21]。如圖1所示，對(duì)某精餾塔進(jìn)行操作時(shí)，當(dāng)塔頂回流量變化時(shí)，第1塊塔板液相發(fā)生變化，通過第1個(gè)降液管，進(jìn)而改變第2塊塔板的液相進(jìn)料，第2塊塔板的液相再經(jīng)過第2個(gè)降液管，到達(dá)第3塊塔板，重復(fù)此過程，直到到達(dá)塔底。若將一個(gè)降液管視為一個(gè)重力作用下的自流式容器，則降液管內(nèi)液相停留時(shí)間為液相體積與流量的比值，即每個(gè)降液管都存在滯后的現(xiàn)象。單個(gè)降液管的停留時(shí)間有限，一般在3～5 s以上，但乙烯精餾塔塔板數(shù)和降液管數(shù)較多，導(dǎo)致當(dāng)塔頂回流量變化時(shí)塔底流量、溫度等參數(shù)的變化出現(xiàn)大滯后，在實(shí)際生產(chǎn)中這種滯后往往能達(dá)到20 min以上。

在建立乙烯精餾塔穩(wěn)態(tài)模擬的過程中，基于平衡級(jí)假設(shè)和全混級(jí)假設(shè)，根據(jù)每個(gè)塔板的質(zhì)量平衡、能量平衡和相平衡方程等計(jì)算得到穩(wěn)態(tài)結(jié)果。動(dòng)態(tài)模擬在穩(wěn)態(tài)模擬的基礎(chǔ)上建立，只考慮各個(gè)理論板上的汽液相條件，忽略了降液管的液相流動(dòng)及傳質(zhì)的滯后，導(dǎo)致仿真得到的精餾塔動(dòng)態(tài)特性無法反映實(shí)際情況。因此，有必要根據(jù)降液管的結(jié)構(gòu)特征和內(nèi)部機(jī)理建立精確的降液管機(jī)理模型，以完善乙烯精餾塔的動(dòng)態(tài)模擬過程。

2 考慮降液管的乙烯精餾塔模型

考慮降液管模型，首先要確定降液管的布置情況，規(guī)定了板上液體流動(dòng)的途徑。其降液管的布置情況一般包括U形流、單溢流、雙溢流和階梯流等類型。本工作針對(duì)某實(shí)際乙烯精餾塔，該精餾塔采用雙溢流的方式，如圖1所示。這種方式占塔板面積較多，一般用于負(fù)荷較大，塔徑2 m以上的精餾塔。如圖2所示，來自上一層塔板的液相分別從兩側(cè)的降液管進(jìn)入受液盤，流過半個(gè)鼓泡區(qū)長(zhǎng)度后進(jìn)入塔板中間的降液管，在下一塔板上液相經(jīng)降液管流入受液盤后流向兩側(cè)，再次通過半個(gè)鼓泡區(qū)長(zhǎng)度后進(jìn)入兩側(cè)的降液管，如此重復(fù)，直到塔底。

2.1 降液管模型

降液管結(jié)構(gòu)示意圖如圖2所示。傳統(tǒng)板式精餾塔的平衡級(jí)計(jì)算，通過引入全塔板效率，將1/塊實(shí)際板和降液管的分離效果等效為1塊理論板的分離效果，忽略了1/個(gè)降液管的質(zhì)量?jī)?chǔ)存容積。本工作引入“理論降液管”的概念，是指理論板對(duì)應(yīng)的降液管。由于理論板滿足平衡級(jí)假設(shè)和全混級(jí)假設(shè)，理論降液管中只存在液相，沒有分離功能，只進(jìn)行液相傳質(zhì)，彌補(bǔ)了由于忽略降液管容積造成的偏差。另外假設(shè)降液管中液相完全混合。

1—downcomer; 2—hole; 3—bubbling zone; 4—spume zone; 5—seal pan

圖2為降液管結(jié)構(gòu)的俯視圖和正視圖。其中，為精餾塔直徑，m；w為堰長(zhǎng)，m；d為降液管寬度，m；s為破沫區(qū)寬度，m；為半個(gè)鼓泡區(qū)寬度，m；w為出口堰高，m；′w為入口堰高，m；ow為堰上液層高度，m；0為降液管底隙高度，m；T為塔板距離，m；d為降液管液位高度，m；為塔板上流入降液管的液相摩爾流量，kmol·h-1；d為降液管底隙出口摩爾流量，kmol·h-1；為塔板上汽相摩爾流量，kmol·h-1；d為降液管底隙液體流速，m·s-1；為塔板壓力，Pa。該乙烯精餾塔采用雙溢流類型，因此存在兩種類型的塔板及降液管結(jié)構(gòu)（Ⅰ和Ⅱ）。對(duì)于Ⅰ類型，受液盤在塔板中間，降液管在塔板兩側(cè)，液相由中間流向兩側(cè)；對(duì)于Ⅱ類型，受液盤在塔板兩側(cè)，降液管在塔板中間，液相由兩側(cè)流向中間。另外，可近似認(rèn)為兩種類型的塔板的降液管面積相等。

對(duì)于由塊理論板組成的精餾塔，第1塊板為冷凝器，第塊板為再沸器，理論降液管個(gè)數(shù)為-3。

傳統(tǒng)板式精餾塔平衡級(jí)模型中只研究了塔板鼓泡區(qū)的分離過程的計(jì)算。本工作主要在原平衡級(jí)模型基礎(chǔ)上對(duì)塔板上的降液管進(jìn)行分析建模，并將降液管動(dòng)態(tài)模型與平衡級(jí)動(dòng)態(tài)模型結(jié)合起來，實(shí)現(xiàn)了動(dòng)態(tài)仿真。

圖2(b)中，對(duì)第塊塔板下的第個(gè)降液管建立模型方程。

（1）物料平衡 對(duì)第個(gè)降液管列寫總物料守恒方程和各組分質(zhì)量守恒方程

(2)

將式(1)代入式(2)，可得

式中，d為降液管截面積，m2；d為降液管內(nèi)液相摩爾密度，kmol·m-3；d為降液管內(nèi)液位高度，m；(,)為第塊塔板上流下液相的第種組分的摩爾分?jǐn)?shù)；d(,)為第個(gè)底隙出口液體的第種組分的摩爾分?jǐn)?shù)。

（2）熱量平衡 對(duì)第個(gè)降液管列寫熱量平衡方程

將式(1)代入式(4)，可得

(5)

式中，dL為降液管內(nèi)液相焓值，kJ·kmol-1；L為塔板上液相焓值，kJ·kmol-1。

（3）其他方程 將降液管等效為重力作用下的自流式容器，則降液管底隙流速d滿足

式中，dM為降液管內(nèi)液相質(zhì)量密度，kg·m-3；L為塔板上液層高度，m；d為降液管底隙液體流速，m·s-1；為塔板壓力，Pa；Δd為降液管底隙壓頭損失，Pa。

而且滿足

式中，w為組分摩爾質(zhì)量，kJ·kmol-1；為常系數(shù)。

降液管底隙面積為

式中，0為降液管底隙出口面積，m2。

降液管底隙高度0為

降液管底隙出口流量滿足

(10)

根據(jù)塔板水力學(xué)方程

由式(6)～式(11)，可得

(12)

對(duì)式(12)求導(dǎo)可得

(14)

(16)

該模型狀態(tài)變量包括d、d、d，輸入變量包括、、。即當(dāng)上一塊板進(jìn)入的液相已知的條件下，降液管液位、底隙流量、組成、溫度等變量唯一確定。

2.2 乙烯精餾塔模型

在全混級(jí)假設(shè)和平衡級(jí)假設(shè)的基礎(chǔ)上，忽略塔板上方的汽相滯留量，并認(rèn)為塔板上的熱量傳遞非常迅速。如圖2(b)所示，分析典型塔板，考慮降液管情況下的MESH方程如下。

（1）組分物料平衡方程

式中，為塔板持液量，kmol；為進(jìn)料量，kmol·h-1；為進(jìn)料中各組分摩爾分?jǐn)?shù)；為液相抽出量，kmol·h-1；為汽相抽出量，kmol·h-1；為汽相各組分摩爾分?jǐn)?shù)。

（2）總物料平衡方程

（3）熱量平衡方程

(19)

式中，F(xiàn)為進(jìn)料焓值，kJ·kmol-1；V為汽相焓值，kJ·kmol-1；為與外界的換熱量，kJ·h-1。

（4）相平衡方程

式中，為相平衡常數(shù)。

（5）摩爾分?jǐn)?shù)加和式

（6）水力學(xué)方程

(22)

式中，為塔板鼓泡區(qū)面積，m2。

（7）壓力降方程

(23)

(25)

(26)

綜上，式(14)～式(16)和式(24)～式(26)構(gòu)成了乙烯精餾塔的動(dòng)態(tài)模型，而且該模型考慮了降液管的機(jī)理過程。

3 乙烯精餾塔動(dòng)態(tài)特性分析

降液管的時(shí)滯效應(yīng)會(huì)對(duì)乙烯精餾塔的動(dòng)態(tài)特性造成影響，這是因?yàn)閷?duì)塔頂回流量進(jìn)行操作時(shí)變化的液相需要通過3個(gè)降液管后才能夠到達(dá)塔底，進(jìn)而影響塔底的溫度、組成等。塔頂對(duì)塔底操作的滯后時(shí)間為全塔降液管液相儲(chǔ)存量與全塔平均液相負(fù)荷（流量）的比值，對(duì)于塔板數(shù)較多的精餾塔，如乙烯精餾塔，滯后時(shí)間能達(dá)到20 min。

根據(jù)考慮降液管的乙烯精餾塔動(dòng)態(tài)模型，利用流程模擬軟件gPROMS對(duì)乙烯精餾塔模型式(14)～式(16)、式(24)～式(26)進(jìn)行仿真。另外，若令d()=()、d(,)=(,)且d()=()，則相當(dāng)于忽略了降液管的時(shí)滯效應(yīng)，其仿真結(jié)果可作為對(duì)比實(shí)驗(yàn)，與考慮降液管模型的仿真結(jié)果進(jìn)行比較。

圖3中，塔頂回流量增加15 kmol·h-1時(shí)，塔頂冷量增加，塔頂溫度降低，產(chǎn)品（乙烯）質(zhì)量升高。考慮降液管時(shí)，延緩了全塔響應(yīng)速度，響應(yīng)時(shí)間加長(zhǎng)。

圖4中塔頂回流量增加15 kmol·h-1時(shí)，塔底溫度降低，塔底產(chǎn)品（乙烷）質(zhì)量降低?？紤]降液管時(shí)，塔頂液相操作不能直接作用于塔底，會(huì)出現(xiàn)10～20 min的滯后，同時(shí)響應(yīng)時(shí)間變長(zhǎng)。

圖5中，再沸蒸汽量增加15 kmol·h-1時(shí)，塔頂溫度升高，產(chǎn)品（乙烯）質(zhì)量降低。考慮降液管時(shí)，延緩了全塔響應(yīng)速度，響應(yīng)時(shí)間加長(zhǎng)。

圖6中，再沸蒸汽量增加15 kmol·h-1時(shí)，塔底溫度升高，塔底產(chǎn)品（乙烷）質(zhì)量升高?？紤]降液管時(shí)，響應(yīng)時(shí)間變長(zhǎng)。

根據(jù)以上實(shí)驗(yàn)結(jié)果，可以直觀看出降液管對(duì)乙烯精餾塔的動(dòng)態(tài)特性的影響。可引入“作用距離”的概念分析以上實(shí)驗(yàn)結(jié)果。在控制系統(tǒng)中，定義“作用距離”為操作變量與被控變量之間信息傳遞的距離。在乙烯精餾塔中，作用距離可表征為

式中，R為某操作變量對(duì)被控變量的作用距離，S為最短作用距離，為自然數(shù)集，0為相鄰板之間的作用距離。

圖7為乙烯精餾塔作用距離示意圖。圖7(a)中，討論回流量對(duì)塔頂溫度的影響，可將其分解為個(gè)過程：第1個(gè)過程為作用距離Sa=0，即操作信號(hào)（回流量）直接作用于塔頂溫度；第2個(gè)過程作用距離為Sa+0，即回流量作用于第2塊塔板后間接作用于塔頂溫度；第個(gè)過程作用距離為Sa+(-1)0，即回流量作用于第塊塔板后間接作用于塔頂溫度；第-2個(gè)過程作用距離為Sa+(-3)0，即回流量作用于塔底后間接作用于塔頂溫度（除去再沸器和冷凝器后塔板數(shù)為-2）。塔頂溫度的響應(yīng)過程為這-2個(gè)過程的疊加。

對(duì)于不考慮降液管的模型，忽略了全塔氣液相流動(dòng)和傳質(zhì)的動(dòng)態(tài)過程，這-2個(gè)過程的作用距離都可認(rèn)為為零。對(duì)于考慮降液管的模型，液相通過降液管產(chǎn)生時(shí)滯，令單個(gè)降液管時(shí)滯為，則0可通過表征。若將-2個(gè)過程近似為一階慣性加純滯后環(huán)節(jié)，則圖7中4個(gè)響應(yīng)過程的傳遞函數(shù)為

(29)

(30)

降液管對(duì)精餾塔動(dòng)態(tài)特性的影響可通過式(28)～式(31)體現(xiàn)，不存在降液管時(shí)傳遞函數(shù)中所有純滯后項(xiàng)均為零，考慮降液管的時(shí)滯效應(yīng)后純滯后項(xiàng)不為零。對(duì)比式(28)～式(31)發(fā)現(xiàn)，式(29)與其他3個(gè)公式滯后項(xiàng)存在明顯差異，這是由于其最短距離不同產(chǎn)生的。對(duì)于模式(b)，塔頂回流量至少需要經(jīng)過-2塊板后才能對(duì)塔底溫度產(chǎn)生影響，即Sb=(-3)0；對(duì)于其他3種模式，滿足Sa=Sc=Sd=0。

對(duì)于式(28)～式(31)，響應(yīng)過程的調(diào)節(jié)時(shí)間見表1。表1中，取允許誤差帶寬為5%，則對(duì)于一階慣性環(huán)節(jié)其響應(yīng)時(shí)間為時(shí)間常數(shù)的3倍。分析表1可知，考慮降液管對(duì)精餾塔的動(dòng)態(tài)特性（響應(yīng)時(shí)間）產(chǎn)生了很大的影響。對(duì)于情形(a)，降液管導(dǎo)致-2個(gè)過程時(shí)滯，則總體體現(xiàn)出來的響應(yīng)過程延緩，響應(yīng)速度降低，與圖3實(shí)驗(yàn)結(jié)果匹配；對(duì)于情形(b)，降液管導(dǎo)致-2個(gè)過程時(shí)滯，而且最小時(shí)滯為(-3)，則總體體現(xiàn)出來的響應(yīng)過程初始滯后一段較長(zhǎng)時(shí)間，而且后續(xù)響應(yīng)速度降低，與圖4實(shí)驗(yàn)結(jié)果匹配；情形(c)、(d)與(a)類似。因此，在動(dòng)態(tài)建模中不能忽略降液管的時(shí)滯效應(yīng)。

表1 乙烯精餾塔動(dòng)態(tài)響應(yīng)時(shí)間

4 乙烯精餾塔控制器設(shè)計(jì)比較

根據(jù)對(duì)有、無降液管的乙烯精餾塔仿真結(jié)果的比較可知，降液管的存在使得塔頂液相的變化無法直接作用于塔底，要經(jīng)過多個(gè)降液管后才能到達(dá)塔底，這種降液管的時(shí)滯效應(yīng)直接表現(xiàn)為：操作塔頂液相時(shí)，塔底參數(shù)狀態(tài)保持一段時(shí)間（10～20 min）后對(duì)該操作做出響應(yīng)，而且由于響應(yīng)滯后的存在，整個(gè)精餾塔的動(dòng)態(tài)響應(yīng)時(shí)間變長(zhǎng)。因此，如果在建立精餾塔動(dòng)態(tài)模型的過程中忽略降液管的時(shí)滯效應(yīng)，會(huì)導(dǎo)致模型出現(xiàn)偏差，進(jìn)而降低在此基礎(chǔ)上設(shè)計(jì)的控制系統(tǒng)的可用性。

對(duì)系統(tǒng)進(jìn)行控制系統(tǒng)設(shè)計(jì)，該控制系統(tǒng)的目標(biāo)為控制塔頂乙烯和塔底乙烷產(chǎn)品質(zhì)量，由于在實(shí)際裝置中通過間接控制溫度實(shí)現(xiàn)對(duì)產(chǎn)品質(zhì)量的控制，操縱變量為回流量和再沸蒸汽量。根據(jù)圖7可知存在兩種情況的變量配對(duì)，一為模式(b)和(c)，二為模式(a)和(d)。由于模式(b)的最短作用距離最長(zhǎng)，不利于控制實(shí)現(xiàn)，第1種變量配對(duì)不可取。因此只能選擇模式(a)和(d)完成被控變量與操作變量的配對(duì)。具體控制方案見表2。

表2 乙烯精餾塔產(chǎn)品質(zhì)量控制方案

比較式(28)和式(31)中考慮與不考慮降液管的參數(shù)區(qū)別，理論上可說明控制器參數(shù)應(yīng)該存在很大的差異。另外根據(jù)表2中的控制方案，對(duì)由gPROMS實(shí)現(xiàn)的仿真過程添加控制器，并對(duì)PID參數(shù)進(jìn)行仿真整定，PID參數(shù)見表3和表4，控制效果如圖8和圖9所示。

表3 塔頂溫度控制器參數(shù)整定結(jié)果比較

表4 塔底溫度控制器參數(shù)整定結(jié)果比較

不考慮降液管的情況下，采用表3和表4中對(duì)應(yīng)參數(shù)設(shè)定塔頂和塔底溫度控制器。當(dāng)塔頂溫度出現(xiàn)0.1 K的波動(dòng)時(shí)，乙烷和乙烯產(chǎn)品質(zhì)量的調(diào)節(jié)曲線如圖8和圖9中實(shí)線所示。

考慮降液管的情況下，若采用相同的PID參數(shù)，模型計(jì)算至20 min時(shí)，由于變量超限，計(jì)算終止。調(diào)整PID參數(shù)，乙烷和乙烯產(chǎn)品質(zhì)量的調(diào)節(jié)曲線如圖8和圖9中虛線所示。

通過比較表3和表4可知，在無降液管模型和有降液管模型的基礎(chǔ)上設(shè)計(jì)的PID控制器參數(shù)差異很大，在實(shí)際應(yīng)用中，若在忽略降液管時(shí)滯效應(yīng)的條件下對(duì)實(shí)際精餾塔建模，仿真驗(yàn)證通過的控制器設(shè)計(jì)方案很可能不適用于實(shí)際生產(chǎn)。因此，在動(dòng)態(tài)建模中考慮降液管的時(shí)滯效應(yīng)具有一定的現(xiàn)實(shí)意義。

5 結(jié) 論

在精餾塔平衡級(jí)動(dòng)態(tài)建模中忽略了降液管液相流動(dòng)及傳質(zhì)的滯后，導(dǎo)致建立的模型與實(shí)際裝置存在明顯的差異。本工作針對(duì)某實(shí)際乙烯精餾塔，通過機(jī)理分析建立理論降液管模型，并將其引入到該塔的原平衡級(jí)模型中，構(gòu)建出了考慮降液管時(shí)滯效應(yīng)的乙烯精餾塔動(dòng)態(tài)模型。并對(duì)模型進(jìn)行解算，比較討論有、無降液管情況下乙烯精餾塔的動(dòng)態(tài)特性，說明忽略降液管降低了模型的準(zhǔn)確性。由于有、無降液管的差異性，對(duì)兩種動(dòng)態(tài)模型分別進(jìn)行了控制器設(shè)計(jì)，發(fā)現(xiàn)分別基于兩者得到的控制器參數(shù)差異很大。因此，在動(dòng)態(tài)建模中考慮降液管能夠更加準(zhǔn)確地對(duì)乙烯精餾塔進(jìn)行分析、控制和優(yōu)化，具有一定的現(xiàn)實(shí)意義。

References

[1] Abou-Jeyab R A, Gupta Y P, Gervais J R,. Constrained multivariable control of a distillation column using a simplified model predictive control algorithm[J]. J. Process Contr., 2001, 11(5): 509-517.

[2] Hiller C, Meier R, Niggemann G,. Distillation in special chemistry[J]. Chem. Eng. Res. Des., 2015, 99: 220-227.

[3] 黃克謹(jǐn), 錢積新, 孫優(yōu)賢, 等. 精餾塔分段正交配置建模方法[J]. 化學(xué)工程, 1994, 22(2): 18-23.Huang K J, Qian J X, Sun Y X,. Compartmental orthogonal collocation methodology for distillation modeling[J]. Chemical Engineering(China), 1994, 22(2): 18-23.

[4] 孫軍, 張貝克, 馬昕, 等. 改進(jìn)的信息流圖及其在精餾塔建模與仿真中的應(yīng)用[J]. 化工學(xué)報(bào), 2013, 64(12): 4391-4395.Sun J, Zhang B K, Ma X,. Improved information flow diagram and its application in distillation column modeling and simulation process[J]. CIESC Journal, 2013, 64(12): 4391-4395.

[5] 高翔, 劉偉, 陳海勝, 等. 一種外部熱耦合反應(yīng)蒸餾系統(tǒng)的模擬研究[J]. 化工學(xué)報(bào), 2012, 63(2): 538-544.Gao X, Liu W, Chen H S,. A simulation-based study of externally heat-integrated reactive distillation system[J]. CIESC Journal, 2012, 63(2): 538-544.

[6] de Canete J F, Gonzalez S, del Saz-Orozco P,. A harmonic balance approach to robust neural control of MIMO nonlinear processes applied to a distillation column[J]. J. Process Contr., 2010, 20(10): 1270-1277.

[7] 江愛朋, 邵之江. 乙烯生產(chǎn)流程中聯(lián)塔模擬與優(yōu)化[J]. 化工學(xué)報(bào), 2006, 57(9): 2128-2134.Jiang A P, Shao Z J. Simulation and optimization of distillation column sequence in large-scale ethylene production[J]. Journal of Chemical Industry and Engineering(China), 2006, 57(9): 2128-2134.

[8] 龔超, 余愛平, 羅祎青, 等. 完全能量耦合精餾塔的設(shè)計(jì)、模擬與優(yōu)化[J]. 化工學(xué)報(bào), 2012, 63(1): 177-184.Gong C, Yu A P, Luo Y Q,. Design, simulation and optimization of fully thermally coupled distillation column[J]. CIESC Journal, 2012, 63(1): 177-184.

[9] 薛美盛, 祁飛, 吳剛, 等. 精餾塔控制與節(jié)能優(yōu)化研究綜述[J]. 化工自動(dòng)化及儀表, 2006, 33(6): 1-6.Xue M S, Qi F, Wu G,. Survey on control and energy saving optimization of distillation column[J]. Control and Instruments in Chemical Industry, 2006, 33(6): 1-6.

[10] Gubitoso F, Pinto J M. A planning model for the optimal production of a real-world ethylene plant[J]. Chem. Eng. Process. Process Intensification, 2007, 46(11): 1141-1150.

[11] Kirshnamurthy R, Martinez H, Martinez R. Murphree point efficiencies in multicomponent system[J]. Trans. Inst. Chem. Eng., 1977, 55: 178-184.

[12] Kirshnamurthy R, Taylor R. A nonequilibrium stage model of multicomponent separation processes[J]. AIChE Journal, 1985, 31(3): 449-465.

[13] 李紅海, 姜奕. 精餾塔設(shè)備的設(shè)計(jì)與節(jié)能研究進(jìn)展[J]. 化工進(jìn)展, 2014, 33(S1): 14-18.Li H H, Jiang Y. Research progress of design and energy-saving of distillation column[J]. Chemical Industry and Engineering Progress, 2014, 33(S1): 14-18.

[14] 許鋒, 羅雄麟. 控制與工藝集成優(yōu)化設(shè)計(jì)研究進(jìn)展[J]. 化工進(jìn)展, 2005, 24(5): 483-488.Xu F, Luo X L. Recent advances in simultaneous design and control optimization[J]. Chemical Industry and Engineering Progress, 2005, 24(5): 483-488.

[15] 何仁初, 羅雄麟. 考慮降液管影響的多元精餾塔動(dòng)態(tài)機(jī)理建模. (Ⅰ):動(dòng)態(tài)機(jī)理模型的建立[J]. 石油學(xué)報(bào)(石油加工), 2007, 23(3): 58-62.He R C, Luo X L. Dynamic mechanism modeling of multicomponent distillation column by considering influence of downcomer(Ⅰ): Building of dynamic mechanism model[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2007, 23(3): 58-62.

[16] Loockett M J, Gharani A W. Downcomer hydraulics at high liquid flowrates[C]// IChemE Symp Ser No. 56. London: IChemE, 1979, 2: 343-364.

[17] Zhang M Q, Wang Z C, Zhang Y Z,. Downcomer mixing and its effect on distillation tray efficiency[J]. Transactions of Tianjin University, 1998, 4(1): 17-23.

[18] Weiler D W, Kirkpatrick R D, Lockett M J. Effect of downcomer mixing on distillation tray efficiency[J]. Chem. Eng. Progr., 1981, 77(1): 63-69.

[19] 黃雪雷, 李育敏. 板式塔弓形降液管液相流場(chǎng)CFD數(shù)值模擬[J]. 浙江工業(yè)大學(xué)學(xué)報(bào), 2010, 38(5):518-521.Huang X L, Li Y M. CFD simulation of liquid flow field in downcomer of trayed tower[J]. Journal of Zhejiang University of Technology, 2010, 38(5): 518-521.

[20] 曹振恒, 王彩琴, 趙立功, 等. 分散降液篩板塔傳質(zhì)性能[J]. 化工學(xué)報(bào), 2015, 66(10): 4061-4066.Cao Z H, Wang C Q, Zhao L G,. Mass transfer characteristics in decentralized downcomer sieve plate column[J]. CIESC Journal, 2015, 66(10): 4061-4066.

[21] Yang X, Xu Q, Li K. Safety-considered proactive flare minimization strategy under ethylene plant upsets[J]. Chem. Eng. Technol., 2011, 34(6): 893-904.

Analyzing time-lag effect of downcomer for control of ethylene column

HUANG Dong, ZHAO Minshuai, LUO Xionglin

(Department of Automation, China University of Petroleum, Beijing 102249, China)

The time-lag of downcomer caused by the influence of volume on liquid flow and mass transfer has been ignored for establishing the dynamic model of distillation column succinctly. And there is some discrepancies with commercial column. Based on the information of a commercial ethylene column, the paper set up the mechanism model of the theoretical downcomer. In addition, the improved model is the integration of the downcomer model and equilibrium stage model. The simulations of the improved model and equilibrium stage model were achieved by gPROMS, and the differences between them were shown by the comparison of the simulation results. The differences indicated that the accuracy of dynamic model by considering the influence of downcomer. Furthermore, a controller was designed for the improved model and equilibrium stage model, and there was a difference between the control parameters. It was to say, the controller designed on the foundation of ignoring downcomer may be not suit for commercial columns. Thus, considering the influence of downcomer was significant for studying the method of distillation column modeling which was the foundation of the analysis, control and optimization for distillation column.

process systems; dynamic simulation; computer simulation; distillation; downcomer

2016-08-24.

Prof. LUO Xionglin, luoxl@cup.edu.cn

10.11949/j.issn.0438-1157.20161177

TE 624

A

0438—1157（2016）11—4696—09

黃冬（1990—），男，博士研究生。

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2012CB720500）。

2016-08-24收到初稿，2016-08-31收到修改稿。

聯(lián)系人：羅雄麟。

supported by the National Basic Research Program of China (2012CB720500).