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, urea reaction was modeled as the CO2 generating chemical decomposing to CO2 and ammonia under 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 oil components to generate surfactant, thus decreasing 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 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. I matched previous laboratory data for three different oils including dodecane, Earlsboro oil, and DeepStar oil, the same oils used in laboratory studies. The phase behavior was modeled with the Equation of State (EOS) under different mole fractions of CO2. The estimated reduction of oil viscosity was calculated, 79% for Earlsboro oil, 91% in DeepStar oil, and 76% in dodecane oil. The oil swelling factors ranged from 10% to 50% in the three lab models, which translate to the recovery factor of oil. Then I modified the endpoints of relative permeability 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 which were used in the lab cases. The study concluded that the incremental recovery due to oil swelling ranges between 6.4% and 18.0%, and the from 24% to 38% is due to IFT and viscosity 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. Later on, I upscaled the lab-scale model to a 3D sector model which features reaction kinetics, multi-components, and sensitivity study on selected parameters. To my knowledge, this is the first in modeling and differentiating the individual contributions on recovery from the synergetic mechanisms of ICE. The success of this model will tremendously reduce the potential laboratory experiment efforts and significantly improve the modeling capability in field application of the ICE technology. The developed model can be easily upscaled and retrofitted for other applications such as development planning and production forecast for ICE.