The void space properties of a set of gas reservoir sandstone samples have been measured. The properties include porosity, absolute gas permeability, electrical resistivity formation factor and tortuosity. The mineralogy of each sandstone was determined by scanning electron microscopy and energy dispersive x-ray analysis. Mercury intrusion and extrusion data have been measured for most of the sandstone samples. A new procedure for measuring the degree and range of void size correlations within resin-filled sandstones has been developed. Image analysis of backscattered electron micrographs of these samples supplies void size and positional information. A "semi-variogram" study of void size and coordinate data ascertains the degree and range of void size correlation. Measurable correlation has been found in two sandstone samples, but was absent from four others. Diffusion coefficients of methane, iso-butane and n-butane through dry sandstones have been measured using an adaptation of a non-steady state method, using a redesigned apparatus. A repeatability and error analysis of diffusion coefficient measurement has also been performed. A correlation between diffusion coefficients, absolute permeability, porosity and formation factor was detected for sandstones containing little clay. The diffusion coefficients measured for clay affected sandstones did not correlate with any petrophysical properties of these samples. A computer model capable of simulating porous media has been previously developed. It consists of a 10x10x10 network of cubic pores and cylindrical throats, and simulates die mercury intrusion curve. The void size distribution is modified until both simulated and experimental curves closely match. New void size distribution input and curve fit algorithms have been developed to increase the speed and accuracy of die simulations and a new modelling procedure allows the modelling of samples with void size correlation. The model is capable of simulating porosity, permeability, tortuosity and mercury extrusion. Each of the reservoir sandstones has been modelled and their characteristic properties simulated. Successful simulations were obtained for all relatively clay-free reservoir sandstones. Clay affected sandstone simulations were less successful due to the high complexity of these samples. A study into formation damage witiiin reservoir sandstones was also undertaken. The effect of colloidal particulate void space penetration is measured and simulated.

Accurately estimating hydrocarbon reserves is important, because it affects every phase of the oil and gas business. Unfortunately, reserves estimation is always uncertain, since perfect information is seldom available from the reservoir, and uncertainty can complicate the decision-making process. Many important decisions have to be made without knowing exactly what the ultimate outcome will be from a decision made today. Thus, quantifying the uncertainty is extremely important. Two methods for estimating original hydrocarbons in place (OHIP) are volumetric and material balance methods. The volumetric method is convenient to calculate OHIP during the early development period, while the material balance method can be used later, after performance data, such as pressure and production data, are available. In this work, I propose a methodology for using a Bayesian approach to quantify the uncertainty of original gas in place (G), aquifer productivity index (J), and the volume of the aquifer (Wi) as a result of combining volumetric and material balance analysis in a water-driven gas reservoir. The results show that we potentially have significant non-uniqueness (i.e., large uncertainty) when we consider only volumetric analyses or material balance analyses. By combining the results from both analyses, the non-uniqueness can be reduced, resulting in OGIP and aquifer parameter estimates with lower uncertainty. By understanding the uncertainty, we can expect better management decision making.

he feasibility of sequestering supercritical CO2 in depleted gas reservoirs. The experimental runs involved the following steps. First, the 1 ft long by 1 in. diameter carbonate core is inserted into a viton Hassler sleeve and placed inside an aluminum coreholder that is then evacuated. Second, with or without connate water, the carbonate core is saturated with methane. Third, supercritical CO2 is injected into the core with 300 psi overburden pressure. From the volume and composition of the produced gas measured by a wet test meter and a gas chromatograph, the recovery of methane at CO2 breakthrough is determined. The core is scanned three times during an experimental run to determine core porosity and fluid saturation profile: at start of the run, at CO2 breakthrough, and at the end of the run. Runs were made with various temperatures, 20°C (68°F) to 80°C (176°F), while the cell pressure is varied, from 500 psig (3.55 MPa) to 3000 psig (20.79 MPa) for each temperature. An analytical study of the experimental results has been also conducted to determine the dispersion coefficient of CO2 using the convection-dispersion equation. The dispersion coefficient of CO2 in methane is found to be relatively low, 0.01-0.3 cm2/min.. Based on experimental and analytical results, a 3D simulation model of one eighth of a 5-spot pattern was constructed to evaluate injection of supercritical CO2 under typical field conditions. The depleted gas reservoir is repressurized by CO2 injection from 500 psi to its initial pressure 3,045 psi. Simulation results for 400 bbl/d CO2 injection may be summarized as follows. First, a large amount of CO2 is sequestered: (i) about 1.2 million tons in 29 years (0 % initial water saturation) to 0.78 million tons in 19 years (35 % initial water saturation) for 40-acre pattern, (ii) about 4.8 million tons in 112 years (0 % initial water saturation) to 3.1 million tons in 73 years (35 % initial water saturation) for 80-acre pattern. Second, a significant amount of natural gas is also produced: (i) about 1.2 BSCF or 74 % remaining GIP (0 % initial water saturation) to 0.78 BSCF or 66 % remaining GIP (35 % initial water saturation) for 40-acre pattern, (ii) about 4.5 BSCF or 64 % remaining GIP (0 % initial water saturation) to 2.97 BSCF or 62 % remaining GIP (35 % initial water saturation) for 80-acre pattern. This produced gas revenue could help defray the cost of CO2 sequestration. In short, CO2 sequestration in depleted gas reservoirs appears to be a win-win technology.

In the coming decades, the world will require additional supplies of natural gas to meet the demand for energy. Tight gas reservoirs can be defined as reservoirs where the formation permeability is so low (< 0.1 md) that advanced stimulation technologies, such as large volume fracture treatments, are required before a reasonable profit can be made. Hydraulic fracturing is one of the best methods to stimulate a tight gas well. Most fracture treatments result in 3-6 fold increases in the productivity index. However, if one computes the effective fracture length of most wells, we usually find that the effective length is less than the designed propped fracture length. The “propped length” is the distance down the fracture from the wellbore where proppants have been placed at a high enough concentration to “prop open” the fracture. The “effective length” is the portion of the propped fracture that cleans up and allows gas flow from the reservoir into the fracture then down the fracture to the wellbore. Whenever the effective length is much shorter than the designed propped length, several reasons must be evaluated to determine what might have occurred. For example, the difference could be caused by one or more of the following issues: insufficient fracture fluid cleanup, proppant settling, proppant embedment, proppant crushing, or poor reservoir continuity. Although all these causes are possible, we believe that fracture fluid cleanup issues may be the most common reason the industry fails to achieve the designed propped fracture length in most cases. In this research, we have investigated fracture fluid cleanup problems and developed a better understanding of the issues involved which hopefully will lead to ways to improve cleanup. Fracture fluid cleanup is a complex problem, that can be influenced by many parameters such as the fluid system used, treatment design, flowback procedures, production strategy, and reservoir conditions. Residual polymer in the fracture can reduce the effective fracture permeability and porosity, reduce the effective fracture half-length, and limit the well productivity. Our ability to mathematically model the fundamental physical processes governing fluid recovery in hydraulic fractures in the past has been limited. In this research, fracture fluid damage mechanisms have been investigated, and mathematical models and computer codes have been developed to better characterize the cleanup process. The codes have been linked to a 3D, 3-phase simulator to model and quantify the fracture fluid cleanup process and its effect on long-term gas production performances. Then, a comprehensive systematic simulation study has been carried out by varying formation permeability, reservoir pressure, fracture length, fracture conductivity, yield stress, and pressure drawdown. On the basis of simulation results and analyses, new ways to improve fracture fluid cleanup have been provided. This new progress help engineers better understand fracture fluid cleanup, improve fracture treatment design, and increase gas recovery from tight sand reservoirs, which can be extremely important as more tight gas reservoirs are developed around the world.

Optimizing well spacing in unconventional gas reservoirs is difficult due to complex heterogeneity, large variability and uncertainty in reservoir properties, and lack of data that increase the production uncertainty. Previous methods are either suboptimal because they do not consider subsurface uncertainty (e.g., statistical moving-window methods) or they are too time-consuming and expensive for many operators (e.g., integrated reservoir characterization and simulation studies). This research has focused on developing and extending a new technology for determining optimal well spacing in tight gas reservoirs that maximize profitability. To achieve the research objectives, an integrated multi-well reservoir and decision model that fully incorporates uncertainty was developed. The reservoir model is based on reservoir simulation technology coupled with geostatistical and Monte Carlo methods to predict production performance in unconventional gas reservoirs as a function of well spacing and different development scenarios. The variability in discounted cumulative production was used for direct integration of the reservoir model with a Bayesian decision model (developed by other members of the research team) that determines the optimal well spacing and hence the optimal development strategy. The integrated model includes two development stages with a varying Stage-1 time span. The integrated tools were applied to an illustrative example in Deep Basin (Gething D) tight gas sands in Alberta, Canada, to determine optimal development strategies. The results showed that a Stage-1 length of 1 year starting at 160-acre spacing with no further downspacing is the optimal development policy. It also showed that extending the duration of Stage 1 beyond one year does not represent an economic benefit. These results are specific to the Berland River (Gething) area and should not be generalized to other unconventional gas reservoirs. However, the proposed technology provides insight into both the value of information and the ability to incorporate learning in a dynamic development strategy. The new technology is expected to help operators determine the combination of primary and secondary development policies early in the reservoir life that profitably maximize production and minimize the number of uneconomical wells. I anticipate that this methodology will be applicable to other tight and shale gas reservoirs.

One of the major challenges in fracturing low permeability/tight/unconventional gas formations is the loss of frac water and well productivity due to fluid entrapment in the matrix or fracture. Field results have indicated that only 15-30% of the frac fluid is recovered at the surface after flow back is initiated. Past studies have suggested that this water is trapped in the rock matrix near the fracture face and remains trapped due to the high capillary pressure in the matrix. Significant efforts have been made in the past to understand the impact of liquid blocking in hydraulically fractured conventional gas wells. Numerous remediation measures such as huff and puff gas cycling, alcohol or surfactant based chemical treatments have been proposed to reduce fracture face damage. However, when considering hydraulic fractures in unconventional reservoirs horizontal wells, the fluid may also be trapped within the fracture itself and may impact the cleanup as well as productivity. This study shows that under typical gas flow rates in tight / shale gas formations, liquid loading within the fractures is likely to occur. Most of the previous simulation studies consider a 2D reservoir model and ignore gravity, considering the high vertical anisotropy (or extremely low vertical permeability) in these tight reservoirs matrix. However, this study presents the results of 3D simulations of liquid loading in hydraulic fractures in horizontal wells, including gravity and capillary pressure effects. Both CMG IMEX and GEM have been used to study this phenomenon in dry and wet gas cases. The impact of drawdown, fracture and reservoir properties on liquid loading and well productivity is presented. Results show that low drawdown, low matrix permeability or low initial gas rates aggravate the liquid loading problem inside the fracture and thereby impact the cleanup and gas productivity during initial production. A clear understanding of the phenomena could help in selection of optimal production facilities and well profile. / text

he feasibility of sequestering supercritical CO2 in depleted gas reservoirs. The experimental runs involved the following steps. First, the 1 ft long by 1 in. diameter carbonate core is inserted into a viton Hassler sleeve and placed inside an aluminum coreholder that is then evacuated. Second, with or without connate water, the carbonate core is saturated with methane. Third, supercritical CO2 is injected into the core with 300 psi overburden pressure. From the volume and composition of the produced gas measured by a wet test meter and a gas chromatograph, the recovery of methane at CO2 breakthrough is determined. The core is scanned three times during an experimental run to determine core porosity and fluid saturation profile: at start of the run, at CO2 breakthrough, and at the end of the run. Runs were made with various temperatures, 20°C (68°F) to 80°C (176°F), while the cell pressure is varied, from 500 psig (3.55 MPa) to 3000 psig (20.79 MPa) for each temperature. An analytical study of the experimental results has been also conducted to determine the dispersion coefficient of CO2 using the convection-dispersion equation. The dispersion coefficient of CO2 in methane is found to be relatively low, 0.01-0.3 cm2/min.. Based on experimental and analytical results, a 3D simulation model of one eighth of a 5-spot pattern was constructed to evaluate injection of supercritical CO2 under typical field conditions. The depleted gas reservoir is repressurized by CO2 injection from 500 psi to its initial pressure 3,045 psi. Simulation results for 400 bbl/d CO2 injection may be summarized as follows. First, a large amount of CO2 is sequestered: (i) about 1.2 million tons in 29 years (0 % initial water saturation) to 0.78 million tons in 19 years (35 % initial water saturation) for 40-acre pattern, (ii) about 4.8 million tons in 112 years (0 % initial water saturation) to 3.1 million tons in 73 years (35 % initial water saturation) for 80-acre pattern. Second, a significant amount of natural gas is also produced: (i) about 1.2 BSCF or 74 % remaining GIP (0 % initial water saturation) to 0.78 BSCF or 66 % remaining GIP (35 % initial water saturation) for 40-acre pattern, (ii) about 4.5 BSCF or 64 % remaining GIP (0 % initial water saturation) to 2.97 BSCF or 62 % remaining GIP (35 % initial water saturation) for 80-acre pattern. This produced gas revenue could help defray the cost of CO2 sequestration. In short, CO2 sequestration in depleted gas reservoirs appears to be a win-win technology.

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.

Shale gas reservoirs have become a major source of energy in recent years. Developments in hydraulic fracturing technology have made these reservoirs more accessible and productive. Apart from other dissimilarities from conventional gas reservoirs, one major difference is that a considerable amount of gas produced from these reservoirs comes from desorption. Ignoring a major component of production, such as desorption, could result in significant errors in analysis of these wells. Therefore it is important to understand the adsorption phenomenon and to include its effect in order to avoid erroneous analysis. The objective of this work was to imbed the adsorbed gas in the techniques used previously for the analysis of tight gas reservoirs. Most of the desorption from shale gas reservoirs takes place in later time when there is considerable depletion of free gas and the well is undergoing boundary dominated flow (BDF). For that matter BDF methods, to estimate original gas in place (OGIP), that are presented in previous literature are reviewed to include adsorbed gas in them. More over end of the transient time data can also be used to estimate OGIP. Kings modified z* and Bumb and McKee’s adsorption compressibility factor for adsorbed gas are used in this work to include adsorption in the BDF and end of transient time methods. Employing a mass balance, including adsorbed gas, and the productivity index equation for BDF, a procedure is presented to analyze the decline trend when adsorbed gas is included. This procedure was programmed in EXCEL VBA named as shale gas PSS with adsorption (SGPA). SGPA is used for field data analysis to show the contribution of adsorbed gas during the life of the well and to apply the BDF methods to estimate OGIP with and without adsorbed gas. The estimated OGIP’s were than used to forecast future performance of wells with and without adsorption. OGIP estimation methods when applied on field data from selected wells showed that inclusion of adsorbed gas resulted in approximately 30 percent increase in OGIP estimates and 17 percent decrease in recovery factor (RF) estimates. This work also demonstrates that including adsorbed gas results in approximately 5percent less stimulated reservoir volume estimate.

It is essential to understand natural fracture systems embedded in shale-gas reservoirs and the stress fields that influence how induced fractures form in targeted shale units. Multicomponent seismic technology and elastic seismic stratigraphy allow geologic formations to be better images through analysis of different S-wave modes as well as the P-wave mode. Significant amounts of energy produced by P-wave sources radiate through the Earth as downgoing SV-wave energy. A vertical-force source is an effective source for direct SV radiation and provides a pure shear-wave mode (SV-SV) that should reveal crucial information about geologic surfaces located in anisotropic media. SV-SV shear wave modes should carry important information about petrophysical characteristics of hydrocarbon systems that cannot be obtained using other elastic-wave modes. Regardless of the difficulties of extracting good-quality SV-SV signal, direct shear waves as well as direct P and converted S energy should be accounted for in 3C seismic studies. Acquisition of full-azimuth seismic data and sampling data at small intervals over long offsets are required for detailed anisotropy analysis. If 3C3D data can be acquired with improved signal-to-noise ratio, more uniform illumination of targets, increased lateral resolution, more accurate amplitude attributes, and better multiple attenuation, such data will have strong interest by the industry. The objectives of this research are: (1) determine the feasibility of extracting direct SV-SV common-mid-point sections from 3-C seismic surveys, (2) improve the exploration for stratigraphic traps by developing systematic relationship between petrophysical properties and combinations of P and S wave modes, (3) create compelling examples illustrating how hydrocarbon-bearing reservoirs in low-permeable rocks (particularly anisotropic shale formations) can be better characterized using different S-wave modes (P-SV, SV-SV) in addition to the conventional P-P modes, and (4) analyze P and S radiation patterns produced by a variety of seismic sources. The research done in this study has contributed to understanding the physics involved in direct-S radiation from vertical-force source stations. A U.S. Patent issued to the Board of Regents of the University of Texas System now protects the intellectual property the Exploration Geophysics Laboratory has developed related to S-wave generation by vertical-force sources. The University’s Office of Technology Commercialization is actively engaged in commercializing this new S-wave reflection seismic technology on behalf of the Board of Regents. / text