《精馏塔的动态模拟》由会员分享，可在线阅读，更多相关《精馏塔的动态模拟（18页珍藏版）》请在金锄头文库上搜索。
1、精馏塔的动态模拟目 录一、数学建模二、分析与讨论三、优化四、程序清单（）分析和讨论（）分析和讨论（）最优化部分前 言化学工业中，精馏过程是能量消耗最大的单元操作之一，自从发生了世界性的能源问题以来，精馏过程的节能问题已广泛引起了人们的重视。近年来，已经开发了多种精馏节能的工艺流程，如多效精馏.热泵精馏、热偶精馏等。多效精馏作为一种精馏节能新工艺近几年来其理论研究不断深入，在工业生产中的应用日益广泛。工业上普遍存在非稳定状态或动态的精馏过程.实际的生产过程不可避免地受到各种人为或非人为因素的扰动，使一些操作参数和过程变量随时间发生变化，因此对于连续精馏过程的动态特性的研究和模拟具有重要的实际意义
2、.研究精馏塔的动态特性时通常使用数学模拟方法，这首先需要建立 1 个模型精馏塔，然后对模型塔中的各个塔板作非稳态物料衡算、热量衡算等，得到动态精馏过程的数学模型，在一定的初始条件下经过求解，可得到操作条件发生扰动时各个精馏过程参数随时间变化的历程，即动态响应.若精馏系统内无约束某些参数变化速度和变化范围的控制器，其动态响应称为开环响应，否则为闭环响应.反应精馏技术将反应与分离过程在一个塔内实现, 相对传统的先反应后分离过程具有转化率高、选择性好、操作易、投资省等一系列优点, 但反应精馏过程需同时遵循质量作用定律和精馏分离原理, 其过程影响因素复杂 . 自 20 世纪 70 年代以来, 有关反应
3、精馏的研究重点从工艺转向数学模拟. 现已开发有灵活可靠的过程模拟计算软件. 另一方面, 自 20 世纪 70 年代末催化精馏技术成功地应用于甲基叔丁基醚(M TBE) 生产以后, 该技术的应用受到了学者们的关注. 一、精馏塔数学模型的建立：根据对过程系统中状态变量分布特征的不同描述方法，一般可以把数学模型分为集中参数模型，分布参数模型和多级集中参数模型。本次大作业利用多级集中参数模型对精馏塔动态特性进行分析和模拟。对于控制的动态数学模型，我们希望用最简单的形式，最大限度地概括出过程的特性。所以为了简化数学模型，我们必须做出以下必要的假定：每块塔板上气相和液相分别为理想混合；1、组分 A 和组分
4、 B 的摩尔汽化热近似相等；汽相和液相在沿塔轴向运动过程中，显热变化对热量衡算和热损失的影响均可忽略不计；2、整个塔分离过程 A 对 B 的相对挥发度 保持恒定；3、每块踏板上都达到了汽液相平衡；4、冷凝器和再沸器均处于稳态操作；5、采用泡点进料和常压进料；6、每块踏板上持液量远大于持汽量，后者及其变化可以忽略；7、塔顶冷凝器为全冷凝器。在这些假设的基础上，以物料衡算及能量衡算为基础列出精馏塔动态数学模型：全凝器及馏出塔总物料衡算：DLVtMRD d全凝器及馏出液易挥发组分衡算：)( d1DDXYVtXM 第 n 块塔板总物料衡算：nnLtM1nd 第 n 块板易挥发组分衡算：)()( d11
5、n nnnn YVXLtXM 离开第 n 块板汽液相浓度关系：nnn11XY） （ 加料板物料衡算及易挥发组分汽液相浓度平衡关系：FLMFt 1 d)()()( d 1f11 FFFFF YVXXLtXM FFXY)1( 再沸器及塔底总物料衡算： BVLtMNB d 再沸器及塔底易挥发组分衡算：)()( dBBNB YXVXLtM 离开再沸器及塔底的汽液相浓度关系：BBXY)1( 塔板水力学关系：（12）51n028.( .AhMLLWN ） 混 混混混 混 （13）iX混（14）iWrr混查得苯和甲苯的物性数据如下表：T/苯0PKpa甲 苯0KpaK苯g/m3甲 苯Kg/m3甲KJ/Kg甲K
6、J/Kg苯x%苯y%80.2 101.7941 39.3471 814.778 809.98 393.956 379.778 100.0 100.080.66 103.243 39.9743 814.267 809.35 393.625 379.497 97.0 98.8 2.5463981.2 104.9646 40.7209 813.668 808.824 393.236 379.168 95.0 97.9 2.4536382.3 108.5409 42.2765 812.447 807.746 392.444 378.497 90.3 95.7 2.3907183.33 111.9754
7、 43.7761 811.303 806.74 391.702 377.869 85.3 93.6 2.5203784.4 115.6325 45.3789 810.116 805.688 390.932 377.216 80.3 91.4 2.6073585.44 119.2759 46.9817 808.962 804.669 390.183 376.582 75.0 88.5 2.5652286.8 124.1747 49.1458 807.452 803.336 389.204 375.752 70.0 85.3 2.4868887.63 127.2405 50.5054 806.53
8、1 802.523 388.606 375.246 65.0 82.5 2.5384689.4 133.9751 53.5055 799.566 800.788 387.332 374.166 59.2 78.9 2.5771192.1 144.778 58.3549 801.506 798.121 385.304 372.498 48.9 71.0 2.5584295.2 158.0008 64.3490 795.052 795.052 382.948 370.576 39.7 61.8 2.4572798.6 173.5565 71.4767 794.096 791.686 380.364
9、 368.468 30.0 50.0 2.33333102.2 191.2863 79.6929 789.816 788.1 377.584 366.192 20.0 37.0 2.34921106.1 212.029 89.4199 785.058 784.2 374.542 363.696 8.0 21.2 2.78817110.6 238.054 101.783 779.616 779.682 371.002 360.804 0.0 0.0 0.00000塔的示意图和设备尺寸如下：溢流管形式 堰长 Wl堰高 孔数 n 管底与受液盘距离 塔径 塔板 数 孔径弓形 0.66m 0.04m 2
10、738 个 0.032m 1.0m 10 5mm设进料温度为 90.8， 2.4方程，简化为：（15） nnXY4.12（16） FF. （17）BBXY4.12x50，y50： g/mol8592.0785混 MKJ/Kg.37.395.0混 WrKg/m32.0.284混（18）QEV61. 5.125.8039.028.475.08. .5.3n nML（19）.7.16412.43nE二模型微分方程的仿真计算：计算流程如下：分析和讨论总结：开 始给定 Q、F、D、B;赋 M、X 初值计算：V, L 计算：Y求解常微分方程初值问题XDX D,设计值XB XB,设计值 结束赋t从模型数据结
11、果分析来看，与实际情况即理论计算都不吻合。这主要是与所建立的模型的精确性以及程序算法有关。由于平衡级动态模型采用简化假设与实际有出入，没有考虑精馏塔复杂的流体力学特性和非平衡的传质过程，造成模型的不准确。可选用适用于多股进料及中间出料的多元物系的分段正交配置方法建立精馏塔简化数学模型，或选用考虑非平衡效应和惯性效应的三维非平衡混合动态精馏塔模型，效果应该会好些。在此模型中有很多的求导微分算式。目前应用广泛的求导算法有手工推导法、差商法（finite difference，FD） 、符号微分法（symbolic differentiation,SD） 、逆波兰表示计算法（reverse poli
12、sh notation evalution,RPN）以及自动微分法（automatic differentiation,AD）,其中手工推导法效率低对大规模计算不适用；差商法使用条件宽乏，占用内存少，但精度低，速度也一般；符号微分法精度高，速度快，但是在实际运用中会遇到无法处理的复杂程序代码，符号膨胀，导函数不连续等问题；逆波兰表示计算法是符号微分法的一种特殊情况，其实质是运用堆栈而不是树的数据结构来表示被求导对象；而自动微分法是进来发展迅速的一种求导方法，可以灵活地对复杂程序代码进行求导。应该比本次计算所采用的算法更适合。在采用软件模拟时，也可用选用 ASPEN PLUS 软件中的Mulit
13、ifral 模型。ASPEN PLUS 软件的灵敏度分析也可应用在本次分析中。选用方程中要使其更加严格，汽液相平衡可先用 PRSV 方程，模型求解可采用 Wegestein 法。从模型分析所采取的物系来看，由于是高度理想的二元物系，使得在原来平衡级假设基础上又省去对每个塔板的能量衡算方程，设相对挥发度不随温度改变，省去了温度迭代，而一般情况要对温度进行迭代。因此此模型分析结果很特殊没有代表性。但是，此次是建立动态精馏塔模型，因此基本符合要求。四、优化概述：精馏塔是广泛应用于化工、石油、轻工等行业的汽液传质设备，也是化学工业生产中能源消耗最大的单元操作之一。随着能源价格的不断上涨，精馏系统的节能
14、在化学工业的节能研究和开发中占有越来越重要的地位，有关精馏过程的节能措施一直是人们普遍关注的问题，研究比较多的主要又一下几个方面：对于塔顶塔釜温差较大的精馏塔，在精馏塔中间设冷凝器或再沸器，能节省能量 20左右。采用热泵系统，使塔顶低温处的热量传递给塔釜高温处，也能达到较好的节能效果。采用减压操作，可增大组分之间的相对挥发度，从而减小了塔底再沸器的热负荷。降低精馏塔消耗的途径是多种多样的，无论采取那一种措施，均能获得一定程度的节能效果，但有些措施操作复杂，实施起来有一定难度。实际上，在精馏塔设计使如果能够注意回流比的选择，也能达到较好的分离效果。在精馏塔的设计和操作中，回流比是一个很重要的控制
15、参数，它的选择合理与否，直接影响着精馏塔的操作费用和设备折旧费用。按照本文前面的讨论，精馏塔的塔板数、塔板直径已经给定，也就是该精馏塔已经投产，设备已安装好，塔板数和上升蒸汽量已经给定的情况下，本应从最优控制方面来进行优化，但最优控制过程复杂，计算量大，本文只作定性分析，然后针对精馏塔的设计方面，从系统工程的观点出发，建立经济模型，通过优化方法找出适宜的回流比。经济模型的建立及求解：问题的分析：实用的回流比应在 Rmin 和 R之间选取。可以从经济核算的角度来选定最优回流比 Ropt。精馏费用可从设备投资费和操作费两方面来看。当 R=Rmin，达成分离要求所需的理论板数 N= ，相应的设备费亦为无穷大；当 R 稍稍增大，N 即从无限大急剧减少，随 R 继续增大，R 对 N 的影响逐渐减弱。另方面，随 R的增大，为得到同样的产品 D，精馏段上升蒸汽量 V=(R+1)D 随 R 线性增加，使得再沸器，冷凝器的负荷随之增加，而且塔径也要相应增大。当这些增加的费用超过塔板数减少的费用时，设备费将随 R的增大而增大。因此，随着 R 从 Rmin 起逐渐增大，设备费先是由急剧减小，经过一最小值后又重新增大，如图一中的线 1 示意。操作费主要有：加热蒸汽和冷却水的费用，可称为能源费，它随着热负荷，
