Application of the extended LJ potential-based equation of state to predict the density of five different classes of refrigerant systems (HCFCs,HFCs, HFEs,PFAs, and PFAAs)

The major goal of this work is to apply the extended Lennard-Jones potential-based equation of state (ELJ-based EoS) to predict the density of five different classes of refrigerant systems including chlorofluorocarbons, hydrochlorofluorocarbons, hydrofluorocarbons, hydrofluoroethers, perfluoroalkanes, and perfluoroalkylalkane. This EoS is based on an effective near-neighbor pair potential of the Lennard-Jones (12,6,3) type. The temperature dependencies of the parameters of the equation of state can be calculated at any temperature for each refrigerant. The calculated parameters along with the ELJ-based EoS have been used to calculate the density and isothermal compressibility coefficient of different refrigerants. A comparison between the predicted results and experimental data shows that the agreement is good. The total absolute average deviation of density for 14,871 data was found to be 0.34 compared with experimental data. Comparisons with the other EoSs show that the ELJ-based EoS is more accurate than other EoSs for most of the studied refrigerants.     

An accurate empirical correlation for predicting natural gas compressibility factors

The compressibility factor of natural gas is an important parameter in many gas and petroleum engineering calculations. This study presents a new empirical model for quick calculation of natural gas compressibility factors. The model was derived from 5844 experimental data of compressibility factors for a range of pseudo reduced pressures from 0.01 to 15 and pseudo reduced temperatures from 1 to 3. The accuracy of the new empirical correlation has been compared with commonly used existing methods. The comparison indicates the superiority of the new empirical model over the other methods used to calculate compressibility factor of natural gas with average absolute relative deviation percent (AARD%) of 0.6535.     

Rapidly Estimating Natural Gas Compressibility Factor

Natural gases containing sour components exhibit different gas compressibility factor(Z) behavior than do sweet gases.Therefore,a new accurate method should be developed to account for these differences.Several methods are available today for calculating the Z-factor from an cquation of state. However,these equations are more complex than the foregoing correlations,involving a large number of parameters,which require more complicated and longer computations.The aim of this study is to develop a simplified calculation method for a rapid estimating Z-factor for sour natural gases containing as much as 90% total acid gas.In this article,two new correlations are first presented for calculating the pseudo- critical pressure and temperature of the gas mixture as a function of the gas specific gravity.Then,a simple correlation on the basis of the standard gas compressibility factor chart is introduced for a quick estimation of sweet gases' compressibility factor as a function of reduced pressure and temperature.Finally,a new corrective term related to the mole fractions of carbon dioxide and hydrogen sulfide is developed.     

Calculation of isentropic compressibility and sound velocity in two-phase fluids

Derivative properties from equations of state (EoS) are well defined for homogeneous fluid systems. However, some of these properties, such as isothermal and isentropic (or adiabatic) compressibilities and sound velocity need to be calculated at conditions for which a homogeneous fluid splits into two (or more) phases, liquid or vapor. The isentropic compressibility and sound velocity of thermodynamically equilibrated fluids exhibit important discontinuities at phase boundaries, as noticed long ago by Landau and Lifschitz in the case of pure fluids. In this work, the two-phase isentropic compressibility (or inverse bulk modulus) is expressed in terms of the two-phase isothermal compressibility, two-phase thermal expansivity and an apparent heat capacity, defined as the partial derivative of total enthalpy with respect to temperature at constant pressure and composition. The proposed method is simple (simpler than previous approaches), easy to implement and versatile; it is not EoS-dependent and it requires only a flash routine and the expression of total enthalpy at given pressure, temperature and composition. Our approach is applied to a variety of fluid systems representative of reservoir applications and geophysical situations, including petroleum fluids (oil and gas condensate) and mixtures of water and gas (methane or CO₂). For low gas content in the two-phase fluid, i.e., near bubble point conditions, we obtain significantly lower bulk moduli and sound velocities than predicted within Wood's conventional approach, in which the liquid and gas phases are considered to be “frozen” at the passage of the acoustic wave.     

Bridging continuum and statistical thermodynamics via equations of state and the density of states

The connection between molecular force fields and equations of state (EoS) is typically established at the level of predicted quantities, e.g., by comparing simulation data and EoS data. In this paper we show how an EoS can be used to extract the density of states (Omega) of a system thus establishing a deeper connection between EoSs and statistical thermodynamics. We also show how such an EoS Omega can be used to aid molecular simulation methods designed to map out Omega (like the multicanonical approach). Central to the implementation of these ideas is the fact that the configurational Omega is related to thermodynamic properties accessible by an EoS via Boltzmann's equation. Sample calculations are presented for the Omega relevant to isothermal-isobaric and grand canonical ensemble simulations using the hard-sphere system and the Lennard-Jones system as model fluids, and the Carnahan-Starling EoS and a cubic EoS, respectively, as thermodynamic models.     

Vapor liquid equilibrium modeling of alkane systems with Equations of State: “Simplicity versus complexity”

Two types of Equations of State (EoS), which are characterized here as “simple” and “complex” EoS, are evaluated in this study. The “simple” type involves two versions of the Peng–Robinson (PR) EoS: the traditional one that utilizes the experimental critical properties and the acentric factor and the other, referred to as PR-fitted (PR-f), where these parameters are determined by fitting pure compound vapor pressure and saturated liquid volume data. As “complex” EoS in this study are characterized the EoS derived from statistical mechanics considerations and involve the Sanchez–Lacombe (SL) EoS and two versions of the Statistical Associating Fluid Theory (SAFT) EoS, the original and the Perturbed-Chain SAFT (PC-SAFT).

The evaluation of these two types of EoS is carried out with respect to their performance in the prediction and correlation of vapor liquid equilibria in binary and multicomponent mixtures of methane or ethane with alkanes of various degree of asymmetry. It is concluded that for this kind of systems complexity offers no significant advantages over simplicity. Furthermore, the results obtained with the PR-f EoS, especially those for multicomponent systems that are encountered in practice, even with the use of zero binary interaction parameters, indicate that this EoS may become a powerful tool for reservoir fluid phase equilibria modeling.     

环脲硝胺化合物的性能和电子结构之间的关系I. 羰基伸缩频率和键级、双原子作用能的近似线性相关

本文通过CNDO/2计算研究环脲硝胺同系化合C=O基伸缩频率与所得计算参量的关系。获得的结果是: 羰基伸缩频率实验值与C=O的Mulliken键级、双原子作用能等计算值之间, 存在近似的线性关系, 对现有的9个文献化合物和新合成的2,4,6,8-四硝基半甘脲的~γCO和计算参量作了线性回归处理, 相关系数为0.95。讨论了一些理论有机的传统观点与本文结果的关系。     

Measurement and calculation of gas compressibility factor for condensate gas and natural gas under pressure up to 116 MPa

The volumetric properties of two reservoir fluid samples collected from one condensate gas well and one natural gas well were measured under four groups of temperatures, respectively, with pressure up to 116 MPa. For the two samples examined, the experimental results show that the gas compressibility factor increases with the increase of pressure. But the influence of the temperature is related to the range of the experimental pressure. It approximately decreases with the increase of temperature when the pressure is larger than (45 to 50) MPa, while there is the opposite trend when the pressure is lower than (45 to 50) MPa. The dew point pressure was also determined for the condensate gas sample, which decreases with the increase of temperature. The capabilities of four empirical correlations and a thermodynamic model based on equation of state for describing gas compressibility factor of reservoir fluids under high pressure were investigated. The comparison results show that the thermodynamic model recommended is the most suitable for fluids whatever produced from high-pressure reservoirs or conventional mild-pressure reservoirs.     

Efficient methods for calculations of compressibility,density and viscosity of natural gases

This study presents two new methods for calculating properties of natural gases. The first is an efficient empirical model to calculate compressibility and density of natural gases containing high amount of heptane plus and none-hydrocarbon components. The model is derived from 2400 measurements of compressibility and density of various gases presented in this study. Accuracy of the model is compared to various equations of state (EOS), corresponding state, and empirical methods. The study shows that the new model is simpler and more efficient than EOS. It eliminates the numerous computations involved in EOS calculations. The new method also eliminates the characterization of the heptane plus fraction and estimation of binary interaction parameters needed for EOS calculations. Experimentally measured density of several gases has been used to study the validity of the proposed method. These measurements indicate that the new method successfully capture the physical trend of changing gas density as a function of pressure, temperature, and composition.The second method is a modification of Lee–Gonzalez–Eakin gas viscosity correlation. The new method accounts for the presence of heptane plus, hydrogen sulfide, and carbon dioxide in natural gases. The proposed method is compared to other EOS-based viscosity model, corresponding state methods, and correlations. The comparison indicates the superiority of the new method over the other methods used to calculate viscosity of natural gases.     

Thermodynamic derivative properties and densities for hyperbaric gas condensates: SRK equation of state predictions versus Monte Carlo data

The ability of Soave–Redlich–Kwong cubic equation of state (SRK EoS) to predict densities and thermodynamic derivative properties such as thermal expansivity, isothermal compressibility, calorific capacity, and Joule–Thompson coefficients, for two gas condensates over a wide range of pressures (up to 110 MPa) was studied. The predictions of the EoS were compared to Monte Carlo simulation data obtained by Lagache et al. [M.H. Lagache, P. Ungerer, A. Boutin, Fluid Phase Equilibr. 220 (2004) 221]. Two completely different alpha functions for the SRK EoS attractive term were used and their respective effects on the predictions of such properties were analyzed. Also, two different forms of the crossed terms of the attractive parameter, a_ij, and three expressions of the crossed terms of the repulsive parameter, b_ij, were combined in different ways, and predictions were carried out. Little sensitivity of the properties on the chosen alpha function, except for the calorific capacities, was found in the systems studied. The most commonly used combination rules to model phase behavior of reservoir fluids, i.e. geometric and arithmetic forms of a_ij and b_ij, respectively, predicted very deficient results for these fluids at extreme conditions, specially for density calculations.     

An efficient correlation for calculating compressibility factor of natural gases

Compressibility factor (z-factor) values of natural gases are necessary in most petroleum engineering calculations. Necessity arises when there are few available experimental data for the required composition, pressure and temperature conditions. One of the most common methods of calculating z-factor values is empirical correlation. Firstly, a new correlation based on the famous Standing-Katz (S-K) Chart is presented to predict z-factor values. The advantage of this correlation is that it is explicit inzand thus does not require an iterative solution as is required by other methods. Secondly, the comparison between new one and other correlations is carried out and the results indicate the superiority of the new correlation over the other correlations used to calculate z-factor.     

A new approach in modelling phase equilibria and gas solubility in electrolyte solutions and its applications to gas hydrates

This paper presents a new predictive model for phase equilibria and gas solubility calculations in the presence of electrolyte solutions. It treats salts as pseudo-components in an equation of state (EoS) by defining the critical properties and acentric factor for each salt. The water–salt, gas–salt and salt–salt binary interaction parameters (BIP) have been determined by using available experimental data on freezing point depression and boiling point elevation as well as gas solubility and salt solubility data in saline solutions.The methodology has been applied in modelling sodium chloride, potassium chloride and their mixtures, as well as solubility of methane and carbon dioxide in aqueous single and mixed electrolyte solutions.The developed model is capable of accurately predicting the phase behaviour, gas hydrate stability zone and potential salt precipitation in single and mixed electrolyte solutions. The model predictions are compared with available independent experimental data, including hydrate inhibition characteristics of single and mixed electrolyte solutions, and good agreement is demonstrated.     

Calculations of nuclear quadrupole coupling in noble gas-noble metal fluorides: interplay of relativistic and electron correlation effects

The nuclear quadrupole coupling constants (NQCCs) of noble gas and noble metal nuclei in the recently found noble gas-noble metal fluorides (NgMF, where Ng=Ar,Kr,Xe and M=Cu,Ag,Au) are obtained theoretically by high-level ab initio calculations, where both relativistic and electron correlation effects are included, and compared to experimental results. Fully relativistic four-component Dirac-Hartree-Fock (DHF) calculations are carried out at the basis set limit for electric field gradient that couples with the electric quadrupole moment of the nucleus, and uncorrelated relativistic effects are extracted by comparing DHF results to nonrelativistic (NR) HF calculations. Electron correlation effects are investigated both at fully relativistic second-order Moller-Plesset (DMP2) and at NR MP2 levels of theory, as well as at the NR coupled-cluster singles and doubles with perturbational triples [CCSD(T)] level. The validity of the approximation where relativistic effects, on the one hand, and nonrelativistically obtained correlation effects, on the other hand, are evaluated separately and assumed to be additive, is investigated by comparison with the DMP2 results. Inclusion of relativistic effects is shown to be necessary for obtaining the correct NQCC trends as the nucleus of interest and/or its neighbors become heavier. Electron correlation treatment is needed for approaching quantitative agreement with the experimental NQCCs. The assumption of additive electron correlation and relativistic effects, corresponding to the NR correlation treatment added on top of relativistic DHF data, gives qualitatively correct noble gas NQCCs. For noble metal NQCCs, correlation treatment at the relativistic level of theory is mandatory for reaching agreement with experimental results. Current work also confirms the experimental trends of NQCCs, which have been taken as an indication of nearly covalent interaction between noble gas and noble metal in the heaviest present systems, especially in XeAuF.     

Density functional theory-based prediction of the formation constants of complexes of ammonia in aqueous solution: indications of the role of relativistic effects in the solution chemistry of gold(I)

A prediction of the formation constants (log K1) for complexes of metal ions with a single NH3 ligand in aqueous solution, using quantum mechanical calculations, is reported. DeltaG values at 298 K in the gas phase for eq 1 (DeltaG(DFT)) were calculated for 34 metal ions using density functional theory (DFT), with the expectation that these would correlate with the free energy of complex formation in aqueous solution (DeltaG(aq)). [M(H2O)6]n+(g) + NH(3)(g) = [M(H2O)5NH3]n+(g) + H2O(g) (eq 1). The DeltaG(aq) values include the effects of complex changes in solvation on complex formation, which are not included in eq 1. It was anticipated that such changes in solvation would be constant or vary systematically with changes in the log K(1) value for different metal ions; therefore, simple correlations between DeltaG(DFT) and DeltaG(aq) were sought. The bulk of the log K1(NH3) values used to calculate DeltaG(aq) were not experimental, but estimated previously (Hancock 1978, 1980) from a variety of empirical correlations. Separate linear correlations between DeltaG(DFT) and DeltaG(aq) for metal ions of different charges (M2+, M3+, and M4+) were found. In plots of DeltaG(DFT) versus DeltaG(aq), the slopes ranged from 2.201 for M2+ ions down to 1.076 for M4+ ions, with intercepts increasing from M2+ to M4+ ions. Two separate correlations occurred for the M3+ ions, which appeared to correspond to small metal ions with a coordination number (CN) of 6 and to large metal ions with a higher CN in the vicinity of 7-9. The good correlation coefficients (R) in the range of 0.97-0.99 for all these separate correlations suggest that the approach used here may be the basis for future predictions of aqueous phase chemistry that would otherwise be experimentally inaccessible. Thus, the log K1(NH3) value for the transuranic Lr3+, which has a half-life of 3.6 h in its most stable isotope, is predicted to be 1.46. These calculations should also lead to a greater insight into the factors governing complex formation in aqueous solution. All of the above DFT calculations involved corrections for scalar relativistic effects (RE). Au has been described (Koltsoyannis 1997) as a "relativistic element". The chief effect of RE for group 11 ions is to favor linear coordination geometry and greatly increase covalence in the M-L bond. The correlation for M+ ions (H+, Cu+, Ag+, Au+) involved the preferred linear coordination of the [M(H2O)2]+ complexes, so that the DFT calculations of DeltaG for the gas-phase reaction in eq 2 were carried out for M = H+, Cu+, Ag+, and Au+. [M(H2O)2]+(g) + NH3(g) = [M(H2O)NH3]+(g) + H2O(g) (eq 2). Additional DFT calculations for eq 2 were carried out omitting corrections for RE. These indicated, in the absence of RE, virtually no change in the log K1(NH3) value for H+, a small decrease for Cu+, and a larger decrease for Ag+. There would, however, be a very large decrease in the log K1(NH3) value for Au(I) from 9.8 (RE included) to 1.6 (RE omitted). These results suggest that much of "soft" acid behavior in aqueous solution in the hard and soft acid-base classification of Pearson may be the result of RE in the elements close to Au in the periodic table.     

Initial estimation of hydrate formation temperature of sweet natural gases based on new empirical correlation

Production,processing and transportation of natural gases can be significantly affected by clathrate hydrates.Knowing the gas analysis is crucial to predict the right conditions for hydrate formation.Nevertheless,Katz gas gravity method can be used for initial estimation of hydrate formation temperature (HFT) under the circumstances of indeterminate gas composition.So far several correlations have been proposed for gas gravity method,in which the most accurate and reliable one has belonged to Bahadori and Vuthaluru.The main objective of this study is to present a simple and yet accurate correlation for fast prediction of sweet natural gases HFT based on the fit to Katz gravity chart.By reviewing the error analysis results,one can discover that the new proposed correlation has the best estimation capability among the widely accepted existing correlations within the investigated range.     

量子化学MNDO法中键强度的描述

1 前言以CNDO 法为基础建立的键函数方法,把量子化学的成键强度与传统的成键和非成键相互作用的有关项直接联系了起来.它对σ、π体系都能作定量描述,并且有旋转不变性,满足对称性要求,能很好地随键长变化反映键强度情况,并作了很好的尝试.按照同样原理,本文用理论上已作了许多改进的MNDO 法建立原子对作用能E_(AB)来量度成键强度.理论上可以证明,MNDO 法的原子     

High-pressure vapor–liquid equilibria of systems containing ethylene glycol,water and methane: Experimental measurements and modeling

This work presents new experimental phase equilibrium measurements of the binary MEG–methane and the ternary MEG–water–methane system at low temperatures and high pressures which are of interest to applications related to natural gas processing. Emphasis is given to MEG and water solubility measurements in the gas phase. The CPA and SRK EoS, the latter using either conventional or EoS/G^E mixing rules are used to predict the solubility of the heavy components in the gas phase. It is concluded that CPA and SRK using the Huron–Vidal mixing rule perform equally satisfactory, while CPA requires fewer interaction parameters.     

Prediction of gas condensate properties by Esmaeilzadeh–Roshanfekr equation of state

The Esmaeilzadeh–Roshanfekr (ER) equation of state (EOS) is used to predict the PVT properties of gas condensate reservoir fluids. Three gas condensate fluid samples taken from three wells in a real field in Iran, referred here as SA1, SA4 and SA8, as well as five samples from literature have been used to check the validity of the ER EOS in calculating the PVT properties of gas condensate mixtures. Some experiments such as constant composition expansion (CCE), constant volume depletion (CVD) and dew point pressures are carried out on these samples. In order to have an unbiased comparison between the ER and the Peng–Robinson (PR) equation of state, van der Waals mixing rules are used without using any adjustable parameters (k_ij = 0). Also, no pure component parameters are adjusted. The critical properties and acentric factor for plus-fraction are estimated by the Kesler–Lee, Pedersen et al. and Riazi–Daubert characterization methods. The results of dew point pressure calculations show that the ER EOS has smaller error than the PR EOS. For some mixtures, relative volume, gas compressibility factor and condensate drop-out in CVD and CCE test were also predicted. Comparison results between experimental and calculated data indicate that the ER EOS has smaller error than the PR EOS. The total average absolute deviation was found to be 0.82% and 2.97% for calculating gas compressibility factor and gas specific gravity in CVD test. Also, the total average absolute deviation was found to be 2.06% and 3.42% for calculating gas compressibility factor and relative volume in CCE test.