Kaibel分壁精馏塔分离芳烃的稳态和动态模拟

2015年4月 CIESC Journal

April 2015

第66卷 第4期 化 工 学 报 V ol.66 No.4

Kaibel 分壁精馏塔分离芳烃的稳态和动态模拟

蔺锡钰，吴昊，沈本贤，凌昊

（华东理工大学化学工程联合国家重点实验室，上海200237）

摘要：Kaibel 分壁精馏塔（Kaibel divided-wall column ，KDWC ）可在一个塔内实现四组分混合物的高纯度分离。本文以分离苯、甲苯、二甲苯和均三甲苯为研究对象，建立了KDWC 严格稳态模型，获得了优化的塔体结构。通过研究KDWC 分离17组进料组成的基础上，获得了KDWC 的稳态分离的初步规律：实现KDWC 的高纯度分离，既需要通过分液比控制预分馏段顶部馏出气相中的二甲苯含量，又需要通过分气比控制预分馏底部馏出液相中甲苯含量；两个侧线的组成中，重组分杂质的含量要远多于轻组分杂质的含量；中间组分甲苯在预分馏段仍有返混。随后，在Aspen Dynamic 环境下建立了KDWC 的组分控制模型，控制结果表明该模型可以应对±10%的流量和进料组成波动，但二甲苯产品纯度会出现少量偏差。

关键词：Kaibel 分壁精馏塔；芳烃；蒸馏；分离；稳态；动态控制；模拟 DOI ：10.11949/j.issn.0438-1157.20141611

中图分类号：TQ 202 文献标志码：A 文章编号：0438—1157（2015）04—1353—10

Steady-state behavior and control of Kaibel divided-wall column

for aromatics separation

LIN Xiyu, WU Hao, SHEN Benxian, LING Hao

(State Key Laboratory of Chemical Engineering , East China University of Science and Technology , Shanghai 200237, China )

Abstract : Kaibel divided-wall column (KDWC) permits separation of a four-component mixture into four pure fractions in a divided-wall column. This work studied separation of a four-component mixture of benzene, toluene, o -xylene and tri-methyl-benzene in a KDWC. The optimum economic design of the Kaibel divided-wall column was obtained. Then seventeen cases of various feed compositions were investigated for obtaining the steady state behavior of the KDWC. Xylene should be controlled at the top of the prefractionator, while toluene should be controlled at the bottom of the prefractionator; the content of heavy impurities was much larger than that of the light impurities in the two side streams; the intermediate component, toluene re-mixed in the prefractionator. At last, a control structure with five composition control loops was established. Simulation showed that it could handle ±10% disturbances of flow rate and feed composition, but the purity of xylene sometimes had small deviations.

Key words : Kaibel divided-wall column; aromatics; distillation; separation; steady-state; dynamic control; simulation

2014-10-27收到初稿，2014-12-31收到修改稿。

联系人：凌昊。第一作者：蔺锡钰（1989—），男，硕士研究生。基金项目：国家自然科学基金项目（21476081)；中央高校基本科研业务费专项资助项目；上海市教育委员会科研创新重点项目（14ZZ058）。

Received date : 2014-10-27.

Corresponding author : LING Hao, linghao@https://www.360docs.net/doc/ab11864710.html,

Foundation item : supported by the National Natural Science Foundation of China (21476081), the Fundamental Research Funds for the Central Universities of Ministry of China and Shanghai Municipal Education Commission (14ZZ058).

化工学报第66卷·1354·

引言

分壁式精馏塔（divided-wall column, DWC）也叫全热耦合精馏塔，可以高纯度地分离三组分和四组分混合物，相比常规两塔和三塔的分离方式可以节省能耗30%以上[1]。在常规的两塔分离三组分混合物（A, B, C）的直接序列（DS）过程中，轻组分A在第1个精馏塔塔顶馏出，中间组分B和重组分C在塔底馏出。在多数三组分混合物的分离过程中，中间组分B会在第1个精馏塔靠近塔釜的塔板上浓度达到最大值，然后逐渐下降。这一热力学损失现象在常规的精馏塔分离三组分混合物的过程中无法避免，DWC是解决这一问题的良好途径[2-4]。

目前，关于分离三组分的分壁式精馏塔的研究非常广泛，国外工业应用也趋于成熟。Skogestad等首次对DWC的动态控制进行了研究，提出了以分液比作为调控变量的四点控制法[5]；通过研究恒定相对挥发度的虚拟系统，发现在预分馏塔中应尽量避免重组分从塔顶移出[6]。Abdul等[7]在研究相同体系中首次使用温度代替组分作为控制变量，提出温度控制结构。Serra等[8]研究了3个具有不同相对挥发度的系统，比较PID和动态矩阵控制，发现PID 控制能提供更好的动态控制。Ling等[9-11]以苯、甲苯和二甲苯三元体系作为研究对象，使用Aspen Plus构建严格精馏模型，利用4个精馏模块模拟DWC，分别采用组分控制、温度控制和温差控制结构对DWC进行测试，发现组分控制可以很好地应对组分组成和进料流量组成波动。许锡恩等[12]对三相热耦合精馏过程进行了模拟，并对具有双液相的精馏过程的数学模型予以修正。袁希钢等[13-14]针对热耦合精馏所需塔段数目以及冷凝器和再沸器数目的不确定性，提出了分解求解策略和编码表达法，使流程结构的描述和优化变得更加简便。曲平等[15]以丁二烯分离装置为对象，对热耦合精馏用于非理想的操作特性和节能效果进行了模拟和研究。孙兰义等对分壁精馏塔反应精馏水解醋酸甲酯的控制及空分应用性能做了详细的研究[16-18]，最近在反应精馏塔的设计和应用取得了很多成绩[19-20]。裘兆蓉等[21-23]对分壁精馏塔开展了大量的稳态实验和模拟研究，采用分壁精馏塔分离裂解汽油、萃取精馏制无水乙醇和无水叔丁醇均取得了满意的效果。Huang 等[24-26]近几年对分壁精馏塔操作过程的“黑洞”现象、温度和温差过程控制也开展了大量研究，取得了诸多成果。杨剑等[27]通过串联常规塔和分壁精馏塔实现了芳烃四组分混合物的高纯度分离，并实现了组分控制结构对进料流量和组成10%波动的平稳控制。

虽然上述方法很好地解决了三组分分壁精馏塔分离过程的稳态和动态问题，然而针对更为复杂的四元混合物的分壁精馏塔精馏体系的工作尚需进一步研究。1987年Kaibel[28]提出在单个分壁精馏塔（Kaibel divided-wall column, KDWC）中可以实现四组分混合物的清晰分割，但在随后的20多年来很少有文献报道。Skogestad课题组[29-32]首次采用Matlab研究了非严格精馏模型的四组分分壁精馏塔的分离过程的稳态和动态过程，开创了此方面学术研究的先河。KDWC用于设计和操作使用的有5个自由度，分别是：回流量R，再沸器负荷Q R，上侧线流量S1，下侧线流量S2和分液比βL。这5个可以控制的自由度可以用来操控5个变量，其中4个产品的浓度是必须要被控制的，而第5个自由度βL 可以实现能耗优化目的。通过对甲醇、乙醇、丙醇和正丁醇的模拟和24块理论板（高8 m，直径50 mm）的小试实验，他们发现4×4的PID温度控制回路的各温度点和Matlab模拟温度接近，实现了KDWC的平稳控制，但侧线产品纯度最高仅达到88%（摩尔分数），无法达到高纯度分离的要求。在模拟方面，该课题组研究发现分气比βV对于能耗的影响是十分敏感的，在他们的研究中，常用βV代替βL作为一个潜在的自由度。Dejanovi?等[33]对KDWC分离多组分混合物进行模拟研究，将十五组分芳烃混合物分离成4种给定的产品组成，模拟结果显示可以节省年总投资（TAC）40.1%。此外，Akzo Nobel公司的Kiss等[34]进行了正戊烷-苯-甲苯-二甲苯-均三甲苯体系的分离研究，分离目标为99.5%（质量分数）以上。

本文研究KDWC分离苯-甲苯-二甲苯-均三甲苯（BTXH；benzene，B；toluene，T；o-xylene，X；1,3,5-trimethylbenzene，H）四组分混合物的稳态和动态分离过程。首先，通过Aspen Plus中的RadFrac 模块建立了KDWC的稳态严格精馏模型，通过经济优化获得了优化的塔体结构；随后考察了17种进料组成、不同βV和βL以及KDWC内部关键塔板的气液相组成的关系，获得了KDWC的稳态分离的初步规律；最后建立了KDWC的浓度控制模型，考察±10%的流量和进料组成波动时控制效果。

第4期蔺锡钰等：Kaibel分壁精馏塔分离芳烃的稳态和动态模拟·1355·

1 稳态模拟

1.1稳态设计

以BTXH四组分混合物作为研究对象，其沸点分别为353、385、419和437 K，相对挥发度αB/αT/αX/αH约为17.0/6.51/2.76/1。进料流量为 1 kmol·s?1，进料温度为351 K，四组分混合物等摩尔进料。产品纯度要求为99%（摩尔分数）。物性方法为Chao-Seader。采用Aspen Plus中的RadFrac 模块，KDWC是利用一个提馏塔（只有一个再沸器），两个平行的吸收塔（无再沸器和冷凝器）和一个精馏塔（只有一个冷凝器）搭建而成。与以往分离三组分混合物的DWC相比，KDWC多了一个侧线抽出。

图1是以最小年总投资（TAC）为目标优化得到的最优塔体结构。年总投资（TAC）同时考虑到了能耗与设备投资这两方面因素[35]。可以用来优化年总投资（TAC）的变量有：四段塔的塔板数、进料位置、两条侧线的出料位置、分液比βL和分气比βV。在KDWC中，隔板左侧的预分馏段塔板数与隔板右侧的侧线出料段塔板数相同。

图1中模拟塔顶压力为101325 Pa（1 atm），板压降设定为689 Pa（0.0068 atm）。优化后的KDWC 共83块理论板。在Aspen Plus中，冷凝器为第1块塔板，再沸器为第83块塔板。KDWC精馏段是从第2块塔板到第11块塔板。预分馏段和侧线出料段是从第12块塔板到第59块塔板。提馏段是从第60块塔板到第82块塔板。进料位置在隔板预分馏侧第30块板，上下侧线的出料位置分别在隔板侧线出料侧第19、第34块板。操作参数为：回流比6.67，再沸器负荷69.93 MW，βL为0.257，βV为0.586。全塔直径为10.8 m。而用常规精馏方法分离四组分至少需要3个精馏塔，有两种序列，即直接序列（DS）和间接序列（IS），本文使用直接序列分离流程与KDWC进行比较，流程如图2所示。表1给出KDWC结构与常规三塔结构能耗和经济对比

图1 KDWC的稳态设计Fig.1 Steady state design of KDWC

化 工 学 报 第66卷

·1356·

图2 直接序列常规三塔结构

Fig.2 Direct-separation sequence of conventional

three-column structure

数据。从表1中可以看出，KDWC 较常规三塔流程

可减小再沸器负荷Q R 28.5%，可减小冷凝器负荷Q D 30.6%，可降低TAC 约33.7%。 1.2 稳态分析

图3为KDWC 的全塔塔板温度变化曲线，其中虚线为主塔塔板温度，实线为预分馏段塔板温度。从图3中可以看出，主塔的塔板温度在第11板至第48板时温度明显高于预分馏段塔板温度。其中在第30块板附近达到23℃左右，这一现象说明分隔壁两端的温差明显，分隔壁在实际设计时应考虑避免温差过大导致分隔壁两端的传热。

图4为KDWC 的液相组成分布。从图4 (a)可以看出：①甲苯的浓度变化较为复杂，实线的预分馏段中甲苯浓度x P,T 从第12块塔板（预分馏段的

图3 KDWC 的温度分布曲线 Fig.3

Temperature profiles in KDWC

图4 KDWC 液相分布 Fig.4 Composition profile in KDWC

第1块塔板）到第29块塔板浓度逐渐降低，随后到

第33块塔板逐渐升高，然后逐渐降低，在第59块板时基本消失。在预分馏塔中x P,T 浓度的变化自上而下呈现了一个从高→低→次高→消失的趋势，这一趋势为典型的返混现象，热力学角度上导致了能量损失。此外，在预分馏塔底部（第59块板）甲苯基本消失，这就避免了甲苯从主塔的提馏段（60～82）进入分隔壁的另一侧，保证了二甲苯的产品纯度。②苯的浓度从塔顶到第40块板时呈持续下降趋势，在主塔中苯浓度x M,B 接近第20块板时已很低，预分馏段苯浓度x P,B 接近第40块板时接近0，无返混现象。

从图4 (b)可以看出：①二甲苯的浓度变化也较为复杂，预分馏段中二甲苯浓度x P,X 从第59块塔板（预分馏段的最后1块塔板）到第29块塔板浓度逐渐降低，随后到第24块塔板逐渐升高，然后逐渐降低，在第12块板时基本消失。在预分馏塔中x P,X 浓度的变化自下而上呈现了一个从高→低→次高→消

表1 Kaibel 分壁精馏塔结构与常规三塔结构比较

Table 1 Comparison of Kaibel divided-wall column structure and conventional DS three-column structure

Conventional three-column structure

Item C1 C2 C3 Total Kaibel divided-wall column

structure Relative gain/%

N

30 49 47 126

83 — D /m 6.96 6.34 7.72 — 10.8 — Q R /MW 31.69 26.70 39.41 97.80 69.93 28.5 Q D /MW 20.89 25.13 39.89 85.91 59.62 30.6 TAC ×10?6/$·a ?1

5.83 5.29 7.51 18.63

12.36

33.7

Note: relative gain =1?(Kaibel divided-wall column value/conventional three-column value).

第4期蔺锡钰等：Kaibel分壁精馏塔分离芳烃的稳态和动态模拟·1357·

失的趋势，也存在导致能量损失的返混现象。此外，在预分馏塔顶部（第12块板）x P,X含量很低，这就避免了二甲苯从主塔的精馏段（2～11）进入分隔壁的另一侧，保证了甲苯的产品纯度。②均三甲苯的浓度从塔底到第20块板时呈持续下降趋势，在主塔中均三甲苯浓度x M,H接近第40块板时已很低，预分馏段均三甲苯浓度x P,H接近第20块板时接近0，无返混现象。

从上述现象可以看出：KDWC预分馏段塔顶气相的重组分（均三甲苯）含量与在预分馏段塔底液相的轻组分（苯）含量均可忽略不计，因此，预分馏段的作用是防止二甲苯通过预分馏段塔顶进入侧线出料段，除此以外，防止甲苯通过预分馏段塔底进入侧线出料段。而在模拟过程中发现，预分馏段塔底甲苯的含量也非常低（摩尔分数2.88×10?5%），这意味着只需利用分液比控制二甲苯在预分馏段塔顶的含量，就可以同时实现预分馏段的两个主要作用。此外，预分馏段存在甲苯和二甲苯的返混现象，通过增加全塔的塔板数进一步研究发现甲苯的返混现象无法避免，而二甲苯的返混现象是可以避免的。其原因尚需进一步考察。

在Kaibel精馏塔中，分液比对能耗有很显著的影响，合理利用和控制分液比非常重要。图5为βL 与能耗及预分馏段塔顶气相中二甲苯含量y P12（摩尔分数）的组成关系。从图5可以看出，控制y P12（摩尔分数）可以实现进料组成波动时的能耗最小化。即当进料组成发生波动时，为了使能耗最小化，βL需要自发进行调整以保证y P12（摩尔分数）在0.008%附近。图5 (a)为进料中苯浓度波动时，能耗随分液比的变化情况。在进料组成未发生波动前，苯、甲苯、二甲苯和均三甲苯的含量分别为25%（摩尔分数）。进料组成变化以苯为例，苯含量增加10%意味着苯含量从25%提升至27.5%。苯含量增加20%意味着苯含量从25%提升至30%。而其他3个组分的含量随之发生变化，但相互之间的比例保持不变，最终使四者总和为100%。图5 (a)上图中，用星号标注的为最优βL对应的最小能耗，从图中可以发现，能耗随着进料组成的变化而变化，苯含量增加分离能耗逐渐降低。图5 (a)下图显示了进料中苯浓度波动时，y P12（摩尔分数）随βL的变化情况。从上下两图中可以明显地看出，预分馏段塔顶二甲苯含量为

P12,X

y=0.0081所对应的βL，基本上与最优能耗所对应的βL相同。图5 (b)为进料甲苯浓度发生变化的情况，y P12,X的变化规律与进料组成中苯含量

图5 分液比对能耗及预分馏段塔顶气相组成的影响Fig.5 Effect of liquid split ratio on energy consumption and prefractionator vapor composition

变化类似，但体系的能耗随着甲苯含量的增加呈现先下降后增加的趋势。二甲苯和均三甲苯的变化也与苯和甲苯变化相似，文中省略相关图表。

2 组分浓度控制结构

当进料组成和进料流量发生波动时，控制的目的是保持产品纯度的同时实现最优能耗操作。本文采用传统PID控制结构，选择控制杂质的含量，即控制甲苯在塔顶的含量为

DT

x=0.01，控制二甲苯在

塔底的含量为

BX

x=0.01，控制二甲苯在上侧线的含

量为

S1X

x=0.008，控制均三甲苯在下侧线的含量为S2H

x=0.009。

4个被控变量需要4个控制参数，有5个可用的控制参数：R、Q R、S1、S2和βL。如前所述，分液比

化 工 学 报 第66卷 ·1358· 可用来控制y P12（摩尔分数）从而达到能量优化的目的，因此需要建立βL →y P12,X 的控制回路，即通过调节βL 控制y P12保持在0.008附近。其余4个控制参数可以用来控制x DT 、x BX 、x S1X 和x S2H 。根据分壁

精馏塔分离三组分的经验，采用R →x DT 、Q R →x BX 、S 1→x S1X 和S 2→x S2H 进行控制参数与控制变量配对。

按照上述控制思想在Aspen Dynamic 环境中进行控制结构设计，图6和图7分别为控制流程图和

图6 组分控制回路

Fig.6

Aspen Dynamic implementation of composition control

图7 组分控制控制器面板

Fig.7 Controller faceplates of composition control

第4期 蔺锡钰等：Kaibel 分壁精馏塔分离芳烃的稳态和动态模拟

·1359·

表2 组分控制回路中控制器调谐参数 Table 2 Controller tuning parameters of

composition control

Control loop Controlled variable Manipulated variable

Controller gain K C Controller integral time τI /min CC0 y P12X

βL

0.18 47.52

CC1 x DT R 0.07 89.76 CC2 x S1X S 1 0.52 69.96 CC3 x S2H S 2 0.44 50.16 CC4

x BX

Q R

0.15 76.56

控制面板。5个组分控制回路均有一个5 min 时滞，

通过以下的调谐方法进行调谐。首先，进行Q R →x BX 回路的调谐，其余4个控制器调为手动模式，利用中继-反馈测试找到最终的增益和周期，使用Tyreus-Luyben 的调谐方法进行计算。随后，将Q R →x BX 回路调为自动模式，其余仍为手动模式，用相同的方法进行R →x DT 回路的调谐。然后，将R →x DT

回路和Q R →x BX 回路调为自动，对S 1→x S1X 回路和

S 2→x S2X 回路的调谐，之前调谐过的3个回路调为自动模式。最后，完成βL →y P12,x 回路调谐，表2为5个控制回路的调谐参数。

3 控制结果

图8～图12为组分控制结构在进料流量和进料组成发生±10%波动时的动态控制结果。以图8为例，进料流量增加或减少10%从2 h 开始，从左到右分别为B 、T 、X 和H 的变化情况，图中实线和虚线分别为流量减少10%和增加10%的控制结果。从图中可以看出，采用的控制方案可以实现产物纯度的平稳控制，与其他3个组分相比，二甲苯X 的浓度波动范围最大，重组分H 的波动范围最小。所有产品在10 h 左右纯度达到设定值，说明组分控制结构可以很好地控制进料流量±10%的波动。

图9为进料组成中苯发生±10%波动时的动态

图8 进料流量波动±10%的控制结果

Fig.8 Control results of ±

10% feed flow rate disturbance

图9 苯组分含量波动±10%的控制结果 Fig.9 Control results of ±10% benzene disturbances

化 工 学 报 第66卷

·1360· 控制结果。从图中可以看出，B 、T 和H 均得到很好的控制，产品纯度可较快回归到99%附近；S 2中的二甲苯的纯度经过50 h ，最终稳定在98.7%（摩尔分数），略有偏差。图10为进料组成中甲苯发 生±10%波动时的动态控制结果，与图9类似，当T ?10%时4

个产品的纯度均能得到很好的控制；当

图10 甲苯组分含量波动±10%的控制结果

Fig.10 Control results of ±

10% toluene disturbances

图11 二甲苯组分含量波动±

10%的控制结果

Fig.11 Control results of ±10% o -xylene disturbances

图12 均三甲苯组分含量波动±10%的控制结果 Fig.12 Control results of ±10% heavy composition disturbances

第4期蔺锡钰等：Kaibel分壁精馏塔分离芳烃的稳态和动态模拟·1361·

T+10%时S2中的二甲苯经过20 h波动，最终稳定在98.1%（摩尔分数）。图11为进料组成中二甲苯发生±10%波动时的动态控制结果，当进料X-10%时S2中的二甲苯含量经过16 h波动，最终稳定在97.5%（摩尔分数）。图12为进料组成中重组分发生±10%波动时的动态控制结果，所有产品的纯度均能得到良好控制。之所以二甲苯的含量在控制过程中出现少量偏差，其原因在于KDWC与传统的三组分DWC相比自由度更多，两个侧线相互耦合，控制难度增大。若要实现更好的控制效果，还需进一步深入研究。

4 结论

分离芳烃四组分BTXH的KDWC稳态模拟结果表明：优化的塔体结构下预分馏段存在二甲苯组分返混；实现四组分的清晰分割必须控制预分馏段顶部馏出的二甲苯含量和底部的甲苯含量。不同进料组成稳态模拟结果表明：控制预分馏段顶部馏出的二甲苯含量可实现KDWC的能耗优化的目标。基于上述分析，在Aspen Dynamic环境下建立了KDWC的5×5组分浓度控制结构，控制结果表明该结构可很好地应对进料流量和组成发生±10%波动，二甲苯产品纯度有时会略有偏差。

符 号 说 明

A C——塔顶冷凝器换热面积，m2

A R——塔釜再沸器换热面积，m2

D——塔的直径，m

N——塔板数

N F——进料板位置

N S——出料板位置

Q D——塔顶冷凝器负荷，MW

Q R——塔釜再沸器负荷，MW

R——分壁精馏塔回流比

S1——上侧线流量，kmol·s?1

S2——下侧线流量，kmol·s?1

α——相对挥发度

βL——分壁精馏塔分液比

βV——分壁精馏塔分气比

下角标

B ——苯

H ——均三甲苯

M ——主塔

P ——预分馏段

T——甲苯

X——二甲苯

References

[1] Stankiewicz A, Moulijn J A. Process intensification [J]. Ind. Eng.

Chem. Res., 2002, 41: 1920-1924

[2] Hernandez S, Pereira-Pech S, Jimenez A, Rico-Ramírez V. Energy

efficiency of an indirect, thermally coupled distillation sequence [J].

Can. J. Chem. Eng., 2003, 81: 1087-1091

[3] Lee J Y, Kim Y H, Hwang K S. Application of a fully thermally

coupled distillation column for fractionation process in naphtha reforming plant [J]. Chem. Eng. Process.: Process Intensif., 2004, 43:

495-501

[4] Ling H, Cai Z, Wu H, Wang J, Shen B. Remixing control for

divided-wall column [J]. Ind. Eng. Chem. Res., 2011, 50: 12694-12705 [5] Wollf E A, Skogestad S. Operation of integrated three-product

distillation columns [J]. Ind. Eng. Chem. Res., 1995, 34: 2094-2103 [6] Halvorsen I J, Skogestad S. Optimal operation of Petlyuk distillation:

steady-state behavior [J]. Process Control, 1999, 9: 407–424

[7] Abdul Mutalib M I, Zeglam A O, Smith R. Operation and control of

dividing wall columns (Ⅱ): Simulation and pilot plant studies using temperature control [J]. Trans. Inst. Chem. Eng. Part A, 1998, 76:

319-334

[8] Serra M, Perrier, Esuna A, Puigjaner L. Analysis of different control

possibilities for the divided wall column: feedback diagonal and dynamic matrix control [J]. Comput. Chem. Eng., 2001, 25: 859–866 [9] Ling H, Luyben W L. New control structure for divided-wall columns

[J]. Ind. Eng. Chem. Res., 2009, 48: 6034-6049

[10] Ling H, Luyben W L. Temperature control of the BTX divided-wall

column [J]. Ind. Eng. Chem. Res., 2010, 49: 189-203

[11] Ling H, Cai Z, Wu H, Wang J, Shen B J. Remixing control for

divided-wall columns [J]. Ind. Eng. Chem. Res., 2011, 50 (22): 12694-12705

[12] Xu Xi’en (许锡恩), Li Shouchun (李守春), Wang Baoguo (王保国),

Ma Haihong (马海洪). Simulation of thermally coupled distillation processes with two liquid phases [J]. Journal of Chemical Industry and Engineering (China) (化工学报), 1989, 40 (1): 28-37

[13] An Weizhong (安维中), Y uan Xigang (袁希钢). Synthesis of

thermally coupled complex distillation system based on stochastic optimization (Ⅰ): Problem formulation [J]. Journal of Chemical Industry and Engineering(China) (化工学报), 2006, 57 (7): 1591-1598

[14] An Weizhong (安维中), Y uan Xigang (袁希钢). Synthesis of

thermally coupled complex distillation system based on stochastic optimization (Ⅱ): Example illustration and analysis [J]. Journal of Chemical Industry and Engineering (China) (化工学报), 2006, 57 (7): 1599-1604

[15] Qu Ping (曲平), Yan Rui (闫蕊), Yu Y uguo (俞裕国), Yao Pingjing

(姚平经). Operation performance of thermally coupled distillation of the butadiene separation unit [J]. J. Chem. Eng. Chinese Universities

(高校化学工程学报), 1996, 10 (2): 146-151

[16] Sun Lanyi (孙兰义), Zhang Mingyue (张明月), Zhou Hui (周辉)，

Yang Delian (杨德连), Li Jun (李军), Li Qingsong (李青松). Control of reactive distillation in dividing wall column for hydrolysis of methyl acetate [J]. Petrochemical Technology (石油化工), 2009, 11:

化工学报第66卷·1362·

1194-1200

[17] Sun L, Qi C, Li J, Li Q. Control of catalytic divided wall column [J].

Adv. Mater. Res., 2011, 225/226: 496-499

[18] Ma Zhanhua (马占华), Zhai Cheng (翟诚), Li Jun (李军)，Zhang Hui

(张慧)，Sun Lanyi (孙兰义). Application of air separation dividing wall column [J]. Chemical Engineering (化学工程), 2011, 40 (2): 1-6 [19] Sun L Y, Chang X W, Qi C X, Lia Q S. Implementation of ethanol

dehydration using dividing-wall heterogeneous azeotropic distillation column [J]. Separ. Sci. Technol., 2011, 46 (8): 1365-1375

[20] Sun L Y, Bi X X. Shortcut method for the design of reactive dividing

wall column [J]. Ind. Eng. Chem. R es., 2014, 53 (6): 2340-2347 [21] Ye Qing (叶青), Qian Chunjian (钱春健), Qiu Zhaorong (裘兆蓉).

Simulation of separation of aromatic hydrocarbon mixture in distillation column with diving wall [J]. Petrochemical Technology

(石油化工), 2007, 36 (11): 1134-1138

[22] Qian Chunjian (钱春健), Ye Qing (叶青)，Qiu Zhaorong (裘兆蓉),

Wen Yilu (闻以露). Separation of aromatic mixture by divided wall distillation column [J]. Petroleum Processing and Petrochemicals (石

油炼制与化工), 2007, 38 (5): 40-43

[23] Ye Qing (叶青), Li Langtao (李浪涛), Qiu Zhaorong (裘兆蓉).

Separation of pyrolysis gasoline by divided wall distillation column [J]. Acta Petrolei Sinica: Petroleum Processing Section (石油学报：石

油加工), 2010, 26 (4): 617-621

[24] Gao J, Huang K, Luan S, Jiao Y, Li S. Avoiding the black-hole

problem by the arrangements of multiple intermediate products to dividing-wall distillation columns [J]. Ind. Eng. Chem. Res., 2013, 52

(11): 4178-4201

[25] Luan S, Huang K, Wu N. Operation of dividing-wall columns (Ⅰ): A

simplified temperature difference control scheme [J]. Ind. Eng. Chem.

Res., 2013, 52 (7): 2642-2660

[26] Wu N, Huang K, Luan S. Operation of dividing-wall distillation

column (Ⅱ): A double temperature difference control scheme [J]. Ind.

Eng. Chem. Res., 2013, 52 (15): 5365-5383

[27] Yang Jian (杨剑), Shen Benqiang (沈本强), Lin Xiyu (蔺锡钰), Wu

Hao (吴昊), Ling Hao (凌昊). Steady state and dynamic control of divided-wall column for separating aromatics [J]. CIESC Journal (化

工学报), 2014, 65 (10): 3993-4003

[28] Kaibel G. Distillation columns with vertical partitions [J]. Chem. Eng.

Technol., 1987, 10: 92-98

[29] Dwivedi D, Strandberg J P, Halvorsen I J, Skogestad S. Steady state

and dynamic operation of four-product dividing-wall (Kaibel) columns: experimental verification [J]. Ind. Eng. Chem. Res., 2012,

51 (48): 15696-15709

[30] Dwivedi D, Strandberg J P, Halvorsen I J, Preisig H A, Skogestad S.

Active vapor split control for dividing-wall columns [J]. Ind. Eng.

Chem. Res., 2012, 51 (46): 15176-15183

[31] Ghadrdan M, Halvorsen I J, Skogestad S. Optimal operation of Kaibel

distillation columns [J]. Chem. Eng. Res. Des., 2011, 89 (8): 1382-1391

[32] Maryam Ghadrdan, Ivar J Halvorsen, Sigurd Skogestad. Manipulation

of vapor split in Kaibel distillation arrangements [J]. Chem. Eng.

Process. Process. Intensif., 2013, 72: 10-23

[33] Dejanovi? I, Matija?evi? L J, Halvorsen I J, Skogestad S, Jansen H,

Kaibel B, Oluji? ?. Designing four-product dividing wall columns for

separation of a multicomponent aromatics mixture [J]. Chem. Eng.

Res. Des., 2011, 89 (8): 1155-1167

[34] Kiss Anton A, Ignat Radu M, Flores Landaeta Servando J, de Haan

André B. Intensified process for aromatics separation powered by Kaibel and dividing-wall columns [J]. Chem. Eng. Process.: Process Intensif., 2013, 67: 39-48

[35] Luyben W L. Distillation Design and Control Using Aspen

Simulation [M]. New York: John Wiley & Sons Inc.,2006: 85-97