论坛首页 › 『Aspen Plus交流』 › Aspen入门实践 实例1 苯-噻吩体系的NRTL模型二元交互作用参数的回归

Mellapak_tju 发表于 2013-5-29 05:55:10

Aspen入门实践 实例1 苯-噻吩体系的NRTL模型二元交互作用参数的回归

本贴以348.15K的苯-噻吩体系pxy数据回归NRTL模型交互作用参数为例，说明使用Aspen Plus回归模型参数的个人看法与观点，欢迎大家交流
计划长期发关于Aspen的帖子，和大家分享我过去在天大和目前在德国埃尔朗根纽伦堡大学学习Aspen的经验，希望多多交流，也算是为我这学期的这门课复习啦~
附件中有关于该实例的文本和bkp文档，还有Aspen Plus帮助中对NRTL模型的说明，以及我所使用的VLE数据的excel文档。
本科做毕业设计粗苯精制的时候，有尝试过对实验数据进行回归，但现在看来当时的问题还真不少，但估计现在的这篇也会存在一些问题，欢迎指正。

zy2013207188 发表于 2013-5-29 05:55:10

bkp文件下载下来 打开发现，properties-parameters-binary interaction-NRTL-1页面的NRTL参数的TLOWER和TUPPER明显不对。是否应该在properties-regression-R-1-input页面的parameter选项卡中设置Name/Element为NRTL 7和NRTL8的Usage为fix，initial value为348.15K。运行了一下，交互参数结果是一样的，只是改一下看起来会更规范吧。

tdl522 发表于 2013-5-29 05:55:10

1. aji = 0.20 for hydrocarbons and polar, nonassociated species (e.g., n-heptane/acetone).
2. aji = 0.30 for nonpolar compounds (e.g., benzene/ n-heptane), except fluorocarbons and paraffins; nonpolar and polar, nonassociated species (e.g., benzene/acetone); polar species that exhibit negative deviations from Raoult’s law (e.g.,acetone/chloroform) and moderate
positive deviations (e.g., ethanol/water); mixtures of water and polar nonassociated species (e.g., water/acetone).
3. aji = 0.40 for saturated hydrocarbons and homolog perfluorocarbons (e.g., n-hexane/perfluoro-n-hexane).
4. aji = 0.47 for alcohols or other strongly self-associated species with nonpolar species (e.g., ethanol/benzene);carbon tetrachloride with either acetonitrile or nitromethane;
water with either butyl glycol or pyridine.

zy2013207188 发表于 2013-5-29 05:55:10

你好，不管温度设置或不设置为恒定，得到的交互参数结果都是一样的（至少对于这个例子是这样的）。其实，我也不清楚参数回归具体的计算过程，你说的不等式约束优化什么的我也没太听懂。但如果Aspen参数回归计算过程中涉及到温度的话，既然实验数据来自一个温度下，我觉得就应该用这个温度来计算吧。在不设置温度恒定的情况下，Aspen有点不太智能，它在properties-parameters-binary interaction-NRTL-1页面的NRTL参数的TLOWER和TUPPER显示的不是实验数据的温度，好像是默认情况下的温度范围0-1000，看起来有点别扭哦。哈哈，说了这么多，不知道你听懂我在说什么了吗。其实这个问题我觉得没什么意思。

Mellapak_tju 发表于 2013-5-29 05:55:10

zy2013207188 发表于 2014-12-29 11:01
bkp文件下载下来 打开发现，properties-parameters-binary interaction-NRTL-1页面的NRTL参数的TLOWER和TUP ...

然而我一直觉得温度上下限的设置没有必要规定为定值，虽然初值的好坏是对回归和优化有影响，但是我们常规的情况下应该不会出现太大的问题。即使不规定Tupper=Tlower=Texp，但是基于实验数据的回归，所回归的参数势必都会适应于实验温度下的pxy，我是这样认为的。我觉得温度上下限的设置类似于优化问题中的不等式约束，影响应该不大。此外我觉得我这个例子中使用pxy+给定constant的温度值的set-up我认为也是合理的。回头我再试试几个例子

竹影清风 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2014-11-7 22:03
哦，关于这个呀，我可能当时没有给出这个图例来，很多人都问过我，这是aspen自带的一个热力学数据库，在 ...

您好，非常感谢您的耐心解答，您真的对Aspen研究的很细致，我研究方向也是ASPEN方向，希望可以多多向您学习。
我按照您的方法尝试了一下，可是我好像只能选择一个组分，并不能选择两个组分呢，是不是因为版本问题呀？我是Aspen V7.2。
我很想了解二元系的气液平衡数据如何得到的，因为我课题中也是要回归二元相互作用参数，文献中我几乎没找到我需要组分的气液平衡数据，还需要您再解答一下呀！{:1106_398:}

Mellapak_tju 发表于 2013-5-29 05:55:10

竹影清风 发表于 2014-11-7 20:41
请教楼主你编写的说明中“选择Retrieve TDE binary data（搜索热力学数据引擎的二元数据），会发现Aspen Pl ...

哦，关于这个呀，我可能当时没有给出这个图例来，很多人都问过我，这是aspen自带的一个热力学数据库，在工具栏上有个写着TDE的上红下蓝的按钮就能直接打开，这是免费的，你点开后会出现一个对话框，可以选择纯组分或是二元体系，和他挨着的是一个烧瓶状按钮，那是DECHMA的数据库，要花钱的。

Mellapak_tju 发表于 2013-5-29 05:55:10

zy2013207188 发表于 2014-12-29 11:01
bkp文件下载下来 打开发现，properties-parameters-binary interaction-NRTL-1页面的NRTL参数的TLOWER和TUP ...

但我一直觉得温度上下限的设置没有必要规定为定制，因为这是基于实验数据的回归，所以如果不出太大问题，所回归的参数势必都会适应于实验温度下的pxy，我是这样认为的，我觉得温度上下限的设置类似于优化问题中的不等式约束，影响应该不大。回头我再试试几个例子

Mellapak_tju 发表于 2013-5-29 05:55:10

sh_ssp 发表于 2013-7-28 11:29
请问高手，“c. 选择Retrieve TDE binary data（搜索热力学数据引擎的二元数据）”从何处找啊？

有的破解版的比如我的V7.1就没有Detherm和NIST的二元数据库，我这里的截图是在学校机房做的，他的版本里这些数据是全的。如果你的版本中找不到这一步，你可以试试点击Detherm或是NIST的那个按钮查找，在工具栏上，Detherm的按钮像个小烧瓶，NIST的一般在他旁边

Mellapak_tju 发表于 2013-5-29 05:55:10

竹影清风 发表于 2014-11-8 19:56
您好，非常感谢您的耐心解答，您真的对Aspen研究的很细致，我研究方向也是ASPEN方向，希望可以多多向您学 ...

实在找不到数据，或是现有的数据很贵，自己实验测试VLE数据也不失为一个很好的办法，当然貌似实验测数据也并非一件容易的实习，虽然有现成的仪器，但是一组可靠的实验数据需要很好的重复性，因此很费时间。总之，为了研究，时间，金钱，体力，总得支出一项~

tdl522 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2013-6-4 01:43
非常感谢您，我已将其添加到我的文档里了，我想问一下像这个nonassociated species是不是翻译成缔合分子 ...

Renon, H., and J.M. Prausnitz, AIChE J., 14, 135–144 (1968).
http://onlinelibrary.wiley.com/doi/10.1002/aic.690140124/abstract

Renon, H., and J.M. Prausnitz, Ind. Eng. Chem. Process Des. Dev., 8,413–419 (1969).

Mellapak_tju 发表于 2013-5-29 05:55:10

璇雅 发表于 2014-2-18 11:01
您好，我看了其中定义组分和物性方法以后点击运行会出现一组NRTL的交互作用参数，这个与后面提供实验数据后 ...

首先是要通过热力学的一致性校验，可以是点检测也可以是面检测，这个在regression的结果中可以看到，直观的判断方式我觉得应该是对比基于拟合的数据的VLE曲线和你的实验数据，看相关性

Mellapak_tju 发表于 2013-5-29 05:55:10

tdl522 发表于 2013-6-3 07:40
1. aji = 0.20 for hydrocarbons and polar, nonassociated species (e.g., n-heptane/acetone).
2. aji = ...

非常感谢您，我已将其添加到我的文档里了，我想问一下像这个nonassociated species是不是翻译成缔合分子好一点，还有能不能告诉我这些数据的出处，谢啦~

454449711 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2017-9-22 21:59
真是抱歉，我不记得UNIFAC可以&#36 ...

这句话说得什么呢？我这边显示乱码了，难道当时涉及敏感词汇？{:1106_388:}

竹影清风 发表于 2013-5-29 05:55:10

本帖最后由 竹影清风 于 2014-11-7 20:44 编辑

请教楼主你编写的说明中“选择Retrieve TDE binary data（搜索热力学数据引擎的二元数据），会发现Aspen Plus中存在很多数据，我们选择VLE007”，这个是怎么操作的，是在Aspen中存在这些气液平衡数据吗？{:1106_368:}期待回复，谢谢@Mellapak_tju

Mellapak_tju 发表于 2013-5-29 05:55:10

竹影清风 发表于 2014-11-8 19:56
您好，非常感谢您的耐心解答，您真的对Aspen研究的很细致，我研究方向也是ASPEN方向，希望可以多多向您学 ...

这是你aspen本身的问题，实在是帮不了你了，你还是自己在网上找找关于气液平衡数据的文献吧

Mellapak_tju 发表于 2013-5-29 05:55:10

zy2013207188 发表于 2014-12-29 11:01
bkp文件下载下来 打开发现，properties-parameters-binary interaction-NRTL-1页面的NRTL参数的TLOWER和TUP ...

好的，感谢你的建议，话说我还真从没有注意过这些问题，我回头好好再看看，谢谢~

Mellapak_tju 发表于 2013-5-29 05:55:10

imeagor 发表于 2013-6-4 08:53
支持并学习，希望师弟在他乡学成归来，为你高兴！
11年的时候天大开会见过曾老师，祝愿他健康开心

谢谢师兄，才来德国不久，虽然觉得在这里有很多新知识，但大多都很麻烦，我觉得我的硕士都很难毕业……

Mellapak_tju 发表于 2013-5-29 05:55:10

矿大 发表于 2017-7-3 15:57
楼主，做过三元液液相平衡回归二元交互参数的么，或者我怎么能联系上您呢，请教几下问题

我自己还没做过三元LLE的参数回归，但我觉得原理类似，我们可以讨论。我微信machenhaoxa

璇雅 发表于 2013-5-29 05:55:10

本帖最后由 璇雅 于 2014-2-18 11:04 编辑

您好，我看了其中定义组分和物性方法以后点击运行会出现一组NRTL的交互作用参数，这个与后面提供实验数据后模拟出来的结果哪个更准确啊，如何去判断，谢谢，请不吝赐教@Mellapak_tju

Mellapak_tju 发表于 2013-5-29 05:55:10

sh_ssp 发表于 2013-7-29 15:17
找到了，谢了。好像是要欧元支付80欧元，天价啊
- 本文出自马后炮化工论坛，原文地址：http: ...

有的破解版中是包含NIST数据的，你可以找人帮你查查嘛

竹影清风 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2014-11-9 00:13
实在找不到数据，或是现有的数据很贵，自己实验测试VLE数据也不失为一个很好的办法，当然貌似实验测数据 ...

嗯嗯，还是学到了很多，谢谢啦！

snay 发表于 2013-5-29 05:55:10

zy2013207188 发表于 2014-12-29 11:01
bkp文件下载下来 打开发现，properties-parameters-binary interaction-NRTL-1页面的NRTL参数的TLOWER和TUP ...

楼主厉害哇。

hide32 发表于 2013-5-29 05:55:10

楼主对于二元交互参数的详细介绍，解答了这两天一直困扰我的问题：为什么包宗宏老师的书上写着NRTL方程只有三组参数。找了两天资料，终于明白了。谢谢~！

Mellapak_tju 发表于 2013-5-29 05:55:10

454449711 发表于 2020-7-10 13:13
您好，请问能否加您个联系方式呢？准备开始做硕士论文乐，aspen方面还希望能向您学习学习！

527475666，我的微信

Mellapak_tju 发表于 2013-5-29 05:55:10

yingxiongmolu 发表于 2015-5-21 08:27
@Mellapak_tju这位马兄请问一下推荐本关于Aspen学习方面的书籍！

我觉得鲁本的书就很赞！

Mellapak_tju 发表于 2013-5-29 05:55:10

tdl522 发表于 2013-6-4 08:05
Renon, H., and J.M. Prausnitz, AIChE J., 14, 135–144 (1968).
http://onlinelibrary.wiley.com/doi/ ...

我记下来了，非常感谢~

18562317750 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2013-7-29 16:11
有的破解版中是包含NIST数据的，你可以找人帮你查查嘛

打开没有Retrievebinary

454449711 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2020-7-10 11:05
谢谢！

您好，请问能否加您个联系方式呢？准备开始做硕士论文乐，aspen方面还希望能向您学习学习！{:1106_398:}

sh_ssp 发表于 2013-5-29 05:55:10

找到了，谢了。好像是要欧元支付80欧元，天价啊{:1106_361:}
- 本文出自马后炮化工论坛，原文地址：http://meng.horse/thread-67088-4-1.html

zy2013207188 发表于 2013-5-29 05:55:10

或者data的类型设为tpxy,t都输入为348.15K.其实如何进行数据回归的例子在Aspen用户手册和Aspen帮助里边都有，可以去看一下。

sanernim 发表于 2013-5-29 05:55:10

楼主大人，请问你知道unifac模型二元交互作用参数的回归该怎么弄吗？ 或者怎样修改unifac基团啊？ 这个问题困扰我很久了

Mellapak_tju 发表于 2013-5-29 05:55:10

muyue3 发表于 2016-5-25 09:42
请教一下楼主，最后的那个残差图是什么意思？

比较的是实验数据和预测数据之间的残差呀

Mellapak_tju 发表于 2013-5-29 05:55:10

454449711 发表于 2020-7-6 20:01
看完了作者自己写的回归说明书，实在厉害，要学透aspen看来还有很长一段路要走。

谢谢！

imeagor 发表于 2013-5-29 05:55:10

师弟的深入探讨果然精彩，希望你有更多的原创资料，qq已加你
另外nonassociated species翻译为非缔合物质，没有问题

tdl522 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2013-6-4 17:28
我记下来了，非常感谢~
哈哈，客气，期待与imeagor老师和Mellapak_tju同学交流。

Mellapak_tju 发表于 2013-5-29 05:55:10

18562317750 发表于 2019-2-19 16:50
打开没有Retrievebinary

那你重新下载一个啊

imeagor 发表于 2013-5-29 05:55:10

支持并学习，希望师弟在他乡学成归来，为你高兴！
11年的时候天大开会见过曾老师，祝愿他健康开心

sh_ssp 发表于 2013-5-29 05:55:10

请问高手，“c. 选择Retrieve TDE binary data（搜索热力学数据引擎的二元数据）”从何处找啊？

Mellapak_tju 发表于 2013-5-29 05:55:10

gdszlwmg 发表于 2013-6-3 08:09
看帖看完了至少要顶一下，还可以加入到淘帖哦！

什么是淘贴呢？

矿大 发表于 2013-5-29 05:55:10

楼主，做过三元液液相平衡回归二元交互参数的么，或者我怎么能联系上您呢，请教几下问题

Mellapak_tju 发表于 2013-5-29 05:55:10

有错觉 发表于 2019-2-25 10:38
mark11111111111111111

这么老的帖子还有人看，真是感动

Mellapak_tju 发表于 2013-5-29 05:55:10

leihaos 发表于 2018-4-19 10:04
你好 能在聚合物模拟给我点意见吗

没做过聚合物模拟

yingxiongmolu 发表于 2013-5-29 05:55:10

@Mellapak_tju这位马兄请问一下推荐本关于Aspen学习方面的书籍！

454449711 发表于 2013-5-29 05:55:10

看完了作者自己写的回归说明书，实在厉害，要学透aspen看来还有很长一段路要走。

Mellapak_tju 发表于 2013-5-29 05:55:10

hanhuiziling 发表于 2015-4-18 18:34
楼主可以分享一下亚硝酸甲酯么

没做过呀

Mellapak_tju 发表于 2013-5-29 05:55:10

yingxiongmolu 发表于 2015-5-25 22:19
可是他的书，主要是写控制方面的！

都有啊

454449711 发表于 2013-5-29 05:55:10

Mellapak_tju 发表于 2020-7-27 12:33
527475666，我的微信

已加，请您通过

Mellapak_tju 发表于 2013-5-29 05:55:10

wnw3248 发表于 2020-8-9 07:09
这个真的是不错，下载了！

谢谢！

gdszlwmg 发表于 2013-5-29 05:55:10

看帖看完了至少要顶一下，还可以加入到淘帖哦！

zhdgzhdg 发表于 2013-5-29 05:55:10

支持原创，建议设为经典。

pklwz0608 发表于 2013-5-29 05:55:10

啥也不说了，楼主就是给力！

chenhaijunsjy 发表于 2013-5-29 05:55:10

例子非常好,再有过程说明书!{:1106_362:}

xtu_xz 发表于 2013-5-29 05:55:10

楼主真强大啊，好好学习了{:1106_383:}

zhangyili 发表于 2013-5-29 05:55:10

啥也不说了，楼主就是给力！

罗修绝 发表于 2013-5-29 05:55:10

膜拜神贴，后面的请保持队形~

zy7181 发表于 2013-5-29 05:55:10

看了LZ的帖子，我只想说一句很好很强大！

bigchenrui 发表于 2013-5-29 05:55:10

这是一个好贴，谢谢分享！

kaitian007 发表于 2013-5-29 05:55:10

此乃严谨之学风！赞一个

秋古风 发表于 2013-5-29 05:55:10

不管你信不信，反正我是信了。

sdlg348 发表于 2013-5-29 05:55:10

不错，又占了一个沙发！

appleliangzi 发表于 2013-5-29 05:55:10

谢谢分享，灰常感谢

山中来客 发表于 2013-5-29 05:55:10

lz强大，支持原创！

苦海女神龙 发表于 2013-5-29 05:55:10

咱还是新手

苦海女神龙 发表于 2013-5-29 05:55:10

就需要这个

查看完整版本: Aspen入门实践 实例1 苯-噻吩体系的NRTL模型二元交互作用参数的回归