Foam has been used as a mobility control technique in CO₂ flooding to improve volumetric sweep efficiency. Stabilizing CO₂ foam with nanoparticle instead of surfactant has some notable advantages. Nanoparticle-stabilized foam is very stable because a large adsorption energy is required to bring nanoparticles to the bubble interfaces. As a solid, nanoparticle can potentially withstand the high temperature in the reservoir, providing a robust foam stability for an extended period of time. The ability of nanoparticles to generate foam only above a threshold shear rate is promising as foam can be engineered to form only in the high permeability zone. These nanoparticles are hundreds of times smaller than pore throats and thus can travel in the reservoir without plugging the pore throats. Surface-modified silica nanoparticle was found to stabilize CO₂ -in-water foam at temperature up to 80 ˚C and salinity as high as 7.2 wt%. The foam was generated through the co-injection of aqueous nanoparticle dispersion and CO₂ into consolidated rock cores, primarily sandstones, with and without an induced fracture in the core. A critical shear rate for foam generation was found to exist in both matrix and fracture, however, this critical rate varied with the experiment conditions. The effects of experimental parameters on the critical shear rate and foam apparent viscosity were also investigated. Additionally, the flow distribution calculation in fractured sandstone cores revealed a diversion of flow from fracture toward matrix once foam was generated, suggesting conformance control potential in fractured reservoirs. In order to study foam rheology, high-permeability beadpack was installed upstream of the core to serve as a foam generator. This allows the foam mobility to be measured solely while being transported through the core, without the complicating effect of transient foam generation in the core. The injection of the pre-generated foam into the core at residual oil condition was found to reduce the residual oil saturation to the same level as CO₂ flood, however, with the advantage of mobility control. The 'coalescence-regeneration' mechanism of foam transport in porous media possibly allowed the foam's CO₂ to contact and mobilize the residual oil. The injection of the foam slug followed by a slug of only CO₂ was also tested, showing similar viscosification as the continuous foam injection, however, required less nanoparticles. / text

Surfactant-polymer (SP) and alkali-surfactant-polymer (ASP) flooding has gained little traction among different tertiary recovery strategies such as thermal and miscible gas flooding; however, many mature onshore reservoirs could be potential candidates. More than four decades of research has detailed technical challenges and successes through laboratory experimentation, chemical flood simulation, and some pilot projects, which have provided technical screening procedures to efficiently filter unfeasible projects. Therefore, technical understanding seems sufficient to advance projects through early development stages; however, a project value identification and realization process ultimately dictates project implementation in the oil and gas industry, with technical feasibility merely supporting overall valuation and project feasibility. A quick early screening methodology integrating important project valuation criteria can efficiently assess large numbers of projects. The relatively few studies detailing chemical flooding valuation from just an economic standpoint reflects the need for an integrated process-oriented framework for quick early screening valuation of chemical flooding opportunities. This study develops an integrated process-oriented framework for early screening and valuation, with an overall objective to quickly filter unfeasible projects based on valuation criteria, rather than technical feasibility alone. A reservoir-to-market model was developed, integrating information from laboratory experiments (phase behavior, core flood), field analogues (well performance and layout), facilities, rigs, costs, scheduling, and economics. Recently published ASP flood data of the central Xing2 area in Daqing, China was used for model inputs. A reservoir-to-market benchmark model for a typical mature onshore field was successfully built and tested, and could value projects using standard economic metrics (net present value, internal rate of return, value investment ratio, unit technical cost, and payback period). Model simplification was achieved through global sensitivity analysis. Using a mean-reversion oil price model, the oil price accounted for 98% of the total sensitivity. . Model efficiency was achieved through discretization of input parameter uncertainties, which sped the screening process. Decision-making between model alternatives given information and different states of nature was performed through decision-tree techniques based on overall project valuation. Overall, this study was novel and provided benefit as a robust, integrated process-oriented framework for chemical EOR project screening, valuation, and decision-making. / text

The efficiency of one fluid displacing another in permeable media depends greatly on the pore-scale dynamics at the main wetting front. Experiments have shown that the frontal dynamics can result in two different flow regimes: a stable and an unstable front. In stable displacements, any perturbation of the front will diminish with time and the effect of variation in permeability will be lessened. In contrast, in unstable displacements any perturbation of the front will grow with time and any variation in permeability will be magnified. In this dissertation, the stability of two different displacement processes are contemplated; a) vertical infiltration of dense liquid into dry sand from above and b) horizontal displacement of nanoparticle suspension with high pressure liquid CO₂. Significant insights are obtained by measuring the in-situ flow patterns in real time with a light transmission method and CT scanning. Vertical infiltration of dense fluid into dry sands from above is often observed to be unstable and produce gravity driven fingers. The formation of gravity fingers can have large consequences on the sweep efficiency of a displacement. Infiltration experiments showed that gravity driven fingers have a unique saturation profile known as saturation overshoot with a higher saturation at the finger tips than the saturation at the finger tail. Despite the vast number of theoretical and experimental investigations, conditions under which the front is unstable, remain unclear. To determine what controls the saturation overshoot and how it relates to the dynamics at the initial wetting front, saturation overshoot was measured as a function of flux for seven different liquids. These liquids gave a range of molecular weights, viscosities, and vapor pressures. It is found that for each fluid there is a flux (called overshoot flux) below which saturation overshoot ceases and the front is diffuse. The magnitude of the overshoot flux depends inversely on the invading fluid's viscosity and shows little or no dependence on the invading fluid's surface tension, vapor pressure, or miscibility with water. Since the saturation overshoot is not described by the continuum multi-phase flow models, the experimental results are used to develop a semi-continuum model that bridges the continuum-scale and pore-scale physics. The proposed model predicts the observed dependence of overshoot on media permeability and invading fluid properties. At the planned depth for CO₂ injection, either as an enhanced oil recovery technique or for CO₂ storage, CO₂ is typically less dense and less viscous than the in-situ fluid. Therefore, CO₂ injection is unstable and produces viscous fingers. This can greatly reduce the efficiency of a CO₂ flood or CO₂ storage capacity of an aquifer. To remedy this behavior, surface treated nanoparticles were used to reduce the mobility of injected CO₂. Displacement experiments were performed at low pressure with a CO₂ analogue (n-octane) fluid and at high pressure with liquid CO₂. Saturation distributions and pressure drops were measured in real time with the CT scanner when high pressure liquid CO₂ or n-octane was used to displace brine in different cores with and without suspended nanoparticles. In the presence of nanoparticles, the displacement front is more spatially uniform with a later breakthrough compared to the same experiment with no suspended nanoparticles. These observations suggest that nanoparticle stabilized foam, which forms during the displacement, acts to suppress the instability. It is argued that the generation of droplets occurs at the leading front of all drainage displacements. In the presence of nanoparticles, these droplets are preserved when nanoparticle adhere at the fluid-fluid interface. The new mechanism for foam generation described here, provides an interesting alternative for mobility control in CO₂ floods. Moreover, the same mechanism can potentially a) increase the CO₂ storage capacity of an aquifer, b) enhance the CO₂ capillary trapping, and c) provide an engineered barrier to CO₂ leakage from a storage sites, thereby alleviating the risk of contaminating the overlying fresh groundwater resources for CO₂ storage projects. / text

Operators can improve a reservoir’s value by optimizing it in a more holistic manner, or over its entire life cycle. This thesis developed approaches to life cycle optimization, with emphasis on accessible technical and economic modeling techniques for production. The challenges of life cycle optimization are properly scheduling the times at which the operator should switch from one recovery phase to the next, along with determining other field design parameters such as well spacing and injection pressures for waterflooding and enhanced oil recovery processes. To deliver the most value, the operator needs to produce from a reservoir the greatest quantity of oil, at a relatively low cost, reasonably soon, and ideally at a time when the oil price is high. This is quite a tall order, as these goals are often in conflict. This thesis extended existing research regarding lifecycle optimization, first modeling production from a reservoir using an exponential decline model and assuming the oil price’s behavior can be approximated with mean-reverting processes. Implications of operating and capital costs potentially being correlated with the oil price were also examined. Finally, a mean-reverting price model that forecasts the mean oil price as increasing and described by a logistic model was proposed to accommodate both recent price forecasts and economic reality. As exponential decline models are more appropriate for characterizing existing production history rather than making a priori predictions, a geologic-parameter-based model was developed using a tank model for primary recovery and a model based on Koval theory and parameterizing a reservoir in terms of flow capacity and storage capacity for waterflooding and CO2 flooding. This model was adapted from existing theory to account for situations where a waterflood has incompletely swept a reservoir at the start of CO2 flooding. Analytical expressions were also derived for estimating injection rates into a formation parameterized by flow capacity and storage capacity. The geologic-parameter-based model was combined with economic assumptions and optimized using a genetic algorithm. This optimization suggested an operator should switch from primary recovery to a CO2 flood with a large WAG ratio relatively early in the reservoir’s life. / text

Developing hybrid processes for heavy oil recovery is a major area of interest in recent years. The need for such processes originates from the challenges of heavy oil recovery relating to fluid injectivity, reservoir heating, and oil displacement and production. These challenges are particularly profound in shaley thin oil deposits where steam injection is not feasible and other recovery methods should be employed. In this work, we aim to develop and optimize a hybrid process that involves moderate reservoir heating and chemical enhanced oil recovery (EOR). This process, in its basic form, is a three-stage scheme. The first stage is a short electrical heating, in which the reservoir temperature is raised just enough to create fluid injectivity. After electrical heating has created sufficient fluid injectivity, high-rate high-pressure hot water injection accelerates the raise in temperature of the reservoir and assists oil production. At the end of hot waterflooding the oil viscosities are low enough for an Alkali-Co-solvent-Polymer (ACP) chemical flood to be performed where oil can efficiently be mobilized and displaced at low pressure gradients. A key aspect of ultra-low IFT chemical flood, such as ACP, is the rheology of the microemulsions that form in the reservoir. Undesirable rheology impedes the displacement of the chemical slug in the reservoir and results in poor process performance or even failure. The viscosity of microemulsions can be altered by the addition of co-solvents and branched or twin-tailed co-surfactants and by an increase in temperature. To reveal the underlying mechanisms, a consistent theoretical framework was developed. Employing the membrane theory and electrostatics, the significance of charge and/or composition heterogeneity in the interface membrane and the relevance of each to the above-mentioned alteration methods was demonstrated. It was observed that branched co-surfactants (in mixed surfactant formulations) and temperature only modify the saddle-splay modulus (k ̅) and bending modulus (k) respectively, whereas co-solvent changes both moduli. The observed rheological behavior agrees with our findings. To describe the behavior of microemulsions in flow simulations, a rheological model was developed. A key feature of this model is the treatment of the microemulsion as a bi-network. This provides accuracy and consistency in the calculation of the zero-shear viscosity of a microemulsion regardless of its type and microstructure. Once model parameters are set, the model can be used at any concentration and shear rate. A link between the microemulsion rheological behavior and its microstructure was demonstrated. The bending modulus determines the magnitude of the viscous dissipations and the steady-shear behavior. The new model, additionally, includes components describing the effects of rheology alteration methods. Experimental viscosity data were used to validate the new microemulsion viscosity model. Several ACP corefloods showing the large impact of microemulsion viscosity on process performance were matched using the UTCHEM simulator with the new microemulsion rheology model added to the code. Finally, numerical simulations based on Peace River field data were performed to investigate the performance of the proposed hybrid thermal-chemical process. Key design parameters were identified to be the method of heating, duration of the heating, ACP slug size and composition, polymer drive size, and polymer concentration in the polymer drive. An optimization study was done to demonstrate the economic feasibility of the process. The optimization revealed that short electrical heating and high-rate high-pressure waterflooding are necessary to minimize the energy use and operational expenses. The optimum slug and polymer drive sizes were found to be ~0.25 PV and ~1 PV, respectively. It was shown that the well costs dominate the expenditure and the overall cost of the optimized process is in the range of 20-30 $⁄bbl of incremental oil production. / text

Polymers used in enhanced oil recovery (EOR) help to control the mobility ratio between oil and aqueous phases and as a result, polymer flooding improves sweep efficiency in reservoirs. However, the conventional wisdom is that polymer flooding does not have considerable effect on pore-level displacement because pressure forces would not be enough to overcome trapping caused by capillary forces. Recently, both coreflood experiments and field data suggest that injecting viscoelastic polymers, such as hydrolyzed polyacrylamide (HPAM), can result in lower residual oil saturation. The hypothesis is that the polymer elasticity provides several pore-level mechanisms for oil mobilization that are generally not significant for purely-viscous fluids. Both experiments and modeling need to be performed to investigate the effect of polymer elasticity on residual oil saturation. Pore-scale modeling and micro-fluidic experiments can be used to investigate pore-level physics, and then used to upscale to the macro-scale. The objective of this work is to understand the effect of polymer elasticity on apparent viscosity and residual oil saturation in porous media. Single- and multi-phase pore-level computational fluid dynamics (CFD) modeling for viscoelastic polymer flow is performed to investigate the dominant mechanisms at the pore level to mobilize trapped oil. Several interesting results are found from the CFD results. First, the elasticity of the polymer results in an increase in normal stress at the pore-level; therefore, the normal stresses exerted on a static oil droplet are significant and not negligible as for a purely-viscous fluid. The CFD results show that viscoelastic fluid exerts additional forces on the oil-phase which may help mobilize trapped oil out of the porous medium. Second, due to the elasticity of polymer, the viscoelastic polymer has some level of pulling effect; while passing above a dead-end pore it can pull out the trapped oil phase and then mobilize it. However, both CFD modeling and micro-fluidic experiments show the pulling-effect is not likely the main mechanism to reduce oil saturation at pore-level. Third, dynamic CFD simulations show less deformation of the oil phase while viscoelastic polymer is displacing fluid compared to purely viscous fluid. It may justify the hypothesis that polymer elasticity resists against snap-off mechanism. As a result, when viscoelastic polymer displaces the oil ganglia, the oil phase does not snap off, and the oil phase remains connected, and therefore easier to move in porous media compared to disconnected oil. For single phase flow, a closed-form flow equation has been developed based on CFD modeling in converging/diverging ducts representative of pore throats. The pore-level equations were substituted into a pore-network model and validated against experimental data. Good agreement is observed. This study reveals important findings about the effect of polymer elasticity to reduce the residual oil saturation; however, more experiments and simulations are recommended to fully-understand the mobilization mechanisms and take advantage of them to optimize the polymer-flooding process in the field. / text

Polymer flooding is a commercial enhanced oil recovery (EOR) method used to increase the sweep efficiency of water floods. Hydrolyzed polyacrylamide (HPAM), a synthetic commercial polymer, is widely used in commercial polymer floods and it is also used for mobility control of chemical floods using surfactants such as surfactant-polymer flooding and alkaline-surfactant-polymer flooding. The increase in the degree of hydrolysis of HPAM at elevated temperature or pH with time affects the polymer solution viscosity and its adsorption on rock surfaces. A polymer hydrolysis model based on published laboratory data was implemented in UTCHEM, a chemical EOR simulator, in order to assess the effect of hydrolysis on reservoir performance. Both 1D and 3D simulations were performed to validate the implementation of the model. The simulation results are consistent with the laboratory observations that show an increase in polymer solution viscosity as hydrolysis progresses. The numerical results indicate that hydrolysis occurs very rapidly and impacts the near wellbore region polymer injectivity. / text

Studying chemical enhanced oil recovery (EOR) in a high-temperature/high-salinity (HT/HS) reservoir will help expand the application of chemical EOR to more challenging environments. Until recently, chemical EOR was not recommended at reservoirs that contain high concentrations of divalent cations without the need to recondition the reservoir by flooding it with less saline/ less hardness brines. This strategy was found ineffective in preparing the reservoir for chemical flooding. Surfactants used for chemical flooding operating in high temperatures tend to precipitate when exposed to high concentrations of divalent cations and will partition to the oil phase at high salinities. In this study amphoteric surfactant was used to replace the traditionally used anionic surfactants. Amphoteric surfactants show higher multivalent cations tolerance with better thermal stability. A modified amphoteric surfactant with lower adsorption properties was evaluated for oil recovery. Organic alkali was used to eliminate the water softening process when preparing the chemical solution and reduce potential scale problems caused by precipitation due to incompatibility between chemical slug containing alkali and formation brine. Using organic alkali helped in minimizing softening required when preparing an alkali-surfactant-polymer (ASP) solution using seawater. Solution prepared with organic alkali showed the least injectivity decline when compared to traditional alkalis (NaOH and Na2CO3) and sodium metaborate. Adding organic alkali helped further reduce IFT values when added to surfactant solution. Amphoteric surfactant was found to produce low IFT values at low concentrations and can operate at high salinity / high hardness conditions. When mixed with polymer it improved the viscosity of the surfactant-polymer (SP) solution when prepared in high salinity mixing water (6% NaCl). When prepared in seawater and tested in reservoir temperature (95°C) no reduction in viscosity was found. Unlike the anionic surfactant that causes reduction in viscosity of the SP solution at reservoir temperature. This will not require increasing the polymer concentration in the chemical slug. Unlike the case when anionic surfactant was used and more polymer need to be added to compensate the reduction in viscosity. Berea sandstone cores show lower recovery compared to dolomite cores. It was also found that Berea cores were more sensitive to polymer concentration and type and injectivity decline can be a serious issue during chemical and polymer injection. Dolomite did not show injectivity decline during chemical and polymer flooding and was not sensitive to the polymer concentration when a polymer with low molecular weight was used. CT scan was employed to study the displacement of oil during ASP, SP, polymer and surfactant flooding. The formation and propagation oil bank was observed during these core flood experiments. ASP and SP flooding showed the highest recovery, and formation and propagation of oil bank was clearer in these experiments compared to surfactant flooding. It was found that in Berea sandstone with a permeability range of 50 to 80 md that the recovery and fluid flow was through some dominating and some smaller channels. This explained the deviation from piston-like displacement, where a sharp change in saturation in part of the flood related to the dominated channels and tapered front with late arrival when oil is recovered from the smaller channels. It was concluded that the recovery in the case of sandstone was dominated by the fluid flow and chemical propagation in the porous media not by the effectiveness of the chemical slug to lower the IFT between the displacing fluid and oil.

A model for the evaluation of CO2-EOR projects has been developed. This model includes both reservoir simulation to handle reservoir properties, fluid flow and injection and production schedules, and a numerical economic model that generates a monthly cash flow stream from the outputs of the reservoir model. This model is general enough to be used with any project and provide a solid common basis to all of them. This model was used to evaluate CO2-EOR injection and production strategies and develop an optimization workflow. Producer constraints (maximum oil and gas production rates) should be optimized first to generate a reference case. Further improvements can then be obtained by optimizing the injection starting date and the injection plateau rate. Investigation of sensitivity of CO2-EOR to the presence of an aquifer showed that CO2 injection can limit water influx in the reservoir and is beneficial to recovery, even with a strong water drive. The influence of some key parameters was evaluated: the producer should be completed in the top part of the reservoir, while the injector should be completed over the entire thickness; it is recommended but not mandatory that the injection should start as early as possible to allow for lower water cut limit. Finally, the sensitivity of the economics of the projects to some key parameters was evaluated. The most influent parameter is by far the oil price, but other parameters such as the CO2 source to field distance, the pipeline cost scenario, the CO2 source type or the CO2 market price have roughly the same influence. It is therefore possible to offset an increase of one of them by reducing another.

The southern end of the Black Warrior Basin has been the site of limited drilling operations, but a critical need now exists to establish a greater understanding of the regional stratigraphy. The objectives of this study were to define a sequence stratigraphic framework for the southernmost Black Warrior Basin, to identify chronostratigraphic timelines within depositional environments, identify regional transgressive and high stand systems tracts. This information was used to identify three target reservoirs, characterize petrophysical properties, and confirm integrity of reservoir and seal formations for geologic storage. Methods include correlation of petrophysical well logs in the study area, well log analysis, as well as petrographic and core analyses. Five cycles were identified in well log cross sections. Sequence boundaries will be identified in both cross sections seismic data. Cretaceous sediments deposited above a regional sequence boundary above the Paleozoic that may represent as much as ~141 ma of erosion or non-deposition. The results of this study will contribute to development of a proposed geologic carbon sequestration facility in Kemper County, Mississippi.