03 Geochemistry
03 Geochemistry
UNCONVENTIONAL RESERVOIRS
Hamid Emami-Meybodi
June 2023
Outline 2
Organic Matter
• Vitrinite Reflectance
Organic Content
• TOC determination
• Rock-Eval Pyrolysis
• Geochemical calculations
Petroleum Geochemistry 3
Why Do We Care?
• Soluble fraction contains free HC, whereas insoluble fraction can be converted to HC.
Hydrocarbon Source
Total rock
Organic matter Minerals
Organic matter
Other
(Pyrobitumen)
Bitumen
Kerogen
Kerogen Type
The van Krevelen diagram is commonly used
• Type I (H/C>1.2; O/C~0.5) to assess origin and maturity of kerogen.
hydrogen rich derived from zooplankton,
phytoplankton, microorganisms (bacteria) and
the lipid rich components of higher plants.
Type I
• Type II (1.2<H/C<1.4; O/C~0.1)
Increasing
intermediate composition, derived from Type II burial
mixtures of degraded/partly oxidized remnants
of higher plants and marine phytoplankton.
Type III
• Type III (0.7<H/C<1.0; O/C~0.2)
hydrogen poor but oxygen rich, derived from
Type IV
lignin-cellulose from higher plants.
Temperature (oC)
2) Catagenesis: progressive diagenesis in response
to thermal exposure (time and temperature).
Occurs at 40-150°C. Kerogen increases in carbon
content and decreases in volatiles (H and O).
Significant amounts of thermogenic CH4 and CO2.
Overall reaction example: 𝑅𝑅𝐶𝐶𝐶𝐶2 𝐶𝐶𝐶𝐶2 𝐶𝐶𝐶𝐶3 + 4𝐻𝐻2 𝑂𝑂 → 𝑅𝑅 + 2𝐶𝐶𝐶𝐶2 + 2𝐶𝐶𝐶𝐶4 + 5𝐻𝐻2
Maturation Process and Thermal Maturity 9
Vitrinite Reflectance
• Vitrinite reflectance index is the most commonly used method for determination of
coal rank and shale thermal maturity.
Vitrinite Reflectance
• Vitrinite reflectance (%Ro) is measured as the percentage of light reflected from the
sample, calibrated against a material which has ~ 100% reflectance (e.g., a mirror):
• Limitations
o is vitrinite present in the rock (for shales)?
o is vitrinite correctly identified (for shales)?
o is the sample polished properly?
o is the photomultiplier correctly calibrated?
o are any extraneous light sources affecting the reading by the photomultiplier?
Maturation Process and Thermal Maturity 11
• Proximate analysis
o Determination of the overall composition (i.e., moisture, volatile matter, ash, and
fixed carbon content).
• Ultimate analysis
o Absolute measurement of the elemental composition (i.e., carbon, hydrogen,
sulfur, nitrogen, and oxygen content).
Maturation Process and Thermal Maturity 13
• Proximate analysis provides relative abundance of moisture, ash volatiles and fixed
carbon content in coal (ASTM procedure).
As received
Surface moisture
Inherent moisture
Air dried
Ambient conditions
Dry
Volatile matter
105 -110oC Solid fuel
(CH4, H2, CO, CO2, N2) Dry, (C, H, O, S, N)
volatile
free
930 - 970oC Ash
Inert atmosphere (Nitrogen)
Combusted at 750oC
Oxygen atmosphere
wood 50 > 65 - -
sub-bituminous 80 40 33.5 5
semi-anthracite 92 8 36 1
anthracite 95 2 35.2 2
Maturation Process and Thermal Maturity 16
o Hydrocarbon migration
Gas
o OM type prone
Biodegeredation
Biogenic
Kerogen Bitumen Oil (Dry) Gas
Secondary
cracking
Source of gas
Thermogenic 1) Thermogenic (Kerogen
(Wet/Dry) Gas cracking, bitumen
Source rock cracking, oil cracking)
processes 2) Biogenic
Dead Carbon
(Modified from Jarvie et al., 2007)
Maturation Process and Thermal Maturity 17
• Used to evaluate the quality of source rock and to estimate adsorbed gas.
Question: In an exploration project area, the average annual temperature is 67°F, and
geothermal gradient is 1.9°F/100 ft. The project is in a structurally simple basin that, to date,
has experienced slow, continuous subsidence since basin initiation and only minor
deformation that resulted in simple fold structures. You have drilled a vertical exploration well
for an anticlinal trap in a sandstone reservoir identified from seismic data. In the exploration
well on the crest of the anticline, top of the reservoir sandstone was encountered at 21,113-ft
drilling depth, base was at 21,317-ft depth, and net sand thickness was 180 ft. Source rock
for the reservoir was a 300-ft thick shale that direct underlies the reservoir and has 3.5%
Type I kerogen. When drilling the well, you also encountered three coal beds having net
thickness of 45 ft between the depths of 2,230 and 2,297 ft.
• Knowing kerogen type, restoration of TOC at full conversion (only TOC remaining is
the spent organic carbon that has no remaining potential to generate petroleum):
Question: What is TOCpd for a fully converted sample with TOCo = 10 wt% for
a) Type I:
b) Type II:
c) Type III:
Maturation Process and Thermal Maturity 21
TOC (Original)
• When kerogen matures and produces gas and oil, the amount of kerogen will
decrease, whereas TOC remains constant until the generated hydrocarbons
expelled to other reservoirs.
TOC (Original)
6
Hybrid shale containing organic-rich and
Matrix & Organic Porosity (%)
5 organic-lean intervals
1
Organic-rich shale
0
0 1 2 3 4 5 6 7 8
TOC (wt%)
2.5
Kerogen is 44X more adsorptive than quartz
2.0
Illite 0.160 0
Younger shale
Kerogen 1.293 0 2 4 6 8 10 12 14
TOC (wt%)
(Data from Schettler and Parmely, 1991)
(Modified from Ross and Bustin, 2007)
Maturation Process and Thermal Maturity 26
3000
Pressure (psia)
2000
1000
0
0 20 40 60 80 140
Temperature (oC) (Modified from Bustin et al., 2009)
Maturation Process and Thermal Maturity 28
1.6
1.2
RO max (%)
0.8
0.4
0
0 4 8 12 16 20 24
TOC (wt%) (Modified from Bustin et al., 2009)
Organic Content 29
TOC Determination
1) Well log interpretations: Passey methodology
o Using sonic and resistivity logging information.
o Sonic log measures of interval travel time in formation.
o Resistivity log measures resistivity of formation.
Sonic OM increases the apparent transit time of acoustic logs (lower velocity).
TOC Determination
1) Well log interpretations: Passey methodology
TOC Determination
1) Well log interpretations: Passey methodology
𝑅𝑅
Δ𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 = 𝑙𝑙𝑙𝑙𝑙𝑙10 + 0.02 ∆𝑡𝑡 − ∆𝑡𝑡𝑏𝑏𝑏𝑏𝑏𝑏𝑏𝑏
𝑅𝑅𝑏𝑏𝑏𝑏𝑏𝑏𝑏𝑏
o The LOM could be calculated from vitrinite reflectance measurement (but it may not be
known in exploration wells!).
TOC Determination
Free hydrocarbons
(mg HC/g rock)
Generated hydrocarbon
assuming HC contain (mg HC/g rock)
83.33 wt% carbon
Rock-Eval Pyrolysis
• Based on the pyrolysis of the organic matter, four basic parameters are determined:
Temperature
generated by thermal cracking of kerogen.
Rock-Eval Pyrolysis
𝑆𝑆2
• Hydrogen index (mg HC/g C): related to H/C ratio 𝐻𝐻𝐻𝐻 =
𝑤𝑤𝑇𝑇𝑇𝑇𝑇𝑇
and origin/maturity of kerogen.
𝑆𝑆3
• Oxygen index (mg CO2/g C): related to O/C ratio 𝑂𝑂𝑂𝑂 =
𝑤𝑤𝑇𝑇𝑇𝑇𝑇𝑇
and is higher for land plants/inert OM than for
marine/algal OM.
𝑆𝑆1
• Productivity Index: ratio of already generated 𝑃𝑃𝑃𝑃 =
𝑆𝑆1 + 𝑆𝑆2
hydrocarbon to potential hydrocarbon. Also called
transformation ratio.
Organic Content 35
12.0
Rock-Eval Pyrolysis
Type I Type II
10.0
• Kerogen typing may be based on
S2 (mgHC/gRock)
o Visual examination of particulate matter. 8.0
Type III
o Its chemical composition: 6.0
HI vs. OI plots
4.0
S2 vs. TOC plots
2.0 Type IV
HI vs. Tmax
0 0 1.0 2.0 3.0 4.0 5.0 6.0
TOC (w%)
900 900
Type I Type I
800 Oil prone 800
Oil
Condensate
700 700
Type II
HI (mgHC/gTOC)
HI (mgHC/gTOC)
Type II
600 600
Oil/Gas prone
500 500
Gas
400 400 Type II & III
300 300
Type III Gas prone
200 200 Type III
100 Type IV No potential 100
Type IV
0 0
0 50 100 150 200 350 400 450 500 550
OI (mgCO2/gTOC) Tmax (oC)
Organic Content 36
Ro (%) Tmax (oC) TR (%) HIpd (mg HC/g TOC) Dry gas (%) C20+ (%)
Minimum 1.00 455 80 76 – 100 80 5
Best 1.40 475 95 <50 95 1
Highly converted
shale, potentially high-
Productive well
flow-rate gas system
Non-productive well
Low-level conversion,
low-flow-rate gas system
Geochemical Log
• Oil crossover effect: Oil saturation index, OSI = S1 x 100/TOC > 100 mg oil/g TOC.
(Jarvie, 2012)
Organic Content 38
(Jarvie, 2014)
Organic Content 39
Original HI
There are multiple ways to derive the original TOC value in shales from Rock-Eval
data. One way is from HIo. Two ways to estimate HIo:
Geochemical Calculations
For any given HIo, the percentage of generative organic carbon is calculated by:
Geochemical Calculations
• Estimating Original TOC from
o HIo: estimated from immature rocks or correlations.
o TOCpd: obtained from Rock-Eval, or Passey method, or etc.
o S1pd (present-day retained free HC): obtained from Rock-Eval.
o S2pd (HC equivalent of present-day GOC): obtained from Rock-Eval.
TOCo
Original
GOCo NGOCo
NGOCpd
Present-day
TOCpd
Geochemical Calculations
%𝐺𝐺𝐺𝐺𝐺𝐺𝐺𝐺𝑜𝑜 = 0.085 𝐻𝐻𝐻𝐻𝑜𝑜
1) Original GOC in TOCo
%𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑜𝑜
5) Original TOC %𝑇𝑇𝑇𝑇𝑇𝑇𝑜𝑜 = 100 ×
100 − %𝐺𝐺𝑇𝑇𝑂𝑂𝑂𝑂𝑜𝑜
%𝐺𝐺𝑂𝑂𝑂𝑂𝑜𝑜
7) Original hydrocarbon (mg HC/g rock) 𝑆𝑆𝑆𝑜𝑜 =
0.085
Organic Content 43
Geochemical Calculations
𝜌𝜌𝑠𝑠𝑠
𝜙𝜙𝐺𝐺𝐺𝐺𝐺𝐺 = %𝐺𝐺𝐺𝐺𝐺𝐺𝑜𝑜 − %𝐺𝐺𝐺𝐺𝐺𝐺𝑝𝑝𝑝𝑝 − %𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑎𝑎𝑎𝑎𝑎𝑎
𝜌𝜌𝑘𝑘𝑘𝑘
𝜌𝜌𝑠𝑠𝑠
𝜙𝜙𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 = %𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑎𝑎𝑎𝑎𝑎𝑎
𝜌𝜌𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑
9) The bbl equivalent of generated HC per one acre-ft of the source rock:
Geochemical Calculations
Question: The following data are available for a shale gas formation:
o HIo = 400 mg HC/g TOC o S1pd = 0.86 mg HC/g rock
o %TOCpd = 9.21 wt% o S2pd = 2.94 mg HC/g rock
o 𝜌𝜌𝑠𝑠𝑠 = 2.38 g/cc o 𝜌𝜌𝑘𝑘𝑘𝑘𝑘𝑘 = 1.18 g/cc
o 𝜌𝜌𝐻𝐻𝐻𝐻 = 0.85 g/cc o 𝜌𝜌𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 = 1.35 g/cc
Calculate original TOC, organic porosity, and the bbl equivalent of generated HC per
one acre-ft of the source rock.