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

Unconventional reservoirs are typically characterized by very low permeabilities, and thus, the pressure depletion from a producing well may not propagate far from the well during the life of a development. Currently, two approaches are widely utilized to perform unconventional reservoir analysis: analytical techniques, including the decline curve analysis and the pressure/rate transient analysis, and numerical simulation. The numerical simulation can rigorously account for complex well geometry and reservoir heterogeneity but also is time consuming. In this thesis, we propose and apply an efficient technique, fast marching method (FMM), to analyze the shale gas reservoirs. Our proposed approach stands midway between analytic techniques and numerical simulation. In contrast to analytical techniques, it takes into account complex well geometry and reservoir heterogeneity, and it is less time consuming compared to numerical simulation. The fast marching method can efficiently provide us with the solution of the pressure front propagation equation, which can be expressed as an Eikonal equation. Our approach is based on the generalization of the concept of depth of investigation. Its application to unconventional reservoirs can provide the understanding necessary to describe and optimize the interaction between complex multi-stage fractured wells, reservoir heterogeneity, drainage volumes, pressure depletion, and well rates. The proposed method allows rapid approximation of reservoir simulation results without resorting to detailed flow simulation, and also provides the time-evolution of the well drainage volume for visualization. Calibration of reservoir models to match historical dynamic data is necessary to increase confidence in simulation models and also minimize risks in decision making. In this thesis, we propose an integrated workflow: applying the genetic algorithm (GA) to calibrate the model parameters, and utilizing the fast marching based approach for forward simulation. This workflow takes advantages of both the derivative free characteristics of GA and the speed of FMM. In addition, we also provide a novel approach to incorporate the micro-seismic events (if available) into our history matching workflow so as to further constrain and better calibrate our models.

Modern reservoir management typically involves simulations of geological models to predict future recovery estimates, providing the economic assessment of different field development strategies. Integrating reservoir data is a vital step in developing reliable reservoir performance models. Currently, most effective strategies for traditional manual history matching commonly follow a structured approach with a sequence of adjustments from global to regional parameters, followed by local changes in model properties. In contrast, many of the recent automatic history matching methods utilize parameter sensitivities or gradients to directly update the fine-scale reservoir properties, often ignoring geological inconsistency. Therefore, there is need for combining elements of all of these scales in a seamless manner. We present a hierarchical streamline-assisted history matching, with a framework of global-local updates. A probabilistic approach, consisting of design of experiments, response surface methodology and the genetic algorithm, is used to understand the uncertainty in the large-scale static and dynamic parameters. This global update step is followed by a streamline-based model calibration for high resolution reservoir heterogeneity. This local update step assimilates dynamic production data. We apply the genetic global calibration to unconventional shale gas reservoir specifically we include stimulated reservoir volume as a constraint term in the data integration to improve history matching and reduce prediction uncertainty. We introduce a novel approach for efficiently computing well drainage volumes for shale gas wells with multistage fractures and fracture clusters, and we will filter stochastic shale gas reservoir models by comparing the computed drainage volume with the measured SRV within specified confidence limits. Finally, we demonstrate the value of integrating downhole temperature measurements as coarse-scale constraint during streamline-based history matching of dynamic production data. We first derive coarse-scale permeability trends in the reservoir from temperature data. The coarse information are then downscaled into fine scale permeability by sequential Gaussian simulation with block kriging, and updated by local-scale streamline-based history matching. he power and utility of our approaches have been demonstrated using both synthetic and field examples.

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.