精馏塔设计计算书

2/2精馏塔设计计算书项目名称：苯-甲苯二元物系常压连续精馏塔设计设计阶段：工艺设计计算编制依据：《化工工艺设计手册》、《精馏设计手册》、《化学工程手册》计算目的：完成精馏塔的物料衡算、最小回流比及操作回流比确定、理论塔板数计算（逐板法及简捷法）、实际塔板数估算、塔径计算、塔高计算、水力学校核，保证塔在操作条件下满足分离要求。一、设计条件与物性数据1.1设计任务项目符号数值单位备注处理量（进料）F5000kg/h苯-甲苯混合物进料组成（苯质量分数）x0.40-换算摩尔分数塔顶馏出液组成（苯质量分数）x0.98-塔底残液组成（苯质量分数）x0.02-操作压力（塔顶）p101.325kPa常压进料热状态-饱和液体进料-泡点进料(q=1)1.2物性数据（苯-甲苯，常压）组分分子量M(g/mol)密度ρ(kg/m³)表面张力σ(mN/m)粘度μ(mPa·s)苯78.1187828.80.56甲苯92.1486628.50.59操作温度：塔顶约80.5℃，塔底约110℃，进料约93℃。1.3相平衡数据（苯-甲苯，101.3kPa）使用Antoine方程计算。Antoine常数（苯：A=6.01907，B=1204.682，C=223.466；甲苯：A=6.08426，B=1347.110，C=219.622，压力单位kPa，温度℃）。本设计采用平均相对挥发度法简化。根据文献，苯-甲苯在常压下平均相对挥发度α≈2.47（取80~110℃平均值）。相平衡方程：y二、物料衡算2.1进料、塔顶、塔底流量的摩尔换算进料总摩尔流量：M其中M=78.11，M=92.14g/mol。分步计算：1.进料中苯的质量分数0.40，甲苯0.60。2.平均摩尔质量：M=1/(0.40/78.11+0.60/92.14)=1/(0.005121+0.006511)=1/0.011632=85.97g/mol。3.进料摩尔流量：F=5000/85.97=58.16kmol/h。进料组成（摩尔分数）：x塔顶组成（摩尔分数）：x塔底组成（摩尔分数）：x2.2全塔物料衡算总物料衡算：F=D+W轻组分衡算：Fx=Dx+Wx代入数值：58.16=58.16×0.440=25.590=解方程组：由第一式：W=58.16-D代入第二式：25.590=0.983D+0.0235(58.16-D)=0.983D+1.366-0.0235D=0.9595D+1.3660.9595D=25.590-1.366=24.224D=24.224/0.9595=25.25kmol/hW=58.16-25.25=32.91kmol/h结论：-塔顶产品流量D=25.25kmol/h，摩尔分数x=0.983-塔底产品流量W=32.91kmol/h，摩尔分数x=0.0235-进料F=58.16kmol/h，x=0.440三、最小回流比与操作回流比3.1泡点进料（q=1）时最小回流比精馏段操作线在q线与平衡线交点处斜率最小。对于泡点进料，交点坐标为(x,y)，其中y由相平衡方程计算：y最小回流比公式（q=1）：R代入：R3.2操作回流比取操作回流比R=1.5R=1.5×1.468=2.20。结论：R=1.47，设计回流比R=2.20。四、理论塔板数计算（逐板法）4.1精馏段操作线方程y代入R=2.20：Rx精馏段方程：y=0.6875x+0.30724.2提馏段操作线方程泡点进料，q=1，提馏段液相流量L'=L+qF=L+F，汽相流量V'=V=(R+1)D。首先计算L：L=RD=2.20×25.25=55.55kmol/h。V=(R+1)D=3.20×25.25=80.80kmol/h。提馏段斜率：L截距（以x表达）：y其中W/V'=32.91/80.80=0.4073。提馏段方程：y=1.408x-0.4073×0.0235=1.408x-0.009574.3逐板计算平衡关系：y=2.47x/(1+1.47x)，反函数x=y/(2.47-1.47y)。从塔顶开始（全凝器，y=x=0.983）。第1块板（塔顶）：x=y/(2.47-1.47×y)=0.983/(2.47-1.47×0.983)=0.983/(2.47-1.445)=0.983/1.025=0.959。第2块板（精馏段）：y=0.6875×x+0.3072=0.6875×0.959+0.3072=0.659+0.3072=0.9662。x=0.9662/(2.47-1.47×0.9662)=0.9662/(2.47-1.420)=0.9662/1.050=0.920。第3块板：y=0.6875×0.920+0.3072=0.6325+0.3072=0.9397。x=0.9397/(2.47-1.47×0.9397)=0.9397/(2.47-1.382)=0.9397/1.088=0.864。第4块板：y=0.6875×0.864+0.3072=0.594+0.3072=0.9012。x=0.9012/(2.47-1.47×0.9012)=0.9012/(2.47-1.325)=0.9012/1.145=0.787。第5块板：y=0.6875×0.787+0.3072=0.541+0.3072=0.8482。x=0.8482/(2.47-1.47×0.8482)=0.8482/(2.47-1.247)=0.8482/1.223=0.694。第6块板：y=0.6875×0.694+0.3072=0.477+0.3072=0.7842。x=0.7842/(2.47-1.47×0.7842)=0.7842/(2.47-1.153)=0.7842/1.317=0.595。第7块板：y=0.6875×0.595+0.3072=0.409+0.3072=0.7162。x=0.7162/(2.47-1.47×0.7162)=0.7162/(2.47-1.053)=0.7162/1.417=0.505。因x=0.505>x=0.440，继续精馏段。第8块板：y=0.6875×0.505+0.3072=0.347+0.3072=0.6542。x=0.6542/(2.47-1.47×0.6542)=0.6542/(2.47-0.9617)=0.6542/1.5083=0.434。此时x=0.434<x=0.440，可认为进料板在第8块（从上往下数），第9块开始换为提馏段操作线。转入提馏段：从x=0.434（精馏段最后一块板液相组成）开始计算提馏段汽相组成。提馏段第一块板（即第9块）上升汽相y对应x由提馏段方程计算：y=1.408×x-0.00957=1.408×0.434-0.00957=0.611-0.00957=0.6014。然后由平衡求x：x=0.6014/(2.47-1.47×0.6014)=0.6014/(2.47-0.884)=0.6014/1.586=0.379。第10块板：y=1.408×0.379-0.00957=0.534-0.00957=0.5244。x=0.5244/(2.47-1.47×0.5244)=0.5244/(2.47-0.771)=0.5244/1.699=0.309。第11块板：y=1.408×0.309-0.00957=0.435-0.00957=0.4254。x=0.4254/(2.47-1.47×0.4254)=0.4254/(2.47-0.625)=0.4254/1.845=0.231。第12块板：y=1.408×0.231-0.00957=0.325-0.00957=0.3154。x=0.3154/(2.47-1.47×0.3154)=0.3154/(2.47-0.464)=0.3154/2.006=0.157。第13块板：y=1.408×0.157-0.00957=0.221-0.00957=0.2115。x=0.2115/(2.47-1.47×0.2115)=0.2115/(2.47-0.311)=0.2115/2.159=0.0980。第14块板：y=1.408×0.0980-0.00957=0.138-0.00957=0.1284。x=0.1284/(2.47-1.47×0.1284)=0.1284/(2.47-0.189)=0.1284/2.281=0.0563。第15块板：y=1.408×0.0563-0.00957=0.0793-0.00957=0.0697。x=0.0697/(2.47-1.47×0.0697)=0.0697/(2.47-0.1025)=0.0697/2.3675=0.0295。第16块板：y=1.408×0.0295-0.00957=0.0415-0.00957=0.0319。x=0.0319/(2.47-1.47×0.0319)=0.0319/(2.47-0.0469)=0.0319/2.4231=0.0132。此时x=0.0132<x=0.0235，说明塔釜已满足要求。计算中塔底再沸器相当于一块理论板，所以总理论板数包括再沸器。理论板数汇总：-精馏段理论板数：7块（从第1到第7板，因第8板已低于进料组成，通常进料板算作精馏段最后一块或提馏段第一块。通常将进料板作为提馏段第一块，则精馏段为第1~7块，共7块）-提馏段理论板数：从进料板（第8块）到塔底，共16-7=9块（含再沸器）-总理论塔板数（不含再沸器）：15块-含再沸器：16块结论：理论塔板数（含再沸器）N=16，其中精馏段7板，提馏段9板（含再沸器）。五、实际塔板数估算采用奥康奈尔全塔效率关联式：E其中α为相对挥发度（取2.47），μ为进料液相平均粘度（mPa·s）。计算进料粘度（摩尔分数基准）：进料x=0.440，苯粘度0.56，甲苯0.59，则：μαE求1.422=e=e=e=0.917。E实际塔板数：N或含再沸器不计，塔板数为35块。结论：实际塔板数N=35块（不含再沸器）。六、塔径计算6.1气液相体积流量精馏段：汽相摩尔流率V=(R+1)D=3.20×25.25=80.80kmol/h。液相摩尔流率L=RD=2.20×25.25=55.55kmol/h。操作温度（塔顶~塔底平均约90℃），压力常压。计算汽相密度：平均分子量M=xM+(1-x)M=0.983×78.11+0.017×92.14=76.78+1.57=78.35g/mol。汽相密度ρ=(pM)/(RT)，p=101.3kPa，T=273+90=363K，R=8.314kJ/(kmol·K)ρ液相密度：90℃时苯密度约813kg/m³，甲苯密度约794kg/m³（近似，因温度影响）。按塔顶组成取液相密度近似：苯质量分数0.98，则ρ=1/(0.98/813+0.02/794)=1/(0.001205+0.0000252)=1/0.0012302=813kg/m³。汽相体积流量：V液相体积流量：L液相平均分子量M=xM+(1-x)M=0.983×78.11+0.017×92.14=78.35g/mol（相同）。L6.2泛点气速计算采用Smith关联图法。先计算L/V·(ρ/ρ)：L√309.1=17.58，乘积0.002226×17.58=0.0391。查Smith图（或经验公式），取板间距H=0.45m，液泛因子C=0.082（根据图估）。表面张力σ约28.6mN/m，校正：C1.43=e=e=e=1.074，C=0.082×1.074=0.0881。泛点气速：u√308.1=17.55，u=0.0881×17.55=1.546m/s。取操作气速为泛点气速的70%：u6.3塔截面积与塔径所需截面积：A塔径：D圆整至D=0.9m。校核：实际截面积A=π×(0.9)/4=0.636m²，实际空塔气速u=0.668/0.636=1.050m/s，为泛点气速的68%，合理。结论：塔径D=900mm。七、塔高计算7.1塔体有效高度板间距H=0.45m，实际塔板数35块，则有效塔高（不含裙座及开孔）：H7.2附加高度-进料段高度：H=1.0m-塔顶空间（含除沫器）：H=1.2m-塔底空间（液位+再沸器入口）：H=2.5m-裙座高度：H=2.0m总塔高：H八、塔盘水力学校核（F1型浮阀塔为例）8.1堰高与堰上液层高度取堰长l=0.7D=0.63m。液相流量L=0.001487m³/s，L/l=0.001487/0.63。0.63=0.63×0.63=0.3969×0.7937=0.315，则L/l=0.001487/0.315=0.00472。堰上液层高度（Francis公式）：hL国际上常用公式为h=0.00284×(L/l)，但单位需注意。采用标准单位（m³/s，m）：h(0.002361)=e=e=e=0.0178。h=0.00284×0.0178=0.0000506m=0.05mm，明显错误。更正：Francis公式通常采用h=0.00284(3600L/l)其中L为m³/h。L(m³/h)=0.001487×3600=5.353m³/hLhln8.50=2.140,2.140×0.6667=1.427,e=4.17，则h=0.00284×4.17=0.01184m=11.8mm。取堰高h=50mm，则清液层高度h=h+h=50+11.8=61.8mm。8.2雾沫夹带与降液管负荷计算泛点率（已控制68%），雾沫夹带可接受。降液管停留时间：降液管停留时间校核：取降液管宽度W_d=0.12D=0.12×0.9=0.108m，降液管截面积A_d=W_d×H_T（板间距）=0.108×0.45=0.0486m²。降液管内液体体积流量与精馏段液相体积流量一致，即L_l=0.001487m³/s。降液管停留时间τ=A_d/L_l=0.0486/0.001487≈32.7s。工业上要求降液管停留时间τ≥3~5s（浮阀塔通常要求≥5s），此处τ=32.7s，满足设计要求，可有效避免液体在降液管内夹带气泡，保证液流稳定。雾沫夹带校核：雾沫夹带量e_V是衡量塔盘操作稳定性的重要指标，工业上要求e_V≤0.1kg（液）/kg（气）。采用经验公式计算雾沫夹带量：e其中，u为空塔气速（1.050m/s），H_T为板间距（0.45m），h_{cl}为清液层高度（0.0618m）。代入数值计算：H_T-h_{cl}=0.45-0.0618=0.3882mfracuee_V=0.0915kg/kg<0.1kg/kg，满足设计要求，雾沫夹带量在允许范围内，不会因雾沫夹带导致分离效率下降。8.3气体通过浮阀的压力降F1型浮阀塔的气体通过浮阀的压力降Δp由干阀压力降Δp_d和液层压力降Δp_l两部分组成，总压力降Δp_total=Δp_d+Δp_l。干阀压力降计算：采用经验公式：Delta其中，u_0为气体通过浮阀孔的气速（m/s），ρ_V为汽相密度（2.63kg/m³），g为重力加速度（9.81m/s²）。取浮阀孔径d_0=39mm（标准F1浮阀），每块塔盘浮阀数N=120个（根据塔径0.9m估算，浮阀中心距约100mm），浮阀总开孔面积A_0=N×(πd_0²/4)=120×(3.1416×0.039²/4)≈120×0.00119≈0.1428m²。气体通过浮阀孔的气速u_0=V_g/A_0=0.668/0.1428≈4.68m/s。代入公式计算干阀压力降：Delta液层压力降计算：液层压力降与清液层高度h_{cl}相关，经验公式为：Delta总压力降：Delta浮阀塔的总压力降通常控制在0.1~0.3m液柱，此处Δp_total≈0.198m液柱，符合设计要求，不会因压力降过大导致塔内操作不稳定或能耗过高。8.4液体在塔盘上的停留时间液体在塔盘上的停留时间t_L反映液体与气体的接触时间，接触时间不足会影响传质效率，通常要求t_L≥2~5s。塔盘有效面积A_active=A-2A_d=0.636-2×0.0486=0.5388m²（扣除两侧降液管面积）。液体在塔盘上的平均流速u_L=L_l/A_active=0.001487/0.53