Carbonate reservoirs may have a variety of porosity types created by depositional, diagenetic, and fracture processes. This leads to the formation of complex pore systems, and in turn creates heterogeneities in reservoir performance and quality. In carbonate reservoirs affected by diagenesis and fracturing, porosity and peremeability can be independent of depositional facies or formation boundaries; consequently, conventional reservoir characterization methods are unreliable for predicting reservoir flow characteristics. This thesis provides an integrated petrographic, stratigraphic, and petrophysical study of the 'Canyon Reef' reservoir, a Pennsylvanian phylloid algal mound complex in the Horseshoe atoll. Core descriptions on three full-diameter cores led to the identification of 5 distinct depositional facies based on fundamental rock properties and biota. Fifty-four thin sections taken from the core were described are pores were classified using the Humbolt modification of the Ahr porosity classification. In order to rank reservoir quality, flow units were established on the basis of combined porosity and permeability values from core analysis. A cut off criterion for porosity and permeability was established to separate good and poor flow units. Ultimately cross sections were created to show the spatial distribution of flow units in the field.

Application of optimization techniques for determining the optimal operating policy for reservoirs is a major title in water resources planning and management. Genetic algorithms, ruled by evolution techniques, have become popular for solving optimization problems in diversified fields of science. The main aim of this research was to explore the efficiency and effectiveness of the applicability of genetic algorithm in optimization of multi-reservoirs. A computer code has been constructed for this purpose and verified by means of a reference problem with a known global optimum. Three reservoirs in the Colorado River Storage Project were optimized for maximization of energy production. Besides, a real-time approach utilizing a blend of online and a posteriori data was proposed. The results achieved were compared to the real operational data and genetic algorithms were found to be effective, competitive and can be utilized as an alternative technique to other traditional optimization techniques.

The purpose of this study was first to develop a method by which a detailed porosity classification system could be utilized to understand the relationship between pore/pore-throat geometry, genetic porosity type, and facies. Additionally, this study investigated the relationships between pore/pore-throat geometry, petrophysical parameters, and reservoir performance characteristics. This study focused on the Jurassic Smackover reservoir rocks of Grayson field, Columbia County, Arkansas. This three part study developed an adapted genetic carbonate pore type classification system, through which the Grayson reservoir rocks were uniquely categorized by a percent-factor, describing the effect of diagenetic events on the preservation of original depositional texture, and a second factor describing if the most significant diagenetic event resulted in porosity enhancement or reduction. The second part used petrographic image analysis and mercury-injection capillary pressure tests to calculate pore/pore-throat sizes. From these data sets pore/pore-throat sizes were compared to facies, pore type, and each other showing that pore-throat size is controlled by pore type and that pore size is controlled primarily by facies. When compared with each other, a pore size range can be estimated if the pore type and the median pore-throat aperture are known. Capillary pressure data was also used to understand the behavior of the dependent rock properties (porosity, permeability, and wettability), and it was determined that size-reduced samples, regardless of facies, tend to show similar dependent rock property behavior, but size-enhanced samples show dispersion. Finally, capillary pressure data was used to understand fluid flow behavior of pore types and facies. Oncolitic grainstone samples show unpredictable fluid flow behavior compared to oolitic grainstone samples, yet oncolitic grainstone samples will move a higher percentage of fluid. Size-enhanced samples showed heterogeneous fluid flow behavior while the size-reduced samples could be grouped by the number of modes of pore-throat sizes. Finally, this study utilized petrographic image analysis to determine if 2- dimensional porosity values could be calculated and compared to porosity values from 3-dimensional porosity techniques. The complex, heterogeneous pore network found in the Grayson reservoir rocks prevents the use of petrographic image analysis as a porosity calculation technique.

With the rapidly increasing demand for energy and the increasing prices for oil and gas, the role of unconventional gas reservoirs (UGRs) as energy sources is becoming more important throughout the world. Because of high risks and uncertainties associated with UGRs, their profitable development requires experts to be involved in the most critical development stages, such as drilling, completion, stimulation, and production. However, many companies operating UGRs lack this expertise. The advisory system we developed will help them make efficient decisions by providing insight from analogous basins that can be applied to the wells drilled in target basins. In North America, UGRs have been in development for more than 50 years. The petroleum literature has thousands of papers describing best practices in management of these resources. If we can define the characteristics of the target basin anywhere in the world and find an analogous basin in North America, we should be able to study the best practices in the analogous basin or formation and provide the best practices to the operators. In this research, we have built an advisory system that we call the Unconventional Gas Reservoir (UGR) Advisor. UGR Advisor incorporates three major modules: BASIN, PRISE and Drilling & Completion (D&C) Advisor. BASIN is used to identify the reference basin and formations in North America that are the best analogs to the target basin or formation. With these data, PRISE is used to estimate the technically recoverable gas volume in the target basin. Finally, by analogy with data from the reference formation, we use D&C Advisor to find the best practice for drilling and producing the target reservoir. To create this module, we reviewed the literature and interviewed experts to gather the information required to determine best completion and stimulation practices as a function of reservoir properties. We used these best practices to build decision trees that allow the user to take an elementary data set and end up with a decision that honors the best practices. From the decision trees, we developed simple computer algorithms that streamline the process.

Thesis (M.S.)--Eastern Illinois University, 2009. / Includes bibliographical references (leaves [117]-125).

Thesis (M. Phil.)--Hong Kong University of Science and Technology, 2003. / Includes bibliographical references (leaves 208-219). Also available in electronic version. Access restricted to campus users.

Concerns about elevated atmospheric carbon dioxide (CO 2 ) and the effect on global climate have created proposals for the reduction of carbon emissions from large stationary sources, such as power plants. Carbon dioxide capture and sequestration (CCS) in deep geologic units is being considered by Florida electric-utilities. Carbon dioxide-enhanced oil recovery (CO 2 -EOR) is a form of CCS that could offset some of the costs associated with geologic sequestration. Two potential reservoirs for geologic sequestration were evaluated in south-central and southern Florida: the Paleocene Cedar Keys Formation/Upper Cretaceous Lawson Formation (CKLIZ) and the Lower Cretaceous Sunniland Formation along the Sunniland Trend (Trend). The Trend is a slightly arcuate band in southwest Florida that is about 233 kilometers long and 32 kilometers wide, and contains oil plays within the Sunniland Formation at depths starting around 3,414 meters below land surface, which are confined to mound-like structures made of coarse fossil fragments, mostly rudistids. The Trend commercial oil fields of the South Florida Basin have an average porosity of 16% within the oil-producing Sunniland Formation, and collectively have an estimated storage capacity of around 26 million tons of CO 2 . The Sunniland Formation throughout the entire Trend has an average porosity of 14% and an estimated storage capacity of about 1.2 billion tons of CO 2 (BtCO2 ). The CKLIZ has an average porosity of 23% and an estimated storage capacity of approximately 79 BtCO 2 . Porous intervals within the CKLIZ and Sunniland Formation are laterally homogeneous, and low-permeability layers throughout the units provide significant vertical heterogeneity. The CKLIZ and Sunniland Formation are considered potentially suitable for CCS operations because of their geographic locations, appropriate depths, high porosities, estimated storage capacities, and potentiallyeffective seals. The Trend oil fields are suitable for CO 2 -EOR in the Sunniland Formation due to appropriate injected-CO 2 density, uniform intergranular porosity, suitable API density of formation-oil, sufficient production zones, and adequate remaining oil-in-place following secondary recovery. In addition to these in-depth investigations of the CKLIZ and Sunniland Formation, a more-cursory assessment of deep geologic units throughout the state of Florida, which includes rocks of Paleocene and Upper Cretaceous age through to rocks of Ordovician age, shows additional units in Florida that may be suitable for CO 2 -EOR and CCS operations. Furthermore, this study shows that deep geologic units throughout Florida potentially have the capacity to sequester billions of tons of CO 2 for hundreds of fossil-fuel-fired power plants. Geologic sequestration has not yet been conducted in Florida, and its implementation could prove useful to Florida utility companies, as well as to other energy-utilities in the southeastern United States.

This dissertation advances our understanding of how turbidity currents interact with three-dimensional (3-D) minibasin topography and the resulting deposits that form. Conceptual Gulf of Mexico-centric models of minibasin fill development have become the foundation for exploring and identifying strategic deep-water hydrocarbon reserves on continental slopes around the world. Despite the abundance of subsurface data, significant questions remain about the 3-D physical processes through which minibasins fill and the relationship between these processes and the topography of the basin. To overcome this problem, I utilize techniques in physical laboratory modeling to query established models of the role that turbidity currents play in minibasin fill development, and observe the relationships between fill from the Lobster minibasin located in a proximal continental slope position in the Gulf of Mexico and from the Safi Haute Mer (SHM) minibasin located in the distal continental slope of offshore western Morocco. First, existing published literature are reviewed and assessed for the known state of minibasin development and fill processes, and the strengths and weaknesses of our current knowledge base. Second, results are presented from two series of experiments that document the interaction between steady, depletive turbidity currents and 3-D minibasin topography. Experimental results suggest that turbidity currents produce deposits that are more likely to drape pre-flow topography than pond within it. Turbidity current velocity data show a strong 3-D physical component in minibasin fill sedimentation that also influences extra-basinal sedimentation patterns. Details of these results provide insight into processes that have not been previously considered in published conceptual models of minibasin fill. Third, a comparison of the two subsurface datasets show that the types and abundance of architectural elements vary depending on the location of the minibasin on the continental slope (i.e. proximal vs. distal), and suggests key differences in the processes responsible for their infilling. Finally, a comparison of experimental results to preserved deposit architectures in the Lobster and SHM datasets suggest a more complex relationship of process-driven sedimentation than that derived primarily from suspension fallout. This improved understanding of minibasin fill is applicable to industry for increasing confidence in subsurface interpretations and reducing risk while exploring for quality reservoirs in deepwater regions. / text

The growing demand for energy provides an incentive to pursue unconventional resources. Among these resources, tight gas and shale gas reservoirs have gained significant momentum because recent advances in technology allowed us to produce them at an economical rate. More importantly, they seem likely to contain a significant volume of hydrocarbon. There are, however, many questions concerning hydrocarbon production from these unconventional resources. For instance, in tight gas sandstone, we observe a significant variability in the producibilities of wells in the same field. The heterogeneity is even present in a single well with changes in depth. It is not clear what controls this heterogeneity. In shale gas, the pore connectivity inside the void space is not well explored and hence, a representative pore model is not available. Further, the effects of an adsorbed layer of gas and gas slippage on shale permeability are poorly understood. These effects play a crucial role in assigning a realistic permeability for shale in-situ from a laboratory measurement. In the laboratory, in contrast to in-situ, the core sample lacks the adsorbed layer because the permeability measurements are typically conducted at small pore pressures. Moreover, the gas slippages in laboratory and in-situ conditions are not identical. The present study seeks to investigate these discrepancies. Drainage and imbibition are sensitive to pore connectivity and unconventional gas transport is strongly affected by the connectivity. Hence, there is a strong interest in modeling mercury intrusion capillary pressure (MICP) test because it provides valuable information regarding the pore connectivity. In tight gas sandstone, the main objective of this research is to find a relationship between the estimated ultimate recovery (EUR) and the petrophysical properties measured by drainage/imbibition tests (mercury intrusion, withdrawal, and porous plate) and by resistivity analyses. As a measure of gas likely to be trapped in the matrix during production---and hence a proxy for EUR---we use the ratio of residual mercury saturation after mercury withdrawal (S[subscript gr]) to initial mercury saturation (S[subscript gi]), which is the saturation at the start of withdrawal. Crucially, a multiscale pore-level model is required to explain mercury intrusion capillary pressure measurements in these rocks. The multiscale model comprises a conventional network model and a tree-like pore structure (an acyclic network) that mimic the intergranular (macroporosity) and intragranular (microporosity) void spaces, respectively. Applying the multiscale model to porous plate data, we classify the pore spaces of rocks into macro-dominant, intermediate, and micro-dominant. These classes have progressively less drainage/imbibition hysteresis, which leads to the prediction that significantly more hydrocarbon is recoverable from microporosity than macroporosity. Available field data (production logs) corroborate the higher producibility of the microporosity. The recovery of hydrocarbon from micro-dominant pore structure is superior despite its inferior initial production (IP). Thus, a reservoir or a region in which the fraction of microporosity varies spatially may show only a weak correlation between IP and EUR. In shale gas, we analyze the pore structure of the matrix using mercury intrusion data to provide a more realistic model of pore connectivity. In the present study, we propose two pore models: dead-end pores and Nooks and Crannies. In the first model, the void space consists of many dead-end pores with circular pore throats. The second model supposes that the void space contains pore throats with large aspect ratios that are connected through the rock. We analyze both the scanning electron microscope (SEM) images of the shale and the effect of confining stress on the pore size distribution obtained from the mercury intrusion test to decide which pore model is representative of the in-situ condition. We conclude that the dead-end pores model is more representative. In addition, we study the effects of adsorbed layers of CH₄ and of gas slippage in pore walls on the flow behavior in individual conduits of simple geometry and in networks of such conduits. The network is based on the SEM image and drainage experiment in shale. To represent the effect of adsorbed gas, the effective size of each throat in the network depends on the pressure. The hydraulic conductance of each throat is determined based on the Knudsen number (Kn) criterion. The results indicate that laboratory measurements made with N₂ at ambient temperature and 5-MPa pressure, which is typical for the transient pulse decay method, overestimate the gas permeability in the early life of production by a factor of 4. This ratio increases if the measurement is run at ambient conditions because the low pressure enhances the slippage and reduces the thickness of the adsorbed layer. Moreover, the permeability increases nonlinearly as the in-situ pressure decreases during production. This effect contributes to mitigating the decline in production rates of shale gas wells. Laboratory data available in the literature for methane permeability at pressures below 7 MPa agree with model predictions of the effect of pressure. / text

Successful creations of multiple hydraulic fractures in horizontal wells are critical for economic development of unconventional reservoirs. The recent advances in diagnostic techniques suggest that multi-fracturing stimulation in unconventional reservoirs has often caused complex fracture geometry. The most important factors that might be responsible for the fracture complexity are fracture interaction and the intersection of the hydraulic and natural fracture. The complexity of fracture geometry results in significant uncertainty in fracturing treatment designs and production optimization. Modeling complex fracture propagation can provide a vital link between fracture geometry and stimulation treatments and play a significant role in economically developing unconventional reservoirs. In this research, a novel fracture propagation model was developed to simulate complex hydraulic fracture propagation in unconventional reservoirs. The model coupled rock deformation with fluid flow in the fractures and the horizontal wellbore. A Simplified Three Dimensional Displacement Discontinuity Method (S3D DDM) was proposed to describe rock deformation, calculating fracture opening and shearing as well as fracture interaction. This simplified 3D method is much more accurate than faster pseudo-3D methods for describing multiple fracture propagation but requires significantly less computational effort than fully three-dimensional methods. The mechanical interaction can enhance opening or induce closing of certain crack elements or non-planar propagation. Fluid flow in the fracture and the associated pressure drop were based on the lubrication theory. Fluid flow in the horizontal wellbore was treated as an electrical circuit network to compute the partition of flow rate between multiple fractures and maintain pressure compatibility between the horizontal wellbore and multiple fractures. Iteratively and fully coupled procedures were employed to couple rock deformation and fluid flow by the Newton-Raphson method and the Picard iteration method. The numerical model was applied to understand physical mechanisms of complex fracture geometry and offer insights for operators to design fracturing treatments and optimize the production. Modeling results suggested that non-planar fracture geometry could be generated by an initial fracture with an angle deviating from the direction of the maximum horizontal stress, or by multiple fracture propagation in closed spacing. Stress shadow effects are induced by opening fractures and affect multiple fracture propagation. For closely spaced multiple fractures growing simultaneously, width of the interior fractures are usually significantly restricted, and length of the exterior fractures are much longer than that of the interior fractures. The exterior fractures receive most of fluid and dominate propagation, resulting in immature development of the interior fractures. Natural fractures could further complicate fracture geometry. When a hydraulic fracture encounters a natural fracture and propagates along the pre-existing path of the natural fracture, fracture width on the natural fracture segment will be restricted and injection pressure will increase, as a result of stress shadow effects from hydraulic fracture segments and additional closing stresses from in-situ stress field. When multiple fractures propagate in naturally fracture reservoirs, complex fracture networks could be induced, which are affected by perforation cluster spacing, differential stress and natural fracture patterns. Combination of our numerical model and diagnostic methods (e.g. Microseismicity, DTS and DAS) is an effective approach to accurately characterize the complex fracture geometry. Furthermore, the physics-based complex fracture geometry provided by our model can be imported into reservoir simulation models for production analysis. / text