Numerical Mechanistic Study of In-Situ CO2 EOR – Kinetics and Recovery Performance Analysis

The success of supercritical CO2 Enhanced Oil Recovery (EOR) cannot be duplicated if the cost of CO2 transposition and processing becomes prohibitive. Research results of the in-situ CO2 EOR (ICE) approach offered a potential technology for many waterflooded stripper wells that lack access to affordable CO2 sources. Previously the ICE synergetic mechanisms were only qualitatively attributed to oil swelling and viscosity reduction due to the preferential partition of CO2 into the oleic phase. This study aims to quantify the contributions to recovery factors from several plausible mechanisms with numerical modeling and simulation. First, the urea reaction was modeled as the CO2 generating chemical decomposing to CO2 and ammonia in the reservoir conditions. The CO2 partitions into oil, which leads to the reaction continuation to generate more CO2. The resulting ammonia largely left in water may further react with certain crudes to generate surfactants, thus, decrease the oil/water interfacial tension (IFT). It is expected that the oil containing CO2 also has a lower IFT with water. The reaction kinetics under different temperatures were incorporated into the numerical model. A numerical model featuring the synergetic mechanisms was built including stoichiometry and kinetics of urea reaction, oil swelling effect, oil viscosity reduction, and IFT reduction effect on the relative permeabilities. The laboratory experiments, pore volume injection versus oil saturation were history matched for three different oils including Dodecane, Earlsboro oil, and DeepStar oil. The phase behavior was modeled with the equation of state (EOS) under different mole fractions of CO2. The reaction kinetics were also modified to history match the laboratory experiment. The estimated reduction of oil viscosity was calculated, 76% for Earlsboro oil, 91% in DeepStar oil, and 75% in dodecane oil. The oil swelling factors ranged from 1.60% to 19% in the three lab models, which translates to the recovery factor of oil. The endpoints of relative permeability were modified to account for the recovery contribution to the IFT and viscosity reduction. The impact of reaction kinetics on oil swelling and recovery factor was also determined, and they are not numerically close to reaction kinetics used in the lab cases. The matched reaction kinetics, activation energy and reaction frequency factor for the dodecane laboratory experiment were 91.80 kJ/gmol and 6.5E+09 min−1. The study concluded that the incremental recovery due to oil swelling ranges between 3.16% and 18.30%, and then from 12.91% to 41.59% is due to IFT reduction for all the cases. The relative permeability and urea reaction kinetics remained the most uncertain parameters during history matching and modeling the ICE synergetic mechanisms.