Unconventional natural gas extraction from the Marcellus Shale requires millions of gallons of water to fracture shale and release natural gas from the formation. This process produces water with high levels of total dissolved solids (TDS); and, efforts to recycle these fluids has stimulated microbial growth in produced water. The objective of this study was to analyze the ionic composition of and characterize microorganisms from Marcellus produced water samples. A semi-synthetic culture medium was designed with high TDS to enrich for halophilic microbes, which yielded robust cultures that were able to grow over a wide range of salinities. DNA extracted from aerobic cultures was used for 16s rDNA clone libraries and Automated Ribosomal Intergenic Spacer Analysis (ARISA). ARISA and 16S gene sequencing revealed differences in bacterial composition between Marcellus and freshwater samples. Sequencing of 16S gene indicated the presence of Halomonas, Thalassospira and other genera related to halophilic and petroleum degrading species. / Bayer School of Natural and Environmental Sciences / Environmental Science and Management (ESM) / MS / Thesis

Time-rate analysis and time-rate-pressure analysis methods are available to estimate reserves and study flow performance of wells in unconventional gas reservoirs. However, these tools are often incorrectly used or the analysis can become difficult because of the complex nature of the reservoir system. Conventional methods (e.g., Arps' time-rate relations) are often used incorrectly to estimate reserves from such reservoirs. It was only recently that a serious study was conducted to outline the limitations of these relations and to set guidelines for their correct application. New time-rate relations, particularly the Duong and logistic growth model, were introduced to estimate reserves and forecast production from unconventional reservoirs. These new models are being used with limited understanding of their characteristics and limitations. Moreover, well performance analyses using analytical/semi-analytical solutions (time-rate-pressure) are often complicated from non-uniqueness that arises when estimating well/formation properties. In this work, we present a detailed study of the Duong model and logistic growth model to investigate the behaviors and limitations of these models when analyzing production data from unconventional reservoirs. We consider production data generated from numerical simulation cases and data obtained from unconventional gas reservoirs to study the quality of match to specific flow regimes and compare accuracy of the reserve estimates. We use the power-law exponential model (PLE), which has been shown to model transient, transition and boundary-dominated flow regimes reliably, as a benchmark to study performance of Duong and logistic growth models. Moreover, we use the "continuous EUR" approach to compare these models during reserve estimation. Finally, we develop four new time-rate relations, based on characteristics of the time-rate data on diagnostic plots. Using diagnostic plots we show that the new time-rate relations provide a quality match to the production data across all flow regimes, leading to a reliable reserve estimate. In a preliminary study, we integrated time-rate model parameters with fundamental reservoir properties (i.e., fracture conductivity (Fc) and 30 year EUR (EUR30yr)), by studying 15 numerical simulation cases to yield parametric correlations. We have demonstrated a methodology to integrate time-rate model parameters and reservoir properties. This method avoids the non-uniqueness issues often associated with model-based production data analysis. This study provides theoretical basis for further demonstration of the methodology using field cases.

Fractured reservoirs have been traditionally idealized using dual-porosity models. In these models, all matrix and fractures systems have identical properties. However, it is not uncommon for naturally fractured reservoirs to have orthogonal fractures with different properties. In addition, for hydraulically fractured reservoirs that have preexisting natural fractures such as shale gas reservoirs, it is almost certain that these types of fractures are present. Therefore, a triple-porosity (dual-fracture) model is developed in this work for characterizing fractured reservoirs with different fractures properties. The model consists of three contiguous porous media: the matrix, less permeable micro-fractures and more permeable macro-fractures. Only the macro-fractures produce to the well while they are fed by the micro-fractures only. Consequently, the matrix feeds the micro-fractures only. Therefore, the flow is sequential from one medium to the other. Four sub-models are derived based on the interporosity flow assumption between adjacent media, i.e., pseudosteady state or transient flow assumption. These are fully transient flow model (Model 1), fully pseudosteady state flow model (Model 4) and two mixed flow models (Model 2 and 3). The solutions were mainly derived for linear flow which makes this model the first triple-porosity model for linear reservoirs. In addition, the Laplace domain solutions are also new and have not been presented in the literature before in this form. Model 1 is used to analyze fractured shale gas horizontal wells. Non-linear regression using least absolute value method is used to match field data, mainly gas rate. Once a match is achieved, the well model is completely described. Consequently, original gas in place (OGIP) can be estimated and well future performance can be forecasted.

Production forecasting in shale (ultra-low permeability) gas reservoirs is of great interest due to the advent of multi-stage fracturing and horizontal drilling. The well renowned production forecasting model, Arps? Hyperbolic Decline Model, is widely used in industry to forecast shale gas wells. Left unconstrained, the model often overestimates reserves by a great deal. A minimum decline rate is imposed to prevent overestimation of reserves but with less than ten years of production history available to analyze, an accurate minimum decline rate is currently unknown; an educated guess of 5% minimum decline is often imposed. Other decline curve models have been proposed with the theoretical advantage of being able to match linear flow followed by a transition to boundary dominated flow. This thesis investigates the applicability of the Stretched Exponential Production Decline Model (SEPD) and compares it to the industry standard, Arps' with a minimum decline rate. When possible, we investigate an SEPD type curve. Simulated data is analyzed to show advantages of the SEPD model and provide a comparison to Arps' model with an imposed minimum decline rate of 5% where the full production history is known. Long-term production behavior is provided by an analytical solution for a homogenous reservoir with homogenous hydraulic fractures. Various simulations from short-term linear flow (~1 year) to long-term linear flow (~20 years) show the ability of the models to handle onset of boundary dominated flow at various times during production history. SEPD provides more accurate reserves estimates when linear flow ends at 5 years or earlier. Both models provide sufficient reserves estimates for longer-term linear flow scenarios. Barnett Shale production data demonstrates the ability of the models to forecast field data. Denton and Tarrant County wells are analyzed as groups and individually. SEPD type curves generated with 2004 well groups provide forecasts for wells drilled in subsequent years. This study suggests a type curve is most useful when 24 months or less is available to forecast. The SEPD model generally provides more conservative forecasts and EUR estimates than Arps' model with a minimum decline rate of 5%.

Horizontal drilling and hydraulic stimulation make hydrocarbon production from organic-rich shales economically viable. Identification of suitable zones to drill a horizontal well and to initiate or contain hydraulic fractures requires detection and quantification of many factors, including elastic mechanical properties. Elastic behavior of rocks is affected by rock composition and fabric, pore pressure, confining stress, and other factors. Rock fabric refers to the arrangement of the rock’s solid and fluid constituents. The objective of this thesis is to quantify rock fabric properties of hydrocarbon-bearing shales affecting elastic properties, including load-bearing matrix, anisotropic cracks, and shape of rock components. Once rock fabric is validated with sonic logs, results can be used to identify suitable zones to drill a horizontal well, initiate hydraulic stimulation, and contain fracture propagation. We develop a method to estimate elastic properties based on rock composition. The differential effective medium (DEM) theory is invoked to model rock elastic properties with a load-bearing component in which remaining minerals and pores are added as spheres or ellipsoids. The method can be combined with the self-consistent approximation (SCA) to construct a load-bearing matrix made of two solid phases. Anisotropic inclusions are added via Hudson’s model. Subsequently, Gassmann’s theory is invoked to saturate the rock with fluids and determine low-frequency elastic properties for comparison to sonic logs. Rock fabric properties remain constant in a vertically homogeneous formation. In vertically heterogeneous strata, the depth interval of interest is divided into rock types, based on rock solid composition, and each rock type is associated with a specific fabric. Quantification of rock fabric properties is a non-unique process, and one should take into account as much petrophysical and geological information as possible to ensure physically viable results. Our simulation and interpretation method is implemented in two wells in both the Haynesville and Barnett shales. Averages of relative errors between estimated velocities and sonic logs are less than 4% in the four wells. Simulations in the Haynesville shale are isotropic, and therefore indicate that rock fabric may not be the main cause of mechanical anisotropy in cases where such behavior is inferred from field data. Rock fabric properties are constant with depth in both wells. Consequently, identification of suitable zones to drill a horizontal well or to contain fracture propagation is not based on rock fabric; it is deduced from Young’s modulus. Simulated Poisson’s ratio is shown to be more sensitive to errors in velocities than Young’s modulus and is therefore not used in the interpretation. Favorable depth intervals for gas production exhibit sizeable volumes of gas and organic content. In the Barnett shale, the two wells exhibit different rock fabrics. Such a behavior indicates that the formation is laterally heterogeneous. Rock physics models should therefore be extrapolated from one well to another with caution. Simulations assume anisotropic elastic behavior and suggest the presence of compliant horizontal pores in one case. Natural vertical fractures are observed on electric image logs in the remaining case and are modeled with Hudson’s theory. This behavior suggests that rock fabric causes mechanical anisotropy in the formation. Heterogeneity of the Barnett shale rock fabric is inferred from the necessary use of rock typing to adequately reproduce sonic logs in both wells. Intervals with large porosity and high gas saturation identify suitable zones to perform hydraulic stimulation. Among such zones, rock types that exhibit stiff load-bearing matrices (comprising mostly calcite, for example) indicate suitable depths to drill horizontal wells or to contain hydraulic fractures. Intervals with dense layering of different rock types are unsuitable for fracture propagation and should be avoided during hydraulic-fracturing operations. / text

The vast shale gas and tight oil reservoirs in North America cannot be economically developed without multi-stage hydraulic fracture treatments. Owing to the disparity in the density of natural fractures in addition to the disparate in-situ stress conditions in these kinds of formations, microseismic fracture mapping has shown that hydraulic fracture treatments develop a range of large-scale fracture networks in the shale plays. In this thesis, an approach is presented, where the fracture networks approximated with microseismic mapping are integrated with a commercial numerical production simulator that discretely models the network structure in both vertical and horizontal wells. A novel approach for reservoir simulation is used, where porosity (instead of permeability) is used as a scaling parameter for the fracture width. Two different fracture geometries have been broadly proposed for a multi stage horizontal well, orthogonal and transverse. The orthogonal pattern represents a complex network with cross cutting fractures orthogonal to each other; whereas transverse pattern maps uninterrupted fractures achieving maximum depth of penetration into the reservoir. The response for a vii single-stage fracture is further investigated by comparing the propagation of the stage to be dendritic versus planar. A dendritic propagation is bifurcation of the hydraulic fracture due to intersection with the natural fracture (failure along the plane of weakness). The impact of fracture spacing to optimize these fracture geometries is studied. A systematic optimization for designing the fracture length and width is also presented. The simulation is motivated by the oil window of Eagle Ford shale formation and the results of this work illustrate how different fracture network geometries impact well performance, which is critical for improving future horizontal well completions and fracturing strategies in low permeability shale and tight oil reservoirs. A rate transient analysis (RTA) technique employing a rate normalized pressure (RNP) vs. superposition time function (STF) plot is used for the linear flow analysis. The parameters that influence linear flow are analytically derived. It is found that picking a straight line on this curve can lead to erroneous results because multiple solutions exist. A new technique for linear flow analysis is used. The ratio of derivative of inverse production and derivative of square root time is plotted against square root time and the constant derivative region is seen to be indicative of linear flow. The analysis is found to be robust because different simulation cases are modeled and permeability and fracture half-length are estimated. / text

After a fracturing job in a shale reservoir, only a fraction of injected water is recovered. Water is trapped inside the reservoir and reduces the relative permeability of gas. By reducing the interfacial tension between water and hydrocarbon, more water can be recovered thus increasing overall gas production. By adding surfactants into the fracturing fluid, the IFT can be reduced and will help mobilize trapped water. From previous research, two types of surfactant have been identified to be CO₂ soluble. These are the ethoxylated tallow amine and ethoxylated coco amine with varying ethoxylate length. Experiments were performed to test the solubility of these surfactants in water, observe how they change the interaction between HC and water, and measure the IFT reduction between HC and water. Surfactants with more than 10 EO groups were soluble at all salinities, temperature and pH. They also form a non-typical water-in-oil emulsion at all salinities. The surfactants, Ethomeen T/25, T/30, C/15, and C/25 were used in the IFT measurements. They showed interesting trends that exhibit their hydrophilic/hydrophobic nature. These surfactants reduce the IFT between pentane and water to approximately 5 mN/m. The results show that these surfactants do reduce the IFT between water and hydrocarbon, but not as well as conventional EOR surfactants. They do have other added benefits such as being CO₂ soluble, form water in oil emulsions, and tolerant to high temperature and salinity. / text

Recent optimistic revisions to projections for recoverable natural gas resources in the United States have generated renewed interest in the possibility of greater utilization of natural gas as a transportation fuel. Against a backdrop of significant policy challenges for the United States, including air quality concerns in urban areas, slow economic growth and high unemployment, and a rising unease with regard to an increasing dependence on foreign oil; natural gas offers the nation’s transportation sector an opportunity to reduce mobile emissions, lower fuel costs, create jobs and reduce dependence on imported oil. While the current focus for expanded use of natural gas in the transportation sector emphasizes heavy duty and fleet vehicles, there may also be potential for increased use for passenger vehicles. Inconvenience, with regard to refueling, and high incremental vehicle costs, however, are seen as major obstacles to greater adaptation. This analysis examines the benefits and drawbacks of natural gas vehicles from the passenger vehicle perspective and includes data from a cross-country road trip. The report includes a review of market trends and possible development scenarios and concludes with recommendations to minimize the potential challenges of greater adaptation of natural gas vehicles in the passenger vehicle market. / text

An investigation of shale pore structure and compositional/geochemical heterogeneities has been undertaken to elucidate the controls upon gas capacities of potential shale gas reservoirs in northeastern British Columbia, western Canada. Methane sorption isotherms, pore structure and surface area data indicate a complex interrelationship of total organic carbon (TOC) content, mineral matter and thermal maturity affect gas sorption characteristics of Devonian-€”Mississippian (D-€”M) and Jurassic strata. Methane and carbon dioxide sorption capacities of D-€”M shales increase with TOC content, due to the microporous nature of the organic matter. Clay mineral phases arealso capable of sorbing gas to their internal structure; hence D-€”M shales which are both TOC- and clay-rich have the largest micropore volumes and sorption capacities on a dry basis. Jurassic shales, which are invariably less thermally mature than D-€”M shales, do not have micropore volumes which correlate with TOC. The covariance of methane sorption capacity with TOC, independent of micropore volume, indicates a solute gas contribution (within matrix bituminite) to the total gas capacity. On a wt% TOC basis, D-M shales sorb more gas than Jurassic shales: a result of thermal-maturation induced, structural transformation of the D-€”M organic fraction. Organic-rich D-€”M strata are considered to be excellent candidates for gas shales in Western Canada. These strata have TOC contents ranging between 1-5.7 wt%, thermal maturities into the dry-gas region, and thicknesses in places of over 1000 m. Total gas capacity estimates range between 60 and 600 bcf/section where a substantial percentage of the gas capacity is free gas, due to high reservoir temperatures and pressures. Inorganic material influences modal pore size, total porosity and sorption characteristics of D-M shales. Carbonate-rich samples often have lower organic carbon contents (oxic deposition) and porosity, hence potentially lower sorbed and free-gas capacities. Highly mature Devonian shales are both silica and TOC-rich (up to 85% quartz and 5 wt% TOC) and as such, deemed excellent potential shale gas reservoirs because they are both brittle (fracable), and gas-charged. However, quartz-rich Devonian shales display tight-rock characteristics, with poorly developed fabric, small median pore diameters and low permeabilities. Hence potential `frac-zones' will require an increased density of hydraulic fracture networks for optimum gas production.

Excess water production is a major economic and environmental problem for the oil and gas industry. The cost of processing excess water runs into billions of dollars. Polymer gel technology has been successfully used in controlling water influx without damaging hydrocarbon production in conventional naturally fractured or hydraulically fractured reservoirs. However, there has been no systematic investigation on effectiveness and placement conditions of polymer gels for shutting off water flow from fractures with narrow apertures in shale and tight gas reservoirs. The existing polymer gels, like those based on Chromium(III) Acetate, as a crosslinker will exert very high extrusion pressure to effectively penetrate the narrow aperture fractures present in shale and tight gas reservoirs. This gives rise to a need for a new polymer gel system that can be used for selectively shutting off water flow from narrow aperture fractures in shale and tight gas reservoirs. The new gel system will have a longer gelation time than the existing polymer gels; this ensures minimum crosslinking of the gel by the time it reaches bottom hole. The gelant solution will be pumped at low pressure so that, it penetrates only pre-existing fractures in the formation with ease. This study for the first time focuses on developing an environmentally benign polymer gel system based on high molecular weight HPAM, as a base polymer and a commercial grade PEI as an organic crosslinker. Gel samples of different concentration ratios of the polymer and crosslinker were prepared and classified under Sydansk code of gel strength to find optimum concentration ratios that gave good gels. The gel system was characterized using Brookfield DV-III Ultra Rheometer and Fann-35 Viscometer.