Miscible gas flooding is a common method for enhanced oil recovery. Reliable design of miscible gas flooding requires compositional reservoir simulation that can accurately predict the fluid properties resulting from mass transfer between reservoir oil and injection gas. Drawbacks of compositional simulation are the efficiency and robustness of phase equilibrium calculations consisting of flash calculations and phase stability analysis. Simulation of multicontact miscible gas flooding involves a large number of phase equilibrium calculations in a near-critical region, where the calculations are time-consuming and difficult. Also, mixtures of reservoir oil and solvent such as CO₂ and rich gas can exhibit complex phase behavior at temperatures typically below 120°F, where three hydrocarbon-phases can coexist. However, most compositional simulators do not attempt to solve for three hydrocarbon-phases because three-phase equilibrium calculations are more complicated, difficult, and time-consuming than traditional two-phase equilibrium calculations. Due to the lack of robust algorithms for three-phase equilibrium calculations, the effect of a third hydrocarbon-phase on low-temperature oil displacement is little known. We develop robust and efficient algorithms for phase equilibrium calculations for two and three phases. The algorithms are implemented in a compositional reservoir simulator. Simulation case studies show that our algorithms can significantly decrease the computational time without loss of accuracy. Speed-up of 40% is achieved for a reservoir simulation using 20 components, compared to standard algorithms. Speed-up occurs not only because of improved computational efficiency but also because of increased robustness resulting in longer time-step sizes. We demonstrate the importance of three-phase equilibrium calculations, where simulations with two-phase equilibrium approximations proposed in the literature can result in complete failure or erroneous simulation results. Using the robust phase equilibrium algorithms developed, the mechanism is investigated for high efficiency of low-temperature oil displacements by CO₂ involving three hydrocarbon-phases. Results show that high displacement efficiency can be achieved when the composition path goes near the critical endpoint where the gaseous and CO₂-rich liquid phases merge in the presence of the oleic phase. Complete miscibility may not be developed for three-phase flow without considering the existence of a tricritical point. / text

Minimum miscibility pressure (MMP) is a key parameter in the design of gas flooding. There are experimental and computational methods to determine MMP. Computational methods are fast and convenient alternatives to otherwise slow and expensive experimental procedures. This research focuses on the computational aspects of MMP estimation. It investigates the shortcomings of the current computational models and offers ways to improve the robustness of MMP estimation. First, we develop a new mixing cell method of estimating MMP that, unlike previous "mixing cell" methods, uses a variable number of cells and is independent of gas-oil ratio, volume of the cells, excess oil volumes, and the amount of gas injected. The new method relies entirely on robust P-T flash calculations using any cubic equation-of-state (EOS). We show that mixing cell MMPs are comparable with those of other analytical and experimental methods, and that our mixing cell method finds all the key tie lines predicted by MOC; however, the method proved to be more robust and reliable than current analytical methods. Second, we identify a number of problems with analytical methods of MMP estimation, and demonstrate them using real oil characterization examples. We show that the current MOC results, which assume that shocks exist from one key tie line to the next may not be reliable and may lead to large errors in MMP estimation. In such cases, the key tie lines determined using the MOC method do not control miscibility, likely as a result of the onset of L₁-L₂-V behavior. We explain the problem with a simplified pseudo-ternary model and offer a procedure for determining when an error exists and for improving the results. Finally, we present a simple mathematical model for predicting the MMP of contaminated gas. Injection-gas compositions often vary during the life of a gasflood because of reinjection and mixing of fluids in situ. Determining the MMP by slim-tube or other methods for each possible variation in the gas-mixture composition is impractical. Our method gives an easy and accurate way to determine impure CO₂ MMPs for variable field solvent compositions on the basis of just a few MMPs. Alternatively, the approach could be used to estimate the enrichment level required to lower the MMP to a desired pressure. / text

Tight gas and shale oil play an important role in energy security and in meeting an increasing energy demand. Hydraulic fracturing is a widely used technology for recovering these resources. The design and evaluation of hydraulic fracture operation is critical for efficient production from tight gas and shale plays. The efficiency of fracturing jobs depends on the interaction between hydraulic (induced) and naturally occurring discrete fractures. In this work, a coupled reservoir-fracture flow model is described which accounts for varying reservoir geometries and complexities including non-planar fractures. Different flow models such as Darcy flow and Reynold's lubrication equation for fractures and reservoir, respectively are utilized to capture flow physics accurately. Furthermore, the geomechanics effects have been included by considering a multiphase Biot's model. An accurate modeling of solid deformations necessitates a better estimation of fluid pressure inside the fracture. The fractures and reservoir are modeled explicitly allowing accurate representation of contrasting physical descriptions associated with each of the two. The approach presented here is in contrast with existing averaging approaches such as dual and discrete-dual porosity models where the effects of fractures are averaged out. A fracture connected to an injection well shows significant width variations as compared to natural fractures where these changes are negligible. The capillary pressure contrast between the fracture and the reservoir is accounted for by utilizing different capillary pressure curves for the two features. Additionally, a quantitative assessment of hydraulic fracturing jobs relies upon accurate predictions of fracture growth during slick water injection for single and multistage fracturing scenarios. It is also important to consistently model the underlying physical processes from hydraulic fracturing to long-term production. A recently introduced thermodynamically consistent phase-field approach for pressurized fractures in porous medium is utilized which captures several characteristic features of crack propagation such as joining, branching and non-planar propagation in heterogeneous porous media. The phase-field approach captures both the fracture-width evolution and the fracture-length propagation. In this work, the phase-field fracture propagation model is briefly discussed followed by a technique for coupling this to a fractured poroelastic reservoir simulator. We also present a general compositional formulation using multipoint flux mixed finite element (MFMFE) method on general hexahedral grids with a future prospect of treating energized fractures. The mixed finite element framework allows for local mass conservation, accurate flux approximation and a more general treatment of boundary conditions. The multipoint flux inherent in MFMFE scheme allows the usage of a full permeability tensor. An accurate treatment of diffusive/dispersive fluxes owing to additional velocity degrees of freedom is also presented. The applications areas of interest include gas flooding, CO₂ sequestration, contaminant removal and groundwater remediation. / text

Hybrid gas/chemical Enhanced Oil Recovery (EOR) methods are such novel techniques to increase oil production and oil recovery efficiency. Gas flooding using carbon dioxide, nitrogen, flue gas, and enriched natural gas produce more oil from the reservoirs by channeling gas into previously by-passed areas. Surfactant flooding can recover trapped oil by reducing the interfacial tension between oil and water phases. Hybrid gas/chemical EOR methods benefit from using both chemical and gas flooding. In hybrid gas/chemical EOR processes, surfactant solution is injected with gas during low-tension-gas or foam flooding. Polymer solution can also be injected alternatively with gas to improve the gas volumetric sweep efficiency. Most fundamentally, wide applications of hybrid gas/chemical processes are limited due to uncertainties in reservoir characterization and heterogeneity, due to the lack of understanding of the process and consequently lack of a predictive reservoir simulator to mechanistically model the process. Without a reliable simulator, built on mechanisms determined in the laboratory, promising field candidates cannot be identified in advance nor can process performance be optimized. In this research, UTCHEM was modified to model four-phase water, oil, microemulsion, and gas phases to simulate and interpret chemical EOR processes including free and/or solution gas. We coupled the black-oil model for water/oil/gas equilibrium with microemulsion phase behavior model through a new approach. Four-phase fluid properties, relative permeability, and capillary pressure were developed and implemented. The mass conservation equation was solved for total volumetric concentration of each component at standard conditions and pressure equation was derived for both saturated and undersaturated PVT conditions. To model foam flow in porous media, comprehensive research was performed comparing capabilities and limitations of implicit texture (IT) and population-balance (PB) foam models. Dimensionless foam bubble density was defined in IT models to derive explicitly the foam-coalescence-rate function in these models. Results showed that each of the IT models examined was equivalent to the LE formulation of a population-balance model with a lamella-destruction function that increased abruptly in the vicinity of the limiting capillary pressure, as in current population-balance models. Foam models were incorporated in UTCHEM to model low-tension-gas and foam flow processes in laboratory and field scales. The modified UTCEM reservoir simulator was used to history match published low-tension-gas and foam coreflood experiments. The simulations were also extended to model and evaluate hybrid gas/chemical EOR methods in field scales. Simulation results indicated a well-designed low-tension-gas flooding has the potential to recover the trapped oil where foam provides mobility control during surfactant and surfactant-alkaline flooding in reservoirs with very low permeability. / text