Experiments and Modeling of Volumetric Properties and Phase

May 8, 2012 - In the Tarim region of. China, ultra-high-pressure gas reservoirs/condensate gas reservoirs were also discovered with buried at great de...
0 downloads 0 Views 1MB Size
Article pubs.acs.org/IECR

Experiments and Modeling of Volumetric Properties and Phase Behavior for Condensate Gas under Ultra-High-Pressure Conditions Chang-Yu Sun,* Huang Liu, Ke-Le Yan, Qing-Lan Ma, Bei Liu, and Guang-Jin Chen* State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249, China

Xiang-Jiao Xiao and Hai-Ying Wang Research Institute of Exploration and Development, PetroChina Tarim Oilfield Company, Korla 841000, China

Xi-Tan Zheng and Shi Li Research Institute of Petroleum Exploration and Development, PetroChina, Beijing 100083, China S Supporting Information *

ABSTRACT: Four reservoir samples under ultra-high-pressure and high-temperature conditions were collected from condensate gas fields in China. Constant-composition expansion tests were performed to determine the phase behavior and volumetric properties of reservoir fluid using an ultra-high-pressure fluid PVT test system. The compressibility factor and dewpoint pressure were obtained at four temperatures for four samples. The range of pressure was from 22.03 to 118.89 MPa. For the samples studied, the experimental results showed that the dew-point pressure decreased with increasing temperature and the compressibility factors increased with increasing pressure but decreased with increasing temperature at a given high reduced pressure. A thermodynamic model based on an equation of state was developed to describe the volumetric properties and phase behavior of the condensate gas under ultra-high-pressure conditions. The calculated results are in good accordance with the experimental data, which is important for the development of condensate gas reservoirs in ultra-high-pressure environments.

1. INTRODUCTION Although most producing oil and gas fields exhibit pressures lower than 50 MPa and temperatures lower than 423 K, oil and gas fields have been found under more severe conditions in several regions as a result of improvements in drilling techniques.1 These kinds of complex reservoirs are called ultra-high-pressure gas reservoirs with high pressures (generally higher than 100 MPa) and high temperatures. For instance, several condensate gas reservoirs at temperatures as high as 463 K and pressures above 110 MPa have been found in the North Sea.1 Ultra-high-pressure gas reservoirs were also discovered in the Gulf of Mexico, Japan, and Australia. In the Tarim region of China, ultra-high-pressure gas reservoirs/condensate gas reservoirs were also discovered with buried at great depths with high reservoir pressures and temperatures. Condensate gas is considered to be a very complex reservoir fluid because of its high gas/oil ratio. Knowledge of the phase behavior of condensate gas systems is important for predicting reservoir performance and future processing needs.2 Two of the most important factors to be evaluated by engineers in planning the development of a condensate gas reservoir are the fluid compressibility factor, which determines its volumetric behavior, and the dew-point pressure of the original fluid, which determines the retrograde dew-point curve below which heavier components form a liquid phase.3−6 For reservoirs under high pressure and high temperature, it becomes a very interesting challenge not only to sample fluids and perform © 2012 American Chemical Society

PVT experiments, but also to predict the phase behavior using typical thermodynamic models.5 Three types of methods are used to determine dew-point pressures and volumetric behavior: experimental determination, empirical correlations (with some accuracy limitations), and thermodynamic models mostly based on equations of state. As examples of experimental studies, Ungerer et al.1 studied the phase behavior of four synthetic gas condensates at 298−426 K and pressures up to 60 MPa. Guo et al.7 measured the compressibility factors of three natural gas samples in the ranges of 322.5−413.4 K and 20.1−95.9 MPa. Gozalpour et al.8 measured the densities and viscosities of two North Sea-based gas condensates at temperature up to 423 K and pressures up to 140 MPa. They also measured the vapor−liquid equilibrium data of condensate gas at five temperatures in the range of 278.2−383.2 K and pressures of 10.45−41.47 MPa.9 Only a small number of experiences with reservoirs where the pressure and temperature of formation can reach up 110 MPa have been reported in the open literature.10,11 As examples of empirical correlations, Standing and Katz12 developed a chart according to the theory of corresponding states for the compressibility factor. Since their work, more than 20 correlations have become available for calculating the Received: Revised: Accepted: Published: 6916

November 9, 2011 April 24, 2012 April 27, 2012 May 8, 2012 dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

Figure 1. Schematic diagram of the ultra-high-pressure fluid PVT system.

compressibility factor, including those of Papay,13 Hall and Yarborogh,14 Brill and Beggs,15 Dranchuk et al.,16 Dranchuk and Abou-Kassem,17 Porshakov et al.,18 and Heidaryan et al.19 Many researchers, such as Kurata and Katz,20 Eilerts and Smith,21 Olds et al.,22,23 Reamer and Sage,24 Nemeth and Kennedy,25 Carlson and Cawston,26 Potsch and Braeuer,27 Marruffo et al.,28 Humoud and Al-Marhoun,29 and Elsharkawy,3 have correlated the dew-point pressures of condensate gas reservoirs with gas/oil ratio, temperature, composition, boiling point, oil gravity, or characteristics of the heptanes-plus (C7+) fraction. Such correlations are limited to the complexity and range of condensate gas parameters used. To overcome the limitations of empirical correlations, thermodynamic models such as equations of state are used extensively in the petroleum industry to study the volumetric and phase behavior of petroleum reservoir fluids. However, such an approach can deteriorate for phase-behavior modeling of complex hydrocarbons, such as volatile oils and condensate gas, especially in the retrograde region.6,30 Choosing a suitable equation of state (EOS)5,31−36 is important. In addition, it is necessary to fit some parameters to predict the thermodynamic behavior of condensate gas for thermodynamic models.37 The following tuning methods have been employed to improve EOS capabilities for predicting the phase behavior of condensate gas: splitting and characterization of the C7+ fraction,38−44 adjustment of the methane interaction coefficient,45−47 adjustment or regression of the EOS parameters,48 adjustment of the critical properties,39,49 adjustment of the molecular weight of the heavy fraction,50 and volume translation.51−53 As an alternative to EOS methods, Elsharkawy and Foda2 developed a neural network model capable of predicting the dew-point pressure and simulating the depletion performance of condensate gas reservoirs. Nowroozi et al.54 proposed a fuzzy neural system to estimate dew-point pressures. Shokir6 applied genetic programming with an orthogonal least-squares algorithm to generate a dew-point-pressure prediction model represented by tree structures. In addition, Monte Carlo simulations have also been used to predict the thermodynamic properties of condensate gases.4,55,56 Although many correlations and models for describing the phase behavior of condensate gas exist, their suitable application range generally requires a reduced pressure of less than 15, as high errors occur when they are applied at reduced pressures higher than 15 in view that the lack of corresponding experimental data under ultra-high-pressure conditions for this type of condensate gas reservoirs. Therefore, in this work, the compressibility factors and dew-point pressures of four condensate gas samples at high pressure and temperature were measured. A thermodynamic model based on an EOS was

developed to describe the volumetric properties and phase behavior of condensate gas under ultra-high-pressure conditions. The obtained results could be of importance for the development of condensate gas reservoirs under ultra-highpressure environments.

2. EXPERIMENTAL SECTION 2.1. Sample Preparation. Four groups of condensate gas fluids from China gas reservoirs were sampled using the separator sampling method. The separator gas composition was analyzed with a HP6890 gas chromatograph. The gas chromatograph was calibrated by injecting standard gases in advance, and its minimum detection limit was found to be 200 ppmv. Each gas sample was measured at least three times, and the average value was used as the composition of the gas mixture. The separator oil was flashed, the tank gas was analyzed with the HP6890 instrument, and the tank oil was analyzed with an HP5890A gas chromatograph. The compositions of the separator oil and in situ reservoir fluids were then obtained by recombining the oil and gas according to the gas/oil ratio. The compositions of the separator gas, separator oil, and reservoir fluid; the gas/oil ratio; and the molecular weight and relative densities of the C11+ fraction for four groups of condensate gas samples are listed in Table S1 in the Supporting Information. The reservoir fluids were then prepared in the laboratory by recombining the separator oil with the separator gas according to the gas/oil ratio data based on a standard analytical method for condensate gas properties (Chinese standard SY/T 55432002). An ultra-high-pressure fluid PVT system was used for preparation of the reservoir fluid. To check the reliability of the prepared reservoir fluid samples, a sample was pressed to a single phase at the corresponding reservoir temperature, and a flash test was performed. Only when the compositions of the fluid recombined from the flashed gas and flashed oil were in agreement with those of the reservoir fluid could the sample prepared be used for the subsequent phase behavior measurements. 2.2. Constant-Composition Expansion Test. An ultrahigh-pressure fluid PVT system manufactured by Sanchez Technologies Co., Viarmes, France, was used to carry out phase behavior measurements on the condensate gas in this work. A schematic diagram of the experimental device is shown in Figure 1. The maximum operating pressure was up to 150 MPa, and the operating temperature was within the range from room temperature to 473 K. The accuracies of pressure and temperature measurements were ±0.01 MPa and ±0.1 K, respectively. The maximum working volume of the cell was 240 mL. 6917

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

The experimental measurements were processed by loading reservoir fluid into the condensate gas PVT cell. A constantcomposition expansion (CCE) test was performed to determine the compressibility factor, dew-point pressure, and liquid dropout behavior of the fluid. Thereafter, the dew point was determined from a dew-point detection system. The pressure was dropped in steps, and the liquid volume was recorded at each pressure stage after the experimental pressure had decreased below the dew-point pressure. For each group of sample, volumetric properties, liquid dropout curves at constant mass, and dew-point pressures at the reservoir temperature and three other temperatures below reservoir conditions were measured. 2.3. Experimental Results. In this work, similar experimental phenomena for the CCE test were observed through the windows of the PVT device for condensate gas fluid sampled from China gas fields. Sample 1, at a reservoir temperature of 407.2 K, is described here as an example to introduce the phenomena observed during the CCE test. The prepared sample was pressurized to 106.30 MPa at a reservoir temperature of 407.2 K. The reservoir fluid was a single phase that was transparent and pale yellow in color. As the CCE test proceeded and the pressure dropped, the color of the reservoir fluid also changed gradually from pale yellow to pale red. When the pressure was decreased to below 45 MPa, which is close to the dew point under reservoir conditions, the fluid began to fog up, and the fog thickened with decreasing pressure. The convective phenomena of the thick fog could also be observed by the naked eye, because of differences in the density of the fluid at different heights of the PVT cell and the adjustment of the fluid to maintain an equilibrium state. With a continuous drop in pressure, a tiny liquid drop precipitated, and the corresponding pressure was then assumed to be the dewpoint pressure of the reservoir fluid. In this group of experiments, the dew-point pressure was determined to be 40.2 MPa at a reservoir temperature of 407.2 K after repetitive measurements. Below the dew-point pressure, many tiny anomalous particles were observed to fall in the PVT cell with a further decrease of the pressure, whereas some of them adhered to the cell wall and others accumulated on the bottom of the cell and formed an unflowing liquid interface. The amount and size of the particles and the small liquid drop adhered to the cell wall also increased with decreasing pressure and then slid down and fell to the bottom of the cell after they attained a certain size. These phenomena continued until a new gas−liquid equilibrium was reached at the current pressure. For sample 1, CCE tests were performed at 313.2, 343.2, 373.2, and 407.2 K to determine the phase behavior at different temperatures from the stratum zone to the wellhead. Experimental phenomena similar to that occurring at 407.2 K were observed at the three other temperatures. Similarly, the constant-composition expansion tests for the other three samples, namely, samples 2−4, were also performed at four temperatures. The compressibility factors at different pressures above the dew-point pressure at four temperatures for samples 1−4 are listed in Tables S2−S5, respectively, of theSupporting Information. The four compressibility charts for the four samples at different reduced temperatures and reduced pressures are presented in Figures 2−5, where the calculation of the critical properties of the condensate gas will be introduced in a later section of this article. It was found that, in the single-phase zone, the compressibility factor increased with increasing pressure. Higher compressibility factor values

Figure 2. Compressibility chart at different reduced temperatures and reduced pressures for sample 1.

Figure 3. Compressibility chart at different reduced temperatures and reduced pressures for sample 2.

Figure 4. Compressibility chart at different reduced temperatures and reduced pressures for sample 3.

were obtained at lower temperatures for a given high reduced pressure. The variations of the retrograded liquid amount with pressure at four temperatures for samples 1−3 are shown in Figures 6−8, respectively. It was found that, in the initial stage, the retrograded liquid amount was not obvious. It increased with decreasing pressure and reached a maximum at the maximum retrograded pressure. After that, the retrograded liquid amount decreased with a continuous decrease of the pressure. The maximum retrograded liquid amounts for sample 1 at 313.2, 343.2, 373.2, and 407.2 K were 1.79%, 0.99%, 0.69%, and 0.52% at the corresponding maximum retrograded pressures of 17.47, 17.33, 16.07, and 17.13 MPa, respectively. 6918

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

Figure 8. Variation of the retrograded liquid percent of sample 3 with pressure during constant-composition expansion experiments at four different temperatures.

Figure 5. Compressibility chart at different reduced temperatures and reduced pressures for sample 4.

maximum retrograded liquid amount at stratum conditions was only 0.033% at a pressure of 17.04 MPa. The retrograded liquid amount for sample 4, with a gas/oil ratio of 112193 m3/m3, was very small and could not even be measured accurately. Table 1 lists the dew-point pressures at four temperatures for the four groups of condensate gas samples. It can be seen that, Table 1. Experimental Dew-Point Pressures (DPPs) for Samples 1−4 at Different Temperatures T (K)

DPP (MPa)

T (K)

sample 1 313.2 343.2 373.2 407.2

Figure 6. Variation of the retrograded liquid percent of sample 1 with pressure during constant-composition expansion experiments at four different temperatures.

46.13 42.09 41.02 40.20

313.2 343.2 373.2 411.1

29.95 29.20 27.68 24.65

323.2 353.2 383.2 418.6

sample 3 303.2 333.2 363.2 397.4

DPP (MPa) sample 2 46.55 43.84 42.18 40.99 sample 4 28.88 24.07 23.44 22.03

for each group of samples, the dew-point pressure decreased with increasing temperature. For samples 3 and 4, which had higher gas/oil ratios, lower dew-point pressure values were obtained than for samples 1 and 2.

3. THERMODYNAMIC MODEL Phase equilibrium and volumetric properties of condensate gas in the pressure and temperature range encountered during production operations are expected to be well-predicted. For this purpose, cubic equations of state are generally used to describe vapor−liquid equilibria of these fluids. However, these equations of states do not generally lead to good predictions under all PVT conditions.57 In particular, for condensate gas fluids, the calculated PVT properties do not match the experimental data.53,58 As reviewed in the Introduction, many tuning methods have been applied to improve EOS capabilities for describing the thermodynamic behavior of condensate gas. In this work, a thermodynamic model was used to describe the vapor−liquid equilibria and PVT properties for condensate gas under high-pressure and high-temperature conditions by introducing suitable methods for plus-fraction splitting and characterization of pseudocomponents. A transition function wais also introduced to improve the description of the

Figure 7. Variation of the retrograded liquid percent of sample 2 with pressure during constant-composition expansion experiments at four different temperatures.

The maximum retrograded liquid amount increased with decreasing temperature, meaning that more retrograded liquid will precipitate during condensate gas exploitation from the bottom of the well to the ground. A comparison of Figures 6−8 shows that the retrograded liquid amount decreased to a large extent with increasing gas/oil ratio. For sample 2, with a gas/oil ratio of 10683 m3/m3, the maximum retrograded liquid amount at stratum conditions was 0.48% at a pressure of 18.08 MPa. However, for sample 3, with a gas/oil ratio of 96424 m3/m3, the 6919

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

volumetric properties of condensate gas throughout the entire pressure range. 3.1. Development of the Model. A suitable equation of state is important for describing the phase behavior of condensate gas. In this work, the Soave−Redlich−Kwong (SRK) EOS31 was used, which is given by p=

a(T ) RT − v−b v(v + b)

a(T ) = 0.42747

b = 0.08664

R2Tc 2 α(Tr) Pc

RTc Pc

(1)

(2)

(3)

Tr )]2

m = 0.480 + 1.574ω − 0.176ω 2

(4) (5)

van der Waals-type single-fluid mixing rules for EOS parameters a and b were used a=

∑ ∑ xixj(aiaj)0.5 (1 − kij) i

j

(6)

n

b=

∑ xibi i=1

(7)

(13)

Pr 2 Tr

(14)

where β(T r ,P r ) is the transition function; Z 0 is the compressibility factor before the introduction of the transition function; and C1, C2, and C3, are the transition function parameters, which are 0.94134, 0.19822 × 10−2, and 0.34748 × 10−4, respectively, determined from the measured PVT data. The critical temperature and critical pressure of condensate gas were calculated by the following expressions according to the critical properties of the definite components and pseudocomponents

(8)

(9)

where, in this case, j refers to components other than CO2 and CH4 in the condensate gas sample. For ωj values larger than 0.35, the values of kCO2−j were set equal to 0.105. The methods of plus-fraction splitting and pseudocomponent characterization are also important for describing the phase behavior of condensate gas. In this work, because the C11+ component makes up only a small fraction of the total mixture, the extended analysis from C7 to C11+ was first lumped into a new pseudocomponent (i.e., C7+). Afterward, a suitable plus-fraction splitting method for the C7+ fraction was chosen to improve the capacity of the EOS under high-pressure and high-temperature conditions. Here, the following formula proposed by Ahmed et al.59 was used to split the C7+ fraction ⎛ MW(n + 1) + − MWn + ⎞ ⎟⎟ zn = zn +⎜⎜ ⎝ MW(n + 1) + − MWn ⎠

(12)

β(Tr , Pr) = C1 + C2Pr + C3

where j refers to components other than methane in the condensate gas sample. For kij between CO2 and other components (except CH4), the following correlation was used k CO2 − j = 0.12443 + 0.083828ωj − 0.093302ωj 2

MWn + = MWC7+ + S(n − 7)

Z = Z0β(Tr , Pr)

The binary interaction parameters between components i and j, kij, were considered to be equal to zero for hydrocarbon− hydrocarbon interactions, except for those between methane or CO2 and other components. For interactions between CH4 and other components, the interaction parameters were calculated by the following correlations k CH4 − j = 0.0575(ωj 0.5 − ωCH4 0.5)

(11)

where n is the single carbon number, which is greater than or equal to 7. For n = 7 or 8, S is 15.5; for n > 8, S is 17.0. For a given plus fraction, constraint equations were used to ensure that the distributions of mole fractions, molecular weights, and densities were consistent with the measured data for the plus fraction. The number of components and C7+ fractions in the final distribution might be too large. To reduce the overall number of components, a common procedure of lumping several components into a pseudocomponent, as described by Pedersen et al.,60 such that the weight fractions of all groups are approximately equal was used. In this work, 10 hydrocarbon groups were used to represent the hydrocarbon part of a reservoir fluid. The method proposed by Kesler and Lee61 was used to evaluate the critical properties and acentric factors for pseudocomponents. It is known that the main role of an EOS is to reproduce the PVT data accurately, whereas the correct vapor−liquid equilibrium (VLE) description can be achieved by devising an integration path.62 The calculations of PVT properties and VLE behavior using an EOS are relatively independent.63 Volume translation methods51−53 have been used to improve the description of the volume of the liquid phase by the equation of state. Similarly, in this work, a transition function was introduced to improve the description of the compressibility factor throughout the entire pressure range, without modifying the description of vapor−liquid equilibrium, that is

The form of α(Tr) in the original SRK EOS was employed for the whole Tr range α(Tr) = [1 + m(1 −

zn + = z(n − 1) + − z(n − 1)

Pc =

∑ xiPci

Tc =

∑ xiTci

(15)

n i=1

(16)

3.2. Comparison with Different α Functions in the Attractive Term. In view that the modification of the attraction term of an equation of state will favor the description of the phase behavior of condensate gas, the SRK EOS31 and two other models proposed by Twu et al.64 and Souahi et al.65 that are based on the modification of the α function of the attractive term of the Redlich−Kwong equation of state66 were used to examine the influence of the EOS form. Table 2 lists the corresponding results of calculating the dew-point pressure for sample 1. The absolute average deviations (AADs) for four temperatures were 11.02%, 4.44%, and 4.06% for the SRK EOS and the modifications by Twu et al.64 and Souahi et al.,65 respectively. The plus-fraction splitting method and the

(10)

where 6920

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

calculating the dew-point pressure for sample 2 are listed in Table 4. The absolute average deviations at four temperatures were 2.44%, 16.86%, 17.15%, 14.76%, 19.19%, and 18.33%, respectively, for the six combination methods. It was found that the model developed in this work, in which the SRK EOS was combined with the plus-fraction splitting method proposed by Ahmed et al.59 and the pseudocomponent characterization method proposed by Kesler and Lee,61 provided the best results for describing the phase behavior of condensate gas under highpressure and high-temperature conditions. 3.4. Prediction of Dew-Point Pressure for Four Samples. Table 5 lists the dew-point pressures predicted

Table 2. Comparison of Dew-Point-Pressure Calculation Results for Sample 1 with Different α Functions in the Attractive Terms of Three Equations of State RDa (%) T (K)

DPP (MPa)

SRK

Twu et al.64

Souahi et al.65

46.13 42.09 41.02 40.20

−11.67 −6.34 −9.32 −16.74 11.02

0.65 3.45 −2.01 −11.68 4.44

0.96 7.38 4.62 −3.29 4.06

313.2 343.2 373.2 407.2 AADa (%) a

RD = relative deviation, AAD = absolute average deviation. The same definitions of RD and AAD apply to Tables 4−6

Table 5. Comparison of Experimental and Calculated DewPoint Pressures

pseudocomponent characterization method used by these three models are the same and can be referred to as the methods of Whitson38 and Kesler and Lee,61 respectively. The binary interaction parameters kij between CH4 and other components and between CO2 and other components were set according to eqs 8 and 9, respectively. Other interactions between hydrocarbons were set to zero. It was found that the modification of the α function of the attractive term improved the description of the phase behavior of condensate gas under high-pressure and high-temperature conditions. Therefore, in the subsequent discussion, the results calculated with the models of Twu et al.64 and Souahi et al.65 are also listed for comparison. 3.3. Comparison with Different Characterization Methods. To examine the influence of the plus-fraction splitting and pseudocomponent characterization methods on describing the phase behavior of condensate gas, the three plusfraction splitting methods proposed by Ahmed et al.,59 Pederson et al.,67 and Whitson38 and the two characterization methods proposed by Kesler and Lee61 and Pederson et al.67 were investigated with the SRK EOS. The combinations of methods are listed in Table 3. The corresponding results of

RD (%) T (K)

Table 3. Methods of Plus-Fraction Splitting and Characterizationa

a

overall method

method for plus-fraction splitting

method for pseudocomponent characterization

this work method 2 method 3 method 4 method 5 method 6

Ahmed et al.59 Ahmed et al.59 Pederson et al.67 Pederson et al.67 Whitson38 Whitson38

Kesler and Lee61 Pederson et al.67 Kesler and Lee61 Pederson et al.67 Kesler and Lee61 Pederson et al.67

DPP (MPa)

313.2 343.2 373.2 407.2 AAD (%)

46.13 42.09 41.02 40.20

313.2 343.2 373.2 411.1 AAD (%)

46.55 43.84 42.18 40.99

303.2 333.2 363.2 397.4 AAD (%)

29.95 29.20 27.68 24.65

323.2 353.2 383.2 418.6 AAD (%)

28.88 24.07 23.44 22.03

this work sample −2.94 3.95 2.09 −3.88 3.22 sample 0.01 3.46 2.77 −3.54 2.44 sample 1.87 −1.90 −8.30 −14.66 6.68 sample 3.59 14.43 1.59 −20.77 10.10

AAD (%) for four samples

5.61

Twu et al.64

Souahi et al.65

0.65 3.45 −2.00 −11.68 4.45

0.96 7.38 4.62 −3.29 4.06

2.99 2.85 −0.78 −10.48 4.27

3.33 6.95 6.21 −1.19 4.42

13.47 6.27 −2.67 −16.39 9.69

12.92 8.87 2.20 −9.90 8.47

0.049 3.96 10.96 −50.07 16.26

1.43 7.94 −11.09 −50.07 17.63

1

2

3

4

8.67

8.65

using the model developed in this work for four samples at different temperatures. For comparison, the results calculated using two other models (those of Twu et al.64 and Souahi et al.65) based on the modification of the α function of the

Used with the SRK EOS.

Table 4. Comparison of Dew-Point-Pressure Calculation Results for Sample 2 Using Different Methods of Plus-Fraction Splitting and Characterization RD (%) T (K) 313.2 343.2 373.2 411.1 AAD (%)

DPP (MPa)

this work

method 2

method 3

method 4

method 5

method 6

46.55 43.84 42.18 40.99

0.01 3.46 2.77 −3.54 2.44

30.40 21.40 11.40 −4.25 16.86

−16.85 −13.82 −15.18 −22.73 17.15

−6.45 −9.11 −15.20 −28.27 14.76

−18.50 −15.67 −17.29 −25.31 19.19

−10.82 −12.72 −18.42 −31.39 18.33

6921

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

three models for sample 1 at four temperatures are shown in Figures 9−12, which clearly demonstrate the beneficial effect of

attractive term of the Redlich−Kwong equation of state are also included because of their high accuracy. As shown in Table 5, for sample 1, the predicted precisions were 3.22%, 4.45%, and 4.06% for this work and the models of Twu et al.64 and Souahi et al.,65 respectively. For sample 2, the corresponding predicted precisions were 2.44%, 4.27%, and 4.42%. It was thus found that the model developed in this work can be successfully used to describe the phase equilibrium of condensate gas at the low gas/oil ratio of about 10000 m3/m3. However, the errors increased with increasing gas/oil ratio. The predicted precisions for sample 3 were 6.68%, 9.69%, and 8.47% for this work and the models of Twu et al.64 and Souahi et al.,65 respectively. The corresponding predicted precisions for sample 4 were 10.10%, 16.26%, and 17.63%. For samples 3 and 4, the gas/oil ratio was about 100000 m3/m3, or about 10 times that for samples 1 and 2. The concentration of C11+ in condensate gas is low, and the measurement deviation for the C11+ composition is relatively high. Because the measured value of the dew-point pressure is sensitive to the C11+ content, the experimental data at higher gas/oil ratios have higher experimental deviations, which results in higher deviations of dew-point pressure predictions. Compared with those of the Twu et al.64 and Souahi et al.65 models, the predictive capabilities of the model developed in this work were always higher for the four groups of samples investigated. 3.5. Prediction of Compressibility Factor for Four Samples. Table 6 lists the compressibility factor values for four

Figure 9. Experimental and calculated compressibility factors of sample 1 at 313.2 K.

Table 6. Comparison of Experimental and Calculated Compressibility Factors AAD (%) T (K)

this work

Twu et al.64

Souahi et al.65

sample 1 313.2 343.2 373.2 407.2

2.07 1.59 1.06 0.69

3.20 3.63 4.18 4.62

3.69 4.58 5.01 5.54

313.2 343.2 373.2 411.1

0.54 1.53 1.76 1.71

2.85 4.05 4.37 4.32

2.83 4.47 5.11 5.74

303.2 333.2 363.2 397.4

0.89 1.16 1.60 1.95

3.38 3.45 3.71 3.82

3.26 3.63 4.22 4.72

323.2 353.2 383.2 418.6

2.05 1.54 2.15 1.15

2.28 1.93 1.28 1.84

2.75 2.23 2.72 3.51

total AAD (%)

1.46

3.31

4.00

Figure 10. Experimental and calculated compressibility factors of sample 1 at 343.2 K.

sample 2

sample 3

sample 4

Figure 11. Experimental and calculated compressibility factors of sample 1 at 373.2 K.

introducing the transition function for describing the compressibility factor. In addition, it can be seen that similar precisions for the four samples were obtained for the description of the compressibility factor, whereas the precision for the description of the dew-point pressure is related to the gas/oil ratio and the composition of plus fraction. It can also be seen that, compared with the calculation of the dew-point pressure, the calculation of the compressibility factor is significantly less sensitive to the characterization of the petroleum fraction.

samples at different temperatures predicted by the model developed in this work and by the models proposed by Twu et al.64 and Souahi et al.65 It was found that the total absolute average deviation was only about 1.46% using the model developed in this work. In comparison, the absolute average deviations were 3.31% and 4.00% using the models of Twu et al.64 and Souahi et al.,65 respectively. The comparisons of the 6922

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

equations of state. The model developed in this work favors the description of the phase behavior and volumetric properties of condensate gas under high-pressure and high-temperature conditions.



ASSOCIATED CONTENT

S Supporting Information *

Properties of reservoir fluid for four groups of condensate gas samples (Table S1) and compressibility factors at different pressures and four temperatures for four samples (Tables S2− S5). This material is available free of charge via the Internet at http://pubs.acs.org.



Figure 12. Experimental and calculated compressibility factors of sample 1 at 407.2 K.

AUTHOR INFORMATION

Corresponding Author

*Fax: +86 10 89732126. E-mail: [email protected] (C.-Y.S.), [email protected] (G.-J.C.).

In addition, the liquid drop-out amount was also predicted using the model developed in this work. Figure 13 shows the

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The financial support received from PetroChina Tarim Oilfield Company; the National Natural Science Foundation of China (Nos. 20925623, U1162205, 21076225); and the Science Foundation of China University of Petroleum, Beijing (LLYJ2011-63), is gratefully acknowledged.



REFERENCES

(1) Ungerer, P; Faissat, B.; Leibovici, C.; Zhou, H.; Behar, E.; Moracchini, G.; Courcy, J. P. High Pressure−High Temperature Reservoir Fluids: Investigation of Synthetic Condensate Gases Containing a Solid Hydrocarbon. Fluid Phase Equilib. 1995, 111, 287−311. (2) Elsharkawy, A. M.; Foda, S. G. EOS Simulation and GRNN Modeling of the Constant Volume Depletion Behavior of Gas Condensate Reservoirs. Energy Fuels 1988, 12, 353−364. (3) Elsharkawy, A. M. Predicting the Dew Point Pressure for Gas Condensate Reservoirs: Empirical Models and Equations of State. Fluid Phase Equilib. 2002, 193, 147−165. (4) Lagache, M. H.; Ungerer, P.; Boutin, A. Prediction of Thermodynamic Derivative Properties of Natural Condensate Gases at High Pressure by Monte Carlo Simulation. Fluid Phase Equilib. 2004, 220, 211−223. (5) Cañas-Marín, W. A.; Ortiz-Arango, J. D.; Guerrero-Aconcha, U. E.; Soto-Tavera, C. P. Thermodynamic Derivative Properties and Densities for Hyperbaric Gas Condensates: SRK Equation of State Predictions Versus Monte Carlo Data. Fluid Phase Equilib. 2007, 253, 147−154. (6) Shokir, E. M. El-M. Dewpoint Pressure Model for Gas Condensate Reservoirs Based on Genetic Programming. Energy Fuels 2008, 22, 3194−3200. (7) Guo, X. Q.; Wang, F.; Chen, G. J. Measurement of the Compressibility Factor of Natural Gases at Super High Pressure. J. Chem. Eng. Chin. Univ. 1999, 13, 393−397 (in Chinese). (8) Gozalpour, F.; Danesh, A.; Todd, A. C.; Tohidi, B. Viscosity and Density Data of Two North Sea Gas Condensate Samples at Temperatures to 423 K and Pressures to 140 MPa. J. Chem. Eng. Data 2001, 46, 1305−1308. (9) Gozalpour, F.; Danesh, A.; Todd, A. C.; Tehrani, D. H.; Tohidi, B. Vapour−Liquid Equilibrium Volume and Density Measurements of a Five-Component Gas Condensate at 278.15−383.15 K. Fluid Phase Equilib. 2003, 206, 95−104. (10) Neuhaus, D.; Browning, G.; Estensen, O. Managing Seismic Positioning Uncertainty in HPHT Developments. In Offshore Europe Oil and Gas Exhibition and Conference, 7−10 September 1999, Aberdeen,

Figure 13. Model-predicted results for the condensate liquid amount of sample 1.

results predicted for sample 1. It was found that the calculated retrograde liquid amount was approximatively in agreement with the experimental data, although a certain deviation still existed. One reason for the deviation is that the amount of retrograde liquid is very low for condensate gas at such high gas/oil ratios as used in this work, which could result in a relatively large measurement error.

4. CONCLUSIONS Volumetric properties, liquid drop-out curves at constant mass, and dew-point pressures of four groups of condensate gas fluid from China natural gas reservoirs at the reservoir temperature and three other temperatures below reservoir conditions were measured using an ultra-high-pressure fluid PVT system. A model based on the original SRK equation of state was presented in this study to describe the vapor−liquid equilibria and PVT properties of condensate gas under high-pressure and high-temperature conditions by introducing suitable plusfraction splitting and pseudocomponent characterization methods. A transition function was also introduced to improve the prediction of volumetric properties of condensate gas throughout the entire pressure range. The accuracies of the proposed model were verified by comparison with other 6923

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

United Kingdom; Society of Petroleum Engineers: Richardson, TX, 1999; Paper SPE 56950. (11) Fort, J. The Elgin/Franklin Project: Developing the Largest High Pressure/High Temperature Fields in the World. In Offshore Technology Conference 2000; Curran Associates, Inc.: Richardson, TX, 2000; Paper OTC 12117. (12) Standing, M. B.; Katz, D. L. Density of Natural Gases. Trans AIME 1942, 142, 140−149. (13) Papay, J. A Termelestechnologiai Parameterek Val-tozasa a Gastelepek Muvelese Soran. Ogil Musz, Tud, Kuzl., Budapest 1968, 267−273. (14) Hall, K. R.; Yarborough, L. A New Equation of State for ZFactor Calculations. Oil Gas J. 1973, No. June 18, 82−85. (15) Brill, J. P.; Beggs, H. D. Two Phase Flow in Pipes; University of Tulsa: Tulsa, OK, 1974. (16) Dranchuk, P. M.; Purvis, R. A.; Robinson, D. B. Computer Calculation of Natural Gas Compressibility Factors Using the Standing and Katz Correlation. In Institute of Petroleum Technical Series; Petroleum Society of Canada: Calgary, Alberta, Canada, 1974; Paper IP74-008. (17) Dranchuk, P. M.; Abou-Kassem, J. H. Calculation of Z Factors for Natural Gases Using Equations of State. J. Can. Petrol. Technol. 1975, No. July−September, 34−36. (18) Porshakov, Yu. P.; Kosachenko, A. N.; Nikishin, V. I. Power of Pipeline Gas Transportation; Oil and Gas: Moscow, Russia, 2001 (in Russian). (19) Heidaryan, E.; Salarabadi, A.; Moghadasi, J. A Novel Correlation Approach for Prediction of Natural Gas Compressibility Factor. J. Nat. Gas Chem. 2010, 19, 189−192. (20) Kurata, F.; Katz, D. L. Critical Properties of Volatile Hydrocarbon Mixtures. Trans. AIChE 1942, 38, 995−1021. (21) Eilerts, K.; Smith, R. V. Specific Volumes and Phase-Boundary Properties of Separator-Gas and Liquid-Hydrocarbon Mixtures; Report BM-RI-3642; Bureau of Mines: Bartlesville, OK, 1942. (22) Olds, R. H.; Sage, B. H.; Lacey, W. N. Volumetric and Phase Behavior of Oil and Gas from Paloma Field. Trans. AIME 1945, 160, 77−99. (23) Olds, R. H.; Sage, B. H.; Lacey, W. N. Volumetric and Viscosity Studies of Oil and Gas from a San Joaquin Valley Field. Trans. AIME 1949, 179, 287−302. (24) Reamer, H. H.; Sage, B. H. Volumetric Behavior of Oil and Gas from a Louisiana Field. Trans. AIME 1950, 189, 261−268. (25) Nemeth, L. K.; Kennedy, H. T. A Correlation of Dew Point Pressure with Fluid Composition and Temperature. SPE J. 1967, 7, 99−104. (26) Carlson, M. R.; Cawston, W. B. Obtaining PVT Data for Very Sour Retrograde Gas and Volatile Oil Reservoirs: A Muti-Disciplinary Approach. In SPE Gas Technology Symposium, 28 April−1 May 1996, Calgary, Alberta, Canada; Society of Petroleum Engineers: Richardson, TX, 1996; Paper SPE 35653. (27) Potsch, K. T.; Braeuer, L. A Novel Graphical Method for Determining Dew Point Pressure of Gas Condensates. In European Petroleum Conference, 22−24 October 1996, Milan, Italy; Society of Petroleum Engineers: Richardson, TX, 1996; Paper SPE 36919. (28) Marruffo, I.; Maita, J.; Him, J.; Rojas, G. Statistical Forecast Models to Determine Retrograde Dew Pressure and C7+ Percentage of Gas Condensates on Basis of Production Test Date of Eastern Venezuelan Reservoirs. In SPE Latin American and Caribbean Petroleum Engineering Conference, 25−28 March 2001, Buenos Aires, Argentina; Society of Petroleum Engineers: Richardson, TX, 2001; Paper SPE 69393. (29) Humoud, A. A.; Al-Marhoun, M. A. A New Correlation for Gas Condensate Dew Point Pressure Prediction. In Proceedings of Society of Petroleum Engineers (SPE) Middle East Oil Show Conference, Bahrain, March 17−20, 2001; Society of Petroleum Engineers: Richardson, TX, 2001; Paper SPE 68230. (30) Saker, R.; Danesh, A. S.; Todd, A. C. Phase Behavior Modeling of Gas Condensate Fluids Using an Equation of State. In Proceedings of Society of Petroleum Engineers (SPE) Annual Technical Conference and

Exhibition, Dallas, TX, Oct 6−9, 1991; Society of Petroleum Engineers: Richardson, TX, 1991; Paper SPE 22714. (31) Soave, G. Equilibrium Constants from Modified Redlich− Kwong Equation of State. Chem. Eng. Sci. 1972, 27, 1197−1203. (32) Peng, D. Y.; Robinson, D. B. A New Two Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59−64. (33) Martine, J. J. Cubic Equations of States: Which? Ind. Eng. Chem. Fundam. 1979, 18, 81−97. (34) Yang, T.; Chen, G. J.; Yan, W.; Guo, T. M. Extension of the Wong−Sandler Mixing Rule to the Three-Parameter Patel−Teja Equation of State: Application up to the Near-Critical Region. Chem. Eng. J. 1997, 67, 27−36. (35) Aasberg-Petersen, K.; Stenby, E. Prediction of Thermodynamic Properties of Oil and Gas Condensate Mixtures. Ind. Eng. Chem. Res. 1991, 30, 248−254. (36) Bonyadi, M.; Esmaeilzadeh, F. Prediction of Gas Condensate Properties by Esmaeilzadeh−Roshanfekr Equation of State. Fluid Phase Equilib. 2007, 260, 326−334. (37) Barreau, A.; Gaillard, K.; Béhar, E.; Daridon, J. L.; Lagourette, B.; Xans, P. Volumetric Properties, Isobaric Heat Capacity and Sound Velocity of Condensate Gas Measurements and Modeling. Fluid Phase Equilib. 1997, 127, 155−171. (38) Whitson, C. H. Characterizing Hydrocarbon-Plus Fractions. SPE J. 1983, No. August, 683−694. (39) Whitson, C. H. Effect of Physical Properties Estimation on Equation-of-State Predictions. SPE J. 1984, No. December, 685−696. (40) Whitson, C. H.; Anderson, T. F.; Soreide, I. C7+ Characterization of Related Equilibrium Fluids Using the Gamma Distribution. Presented at the AIChE Spring National Meeting, New Orleans, LA, Mar 6−10, 1988. (41) Barreau, A.; Braunschweig, B.; Emami, E.; Behar, E. A Knowledge-Based System for the Automation of Thermodynamic Models Adjustment. Presented at the Conference on Artificial Intelligence in Petroleum Exploration and Production, College Station, TX, May 15−17, 1991. (42) Braunschweig, B.; Barreau, A.; Emami, E. Genetic Algorithms for the Automatic Adjustment of Thermodynamic Models. In Artificial Intelligence in the Petroleum Industry; Technip: Paris, 1995. (43) Sportisse, M.; Barreau, A.; Ungerer, P. Modeling of Gas Condensates Properties Using Continuous Distribution Functions for the Characterization of the Heavy Fraction. Fluid Phase Equilib. 1997, 139, 255−276. (44) Jaubert, J. N.; Solimando, R.; Paricaud, P.; Barreau, A. Use of Distribution Functions: A Useful Tool to Calculate the Properties of Condensate Gases. Ind. Eng. Chem. Res. 2000, 39, 5029−5036. (45) Katz, D. L.; Firoozabadi, A. Predicting Phase Behavior of Condensate/Crude Oil Systems Using Methane Interaction Coefficients. J. Petrol. Technol. 1978, No. November, 1649−1655. (46) Gani, R.; Fredenslund, Aa. Prediction of Properties of Petroleum Mixtures: Improvement in Reliability and Validity of Tuned Parameters. Fluid Phase Equilib. 1986, 29, 575−582. (47) Gani, R.; Fredenslund, Aa. Thermodynamics of Petroleum Mixtures Containing Heavy Hydrocarbons: An Expert Tuning System. Ind. Eng. Chem. Res. 1987, 26, 1304−1312. (48) Coats, K. H.; Smart, G. T. Application of a Regression-Based EOS PVT Program to Laboratory Data. SPE Reservoir Eng. 1986, No. May, 277−299. (49) Stewart, G.; Stamataki, S.; Clancy, M.; Todd, A. C. Modelling of Gas Condensate Phase Behaviour Using Cubic Equations of State. North Sea Oil Gas Reservoirs 1987, 337−347. (50) Pedersen, K. S.; Thomassen, P.; Fredenslund, Aa. Characterization of Gas Condensate Mixtures. Presented at the AIChE Spring National Meeting, New Orleans, LA, Mar 6−10, 1988; Paper SEP 8804. (51) Peneloux, A.; Rauzy, E.; Freze, R. A Consistent Correction for Redlich−Kwong−Soave Volumes. Fluid Phase Equilib. 1982, 8, 7−23. (52) Jhaveri, B. S.; Youngren, G. K. Three-Parameter Modification of the Peng−Robinson Equation of State to Improve Volumetric Predictions. SPE Reservoir Eng. 1988, 3, 1033−1040. 6924

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925

Industrial & Engineering Chemistry Research

Article

(53) Esmaeilzadeh, F.; Samadi, F. Modification of Esmaeilzadeh− Roshanfekr Equation of State to Improve Volumetric Predictions of Gas Condensate Reservoir. Fluid Phase Equilib. 2008, 267, 113−118. (54) Nowroozi, S.; Ranjbar, M.; Hashemipour, H.; Schaffie, M. Development of a Neural Fuzzy System for Advanced Prediction of Dew Point Pressure in Gas Condensate Reservoirs. Fuel Process. Technol. 2009, 90, 452−457. (55) Neubauer, B.; Boutin, A.; Tavitian, B.; Fuchs, A. H. Gibbs Ensemble Simulations of VapourLiquid Phase Equilibria of Cyclic Alkanes. Mol. Phys. 1999, 97, 769−776. (56) Escobedo, F. A.; Chen, Z. Simulation of Isoenthalps and Joule− Thomson Inversion Curves of Pure Fluids and Mixtures. Mol. Simul. 2001, 26, 395−416. (57) Marteau, Ph.; Tobaly, P.; Ruffier-Meray, V.; Barreau, A. In Situ Determination of High Pressure Phase Diagrams of Methane-Heavy Hydrocarbon Mixtures Using an Infrared Absorption Method. Fluid Phase Equilib. 1996, 119, 213−230. (58) Anastasiades, K.; Stamataki, S.; Tassios, D. Prediction/ Correlation of Vapor−Liquid Equilibria in Synthetic Gas and Oil Systems. Fluid Phase Equilib. 1994, 93, 23−54. (59) Ahmed, T.; Cady, G.; Story, A. A Generalized Correlation for Characterizing the Hydrocarbon Heavy Fractions. In SPE Annual Technical Conference and Exhibition, 22−26 September 1985, Las Vegas, Nevada; Society of Petroleum Engineers: Richardson, TX, 1985; Paper SPE 14266. (60) Pedersen, K. S.; Thomassen, P.; Fredenslund, Aa. Thermodynamics of Petroleum Mixtures Containing Heavy Hydrocarbons III: Efficient Flash Calculations Using Cubic Equations of State. Ind Eng. Chem. Process Des. Dev. 1985, 24, 948−954. (61) Kesler, M. G.; Lee, B. I. Improve Prediction of Enthalpy of Fractions. Hydrocarbon Process. 1976, 55, 153−158. (62) Nitsche, J. M. New Applications of Kahl’s VLE Analysis to Engineering Phase Behavior Calculations. Fluid Phase Equilib. 1992, 78, 157−190. (63) Chen, G. J.; Sun, C. Y.; Guo, T. M. A Theoretical Revision of The Derivation of Liquid Property Expressions from an Equation of State and Its Application. Chem. Eng. Sci. 2000, 55, 4913−4923. (64) Twu, C. H.; Coon, J. E.; Cunningham, J. R. A New Generalized Alpha Function for a Cubic Equation of State Part 2. Redlich−Kwong Equation. Fluid Phase Equilib. 1995, 105, 61−69. (65) Souahi, F.; Sator, S.; Albane, S. A.; Kies, F. K.; Chitour, C. E. Development of a New Form for the Alpha Function of the Redlich− Kwong Cubic Equation of State. Fluid Phase Equilib. 1998, 153, 73− 80. (66) Redlich, O.; Kwong, J. N. S. The Thermodynamics of Solutions. V. An Equation of State. Fugacities of Gaseous Solutions. Chem. Rev. 1949, 44, 233−244. (67) Pederson, K. S.; Thomassen, P.; Fredenslund, Aa. Thermodynamics of Petroleum Mixtures Containing Heavy Hydrocarbons. 1. Phase Envelope Calculations by Use of the SRK Equation of State. Ind. Eng. Chem. Process Des. Dev. 1984, 23, 163−170.

6925

dx.doi.org/10.1021/ie2025757 | Ind. Eng. Chem. Res. 2012, 51, 6916−6925