**Current Filters**

**Peer Reviewed**

**Source**

**Conference**

- 8th European Formation Damage Conference (1)
- Asia Pacific Oil and Gas Conference and Exhibition (1)
- Canadian Unconventional Resources and International Petroleum Conference (2)
- EUROPEC/EAGE Conference and Exhibition (4)
- SPE Annual Technical Conference and Exhibition (9)
- SPE Asia Pacific Improved Oil Recovery Conference (1)
- SPE European Formation Damage Conference (1)
- SPE Europec/EAGE Annual Conference (2)
- SPE EUROPEC/EAGE Annual Conference and Exhibition (3)
- SPE Improved Oil Recovery Symposium (3)
- SPE International Oil and Gas Conference and Exhibition in China (1)
- SPE International Symposium and Exhibiton on Formation Damage Control (1)
- SPE Latin American and Caribbean Petroleum Engineering Conference (1)
- SPE Reservoir Simulation Symposium (1)
- SPE Western Regional Meeting (2)
- SPE/DOE Enhanced Oil Recovery Symposium (1)
- SPE/EAGE Reservoir Characterization and Simulation Conference (1)

to

Go **Publisher**

**Theme**

**Author**

- Ameri Ghasrodashti, Amin (1)
- Andrianov, Alexey (1)
- Battistutta, Elisa (1)
- Betting, Marco (1)
- Biezen, E.N. (1)
- Botermans, C. Wouter (1)
- Bottero, Simona (1)
**Bruining, Johannes (47)**- Buijzen, Floris (1)
- de Ridder, Maarten (1)
- Delil, Harm Ad (1)
- Eftekhari, Ali Akbar (2)
- Elsinga, Gerrit (1)
- Enzien, Michael V. (1)
- Farajzadeh, Rouhollah (9)
- Geels, A.M. (1)
- Geels, Amanda M. (1)
- Godderij, R.R. (2)
- Godderij, R.R.G.G. (1)
- Gumrah, Fevzi (1)
- Heimovaara, Timo (1)
- Hopman, Victor (2)
- Kahrobaei, Siavash (1)
- Kalikmanov, Vitaly (1)
- Khosrokhavar, Roozbeh (2)
- Laros, Mario (1)
- Lomans, Bart (1)
- Mailybaev, Alexei (1)
- Mailybaev, Alexei A. (1)
- Marchesin, Dan (4)
- Mazumder, Saikat (3)
- Mojaddam, Ali (1)
- Molenaar, J. (1)
- Molenaar, Johannes (1)
- Muyzer, Gerard (1)
- Nguyen, Buu-Long (8)
- Palmgren, C.T.S. (2)
- Picioreanu, Cristian (1)
- Plug, Willem-Jan (1)
- Poelma, Christian (1)
- Ranganathan, Panneer (1)
- Rossen, William Richard (1)
- Rudolph-Floter, Susanne (1)
- Salimi, Hamidreza (11)
- Schoemaker, Christiaan (1)
- Slob, E.C. (2)
- Slob, Evert (1)
- Slob, Evert C. (5)
- Smeulders, David M. (1)
- Suicmez, Vural Sander (2)
- Thuss, Barbara (1)
- van Batenburg, Diederik W. (1)
- Van Der Kraan, Geert Marinus (1)
- van Hemert, Patrick (2)
- van Loosdrecht, Mark (1)
- Van Loosdrecht, Mark (1)
- Vermolen, Fred (1)
- Wahanik, Helmut (2)
- Wolf, Karl-Heinz (1)
- Wolf, Karl-Heinz A.A. (6)
- Zanganeh, Maryam Namdar (1)
- Zitha, Pacelli Lidio Jose (5)

to

Go **Concept Tag**

- Behavior (4)
- breakthrough (3)
- carbon dioxide (5)
- case (4)
- Case II (2)
- cell (6)
- coal (3)
- coal bed methane (4)
- coal seam gas (4)
- coalbed methane (4)
- coefficient (9)
- complex reservoir (11)
- component (3)
- concentration (11)
- condition (5)
- core (4)
- Csug Spe (2)
- cumulative oil production (4)
- density-driven natural convection (3)
- dielectric permittivity (3)
- dielectric property (4)
- diffusion (6)
- dissolution (3)
- drainage (3)
- Eclipse (2)
- effect (11)
- enhanced recovery (24)
- equation (18)
- equation of state (3)
- equilibrium (3)
- experiment (16)
- Farajzadeh (3)
- field (4)
- flow in porous media (26)
- Fluid Dynamics (26)
- fluid modeling (3)
- formation evaluation (8)
- fraction (6)
- fracture (10)
- fracture system (3)
- frequency (6)
- function (6)
- gas saturation (3)
- geothermal reservoir (4)
- GPR (2)
- gravity (9)
- health safety security environment and social responsibility (3)
- heat (4)
- heat loss (3)
- homogenization (4)
- HTO (2)
- imbibition (6)
- increase (4)
- Inj (2)
- injection (12)
- interface (9)
- mass transfer (3)
- matrix (9)
- matrix block (5)
- model (16)
- mole fraction (4)
- natural convection (3)
- paper (4)
- permeability (10)
- polymer (3)
- process (6)
- production (11)
- profile (3)
- Rayleigh (3)
- reflection (4)
- regime (4)
- region (4)
- Reservoir Characterization (5)
- reservoir description and dynamics (47)
- reservoir simulation (10)
- SAGD (6)
- sample (5)
- sample cell (3)
- saturation (14)
- scaling method (3)
- scatter function (4)
- shock (3)
- signal (3)
- Simulation (5)
- society of petroleum engineers (5)
- sorption (4)
- SPE (10)
- steam (3)
- steam-assisted gravity drainage (6)
- structure (3)
- subsurface storage (7)
- system (6)
- thermal methods (6)
- transfer (5)
- transmission (3)
- Upstream Oil & Gas (47)
- viscosity (3)
- water (27)
- water saturation (4)
- Wave (5)

to

Go

In a transition period from a fossil fuel based society to a sustainable energy society it is expected that CO2 capture and subsequent sequestration (CCS) in geological formations will play a major role in reducing greenhouse gas emissions. Possibilities of sequestration include storage in aquifers and depleted gas reservoir. The storage capacity of gas reservoirs for CO2 depends also on the sorption in the omnipresent minerals and shales. It is important to investigate whether adsorption on shales gives an important contribution to the storage capacity. It is also important to relate the adsorption to the carbon content in the shale. Only a few measurements have been reported in the literature for high-pressure gas sorption on shales, and interest is largely focused on shales occurring outside Europe. We present results using a high pressure manometric set-up on a dried black shale sample from Belgium. It consists of more than 57% of clay minerals and 6.58% organic matter. The excess sorption isotherm shows an initial increase to a maximum value of 0.19 mmol/gram and then starts to decrease until it becomes zero at 82 bar and subsequently the excess sorption becomes negative. Similar behavior was also observed for other shales and coal reported in the literature.

We derive the equation for excess sorption in the manometric set-up allowing for a changing void volume. This equation is based on the finite density of the adsorbed phase. However, this is not the only mechanism causing a maximum in the sorption curve. Other reasons for void volume change are swelling of the shale and volume changes due to chemical reactions excluding sorption. Further research is necessary to investigate reasons for void volume changes in shales.

adsorption, air emission, capacity, carbon dioxide, cell, climate change, CO2 capture, complex reservoir, determination, excess sorption, experiment, health safety security environment and social responsibility, helium, manometric setup, reservoir description and dynamics, sample, sample cell, shale gas, sorption, sorption capacity, SPE, subsurface storage, Upstream Oil & Gas, void

Oilfield Places:

- North America > United States > West Virginia > Big Sandy Field (0.99)
- North America > United States > Kentucky > Big Sandy Field (0.99)
- Europe > Germany > Saar Basin (0.98)
- (2 more...)

SPE Disciplines:

- Health, Safety, Security, Environment and Social Responsibility > Environment > Climate change (1.00)
- Facilities Design, Construction and Operation > Unconventional Production Facilities > CO2 capture and management (1.00)
- Reservoir Description and Dynamics > Storage Reservoir Engineering > CO2 capture and sequestration (0.93)
- (2 more...)

We quantify the capillary pressure effect on the phase equilibrium of the CO2-water system. Our interest is in the capillarypressure range between 0 and 100 bars for temperatures between 293 and 372 K and bulk (wetting-phase) pressures between 25 and 255 bars. For this purpose, we have implemented the capillary pressure effect in the PRSV equation of state. Inclusion of capillary pressure in the phase equilibrium of the CO2-water system makes it possible to determine the capillary-pressure effect on the CO2 storage capacity and heat-energy recovery for CO2-water injection into geothermal reservoirs. We illustrate the process using a 2D model of the geothermal reservoir in the Delft Sandstone Member, below the city of Delft (The Netherlands). The process involves phase transitions between single-phase and two-phase regions. To deal with phase appearance and disappearance, we have applied a new and effective solution approach, the so-called "nonisothermal negative saturation?? (NegSat) solution approach.

The results show that the capillary pressure promotes evaporation. In the pressure and temperature range of our interest, capillary pressure reduces the CO2 solubility in water and the aqueous-phase density up to 64% and 1.3%, respectively, whereas it increases the water solubility in the CO2-rich phase and the CO2-rich-phase density up to 3,945% (40.5 times) and 1,544%, respectively. Capillary pressure shifts the CO2 liquid-vapor transition and consequently the upper critical point of the CO2-water system to a lower pressure. The intensity of the shift depends on the value of the capillary pressure and the bulk (wetting-phase) pressure. For instance, the CO2 liquid-vapor transition at T = 293 K occurs approximately at 60 bars for Pc = 0 bars, whereas it occurs at 15 bars for Pc = 45 bars.

For mixed CO2-water injection into the geothermal reservoir (200 bars < P < 260 bars, 290 K < T < 360 K), inclusion of the capillary pressure effect in the phase-equilibrium behavior does not significantly alter the capillary CO2-trapping mechanism. In other words, CO2 banks are mainly formed in the highly permeable zones that are surrounded by less permeable zones. However, for injected CO2 concentrations close to the bubble point, the effect of capillary pressure on the

phase equilibrium reduces the heat recovery by 37% and the CO2-storage capacity also by 37%. For overall injected CO2 mole fractions between 4% and 13%, the reduction in the heat recovery and CO2-storage capacity is 10%. Based on simulations, we construct a plot of the recuperated heat energy versus the maximally stored CO2 for a variety of conditions; we compare the results including and excluding the effect of capillary pressure in the phase-equilibrium calculations.

Salimi, Hamidreza (Delft University of Technology) | Wolf, Karl-Heinz (Delft University of Technology) | Bruining, Johannes (Delft University of Technology)

Cold mixed CO_{2}/water injection into hot-water reservoirs can be used for simultaneous geothermal-energy (heat) production and subsurface CO_{2} storage. This paper studies this process in a 2D geothermal homogeneous reservoir, a layered reservoir, and a heterogeneous reservoir represented by a stochastic-random field. We give a set of simulations for a variety of CO_{2}/water-injection ratios. In this process, often regions of two-phase flow are connected to regions of single-phase flow. Different systems of equations apply for single-phase and two-phase regions. We develop a solution approach, called the nonisothermal-negative-saturation (NegSat) solution approach, to solve efficiently nonisothermal compositional flow problems (e.g., CO_{2}/water injection into geothermal reservoirs) that involve phase appearance, phase disappearance, and phase transitions. The advantage of this solution approach is that it circumvents using different equations for single-phase and two-phase regions and the ensuing unstable switching procedure. In the NegSat approach, a single-phase multicomponent fluid is replaced by an equivalent fictitious two-phase fluid with specific properties. The equivalent properties are such that the extended saturation of a fictitious gas is negative in the single-phase aqueous region.

We discuss the salient features of the simulations in detail. When two phases are present at the injection side, heterogeneity and layering lead to more CO_{2} storage compared with the homogeneous case because of capillary trapping. In addition, layering avoids movement of the CO_{2} to the upper part of the reservoir and thus reduces the risk of leakage. Our results also show that heterogeneity and layering change the character of the solution in terms of useful-energy production and CO_{2} storage. The simulations can be used to construct a plot of the recovered useful energy vs. maximally stored CO_{2}. Increasing the amount of CO_{2} in the injection mixture leads to bifurcation points at which the character of the solution in terms of energy production and CO_{2} storage changes. For overall injected-CO_{2} mole fractions less than 0.04, the result with gravity is the same as the result without gravity. For larger overall injected-CO_{2} mole fractions, however, the plot without gravity differs from the plot with gravity because of early breakthrough of a supercritical-CO_{2} tongue near the caprock. The plot of the useful energy (exergy) vs. the CO_{2}-storage capacity in the presence of gravity shows a Z-shape. The top horizontal part represents a branch of high exergy recovery and a relatively lower storage capacity, whereas the bottom part represents a branch of lower exergy recovery and a higher storage capacity.

enhanced recovery, fraction, gas saturation, geothermal reservoir, gravity, homogeneous case, homogeneous-permeability field, increase, injection, mole fraction, production, region, reservoir, reservoir description and dynamics, saturation, single-phase region, two-phase region, Upstream Oil & Gas, water

Solvent injection has been considered as an efficient method for enhancing oil recovery from fractured reservoirs. If the mass transfer would be solely based on diffusion, oil recovery would be unacceptably slow. The success of this method therefore depends on the degree of enhancement of the mass exchange rate between the solvent residing in the fracture and the oil residing in the matrix.

A series of soak experiments have been conducted to investigate the mass transfer rate between the fracture and the matrix. In a soak experiment, a porous medium containing oil is immersed in an open space containing the solvent to simulate the matrix and the fracture respectively. We use a CT scanner to visualize the process. The experimental data are compared with a simulation model that takes diffusive and gravitational forces into account.

We find that the initial stage of all experiments can be described by a diffusion-based model with an enhanced "effective diffusion coefficient??. In the second stage enhancement of the transfer rate occurs due to the natural convection of solvent in the fracture. The experiments are quantitatively modeled by numerical simulations. We find that transfer rates depend on the properties of the rock permeability, the viscosity and the density of solvent and oil. The gravity enhanced transfer is quantified by comparison of experimental and simulated results.

bottom, coefficient, concentration, core, decalin, diffusion, enhanced recovery, experiment, flow in porous media, Fluid Dynamics, fracture, gravity, hexadecane, injection, matrix, mechanism, pentane, process, reservoir description and dynamics, Simulation, society of petroleum engineers, transfer, Upstream Oil & Gas, viscosity

SPE Disciplines:

Ameri Ghasrodashti, Amin (Tarbiat Modares University) | Farajzadeh, Rouhollah (Shell Intl. E&P BV) | Suicmez, Vural Sander (Brunei Shell Petr Sdn Bhd.) | Bruining, Johannes (Delft U. of Technology)

Gas oil gravity drainage is an effective oil recovery process, which has been proven in the field. Under favorable conditions the displacement is stable and for the right surface tension combinations the residual oil saturation is low. In the absence of gas dissolution, the recovery after gas injection is usually low as a large amount of oil remains capillary trapped in the matrix blocks. However, when the main gas constitutes is soluble in the oil, the dissolution leads to mixing and interfacial tension (IFT) reduction, which cause gravity enhanced transfer between matrix and fracture. Therefore, a study of the mechanisms that control the interactions between fracture and matrix (e.g. capillarity, gravity, phase behavior and flow behavior) can help to optimize recovery. This paper concerns an experimental study to investigate whether gravity drainage is also an effective recovery process in fractured reservoirs. In this study, we describe six gas injection experiments conducted at different miscibility conditions, i.e., immiscible, developed miscible and first contact miscible (FCM), using CO2, nitrogen and flue gas.

In addition, the impact of switching from an immiscible (Nitrogen, Flue gas) injection gas to non-equilibrium and fully miscible CO2 injection is investigated. In one of the experiments, we study the effect of a permeability barrier on the recovery efficiency from the matrix block when CO2 is injected in the fracture at immiscible and miscible conditions. Accurate modeling for the transfer between fracture and matrix is also essential for accurate recovery predictions. In this study, a numerical model is developed to perform compositional simulations of gas injection for different miscibility scenarios. Results revealed that ultimate oil recovery increases considerably once miscibility is reached. Miscibility can usually be achieved at high pressures only. High pressure gas injection has two disadvantages, viz., (1) one may need a larger mass of gas to fill the pore space from where the oil is recovered and (2) the density of injected gas increases significantly, which reduces the density difference between the gas and oil. This leads to less effective gravity mediated recovery. Even if the impermeable layer

impairs the performance of the gas oil gravity drainage (GOGD) process for immiscible gas injection, it improves the recoveries for first contact miscible gas injection.

SPE Disciplines:

This paper gives an analysis of the Thomas and Windle model (Thomas and Windle 1982) to determine its usefulness for describing anomalous diffusion of CO_{2} in coal and its relation to matrix swelling. In addition, a finite-element description for this model is derived. For reasons of easy reference, a shortened derivation of the Thomas and Windle model is presented, which was originally derived to describe diffusion in polymers. proposed by Hui et al. (1987a, 1987b). Because the cumulative sorption showed *t ^{a}* behavior with

Industry:

- Energy > Oil & Gas > Upstream (1.00)
- Materials > Chemicals > Commodity Chemicals > Petrochemicals (0.46)

SPE Disciplines:

There is a renewed interest in using combustion to recover medium- or high-viscosity oil. Despite numerous experimental, numerical, and analytical studies, the mechanisms for incomplete fuel combustion or oxygen consumption are not fully understood. Incomplete oxygen consumption may lead to low-temperature oxidation reactions downstream. This paper shows that these features emerge in a relatively simple 1D model, where air is injected in a porous medium filled with inert gas, water, and an oil mixture consisting of precoke, medium oil, and light oil. Precoke is a component that is dissolved in the oil but has essentially the same composition as coke. At high temperatures, precoke is converted to coke, which participates in high-temperature oxidation. At high temperatures, medium-oil components are cracked, releasing gaseous oil. Light-oil components and water are vaporized. The model possesses an analytical solution, which was obtained by a concept introduced by Zeldovich et al. (1985). This concept, which underlies most analytical approaches such as the reaction-sheet approximation and large-activation-energy asymptotics, entails that reaction can occur only in a very small temperature range because of the highly nonlinear nature of the Arrhenius factor. For a temperature below this range, the reaction rate is too slow, and for temperatures above this range, the reaction rate is so fast that either the fuel or oxygen concentrations become zero. The model results, in the absence of external heat losses, show that there are two combustion regimes in which coke or oxygen is partially consumed. In one regime, the reaction zone moves in front of the heat wave; whereas, in the other regime, the order of the waves is reversed. There are also two combustion regimes in which the coke and oxygen are completely consumed. Also, here the reaction zone can move in front of or at the back of the heat wave. Each combustion regime is described by a sequence of waves; we derive formulas for parameters in these waves. We analyze our formulas for typical in-situ-combustion data and compare the results with numerical simulation. The main conclusion is that mainly two key parameters (i.e., the injected oxygen mole fraction and the fuel concentration) determine the combustion-front structure and when either incomplete oxygen consumption or incomplete fuel consumption occurs in the high-temperature oxidation zone.

coke, combustion, concentration, consumption, enhanced recovery, heat, HTO, HTO reaction, hto region, in-situ combustion, Inj, occur, oxygen, reaction, reaction wave, reaction-leading structure, regime, region, reservoir description and dynamics, SAGD, steam-assisted gravity drainage, thermal methods, Upstream Oil & Gas, vaporization, Wave

SPE Disciplines: Reservoir Description and Dynamics > Improved and Enhanced Recovery > Thermal methods (1.00)

We use upscaling through homogenization to predict oil recovery from fractured reservoirs consisting of matrix columns, also called vertically fractured reservoirs (VFRs), for a variety of conditions. The upscaled VFR model overcomes limitations of the dual-porosity model, including the use of a shape factor. The purpose of this paper is to investigate three main physical aspects of multiphase flow in fractured reservoirs: reservoir wettability, viscosity ratio, and heterogeneity in rock/fluid properties. The main characteristic that determines reservoir behavior is the Péclet number that expresses the ratio of the average imbibition time divided by the residence time of the fluids in the fractures. The second characteristic dimensionless number is the gravity number.

Upscaled VFR simulations, aimed at studying the mentioned features, add new insights. First, we discuss the results at low Péclet numbers. For only small gravity numbers, the effect of contact angle, delay time for the nonequilibrium capillary effect, the heterogeneity of the matrix-column size, and the matrix permeability can be ignored without appreciable loss of accuracy. The ultimate oil recovery for mixed-wet VFRs is approximately equal to the Amott index, and the oil production does not depend on the absolute value of the phase viscosity but on viscosity ratio. However, large gravity numbers enhance underriding, aggravated by large contact angles, longer delay times, and higher viscosity ratios. Layering can lead to an improvement or deterioration, depending on the fracture aperture and permeability distribution. At low Péclet numbers, the fractured reservoir behaves very similarly to a conventional reservoir and depends largely on the viscosity ratio and the gravity number. At high Péclet numbers, after water breakthrough, the oil recovery appears to be proportional to the cosine of the contact angle and inversely proportional to the sum of the oil and water viscosity. In addition, the mixed-wetting effect is more pronounced; there are significant influences of delay time (nonequilibrium effects), matrix permeability, matrix-column size, and the column-size distribution on oil recovery. At low gravity numbers and an effective length/thickness ratio larger than 10, the oil recovery is independent of the vertical-fracture-aperture distribution. For the same amount of injected water, the recovery at low Péclet numbers is larger than the recovery at high Péclet numbers.

case, complex reservoir, cumulative oil production, effect, enhanced recovery, flow in porous media, Fluid Dynamics, fracture, fractured reservoir, gravity, imbibition, lateral matrix-column, matrix, matrix column, matrix permeability, Oil Recovery, permeability, production, recovery, reservoir description and dynamics, reservoir simulation, Upstream Oil & Gas, viscosity, water

SPE Disciplines:

- Reservoir Description and Dynamics > Unconventional and Complex Reservoirs > Naturally-fractured reservoirs (1.00)
- Reservoir Description and Dynamics > Reservoir Simulation (1.00)
- Reservoir Description and Dynamics > Reservoir Fluid Dynamics > Flow in porous media (1.00)
- Reservoir Description and Dynamics > Improved and Enhanced Recovery (1.00)

Flow modeling in fractured reservoirs is largely confined to the so-called sugar-cube model. Here, we consider a situation where matrix blocks are connected to neighboring blocks so that part of the global flow occurs only in the matrix domain. We call this a partially fractured reservoir (PFR). As opposed to the sugar-cube model, global flow in the matrix blocks plays an important role in the PFR when the interconnections between the matrix blocks are sufficiently large. We apply homogenization to derive an upscaled model for PFRs that combines dual-porosity and dual-permeability concepts simultaneously. We formulate a well-posed fully implicit 3D upscaled numerical model and investigate oil-recovery mechanisms for different dimensionless characteristic numbers. As we found previously for the sugar-cube model, the Péclet number, defined here as the ratio of the capillary diffusion time in the matrix to the residence time of the fluids in the fracture, plays a crucial role. The gravity number and specific fracture/matrix-interface area play a secondary role. For low Péclet numbers and high gravity numbers, the results are sensitive to gravity and water-injection rates, but relatively insensitive to the specific fracture/matrix-interface area, matrix-block size, and reservoir geometry (i.e., sugar cube vs. PFR). At low Péclet numbers and high gravity numbers, ECLIPSE simulations using the Barenblatt or Warren and Root (BWR) approach give poor predictions and overestimate the oil recovery, but, at short injection times, show good agreement with the solution of the PFR model at intermediate Péclet numbers. At high Péclet numbers, the results are relatively insensitive to gravity, but sensitive to the other conditions mentioned. In particular, when the specific fracture/matrix-interface area is large, it enhances the imbibition and, consequently, leads to a higher oil production. If this specific interface area is small, it leads to a considerable retardation of the imbibition process, which leads to an earlier water breakthrough and lower oil recovery. The BWR (commercial simulator) simulations and the sugar-cube model result in inaccurate predictions of the oil-production rate at high Péclet numbers. This can be inferred from the discrepancy with respect to the PFR model for which we assert that it accurately predicts the oil recovery. We conclude that, at low Péclet numbers and large gravity numbers, it is advantageous to increase the water-injection rate to improve the net present value. However, at high Péclet numbers, increasing the flow rate may lead to uneconomical water cuts.

cumulative oil production, enhanced recovery, equation, flow in porous media, Fluid Dynamics, fracture, gravity, homogenization, hydraulic fracturing, imbibition, matrix, matrix block, model, permeability, PFR model, production, reservoir description and dynamics, reservoir simulation, specific fracture matrix-interface, system, Upstream Oil & Gas, water, well completion

SPE Disciplines:

In this paper we follow a similar procedure as proposed by Koval (1963) to analytically model the performance of gravitationally unstable flow in porous media. The Koval model is analogous to the Buckley-Leverett method and multiplies the heterogeneity index of the system as an input (H-factor) with the fluid-flow (here gravity) induced instability factor, E to obtain the Koval factor KG = HE. This paper only considers the gravity induced instability factor E (H=1). The Koval factor is implemented in a modified fractional flow function that includes a dilution effect when the CO2 moves away from the interface to describe countercurrent gravity flow. The pseudo two-phase flow problem provides the average concentration of CO2 in the brine as a function of distance. The KG-factor can be used in commercial simulators to account for the density-driven natural convection, which cannot be currently captured because the grid cells are typically orders of magnitude larger than the wavelength of the initial fingers. Such natural convection effects occur in storage of greenhouse gases in aquifers and EOR processes using carbon dioxide or other solvents.

A comparison of the analytical model with the horizontally-averaged concentrations obtained from 2-D numerical simulations provides a correlation for calculation of the KG-factor for different Rayleigh numbers. The model shows a rarefaction followed by shock-like behavior because the CO2 concentration decreases away from the gaseous CO2-liquid interface. The agreement between the analytical model and full numerical simulation is practically acceptable. We leave the introduction of the heterogeneity factor for future work.

analytical model, concentration, condition, density-driven natural convection, diffusion, enhanced recovery, Farajzadeh, flow in porous media, Fluid Dynamics, heterogeneity, inner solution, instability, Koval, mixture, model, rarefaction, Rayleigh, reservoir description and dynamics, reservoir simulation, saturation, solution, SPE, Upstream Oil & Gas, viscosity, water

SPE Disciplines:

Thank you!