Research Article  Open Access
Prediction of Molar Volumes of the Sudanese Reservoir Fluids
Abstract
This paper provided important experimental PVT data of the Sudanese reservoir fluids. It includes composition analysis of 11 mixtures and about 148 PVT data points of constant mass expansion (CME) tests at pressures below the bubble point. The datasets are compared with eight equations of state (EOS), namely, Peng Robinson (PR), SoaveRedlichKwong (SRK), LawalLakeSilberberg (LLS), AdachiLuSugie (ALS), SchmidtWenzel (SW), PatelTeja (PT), ModifiedNasrifarMoshfeghian (MNM), and HarmensKnapp (HK). The results of comparison reveals that, with the exception of PR and ALS EOSs, all other EOSs yield consistently a higher average absolute percent deviation (AAPD) in the prediction of molar volume; it exceeds 20% by all mixtures. The grand average AAPD of all mixtures is 17 and 16 for PR and ALS, respectively. ALS is selected to represents the mixtures. It is modified by replacing the coefficient () of the parameter () in the dominator of repulsive term by that of PR. This procedure enhanced the accuracy of ALS by 30 to 90% for individual mixtures and the grand average AAPD is significantly reduced from 16 to about 7.
1. Introduction
In the absence of the experimental PVT study, properties such as isothermal compressibility factor, factor, and formation volume factor, are determined from empirically derived correlations or equations of state (EOSs). The correlations are basically developed for crude from certain geographical region with certain hydrocarbon and nonhydrocarbon contents and API. Hence such correlations may not be valid to crude oils of geographical regions other than those for which these correlations have been developed. Although EOSs are generalized correlations, their validity to different crudes varies.
Adepoju [1] has made extensive study on Texas oil and found that Peng Robison (PR) [2], and SoaveRedlichKwong (SRK) [3] give a higher average absolute percent deviation (AAPD) in the prediction of the total volume of reservoir fluids. He obtained a good result when PR and SRK are modified by replacing the repulsion and attraction terms by that of LawalLakeSilberberg (LLS) EOS [4]. Akberzadeh et al. [5] have investigated the ModifiedNasrifarMoshfeghian (MNM) EOSs, PR, and SRK for Western Canadian heavy oils. They have shown that MNM without any volume correction predicted the densities with accuracy similar to SRK EOS with volume correction.
Jensen [6] found that AdachiLuSugie (ALS) EOS [7] is the most accurate for prediction of the phase behavior of welldefined hydrocarbon mixtures with and without a considerable content of CO_{2} or N_{2}. The ALS EOS seems to be well suited for calculation of the phase equilibrium of reservoir fluids but often proves to give inaccurate predictions of the densities of hydrocarbon mixtures [8]. By incorporating the volume translation principle of Peneloux et al. [9], ALS equation was found to give good results for hydrocarbon mixtures with and without a considerable content of CO_{2} or N_{2} [10].
Pedersen et al. [10] developed a characterization procedure for SRK coupled with the volume correction term of Peneloux et al. [9]. This procedure does not need experimental data and generally gives good prediction of saturation points and vaporliquid equilibrium. However, the model frequently calculates a too large liquid precipitation for gas condensate when simulating constant composition expansion experiments. In addition, prediction of liquid density is sometimes inaccurate, they added.
Almehaideb et al. [11] tested crude and gas samples from 17 UAE reservoirs and found that PR EOS predicted the density and bubble point pressure of UAE petroleum reservoir with an error of 9.28%. Yu and Chen [12] have evaluated PR and PT for binary and tertiary nonpolar and polar mixtures. For binary mixtures, the grand average AAPD for PR and RT is 9.69 and 10.39, respectively. For tertiary mixtures AAPD of PR and PT is 17.20 and 17.04, respectively. Kumer [13] has used 3100 data points of reservoir fluids mainly sweet and sour dry gases from various sources to evaluate the compressibility factor using eight EOSs, namely LLS, VDW, PR, RK, SRK, SW, PT, and TrebbleBishnoi (TB). He concluded that LLS is superior to other EOSs in the prediction of factor.
Sudanese crude oil has come to surface on a commercial scale in the mid of 1990s, and there is little known in open literature about its PVT properties. Therefore the purpose of this work is to select an EOS that best represents the PVT data of the Sudanese reservoir fluids. The candidate EOSs includes PR, SRK, LLS, ALS, SW, PT, MNM, and HK.
2. PVT Study
Experimental PVT data were supplied by the Ministry of Energy and Mining, Sudan, for a number of wells representing different reservoirs. The data include compositional analysis of single carbon numbers of up to Eicosanes plus (C20^{+}) and Hexatriacontanes plus (C36^{+}) and PVT of CME of bottomhole samples. Table 1 shows the composition analysis of eleven reservoir fluids lumped up to C7^{+}. Samples 1 and 2 are taken from [14], samples 3 to 11 are presented for the first time. The data also include bubble point pressure, reservoir temperature and molecular weight, and specific gravity of C7^{+}. In some PVT reports, the molecular weight of C7^{+} is not available; under such condition it is obtained using the material balance as where and are the component molecular weight and the apparent molecular weight, respectively, and is the mole fraction. The plus fraction specific gravity is calculated as The experimental error in pressure is 5 psi, temperature is 0.5°F, and cell volume is 0.3 cc as reported in PVT studies.

3. Flash Algorithm
Although flash calculation procedure is well documented, Figure 1 shows the Flash algorithm used in this work. The input data include the reservoir pressure (), bubble point pressure (), reservoir temperature (), reservoir fluid composition (), and the density and the molecular weight of C7^{+}. For pure hydrocarbon and nonhydrocarbon, the critical properties (), the molecular weight (MW), and the acentric factor () are obtained from the generalized properties tables [15]. For C7^{+} the critical properties and the acentric factor are estimated for a given molecular weight and a specific gravity from LawalTododoHeinze [16] correlations. Binary interaction parameters (BIPs) for hydrocarbonhydrocarbon, nonhydrocarbonhydrocarbon, and nonhydrocarbonnonhydrocarbon systems are taken as zero because all samples contains small amount of CO_{2} and N_{2}. The vapor and liquid molar volumes are calculated using the flash algorithm. Mathematically, the twophase flash calculation is the solution of RachfordRice equation that satisfies the equal fugacity constraint [17]. NewtonRaphson iteration scheme was employed. The calculation is initiated with value obtained using Wilson correlation. If a convergence is not obtained value is modified as [18]. The flash program is executed for eight EOSs. These are PR, SRK, LLS, ALS, SW, PT, MNM, and HK EOSs.
4. Results and Discussions
To compare EOSs to experimental data, numerous quality measurements based on statistical error analysis are computed. These include the percent deviation (PD), the average absolute percent deviation (AAPD), the minimum absolute percent deviation (), the maximum absolute percent deviation (), and the grand average AAPD. The percent deviation is defined as and the average absolute percent deviation (AAPD) is defined as whereis experimental molar volume (ft^{3}/lbmole),is the calculated molar volume, and is the number of the data points.
Table 2 shows a calculation sample of mixture number 9. The rest of the results are given in Tables 6, 7, 8, and 9. Table 3 shows the summary of statistical parameters for all mixtures. It is also shown in Table 3 the reservoir temperature, bubble point pressure, and C1 and C7^{+} content of each mixture. It should be noted that none of the mixtures at hand contains H_{2}S and they contain a little amount of CO_{2} and N_{2}. The result of comparison reveals that none of EOSs has a grand average AAPD of less than 16. EOSs such as LLS, HK, MNM, PT and SW yield consistently high AAPD of all mixtures irrespective to their bubble point pressures and C1 and C7^{+} contents. The rest of EOSs (SRK, PR, and ALS) perform better for mixtures with a higher C1 content than that of a lower C1 however, with few exceptions (e.g., mixture no. 2). In the overall evaluation ALS has the least grand average AAPD (=16) among all tested EOSs.

 
${}_{}{}^{*}\text{Total}$ number of data points. 
The reported inaccuracy is not peculiar to Sudanese reservoir fluids but it is rather a known problem associated with EOSs. Many investigators such as Coats and Smart [19] and Ahmed [20], to mention a few, have reported the inaccuracy of EOSs to reservoir fluids. The inaccuracy can be attributed to the following EOSs limitations and plus fraction properties. ()High level of uncertainty in the prediction of the critical properties and the accentric factor of plus fraction which are not measured in the laboratory. Bearing in mind that plus fraction is the main constituent of mixtures. It constitutes more than 75% in most of the mixtures at hand (mixtures no. 3 to no. 7 and mixture no. 9). The inaccuracy in prediction of critical properties and acentric factor is due to the fact that the plus fraction lumps millions of compounds that only few of them (C7 to C36) are known by measurement. Hence it is expected that as the content of the plus fraction in the mixture increases the inaccuracy of EOS increases. This may justify the relatively good performance of SRK, PR and ALS for mixtures of a lower percentage of plus fraction. ()The parameters of the attraction term , and covolume of EOS are determined based on van der Waal critical point assumption while the reservoirs temperatures in this work (cf. Table 1) are higher than the critical temperatures of N_{2}, CO_{2}, C1, and C2. For mixtures numbers 1, 2, 8, and 11, these components together constitute more than 20% of the said mixture. This means that the applicability limit of EOS is violated and hence the inaccuracy of EOSs is not a surprise. ()Lack of information on BIPs of N_{2}–CO_{2}, N_{2}hydrocarbon, and CO_{2}hydrocarbon, However, the presence of the nonhydrocarbon components (N_{2} and CO_{2}) is small in the investigated mixtures. ()All EOSs used in this work contain three [] and four [] parameters, hence besides the van der Waal critical point conditions, the parameters were determined by regressing experimental data for pure components. The commonly used experimental data for such a purpose include vapor pressure, normal boiling point and density at standard conditions ( = 15°C and = 1 atm). These data are generally for lower molecular weight components. Hence an EOS that developed on these data unlikely will suffice for reservoir fluid which contains a higher molecular weight plus fraction.
Because the inaccuracy of EOSs rests on the four reasons outlined above, a number of methods have been proposed over years to enhance the capability of EOS yet maintaining its original characteristics. These accuracy enhancement methods include the following.()Development of accurate models for the prediction of critical properties and acentric factor of plus fraction. ()Application of a volumetranslation technique such as Peneloux shift factor to the EOS.()Tuning EOS to experimental data. However, Pedersen [21] warned that “using equation of state parameters “tuned” to one specific property yields unreliable predictions of other thermodynamic properties.” ()Modification of the expression used in the denominator of the attractive term Yu et al. [22]. This method as mentioned earlier is employed by Adepoju [1]. He has made a significant improvement in the accuracy of PR and SRK by replacing their parameters by those of LLS.
In this work ALS EOS which produced the least grand average AAPD (=16) of all mixtures is considered as the candidate to predict PVT data of the Sudanese reservoir fluids. Enhancement procedure, using the technique number 4 listed above, is considered.
Prior to the employment of the modification, ALS is described asThe modification includes the replacement of the first term of the coefficient (cf. (6b)) by that of PR EOS (= 0.07780). Hence the modified form of is Remember that the modification of will automatically modify the coefficients and.
Table 4 shows calculation samples of mixtures numbers 1 and 6. Table 5 shows the summary of AAPD for all mixtures after the modification of the ALS parameters. It can be seen that the accuracy of ALS enhanced by a factor of 30 to 90%. The grand average AAPD is reduced significantly from 16 to 7.






5. Conclusion
The work provided important information on PVT data on Sudanese reservoir fluids. It includes composition analysis of fraction plus up to C7^{+} and about 148 data points of CME test (pressurevolume data) at pressures below the bubble point. The paper presents also a modified form of ALS that describes the Sudanese reservoir fluids with a good level of accuracy.
Acknowledgment
The authors acknowledge the support of the Ministry of Energy and Mining, Sudan.
References
 O. O. Adepoju, Coefficient of isothermal oil compressibility for reservoir fluids by cubic equation of state, M.Sc. thesis, University of Texas, Austin, Tex, USA, 2006.
 D.Y. Peng and D. B. Robinson, “A new twoconstant equation of state,” Industrial and Engineering Chemistry Fundamentals, vol. 15, no. 1, pp. 59–64, 1976. View at: Google Scholar
 G. Soave, “Equilibrium constants from a modified RedlichKwong equation of state,” Chemical Engineering Science, vol. 27, no. 6, pp. 1197–1203, 1972. View at: Google Scholar
 A. S. Lawal, E. T. Van der Laan, and R. K. M. Thambynayagam, “Fourparameter modification of the LawalLakeSilberberg equation of state for calculating gascondensate phase equilibria,” in Proceedings of the Annual Technical Conference and Exhibition, Las Vegas, Nev, USA, September 1985, paper SPE 14269. View at: Google Scholar
 K. Akbarzadeh, Sh. Ayatollahi, Kh. Nasrifar, H. W. Yarranton, and M. Moshfeghian, “Prediction of the densities of Western Canadian heavy oils and their SARA fractions from the cubic equations of state,” Iranian Journal of Science and Technology, Transaction B, vol. 28, no. B6, pp. 695–699, 2004. View at: Google Scholar
 B. H. Jensen, Densities, viscosities and phase equilibria in enhanced oil recovery, Ph.D. thesis, Department of Chemical Engineering, the Technical University of Denmark, Lyngby, Denmark, 1987.
 Y. Adachi, B. C.Y. Lu, and H. Sugie, “A fourparameter equation of state,” Fluid Phase Equilibria, vol. 11, no. 1, pp. 29–48, 1983. View at: Google Scholar
 K. AasbergPetersen, Bulk phase properties and phase equilibria for miscible and immiscible oil displacement processes, Ph.D. thesis progress report, Department of Chemical Engineering, the Technical University of Denmark, Lyngby, Denmark, 1989.
 A. Peneloux, E. Rauzy, and R. Freze, “A consistent correction for RedlichKwongSoave volumes,” Fluid Phase Equilibria, vol. 8, no. 1, pp. 7–23, 1982. View at: Google Scholar
 K. AasbergPetersen and E. Stenby, “Prediction of thermodynamic properties of oil and gas condensate mixtures,” Industrial and Engineering Chemistry Research, vol. 30, no. 1, pp. 248–254, 1991. View at: Google Scholar
 R. A. Almehaideb, I. Ashour, and K. A. ElFattah, “Improved Kvalue correlation for UAE crude oil components at high pressures using PVT laboratory data,” 2003. View at: Google Scholar
 M.L. Yu and Y.P. Chen, “VLE calculations by applying a modified perturbed hard sphere EOS,” Fluid Phase Equilibria, vol. 129, no. 12, pp. 21–35, 1997. View at: Google Scholar
 N. Kumer, Compressibility factor for natural and sour reservoir gases by correlations and cubic equations of state, M.Sc. thesis, Texas Tech University, Lubbock, Tex, USA, 2004.
 A. A. Rabah and S. A. Mohamed, “Prediction of molar volumes of undersaturated Sudanese reservoir fluids,” submitted to Journal of Petroleum Science and Engineering. View at: Google Scholar
 A. Danesh, PVT and Phase Behavior of Petroleum Reservoir Fluids, Elsevier Science, Amsterdam, The Netherlands, 2nd edition, 1998.
 A. O. Tododo, Thermodynamically equivalent pseudo components for compositional reservoir simulation models, Ph.D. dissertation, Texas Tech University, Lubbock, Tex, USA, 2005.
 A. Rabah and S. Kabelac, “Flow boiling of R134a and R134a/propane mixtures at low saturation temperatures inside a plain horizontal tube,” Journal of Heat Transfer, vol. 130, no. 6, Article ID 061501, 9 pages, 2008. View at: Publisher Site  Google Scholar
 C. H. Whitson and M. R. Brule', Phase Behavior, vol. 20 of Monograph: SPE Henry L. Doherty Series, SPE, Richardson, Tex, USA, 2000.
 K. H. Coats and G. T. Smart, “Application of a regressionbased EOS PVT program to laboratory data,” SPE Reservoir Engineering, vol. 1, no. 3, pp. 277–299, 1986. View at: Google Scholar
 T. Ahmed, “On equation of state,” in Proceedings of the SPE Latin American and Caribbean Petroleum Engineering Conference, vol. 1, pp. 539–559, Buenos Aires, Argentina, 2007. View at: Google Scholar
 K. S. Pedersen, “On the dangers of “Tuning” equation of state parameter,” Society of Petroleum Engineers of AIME, p. 14487, 1985. View at: Google Scholar
 J. M. Yu, S. H. Huang, and M. Radosz, “Phase behavior of reservoir fluids: supercritical carbon dioxide and cold lake bitumen,” Fluid Phase Equilibria, vol. 53, pp. 429–438, 1989. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2010 A. A. Rabah and S. A. Mohamed. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.