Impact of Lumping Techniques for Fluid ...

33 downloads 0 Views 774KB Size Report
To using this rule, normalized mole fraction of the component i in ..... Danesh, Ali, Dong-hai Xu, Adrian C. Todd, and U. Heriot-Watt, 1992. A Grouping Method To ...
Australian Journal of Basic and Applied Sciences, 7(1): 320-333, 2013 ISSN 1991-8178

Impact of Lumping Techniques for Fluid Characterization in Gas Condensate Reservoir 1

Hamidreza moghadamzadeh, 2Hojjatolla Maghsoodloorad, 3Adel zarabpour, 2,4Ali Hemmati, 5Shahin Shahsavari, 6Shadab Shahsavari 1

Islamic Azad University South Tehran Branch, P.O. Box 19585-466, Tehran, Iran School of Chemical Engineering, Iran University of Science and Technology (IUST), Narmak, Tehran, Iran 3 Islamic Azad University Tabas Branch, 4-police University 5 Department of Chemical Engineering, University of Sistan and Baluchestan, Zahedan, Iran 6 Department of Chemical Engineering, Varamin (Pishva) Branch, Islamic Azad University, Tehran, Iran 2

Abstract: Several problems are associated with “regrouping” the original components into a smaller number without losing the predicting power of the equation of state. These problems include (A) How to select the groups of pure components to be represented by one pseudo-component. (B) What mixing rules should be used for determining the EOS constants (pc, Tc, and ω) for the new lumped pseudocomponents? In this paper the results of the lumping methods are compared with the data for gas condensate systems to find accurate lumping method for accurate fluid characterization. For this purpose, Whitson’s lumping scheme, Pedersen lumping scheme, Danesh et al. lumping scheme, Lee et al. lumping scheme, Behras and Stanndler lumping scheme are used as lumping scheme and Lee mixing rules was used for calculation of thermodynamics properties of lumped group. To compare the Impact of each Lumping method, its effect was examined on the phase diagram. The results show that the Whitson lumping scheme (it is grouping the hydrocarbon from C7 to C16+ into 4 groups) can predict the phase diagram of this gas condensate systems fluid with good approximation, the Pedersen lumping scheme (it is grouping the hydrocarbon from C7 to C16+ into 3 groups ) can predict the phase diagram of this gas condensate systems fluid somewhat and the Pedersen lumping schemes and Danesh et al. lumping schemes (it is grouping the hydrocarbon from C7 to C16+ into 4 groups) approximately predicted the same phase diagram and the Lee et al. lumping scheme (it is grouping the hydrocarbon from C7 to C16+ into 3 groups) can predict the phase diagram of this gas condensate systems fluid exactly. Therefore, the Lee at al. lumping scheme has very good accuracy compared to the other methods. Also, All Lumping Scheme has an effect on dew point line of phase diagram. Key words: Lumping, Groping, Pseudoization, Pseudo component, Mixing rule. INTRODUCTION Petroleum is composed of various hydrocarbon and non-hydrocarbon components. For characterizing and simulating the oil reservoir, the reservoir fluid must be characterized. The large number of components is necessary to describe the hydrocarbon mixture for accurate phase behavior modeling frequently burdens EOS calculations. Generally, with a sufficiently large number of pseudo-components used in characterizing the heavy fraction of a hydrocarbon mixture, a satisfactory prediction of the PVT behavior by the equation of state can be obtained. However the cost and computing time can be increased significantly with the increased number of components in the system. Therefore, strict limitations are placed on the maximum number of components that can be used in compositional models and the original components have to be lumped into a smaller number of pseudo components (Ahmed 2010, 1989). Several problems are associated with “regrouping” the original components into a smaller number without losing the predicting power of the equation of state. These problems include How to select the groups of pure components to be represented by one pseudo-component each. Several unique techniques have been published that can be used to address the above lumping problems; notably the methods proposed by Mehra et al, Schlijper, Gonzalez, Colonomos, and Rusinek, Aguilar and McCain, Pedersen et al. (Ahmed 2010, 1989; Hong 1982; Guo and Du 1989; Lee et al. 1981; Montel and Gouel 1984; Al-Meshari 2004), Whitson (Whitson 1983), Lee et al. (Lee et al. 1981), Montel and Gouel (Montel and Gouel 1984), Behrens and Sandler (Behrens and Sandler 1988), Hong (Hong 1982) and Danesh et al. (Danesh et al. 1992). Corresponding Author: Hamidreza moghadamzadeh, Islamic Azad University South Tehran Branch, P.O. Box 19585-466, Tehran, Iran E-mail: [email protected]

320

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

For describing the phase behavior of a fluid mixture, critical pressure, critical temperature, critical volume, acentric factor, molecular weight and binary interaction parameters of component must be specified. Different mixing rules were evaluated for, characterizing the lumped C7+ fraction from its constituent component properties (Ahmed 2010; Hong 1982; Adepoju 2006; Joergensen and Stenby 1995). Impact of these lumping schemes on one system, especially gas condensate systems, isn’t studied. In this work, the accuracy and impact of Whitson, Pedersen, Danesh et al., Lee et al., Behras and stanndler lumping scheme on phase diagram of the gas condensate systems are researched. Lee’s mixing rule is one of these mixing rules which are used in this research. 2. Methods: 2.1 Experiment: A sample of reservoir fluid is prepared from one of gas condensate oil fields. In the laboratory, it composition is determined by using gas chromatography. Also, the properties of heavy oil fraction were measured with the existing standards methods in the laboratory. 2.2 Lee’s mixing rules: There are many mixing rules which obtained by researchers. Lee’s mixing rules is one of useful rules. Lee et al. (1979) (Ahmed 2010) used Kay’s mixing rules as the characterizing approach for determining the properties of the lumped fractions in their suggested regrouping model. To using this rule, normalized mole fraction of the component i in the lumped fraction must be defined. It was defined by:

φi =

zi L

∑z

(1) i

After that, Lee et al. are proposed the following rules for calculation of thermodynamic property of multi component mixture: (2)

γL =

MwL L

∑ φ Mw i

VcL = ∑

(3)

i

φi MwiVci

(4)

MwL

L

Pcl = ∑ φi Pci

(5)

L

TcL = ∑ φi Tci

(6)

L

ω l = ∑ φi ω i

(7)

3. Results: In this paper, Whitson’s lumping scheme, Pedersen lumping scheme, Danesh et al. lumping scheme, Lee et al. lumping scheme, Behras and Stanndler lumping scheme are used. Calculation of the method described in the appendix. A compositional analysis of gas condensate systems fluid described in Table 1 and the molecular weight and the specific gravity of C16+ is equal to 259 and 0.908 respectively. The phase diagram of this gas condensate system is shown in Figure 1. 3.1 Whitson’s lumping scheme: Whitson (Whitson 1983) proposed a regrouping scheme whereby the compositional distribution of the C7+ fraction is reduced to only a few multiple carbon- number (MCN) groups. Whitson suggested that the number of MCN groups is necessary to describe the plus fraction is given by the following empirical rule:

N G = Int [1 + 3.3 log( N − n )]

(8)

321

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

Not that, for black oil system this number probably can be reduced by one.6 The molecular weights separating each MCN group are calculated from the following expression:

 MwN MwI = Mwn   Mwn

I

 NG  

(9)

Table 1: Composition of reservoir fluid Component C1 C2 C3 i-C4 n-C4 i-C5 n-C5 C6 C7

zi 0.413506 0.0403 0.2153 0.0539 0.0543 0.0515 0.0519 0.1039 0.00347

Component C8 C9 C10 C11 C12 C13 C14 C15 C16+

zi 0.00268 0.00207 0.00159 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

Fig. 1: Phase diagram of reservoir fluid Components with molecular weight falling within the boundaries of MI–1 to MI are included in the I’th MCN group. By applying the Whitson lumping scheme can see that it is grouping the hydrocarbon from C7 to C16+ into 4 groups. The results of Whitson’s lumping scheme are shown in Table 2. The impact of this scheme on the phase diagram is shown in Figure 2. Thus, this figure shows that the Whitson lumping scheme can predict the phase diagram of this gas condensate systems fluid with good approximation. Table 2: Whitson's lumping scheme results Group I 1 2 3

4

component

zi

C7 C8 C9 C10 C11 C12 C13 C14 C15

0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.0016 71

C16+

zI

MwI

γI

Vcl (ft3/lb)

Pcl (psia)

Tcl (R)

ωl

0.00822

105.9

0.746

0.0627

424

1020

0.3076

0.002826

139.7

0.787

0.0628

339.7

1144.5

0.4000

0.002246

172.9

0.814

0.0631

288

1230.6

0.4794

0.002108

248

0.892

0.0637

223.3

1433

0.6531

Whitson's lumping scheme accuracy can be attributed to function of the number of MCN groups (equation 8) and iterative method to find the molecular weight of each group Determination of the number of groups is the excellent characteristics of this approach.

322

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

Fig. 2: Impact of Whitson’s lumping scheme on phase diagram 3.2 Pedersen Lumping Scheme: Pedersen et al. (Ahmed 2010, 1989; Hong 1982; Guo and Du 1989; Lee et al. 1981; Montel and Gouel 1984; Al-Meshari 2004) suggested the process of grouping the components on the basis of each group containing approximately the same weight fraction (equal weight fraction), which will give all hydrocarbon segments of the C7+ fractions equal importance. The C7+ fraction is divided into three or more groups that, by weight, are equal size approximately. The weight for each pseudo component, Wj, can be calculated as: nc

w j = ∑ zi × Mwi

(10)

i =1

By applying the Pedersen lumping scheme can see that it is grouping the hydrocarbon from C7 to C16+ into 3 groups (3 groups is selected according to some restrictions). The results of Pedersen lumping scheme are shown in Table 3. The phase diagram of this system is under the impact of Pedersen lumping scheme is shown in Figure 3. Therefore, this figure shows that the Pedersen lumping scheme can predict the phase diagram of this gas condensate systems fluid somewhat. Table 3: Pedersen lumping scheme results Group I 1

2

3

component

zi

C7 C8 C9 C10 C11 C12 C13 C14 C15 C16+

0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

zI

MwI

γI

Vcl (ft3/lb)

Pcl (psia)

Tcl (R)

ωl

0.00615

100.794

0.73657

0.06277

438.184

1007.2

0.29395

0.00584

136.522

0.78409

0.06279

348.896

1123.4

0.39027

0.003404

222.709

0.86866

0.06353

244.507

1363.8

0.59401

As you can see, a simple relationship exists between the concentration and molecular weight in Pedersen lumping scheme. Choosing the number of groups is based on the limitations and experiences, so, the accuracy of this method would not be high. 3.3 Danesh et al. lumping scheme: Danesh et al. (Danesh et al. 1992) proposed a lumping method based on the concentration and molecular weight of compounds in a mixture. This grouping method arranged the original components in order of their normal boiling point temperatures and grouped together in ascending order to form Np groups so that the values of Σ (zilnMi) for all the groups become nearly equal, (Quasi-equal-weight criterion). That is:

  l  1  n  × ∑ z i × ln Mwi  ≤ 0 ∑ z i × ln Mwi −     i =1  N p  i =1 And

323

(11)

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

  l +1  1 n ∑ z i × ln Mwi  ≥ 0 ∑ z i ln Mwi −    N p  i =1   i =1 i = 1,2,..., n, l = 1,2,...N p

(12)

Fig. 3: Impact of Pedersen lumping scheme on phase diagram According to the Danesh et al. lumping scheme can see that it is grouping the hydrocarbon from C7 to C16+ into 4 groups. The results of Danesh et al. lumping scheme are shown in Table 4. The phase diagram of this system is under the impact of Danesh et al. lumping scheme is shown in Figure 4. The predicted phase diagram of this gas condensate systems fluid by the Danesh et al. lumping scheme is look like the phase diagram was predicted by Pedersen lumping scheme. Figure 4 show that this lumping scheme can predict the phase diagram somewhat. In this method, choosing the number of groups was done similar to the Pedersen lumping scheme. In Danesh et al. lumping scheme some limitation is defined for choosing the molecular weight of groups, so it would be more accurate than the Pedersen lumping scheme. Table 4: Danesh et al. lumping scheme results Group zI component zi I 1 C7 0.00347 0.00347

Pcl (psia) 453

Tcl (R) 985

ωl

0.727

Vcl (ft3/lb) 0.06289

MwI

γI

96

96

2

C8 C9

0.00268 0.00207

0.00475

113.101

0.75719

0.06261

403.312

1045.59

113.101

3

C10 C11 C12

0.001596 0.00123 0.00095

0.00377

145.045

0.79169

0.06288

330.17

1159.3

145.045

4

C13 C14 C15 C16+

0.00073 0.000566 0.000437 0.001671

0.003404

222.709

0.86866

0.06353

244.507

1363.78

222.709

3.4 Lee et al. lumping scheme: Lee et al. (Lee et al. 1981) devised a simple procedure that can be served as a guideline for lumping the oil fractions. The idea for the procedure lies in the physical reasoning that crude-oil fractious having relatively close physicochemical properties can be represented fairly accurately by a single fraction, the closeness of these properties can be reflected by the slopes of curve when these properties are plotted against some characteristic independent variables. According to the Lee et al. lumping scheme can see that it is grouping the hydrocarbon from C7 to C16+ into 3 groups. The results of Lee et al. lumping scheme are shown in Table 5. The accuracy of this scheme can be seen in the phase diagram in Figure 5. Hence, Figure 5 shows that the Lee et al. lumping scheme can predict the phase diagram of this gas condensate systems fluid exactly.

324

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

Fig. 4: Impact of Danesh et al. lumping scheme on phase diagram Lee et al. lumping scheme is more complex than lumping schemes which is discussed in this paper but it gives highly accurate answer. On the other hand, this lumping scheme is determined number of groups which is the key points of this method. 3.5 Behras And Stanndler Lumping Scheme: Behrens and Standler (Behrens and Sandler 1988) have been used a semi continuous thermodynamic description to model the C7+ fraction for equation of state (EOS) calculations. The step-by-step procedure will be illustrated. Step 1: Finding the ends of the distribution the ends of the distribution can be calculated by: (13)

A= starting carbon number - 1/2 B= ending carbon number + 1/2 Table 5: Lee et al. lumping scheme results Group component zi I C7 0.00347 1 C8 0.00268 C9 0.00207 C10 0.001596 C11 0.00123 2 C12 0.00095 C13 0.00073 C14 0.000566 C15 0.000437 3 C16+ 0.001671

(14)

zI

MwI

γI

Vcl (ft3/lb)

Pcl (psia)

Tcl (R)

ωl

0.00615

100.794

0.73657

0.06277

438.184

1007.2

0.29395

0.007573

148.237

0.79453

0.0629

331.518

1155.6

0.41792

0.001671

259

0.908

0.0638

215

1467

0.68

The parameters of the semi continuous distribution in this system are A = 7 - 1/2 = 6.5 B =16 + 1/2 =16.5. Step 2 – determining the slope of the distribution The slope of the distribution is found from the average carbon number Cn which is related to the where average MW approximately as −

c

n

=

( MW + 4) 14

(15)

The slope of the distribution, α, is then found from the following equation:

−  (B − A)e − Bα  = c − A +  − Aα  n α − e − Bα  e

1

(16)

325

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

Fig. 5: Impact of Lee et al. lumping scheme on phase diagram If the upper endpoint were infinity, the term in the brackets would vanish and l/α could be found explicitly. Since the upper endpoint is finite, α is found by successive substitutions. Slope of the distribution: −

c

n

=

(141.25 + 4) = 10.375 14

Then

 6 × e −12.5*α = 10.375 − 6.5 +  −6.5*α α − e −12.5*α e 1

 → α = 0.0011  

Step 3 - Finding integration points and weights The integration range can be calculated by Equation (17).

c = ( B − A)α

(17)

Integration points and weights C = (16.5-6.5)*0.0011 = 0.011 Because C< 0.3 so with this method cannot lump this fluid. Conclusions: By comparing the phase diagram for reservoir fluid (Figure 1), and data that is lumped with some lumping schemes like Whitson’s Lumping Scheme (Figure 2), Pedersen lumping scheme (Figure 3), Danesh lumping scheme (Figure 4), Lee et al. lumping scheme (Figure 5) and also plot all lumping scheme and the phase diagram of reservoir fluid on one diagram (Figure 6) can be seen that:

Fig. 6: All lumping scheme

326

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

1- With the Pedersen Lumping Scheme and Danesh et al. Lumping Scheme can reach to the approximately same the phase diagram. 2- Lee at al. Lumping Scheme has very good accuracy related to the other methods in characterization of gas condensate reservoir. 3- All Lumping Scheme has an effect on dew point line. REFERENCES Adepoju, O.O., 2007. Coefficient of Isothermal Oil Compressibility- A Study for Reservoir Fluids by Cubic Equation-of-State. VDM Verlag Dr. Mueller e.K Press. Ahmed, T.H., 1989. Hydrocarbon phase behavior, Contributions in petroleum geology & engineering 7. Gulf Pub. Co Press. Ahmed, Tarek H., 2010. Reservoir engineering handbook. Gulf Professional Pub Press. Al-Meshari, A.A., 2004. New strategic method to tune equation-of-state to match experimental data for compositional simulation, Ph.D. thesis, Texas A&M University, Texas Behrens, R.A., and S.I. Sandler, 1988. The Use of Semicontinuous Description To Model the C7+ Fraction in Equation of State Calculations. SPE Reservoir Engineering, 3(3): 1041-1047. Danesh, Ali, Dong-hai Xu, Adrian C. Todd, and U. Heriot-Watt, 1992. A Grouping Method To Optimize Oil Description for Compositional Simulation of Gas-Injection Processes. SPE Reservoir Engineering, 7(3): 343-348. Guo, T.M., and L. and Du., 1989. A New Three Parameter Cubic Equation of State for Reservoir Fluids– Part II. Application to Reservoir Crude Oils. SPE pp: 19373. Hong, K.C., 1982. Lumped-Component Characterization of Crude Oils for Compositional Simulation. In SPE Enhanced Oil Recovery Symposium, Tulsa, Oklahoma Joergensen, M., and E.H. Stenby, 1995. Optimization of Pseudo-component Selection for Compositional Studies of Reservoir Fluids. In 70th Annul Technical Conference and Exhibition. Daltas, Taxas. Lee, S.T., R.H. Jacoby, W.H. Chen, W.E. Culham, and Gulf Research and Development Co., 1981. Experimental and Theoretical Studies on the Fluid Properties Required for Simulation of Thermal Processes. SPE Journal., 21(5): 535-550. Montel, F., and P.L. Gouel., 1984. A New Lumping Scheme of Analytical Data for Compositional Studies. In SPE Annual Technical Conference and Exhibition. Houston, Texas. Whitson, C., 1983. Characterizing Hydrocarbon Plus Fractions. SPE Journal, 23(4): 683-694. Nomanclature: I = 1, 2. . . Ng i = component i Int = integer L = lumped fraction MwC7 = molecular weight of C7 Mwi = Molecular weight MwN+ = molecular weight of the last reported component in the extended analysis of the hydrocarbon system n = number of carbon atoms of the first component in the plus fraction, i.e., n = 7 for C7+ N = number of carbon atoms of the last component in the hydrocarbon system Ng = number of MCN groups Pci = Critical pressure Tci = Critical temperature vci = Critical volume Wj = weight for each pseudo component Zi = mole fraction Фi = weighting factor ωi = Acentric factor Appendix A: A.1 Whitson’s lumping scheme: Determining the appropriate number of pseudo-components forming in the C7+ by using the Whitson’s Lumping Scheme Step 1. Determining the molecular weight of each component in the system:

327

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

TableA-1: Molecular weight of components component C7 C8 C9 C10 C11 C12 C13 C14 C15 C16+

zi 0.00347 0.00268 0.00207 0.00159 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

Mwi 96 107 121 134 147 161 175 190 206 259

Step 2. Calculating the number of pseudo-components: For calculation of pseudo-components having N=16 n=7

N g = Int [1 + 3.3 log( N − n)]

Therefore,

N g = Int [1 + 3.3 log(16 − 7)]

N g = Int [4.15] Ng = 4

Step 3. Determining the molecular weights separating the hydrocarbon groups:

 MwN + MwI = MwC 7   MwC 7  259  MwI = 96   96 

I

  

I

Ng

4

Therefore,

MwI = 96(2.698)

I

4

Table A-2: Molecular weight of each group I MWI 1 123.0359 2 157.6857 3 202.0937 4 259.008

• First pseudo-component: The first pseudo-component includes all components with a molecular weight in the range of 96 to 123. This group then includes C7, C8, and C9. • Second pseudo-component: The second pseudo-component contains all components with a molecular weight higher than 123 to a molecular weight of 158. This group includes C10 and C11. • Third pseudo-component: The third pseudo-component includes components with a molecular weight higher than 158 to a molecular weight of 202. Therefore, this group includes C12, C13, and C14. • Fourth pseudo-component: This pseudo-component includes all the remaining components, i.e., C15 and C16+. Table A-3: Molecular percentofeach group Group I component C7 1 C8 C9 2 C10

zi 0.00347 0.00268 0.00207 0.001596

328

zI 0.00822 0.002826

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

C11 C12 C13 C14 C15

3 4

0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

C16+

0.002246 0.002108

Zl is calculated by the equation of: l

∑z

i

= zl

there are numerous ways to mix the properties of the individual components, all giving different properties for the pseudo-components, the choice of a correct mixing rule is as important as the lumping scheme. The mixing rules can be employed to characterize the pseudo-component in terms of its pseudo physical and pseudo-critical properties. In this Example Using Lee’s mixing rules to determine the physical and critical properties of the four pseudo-components TableA-4: Thermodynamics property of components Group I 1 2 3 4

component

zi

C7 C8 C9 C10 C11 C12 C13 C14 C15 C16+

0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

zI

0.00822 0.002826 0.002246 0.002108

γi from appendix 0.727 0.748 0.768 0.782 0.793 0.804 0.815 0.825 0.826 0.908

Mwi from appendix 96 107 121 134 147 161 175 190 206 259

vci from appendix 0.06289 0.06264 0.06258 0.06273 0.06291 0.06306 0.06311 0.06316 0.06325 0.0638*

Pci from appendix 453 419 383 351 325 302 286 270 255 215*

Tci from appendix 985 1036 1058 1128 1166 1203 1236 1270 1304 1467

* Calculated Step 2. Calculate the physical and critical properties of each group by applying Equations (4-23) through (4-29) to give: Table A-5: φi of components Group I component C7 1 C8 C9 C10 2 C11 C12 3 C13 C14 C15 4 C16+

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

zl 0.00822 0.002826 0.002246 0.002108

Table A-6: Thermodynamics property of each group by Whitson’s lumping scheme Group 1 2 3 zI 0.00822 0.002826 0.002246 MI 105.9 139.7 172.9 0.746 0.787 0.814 γI Vcl 0.0627 0.0628 0.0631 Pcl 424 339.7 288 Tcl 1020 1144.5 1230.6 ωl 0.3076 0.4000 0.4794

4 0.002108 248 0.892 0.0637 223.3 1433 0.6531

A.2 Pedersen Lumping Scheme: Table A-7: Molecular weight of components component zi C7 0.00347 C8 0.00268 C9 0.00207 C10 0.00159 C11 0.00123 C12 0.00095 C13 0.00073

Mwi 96 107 121 134 147 161 175

329

φi=zi/zl 0.423171 0.326829 0.252439 0.57 0.439286 0.431818 0.331818 0.257273 0.208095 0.795714

zi Mwi 0.3331 0.2868 0.2505 0.2131 0.1808 0.153 0.1278

ωi from appendix 0.280 0.312 0.348 0.385 0.419 0.454 0.484 0.516 0.550 0.68*

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

C14 C15 C16+

0.000566 0.000437 0.01671

190 206 259

0.1075 0.09 4.3279

∑ z Mw =2.1753 i

i

Let take 3 groups so wl for each groups is equal to:

1 2.1753 z i Mwi = = 0.7254 ∑ 3 Wl = 3 According to this value of wl groups are Table A-8: Molecular percent of each group Group I component C7 1 C8 C9 C10 2 C11 C12 C13 C14 3 C15 C16+ By using Lee’s mixing rules: Table A-9: φi of components Group I component C7 1 C8

2

3

C9 C10 C11 C12 C13 C14 C15 C16+

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

zI 0.00615

0.00584

0.003404

φi=zi/zl 0.564228 0.435772 0.354452 0.27226 0.210616 0.162671 0.214454 0.166275 0.128378 0.490893 0.564228

zI 0.00615

0.00584

0.003404

Table A-10: Thermodynamics property of each group by Behras and Stanndler lumping scheme Group 1 2 zI 0.00615 0.00584 MI 100.7935 136.5223 0.736566 0.784092 γI Vcl 0.062774 0.062787 Pcl 438.1837 348.8955 Tcl 1007.2 1123.4 ωl 0.293945 0.390271

3 0.003404 222.7089 0.868659 0.063528 244.5065 1363.8 0.594009

A.3 Danesh et al. Lumping scheme: Table A-11: Molecular weight of components component zi C7 0.00347 C8 0.00268 C9 0.00207 C10 0.00159 C11 0.00123 C12 0.00095 C13 0.00073 C14 0.000566 C15 0.000437 C16+ 0.001671

Ln(Mwi) 4.5643 4.6728 4.7958 4.8978 4.9904 5.0814 5.1648 5.247 5.3279 5.5568

zi *Ln(Mwi) 0.015838 0.012523 0.009927 0.007788 0.006138 0.004827 0.00377 0.00297 0.002328 0.009285

∑ z ln(Mw ) =0.075396 i

Let take 4 groups, so groups are:

330

i

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

Table A-12: Molecular percent of each group Group I component 1 C7 C8 2 C9 C10 3 C11 C12 C13 C14 4 C15 C16+ By using Lee’s mixing rules:

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

Table A-13: φi of components Group I component 1 C7 C8 2 C9 C10 C11 3 C12 C13 C14 4 C15 C16+

zI 0.00347 0.00475 0.00377

0.003404

φi=zi/zl 1 0.564211 0.435789 0.421751 0.32626 0.251989 0.214454 0.166275 0.128378 0.490893

zI 0.00347

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

0.00475 0.00377

0.003404

Table A-14: Thermodynamics property of each group by Danesh et al. lumping scheme Group 1 2 3 zI 0.00347 0.00475 0.00377 MI 96 113.1011 145.0451 0.727 0.757193 0.791691 γI Vcl 0.06289 0.062612 0.062882 Pcl 453 403.3116 330.1698 Tcl 985 1045.587 1159.297 ωl 0.28 0.327688 0.41348

4 0.003404 222.7089 0.868659 0.063528 244.5065 1363.779 0.594009

A.4 Lee et al. Lumping scheme: First plot zi vs. Tb and Mwi vs. Tb Zi = 0.2926e-0.012*Tb

Zi vs. Tb

0.004 0.0035 0.003

Zi

0.0025 0.002 0.0015 0.001 0.0005 0 0

100

200

300

400

500

Tb

Fig. A -1: Molecular percent of each components with bubble point temperature

331

600

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

Mwi vs. Tb 250 200

Mwi

Mwi = 0.7571*Tb - 203.2 150 Mwi = 0.5224*Tb - 95.878 100 50 0 0

100

200

300

400

500

600

Tb

Fig. A -2: Molecular weight of components with bubble point temperature Set property of zi number 1 and property of Mwi number 2

wi =

1 2 −



component

z

Mwi

Tb

m i1

m i2

mi1

mi 2

Mi

C7

0.00347

96

366

-4.3454E-05

0.5224

1

0.69

0.845

C8

0.00268

107

390

-3.258E-05

0.5224

0.7498

0.69

0.7199

C9

0.00207

121

416

-2.3848E-05

0.2524

0.5488

0.69

0.6194

C10

0.00159

134

439

-1.8096E-05

0.2524

0.4164

0.69

0.5532

C11

0.00123

147

461

-1.3898E-05

0.7571

0.3198

1

0.6599

C12

0.00095

161

482

-1.0802E-05

0.7571

0.2486

1

0.6243

C13

0.00073

175

501

-8.5996E-06

0.7571

0.1979

1

0.599

C14

0.000566

190

520

-6.8463E-06

0.7571

0.1576

1

0.5788

C15

0.000437

206

539

-5.4505E-06

0.7571

0.1254

1

0.5627

Table A-15: Molecular percent of each group Group I component C7 1 C8 C9 C10 C11 2 C12 C13 C14 C15 3 C16+ By using Lee’s mixing rules: Table A-16: φi of components Group I component C7 1 C8 C9 C10 C11 2 C12 C13 C14

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566 0.000437 0.001671

zi 0.00347 0.00268 0.00207 0.001596 0.00123 0.00095 0.00073 0.000566

332

zI 0.00615

0.007573

0.001671

zI 0.00615

0.007573

φi=zi/zl 0.564228 0.435772 0.273339 0.209956 0.162419 0.125446 0.096395 0.074739

Aust. J. Basic & Appl. Sci., 7(1): 320-333, 2013

C15 C16+

3

0.000437 0.001671

0.057705 1

0.001671

Table A-17: Thermodynamics property of components by Lee lumping scheme Group 1 2 zI 0.00615 0.007573 MI 100.7935 148.2374 0.736566 0.794528 γI Vcl 0.062774 0.062897 Pcl 438.1837 331.5179 Tcl 1007.2 1155.6 ωl 0.293945 0.41792 A.5 Behras and Stanndler lumping scheme Staring component Mole fraction (xD) MW

C7 0.0154 141.25

Step 1: The parameters of the semi continuous distribution in this example are A = 7 - 1/2 = 6.5 B =16 + 1/2 =16.5. Step 2: determine the slope of the distribution −

c

n

=

3 0.001671 259 0.908 0.0638 215 1467 0.68

(141.25 + 4) = 10.375 14

Then

 6 × e −12.5*α  = 10.375 − 6.5 +  −6.5*α α − e −12.5*α  e 1

α = 0.0011 Step 3 - Find integration points and weights C = (16.5-6.5)*0.0011 = 0.011 Because C< 0.3 so with this method cannot lump this oil.

333