No abstract

Computer simulation of interphase mass transfer processes in natural hydrocarbon mixtures, as well as their qualitative laboratory studies, allow the specialist to study the processes in the oil reservoir more deeply. Today, the industry uses a number of software solutions (PVTSim, WinProp CMG, Aspen HYSYS, etc.), simulating interphase mass transfer in oils, gas condensate mixtures, etc. In most of these industrial software products, a method for calculating the phase equilibrium based on the equation of state of a multicomponent hydrocarbon system is implemented (flash calculations). For this purpose, as a rule, various versions of the Peng-Robinson (PR) and Soave-Redlich-Kwong (SRK) cubic equations of state are used. At the same time, in the last decades the search for a reliable and simultaneously simple equation of state is one of the most important directions of research (Amao 2014; Yan 2011; Ahmed 2010; Brusilovsky 2002; Gross 2001; Brusilovsky 1992). One of the modifications of cubic equations of state is the Brusilovsky equation (Brusilovsky 2002; Brusilovsky 1992). When developing this equation, the task implied the formulation of a cubic equation of state describing the PVT properties of natural hydrocarbon mixtures at pressures up to 100 MPa and temperatures up to 200° with sufficient accuracy for engineering purposes. The existing equations of state were reduced to a single form and a new generalized equation of state was proposed (Brusilovsky 2002; Brusilovsky 1992). The aim of the present paper is to identify the advantages of using a particular equation of state (PR, SRK, Brusilovsky) to describe the vapor-liquid equilibrium of hydrocarbon mixtures under various thermobaric conditions. To achieve this goal, a software in Python was designed that implements an iterative algorithm for calculating the equilibrium ratios (K-values) of components for the multicomponent hydrocarbon mixtures (flash calculations) using the Brusilovsky equation of state. The designed software is available for public access at the Internet under the link: https://gist.github.com/GrushnikovIU/982ba4dff03c0fa1fa5753b590a477d2. Since the PR and SRK equations of state can be considered as a special case of the Brusilovsky equation of state, the designed software made it possible to compare the K-values calculated using different equations of state. For some mixtures, the calculated K-values were compared with the known experimental K-values (Kogan 1966). It has been revealed in which cases the use of the generalized Brusilovsky equation of state can give advantages in comparison with the application of the PR and SRK equations of state. The results obtained in the present paper made it possible to formulate recommendations on the use of the considered equations of state (PR, SRK, Brusilovsky) in modeling the vapor-liquid equilibrium of hydrocarbon mixtures of various compositions. With the help of the software designed, additional studies were carried out (calculations of the mixtures’ composition of the vapor and liquid phases, the dew-point curve and the boiling curve) in order to demonstrate its capabilities.

<p>This work is devoted to developing the self-consistent thermo-hydro-chemo-mechanical reactive transport model to predict and describe natural and industrial petroleum processes at different scales.</p><p>We develop a version of the front tracking approach for multicomponent multiphase flow in order to treat spontaneous splitting of discontinuities. We revisit the solution for the Riemann problem and systematically classify all possible configurations as functions of initial concentrations on both sides of the discontinuity. We validate the algorithm against finite volume high-resolution technics and high-order spectral finite elements.</p><p>To calculate the parameters of phase equilibria, we utilize an approach based on the direct minimization of the Gibbs energy of a multicomponent mixture. This method ensures the consistency of the thermodynamic lookup tables. The core of the algorithm is the non-linear free-energy constrained minimization problem, formulated in the form of a linear programming problem by discretization in compositional space.</p><p>The impact of the complex rheological response of porous matrix on the morphology of fluid flow and shear deformation localization is considered. Channeling of porosity waves and shear bands morphology and their orientation is investigated for viscoelastoplastic both shear and bulk rheologies.</p>

We discuss numerical simulation of carbon dioxide injection considered by oil and gas companies. Complex behavior of multicomponent reservoir fluids mixed with carbon dioxide may cause the occurrence of vapor-liquid-liquid equilibria (VLLE), salt precipitation in aquifers, pore-clogging, etc. We propose a simple algorithm for phase equilibria calculations based on the minimization of the multicomponent system free energy. This algorithm can be used to calculate phase separations and component partitioning between the phases under various conditions (critical region, two- and three-phase equilibria, etc.). We demonstrate the applicability of the proposed algorithm in a series of calculations. We consider binary and ternary mixtures that include carbon dioxide and hydrocarbons. We examine the algorithm in two- and three-phase equilibrium calculations and compare its performance with the popular iterative fugacity equilibration technique. We show that both calculation techniques give near-identical results for the considered mixtures. Thus, we show that the free energy minimization algorithm can be used interchangeably with the fugacity equilibration technique for calculating phase equilibria. This algorithm is applicable for VLLE calculations, which is important when considering multicomponent reservoir fluids that include carbon dioxide.