利用AspenProperties作氣液平衡數(shù)據(jù)多模型多參數(shù)多套試驗(yàn)數(shù)據(jù)回歸

氣液平衡實(shí)驗(yàn)數(shù)據(jù)二元相互作用參數(shù)回歸謹(jǐn)以此文獻(xiàn)給曾經(jīng)一起工作生活的同事朋友們-DZ,YSW,XB,CW,ZGW,LCH,HQ,ZXA,QHX,CGC,DYL,WSY,ZMM,

我們?cè)鴵碛泄餐拿?2S，希望將來的某一天這個(gè)logo能夠在藍(lán)天白云下高高飄揚(yáng)。

WJH20110403簡介此文完全依照AspenProperties幫助文檔中的例子翻譯而來，配以必要的補(bǔ)充說明，這個(gè)例子演示了如何利用該軟件作乙酸乙酯-乙醇系統(tǒng)的多套實(shí)驗(yàn)數(shù)據(jù)的多模型的(Wilson、NRTL、UNIQUAC)二元相互作用參數(shù)回歸數(shù)據(jù)來源：40Cand70CdataofMartl,Collect.

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

Data31,339(1986):數(shù)據(jù)在相同文件夾下的Word文檔里，可以用來向AspenProperties拷貝步驟1：選擇單位步驟2：添加組分，注意順序要和我的一致，否則后面就沒法作了步驟3：選取熱力學(xué)方法由于要作多套模型參數(shù)的回歸，因此在referenced頁面將需要回歸的其他模型選上步驟4：數(shù)據(jù)實(shí)驗(yàn)輸入，點(diǎn)New選擇MIXTURE后點(diǎn)擊OK由于我們要同時(shí)回歸等溫等壓數(shù)據(jù)，因此在Setup頁面作如下選擇在Data頁面輸入數(shù)據(jù)，要用拷貝粘貼功能，否則將是一個(gè)痛苦的過程，數(shù)據(jù)在附件中再次點(diǎn)擊Data輸入第二、三套實(shí)驗(yàn)數(shù)據(jù)輸入第二套實(shí)驗(yàn)數(shù)據(jù)輸入第三套實(shí)驗(yàn)數(shù)據(jù)步驟5：建立回歸，點(diǎn)擊Regression->New為Wilson模型回歸添加三套實(shí)驗(yàn)數(shù)據(jù)設(shè)置回歸Wilson的模型參數(shù)，有一些人不明白這里的Element中的1、2到底是什么意思，這里解釋下如果你對(duì)教科書上的公式情有獨(dú)鐘，扯點(diǎn)題外的Aspen也是支持的，例如你想輸入老式參數(shù)，選中你要輸入的列，使之變蘭，然后點(diǎn)擊右上角那個(gè)不起眼的Dechema按鈕，在跳出的窗口中就可以輸入了，做這個(gè)頁面的程序員不知道怎么考慮的，是不是就怕用戶看見呀，放那么偏的位置，我很久都沒有注意到！建立第二套回歸計(jì)算并為NRTL模型回歸添加三套實(shí)驗(yàn)數(shù)據(jù)設(shè)置NRTL模型被回歸的參數(shù)，為了說明圖中Element的作用，我們選1、2、5作為被回歸的參數(shù)，以和Wilosn模型作對(duì)比，當(dāng)然你可以去回歸3、4同樣的方法設(shè)置UNIQUAC模型的回歸參數(shù)UNIQUAC模型我們只回歸Element2選取Regression后就可以回歸了點(diǎn)擊N->開始回歸回歸結(jié)束，選擇YestoAll就用回歸得到的參數(shù)取代了數(shù)據(jù)庫中的參數(shù)回歸結(jié)果點(diǎn)擊Properties->Parameters->Wilson-1、NRTL-1及UNIQ-1頁面可以看到回歸得到的參數(shù)已經(jīng)取代了舊參數(shù)，且Wilson模型回歸了Element1、2，NRTL模型回歸了1、2、5，UNIQUAC模型只回歸了Element2，有一點(diǎn)要注意回歸UNIQUAC模型時(shí)，由于我們未回歸Aij，因此它的值還是原來的數(shù)值，并且Source也是R-R-1,2,3，要用回歸得到UNIQUAC的BIJ參數(shù)，我認(rèn)為必須和原有的AIJ一起使用，因此該參數(shù)也參與了回歸計(jì)算，必須一起使用。Wilson模型的熱力學(xué)一致性檢驗(yàn)結(jié)果殘差分析及作圖

這是Wilson模型第一套實(shí)驗(yàn)數(shù)據(jù)的殘差數(shù)據(jù)點(diǎn)擊菜單中Plot->PlotWizard…選擇殘差分析作圖由于D-1是等溫?cái)?shù)據(jù)，因此我們選壓力作為殘差分析項(xiàng)從結(jié)果可以看出殘差并未均勻分布在零值附近，這有由于回歸是三套實(shí)驗(yàn)同時(shí)回歸，需要照顧三套實(shí)驗(yàn)數(shù)據(jù)現(xiàn)在我們把第一、二套的殘差放在一起分析，不要關(guān)閉剛才的圖，繼續(xù)點(diǎn)擊Plot->PlotWizard…只是在這一步，我們做如下選擇，這樣就把第一、二套的殘差數(shù)據(jù)放在一起了第一、二套實(shí)驗(yàn)數(shù)據(jù)殘差分析，可以看出雖然有部分?jǐn)?shù)據(jù)位于零值上方，但大部分?jǐn)?shù)據(jù)位于零值下方由于第三套實(shí)驗(yàn)數(shù)據(jù)是等壓實(shí)驗(yàn)數(shù)據(jù)，因此需要用溫度作為殘差分析選擇可以看出第三套實(shí)驗(yàn)數(shù)據(jù)的偏差大部分位于零值下方，這和第一、二套的數(shù)據(jù)是相同的，但是你要知道溫度與壓力應(yīng)該具有相反的方向，壓力偏高，通常意味著溫度偏低，因