Basin Modeling: The Key For Unconventional Shale Gas Assessment
Basin Modeling: The Key For Unconventional Shale Gas Assessment
assessment
J.M. Laigle*, BeicipFranlab, 232, Av. Napolon Bonaparte, 92502 Rueil-Malmaison, France
[email protected]
and
F. Schneider, L. Kuhfuss Monval, P. Lemouzy, BeicipFranlab, 232, Av. Napolon Bonaparte, 92502 Rueil-
Malmaison, France
Summary
The objective of this paper is to demonstrate how a basin modeling study can help in the evaluation of
unconventional shale gas / oil resources. We present an improved expulsion model that accounts for
the evolution of retention capacity as a function of maturity. Basin modeling is indeed key tools to help
in 1) the evaluation of the initial organic matter distribution, type and quality inside a formation, and 2)
the assessment of the present day total organic content (TOC) and maturity level of the formation.
These two parameters have a first order control on the volume of hydrocarbons generated and still
retained in the formation.
Introduction
Unconventional shale gas / shale oil resources refer to the gas/oil retained in source rocks at present
day.
Hydrocarbons are generated from kerogen during the source rock burial thanks to an elevation of the
temperature. As for the gas, it can be either of biogenic origin (generated at low maturity), or of
thermogenic origin (generated from the cracking of kerogen and oil at higher maturity level).
Part of the gas / oil generated is expelled from the source rock and then migrates to conventional
reservoir through carrier beds. The other part is retained by the source rock. Hydrocarbons composition
and volumes generated and retained depend on the initial organic matter properties and maturity.
Expulsion model in the past have been designed in order to quantify the expelled amounts. These
models have to be revisited and improved in order to better quantify the retained gas/oil.
The source rock retention capacity is tightly linked to the TOC content, the organic matter type and the
maturity. Hydrocarbons can be stored in source rocks by adsorption (this process concerns mainly the
gas) and as free hydrocarbons inside the porosity.
The porosity in mature source rock can be of two types: 1/ an effective porosity inside the non-organic
fraction of the rock. The evolution of this porosity vs. depth is controlled by burial and possible
diagenetic phenomena and 2/ porosity inside the organic matter itself. This second porosity is created
consequently to the loss of kerogen mass resulting of the transformation of organic matter into
hydrocarbons. It is assumed that the organic porosity is poorly connected to the non-organic effective
porosity as the pore threshold size in this organic porosity is believed to be of the nanometer order of
Retention Model
The adsorption represents also a a very efficient way to retain important volume of hydrocarbons inside
the source rocks, especially gas, as the density of the adsorbed layer is much higher than the density of
the free hydrocarbons phase in bottom conditions. The adsorption capacity of the rock is modeled as a
function of TOC, pressure and temperature through Langmuir isotherms that are calibrated from lab
experiments.
Workflow
To assess the potential for shale gas/oil in a formation, and as the volume and quality of hydrocarbons
generated depends on the initial geochemical parameters of the formation and the maturity evolution
through time, it is so necessary to evaluate:
The initial TOC, the organic matter distribution, type and quality;
The source rock present day maturity;
The storage capacity in the non-organic porosity, the organic porosity and the adsorption;
The generated HC quantity and quality, in order to quantify how much is retained in the source
rock intervals, in the porosity or by adsorption.
The initial TOC and kerogen initial HI are inverted from present day values. The present day TOC is
deduced from logs using an improved Carbolog methodology, which combines sonic and resistivity
logs (Carpentier et al, 1991; Schneider et al, 2010). For thick formations, or if the source rock is not well
known, a forward stratigraphic model, such as Dionisos, can help in assessing the source rock extent
and its internal sedimentary architecture.
The geochemical study, and especially the analysis of the Rock Eval data, enables to define the
kerogen type and quality and helps in selecting a relevant kinetic scheme. The basin model, calibrated
in maturity, allows estimating the transformation ratio of the kerogen, which is used to calculate the
initial TOC from the present day TOC.
Accurate compositional kinetic parameters for kerogen thermal decomposition are key elements in this
kind of study to predict the hydrocarbons fluids quality and type as a function of maturity. IFP kinetics
schemes have proven to be effective for representing accurately kerogen transformation into oil and
gas (Behar et al., 1997, 2008). They provide accurate and proven reactions scheme, able to predict
early or late gas. They have been successfully applied to Barnett and Posedonia shales (Behar &
Jarvie 2011). A solution for modeling biogenic gas has also been proposed (Schneider et al, 2010).
The modeling of hydrocarbons generation allows estimating the mass loss in the kerogen and the
organic porosity consequently created as a function of maturity.
The effective porosity in the shale matrix is estimated from log analysis, using a methodology that
allows correcting the log response from the effect of organic matter.
The integration of these tasks allows estimating the retention capacity of the formation, accounting for
shale effective porosity, organic porosity and adsorption capacity. This retention capacity depends
primarily of the TOC distribution and the TOC decreases with maturity.
The Gas Initially In Place (GIIP) can then be directly deduced from the basin modeling simulation
results, as a function of this retention capacity. The numerical model computes a GIIP value at play
scale and gives also access to lateral and vertical variations of gas distribution, depending on the
resolution of the source rock architecture implemented in the model.
Using experimental design and response surface-based methodology, the calculated GIIP may be
risked as a function of the uncertainty associated to the values of some key input parameters. GIIP
percentiles are then computed (P10, P50 and P90).
Examples
The methodology has been successfully applied in a number of formations over the world, including
Silurian in Algeria, Abu Gabra formation in Sudan, Los Monos formation in Argentina, or Lower Barnett
unit in the Fort Worth Basin (Romero-Sarmiento et al., 2012). Most of the formations studied present
significant differences with US plays, in terms of thickness, organic matter type, present day TOC,
adsorbed gas proportion, volume of gas per ton of rock.
Conclusions
The methodology presented allows estimating the hydrocarbons retention capacity of the formation and
the composition of hydrocarbons still retained. Our experience proves that these properties can vary
significantly laterally and vertically, and that the shale gas potential of a formation can then not be
simply deduced from a local evaluation that would be extrapolated throughout the basin.
Shale development presents huge opportunity and faces large geological uncertainty. Innovative tools
and methods are key to diminish the geological risk from resource evaluation to optimal well
development. Basin modeling can play a key role in the shale evaluation process as it allows
integrating all available knowledge into a consistent and comprehensive model of the studied area,
from early stages of exploration to pilot zone identification
Large research effort is still required to understand and to quantify the physical processes that occur in
shales and organic matter, such as: retention, adsorption, diffusion.
Acknowledgements
References
Behar, F. and D. Jarvie, 2011, Compositional modeling of gas generation from two shale gas resource systems: Mississippian
Barnett Shale (USA) and Lower Jurassic Posidonia Shale (Germany), Hedberg Research Conf., submitted
Behar, F., M. Vandenbroucke, Y. Tang, F. Marquis, and J. Espitalie, 1997, Experimental cracking of kerogen in open and
closed systems: determination of kinetic parameters and stoichiometric coefficients for oil and gas generation: Organic
Geochemistry, v. 26, p. 321339
Bhar F., F. Lorant, and L. Mazeas, 2008, Elaboration of a new compositional kinetic scheme for oil cracking: Organic
Geochemistry, v. 39, p. 764-782.
Carpentier, B., A. Huc, and G. Bessereau, 1991, Wireline 550 logging and source rocks. Estimation of organic carbon by the
CARBOLOG method, The Log Analyst, MayJune, 1991, p. 279297.
Loucks R., R. M. Reed, S. C. Ruppel, and D. M. Jarvie, 2009, Morphology, Genesis, and Distribution of Nanometer-Scale
Pores in Siliceous Mudstones of the Mississippian Barnett Shale: Journal of Sedimentary Research, v. 79, p. 848-861.
Romero-Sarmiento M.F., M. Ducros, B. Carpentier, F. Lorant, S. Rohais, I. Moretti, and A.Y. Huc, 2012, Predicting TOC,
Organic Porosity and Gas Retention Distribution in a Shale-Gas System Using Petroleum Basin Modelling: Search and
Discovery Article #40973
Schneider, F. J. S., J. A. Noya, and C. Magnier, 2010, Model of low-maturity generation of hydrocarbons applied to the
Carupano Basin, Offshore Venezuela, in K. Peters, D. Curry, and M. Kacewicz, eds., Basin and Petroleum System Modeling:
2009 Napa Hedberg Conference: Hedberg Series No. 4, p. 1 19.