一种计算烃类化合物密度的方法与流程

本发明涉及一种计算烃类化合物密度的方法。
背景技术：
：石油中包含成千上万种化合物，在炼油工业中，原油及油品的密度在储运中也具有重要的意义，在原料及产品的计量以及炼油装置的设计等方面都是必不可少的。获取石油中各种类型化合物的密度，可为认识油品混合物的密度提供依据。获取化合物密度最直接的手段是通过实验方法测定，但有些化合物因分离提纯的困难，很难通过实验的手段获取其密度数据，有必要对其密度进行估算或预测。化合物密度与其分子结构密切相关，所以寻求化合物的结构-物性关联具有重要意义。烃类化合物是由碳和氢两种元素所组成的结构复杂的化合物，通常可分为烷烃、烯烃、炔烃、环烷烃、芳烃等几大类。在计算烃类化合物的基本物性时，如果仅将烃类化合物分成这几类物质来计算的话，不足以体现烃分子结构的复杂性，也得不到理想的计算效果。目前，建立结构-物性关联的常见方法是基团贡献法和拓扑方法，基团贡献法不能合理的外推至高碳数化合物的密度，对同分异构体的区分能力有限。拓扑方法对同分异构体有较好的区分能力，计算结果可靠，但使用范围相对较小(比如适合的碳数范围相对较小，不同类型的化合物需要不同的拓扑指数)。技术实现要素：本发明提出了一种计算烃类化合物密度的方法。本发明计算烃类化合物密度的方法包括以下步骤：(1)将烃类化合物中分子骨架相同而其骨架所连接的烷基不同的化合物归类；(2)按照步骤(1)的归类方法确定待计算化合物所属的类，选定此类化合物中一个已知其密度的化合物；(3)获取正构烷烃的密度D与表征分子结构的变量A的关联公式D＝f(A)；(4)利用与待计算化合物属于同类的已知化合物的密度数据D0及其表征分子结构的变量A的值A0，对步骤(3)中正构烷烃的密度D-表征分子结构的变量A的关联公式进行校正，得到该类化合物的密度D-表征分子结构的变量A的关联公式D＝f(A+A′)，其中A′为常数；(5)确定待计算化合物的表征分子结构的变量A的值，代入步骤(4)中得到的该类化合物的密度D-表征分子结构的变量A的关联公式D＝f(A+A′)，得到待计算化合物的密度。步骤(1)中优选将烃类化合物中分子骨架相同而其骨架所连接的直链烷基不同的化合物归类。例如，对于2-甲基丙烷、2-甲基丁烷、2-甲基戊烷、2-甲基己烷、…这一类化合物，可以将其归类为2-甲基烷烃类的化合物，显然归入此类的化合物均具有2-甲基烷烃的骨架，只是骨架上所连接的烷基不同，连接甲基的为2-甲基丙烷，连接乙基的为2-甲基丁烷，连接正丙基的为2-甲基戊烷，连接正丁基的为2-甲基己烷，…；对于甲苯、乙苯、丙苯、…这一类物质，可以将其归类为单烷基苯类的化合物；其它分子结构的烃类化合物，也按此方法归类。有的烃类化合物可能会归入多个表示不同分子骨架的类中。步骤(3)中所述的表征分子结构的变量A是体现分子结构特征的一个变量，例如可以为分子量、碳原子数或拓扑指数，优选采用分子量作为表征分子结构的变量。步骤(3)中所述的正构烷烃的密度D-表征分子结构的变量A的关联公式D＝f(A)可以采用已知的关联公式；也可以根据已知的正构烷烃的密度数据，选择一个表征分子结构的变量A，将二者关联而得到其关联公式D＝f(A)。下面举例说明：已知的沸点-碳原子数的关联公式(出自文献Kreglewski,A.andB.J.Zwolinski,ANewRelationforPhysicalPropertiesofn-Alkanesandn-AlkylCompounds.J.Phys.Chem.1961(65):1050-1052)的函数表达式为：ln(BP∞-BP)＝a-bCN2/3，式中BP∞、a、b为恒定参数，BP为沸点，CN为碳原子数；由于同一类型的烃类化合物随着碳原子数的递增，密度呈现较好的规律性，即密度趋近于正构烷烃类化合物，所以在计算烃类化合物的密度时，可以由上述的沸点-碳原子数的关联公式修改得到，即用因变量密度D取代上述关联公式中的因变量沸点BP，用自变量(CN+CN0)取代上述关联公式中的自变量CN，用参数变量γ来取代恒定常数2/3，并沿用其函数表达式，修改为：ln(D∞-D)＝p-n(CN+CN0)γ，式中D∞、p、n、CN0、γ为恒定参数，D为密度，CN为碳原子数；式中D∞、p、n、CN0、γ这些恒定参数可以由已知密度数值的一组烃类化合物关联得到，通常可以由一组已知其密度值的正构烷烃关联得到。由修正后的公式可以看出CN的修正CN0，可以保证侧链取代基对密度的影响，引入的参数变量γ可以保证同一类化合物随着碳原子数递增其密度呈现的规律性的影响，所以采用该关联式保证了其准确性。由于以拓扑指数或相对分子质量代替碳原子数作为表征分子结构的变量会更为准确，因此也可以将上述沸点-碳原子数的关联公式修改为密度-w指数的关联公式或密度-相对分子质量的关联公式：ln(D∞-D)＝p-n(W+W0)γ，式中D∞、p、n、W0、γ为恒定的参数，D为密度，w为wiener指数；或ln(D∞-D)＝p-n(MW+MW0)γ，式中D∞、p、n、MW0、γ为恒定的参数，D为密度，MW为分子量；上面三个关联公式中的恒定参数的值在每个公式中可能是不同的。同样地，关联方法可以选用回归分析方法、拟合方法、遗传算法、马夸尔特法，优选采用回归分析方法，所述回归分析方法包括最小二乘法、牛顿法、马夸特法和偏最小二乘法中的一种。当然这些关联公式在其关系表达式上会存在差异，比如函数关系及参数不同。由于正构烷烃的密度数据较为齐全，因此优选将正构烷烃化合物的密度与表征分子结构的变量相关联，得到密度与分子结构变量的关联公式。在步骤(4)中，首先选择一个与待计算化合物属于同类的已知化合物，确定其表征分子结构的变量A0、获取其密度值D0，可以通过已知的文献或通过实验方法或计算方法获得该值，然后引入校正值A′，将A0代替A、D0代替D代入正构烷烃的关联公式D＝f(A+A′)，对该公式进行校正，得 到计算该类化合物时表征分子结构的变量A的校正值A′，该值为一常数，将A′回代即得到该类化合物的密度D-表征分子结构的变量A的关联公式D＝f(A+A′)。在步骤(5)中，确定待计算化合物的表征分子结构的变量A的值，代入步骤(4)中得到的关联公式D＝f(A+A′)，得到待计算化合物的密度。本发明方法简单可靠，计算精度高，克服了基团贡献法和拓扑方法的不足，对同分异构体以及高碳数烃类化合物的计算较为准确。附图说明图1为本发明方法的逻辑框图。图2为本发明的3-甲基烷烃类化合物密度推测数据与文献数据对比图。图中正方形黑点为3-甲基烷烃类化合物的文献数据，圆形黑点为3-甲基辛烷的已知数据，平滑曲线为计算得到的3-甲基烷烃类化合物密度的关联公式。图3为本发明的2-甲基-1-烯烃类化合物密度推测数据与文献数据对比图。图中正方形黑点为2-甲基-1-烯烃类化合物的文献数据，圆形黑点为2-甲基-1-戊烯的已知数据，平滑曲线为计算得到的2-甲基-1-烯烃类化合物密度的关联公式。图4为本发明的2-甲基-1-烯烃类化合物密度误差分析。具体实施方式实施例1计算2-甲基烷烃类化合物的密度将2-甲基丙烷、2-甲基丁烷、2-甲基戊烷、2-甲基己烷等具有相同分子骨架而所连接烷基不同的化合物归入2-甲基烷烃类，由于2-甲基丁烷至2-甲基十三烷这些化合物的密度可以从已有文献中获取，因此利用本发明方法在计算这些化合物的密度后可以与文献数据相比较，以检验本发明方法的准确性。以2-甲基戊烷作为已知密度的同类化合物作校正，来计算此类化合物中其它物质的密度。计算方法如下：(1)将2-甲基丙烷、2-甲基丁烷、2-甲基戊烷、2-甲基己烷等其分子骨架相同而其骨架连接的烷基不同的化合物归类；(2)按照步骤(1)的归类方法确定待计算化合物属于2-甲基烷烃类，选定2-甲基戊烷作为已知密度的化合物，计算此类化合物中其它化合物的密度；(3)将已知的沸点-碳原子数的关联公式修改为密度-w指数的关联公式 ln(D∞-D)＝p-n(W+W0)γ，根据已知正构烷烃的密度数据，采用回归分析方法中的最小二乘法，对此公式进行回归分析，确定D∞、p、n、W0、γ的参数值，得到密度-w的关联公式为ln(0.83592-D)＝1.6351-2.60621(w+1.15777)0.08127；(4)将2-甲基戊烷作为2-甲基烷烃类已知密度的化合物，2-甲基戊烷的密度为0.6532g/cm3，其w指数为32，将这些值代入正构烷烃的密度-w指数的关联公式进行校正，此关联公式引入了校正参数w′：ln(0.83592-D)＝1.6351-2.60621(w+1.15777+w′)0.08127确定w′的值为-13.83277，将其回代入上式，得到计算2-甲基烷烃类化合物密度的关联公式ln(0.83592-D)＝1.6351-2.60621(w-12.675)0.08127(5)根据步骤(4)中得到的2-甲基烷烃类化合物密度的关联公式，计算2-甲基烷烃类化合物中除了2-甲基戊烷以外其它化合物的密度。对于2-甲基烷烃类化合物密度的计算结果见表1。计算结果表明本发明方法具有较高的精度。表12-甲基烷烃类化合物密度数据及误差分析名称wiener指数密度/(g/cm3)密度计算值/g/cm3)绝对误差2-甲基丁烷180.619720.613123-0.00662-甲基戊烷320.65320.649653-0.003552-甲基己烷520.678690.67655-0.002142-甲基庚烷790.69870.69718-0.001522-甲基辛烷1140.71350.7135-4.2E-092-甲基壬烷1580.72810.726729-0.001372-甲基癸烷2120.73690.7376650.0007652-甲基十一烷2770.74580.7468530.0010532-甲基十二烷3540.7540.7546770.0006772-甲基十三烷4440.76030.7614180.001118、、、、、、、、、、、、、、、、、、、、、实施例2计算3-甲基烷烃类化合物的密度将3-甲基戊烷、3-甲基己烷、3-甲基庚烷、…、3-甲基二十一烷等具有相同分子骨架而所连接烷基不同的化合物归入3-甲基烷烃类，由于3-甲基戊烷至3-甲基二十一烷这些化合物的密度可以从已有文献中获取，因此利用本发明方法在计算这些化合物的密度后可以与文献数据相比较，以检验本发明方法的准确性。此类化合物的分子骨架上所连接烷基的碳数大于二十一以上的化合物也可以采用本发明方法计算，只是没有相应的数据用来检验计算结果的准确性。以3-甲基辛烷作为已知密度的同类化合物作校正，来计算此类化合物中其它物质的密度。计算方法如下：(1)将3-甲基戊烷、3-甲基己烷、3-甲基庚烷、…、3-甲基二十一烷等其分子骨架相同而其骨架连接的烷基不同的化合物归类；(2)按照步骤(1)的归类方法确定待计算化合物属于3-甲基烷烃类，选定3-甲基辛烷作为已知密度的化合物，计算此类化合物中其它化合物的密度；(3)将已知的沸点-碳原子数的关联公式修改为密度-相对分子质量的关联公式ln(D∞-D)＝p-n(MW+MW0)γ，根据已知正构烷烃的密度数据，采用最小二乘法对此公式进行回归分析，确定D∞、p、n、MW0、γ的参数值，得到密度-w的关联公式为：ln(0.83952-D)＝2.56254-1.67274(MW-4.71336)0.21298(4)将3-甲基辛烷作为3-甲基烷烃类已知密度的化合物，3-甲基辛烷的密度为0.7209g/cm3，其相对分子质量为128.26，将这些值代入正构烷烃的密度-相对分子质量的关联公式进行校正，此关联公式引入了校正参数MW′：ln(0.83952-D)＝2.56254-1.67274(MW-4.71336+MW′)0.21298确定MW′的值为4.20096，将其回代入上式，得到计算3-甲基烷烃类化合物密度的关联公式ln(0.83952-D)＝2.56254-1.67274(MW-0.5124)0.21298(5)根据步骤(4)中得到的3-甲基烷烃类化合物密度的关联公式，计算3-甲基烷烃类化合物中除了3-甲基辛烷以外其它化合物的密度。对于3-甲基烷烃类化合物密度的计算结果见表2及图2。计算结果表明本发明方法具有较高的精度。表23-甲基烷烃类化合物的密度理论值与计算值及误差实施例3计算2-甲基-1-烯烃类化合物的密度将2-甲基-1-丁烯、2-甲基-1-戊烯、2-甲基-1-己烯等具有相同分子骨架而所连接烷基不同的化合物归入2-甲基-1-烯烃类，由于2-甲基-1-丁烯至2-甲基-1- 十五烯这些化合物的密度可以从已有文献中获取，因此利用本发明方法在计算这些化合物的密度后可以与文献数据相比较，以检验本发明方法的准确性。以2-甲基-1-戊烯作为已知密度的同类化合物作校正，来计算此类化合物中其它物质的密度。计算方法如下：(1)将2-甲基-1-丁烯、2-甲基-1-戊烯、2-甲基-1-己烯等其分子骨架相同而其骨架连接的烷基不同的化合物归类；(2)按照步骤(1)的归类方法确定待计算化合物属于2-甲基-1-烯烃类，选定2-甲基-1-丁烯作为已知密度的化合物，计算此类化合物中其它化合物的密度；(3)将已知的沸点-碳原子数的关联公式修改为密度-碳原子数的关联公式ln(D∞-D)＝p-n(CN+CN0)γ，根据已知正构烷烃的密度数据，采用最小二乘法对此公式进行回归分析，确定D∞、p、n、CN0、γ的参数值，得到密度-碳原子数的关联公式为：ln(0.83952-D)＝2.566-2.9377(CN-0.18509)0.21298(4)将2-甲基-1-戊烯作为2-甲基-1-烯烃类已知密度的化合物，2-甲基-1-戊烯的密度为0.941599，其相对碳原子数为6，将这些值代入正构烷烃的密度-碳原子数的关联公式进行校正，此关联公式引入了校正参数CN′ln(0.83952-D)＝2.566-2.9377(CN-0.18509+CN′)0.21298确定CN′的值为22.5056，将其回代入上式，得到计算2-甲基-1-烯烃类类化合物密度的关联公式ln(0.83952-D)＝2.566-2.9377(CN+0.756509)0.21298；(5)根据步骤(4)中得到的2-甲基-1-烯烃类化合物密度的关联公式，计算2-甲基-1-烯烃类化合物中除了2-甲基-1-戊烯以外其它化合物的密度。对于2-甲基-1-烯烃类化合物密度的计算结果见图3及图4。计算结果表明本发明方法具有较高的精度。当前第1页1 2 3