Numerical reservoir simulation has been used to describe mechanism of methane gas desorption process, diffusion process, and fluid flow in a coalbed methane reservoir. The reservoir simulation model reflects the response of a reservoir system and the relationship among coalbed methane reservoir properties, operation procedures, and gas production. This work presents a procedure to select the optimum well spacing scenario by using a reservoir simulation. This work uses a two-phase compositional simulator with a dual porosity model to investigate well-spacing effects on coalbed methane production performance and methane recovery. Because of reservoir parameters uncertainty, a sensitivity and parametric study are required to investigate the effects of parameter variability on coalbed methane reservoir production performance and methane recovery. This thesis includes a reservoir parameter screening procedures based on a sensitivity and parametric study. Considering the tremendous amounts of simulation runs required, this work uses a regression analysis to replace the numerical simulation model for each wellspacing scenario. A Monte Carlo simulation has been applied to present the probability function. Incorporated with the Monte Carlo simulation approach, this thesis proposes a well-spacing study procedure to determine the optimum coalbed methane development scenario. The study workflow is applied in a North America basin resulting in distinct Net Present Value predictions between each well-spacing design and an optimum range of well-spacing for a particular basin area.

The main objective of this research is to devise a risk-based methodology for developing emergency operation schedules (EOS). EOS are decision tools that provide guidance to reservoir operators in charge of making real-time release decisions during major flood events. A computer program named REOS was created to perform the computations to develop risk-based EOS. The computational algorithm in REOS is divided in three major components: (1) synthetic streamflow generation, (2) mass balance computations, and (3) frequency analysis. The methodology computes the required releases to limit storage to the capacity available based on the probabilistic properties of future flows, conditional to current streamflow conditions. The final product is a series of alternative risk-based EOS in which releases, specified as a function of reservoir storage level, current and past inflows, and time of year, are associated with a certain risk of failing to attain the emergency operations objectives. The assumption is that once emergency operations are triggered by a flood event, the risk associated with a particular EOS reflects the probability of exceeding a pre-established critical storage level given that the same EOS is followed throughout the event. This provides reservoir operators with a mechanism for evaluating the tradeoffs and potential consequences of release decisions. The methodology was applied and tested using the Addicks and Barker Reservoir system in Houston, TX as a case study. Upstream flooding is also a major concern for these reservoirs. Modifications to the current emergency policies that would allow emergency releases based on the probability of upstream flooding are evaluated. Riskbased EOS were tested through a series of flood control simulations. The simulations were performed using the HEC-ResSim reservoir simulation model. Rainfall data recorded from Tropical Storm Allison was transposed over the Addicks and Barker watersheds to compute hypothetical hydrographs using HEC-HMS. Repeated runs of the HEC-ResSim model were made using different flooding and residual storage scenarios to compare regulation of the floods under alternative operating policies. An alternative application of the risk-based EOS in which their associated risk was used to help quantify the actual probability of upstream flooding in Addicks and Barker was also presented.

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

Slickwater with sand is the most commonly used hydraulic fracturing treatment for shale reservoirs. The slickwater treatment produces long skinny fractures, but only the near wellbore region is propped due to fast settling of sand. Adding gel into water can prevent the fast settling of sand, but gel may damage the fracture surface and proppant pack. Moreover, current water-based fracturing consumes a large amount of water, has high water leakage, and imposes high water disposal costs. The goal of this project is to develop non-damaging, less water-intensive fracturing treatments for shale gas reservoirs with improved proppant placement efficiency. Earlier studies have proposed to replace sand with ultra-light weight proppants (ULWP) to enhance proppant transport, but it is not used commonly in field. This study evaluates the performance of three kinds of ULWPs covering a wide range of specific gravity and representing the three typical manufacturing methods. In addition to replacing sand with ULWPs, replacing water with foams can be an alternative treatment that reduces water usage and decreases proppant settling. Polymer-added foams have been used in conventional reservoirs to improve proppant placement efficiency. However, polymers can damage shale permeability in unconventional reservoirs. This dissertation studies polymer-free foams (PFF) and evaluates their performance. This study uses both experiments and simulations to assess the productivity and profitability of the ULWP treatment and PFF treatment. First, a reservoir simulation model is built in CMG to study the impact of fracture conductivity and propped length on fracture productivity. This model assumes a single fracture intersecting a few reactivated natural fractures. Second, a 2D fracturing model is used to simulate the fracture propagation and proppant transport. Third, strength, API conductivity and gravity settling rates are measured for three ULWPs. Fourth, foam stability tests are conducted to screen the best PFF agents and the selected foams are put into a circulating loop to study their rheology. Finally, empirical correlations from the experiments are applied in the fracturing model and reservoir model to predict productivity by using the ULWPs with slickwater or using the PFFs with sand. Experimental results suggest that, at 4000 psi with concentrations varying from partial monolayer (0.05 lb/ft²) to multilayer (1 lb/ft²), ULW-1 (polymeric) is the most deformable with conductivity of 1-10 md-ft. ULW-2 (resin coated and impregnated ground walnut hull) is the second most deformable with similar conductivity. ULW-3 (resin coated porous ceramic) is the least deformable with conductivity of 20-1000 md-ft, which is comparable to sand. Three foam formulations (A, B: regular surfactant foam, C: viscoelastic surfactant foam) are selected based on the stability results of fourteen surfactants. All PFFs exhibit power-law rheological behavior in a laminar flow regime. The power law parameters of the regular surfactant PFF depend on both quality and pressure when quality is higher than 60% but depend on quality only when quality is lower than 60%. Simulation results suggest that under the optimal concentration of 0.04-0.06 v/v (0.37-0.55 lb/gal) for both ULW-1 and ULW-2, and 0.1 v/v (1.46 lb/gal) for ULW-3, 1-year cumulative production for 0.1 µD shale reservoir is higher than sand by 127% for ULW-1, 28% for ULW-2, and 38% for ULW-3. The productivity benefits decrease as shale permeability increases for all three ULWPs. ULW-1 and ULW-2 have higher productivity benefits for longer production time, while ULW-3 has relatively constant productivity benefits over time. The economic profit of ULW-1 when priced at $5/lb is 2.2 times larger than that of sand for 1-year production in 0.1 µD shale reservoirs; the acceptable maximum price is $10/lb for ULW-1, $6/lb for ULW-2, and $2.5/lb for ULW-3. The maximum price increases as production time increases. The PFFs with a quality of 60% carrying mesh 40 sand at a partial monolayer concentration of 0.04 v/v (0.88 lb/gal) can generate 50% higher productivity, 74% higher economic profit, and over 300% higher water efficiency than the best slickwater-sand case (mesh 40 sand at 0.1 v/v) for 1-year production in 0.1µD shale reservoirs. The benefits of using the PFFs decrease with increasing shale permeability, increasing production time, or decreasing pumping time. This dissertation gives a range of field conditions where the ULWP and PFF may be more effective than slickwater-sand fracturing. / text

The primary controls on subsurface reservoir heterogeneities and fluid flow characteristics are sedimentary facies architecture and petrophysical rock fabric distribution in clastic reservoirs and in carbonate reservoirs, respectively. Facies models are critical and fundamental for summarizing facies and facies architecture in data-rich areas. Facies models also assist in predicting the spatial architectural trend of sedimentary facies in other areas where subsurface information is lacking. The method for transferring geological information from different facies models into digital data and then generating associated numerical models is called facies modeling or geological modeling. Facies modeling is also vital to reservoir simulation and reservoir characterization analysis. By extensively studying and reviewing the relevant research in the published literature, this report identifies and analyzes the best and most detailed geologic data that can be used in facies modeling, and the most current geostatistical and stochastic methods applicable to facies modeling. Through intensive study of recent literature, the author (1) summarizes the basic concepts and their applications to facies and facies models, and discusses a variety of numerical modeling methods, including geostatistics and stochastic facies modeling, such as variogram-based geostatistics modeling, object-based stochastic modeling, and multiple-point geostatistics modeling; and (2) recognizes that the most effective way to characterize reservoir is to integrate data from multiple sources, such as well data, outcrop data, modern analogs, and seismic interpretation. Detailed and more accurate parameters using in facies modeling, including grain size, grain type, grain sorting, sedimentary structures, and diagenesis, are gained through this multidisciplinary analysis. The report concludes that facies and facies models are scale dependent, and that attention should be paid to scale-related issues in order to choose appropriate methods and parameters to meet facies modeling requirements. / text

No description available.

In Santa Barbara County, Santa Maria Pacific (an exploration and production company) is expanding their cyclic steam project in a diatomite reservoir. The hilly or mountainous topography and cut and fill restrictions have interfered with the company's ideal development plan. The steep hillsides prevent well pad development for about 22 vertical well locations in the 110 well expansion plan. Conventional production performs poorly in the area because the combination of relatively low permeability (1-10 md) and high viscosity (~220 cp) at the reservoir temperature. Cyclic steam injection has been widely used in diatomite reservoirs to take advantage of the diatomite rocks unique properties and lower the viscosity of the oil. Some companies used deviated wells for cyclic steam injection, but Santa Maria Pacific prefers the use only vertical wells for the expansion. Currently, the inability to create well pads above 22 vertical well target locations will result in an estimated $60,000,000 of lost revenue over a five year period. The target locations could be developed with unstimulated deviated or horizontal wells, but expected well rates and expenses have not been estimated. In this work, I use a thermal reservoir simulator to estimate production based on five potential development cases. The first case represents no development other than the cyclic wells. This case is used to calibrate the model based on the pilot program performance and serves as a reference point for the other cases. Two of the cases simulate a deviated well with and without artificial lift next to a cyclic well, and the final two cases simulate a horizontal well segment with and without artificial lift next to a cyclic well. The deviated well with artificial lift results in the highest NPV and profit after five years. The well experienced pressure support from the neighboring cyclic well and performed better with the cyclic well than without it. Adding 22 deviated wells with artificial lift will increase the project's net profit by an estimated $7,326,000 and NPV by $2,838,000 after five years.

No description available.

Various models featuring horizontal wells with multiple induced fractures have been proposed to characterize flow behavior over time in tight gas and shale gas systems. Currently, there is little consensus regarding the effects of non-ideal fracture geometries and coupled primary-secondary fracture interactions on reservoir performance in these unconventional gas reservoirs. This thesis provides a grid construction tool to generate high-resolution unstructured meshes using Voronoi grids, which provides the flexibility required to accurately represent complex geologic domains and fractures in three dimensions. Using these Voronoi grids, the interaction between propped hydraulic fractures and secondary "stress-release" fractures were evaluated. Additionally, various primary fracture configurations were examined, where the fractures may be non-planar or non-orthogonal. For this study, a numerical model was developed to assess the potential performance of tight gas and shale gas reservoirs. These simulations utilized up to a half-million grid-blocks and consider a period of up to 3,000 years in some cases. The aim is to provide very high-definition reference numerical solutions that will exhibit virtually all flow regimes we can expect in these unconventional gas reservoirs. The simulation results are analyzed to identify production signatures and flow regimes using diagnostic plots, and these interpretations are confirmed using pressure maps where useful. The coupled primary-secondary fracture systems with the largest fracture surface areas are shown to give the highest production in the traditional "linear flow" regime (which occurs for very high conductivity vertical fracture cases). The non-ideal hydraulic fracture geometries are shown to yield progressively lower production as the angularity of these fractures increases. Hence, to design optimum fracture completions, we should endeavor to keep the fractures as orthogonal to the horizontal well as possible. This work expands the current understanding of flow behavior in fractured tight-gas and shale-gas systems and may be used to optimize fracture and completion design, to validate analytical models and to facilitate more accurate reserves estimation.

For a stress-sensitive or stress-dependent reservoir, the interactions between its seepage field and in situ stress field are complex and affect hydrocarbon recovery. A coupled geomechanics and fluid-flow model can capture these relations between the fluid and solid, thereby presenting more precise history matchings and predictions for better well planning and reservoir management decisions. A traditional reservoir simulator cannot adequately or fully represent the ongoing coupled fluid-solid interactions during the production because of using the simplified update-formulation for porosity and the static absolute permeability during simulations. Many researchers have studied multiphase fluid-flow models coupled with geomechanics models during the past fifteen years. The purpose of this research is to develop a coupled geomechanics and compositional model and apply it to problems in the oil recovery processes. An equation of state compositional simulator called the General Purpose Adaptive Simulator (GPAS) is developed at The University of Texas at Austin and uses finite difference / finite control volume methods for the solution of its governing partial differential equations (PDEs). GPAS was coupled with a geomechanics model developed in this research, which uses a finite element method for discretization of the associated PDEs. Both the iteratively coupled solution procedure and the fully coupled solution procedure were implemented to couple the geomechanics and reservoir simulation modules in this work. Parallelization, testing, and verification for the coupled model were performed on parallel clusters of high-performance workstations. MPI was used for the data exchange in the iteratively coupled procedure. Different constitutive models were coded into GPAS to describe complicated behaviors of linear or nonlinear deformation in the geomechanics model. In addition, the geomechanics module was coupled with the dual porosity model in GPAS to simulate naturally fractured reservoirs. The developed coupled reservoir and geomechanics simulator was verified using analytical solutions. Various reservoir simulation case studies were carried out using the coupled geomechanics and GPAS modules. / text