## Confined Phase Envelope of Gas-Condensate Systems in Shale Rocks

#### STANISLAW NAGY*, JAKUB SIEMEK*

*AGH UNIVERSITY OF SCIENCE & TECHNOLOGY, DRILLING & OIL-GAS FACULTY, AL. A. MICKIEWICZA 30, 30-059 KRAKOW, POLAND.*

### Abstract

Natural gas from shales (NGS) and from tight rocks are one of the most important fossil energy resource in this and next decade. Significant increase in gas consumption, in all world regions, will be marked in the energy sector. The exploration of unconventional natural gas & oil reservoirs has been discussed recently in many conferences. This paper describes the complex phenomena related to the impact of adsorption and capillary condensation of gas-condensate systems in nanopores. New two phase saturation model and new algorithm for search capillary condensation area is discussed. The algorithm is based on the Modified Tangent Plane Criterion for Capillary Condensation (MTPCCC) is presented. The examples of shift of phase envelopes are presented for selected composition of gas-condensate systems.

## 1. Introduction

Extraction of natural gas from shale rocks, tight rocks (also production of coal bed methane) successfully implemented in the USA in the last two decades stimulate the American economy (IEA (2011, 2012), EIA 2013). The low cost of fuel and low carbon emission during combustion of natural gas induce high gas demand in the World in the next 30 years (IEA (2011, 2012)). Natural gas is also an important energy source for the transition time towards a low-carbon economy in the World (IEA, 2011, 2012; MIT, 2011; Siemek et al., 2009; Nagy & Siemek, 2011; Siemek & Nagy 2012). The small European Union gas production and the high gas import dependency of Poland (Siemek et al., 2009; Nagy et al., 2009) is a motivation to search own way to increase country energy safety in next twenty years. The first results of exploration suggest wet/conden- sate or light oil systems in Baltic Basin as well as in the Lublin basin (Poprawa, 2010; Siemek & Nagy, 2012; PGI, 2012). This is also motivation to research potential efficiency of extraction dry shale gas (Klimkowski & Nagy, 2014) and to know thermodynamics and flow mechanism in two phase shale/tight systems. Large part of this paper has been presented in last year (Nagy & Siemek, 2013)1. The entire thermodynamic analysis of the reservoir hydrocarbon system is impossible because of the large influence of factors (see e.g. Nagy, 2002, 2003). The reservior engineering models of shale/light gas extraction usually omitt the condensation in parous media, focuset on “dry” gas (e.g. Rubin, 2010; Klimkowski & Nagy, 2014) The traditional attitude to flat vapor-liquid thermodynamics may be extended to curved surfaces with gravity and adsorption/ capillary condensation components and with new nano-effects in case of mezzo- and micropore reservoirs. The new effects related to critical phenomena in narrow pores will be significant for vapor- liquid equilibrium in tight and shale rocks.

## 2. State of the art of research of VLE in the porous media

The interface between phases (i.e. gas, condensate, brine) in a porous media is not flat. The capillarity and gravity segregation effects on the PVT and VLE properties of reservoir fluid are omitted In the classical thermodynamic analysis, which one may find in many works – see e.g. Ahmed (1989), Pedersen et al. (1989), Firoozabadi (1999), Whitson & Brule, 2000). In the range of laboratory PVT-VLE research an additional effect of capillarity is omitted in the case of saturation phase envelope. The pioneered Russian authors Trebin and Zadora (1968) reported a strong influence of the porous media on the dew point pressure and VLE of condensate system. The observed effect was about 10-15% increase of dew point pressure in this system. A similar observation was recorded by Sadyk and Zade (1963) and Tindy and Raynal (1966). The other American and Canadian research (Smith, Yarborough (1968), Sigmund et al. (1973)) suggest that there is no evident impact of capillarity on the saturation pressure. Smith and Yarborough (1968) indicate in their paper about wettability that the porous structure has little impact on the vaporization process of the liquid phase in the rock.

Sigmund and al. (1973) have analysed of laboratory works and concluded that the main reason of the hypothetical impact of porous structure on the saturation pressure was the lack of fluid circulation in the cell. More interesting discussion related to impact of porous media properties on VLE is presented in a paper of Ping et. al. (1996), Brusilovsky (1990), Firoozabadi (1999), Shapiro & Stanby (1996), Nagy (2002), Didar & Akkutlu (2013), Pang et al. (2012), Devegowda et al. (2012), Firincioglu et al. (2012).

In authors opinion the condition of necessary to form continuous liquid phase is oil wettability. The other papers see e.g. Cao Minh et al. (2012) confirm the existence of such type of wettability in the unconventional reservoirs, where up to 50% of pores are located in organic matrix (kerogen). The opinion of Lee (1989), Brusilovsky (1990), and Guo (1996), Firoozabadi (1999), Nagy (2002, 2003), Elamin et al. (2013), Didar & Akkutlu (2013), Whitson & Sunjerga (2012), Orangi (2011), Hartman et al. (2011), Pang et al (2012), Honarpour et. al. (2012), Kuppe F. et al., 2012. Some of the confirm necessary to include the capillary effects in the deep reservoir with the low permeability and oil-wet wettability. Additionally new effect is considered – shift of critical parameters (Didar & Akkutlu (2013), Elamin et al. (2013)). The range of pore diam- eters in the shale gas reservoir is between 1.5-100 nm (Clarkson & Haghshenas (2012), Kang (2011), Diaz eta al. (2009). Kang (2011), Clarkson et al. (2012) Diaz et al. (2009), Zhang et al. 2013, Sanaei et al. 2013. Firincioglu et al. (2012) have analysed porosity and pore structure of several shale reservoir. A molecular layer density for methane at 80°C across of a 3.6 nm organic slit-pore is given based on Diaz et al. (2009) work (Fig. 1). The Figure 1 has been obtained by Diaz et al. using molecular dynamics simulation carried out in the canonical (NVT) ensemble.

The estimated pore pressure for the bulk phase measurements of methane is 4,413 psi. Similar analysis is presented in Ma & Jamili (2014). The example of histogram of pore diameter from MICP experiment in the shale rock is given also in the Figure 2 (Clarkson & Haghshenas (2012)), similar histogram associated to different shale rock is given in Figure 3 (Dhanapal et al., 2014). The mentioned papers presented more detailed possible distribution of porous structures in shale and tight rocks. Based on IUPAC pore size distribution Alharthy et al (2013) proposed three different thermodynamic phase behaviours’ paths related to VLE conditions in porous media:

- Unconfined pore phase behaviour, which is mainly the unshifted phase behaviour in macropores or what would be in a PVT cell with no pore confinement effects. These occupy majority of the macro and mesopores (60-80%).
- Mid-confined pore phase behaviour, which is the partially shifted phase behavior in mesopores and is considered to be between unconfined and confined pore phase behaviors. This category occupies about 10-15% of the pores.
- Confined pore phase behavior, which is the shifted phase behavior in nanopores where the pore size is less than 3 nm and these occupy about 3-5% of pores and is dependent on the mineralogy (clay content).

The paper of Alharthy et. la (2013) also discussed two possible scenario relevant to condensation of hydrocarbon in dew phase envelope (Figs 4 & 5):

1. scenario without phase envelope shift: inside the nanopores, the capillary pressure is high and the gas phase pressure is higher than the original dew point pressure which causes the fluid condense sooner (capillary condensation (see Nagy, 2002)). This is in contrast to the meso and macro pores where at high pressures the supercritical fluid condenses the normal condensation path.

2. scenario with shift of phase envelope (shift of dew pressure): in case of favorable phase envelope shift, the process of condensation does not occur. The fluid in nanopores stays in the supercritical state with minimum liquid dropout. In the meso and macropores the fluid will condense sooner.

### 2.1. Hydrocarbon adsorption phenomena in the porous rocks in the near of saturation curve

The adsorption phenomena in the porous media may have a significant impact on the reserve distribution of tight, shale and coal-bed methane reservoirs. The adsorption process may largely distinguish from surface adsorption observed in the chemical labs. The main two differentiating reasons are: existence of capillary condensation phenomena in the narrow pores and possibility of flow access blocking in the porous network. The progress in the fundamentals of adsorption theory one may find in the Dąbrowski (2001) paper. The adsorption phenomena related to the porous media are discussed in many textbooks (ie. Defay, Prigogine (1966), Adamson (1990), Dullien (1992)).

The advances in the adsorption process in the high pressure porous media may be found in Shapiro, Stenby (1996, 2000, 2001), Guo et al., (1996) and Satik, Horne, Yortsos, (1995), Kang (2011) Altman et al. (2014), Travalloni et al. (2013) papers. In this paper the conservative model of adsorption has been used with the single molecular layer. The maximum thickness of adsorption film found in the real mesopore system is 23 Ǻ (Adamson, 1990), but in this work very conservative film of adsorption 4 Ǻ is assumed. The possible selective molecule adsorption to the kerogene in a shale gas condensate system is presented in Fig. 6 after Altman et al. (2014).

### 2.2. Capillary-adsorption vapor-liquid equilibrium (VLE) model

On the curved surface inside of porous media for the stationary state of the second kind (based upon Prigogine definition) following set of equation is valid (neglecting gravity force):

where T — temperature and m — chemical potential of i — component can be varied indepen- dently, that r is allowed to vary at the same time, t — i the critical adsorption film. At the dew point of capillary condensation phenomena exist three types of chemical potential:

Laplace equation (another proposal for modification of Laplace equation is presented in paper of Firincinoglu et al. (2012), in this work only classic equation has been used):

Interfacial multicomponent surface tension is calculated using (8) equation:

where σ — interfacial tension (IFT), Z V, Z L — vapor and liquid phase compressibility factors, E — empirical dimensionless constant, π_{i} — parachor of i-component (Danesh, 1991; Weinaug & Katz, 1943)).

The set of nonlinear equations may be modified by expanding using Taylor series and taking into consideration only first term, the chemical potential of i-component in the liquid phase may be calculated:

or using fugacity of i-component for isothermal conditions:

The derivative on the right side is a molar partial component volume:

where Vi is partial liquid molar volume o i-component. Equation (10) may be rewritten to the form:

Based upon the above derivation difference between bulk and capillary condensate dew point equations may be formulated:

Based upon equation the new capillary condensation equilibrium constant may be introduced:

or after simplification based of constant partial molar volume:

The proposed algorithm for calculation of the capillary condensation dew pressure is given in the paper of Nagy (2002) together with new form of Rachford_Rice equation and modified modified tangent plane criterion for capillary condensation (MTPCCC) based upon the fundamental work of Michelsen (1982a, 1982b). Generally, our use of stability proposed test (MTPCCC) to direct computation of vapor liquid equilibrium in the capillary condensation area without necessary to compute saturation pressure of flat and curved system.

### 2.3. The nano impact on the critical parameters of natural gas components

A very interesting discussion of the phase nanopore effects is presented in the paper of Devogowda et al. (2012). This paper discusses such change of critical parameters in micro-and meso-pores – so called “confinement effect”. The authors discuss the Monte Carlo simulations Hamach et al. (2007), Sing et al. (2009), Travalloni et al. (2012), Ortiz et al. (2005). Numerical study of Sing et al. (2009) experiment being used by the Devogowda et al. in the construction of function of the critical pressure and temperature as a function of molecular weight. Devogowda et al. (2012) shows the possible impact of the curvature coefficient of porous media to compressibility factor Z, phase PVT properties (without capillary effects, and without adsorption). This work discusses the influence of parameters on the viscosity of the pore gas as a function of the Knudsen number. The present work uses data of Sing et al. (2009) for correlating the critical temperature as a function of the effective pore radius (or diameter).

where T_{cc} — confined critical temperature, T_{c} — critical temperature, r_{0} — reference pore diameter (1.5 nm), r — actual pore diameter (nm). The critical pressure has been estimated on the basis of Sing et al. (2009) as a function similar to the correlation equation Antoione for vapor pressure:

where Pcc — confined critical pressure, Pc — critical pressure. Both equations are valid only for 1.5-5 nm radius (or 3-10 nm pore diameter).

A similar correlation functions for the critical pressure also suggested Didar & Akkutlu (2013). The relative shift of critical parameters (after Sing(2009) may be defined as:

The second correlation relevance to dimensionless shift of critical parameters pressure and temperature is presented in Ma et al. (2013), Ma & Jamili (2014) paper:

where ΔTc, ΔPc — dimensionless shift of critical temperature and pressure relevant to confinement phenomena.

Dp = 2 rp — porous diameter, σ – effective particle diameter based on the formula:

### 2.4. Phase envelope algorithm (classical bulk – hydrocarbon dew pressure case)

The algorithm of calculation of bulk phase envelope in gas condensate system is shown in figure 9 (Nagy, 1992). standard Tangent Criterion (Michelsen, 1982a, 1982b) is used in calcula- tion. The presented algorithm is not sensitive to selection of Equation of State. In this scheme two solution method are preferable, simple successive substitution (λ = 1) or GDEM (General Dominant Eigenvalue Method) where λ-acceleration exponent (Crowe & Nishio, 1975). The method has been tested in several papers (e.g. Nagy, 1992a, 1992b, 2002, 2003).

### 2.5. An algorithm for the calculation of the capillary condensation dew pressure

In the presence of transition region (gas-oil) in the tight/shale rock (where average capillary radius is below 1 μm – see Nagy 2002, 2003) we should include the additional terms related to adsorption and capillary condensation process. The proposed algorithm for calculation of the capillary condensation dew pressure consists of the following steps:

- Compute set of initial equilibrium constants K
^{∞}, - Compute thickness (t) of the adsorption film using any meso-pores adsorption model (i.e. Frankel-Halsey-Hill, simplified de Boer et al., etc.),
- Compute the interfacial tension (IFT) (σ) using scale or parachor model,
- Compute capillary pressure p
_{c}, - Compute the set of liquid partial volume of components V
_{iL}and Pointing factors (ε_{i}), - Compute new-corrected equilibrium constants using equation
- Compute the newly-corrected condensate capillary equilibrium composition of the liquid phase,
- Perform saturation pressure correction for capillary/bulk vapour/liquid equilibria using,
- Repeat steps 2-7 until convergence is reached.

## 3. Application of theory to possible reservoir gas-condensate systems

The selection of proper Equation of State for VLE is important. Although the is hot discussion related to new EOS designated to description of nano confinement phenomena (Zarragoicoechea & Kuz, 2004; Derouane, 2007; Travalloni et al., 2010, 2013) – this paper is focused on classical EOS. The computation procedure in this work has been done using volume translated version of Peng-Robinson equation. From among a lot of versions of PR CEOS (Peng-Robinson, 1976; Peneloux et al., 1982; Stryjek & Vera, 1986; Jhaveri & Youngren, 1988; Patel & Teja, 1982, Stamataki & Magoulas, 2001; Lin & Duan, 2005). The Tsai and Chen (1998) VTPR EOS version have been chosen. The mentioned Tsai Chen VTPR EOS and VTPR in the form of being in the opinion of the authors most versatile EOS for high pressure vapor liquid equilibria. The original Tsai and Chen (1998) form of the equation is given below:

where t — the translated volume parameter, with the volume calculated by the CEOS approaches the experimental value (sees Nagy, 2002; Tsai & Chen, 1998).

The most important parameter which is needed in this algorithm is proper evaluation of surface tension oil/condensate-gas in the high pressure VLE. The presented computation is related to new papers Fanchi (1990),

Danesh et al. (1991). The data from this paper has been used in computations. The models of Fanchi and Danesh et al. are related to critical parameters. The solution procedure is similar as presented by Nagy (1992, 2002) using classical critical properties taken from Reid, Prausnitz (1977) and using confined parameters (eq. 16-17 or 20-21). The confined critical parameters have been used in this calculation for case of 5 and 2 nm pore diameters. The mole composition of gas condensate mixture has been taken from Brusilovsky (1990): methane 0.7643; ethane 0.0746; propane 0.0312; i-butane 0.0059; n-butane 0.0121; i-pentane 0.005; n-pentane 0.0059; n-hexane 0.0079; n-heptane 0.01; n-octane 0.01; n-nonane 0.01; n-decane 0.037; CO_{2} 0.0249; N_{2} 0.0012).

Results of computation are given in Figure 8. The capillary-adsorption impact on the vapor phase fraction near flat upper part dew point line is small except near cricondentherm point – above 420 K. The cricondentherm pressure calculated with the set of equation (4-8) is located 30 bars above “flat surface” cricondentherm pressure. The largest impact of capillary effects is observed in the lower part of the curve, where the difference in increase of pressure is close to 40 bars. The capillary effects are much smaller in case of use of confined critical parameters, but are still important in the lower part of the phase envelope. The confined parameters (eq. (16-17)) are used for calculation of IFT. The effect of capillarity is very small for pore diameter 2 nm in the case of use of confined critical parameters. The impact of new confined critical parameters for calculation of phase equilibria of rocks with 2, 5 20+ nm pore diameter is large, as you may observe in the Figure 8.

In case of the reservoir temperature higher than 385°K, the process of condensation of the heavy component (C_{2}-C_{10}) is stopped for 2 nm pore diameter. In this case the saturation pressure of 5 nm pores is lower up to 150 bars – difference (80 bars). In such example the condensation process is expected in the pores of diameter 100-20 nm. This may be important in the deep shale gas structures, below 3000 m, where the shift of average pore diameters of structure may be valid.

Similar impact of nano-porosity on bubble and dew curve is presented on figure 9 for ternary oil-like composition (CH_{4} – 10%) – nC_{4}H_{10} (25%) – nC_{8}H_{18} (65%) using TC TVPR EOS. This calculation is done without adsorption and capillary pressure modification, but with the shift of critical parameters for confined conditions (eq. 20-21). In this situation we observe lowering of bubble presssure in confinement pores (Dp = 4 nm). An important problem relevant to calculation of the compressibility factor of single phase – like methane – has to be additional research (see Fig. 10). The observed isotherm at 85°C shows different shapes and values in case of confinement pore structure. This may influence significant on estimation Gas In Place Resources.

## 4. Conclusions

The porous media dew pressure of gas condensate in this paper shows mixed impact in case of capillary/adsorption effects. The impact of porous media in near critical condition is small. It is a low probability to expect of any capillarity effect in the high permeability rock. The new numerical experiments may be very useful in understanding of VLE phenomena, but should be verified in the future. The shift of phase envelope in the nano-porous structure should be still investigated. Because only part of the shale gas reservoir has a nanopore structure of 2-5nm diameter, probably the shift of phase envelope parameters must be addressed only to those elements of structure. The rest of the reservoir, which have a pore diameter larger than 10 nm should be described using classical tools with capillary correction as presented in this paper.

#### Acknowledgments

The authors are grateful to the National Centre for Research & Development for providing research fund through “Blue Gas Program” (projects “IRES”) and statutory research program (11.11.190.01/2014) performed at the AGH University of Science and Technology (Faculty of Drilling, Oil and Gas).

## References

- Adamson A.W., 1990. Physical Chemistry of Surfaces. J. Wiley & Sons Inc.
- Ahmed T., 1989. Hydrocarbon Phase Behavior (Contributions in Petroleum Geology and Engineering). Gulf Publishing Co.
- Altman R.M. et al., 2014. Understanding Mechanisms for Liquid Dropout from Horizontal Shale Gas Condensate Wells. Society of Petroleum Engineers. doi:10.2118/170983-MS.
- Ambrose R.J., 2011. Micro-Structure of Gas Shales and Its Effects on Gas Storage and Production Performance. Master’s Thesis, University of Oklahoma, Norman, OK.
- Brusilovsky A.I., 1990. Mathematical Simulation of Phase Behavior of Natural Multicomponent Systems at High Pressures with an Equation of State. SPE 20180.
- Campos M.D., 2010. Uncertainties in Shale Gas-_in-Place Calculations: Molecular Simulation Approach. Ph.D. Dis- sertation, University of Oklahoma, Norman, OK, 2011.
- Campos M.D., Akkutlu I.Y., Sigal R.F., 2009. A Molecular Dynamics Study on Natural Gas Solubility Enhancement in Water Confined to Small Pores, Paper SPE124491-MS.
- Cao Minh et al., 2012. 2D-NMR Applications in Unconventional Reservoirs. SPE 161578.
- Clarkson C.R., Haghshenas B., 2012. Modeling of Supercritical Fluid Adsorption on Organic-Rich Shales and Coal. SPE 164532.
- Clarkson C.R., Wood J.M., Burgis S.E., Aquino S.D., Freeman M., Birss V., 2012. Nanopore Structure Analysis and Permeability Predictions for a Tight Gas/Shale Reservoir Using Low-Pressure Adsorption and Mercury Intrusion Techniques. SPE 155537.
- Crowe C.M., Nishio M., 1975. Convergence Promotion in the Simulation of Chemical Process – the General Dominant Eigenvalue Method. AIChE J.,Vol. 21, No. 3, 528-533.
- Danesh A.S., Dandekar A.Y., Todal A.C., Sarkar R., 1991. A Modified Scaling Law andParachor Approach for Improved Prediction of Interfacial Tension of Gas-Condensate System. SPE 2270.
- Defay R., Prigogine I., 1966. Surface Tension and Adsorption. Longmans, London.
- Derouane E.G., 2007. On the physical state of molecules in microporous solids. Microporous and Mesoporous Materials. 104. 46-51.
- Devegowda D., Sapmanee K., Civan F., Sigal R.F., 2012. Phase Behavior of Gas Condensates in Shales Due to Pore Proximity Effects: Implications for Transport, Reserves and Well Productivity. Society of Petroleum Engineers. doi:10.2118/160099-MS.
- Dhanapal K. et al., 2014. Phase Behavior and Storage in Organic Shale Nanopores: Modeling of Multicomponent Hydrocarbons in Connected Pore Systems and Implications for Fluids-in-place Estimates in Shale Oil and Gas Reservoirs. SPE-169008-MS.
- Didar B.R., Akkutlu I.Y., 2013. Pore-Size Dependence of Fluid Phase Behavior and Properties in Organic-Rich Shale Reservoirs. SPE 164099.
- Elamin A. et al., 2013. Simulation of Multicomponent Gas Flow and Condensation in Marcellus Shale Reservoir. SPE 159869.
- Everett D.H., 1972. IUPAC manual of symbols and terminology for physicochemical quantities and units. App. II, Part I. Pure Appl. Chem. Vol. 21, p. 584-594.
- Fanchi J., 1990. Calculation of Parachors for Compositional Simulation: An Updated. SPE Res. Eng., p. 43 SPE 19453. Firincioglu T. et. al. 2012. Thermodynamics of Multiphase Flow in Unconventional Liquids rich Reservoirs. SPE 123455. Guo P. et al., 1996. A Theoretical Study of the Effect of Porous Media on the Dew Point Pressure of a Gas Condensate. SPE 25644.
- Hamada Y., Koga K., Tanaka H., 2007. Phase equilibria and interfacial tension of fluids confined in narrow pores. J. Chem. Phys., 127 (8): 084908-1-084908-9.
- Hartman R.C. et al., 2011. Shale Gas-in-Place Calculations Part 2 – Multicomponent Gas Adsorption Effects. SPE 144097-MS.
- Honarpour et al., 2012. Characterization of Critical Fluid, Rock, and Rock-Fluid Properties – Impact on Reservoir Performance of Liquid-Rich Shales. SPE 158042.
- IEA, 2012. World Energy Outlook: Golden Rules in the Golden Age of Natural Gas. Report IEA. IEA, 2011. World Energy Outlook 2011: Golden Age of Natural Gas? Report IEA.
- Jhaveri B.S., Youngren G.K., 1988. Three Parameter Modification of the Peng Robinson Equation of State to Improve Volumetric Predictions. SPE Res. Eng., (Aug. 1988), p. 1033.
- Kaliski M., Nagy S., Siemek J., Sikora A., Szurlej A., 2012. Unconventional natural gas-USA, the European Union, Poland. Archiwum Energetyki, 42, 109-122 (In Polish).
- Kang S., Fathi E., Ambrose R.J., Akkutlu I.Y., Sigal R.F., 2011. Carbon Dioxide Storage Capacity of Organicrich Shales. SPE Journal, Vol. 16 (4), 842-855.
- Kang S.M., 2011. Carbon Dioxide Storage Capacity of Barnett Shale. MSc, Oklahoma Univ.
- Klimkowski Ł., Nagy S., 2014. Key Factors in Shale Gas Modeling and Simulation. Arch. Min. Sci., Vol. 59, No 4, p. 987-1004.
- Lee S-T., 1999. Capillary-Gravity Equilibria for Hydrocarbon Fluids in Porous Media. SPE 19650.
- Lin H., Duan Y-Y., 2005. Empirical correction to the Peng-Robinson equation of state for the saturated region. Fluid Phase Equilibria, 233 (2005) 194-203.
- Ma Y., Jamili A., 2014. Modeling the Effects of Porous Media in Dry Gas and Liquid Rich Shale on Phase Behavior. Society of Petroleum Engineers. doi:10.2118/169128-MS.
- Ma Y., Jin L., Jamili A., 2013. Modifying van der Waals Equation of State to Consider Influence of Confinement on Phase Behavior. Society of Petroleum Engineers. doi:10.2118/166476-MS.
- Michelsen M.L., 1982a. The Isothermal Flash Problem I. Stability analysis. Fluid Phase Equilibria Vol. 8, p. 1-19. Michelsen M.L., 1982b. The Isothermal Flash Problem I. Phase Split Calculation. Fluid Phase Equilibria Vol. 8, 21-40. MIT, 2011. Future of Natural Gas, An Interdisciplinary MIT Study. Massachusetts Inst. of Tech., Boston.
- Nagy S., 1992. Isenthalpic throttling effect in multiphase and multicomponent systems. Archiwum Termodynamiki, Vol. 12, No. 1-4, p.116-128.
- Nagy S., 2002. Capillary adsorption effects in gas condensate systems in tight rocks. Arch. Min. Sci., Vol. 47, No 2, p. 205-253
- Nagy S., 2003. Capillary adsorption effects in gas condensate systems in tight rocks – Vertical variation of hydrocarbon composition. Arch. Min. Sci., Vol. 48, No 3, p. 355-402.
- Nagy S., Siemek J., 2011. Shale Gas in Europe: the state of the technology - challenges opportunities. Arch. Min. Sci., Vol. 56, No. 4.
- Nagy S., Siemek J., 2013. Phase envelope of gas condensate systems in nanopores. [In:] ICEE/ICIT-2013 conference: joint International Conference on Engineering Education and Research and International Conference Information Technology: 8-12 December Cape Town, South Africa. – Cape Town: Cape Peninsula University of Technology, p. 111-119.
- Orangi A. et al., 2011. Unconventional Shale Oil and Gas-Condensate Reservoir Production, Impact of Rock, Fluid and Hydraulic Fractures. SPE 140536.
- Ortiz V., Lopez-Álvarez Y.M., Lopez G.E., 2005. Phase diagrams and capillarity condensation of methane confined in single- and multi-layer nanotubes. Molecular Physics, 103 (19): 2587-2592.
- Pang J. et al., 2012. Impact of Porous Media on Saturation Pressures of Gas and Oil in Tight Reservoirs. SPE 161143. Patel N.C., Teja A.S., 1982. A new cubic equation of state for fluids and fluids mixtures. Chem. Eng. Sci., Vol. 37, p. 463-473.
- Pedersen K.S., Fredenslund A., Thomassen P., 1989. Properties of Oils and Natural Gases. Gulf Publishing Co., Houston, Tx.
- Peng D.-Y., Robinson D.B., 1976. A New Two-Constant Equation of State. Ind. Eng. Chem. Fund., 15, No. 1, 59-64. PGI, 2012. Assessment of shale gas and shale oil resources of the Lower Paleozoic Baltic-Podlasie-Lublin Basin in Poland. Report, Warsaw 2012.
- Ping S. et al., 1996. A Theoretical Study of the Effect of Porous Media on the Dew Point Pressure of a Gas Condensate. SPE 35644.
- Poprawa P., 2010. Poszukiwanie złóż gazu ziemnego w łupkach (shale gas) w Polsce. Wiadomości Naftowe I Gazownicze, 2/2010 (In Polish).
- Reid R.C., Prausnitz I.M., Sherwood T., 1977. The Properties of Gases and Liquids. third edition, McGraw-Hill Book Co. Inc., New York City.
- Rubin B., 2010. Accurate simulation of non Darcy flow in stimulated fractured shale reservoirs. SPE paper 132093. Sadyk K., Zade E.S., 1963. Determination of the Beginning of Condensation in the Presence of Porous Media. Neft i Gaz. Sadyk K., Zade E.S., 1968. A Study of the Process at Equilibrium Rate Attainment During Condensate System.
- Neft i Gaz. Shapiro A.A., Stanby E.H., 1996. Effect of Capillary Forces & Adsorption on Reserves Distributions. SPE 36922. Siemek J., Nagy S., 2012. Energy carriers use in the world: natural gas – conventional and unconventional gas resources, Arch. Min. Sci., Vol. 57, No 2, p. 283-312.
- Siemek J., Nagy S., Siemek P., 2013. Challenges for sustainable development: the case of shale gas exploitation in Poland. Problems of Sustainable Development, vol. 8, no. 1, p. 91-104.
- Siemek J., Rychlicki S., Kolenda Z., Stryczek S., Nagy S., 2009. Gaz ziemny i ropa naftowa–Polska i Unia Europejska, potrzeby, własne wydobycie i dostawy. Zarys stanu i perspektywy energetyki polskiej. Studium AGH, AGH Kraków (In Polish).
- Sigmund P.M. et al., 1973. Retrograde Condensation in Porous Media. SPE 3476.
- Singh S.K. at al., 2009. Vapor-Liquid Phase Coexistence, Critical Properties, and Surface Tension of Confined Alkanes. J. of Phys. Chem. C, 113(17).
- Smith L.R., Yarborough L., 1968. Equilibrium Revaporization of Retrograde Condensate by Dry Gas Injection. SPEJ (March 1968), p. 87-94.
- Stryjek R., Vera J.H., 1986. PRSV: an improved Peng-Robinson equation of state for pure compounds and mixtures. Can. J. Chem. Eng., Vol. 64, p. 323-333.
- Tindy R., Raynal M., , 1966. Are Test-Cell Saturation Pressures Accurate Enough? Oil and Gas Journal, 64(126). Travalloni L., Castierb M., Tavaresa F.W., Sandler S.I., 2010. Critical behavior of pure confined fluids from an extension of the van der Waals equation of state. Journal of Supercritical Fluids, 55 (2): 455-461.
- Trebin F.A., Zadora, 1968. Experimental Study of the Effect of a Porous Media on Phase Change in Gas Condensate System. Neft’ i Gaz, 8(37).
- Warowny W., Nagy S., Śniatała K., Roszak R., Szołkowski S., 2004. Wyznaczanie punktu rosy węglowodorów w gazie ziemnym. Gaz, Woda i Technika Sanitarna, 114-122 [In Polish].
- Whitson C., Sunjerga S., 2012. PVT in Liquid-Rich Shale Reservoirs. SPE 155499. Whitson C.H., Brule M R., 2000. Phase Behavior. SPE Monograph Series. Vol. 20.
- Zhang Y., Civan F., Devegowda D., Jamili A., Sigal R. F., 2013. Critical Evaluation of Equations of State for Multi- component Hydrocarbon Fluids in Organic Rich Shale Reservoirs. Society of Petroleum Engineers. doi:10.1190/ URTEC2013-182.
- Zuo Y., Stenby E.H., 1997. Corresponding-States and Parachor Models for the Calculation of Interfacial Tensions. Can. J. Chem. Eng., 75, 1130-1137.

Received: 22 April 2014

#### Acronyms

PVT – Pressure-Volume-Temperature,

VLE – Vapor Liquid Equilibria,

PR CEOS – Peng-Robinson Cubic Equation of State,

VTPR EOS – Volume Translated Cubic Equation of State,

CEOS – Cubic Equation of State,

IFT – Interfacial Tension.

#### Nomenclature

A, B, C – coefficients in Peng-Robinson EOS,

a, b, t – coefficients in Peng-Robinson EOS,

f – fugacity,

k – binary interaction constant,

n – number of components,

p – pressure,

T – temperature,

R – universal gas constant,

x – mole fraction of component in liquid phase,

y – mole fraction of component in vapor phase,

x – mole fraction of component in liquid or vapor phase,

v – molar volume,

w – Pitzer’s acentric factor,

z – mole fraction of component in system,

Ŕ – fugacity coefficient,

ß – phase mole fraction,

#### Indexes

c – critical parameter

i – i-component,

j – j-component

l – liquid phase

n – standard condition

v – vapor phase.

##### Authors: STANISLAW NAGY*, JAKUB SIEMEK*

*AGH UNIVERSITY OF SCIENCE & TECHNOLOGY, DRILLING & OIL-GAS FACULTY, AL. A. MICKIEWICZA 30, 30-059 KRAKOW, POLAND. E-mail: stanislaw.nagy@agh.edu.pl, siemek@agh.edu.pl*