A numerical reservoir model was created to optimize CO₂ Huff-and-Puff operations in the Bakken Shale. Huff-and-Puff is an enhanced oil recovery treatment in which a well alternates between injection, soaking, and production. Injecting CO₂ into the formation and allowing it to “soak” re-pressurizes the reservoir and improves oil mobility, boosting production from the well. A compositional reservoir simulator was used to study the various design components of the Huff-and-Puff process in order to identify the parameters with the largest impact on recovery and understand the reservoir’s response to cyclical CO₂ injection. It was found that starting Huff-and-Puff too early in the life of the well diminishes its effectiveness, and that shorter soaking periods are preferable over longer waiting times. Huff-and-Puff works best in reservoirs with highly-conductive natural fracture networks, which allow CO₂ to migrate deep into the formation and mix with the reservoir fluids. The discretization of the computational domain has a large impact on the simulation results, with coarser gridding corresponding to larger projected recoveries. Doubling the number of hydraulic fractures per stage results in considerably greater CO₂ injection requirements without proportionally larger incremental recovery factors. Incremental recovery from CO₂ Huff-and-Puff appears to be insufficient to make the process commercially feasible under current economic conditions. However, re-injecting mixtures of CO₂ and produced hydrocarbon gases was proven to be technically and economically viable, which could significantly improve profit margins of Huff-and-Puff operations. A substantial portion of this project involved studying alternative numerical methods for modeling hydraulically-fractured reservoir models. A domain decomposition technique known as mortar coupling was used to model the reservoir system as two individually-solved subdomains: fracture and matrix. A mortar-based numerical reservoir simulator was developed and its results compared to a tradition full-domain finite difference model for the Cinco-Ley et al. (1978) finite-conductivity vertical fracture problem. Despite some numerical issues, mortar coupling closely matched Cinco-Ley et al.'s (1978) solution and has potential applications in complex problems where decoupling the fracture-matrix system might be advantageous. / text

[pt] A injeção de emulsão óleo-água vem sendo estudada como método de recuperação avançada, levando a um incremento no volume de óleo recuperando através de um controle de mobilidade do fluido injetado e aumento da eficiência de deslocamento de deslocamento do óleo. A aplicação deste ,étodo requer o entendimento do escoamento tanto na escala de poros como na escala de reservatórios. Neste trabalho, o efeito de emulsão no escoamento bifásico é modelado através de curvas de permeabilidades relativas que variam com a concentração de gotas da fase dispersa da emulsão. Estas curvas descrevem dois fenômenos físicos observados nestes escoamentos, o controle de mobilidade, representado pela diminuição da permeabilidade relativa da fase aquosa devido à presença de gostas, e a melhora na eficiência de deslocamento, representada pela diminuição da saturação residual da fase oleosa. O modelo proposto foi implementado em um simulador comercial (STARS-CMG) em uma geometria bi-dimensional com dois poços verticais ( injetor e produtor). A produção de óleo durante o processo de injeção alterada de água/ emulsão/ água é estudado em diversas condições do escoamento com o objetivo de analisar separadamente o mecanismo de controle de mobilidade e da melhoria do deslocamento no processo, determinando as condições ótimas de operações. Os resultados obtidos mostram que a injeção de emulsão óleo-água pode representar uma importante alternativa na produção de óleo como processo de recuperação avançada. / [en] Emulsion injection has been studied as an enhanced oil recovery method, leading to an increase on the recovered oil volume by mobility control of injected fluid and increase of oil displacement efficiency. The use of this method requires understanding the flow both in pore and reservoir scales. The effect of emulsion on the two phase flow is modeled by using relative permeability curves that vary with the concentration of the dispersed phase of the emulsion. Those curves describe two physical phenomena observed on these flows: control of mobility, represented by a lower relative permeability of the aqueous phase due to the presence of droplets and improved displacement efficiency, represented by the decrease in the residual oil saturation. The proposed model was implemented in a commercial simulator (STARS-CMG) in a two-dimensional geometry with two vertical wells (injector and producer). The oil production during alternate injection of water/emulsion/water process is studied with several flows condition with the purpose to analyze separately the effect of mobility control and improvement of the oil displacement in the process and finding the optimum operating conditions. The results obtained show that the oil-water emulsion injection can represent an important alternative in oil production as an enhanced oil recovery method.

Various analytical, semi-analytical, and empirical models have been proposed to characterize rate and pressure behavior as a function of time in tight/shale gas systems featuring a horizontal well with multiple hydraulic fractures. Despite a small number of analytical models and published numerical studies there is currently little consensus regarding the large-scale flow behavior over time in such systems. The purpose of this work is to construct a fit-for-purpose numerical simulator which will account for a variety of production features pertinent to these systems, and to use this model to study the effects of various parameters on flow behavior. Specific features examined in this work include hydraulically fractured horizontal wells, multiple porosity and permeability fields, desorption, and micro-scale flow effects. The theoretical basis of the model is described in Chapter I, along with a validation of the model. We employ the numerical simulator to examine various tight gas and shale gas systems and to illustrate and define the various flow regimes which progressively occur over time. We visualize the flow regimes using both specialized plots of rate and pressure functions, as well as high-resolution maps of pressure distributions. The results of this study are described in Chapter II. We use pressure maps to illustrate the initial linear flow into the hydraulic fractures in a tight gas system, transitioning to compound formation linear flow, and then into elliptical flow. We show that flow behavior is dominated by the fracture configuration due to the extremely low permeability of shale. We also explore the possible effect of microscale flow effects on gas effective permeability and subsequent gas species fractionation. We examine the interaction of sorptive diffusion and Knudsen diffusion. We show that microscale porous media can result in a compositional shift in produced gas concentration without the presence of adsorbed gas. The development and implementation of the micro-flow model is documented in Chapter III. This work expands our understanding of flow behavior in tight gas and shale gas systems, where such an understanding may ultimately be used to estimate reservoir properties and reserves in these types of reservoirs.

Various analytical, semi-analytical, and empirical models have been proposed to characterize rate and pressure behavior as a function of time in tight/shale gas systems featuring a horizontal well with multiple hydraulic fractures. Despite a small number of analytical models and published numerical studies there is currently little consensus regarding the large-scale flow behavior over time in such systems. The purpose of this work is to construct a fit-for-purpose numerical simulator which will account for a variety of production features pertinent to these systems, and to use this model to study the effects of various parameters on flow behavior. Specific features examined in this work include hydraulically fractured horizontal wells, multiple porosity and permeability fields, desorption, and micro-scale flow effects. The theoretical basis of the model is described in Chapter I, along with a validation of the model. We employ the numerical simulator to examine various tight gas and shale gas systems and to illustrate and define the various flow regimes which progressively occur over time. We visualize the flow regimes using both specialized plots of rate and pressure functions, as well as high-resolution maps of pressure distributions. The results of this study are described in Chapter II. We use pressure maps to illustrate the initial linear flow into the hydraulic fractures in a tight gas system, transitioning to compound formation linear flow, and then into elliptical flow. We show that flow behavior is dominated by the fracture configuration due to the extremely low permeability of shale. We also explore the possible effect of microscale flow effects on gas effective permeability and subsequent gas species fractionation. We examine the interaction of sorptive diffusion and Knudsen diffusion. We show that microscale porous media can result in a compositional shift in produced gas concentration without the presence of adsorbed gas. The development and implementation of the micro-flow model is documented in Chapter III. This work expands our understanding of flow behavior in tight gas and shale gas systems, where such an understanding may ultimately be used to estimate reservoir properties and reserves in these types of reservoirs.

The ever increasing energy demand brings about widespread interest to rapidly, profitably and efficiently develop unconventional resources, among which tight gas sands hold a significant portion. However, optimization of development strategies in tight gas fields is challenging, not only because of the wide range of depositional environments and large variability in reservoir properties, but also because the evaluation often has to deal with a multitude of wells, limited reservoir information, and time and budget constraints. Unfortunately, classical full-scale reservoir evaluation cannot be routinely employed by small- to medium-sized operators, given its timeconsuming and expensive nature. In addition, the full-scale evaluation is generally built on deterministic principles and produces a single realization of the reservoir, despite the significant uncertainty faced by operators. This work addresses the need for rapid and cost-efficient technologies to help operators determine optimal well spacing in highly uncertain and risky unconventional gas reservoirs. To achieve the research objectives, an integrated reservoir and decision modeling tool that fully incorporates uncertainty was developed. Monte Carlo simulation was used with a fast, approximate reservoir simulation model to match and predict production performance in unconventional gas reservoirs. Simulation results were then fit with decline curves to enable direct integration of the reservoir model into a Bayesian decision model. These integrated tools were applied to the tight gas assets of Unconventional Gas Resources Inc. in the Berland River area, Alberta, Canada.

In this study, a three-dimensional, three-phase dynamic simulation model based on geological investigations of Korolev oilfield in Kazakhstan was used as a development planning tool in order to improve performance of three new wells. The model, developed previously by means of a seismic study, well log and core data, incorporating with characteristics of oilfield productivity, properties of reservoir, liquids and gases that are saturating the hydrocarbon-bearing horizon can be used to calculate development parameters for Korolev field, including production well locations, drilling schedules, and to facilitate both long-term and short-term forecasting for the purposes of optimizing the hydrocarbon recovery from the field. The objective of this work is to assess the impact of adding 3 extra producing wells and find ways to optimize cumulative production with the least impact on the existing development plan by means of deeper understanding subsurface dynamic processes of the Korolev field. The challenge is a high degree of connectivity between wells in the productive formation throughout the field so that any change of production parameters affects the whole field&rsquo / s cumulative production. Trying to find a solution to optimum production of the reservoir forecast studies were carried out, the impact of each new well on development parameters was defined, sub-surface processes changes due to extra producers lead-in were explained and as a result of this thesis two optimization models were proposed, one of which will bring nearly 9.7 million barrels more oil.

Nowadays, field development strategy has become increasingly dependent on the results of reservoir simulation models. Reservoir studies demand fast and efficient results to make investment decisions that require a reasonable trade off between accuracy and simulation time. One of the suitable options to fulfill this requirement is streamline reservoir simulation technology, which has become very popular in the last few years. Streamline (SL) simulation provides an attractive alternative to conventional reservoir simulation because SL offers high computational efficiency and minimizes numerical diffusion and grid orientation effects. However, streamline methods have weaknesses incorporating complex physical processes and can also suffer numerical accuracy problems. The main objective of this research is to evaluate the numerical accuracy of the latest SL technology, and examine the influence of different factors that may impact the solution of SL simulation models. An extensive number of numerical experiments based on sensitivity analysis were performed to determine the effects of various influential elements on the stability and results of the solution. Those experiments were applied to various models to identify the impact of factors such as mobility ratios, mapping of saturation methods, number of streamlines, time step sizes, and gravity effects. This study provides a detailed investigation of some fundamental issues that are currently unresolved in streamline simulation.

Reconciling high-resolution geologic models to field production history is still a very time-consuming procedure. Recently streamline-based assisted and automatic history matching techniques, especially production data integration by “travel-time matching,” have shown great potential in this regard. But no systematic study was done to examine the merits of travel-time matching compared to more traditional amplitude matching for field-scale application. Besides, most applications were limited to two-phase water-oil flow because current streamline models are limited in their ability to incorporate highly compressible flow in a rigorous and computationally efficient manner. The purpose of this work is fourfold. First, we quantitatively investigated the nonlinearities in the inverse problems related to travel time, generalized travel time, and amplitude matching during production data integration and their impact on the solution and its convergence. Results show that the commonly used amplitude inversion can be orders of magnitude more nonlinear compared to the travel-time inversion. Both the travel-time and generalized travel time inversion (GTTI) are shown to be more robust and exhibit superior convergence characteristics. Second, the streamline-based assisted history matching was enhanced in two important aspects that significantly improve its efficiency and effectiveness. We utilize streamline-derived analytic sensitivities to determine the location and magnitude of the changes to improve the history match, and we use the iterative GTTI for model updating. Our approach leads to significant savings in time and manpower. Third, a novel approach to history matching finite-difference models that combines the efficiency of analytical sensitivity computation of the streamline models with the versatility of finite-difference simulation was developed. Use of finite-difference simulation can account for complex physics. Finally, we developed an approach to history matching three-phase flow using a novel compressible streamline formulation and streamline-derived analytic sensitivities. Streamline models were generalized to account for compressible flow by introducing a relative density of total fluids along streamlines and a density-dependent source term in the saturation equation. The analytical sensitivities are calculated based on the rigorous streamline formulation. The power and utility of our approaches have been demonstrated using both synthetic and field examples.

The demand for oil and gas is increasing yearly, whereas proven oil and gas reserves are being depleted. The potential of stripper oil and gas fields to supplement the national energy supply is large. In 2006, stripper wells accounted for 15% and 8% of US oil and gas production, respectively. With increasing energy demand and current high oil and gas prices, integrated reservoir studies, secondary and tertiary recovery methods, and infill drilling are becoming more common as operators strive to increase recovery from stripper oil and gas fields. The primary objective of this research was to support optimized production of oil and gas from stripper well fields by evaluating one stripper gas field and one stripper oil field. For the stripper gas field, I integrated geologic and engineering data to build a detailed reservoir characterization model of the Second White Specks (SSPK) reservoir in Garden Plains field, Alberta, Canada. The objectives of this model were to provide insights to controls on gas production and to validate a simulation-based method of infill drilling assessment. SSPK was subdivided into Units A ? D using well-log facies. Units A and B are the main producing units. Unit A has better reservoir quality and lateral continuity than Unit B. Gas production is related primarily to porosity-netthickness product and permeability and secondarily to structural position, minor structural features, and initial reservoir pressure. For the stripper oil field, I evaluated the Green River formation in the Wells Draw area of Monument Butte field, Utah, to determine interwell connectivity and to assess optimal recovery strategies. A 3D geostatistical model was built, and geological realizations were ranked using production history matching with streamline simulation. Interwell connectivity was demonstrated for only major sands and it increases as well spacing decreases. Overall connectivity is low for the 22 reservoir zones in the study area. A water-flood-only strategy provides more oil recovery than a primary-then-waterflood strategy over the life of the field. For new development areas, water flooding or converting producers to injectors should start within 6 months of initial production. Infill drilling may effectively produce unswept oil and double oil recovery. CO2 injection is much more efficient than N2 and CH4 injection. Water-alternating-CO2 injection is superior to continuous CO2 injection in oil recovery. The results of this study can be used to optimize production from Garden Plains and Monument Butte fields. Moreover, these results should be applicable to similar stripper gas and oil field fields. Together, the two studies demonstrate the utility of integrated reservoir studies (from geology to engineering) for improving oil and gas recovery.

Hydraulic fracturing is an important stimulation method for producing unconventional gas reserves. Natural fractures are present in many low-permeability gas environments and often provide important production pathways for natural gas. The production benefit from natural fractures can be immense, but it is difficult to quantify. The Mesaverde Group in the Piceance Basin in Colorado is a gas producing reservoir that has low matrix permeability but is also highly naturally fractured. Wells in the Piceance Basin are hydraulically fractured, so the production enhancements due to natural fracturing and hydraulic fracturing are difficult to decouple. In this thesis, dipole sonic logs were used to quantify geomechanical properties by combining stress equations with critically-stressed faulting theory. The properties derived from this log-based evaluation were used to numerically model hydraulic fracture treatments that had previously been pumped in the basin. The results from these hydraulic fracture models, in addition to the log-derived reservoir properties were used to develop reservoir models. Several methods for simulating the reservoir were compared and evaluated, including layer cake models, geostatistical models, and models simulating the fracture treatment using water injection. The results from the reservoir models were compared to actual production data to quantify the effect of both hydraulic fractures and natural fractures on production. This modeling also provided a framework upon which completion techniques were economically evaluated. / text