A New Way of Compositional Simulation without Phase Labeling

SPE Journal ◽  
2021 ◽  
Vol 26 (02) ◽  
pp. 940-958
Author(s):  
Saeid Khorsandi ◽  
Liwei Li ◽  
Russell T. Johns

Summary Current relative permeability models rely on labeling a phase as “oil” and “gas” and cannot therefore capture accurately the effect of compositional variations on relative permeabilities and capillary pressures in enhanced oil recovery processes. Discontinuities in flux calculations caused by phase labeling problems not only cause serious convergence and stability problems but also affect the estimated recovery factor owing to incorrect phase mobilities. We developed a fully compositional simulation model using an equation of state (EoS) for relative permeabilities (kr) to eliminate the unphysical discontinuities in flux functions caused by phase labeling issues. The model can capture complex compositional and hysteresis effects for three-phase relative permeability. Each phase is modeled separately based on physical inputs that, in part, are proxies to composition. Phase flux calculations from one gridblock to another are also updated without phase labels. The tuned kr-EoS model and updated compositional simulator are demonstrated for simple ternary cases, multicycle three-phase water-alternating-gas (WAG) injection, and three-hydrocarbon-phase displacement with complex heterogeneity. The approach improves the initial estimates and convergence of flash calculations and stability analyses, as well as the convergence in the pressure solvers. The new compositional simulator allows for high-resolution simulation that gives improved accuracy in recovery estimates at significantly reduced computational time.

SPE Journal ◽  
2013 ◽  
Vol 18 (05) ◽  
pp. 841-850 ◽  
Author(s):  
H.. Shahverdi ◽  
M.. Sohrabi

Summary Water-alternating-gas (WAG) injection in waterflooded reservoirs can increase oil recovery and extend the life of these reservoirs. Reliable reservoir simulations are needed to predict the performance of WAG injection before field implementation. This requires accurate sets of relative permeability (kr) and capillary pressure (Pc) functions for each fluid phase, in a three-phase-flow regime. The WAG process also involves another major complication, hysteresis, which is caused by flow reversal happening during WAG injection. Hysteresis is one of the most important phenomena manipulating the performance of WAG injection, and hence, it has to be carefully accounted for. In this study, we have benefited from the results of a series of coreflood experiments that we have been performing since 1997 as a part of the Characterization of Three-Phase Flow and WAG Injection JIP (joint industry project) at Heriot-Watt University. In particular, we focus on a WAG experiment carried out on a water-wet core to obtain three-phase relative permeability values for oil, water, and gas. The relative permeabilities exhibit significant and irreversible hysteresis for oil, water, and gas. The observed hysteresis, which is a result of the cyclic injection of water and gas during WAG injection, is not predicted by the existing hysteresis models. We present a new three-phase relative permeability model coupled with hysteresis effects for the modeling of the observed cycle-dependent relative permeabilities taking place during WAG injection. The approach has been successfully tested and verified with measured three-phase relative permeability values obtained from a WAG experiment. In line with our laboratory observations, the new model predicts the reduction of the gas relative permeability during consecutive water-and-gas-injection cycles as well as the increase in oil relative permeability happening in consecutive water-injection cycles.


SPE Journal ◽  
2016 ◽  
Vol 21 (06) ◽  
pp. 1916-1929 ◽  
Author(s):  
Stefan Iglauer ◽  
Taufiq Rahman ◽  
Mohammad Sarmadivaleh ◽  
Adnan Al-Hinai ◽  
Martin A. Fernø ◽  
...  

Summary We imaged an intermediate-wet sandstone in three dimensions at high resolution (1–3.4 µm3) with X-ray microcomputed tomography (micro-CT) at various saturation states. Initially the core was at connate-water saturation and contained a large amount of oil (94%), which was produced by a waterflood [recovery factor Rf = 52% of original oil in place (OOIP)] or a direct gas flood (Rf = 66% of OOIP). Subsequent waterflooding and/or gasflooding (water-alternating-gas process) resulted in significant incremental-oil recovery (Rf = 71% of OOIP), whereas a substantial amount of gas could be stored (approximately 50%)—significantly more than in an analog water-wet plug. The oil- and gas-cluster-size distributions were measured and followed a power-law correlation N ∝ V−τ , where N is the frequency with which clusters of volume V are counted, and with decays exponents τ between 0.7 and 1.7. Furthermore, the cluster volume V plotted against cluster surface area A also correlated with a power-law correlation A ∝ Vp, and p was always ≈ 0.75. The measured τ- and p-values are significantly smaller than predicted by percolation theory, which predicts p ≈ 1 and τ = 2.189; this raises increasing doubts regarding the applicability of simple percolation models. In addition, we measured curvatures and capillary pressures of the oil and gas bubbles in situ, and analyzed the detailed pore-scale fluid configurations. The complex variations in fluid curvatures, capillary pressures, and the fluid/fluid or fluid/fluid/fluid pore-scale configurations (exact spatial locations also in relation to each other and the rock surface) are the origin of the well-known complexity of three-phase flow through rock.


SPE Journal ◽  
2010 ◽  
Vol 15 (03) ◽  
pp. 689-703 ◽  
Author(s):  
R.. Okuno ◽  
R.T.. T. Johns ◽  
K.. Sepehrnoori

Summary CO2 flooding at low temperatures often results in three or more hydrocarbon phases. Multiphase compositional simulation must simulate such gasfloods accurately. Drawbacks of modeling three hydrocarbon phases are the increased computational time and convergence problems associated with flash calculations. Use of a reduced method is a potential solution to these problems. We first demonstrate the importance of using three-phase flash calculations in compositional simulation by investigating difficulties with two-phase equilibrium approximations proposed in the literature. We then extend an algorithm for reduced two-phase flash calculations to three-phase calculations and show the efficiency and robustness of our algorithm. The reduced three-phase flash algorithm is implemented in a multiphase compositional simulator to demonstrate the speed-up and increased robustness of simulations in various case studies. Results show that use of a two-phase equilibrium approximation in reservoir simulation can result in a complete failure or erroneous simulation results. Simulation case studies show that our reduced method can decrease computational time significantly without loss of accuracy. Computational time is reduced using our reduced method because of the smaller number of equations to be solved and increased timestep sizes. We show that a failure of a flash calculation leads directly to reduced timestep sizes using the UTCOMP simulator.


2005 ◽  
Vol 8 (01) ◽  
pp. 33-43 ◽  
Author(s):  
Yildiray Cinar ◽  
Franklin M. Orr

Summary In this paper, we present results of an experimental investigation of the effects of variations in interfacial tension (IFT) on three-phase relative permeability. We report results that demonstrate the effect of low IFT between two of three phases on the three-phase relative permeabilities. To create three-phase systems in which IFT can be con-trolled systematically, we used a quaternary liquid system composed of hexadecane(C16), n-butanol (NBA), water (H2O), and isopropanol (IPA). Measured equilibrium phase compositions and IFTs are reported. The reported phase behavior of the quaternary system shows that the H2O-rich phase should represent the "gas" phase, the NBA-rich phase should represent the "oil" phase, and the C16-rich phase should represent the "aqueous" phase. Therefore, we used oil-wet Teflon (PTFE) bead packs to simulate the fluid flow in a water-wet oil reservoir. We determined phase saturations and three-phase relative permeabilities from recovery and pressure-drop data using an extension of the combined Welge/Johnson-Bossler-Naumann (JBN) method to three-phase flow. Measured three-phase relative permeabilities are reported. The experimental results indicate that the wetting-phase relative permeability was not affected by IFT variation, whereas the other two-phase relative permeabilities were clearly affected. As IFT decreases, the oil and gas phases become more mobile at the same phase saturations. For gas/oil IFTs in the range of 0.03 to 2.3 mN/m, we observed an approximately 10-fold increase in the oil and gas relative permeabilities against an approximately 100-fold decrease in the IFT. Introduction Variations in gas and oil relative permeabilities as a function of IFT are of particular importance in the area of compositional processes such as high-pressure gas injection, where oil and gas compositions can vary significantly both spatially and temporally. Because gas-injection processes routinely include three-phase flow (either because the reservoir has been water-flooded previously or because water is injected alternately with gas to improve overall reservoir sweep efficiency), the effect of IFT variations on three-phase relative permeabilities must be delineated if the performance of the gas-injection process is to be predicted accurately. The development of multicontact miscibility in a gas-injection process will create zones of low IFT between gas and oil phases in the presence of water. Although there have been studies of the effect of low IFT on two-phase relative permeability,1–14 there are limited experimental data published so far analyzing the effect of low IFT on three-phase relative permeabilities.15,16 Most authors have focused on the effect of IFT on oil and solvent relative permeabilities.17 Experimental results show that residual oil saturation and relative permeability are strongly affected by IFT, especially when the IFT is lower than approximately 0.1 mN/m (corresponding to a range of capillary number of 10–2 to 10–3). Bardon and Longeron3 observed that oil relative permeability increased linearly as IFT was reduced from approximately 12.5 mN/m to 0.04 mN/m and that for IFT below 0.04, the oil relative permeability curves shifted more rapidly with further reductions in IFT. Later, Asar and Handy6 showed that oil relative permeability curves began to shift as IFT was reduced below 0.18 mN/m for a gas/condensate system near the critical point. Delshad et al.15 presented experimental data for low-IFT three-phase relative permeabilities in Berea sandstone cores. They used a brine/oil/surfactant/alcohol mixture that included a microemulsion and excess oil and brine. The measurements were done at steady-state conditions with a constant capillary number of 10–2 between the microemulsion and other phases. The IFTs of microemulsion/oil and microemulsion/brine were low, whereas the IFT between oil and brine was high. They concluded that low-IFT three-phase relative permeabilities are functions of their own saturations only. Amin and Smith18 recently have published experimental data showing that the IFTs for each binary mixture of brine, oil, and gas phases vary as pressure increases(Fig. 1). Fig. 1 shows that the IFT of a gas/oil pair decreases as the pressure increases, whereas the IFTs of the gas/brine and oil/brine pairs approach each other.


1974 ◽  
Vol 14 (06) ◽  
pp. 573-592 ◽  
Author(s):  
K.H. Coats ◽  
W.D. George ◽  
Chieh Chu ◽  
B.E. Marcum

Coats, K.H., Member SPE-AIME, Intercomp Resource Development and Engineering, Inc., Houston, Texas George, W.D., Chu, Chieh, Member SPE-AIME, Getty Oil Co., Houston, Tx. Marcum, B.E., Member SPE-AIME, Getty Oil Co., Los Angeles, Calif. Abstract This paper describes a three-dimensional model for numerical simulation of steam injection processes. The model describes three-phase flow processes. The model describes three-phase flow of water, oil, and steam and heat flow in the reservoir and overburden. The method of solution simultaneously solves for the mass and energy balances and eliminates the need for iterating on the mass transfer (condensation) term.Laboratory data are reported for steamfloods of 5,780-cp oil in a 1/4 five-spot sand pack exhibiting three-dimensional flow effects. These experiments provide additional data for checking accuracy and provide additional data for checking accuracy and assumptions in numerical models. Comparisons of model results with several sets of experimental data indicate a need to account for effects of temperature on relative permeability. Calculated areal conformance of a steamflood in a confined five-spot depends strongly upon the alignment of the x-y grid axes relative to the diagonal joining injection and production wells. It has not been determined which, if either, of the two grid types yields the correct areal conformance.Model calculations indicate that steamflood pressure level strongly affects oil recovery. pressure level strongly affects oil recovery. Calculated oil recovery increases with decreasing pressure level. An example application illustrates pressure level. An example application illustrates the ability of the model formulation to efficiently simulate the single-well, cyclic steam stimulation problem. problem Introduction The literature includes many papers treating various aspects of oil recovery by steamflooding, hot waterflooding, and steam stimulation. The papers present laboratory experimental data, field papers present laboratory experimental data, field performance results, models for calculating fluid performance results, models for calculating fluid and heat flow, and experimental data regarding effects of temperature on relative permeability. The ultimate goal of all this work is a reliable engineering analysis to estimate oil recovery for a given mode of operation and to determine alternative operating conditions to maximize oil recovery.Toward that end, our study proposed to develop and validate an efficient, three-dimensional numerical model for simulating steamflooding, hot waterflooding, and steam stimulation. Laboratory steamflood experiments were conducted to provide additional data for validation. Desired model specifications included three-dimensional capability and greater efficiency than reported for previous models. Omitted from the specifications were temperature-dependent relative permeability and steam distillation effects.This paper describes the main features of the three-dimensional, steamflood model developed. Those features include a new method of solution that includes implicit water transmissibilities, that simultaneously solves for mass and energy balances, and that eliminates the need for iteration on the condensation term. Laboratory data are reported for steamfloods in a 1/4 five-spot model exhibiting three-dimensional flow effects. Numerical model applications described include comparisons with experimental data, a representative field-scale steamflood, and a cyclic steam stimulation example. REVIEW OF PREVIOUS WORK Early efforts in mathematical modeling of thermal methods concentrated on simulation of the heat flow and heat loss. Gottfried, in his analysis of in-situ combustion, initiated a series of models that solve fluid mass balances along with the energy balance. Davidson et al. presented an analysis for well performance during cyclic steam injection. Spillette and Nielsen treated hot waterflooding in two dimensions. Shutler described three-phase models for linears and two-dimensional steamflooding, and Abdalla and Coats treated a two-dimensional steamflood model using the IMPES method of solution. SPEJ P. 573


1999 ◽  
Vol 5 (4) ◽  
pp. 339-346 ◽  
Author(s):  
E. F. Balbinski ◽  
T. P. Fishlock ◽  
S. G. Goodyear ◽  
P. I. R. Jones

2017 ◽  
Vol 139 (3) ◽  
Author(s):  
Bander N. Al Ghamdi ◽  
Luis F. Ayala H.

Gas-condensate productivity is highly dependent on the thermodynamic behavior of the fluids-in-place. The condensation attendant with the depletion of gas-condensate reservoirs leads to a deficiency in the flow of fluids moving toward the production channels. The impairment is a result of condensate accumulation near the production channels in an immobility state until reaching a critical saturation point. Considering the flow phenomenon of gas-condensate reservoirs, tight formations can be inevitably complex hosting environments in which to achieve economical production. This work is aimed to assess the productivity gas-condensate reservoirs in a naturally fractured setting against the effect of capillary pressure and relative permeability constraints. The severity of condensate coating and magnitude of impairment was evaluated in a system with a permeability of 0.001 mD using an in-house compositional simulator. Several composition combinations were considered to portray mixtures ascending in complexity from light to heavy. The examination showed that thicker walls of condensate and greater impairment are attained with mixture containing higher nonvolatile concentrations. In addition, the influence of different capillary curves was insignificant to the overall behavior of fluids-in-place and movement within the pores medium. A greater impact on the transport of fluids was owed to relative permeability curves, which showed dependency on the extent of condensate content. Activating diffusion was found to diminish flow constraints due to the capturing of additional extractions that were not accounted for under Darcy's law alone.


SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1236-1253 ◽  
Author(s):  
Tae Wook Kim ◽  
E.. Vittoratos ◽  
A. R. Kovscek

Summary Recovery processes with a voidage-replacement ratio (VRR) (VRR = injected volume/produced volume) of unity rely solely on viscous forces to displace oil, whereas a VRR of zero relies on solution-gas drive. Activating a solution-gas-drive mechanism in combination with waterflooding with periods of VRR less than unity (VRR < 1) may be optimal for recovery. Laboratory evidence suggests that recovery for VRR < 1 is enhanced by emulsion flow and foamy (i.e., bubbly) crude oil at pressures under bubblepoint for some crude oils. This paper investigates the effect of VRR for two crude oils referred to as A1 (88 cp and 6.2 wt% asphaltene) and A2 (600 cp and 2.5 wt% asphaltene) in a sandpack system (18-in. length and 2-in. diameter). The crude oils are characterized with viscosity, asphaltene fraction, and acid/base numbers. A high-pressure experimental sandpack system (1 darcy and Swi = 0) was used to conduct experiments with VRRs of 1.0, 0.7, and 0 for both oils. During waterflood experiments, we controlled and monitored the rate of fluid injection and production to obtain well-characterized VRR. On the basis of the production ratio of fluids, the gas/oil and /water relative permeabilities were estimated under two-phase-flow conditions. For a VRR of zero, the gas relative permeability of both oils exhibited extremely low values (10−6−10−4) caused by internal gas drive. Waterfloods with VRR < 1 displayed encouraging recovery results. In particular, the final oil recovery with VRR = 0.7 [66.2% original oil in place (OOIP)] is more than 15% greater than that with VRR = 1 (55.6% OOIP) with A1 crude oil. Recovery for A2 with VRR = 0.7 (60.5% OOIP) was identical to the sum of oil recovery for solution-gas drive (19.1% OOIP) plus waterflooding (40.1% OOIP). An in-line viewing cell permitted inspection of produced fluid morphology. For A1 and VRR = 0.7, produced oil was emulsified, and gas was dispersed as bubbles, as expected for a foamy oil. For A2 and VRR < 1, foamy oil was not clearly observed in the viewing cell. In all cases, the water cut of VRR = 1 is clearly greater than that of VRR = 0.7. Finally, three-phase relative permeability was explored on the basis of the experimentally determined two-phase oil/water and liquid/gas relative permeability curves. Well-known algorithms for three-phase relative permeability, however, did not result in good history matches to the experimental data. Numerical simulations matched the experimental recovery vs. production time acceptably after modification of the measured krg and krow relationships. A concave shape for oil relative permeability that is suggestive of emulsified oil in situ was noted for both systems. The degree of agreement with experimental data is sensitive to the details of gas (gas/oil system) and oil (oil/water system) mobility.


Sign in / Sign up

Export Citation Format

Share Document