《利用enProperties作气液平衡数据多模型多参数多套试验数据回归.ppt》由会员分享,可在线阅读,更多相关《利用enProperties作气液平衡数据多模型多参数多套试验数据回归.ppt(42页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、气液平衡实验数据二元相互作用参数回归谨以此文献给曾经一起工作生活的同事朋友们-DZ,YSW,XB,CW,ZGW,LCH,HQ,ZXA,QHX,CGC,DYL,WSY,ZMM,我们曾拥有共同的名字-2S,希望将来的某一天这个logo能够在蓝天白云下高高飘扬。WJH 2011 04 03简介此文完全依照AspenProperties帮助文档中的例子翻译而来,配以必要的补充说明,这个例子演示了如何利用该软件作乙酸乙酯-乙醇系统的多套实验数据的多模型的(Wilson、NRTL、UNIQUAC)二元相互作用参数回归数据来源:40C and 70C data of Martl,Collect.Czech.
2、Chem.Commun.37,266(1972):Atmospheric data of Ortega J.and Pena J.A.,J.Chem.Eng.Data 31,339(1986):数据在相同文件夹下的Word文档里,可以用来向AspenProperties拷贝步骤1:选择单位步骤2:添加组分,注意顺序要和我的一致,否则后面就没法作了步骤3:选取热力学方法由于要作多套模型参数的回归,因此在referenced页面将需要回归的其他模型选上步骤4:数据实验输入,点New选择MIXTURE后点击OK由于我们要同时回归等温等压数据,因此在Setup页面作如下选择在Data页面输入数据,要用
3、拷贝粘贴功能,否则将是一个痛苦的过程,数据在附件中再次点击Data输入第二、三套实验数据输入第二套实验数据输入第三套实验数据步骤5:建立回归,点击Regression-New为Wilson模型回归添加三套实验数据设置回归Wilson的模型参数,有一些人不明白这里的Element中的1、2到底是什么意思,这里解释下 如果你对教科书上的公式情有独钟,扯点题外的如果你对教科书上的公式情有独钟,扯点题外的AspenAspen也是支持的,例如你想输入老式也是支持的,例如你想输入老式参数,选中你要输入的列,使之变兰,然后点击右上角那个不起眼的参数,选中你要输入的列,使之变兰,然后点击右上角那个不起眼的De
4、chemaDechema按钮,按钮,在跳出的窗口中就可以输入了,做这个页面的程序员不知道怎么考虑的,是不是就怕在跳出的窗口中就可以输入了,做这个页面的程序员不知道怎么考虑的,是不是就怕用户看见呀,放那么偏的位置,我很久都没有注意到!用户看见呀,放那么偏的位置,我很久都没有注意到!建立第二套回归计算并为NRTL模型回归添加三套实验数据设置设置NRTLNRTL模型被回归的参数,为了说明图中模型被回归的参数,为了说明图中ElementElement的作用,的作用,我们选我们选1 1、2 2、5 5作为被回归的参数,以和作为被回归的参数,以和WilosnWilosn模型作对比,模型作对比,当然你可以去
5、回归当然你可以去回归3 3、4 4同样的方法设置UNIQUAC模型的回归参数UNIQUAC模型我们只回归Element2选取Regression后就可以回归了点击N-开始回归回归结束,选择Yes to All就用回归得到的参数取代了数据库中的参数回归结果点击点击Properties-Parameters-Wilson-1Properties-Parameters-Wilson-1、NRTL-1NRTL-1及及UNIQ-1UNIQ-1页面可以看到回归得到的页面可以看到回归得到的参数已经取代了旧参数,且参数已经取代了旧参数,且WilsonWilson模型回归了模型回归了Element 1Eleme
6、nt 1、2 2,NRTLNRTL模型回归了模型回归了1 1、2 2、5 5,UNIQUACUNIQUAC模型只回归了模型只回归了Element2Element2,有一点要注意回归,有一点要注意回归UNIQUACUNIQUAC模型时,由于我们模型时,由于我们未回归未回归AijAij,因此它的值还是原来的数值,并且,因此它的值还是原来的数值,并且SourceSource也是也是R-R-1,2,3R-R-1,2,3,要用,要用回归得到回归得到回归得到回归得到UNIQUACUNIQUAC的的BIJBIJ参数,我认为必须和原有的参数,我认为必须和原有的AIJAIJ一起使用,因此该参数也参与了回归计一起
7、使用,因此该参数也参与了回归计算,必须一起使用。算,必须一起使用。Wilson模型的热力学一致性检验结果残差分析及作图这是这是WilsonWilson模型第一套实验数据的残差数据模型第一套实验数据的残差数据点击菜单中Plot-Plot Wizard选择残差分析作图由于D-1是等温数据,因此我们选压力作为残差分析项从结果可以看出残差并未均匀分布在零值附近,这从结果可以看出残差并未均匀分布在零值附近,这有由于回归是三套实验同时回归,需要照顾三套实有由于回归是三套实验同时回归,需要照顾三套实验数据验数据现在我们把第一、二套的残差放在一起分析,不要关闭刚才现在我们把第一、二套的残差放在一起分析,不要关
8、闭刚才的图,继续点击的图,继续点击Plot-Plot WizardPlot-Plot Wizard只是在这一步,我们做只是在这一步,我们做如下选择,这样就把第一、二套的残差数据放在一起了如下选择,这样就把第一、二套的残差数据放在一起了第一、二套实验数据残差分析,可以看出虽然有部分数据位于零值上方,但大部分数据位于零值下方由于第三套实验数据是等压实验数据,因此需要用温度作为残差分析选择可以看出第三套实验数据的偏差大部分位于零值下方,这和可以看出第三套实验数据的偏差大部分位于零值下方,这和第一、二套的数据是相同的,但是你要知道温度与压力应该第一、二套的数据是相同的,但是你要知道温度与压力应该具有相反的方向,压力偏高,通常意味着温度偏低,因此如具有相反的方向,压力偏高,通常意味着温度偏低,因此如果这三套结合在一起,大概会接近于零吧(这一点是我的猜果这三套结合在一起,大概会接近于零吧(这一点是我的猜测)测)用同样的方法做P-xy图,即泡点露点压力图选择第一套数据作图可以看出此物系存在共沸点,有最高共沸压力,同样也有最低共沸温度The End