1、实验二二元气液平衡数据的测定实验二二元气液平衡数据的测定一、 实验目的1.测定苯一正庚烷二元体系在常压下的气液平衡数据。2.通过实验了解平衡釜的结构，掌握气液平衡数据的测定方法和技能。3.应用Wilson方程关联实验数据。二、 实验原理气液平衡数据是化学工业发展新产品、开发新工艺、减少能耗、进行 三废处理的重要基础数据之一。化工生产中的蒸馏和吸收等分离过程设备 的改造与设计、挖潜与革新以及对最佳工艺条件的选择，都需要精确可靠 的气液平衡数据。这是因为化工生产过程都要涉及相间的物质传递， 故这种数据的重要性是显而易见的。随着化工生产的不断发展，现有气液平衡 数据远不能满足需要。许多物系的平衡数据
2、，很难由理论直接计算得到， 必须由实验测定。平衡数据实验测定方法有两类，即间接法和直接法。直接法中又有静 态法、流动法和循环法等。其中循环法应用最为广泛。若要测得准确的气 液平衡数据，平衡釜是关键。现已采用的平衡釜形式有多种， 而且各有特 点，应根据待测物系的特征，选择适当的釜型。用常规的平衡釜测定平衡数据，需样品量多，测定时间长。本实验用 的小型平衡釜主要特点是釜外有真空夹套保温，釜内液体和气体分别形成 循环系统，可观察釜内的实验现象，且样品用量少，达到平衡速度快，因 而实验时间短。以循环法测定气液平衡数据的平衡釜类型虽多，但基本原理相同，如图 1 所示。当体系达到平衡时，两个容器的组成不随
3、时间变化，这时从A和B 两容器中取样分析，即可得到一组平衡数据。当达到平衡时，除了两相的压力和温度分别相等外， 每一组分的化学位也相等，即逸度相等，其热力学基本关系为：ft = /! (1)e 2刃= f凶常压下，气相可视为理想气体，巾=1 ；再忽略压力对液体逸度的影响，fi=Pi0从而得出低压下气液平衡关系式为：pyt = rifiQxl 式中，p体系压力(总压)：PL纯组分i在平衡温度下饱和蒸气压，可用安托尼(itoine)公式计算:忌、系分别为组分i在液相和气相中的摩尔分率：组分I的活度系数。由实验测得等压下气液平衡数据，则可用: 计算出不同组成下的活度系数。In斤=一陨比+ A12.v
4、2 ） + x2 （屮儿?2 _七+；珂 本实验中活度系数和组成关系采用 Wilson方程关联。Wilson方程为：血尸2 =1 hl（ .丫2 + A21.Y| ） + 比（*2+1；円-X+J二,） （5）Wilson方程二元配偶参数 A12和 筠 采用非线性最小二乘法，由二元 气液平衡数据回归而得。目标函数选为气相组成误差的平方和，即尸=工（片实_加）;+（旳实一厂计）; 尸1三、实验装置和试剂1.平衡釜一台。平衡釜的选择原则：易于建立平衡、样品用量少、平衡温度测定准确、气相中不夹带液滴、液相不返混及不易爆沸等。本实验用 气液双循环的小平衡釜，其结构如图 2所示。2.阿贝折射仪一台。3.
5、超级恒温槽一台4.50100十分之一的标准温度计一支、050十分之一的标准温度计一支。5.所用试剂（苯、正庚烷）为优级品。6.1mL注射器十支、5mL注射器三支。n图2小P液平衡釜示总图1磨口； 2气扣取样口： 3气相贮液槽；斗一连通管：5缓冲球；6回流管；7平I8钟罩；温度计套管s 10液根取样口； 11液和贮液僧；H提升管；13沸腾室；1斗一加热套管：15真空变会；16加料液面四、实验步骤1. 首先开启超级恒温槽，调节温度至测定折射率所需温度 30C。2.测温套管中道入甘油，将标准温度计插入套管中，并在温度计露出部 分中间固定一支温度计。因本实验对温度要求较严，需对温度进行校正。3.检查整
6、个系统的气密性，因为我们这个系统都是密闭的。检测方法是将100毫升针筒与系统相连，并使系统与大气隔绝，针筒缓缓抽出一点压 力，发现硅油U型管的两个液柱差不变时（说明系统是密闭的），然后 再通大气。4. 我们做的是常压下的气液平衡，需读出当天实验时的大气压值 （一个标 准大气压）。5.在平衡釜内加入一定浓度的苯 一正庚烷混合液约2030ml。打开冷却 水，安放好加热器，接通电源。控制加热电流，开始时给 0.1A, 5min后 给0.2A，慢慢调到0.25A左右即可，以平衡釜内液体沸腾为准。冷凝回流液控制在每秒23滴。观察温度计示数，两次度数，示数不变即达到平衡。6.达到平衡后，需要记录下两个温度
7、计的读数，此温度为平衡温度，并 用微量注射器分别取两相样品，通过阿贝折射仪测其组成。将加热电流调 为0,釜液停止沸腾后拿下加热器，待温度降下来，换液。7.根据第一次平衡时测得的温度和混合液组成确定应该加入的组分为苯。 用注射器从釜中取出3ml的混合液，然后加入4ml苯，重新建立平衡。均 匀测量4组平衡数据。8.实验完毕，关掉电源和水源。五、实验数据记录和处理(1 )实验数据记录：1号装置日期：2015年5月16日室内温度：29.9C大气压：101.3Kpa1.实验用纯组分在常压下的沸点：苯：80.1 C正庚烷：984C2.表1实验数据记录表t主/Ct辅/C气相折光率气相组成/%液相折光率液相组
8、成/%192.7932.01.397625.601.395623.00287.7931.81.417548.831.405635.65385.9831.71.424454.671.412543.68484.2131.91.435663.411.421652.483.表2温度计示数修正表温度示值/ c505560657075修正值/c-0.02-0.04+0.01-0.03-0.04-0.01温度示值/ c80859095100修正值/c+0.010.00+0.020.00-0.04根据示值计算实际温度的公式：实际温度 二示值+修正值4.表3纯物料的安托尼常数表组分ABCt范围（C）苯6.879
9、871196.76219.1611584正庚烷6.893861264.37216.643 127（2）实验数据处理1平衡温度的计算平衡温度计算方法如下：t实际= t主+1 修正 +1 校正值t修正：通过精查温度计的修正值得到，见表 2t校正值：t校正值 =kn （t 主-ts）k 水银在玻璃中间的膨胀系数取 0.00016;n露出部分的温度系数取 60Ct主在介质中的温度； ts露出水银柱的平均温度（即辅助温度计的读数）t修正 由Matlab2013程序得到：t0=84.21 85.98 87.79 92.79;x=50 55 60 65 70 75 80 85 90 95 100y二-0.0
10、2 -0.04 0.01 -0.03 -0.04 -0.01 0.01 0.00 0.02 0.00 -0.04z=spli ne(x,y,to)z(tx)二-0.0001 0.0025 0.0113 0.0146ts=31.9 31.7 31.8 32.0tj = kn (t0-ts)=0.00016*60*(t0-ts)t=t0+tx+tj输出 t=84.7121 86.5063 88.3388 93.3881以第一组数据为例：t实际= t主+1 修正 +1 校正值t 主=92.79C t 辅=32.0C t 修正=0.0146t 校正值二kn(t 主-t 辅)=0.00016*60*(9
11、2.79-32)=0.5836 Ct 实际=t 主 +t 修正 +t 校正值=92.79C +0.5836C +0.0146C =93.3881 C表4平衡温度计算表t主/ Ct辅/Ct修正/ Ct实际/ C192.7932.0-0.000184.7121287.7931.80.002586.5063385.9831.70.011388.3388484.2131.90.014693.3881注：表中修正温度由主温度在表1所给温度修正值在三次样条插值而得，最终实际温 度计算结果保留一位小数，是因为主温度测定的精度决定的2.苯和正庚烷纯组分在本实验不同温度下的饱和蒸汽压的计算Antoine (安托
12、尼)公式：lgPi0=Ai Bi/(Ci+t); Pi=10(Ai-Bi/(Ci+t)式中：Po为饱和蒸汽压，mmHg； t为实际温度，见表4; A、B、C均为常 数；纯物料的Antoine常数见表3以第一组数据为例：t=934C苯的安托尼常数：A=6.87987 B=1196.76 C=219.161Pa=10八(A-B/(C+t)=10八(6.87987-1196.76/(219.161+93.4)=1124.2000mmHg 正庚烷的安托尼常数：A=6.89386 B=1264.37 C=216.64Pb0=0A(A-B/(C+t)=1QA(6.89386-1264.37/(216.64
13、+93.4)=654.0588mmHg表5纯组分饱和蒸汽压计算结果847C86.5 C88.3C93.4C苯 Pa0/ mmHg874.0081921.7182972.66891124.2000正庚烷Pb0/ mmHg499.1179528.4433559.8819654.05883用非线形最小二乘法回归配偶参数，并求液相组成1活度系数和组成关系采用Wilson方程关联In iln(x上 12X2) X2(12 一 21 )为 + 馆2x2 xA21x1In 2 =Tn(X2 上21X1)洛( 21 )x2 +A21x-! x A12x2目标函数选为气相组成误差的平方和，即m f八(%实-%计
14、)2 厲实-y2计)2j丄用matlab非线性最小二乘法拟合：fun cti on xLs qnon li nbb0=1,1bb,res no rm,residual=ls qnon li n(Funl v,bbO)fun cti on F=F unlv (bb)x1= 0.230 0.3565 0.4368 0.5248;x2=1-x1;y1= 0.256 0.4883 0.5467 0.6341;y2=1-y1;Pa0=1124.2 972.6689 921.7182 874.0081;Pb0二654.0588 559.8899 528.4433 499.1179;for i=1:4F(i)二y1(i)-Pa0(i)*x1(i)/760*exp(-log(x1(i)+bb(1)*x2(i)+x2(i)*(bb(1)/(x1(i)+bb(1)*x2(i)-bb(2)/(x2(i)+bb(2)*x1(i)endfor i=5:8j=i-4F(i)=y2(j)-Pb0(j)*x2(
