CN109839495B - 一种根据分子组成计算轻质石油馏分馏程的方法 - Google Patents
一种根据分子组成计算轻质石油馏分馏程的方法 Download PDFInfo
- Publication number
- CN109839495B CN109839495B CN201910221942.5A CN201910221942A CN109839495B CN 109839495 B CN109839495 B CN 109839495B CN 201910221942 A CN201910221942 A CN 201910221942A CN 109839495 B CN109839495 B CN 109839495B
- Authority
- CN
- China
- Prior art keywords
- distillation
- mixed
- molecular
- gas
- substance
- 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.)
- Active
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明涉及一种根据分子组成计算轻质石油馏分馏程的方法，包括如下步骤：(1)汇合石油生产中所有轻质馏分的分子种类，建立涵盖轻质石油馏分和产品的分子组成数据库；(2)基于分子组成数据库建立分子性质及温度参数数据库：(3)基于分子组成数据库建立分子混合性质计算模型数据库；(4)获取被检测油品真实的单体烃和含氧化合物组成数据，该数据可借助实验室分析获取，或者通过数学模型预测获取；(5)基于步骤(4)获得的单体烃和含氧化合物，计算得到轻质石油馏分的馏程。本发明的通用性好，可以适用于轻质石油馏分和产品的馏程计算，并且准确度高；解决了现有色谱模拟蒸馏与恩氏蒸馏关联转换误差大，难度大的问题。
Description
技术领域
本发明涉及石油炼制及石油化工生产技术领域，尤其涉及一种根据分子组成计算轻质石油馏分馏程的方法。
背景技术
轻质石油馏分和产品的蒸馏特性，是表征燃料或石油产品挥发性能的重要特性。它提供了燃料或者石油产品的组成、性能、使用性能等信息。蒸馏特性对于车用汽油极为重要，它会影响发动机的启动、升温性能、产生气阻的倾向。在这些燃料中存在高沸点的组分可显著影响固体燃烧沉积物的生成程度。因此，在石油产品规格、商业合同协议以及炼化企业生产控制中均会对蒸馏特性的限值有所规定。国际通用检测汽油馏程温度的标准是采用美国试验与材料协会标准ASTM D86《石油产品常压蒸馏试验法》。我国在此标准基础上重新起草了《石油产品常压蒸馏特性测定法》GB/T6536。以上标准规定采用实验室间歇蒸馏仪器定量测定。但是采用间歇蒸馏仪器测定馏程特性，不仅需要大量的样品，需要消耗100～200ml试样，而且费时，通常要进行2～4小时时间。
在气相色谱技术基础上，发展了色谱模拟蒸馏分析方法。色谱模拟蒸馏具有快速、简便、重复性高等优势。在色谱条件下，将试样按组分沸点次序分离，同时进行切片积分，获得对应的累加面积，经过温度-时间的内插校正，就可以得到对用百分收率的温度。色谱模拟蒸馏结果接近实沸点蒸馏(ASTM D2892)。ASTM D2887规定了终馏点不超过260℃的石油产品和馏分的色谱模拟馏程测定方法。但是该方法不适用于包含含氧化合物的石油产品和馏分。
ASTM D7096规定了采用大孔径毛细管气相色谱测定终馏点不超过280℃的石油产品和馏分的色谱模拟馏程测定方法。可以测量乙醇汽油的馏程，其余含氧汽油也不排除，但是预测误差较大。该方法还可以获取汽油中C5以内的饱和烃浓度信息。从色谱模拟蒸馏转换为D86馏程温度采用多元线性回归的办法建立关联式。API technical data book中总结了ASTM D86与色谱模拟蒸馏数据的转换方法。建立这种关联关系，首先需要进行大量的实验。其次，不同类型的馏分和石油产品，关联方程不能通用。或者误差很大。第三，对于不同含量乙醇汽油等需要重新建立关联方程。第四，只能测算标准规定中的初馏点、10％蒸馏温度、30％蒸馏温度、50％蒸馏温度、70％蒸馏温度、90％蒸馏温度、终馏点。不能计算任意馏出体积的对应温度。第五、误差很大。关联式受到环境温度、压力的影响，随着时间的推移，可能越来越不准。或者组分改变，模型也不通用。需要重新建模。
随着仪器分析技术的进展，目前已经能分析出C1-C13的单体烃浓度信息。基于单体烃浓度信息，可以建立间歇蒸馏模型预测D86恩氏馏程特性。荷兰AC分析公司Spieksma(US Patent US6,711,532B1)，在气相色谱检测获得详细烃数据的基础上，模拟了简单蒸馏装置。该方法需要利用保留因子计算各物质的物性，但由于现有技术不能分析含氧的组分的详细单体组成，因此对含氧馏分，如乙醇汽油还无法预测。Greenfield等将蒸馏过程简化为二级非稳态精馏过程，建立质量守恒和能量守恒模型，并考虑了初始时刻蒸馏烧瓶内的空气对温度测量的影响。
现有色谱模拟蒸馏与恩氏蒸馏关联转换误差大，难度大。色谱模拟蒸馏获得的结果接近是实沸点蒸馏，需要采用关联方法转换为D86恩氏蒸馏。关联方法多是回归分析。对于一些组分复杂的石油产品和馏分预测误差较大。现有蒸馏模型计算不同的石油馏分通用性不足。石油产品蒸馏是一个复杂的非稳态传递过程，具体涉及到动量传递、热量传递和质量传递。现有的基于单体烃计算的蒸馏模型，做了过多简化，导致通用性不足。例如将蒸馏过程简化为气液平衡过程，但事实上，由于蒸馏过程的非稳态特性和传递效率不同，一些复杂组分的蒸馏过程是偏离平衡状态的。汽油中添加含氧化合物后性质改变，因此按照烃类组成建立的模型，不能同时预测含氧汽油的馏程温度。现有基于单体烃预测蒸馏曲线的模型误差大。第一，现有模型无法衡量动态滞留量的影响。动态滞留量是指在蒸馏过程中出现在蒸馏烧瓶的瓶颈、支管和冷凝管中的物料。现有模型要么假设蒸气冷凝非常快速，要么假设这部分滞留量是恒定的。从而导致了预测偏差。第二，现有蒸馏模型对轻组分损失、初馏点以及终馏点判断不准。初馏点的判断与动态滞留量以及轻组分损失有关。现有模型无法精确描述动态滞留量以及轻组分损失，导致初馏点判断不准。终馏点的情况更加复杂，不同的组分在蒸馏到终馏点附近时表现不同，现有模型则简单的设置一个残留量，以残留量的沸点做为终馏点。第三，现有模型在计算过程中的物性参数时，假设为固定参数。实际上，这些物性均是温度和组成的函数，不能简单按照常数对待。第四，含氧汽油如乙醇汽油、甲醇汽油的馏程温度预测不准。
发明内容
本发明为克服上述的不足之处，目的在于提供一种根据分子组成计算轻质石油馏分馏程的方法，本发明的通用性好，可以适用于轻质石油馏分和产品的馏程计算，并且准确度高；解决了现有色谱模拟蒸馏与恩氏蒸馏关联转换误差大，难度大的问题。
本发明是通过以下技术方案达到上述目的：一种根据分子组成计算轻质石油馏分馏程的方法，包括如下步骤：
(1)汇合石油生产中所有轻质馏分的分子种类，建立涵盖轻质石油馏分和产品的分子组成数据库；
(2)基于分子组成数据库建立分子性质及温度参数数据库：
(3)基于分子组成数据库建立分子混合性质计算模型数据库；
(4)获取被检测油品真实的单体烃和含氧化合物组成数据，该数据可借助实验室分析获取，或者通过数学模型预测获取；
(5)基于步骤(4)获得的单体烃和含氧化合物，计算得到轻质石油馏分的馏程。
作为优选，所述的分子性质及温度参数数据库包括16项性质：分子气体密度、液体密度、摩尔质量、摩尔体积、表面张力、运动粘度、气体比热、液体比热、气体导热系数、液体导热系数、蒸发热、临界温度、临界压力、偏心因子、气相PVT状态方程参数、活度系数模型参数；其中，上述性质采用模型计算或查阅文献资料两种方法结合获得。
作为优选，所述的分子混合性质计算模型数据库包括混合气相密度、混合液相密度、混合气摩尔质量、混合液相摩尔体积、混合液相表面张力、混合气相运动粘度、混合液相运动粘度、混合气比热、混合液比热、混合气导热系数、混合液导热系数、混合液蒸发热、混合气PVT状态方程参数、混合液活度系数。
作为优选，所述的混合气相密度、混合液相密度、混合气摩尔质量、混合液相摩尔体积、混合气比热、混合液比热、混合气导热系数、混合液导热系数采用以下公式：
其中，xi表示混合体系中物质i的摩尔分数，fi表示混合体系中物质i的纯物质性质。
作为优选，所述的混合液相表面张力采用以下公式：
其中，[Pi]为混合体系中物质i等张比容，ρL，i为混合体系中物质i液相密度，ρV，i为混合体系中物质i气相密度。
作为优选，所述的混合气相运动粘度、混合液相运动粘度采用以下公式：
其中，xi表示混合体系中物质i的摩尔分数，ηi表示混合体系中物质i的纯物质性质。
作为优选，所述的混合液蒸发热采用以下公式：
其中，γi为混合液相体系中物质i的活度系数。
作为优选，所述的混合气PVT状态方程参数a,b采用以下公式计算：
其中，yi表示混合气中物质i的摩尔分数，ai，bi表示混合体系中物质i的纯物质PVT状态方程参数。
作为优选，所述步骤(5)具体如下：
(5.1)根据真实的馏程分析设备，建立严格的蒸馏模型——液膜表面自动更新蒸馏模型，得到模型的方程组；
(5.2)建立一种求解液膜表面自动更新蒸馏模型的数值算法——参数耦合迭代有限差分法；在求解模型过程中，采用内外循环嵌套程序，外循环求解在特定时间步长内的蒸发量以及蒸发气的组成和温度；内循环则联立求解步骤(5.1)中的方程组，方程组初始参数根据外循环获得；根据求解的结果，不断迭代更新方程组参数，直到步骤(5.1)中的方程收敛。
作为优选，所述步骤(5.1)包括如下方程组：
(i)描述蒸馏瓶球部液相质量守恒和能量守恒微分方程组以及初始条件；
(ii)描述蒸馏中所产生气相沿着蒸馏瓶轴向运动的一维动态蒸馏动量传递、质量传递、热量传递方程以及定解条件；
(iii)描述蒸馏中气相沿着瓶壁冷凝的一维动态蒸馏动量传递、质量传递、热量传递方程以及定解条件；
(iv)描述蒸馏中气相沿着瓶壁冷凝的质量守恒、能量守恒方程；
(v)描述蒸馏过程中的气液平衡方程。
本发明的有益效果在于：(1)本发明的通用性好，可以适用轻质石油馏分和产品的馏程计算，包括溶剂油、汽油调和组分和成品汽油、乙醇汽油、甲醇汽油；(2)本发明准确度高，预测轻质石油馏分和产品馏程——包括常用的初馏点、5％馏出温度、10％馏出温度、30％馏出温度、50％馏出温度、70％馏出温度、90％馏出温度、95％馏出温度、终馏点，以及其他任意馏出体积所对应的温度，或任意温度所对应的馏出体积。
附图说明
图1是本发明的流程示意图；
图2是本发明的参数耦合迭代有限差分法流程示意图；
图3是本发明实施例的二甲苯溶剂馏程预测示意图；
图4是本发明实施例的95#汽油馏程预测示意图；
图5是本发明实施例的E10乙醇汽油馏程预测示意图；
图6是本发明实施例的M15甲醇汽油馏程预测示意图。
具体实施方式
下面结合具体实施例对本发明进行进一步描述，但本发明的保护范围并不仅限于此：
实施例：以下通过一个预测乙醇汽油的具体实施方式对本发明作进一步描述；如图1所示，一种根据分子组成计算轻质石油馏分馏程的方法包括如下步骤：
(1)汇总石油生产中所有轻质馏分的分子种类，包括溶剂油、汽油调和组分和成品汽油、乙醇汽油、甲醇汽油，形成一个涵盖碳数C3-C13烃类化合物和含氧化合物的总分子库。总分子库包含超过300个分子，其中烃类化合物包括：
异构烷烃化合物，包括C3-C13单取代、多取代异构烷烃至少100个分子。特征分子包括异丁烷、2,2-二甲基丙烷、2,2,3-三甲基丁烷、2,2,4-三甲基戊烷、2,6-二甲基辛烷、C13异构烷烃等。
正构烷烃化合物，包括C3-C13正构烷烃11个分子。
环烷烃化合物，包括环戊烷、环己烷、1,2-二甲基环戊烷、1,4-二甲基环己烷等至少60个分子。
烯烃化合物，包括3-甲基-1-丁烯、4-甲基-1-戊烯、4,4-二甲基-反-2-戊烯等共等至少60个分子。
芳烃化合物，包括苯、甲苯、对二甲苯、邻二甲苯、间二甲苯、乙苯、丙苯、丁基苯等至少70个分子。
含氧化合物包括：醚类化合物(含甲基叔丁醚、乙基叔丁基醚、叔戊基甲基醚)、醇类化合物(乙醇、甲醇)。
(2)针对步骤(1)所建分子组成数据库，建立分子性质及温度参数数据库：包括16项性质(详见表1)：分子气体密度、液体密度、摩尔质量、摩尔体积、表面张力、运动粘度、气体比热、液体比热、气体导热系数、液体导热系数、蒸发热、临界温度、临界压力、偏心因子、气相PVT状态方程参数、活度系数模型参数。上述性质采用模型计算或可查阅文献资料两种方法结合获得。
| 气相密度，g/cm3 | 气体导热系数，W/(cm.K) |
| 液相密度，g/cm3 | 液体导热系数，W/(cm.K) |
| 摩尔质量，g/mol | 蒸发热，J/mol |
| 摩尔体积，cm3/mol | 临界压力，kpa |
| 表面张力，N/m | 临界温度，kpa |
| 运动粘度，g/(cm.s) | 偏心因子 |
| 气体比热，J/(mol.K) | PVT状态方程参数(PR方程) |
| 液体比热，J/(mol.K) | 活度系数模型参数(UNIFAC) |
表1
(3)针对步骤(1)所建分子组成数据库，建立分子混合性质计算模型数据库，包括混合气相密度、混合液相密度、混合气摩尔质量、混合液相摩尔体积、混合液相表面张力、混合气相运动粘度、混合液相运动粘度、混合气比热、混合液比热、混合气导热系数、混合液导热系数、混合液蒸发热、混合气PVT状态方程参数、混合液活度系数。其中，混合气相密度、混合液相密度、混合气摩尔质量、混合液相摩尔体积、混合气比热、混合液比热、混合气导热系数、混合液导热系数采用以下公式：
其中xi表示混合体系中物质i的摩尔分数，fi表示混合体系中物质i的纯物质性质。
混合液相表面张力采用以下公式：
其中[Pi]为混合体系中物质i等张比容，ρL，i为混合体系中物质i液相密度，ρV，i为混合体系中物质i气相密度。
混合气相运动粘度、混合液相运动粘度采用以下公式：
其中xi表示混合体系中物质i的摩尔分数，ηi表示混合体系中物质i的纯物质性质。
混合液蒸发热采用
其中γi为混合液相体系中物质i的活度系数，本实施例采用UNIFAC活度系数模型。UNIFAC活度系数公式可从已公开的文献获取。
混合气PVT状态方程参数a,b采用以下公式
其中yi表示混合气中物质i的摩尔分数，ai，bi表示混合体系中物质i的纯物质PVT状态方程参数。
(4)获取被检测油品真实的单体烃和含氧化合物组成数据。本实施例中乙醇汽油可借助气相色谱分析获得单体烃组成数据，借助多维气相色谱分析获得含氧化合物组成数据。部分组成数据见表2(二甲苯分子组成数据)、表3(95#汽油分子组成数据)、表4(E10乙醇汽油分子组成数据)、表5(M15甲醇汽油分子组成数据)所示。
| 分子 | 质量百分数 |
| 甲苯 | 0.509 |
| 乙基苯 | 17.171 |
| 间二甲苯 | 51.133 |
| 对二甲苯 | 9.66 |
| 邻二甲苯 | 18.681 |
表2
| 分子 | 质量百分数 |
| 正戊烷 | 5.742 |
| 异戊烷 | 11.281 |
| 反-2-戊烯 | 1.428 |
| 环戊烷 | 2.747 |
| 正己烷 | 1.162 |
| 2-甲基戊烷 | 2.443 |
| 3-甲基戊烷 | 1.231 |
| 甲苯 | 16.105 |
| 正辛烷 | 0.522 |
| 2,2,4-三甲基戊烷 | 4.123 |
| 2,2,4-三甲基戊烷 | 0.563 |
| 2,2-二甲基己烷 | 1.461 |
| 2,5-二甲基己烷 | 1.003 |
| 2,2,3-三甲基戊烷 | 0.876 |
| 3,3-二甲基己烷 | 0.982 |
| 2,3,4-三甲基戊烷 | 3.473 |
| 乙基苯 | 1.438 |
| 间二甲苯 | 4.246 |
| 对二甲苯 | 1.213 |
| 邻二甲苯 | 2.253 |
| 1,2,4-三甲基苯 | 3.708 |
| 甲基叔丁基醚 | 6.469 |
表3
表4
表5
(5)基于步骤(4)获得的单体烃和含氧化合物，计算轻质石油馏分的馏程。具体包括以下两个步骤：
(5.1)根据真实的馏程分析设备，建立严格的蒸馏模型——液膜表面自动更新蒸馏模型。具体包括如下步骤：
(i)描述蒸馏瓶球部液相质量守恒和能量守恒以及初始条件。
(ii)描述蒸馏中所产生气相沿着蒸馏瓶轴向运动的一维动态蒸馏动量传递、质量传递、热量传递以及定解条件。
(iii)描述蒸馏中气相沿着瓶壁冷凝的一维动态蒸馏动量传递、质量传递、热量传递以及定解条件。
(iv)描述蒸馏中气相沿着瓶壁冷凝的质量守恒、能量守恒。其中质量守恒使得蒸馏中气相冷凝量等于冷凝液增加量。能量守恒使得冷凝释放潜热等于热量损失。
(v)描述蒸馏过程中的气液平衡。
(5.2)建立一种求解液膜表面自动更新蒸馏模型的数值算法——参数耦合迭代有限差分法。本实施例在求解模型过程中，采用内外循环嵌套程序，外循环求解在特定时间步长内的蒸发量以及蒸发气的组成和温度。内循环则联立求解步骤(5.1)中的各个步骤，初始参数根据外循环获得。内循环获得气相和冷凝液膜速度、组成、温度等状态变量随着时间变化的轴向分布。根据求解得组成和温度分布，不断迭代更新方程组参数，直到步骤(5.1)中方程组收敛。本实施例具体算法如图2所示。
本实施例最终得到轻质石油馏分的馏程曲线。其中，二甲苯馏程计算的平均绝对误差为0.4℃，95#汽油馏程计算的平均绝对误差为1.77℃，E10乙醇汽油馏程计算的平均绝对误差为2.33℃，M15甲醇汽油馏程计算的平均绝对误差为2.14℃。详细结果见表6(二甲苯馏程预测结果)、表7(95#汽油馏程预测结果)、表8(E10乙醇汽油馏程预测结果)表9(M15甲醇汽油馏程预测结果)。馏程曲线见图3、图4、图5、图6所示。模型计算误差小，均在国家标准试验方法所规定的再现性精度范围内，可在实际应用中使用。
| 馏程指标 | 初馏点 | T10 | T50 | T90 | 终馏点 |
| 预测值/℃ | 138.2 | 138.3 | 138.6 | 139.3 | 141.8 |
| 实测值/℃ | 138.1 | 138.5 | 139 | 139.9 | 141.1 |
| 误差/℃ | 0.1 | -0.2 | -0.4 | -0.6 | 0.7 |
表6
| 馏程指标 | 初馏点 | T10 | T20 | T30 | T40 | T50 | T60 | T70 | T80 | T90 | T95 | 终馏点 |
| 预测值/℃ | 36.5 | 45.5 | 55 | 65 | 76 | 86.4 | 96 | 110 | 132 | 165 | 180 | 195.6 |
| 实测值/℃ | 39.2 | 47.5 | 55.2 | 67.6 | 77.7 | 83.9 | 94.8 | 113.4 | 133.1 | 166.7 | 181.3 | 196.4 |
| 误差/℃ | -2.7 | -2.0 | -0.2 | -2.6 | -1.7 | 2.5 | 1.2 | -3.4 | -1.1 | -1.7 | -1.3 | -0.8 |
表7
| 馏程指标 | 初馏点 | T10 | T20 | T30 | T40 | T50 | T60 | T70 | T80 | T90 | T95 | 终馏点 |
| 预测值/℃ | 34.6 | 45.9 | 55.1 | 59.8 | 64.3 | 70.5 | 98.4 | 121.6 | 141.7 | 164.1 | 178 | 197.9 |
| 实测值/℃ | 37.3 | 48.7 | 57.2 | 62.8 | 67.8 | 70.3 | 101.0 | 118.7 | 143.2 | 162.5 | 177.0 | 201.9 |
| 误差/℃ | -2.7 | -2.8 | -2.1 | -3.0 | -3.5 | 0.2 | -2.6 | 2.9 | -1.5 | 1.6 | 1.0 | -4.0 |
表8
| 馏程指标 | 初馏点 | T5 | T10 | T20 | T30 | T40 | T50 | T60 | T70 | T80 | T90 | 终馏点 |
| 预测值/℃ | 32.3 | 38 | 42.4 | 44.4 | 46.2 | 49.7 | 53.5 | 57 | 80.6 | 120.1 | 181.3 | 204.4 |
| 实测值/℃ | 35.0 | 40.3 | 42.2 | 43.1 | 49.6 | 46.8 | 52.6 | 57.5 | 83.6 | 124.7 | 183.2 | 206.4 |
| 误差/℃ | 2.7 | 2.3 | -0.2 | -1.3 | 3.4 | -2.9 | -0.9 | 0.5 | 3.0 | 4.6 | 1.9 | 2.0 |
表9
以上的所述乃是本发明的具体实施例及所运用的技术原理，若依本发明的构想所作的改变，其所产生的功能作用仍未超出说明书及附图所涵盖的精神时，仍应属本发明的保护范围。
Claims (8)
1.一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于，包括如下步骤：
(1)汇合石油生产中所有轻质馏分的分子种类，建立涵盖轻质石油馏分和产品的分子组成数据库；
(2)基于分子组成数据库建立分子性质及温度参数数据库：
(3)基于分子组成数据库建立分子混合性质计算模型数据库；
(4)获取被检测油品真实的单体烃和含氧化合物组成数据，该数据可借助实验室分析获取，或者通过数学模型预测获取；
(5)基于步骤(4)获得的单体烃和含氧化合物，计算得到轻质石油馏分的馏程；具体如下：
(5.1)根据真实的馏程分析设备，建立严格的蒸馏模型——液膜表面自动更新蒸馏模型，得到模型的方程组；具体包括如下方程组：
(i)描述蒸馏瓶球部液相质量守恒和能量守恒微分方程组以及初始条件；
(ii)描述蒸馏中所产生气相沿着蒸馏瓶轴向运动的一维动态蒸馏动量传递、质量传递、热量传递方程以及定解条件；
(iii)描述蒸馏中气相沿着瓶壁冷凝的一维动态蒸馏动量传递、质量传递、热量传递方程以及定解条件；
(iv)描述蒸馏中气相沿着瓶壁冷凝的质量守恒、能量守恒方程；(v)描述蒸馏过程中的气液平衡方程；
(5.2)建立一种求解液膜表面自动更新蒸馏模型的数值算法——参数耦合迭代有限差分法；在求解模型过程中，采用内外循环嵌套程序，外循环求解在特定时间步长内的蒸发量以及蒸发气的组成和温度；内循环则联立求解步骤(5.1)中的方程组，方程组初始参数根据外循环获得；根据求解的结果，不断迭代更新方程组参数，直到步骤(5.1)中的方程收敛。
2.根据权利要求1所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的分子性质及温度参数数据库包括16项性质：分子气体密度、液体密度、摩尔质量、摩尔体积、表面张力、运动粘度、气体比热、液体比热、气体导热系数、液体导热系数、蒸发热、临界温度、临界压力、偏心因子、气相PVT状态方程参数、活度系数模型参数；其中，上述性质采用模型计算或查阅文献资料两种方法结合获得。
3.根据权利要求1所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的分子混合性质计算模型数据库包括混合气相密度、混合液相密度、混合气摩尔质量、混合液相摩尔体积、混合液相表面张力、混合气相运动粘度、混合液相运动粘度、混合气比热、混合液比热、混合气导热系数、混合液导热系数、混合液蒸发热、混合气PVT状态方程参数、混合液活度系数。
4.根据权利要求3所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的混合气相密度、混合液相密度、混合气摩尔质量、混合液相摩尔体积、混合气比热、混合液比热、混合气导热系数、混合液导热系数采用以下公式：
其中，xi表示混合体系中物质i的摩尔分数，fi表示混合体系中物质i的纯物质性质。
5.根据权利要求3所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的混合液相表面张力采用以下公式：
其中，[Pi]为混合体系中物质i等张比容，ρL，i为混合体系中物质i液相密度，ρV，i为混合体系中物质i气相密度。
6.根据权利要求3所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的混合气相运动粘度、混合液相运动粘度采用以下公式：
其中，xi表示混合体系中物质i的摩尔分数，ηi表示混合体系中物质i的纯物质性质。
7.根据权利要求3所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的混合液蒸发热采用以下公式：
其中，γi为混合液相体系中物质i的活度系数。
8.根据权利要求3所述的一种根据分子组成计算轻质石油馏分馏程的方法，其特征在于：所述的混合气PVT状态方程参数a,b采用以下公式计算：
其中，yi表示混合气中物质i的摩尔分数，ai，bi表示混合体系中物质i的纯物质PVT状态方程参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910221942.5A CN109839495B (zh) | 2019-03-22 | 2019-03-22 | 一种根据分子组成计算轻质石油馏分馏程的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910221942.5A CN109839495B (zh) | 2019-03-22 | 2019-03-22 | 一种根据分子组成计算轻质石油馏分馏程的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109839495A CN109839495A (zh) | 2019-06-04 |
| CN109839495B true CN109839495B (zh) | 2021-07-23 |
Family
ID=66886201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910221942.5A Active CN109839495B (zh) | 2019-03-22 | 2019-03-22 | 一种根据分子组成计算轻质石油馏分馏程的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109839495B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110459271B (zh) * | 2019-08-16 | 2023-04-18 | 青岛科技大学 | 一种纯组分和石油馏分苯胺点的估算方法 |
| CN111077288B (zh) * | 2019-12-06 | 2022-06-24 | 浙江中控技术股份有限公司 | 用于油品在线调合的馏程计算方法 |
| CN113987897B (zh) * | 2021-10-12 | 2023-04-14 | 广西大学 | 预测变压器用植物绝缘油物理化学性质的复合模型及应用 |
| CN114996636B (zh) * | 2022-05-27 | 2024-07-23 | 国家石油天然气管网集团有限公司 | 一种实时计算lng密度的方法及其系统 |
| CN115662532B (zh) * | 2022-10-26 | 2024-01-30 | 北京理工大学 | 一种宽域石油馏分通用表征燃料组分模板构建方法 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0336716A2 (en) * | 1988-04-06 | 1989-10-11 | Exxon Chemical Patents Inc. | Lactone modified viscosity modifiers useful in oleaginous compositions |
| EP0632267A1 (en) * | 1993-01-12 | 1995-01-04 | Nippon Fine Chemical Co., Ltd. | Method of measuring content of polycyclic aromatic compound in lanolin and removal thereof |
| CN1359001A (zh) * | 2000-10-17 | 2002-07-17 | 实验仪器科技有限公司 | 通过快速微蒸馏测定液体石油产品的蒸馏特性的方法和实施该方法所用的装置 |
| CN103091475A (zh) * | 2011-10-27 | 2013-05-08 | 中国石油化工股份有限公司 | 一种测定石油馏分油馏程的方法 |
| CN103091474A (zh) * | 2011-10-27 | 2013-05-08 | 中国石油化工股份有限公司 | 一种用于测定重质馏分油馏程的装置 |
| CN103509593A (zh) * | 2012-06-26 | 2014-01-15 | 中国石油化工股份有限公司 | 一种轻质烃油催化转化方法 |
| CN104020070A (zh) * | 2014-04-29 | 2014-09-03 | 辽宁万泓激光科技股份有限公司 | 煤焦油成分分析仪 |
| CN203811561U (zh) * | 2014-03-18 | 2014-09-03 | 大连瀚泰科技有限公司 | 馏程测定仪的加热单元自动控制系统 |
| CN204064397U (zh) * | 2014-11-11 | 2014-12-31 | 大庆市日上仪器制造有限公司 | 一种自动馏程液位测量装置 |
| CN105136940A (zh) * | 2015-10-15 | 2015-12-09 | 宣化钢铁集团有限责任公司 | 焦化轻苯初馏点和终馏点的测定方法 |
| CN205643216U (zh) * | 2016-03-08 | 2016-10-12 | 安徽天易金属新材料有限公司 | 一种新型石油产品馏程测定仪 |
| CN106324229A (zh) * | 2016-08-22 | 2017-01-11 | 何恺源 | 一种确定原油及石油流股详细分子组成的方法 |
| CN108920416A (zh) * | 2018-07-12 | 2018-11-30 | 廊坊会徕石油制品有限公司 | 一种油品调和参数计算方法 |
-
2019
- 2019-03-22 CN CN201910221942.5A patent/CN109839495B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0336716A2 (en) * | 1988-04-06 | 1989-10-11 | Exxon Chemical Patents Inc. | Lactone modified viscosity modifiers useful in oleaginous compositions |
| EP0632267A1 (en) * | 1993-01-12 | 1995-01-04 | Nippon Fine Chemical Co., Ltd. | Method of measuring content of polycyclic aromatic compound in lanolin and removal thereof |
| CN1359001A (zh) * | 2000-10-17 | 2002-07-17 | 实验仪器科技有限公司 | 通过快速微蒸馏测定液体石油产品的蒸馏特性的方法和实施该方法所用的装置 |
| CN103091475A (zh) * | 2011-10-27 | 2013-05-08 | 中国石油化工股份有限公司 | 一种测定石油馏分油馏程的方法 |
| CN103091474A (zh) * | 2011-10-27 | 2013-05-08 | 中国石油化工股份有限公司 | 一种用于测定重质馏分油馏程的装置 |
| CN103509593A (zh) * | 2012-06-26 | 2014-01-15 | 中国石油化工股份有限公司 | 一种轻质烃油催化转化方法 |
| CN203811561U (zh) * | 2014-03-18 | 2014-09-03 | 大连瀚泰科技有限公司 | 馏程测定仪的加热单元自动控制系统 |
| CN104020070A (zh) * | 2014-04-29 | 2014-09-03 | 辽宁万泓激光科技股份有限公司 | 煤焦油成分分析仪 |
| CN204064397U (zh) * | 2014-11-11 | 2014-12-31 | 大庆市日上仪器制造有限公司 | 一种自动馏程液位测量装置 |
| CN105136940A (zh) * | 2015-10-15 | 2015-12-09 | 宣化钢铁集团有限责任公司 | 焦化轻苯初馏点和终馏点的测定方法 |
| CN205643216U (zh) * | 2016-03-08 | 2016-10-12 | 安徽天易金属新材料有限公司 | 一种新型石油产品馏程测定仪 |
| CN106324229A (zh) * | 2016-08-22 | 2017-01-11 | 何恺源 | 一种确定原油及石油流股详细分子组成的方法 |
| CN108920416A (zh) * | 2018-07-12 | 2018-11-30 | 廊坊会徕石油制品有限公司 | 一种油品调和参数计算方法 |
Non-Patent Citations (4)
| Title |
|---|
| 炼油装置流程模拟及优化系统设计与应用;杨小健 等;《石油炼制与化工》;20110930;第87-91页 * |
| 用ASPEN PLUS软件模拟优化卡宾达原油常压蒸馏的研究;赵晓军 等;《炼油技术与工程》;20051130;第43-48页 * |
| 石油产品馏程测定法;中华人民共和国国家标准;《中华人民共和国国家标准》;19771108;第1-4页 * |
| 色谱模拟蒸馏方法测定蒽油的馏程和初馏点;马振元 等;《煤化工》;20130630;第26-28页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109839495A (zh) | 2019-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109839495B (zh) | 一种根据分子组成计算轻质石油馏分馏程的方法 | |
| Smith et al. | Improvements in the measurement of distillation curves. 3. Application to gasoline and gasoline+ methanol mixtures | |
| Knop et al. | A linear-by-mole blending rule for octane numbers of n-heptane/iso-octane/toluene mixtures | |
| Takeshita et al. | Influence of solvent addition on the physicochemical properties of Brazilian gasoline | |
| De Menezes et al. | Addition of an azeotropic ETBE/ethanol mixture in eurosuper-type gasolines | |
| Firoozabadi et al. | Reservoir depletion calculations for gas condensates using extended analyses in the Peng‐Robinson equation of state | |
| CN109859805B (zh) | 一种基于分子组成的汽油调和优化方法 | |
| JP2005512051A (ja) | 特定の分析データに一致するように計算された既知物質のブレンドとして未知物質を分析し、計算されたブレンドに基づいて未知物質の特性を予測する方法 | |
| Andersson et al. | On-line gas chromatographic analysis of higher alcohol synthesis products from syngas | |
| US6711532B1 (en) | Method and apparatus for predicting a distillation temperature range of a hydrocarbon-containing compound | |
| Mendes et al. | Reid vapor pressure prediction of automotive gasoline using distillation curves and multivariate calibration | |
| CN110431415B (zh) | 确定含烃的燃烧气体混合物的甲烷指数的方法 | |
| EP0984277A1 (en) | A method and apparatus for predicting a distillation temperature range of a hydrocarbon-containing compound | |
| Lanzer et al. | A predictive thermodynamic model for the Brazilian gasoline | |
| Zheng et al. | Application of the perturbed-chain SAFT to phase equilibria in the Fischer–Tropsch synthesis | |
| Trancoso et al. | Vapor–Liquid Equilibrium Measurements and Cubic-Plus-Association Modeling of Triethylene Glycol+ Water+ Methane Systems at 6.0 and 12.5 MPa | |
| Yakovleva et al. | Analysis of residual hydrocarbons and methanol in n-butane on a porous-layer capillary column with poly (1-trimethylsilyl-1-propyne) | |
| Harries et al. | A distillation approach to phase equilibrium measurements of multicomponent fluid mixtures | |
| Suri et al. | Application of phase-titrations for estimation of adulteration of gasoline and high-speed diesel with kerosene | |
| Vieira et al. | Evaluation of pressure on the crystallization of waxes using microcalorimetry | |
| Eliosa-Jiménez et al. | Vapor–liquid equilibrium data for the nitrogen+ n-octane system from (344.5 to 543.5) K and at pressures up to 50 MPa | |
| US20200103340A1 (en) | Systems and processes for performance property determination using optical spectral data | |
| Twu et al. | Determination of the solubilizing character of 1-methoxyethyl-1-methylpiperidinium tris (pentafluoroethyl) trifluorophosphate based on the Abraham solvation parameter model | |
| Spieksma | Prediction of ASTM method D86 distillation of gasolines and naphthas according to the fugacity-filmmodel from gas chromatographic detailed hydocarbon analysis | |
| CN111650324A (zh) | 一种烃类含量的在线检测方法 |
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 | ||
| CB02 | Change of applicant information |
Address after: Room 402, building F, No.98, Xiangxue No.8, Huangpu District, Guangzhou City, Guangdong Province Applicant after: Guangdong Xinfu Technology Co., Ltd Address before: 310053 room 412, 4 building, 66 Dongxin Avenue, Binjiang District, Hangzhou, Zhejiang. Applicant before: HANGZHOU SYSPETRO ENERGY TECHNOLOGY Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |