Rapid Flash Calculations for Compositional Simulation
- Yinghui Li (U. of Texas at Austin) | Russell T. Johns (U. of Texas Austin)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- October 2006
- Document Type
- Journal Paper
- 521 - 529
- 2006. Society of Petroleum Engineers
- 4.2 Pipelines, Flowlines and Risers, 5.4.9 Miscible Methods, 4.6 Natural Gas, 5.3.2 Multiphase Flow, 5.3.1 Flow in Porous Media, 5.6.4 Drillstem/Well Testing, 5.2.2 Fluid Modeling, Equations of State, 5.5 Reservoir Simulation, 5.2 Reservoir Fluid Dynamics, 4.3.4 Scale, 5.4.2 Gas Injection Methods, 5.2 Fluid Characterization, 5.2.1 Phase Behavior and PVT Measurements
- 2 in the last 30 days
- 1,522 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
The computational time for conventional flash calculations increases significantly with the number of components, making it impractical for use in many fine-grid compositional simulations and other applications. Previous research to increase flash-calculation speed has been limited to those with zero binary interaction parameters (BIPs) or approximate methods based on an eigenvalue analysis of the binary interaction matrix. Practical flash calculations, however, nearly always have some nonzero BIPs. Further, the accuracy and speed of the eigenvalue methods varies depending on the choice and number of the dominant eigenvalues.
This paper presents a new and simple method for significantly increasing the speed of flash calculations for any number of nonzero BIPs. The approach requires the solution of up to six reduced parameters regardless of fluid complexity or the number of components and is based on decomposing the BIPs into two parameters using a simple quadratic expression. The new approach is exact in that the equilibrium-phase compositions for the same BIPs are identical to those with the conventional flash calculation; no eigenvalue analysis is required. Further, the new approach eliminates the Rachford-Rice procedure (1952) and is more robust than the conventional flash-calculation procedure. We demonstrate the new approach for several example fluids and show that speedup by a factor of approximately 3 to 20 is obtained over conventional flash calculations, depending on the number of components.
Gas injection into oil reservoirs results in complex interactions of flow with phase behavior that often are not modeled accurately by black-oil simulation. This is especially true for miscible or nearly miscible floods in which significant mass transfer occurs between the hydrocarbon phases. Such floods are best modeled by compositional simulation.
A significant disadvantage of compositional simulation, however, is that it is much more computationally intensive than black-oil simulation. The primary reason for the increased computational time is the result of solving repeated flash calculations with cubic equations of state (EOS). Research has shown that EOS flash calculations can occupy 50 to 70% of total computational time in compositional simulations (Stenby and Wang 1993; Chang 1990). This is also true for other applications, such as multiphase flow in pipelines.
The use of fewer pseudocomponents can reduce the flash computation time, but fewer components results in less accuracy (Hong 1982; Liu 2001; Egwuenu et al. 2005). This is especially true in multicontact miscible displacements, in which miscibility is developed by a combined condensing/vaporizing drive process (Zick 1986; Johns et al. 1993; Egwuenu et al. 2005). Fluid characterization by pseudocomponent models can be improved by tuning to the analytical minimum miscibility enrichment or minimum miscibility pressure (Johns et al. 1994), but those models still require significant computational time, even for fewer pseudocomponents.
Another way to reduce computation time is to reduce the number of gridblocks. With coarse grids, however, numerical dispersion is large, which may cloud the results (Solano et al. 2001). Ideally, fine grids should be used that better match the level of dispersion found at the field scale.
More recently, methods have been examined to find reduced parameters for flash calculations. Michelsen (1982a, 1982b, 1986) significantly increased flash-calculation speed by finding three reduced parameters, regardless of the number of components. His method, however, assumes zero BIPs, which is too restrictive for real fluid characterization. Michelsen also gave a practical method for stability calculations using the tangent plane distance (TPD) (Michelsen 1982b).
|File Size||1 MB||Number of Pages||9|
Chang, Y. 1990. Development andApplication of an Equation of State Compositional Simulator. PhD dissertation,U. of Texas, Austin, Texas.
Egwuenu, A.M., Johns, R.T., and Li, Y.2005. Improved FluidCharacterization for Miscible Gas Floods. Paper SPE 94034 presented at theSPE Europec/EAGE Annual Conference, Madrid, Spain, 13-16 June. DOI:10.2118/94034-MS.
Firoozabadi, A. 1998. Thermodynamicsof Hydrocarbon Reservoirs. New York City: McGraw-Hill.
Firoozabadi, A. and Pan, H. 2002. Fast and Robust Algorithm forCompositional Modeling: Part I—Stability Analysis Testing. SPEJ7 (1): 78-89. SPE-77299-PA. DOI: 10.2118/77299-PA.
Hearn, C.L. and Whitson, C.H. 1995. Evaluating Miscible and ImmiscibleGas Injection in the Safah Field, Oman. Paper SPE 29115 presented at theSPE Reservoir Simulation Symposium, San Antonio, Texas, 12-15 February. DOI:10.2118/29115-MS.
Hendriks, E.M. 1988. Reduction Theoremfor Phase Equilibrium Problems. Ind. Eng. Chem. Res.27 (9): 1728-1733.
Hendriks, E.M. and van Bergen, A.R.D.1992. Application of a Reduction Method to Phase Equilibria Calculations.Fluid Phase Equilibria 74: 17-34. DOI: http://dx.doi.org/10.1016/0378-3812(92)85050-I.
Hong, K.C. 1982. Lumped-Component Characterization ofCrude Oils for Compositional Simulation. Paper SPE 10691 presented at theSPE Enhanced Oil Recovery Symposium, Tulsa, 4-7 April. DOI:10.2118/10691-MS.
Jensen, B.H. and Fredenslund, A. 1987. ASimplified Flash Procedure for Multicomponent Mixtures Containing Hydrocarbonsand One Non-Hydrocarbon Using Two-Parameter Cubic Equations of State.Ind. Eng. Chem. Res. 26 (10): 2129-2134.
Johns, R.T., Dindoruk, B., and Orr, F.M.Jr. 1993. Analytical Theory ofCombined Condensing/Vaporizing Gas Drives. SPE Advanced TechnologySeries 1 (2): 7-16. SPE-24112-PA. DOI:10.2118/24112-PA.
Johns, R.T., Fayers, F.J., and Orr, F.M.Jr. 1994. Effect of Gas Enrichmentand Dispersion on Nearly Miscible Displacements in Condensing/VaporizingDrives. SPE Advanced Technology Series 2 (2): 26-34.SPE-24938-PA. DOI: 10.2118/24938-PA.
Jutila, H.A., Logmo-Ngog, A.B., Sarkar,R., Killough, J.E., and Ross, F.C. 2001. Use of Parallel CompositionalSimulation to Investigate Productivity Improvement Options for aRetrograde-Gas-Condensate Field: A Case Study. Paper SPE 66397 presented atthe SPE Reservoir Simulation Symposium, Houston, 11-14 February. DOI:10.2118/66397-MS.
Liu, K. 2001. Reduce the Number of Components forCompositional Reservoir Simulation. Paper SPE 66363 presented at the SPEReservoir Simulation Symposium, Houston, 11-14 February. DOI:10.2118/66363-MS.
Michelsen, M. 1982a. The Isothermal FlashProblem. Part I: Stability. Fluid Phase Equilibria 9(1):1-19. DOI: http://dx.doi.org/10.1016/0378-3812(82)85001-2.
Michelsen, M. 1982b. The Isothermal FlashProblem. Part II: Phase Split Calculation. Fluid Phase Equilibria9 (1): 21-40. DOI: http://dx.doi.org/10.1016/0378-3812(82)85002-4.
Michelsen, M. 1986. Simplified FlashCalculation for Cubic Equations of State. Ind. Eng. Chem. Proc. Des.Dev. 25 (1): 184-188.
Nichita, D.V., Broseta, D., andHemptinne, J.C. 2004. MultiphaseEquilibrium Calculation Using Reduced Variables. Paper SPE 89439 presentedat the SPE/DOE Symposium on Improved Oil Recovery, Tulsa, 17-21 April. DOI:10.2118/89439-MS.
Pan, H. and Firoozabadi, A. 2001. Fast and Robust Algorithm forCompositional Modeling: Part II—Two-Phase Flash Computations. Paper SPE71603 presented at the SPE Annual Technical Conference and Exhibition, NewOrleans, 30 September-3 October. DOI: 10.2118/71603-MS.
Peng, D.Y. and Robinson, D.B. 1977. ARigorous Method for Predicting the Critical Properties of MulticomponentSystems from an Equation of State. AIChE J. 23 (2):137-144. DOI: http://dx.doi.org/10.1002/aic.690230202.
PVTSim R14. 2004. Calsep Inc., FluidHandling Software, Houston.
Rachford, H.H. Jr. and Rice, J.D. 1952.Procedure for Use of Electronic Digital Computers in Calculating FlashVaporization Hydrocarbon Equilibrium. Trans., AIME195:327-328.
Solano, R., Johns, R.T., and Lake, L.W.2001. Impact of Reservoir Mixingon Recovery in Enriched-Gas Drives Above the Minimum MiscibilityEnrichment. SPEREE 4 (5): 358-365. SPE-73829-PA. DOI:10.2118/73829-PA.
Stenby, E.H. and Wang, P. 1993. Noniterative Phase EquilibriumCalculation in Compositional Reservoir Simulation. Paper SPE 26641presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6October. DOI: 10.2118/26641-MS.
Varotsis, N. 1989. A Robust Prediction Method for RapidPhase-Behavior Calculations. SPERE 4 (2): 237-243.SPE-16943-PA. DOI: 10.2118/16943-PA.
Wilson, G.M. 1969. A ModifiedRedlich-Kwong Equation of State, Application to General Physical DataCalculation. Paper presented at the AIChE Natl. Meeting, Cleveland, Ohio, 4-7May.
Zick, A.A. 1986. A Combined Condensing/VaporizingMechanism in the Displacement of Oil by Enriched Gases. Paper SPE 15493presented at the SPE Annual Technical Conference and Exhibition, New Orleans,5-8 October. DOI: 10.2118/15493-MS.