高级检索

  • ISSN 1006-3080
  • CN 31-1691/TQ
引用本文:
Citation:

一种改进的分隔壁精馏塔简捷计算方法

    作者简介: 叶贞成(1977—),男,博士,副研究员,主要研究方向为化工过程建模与优化。E-mail: yzc@ecust.edu.cn;
    通讯作者: 叶贞成, yzc@ecust.edu.cn
  • 中图分类号: TP273.1

An Improved Short Cut Design Method for Dividing Wall Column

    Corresponding author: Zhencheng YE, yzc@ecust.edu.cn ;
  • CLC number: TP273.1

  • 摘要: 针对现有分隔壁精馏塔四塔等效模型简捷计算方法的不足,提出了一种改进的计算方法,根据Underwood方程确定分隔壁精馏塔的最小气相流率。该方法可以结合设计要求快速确定精馏塔的塔板数、进出料位置、回流比、气液分配比等操作参数。为了验证设计方法的有效性,分别选取三种不同分离因子的进料物系进行设计并将所获得的设计结果在HYSYS软件中进行严格模拟,结果表明:与其他设计方法相比,该方法可以更准确地计算出预分馏塔与主塔之间连接流股的组成,且所需要的年总成本较小。
  • 图 1  四塔等效模型示意图

    Figure 1.  Schematic diagram of the four tower equivalent model

    图 2  简捷计算流程图

    Figure 2.  Shortcut procedure applied for the design of DWC

    图 3  精馏段逐板衡算示意图

    Figure 3.  By-board accounting of Rectifying section

    图 4  提馏段逐板衡算示意图

    Figure 4.  By-board accounting Stripping section

    图 5  HYSYS模型图

    Figure 5.  HYSYS model figure

    表 1  不同分离指数的进料

    Table 1.  Feeds with different separation indices

    Feed Mixtures(A/B/C) ESI
    F1 n-pentane/hexane/heptane 1.04
    F2 pentane/hexane/octane 0.46
    F3 ethane/propane/isobutane 1.79
    下载: 导出CSV

    表 2  设计规定

    Table 2.  Design Specification

    Feed D S xDA xSB xWC
    F1 33.4 31.67 0.9851 0.9900 0.9670
    F2 33.15 31.67 0.9945 0.9900 0.9583
    F3 33.52 31.67 0.9809 0.9900 0.9710
    下载: 导出CSV

    表 3  简捷计算结果

    Table 3.  Shortcut design results

    Feed N1_1 N1_2 N1 N2 N3_1 N3_2 N3 N4 NT R
    F1 10 11 20 6 8 8 17 6 29 4.32
    F2 6 7 12 3 4 6 11 4 18 3.04
    F3 8 6 13 6 8 4 13 2 21 4.32
    下载: 导出CSV

    表 4  连接流组分摩尔分数

    Table 4.  Connection flow component mole fraction

    Feed ${y_{{V_{1\_1}}i}}$ $ {x_{ {L_{1\_1} }i} } $ $ {y_{ {V_{1\_2} }i} } $ $ {x_{ {L_{1\_2} }i} } $
    F1 0.6201 0.5356 0.0009 0.0011
    0.3791 0.4634 0.4865 0.5757
    0.0008 0.0010 0.5125 0.4232
    F2 0.8021 0.7388 0.0009 0.0007
    0.1970 0.2600 0.6592 0.5307
    0.0009 0.0011 0.3399 0.4686
    F3 0.4936 0.3808 0.0018 0.0015
    0.5058 0.6184 0.1583 0.1372
    0.0006 0.0007 0.8390 0.8613
    下载: 导出CSV

    表 5  简捷计算与严格模拟结果对比

    Table 5.  Comparison between shortcut design and strict simulation results

    Feed $\Delta {y_{{V_{1 - 1}}i}}$ $\Delta {x_{{L_{1 - 1}}i}}$ $\Delta {y_{{V_{1 - 2}}i}}$ $\Delta {x_{{L_{1 - 2}}i}}$
    F1 0.0016 −0.0040 0.0000 0.0001
    0.0000 0.0040 −0.0012 −0.0020
    −0.0016 0.0000 0.0012 0.0019
    F2 0.0031 −0.0056 0.0000 0.0002
    0.0001 0.0056 −0.0035 −0.0027
    −0.0032 0.0001 0.0035 0.0025
    F3 0.0003 −0.0180 0.0000 0.0003
    0.0008 0.0179 −0.0046 −0.0054
    −0.0011 0.0000 0.0047 0.0051
    下载: 导出CSV

    表 6  本文方法与文献结果对比

    Table 6.  Comparison of the method in this paper and the results of literature

    Item N1-1 N1-2 N1 N2 N3-1 N3-2 N4 NT R L11 TAC
    This paper 10 11 20 6 8 8 6 29 4.32 35.47 154.28
    Literature[36] 10 16 25 5 12 13 8 38 3.66 26.57 177.08
    Literature[24] 11 11 21 7 11 10 11 39 3.17 24.44 162.63
    下载: 导出CSV

    表 7  本文方法与文献方法计算得到的连接流组分浓度差异对比

    Table 7.  Comparison of the difference in the component concentration of the connecting flow calculated by the method in this paper and the method in the literature

    Item $\Delta {y_{{V_{1 - 1}}i}}$ $\Delta {x_{{L_{1 - 1}}i}}$ $\Delta {y_{{V_{1 - 2}}i}}$ $\Delta {x_{{L_{1 - 2}}i}}$
    This paper 0.0016 0.0000 -0.0016 -0.0040 0.0040 0.0000 0.0000 -0.0012 0.0012 0.0001 -0.0020 0.0019
    Literature[36] 0.0086 0.0055 0.0031 0.0012 0.0004 0.0016 0.0028 0.0176 0.0148 0.0022 0.0086 0.0063
    Literature[24] 0.1098 0.1001 0.0096 0.2572 0.2566 0.0006 0.0026 0.1510 0.1536 0.0056 0.0818 0.0875
    下载: 导出CSV

    表 8  严格模拟结果产物浓度

    Table 8.  Strict simulation result of product concentration

    Feed xDA xSB xWC
    F1 0.9930 0.9902 0.9714
    F2 0.9993 0.9932 0.9690
    F3 0.9868 0.9899 0.9935
    下载: 导出CSV
  • [1] SHEN W F, BENYOUNES H, GERBAUD V. Extension of thermodynamic insights on batch extractive distillation to continuous operation: 1. Azeotropic mixtures with a heavy entrainer[J]. Industrial & Engineering Chemistry Research, 2013, 52(12): 4606-4622.
    [2] SHEN W F, BENYOUNES H, GERBAUD V. Extension of thermodynamic insights on batch extractive distillation to continuous operation: 2. Azeotropic mixtures with a heavy entrainer[J]. Industrial & Engineering Chemistry Research, 2013, 52(12): 4623-4637.
    [3] 李清元, 朱志亮. 四组分隔板塔热耦合精馏节能技术[J]. 化学工程, 2011(12): 6-10, 16. doi: 10.3969/j.issn.1005-9954.2011.12.002
    [4] 刘洪忠. 精馏耦合工艺研究进展及展望[J]. 化工技术与开发, 2018, 47(1): 5-41.
    [5] ÖMER Y, KISS A A, KENIG E Y. Dividing wall columns in chemical process industry: A review on current activities[J]. Separation & Purification Technology, 2011, 80(3): 403-417.
    [6] 孙兰义, 李军, 李青松. 隔壁塔技术进展[J]. 现代化工, 2008, 28(9): 38-41, 43. doi: 10.3321/j.issn:0253-4320.2008.09.008
    [7] 岳金彩, 闫飞, 邹亮, 等. 精馏过程节能技术[J]. 节能技术, 2008(1): 64-67. doi: 10.3969/j.issn.1002-6339.2008.01.017
    [8] 杨祖杰. 隔壁精馏塔节能应用研究[J]. 化工设计, 2019, 29(1): 25-28. doi: 10.3969/j.issn.1007-6247.2019.01.008
    [9] 高凌云, 齐鸣斋. 分壁式精馏塔模拟及节能研究[J]. 现代化工, 2015, 35(7): 135-138.
    [10] LING H, QIU J, HUA T, <italic>et al</italic>. Remixing analysis of four-product dividing-wall columns[J]. Chemical Engineering & Technology, 2018, 41(7): 1359-1367.
    [11] 朱登磊, 尚书勇, 谭超, 等. 基于分壁精馏塔的乙烯装置顺序分离新工艺及其模拟研究[J]. 石油学报(石油加工), 2014, 30(4): 682-686.
    [12] 李建清, 王二强, 黄金成, 等. 隔板塔应用于氯乙烯精馏过程的模拟分析[J]. 现代化工, 2012, 32(9): 70-72, 74. doi: 10.3969/j.issn.0253-4320.2012.09.020
    [13] 贾玉霞, 李玉安, 周文勇, 等. 分壁式精馏塔精制丁二烯流程模拟[J]. 化工进展, 2015, 34(10): 3563-3568, 3594.
    [14] 郭湘波, 王瑾. 分隔壁塔分离苯和乙烯烷基化产物的模拟[J]. 石油学报(石油加工), 2016, 32(3): 597-604.
    [15] DIGGELEN R C V, KISS A A, HEEMINK A W. Comparison of control strategies for dividing-wall columns[J]. Industrial & Engineering Chemistry Research, 2010, 49(1): 288-307.
    [16] HALVORSEN I J, SKOGESTAD S. Minimum energy eonsumption in multicomponent distillation: 1. Vmin diagram for a two-product column[J]. Industrial & Engineering Chemistry Research, 2003, 42(3): 596-604.
    [17] HALVORSEN I J, SKOGESTAD S. Minimum energy consumption in multicomponent distillation: 2. Three-product petlyuk arrangements[J]. Industrial & Engineering Chemistry Research, 2003, 42(3): 605-615.
    [18] HALVORSEN I J, SKOGESTAD S. Minimum energy consumption in multicomponent distillation: 3. More than three products and generalized petlyuk arrangements[J]. Industrial & Engineering Chemistry Research, 2003, 42(3): 616-629.
    [19] 沈鸳语, 黄泽恩, 石剑等. 隔壁精馏塔结构参数优化研究[J]. 化学工程, 2019, 47(1): 26-31. doi: 10.3969/j.issn.1005-9954.2019.01.006
    [20] 黄国强, 靳权. 隔壁精馏塔的设计、模拟与优化[J]. 天津大学学报(自然科学与工程技术版), 2014, 47(12): 1057-1064.
    [21] 何西涛, 干爱华, 陈宁. 分壁精馏塔分离苯/甲苯/二甲苯的模拟工艺研究[J]. 化学工业与工程, 2012, 29(02): 51-57. doi: 10.3969/j.issn.1004-9533.2012.02.011
    [22] 龚超, 余爱平, 罗祎青, 等. 完全能量耦合精馏塔的设计、模拟与优化[J]. 化工学报, 2012, 63(1): 177-184. doi: 10.3969/j.issn.0438-1157.2012.01.025
    [23] SALVADOR H, JUAN G S, VICENTE R. Thermodynamically equivalent distillation schemes to the petlyuk column for ternary mixtures[J]. Energy, 2006, 31(12): 2176-2183. doi: 10.1016/j.energy.2005.10.007
    [24] TRIANTAFYLLOU C, SMITH R. The design and optimization of dividing wall distillation columns[J]. Energy Efficiency in Process Technology, 1993: 351-360.
    [25] KIM Y H. Structural design and operation of a fully thermally coupled distillation column[J]. Chemical Engineering Journal, 2002, 85(2): 289-301.
    [26] KIM Y H. A new fully thermally coupled distillation column With postfractionator[J]. Chemical Engineering and Processing, 2006, 45(4): 254-263. doi: 10.1016/j.cep.2005.03.013
    [27] NELLY R C, ARTURO J G, ANGEL C A, <italic>et al</italic>. Optimum design of petlyuk and divided-wall distillation systems using a shortcut model[J]. Chemical Engineering Research & Design, 2010, 88(10): 1405-1418.
    [28] CHU K T, CADORET L, YU C C, <italic>et al</italic>. A new shortcut design method and economic analysis of divided wall columns[J]. Industrial & Engineering Chemistry Research, 2011, 50(15): 9221-9235.
    [29] BENYOUNES H, BENYAHIA K, SHEN W F, <italic>et al</italic>. Novel procedure for assessment of feasible design parameters of dividing-wall columns: Application to non-azeotropic mixtures[J]. Industrial & Engineering Chemistry Research, 2015, 54(19): 5307-5318.
    [30] SOTUDEH N, SHAHRAKI B H. [J]. Chemical Engineering & Technology, 2007, 30(9): 1284-1291.
    [31] UWITONZE H, HAN S, KIM S, <italic>et al</italic>. Structural design of fully thermally coupled distillation column using approximate group methods[J]. Chemical Engineering and Processing: Process Intensification, 2014, 85: 55-167.
    [32] UWITONZE H, HAN S, HWANG K S. New design method for fully thermally coupled distillation column using group and approximate methods[J]. Industrial & Engineering Chemistry Research, 2014, 53(30): 1979-11988.
    [33] UWITONZE H, HWANG K S, LEEA I. A new design method and operation of fully thermally coupled distillation column[J]. Chemical Engineering and Processing, 2016: 47-58.
    [34] SEIHOUB F Z, BENYOUNES H, SHEN W F, <italic>et al</italic>. An improved shortcut design method of divided wall columns exemplified by a liquefied petroleum gas process[J]. Industrial & Engineering Chemistry Research, 2017, 56(34): 9710-9720.
    [35] LORENA E R, NELLY R, ANGEL C, <italic>et al</italic>. Shortcut Design of fully thermally coupled distillation systems with postfractionator[J]. Industrial & Engineering Chemistry Research, 2011, 50(10): 6287-6296.
    [36] DEJANOVIC I, MATIJASEVIC L, OLUJIC Z. Dividing wall column-a breakthrough towards sustainable distilling[J]. Chemical Engineering and Processing, 2010, 49(6): 559-580. doi: 10.1016/j.cep.2010.04.001
    [37] DEJANOVIC I, MATIJASEVIC L, OLUJIC Z. An effective method for establishing the stage and reflux requirement of three-product dividing wall columns[J]. Chemical and Biochemical Engineering Quarterly, 2011, 25(2): 147-157.
  • [1] 王俊黄秀辉崔国民 . 碘硫循环制氢中HI浓缩精馏的模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180927003
    [2] 贾翔飞钱嘉澍吴幼青陈剑佩 . 鼓泡塔反应器中两相流动CFD-PBM耦合数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191121001
    [3] 杨祺刘士荣 . 多自主车辆队列跟随控制器设计. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190118001
    [4] 许叶龙刘迎圆惠虎於晔鸿 . 固液混合过程的数值模拟及实验研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180709001
    [5] 全丽娜刁妍妍李巧李文杰朱丽丽李洪林 . 靶向PD-1的多肽阻断剂的设计及评价. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20181126003
    [6] 汪伟亮邓卫平 . ChemCloser:一个基于路径与片段匹配的药物设计软件. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200310004
    [7] 高天阳虞慧群范贵生 . 基于模拟退火遗传算法的云资源调度方法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180416001
    [8] 罗安王汉奎王建文 . 基于小冲杆试验数据的力学性能的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180609002
    [9] 胡志梁于新海赖焕新 . 拟塑性流体在六直叶搅拌器中的模拟与比较. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180517006
    [10] 张丹杨敏博冯霄 . 循环流化床甲醇制芳烃分离工艺的模拟与改进. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180911004
    [11] 潘傲谢明辉夏建业庄英萍 . 基于均龄理论模拟搅拌式反应器的混合时间. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180716001
    [12] 李岁王元华 . 油田水套加热炉高温空气燃烧瞬态模拟及最小换向时间. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180723008
    [13] 胡贵华隆建杜文莉 . 考虑汽化影响的乙烯裂解炉对流段数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913002
    [14] 蔡冬莹王跃林段先建胡丹李玉冰胡彦杰 . 预混合高速射流燃烧反应器内温度场的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190221004
    [15] 肖桐王千红沈盈莺惠卿陈晗程振民 . 颗粒床内固液过滤的三维CFD-DEM模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190223001
    [16] 肖卫强周国俊蒋健胡安福詹望成郭杨龙 . 加热卷烟制品传热与烟气流动过程数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191029004
    [17] 陈阳刘锦泽于兆臣王智峰欧阳福生 . 重油催化裂化12集总动力学模型的工业模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200316002
    [18] 汤渊潘伟童梁钦锋许建良代正华于广锁王辅臣 . 管道内高温合成气喷雾激冷过程数值模拟研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191207001
    [19] 罗雪莉郎美东 . 聚合物前药载药性能的计算机模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190514001
    [20] 于苗钱锋胡贵华隆建李天越 . 工业级MIP提升管反应器气固两相流动特性的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913003
  • 加载中
图(5)表(8)
计量
  • 文章访问数:  2039
  • HTML全文浏览量:  585
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-10
  • 网络出版日期:  2019-12-28

一种改进的分隔壁精馏塔简捷计算方法

    通讯作者: 叶贞成, yzc@ecust.edu.cn
    作者简介:叶贞成(1977—),男,博士,副研究员,主要研究方向为化工过程建模与优化。E-mail: yzc@ecust.edu.cn
  • 华东理工大学化工过程先进控制和优化技术教育部重点实验室,上海 200237

摘要: 针对现有分隔壁精馏塔四塔等效模型简捷计算方法的不足,提出了一种改进的计算方法,根据Underwood方程确定分隔壁精馏塔的最小气相流率。该方法可以结合设计要求快速确定精馏塔的塔板数、进出料位置、回流比、气液分配比等操作参数。为了验证设计方法的有效性,分别选取三种不同分离因子的进料物系进行设计并将所获得的设计结果在HYSYS软件中进行严格模拟,结果表明:与其他设计方法相比,该方法可以更准确地计算出预分馏塔与主塔之间连接流股的组成,且所需要的年总成本较小。

English Abstract

  • 精馏是化学工业中应用最广泛的分离技术之一,其基本原理是在塔釜加热液相物流形成气相回流,在塔顶冷凝气相物流形成液相回流,两者在塔内逆向流动并进行传热和传质,实现轻重组分的分离与提纯。据统计,在一个典型石油化工厂中,精馏的能耗约占全厂总能耗的40%左右,精馏过程中,进入再沸器的热量有95%需要在塔顶冷凝器中移走。因此,精馏塔分离系统的高效设计和运行在化工生产中至关重要。为了减少精馏过程中的能耗,开发可持续的、经济的精馏系统显得尤为必要[1-2]。目前精馏过程的主要节能方式可以分为两类[3]:精馏塔设备的改进,如采用高效规整填料代替塔板或者低效填料等;能量的耦合,如分隔壁精馏、结晶精馏、热泵精馏、反应精馏、多效精馏、共沸精馏等[4-5]。分隔壁精馏塔作为一种新型的精馏方式,已被证实既可以大幅提高精馏塔的热力学效率,又可以减少设备投资[6-9],近年来受到广泛的关注。文献[10-14]分别将分隔壁精馏塔运用到不同的生产工艺中,结果表明较传统工艺可节能30%左右。

    与传统的精馏塔工艺相比,分隔壁精馏塔具有更多的自由度,这给分隔壁精馏塔的设计、优化、控制带来很大的困难[15]。Halvoren等人[16-18]采用Underwood方程在无穷塔板数的极限情况下分别对两组分和三组分分离过程所需要的最小上升气相流率进行分析并将其扩展到多组分分离过程,然后通过Vmin图直观地展示分离不同组分需要的最小上升气相流率,但对精馏塔结构设计却未提及。在对分隔壁精馏塔结构进行设计时国内外的学者均采用了不同形式的简化,文献[19-21]采用三个普通的精馏塔等效替代分隔壁精馏塔,基于Aspen plus的DSTWU模块进行简捷计算得到基本结构参数。文献[22-25]采用一个普通精馏塔一个回流吸收塔一个再沸吸收联合等效替代分隔壁精馏塔,基于Fenske-Underwood-Gilliland-Kirkbride(FUGK)方程进行了简捷设计。用普通的精馏塔等效预分馏塔这种过度简化,会对计算结果的准确性产生一定程度的影响。文献[26-29]采用Petlyk结构基于FUGK方程对分隔壁精馏塔进行简捷计算与设计,取得了不错的效果。Sotudeh[30]认为用Fenske方程计算分隔壁精馏塔的塔板数是不正确的,因为在全回流的情况下分隔壁精馏塔中预分馏塔顶部液相进料和气相出料组分浓度以及底部液相出料和气相进料分浓度是不相等的,不同于传统精馏塔全回流情况下顶部、底部进出料组分浓度相等的情况。因此文献[26-29]中的简捷计算方法都具有Sotudeh指出的不足之处。

    针对采用Fenske方程计算预分馏塔塔板数存在的问题,Uwitonze等[31-33]采用逐板衡算的方法对原计算方法进行改进,提高了计算的准确性。由于预分馏塔与主塔之间是相互连接的,通常的多组分设计程序在未给出相互连接流信息的情况下不适用于主塔的设计。Seihoub等[34]对连接流股的组成进行了研究,结合文献[31]中逐板衡算的优点对各段的塔板数进行了计算,然后通过调节回流比使得预分馏塔和侧线采出塔的塔板相等。由于预分馏塔和侧线采出塔分离任务不同,预分馏塔和侧线采出塔的塔板数可以是不相等的,因此这种算法不仅增加了计算复杂度,同时也大大限制了分隔壁精馏塔的灵活性。

    对于分隔壁精馏塔来说因进料与塔板之间、连接流与主塔之间组分不一致而造成的返混,是导致分隔壁精馏塔热力学效率低下的主要原因[25, 35]。本文提出了一种改进的简捷计算方式,采用四塔结构,结合了上述文献中简捷计算方法的优点,采用Underwood方程对气液流量进行分析,并结合逐板衡算的方法对塔板数进行计算,同时取消预分馏塔与侧线采出塔塔板数目相等的约束条件,在HYSYS中进行严格模拟,并与文献[24, 36]中计算结果进行对比分析。结果表明该方法可以更有效地缩小预分馏塔与主塔之间的浓度差,更准确地计算出所需要的塔板数、进料位置、隔板位置等参数。

    • 文献[36]指出四塔等效模型相对三塔、二塔等效模型能更好地体现分隔壁精馏塔的操作特性,因而采用四塔模型进行设计。图1为用四塔模型等效替代分隔壁精馏塔的示意图,将分隔壁精馏塔划分为4个塔段:顶部精馏塔、预分馏塔、侧线采出塔、底部提馏塔。预分馏塔上升气相与侧线采出塔上升气相混合进入顶部精馏塔;预分馏塔下降液相与侧线采出塔下降液相混合进入底部提馏塔;顶部精馏塔下降液相被分成两股,分别进入预分馏塔和侧线采出塔;底部提馏塔上升气相被分为两股,分别进入预分馏塔和侧线采出塔。本文拟对三组分分离的四塔模型设计进行研究。分隔壁精馏塔的设计步骤如图2所示。

      图  1  四塔等效模型示意图

      Figure 1.  Schematic diagram of the four tower equivalent model

      图  2  简捷计算流程图

      Figure 2.  Shortcut procedure applied for the design of DWC

      (1)定义预分馏塔顶部轻组分A的收率 $ {r_{{{\rm{D}}_{{\rm{1,A}}}}}} $ 和重组分C的收率 $ {r_{{{\rm{D}}_{{\rm{1,C}}}}}} $

      (2)根据插值法求出预分馏塔顶部中间组分B的收率 $ {r_{{{\rm{D}}_{{\rm{1,B}}}}}} $

      (3)基于Underwood方程估算达到规定分离任务时所需的最小气相负荷Vmin和最小回流比Rmin;

      (4)确定回流比R和精馏塔内气液两相流量LiVi

      (5)根据预分馏塔顶部与顶部精馏塔的物料衡算得出顶部连接流处的组分浓度,根据预分馏塔底部与底部提馏塔的物料衡算得出底部连接流处的组分浓度。

      (6)通过逐板衡算的方式计算出各段的塔板数Ni和进出料位置。

      其中:Section1为预分馏塔、S1_1为预分馏塔精馏段、S1_2为预分馏塔提馏段、Section2为顶部精馏塔、Section3为侧线采出塔、S3_1为侧线采出塔提馏段、S3_2为侧线采出塔精馏段、Section4为底部提馏塔。

      其中:D1为预分馏塔净采出量,R1min为预分馏塔最小回流比,Rmin为主塔最小回流比,R1为预分馏塔实际回流比,R为主塔实际回流比, ${x_{{L_{1\_1}}i}} $ ${y_{{V_{1\_1}}i}} $ ${x_{{L_{1\_2}}i}} $ ${y_{{V_{1\_2}}i}} $ 为连接流股组分摩尔分数,Ki为气液平衡常数,Ni为塔板数,ai是与气液相流量和Ki相关的数,详见下文。

    • 组分i向上通过塔板n的净流量wi等于向上流动气相与向下流动液相中组分i物质的量之差:

      在预分馏塔顶部:

      其中: $ {D_{\rm{1}}} = {V_{1\_1}} - {L_{1\_1}} $

      在预分馏塔底部:

      其中: $ {W_{\rm{1}}} = {L_{1\_2}} - {V_{1\_2}} $ 。净组分流动的正方向是向上定义的,此时有 $w_i^{\rm{B}} \le 0$

      根据Underwood方程此时预分馏塔的精馏段上升气相值为:

      解方程(4)可得到3个根,有 ${\alpha _A} > {\theta _1} > {\alpha _B} > {\theta _2} > {\alpha _C} > {\theta _3}$

      预分馏塔的提馏段上升气相值为:

      解方程(5)可得到3个根,有 ${\delta _1} > {\alpha _{\rm{A}}} > {\delta _2} > {\alpha _{\rm{B}}} > {\delta _3} > {\alpha _{\rm{C}}}$

      根据上升气相物料平衡方程有:

      其中:q1为进料热状态,当q1<0时为过热蒸汽进料,当q1=0时为饱和蒸汽进料,当0<q1<1时为气液混合进料,当q1=1时为饱和液体进料,当1<q1时为过冷液体进料,这里取q1=1,即 $ {\theta _1} = {\delta _2} = {\phi _{11}}, {\theta _2} = {\delta _3} = {\phi _{12}}$

      Underwood[16]得出对于任何给定的产品分布,方程(4)和方程(5)至少有一对公共根,它们同时满足方程(6)。对方程(6)根据物料平衡,进一步简化可以得到方程(7)

      求解方程(7)可以得到2个根,有 ${\alpha _{\rm{A}}} > {\phi _{11}} > {\alpha _{\rm{B}}} > {\phi _{12}} > {\alpha _{\rm{C}}}$

      因此可以得到精馏段最小上升气相流率为:

      预分馏塔的任务是尽可能的将组分AC进行分离,因此,对预分馏塔顶部组分的收率进行定义:

      预分馏塔的最小理论板数可以由Fenske方程进行估算:

      中间组分B在有限回流下的收率估计值可以由非清晰分割[39]的物料衡算法计算得出:

      对公式(10)进行变换计算可以得到:

      在已知预分馏塔顶部组分B收率的基础上,可以得出预分馏塔塔顶净采出量:

      此时,预分馏塔的最小回流比可以很容易得到:

      预分馏塔实际回流比取为最小回流比的b倍:

    • 对于顶部精馏塔,将预分馏塔顶部进入的液相和离开的气相结合形成一个净流 ${D_{\rm{P}}}$ ,则进入顶部精馏塔的物料热状态 ${q_2}$ 可以定义为:

      求解公式(16)可以得到两个 ${\phi _{2z}}$ 的值,且αA > ϕ21 > αB > ϕ22 > αC。根据所求得的 ${\phi _{2z}}$ 值结合Underwood方程可以得到顶部精馏塔的最小上升气相流率:

    • 对于底部提馏塔,将预分馏塔底部离开的液相和进入的气相结合形成一个净流 ${W_{\rm{P}}}$ ,则进入底部提馏塔的物料热状态 ${q_4}$ 可以定义为:

      求解公式(19)可以得到两个 ${\phi _{4z}}$ 的值,且αA > ϕ41 > αB > ϕ42 > αC。根据所求得的 ${\phi _{4z}}$ 值结合Underwood方程可以得到底部提馏塔的最小上升气相流率:

    • 底部提馏塔的上升气相分别进入预分馏塔和侧线采出塔段,因此有:

      由于侧线采出塔提馏段和侧线采出塔精馏段形成的侧线采出塔中没有气相的采出,因此进入这两个塔段的蒸汽流率相等, ${V_{3\_1\min }} = {V_{3\_2\min }}$

      通过上述对四塔模型每个塔段最小上升气相流率的讨论可以得出,主塔部分的最小上升气相流率为:

      主塔的最小回流比为:

      实际回流比为最小回流比的a倍:

      至此可以确定每个塔段的气液相流量:

    • 为了得到每个塔段的准确塔板数,需准确计算各部分连接流的组分浓度。

      对顶部精馏塔,预分馏塔顶部进行物料衡算:

      为了减小预分馏塔与主塔之间连接流股的浓度差( $\Delta x$ ),联立方程(31)(32),求解可得顶部各个组分的浓度:

      对底部提馏塔,预分馏塔底部进行物料衡算:

      同理联立方程(35)(36),可得到底部各个组分的浓度:

    • 在对4股连接流组分浓度求解的基础上,进一步求解四塔模型各部分塔板数。此外,为了简化计算,对每个部分中的气液平衡常数 K 取值为常数。由每个塔段中关键组分的物料衡算结合相平衡方程,采用逐板衡算的方式进行计算得出各塔段的塔板数。

      (1)精馏段逐板计算示意图如图3所示。

      图  3  精馏段逐板衡算示意图

      Figure 3.  By-board accounting of Rectifying section

      精馏段塔板数计算方法:

      在塔顶处n=N ,可得:

      其中, $a = \displaystyle\frac{L}{{KV}}$

      (2)提馏段逐板计算示意图如图4所示。

      图  4  提馏段逐板衡算示意图

      Figure 4.  By-board accounting Stripping section

      提馏段塔板数计算方程:

      在塔底处n=N ,可得:

      其中, $a = \displaystyle\frac{{KV}}{L}$ ,经过上文关于塔板数计算公式的推算,可以用公式(42)对不同塔段的精馏段进行计算,用公式(46)对不同塔段的提馏段进行计算。

    • 预分馏塔精馏段塔板数:

      预分馏塔提馏段塔板数:

      则预分馏塔塔板数N1为:

    • 顶部精馏塔塔板数:

    • 侧线采出塔提馏段塔板数:

      侧线采出塔精馏段塔板数:

      则侧线采出塔塔板数N3为:

    • 底部提馏塔塔板数:

    • 为了验证方法的合理性,采用了文献[34]中的数据和设计要求,并将本文的设计结果与其进行对比。选取具有不同分离指数的三组分(如表1所示)作为被分离物系,分离指数(Ease of Separability Index,ESI)定义如下:

      由ESI的定义可知:如果ESI<1,组分A与B的分离比组分B与C分离更加困难;如果ESI>1,组分A与B的分离比组分B与C分离更加容易;如果ESI=1,组分A与B的分离与组分B与C的分离难易程度相当。

      Feed Mixtures(A/B/C) ESI
      F1 n-pentane/hexane/heptane 1.04
      F2 pentane/hexane/octane 0.46
      F3 ethane/propane/isobutane 1.79

      表 1  不同分离指数的进料

      Table 1.  Feeds with different separation indices

      混合物在202.65kPa下以液相进料,进料流量为100 kmol/h,组分(A/B/C)摩尔分数分别为0.33、0.33、0.34,要求中间组分B的收率大于等于0.95且侧线采出组分摩尔分数大于0.99。具体设计规定如表2所示。

      Feed D S xDA xSB xWC
      F1 33.4 31.67 0.9851 0.9900 0.9670
      F2 33.15 31.67 0.9945 0.9900 0.9583
      F3 33.52 31.67 0.9809 0.9900 0.9710

      表 2  设计规定

      Table 2.  Design Specification

      采用本文改进的简捷计算方法,在 $a = b = 1.2$ 时对表1中不同进料情况下的分隔壁精馏塔进行计算,结果如表3所示。

      Feed N1_1 N1_2 N1 N2 N3_1 N3_2 N3 N4 NT R
      F1 10 11 20 6 8 8 17 6 29 4.32
      F2 6 7 12 3 4 6 11 4 18 3.04
      F3 8 6 13 6 8 4 13 2 21 4.32

      表 3  简捷计算结果

      Table 3.  Shortcut design results

      与ESI接近1的物系相比,ESI大于1或ESI小于1的物系需要更少的总塔板数。此外,ESI小于1的物系需要更少的精馏段塔板数,由于组分A与B分离比B与C分离容易,因此在预分馏塔中大部分中间组分从塔底采出;ESI大于1的物系需要更少的提馏段塔板数,由于A与B分离比B与C分离困难,因此在预分馏塔中大部分中间组分从塔顶采出。

      将本文通过简捷计算得到的参数,如塔板数,连接流组分浓度(如表4所示)等,输入Aspen HYSYS软件中进行严格模拟,并将HYSYS模拟结果与简捷计算结果进行对比,如表5所示,可以看出两者结果相差很小,表明该方法适用于不同结构得系统,HYSYS模型如图5所示。

      Feed ${y_{{V_{1\_1}}i}}$ $ {x_{ {L_{1\_1} }i} } $ $ {y_{ {V_{1\_2} }i} } $ $ {x_{ {L_{1\_2} }i} } $
      F1 0.6201 0.5356 0.0009 0.0011
      0.3791 0.4634 0.4865 0.5757
      0.0008 0.0010 0.5125 0.4232
      F2 0.8021 0.7388 0.0009 0.0007
      0.1970 0.2600 0.6592 0.5307
      0.0009 0.0011 0.3399 0.4686
      F3 0.4936 0.3808 0.0018 0.0015
      0.5058 0.6184 0.1583 0.1372
      0.0006 0.0007 0.8390 0.8613

      表 4  连接流组分摩尔分数

      Table 4.  Connection flow component mole fraction

      Feed $\Delta {y_{{V_{1 - 1}}i}}$ $\Delta {x_{{L_{1 - 1}}i}}$ $\Delta {y_{{V_{1 - 2}}i}}$ $\Delta {x_{{L_{1 - 2}}i}}$
      F1 0.0016 −0.0040 0.0000 0.0001
      0.0000 0.0040 −0.0012 −0.0020
      −0.0016 0.0000 0.0012 0.0019
      F2 0.0031 −0.0056 0.0000 0.0002
      0.0001 0.0056 −0.0035 −0.0027
      −0.0032 0.0001 0.0035 0.0025
      F3 0.0003 −0.0180 0.0000 0.0003
      0.0008 0.0179 −0.0046 −0.0054
      −0.0011 0.0000 0.0047 0.0051

      表 5  简捷计算与严格模拟结果对比

      Table 5.  Comparison between shortcut design and strict simulation results

      图  5  HYSYS模型图

      Figure 5.  HYSYS model figure

      其中:T1,T2为分配器,M1,M2为混合器,R1,R2为循环器,Q11和Q12分别表示塔顶冷负荷和塔釜热负荷,其它设备号名称与图1对应。

      在同样的条件下分别采用本文的方法与文献[36, 24]提出的方法,对戊烷/己烷/庚烷进料体系进行设计,具体计算结果如表6所示。

      Item N1-1 N1-2 N1 N2 N3-1 N3-2 N4 NT R L11 TAC
      This paper 10 11 20 6 8 8 6 29 4.32 35.47 154.28
      Literature[36] 10 16 25 5 12 13 8 38 3.66 26.57 177.08
      Literature[24] 11 11 21 7 11 10 11 39 3.17 24.44 162.63

      表 6  本文方法与文献结果对比

      Table 6.  Comparison of the method in this paper and the results of literature

      从表6中的计算结果可以看出在相同的设计要求下,本文的设计方法相比文献[24,36]设计方法需要更少的理论板数(29),虽然比其他计算方法需要更大的回流比(4.32),但是根据Dejanovic[37]提出的使用 $N\left( {R + 1} \right)$ 近似计算精馏塔年总成本,可以计算出本文的简捷计算方法年投资成本要小于文献[24,36]方法。

      分别采用本文方法和文献方法对戊烷/己烷/庚烷进料体系得到的预分馏塔与主塔之间的连接流组成差异进行比较,结果如表7所示。可以看出,本文方法的计算结果误差要小于文献[24,36]的计算结果误差。

      Item $\Delta {y_{{V_{1 - 1}}i}}$ $\Delta {x_{{L_{1 - 1}}i}}$ $\Delta {y_{{V_{1 - 2}}i}}$ $\Delta {x_{{L_{1 - 2}}i}}$
      This paper 0.0016 0.0000 -0.0016 -0.0040 0.0040 0.0000 0.0000 -0.0012 0.0012 0.0001 -0.0020 0.0019
      Literature[36] 0.0086 0.0055 0.0031 0.0012 0.0004 0.0016 0.0028 0.0176 0.0148 0.0022 0.0086 0.0063
      Literature[24] 0.1098 0.1001 0.0096 0.2572 0.2566 0.0006 0.0026 0.1510 0.1536 0.0056 0.0818 0.0875

      表 7  本文方法与文献方法计算得到的连接流组分浓度差异对比

      Table 7.  Comparison of the difference in the component concentration of the connecting flow calculated by the method in this paper and the method in the literature

      严格模拟后得到的产物的关键浓度如表8所示。可以看出产物浓度均满足设计要求。

      Feed xDA xSB xWC
      F1 0.9930 0.9902 0.9714
      F2 0.9993 0.9932 0.9690
      F3 0.9868 0.9899 0.9935

      表 8  严格模拟结果产物浓度

      Table 8.  Strict simulation result of product concentration

    • 针对分离三元组分混合物的分隔壁精馏塔设计提供了一种简捷计算方法,采用逐板衡算的方法计算塔板数,避免了FUGK方程的不足。消除了在半严格设计过程中对预分馏塔与主塔之间连接流组分计算时繁琐的迭代过程,大大减少了计算量。同时针对前人要求预分馏塔塔板数与侧线采出塔塔板数相同的不足进行改进,更准确地计算出预分馏塔与主塔之间的连接流的组成。结果表明,本文对分隔壁精馏塔的简捷计算方法的改进是有效的。

      本文所提出的简捷计算方法基于恒摩尔流假定,同时在计算塔板数时假定三元理想混合物具有恒定的相对挥发度,因此该方法计算速度快易于理解,本文方法为分隔壁精馏塔的应用设计提供了理论基础。

(5)  表(8) 参考文献 (37) 相关文章 (20)

目录

    /

    返回文章