Subscriber access provided by Kaohsiung Medical University
Fossil Fuels
Thermodynamic Characteristic of Methane Sorption on Shales from Oil, Gas and Condensate windows Feng Yang, Boyu Hu, Shang Xu, Qingbang Meng, and Bernhard M. Krooss Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.8b02140 • Publication Date (Web): 31 Aug 2018 Downloaded from http://pubs.acs.org on September 2, 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.
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 32 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
Thermodynamic Characteristic of Methane Sorption on Shales from Oil, Gas and Condensate windows
Feng Yang,a Boyu Hu,a Shang Xu,a,* Qingbang Meng,a Bernhard M. Krooss b
a
Key Laboratory of Tectonics and Petroleum Resources (China University of Geosciences),Ministry of
Education, Wuhan 430074, PR China b
Energy and Mineral Resources Group (EMR), Institute of Geology and Geochemistry of Petroleum and
Coal, Lochnerstr. 4-20, RWTH Aachen University, 52056 Aachen, Germany * Corresponding Author:
[email protected] 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
1
ABSTRACT: High-pressure methane sorption isotherms at 40–101 °C and pressures up to 25 MPa were
2
measured on Jurassic lacustrine and Silurian marine shales from China. Shale samples span a thermal
3
maturity range from low mature (oil window) to overmature (dry gas window). Low pressure CO2 and
4
N2 adsorption techniques were used to quantify specific surface area, pore volume, and pore size
5
distributions. Thermodynamic characteristic of methane sorption on shales was assessed based on the
6
experimental multi-temperature isotherms. The effects of physical and chemical properties of shales on
7
thermodynamic properties were analyzed and discussed. Finally, standard enthalpy of sorption was first
8
introduced to evaluate the sorption affinity of methane on shales, and a general pattern describing the
9
evolution of methane sorption as a function of thermal maturity was proposed.
10
The Langmuir sorption capacity of these shales varies from 0.09 to 0.16 mmol/g. The low-TOC,
11
clay-rich lacustrine shales have comparable methane sorption capacities than organic-rich, high
12
thermal-maturity marine shales, though high thermal-maturity shale samples tend to have larger
13
micropore volume than low mature shales. Clay minerals, especially I/S mixed layer minerals, contribute
14
a lot to methane sorption of lacustrine shales in oil window. The isosteric heat of methane sorption on
15
these shales decreases with increasing absorbed amount. The commonly used Clausius-Clapeyron
16
equation, which neglects the real gas behavior and adsorbed volume, tends to overvalue the isosteric heat.
17
The standard enthalpy of sorption reflects the comprehensive effect of physical and chemical properties
18
of shales on gas sorption, and shows a parabolic-like pattern with thermal maturity. The standard enthalpy
19
of sorption first decreases with increasing thermal maturity up to the condensate window shales (late
20
mature, equilibrated vitrinite reflectance 1.0–1.1%) and subsequently increases toward the overmature
21
shales.
22
KEYWORDS: shale; sorption isotherm; thermodynamics; thermal maturation; heat of sorption
2
ACS Paragon Plus Environment
Page 2 of 32
Page 3 of 32 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
23
1 Introduction
24
The commercial success of recovering gas from unconventional reservoirs such as organic-rich shales has
25
led to increased efforts in studying the sorption behavior of these natural porous materials [1–6]. Understanding
26
gas sorption behavior in shales is crucial to an accurate assessment of gas-in-place (GIP) of shale gas systems. It
27
is roughly estimated that 50−80% of the total amount of gas in several shale reservoirs in United States is
28
considered to be stored as sorbed gas [7]. However, gas sorption on shales is complicated since shale is a kind of
29
natural materials with strong heterogeneity. Shales are made up of organic matter, clay minerals, and other
30
inorganic matters. At the same time, pore size distributions of shales extend over a wide range from nanometers
31
to micrometers [8]. Complicated mineral compositions and variable pore sizes offer a great challenge to
32
accurately predict methane sorption capacity in shales.
33
Gas sorption on shales is a function of type and abundance of organic and inorganic matter, thermal
34
maturation, and other external conditions (pressure, temperature, moisture content etc.) [9–15]. Hydrocarbon
35
generation and storage are closely associated with the thermal maturation of organic matter in shales. Though
36
significant progress has been made in analyzing the factors affecting methane sorption of shales, evolution of
37
sorption capacity with thermal maturation is not fully understood [9–15]. Part of this is due to that the total
38
organic matter carbon (TOC) content typically trumps the effect of other geological controls [16]. Previous
39
studies on shales from North America and Europe have demonstrated the critical role of organic matter on
40
methane sorption [8–15]. In order to extract the effect of thermal maturity, sorption data are commonly reported
41
normalized to TOC content [10–17]. However, it is demonstrated that this approach can bring about significant
42
errors due to the uncertainties in measured TOC values [18]. Furthermore, thermal evolution of organic matter is
43
accompanied by changes in both physical properties (pore structure) and chemical compositions (function
44
groups). Though high-resolution scanning electron microscopy (SEM) techniques have revealed various
45
nano-scaled pores in organic matter [19–23], there seems to be no specific relationship between thermal maturity 3
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
46
and porosity in organic matter considering the abundance and type of organic matter [22,23]. Several
47
investigators found that even in the same shale sample, the organic porosity of different organic matter particles
48
might be totally different [23]. What’s more, investigations on the maceral chemistry of coals have shown that
49
low mature organic matter is rich in aliphatic and oxygen containing functional groups [24–26]. In contrast, high
50
rank coals (anthracite) show the accumulation of condensed aromatic structures [24–26]. Variations in both pore
51
structure and function groups of organic matter share a major challenge in fundamentally understanding the
52
relationship between thermal maturity and methane sorption capacity.
53
In addition to organic matter, clay minerals potentially help to gas sorption on shales. Literatures have
54
documented the influences of clay minerals on sorption capacity of shales [27–29]. The contribution from clay
55
minerals complicates the evolution of sorption behavior of shales with thermal maturity. Compared to marine
56
shales, lacustrine shales are often clay rich and low in thermal maturation [30–33]. Unlike organic-rich shales, it
57
is difficult to fully explain the variations in sorption capacity via TOC content alone in low-TOC, clay-rich
58
shales [32]. It is important to compare gas sorption capacities of lacustrine and marine shales, since GIP from
59
lacustrine shales accounts for roughly one third of total recoverable shale gas resource in China [34].
60
Sorption is a process that adsorbate molecules accumulate on the adsorbent surface to reduce the potential
61
energy. Experimental measured sorption data not only provide the absorbed amount, but also imply sorption heat
62
in the sorption process [35, 36]. Though significant progress has been made on the evaluation of sorption
63
capacity, there is relatively few detailed studies on the thermodynamic aspect of gas sorption on shales [6,10,13].
64
The isosteric heats of bulk shales are reported to fall between the main organic and inorganic constituents of
65
shales [10]. The isosteric heat of methane sorption on type III kerogen is thought to be higher than that of type I
66
and type II kerogens because of more aromatic moieties [13]. To improve understanding of the interaction
67
between shale and methane molecules, a thorough investigation on the thermodynamic process of methane
68
sorption in these complex micro- and mesoporous materials is needed. 4
ACS Paragon Plus Environment
Page 4 of 32
Page 5 of 32 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
69
This paper aims at shedding light on methane sorption on shales from the perspective of sorption
70
thermodynamics. High-pressure sorption isotherms for methane on lacustrine and marine shales spanning
71
different thermal maturity stages were measured up to pressures of 25 MPa, at temperatures from 40 to 101 °C
72
using a manometric method. Thermodynamic parameters of methane sorption on shales were estimated
73
according to the experimental multi-temperature sorption isotherms. The effects of physical and chemical
74
properties of shales on the thermodynamic properties have been discussed. Finally, a general evolution pattern
75
relating sorption heat to thermal maturity was proposed by the combination of our isotherms and those from
76
North America and Europe.
77
2 Thermodynamic frameworks
78
2.1 Excess sorption and absolute sorption
79
The measured quantity in all sorption experiments (manometric, volumetric, gravimetric, chromatographic,
80
etc) is “Gibbsian surface excess adsorption” [37]. With the Gibbs definition of adsorption, the experimentally
81
measured excess sorption can be written as:
82
= − ∙ = 1 −
(1)
83
Here nexc (mmol/g) is the Gibbs excess sorbed amount; nabs (mmol/g) is the absolute sorption; Vabs (g/mmol) is
84
the volume of the adsorbed phase; ρgas and ρads are densities of the free phase and adsorbed phase, respectively,
85
both expressed in the same units (e.g., kg/m3).
86
Rigorous thermodynamic analysis relies on absolute sorption values. Unfortunately, the density and volume
87
of the adsorbed phase cannot be experimentally measured. It is common to assume that the excess sorption
88
approximates the absolute sorption in literatures [13,14,27]. Although this assumption is reasonable at low
89
pressure conditions, it cannot be applicable in high-pressure studies. Under low pressure conditions, the volume
90
occupied by the adsorbed phase is relatively small, thus the Gibbs excess sorption approximates the absolute
91
sorption. However, the adsorbed volume is no longer negligible at high pressures. To determine the absolute 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
Page 6 of 32
92
amount of sorption and avoid the well-documented errors associated with equating excess sorption and absolute
93
sorption, an adapted Langmuir function, which explicitly considering the volume of the adsorbed phase, was
94
applied to fit the measured excess sorption data [10, 16, 38]:
= 1 −
95 96
= 1 −
(2)
Where =
97
√
%
exp ( ) &
(3)
98
Here nL (mmol/g) is the amount of substance adsorbed at infinite pressure, also known as the Langmuir
99
volume; K (1/MPa) is the Langmuir constant; A (K1/2MPa-1) is a prefactor and E (kJ/mol) is a binding energy; R
100
(kJ/(mol·K)) is the universal gas constant. Determining the absolute sorption quantity (nabs) is simple when the
101
adsorbed gas density or the volume of the adsorption layer (Vads) is known. Since it is difficult to experimentally
102
measure the adsorbed phase density, a number of assumptions have been presented to estimate it [39]. One of the
103
general approaches is letting the adsorbed phase density to be an independent parameter in the fitted equation of
104
choice. The adapted Langmuir equation has been shown to be suitable fitting equations for determining absolute
105
sorption from excess isotherms since they are monotonically increasing and contain a relatively small number of
106
fitting parameters to achieve a satisfactory fit to the sorption data [10–12]. The variables that remain unknown in
107
Eq. (2) and Eq. (3) are the Langmuir volume nL, prefactor A, binding energy Ei, and adsorbed phase density ρads.
108
These variables can be obtained by simultaneously performing multivariable nonlinear regression over the
109
sorption data recorded at different temperatures.
110
2.2 Isosteric enthalpy of sorption
111
The isosteric enthalpy (heat) of sorption is one of the key thermodynamic variables for gas sorption systems.
112
The isosteric heat is always captured by differentiating a series of sorption isotherms at constant loading using
113
the popular Clausius-Clapeyron equation [40–42]. However, strictly speaking, the approach based on the
114
Clausius-Clapeyron equation only applies to perfect gas. It has been noted that the extension of this method to 6
ACS Paragon Plus Environment
Page 7 of 32 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
115
the case of real gas will lead to considerable confusion [40, 43]. In this paper, we start from the basic
116
thermodynamic theory to derive the general expressions of sorption thermodynamics, and then compare the
117
results of the general and simplified equations.
118 119 120 121 122 123
From the theory of thermodynamic equilibrium, the change in entropy of sorption (∆) ) can be calculated as a function of sorption uptake,
*+
, -,
=.
-.
=.
∆,
(4)
-.
Furthermore, the Gibbs free energy of the adsorbed phase (/ ) is equal to the Gibbs free energy of the gas phase (/ ) at thermodynamic equilibrium, in which case, the isosteric entropy of adsorption (∆) ) is ∆/ = / − / = ∆0 − 1∆) = 0
(5)
∆0 = 1∆)
124
(6)
125
The isosteric enthalpy of adsorption (∆0 ) is a widely used quantity to characterize the strength of
126
binding interactions at a fixed temperature, pressure, and coverage. Typically it is determined using the isosteric
127
method and reported as a positive value (34 ), the so-called isosteric heat in literatures [40–42],
128
34 = −∆0 = −1∆) = −T
*+
6 − 7 = −T
*+
8
9
−
9
:
(7)
129
Eq. (7) does not include any assumption and is arguably the most important equation of sorption
130
thermodynamics. A number of assumptions are proposed to simplify this expression. The most common one is
131
assuming that the volume of adsorbed phase is negligible compared to the bulk phase, which yields:
132 133 134 135
34 = −∆0 = −T
*+
If the bulk phase is assumed as ideal gas, i.e.,
6− 7 = −T 9
=
&
*+
8−
9
:
(8)
, Eq. (7) can be further simplified as the commonly
used Clausius-Clapeyron equation (where it also neglects the volume of adsorbed phase):
34 = −∆0 = −T
*+
−
&
=
& ;
*+
(9)
136
Optionally, the derivative of pressure with respect to temperature can be determined if the functional
137
relation between the amount of substance adsorbed and pressure and temperature are provided. In this article, the 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
138
139
140 141 142 143 144 145 146 147
Page 8 of 32
derivative is further treated based the Langmuir equation, and can be decomposed as follows:
*+
B A + + >? >C A +
(10)
In the case of Langmuir equation, the respective components of the derivative in Eq. (10) can be written as: