CN108875293B - 一种基于基团贡献理论预测混合物临界性质的方法 - Google Patents

一种基于基团贡献理论预测混合物临界性质的方法 Download PDF

Info

Publication number
CN108875293B
CN108875293B CN201810564142.9A CN201810564142A CN108875293B CN 108875293 B CN108875293 B CN 108875293B CN 201810564142 A CN201810564142 A CN 201810564142A CN 108875293 B CN108875293 B CN 108875293B
Authority
CN
China
Prior art keywords
critical
group
mixture
pure
contribution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201810564142.9A
Other languages
English (en)
Other versions
CN108875293A (zh
Inventor
刘向阳
王成杰
张颖
何茂刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201810564142.9A priority Critical patent/CN108875293B/zh
Publication of CN108875293A publication Critical patent/CN108875293A/zh
Application granted granted Critical
Publication of CN108875293B publication Critical patent/CN108875293B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种基于基团贡献理论预测混合物临界性质的方法,通过纯质分子结构划分出第一水平和第二水平基团库;利用纯质临界性质的实验数据拟合得到各个基团对临界性质的基团贡献值;提出了混合物临界性质和基团贡献值之间的关联模型,并搜集二元混合物的临界性质实验数据,利用各类体系临界性质的实验数据拟合得到对应的混合物模型参数;获得基团贡献值和混合物模型参数后便可预测多元混合物的临界性质,本发明具有不需要纯质临界体积数据作为基础且可以在纯质临界性质未知的情况下精确预测混合物临界性质的优点,在预测临界温度和临界压力时仅分别用到临界温度和临界压力的基团贡献值,并不需要临界体积的基团贡献值,能够区分同分异构体。

Description

一种基于基团贡献理论预测混合物临界性质的方法
技术领域
本发明属于流体热物性理论计算领域,涉及一种基于基团贡献理论预测混合物临界性质的方法。
背景技术
临界性质(Tc、pc、Vc)是流体最重要的热物理性质之一,它不仅是流体状态方程的关键参数,也常用于pVTx和黏度等热力学性质和迁移性质的计算。另外,临界性质还是超临界萃取等过程设计时的基础数据。然而,由于临界性质的实验测定较为困难,混合物临界性质的实验数据极为匮乏,难以满足化工生产和研究的需求,因此,准确的混合物临界性质理论估算方法显得尤为重要。
目前为止,已经有很多理论估算方法被提出,主要可以分为两种。一种是基于Gibbs混合物临界态严格热力学判据的解析方法,该方法尤其复杂且计算精度受混合法则的影响较大。另一种是经验估算方法,该方法因为计算简单而被广泛应用。例如,P.L.Chueh,J.M.Prausnitz,Vapor–liquid equilibria at high pressures:calculationof critical temperatures,volumes,and pressures of nonpolar mixtures,AIChEJ.13(1967)1107–1113中提出了著名的Chueh–Prausnitz方程并获取了87种混合物的二元相互作用参数;A.Kreglewski,W.B.Kay,Critical constants of conformed mixtures,J.Phys.Chem.73(1969)3359–3366中运用保形溶液理论来预测混合物的临界温度和临界压力;C.F.Spencer,T.E.Daubert,R.P.Danner,A critical review of correlations forthe critical properties of defined mixtures,AIChE J.19(2010)522–527中用于预测烃类混合物临界性质的估算方法,发现Kreglewski方法在预测临界压力方面更为准确。然而,运用这些经验方法预测混合物的临界性质时一般都需要用到纯质的临界数据。
L.Li,E.Kiran,Estimation of critical properties of binary mixturesusing group contribution methods,Chem.Eng.Commun.94(1990)131–141中将A.L.Lydersen,Estimation of Critical Properties of Organic Compounds,University of Wisconsin,Coll.Eng.,Engineering Experimental Station Report 3,Madison,Wisconsin,USA(1955)提出的基团贡献法和混合法则结合起来预测了二元混合物的临界性质。在该方法中,利用混合法则将混合物的基团贡献值表示为各组分基团贡献值的函数,再将其代入到Lydersen提出的方程中求得混合物的临界性质。该方法虽然没有用到纯质临界数据,但是在预测混合物的临界温度和临界压力时需要临界体积的基团贡献值,而临界体积的基团贡献值要由临界体积的实验数据拟合得到。由于临界体积实验数据非常缺乏,该方法受到了极大的限制。另外,Lydersen基团贡献法不能区分同分异构体,使得Li方法在估算含有同分异构体的混合物临界性质时出现较大偏差,不过目前已有一些基团贡献法可以通过划分第二水平基团而有效区分同分异构体,如Y.H Liang,P.S.Ma,A newgroup-contribution method for critical properties,Chin.J.Chem.Eng.8(2000)74–79基团贡献法。
综上所述,需要提出一种无需纯质临界数据和临界体积基团贡献值,且能区分同分异构体的混合物临界温度和临界压力预测方法。
发明内容
本发明的目的在于提供一种基于基团贡献理论预测混合物临界性质的方法,以克服现有技术的不足。
为达到上述目的,本发明采用如下技术方案:
一种基于基团贡献理论预测混合物临界性质的方法,包括以下步骤:
将纯质分子结构划分为第一水平基团库和第二水平基团库,分别利用纯质临界温度和纯质临界压力的实验数据拟合得到第一水平基团和第二水平基团对临界温度和临界压力的基团贡献值;构建关联混合物临界性质和基团贡献值之间的关联模型,利用二元混合物临界温度和临界压力的实验数据拟合得到混合物模型参数,获得基团贡献值和混合物模型参数后在未知纯质临界性质的条件下对混合物的临界性质进行预测。
进一步的,其中第一水平基团为简单官能团,第二水平基团为以第一水平基团为连接点的官能团。
进一步的,分别利用纯质临界温度和临界压力拟合得到第一水平基团和第二水平基团对临界温度和临界压力的基团贡献值:
Figure BDA0001684069960000031
Figure BDA0001684069960000032
Figure BDA0001684069960000033
Figure BDA0001684069960000034
Figure BDA0001684069960000035
Figure BDA0001684069960000036
式中:Tc为纯质的临界温度,pc为纯质的临界压力;Tb为纯质的沸点,a1–a8为纯质模型参数,ni为第一水平基团的个数,mj为第二水平基团个数;ΔTci和Δpci为第一水平基团i对临界温度和临界压力的贡献值,ΔTcj和Δpcj为第二水平基团j对临界温度和临界压力的贡献值;Tg和pg为第一水平基团对临界温度和临界压力总的基团贡献值;Ts和ps为第二水平基团对临界温度和临界压力总的基团贡献值;w为是否使用第二水平基团的判断参数。
进一步的,当w=0时,仅使用第一水平基团;当w=1时,同时使用第一水平基团和第二水平基团。
进一步的,二元混合物分为正烷烃、其他碳氢化合物、含氧化合物、含卤代物和含硅基化合物五类二元混合体系。
进一步的,若二元混合物中两组分属于不同类化合物,则其分类以排序靠后化合物的分类为准。
进一步的,模型参数有五组,分别对应于五类二元混合体系,并用于预测不同混合体系的临界温度和临界压力。
进一步的,混合物临界性质和基团贡献值之间的关联模型为:
Figure BDA0001684069960000041
Figure BDA0001684069960000042
Figure BDA0001684069960000043
Figure BDA0001684069960000044
式中:Tcm为混合物的临界温度,pcm为混合物的临界压力,c1–c8为混合物模型参数,xi(i=1,2,…,n)为第i种组分的摩尔分数,Tc,i和pc,i为第i种组分的临界温度和临界压力,Tg,i和pg,i为第i种组分第一水平基团对临界温度和临界压力总的基团贡献值,i<j。
与现有技术相比,本发明具有以下有益的技术效果:
本发明一种基于基团贡献理论预测混合物临界性质的方法,通过纯质分子结构划分出第一水平基团库和第二水平基团库;利用纯质临界性质的实验数据拟合得到各个基团对临界性质的基团贡献值;提出了混合物临界性质和基团贡献值之间的关联模型,并搜集二元混合物的临界性质实验数据,利用各类体系临界性质的实验数据拟合得到对应的模型参数;获得基团贡献值和混合物模型参数后便可预测多元混合物的临界性质,本发明同现有方法相比,具有不需要纯质临界体积数据作为基础且可以在纯质临界性质未知的情况下精确预测混合物临界性质的优点,在预测临界温度和临界压力时仅分别用到临界温度和临界压力的基团贡献值,并不需要临界体积的基团贡献值,而且能够区分同分异构体。
进一步的,将混合物分为正烷烃、其他碳氢化合物、含氧化合物、含卤代物和含硅基化合物五类二元混合体系,将二元混合物分为五类二元混合体系,分别利用大量的临界性质实验数据关联得到各体系对应的模型参数,具有较高的预测精度。
附图说明
图1是二元混合物临界温度拟合效果图;
图2是二元混合物临界压力拟合效果图;
图3是二元和三元混合物临界温度预测结果图。
图4是二元和三元混合物临界压力预测结果图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
一种基于基团贡献理论预测混合物临界性质的方法,包括以下步骤:
将纯质分子结构划分出第一水平基团库和第二水平基团库,分别利用纯质临界温度和纯质临界压力的实验数据拟合得到第一水平基团和第二水平基团对临界温度和临界压力的基团贡献值;构建关联混合物临界性质和基团贡献值之间的关联模型,利用二元混合物临界温度和临界压力的实验数据拟合得到混合物模型参数,获得基团贡献值和混合物模型参数后在未知纯质临界性质的条件下对混合物的临界性质进行预测。
基团划分与基团贡献值获取:
选取N1种纯质,划分出第一水平基团和第二水平基团,并利用纯质临界温度的实验数据和关联式(1–3)分别拟合得到第一水平基团和第二水平基团对临界温度的基团贡献值,同理,选取N2种纯质,划分出第一水平基团和第二水平基团,并利用纯质临界压力的实验数据和关联式(4–6)分别拟合得到第一水平基团和第二水平基团对临界压力的基团贡献值;
第一水平基团为简单官能团,第二水平基团为以第一水平基团为连接点的官能团;
所述的第一水平基团和第二水平基团对临界温度和临界压力的基团贡献值由以下公式获得:
Figure BDA0001684069960000061
Figure BDA0001684069960000062
Figure BDA0001684069960000063
Figure BDA0001684069960000064
Figure BDA0001684069960000065
Figure BDA0001684069960000071
式中:Tc为纯质的临界温度,pc为纯质的临界压力;Tb为纯质的沸点,a1–a8为纯质模型参数,ni为第一水平基团的个数,mj为第二水平基团个数;ΔTci和Δpci为第一水平基团i对临界温度和临界压力的贡献值,ΔTcj和Δpcj为第二水平基团j对临界温度和临界压力的贡献值;Tg和pg为第一水平基团对临界温度和临界压力总的基团贡献值;Ts和ps为第二水平基团对临界温度和临界压力总的基团贡献值。w为是否使用第二水平基团的一个判断参数,当w=0时,仅使用第一水平基团;当w=1时,同时使用第一水平基团和第二水平基团;
二元混合物分为正烷烃、其他碳氢化合物、含氧化合物、含卤代物和含硅基化合物五类二元混合体系,若二元混合物中两组分属于上述不同类化合物,则其分类以排序靠后化合物的分类为准;如,碳氢化合物和含卤代物的混合物属于含卤代物二元混合体系。
混合模型参数获取:
获得基团贡献值和纯质模型参数后能够在未知纯质临界性质的条件下对混合物的临界性质进行预测,混合物临界性质和基团贡献值之间的关联模型为:
Figure BDA0001684069960000072
Figure BDA0001684069960000073
Figure BDA0001684069960000074
Figure BDA0001684069960000081
式中:Tcm为混合物的临界温度,pcm为混合物的临界压力,c1–c8为混合物模型参数,xi(i=1,2,…,n)为第i种组分的摩尔分数,Tc,i和pc,i为第i种组分的临界温度和临界压力,Tg,i和pg,i为第i种组分第一水平基团对临界温度和临界压力总的基团贡献值,i<j。模型参数有五组,分别对应于五类二元混合体系,并用于预测不同混合体系的临界温度和临界压力。
将二元混合物分为正烷烃、其他碳氢化合物、含氧化合物、含卤代物和含硅基化合物五类二元混合体系,利用混合体系的临界性质实验数据和关联式(7–10)拟合得到各体系临界温度和临界压力的模型参数。图1是379种二元混合物3056个临界温度的拟合效果图,拟合的相对偏差最大不超过±11%。图2是188种二元混合物1602个临界压力的拟合效果图,拟合的相对偏差均不超过±21.5%。
选取相应的基团贡献值和模型参数代入关联式(7–10)中并预测了42种二元混合物的281个临界温度和8种二元混合物的58个临界压力,并将预测结果和搜集到的临界性质实验数据进行对比。所有的临界性质实验数据均没有参与基团划分与基团贡献值获取和混合模型参数获取中的拟合计算。42种二元混合物包括6种正烷烃,4种其他碳氢化合物,7种含氮氧化合物和25种含卤代物二元混合体系。8种二元混合物包括2种其他碳氢化合物,4种含氮氧化合物和2种含卤代物二元混合体系。每类二元混合体系临界温度的平均相对绝对预测偏差分别为2.83%,5.76%,1.51%和1.92%,临界压力的平均相对绝对预测偏差分别为5.45%,5.86%和9.49%。
另外,本发明所提出的方法可以拓展到三元混合物临界性质的预测。利用本发明所提出的方法预测了12种三元混合物的289个临界温度和临界压力,并将预测结果与搜集的实验数据进行对比。三元混合体系的分类类似于二元混合体系,12种三元混合物包括2种正烷烃和10种含氮氧化合物三元混合体系。两类三元混合体系临界温度的平均相对绝对预测偏差分别为4.31%和1.29%,临界压力的平均相对绝对预测偏差分别为19.07%和3.59%。图3和图4是二元和三元混合物临界温度和临界压力的预测结果图,可以看出临界温度的相对偏差绝大部分在±15%以内,而临界压力的相对偏差绝大部分在±30%以内。

Claims (7)

1.一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,包括以下步骤:
将纯质分子结构划分为第一水平基团库和第二水平基团库,分别利用纯质临界温度和纯质临界压力数据拟合得到第一水平基团和第二水平基团对临界温度和临界压力的基团贡献值;构建关联混合物临界性质和基团贡献值之间的关联模型,利用二元混合物临界温度和临界压力的实验数据拟合得到混合物模型参数,获得基团贡献值和混合物模型参数后在未知纯质临界性质的条件下对混合物的临界性质进行预测;分别利用纯质临界温度和纯质临界压力拟合得到第一水平基团和第二水平基团对临界温度和临界压力的基团贡献值:
Figure FDA0002456113080000011
Figure FDA0002456113080000012
Figure FDA0002456113080000013
Figure FDA0002456113080000014
Figure FDA0002456113080000015
Figure FDA0002456113080000016
式中:Tc为纯质的临界温度,pc为纯质的临界压力;Tb为纯质的沸点,a1–a8为纯质模型参数,ni为第一水平基团的个数,mj为第二水平基团个数;△Tci和△pci为第一水平基团i对临界温度和临界压力的贡献值,△Tcj和△pcj为第二水平基团j对临界温度和临界压力的贡献值;Tg和pg为第一水平基团对临界温度和临界压力总的基团贡献值;Ts和ps为第二水平基团对临界温度和临界压力总的基团贡献值;w为是否使用第二水平基团的判断参数。
2.根据权利要求1所述的一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,其中第一水平基团为简单官能团,第二水平基团为以第一水平基团为连接点的官能团。
3.根据权利要求2所述的一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,当w=0时,仅使用第一水平基团;当w=1时,同时使用第一水平基团和第二水平基团。
4.根据权利要求1所述的一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,二元混合物分为正烷烃、其他碳氢化合物、含氧化合物、含卤代物和含硅基化合物五类二元混合体系。
5.根据权利要求4所述的一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,若二元混合物中两组分属于不同类化合物,则其分类以排序靠后化合物的分类为准。
6.根据权利要求4所述的一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,模型参数有五组,分别对应于五类二元混合体系,并用于预测不同混合体系的临界温度和临界压力。
7.根据权利要求1所述的一种基于基团贡献理论预测混合物临界性质的方法,其特征在于,获得基团贡献值和混合物模型参数后在未知纯质临界性质的条件下对混合物的临界性质进行预测,混合物临界性质和基团贡献值之间的关联模型为:
Figure FDA0002456113080000021
Figure FDA0002456113080000022
Figure FDA0002456113080000031
Figure FDA0002456113080000032
式中:Tcm为混合物的临界温度,pcm为混合物的临界压力,c1–c8为混合物模型参数,xi(i=1,2,…,n)为第i种组分的摩尔分数,Tc,i和pc,i为第i种组分的临界温度和临界压力,Tg,i和pg,i为第i种组分第一水平基团对临界温度和临界压力总的基团贡献值,i<j。
CN201810564142.9A 2018-06-04 2018-06-04 一种基于基团贡献理论预测混合物临界性质的方法 Expired - Fee Related CN108875293B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810564142.9A CN108875293B (zh) 2018-06-04 2018-06-04 一种基于基团贡献理论预测混合物临界性质的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810564142.9A CN108875293B (zh) 2018-06-04 2018-06-04 一种基于基团贡献理论预测混合物临界性质的方法

Publications (2)

Publication Number Publication Date
CN108875293A CN108875293A (zh) 2018-11-23
CN108875293B true CN108875293B (zh) 2020-07-28

Family

ID=64335294

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810564142.9A Expired - Fee Related CN108875293B (zh) 2018-06-04 2018-06-04 一种基于基团贡献理论预测混合物临界性质的方法

Country Status (1)

Country Link
CN (1) CN108875293B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111892938B (zh) * 2020-06-12 2022-08-30 中国石油天然气股份有限公司 一种原油分子组成的确定方法、系统、装置及存储介质
CN111899814B (zh) * 2020-06-12 2024-05-28 中国石油天然气股份有限公司 一种单分子、混合物物性计算方法、设备及存储介质
CN112382348B (zh) * 2020-11-27 2022-03-29 华南理工大学 一种基于枚举法设计和筛选结晶过程溶剂的方法
CN113380347B (zh) * 2021-06-30 2024-07-16 北京化工大学 复合固体推进剂体系相关化学反应热力学性质的分析方法
CN115482885A (zh) * 2022-09-30 2022-12-16 广东电网有限责任公司 一种二元汽液相平衡预测方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103217513A (zh) * 2012-01-18 2013-07-24 中国石油天然气股份有限公司 一种基于调合指数概念的轻柴油冷滤点的预测方法
CN106096238A (zh) * 2016-06-01 2016-11-09 东南大学 一种预测生物质热解动力学参数的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9026421B2 (en) * 2011-08-03 2015-05-05 Purdue Research Foundation Method of modeling cloud point of a mixture of fatty acid methyl esters using a modified UNIFAC model and a system therefor
WO2013047494A1 (ja) * 2011-09-29 2013-04-04 公立大学法人首都大学東京 相平衡データの高精度相関方法及び相平衡推算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103217513A (zh) * 2012-01-18 2013-07-24 中国石油天然气股份有限公司 一种基于调合指数概念的轻柴油冷滤点的预测方法
CN106096238A (zh) * 2016-06-01 2016-11-09 东南大学 一种预测生物质热解动力学参数的方法

Also Published As

Publication number Publication date
CN108875293A (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
CN108875293B (zh) 一种基于基团贡献理论预测混合物临界性质的方法
Dijkstra et al. Entropic wetting and many-body induced layering in a model colloid-polymer mixture
Ghasemitabar et al. Estimation of the normal boiling point of organic compounds via a new group contribution method
Zamani et al. Effect of porous media properties on the onset of polymer extensional viscosity
Gani et al. Molecular structure based estimation of properties for process design
Papaioannou et al. Dynamic viscosity of multicomponent liquid mixtures
Mansour et al. Predicting PVT properties of Egyptian crude oils by a modified Soave–Redlich–Kowng equation of state
He et al. Prediction of critical temperature and critical pressure of multi-component mixtures
Delhommelle Should “lane formation” occur systematically in driven liquids and colloids?
Tsai et al. A model for predicting the auto-ignition temperature using quantitative structure property relationship approach
Perry et al. Estimation of the change point of the process fraction nonconforming in SPC applications
Panah Modeling binary vapor–liquid equilibrium data containing perfluorocarbons using the Peng–Robinson and the PC-SAFT equations of state
Shene et al. Mathematical modeling of elution curves for a protein mixture in ion exchange chromatography and for the optimal selection of operational conditions
Jin et al. Prediction of the normal boiling point of oxygen containing organic compounds using quantitative structure–property relationship strategy
Khosronezhad Gheshlaghi et al. Prediction of kinematic viscosity of petroleum fractions using artificial neural networks
Toxvaerd et al. Phase-stability lines, spinodals and distribution functions for binary immiscible mixtures
Tamaki et al. Multiparametric virtual metrology model building by job-shop data fusion using a Markov chain Monte Carlo method
Fotouh et al. A comparative study of numerical methods for calculating phase equilibrai in fluid mixtures from an equation of state
Kholmatov et al. Monte Carlo simulations of phase equilibria of binary mixtures containing methane, ethane, n-pentane, nitrogen, and oxygen comparison with experimental measurements and predictions of Peng-Robinson EOS and Perturbed-Chain SAFT
Issaoui et al. Correlation of the binary interaction factor for polar solutes dissolved in supercritical carbon dioxide
KR20140145753A (ko) 기준 분자들의 물성 값 및 분자 표현자에 기반한 분자간 근원적 상관성을 이용하여 혼합물을 구성하는 각 분자들의 물성 값을 예측하는 방법
Lindegren et al. Bottom water formation in the Weddell Sea resolved by principal component analysis and target estimation
Vasquez et al. Effect of data type on thermodynamic model parameter estimation: A Monte Carlo approach
Tham et al. Densities of potassium hydroxide solutions
Kumar et al. Prediction of Liquid Mixture Viscosity for Dry and Water-Saturated TBP/[italic] n [roman]-Dodecane Solutions

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200728