二元恒殘液組成提餾式間歇精餾的捷算法

王為國，段曉玲，曾真，王存文，吳元欣

［1.武漢工程大學化工與制藥學院，湖北武漢430074；2.綠色化工過程教育部重點實驗室（武漢工程大學），湖北武漢430074；3.武漢工程大學機電工程學院,湖北武漢430074］

二元恒殘液組成提餾式間歇精餾的捷算法

王為國1，段曉玲2，曾真3，王存文2，吳元欣2

［1.武漢工程大學化工與制藥學院，湖北武漢430074；2.綠色化工過程教育部重點實驗室（武漢工程大學），湖北武漢430074；3.武漢工程大學機電工程學院,湖北武漢430074］

探討了二元恒殘液組成提餾式間歇精餾在理想操作條件下汽化總量的數值計算方法.采用數值計算方法計算出二元恒殘液組成提餾式間歇精餾，給定分離任務理論板數與汽化總量的對應關系，并通過關聯計算結果得到其理論板數與汽化總量的近似關系式.對于給定的分離任務，如果已知汽化總量，用該近似關系式計算理論板數與精確解的最大正相對偏差為0.027 1，最大負相對偏差為-0.028 8，平均相對偏差為0.014 6；如果已知理論板數，用該近似關系式計算汽化總量與精確解的最大正相對偏差為0.041 8，最大負相對偏差為-0.035 1，平均相對偏差為0.016 5.因此，該近似關系式具有很好的計算精度.

間歇精餾；恒殘液組成；提餾式；捷算法

0 引言

間歇精餾提餾式操作方式最早由Robinson和Gilliland［1］于1950年提出，其流程與操作方式見文獻［2］.Chiotti［3-4］等、S?rensen［5］等研究均得到，當原料中輕組分含量很低的時，較常規間歇精餾，提餾式操作方式能耗相對較低的結論.許松林等［6-7］研究了提餾式間歇精餾塔的設計方法及控制方法.張雪梅［7］等推導了二元恒再沸比提餾式間歇精餾過程的最小再沸比計算方法.筆者［9-10］分別對二元常規間歇精餾恒餾出液組成操作方法的理論板與再沸器汽化總量、恒回流比操作方法的理論板與回流比進行了關聯，得到相應的捷算法.近來筆者［2］推導出二元恒殘液組成提餾式間歇精餾，理想操作條件下最小汽化總量的計算公式，為其捷算法提出提供可能.本研究通過對二元恒殘液組成提餾式間歇精餾的理論板與塔釜汽化總量進行關聯，得到其捷算法.

1 再沸器汽化總量的數值計算方法

二元恒殘液組成提餾式間歇精餾的流程見文獻［2］圖1.操作條件同文獻［2］，本研究稱之為理想操作條件.基本方程見文獻［2］式（1）～（8），式（20）～（23）.理想操作條件下二元恒殘液組成提餾式間歇精餾，在對塔頂貯槽內液體中輕組分濃度由xF變為xDe的精餾過程離散后，可用試差圖解法求解出給定的分離任務理論板數與汽化總量的對應關系，但因其計算量非常大，用試差圖解法大規模的計算幾乎不可能.以下簡要介紹用圖解法原理，編程計算的方法.本研究所介紹的數值計算方法，主要滿足本文關聯N+1與nV的需要，其步驟及相關的（離散化）方程如下：

①輸入：α、nF、xF、xW、ηh、N+1和k.

②計算xDe.

由文獻［2］式（8）、（2）、（3）得

③i=0.

④確定x(i)D.

⑤計算R(i)b.

計算R(i)b的過程較為復雜，后文將詳細說明.

⑥判斷i+1＞k是否成立.

如果不成立，i=i+1，轉到④

如果成立，計算nV，輸出計算結果.

由文獻［2］式（23）得

上述計算過程步驟⑤是計算已知α、xD、xW和N+1時的Rb，具體過程（見圖1示意，作圖1的條件為：α=2.5、xD=0.50、xW=0.10、N+1=6，試差作出的操作線為-ba，此時，y1=0.68， Rb=2.22.）為：假定Rb的數值，從xW開始，交替使用文獻［2］式（1）與式（4）（即：逐板計算法）計算出x0.如果x0＞xD表明假定Rb的數值過大（圖中的線-ba"所示）.如果x0＜xD表明假定Rb的數值過小（圖中的線-ba'所示）.如果x0=xD，表明假定Rb的數值就是所要求的（如：圖中的線-ba所示）.由圖1可見，對于正常的精餾過程，當α、N+1和xW一定時，x0是Rb單調增函數，因此，數值法計算已知α、xD、xW和N+1時的Rb較為合理方法為：采用Newton法修正Rb迭代求解.

圖1 Rb求解示意圖Fig.1 schematic diagram of solving Rb

2 最少理論板數

二元恒殘液組成提餾式間歇精餾過程中，再沸器汽化速率V為有限值，只有在塔底出料液速率W→0時，才有Rb=V W→∞.當Rb=V W→∞時，N+1→Nm+1.由塔頂貯槽內液體中輕組分濃度由xF變為xDe的精餾過程中，xW不變，xD在增大，所需的最少理論板數是變化的，瞬時最少理論板數由芬斯克方程得

精餾起始時，所需的最少理論板數最小為

精餾終了時，所需的最少理論板數最大為

3 最小汽化總量

二元恒殘液組成提餾式間歇精餾，當N+1→∞時，nV→nVm.塔頂貯槽內液體中輕組分濃度由xF變為xDe的精餾過程中，塔釜汽化總量用下式計算，見文獻［2］式（25）.

4 捷算法

獲得二元恒殘液組成提餾式間歇精餾N+1與nV復雜對應關系的近似計算式的具體步驟如下：

（1）確定分離任務對在理想操作條件下，確定二元理想混合物間歇精餾分離任務的4個獨立變量(α,xF,xW,ηh)用5水平正交設計方案取值給定25個分離任務（見表1）.

（2）由式（8）計算的各分離任務的nVm/nF.

（3）由式（6）、（7）計算的各分離任務的Nms+1、Nme+1.

（4）按N+1≈(1.3～1.5)(Nme+1)整數，選定各分離任務的N+1，用數值計算法計算對應的nV/nF.

（5）關聯.對對應的數據進行變量組合，令

關聯對應的X與Y得

由表1、圖2可見，式（11）很好地表達了Y與

X的近似關系，如果已知X，用式（11）計算Ya與精確解的最大正偏差（Y-Ya）為0.018 1［見表（1）編號3分離任務］，最大負偏差為-0.015 9［見表（1）編號21分離任務］，平均偏差0.008 4.對于給定的分離任務，如果已知nVnF，用式（11）計算Na+1與精確解的最大正相對偏差［(Na-N)(N+1)］為0.027 1［見表（1）編號21分離任務］，最大負相對偏差為-0.028 8［見表（1）編號3分離任務］，平均相對偏差為0.014 6（見圖3）；如果已知N+1，用式（11）計算nVa與精確解的最大正相對偏差［(nVa-nV)nV］為0.041 8［見表（1）編號21分離任務］，最大負相對偏差為-0.035 1［見表（1）編號3分離任務］，平均相對偏差為0.016 5（見圖4）.此外，必須說明的是，關聯出式（11）中N+1的取值范圍為N+1≈(1.3～1.5)(Nme+1)，應用式（11）計算，如果N+1的值過大（或過小），均可能引起較大的偏差.

表1 N+1和nV的關聯數據表Table 1 Data of correlation N+1 w ith nV

圖2 Y（Ya）和X的關系Fig.2 Relationship of Y（Ya）and X

圖3 Na+1和N+1的偏差Fig.3 Deviation of Na+1 from N+1

圖4 nVanF和nVnF的偏差Fig.4 Deviation of nVanFfrom nVnF

5 結語

a.介紹了二元恒殘液組成提餾式間歇精餾，理想操作條件下，汽化總量的數值計算方法.

b.關聯了恒殘液組成二元提餾式間歇精餾時，給定的分離任務的理論板數與汽化總量，得其近似關系式（式（11）），該近似關系式具有很好的計算精度.

符號說明

Nm+1——瞬時最少理論板數（包括塔釜）；

Nms+1——精餾起始時的最少理論板數（包括塔釜）；

Nme+1——精餾終了時的最少理論板數（包括塔釜）；

nF——原料量，kmol；

nDe——塔頂產品量，kmol；

nV——再沸器汽化總量，kmol；

nVm——再沸器最小汽化總量，kmol；

nWe——塔底產品量，kmol；

Rb——再沸比；

V——塔內汽相流率，kmol·h-1；

W——塔底出料液速率，kmol·h-1；

x——塔板上輕組分的摩爾分數；

xD——塔頂貯槽中輕組分的瞬時摩爾分數；

xF——原料中輕組分的摩爾分數；

xW——塔底出料中輕組分的摩爾分數；

α——相對揮發度；

ηh——塔底產品中重組分的收率；

上角標

a——近似計算結果；

(i)——數值計算時，各物理量第(i)段對應的數值.

致謝

感謝武漢工程大學化工與制藥學院提供的實驗平臺.

［1］Robinson C S，Gilliland E R.Elements of fractional Distillation［M］.4th ed.New York：M cGraw-Hill，1950.

［2］王為國，曾真，覃遠航，等.恒再沸比提餾式間歇精餾的最小再沸流比與能耗分析［J］.化工學報，2012，63（7）：2106-2112.

WANGWei-guo，ZENG Zhen，QIN Yuan-hang，etal. Energy consumption analysis and m inimum reboil ratio of stripping batch distillation under constant reboil ratio［J］.Journal of Chem ical Industry and Engineering，2012，63（7）：2106-2112.（in Chinese）

［3］Chiotti O J，Iribarren O A.Simplified models for binary batch distillation［J］.Computers&Chem ical Engineering，1991，15（1）：1-5.

［4］Chiotti O J，Salomone H E，Iribarren O A.Selection of multicomponent batch distillation sequences［J］. Chem ical Engineering Communications 1993，119（1）：1-21.

［5］S?rensen E，Skogestad S．Com parison of regular and inverted batch distillation［J］.Chem icalEngineering Science，1996，51（22）：4949-4952.

［6］XU Song-lin，Salomone H E，Iribarren O A．Shortcut procedure for inverted batch distillation column（I）multicom ponent ideal system［J］.Chinese Journal of Chem ical Engineering，2001，9（1）：28-33.

［7］Xu S L，Espinosa J，Salomone H E，et al.Operation of a batch stripping distillation column［J］.Chinese Journal of Chem ical Engineering，2001，9（2）：141-144.

［8］張雪梅，郭錦棠，張衛江.二元物系提餾式分批精餾過程［J］.大津大學學報，2005，38（11）：950-954.

ZHANG Xue-mei，GUO Jin-tang，ZHANG Wei-jian．Binary Stripping Batch Distillation［J］.Journal of Tianjin University，2005，38（11）：950-954.（in Chinese）

［9］王為國，曾真，張所信，等.二元混合物恒餾出液組成間歇精餾的捷算法［J］.武漢化工學院學報，1999，21（3）：5-9.

WANGWei-guo，ZENG Zhen，ZHANG Suo-xin，et al. Shortcut method for batch distillation of binary m ixture under constant distillate com position［J］. Journal of Wuhan Institute of Chem ical Technology，1999，21（3）：5-9.（in Chinese）

［10］王為國，曾真，畢亞凡，等.二元混合物恒回流比間歇精餾的捷算法［J］.武漢化工學院學報，2000，22（1）：16-19.

WANG Wei-guo，ZENG Zhen，BI Ya-fan，et al. Shortcut method for batch distillation of binary m ixture under constant reflex ratio［J］.Journal of Wuhan Institute of Chem ical Technology，2000，22（1）：16-19.（in Chinese）

Shortcut method of binary stripping batch distillation under constantresidual fraction composition

WANG Wei-guo1，DUAN Xiao-ling2，ZENG Zhen3，WANG Cun-wen2，WU Yuan-xin2
［1.School of Chemical Engineering and Pharmacy，Wuhan Institute of Technology，Wuhan 430074，China；2. Key Laboratory for Green Chemical Process of Ministry of Education(Wuhan Institute of Technology)，Wuhan 430074，China；3.School of Mechanical and Electrical Engineering，Wuhan Institute of Technology，Wuhan 430074，China］

The numerical method for calculating total vapor capacity of binary stripping batch distillationunder constant residual fraction composition and ideal operation conditions was explored. The correspondingrelationship between the theoretical trays and the total vapor capacity of binary stripping batch distillationunder constant residual fraction composition was figured out by the numerical method for given separationtasks，and an approximate relationship between the theoretical trays and the total vapor capacity was set upthrough correlating the calculated results. For a given separation task，if the total vapor capacity was known，the maximum positive relative deviation and the maximum negative relative deviation between the calculatedtheoretical trays based on the approximate relationship and the exact solution were 0.027 1 and -0.028 8，respectively，and the average relative deviation was 0.014 6. If the theoretical trays were known，themaximum positive relative deviation and the maximum negative relative deviation between the calculatedtheoretical trays based on the approximate relationship and the exact solution were 0.041 8 and -0.035 1，respectively，and the average relative deviation was 0.016 5. Therefore，the approximate relationship has agood calculation accuracy.

batch distillation；constant residual fraction composition；stripping；shortcut method

TQ015

A

10.3969/j.issn.1674-2869.2014.01.003

1674－2869（2014）01－0009－05

本文編輯：張瑞

2013-06-03

王為國（1963-），男，江蘇建湖人，副教授，碩士.研究方向：間歇精餾基礎理論.