Shale Properties
Tight/shale gas reservoirs have very low matrix permeability and in some instances contain organic-rich rocks where gas can also be adsorbed. (See Shale Adsorption Correction for additional information.) Production from these reservoirs requires long horizontal wells with multiple staged hydraulic fracture treatments to create extensive artificial fracture networks near the wellbore.
Estimation of original gas-in-place in tight/shale gas reservoirs requires:
- Conventional volumetric calculation using rock volume, porosity, fluid saturations, and pressure.
- Use of Langmuir Isotherms to estimate the volume of adsorbed gas.
Typically, all water in a tight/shale gas reservoir is not moveable and so reported water production is almost always attributable to production of load fluids (completion/fracture treatment fluid).

Conventional vs. Tight / Shale Gas Reservoirs
A comparison of conventional and tight / shale gas reservoirs is shown below:
| Characteristic |
Conventional |
Shale |
|---|---|---|
|
Gas Generation |
Gas is generated in the source rock and then migrates into the reservoir. |
Gas is generated and trapped within the source rock. |
|
Gas Storage Mechanism |
Compression. |
Compression and adsorption. |
|
Gas Produced |
Free gas only. |
Free and adsorbed gas. |
|
Production Performance |
Minimal transient period followed by a long boundary-dominated flow period. Production rates are mainly related to permeability and declining reservoir pressure. From a traditional decline perspective, ‘b’ values typically range from 0 to 0.5 but can be somewhat larger if there are commingled layers. Recovery factor = 50% – 90% |
Very long transient (linear) flow period that can extend many years. In some cases, it is debatable if boundary-dominated flow will ever be fully realized. Production rates are mainly related to the success of creating a large fracture network around a long horizontal wellbore and to the matrix permeability. From a traditional decline perspective, ‘b’ values initially equal 2.0 (indicating linear flow), and then transition to <1.0 as boundary-dominated flow becomes prevalent. Recovery factor = 5% – 20% |
A comparison of conventional vs. tight/shale gas rate-time plots is shown below:
Shale Adsorption Correction
The shale adsorption correction is based on SPE paper 131772 by R. Ambrose et al.
In shale reservoirs that display adsorption, it has traditionally been thought that the adsorbed gas resides in the matrix of the organic-rich portion of the rock, which does not contribute significantly to the porosity of the rock. New evidence (Ambrose et al.) is suggesting that this organic-rich portion of the shale does, in fact, contain the majority of the porosity, as well as the adsorbed gas. If this assumption is valid, it is possible that the adsorbed gas is occupying a portion of the measured porosity, thus decreasing the pore space available to free gas. The volume occupied by the adsorbed gas must therefore be subtracted from the measured pore volume when calculating the volume available for free gas storage. The Ads. Sat. (Adsorption Saturation) option in Harmony performs this correction.
When gas is adsorbed, the molecules are so close to each other that they exist in a near liquid phase. The adsorption volume correction in Harmony converts the measured adsorbed gas content into a liquid, and subsequently treats this phase as an extra saturation term, similar to water saturation.
When the correction is enabled, an Adsorbed Saturation, “Sads” is calculated, and gas saturation is calculated as follows:
To calculate adsorbed saturation, define:
The adsorbed saturation can therefore be calculated as follows:
From SPE 131772, is defined as:
where
= Fractional volume occupied by adsorbed phase
M = Molecular weight of the gas (lbm/lbmole)
VL = Langmuir Volume (scf/ton)
pL = Langmuir Pressure (psi(a))
pi = Initial reservoir pressure (psi(a))
= Bulk density of the rock (g/cm3)
= Density of the adsorbed phase (g/cm3)
Note:
is the density of the adsorbed “liquid” phase. The default in Harmony for this value corresponds to the liquid methane density of 0.42 g/cm3.
The adsorption saturation correction is available for Gas Volumetric and Gas Analytical models. In Gas Volumetric analyses, the correction is available under Adsorption Properties. To enable the correction, click the Enable checkbox and then the Ads. Sat. checkbox.
In Gas Analytical models, click the Adsorption checkbox, then click the Ads. Sat. checkbox.
Enabling the Ads. Sat. correction will result in an adsorbed saturation (Sads) being calculated, which will in turn be subtracted from the original gas saturation (Sg). This corrected gas saturation will now be used in the volumetric calculations to calculate a new free-gas OGIP.
Note: The adsorption saturation correction will also affect the desorption compressibility (cd) term in the total compressibility used in the advanced material balance equation when calculating flow rates and pressures from the analytical model.