Abstract This study presents a new multipurpose compositional model for petroleum reservoirs. The proposed formulation allows the simulation of the isothermal and non-isothermal secondary or tertiary recovery processes, as well as condensate reservoirs where mass transfer of components between phases is considered to be significant. Although the current study limits the number of phases to four (4), the formulation presented in this study does not impose any limitation to the number of phases to be considered. The present compositional model utilizes either Soave–Redlich–Kwong (1972) or Peng and Robinson [Ind. Eng. Chem. Fundam. 15 (1976) 59; Can. J. Chem. Eng. (1976) 595] equations of state to describe the phase behavior of the system in lieu of predetermined equilibrium constants. The phase equilibrium package is capable of handling co-existing liquid phases and is equipped with a robust critical point determination technique, which is used in identifying the phases obtained under the existing temperature and pressure conditions. A linear programming algorithmic protocol utilized to determine the interaction coefficients and/or C7+ component(s) critical properties is also presented. The partial differential equations describing the transport phenomena in porous media are approximated using the finite-difference technique and are coupled with the flow equations. The resulting system of equations is solved simultaneously using the generalized Newton–Raphson procedure. The simulator has been tested using published data for several options. Results are found to be in good agreement. Also, various sensitivity checks were conducted by varying different parameters in order to test the performance of the model under different conditions. The fully implicit nature of the proposed model is found to be an advantage so that no stability problems have been encountered when large time-step sizes are utilized.
[1]
K. Aziz,et al.
Computation of Multiphase Equilibrium For Compositional Simulation
,
1982
.
[2]
K. H. Coats,et al.
Three-Dimensional Simulation of Steamflooding
,
1974
.
[3]
L. X. Nghiem,et al.
A Robust Iterative Method for Flash Calculations Using the Soave-Redlich-Kwong or the Peng-Robinson Equation of State
,
1983
.
[4]
K. Aziz,et al.
Compositional modeling with an equation of state
,
1981
.
[5]
R. Risnes,et al.
Equilibrium Calculations for Coexisting Liquid Phases
,
1984
.
[6]
D. E. Kenyon,et al.
Third SPE comparative solution project: Gas cycling of retrograde condensate reservoirs
,
1987
.
[7]
J. H. Abou-Kassem.
Investigation of grid orientation in a two-dimensional, compositional, three-phase steam model
,
1981
.
[8]
S. T. Obut,et al.
Determination of Interaction Coefficients and Critical Properties of Heavy Fractions for Equation of State Computations: An Automated Algorithmic Approach
,
1987
.
[9]
J. Douglas Faires,et al.
Numerical Analysis
,
1981
.
[10]
K. H. Coats,et al.
A New Technique for Determining Reservoir Description from Field Performance Data
,
1970
.
[11]
K. Aziz,et al.
Petroleum Reservoir Simulation
,
1979
.
[12]
Khalid Aziz,et al.
Fourth SPE comparative solution project: Comparison of steam injection simulators
,
1985
.
[13]
Abbas Firoozabadi,et al.
Predicting Phase Behavior of Condensate/Crude-Oil Systems Using Methane Interaction Coefficients
,
1978
.