用Lennard—Jones和Kihara位能模型预测真实流体的PVT性质

分别采用Ｌｅｎｎａｒｄ－Ｊｏｎｅｓ（Ｌ－Ｊ）和Ｋｉｈａｒａ位能模型并结合顺序解析法计算了径向分布函数，然后代入统计热力学所给出的严格理论状态方程来预测真实流体的ＰＶＴ性质。位能模型参数由纯流体的临界温度性质确定。计算结果表明，该方法可用于预测较度广的压力、温度范围内的非极性和弱极性纯流体的ＰＶＴ性质。讨还上于正己烷时，近临界区除外，由Ｌ－Ｊ模型预测的饱和液体体积的相对平均误差一般小于５％。对于分子     

纯组分高压流体的粘度模型

基于ｐＶＴ和Ｔμｐ图形的相似性,并结合两参数ＰｅｎｇＲｏｂｉｎｓｏｎ状态方程,建立了一个能够同时预测纯流体气、液相粘度的统一模型。该模型能够同时描述气、液相及超临界流体的粘度随温度、压力和组成的变化情况。模型的纯组分参数已普遍化为对比温度、对比压力和偏心因子的函数。在较宽的温度、压力范围内（Ｔ＝２８０～６００Ｋ,ｐ＝０．１～１００ＭＰａ）,对２２种烷烃、二氧化碳和氮气的粘度进行了计算,其相对平均误差为７．０１％．     

非极性及弱极性流体汽液平衡的微扰理论模型

建立了一个微扰理论模型来计算流体的汽液平衡。该模型将分子分为等直径的链节,以硬球流体为参考流体,以成链项、Ｌ－Ｊ作用项和偶极－偶极作用项为微扰项。计算了链烷烃、环烷烃、烯烃、芳烃和酮等共４８种非极性及弱极性纯流体的汽液平衡,并计算了１０个二元混合流体的汽液平衡,结果令人满意。     

纯组分高压流体的粘度模型

基于ｐ－Ｖ－Ｔ和Ｔ－μ－ｐ图形有相似性,并结合两参数Ｐｅｎｇ－Ｒｏｂｉｎｓｏｎ状态方程,建立了一个能够同时预测纯流体气、液相粘度的统一模型。该模型能够同时描述气、液相及超临界流体的粘度随温度、压力和组成的变化情况。模型的纯组分参数已普遍化为对比温度、对比压力和偏心因子的函数。在较宽的温度、压力范围内（Ｔ＝２８０￣６００Ｋ,ｐ＝０．１￣１００ＭＰａ）,对２２种烷烃、二氧化碳和氮气的粘度进行了计算,其     

Lennard-Jones 纯流体液汽界面的分子动力学模拟

在表面张力长程校正公式的基础上,应用正则系综分子动力学模拟方法较为系统地对Ｌｅｎｎａｒｄ－Ｊｏｎｅｓ（ＬＪ）纯流体的液汽平衡及界面性质进行了模拟研究。在温度Ｔ＊＝ｋＢＴ／ε为０．７５～１．２０范围内,分别计算了７个温度条件下ＬＪ纯流体的汽液相共存平衡密度、密度分布、表面张力和界面厚度。研究结果表明：文中得到的ＬＪ纯流体汽液相平衡密度与采用Ｇｉｂｂｓ系综ＭｏｎｔｅＣａｒｌｏ模拟计算结果基本一致;表面张力的模拟计算值与实验数据符合较好;另外,随着温度的升高,流体的表面张力降低而界面厚度变大。     

由饱和容积性质获取三参数范德华型状态方程的参数

本文认为使用由纯流体实际的饱和容积数据拟合的三个（温度）可调参数的立方型方程，可 以提高推算流体临界点附近区域的容积精度．文章提出了由纯流体饱和容积性质获取范德华型方 程参数的普遍化计算方法．所提出的计算方法物理意义明了，计算简便．文章还用纯ＣＯ２饱和容 积性质数据处理Ｐ－Ｔ和Ｍａｒｔｉｎ的方程．进一步说明此法获取状态方程参数的过程．     

极性流体的微扰理论研究

建立了一个用以计算极性流体热力学性质的微扰理论模型。模型中考虑了粒子间硬球排斥、色散、静电和诱导等相互作用能，由同时拟合饱和蒸汽压和液相密度获得分子微观参数。对几种常见的极性流体的计算结果表明，这一模型具有比较简单、准确、参数的物理意义明确等优点。     

用基团贡献法预测原油平衡闪蒸气化率

用正构烷烃模型分子法确定假组分的基团组成,根据ＵＮＩＦＡＣ基团贡献法饱和蒸气压模型与活度系数模型,以及修正的活度系数模型,计算了假组分的相平衡常数Ｋ,并预测了大庆原油平衡闪蒸气化率。结果表明,用平均沸点确定基团组成,ＲｉａｚｉＤａｕｂｅｒｔ关联式计算相对分子质量的ＲａｕＫ模型的预测值与实验值的偏差最小,为３．３６％;用ＡＰＩ式计算相对分子质量的ＫｉｋＫ模型和ＬａｒＫ模型次之,预测偏差低于５．５％,均可满足工程计算的要求。而且运用正构烷烃模型分子ＵＮＩＦＡＣ基团贡献法时,仅需平均沸点或相对分子质量之中的一个物性参数即可确定假组分的基团组成,计算过程更加简单。     

加压下混合流体汽液界面张力的分子热力学模型

建立了两个计算近临界点区域混合流体的汽液界面张力统计力学模型：一是将范德华理论１混合规则引入Ｋｉｒｋｗｏｏｄ－Ｂｕｆ纯流体界面张力模型;另一模型是建立在Ｄａｖｉｓ混合流体汽液界面张力模型。Ｄａｖｉｓ模型对单一纯组分流体可还原为Ｋｉｒｋｗｏｏｄ－Ｂｕｆ模型。对２１个纯组分流体的界面张力进行了关联,得到模型的分子参数,对纯流体平均标准差为１２０μＮ／ｍ,温度范围从常温到临界点,压力范围从常压到临界压力。应用两种模型计算了二氧化碳／丁烷混合体系的界面张力。     

链状流体的分子热力学模型：Ⅲ.混合物pVT和三元汽液平衡计算

将链状流体混合物的分子热力学模型应用于二元混合物ｐＶＴ数据关联和三元汽淮平衡数据的预测。结果表明，本模型用一个可调相互作用参数可以拟合二元混合物ｐＶ６数据。利用实验数据拟合得到二元相互作用参数，能令人满意地预测三元混合物的汽淮平衡数据。     

用第二音速维里系数确定分子间势能模型参数研究

利用工质的音速与第二维里系数的关系,推导出第二音速维里系数与分子间势能函数的关系式,提出了一种新的拟合分子间势能模型参数方法.并利用此方法,拟合得到了一些简单流体和制冷剂的的方阱、LJ(6,12)及Kihara等势能模型的参数.与传统方法相比,该法拟合势能参数具有更高的精度.     

GE型状态方程混合规则计算烃-水相平衡的结果对比

将PRSV方程分别与 5种有代表性的GE 型状态方程混合规则相结合 ,计算了烃水体系相平衡参数 ,并对计算结果进行了对比。轻质烃类水体系相平衡的计算结果表明 ,以范德华流体为参考态推导的Kurihara和Twu Coon混合规则无法正确描述烃水两相的组成 ,而Huron Vidal与Wong Sandler混合规则的计算偏差也较大 ,只有MHV2的计算精度与传统的使用组成依赖的交互作用参数的Kabadi混合规则相当。对正己烷水体系的三相平衡和临界轨迹线的计算表明 ,MHV2混合规则虽然在烃水互溶度的计算上具有优点 ,但在靠近临界点附近的计算还有改进的必要     

用音速确定HFC125的第二维里系数和分子间势能参数

利用ＨＦＣ１２５的音速和非线性最小二乘法中的Ｌｅｖｅｎｂｅｒｇ－Ｍａｒｑｕａｒｄｔ方法，拟合得到一个新的具有更高精度和更广泛用范围的ＨＣＦ１２５第二维里系数计算关联式，可以用于流体其他热物性参数的计算，从第二维里系数和结一音速维里系数及其与分子间势能函数的关系，提出了用第二音速维里系数确定分子间势能模型参数的新方法，并推导出用第二音速维里系数拟合方阱，ＬＪ（６，１２）及Ｋｉｌｈａｒａ等势能模型参数     

HFC-227ea热力学物性的分子动力学模拟

利用分子动力学模拟算法研究了七氟丙烷 (HF C-2 2 7ea)的热力学物性。通过将 HFC-2 2 7ea作为极性的两质点 L ennard-Jones流体处理建立了相应的势能模型。利用等温等压系综 +测试粒子 (NPT+ Test Particle)算法对 HF C-2 2 7ea的气液相平衡性质、过冷液和过热蒸气性质进行了模拟。模拟结果与美国国家标准研究所 (NIST)的数据库的最大相对偏差分别在 1.3 % (单相区 )和 1.6% (气液平衡区 )以内 ,表明该模拟方法已可用于物性的预测和工程实际     

真实流体的新状态方程

经统计力学椎证,建立新的流体状态方程。确定了新方程的气体形式,实现了气体参数普遍化。结果表明,该气体方程用于PVT性质计算,压力偏差小于3%、密度偏差小于5%的置信水平在91%以上;用于纯气体逸度系数计算,偏差小于5%的置信水平在97%以上。     

多维正态总体零均值的假设检验

讨论了在协方差已知或未知情况下多维正态总体零均值假设检验的问题。首先利用降维法 ,将多维正态总体转化为一维正态总体 ,然后利用T 统计量或U 统计量来判断拒绝或接受原假设。     

正庚烷临氢异构化研究进展

直馏汽油中一般含有较多的正庚烷 ,而其辛烷值为零 ,如能将其异构成带多个侧链的异庚烷 ,则可拓宽高辛烷值汽油理想组分的来源。正庚烷异构化的难点在于高选择性地生成多侧链异庚烷和抑制异庚烷的裂化反应。对正庚烷异构化的热力学、动力学、催化剂及反应的主要影响因素进行了分析 ,提出了结合反应器技术来研究正庚烷异构化的方法。从理论上来说 ,采用浆态床反应器优于传统的固定床     

从fcc到sc的相变模型：压力诱导相变

提出了一个从面心立方（ｆｃｃ）到简立方（ｓｃ）的相变模型，模型中的粒子之间的相互作用是势是各向异性的Ｌｅｎｎａｒｄ－Ｊｏｎｅｓ势，用ＭｏｎｔｅｔＣａｒｌｏ方法研究了该模型的有序无序转变以及相变温度与压力的关系，结果表明，序参量和晶格常数随温度而增大并在相变点有一个突变，亚稳态的存在和两相共存表明该相变是一个典型的一级相变，该模型可以唯象地描述固体Ｃ５０的相变特征。     

油气藏流体-CO2体系相行为的实验测定与计算

利用高压相平衡实验装置 ,对某油田F6 6A井和F30井的原油在油藏温度下注入不同量的CO2 (注入CO2 与油藏流体的摩尔比为 0～ 2 .0 0 )时的PVT关系进行了测定 ,共取得PVT数据点 16 5个 ,泡点压力点和体积系数点各 11个。实验结果表明 ,由于原始油样的组分结构和组成不同 ,其注气体系相行为有较大的不同。F6 6A井油藏原油的饱和 (泡点 )压力随着CO2 注入量的增加近似呈线性增加 ,而F30井的饱和压力和CO2 注入量之间呈二次曲线关系 ,上升速度比F6 6A井缓慢 ,注气体系的体积系数的增加速度与CO2 注入量近似呈直线关系。对油藏流体中各组分进行适当组合合并为假组分后 ,用PR和SRK状态方程对所测注CO2 原油体系的相行为进行了预测。结果表明 ,PR状态方程的预测精度优于SRK方程 ,但两者均不能用于计算高CO2 注入量时油藏流体的相行为     

一种新的二元Lennard-Jones链式流体互扩散系数计算公式

在前期工作提出的Lennard-Jones链式(LJC)流体自扩散系数计算公式的基础上,将计算平衡性质物性参数的范得瓦尔斯混合法则尝试应用于迁移性质的计算中,提出了一种新的计算二元LJC流体互扩散系数的计算公式.公式不包含二元可调参数,物理意义更加明确.通过不同温度、组分范围的12对二元LJC流体的270个实验点对新的方法进行了验证,新公式计算结果与实验值的平均绝对偏差为6.98%.与其他计算方法的比较结果显示,该计算方法具有更高的精度.