Predicting long-term production from gas shale reservoirs has been a major challenge for the petroleum industry. To better understand how production profiles are likely to evolve with time, we have conducted laboratory experiments examining the effects of confining stress and pore pressure on permeability. Experiments were conducted on intact core samples from the Barnett, Eagle Ford, Marcellus, and Montney shale reservoirs. The methodology used to measure permeability allows us to separate the reduction of permeability with depletion (because of the resultant increase in effective confining stress) and the increase in permeability associated with Knudsen diffusion and molecular slippage (also known as Klinkenberg) effects at very low pore pressure. By separating these effects, we are able to estimate the relative contribution of both Darcy and diffusive fluxes to total flow in depleted reservoirs. Our data show that the effective permeability of the rock is significantly enhanced at very low pore pressures ( []) because of the slippage effects. We use the magnitude of the Klinkenberg effect to estimate the effective aperture of the flow paths within the samples and compare these estimates to scanning electron microscopy image observations. Our results suggest effective flow paths to be on the order from tens of nanometers in most samples to 100–200 nm, in a relatively high-permeability Eagle Ford sample. Finally, to gain insight on the scale dependence of permeability measurements, the same core plugs were crushed, and permeability was again measured at the particle scale using the so-called Gas Research Institute method. The results show much lower permeability than the intact core samples, with very little correlation to the measurements on the larger scale cores.

Experimental investigation of matrix permeability of gas shales

Geological controls on matrix permeability of Devonian Gas Shales in the Horn River and Liard basins, northeastern British Columbia, Canada

Experimental study and modelling of methane adsorption and diffusion in shale

On the basis of microand mesoscale investigations, a new mathematical formulation is introduced in detail to investigate multiscale gas-transport phenomena in organic-rich-shale core samples. The formulation includes dual-porosity continua, where shale permeability is associated with inorganic matrix with relatively large irregularly shaped pores and fractures, whereas molecular phenomena (diffusive transport and nonlinear sorption) are associated with the kerogen pores. Kerogen is considered a nanoporous organic material finely dispersed within the inorganic matrix. The formulation is used to model and history match gaspermeation measurements in the laboratory using shale core plugs under confining stress. The results indicate significance of molecular transport and strong transient effects caused by gas/solid interactions within the kerogen. In the second part of the paper, we present a novel multiscale perturbation approach to quantify the overall impact of local porosity fluctuations associated with a spatially nonuniform kerogen distribution on the adsorption and transport in shale gas reservoirs. Adopting weak-noise and meanfield approximation, the approach applies a stochastic upscaling technique to the mathematical formulation developed in the first part for the laboratory. It allows us to investigate local kerogenheterogeneity effects in spectral (Fourier-Laplace) domain and to obtain an upscaled “macroscopic” model, which consists of the local heterogeneity effects in the real time—space domain. The new upscaled formulation is compared numerically with the previous homogeneous case using finite-difference approximations to initial/boundary value problems simulating the matrix gas release. We show that macrotransport and macrokinetics effects of kerogen heterogeneity are nontrivial and affect cumulative gas recovery. The work is important and timely for development of newgeneration shale-gas reservoir-flow simulators, and it can be used in the laboratory for organic-rich gas-shale characterization.

Multiscale Gas Transport in Shales With Local Kerogen Heterogeneities

Various attempts have been made to model flow in shale gas systems. However, there is currently little consensus regarding the impact of molecular and Knudsen diffusion on flow behavior over time in such systems. Direct measurement or model-based estimation of matrix permeability for these “ultra-tight” reservoirs has proven unreliable. The composition of gas produced from tight gas and shale gas reservoirs varies with time for a variety of reasons. The cause of flowing gas compositional change typically cited is selective desorption of gases from the surface of the kerogen in the case of shale. However, other drivers for gas fractionation are important when pore throat dimensions are small enough. Pore throat diameters on the order of molecular mean free path lengths will create non-Darcy flow conditions, where permeability becomes a strong function of pressure. At the low permeabilities found in shale gas systems, the dusty-gas model for flow should be used, which couples diffusion to advective flow. In this study we implement the dusty-gas model into a fluid flow modeling tool based on the TOUGH+ family of codes. We examine the effects of Knudsen diffusion on gas composition in ultra-tight rock. We show that for very small average pore throat diameters, lighter gases are preferentially produced at concentrations significantly higher than in situ conditions. Furthermore, we illustrate a methodology which uses measurements of gas composition to more uniquely determine the permeability of tight reservoirs. We also describe how gas composition measurement could be used to identify flow boundaries in these reservoir systems. We discuss how new measurement techniques and data collection practices should be implemented in order to take advantage of this method. Our contributions include a new, fit-for-purpose numerical model based on the TOUGH+ code capable of characterizing transport effects including permeability adjustment and diffusion in micro- and nano-scale porous media.

A numerical study of microscale flow behavior in tight gas and shale gas reservoir systems

In this paper we analyze by means of numerical simulation the mechanisms and processes of flow in two types of fractured tight gas reservoirs: shale and tight-sand systems. The numerical model includes Darcy’s law as the basic equation of multiphase flow and accurately describes the thermophysical properties of the reservoir fluids, but also incorporates other options that cover the spectrum of known physics that may be involved: non-Darcy flow, as described by a multi-phase extension of the Forschheimer equation that accounts for laminar, inertial and turbulent effects; stress-sensitive flow properties of the matrix and of the fractures, i.e., porosity, permeability, relative permeability and capillary pressure; gas slippage (Klinkenberg) effects; and, non-isothermal effects, accounting for the consequences of energy balance and temperature changes in the presence of phenomena such as Joule-Thompson cooling in the course of gas production. The flow and storage behavior of the fractured media (shale or tight sand) is represented by various options of the Multiple Interactive Continua (MINC) conceptual model, in addition to an Effective Continuum Method (ECM) option, and includes a gas sorption term that follows the Langmuir isotherm. Comparison to field data, analysis of the simulation results and parameter determination through history matching indicates that (a) the ECM model is incapable of describing the fractured system behavior, and (b) shale and tight-sand reservoirs exhibit different behavior that can be captured (albeit imperfectly) using some of the more complex options of the multi-continua fractured-system models. The sorption term is necessary to describe the behavior of shale gas reservoirs, and significant deviations from the field data are observed if it is omitted. Conversely, production data from tight-sand reservoirs can be adequately represented without accounting for gas sorption. All the other processes and mechanisms allow refinement of the match between predictions and observations, but appear to have secondorder effects in the description of flow through fractured tight gas reservoirs.

Analysis of Mechanisms of Flow in Fractured Tight-Gas and Shale-Gas Reservoirs

A numerical study of performance for tight gas and shale gas reservoir systems

High-Resolution Numerical Modeling of Complex and Irregular Fracture Patterns in Shale-Gas Reservoirs and Tight Gas Reservoirs

Various analytical, semi-analytical, and empirical models have been proposed to characterize rate and pressure behavior as a function of time in tight gas and shale gas systems featuring horizontal wells with multiple hydraulic fractures. Despite a few analytical models, as well as a small number of published numerical studies, there is currently little consensus regarding the large-scale flow behavior over time in such systems, particularly regarding the dominant flow regimes and whether or not reservoir properties or volumes can be estimated from well performance data. We constructed a fit-for-purpose numerical simulator which accounts for a variety of production features pertinent to these systems—specifically ultra-tight matrix permeability, hydraulically fractured horizontal wells with induced fractures of various configurations, multiple porosity and permeability fields, and desorption. These features cover the production mechanisms which are currently believed to be most relevant in tight gas and shale gas systems. We employ the numerical simulator to examine various tight gas and shale gas systems and to identify and illustrate the various flow regimes which progressively occur over time. We perform this study at fine grid discretization on the order of 1 mm near fractures to accurately capture flow effects at all time periods. We visualize the flow regimes using specialized plots of rate and pressure functions, as well as maps of pressure and sorption distributions. We use pressure maps to visualize the various flow regimes and their transitions in tight gas systems. In a typical tight gas system, we illustrate the initial linear flow into the hydraulic fractures (i.e., formation linear flow ), transitioning to compound formation linear flow , and eventually transforming into elliptical flow . We explore variations of possible shale gas system models. Based on diffusive flow (with and without desorption), we show that due to the extremely low permeability of shale matrix (a few nanodarcies), the flow behavior is dominated by the extent of and configuration of the fractures. This work expands our understanding of flow behavior in tight gas and shale gas systems, where such an understanding may ultimately be used to estimate reservoir properties and reserves in these types of reservoirs.

Craig M. Freeman

Papers

A numerical study of microscale flow behavior in tight gas and shale gas reservoir systems

Analysis of Mechanisms of Flow in Fractured Tight-Gas and Shale-Gas Reservoirs

A numerical study of performance for tight gas and shale gas reservoir systems

High-Resolution Numerical Modeling of Complex and Irregular Fracture Patterns in Shale-Gas Reservoirs and Tight Gas Reservoirs

A Numerical Study of Performance for Tight Gas and Shale Gas Reservoir Systems