1、第五章第五章 多组分多级分离的多组分多级分离的严格计算严格计算 掌握掌握掌握掌握掌握掌握第五章第五章 主要内容主要内容5.1多组分复杂精馏过程的模型多组分复杂精馏过程的模型5.2三对角矩阵法三对角矩阵法5.3逐板计算法逐板计算法5.4松弛法松弛法了解了解本章要求本章要求思考题思考题v1）了解平衡级的基本假定，掌握多级平）了解平衡级的基本假定，掌握多级平衡分离过程的物理模型和数学模型衡分离过程的物理模型和数学模型v2）掌握逐板计算法原理、计算起点和步）掌握逐板计算法原理、计算起点和步骤骤v3）掌握三对角线矩阵法（泡点法）的计）掌握三对角线矩阵法（泡点法）的计算原理和解法，该法的优缺点和改进方法算
2、原理和解法，该法的优缺点和改进方法v4）了解流率加和法（）了解流率加和法（SR法）和松驰法的法）和松驰法的计算原理计算原理本章要求：本章要求：5.1平衡级的理论模型平衡级的理论模型 5.1.1 复杂精馏塔物理模型复杂精馏塔物理模型v对于多于一股进料和多于两股出料对于多于一股进料和多于两股出料的精馏，称为复杂精馏的精馏，称为复杂精馏v采用复杂精馏进行分离是为了采用复杂精馏进行分离是为了节省节省能量能量和和减少设备减少设备的数量的数量 5.1.1.1复杂精馏塔类型复杂精馏塔类型多股进料多股进料侧线采出侧线采出中间冷凝或中间再沸中间冷凝或中间再沸 5.1.1.2模型塔模型塔j平衡级平衡级jijjjj
3、HyVPT,jijjjjhxLPT,11,111,jjijjjhxLPT11,111,jjijjjHyVPTFjjijFjFjHZFPT,jGjUjQ平衡级平衡级NG1NUrQCQ5.1.2平衡级的理论模型 5.1.2.1多级分离过程的平衡级多级分离过程的平衡级 (1)在每一分离级上的每一相流体都是在每一分离级上的每一相流体都是完全完全混合混合的，级上的温度、压力和组成均一致，的，级上的温度、压力和组成均一致，且与离开分离级的该相流体相同且与离开分离级的该相流体相同 (2)离开分离级的两相流体之间成离开分离级的两相流体之间成相平衡相平衡 具备这两个条件的分离级就是平衡级具备这两个条件的分离级就
4、是平衡级 由平衡级假设引起的误差，可以用由平衡级假设引起的误差，可以用级效级效率率进行修正进行修正5.1.2.25.1.2.2多级分离过程的数学模型多级分离过程的数学模型MESH方程组方程组 物料平衡式物料平衡式 Material Balance Equation ijjjijijjjijjjjijzFyVxULyGVxL1,11,1)()(ijijijxKy 1ijx1ijy 热量平衡式热量平衡式 Heat Balance Equation 或或摩尔分率加和式摩尔分率加和式 Summary Equation相平衡关系式相平衡关系式Equilibrium Balance EquationFjj
5、jjjjjjjjjjHFHVhULHGVhL1111)()(jQ 除除MESH模型方程组外模型方程组外, jjijhHK,),(ijijjjijijyxpTKK ),(ijjjjjxpThh ),(ijjjjjypTHH) 32(CNNuc 1)93(CNNuv个变量个变量个方程和个方程和将上述将上述N个平衡级按逆流方式个平衡级按逆流方式串联串联起来，有起来，有必须知道必须知道的关联式的关联式5.1.2.35.1.2.3分离过程的变量分析分离过程的变量分析 v设计变量总数 )32(cNNucvMESH方程v装置变量总数 1)93(cNNuv1)6(cNNuiv固定设计变量 )3(cNv可调设计
6、变量 13Nv其中串级单元数为1，侧线采出单元数为2（N1），传热单元数为N。v设计型，操作型5.1.3模拟计算方法5.1.3.1模拟计算方法模拟计算方法 开发前的准备开发前的准备 构成一个精馏塔模拟计算的算法。必须对构成一个精馏塔模拟计算的算法。必须对下列下列三点作出选择和安排三点作出选择和安排 1）迭代变量的选择）迭代变量的选择 2）迭代变量的组织）迭代变量的组织 3）一些变量的圆整和归一的方法以及迭代）一些变量的圆整和归一的方法以及迭代的加速方法的加速方法 5.1.3.2严格计算法的种类方程解离法方程解离法v一类是逐级计算法，将一类是逐级计算法，将MESH方程组方程组按平衡按平衡级分组级
7、分组，从塔两端，从塔两端逐级逐级求解，即逐级求解相求解，即逐级求解相平衡，物料平衡和热平衡方程式；平衡，物料平衡和热平衡方程式；适用于清适用于清晰分割晰分割。v另一类是逐次替代法，将另一类是逐次替代法，将MESH方程方程按类别按类别组合组合，对其中，对其中一类一类和几类方程组用和几类方程组用矩阵法对矩阵法对各级同时求解。各级同时求解。 同时校正法同时校正法v利用某种迭代技术求解利用某种迭代技术求解全部全部或或大部大部MESH方程。应方程。应先将先将MESH方程组方程组线性化线性化，后由某种同时校正，后由某种同时校正(SC)法法模拟模拟各类多组元各类多组元精馏分离过程精馏分离过程，用迭代方法，用
8、迭代方法(如如Newton法法)对其同时求解。对其同时求解。适用于非理想体系适用于非理想体系。非稳态方程计算方法非稳态方程计算方法v采用采用不稳定不稳定状态的物料平衡方程和热平衡方程，求状态的物料平衡方程和热平衡方程，求解解稳定稳定状态下多级平衡过程状态下多级平衡过程 5.1.3.3计算类型计算类型设计型计算。设计型计算。其目的在于解决完成一预定其目的在于解决完成一预定的分离任务的新过程设计问题的分离任务的新过程设计问题v需知道轻、重关键组分的回收率，求解需知道轻、重关键组分的回收率，求解所需所需的的理论板数理论板数，和最佳，和最佳进料位置进料位置和和侧线采出位侧线采出位置置操作型计算。操作型
9、计算。已知操作条件下，分析和考已知操作条件下，分析和考察已有的分离设备的性能察已有的分离设备的性能 5.1.3.4典型精馏计算法优缺点比较典型精馏计算法优缺点比较 算法类别优点缺点适用范围逐级法可用于设计型计算需估计塔顶、塔釜组成；截断误差和累计误差大；难于在复杂塔中应用三对角矩阵法用于操作型计算，没有逐级法的缺点算法包含泡（露）点计算，因它也是迭代过程，费机时。对高度非理想系难收敛BP法算法简单，对初值要求不高，占用内存单元少当临近解时出现收敛速度慢，迭代次数，多当趋于解时，收敛缓慢精馏SR法算法中不包含泡（露）点计算，收敛速度快级数多时不稳定，不适用于窄沸程精馏系统适用吸收塔及宽沸程精馏系
10、统BP-SR法同BP法只能用于吸收精馏吸收精馏2N牛顿-拉夫森法收敛速度快，稳定性好不包含泡露点计算计算工作量大，占用内存单元多，对初值要求高弱非理想体系的各种场合多元牛顿-拉夫森法收敛速度快，稳定性好， 计算雅可比矩阵耗时长，计算工作量大，占用内存单元多，对初值要求高各种场合松弛法对初值要求不高，稳定性好，适用范围宽收敛速度极慢，不宜用作常规算法各种场合5.2.三对角矩阵法三对角矩阵法 5.2.1计算方法和原理计算方法和原理v在初步假定的沿塔高温度、汽、液流量的情况下，在初步假定的沿塔高温度、汽、液流量的情况下，v逐板地用物料平衡（逐板地用物料平衡（MM）和汽液平衡（）和汽液平衡（E E）方
11、程联立）方程联立求得一组方程，并用矩阵求解各板上求得一组方程，并用矩阵求解各板上组成组成v用用S S方程求各板上新的方程求各板上新的温度温度v用用H H方程方程求各板上新的求各板上新的汽液流量汽液流量v如此循环计算直到如此循环计算直到稳定稳定为止为止 泡点法泡点法(BP法法) v初步假定的沿塔高温度、汽、液流量的情况下初步假定的沿塔高温度、汽、液流量的情况下v逐板地用物料平衡（逐板地用物料平衡（MM）和汽液平衡（）和汽液平衡（E E）方程联）方程联立求得一组方程，并用矩阵求解各板上立求得一组方程，并用矩阵求解各板上组成组成x xi,ji,j。 v用用S S方程方程求各板上新的求各板上新的温度（
12、泡点方程）温度（泡点方程）v用用H H方程方程求各板上新的求各板上新的汽液流量汽液流量v如此循环计算直到如此循环计算直到稳定稳定为止为止 v主要应用于窄沸程混合物的精馏过程。主要应用于窄沸程混合物的精馏过程。 露点法露点法(DP法法) v初步假定的沿塔高温度、汽、液流量的情况下初步假定的沿塔高温度、汽、液流量的情况下v逐板地用物料平衡逐板地用物料平衡（MM）和汽液平衡（和汽液平衡（E E）方程联方程联立求得一组方程，并用矩阵求解各板上立求得一组方程，并用矩阵求解各板上组成组成y yi,ji,j。 v用用S S方程方程求各板上新的求各板上新的温度温度( (露点方程露点方程）v用用H H方程方程求
13、各板上新的求各板上新的汽液流量汽液流量v如此循环计算直到如此循环计算直到稳定稳定为止为止 v分离含氢较多的混合物时常用。分离含氢较多的混合物时常用。 流量加和法流量加和法 (SR法法) v先由MEME方程方程求出组分流率v vi,ji,j（或（或l li,ji,j），），后求出总流率V Vj j，v再由H H方程方程求T Tj j；v应用于吸收塔、解吸塔和萃取塔。矩阵求逆法矩阵求逆法 v初步假定的沿塔高温度、汽、液流量的情况下初步假定的沿塔高温度、汽、液流量的情况下v采用采用矩阵求逆矩阵求逆的方法的方法解解ME方程方程求出新的液相组成求出新的液相组成xi,j或气相组成或气相组成yi,j。 v用
14、用S方程方程求各板上新的求各板上新的温度温度(露点方程）露点方程）v用用H方程方程求各板上新的求各板上新的汽液流量汽液流量v如此循环计算直到如此循环计算直到稳定稳定为止为止 v分离含氢较多的混合物时常用。分离含氢较多的混合物时常用。同时校正法同时校正法(SC法法) v通过某种迭代技术通过某种迭代技术求解全部求解全部或或大部大部分分MESH方程方程或与之等价的方程式。或与之等价的方程式。vSC法适用于法适用于非理想非理想性很强的液体混合物的精性很强的液体混合物的精馏过程，如萃取精馏和共沸精榴。馏过程，如萃取精馏和共沸精榴。vSC法还适用于带有化学反应的分离过程的计法还适用于带有化学反应的分离过程
15、的计算，如反应精馏和催化精馏等。算，如反应精馏和催化精馏等。5.2.2泡点法泡点法(BP法法) 5.2.2.1解三对角矩阵方程解三对角矩阵方程(ME方程方程)求组分分布求组分分布 ME方程方程将将E方程带入方程带入M方程方程消去消去 1jjLA1,1jijjKVCijjjzFDjjijijjjijDxCxBxA1,1,令令ijijjjijijjjjijijjijFzxULxKGVxKVxL)()(1,1,11,1ijjjijjijjjijjjjijzFxKVxULKGVxL1,111,1)()()()(jjijjjjULKGVBijy 当当001 LA)(11111ULKVBi221iKVC
16、01D1j1NNLA)(BKVBiNNN011NNNKVC0NDNj 12111BDxCxiijjijijjjijDxCxBxA1,1,NNiNNiNDxBxA,1,1j12NjNj 时，即塔釜，无下一板来蒸汽时，即塔釜，无下一板来蒸汽VN+10。当当的液体的液体时，即塔顶冷凝器，由于没有上一板来时，即塔顶冷凝器，由于没有上一板来或简写为或简写为 NNjNiNiijiiNNNNNjjjDDDDDxxxxxBACBACBACBACB121,1,2111122211 jijDxCBA,初值的确定初值的确定 jT) 1)(1(jNTTTTDBDj初jVjjjjjjjULGVLVF11jjjjjjjF
17、LULGVV11jjjjULqFL1)1 (1jjjjjqFGVV12) 1(LDDRV其中：RDL 1设为恒摩尔流设为恒摩尔流j板平衡：板平衡：液相平衡：液相平衡：根据塔顶、塔釜的温度线性分布根据塔顶、塔釜的温度线性分布 jLjVDGULFVjkkjkkjjkkj2221jkkkkjjDUGFVL21)(jjjjDCBA, 由由j板与塔顶作物料平衡板与塔顶作物料平衡常数中如果某些物料没有常数中如果某些物料没有可以用零代入可以用零代入求求由由三对角矩阵中求解三对角矩阵中求解 的方法的方法jixNNjiNNiijiiNjqqqqqxxxxxPPPP1211,21121101010101 其中其中
18、 1jjjjjPABCP11jjjjjjjPABqADq12111,111,qxpxqxpxqxpxqxiijjijijNiNNNiNiNjPjq由此求出各由此求出各和和并可以求出某一组分在各块板上的液相组成并可以求出某一组分在各块板上的液相组成同理求出所有组分在各块板的液相组成各块板的液相组成5.2.2.2用用S方程计算新的温度分布方程计算新的温度分布v归一处理归一处理v1）圆整）圆整NTTNjkjkjT01. 0)()(121ijijijxxxv用泡点法求用泡点法求 ，并同时得，并同时得 jTijy 结束满足精度要求否方程 Y1,kjjisjikjTyxTN重设温度重设温度Tj5.2.2.
19、3用用H方程计算各块板的方程计算各块板的 和和 jVjLjjjjjjjFjjjjjjQhULHGVHFHVhL)()(1111)(11jjjjjjjGVFVLULH方程：方程：任一板总物料衡算任一板总物料衡算jjjjjFjjjjjjjjjhHQFhHLhhGVhHV1111)()()(由假定的由假定的DUGFVLKjKKKjj)(211V1jVNVDRLDVRDLUDV) 1(,12111为止为止冷凝器开始，然后随着冷凝器开始，然后随着j的递增而求得的递增而求得，计算顺序从，计算顺序从初始值即可求得初始值即可求得01. 0)()()(1kjkjkjHVVV5.2.2.4 热负荷计算热负荷计算
20、v冷凝器11111122hLhUHVHVQcv再沸器11)(NNNNNNNrhLhLHGVQ5.2.2.5计算步骤计算步骤 确定必要条件和基础数据确定必要条件和基础数据 按塔顶、塔釜的温度假定在塔内温度为线按塔顶、塔釜的温度假定在塔内温度为线性分布的温度初始值性分布的温度初始值jTjVjTijKjjjjjjqPDCBA,ijx1ijx则圆整则圆整 ，若，若 用高斯消去法解矩阵得用高斯消去法解矩阵得中的中的，然后计算，然后计算ME矩阵方程矩阵方程的计算的计算由假设由假设组初始的蒸汽量分布组初始的蒸汽量分布，按恒摩尔流假定一，按恒摩尔流假定一 由计算出的由计算出的ijxjTijyijxijyjTj
21、jhH ,jjLV,NTTkjkjT01. 0)()(2101. 0)()()(1kjkjkjHVVV 判断是否满足收敛条件判断是否满足收敛条件新的汽液相流量新的汽液相流量 用用H方程从冷凝器开始向下计算出各板的方程从冷凝器开始向下计算出各板的计算计算、由由，同时计算，同时计算，用，用S方程试差迭代出新方程试差迭代出新的的温度温度v例5-1 某精馏塔进料流量为100kmol/h,，进料中有三个组分，其摩尔分数分别为0.3、0.4和0.3，进料位置如图(5-5)所示。塔顶馏出液流量为50kmolh，回流比为1，泡点进料。若按恒摩尔流的假定取Vj初值，且假定各级各组分的K值如表5-2所示。试用托马
22、斯法求初次迭代所得的组成断面 。 jix,jix,v解：因进料为泡点，按恒摩尔流的假定有：解：因进料为泡点，按恒摩尔流的假定有： hkmolUVVL505001001121将将Vj代入式（代入式（5-21）同样，可得各级的液相流率同样，可得各级的液相流率Lj如表所示。如表所示。hkmolDRVVV1005021234级序GjUjVjLjK1,jK2,jK3,j1050050210.520010050210.5300100150210.540010050210.5150,50,4321AAALAjjDUGFVLKjKKKjj)(21v由进料流量及组成，求得各级各组分的由进料流量及组成，求得各级各
23、组分的 ，结果见表结果见表 jGjUjVjLjK, 1 jK,2 jK, 3jijjizFD,级序组分的Di,j123100020003-30-40-304000v托马斯法，按步骤计算如下：托马斯法，按步骤计算如下：v对组分对组分1（j1）：）： vj2 1005050111 , 111ULKVB2)100/(2100/12, 12111BKVBCP01000/11 , 11BDQ250502100)(222, 1222ULKGVB3333. 1)2(5025021001123 , 1312222PLBKVPABCP02502500500122122, 12PABQADQ5012 LAj3vj
24、43501502100)(333 , 1333ULKGVB7059. 0)3333. 1(5035021002234, 1423333PLBKVPABCP10588. 0)3333. 1(5035005030233233 , 13PABQADQ250502100)(444, 1444ULKGVB1102. 0)7059. 0(15025010588. 01500344344, 144, 1PABQADQxv回代得：5023 LA15034 LAv同法对二、三两组分进行计算，结果见表同法对二、三两组分进行计算，结果见表 1837. 01102. 0)7059. 0(10588. 04, 1333
25、 , 1xPQx2449. 01837. 0)3333. 1(03 , 1222, 1xPQx4898. 02449. 0)2(02, 1111 , 1xPQx级序组分xi,j合计12310.48980.40.10910.998920.24490.40.21820.863130.18370.40.32730.911040.11020.40.49091.0011NNjiNNiijqqqqqxxxxx1211,12112501500020035050002002505000200100jixNNjiNNiijiiNjqqqqqxxxxxPPPP1211,211211010101015.2.3 法收
26、敛和法收敛和C.M.B.矩阵法矩阵法 为克服三对角矩阵法的缺点，在归1办法上作了改进。v5.2.3.1法收敛法收敛 法归一取代简单归一而进行的计算称为精馏塔的法收敛。法收敛包括三对角矩阵法求xij，对xij的法归一，加和方程求各级温度Tj及焓方程求各级流率。v5.2.3.2 C.M.B.矩阵法矩阵法 v用物料衡算来校正圆整后的液相组成，将顶、底产品和各板的圆整组成予以校正，使之使之不仅满足S方程，也尽量符合M方程，然后进行泡点计算提高矩阵法的稳定性和收敛速度。该法称为C.M.B. 矩阵法。 5.2.4 流量加和法（流量加和法（SR法）和矩阵求逆法法）和矩阵求逆法5.2.4.1 流量加和法流量加
27、和法（SR法法）解解MEME三对角矩阵方程求出组分三对角矩阵方程求出组分流率流率ijlijvjLjVjTijxijlv相应的衡算式的形式也要有所改变相应的衡算式的形式也要有所改变v即用各板之组分流率代替即用各板之组分流率代替v用这一顺序为计算方便，独立变量取用这一顺序为计算方便，独立变量取v再用再用H方程校正温度方程校正温度 或或和和v组分流率加和组分流率加和(S)得到得到5.2.4.2 矩阵求逆法矩阵求逆法 v采用矩阵求逆的方法求解各板液相组成和汽相流量v应用矩阵求逆法的优点是便于处理不按常规流程操作的精馏塔，如设中间冷凝器和中间再沸器的塔。对这类采用特殊操作方案的塔，其ME方程组并不能表示
28、成三对角矩阵的形式。 v5.2.5Tj 和和Vj同时校正法（多元同时校正法（多元NewtonRaphson法）法）5.2.6.讨论讨论（1）由子程序求解）由子程序求解ME三对角矩阵的三对角矩阵的 ，若若为为负值负值均置为零均置为零v（2）对接近）对接近理想理想溶液的物系，采用溶液的物系，采用大循环大循环计计算速度快算速度快v对对非理想非理想溶液物系，采用溶液物系，采用T循环循环有利收敛有利收敛ijx（3）对接近理想溶液体系，对）对接近理想溶液体系，对 初值要初值要求不苛刻求不苛刻 jjVT ,对非理想物系，对非理想物系， 初值初值会影响到是否收会影响到是否收敛和收敛速度。敛和收敛速度。 jjV
29、T ,v大循环大循环vT循环循环（4）收敛判据）收敛判据不唯一不唯一（5）BP法与法与SR法比较法比较 vBP法法是以求得的浓度来决定下一循环所用的是以求得的浓度来决定下一循环所用的温度。即适用于塔板温度。即适用于塔板温度温度主要决定于浓度主要决定于浓度(组组成成），而），而流率流率主要决定于主要决定于热平衡热平衡 vSR法法用在组成对流率的影响大于焓平衡对流用在组成对流率的影响大于焓平衡对流率的影响，率的影响，温度温度决定于决定于热量平衡热量平衡而不是组成，而不是组成，即热平衡中显热的影响较为显著的情况即热平衡中显热的影响较为显著的情况（6）BP法的法的缺陷缺陷v分离程度分离程度要求高要求高
30、时，收敛时，收敛速度速度明显明显，后期收敛速，后期收敛速度变慢；度变慢；vBP法用于非理想性较强的物系，往往不收敛法用于非理想性较强的物系，往往不收敛v原因：原因：对初值有要求对初值有要求v由由M方程得到的液相组成在没有收敛前，不满足方程得到的液相组成在没有收敛前，不满足 的关系式的关系式xxi,ji,j=1=1；v若将这些组成直接返回第二步运算，结果会发散若将这些组成直接返回第二步运算，结果会发散v若将得到的值圆整，若将得到的值圆整， xi,j不再满足不再满足M方程方程 。5.3逐板计算法逐板计算法5.3.1 可调设计变量、计算内容与计算起点可调设计变量、计算内容与计算起点5.3.1.15.
31、3.1.1可调设计变量可调设计变量v按按设计型设计型规定设计变量规定设计变量（1 1）饱和液体回流；）饱和液体回流； （2 2）LKLK的回收率；的回收率； （3 3） HKHK的回收率的回收率（4 4）R R大于大于RminRmin；（5 5）适宜的进料位置）适宜的进料位置 DxLK，D L0 V L F，Zi，F V LN-1 N W xHK,W开发思维：二元精馏开发思维：二元精馏图解梯级法图解梯级法。 5.3.1.25.3.1.2计算内容计算内容v以塔顶或塔釜为计算起点，根据以塔顶或塔釜为计算起点，根据MESHMESH关系，关系，逐板计算出各板的条件及满足关键组分分离逐板计算出各板的条件
32、及满足关键组分分离要求所需的要求所需的N Nv恒摩尔流恒摩尔流：求：求 yxTN,用用MESMES方程联立求解方程联立求解 变摩尔变摩尔流：求流：求 LVyxTN,用用MESHMESH方程联立求解方程联立求解 5.3.1.3计算起点计算起点 起点选择应从起点选择应从组成较精确组成较精确的一端算起的一端算起 分离的混合物中除分离的混合物中除关键组分关键组分外，仅有外，仅有轻组分轻组分存在存在,（两个（两个难挥发难挥发组分作关键组分），则以组分作关键组分），则以塔顶塔顶为起点为起点向下逐板计算，各板的温度由向下逐板计算，各板的温度由露点方程露点方程求取求取。除除关键组分关键组分外仅有外仅有重组分重
33、组分存在时，（两个存在时，（两个易挥发易挥发组组分作关键组分），则以分作关键组分），则以塔釜塔釜作起点向上逐级计算，作起点向上逐级计算，各板的温度由各板的温度由泡点方程泡点方程求取求取关键组分是关键组分是中间组分中间组分，可以从两端同时算起在加料，可以从两端同时算起在加料板处契合板处契合 5.3.2恒摩尔流计算方法恒摩尔流计算方法5.3.2.1从塔从塔顶向下顶向下计算，塔序从塔顶向下计算，塔序从塔顶向下 （1）MES方程方程 MM方程方程：精馏段：精馏段 DininninDxxLyV,1,1DinninnnixVDxVLy,1,11,恒摩尔流：恒摩尔流： ,) 1(,DRVRDLR由由Unde
34、rwood法法求得求得 提馏段提馏段: WimmimmmixVWxVLy,1,11,E E方程方程： jijijijijijijiyyKyx,S S方程方程： 1,jix（2）计算步骤）计算步骤 露点），分凝（SEyDi, 精MtxDDi,Diy,全凝全凝 露点，SEyi 1 ,精，Mxti 1 ,1加料板加料板 niMy,提 塔釜塔釜 5.3.2.2从塔釜向上计算，塔序从塔釜向上从塔釜向上计算，塔序从塔釜向上 (1)MES(1)MES方程方程MM方程方程：提馏段：提馏段 WimmimmmixLWyLVx,1,11,精馏段精馏段 DinninnnixLDyLVx,1,11,E E方程方程： j
35、irijirijijijixxxKy,S S方程方程： 1,jiy(2)(2)计算步骤计算步骤 泡点，SExWi, 提MtyWWi,泡点，SExi 1 ,提，Myti 1 ,1加料板加料板 niMx, 精 塔顶塔顶 5.3.2.35.3.2.3从两头往中间算从两头往中间算 v方法同上，但由塔顶向下计算到加料板的组方法同上，但由塔顶向下计算到加料板的组成，与由塔釜向上计算到加料板处的组成一成，与由塔釜向上计算到加料板处的组成一般是不吻合的般是不吻合的v根据物料分布，根据物料分布，塔顶塔顶产品中往往不含产品中往往不含重组分重组分，由上向下计算时，应在适当部位加入这些组由上向下计算时，应在适当部位加
36、入这些组分，才能保证在进料板处得到吻合。分，才能保证在进料板处得到吻合。v由由塔釜塔釜向上计算时，要在适当部位加入一些向上计算时，要在适当部位加入一些轻组分轻组分 5.3.35.3.3进料板的确定和计算结束的判断进料板的确定和计算结束的判断v适宜的进料位置定义为达到规定分离要适宜的进料位置定义为达到规定分离要求所需总级数求所需总级数最少最少的进料位置的进料位置v近似确定方法是以近似确定方法是以轻、重轻、重关键组分的关键组分的浓浓度之比度之比作为精馏效果的准则作为精馏效果的准则v改变改变操作线操作线效果效果比不改变比不改变更好更好时，换操时，换操作线方程作线方程 5.3.3.1 。从塔釜向上计算
37、。从塔釜向上计算 v为使分离效果好，要求轻、重关键组分液相为使分离效果好，要求轻、重关键组分液相浓度比值浓度比值增加增加越快越好越快越好 v如果如果 则第则第j j极为适宜进料位置，极为适宜进料位置，j+1j+1级应转换成级应转换成用精馏段操作线计算用精馏段操作线计算 计算计算结束判据结束判据： DHLNHLxxxxSjHLRjHLxxxx,SjHLRjHLxxxx, 1, 1且5.3.3.25.3.3.2从塔顶向下计算从塔顶向下计算 v要求轻、重关键组分汽相浓度比值要求轻、重关键组分汽相浓度比值降低降低的的越快越好越快越好 v如果如果 则第则第j j极为适宜进料位置，极为适宜进料位置，j+1
38、j+1级应转换成用级应转换成用提馏段操作线计算提馏段操作线计算 计算计算结束判据结束判据： WHLNHLxxxxSjHLRjHLyyyy,SjHLRjHLyyyy, 1, 1且且5.3.3.3从两头向中间计算（近似法）从两头向中间计算（近似法）nnnHLFHLnHLxxxxxx11取（取（n+1n+1）板为进料）板为进料板 11mmmHLFHLmHLxxxxxx取（取（m+1m+1）板为进料板）板为进料板 mnN5.3.3.45.3.3.4计算结果的校核计算结果的校核 按清晰分割处理往往存在一定误差，如从上按清晰分割处理往往存在一定误差，如从上往下计算时，逐板计算至第往下计算时，逐板计算至第N
39、 N板得板得x xLNKLNK，N N，由此值，由此值结合结合LKLK、HKHK回收率和全塔物料衡算得回收率和全塔物料衡算得x xLNKLNK，D D，将，将此值与前面估计值比较，若不能满足计算要求，此值与前面估计值比较，若不能满足计算要求，则需校核，可将计算得的则需校核，可将计算得的x xLNK,DLNK,D作为下一次的估作为下一次的估计值计值 。01. 0,一般取计算计算估计DLNKDLNKDLNKxxx校核：校核：5.3.4、变摩尔流的逐板法、变摩尔流的逐板法变摩尔流，则各板的变摩尔流，则各板的 nnnnLLVV11,精馏段热量平衡式：精馏段热量平衡式： cDnnnnqDhhLHV11n
40、ncnDnhHqhhDV11)(由塔顶向下计算至由塔顶向下计算至进料板进料板，从而求得各板，从而求得各板上的新流量上的新流量 LV,提馏段提馏段热量衡算式热量衡算式 wmmrmmhWHVqhL11)()(11mmrwmmHhqhhWV由釜向上逐板计算直至由釜向上逐板计算直至加料板加料板为止，即可为止，即可计算提馏段各板的计算提馏段各板的 mmVL 及1111nnnnVHHyMV方程热平操设设调整调整v例例5-4 见书。例题见书。例题5.4非稳态方程计算方法非稳态方程计算方法 （松弛法）松弛法） 开发思维开发思维: 精馏塔的开车精馏塔的开车 5.4.1松弛法的基本方程松弛法的基本方程v由开始的由
41、开始的非稳定态非稳定态向稳定态变化的过程中，对某向稳定态变化的过程中，对某一时间间隔内每块塔板上的物料变化进行衡算一时间间隔内每块塔板上的物料变化进行衡算 v假定塔内为恒摩尔流率，每块板均为理论板，并假定塔内为恒摩尔流率，每块板均为理论板，并假定无侧线采出和热量交换，则在不稳定态时，假定无侧线采出和热量交换，则在不稳定态时，在每个时间间隔内必有累积现象发生在每个时间间隔内必有累积现象发生v累积量累积量=输入量输入量输出量输出量 对第对第j板上的组分进行物料衡算板上的组分进行物料衡算 dyVxLjijjij)(1,11,1avjijjijyVxL)(1,11,1dxLyVijjijj)(avij
42、jijjxLyV)(组分组分i离开离开j板的量板的量=j平衡级平衡级ijjyV ,ijjxL ,1,1,jijxL1,1,jijyV组分组分i进入进入j板的量板的量= 组分组分i在在j板上的累积量板上的累积量=|ijjijjxExEdxEdijj)(ddExddxEjijijjjjjWUEjUjWddExddxExLyVyVxLjijijjavijjijjjijjij1,11,1jE00ddEjddxExLyVyVxLijjijjijjjijjij1,11,1上式变为瞬时值上式变为瞬时值 时，时，不随时间变化，且不随时间变化，且设设为为j板上的存气量折算成存液量板上的存气量折算成存液量 其中其
43、中为为j板上的存液量板上的存液量 当当yxVL,ddxijddxxxijijij)()(ijxddxddxijij关系即关系即的变化率为线的变化率为线性性随时间随时间Rose等假定等假定时的组成时的组成 由拉格朗日中值定理求由拉格朗日中值定理求j板上板上i组分在组分在的变化率的变化率时各板上组分时各板上组分i随时间随时间即可由上式算出在即可由上式算出在均为已知，均为已知，时，全塔各板上的时，全塔各板上的 F为进料量，若各板的存液量相等，各板的为进料量，若各板的存液量相等，各板的松弛系数也可取等值，但对不同组分有时可松弛系数也可取等值，但对不同组分有时可用不同的松弛系数值用不同的松弛系数值kij
44、jijjjijjijjkijkijxLyVyVxLExx1,11,11)()(jjEFEj)105(1称为松弛系数称为松弛系数松弛法基本方程松弛法基本方程5.4.2 计算步骤计算步骤 假定一组初值。假定一组初值。 选选 FEjj)105(1ijxiFijxx ijyjVjL1)(kijx11)(|)()( |kijkijkijxxx 直至稳定为止直至稳定为止判断判断由松弛法基本方程计算新的由松弛法基本方程计算新的（设为恒摩尔流）（设为恒摩尔流） 或或设一组各板的设一组各板的），泡点试差得），泡点试差得（设（设设一组各板的设一组各板的 由由S方程求方程求 11)( ,)(kijkjyT1)(ki
45、jx11)( ,)(kijkjyTjjhH ,11)( ,)(kjkjLVvkjkjkjVVV11)(|)()( |不成立则重复不成立则重复式式 判断判断 由由H方程计算方程计算计算计算，由由松弛法的优缺点松弛法的优缺点 优点优点: 适用于各种复杂的精馏过程，对初值适用于各种复杂的精馏过程，对初值选定没有严格要求，且某一组分的量发生变选定没有严格要求，且某一组分的量发生变化时，对计算影响不大，所以收敛很稳定化时，对计算影响不大，所以收敛很稳定 缺点：缺点：必须使用计算机，且收敛速度也太慢，必须使用计算机，且收敛速度也太慢，当其他方法无法解时，才用松弛法当其他方法无法解时，才用松弛法 改进改进:
46、 石川、平田等人进行了修正，但收敛石川、平田等人进行了修正，但收敛速度仍较慢速度仍较慢v由丁惠华提出的新松弛法则在由丁惠华提出的新松弛法则在Rose松弛法的松弛法的基础上，经过适当处理允许采用较大的松弛基础上，经过适当处理允许采用较大的松弛系数，克服了收敛速度慢的缺点，可用于非系数，克服了收敛速度慢的缺点，可用于非理想物系的精馏计算。理想物系的精馏计算。 5.4.3 模拟计算方法的改进模拟计算方法的改进v5.4.3.1松弛法与松弛法与N-R法的联合算法法的联合算法v5.4.3.2 自由度自由度N(C十十2)同时校正法同时校正法 v5.4.3.3 二对角矩阵法二对角矩阵法 v5.4.3.4 非平
47、衡级模型非平衡级模型 v1 1、三对角矩阵法（泡点法）的计算原理和计、三对角矩阵法（泡点法）的计算原理和计算框图，并简单说明该法在应用上的局限性，算框图，并简单说明该法在应用上的局限性，即如何改进即如何改进v2 2、写出三对角矩阵法的、写出三对角矩阵法的MESHMESH方程，并说明方程，并说明流量加合法和三对角矩阵法有何不同流量加合法和三对角矩阵法有何不同v3 3、开发精馏过程的新算法需对那几点作出选、开发精馏过程的新算法需对那几点作出选择和安排择和安排, ,以三对角矩阵法为例说明以三对角矩阵法为例说明4 4、在严格的精馏计算中，应采用何种方法进行物料、在严格的精馏计算中，应采用何种方法进行物
48、料预分布，其物料衡算的原理是什么？预分布，其物料衡算的原理是什么？5 5、逐板计算法的计算起点如何选择，怎样确定适宜、逐板计算法的计算起点如何选择，怎样确定适宜的加料位置，叙述从塔釜向上逐板计算的步骤。的加料位置，叙述从塔釜向上逐板计算的步骤。6 6、松弛法的开发思路是什么？具有何种物理意义？、松弛法的开发思路是什么？具有何种物理意义？其基本方程式什么？其基本方程式什么？7 7、松弛法的基本思想是什么？列出松弛法的基本方、松弛法的基本思想是什么？列出松弛法的基本方程。程。 学习本章要学习本章要从复杂塔的物理模型推导理论板的物理模型从复杂塔的物理模型推导理论板的物理模型,由此推导由此推导MESH
49、方程即数学模型，逐板计算法重点讲解恒方程即数学模型，逐板计算法重点讲解恒摩尔流的摩尔流的MES方程方程,以及计算步骤，了解进料板的确定方法。以及计算步骤，了解进料板的确定方法。从三对角矩阵法的计算原理出发从三对角矩阵法的计算原理出发,推导推导ME方程方程,说明该法迭说明该法迭代变量的选择和组织的方法，说明该法的缺点和用代变量的选择和组织的方法，说明该法的缺点和用CMB矩矩阵改进的方法，了解松弛法的计算原理和基本方程。阵改进的方法，了解松弛法的计算原理和基本方程。重点：重点：逐级计算方法；托玛斯法解三对角矩阵方程；精馏问逐级计算方法；托玛斯法解三对角矩阵方程；精馏问题的求解。题的求解。 难点：难
50、点：多级分离过程的多级分离过程的MESH数学模型，三对角矩阵方程托数学模型，三对角矩阵方程托玛斯解法计算精馏过程的泡点温度。玛斯解法计算精馏过程的泡点温度。 21. 0, 0 . 1, 5 . 2CTTTBT 习题精馏塔内各组分的相对挥发度 ，进料板上升蒸汽的组成yB=0.35(mol)，yT=0.20，yC=0.45，回流比为1.7，饱和液体回流，进料板上一板下流液体组成为xB=0.24(mol)，xT=0.18，xC=0.58。求进料板以上第二板的上升蒸汽组成。解：对精馏段作物料衡算 总物料总物料 V=L+DV=L+D对对i i组分组分 DxLxVyDijiji,1,因为因为R=1.7R=
