Chemical, Molecular, and Microstructural Evolution of Kerogen during

Mar 21, 2018 - Total organic sulfur and sulfur-XANES data show that sulfur in Woodford kerogens is dominated by aromatic sulfur (thiophene) and that r...
0 downloads 6 Views 2MB Size
Subscriber access provided by TUFTS UNIV

Article

Chemical, Molecular, and Microstructural Evolution of Kerogen during Thermal Maturation: Case study from the Woodford Shale of Oklahoma Paul R. Craddock, Kyle D Bake, and Andrew E Pomerantz Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.8b00189 • Publication Date (Web): 21 Mar 2018 Downloaded from http://pubs.acs.org on March 21, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

1

Chemical, Molecular, and Microstructural Evolution of Kerogen during

2

Thermal Maturation: Case study from the Woodford Shale of Oklahoma

3

Paul R. Craddock*, Kyle D. Bake, Andrew E. Pomerantz

4

Schlumberger-Doll Research Center, Cambridge, MA, USA

5

* Corresponding author: [email protected], +1-617-768-2042

6 7

Manuscript submitted to: Energy & Fuels

8

Original submission: January 15, 2018

9

Revised submission: March 20, 2018

10 11

Abstract

12

Integrated elemental, spectroscopic (infrared spectroscopy, X-ray absorption near edge

13

structure), and gas intrusion (helium pycnometry, nitrogen adsorption) analyses are used to

14

characterize the bulk chemical, molecular, and physical microstructures of kerogen spanning a

15

thermal maturity transect (vitrinite reflectance, Ro, from 0.5% to 2.6%) across the Woodford

16

Shale of the Anadarko Basin, Oklahoma. The integration takes advantage of novel procedures to

17

prepare kerogen isolates that preserve both the chemical and physical properties of the organic

18

matter in the bulk shale. The Woodford kerogens follow the expected trends in H/C and O/C

19

coordinates during thermal maturation for type II kerogen. Infrared spectra show that loss of

20

hydrogen from kerogen is related to cracking of hydrogen-rich aliphatic (alkyl) carbon structures

21

from aromatic carbons. Within the range of Ro values < 1.5%, peripheral aromatic carbons

22

remain highly substituted with alkyl (methyl, methylene) and probably heteroatom functional

23

groups. At Ro values > 1.5%, these substitutions are substantially removed and replaced by

24

hydrogen. The evolution of carbon structures inferred from the IR spectra is supported by known 1 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

25

carbon bond dissociation energies for carbonaceous materials. Total organic sulfur and sulfur-

26

XANES data show that sulfur in Woodford kerogens is dominated by aromatic sulfur

27

(thiophene) and that reactive, aliphatic sulfur (sulfide) is eliminated at low degrees of thermal

28

stress (Ro ≤ 0.9%). At higher thermal stress, sulfur speciation is stable and dominated by

29

thermally stable thiophene (sulfur in aromatic rings). The physical properties of Woodford

30

kerogens evolve during thermal maturation in a manner consistent with their molecular

31

characteristics. Skeletal density of kerogen increases during maturation in a manner that is

32

linearly correlated to its atomic H/C ratio and inferred aromatic carbon content. Specific surface

33

area of kerogen also increases during maturation, reflecting development of internal pores within

34

the kerogen skeletal framework as aliphatic carbon structures are preferentially cracked and

35

expelled from solid kerogen. The quantitative chemical and structural changes expressed by the

36

Woodford kerogens during thermal maturation, including their hydrogen and carbon content,

37

carbon speciation, and skeletal density, are shown by comparison to not be measurably different

38

for other type II kerogens from numerous oil- and pas-producing shale plays, indicating that

39

thermal stress acts to drive maturation of type II kerogen in a similar way globally.

40

2 ACS Paragon Plus Environment

Page 2 of 49

Page 3 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

41

Energy & Fuels

1. Introduction

42

Kerogen is the most abundant form of sedimentary organic matter (OM)1 and a defining

43

component of unconventional petroleum source rocks, commonly termed shale. From a

44

geochemical perspective, kerogen is the solid, insoluble, and nonvolatile sedimentary OM. Its

45

molecular structure (comprised essentially of C, H, N, S, and O) controls the type, amount, and

46

quality of hydrocarbons generated during thermal maturation.2 In a geological context, kerogen

47

is commonly studied using petrographic inspection (e.g., vitrinite reflectance) and thermal

48

decomposition methods (e.g., elemental analysis, programmed pyrolysis). Vitrinite reflectance3-6

49

and programmed pyrolysis7-9 (e.g., Rock-Eval) provide a determination of the thermal maturity

50

of kerogen from which structural characteristics may generally be inferred, but they give no

51

direct insights on structural properties. From a chemical perspective, elemental analysis10 offers

52

bulk chemical concentrations, but also provides no direct molecular structural information.

53

Several solid-state spectroscopic methods are available for the analysis of specific

54

molecular structures in carbonaceous materials such as kerogen, including infrared (IR)

55

spectroscopy, Raman spectroscopy,

56

near-edge structure (XANES), and X-ray photoelectron spectroscopy (XPS). Raman

57

spectroscopy11-18 and

58

environment of carbon comprising aliphatic and aromatic moieties, as well as their absolute

59

abundance. Advanced

60

discrimination and quantification of specific C, CH, CH2, and CH3 functional groups including

61

their connectivity or proximity, from which parameters such as the average carbon number of

62

aromatic clusters and of aliphatic chains, can be estimated.22,

63

routinely identify functions associated with C, H, and O, including aromatic C=C and CH,

13

C nuclear magnetic resonance (NMR), X-ray absorption

13

C NMR19-27 are two techniques capable of probing the bonding

13

C NMR spectral acquisition and editing techniques now offer

3 ACS Paragon Plus Environment

28, 29

Infrared spectroscopy can

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 49

64

aliphatic CH, CH2, and CH3, and oxygenated C–O and C=O groups.26,

65

analysis of organic sulfur and nitrogen species is now established using sulfur- and nitrogen-

66

XANES.25, 48-67 Solid-state XPS is capable of quantifying several functionalities in carbonaceous

67

materials associated with carbon, oxygen, sulfur, and nitrogen.51,

68

several of these techniques and their applications in the analysis of kerogen and other complex

69

carbonaceous materials is given by Pomerantz.73 Separately, methods such as scanning electron

70

microscopy,74-79 gas adsorption,80-84 and small angle neutron and X-ray scattering,84-88 exist for

71

analysis of aspects of the physical microstructure of kerogen, such as its pore size distribution

72

and surface area, performed on either bulk shales or kerogen isolated from the rock.

73

Measurements run on bulk shales can be challenged by presence of volumetrically-dominant

74

inorganic matrix that obscures the signal from the organic matter. Kerogen isolates have

75

generally been prepared using acid demineralization and drying procedures that preserve its

76

chemical structure but alter its physical microstructure, such that microstructural properties of

77

kerogen isolates are not representative of the organic matter in situ.87 More often than not,

78

methods for the analysis of kerogen are used in isolation, providing a far-from-complete picture

79

as to the molecular makeup of these complex materials.

55, 67-72

27, 30-47

Quantitative

A recent review of

80

Kelemen et al.25 carried out one of the most comprehensive surveys of kerogen

81

compositions of different types and from several organic-rich sedimentary facies. These kerogen

82

types [see

83

algal OM in freshwater environments), type II (initial intermediate H/C and O/C ratio, derived

84

predominantly from planktonic OM in marine environments) and type III (initial low H/C and

85

high O/C ratio, derived from terrestrial higher-order plants). Their results from 13C NMR, XPS,

89, 90, 91

] included type I (initial high H/C and low O/C ratios, typical of aliphatic-rich

4 ACS Paragon Plus Environment

Page 5 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

86

and XANES provided a detailed comparison of kerogen chemical structure of different types and

87

organo-facies, and to a limited extent as a function of thermal maturity.

88

Here, we present an integrated analysis of kerogen, with a focus on a single organic-rich

89

sedimentary formation (predominantly type II kerogen of the Woodford Shale, Anadarko Basin,

90

Oklahoma) that presents a natural thermal maturity transect from immature to post-mature

91

(vitrinite reflectance, Ro, range from at least 0.5 to 4%). Samples from a single formation were

92

chosen to minimize differences in original kerogen composition as is encountered when

93

comparing kerogens of different types (i.e., types I, II, and III) and/or from different organo-

94

facies. In this way, measured differences in kerogen characteristics among the samples are

95

predominantly, if not entirely, a result of their level of thermal maturation. The study focuses on

96

type II kerogen because it represents by far the major fraction of sedimentary OM in all paleo-

97

marine petroleum source rocks of North America and elsewhere that are producing substantial

98

and economic quantities of unconventional oil and gas. As such, the results here may be used

99

more broadly as an analog to infer or contrast kerogen characteristics in other shale plays.

100

Using novel procedures for the preparation and analysis of isolated kerogen, here we

101

assess both its chemical and physical microstructural characteristics. The results reveal detailed

102

and complementary trends in chemical composition and physical microstructure of type II

103

kerogen as a function of thermal maturity. This type of information is critical, for example, to

104

optimize the properties of kerogen needed for petrophysical interpretations of shale, as well as to

105

constrain molecular models of kerogen structure92-97 that can provide insights to hydrocarbon

106

storage and transport.

107

2. Materials and Methods

5 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

108

Page 6 of 49

2.1. Nomenclature

109

Kerogen is defined here as all solid, insoluble, and non-volatile fractions of sedimentary

110

organic matter (OM), following established geochemical conventions.1 This convention arises

111

from the methods by which bulk kerogen is chemically isolated from forms of sedimentary OM

112

(e.g., source-rock bitumen) that are soluble in usual organic solvents and from the inorganic

113

mineral matrix of the bulk rock, which makes kerogen amenable to study by numerous analytical

114

techniques. This definition is not precisely the same as the petrographic definition of kerogen as

115

it is studied by optical inspection, in which this term is restricted to “primary” [e.g.,

116

“depositional” OM [e.g.,

117

“secondary”—being formed from thermal maturation of previously soluble source-rock bitumen

118

or liquid oil98—but also insoluble. This secondary insoluble OM is commonly termed

119

“pyrobitumen” or “post-oil solid bitumen” [e.g., 6, 78, 79, 98], and it is distinguished optically from

120

primary insoluble OM on the basis of texture, such pore-filling characteristics.6 These two

121

insoluble fractions are inseparable using chemical techniques for the isolation of bulk insoluble

122

OM, and are together termed kerogen herein. New techniques are being developed to study the

123

chemical and mechanical properties of insoluble OM in situ without the need for its isolation,99

124

such that certain characteristics of primary and secondary insoluble OM may be investigated

125

separately.

126

2.2. Samples

76

78, 79

] or

]. There may exist a fraction of sedimentary OM that is

127

The Woodford Shale is a Late Devonian-Early Mississippian paleo-marine, organic-rich

128

sedimentary facies, containing predominantly kerogen derived from marine planktonic

129

organisms (type II kerogen) at concentrations up to 25 wt% total organic carbon (TOC) and is an

130

important hydrocarbon source rock in Oklahoma.100 The Woodford Shale is widely distributed 6 ACS Paragon Plus Environment

Page 7 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

131

across a geologic province called the Southern Oklahoma Aulacogen,101, 102 with recent drilling

132

and exploration concentrated in four plays: the Anadarko Basin in western Oklahoma, the

133

Arkoma Basin in eastern Oklahoma, the Ardmore Basin in southern Oklahoma, and the

134

Cherokee Platform in northwest Oklahoma.103 The tectonic history of the Southern Oklahoma

135

Aulacogen and general characteristics of the Woodford Shale have been detailed in earlier

136

publications [e.g., 102, 104, 105].

137

The samples studied in this paper are core and cuttings from wells drilled into the

138

Woodford Shale in the Anadarko Basin, which were curated and studied previously for their

139

vitrinite reflectance by the Organic Petrography Laboratory (OPL) of the Oklahoma Geological

140

Survey.106 The samples (Table 1) range from immature (vitrinite reflectance mean Ro = 0.55%)

141

to post-mature (Ro = 4.05%), which can be used to assess trends among kerogen properties as a

142

function of level of thermal maturation within a coherent sedimentary facies. The method used

143

for previous determination of the mean, random vitrinite reflectance values of these samples is

144

described elsewhere.104 For this study, between 15 and 90 g of material was received for each

145

sample, of which ~ 15 g was crushed to fine powder (~ 10 µm) using an auto-mortar. The

146

powders were cleaned of extractable organics (e.g., residual drilling fluid, source-rock bitumen,

147

oil) with dichloromethane in a Soxhlet extractor. After air-drying, the powders were split into

148

homogeneous fractions for subsequent whole-rock screening comprising TOC and programmed

149

pyrolysis (Rock-Eval). This screening was used to select a subset of samples for kerogen

150

isolation that both spanned a range of thermal maturity and contained sufficient mass of kerogen

151

required for a complete chemical and microstructural analysis. Kerogen isolates were obtained

152

from the bulk powders using open-system chemical treatment procedures.107 Briefly, a series of

153

concentrated HF and HCl acid additions were made to remove silicate and carbonate minerals,

7 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

154

followed by addition of acidic CrCl2 to remove metal sulfides (predominantly pyrite). Removal

155

of inorganic minerals was confirmed by high-temperature (~ 1000 °C) ash analysis, with ash

156

concentrations of 0.1–2 wt% relative.

157

The isolated kerogen fractions were suspended in ethanol and subsequently dried using

158

critical point drying (CPD) techniques to preserve their microstructure.108 The procedure

159

employs liquid-gas phase transition without crossing the phase boundary. Kerogen isolates were

160

flushed with liquid CO2 at 5.5 × 106 Pa (800 psi) and 0 °C for 2400 s (40 min) to remove

161

miscible ethanol, followed by vaporization and purging of supercritical CO2 above its critical

162

point (31.1°C, 1071 psi). The cycle was repeated using a 1200 s (20 min) flush. The dry kerogen

163

isolates were then ready for study. This procedure is demonstrated to preserve aspects of kerogen

164

microstructure, such as surface area, which collapse under traditional oven-drying procedures.108

165

2.3. Total organic carbon and Rock-Eval pyrolysis

166

Total organic carbon (TOC) measurements were run on the cleaned, bulk powders using

167

established combustion-coulometry methods.109 Because the samples were solvent-cleaned, the

168

reported TOC values represent only the organic carbon content from the insoluble organic

169

fraction, kerogen.

170

Programmed pyrolysis (Rock-Eval) is a common alternative to vitrinite reflectance for

171

determination of thermal maturity, and was run in this study to evaluate performance against the

172

vitrinite reflectance benchmark. The technique uses temperature-programmed pyrolysis of a

173

powdered sample in an inert atmosphere from which the amounts of free hydrocarbons (S1 peak,

174

units of mg of hydrocarbon per g of sample) and potentially generative hydrocarbons (S2 peak,

175

units of mg of hydrocarbon per g of sample) are quantified from the measured pyrograms. For

8 ACS Paragon Plus Environment

Page 8 of 49

Page 9 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

176

this set of samples, the S1 peak is negligible (< 1 mg/g) because they were solvent-cleaned of

177

existing hydrocarbons. Rock-Eval hydrogen index (H.I.) and Tmax are two metrics commonly

178

calculated from the measured pyrograms. H.I. is defined as equal to 100 × S2/TOC (units of mg

179

hydrocarbon per g of total organic carbon). Tmax is the temperature (°C) at which the maximum

180

amount of hydrocarbon is generated from kerogen decomposition (i.e., temperature of the S2

181

peak maximum), and can be calibrated specifically to equivalent vitrinite reflectance, for

182

example using the correlation: %VRe = 0.018 × Tmax − 7.16.110

183

2.4. Elemental Analysis

184

Kerogen bulk compositions were measured using routine combustion elemental

185

analysis.111 Carbon, hydrogen, nitrogen, and sulfur were quantified by flash combustion in a

186

stream of oxygen at ~ 1000 °C. Oxygen was quantified by flash pyrolysis in an inert atmosphere

187

at ~ 1000 °C.

188

2.5. Infrared spectroscopy

189

Infrared (IR) spectroscopy quantifies the distribution of vibrational frequencies of a

190

sample and has been used extensively to study the functional group abundances associated with

191

carbon, hydrogen, and oxygen in various carbonaceous materials including coals30-37, kerogens26,

192

27, 38-45, 112

, bitumens43, and asphaltenes45-47, 64.

193

Isolated kerogens were prepared for IR analysis using the potassium bromide (KBr) pellet

194

method, which enables the collected IR spectra to be interpreted quantitatively according to the

195

Beer-Lambert Law and thus for the IR spectra of different analytes to be compared directly.

196

Splits of kerogen isolates were gently hand-ground using a clean mortar and pestle, and a

197

precisely known mass of disaggregated kerogen powder (2.000 ± 0.010 mg) was added to a 59 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

198

mL tungsten carbide mixing vial containing two 7-mm diameter agate spheres. Precisely 900.00

199

± 0.50 mg of KBr powder (XL Spectrograde powder form, International Crystal Laboratoroes,

200

Garfield, NJ) was weighed out, and the KBr was added incrementally to the mixing vial, with the

201

mass of KBr added being approximately equal to the total mass of material in the mixing vial.

202

The mixture after each addition was mixed using a Retsch MM400 mixer mill (Verder Scientific

203

Inc., Newtown, PA) run at 23 Hz for 5 min. In total, six additions of KBr were made for a

204

mixing time of 30 min. Next, 200.00 ± 0.50 mg of the homogeneous powdered mixture was

205

placed into a prepared die and pressed into a 13-mm diameter pellet at 10 ton/cm2 for 10 min

206

under vacuum. The prepared kerogen pellets were scanned in transmission mode using a Bruker

207

Vertex70 dual-range Fourier transform infrared (FTIR) spectrometer. The measurement

208

comprised 25 individual scans of 2 s each across the mid-IR range (~ 5200–450 cm-1) with a

209

resolution of 2 cm-1. The collected spectrum was calculated as the mean of the 25 scans in

210

reported in absorbance units. A background (blank) absorption spectrum was scanned at the

211

beginning of each measurement and subtracted automatically from the kerogen absorption

212

spectrum.

213

Several IR structural indices were computed from the acquired spectra, based on band

214

assignments and conventions established by earlier studied.35, 36, 39, 42 The IR intensity of each

215

band of interest was computed as its integrated area following the curve deconvolution

216

procedures described in a previous publication.43 Table 2 defines the IR structural indices. The

217

CH3/CH2 ratio is inversely related to the average length of aliphatic chains and is positively

218

related to the degree of chain branching. Aromaticity is positively correlated to the relative

219

abundance of aromatic versus aliphatic carbons that are bound to one or more hydrogens. The A-

220

factor and C-factor represent, respectively, the amount of aliphatic carbon and oxygenated

10 ACS Paragon Plus Environment

Page 10 of 49

Page 11 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

221

carbon species relative to that of aromatic carbon. These indices were used to quantify changes

222

in the functional group abundances of kerogen during thermal maturation.

223

2.6. Sulfur K-edge X-ray absorption near-edge structure spectroscopy

224

Sulfur K-edge X-ray absorption near-edge structure spectroscopy (XANES) can quantify

225

the abundance of sulfur-containing functional groups in carbonaceous materials such as coal48-52,

226

kerogen25, 53-56, bitumen56-58, asphaltenes59-64, and asphalts.113 The techniques measures electronic

227

transitions from 1s orbitals to vacant molecular orbitals with 3p character. The energy of the

228

transition is a function of sulfur oxidation state (increasing energy with increasing oxidation

229

state), which relates to the bonding environment of the sulfur atom. Sulfur chemistry of kerogen

230

is shown to impact the kinetics of petroleum (oil and gas) generation in source rocks.114. Sulfur

231

content and chemistry can influence the surface polar functionality of carbonaceous materials,

232

including coal, kerogen, and bitumen, that may impact the storage (e.g., adsorption) and

233

transport of petroleum in the subsurface.56, 115, 116

234

Kerogen isolates were prepared for XANES by dilution of ~ 50 mg of analyte in a sulfur-

235

free matrix (boron nitride powder; > 99.5 % from Alfa Aesar, Ward Hill, MA) to ~ 1 wt% total

236

sulfur to prevent self-absorption effects.56 The powder mixture was then manually pressed into a

237

6-mm diameter pellet for analysis. XANES spectra were collected in fluorescence mode at the

238

Advanced Photon Source beamline 9-BM using a Stern-Heald-Lytle detector mounted in a

239

helium-purged sample chamber. The range of photon energies for spectral acquisition was 2450

240

to 2600 eV, and energy was calibrated against a Na-thiosulfate pre-edge feature at 2469.20 eV.

241

Measured XANES spectra were recorded as the ratio of the intensity of the total fluorescence

242

signal, I, to the intensity of the excitation radiation, I0, as a function of the excitation photon

243

energy. Spectra were baseline-corrected by setting I/I0 to zero in the pre-edge region (< 2465 eV) 11 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

244

and normalized by setting I/I0 to unity in the post-edge region (> 2490 eV). The relative

245

concentrations of organic sulfur species in the kerogens were quantified by fitting the resulting

246

spectra to a linear combination of model sulfur compounds (Figure 1). Inorganic sulfur forms

247

such as pyrite and sulfate (anhydrite, barite), if present in the isolate phase, can overlap the

248

XANES spectra of kerogen and so were included in the spectrum fitting, but their abundance was

249

not reported when presenting the molar fractions of organic sulfur species in kerogen. The

250

uncertainty on the analysis is approximately 10 mol%.

251

2.7. Helium pycnometry

252

Kerogen density is known to increase with level of thermal maturation112, 117-119 and is a

253

critical parameter for petrophysical determination of formation volumes and porosity. Helium

254

pycnometry was run to determine the absolute (“skeletal”) density of bulk kerogen isolates.

255

Measurements were made using a Micromeritics AccuPyc II 1340 gas pycnometer that

256

determines the volume of a precisely known mass of analyte based on Boyle’s Law.

257

Approximately 300–400 mg of sample was analyzed inside of a 1 cm3 volumetric chamber. Ten

258

volumetric determinations were made for each sample, and the volume computed as the mean of

259

the replicates. Density was calculated according to mass/volume relationship and is reproducible

260

for a given sample within ± 0.02 g/cm3. Densities were corrected for the presence of residual

261

inorganic minerals in the kerogen isolate assuming the ash was from sulfate (as indicated by

262

XANES sulfur species, here represented by anhydrite, ρ = 3 g/cm3), with the correction being

263

between 0.02 and 0.06 g/cm3.

264

2.8. Nitrogen adsorption surface area analysis

12 ACS Paragon Plus Environment

Page 12 of 49

Page 13 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

265

The preparation of kerogen isolates incorporating novel critical point drying (CPD)

266

methods enabled microstructural analysis of these samples. Specific surface area (SSA)

267

determinations were made on a Micromeritics ASAP 2420 surface area analyzer using multi-

268

stage nitrogen adsorption. The SSA of each sample was calculated from the Brunauer-Emmett-

269

Teller (BET) adsorption isotherm. Several hundred milligrams of sample was precisely weighed

270

into a surface area tube followed by a fill-rod to minimize head space. The pre-measurement

271

mass of the sample plus tube was recorded. The samples were then degassed for 6 h at

272

programmed temperature and pressure (60 °C, < 5 µm Hg) to eliminate surface-adsorbed

273

atmospheric gases. The mass of sample plus tube was recorded again to confirm adsorbed gas

274

loss. The surface area tube was wrapped in a thermal jacket, loaded onto the SSA analyzer inside

275

of a dewar, and immersed in liquid nitrogen. The jacket ensured thermal equilibrium inside of

276

the surface area tube during the measurement. The measurement consisted of sequential

277

quantification of adsorbed nitrogen as a function of relative pressure, P/P0, between 0 and 1

278

where P and P0 are the equilibrium and saturation pressure of the adsorbate (nitrogen). The

279

specific surface area was then computed from the linear BET adsorption isotherm plot between

280

P/P0 equal to 0.012 and 0.286.

281

2.9. Assessing preservation of kerogen properties during acid demineralization

282

Tests were performed to confirm that the acid demineralization procedure used to isolate

283

kerogen from the bulk sample did not oxidize or otherwise measurably alter its original

284

composition or structure. For this purpose, it was necessary to study as a proxy for kerogen, a

285

solid carbonaceous material (asphaltene) that could be obtained in a pure form without acid

286

demineralization, allowing the carbonaceous material to be directly analyzed and compared

287

before and after exposure to concentrated acids. Asphaltenes are macromolecular compounds 13 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

288

representing the heaviest fraction of petroleum and whose composition shares commonalities

289

with that of kerogen.120 Approximately 1 g of solid asphaltene was precipitated from a Middle

290

East black oil (UG8) by diluting in excess volume of n-heptane (40:1 ratio by mass of n-heptane

291

to oil) and stirring for ≥ 12 h. The asphaltene precipitate was filtered under vacuum over a 0.45

292

µm Teflon filter, and then transferred from the filter to a Soxhlet extractor and washed with n-

293

heptane for 24-48 h to eliminate any remaining soluble impurities. The precipitate was split into

294

two equal fractions. One fraction was analyzed by transmission FTIR spectroscopy and XANES

295

in its original form. The second fraction was exposed to the same concentrated acid procedure as

296

used to isolate kerogen from shale and then also analyzed by transmission FTIR spectroscopy

297

and XANES.

298

Figure 2 compares the FTIR and XANES spectra of the original and acid-treated

299

asphaltene splits. The IR spectra are practically identical with respect to the characteristic IR

300

absorption bands (e.g., C—C, C—H, and C—O) for carbonaceous materials. The IR spectra

301

show a 0.1 absorbance unit difference between 1150–1250 cm-1, but this region of the IR

302

spectrum is commonly associated with strongly absorbing inorganic bonds and does not impact

303

the computed IR structural indices. These bands are inferred to represent mineral impurities and

304

are not always observed in the kerogen IR spectra, described below. The XANES spectra and

305

sulfur species abundances are also the same within the reproducibility of the data. The apparent

306

small decrease, albeit within analytical error, in the abundance of oxidized sulfur species

307

(sulfoxide, sulfonate) in the acid-treated asphaltene is inconsistent with oxidation of

308

carbonaceous materials during treatment or storage. These data support the resistance to

309

measurable oxidation or other chemical alteration of carbonaceous materials prepared using the

14 ACS Paragon Plus Environment

Page 14 of 49

Page 15 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

310

same acid treatment as the kerogens analyzed here, in line with the similar observations made by

311

earlier studies of kerogen structure.43, 50, 56, 61

312 313

3. Results

314

3.1. Thermal maturity estimates

315

Figure 3 compares Rock-Eval parameters versus benchmark vitrinite reflectance for the

316

set of Woodford Shale samples. There exists the expected positive correlation between Rock-

317

Eval Tmax and vitrinite reflectance to values up to ~ 1.4% Ro (approximately late-oil window),

318

and within this range Ro-equivalence computed from Tmax [VRe = 0.018 × Tmax – 7.16]110 agrees

319

with that from vitrinite reflectance. Beyond this range (condensate and dry-gas windows), Tmax

320

estimates, and thermal maturity estimates from Tmax, are lower than expected. Consistent with

321

previous studies,8,

322

source rock samples with S2 signals less than ~ 1 mg HC/g TOC.

323

3.2. Elemental compositions

44

Tmax is found here to be an unreliable indicator of thermal maturity for

324

Table 3 reports the bulk elemental (C, H, N, S, O) concentrations of the kerogen isolates

325

from the Woodford Shale. It is known from the Van Krevelen diagram (H/C vs. O/C) that the

326

hydrogen and oxygen content of kerogen decreases during increasing thermal maturation,10, 90

327

Figures 4 and 5 plot, respectively, for the Woodford kerogens, their atomic H/C versus atomic

328

O/C and atomic H/C versus vitrinite reflectance, Ro. The data identify the kerogens as type II,

329

confirm their expected decrease in hydrogen content with increasing level of thermal maturation

330

(H/C ratio decreasing from 1.09 to 0.45), and indicate that they have low oxygen content, even at

331

or prior to the onset of oil generation (O/C ratio ranging from 0.04 to 0.08). The increase in the 15 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

332

oxygen content of the most thermally mature samples is unexpected, and likely reflects

333

contamination from inorganic oxygen from trace amounts of refractory oxides in the kerogen

334

isolates or from water in the pyrolysis measurement. For the Woodford kerogens, total sulfur

335

content sharply decreases from 3.1–3.5 wt% to 1.3–2.1 wt% at a vitrinite reflectance, Ro, value

336

of ~ 0.9%. Total nitrogen content is low (2.0–2.6 wt%) and shows no clear evolution during

337

thermal maturation.

338

3.3. Infrared structural indices

339

Figure 6 plots the FTIR spectra of the Woodford kerogens in order of increasing thermal

340

maturity, with their key structural groups indicated. The IR intensities of structural groups reveal

341

clear trends with increasing thermal maturity, comprising a decrease in the intensities of aliphatic

342

C—H bands and, to a lesser extent, in oxygenated (C—O, C=O) bands, and an increase in the

343

intensities of aromatic C=C and C—H bands. Table 3 reports group relative abundances

344

quantified by IR structural indices: CH3/CH2, Aromaticity, A-factor, and C-factor. Figure 7 plots

345

the computed IR structural indices of the kerogens as a function of their thermal maturity

346

(vitrinite reflectance, Ro) and atomic H/C ratio. The increase in ratio of methyl (—CH3) to

347

methylene (—CH2) groups indicates a shortening of the mean length of aliphatic chains in

348

kerogen during maturation. The A-factor, representing the relative abundance of aliphatic C—H

349

versus aromatic C=C groups, decreases during maturation. Aromaticity as defined from IR

350

spectroscopy reflects the relative amount of hydrogen bonded to aromatic carbons versus

351

aliphatic carbons, and increases during thermal maturation. The C-factor, representing the

352

relative abundance of C—O and C=O versus aromatic C=C groups, also decreases during

353

maturation. This FTIR data supports that the unexpected increase in bulk oxygen content of the

16 ACS Paragon Plus Environment

Page 16 of 49

Page 17 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

354

most mature kerogens reflects contamination in the bulk chemical analysis from inorganic

355

oxygen in mineral or water impurities.

356

3.4. Sulfur XANES

357

The sulfur functional group abundances of kerogen derived from the sulfur XANES

358

analysis are reported in Table 3. The values are expressed as mole percent of total sulfur. Small

359

amounts (2–9 mol%) of inorganic sulfate were identified in the XANES spectra, and excluded

360

when comparing the organic-sulfur distribution between samples. Figure 8 plots the molar

361

abundances of organic-sulfur species in kerogen as a function of their thermal maturity. Organic

362

sulfur in all kerogens is dominated by thiophene (sulfur bonded to aromatic carbon). The two

363

least mature kerogens studied (0.55–0.66% Ro, 3.1–3.5 wt% Stotal) contain 5–6 mol% of organic

364

sulfide (sulfur bonded to aliphatic carbon) that is not identified in the more mature samples (0.9–

365

2.6% Ro, 1.3–2.1 wt% Stotal).

366

3.5. Skeletal density

367

Table 3 reports the skeletal density of the Woodford kerogens derived from helium

368

pycnometry measurements. Figure 9a shows that the skeletal density of kerogen increases

369

systematically with increasing level of thermal maturation from ~ 1.2 g/cm3 (0.55% Ro) to 1.55

370

g/cm3 (2.6% Ro), reflecting an increase of 30% relative. Fig 9b shows that skeletal density plots

371

along an array that is nearly inverse linear to their atomic H/C ratio.

372

3.6. Specific surface area

373

The Woodford kerogen isolates were prepared using novel acid demineralization plus

374

critical point drying procedures that retain the microstructural properties of kerogen.108 Table 3

375

reports the resulting Brunauer-Emmett-Teller (BET) calculated specific surface area (SSA) of 17 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

376

the Woodford kerogens, which range over an order of magnitude from ~ 50 to ~ 300 m2/g.

377

Figure 10 plots SSA as a function of thermal maturity, showing the robust, positive correlation

378

between the two.

379

4. Discussion

380

The compositional evolution of kerogen from the paleo-marine Woodford Shale of the

381

Anadarko Basin, Oklahoma, is assessed here as a function of thermal maturation using a

382

combination of solid-state analytical techniques. The focus on sedimentary organic matter (OM)

383

from a single sedimentary formation aids to minimize compositional heterogeneity among

384

kerogens resulting from processes other than thermal maturation (e.g., differences in kerogen

385

type and depositional environment). Our data show that kerogen in the Woodford Shale evolves

386

dramatically and in a systematic manner during thermal maturation with respect to both its

387

chemical (e.g., carbon, hydrogen, and sulfur abundance and structure) and microstructural

388

properties (e.g., absolute density, specific surface area). The compositional characteristics of

389

kerogen in the Woodford Shale are compared, where available, with those measured for

390

kerogens from other paleo-marine sedimentary formations (type II kerogen), which provides

391

insights to whether thermal maturation of kerogen in other major oil- and gas-producing shales

392

follows compositional trends like those demonstrated here for the Woodford Shale.

393

The bulk elemental and FTIR data demonstrate, as expected, a substantial decrease in the

394

hydrogen content of Woodford kerogens during thermal maturation (Figures 5 and 7). The IR

395

structural indices show a decrease in the abundance of hydrogen-rich aliphatic carbons relative to

396

that of hydrogen-poor aromatic carbons (A-factor; Figure 7), as well as a shortening of the mean

397

length of remaining aliphatic chains by an increasing ratio of methyl to methylene groups

398

(CH3/CH2 ratio; Figure 7). The loss of hydrogen-rich aliphatic structural groups is not linear with 18 ACS Paragon Plus Environment

Page 18 of 49

Page 19 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

399

respect to increasing vitrinite reflectance. Rather, the magnitude is greatest at Ro values less than

400

or equal to 1.2% (immature to oil window [e.g.,

401

from kerogen at low thermal maturity is likely responsible for the generation of bitumen as a

402

precursor to generation of crude oils, as is indicated by laboratory artificial maturation

403

experiments.122, 123 Shown for comparison in Figure 5 are published H/C and vitrinite reflectance

404

data for type II kerogens represented by North American, European, and Middle East paleo-

405

marine organic-rich shales with ages from ~ 480 Ma to 90 Ma. Type II kerogens from these

406

different shales follow a H/C compositional trend with vitrinite reflectance that is

407

indistinguishable from that exhibited by the Woodford kerogens.40, 112, 124 The effect of thermal

408

stress on the cracking of carbon-carbon bonds appears to be largely similar for type II kerogen

409

globally.

121

]). This cracking of alkyl (C—C) structures

410

Insights to sulfur evolution during thermal maturation is provided by the bulk elemental

411

and sulfur XANES data. Sulfur has, despite its generally low abundance in type II kerogen,

412

particular significance for the rates of oil generation because of its proposed role in initiating free

413

radicals that promote cracking of alkyl C—C bonds at faster rate and lower thermal stress.114 The

414

sulfur data for the Woodford kerogens show an early and sharp change in both sulfur content and

415

organic-sulfur speciation as a function of thermal maturity. Total sulfur concentrations decrease

416

by ~ 50% relative from 3.1–3.5 wt% (atomic S/C ~ 0.015) to 1.3–2.1 wt% (atomic S/C ~ 0.01) at

417

vitrinite reflectance, Ro, values between 0.7% and 0.9%, but remain relatively uniform thereafter

418

(Figure 11). These data support the reactive nature of one or more sulfur species in kerogen, and

419

are at least consistent with the hypothesis that cracking of sulfur bonds occurs before and

420

initiates the cracking of alkyl C—C bonds, discussed above. The complementary sulfur XANES

421

data indicates that reactive sulfur is mostly aliphatic sulfide (~ 6 mol%), which is observed in

19 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

422

kerogens with Ro values of less than 0.7%, but not in kerogens of higher thermal maturity

423

(Figure 8). Sulfur species abundances, expressed as the ratio of reactive sulfur (aliphatic sulfide

424

± elemental) to thermally stable sulfur (aromatic thiophene), exhibits the same sharp and early

425

decrease as does the total sulfur content (Figure 11). Considering sulfur mass balance, the factor

426

of two reduction in total sulfur concentration cannot be accommodated only by the removal of

427

reactive sulfide. Thiophenic sulfur dominates the sulfur speciation of all Woodford kerogens

428

(83–89 mol% of organic sulfur), and some aromatic sulfur must be lost of kerogen by direct

429

cracking of aromatic S—C bonds or as collateral loss from cracking of lower-molecular-weight

430

fragments (i.e., bitumen, oil) in which thiophene is present. Indeed, the latter can explain the

431

predominance of thiophene as the major sulfur species in the heaviest molecular fractions of

432

petroleum (i.e., asphaltene).63 The narrow range of thiophene-dominated sulfur distributions of

433

the Woodford kerogens is similar to that documented for kerogen from the Silurian Qusaiba

434

Shale in Saudi Arabia spanning a similar range of thermal maturity (Ro ~ 0.9–2.0%).112 But, the

435

Woodford kerogens do not express as high sulfur concentrations (maximum atomic S/C ~ 0.03),

436

nor the wide range of thiophene abundances (30–90 mol%) measured by Kelemen et al. for other

437

type II kerogens from several shales of mostly lower thermal maturity (mean Rock-Eval Tmax
1.5%.20, 22, 25, 27 These studies have also shown a concomitant decreasing fraction

481

of aromatic carbon bonded to alkyl and heteroatom group attachments. Kelemen et al.25

482

documented for a diverse set of kerogens a robust, negative correlation between their fraction of

483

total carbon in aromatic structures from 13C NMR and XPS measurements and their bulk atomic

484

H/C ratio. The correlation is expected because aromatic carbon structures in kerogen have atomic

485

H/C ratios ≤ 1 whereas aliphatic carbon structures have atomic H/C ratios ≥ 2. Our regression

486

through this correlation [%Caromatic = 110 – 59 × H/C] allows us to estimate an equivalent

487

aromatic carbon fraction in the Woodford kerogens from their measured atomic H/C ratios. The

488

estimated %Caromatic values in the Woodford kerogens increases from ~ 45% at Ro = 0.55% to ~

489

70 % at Ro = 1.55% (Table 3), which is indistinguishable from the values calculated from NMR

490

spectroscopy for other type II kerogens of equal thermal maturity.20,

22 ACS Paragon Plus Environment

22, 25, 27

Extending this

Page 23 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

491

correlation to the most mature Woodford kerogen studied (atomic H/C = 0.45) predicts a

492

%Caromatic value of ~ 83% at Ro = 2.6% (Table 3). The relative trends in the abundance of

493

aromatic carbon and its attachments derived from the FTIR analysis are, for the set of Woodford

494

kerogens, consistent with the interpretations of changes in carbon structures for other type II

495

kerogens derived from NMR spectroscopy.

496

Substantial changes in the abundance of molecular functional groups in kerogen during

497

maturation, described above, are shown here to be complemented by evolution of its physical

498

properties. The skeletal density of the Woodford kerogen increases substantially with increasing

499

vitrinite reflectance, from ~ 1.2 g/cm3 to at least 1.55 g/cm3 (Figure 9). This trend is interpreted

500

to reflect the gross aliphatic versus aromatic character of kerogen, which evolves from aliphatic-

501

rich (analogous to alkyl structures such as waxes that have lower density of < 1 g/cm3) to

502

aromatic-rich (analogous to graphite that has higher density of ~ 2 g/cm3). The density versus

503

maturity trend shown by the Woodford kerogens is indistinguishable from that shown by other

504

type II kerogens from numerous organic-rich shales (Figure 9a).112,

505

chemical and physical measurements of the Woodford kerogens here shows that the absolute

506

density of kerogen is nearly linearly correlated to its atomic H/C content, and by extension its

507

aromatic carbon content (Figure 9b). Assuming this linear correlation holds for any given H/C

508

ratio, the inferred density of kerogen is less than that of pure graphite, reflecting presence of

509

heteroatoms and imperfect ordering of aromatic clusters in catagenic kerogen in contrast to

510

graphite that is composed of essentially infinite and perfectly ordered sheets of aromatic carbon.

117-119

Integrating the

511

Not only does the density of Woodford kerogens measurably change during maturation,

512

so does their internal specific surface area (SSA). Novel and quantitative insights to kerogen

513

pore structure is enabled here using combined acid demineralization and critical-point drying 23 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 49

514

methods that are shown to preserve kerogen microstructure.108 For these samples, SSA increases

515

by nearly a factor of ten between immature kerogen (~ 50 m2/g) and over-mature kerogen in the

516

dry-gas window (~ 300 m2/g; Figure 10). The calculated SSA values for the Woodford kerogens

517

overlap those of other type II kerogens prepared using the same critical-point drying methods,112

518

as well as with those calculated from surface area measurements of bulk shales after correcting

519

for inorganic porosity.130 The development of internal pore surfaces in Woodford kerogens is

520

consistent, for example, with X-ray and electron imaging studies75,

521

qualitatively show abundant nanometer- to micrometer-sized pores formed within kerogen

522

during maturation. This evolution, when interpreted in context of complementary chemical

523

evolution, is explained by the creation of porous volumes within a skeletal framework of

524

imperfectly ordered aromatic carbons, resulting from the cracking and expulsion of aliphatic-

525

dominated structures during maturation.

526

4. Summary and Conclusions

77-79, 95, 131, 132

that

527

A set of kerogens (defined here as the solid, insoluble, and non-volatile fraction of

528

sedimentary organic matter) spanning a maturity transect (vitrinite reflectance, Ro, between 0.5%

529

to 2.6%) across the Woodford Shale of the Anadarko Basin in Oklahoma (paleo-marine, type II

530

kerogen) was studied using a combination of spectroscopic and gas intrusion methods, providing

531

detailed insights to their bulk chemical, molecular, and physical evolution during geological

532

maturation. The following conclusions are made: Woodford kerogens follow the expected trend of substantially decreasing bulk

533



534

H/C content with increasing thermal maturity. IR spectra indicate that the loss of

535

hydrogen from kerogen is related to cracking of hydrogen-rich aliphatic (alkyl) chains

536

from aromatic carbon clusters. Within the oil-generation window (Ro ≤ ~ 1.4%), aromatic 24 ACS Paragon Plus Environment

Page 25 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

537

carbons remain highly substituted with alkyl (methyl, methylene) and probably

538

heteroatom (O, S, N) functional groups. At higher thermal stress (Ro > ~ 1.5%) these

539

aromatic substitutions are substantially removed and eventually replaced by hydrogen.

540

The structural evolution inferred from spectroscopic methods is consistent with known

541

carbon bond dissociation energies for carbonaceous materials.

542



543

have low sulfur concentrations (< 4 wt%), dominated by aromatic moieties (thiophene; ~

544

80–90 mol%). Reactive sulfur (aliphatic sulfide) is eliminated at low thermal maturity

545

(Ro ≤ 0.9%). At higher thermal stress, sulfur content and speciation remains stable.

546



547

manner consistent with their molecular characteristics. Kerogen skeletal density increases

548

by at least 30% relative (1.2 g/cm3 to 1.55 g/cm3) for the range of thermal maturity

549

studied. Density is linearly and negatively correlated with H/C ratio, and by inference is

550

linearly and positively correlated with aromatic carbon content, evolving from aromatic-

551

lean (i.e., alkyl-rich) to aromatic-rich (graphitic).

552



553

procedures to isolate and dry kerogen, and show a nearly order of magnitude increase (~

554

50 m2/g to 300 m2/g) for the range of thermal maturity studied. Increasing surface area is

555

consistent with development of organic-hosted pores within an imperfectly ordered

556

framework of aromatic carbons, resulting from preferential cracking and expulsion of

557

aliphatic carbons on the periphery of aromatic clusters during maturation and petroleum

558

generation.

Organic sulfur contents and sulfur XANES results show that Woodford kerogens

Physical properties of Woodford kerogens evolve during thermal maturation in a

Internal pore surface area measurements of kerogen are enabled by novel

25 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Globally, chemical and physical properties of type II kerogens from numerous

559



560

unconventional oil and gas shale plays follow maturity trends indistinguishable from that

561

exhibited by Woodford kerogen, with respect to carbon and hydrogen content, carbon

562

speciation, and density. Subtle differences may occur among type II kerogens with

563

respect to heteroatom (e.g., S, O) content and speciation, although their evolution with

564

maturity is similar.

565

Acknowledgements

566

We gratefully acknowledge Brian Cardott of the Oklahoma Geological Survey for

567

providing the Woodford Shale samples and their vitrinite reflectance data, for numerous

568

beneficial discussions and a review of an earlier draft of this manuscript, and for support to

569

publish these results. We appreciate the constructive reviews from three reviewers on our

570

original submission. We thank Schlumberger-Doll Research PetroLabs for assistance with

571

sample preparation and analysis.

26 ACS Paragon Plus Environment

Page 26 of 49

Page 27 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

572

Figure Captions

573

Fig. 1. Sulfur K-edge XANES spectra for sulfur standards, arranged top-to-bottom in order of

574

increasing sulfur oxidation state.

575

Fig. 2. Transmission FTIR and sulfur K-edge XANES spectra of asphaltene extracts from a

576

Middle East black oil (UG-8) before and after treatment with concentrated acids. The spectra and

577

the resulting IR and XANES structural indices of the two fractions are the same within

578

uncertainty, indicating that the acid treatment does not measurably impact the chemical

579

composition of solid carbonaceous materials.

580

Fig. 3. Rock-Eval Tmax and VRe maturity parameters plotted against vitrinite reflectance, Ro.

581

Fig. 4. Van Krevelen diagram showing that the Woodford kerogens follow the expected

582

evolution of atomic H/C and O/C ratios during thermal maturation for type II organic matter.

583

Fig. 5. Atomic H/C ratio of Woodford kerogen plotted versus their vitrinite reflectance, Ro.

584

Previously published data for other type II organic matter are plotted for comparison.40, 112, 124

585

Fig. 6. Transmission FTIR spectra of Woodford kerogens, arranged top-to-bottom in order of

586

increasing thermal maturity.

587

Fig. 7. IR structural indices (CH3/CH2, A-factor, Aromaticity, and C-factor) plotted against

588

measured vitrinite reflectance, Ro, and atomic H/C.

589

Fig. 8. Sulfur K-edge XANES spectra and calculated distribution of organic sulfur species for

590

Woodford kerogens, arranged top-to-bottom in order of increasing thermal maturity. Kerogen

591

sulfur species are dominated by thermally stable thiophene (sulfur in aromatic rings) and show

592

the early loss of reactive sulfide (sulfur in aliphatic chains).

593

Fig. 9. Skeletal density of Woodford kerogen plotted versus (a) vitrinite reflectance, Ro, and (b)

594

atomic H/C ratio. Published data for other type II kerogen are shown for comparison.112, 118, 119,

595

133

596

Fig. 10. Specific surface area of Woodford kerogen plotted versus vitrinite reflectance, Ro.

597

Published data for other type II kerogens from bulk shale and kerogen isolate analyses are shown

598

for comparison.112, 130 27 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

599

Fig. 11. Sulfur content (weight percent and atomic S/C) and sulfur speciation (molar ratio of

600

aliphatic sulfur [sulfide + elemental] to aromatic sulfur [thiophene]) of Woodford kerogen

601

plotted versus vitrinite reflectance, Ro, and atomic H/C.

602 603

References

604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643

1. Durand, B., Sedimentary organic matter and kerogen: Definition and quantitative importance of kerogen. In Kerogen—Insoluble Organic Matter from Sedimentary Rocks, Durand, B., Ed. Editions Technip: Paris, 1980; pp 13-34. 2. Dow, W. G., Kerogen studies and geological interpretations. Journal of Geochemical Exploration 1977, 7, 79-99. 3. Dutcher, R. R.; Hacquebard, P. A.; Schopf, J. M.; Simon, J. A., Carbonaceous Materials as Indicators of Metamorphism. In The Geological Society of America, Special Paper 153, 1974; p 116. 4. Corcoran, D. V.; Doré, A. G., A review of techniques for the estimation of magnitude and timing of exhumation in offshore basins. Earth Science Reviews 2006, 72, 129-168. 5. Pawlewicz, M. J.; King, D. J., Vitrinite and solid bitumen reflectance: Some correlations and applications. In The Petroleum System–Status of Research and Methods, Magoon, L. B., Ed. U.S. Geological Survery Bulletin 2007: 1992; pp 58-60. 6. Hackley, P. C.; Cardott, B. J., Application of organic petrography in North American shale petroleum systems: A review. International Journal of Coal Geology 2016, 163, 8-51. 7. Clementz, D. M.; Demaison, G. J.; Daly, R. A. In Well site geochemistry by programmed pyrolysis, Proceedings of the 11th Annual Offshore Technology Conference, 1979; 1979; pp 465-470. 8. Peters, K. E., Guidelines for evaluating petroleum source rock using programmed pyrolysis. AAPG Bulletin 1986, 70, 318-329. 9. Behar, F.; Beaumont, V.; Penteado, H. L. D., Rock-Eval 6 technology: Performances and developments. Oil & Gas Science and Technology - Rev. IFP 2001, 56, 111-134. 10. Van Krevelen, D. W., Graphical-statistical method for the study of structure and reaction processes of coal. Fuel 1950, 29, 269-283. 11. Beyssac, O.; Goffé, B.; Petitet, J.-P.; Froigneux, E.; Moreau, M.; Rouzaud, J.-N., On the characterization of disordered and heterogeneous carbonaceous materials by Raman spectroscopy. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 2003, 59, 2267-2276. 12. Kelemen, S. R.; Fang, H. L., Maturity trends in Raman spectra from kerogen and coal. Energy & Fuels 2001, 15, 653-658. 13. Wopenka, B.; Pasteris, J. D., Structural characterization of kerogens to granulite-facies graphite: Applicability of Raman microprobe spectroscopy. American Mineralogist 1993, 78, 533-557. 14. Liu, Y.; Ferralis, N.; Bryndzia, L. T.; Grossman, J. C., Genome-inspired molecular identification in organic matter via Raman spectroscopy. Carbon 2016, 101, 361-367. 15. Sauerer, B.; Craddock, P. R.; AlJohani, M. D.; Alsamadony, K. L.; Abdallah, W., Fast and accurate shale maturity determination by Raman spectroscopy measurement without sample preparation. International Journal of Coal Geology 2017, 173, 150-157. 16. Schito, A.; Romano, C.; Corrado, S.; Grigo, D.; Poe, B., Diagenetic thermal evolution of organic matter by Raman spectroscopy. Organic Geochemistry 2017, 106, 57-67. 28 ACS Paragon Plus Environment

Page 28 of 49

Page 29 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693

Energy & Fuels

17. Schmidt Mumm, A.; İnan, S., Microscale organic maturity determination of graptolites using Raman spectroscopy. International Journal of Coal Geology 2016, 162, 96-107. 18. Spötl, C.; Houseknecht, D. W.; Jaques, R. C., Kerogen maturation and incipient graphitization of hydrocarbon source rocks in the Arkoma Basin, Oklahoma and Arkansas: A combined petrographic and Raman spectrometric study. Organic Geochemistry 1998, 28, 535542. 19. Burdelnaya, N.; Bushnew, D.; Mokeev, M.; Dobrodumov, A., Experimental study of kerogen maturation by solid-state 13C NMR spectroscopy. Fuel 2014, 118, 308-315. 20. Cao, X.; Birdwell, J. E.; Chappell, M. A.; Li, Y.; Pignatello, J. J.; Mao, J., Characterization of oil shale, isolated kerogen, and postpyrolysis residues using advanced 13C nuclear magnetic resonance spectroscopy. AAPG Bulletin 2013, 97, 421-436. 21. Clough, A.; Sigle, J. L.; Jacobi, D.; Sheremata, J.; White, J. L., Characterization of kerogen and source rock maturation using solid-state NMR spectroscopy. Energy & Fuels 2015, 29, 6370-6382. 22. Mao, J.; Fang, X.; Lan, Y.; Schimmelmann, A.; Mastalerz, M.; Xu, L.; Schmidt-Rohr, K., Chemical and nano-meter schale structure of kerogen and its change during thermal maturation investigated by advanced solid-state 13C NMR spectroscopy. Geochimica et Cosmochimica Acta 2010, 74, 2110-2127. 23. Petsch, S. T.; Smernik, R. J.; Eglinton, T. I.; Oades, J. M., A solid state 13C-NMR study of kerogen degradation during black shale weathering. Geochimica et Cosmochimica Acta 2001, 65, 1867-1882. 24. Solum, M.; Pugmire, R. J.; Grant, D. M., 13C solid-state NMR of Argonne premium coals. Energy & Fuels 1989, 3, 187-193. 25. Kelemen, S. R.; Afeworki, M.; Gorbaty, M. L.; Sansone, M.; Kwiatek, P. J.; Walters, C. C.; Freund, H.; Siskin, M.; Bence, A. E.; Curry, D. J.; Solum, M.; Pugmire, R. J.; Vandenbroucke, M.; Leblond, M.; Behar, F., Direct characterization of kerogen by X-ray and solid-state 13C nuclear magnetic resonance methods. Energy & Fuels 2007, 21, 1548-1561. 26. Tong, J.; Han, X.; Wang, S.; Jiang, X., Evaluation of structural characteristics of Huadian oil shale kerogen using direct techniques (solid-state 13C NMR, XPS, FT-IR, and XRD). Energy & Fuels 2001, 25, 4006-4013. 27. Witte, E. G.; Schenk, H. J.; Müller, P. J.; Schwochau, K., Structural modifications of kerogen during natural evolution as derived from 13C CP/MAS NMR, IR spectroscopy and RockEval pyrolysis of Toarcian shales. Organic Geochemistry 1988, 13, 1039-1044. 28. Cao, X.; Yang, J.; Mao, J., Characterization of kerogen using solid-state nuclear magnetic resonance spectroscopy: A review. International Journal of Coal Geology 2013, 108, 83-90. 29. Mao, J.; Cao, X.; Olk, D. C.; Chu, W.; Schmidt-Rohr, K., Advanced solid-state NMR spectroscopy of natural organic matter. Progress in Nuclear Magnetic Resonance Spectroscopy 2017, 100, 17-51. 30. Fuller, M. P.; Hamadeh, I. M.; Griffiths, P. R.; Lowenhaupt, D. E., Diffuse reflectance infrared spectrometry of powdered coals. Fuel 1982, 61, 529-536. 31. Wang, S.-H.; Griffiths, P. R., Resolution enhancement of diffuse reflectance i.r. spectra of coals by Fourier self-deconvolution. 1. C—H stretching and bending modes. Fuel 1985, 64, 229-236. 32. Guo, Y.; Bustin, R. M., Micro-FTIR spectroscopy of liptinite macerals in coal. International Journal of Coal Geology 1998, 36, 259-275. 33. Ibarra, J. V.; Munoz, E.; Moliner, R., FTIR study of the evolution of coal structure during the coalification process. Organic Geochemistry 1996, 24, 725-735. 34. Iglesias, M.; Jiménez, A.; Laggoun-Défarge, F.; Suárez-Ruiz, I., FTIR study of pure vitrains and associated coals. Energy & Fuels 1995, 9, 458-466.

29 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742

35. Painter, P. C.; Snyder, R. W.; Starsinic, M.; Coleman, M. M.; Kuehn, D. W.; Davis, A., Concerning the application of FTIR to the study of coal: A critical assessment of band assignments and the application of spectral analysis programs. Applied Spectroscopy 1981, 35, 475-485. 36. Lin, R.; Ritz, G. P., Studying individual macerals using i.r. microspectroscopy, and implications on oil versus gas/condensate proneness and “low-rank” generation. Organic Geochemistry 1993, 20, 697-706. 37. Chen, Y.; Mastalerz, M.; Schimmelmann, A., Characterization of chemical functional groups in macerals across different coal ranks via micro-FTIR spectroscopy. International Journal of Coal Geology 2012, 104, 22-33. 38. Durand, B.; Espitalié, J., Geochemical studies on the organic matter from the Douala Basin (Cameroon)—II. Evolution of kerogen. Geochimica et Cosmochimica Acta 1976, 40, 801808. 39. Ganz, H.; Kalkreuth, W., Application of infrared spectroscopy to the classification of kerogen-types and the evaluation of source rock and oil shale potentials. Fuel 1987, 66, 708711. 40. Lis, G. P.; Mastalerz, M.; Schimmelmann, A.; Lewan, M. D.; Stankiewicz, B. A., FTIR absorption indices for thermal maturity in comparison with vitrinite reflectance Ro in type-II kerogen from Devonian black shales. Organic Geochemistry 2005, 36, 1533-1552. 41. Rouxhet, P. G.; Robin, P. L., Infrared study of the evolution of kerogens of different origins during catagenesis and pyrolysis. Fuel 1978, 57, 533-540. 42. Rouxhet, P. G.; Robin, P. L.; Nicaise, G., Characterization of kerogens and of their evolution by infrared spectroscopy. In Kerogen—Insoluble Organic Matter from Sedimentary Rocks, Durand, B., Ed. Editions Technip: Paris, 1980; pp 163-190. 43. Craddock, P. R.; Le Doan, T. V.; Bake, K. D.; Polyakov, M.; Charsky, A. M.; Pomerantz, A. E., Evolution of kerogen and bitumen during thermal maturation by semi-open pyrolysis investigated by infrared spectroscopy. Energy & Fuels 2015, 29, 2197-2210. 44. Craddock, P. R.; Prange, M.; Pomerantz, A. E., Kerogen thermal maturity and content of organic-rich mudrocks determined using stochastic linear regression models applied to diffuse reflectance IR Fourier transform spectroscopy (DRIFTS). Organic Geochemistry 2017, 110, 122-133. 45. Christy, A. A.; Kvalheim, O. M.; Oygard, K.; Dahl, B.; Karstang, T. V., Maturity of kerogen and asphaltenes determine by partial-least-squares (PLS) calibration and target projection of diffuse reflectance Fourier transformed infrared spectra. Organic Geochemistry 1991, 17, 63-74. 46. Christy, A. A.; Hopland, A. L.; Barth, T.; Kvalheim, O. M., Quantitative determination of thermal maturity in sedimentary organic matter by diffuse reflectance infrared spectroscopy of asphaltenes. Organic Geochemistry 1989, 14, 77-81. 47. Barth, T.; Seim, M.; Christy, A. A.; Kvalheim, O. M., Maturity trends in asphaltenes from pyrolyzed source rocks and natural coals-multivariate modelling of diffuse reflectance Fouriertransform infrared spectra. Organic Geochemistry 1995, 23, 139-158. 48. Spiro, C. L.; Wong, J.; Lytle, F. W.; Greegor, R. B.; Maylotte, D. H.; Lamson, S. H., X-ray absorption spectroscopic investigations of sulfur sites in coal: Organic sulfur identification. Science 1984, (226), 48-50. 49. Huffman, G. P.; Mitra-Kirtley, S.; Huggins, F. E.; Shah, N.; Vaidya, S.; Lu, F., Quantitative analysis of all major forms of sulfur in coal by X-ray absorption fine structure spectroscopy. Energy & Fuels 1991, 5, 574-581. 50. Gorbaty, M. L.; Kelemen, S. R.; George, G. N.; Kwiatek, P. J., Characterization and thermal reactivity of oxidized organic sulfur forms in coals. Fuel 1992, 71, 1255-1264.

30 ACS Paragon Plus Environment

Page 30 of 49

Page 31 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793

Energy & Fuels

51. Kelemen, S. R.; Afeworki, M.; Gorbaty, M. L.; Kwiatek, P. J.; Sansone, M.; Walters, C. C.; Cohen, A. D., Thermal transformations of nitrogen and sulfur forms in peat related to coalification. Energy & Fuels 2006, 20, 635-652. 52. Bolin, T. B., Direct determination of pyrite content in Argonne Premium Coals by the use of sulfur X-ray near edge absorption spectroscopy (XANES). Energy & Fuels 2010, 24, 54795482. 53. Wiltfong, R.; Mitra-Kirtley, S.; Mullins, O. C.; Andrews, B.; Fujisawa, G.; Larsen, J. W., Sulfur speciation in different kerogens by XANES spectroscopy. Energy & Fuels 2005, 19, 1971-1976. 54. Kelemen, S. R.; Sansone, M.; Walters, C. C.; Kwiatek, P. J.; Bolin, T. B., Thermal transformations of organic and inorganic sulfur in Type II kerogen quantified by S-XANES. Geochimica et Cosmochimica Acta 2012, 83, 61-78. 55. Wang, Q.; Liu, Q.; Wang, Z.-C.; Liu, H.-P.; Bai , J.-R.; Ye, J.-B., Characterization of organic nitrogen and sulfur in the oil shale kerogens. Fuel Processing Technology 2017, 160, 170-177. 56. Pomerantz, A. E.; Bake, K. D.; Craddock, P. R.; Kurzenhauser, K. W.; Kodalen, B. G.; Mitra-Kirtley, S.; Bolin, T. B., Sulfur speciation in kerogen and bitumen from gas and oil shales. Organic Geochemistry 2014, 68, 5-12. 57. Kelemen, S. R.; Walters, C. C.; Kwiatek, P. J.; Freund, H.; Afeworki, M.; Sansone, M.; Lamberti, W. A.; Pottorf, R. J.; Machel, H. G.; Peters, K. E., Characterization of solid bitumens originating from thermal chemical alteration and thermochemical sulfate reduction. Geochimica et Cosmochimica Acta 2010, 74, 5305-5332. 58. Bolin, T. B.; Birdwell, J. E.; Lewan, M. D.; Hill, R. J.; Grayson, M. B.; Mitra-Kirtley, S.; Bake, K. D.; Craddock, P. R.; Abdallah, W.; Pomerantz, A. E., Sulfur species in source rock bitumen before and after hydrous pyrolysis determined by XANES. Energy & Fuels 2016, 30, 6264-6270. 59. George, G. N.; Gorbaty, M. L., Sulfur K-edge X-ray absorption spectroscopy of petroleum asphaltenes and model compounds. Journal of the American Chemical Society 1989, 111, 3182-3186. 60. Waldo, G. S.; Carlson, R. M. K.; Moldowan, J. M.; Peters, K. E.; Penner-Hahn, J. E., Sulfur speciation in heavy petroleums: Information from X-ray absorption near-edge structure. Geochimica et Cosmochimica Acta 1991, 55, 801-814. 61. Waldo, G. S.; Mullins, O. C.; Penner-Hahn, J. E.; Cramer, S. P., Determination of the chemical environment of sulphur in petroleum asphaltenes by X-ray absorption spectroscopy. Fuel 1992, 71, 53-57. 62. Mitra-Kirtley, S.; Mullins, O. C.; Ralston, C. Y.; Sellis, D.; Pareis, C., Determination of sulfur species in asphaltene, resin, and oil fractions of crude oils. Applied Spectroscopy 1998, 52, 1522-1525. 63. Pomerantz, A. E.; Seifert, D. J.; Bake, K. D.; Craddock, P. R.; Mullins, O. C.; Kodalen, B. G.; Mitra-Kirtley, S.; Bolin, T. B., Sulfur chemistry of asphaltenes from a highly compositionally graded oil column. Energy & Fuels 2013, 27, 4604-4608. 64. Pomerantz, A. E.; Le Doan, T. V.; Craddock, P. R.; Bake, K. D.; Kleinberg, R. L.; Burnham, A. K.; Qinghao, W.; Zare, R.; Brodnik, G.; Lo, C.-H.; Grayson, M.; Mitra-Kirtley, S.; Bolin, T. B.; Wu, T., Impact of laboratory-induced thermal maturity on asphaltene molecular structure. Energy & Fuels 2016, 30, 7025-7036. 65. Mitra-Kirtley, S.; Mullins, O. C.; Branthaver, J. F.; Cramer, S. P., Nitrogen chemistry of kerogens and bitumens from X-ray absorption near-edge structure spectroscopy. Energy & Fuels 1993, 7, 1128-1134. 66. Vairavamurthy, A.; Wang, S., Organic nitrogen in geomacromolecules:  Insights on speciation and transformation with K-edge XANES spectroscopy. Environmental Science and Technology 2002, 36, 3050-3056. 31 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844

67. Kelemen, S. R.; Afeworki, M.; Gorbaty, M. L.; Kwiatek, P. J.; Solum, M. S.; Hu, J. Z.; Pugmire, R. J., XPS and 15N NMR Study of Nitrogen Forms in Carbonaceous Solids. Energy & Fuels 2002, 16, 1507-1515. 68. Kelemen, S. R.; Afeworki, M.; Gorbaty, M. L.; Cohen, A. D., Characterization of organically bound oxygen forms in lignites, peats, and pyrolyzed peats by X-ray photoelectron spectroscopy (XPS) and solid-state 13C NMR methods. Energy & Fuels 2002, 16, 1450-1462. 69. Kelemen, S. R.; Kwiatek, P. J., Quantification of organic oxygen species on the surface of fresh and reacted Argonne Premium coal. Energy & Fuels 1995, 9, 841-848. 70. Kelemen, S. R.; Freund, H.; Gorbaty, M. L.; Kwiatek, P. J., Thermal chemistry of nitrogen in kerogen and low-rank coal. Energy & Fuels 1999, 13, 529-538. 71. Stańczyk, K.; Dziembaj, R.; Piwowarska, Z.; Witkowski, S., Transformation of nitrogen structures in carbonization of model compounds determined by XPS. Carbon 1995, 33, 13831392. 72. Pels, J. R.; Kapteijn, F.; Moulijn, J. A.; Zhu, Q.; Thomas, K. M., Evolution of nitrogen functionalities in carbonaceous materials during pyrolysis. Carbon 1995, 33, 1641-1653. 73. Pomerantz, A. E., Toward molecule-specific geochemistry of heavy ends: Applications to the upstream oil industry. Industrial & Engineering Chemistry Research 2017, 55, 4403-4414. 74. Loucks, R. G.; Reed, R. M.; Ruppel, S. C.; Jarvie, D. M., Morphology, genesis, and distribution of nanometer-scale pores in siliceous mudstones of the Mississippian Barnett Shale. Journal of Sedimentary Research 2009, 79, 848-861. 75. Loucks, R. G.; Reed, R. M.; Ruppel, S. C.; Hammes, U., Spectrum of pore types and networks in mudrocks and a descriptive classification of matrix-related mudrock types. AAPG Bulletin 2012, 96, 1071-1098. 76. Loucks, R. G.; Reed, R. M., Scanning-electron-microscope petrographic evidence for distinguishing organic-matter pores associated with depositional organic matter versus migrated organic matter in mudrocks. GCAGS Journal 2014, 3, 51-60. 77. Curtis, M. E.; Cardott, B. J.; Sondergeld, C. H.; Rai, C. S., Development of organic porosity in the Woodford Shale with increasing thermal maturity. International Journal of Coal Geology 2012, 103, 26-31. 78. Cardott, B. J.; Landis, C. R.; Curtis, M. E., Post-oil solid bitumen network in the Woodford Shale, USA—A potential primary migration pathway. International Journal of Coal Geology 2015, 139, 106-113. 79. Cardott, B. J.; Curtis, M. E., Identification and nanoporosity of macerals in coal by scanning electron microscopy. International Journal of Coal Geology 2017, In press. 80. Ross, D. J. K.; Bustin, R. M., The importance of shale composition and pore structure upon gas storage potential of shale gas reservoirs. Marine and Petroleum Geology 2009, 26, 916-927. 81. Rexer, T. F.; Mathia, E. J.; Aplin, A. C.; Thomas, K. M., High-pressure methane adsorption and characterization of pores in Posidonia shales and isolated kerogens. Energy & Fuels 2014, 28, 2886-2901. 82. Topór, T.; Derkowski, A.; Ziemiański, P.; Szczurowski, J.; McCarty, D. K., The effect of organic matter maturation and porosity evolution on methane storage potential in the Baltic Basin (Poland) shale-gas reservoir. International Journal of Coal Geology 2017, 180, 46-56. 83. Mastalerz, M.; Schimmelmann, A.; Drobniak, A.; Chen, Y., Porosity of Devonian and Mississippian New Albany Shale across a maturation gradient: Insights from organic petrology, gas adsorption, and mercury intrusion. AAPG Bulletin 2013, 97, 1621-1643. 84. Mastalerz, M.; He, L.; Melnichenko, Y. B.; Rupp, J. A., Porosity of coal and shale: Insights from gas adsorption and USANS techniques. Energy & Fuels 2012, 26, 5109-5120. 85. Radlinski, A. P.; Boreham, C. J.; Wignall, G. D.; Lin, J. S., Microstructural evolution of source rocks during hydrocarbon generation: A small-angle-scattering study. Physical Review B 1996, 53, 14152-14160. 32 ACS Paragon Plus Environment

Page 32 of 49

Page 33 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894

Energy & Fuels

86. Ruppert, L. F.; Sakurovs, R.; Blach, T. P.; He, L.; Melnichenko, Y. B.; Mildner, D. F. R.; Alcantar-Lopez, L., A USANS/SANS study of the accessibility of pores in the Barnett Shale to methane and water. Energy & Fuels 2013, 27, 772-779. 87. Thomas, J. J.; Valenza, J. J.; Craddock, P. R.; Bake, K. D.; Pomerantz, A. E., The neutron scattering length density of kerogen and coal as determined by CH3OH/CD3OH exchange. Fuel 2014, 117, 801-808. 88. Clarkson, C. R.; Freeman, M.; He, L.; Agamalian, M.; Melnichenko, Y. B.; Mastalerz, M.; Bustin, R. M.; Radliński, A. P.; Blach, T. P., Characterization of tight gas reservoir pore structure using USANS/SANS and gas adsorption analysis. Fuel 2012, 95, 371-385. 89. Tissot, B. P.; Welte, D. H., Petroleum Formation and Occurrence, 2nd Ed. Springer: Berlin Heidelberg, Germany, 1984. 90. Van Krevelen, D. W., Organic Geochemistry—old and new. Organic Geochemistry 1984, 6, 1-10. 91. Vandenbroucke, M.; Largeau, C., Kerogen origin, evolution and structure. Organic Geochemistry 2007, 38, 719-833. 92. Lille, Ü.; Heinmaa, I.; Pehk, T., Molecular model of Estonian kukersite kerogen evaluated 13 by C MAS NMR spectra. Fuel 2003, 82, 799-804. 93. Collell, J.; Ungerer, P.; Galliero, G.; Yiannourakou, M.; Montel, F. o.; Pujol, M., Molecular simulation of bulk organic matter in type II shales in the middle of the oil formation window. Energy & Fuels 2014, 28, 7457-7466. 94. Ungerer, P.; Collell, J.; Yiannourakou, M., Molecular modeling of the volumetric and thermodynamic properties of kerogen: Influence of organic type and maturity. Energy & Fuels 2014, 29, 91-105. 95. Bousige, C.; Ghimbeu, C. M.; Vix-Guterl, C.; Pomerantz, A. E.; Suleimenova, A.; Vaughan, G.; Garbarino, G.; Feygenson, M.; Wildgruber, C.; Ulm, F.-J.; Pellenq, R. J. M.; Coasne, B., Realistic molecular model of kerogen's nanostructure. Nature Materials 2016, 15, 576-583. 96. Falk, K.; Pellenq, R.; Ulm, F. J.; Coasne, B., Effect of chain length and pore accessibility on alkane adsorption in kerogen. Energy & Fuels 2015, 29, 7889-7896. 97. Behar, F.; Vandenbroucke, M., Chemical modelling of kerogen. Organic Geochemistry 1987, 11, 15-24. 98. Curiale, J. A., Origin of solid bitumens, with emphasis on biological marker results. Organic Geochemistry 1986, 10, (1), 559-580. 99. Yang, J.; Hatcherian, J.; Hackley, P. C.; Pomerantz, A. E., Nanoscale geochemical and geomechanical characterization of organic matter in shale. Nature Communications 2017, 8, 2179. 100. Johnson, K. S.; Cardott, B. J., Geologic framework and hydrocarbon source rocks of Oklahoma. In Source Rocks in the Southern Midcontinent, 1990 Symposium: Oklahoma Geological Survey Circular 93, Johnson, K. S.; Cardott, B. J., Eds. 1992; pp 21-37. 101. Northcutt, R. A.; Campbell, J. A., Geologic provinces of Oklahoma. In Basement Tectonics 12, Central North America and Other Regions. Proceedings of the 12th International Conference on Basement Tectonics, Hogan, J. P.; Gilbert, M. C., Eds. Kluwer Academic Publishers, Boston: 1998; pp 29-37. 102. Hoffman, P.; Dewey, J. F.; Burke, K., Aulacogens and their genetic relation to geosynclines, with a Proterozoic examples from Great Slave lake, Canada. In Moden and Ancient Geosynclinal Sedimentation: SEPM Special Publication 19, 1974; pp 38-55. 103. Cardott, B. J., Thermal maturity of Woodford Shale gas and oil plays, Oklahoma, USA. International Journal of Coal Geology 2012, 103, 109-119. 104. Cardott, B. J.; Lambert, M. W., Thermal maturation by vitrinite reflectance of Woodford Shale, Anadarko Basin, Oklahoma. AAPG Bulletin 1985, 69, 1982-1998.

33 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944

105. Comer, J. B., Organic geochemistry and paleogeography of Upper Devonian formations in Oklahoma and western Arkansas. In Source Rocks in the Southern Midcontinent, 1990 Symposium: Oklahoma Geological Survey Circular 93, Johnson, K. S.; Cardott, B. J., Eds. 1992; pp 70-93. 106. Cardott, B. J., Thermal maturation of the Woodford Shale in the Anadarko basin. In Anadarko Basin Symposium 1988: Oklahoma Geological Survey Circular 90, Johnson, K. S., Ed. 1989; pp 32-46. 107. Ibrahimov, R. A.; Bissada, K. K., Comparative analysis and geological significance of kerogen isolated using open-system (palynological) versus chemically and volumetrically conservative closed-system methods. Organic Geochemistry 2010, 41, 800-811. 108. Suleimenova, A.; Bake, K. D.; Ozkan, A.; Valenza, J. J.; Kleinberg, R. L.; Burnham, A. K.; Ferralis, N.; Pomerantz, A. E., Acid demineralization with critical point drying: A method for kerogen isolation that preserves microstructure. Fuel 2014, 135, 492-497. 109. Jackson, L., L.; Roof, S. R., Determination of the forms of carbon in geological materials. Geostandards Newsletter 1992, 16, 317-323. 110. Jarvie, D. M.; Claxton, B. L.; Henk, F.; Breyer, J. T., Oil and shale gas from the Barnett Shale, Ft. Worth Basin, Texas (abs). AAPG Annual Meeting, Program with Abstracts, p. A100 2001. 111. Durand, B.; Monin, J. C., Elemental analysis of kerogens (C, H, O, N, S, Fe). In Kerogen—Insoluble Organic Matter from Sedimentary Rocks, Durand, B., Ed. Editions Technip: Paris, 1980; pp 113-142. 112. Cheshire, S.; Craddock, P. R.; Xu, G.; Sauerer, B.; Pomerantz, A. E.; McCormick, D.; Abdallah, W., Assessing thermal maturity beyond the reaches of vitrinite reflectance and RockEval pyrolysis: A case study from the Silurian Qusaiba formation. International Journal of Coal Geology 2017, 180, 29-45. 113. Greenfield, M. L.; Byrne, M.; Mitra-Kirtley, S.; Kercher, E. M.; Bolin, T. B.; Wu, T.; Craddock, P. R.; Bake, K. D.; Pomerantz, A. E., XANES measurements of sulfur chemistry during asphalt oxidation. Fuel 2015, 162, 179-185. 114. Lewan, M. D., Sulphur-radical control on petroleum formation rates. Nature 1998, 391, 164-166. 115. Zhang, J.; Clennell, M. B.; Dewhurst, D. N.; Liu, K., Combined Monte Carlo and molecular dynamics simulation of methane adsorption on dry and moist coal. Fuel 2014, 122, 186-197. 116. Zhao, T.; Li, X.; Zhao, H.; Li, M., Molecular simulation of adsorption and thermodynamic properties on type II kerogen: Influence of maturity and moisture content. Fuel 2017, 190, 198207. 117. Guidry, K.; Luffel, D.; Curtis, J. Development of Laboratory and Petrophysical Techniques for Evaluating Shale Reservoirs, Final Report; Gas Research Institute, Chicago, Illinois: 1995; p 304. 118. Alfred, D.; Vernik, L., A new petrophysical model for organic shales. Petrophysics 2013, 54, 240-247. 119. Okiongbo, K. S.; Aplin, A. C.; Larter, S. R., Changes in Type II kerogen density as a function of maturity: Evidence from the Kimmeridge Clay Formation. Energy & Fuels 2005, 19, 2495-2499. 120. Pelet, R.; Behar, F.; Monin, J. C., Resins and asphaltenes in the generation and migration of petroleum. Organic Geochemistry 1986, 10, 481-498. 121. Dembicki Jr., H., Three common source rock evaluation errors made by geologists during prospecct or play appraisals. AAPG Bulletin 2009, 93, 341-356. 122. Lewan, M. D., Experiments on the role of water in petroleum formation. Geochimica et Cosmochimica Acta 1997, 61, 3691-3723.

34 ACS Paragon Plus Environment

Page 34 of 49

Page 35 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972

Energy & Fuels

123. Le Doan, T. V.; Bostrom, N. W.; Burnham, A. K.; Kleinberg, R. L.; Pomerantz, A. E.; Allix, P., Green River oil shale pyrolysis: Semi-open conditions. Energy & Fuels 2013, 27, 74476459. 124. Buchardt, B.; Lewan, M. D., Reflectance of vitrinite-like macerals as a thermal maturity index for Cambrian-Ordovician Alum Shale, southern Scandinavia. AAPG Bulletin 1990, 74, 394-406. 125. Tissot, B.; Deroo, G.; Hood, A., Geochemical study of the Uinta Basin: formation of petroleum from the Green River formation. Geochimica et Cosmochimica Acta 1978, 42, 14691485. 126. Huss, E. B.; Burnham, A. K., Gas evolution during pyrolysis of various Colorado oil shales. Fuel 1982, 61, 1188-1196. 127. Burnham, A. K.; Clarkson, J. E.; Singleton, M. F.; Wong, C. M.; Crawford, R. W., Biological markers from Green River kerogen decomposition. Geochimica et Cosmochimica Acta 1982, 46, 1243-1261. 128. Sanderson, R. T., Chemical Bonds and Bond Energy. Academic Press: New York, N.Y., 1976. 129. Luo, Y.-R., Comprehensive Handbook of Chemical Bond Energies. CRC Press: Boca Raton, FL, 2007. 130. Valenza, J. J.; Drenzek, N.; Marques, F.; Pagels, M.; Mastalerz, M., Geochemical controls on shale microstructre. Geology 2013, 41, 611-614. 131. Bernard, S.; Wirth, R.; Schreiber, A.; Schulz, H.-M.; Horsfield, B., Formation of nanoporous pyrobitumen residues during maturation of the Barnett Shale (Fort Worth Basin). International Journal of Coal Geology 2012, 103, 3-11. 132. Milliken, K. L.; Rudnicki, M.; Awwiller, D. N.; Zhang, T., Organic matter-hosted pore system, Marcellus Formation (Devonian), Pennsylvania. AAPG Bulletin 2013, 97, 177-200. 133. Guidry, F. K.; Luffel, D. L.; Olszewskl, A. J., Devonian shale formation evaluation models based on logs, new core analysis methods, and production tests. In Transactions of the SPWLA 31st Annual Logging Symposium, June 24-27, Lafayette, Louisiana, USA, 1990.

973

35 ACS Paragon Plus Environment

Energy & Fuels

12

FeS2

9

Pyrite

S8

Elemental

S S

I/I0

Sulfide

6

Thiophene

S

O

O

OH H5C2O

Sulfinic

3

OH

S

Sulfoxide

NH2

O

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 49

OC2H5

S

O

O

Sulfite

H3C

S

O

O

O

Sulfone

H3C

Sulfonate

O-

S O S

CH3

CH3

O O-

O

Sulfate

2465

2470

2475

2480

Photon Energy, eV

2485

Fig. 1 1.5-column width (request figure caption right of figure, not below) ACS Paragon Plus Environment

Page 37 of 49

IR absorption, a.u.

0.8

FTIR structural indices

1.0

original acid-treated

1.0

0.8 0.6 0.4

0.6

0.2

0.4

0.0

CH3/CH2

0.2 0.0 500

1000

1500

2000

2500

Wavenumber, cm-1

3000

3.0

Aromaticity

C-Factor

3500

Organic sulfur species, mole percent

Acid-treated

2.0

0%

1.5

20%

elemental

1.0

sulfide

0.5 0.0 2460

A-Factor

Original

2.5

XANES I/I0

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

2465

2470

2475

2480

2485

2490

Photon energy, eV

Fig. 2 Double-column width ACS Paragon Plus Environment

40%

60%

thiophene

80%

sulfoxide

100% sulfone

Energy & Fuels

480

Rock-Eval Tmax, oC

460 440 420 400 380 4

Rock-Eval VRe, %

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 49

3 2 1 0

0

1

2

Ro, %

3

4

Fig. 3 Single-column width ACS Paragon Plus Environment

Page 39 of 49

type I kerogen

(Peters, 1986)

1.5

Atomic H/C

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

type II kerogen 0.55 %Ro 0.66 %Ro

1.0

type III kerogen 0.90 %Ro 1.16 %Ro 1.55 %Ro

0.5

1.82 %Ro 2.59 %Ro

0.0 0.0

0.1

0.2

Atomic O/C

0.3

Fig. 4 Single-column width

ACS Paragon Plus Environment

Energy & Fuels

1.75

Buchardt & Lewan 1990 Lis et al. 2005 Cheshire et al. 2017

1.5

Atomic H/C

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 49

1.25 1.0 0.75 0.5 0.5

1

1.5

2

2.5

Ro, %

Fig. 5 Single-column width

ACS Paragon Plus Environment

Page 41 of 49

aliphatic CH2, CH3

aliphatic CH2, CH3

Absorption, a.u. (spectra offset)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

aliphatic CH3

aromatic CC aromatic CH

CO

0.55 %Ro 0.90 %Ro 1.16 %Ro 1.82 %Ro 2.59 %Ro

500

1000 1500 2000 2500 3000 3500

Wavenumber, cm-1

Fig. 6 Single-column width

ACS Paragon Plus Environment

Energy & Fuels

CH3/CH2

1.25

2.59 %Ro 1.82 %Ro

1

1.16 %Ro 0.90 %Ro

0.75 0.55 %Ro

A-Factor

0.5 1 0.75 0.5 0.25

Aromaticity

0 0.5 0.4 0.3 0.2 0.1 0 0.6 0.5

C-Factor

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 49

0.4 0.3 0.2 0.1

0

1

2

Ro, %

3 1.2

1.0

0.8

0.6 0.4

Atomic H/C

Fig. 7 Single-column width ACS Paragon Plus Environment

Page 43 of 49

elemental

Ro, %

0.55 %Ro

Fluorescence, I/I0

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

0.66 %Ro 0.90 %Ro 1.16 %Ro 1.82 %Ro

sulfide

thiophene

sulfoxide

sulfone

0.55 7 6

83

4

0.66 7 6

83

4

0.90 8

89

3

1.16 8

89

3

1.82 8

88

4

2.59 8

83

0

20

40

7

60

80

100

Organic sulfur species concentration, mol%

2.59 %Ro

2455 2460 2465 2470 2475 2480 2485 2490

Photon energy, eV

Fig. 8 Double-column width ACS Paragon Plus Environment

Energy & Fuels

Kerogen density, g/cm3

1.8

(a)

1.6 1.4 1.2

Guidry et al., 1995 Okiongbo et al., 2005 Alfred & Vernik, 2013 Cheshire et al., 2017

1.0 0.5

1.8

Kerogen density, g/cm3

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 49

1

1.5

Ro, %

2

2.5

(b)

1.6 1.4 1.2 1.0 1.2

1.0

0.8

Atomic H/C

0.6

0.4

Fig. 9 Single-column width ACS Paragon Plus Environment

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Specific surface area, m2/g

Page 45 of 49

Valenza II et al. 2013 Cheshire et al. 2017

600

400

200

0

0.5

1

1.5

2

2.5

Ro, %

Fig. 10 Single-column width ACS Paragon Plus Environment

3

Energy & Fuels

5

S, wt%

4 3 2 1

Atomic S/C x 100

0 2.0 1.5 1.0 0.5 0 0.25

Saliphatic/Saromatic

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 49

0.2 0.15 0.1 0.05 0

0

1

2

Ro, %

3

1.2

1.0

0.8

0.6 0.4

Atomic H/C

Fig. 11 Single-column width ACS Paragon Plus Environment

Page 47 of 49 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Table 1. Woodford shale samples, ordered by increasing vitrinite reflectance VRo OPL # County Type Depth TOC ft % N wt% 549G Garvin Cuttings 8,530-8,540 0.55 41 5.8 536G Carter Cuttings 7,690-7,700 0.57 55 6.0 576D Garvin Cuttings 10,970-10,990 0.63 44 4.3 787G McClain Cuttings 9,960-9,980 0.64 50 3.9 758M Garvin Cuttings 8,780-8,800 0.66 40 5.9 165 Blaine Core 8,414-8,437 0.71 145 6.6 150 Caddo Core 6,129-6,236 0.77 25 5.8 243M Grady Cuttings 12,120-12,130 0.83 40 2.2 569F Grady Cuttings 12,530-12,560 0.90 52 4.1 270H Blaine Cuttings 9,620-9,630 0.94 40 2.8 570E Stephens Cuttings 14,670-14,680 0.94 54 5.6 436U Stephens Cuttings 14,280-14,290 1.07 40 5.1 241 Grady Cuttings 12,960-12,970 1.16 65 8.3 241G Grady Cuttings 12,940-12,950 1.16 65 5.3 246D Canadian Cuttings 12,401-12,408 1.32 45 0.4 353H Blaine Cuttings 13,730-13,740 1.55 40 4.3 803M Beckham Cuttings 16,440-16,450 1.82 42 5.2 271D Roger Mills Cuttings 16,250-16,260 2.48 42 0.8 581E Roger Mills Cuttings 18,150-18,150 2.59 69 2.6 392E Beckham Cuttings 23,920-23,930 4.05 70 2.7

S1 mgHC/g 0.29 0.34 0.51 0.26 0.55 0.05 0.15 0.10 0.28 0.15 0.20 0.21 0.09 0.18 0.05 0.10 0.07 0.10 0.08 0.08

S2 mgHC/g 22.7 24.3 9.1 8.5 18.0 13.1 19.9 1.7 3.4 2.4 2.9 1.5 4.6 1.9 0.07 0.56 0.39 0.10 0.05 0.50

ACS Paragon Plus Environment

S3 mgHC/g 0.24 0.23 0.43 0.27 0.35 0.21 0.16 0.28 0.47 0.62 0.65 0.32 0.69 0.41 0.44 0.33 0.82 0.49 0.33 0.34

T max Calc. Ro H.I. O.I. °C % mgHC/g mgHC/g 434 0.65 392 4 440 0.76 402 4 440 0.76 213 10 437 0.71 215 7 435 0.67 308 6 444 0.83 200 3 441 0.78 343 3 452 0.98 78 13 446 0.87 82 11 438 0.72 83 22 442 0.80 52 12 465 1.21 29 6 462 1.16 55 8 452 0.98 36 8 474 1.37 17 107 456 1.05 13 8 444 0.83 7 16 436 0.69 13 65 414 0.29 2 13 397 0.00 19 13

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 2. Definition of IR structural indices extracted from transmission FTIR spectra Parameter Description Formula Aliphatic CH3/CH2 ratio I 2962/I 2925 CH3/CH2 Aromaticity Aromatic CH/[Aromatic CH + Aliphatic CHx] I 3000-3100/[I 3000-3100 + I 2800-3000] A-factor Aliphatic CH2/[Aliphatic CH2 + Aromatic C=C] [I 2857 + I 2925]/[I 2857 + I 2925 + I 1600-1630] C-factor C=O/[C=O + Aromatic C=C] I 1650-1770/[I 1650-1770 + I 1600-1630]

ACS Paragon Plus Environment

Page 48 of 49

S aliphatic mol% 14 13 8 8 n.d. 8 8

S aromatic mol% 83 83 89 89 n.d. 88 83

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 3. Kerogen density, surface area, elemental analysis, FTIR, and XANES results Skeletal Surface Elemental Analysis Calc. IR Structural Indices XANES Ro Density Area C H N S O CAromatic CH3/CH2 Aromaticity A-Factor C-Factor Elemental Sulfide Thiophene Sulfoxide Sulfone Sulfonate Sulfate g/cm3 m2/g % wt% wt% wt% wt% wt% % mol% mol% mol% mol% mol% mol% mol% 0.55 1.203 54 80.2 7.3 2.0 3.5 5.9 46 0.59 0.01 0.74 0.45 7 6 78 3 0 3 2 0.66 1.228 57 80.0 6.9 2.4 3.1 5.4 49 n.d. n.d. n.d. n.d. 7 5 75 4 0 4 5 0.90 150 81.2 5.4 2.6 1.8 6.8 64 0.81 0.06 0.46 0.44 7 0 78 3 0 6 6 1.16 1.356 114 82.0 5.2 2.6 1.6 4.3 66 0.99 0.10 0.31 0.41 7 0 76 3 0 6 9 1.55 1.385 235 79.7 4.6 2.3 1.3 5.5 69 n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. 1.82 1.380 248 80.1 3.9 2.5 2.1 8.1 76 1.09 0.20 0.22 0.39 7 0 76 4 0 6 8 2.59 1.548 298 81.0 3.1 2.0 1.6 8.7 83 1.11 0.37 0.10 0.28 7 0 70 6 1 10 6 n.d. = no data acquired due to limited sample

Page 49 of 49 Energy & Fuels

ACS Paragon Plus Environment