Accepted Manuscript Pore structure, gas storage and matrix transport characteristics of lacustrine Newark shale R. Fink, A. Amann-Hildenbrand, P. Bertier, R. Littke PII:
S0264-8172(18)30278-2
DOI:
10.1016/j.marpetgeo.2018.06.035
Reference:
JMPG 3404
To appear in:
Marine and Petroleum Geology
Received Date: 13 December 2017 Revised Date:
29 June 2018
Accepted Date: 29 June 2018
Please cite this article as: Fink, R., Amann-Hildenbrand, A., Bertier, P., Littke, R., Pore structure, gas storage and matrix transport characteristics of lacustrine Newark shale, Marine and Petroleum Geology (2018), doi: 10.1016/j.marpetgeo.2018.06.035. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
ACCEPTED MANUSCRIPT 1
Pore structure, gas storage and matrix transport characteristics
2
of lacustrine Newark Shale
3
R. Finka,*, A. Amann-Hildenbranda, P. Bertierb, R. Littkea
4 a
Institute of Geology and Geochemistry of Petroleum and Coal, Energy and Mineral Resources
a
RI PT
Group (EMR), RWTH Aachen University, Lochnerstr. 4-20, D-52056 Aachen, Germany
Department of Clay & Interface Mineralogy, Energy and Mineral Resources Group (EMR), RWTH
Aachen University, Bunsenstr. 8, D-52072 Aachen, Germany *
[email protected]
SC
5 6 7 8 9 10 11
Abstract
13
Shale gas production in the U.S. fuelled research activities in unconventional
14
reservoir rock characterisation. Most studies focused on organic-rich shales of
15
marine origin, while disregarding lacustrine sequences. In this study, thirteen
16
lacustrine shale samples from the Newark Basin, NJ, USA are comprehensively
17
characterised in terms of pore structure, gas storage and matrix transport
18
characteristics. These thermally overmature (VRr 1.4 to 2.7 %) shales have a Na-
19
rich, heterogeneous mineralogy with TOC contents of up to 3.6 %.
20
Methane sorption capacity and pore structure parameters as identified with low-
21
pressure N2 physisorption (microporosity, BET surface area) are neither
22
interrelated with each other nor with any shale components (e.g. clay content,
23
TOC).
24
In contrast, porosity shows a positive correlation with TOC content, which is also
25
typical
26
Correspondingly, porosity and TOC positively correlate to bedding parallel
27
matrix permeability coefficients (between 2 and 80 nD at 40 MPa effective
28
stress). In contrast, permeability coefficients perpendicular to bedding are two
29
to three orders of magnitude lower.
30
Compared to previous studies on marine lithotypes, Newark shales have rather
31
poor gas storage properties with average He-porosities of 2.3 % and average
32
methane sorption capacities of 0.047 mmol g-1.
33
Keywords: Shale Gas; Permeability; Porosity; N2 Physisorption; CH4 Sorption
AC C
EP
TE D
M AN U
12
for
many
thermally
overmature
1
marine
shale
sequences.
ACCEPTED MANUSCRIPT
1 Introduction Knowledge of pore structure is critical for understanding both, fluid flow and
36
storage mechanisms in shale (Chalmers & Bustin 2012, Clarkson et al. 2013, Slatt
37
& O’Brian 2011). Porosity characteristics of shale are very complicated due to
38
the large variability in pore sizes (micro-, meso- and macroporosity),
39
microstructure, mineralogical composition, organic carbon (TOC) content and
40
the geologic history of each system. As a result, advanced pore structure
41
investigation techniques from other fields were applied to gain new insights into
42
microporous (< 2 nm, see IUPAC classification (Sing et al. 1985)) and
43
mesoporous (2 – 50 nm, IUPAC classification (Sing et al. 1985)) characteristics of
44
shale. These methods can be divided into fluid invasion and radiation methods.
45
Fluid invasion methods include high-pressure mercury intrusion (MIP) and low-
46
pressure gas adsorption (N2 and CO2) (e.g. Clarkson et al. 2012, 2013, Kuila &
47
Prasad 2013, Ross & Bustin 2009). Radiation methods include field emission
48
scanning electron microscopy (FE-SEM) and transmission electron microscopy
49
(TEM) after preparation of high-quality surfaces by argon ion beam milling (e.g.,
50
Klaver et al. 2016, 2015, 2012, Loucks et al. 2009, Milliken et al. 2013 and
51
references therein) and small and ultra-small angle neutron scattering
52
(SANS/USANS) (Bahadur et al. 2015, Clarkson et al. 2013, King et al. 2015,
53
Mastalerz et al. 2012, Sun et al. 2016). Application of these methods on marine
54
gas shales revealed that porosity, pore size distributions, specific surface area
55
and pore morphology may be related to TOC, mineralogy, diagenesis and
56
microstructure. The controls are, however, poorly constrained and strongly
57
differ between geologic formations.
58
Natural gas in shale is stored in both free and adsorbed states. The adsorbed gas
59
capacity is strongly controlled by microporosity due to its large internal surface
60
area. This is why several studies found that sorption capacity is often positively
61
correlated with TOC, which in turn often positively correlates to microporosity
62
(e.g. Chalmers & Bustin 2007, Ross & Bustin 2009). However, also clay minerals
63
may significantly contribute to adsorbed gas storage capacity (Gasparik et al.
64
2012, Ji et al. 2012, Merkel et al. 2016, 2015, Ross & Bustin 2009). In contrast,
65
only few studies report experimental gas permeability data and relate the
AC C
EP
TE D
M AN U
SC
RI PT
34 35
2
ACCEPTED MANUSCRIPT transport properties to pore structure characteristics (Ghanizadeh et al. 2015,
67
Ismael & Zoback 2016).
68
Lacustrine shales have not been in the focus of shale gas research, because the
69
currently developed shale gas systems in the U.S. are exclusively of marine
70
origin. In contrast, as lacustrine shale is widely distributed in the sedimentary
71
basins of China (Ordos, Bohai, Songliao and Sichuan), there is an increasing
72
interest in its transport and storage properties for shale gas and shale oil
73
development (e.g. Fu et al. 2015, Guo et al. 2014, Ji et al. 2014, Liu et al. 2015,
74
Wang et al. 2015, Xiong et al. 2016).
75
Lacustrine sequences may significantly differ from marine sequences and it is
76
therefore important to asses their shale gas potential individually. Among the
77
key differences is their greater sensitivity to high frequency climatic variability
78
leading to rapid lake level fluctuations that result in rapid facies changes.
79
Although this significantly reduces the core unconventional target area, the
80
thickness of lacustrine shales may significantly exceed the thickness of marine
81
shales (Katz & Lin 2014). If favourable reservoir conditions are present, the
82
greater thickness and the lack of carrier systems in many lacustrine system may
83
lead to a lower expulsion efficiency increasing the potential viability of these
84
unconventional plays,
85
2 Newark basin Geology and Stratigraphy
86
The Newark basin is located in Eastern North America and developed in a huge
87
rifting zone (~8000 km long), which resulted in the break-up of Pangea and
88
formation of the central part of the Atlantic Ocean from the Late Permian to the
89
Early Jurassic. Normal fault systems with significant hanging wall subsidence
90
(1000s of meters) lead to the development of individual rift valley basins that
91
were most of the time filled by great lakes, far from access to the sea (Van
92
Houten 1964, Olsen 2010, 1997, Olsen et al. 1996, Schlische 2003).
93
The sediment infilling of the Newark half graben system is of continental, largely
94
lacustrine origin ranging from early Triassic (Carnian) to Early Jurassic
95
(Hettangian) age (Olsen 1986) (Figure 1). The fluvial Stockton formation is
96
followed by highly cyclical lacustrine strata with a thickness of > 3000 m. They
AC C
EP
TE D
M AN U
SC
RI PT
66
3
ACCEPTED MANUSCRIPT are subdivided into the mostly gray and black Lockatong Formation, which is
98
overlain by the mostly red Passaic Formation. Both of them preserve the longest
99
unbroken record of orbital forcing that was studied in detail using the Newark
100
Basin Coring Project cores that were also used for this study (e.g. Olsen & Kent
101
1996, 1990, Olsen et al. 1996) (Figure 1).
102
Cyclical lacustrine strata formed due to lake level fluctuations in a tropical
103
climate controlled by orbital parameters. The sedimentary rocks are subdivided
104
based on lithology and sedimentary structures into “depth ranks” from 0 – 5
105
related to the relative water depth of the lake during deposition. Rank 0 consists
106
of red mudstone produced by playas or dry lakes whereas rank 5 is comprised of
107
microlaminated black organic-rich calcareous mudstones formed in deep lakes
108
with perennial chemical stratification. High depth rank (3 – 5) units have
109
elevated organic carbon contents (> 1 %) and are usually gray or black (Olsen
110
1986, Olsen & Kent 1996, Olsen et al. 1996).
111
3 Samples
112
For this study, we selected lacustrine mudstone samples from the Lockatong and
113
Passaic formations. Overall thirteen locations were sampled, one of them from
114
outcrop (Palmyra, Pennsylvania) and twelve from the Newark Basin Coring
115
Project cores (NBCP) (Figure 1 and Table 1). The NBCP recovered over 6770 m
116
(22200 ft) of continuous core of excellent quality by offset drilling between
117
seven carefully selected sites. Details on the NBCP are published in Olsen et al.
118
(1996) and Olsen & Kent (1996).
119
Sampling was focused on black mudstones to cover a highly variable TOC
120
sequence and on a maximal variation in composite depth of > 2500 m (8600 ft)
121
to cover maturity changes. To analyse the effect of facies changes, we
122
furthermore sampled different facies at location 8, 9 and 10 in a very narrow
123
depth interval of around 3.5 m (11.4 ft) (Figure 1).
124
Samples 1,3, 4, 5, 7 and 11 are finely-microlaminated darkgrey-black mudstones
125
with depth rank of 4-5 deposited during lake high-stands (Table 1). In contrast,
126
samples 6, 12 and 13 with depth rank 3 are thin-bedded dark-grey mudstone
127
deposited in a relatively shallow perennial lake and samples 2, 8 and 10 of depth
AC C
EP
TE D
M AN U
SC
RI PT
97
4
ACCEPTED MANUSCRIPT rank 1-2 are gray massiv mudstones (sample ten with mudcracks) typical for
129
lakes with saline and drying out episodes. Sample 9 is the only red massive
130
mudstone of depth rank 0 which is typical for dry playa mudflat deposition
131
(Olsen et al. 1996, 1986, Olsen & Kent 1996).
132
From each location cylindrical sample plugs oriented perpendicular to bedding
133
and a slab (length: 10 – 30 cm, width: 2 – 3 cm) were cut from the whole core
134
using water as lubricant. Additionally, on three locations (3, 11 and 13) plugs
135
were drilled parallel to bedding.
RI PT
128
SC
136
4 Methods
138
4.1 Sample preparation
139
Cylindrical plugs of 28.5 mm diameter and 15 to 47 mm length were drilled
140
either perpendicular or parallel to bedding. Plugs were dried in a vacuum oven
141
for at least 48 hr at 105 °C until weight constancy and used for transport and
142
porosity experiments.
143
For high-pressure methane sorption experiments, total sulphur and Rock-Eval
144
Pyrolysis, a subsample (20 to 40 g) of the core slab material was ground to a size
145
of < 0.1 mm and dried for at least 48 hr at 105 °C in a vacuum oven.
146
X-ray diffraction (XRD) and low-pressure physisorption experiments were
147
performed on sieved size fractions prepared after manual crushing of another
148
core slab subsample (10 to 20 g) in a mortar.
149
4.2 Petrographic analysis and VRr
150
Polished blocks used for petrographic analysis and vitrinite reflectance
151
measurements (VRr) were produced following the procedure and using
152
instruments described in Sachse et al. (2012).
153
Vitrinite reflectance measurements were performed under oil immersion (ne =
154
1.518) on a Zeiss Axio Imager microscope equipped with a tungsten-halogen
155
lamp, a 50/0.85 Epiplan-NEOFLUAR oil immersion objective and a 546 nm filter.
156
A cubic-zirconia standard with 3.125 % reflectance was used for calibration.
AC C
EP
TE D
M AN U
137
5
ACCEPTED MANUSCRIPT Details of the analytical procedure and instrumentation are described in Littke et
158
al. (2012).
159
4.3 Geochemical and Mineralogical analysis
160
4.3.1 TOC, TIC, TS and Rock-Eval Pyrolysis
161
Total organic carbon (TOC) and total inorganic carbon (TIC) were measured on
162
powdered samples with a liquiTOC II analyser. The instrument operates in a
163
non-isothermal mode and analyses the released CO2 during heating with a non-
164
dispersive infra-red detector (NDIR). TOC and TIC can be analysed during one
165
temperature ramp without previous acidification.
166
Total sulphur (TS) content was measured with a LECO S200 analyser and Rock-
167
Eval pyrolysis measurements were performed on a Rock-Eval VI instrument
168
following the procedure initially described in Espitalié (1977).
169
4.3.2 Quantitative X-ray diffraction analysis
170
Bulk mineralogical compositions were derived from diffraction patterns of
171
randomly oriented powder preparates. Core slab material was crushed manually
172
in a mortar to avoid strain damage. To ensure uniform crystallite sizes, the 400 –
173
800 µm size fraction with an internal standard (Corundum, 20 wt %) was milled
174
in ethanol for 15 min with a McCrone Micronising mill. The standard was added
175
to evaluate the accuracy of the analysis.
176
The measurement was done on a BrukerAXS D8 Advance diffractometer using a
177
CuKα-radiation produced at 40 kV and 40 mA. Diffractograms were recorded
178
from 2° to 92° 2θ in 0.02° steps.
179
Quantitative phase analysis was performed by Rietveld refinement usig the
180
BSGM software with customized clay mineral structure models (Ufer et al. 2008).
181
Mineralogical compositions were related to the crystallite content of the
182
analysed samples.
183
4.4 Porosity and pore size distribution measurements
184
4.4.1 Low-pressure N2 physisorption analysis
185
Measurements were performed using a Micrometrics Gemini VII 2390t apparatus
186
on the 63 – 200 µm grain size fraction. Approximately 1 g of sample material was
AC C
EP
TE D
M AN U
SC
RI PT
157
6
ACCEPTED MANUSCRIPT degassed at a vacuum below 2 Pa in a Micromeritics VacPrep 061 at 110 °C for at
188
least 24h.
189
For nitrogen physisorption analysis, adsorption and desorption measurements
190
were performed in a nitrogen bath (77.3 K) at 90 relative pressure steps (p/p0)
191
between