利用AspenProperties作气液平衡数据多模型多参数多套试验数据回归

气液平衡实验数据二元相互作用参数回归谨以此文献给曾经一起工作生活的同事朋友们-DZ,YSW,XB,CW,ZGW,LCH,HQ,ZXA,QHX,CGC,DYL,WSY,ZMM,

我们曾拥有共同的名字-2S，希望将来的某一天这个logo能够在蓝天白云下高高飘扬。

WJH20110403简介此文完全依照AspenProperties帮助文档中的例子翻译而来，配以必要的补充说明，这个例子演示了如何利用该软件作乙酸乙酯-乙醇系统的多套实验数据的多模型的(Wilson、NRTL、UNIQUAC)二元相互作用参数回归数据来源：40Cand70CdataofMartl,Collect.

Czech.Chem.Commun.37,266(1972):AtmosphericdataofOrtegaJ.andPenaJ.A.,J.Chem.Eng.

Data31,339(1986):数据在相同文件夹下的Word文档里，可以用来向AspenProperties拷贝步骤1：选择单位步骤2：添加组分，注意顺序要和我的一致，否则后面就没法作了步骤3：选取热力学方法由于要作多套模型参数的回归，因此在referenced页面将需要回归的其他模型选上步骤4：数据实验输入，点New选择MIXTURE后点击OK由于我们要同时回归等温等压数据，因此在Setup页面作如下选择在Data页面输入数据，要用拷贝粘贴功能，否则将是一个痛苦的过程，数据在附件中再次点击Data输入第二、三套实验数据输入第二套实验数据输入第三套实验数据步骤5：建立回归，点击Regression->New为Wilson模型回归添加三套实验数据设置回归Wilson的模型参数，有一些人不明白这里的Element中的1、2到底是什么意思，这里解释下如果你对教科书上的公式情有独钟，扯点题外的Aspen也是支持的，例如你想输入老式参数，选中你要输入的列，使之变兰，然后点击右上角那个不起眼的Dechema按钮，在跳出的窗口中就可以输入了，做这个页面的程序员不知道怎么考虑的，是不是就怕用户看见呀，放那么偏的位置，我很久都没有注意到！建立第二套回归计算并为NRTL模型回归添加三套实验数据设置NRTL模型被回归的参数，为了说明图中Element的作用，我们选1、2、5作为被回归的参数，以和Wilosn模型作对比，当然你可以去回归3、4同样的方法设置UNIQUAC模型的回归参数UNIQUAC模型我们只回归Element2选取Regression后就可以回归了点击N->开始回归回归结束，选择YestoAll就用回归得到的参数取代了数据库中的参数回归结果点击Properties->Parameters->Wilson-1、NRTL-1及UNIQ-1页面可以看到回归得到的参数已经取代了旧参数，且Wilson模型回归了Element1、2，NRTL模型回归了1、2、5，UNIQUAC模型只回归了Element2，有一点要注意回归UNIQUAC模型时，由于我们未回归Aij，因此它的值还是原来的数值，并且Source也是R-R-1,2,3，要用回归得到UNIQUAC的BIJ参数，我认为必须和原有的AIJ一起使用，因此该参数也参与了回归计算，必须一起使用。Wilson模型的热力学一致性检验结果残差分析及作图

这是Wilson模型第一套实验数据的残差数据点击菜单中Plot->PlotWizard…选择残差分析作图由于D-1是等温数据，因此我们选压力作为残差分析项从结果可以看出残差并未均匀分布在零值附近，这有由于回归是三套实验同时回归，需要照顾三套实验数据现在我们把第一、二套的残差放在一起分析，不要关闭刚才的图，继续点击Plot->PlotWizard…只是在这一步，我们做如下选择，这样就把第一、二套的残差数据放在一起了第一、二套实验数据残差分析，可以看出虽然有部分数据位于零值上方，但大部分数据位于零值下方由于第三套实验数据是等压实验数据，因此需要用温度作为残差分析选择可以看出第三套实验数据的偏差大部分位于零值下方，这和第一、二套的数据是相同的，但是你要知道温度与压力应该具有相反的方向，压力偏高，通常意味着温度偏低，因