High-Pressure Methane Sorption on Dry and Moisture-Equilibrated

Dec 13, 2016 - Energy and Mineral Resources Group (EMR), Institute of Geology and Geochemistry of Petroleum and Coal, RWTH Aachen University, ...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Newcastle, Australia

Article

High-Pressure Methane Sorption on Dry and Moisture-Equilibrated Shales Feng Yang, Congjiao Xie, Zhengfu Ning, and Bernhard M. Krooss Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.6b02999 • Publication Date (Web): 13 Dec 2016 Downloaded from http://pubs.acs.org on December 18, 2016

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 free 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 accessible to all readers and 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 28

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

High-Pressure Methane Sorption on Dry and Moisture-Equilibrated Shales Feng Yang a,b,*, Congjiao Xie a, Zhengfu Ning b, Bernhard M. Krooss c

a

Key Laboratory of Tectonics and Petroleum Resources (China University of Geosciences), Ministry of

Education, Wuhan 430074, PR China b

State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum (Beijing),

#18, Fuxue Rd, Changping, Beijing 102249, PR China c

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 were collected on selected Paleozoic shales

2

from Sichuan Basin. Excess sorption measurements were performed on shales with varied water content

3

(dry, moisture equilibrated at 33, 53, 75, and 97% relative humidities) at 39 °C and up to 25 MPa. Water

4

uptake isotherms were collected at 24 °C and parameterized by the Guggenheim-Anderson-de Boer

5

(GAB) model. The effect of organic richness, mineral compositions, and pore structure characteristics on

6

water uptake and methane sorption behavior has been investigated. The mechanism responsible for the

7

decrease in methane sorption capacity of moisture-equilibrated shales is discussed.

8

Water uptake of shales is primarily controlled by clay minerals, and shows a positive correlation with

9

clay mineral content. Water sorption isotherms of shales can be approximately expressed as the sum of

10

the isotherms of individual clay minerals on a mass-fraction base. Methane sorption capacity of these

11

shales is controlled by TOC content. The maximum Langmuir sorption capacity of shales under both dry

12

and 97% RH conditions correlates positively with TOC content. Compared to dry condition, methane

13

sorption capacity of shales moisture-equilibrated at 97% RH is reduced by 44 to 63%. The experimental

14

results indicate a step-wise decline in methane sorption with increasing water content. Evolution of

15

sorption capacity as a function of water content can be divided into three stages: (1) initial decline stage:

16

the decrease of methane sorption capacity is mainly due to competitive sorption of methane and water

17

on hydrophilic clay minerals. (2) steep decline stage: clusters of water molecules block pore space and

18

reduce the sorption capacity significantly. (3) slow decline stage: a contiguous water phase successively

19

fills the macropores and slightly reduces methane sorption by volume displacement.

20

KEYWORDS: shale gas; methane; sorption isotherm; water sorption; Langmuir

2

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28

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

21

1. Introduction

22

Advances in horizontal drilling and hydraulic fracturing have permitted to extract hydrocarbon from

23

unconventional reservoirs. Shale gas and tight oil formations provided 50% of the US dry gas production

24

in 2015, and the percentage is estimated to increase to 69% by 2040.1 Sorption is an important storage

25

mechanism of natural gas in organic-rich shale reservoirs. More than 50% of the total gas capacity of the

26

Devonian shales in the Appalachian basin2 is associated with gas sorption, and the ratio of sorbed gas to

27

total stored gas is up to 70% in biogenic Antrim shales 3.

28

Gas storage in unconventional shale reservoirs is a complex multi-parameter problem. Studies to

29

date identify organic matter characteristics (total organic carbon (TOC) content, thermal maturation, and

30

kerogen type), mineralogy, pore system, moisture, pressure and temperature as important parameters in

31

assessment of gas sorption capacity of shales.4–17 Organic matter in shales has been considered as the first

32

contributor to gas sorption capacity of shales. Previous studies have reported the positive correlation

33

between TOC content and methane sorption on black shales of North America7 and Europe9,15. Thermal

34

maturity of organic matter has positive effect on the methane sorption capacity of shales, and this is

35

ascribed to the increase in the organic micropores and/or change in chemistry characteristics of organic

36

matter during thermal maturation.6,7 When it comes to the type of organic matter of shales, sorption

37

capacity on TOC basis was reported to increase in order: type I < type II < type III, and this was attributed

38

to the larger micropore surface area of humic kerogen compared to other maceral composition when the

39

organic matter is highly mature and the increasing kerogen aromaticity of organic matter in the

40

progression from type I to type III kerogen.5–7 In addition to organic matter, pores of inorganic

41

constituents (clay minerals) accommodate additional sorbed gas due to their developed internal surface

42

area, and contribute appreciable amounts of gas sorption in clay-rich shales under dry condition. 8, 15, 18, 19

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

43

Moisture is known to have a strong reducing effect on methane sorption in coals and shales. Though

44

researchers realize the negative influence of water on sorption behavior of shales, many of them do not

45

address the issue of water in much detail.9, 12, 13 Day et al.20 described a linear decrease in sorption

46

capacity of coals up to certain critical value above which gas sorption capacity is unaffected. However,

47

the interaction between natural carbonaceous materials (such as coal and shales) and water is complex

48

and the mechanisms of sorption capacity reduction with moisture are not fully understood. Sorption of

49

water on coal is connected with physisorption on functional groups and chemisorption onto mineral

50

surfaces.21 Typically, the chemisorbed water cannot be desorbed completely under vacuum at low

51

temperatures. In addition, water tends to form large clusters via hydrogen bonds and block large pores and

52

interparticle voids.20, 21

53

Although significant progress has been achieved in the experimental methodology and technique for

54

measuring gas sorption isotherms, data on the high-pressure sorption isotherms on shales in the presence

55

of moisture/water are scarce, and the role of water on gas sorption of shales is still poorly constrained.

56

The present paper summarizes the experimental results of recent investigations on methane sorption

57

capacity of dry and moisture-equilibrated marine shales from Sichuan Basin that are being considered as

58

the main target for shale gas exploration in China. High-pressure methane sorption isotherms at 39 °C and

59

pressures up to 25 MPa were determined for four Paleozoic shales at five different relative humidity

60

conditions. The main objectives of the study were to assess the sorption capacity of these shales, with

61

special reference to the mechanism of the reduction of their methane sorption capacity caused by the

62

presence of water. Besides, the effects of organic richness, mineral compositions, and pore structure

63

characteristics on water uptake of shales are also discussed.

64

2. Experimental section 4

ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28

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

65

2.1 Samples

66

Marine shales from the Upper Ordovician Wufeng Formation and the Lower Silurian Longmaxi

67

Formation in Sichuan Basin were collected for this experimental study (Figure 1). The Wufeng-Longmaxi

68

shales in Sichuan Basin are currently considered as the main targets for shale gas exploration in China.

69

Several shale gas fields operated by PetroChina and Sinopec have been developed in Sichuan Basin. One

70

upper Ordovician (CN_11) and two lower Silurian shale samples (the bottom of Longmaxi Formation:

71

CN_22, the top of Longmaxi Formation: CN_33) were collected from the field standard stratotype section

72

in Changning area of Southern Sichuan Basin. The Wufeng Formation (CN_11) and the Lower part of

73

Longmaxi Formation (CN_22) are mainly black graptolite shales, while the Upper part of Longmaxi

74

Formation is mainly silty shale (CN_33). TOC content is high at the bottom of the stratotype section and

75

then decreases upwards, which indicates that the sea level decreased, and the hydrodynamic force

76

increased from the Late Ordovician to Late Silurian.22 The fall of sea level was unfavorable for the

77

preservation of organic matter. In addition, one Lower Silurian Longmaxi Formation shale sample

78

(CQ_14) from Chongqing area was also added to the sample list. General geological descriptions of the

79

sedimentary environment, petrofacies, and potential evaluations of the shales in Sichuan Basin have

80

recently been published.15,22–24 Results from organic geochemistry and XRD measurements of the samples

81

are presented in Table 1. The samples investigated exhibit TOC values of 0.96–9.40%, with helium

82

porosity values ranging between 4.4 and 13% (Table 2). The “equivalent” vitrinite reflectance (Req),

83

derived from the bitumen reflectance, ranges from 2.4% to 2.8%, indicating that these samples having

84

reached a very high thermal maturity.

85

2.2 Pore size distribution and specific surface area

86

The application of N2 physisorption analysis, originally developed for material sciences,25 has been

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 28

87

addressed in shale characterization previously.26,27 Low pressure nitrogen adsorption-desorption

88

measurements at 77 K (–196 °C) were conducted to quantitatively characterize the pore structure of

89

shale samples. Sorption measurements on shales are commonly performed with powdered samples (60–

90

100 mesh particle size). 4,6,7–9,15,18 In this study, powder samples (~100 mesh, 0.149 mm) were outgassed

91

overnight at 105 °C under vacuum. The adsorption-desorption isotherms were collected on Micromeritics

92

Surface Area and Pore Size Analyzer (Gemini VII 2390t). Adsorbed amounts of N2 were measured for

93

each relative pressure (absolute/saturation pressure) point, and dynamic equilibrium was assumed when

94

the pressure change over an interval (10 s) was < 0.01% of the average pressure in the latter interval.

95

Multi-point Brunauer-Emmett-Teller (BET) analysis was used to quantify the specific surface area, and

96

the Barrett-Joyner-Halenda (BJH) approach was applied to calculate pore size distributions.

97

2.3 Water uptake measurements

98

Water uptake as a function of water vapor pressure was determined by moisture-equilibrating the

99

samples at different relative humidities (RH) at room temperature (~24 °C). Powdered samples (~100

100

mesh, 0.149 mm) were first dried in a 105 °C vacuum oven overnight before determining their dry weight

101

(mdry). Then the dried samples were placed in a desiccator where containing a saturated salt solution to

102

establish a defined water vapour partial pressure (relative humidity, RH). Five different saturated salts

103

solutions were used for humidity control: LiCl (11% RH), MgCl2 · 6H2O (33% RH), Mg(NO3)2 · 6H2O

104

(53% RH), NaCl (75% RH), and K2SO4 (97% RH).28 Moisture equilibration on shale samples was carried

105

out at different humidity levels until the sample weights remained constant. The moisture uptake (in g/g)

106

was calculated as:

107 108

(/ ) =

 (/ ) 

(Eq. 1)

Water sorption isotherms were obtained by plotting the water uptake against the relative pressure

6

ACS Paragon Plus Environment

Page 7 of 28

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

109

(p/p0). The water uptake isotherms were parameterized by the 3-parameter Guggenheim-Anderson-de

110

Boer (GAB) model: 29

111

(/ ) =

  ( )(() )

(Eq. 2)

112

Here wm is the GAB monolayer capacity; aw is the sorbate (water) activity and equivalent to relative

113

pressure (aw = p/p0); c and k are GAB constants, which are correlated with the heat of adsorption. The

114

GAB equation has a similar form as the well-known BET equation, and reduces to the BET equation

115

when k = 1. The unknown parameters in the GAB equation can be determined by fitting the measured

116

data using a non-linear regression.

117

2.4 High-pressure methane sorption experiments

118

High-pressure sorption experiments were conducted on a manometric apparatus at 312 K (39 °C)

119

and up to 25 MPa. The experimental set-up consists of a stainless steel sample cell (SC), two

120

high-pressure electric-pneumatic valves and a high-precision Keller pressure transducer (PT maximum

121

pressure 30 MPa with an accuracy of 0.05% at full scale value). The reference cell (RC) is the volume of

122

the tubing connecting the two valves and the dead volume of the PT. Volumes of the RC and SC were

123

calibrated by helium expansion using stainless steel cylinders of known dimensions placed in the SC.

124

Dry/moisture-equilibrated powder samples (~100 mesh) were placed into the SC and the void volume

125

(Vvoid) of the SC was determined by helium expansion. The measuring procedure starts by charging RC

126

with a certain amount of methane and allowing certain time for equilibration in RC. Then methane

127

expands to the SC and sorption begins. The sorption system reaches equilibrium in 1–3 h, depending on

128

the moisture condition of the shale sample. Furthermore, blank methane expansion tests were performed

129

on steel cylinders of different sizes placed in the SC at target temperature. The blank sorption isotherms

130

of steel cylinders with different “void volumes” were interpolated to acquire the blank expansion value of

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

Page 8 of 28

131

a shale sample at an equivalent void volume. From the “raw” expansion data measured on a shale sample,

132

the blank expansion value was subtracted to obtain the final corrected excess sorption isotherm. A detailed

133

description about the experimental setup and measuring method has been published previously.8, 9, 30, 31

134 135

The excess sorbed mass (mexc), also denoted “Gibbs surface excess”, is calculated according to the following mass balance:

136

 = 

!

− #$% (&, ())*

+,

(Eq. 3)

137

Here Vvoid is the pore space that is not occupied by the powered sample in SC and determined by helium

138

expansion; mtotal is the total mass of adsorbate (methane) transferred into the SC; ρgas(T, P) is free gas

139

density determined from equations of state provided by Kunz et al.32 Excess adsorbed mass was

140

converted to excess adsorbed amount of substance (nexc, mmol/g). The experimental excess sorption data

141

was parameterized by an adapted Langmuir function:8,9,15, 30,33 345

/

142

- = -. //

0 (1)

(1 − 3

)

(Eq. 4)

5

143

Where nL (mmol/g) is the Langmuir sorption capacity (also denoted as the Langmuir volume when

144

expressed in units of volume per mass); P (MPa) is the pressure of the free gas; PL (MPa) is the Langmuir

145

pressure; ρads (kg/m3) is the density of the adsorbed phase.

146

3. Results

147

3.1 Pore structure characteristics by nitrogen adsorption

148

All the low-pressure nitrogen adsorption isotherms of the shales studied are of type IV (Figure 2).

149

The low-pressure nitrogen adsorption isotherms of these shale samples exhibit hysteresis in terms of the

150

divergence of the adsorption-desorption branches, which indicates that mesopores (2 nm < pore widths
33% for CN_22 and CQ_14. The k values obtained from the fitting procedure range

193

from 0.7 to 0.9, which is in accordance with the limit values (0.24< K 53%, and shows remarkable water uptake between 75% and 97% RH. The

246

intersection of water uptake isotherms at high relative pressure is probably caused by larger

247

micro-porosity (Figure 2b), hence greater amount of capillary condensate water. PSD derived from N2

248

adsorption show that the micropore volume of sample CQ_14 is significantly larger than those for other

249

three samples (Figure 2b). Micropores are apt to be filled by water molecules. Water molecules firstly

250

sorb on primary sorption centers of clay minerals (functional groups containing oxygen). Table 4 shows

251

that shale samples with higher clay mineral contents have higher relative loss in methane sorption

252

capacity of shales at the initial stage of water uptake. Monolayer coverage finish when the polar sites are

253

occupied by water molecules. Then the adsorbed water molecules act as secondary sorption sites

254

(nucleation sites) for further water sorption by hydrogen bonding, which results in the formation of

255

three-dimensional water clusters. Water sorption will be remarkably enhanced when the clusters of water

256

molecules connect together by “bridges”.40 Eventually, pore filling will occur at high relative pressure.

257

There is a positive relationship between BET surface area and water content of the 97% RH

258

moisture-equilibrated samples (Figure 6c). The intercept value is about 0.01 g/g, indirectly indicating that

259

water in shales is likely to be in the form of condensed phase apart from sorbed state.

260

4.2 Quantification of water uptake as related to mineralogy

261

Since water sorption on shales is mostly attributed to clays, a linear combination approach is applied

262

for predicting sorption isotherms of water on shales based on the mass fractions of mineral components.

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

263

This concept has been successfully used to evaluate several rock physical parameters in our previous

264

studies.15,16 It was assumed that water sorption on framework minerals (quartz, calcite, etc.) is negligible,

265

and the water uptake of the shales is expressed as the sum of amounts of adsorbed water on individual

266

mineral components in mass fraction. Combined with the available water isotherms for typical clay

267

minerals (Figure 7a),35 the constructed water sorption isotherms of shale samples are obtained. The

268

calculated water sorption isotherms reasonably fit the measured data at most of the relative humilities

269

range, and show typical inversed S-shaped Type II isotherm (Figure 7b). This demonstrates the

270

applicability of the developed method, and indicates that the water uptake of shale samples is primarily

271

associated with clay minerals. The calculated water uptake overestimates the measured value of CN_22

272

by a factor of 1.4 at 97% RH. The difference between the calculated water isotherms and the measured

273

data in high relative humidities could be attributed to the fact that the milling process may increase the

274

accessible internal surface and finely ground (< 2 µm in Martin 35) provides better accessibility to water

275

than the 150 µm particles in present study. It is also possible that the water uptake behavior of pure clay

276

minerals, not constrained by the rock fabric, is different from clays components of rocks. The Paleozoic

277

marine shales in Sichuan Basin experienced deep burial before an intense uplift, where the high

278

temperature and pressure changed the pore structure of the clay minerals. Many pores of soft clay

279

minerals in shales are lost because of strong post-compaction, especially when the rigid minerals of the

280

sample are low,16 which will significantly reduce the water content of samples (CN_33 in Figure 6b). In

281

addition, the sorbed water of the sample will lose at a certain degree when the moisture-equilibrated

282

sample was transferred from the desiccator to the electronic balance for moisture measurement. The

283

moisture loss in the transfer process in our experiments is up to 0.5 wt.%. Nevertheless, Figure 7

284

demonstrates that water uptake of these shales is primarily associated with clay minerals, and the linear

14

ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28

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

285

combination scheme provides a convenient method to determine the water content of shales at any given

286

relative humidity.

287

4.3 Effect of water on methane sorption

288

The detrimental impact of water on gas sorption capacity has been previously reported in

289

coals.20,29,38–43 However, there are relatively few detailed studies of methane sorption on

290

moisture-equilibrated shales.9,12,13 Measurements of methane sorption capacity of shale samples from

291

Sichuan Basin under dry and moisture-equilibrated conditions show that sorption capacity of dry samples

292

is substantially greater than that under fully moisture-equilibrated condition (97% RH). There are obvious

293

positive relationships between TOC content and maximum methane sorption capacity of shales under

294

both dry and moisture-equilibrated conditions at 97% RH (Figure 8), demonstrating the first-order effect

295

of organic matter rather than other minerals on methane sorption. This agrees well with the observations

296

in Posidonia and Alum shales

297

between clay minerals content and sorption capacity of these shales (Figure 8a). Clay minerals generally

298

play a role in gas sorption of low-TOC, clay-rich shales in the dry state.8,15,18

9

and Devonian-Mississippian shales 6. No obvious relationship exists

299

The impact of water on methane sorption capacity of shales is further illustrated in Figure 9, which

300

demonstrates the reduction in methane sorption because of pre-sorbed water. The marine shales studied

301

show a variation of decreasing sorption capacities owing to water uptake ranging from 44 to 63%

302

depending on mineral components, type, abundance, distribution, and PSD. Experimental results

303

demonstrate a step-wise decline in methane capacity with increasing water content. Combined with the

304

water uptake process (Figure 4), the evolution of methane sorption capacity of shales as a function of

305

water content can be roughly divided into three stages: (1) initial decline stage: At the beginning of water

306

uptake, water molecules preferentially attach to the oxygen-containing functional groups, and occupy 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

307

polar sites of hydrophilic clay minerals. This process is also referred to as competitive sorption of water

308

and methane to the surface sites of shales. Monolayer coverage is complete when all the polar sites are

309

covered by water molecules. The monolayer moisture capacities (wm) calculated from the GAB model are

310

reached at 33% RH for CN_11 and CN_33, while they are reached at relative humidities between 33%

311

and 53% for CN_22 and CQ_14. The loss of sorption capacity at this stage is mainly related to the type,

312

nature, and abundance of hydrophilic clay minerals. Shale samples with higher clay mineral contents have

313

a greater reduction of methane sorption capacity. The reduction in sorption capacity of clay-rich CN_33 is

314

about 26% after the monolayer coverage of water, while the lost sorption capacities are less than 16% for

315

other samples. The type of clay minerals also play a role on the water uptake and thus reduction in

316

methane sorption. Water uptake of clay minerals generally decreases in the order: montmorillonite > illite >

317

kaolinite (Figure 7a).35 The composition of clay minerals in our shales is dominated by illite, which

318

adsorbs a lot of water. Conversely, water sorption capacity of organic matter is negligible in our highly

319

mature shale samples because many oxygen-containing functional groups (such as O-H) in the organic

320

matter are lost coupled with the generation of hydrocarbon,36,37 and thus the higher TOC content, the

321

lower relative loss in methane sorption of shales at the initial stage of water uptake (Table 4). CQ_14

322

shows relatively lower loss in methane sorption capacity after the monolayer coverage of water, due to

323

the low clay minerals content and high TOC content. (2) steep decline stage: At this stage, clusters of

324

water molecules form around the polar sites via hydrogen bonding among adjacent water molecules

325

following the monolayer moisture sorption. These water clusters block pore volume that would otherwise

326

be available for methane and thus remarkably reduce the methane sorption. The blockage of pore space

327

due to water clusters is dependent on the pore structure characteristics. Capillary condensation

328

preferentially occurs in micropores, which are associated with organic kerogen of shales. TOC-rich

16

ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28

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

329

CQ_14 has high micropore volume, and shows considerable water uptake at > 53% RH, which

330

corresponds to capillary condensation of water in pores. Because many micropores are blocked by water

331

molecules, methane sorption capacities of the TOC-rich sample reduce remarkably. Methane sorption

332

capacity of CQ_14 decreases from 0.21 (53% RH) to 0.12 mmol/g (75% RH), meaning that about 42% of

333

initial dry methane sorption capacity is lost in this stage. The loss of sorption capacity is about 30% for

334

CN_11 (TOC=4.58%) from 33% RH (0.15 mmol/g) to 75% RH (0.10 mmol/g). The loss in methane

335

sorption capacity of low-TOC samples is relatively less than for the TOC-rich samples at this stage,

336

possibly because capillary condensation is of lesser importance here. About 20% of dry methane sorption

337

capacity is lost for CN_33 (TOC=0.96%) from 33% RH (0.09 mmol/g) to 55% RH (0.06 mmol/g). (3)

338

slow decline stage: The loss of sorption capacity and the amount of water uptake correlates until a certain

339

moisture content, above which there is no significant change in sorption capacity with further increase in

340

moisture. At this stage, many micropores have been filled by water clusters, and the capillary

341

condensation continues in macropores and interparticle voids, successively filling the pore volume with a

342

contiguous water phase. Because water molecules aggregate in macropores and capillaries, which are

343

unavailable for methane sorption, methane sorption capacities would not decrease significantly. The slow

344

decline stage is also referred to as volumetric displacement of methane molecules by water and has been

345

reported for coal but investigated to a lesser extent in shales.20 The critical moisture content for these

346

samples is achieved between 75% and 97% RH except that clay-rich CN_33 show critical moisture

347

characteristic between 53% and 75% RH. CN_33 is rich in clay minerals and low in TOC content, thus

348

the main loss in methane sorption of CN_33 occurs in the early stage of water uptake. Furthermore, it is

349

not easy to form contiguous water phase in CN_33 and the volumetric displacement has relatively less

350

effect on methane sorption.

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

351

Page 18 of 28

Our step-wise decrease in sorption capacity with increasing water uptake in shales differs partly

352

from (but is compatible with) previous studies on gas sorption on moist coals.20 Day et al.

20

353

linear decrease in sorption capacity of coals up to a certain critical moisture content, in which the initial

354

decline stage is unconspicuous. This could be attributed to the high content of organic matter (TOC >

355

75%) and relatively low clay minerals (clay minerals content (g/cm3) is the grain density.

374

Considering the skeletal density and porosity values derived from helium pycnometry (Table 2), the

375

amount of water at 97% RH would yield somewhat high water saturation values (>70%) in shales. Taking

376

the water contents at 75% RH, the equivalent water saturation are calculated as 27%–67% for these shales,

377

which are in the range of actual water saturation vales (10%–70%) reported for the Silurian Longmaxi

378

shales from Sichuan Basin.44 Consequently, sorption isotherms measured on dry samples will

379

overestimate sorbed gas in place, especially in clay-rich formation. Sorption capacity provided by organic

380

matter might be partly decreased because of pore blockage and volumetric displacement of water as an

381

additional effect. It is highly advisable to perform methane sorption on moisture-equilibrated samples.

382

Methane sorption capacity derived from dry sample can be considered as a maximum case, and sorption

383

data of partially moisture-equilibrated sample (75% RH) is a more realistic estimation, and that for fully

384

moisture-equilibrated samples (97% RH) is regarded as a minimum value.

385

Gas transport in unconventional shales is a combination of several flow mechanisms: desorption,

386

diffusion, and slip flow in micropores, Darcy flow (advection) in macropores of the matrix and fractures

387

networks. Capillary condensation of water in pores and capillaries not only occupies the pore spaces for

388

gas storage, but also restricts the gas transport in interconnected pores (transport pores). Experiments on

389

methane sorption and diffusion on shales have shown that there is a considerable decrease in the sorption

390

rates for moisture-equilibrated shales while the uptake rate of corresponding dry shales is relatively

391

rapid.9,33 Moreover, surface diffusion, sorbed gas molecules hopping from sorption sites to the

392

neighboring sorption sites, plays a key role in gas mass transfer of shale reservoirs.45 The presence of

393

water clusters in the pores and microfractures of shales certainly has a detrimental impact on the

394

transportation of both sorbed and free gas. Detailed investigations of kinetics of desorption and gas

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

395

transport mechanism within shales will be discussed in future work.

396

5. Conclusions

397

Methane sorption isotherms on shales of Sichuan Basin at five different moisture levels were

398

measured up to pressures of 25 MPa using a manometric method. Water sorption isotherms of these

399

marine shales were also recorded at 24 °C and relative humidities (RH) from 11 to 97%. The purpose of

400

the present study was to investigate the effect of organic richness, mineral composition, and pore system

401

structure on water uptake and methane sorption characteristics of shales. The following conclusions can

402

be drawn:

403

(1) Small amounts of water have a detrimental effect on methane sorption of shales. Sorption

404

capacities of the moisture-equilibrated shales at 97% RH (0.05–0.08 mmol/g) were 44%–63% of those for

405

dry samples (0.11–0.22 mmol/g).

406

(2) The water sorption capacity of shales is controlled by clay mineralogy and a linear combination

407

scheme on a mass-fraction base provides reasonable predictions of water uptake isotherms for these

408

shales. Total clay mineral contents show a positive relationship with water uptake at 97% RH although

409

deviations from this trend are observed for clay-rich samples with low quartz contents.

410

(3) Methane sorption capacity of the marine shales studied is mainly controlled by the organic matter

411

content. Maximum sorption capacities of these shales under both dry and fully moisture-equilibrated

412

conditions correlate positively with TOC content.

413

(4) Experimental results demonstrate a step-wise decline in sorption capacity with increasing water

414

content. Evolution of sorption capacity as a function of water content of shales can be roughly

415

subdivided into three stages: (1) initial decline stage: loss in sorption capacity is attributed to

416

competitive sorption of gas and water molecules on the polar sites of hydrophilic clay minerals. (2) steep

20

ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28

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

417

decline stage: clusters of water molecules block pore space and remarkably reduce methane sorption. (3)

418

slow decline stage: successive filling of macropores with a contiguous water phase and slight reduction

419

of methane sorption capacity by volumetric displacement.

420

(5) Water saturation calculated based on the porosity values and water contents at 75% RH agree

421

with actual water saturation data reported for the Silurian Longmaxi shales. This suggests that

422

experimental sorption measurements performed on partially moisture-equilibrated samples (75% RH)

423

probably provide more realistic estimations for the in-situ sorption of shale reservoirs. Methane sorption

424

capacities measured on dry samples can be considered as maximum values, and those obtained for fully

425

moisture-equilibrated samples (97% RH) are regarded as a minimum values.

426

Acknowledgments

427 428 429

Foundation of China (Grant No. 51604249, 41572109), and State Key Laboratory of Petroleum

430

References

431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450

(1) Energy Information Administration (EIA). Annual Energy Outlook 2016: With Projections to 2040,

The authors would like to acknowledge the financial support of the National Natural Science Resources and Prospecting Independent Research Subject (Grant No. PRP/open-1606).

U.S. Department of Energy, DOE/EIA–0383, released on August, 2016. (2) Lu, X.C.; Li, F.C.; Watson, A.T. Adsorption measurements in Devonian shales. Fuel 1995,74,599– 603. (3) Curtis, J.B. Fractured shale-gas systems. AAPG Bull. 2002, 86, 1921–1938. (4) Chalmers, G.R.L.; Bustin, R.M. The organic matter distribution and methane capacity of the Lower Cretaceous strata of Northeastern British Columbia, Canada. Int. J. Coal Geol. 2007, 70, 223–239. (5) Chalmers, G.R.L.; Bustin, R.M. Lower Cretaceous gas shales in northeastern British Columbia, part 1: geological controls on methane sorption capacity. Bull. Can. Petrol. Geol. 2008, 56,1–21. (6) Ross, D.J.K.; Bustin, R.M. The importance of shale composition and pore structure upon gas storage potential of shale gas reservoirs. Mar. Petrol. Geol. 2009, 26, 916–927. (7) Zhang, T.W.; Ellis, G.S.; Ruppel, S.C.; Milliken, K.; Yang, R.S. Effect of organic-matter type and thermal maturity on methane adsorption in shale-gas systems. Org. Geochem. 2012, 47, 120–131. (8) Gasparik, M.; Ghanizadeh, A.; Bertier, P.; Gensterblum, Y.; Bouw, S.; Krooss, B.M. High-pressure methane sorption isotherms of black shales from the Netherlands. Energy Fuels 2012, 26, 4995– 5004. (9) Gasparik, M.; Bertier, P.; Gensterblum, Y.; Ghanizadeh, A.; Krooss, B.M.; Littke, R. Geological controls on the methane storage capacity in organic-rich shales. Int. J. Coal Geol. 2014, 123, 34–51. (10) Rexer, T.F.T.; Benham, M.J.; Aplin, A.C.; Thomas, K.M. Methane adsorption on shale under simulated geological temperature and pressure conditions. Energy Fuels 2013, 27, 3099–3109. 21

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

451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495

(11) Rexer, T.F.T.; 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. (12) Tan, J.Q.; Weniger, P.; Krooss, B.; Merkel, A.; Horsfield, B.; Zhang, J.C.; et al. Shale gas potential of the major marine shale formations in the Upper Yangtze Platform, South China, Part II: Methane sorption capacity. Fuel 2014, 129, 204–218. (13) Merkel, A.; Fink, R.; Littke, R. High pressure methane sorption characteristics of lacustrine shales from the Midland Valley Basin, Scotland. Fuel 2016,182,361–372. (14) Vandewijngaerde, W.; Piessens, K.; Dusar, M.; Bertier, P.; Krooss, B.M.; Littke, R.; et al. Investigations on the shale oil and gas potential of Westphalian mudstone successions in the Campine Basin, NE Belgium (well KB174): Palaeoenvironmental and palaeogeographical controls. Geol. Belg. 2016, 19, 1–11. (15) Yang, F.; Ning, Z.F.; Zhang, R.; Zhao, H.W.; Krooss, B.M. Investigations on the methane sorption capacity of marine shales from Sichuan Basin, China. Int. J. Coal Geol. 2015,146, 104–117. (16) Yang, F.; Ning, Z.F.; Wang, Q.; Zhang, R.; Krooss, B.M. Pore structure characteristics of lower Silurian shales in the southern Sichuan Basin, China: Insights to pore development and gas storage mechanism. Int. J. Coal Geol. 2016, 156, 12–24. (17) Xiong, W.; Zuo, L.; Luo, L.T.; Hu, Z.M.; Cui, Y.X. Methane adsorption on shale under high temperature and high pressure of reservoir condition: Experiments and supercritical adsorption modeling. Adsorpt. Sci. Technol. 2016, 34, 193–211. (18) Ji, L.M.; Zhang, T.W.; Milliken, K.; Qu, J.L.; Zhang, X.L. Experimental investigation of main controls to methane adsorption in clay-rich rocks. Appl. Geochem. 2012, 27, 2533–2545. (19) Kuila, U.; McCarty, D.K.; Derkowski, A.D.; Fischer, T.B.; Topór, T.; Prasad, M. Nanoscale texture and porosity of organic matter and clay minerals in organic-rich mudrocks. Fuel 2014, 135, 359–373. (20) Day, S.; Sakurovs, R.; Weir, S. Supercritical gas sorption on moist coals. Int. J. Coal Geol. 2008,74,203–214. (21) Busch, A.; Gensterblum, Y. CBM and CO2-ECBM related sorption processes in coal: a review. Int. J. Coal Geol. 2011,87,49–71. (22) Li, Y.F.; Shao, D.Y.; Lv, H.G.; Zhang, Y.; Zhang, X.L.; Zhang, T.W. A relationship between elemental geochemical characteristics and organic matter enrichment in marine shale of Wufeng Formation-Longmaxi Formation, Sichuan Basin. Acta Pet. Sin. 2015, 36, 1470–1483. (23) Liang, D.G.; Guo, T.L.; Chen, J.P.; Bian, L.Z.; Zhao, H. Some progress on studies of hydrocarbon generation and accumulation in marine sedimentary regions, Southern China (Part I): distribution of four suits of regional marine source rocks. Mar. Orig. Pet. Geol. 2008, 13, 1–16. (24) Wang, T.; Yang, K.M.; Xiong, L.; Shi, H.L.; Zhang, Q.L.; Wei, L.M.; et al. Shale sequence stratigraphy of Wufeng-Longmaxi Formation in sourthern Sichuan and their control on reservoirs. Acta Pet. Sin. 2015, 36, 915–925. (25) Gregg, S.J.; Sing, K.S.W. Adsorption, surface area and porosity, 2nd ed, Academic Press, New York; 1982. (26) Hinai, A.A.; Rezaee, R.; Esteban, L.; Labani, M. Comparisons of pore size distribution: A case from the Western Australian gas shale formations. J. Unconv. Oil Gas Resour. 2014, 8, 1–13. (27) Bertier, P.; Schweinar, K.; Stanjek, H.; Ghanizadeh, A.; Clarkson, C.R.; Busch, A.; et al. On the use and abuse of N2 physisorption for the characterization of the pore structure of shales. The Clay Minerals Society Workshop Lectures Series. 2016, 21, 151–161. (28) N.N. DIN EN ISO 483:2006–02, Norm-Beuth.de n.d. 22

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28

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

496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536

(29) Lewicki, P.P. The applicability of the GAB model to food water sorption isotherms. Int. J. Food Sci. Tech. 1997, 32, 553–557. (30) Krooss, B.M.; Van Bergen, F.; Gensterblum, Y.; Siemons, N.; Pagnier, H.J.; David, P. High-pressure methane and carbon dioxide adsorption on dry and moisture-equilibrated Pennsylvanian coals. Int. J. Coal Geol. 2002, 51, 69–92. (31) Gasparik, M.; Rexer, T.F.T.; Aplin, A.C.; Billemont, P.; De Weireld, G.; Gensterblum, Y.; et al. First international inter-laboratory comparison of high-pressure CH4, CO2 and C2H6 sorption isotherms on carbonaceous shales. Int. J. Coal Geol. 2014, 132,131–146. (32) Kunz, O.; Klimeck, R.; Wagner, W.; Jaeschke, M. The GERG-2004 Wide-Range Equation of State for Natural Gases and Other Mixtures. VDI Verlag, Düsseldorf 2007; 978–3–18–355706–6. (33) Yuan, W.N.; Pan, Z.J.; Li, X.; Yang, Y.X.; Zhao, C.X.; Connell, L.D.; et al. Experimental study and modelling of methane adsorption and diffusion in shale. Fuel 2014, 117, 509–519. (34) Passey, Q.R.; Bohacs, K.M.; Esch, W.L.; Klimentidis, R.; Sinha, S. From oil-prone source rock to gas-producing shale reservoir-geologic and petrophysical characterization of unconventional shale-gas reservoirs. SPE 131350 presented at the CPS/SPE International Oil & Gas Conference and Exhibition in Beijing, China, 8–10 June, 2010. (35) Martin, R.T. Adsorbed water on clay: A review. Clays and Clay Minerals: Proceedings of the Ninth National Conference on Clays and Clay Minerals, Lafayette, Indiana, October 5–8, 1960, 28–70. (36) Mahajan, O.P.; Walker, Jr. P.L. Water adsorption on coals. Fuel 1971,50,308–317. (37) Merkel, A.; Fink, R.; Littke, R. The role of pre-adsorbed water on methane sorption capacity of Bossier and Haynesville shales. Int. J. Coal Geol. 2015,147–148,1–8. (38) Ryan, B. A discussion on moisture in coal-implications for coalbed gas and coal utilization. BC Ministry of Energy. Mines Pet. Resour. 2006, 1, 139–149. (39) Gensterblum, Y., Merkel, A., Busch, A., Krooss, B.M. High-pressure CH4 and CO2 sorption isotherms as a function of coal maturity and the influence of moisture. Int. J. Coal Geol. 2013, 118, 45–57. (40) Müller, E.A.; Rull, L.F.; Vega, L.F.; Gubbins, K.E. Adsorption of water on activated carbons: A molecular simulation study. J. Phys. Chem. 1996, 100, 1189–1196. (41) Joubert, J.I.; Grein, C.T.; Bienstock, D. Effect of moisture on the methane capacity of American coals. Fuel 1974,53,186–191. (42) McCutcheon, A.L.; Barton, W.A. Contribution of mineral matter to water associated with bituminous coals. Energy Fuels 1999,13, 160–165. (43) Gensterblum, Y.; Merkel, A.; Busch, A.; Krooss, B.M.; Littke, R. Gas saturation and CO2 enhancement potential of coalbed methane reservoirs as a function of depth. AAPG Bull. 2014, 98, 395–420. (44) Wang, F.Y.; Guan, J.; Feng, W.P.; Bao, L.Y. Evolution of overmature marine shale porosity and implication to the free gas volume. Pet. Explor. Dev. 2013, 40, 819–824. (45) Chen, Y.D.; Yang, R.T. Concentration dependence of surface diffusion and zeolitic diffusion. AIChE J. 1991, 37, 1579–1582. (46) Feng, G.X.; Chen, S.J. Relationship between the reflectance of bitumen and vitrinite in rock. Nat. Gas Ind. 1988, 8, 20–24.

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

537

Figures

538 539 540

Figure 1. Map showing geological structure of the Sichuan Basin and sampling locations (modified after Yang et al. 15). The roman numbers represent different fold belts in Sichuan Basin.

541

542 543 544

Figure 2. (a) Low-pressure N2 adsorption-desorption isotherms for the shale samples; (b) Pore size distributions of the samples based on the BJH method applied to the adsorption branch of the isotherms.

545

546 547 548

Figure 3. Relationship between TOC content and (a) BET specific surface area and (b) average pore diameters of the samples.

549

24

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28

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

550 551 552

Figure 4. Water sorption isotherms for the samples at 297 K (24 °C) (The lines represent GAB fitting of the measured water sorption isotherms).

553

554 555 556 557

Figure 5. Methane excess sorption isotherms of the shale samples measured at varying relative humidity with corresponding fitted excess sorption functions (It should be noted that the evolution of the isotherms of CQ_14 is different from others).

558

559 560 561 562 563

Figure 6. Effect of (a) TOC content and (b) total clay minerals and (c) BET surface area on water content of the 97% RH moisture-equilibrated shale samples (In Figure 6a: The dote line represents the linear fitting of samples in this study (R2=0.61); the dash line represents the linear fitting of European samples measured by Gasparik et al.9 (R2=0.30)). 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

564 565 566 567

Figure 7. (a) Water sorption isotherms of clays from Martin35; (b) Measured (solid symbols) and estimated (lines) of water uptake on selected shales samples. The water uptake isotherms are estimated based on individual mineral components on a mass-fraction base.

568

569 570 571 572 573

Figure 8. Methane sorption capacity of dry and fully moisture-equilibrated (97% RH) shale samples as a function of (a) total clay minerals and (b) TOC content. Squares and circles denote samples measured by Gasparik et al.9, and Ross and Bustin6, respectively. Gas sorption capacity of Devonian-Mississippian shales in the WCSB are measured by Ross and Bustin6 at 30 °C and 6 MPa pressure.

574

575 576 577

Figure 9. Methane sorption capacity as a function of water content of the shale samples. ①, ②, and ③ denote initial decline, steep decline, and slow decline stage, respectively.

26

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

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

Tables Table 1: Geochemical characteristics and mineralogy composition of the shale samples. Reqb

TOC a Sample

Stratigraphy

Quartz +

Total Carbonates

feldspars (wt. %)

Illite

Kaolinite

Chlorite

(wt. %)

(wt. %)

(wt. %)

clays

(%)

(wt. %) (wt. %)

(wt. %)

CN_11

U. Ordovician

4.83

2.8

74.0

0.0

26.0

25.0

0.3

0.8

CN_22

L. Silurian

2.87

2.8

55.4

25.4

17.2

16.2

0.0

1.0

CN_33

L. Silurian

0.96

2.8

25.3

27.7

45.5

35.0

0.0

10.5

CQ_14

L. Silurian

9.40

2.4

66.7

4.1

22.6

22.6

0.0

0.0

a

The TOC content were determined by a LECO CS230 carbon/sulfur analyzer.

b

The equivalent vitrinite reflectance (Req) were derived from bitumen reflectance according to the equation of Feng

and Chen 46.

Table 2: Pore structure parameters of the shale samples. a

b

Porosity

BET surface area

Pore volume

Average pore diameter

Equivalent water

(g/cm )

(g/cm )

(%)

(m2/g)

(cm3/g)

(nm)

saturation (%)

CN_11

2.22

2.56

13.0

28.75

0.043

6.0

27.0

CN_22

2.54

2.66

4.4

14.96

0.018

4.9

65.3

CN_33

2.51

2.74

8.4

14.95

0.027

7.1

38.0

CQ_14

2.36

2.52

6.4

31.58

0.051

2.4

67.3

Bulk density Sample

3

Skeletal density 3

a

The bulk density was determined by caliper measurements on cylindrical plug samples.

b

The skeletal density were measured using a helium pycnometer based on Boyle's Law.

c

The equivalent water saturation of shales were calculated using water contents at 75% RH according to Eq. (5).

27

ACS Paragon Plus Environment

c

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 28 of 28

Table 3. Results from GAB Adsorption Model Analysis on water uptake isotherms

Shale

Clay mineral a

a

Sample

wm (mmol/g)

c

k

CN_11

0.31

17.67

0.89

CN_22

0.31

8.32

0.72

CN_33

0.25

9.39

0.87

CQ_14

0.34

7.95

0.90

Montmorillonite

8.08

7.70

0.70

Illite

1.40

180.5

0.87

Kaolinite

0.16

844.67

0.95 35

The water uptake data of clay minerals were collected from Martin .

Table 4. Moisture content, and methane sorption capacities and Langmuir parameters for shale samples. Sample

CN_11

CN_22

CN_33

CQ_14

Relative humidity

Water content

Water content

nL

PL

ρads

Lost sorption

(%)

(wt. %)

(mmol/g)

(mmol/g)

(MPa)

(kg/m3)

capacity (%)

0

0

0

0.181

3.3

433.4

0

11

0.38

0.21

-

-

-

-

33

0.72

0.40

0.153

7.1

566.7

15.5

53

1.05

0.58

0.143

10.0

648.5

21.0

75

1.58

0.88

0.099

7.7

744.9

45.3

97

4.06

2.26

0.083

7.1

916.9

54.1

0

0

0

0.113

2.7

471.4

0

11

0.29

0.16

-

-

-

-

33

0.50

0.28

0.099

6.1

551.9

12.4

53

0.72

0.40

0.095

13.2

671.8

15.9

75

1.13

0.63

0.069

9.1

778.4

38.9

97

1.76

0.98

0.063

8.7

746.9

44.2

0

0

0

0.118

7.8

508.2

0

11

0.29

0.16

-

-

-

-

33

0.46

0.26

0.087

10.0

673.1

26.3

53

0.67

0.37

0.063

10.5

-

46.6

75

1.27

0.71

0.058

11.4

-

50.8

97

2.75

1.53

0.051

11.0

-

56.8

0

0

0

0.222

2.5

405.4

0

11

0.38

0.21

-

-

-

-

33

0.53

0.30

0.215

4.2

428.0

3.2

53

0.72

0.40

0.214

5.3

376.7

3.6

75

1.81

1.01

0.123

9.0

396.9

44.6

97

4.50

2.50

0.082

17.3

-

63.1

28

ACS Paragon Plus Environment