Investigation of accessible pore structure evolution under

Jan 7, 2019 - ... could damage the nanopores in the Hazleton coal and Marcellus shale outcrop. Thus, the pore structure stability may be sample depend...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF BARCELONA

Fossil Fuels

Investigation of accessible pore structure evolution under pressurization and adsorption for coal and shale using small-angle neutron scattering Shimin Liu, Rui Zhang, Zuleima T. Karpyn, Hongkyu Yoon, and Thomas Dewers Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.8b03672 • Publication Date (Web): 07 Jan 2019 Downloaded from http://pubs.acs.org on January 17, 2019

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 31 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 2 3 4 5 6 7 8 9 10 11 12 13 14 15

Energy & Fuels

Investigation of accessible pore structure evolution under pressurization and adsorption for coal and shale using small-angle neutron scattering Shimin Liua, *, Rui Zhanga, **, Zuleima Karpyna, Hongkyu Yoonb, Thomas Dewersb aDepartment

of Energy and Mineral Engineering, G3 Center and Energy Institute, The Pennsylvania State University, University Park, PA 16802, USA bGeoscience Research and Applications, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA Corresponding Authors: *Tel. No.: +1 8148634491; Fax No.: +1 8148653248; Email: [email protected] **Emails: [email protected] Abstract

16

Pore structure is an important parameter to quantify the reservoir rock adsorption capability

17

and diffusivity, both of which are fundamental reservoir properties to evaluate the gas production

18

and carbon sequestration potential for coalbed methane (CBM) and shale gas reservoirs. In this

19

study, we applied small-angle neutron scattering (SANS) to characterize the total and accessible

20

pore structures for two coal and two shale samples. We carried out in situ SANS measurements to

21

probe the accessible pore structure differences under argon, deuterated methane (CD4) and CO2

22

penetrations. The results show that the total porosity ranges between 0.25 and 5.8 % for the four

23

samples. Less than 50 % of the total pores are accessible to CD4 for the two coals, while more than

24

75 % of the pores were found to be accessible for the two shales. This result suggests that organic

25

matter pores tend to be disconnected compared to mineral matter pores. Argon pressurization can

26

induce pore contraction because of the mechanical compression of the solid skeleton in both the

27

coal and shale samples. Hydrostatic compression has a higher effect on the nanopores of coal and

28

shale with a higher accessible porosity. Both methane and CO2 injection can reduce the accessible

29

nanopore volume due to a combination of mechanical compression, sorption-induced matrix

30

swelling and adsorbed molecule occupation. CO2 has a higher effect on sorption-induced matrix

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

Page 2 of 31

31

swelling compared to methane for both the coal and shale samples. Gas densification and pore

32

filling could occur at higher pressures and smaller pore sizes. In addition, the compression and

33

adsorption could create nanopores in the San Juan coal and Marcellus shale drilled cores but could

34

have an opposite effect in the other samples, namely, the processes could damage the nanopores

35

in the Hazleton coal and Marcellus shale outcrop. Thus, the pore structure stability may be sample

36

dependent.

37

1. Introduction

38

Coalbed methane (CBM) and shale gas are two important unconventional natural gas

39

resources worldwide.1 Coal and gas shale are both self-sourced reservoir rocks where the gas is

40

stored in both free phase and adsorbed phase. The pore structure of the coal and gas shale matrix

41

is known to be heterogeneous and anisotropic. The pore structure change will depend on not only

42

the pressure variation but also the gas-rock interaction effect through the sorption process during

43

gas production in both CBM and shale gas reservoirs. Therefore, it is crucial to investigate the

44

evolution of pore structure alteration under in situ gas depletion and injection.

45

Previous studies in the literature primarily focused on the gas-induced matrix deformation

46

of coal and shale. It was found that the CO2-induced swelling and strain distribution in coal are

47

heterogeneous, depending on the lithotypes and microstructures.2,

48

microlithotypes exhibited CO2-induced swelling, while the clay and inertite regions showed

49

compression in coal.4 The volatile matter will affect the pore structure evolution during the

50

carbonization process.5 Methane adsorption and expansion mainly occurred in the pore structure

51

filled with clay minerals.6 The investigation of competitive CO2-CH4 adsorption in coal showed

52

that the sorption-induced swelling appeared in the low CO2 mole fraction region, which is

53

insensitive to temperature and pressure conditions.7 The internal swelling due to supercritical CO2

3

Vitrite, liptite and clarite

2 ACS Paragon Plus Environment

Page 3 of 31 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

54

was much greater than that for gaseous CO2, and this difference increased with increasing pore

55

pressure in coal.8 Compared with the different sorptive gases and different types of coal, the

56

sorption-induced swelling was greater with CO2 compared to CH4.9 Moreover, CH4 caused

57

swelling more slowly than CO2 at the same pressure.10 Lower rank coals exhibit a higher swelling

58

effect than higher rank coals,9 where the presence of moisture significantly reduced the gas-

59

induced swelling compared to dry coal samples. The deformation was larger in the direction

60

perpendicular to the bedding plane compared to the direction parallel to the bedding plane for both

61

water and gas adsorption under hydrostatic conditions.11, 12 It was also found that the mechanical

62

deformation could be fully reversible in coal.13 From the simulation perspective, 14, 15 the smallest

63

micropores (0.5 nm) tended to be compressed at low pressure (1-10 MPa) initially, followed by

65

expansion with increasing pressure. However, only a few studies have focused on the deformation

66

of gas shale or pure clay minerals. It was found that both water and CO2 can occupy the

67

montmorillonite interlayers, causing strength reduction.16 With increasing CO2 pressure, the

68

swelling of shale initially increased and then decreased.17 Moreover, the maximum swelling

69

gradually decreased with increasing temperature for the shale sample. Compared with the coal

70

sample, the gas adsorption-induced volumetric swelling strain for shale was approximately one

71

magnitude lower than that for coal.18 For practical purposes, the sorption-induced anisotropic

72

swelling of shale matrix has a significant impact on the caprock sealing efficiency for carbon

73

sequestration.19

74

It is notable that most of the deformation studies are based on bulk coal and shale samples,

75

where the detection scale is in the macroscopic range. However, few studies have considered the

76

deformation at nanopore scales in the rock matrix, especially in meso-/micropores. Since ad-

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

Page 4 of 31

77

/desorption and diffusion mainly occur in the meso-/micropore size range in the coal and gas shale

78

matrix, it is critical to investigate the evolution of the pore structure as a function of the pressure

79

for both inert and sorptive gases. In this study, four samples, including San Juan coal, Hazleton

80

coal, Marcellus drilled core and outcrop shales, were used, where the results of the pore size-

81

dependent pore accessibility determination of the two coals and the evolution of the fractal

82

dimension under pressurization and adsorption of the four tested samples are shown in our

83

previous studies.20, 21 The total and accessible pore structures were characterized based on small-

84

angle neutron scattering (SANS) experiments. The in situ SANS measurements were conducted to

85

evaluate the evolution of the pore size distribution (PSD) and porosity of the accessible pores with

86

penetration of argon, deuterated methane (CD4) and CO2. The outcomes of this study will have

87

important implications for gas storage estimates, transport modeling and gas production

88

optimizations in CBM and shale gas reservoirs.

89

2. Experimental methods

90

2.1 Sample preparation and characterization

91

The subbituminous coal (San Juan coal) and anthracite (Hazleton coal) were collected from

92

the northern San Juan Basin in New Mexico and the town of Hazleton in Pennsylvania,

93

respectively. Two shale samples were collected from the Marcellus shale formation in

94

Pennsylvania. One was a drilled core sample, and the other was an outcrop sample. Both the coals

95

and shales were pulverized to a particle size of ~0.5 mm for the SANS measurements and were

96

comprehensively characterized to obtain the necessary input parameters for the modeling and pore

97

structure evolution interpretation. The chemical composition of each rock sample was obtained by

98

X-ray diffraction (XRD), as shown in Table S1 in the supplementary section. The chemical

99

composition was used to estimate the effective scattering length density (SLD) and electron density

4 ACS Paragon Plus Environment

Page 5 of 31 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

100

for the scattering data fitting for each sample (Table S1). We found that the two coal samples had

101

a lower percentage of mineral content and a higher percentage of total organic carbon (TOC)

102

content compared to the two shale samples. The San Juan coal contained minerals such as quartz

103

and kaolinite, whereas only tobelite was encountered in the Hazleton coal. The TOC content was

104

63.95 % for the San Juan coal and 88.10 % for the Hazleton coal. As for the two Marcellus shale

105

samples, the drilled core sample had the most diverse mineral content, whereas only quartz and

106

muscovite were found in the outcrop sample. The drilled core sample had the lowest TOC content

107

(2.06 %), and there was 8.48 % organic matter in the Marcellus outcrop shale. The outcrop shale

108

sample had the highest quartz content (75.41 %), which could be due to the strong effect of

109

weathering.

110

2.2 SANS experiments

111

The SANS measurements were conducted using the general-purpose small-angle neutron

112

scattering (GP-SANS) instrument in High Flux Isotope Reactor (HFIR) at Oak Ridge National

113

Laboratory (ORNL).22 The pulverized powder sample with a particle size of ~0.5 mm was

114

uniformly loaded in the aluminum sample cell with an average thickness of ~1.65 mm for the

115

testing of the samples. The median particle size could optimize the average scattering for all

116

orientations of the pores and avoid scattering from the interparticle voids for the coal and shale

117

samples. The effect of multiple scattering could be maximally reduced in samples with small

118

thicknesses.23 The neutron wavelength was set to 6 Å, and the neutron wavelength spread 𝛥𝜆/𝜆

119

was 0.13. The detector distances were chosen at 18.5 and 0.3 m to cover an overall range of

120

scattering vectors between 3×10-3 and 0.5 Å-1, which corresponded to the pore size range between

121

0.5 and 82 nm.24 All the 2D scattering profiles were radially averaged during data reduction using

122

Igor macros.25 The representative profiles are shown in Figure S1 in the supplementary section.

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 31

123

All the scattering intensities were normalized to absolute intensities by using the secondary

124

standard.26

125

Three gases, including inert and sorptive gases (argon, CD4 and CO2), were used for the in

126

situ SANS measurements. The scattering at vacuum was determined for each sample as a reference.

127

First, CD4 was injected in incremental pressure steps at a pressure of 20, 40 or 68 bar. Then, the

128

sample cell was evacuated, and argon was injected through incremental pressure steps at a pressure

129

of 68, 340 or 476 bar for each sample. Finally, CO2 was injected, and evacuation was also

130

conducted before the injection procedure. The incremental injecting pressure steps for CO2 were

131

the same as for CD4 at an injection pressure of 20, 40 or 68 bar. Notably, a contrast-matched

132

condition was achieved during CD4 injection at 340 bar for the San Juan coal and at 476 bar for

133

the Hazleton coal, Marcellus shale drilled core and outcrop samples.

134

3. Results and discussion

135

3.1 Experimental scattering data

136

Figure 1 shows the scattering intensities 𝐼(𝑄) of the San Juan coal exhibiting the

137

experimental 𝐼(𝑄) values under vacuum and contrast-matched conditions, as well as the in situ

138

experimental 𝐼(𝑄) values under argon, CD4 and CO2 injection conditions. The 𝐼(𝑄) decreased with

139

increasing scattering of vector 𝑄 in the entire 𝑄 range for all the scattering profiles of the San Juan

140

coal. In the contrast-matched condition, the scattering counts at the detector reached the lowest

141

value. However, for the 1D scattering profile, the 𝐼(𝑄) under contrast-matched conditions was

142

smaller than that under vacuum conditions only in the low 𝑄 region. When 𝑄 was above 0.06 Å-1,

143

the 𝐼(𝑄) at vacuum decreased. One reason could be the densification of CD4 in the small pores23,

144

where the scattering contrast is relatively high. Another possibility is that the 𝑄-independent

145

scattering background increased when CD4 injection occurred at high pressure.

6 ACS Paragon Plus Environment

Page 7 of 31 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

146 147 148 149 150

Figure 1. Representative scattering intensities (San Juan coal) (a) in vacuum and contrast-matched conditions; (b) argon injection; (c) CD4 injection and (d) CO2 injection.

151

range for all the aforementioned conditions (Figures 1b, c and d). This decrease in the scattering

152

intensity for argon injection (Figure 1b) could be primarily attributed to the shrinkage of pores due

153

to the solid skeleton contraction, where the pressure-induced mechanical compression on the solid

154

skeleton resulted in an increase in the solid density due to the grain contraction.27 In addition to

155

the mechanical compression effect due to argon injection, there were sorption-induced effects on

156

the microstructure of the rocks with regards to methane and CO2 adsorptions. As shown in Figures

157

1c and d, the 𝐼(𝑄) decreased more severely compared to argon in the low 𝑄 range. This is due to

158

the increased bulk density of methane or CO2 resulted in a reduction in the scattering contrast

For the in situ data, the 𝐼(𝑄) decreased with increasing injection pressure in the low 𝑄

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 31

159

between the pores and matrix. And also, there was a sorption-induced gas-solid interface

160

densification at which the sorption layer had a higher density compared to the bulk gas.23 In

161

addition, the methane and CO2 sorption-induced pore deformation may have affected the scattering

162

intensity due to the shrinkage/swelling effects.12 Interestingly, a large decrease in the scattering

163

intensity for CO2 at 68 bar was observed (Figure 1d). One possible reason is that CO2 entered the

164

liquid phase at 68 bar at room temperature, where the scattering contrast between the pores and

165

rock matrix experienced a severe drop compared to gaseous CO2. All the results for the other three

166

samples (the Hazleton coal, Marcellus shale drilled core and outcrop shales) were similar to the

167

abovementioned findings of the San Juan coal; the experimental scattering intensities are shown

168

in Figures S2-S4 in the supplementary section.

169

3.2 Pore structure of the total and accessible pores

170

Figure 2 shows the pore size distribution (PSD) of the total and accessible pores for the

171

tested four samples. The PSD of the total pores was estimated using the scattering intensity under

172

vacuum. The PSD of the accessible pores was determined by the scattering intensity difference

173

between the vacuum and contrast-matched conditions. A model regression method was used to

174

estimate the experimental scattering intensity to extract the PSD information.28 Based on this

175

method, a polydisperse spherical pore (PDSP) model was used to fit the experimental scattering

176

data of the four tested rock samples. The scattering intensity 𝐼(𝑄) can be expressed as:29

177

𝐼(𝑄) = (𝜌s∗ ― 𝜌f∗ ) ∫𝑄max𝑉2(𝑟)𝑓(𝑟)𝑃(𝑄,𝑟)𝑑𝑟 +𝐵𝑘𝑔

178

where (𝜌s∗ ― 𝜌f∗ ) is the scattering contrast between the pores and surrounding rock matrix where

179

𝜌s∗ and 𝜌f∗ are the SLD of the solid matrix and pores, respectively; 𝑄min and 𝑄max are the lower

180

and upper limits of 𝑄, respectively; 𝑉(𝑟) is the spherical pore volume with a radius 𝑟; 𝑓(𝑟) is the

181

pore size distribution; 𝑃(𝑄, 𝑟) is the spherical form factor; and 𝐵𝑘𝑔 is the scattering background.

2 𝑄

min

(1)

2

8 ACS Paragon Plus Environment

Page 9 of 31 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

182

We defined the lower and upper limits of the pore diameter to be 0.5 and 200 nm, respectively,

183

during the model fitting. In general, the pore size distribution 𝑓(𝑟) was multimodal for the tested

184

samples, as shown in Figure 2. Notably, the pore range for the 𝑓(𝑟) of the accessible pores was

185

smaller than that of the total pores over the entire detected pore range. This result is because the

186

subtraction of the scattering intensity measured under vacuum from the scattering intensity

187

measured under contrast-matched conditions led to a shorter pore range (𝑄 range) with a greater

188

pore size (smaller 𝑄), which is mainly due to the incoherent scattering of hydrogen and/or the

189

condensation effect of high pressure methane at a high level of 𝑄.23 In general, the 𝑓(𝑟) of the total

190

pores was higher than that of the accessible pores for all the tested samples (Figure 2), indicating

191

a certain percentage of inaccessible pores inside the rock samples.21, 23, 30-33 The difference of 𝑓(𝑟)

192

between the total and accessible pores for the two shale samples (Figures 2c and d) was relatively

193

smaller than that for the two coal samples (Figures 2a and b), indicating that a significant number

194

of inaccessible pores were present in the tested coals, while only a small portion was present in the

195

tested shales.

9 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 10 of 31

196 197 198 199 200

Figure 2. Pore size distribution of total and accessible pores for (a) San Juan coal; (b) Hazleton coal; (c) Marcellus drilled core shale and (d) Marcellus outcrop shale.

201

method, which can be expressed as:34

202

𝑄INV = ∫𝑄max𝑄2𝐼(𝑄)𝑑𝑄 = 2𝜋2(𝜌s∗ ― 𝜌f∗ ) ∅(1 ― ∅)

203

where 𝑄INV is the Porod invariant and ∅ is the porosity. The pore size range for the estimation of

204

the total and accessible porosities is shown in Table S2. As shown in Figure 3a, the Hazleton coal

205

had the lowest total porosity of less than 0.5 %. The San Juan coal had a relatively higher value of

206

~4.5 %. Between the two Marcellus shale samples, the Marcellus shale drilled core sample had a

207

smaller porosity of ~2 %, while the Marcellus shale outcrop sample had a porosity that was higher

208

than 5.5 %. Furthermore, the fraction of accessible porosity was quantitatively estimated by the

Quantitatively, the total and accessible porosities can be estimated by the Porod invariant

𝑄

2

min

(2)

10 ACS Paragon Plus Environment

Page 11 of 31 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

209

total and accessible porosities, as shown in Figure 3b. Both coals had much lower fractions of

210

accessible porosities (lower than 50 %) compared to the shale samples (higher than 75 %). This

211

result suggests that the nanopores in the organic matter tended to be disconnected compared to the

212

mineral matter since the coal had a much higher organic carbon content compared to the shale.

213 214 215 216

3.3 Evolution of the accessible pore structure by gas injection

217

3.3.1 Coal samples under gas injection

218

3.3.1.1 Argon injection

Figure 3. (a) Total porosity and (b) fraction of accessible porosity for the tested samples.

219

Figure 4 shows the PSD and porosity evolution of the accessible pores under Ar injection

220

for the tested coal samples. The SLD contrast for the estimation of the PSD and porosity at each

221

pressure for each sample is shown in Table S3. The pore size range for the porosity estimation at

222

each pressure for each sample is shown in Table S4. The 𝑓(𝑟) slightly decreased with increasing

223

Ar pressure in pore sizes larger than ~10 nm for the San Juan coal, as shown in Figure 4a. The

224

𝑓(𝑟) peak at ~9 nm slightly shifted to a smaller pore size with increasing pressure, although there

225

were variations in the 𝑓(𝑟) peak at ~5 nm (the inserted figure in Figure 4a). The decreasing trend 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

Page 12 of 31

226

of the 𝑓(𝑟) was consistent with that of the porosity results, where the accessible porosity generally

227

decreased with increasing pressure, as shown in Figure 4c. This result suggests that the accessible

228

pore volume decreased during the hydrostatic gas pressurization in the detected pore range for the

229

San Juan coal. The results are different compared to those of a previous study in which there was

230

no obvious change in the pore structure during helium injection at 155 bar for bituminous coals.35

231

The latter suggests that the in situ high pressure (greater than ~150 bar) hydrostatic gas injection

232

could lead to nanopore compression, but low pressure injection (less than ~150 bar) would have a

233

negligible effect on the pore structure for bituminous coals.

234 235 236 237

Figure 4. Evolution of accessible pore size distribution under argon injection for (a) San Juan coal and (b) Hazleton coal; Evolution of accessible porosity under argon injection for (c) San Juan coal and (d) Hazleton coal.

12 ACS Paragon Plus Environment

Page 13 of 31 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

238

However, there were huge variations in the 𝑓(𝑟) for the Hazleton coal, and there was no

239

obvious trend in the 𝑓(𝑟) as a function of pressure (Figure 4b). Consistently, there was no obvious

240

change in the accessible porosity for the Hazleton coal (Figure 4d). This result means that there

241

could be a negligible effect of hydrostatic pressurization on semi-anthracite coal, which is different

242

from the subbituminous coal mentioned above. This phenomenon may be caused by the

243

significantly smaller accessible porosity (less than 0.08 %) for this high rank coal. Thus, the value

244

of the accessible porosity could be an important factor for the degree of nanopore compression by

245

hydrostatic gas pressurization. Notably, the detectable pore size greater than ~20 nm was due to

246

the relatively high scattering background for this typical high rank sample (Figure S2). We could

247

not estimate the 𝑓(𝑟) and accessible porosity of the Hazleton coal at 68 bar of argon injection due

248

to the large error in the obtained scattering intensity.

249

3.3.1.2 Methane injection

250

Figure 5 shows the PSD and porosity evolution of the accessible pores under methane

251

injection for the two coal samples. The SLD contrast for the estimation of the PSD and porosity at

252

each pressure for each sample is shown in Table S5. The pore size range for the porosity estimation

253

at each pressure for each sample is shown in Table S6. There were small changes in the 𝑓(𝑟) as a

254

function of the methane pressure at pore sizes greater than ~20 nm, while the 𝑓(𝑟) peak at a pore

255

size of ~15 nm decreased and the 𝑓(𝑟) peaks at pore sizes of ~10 and ~6 nm shifted to smaller

256

sizes with increasing pressure when the pressure was below 40 bar for the San Juan coal, as shown

257

in Figure 5a. This result indicates that, at a relatively low pressure (< 40 bar), the nanopores with

258

a size < ~20 nm were mainly affected by the methane injection. However, when the pressure was

259

68 bar, a huge variation in the 𝑓(𝑟) occurred, and the detectable pore size range shrunk. The

260

vanished nanopores with pore sizes smaller than ~4 nm and pore sizes ranging between ~6 and

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

Page 14 of 31

261

~20 nm could be due to methane densification in these nanopores.21, 23 The nanopores could be

262

filled with highly densified methane, where the SLD inside the nanopores approached that of the

263

coal matrix. Thus, these nanopores became “invisible” to the neutron beam. However, another

264

reason could be the increase in the incoherent scattering intensity with increasing methane pressure,

265

where the information of the nanopores with a smaller size was smeared. Consistent with 𝑓(𝑟), the

266

accessible porosity generally decreased with increasing methane pressure, as shown in Figure 5c,

267

which is consistent with the macroscopic matrix swelling of bituminous coals under constant

268

effective stresses.8 The accessible porosity significantly decreased at 68 bar, which could be caused

269

by methane densification in the nanopores with pore sizes smaller than 20 nm as described above.

270

Compared with the compressive effect of argon, the methane adsorption additionally showed not

271

only sorption-induced matrix swelling but also filling of the accessible pores at pore sizes < 20 nm

272

for the subbituminous coal sample.

14 ACS Paragon Plus Environment

Page 15 of 31 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

273 274 275 276 277

Figure 5. Evolution of accessible pore size distribution under CD4 injection for (a) San Juan coal and (b) Hazleton coal; Evolution of accessible porosity under CD4 injection for (c) San Juan coal and (d) Hazleton coal.

278

under all pressure conditions were observed as shown in Figure 5b; we could not find a trend for

279

𝑓(𝑟) as a function of the methane pressure. However, the accessible porosity slightly increased

280

with increasing pressure, as shown in Figure 5d, which is unexpected and different from the results

281

of the San Juan coal. Higher rank coal usually has a lower degree of sorption-induced matrix

282

swelling compared to lower rank coal.9 The created accessible pores when the pressure increased

283

for the Hazleton coal could be caused by the accessible pore volume slightly increasing with

284

increasing pressure for this tight coal with a significantly low porosity (< 0.5 %). The result may

For the Hazleton coal, relatively large variations in the 𝑓(𝑟) as a function of the pore size

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

Page 16 of 31

285

indirectly suggest that the fraction of accessible pores may be pressure dependent during gas

286

adsorption and desorption.36

287

3.3.1.3 CO2 injection

288

Figure 6 shows the PSD and porosity evolution of the accessible pores under CO2 injection

289

for the tested coal samples. The SLD contrast for the estimation of the PSD and porosity at each

290

pressure for each sample is shown in Table S7. The pore size range for the porosity estimation at

291

each pressure for each sample is shown in Table S8. Notably, we could not obtain the 𝑓(𝑟) at 68

292

bar for the San Juan coal due to the negative subtracted scattering intensity. There were variations

293

in the 𝑓(𝑟) as a function of the pore size with increasing pressure at pore sizes greater than ~7 nm

294

for the San Juan coal, as shown in Figure 6a. However, the 𝑓(𝑟) peak at a pore size of ~6 nm

295

decreased with increasing CO2 pressure and vanished at 40 bar (the inserted figure in Figure 6a).

296

CO2 densification could occur at 40 bar when the pore size was smaller than ~7 nm for the San

297

Juan coal. Additionally, rather than CO2 densification, the increase in the incoherent scattering

298

could be another reason as was discussed in the previous section. Compared with methane injection,

299

the densification of CO2 occurred at lower pressures. This phenomenon is consistent with the

300

porosity results in that, as shown in Figure 6c, the accessible porosity decreased with increasing

301

CO2 pressure, where the decrease rate was higher than that of methane in the same pressure range.

302

The CO2 injection resulted in a higher degree of sorption-induced matrix swelling compared to

303

methane.9 Also, the CO2 swelling rate was higher than that of methane at the same pressure for the

304

bituminous coals.10 In addition, the pore filling with CO2 could occur at a relatively lower pressure

305

(40 bar) compared to methane for the San Juan coal.

16 ACS Paragon Plus Environment

Page 17 of 31 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

306 307 308 309 310

Figure 6. Evolution of accessible pore size distribution under CO2 injection for (a) San Juan coal and (b) Hazleton coal; Evolution of accessible porosity under CO2 injection for (a) San Juan coal and (b) Hazleton coal.

311

the pressure was below 40 bar, as shown in Figure 6b. We observed a slight increase in the 𝑓(𝑟)

312

with increasing CO2 pressure. However, there was a significant increase in the 𝑓(𝑟) at 68 bar when

313

CO2 was injected. This result is consistent with the porosity results, where the accessible porosity

314

slightly increased with increasing pressure and significantly increased at 68 bar, as shown in Figure

315

6d. The pore accessibility could also increase with increasing pressure during CO2 injection. The

316

result at 68 bar (the liquid phase) was unexpected since the accessible porosity was nearly eight

317

times greater than the values at lower CO2 pressure (the gaseous phase). From Figure S1, we found

318

that the scattering intensity at 68 bar decreased in the low 𝑄 region but increased when 𝑄 was

For the Hazleton coal, there were variations in the 𝑓(𝑟) as a function of the pore size when

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

Page 18 of 31

319

greater than 0.02 Å-1 compared to the injection results at a lower CO2 pressure. One possible reason

320

could be that the estimated SLD inside the accessible pores may not be the SLD of the bulk CO2

321

at 68 bar for the Hazleton coal due to the very small accessible porosity. The SLD contrast between

322

the rock matrix and accessible pores may be underestimated, and thus, the PSD and porosity could

323

be overestimated. However, another possible reason could be that the accessible porosity in the

324

pore size range between 5.4 and 82 nm did significantly increase at 68 bar, at which pressure the

325

CO2 became liquid at room temperature. The liquid CO2 could partially fill the larger accessible

326

pores that were beyond the SANS detectable range, and a significant number of the smaller

327

accessible pores (perhaps the “bubble-like” pores inside the liquid CO2), which became “visible”

328

in the detectable pore range, increased for this particular high rank coal sample. Further analysis

329

needs to confirm the hypothesis.

330

3.3.1.4 Stability

331

Figure 7 shows the PSD and porosity of the total pores before and after the gas injection

332

for the tested coal samples. There were only small variations in the 𝑓(𝑟) of the total pores

333

(including both the accessible and inaccessible pores) when the pore size was > ~20 nm initially,

334

after the methane and argon injection, and before the CO2 injection, as shown in Figures 7a and b.

335

These results indicate that the nanopores with a size greater than ~20 nm were stable after

336

pressurization and adsorption for the two coal samples. However, greater variations in the 𝑓(𝑟)

337

existed for pore sizes < ~20 nm, which was different between the coal samples. For the San Juan

338

coal, the 𝑓(𝑟) peak at ~1 nm increased after methane and argon injection and before CO2 injection

339

(Figure 7a). This result is consistent with the porosity results, where total porosity also increased

340

(Figure 7c). For the Hazleton coal, the 𝑓(𝑟) slightly decreased at a pore size of ~8 nm, and the

341

𝑓(𝑟) peak at ~1 nm shifted to a smaller pore size (Figure 7b). Consistently, the total porosity

18 ACS Paragon Plus Environment

Page 19 of 31 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

342

decreased after methane and argon injection and before CO2 injection, as shown in Figure 7d. The

343

shrinkage of the nanopore structure for the Hazleton coal was consistent with a previous study on

344

anthracite6, where the deformation was related to the clay minerals. However, the macroscopic

345

coal deformation was reversible after helium and methane injection for coal samples from the

346

Illinois Basin12, which is different from the microscopic results in this study. Irreversible changes

347

in the pore structure did occur during the cyclic adsorption and desorption.37, 38 The results of the

348

two tested coals suggest that methane and argon penetration could permanently create nanopores

349

with a size of ~1 nm for the San Juan coal but could damage the nanopores with a size smaller

350

than 10 nm for the Hazleton coal.

351 352 353 354

Figure 7. The stability of total pore size distribution during gas injection for (a) San Juan coal and (b) Hazleton coal; The stability of total porosity before and after gas injection for (c) San Juan coal and (d) Hazleton coal.

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

355

3.3.2 Shale samples under gas injection

356

3.3.2.1 Argon injection

Page 20 of 31

357

Figure 8 shows the PSD and porosity evolution of the accessible pores under Ar injection

358

for the tested Marcellus shale samples. The 𝑓(𝑟) slightly decreased with increasing Ar pressure

359

for pores with a size greater than ~10 nm for the Marcellus shale drilled core sample, as shown in

360

Figure 8a, although there were small variations in the 𝑓(𝑟) as a function of pressure. The 𝑓(𝑟) for

361

pore sizes between ~10 and ~20 nm slightly decreased, and the 𝑓(𝑟) peaks at pore sizes of ~3 and

362

~6 nm shifted to smaller pore sizes (inserted figure in Figure 8a). Similar to the 𝑓(𝑟), the accessible

363

porosity, in general, slightly decreased with increasing pressure (Figure 8c). Thus, the results

364

indicate that the nanopores in the Marcellus shale drilled core were compressed by hydrostatic

365

argon injection. Compared to the coal samples, the compressive effect for the Marcellus shale

366

drilled core with an accessible porosity that was even twice as large was smaller compared to the

367

San Juan coal.

20 ACS Paragon Plus Environment

Page 21 of 31 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

368 369 370 371 372 373

Figure 8. Evolution of accessible pore size distribution under argon injection for (a) Marcellus drilled core shale and (b) Marcellus outcrop shale; Evolution of accessible porosity under argon injection for (c) Marcellus drilled core shale and (d) Marcellus outcrop shale.

374

increasing argon pressure at pore sizes greater than ~10 nm for the Marcellus shale outcrop sample,

375

as shown in Figure 8b. There were also small variations in the 𝑓(𝑟) as a function of the gas pressure.

376

The 𝑓(𝑟) peak at a pore size of ~7 nm decreased and slightly shifted to a smaller pore size. The

377

accessible porosity also showed a decrease with increasing argon pressure (Figure 8d), which was

378

consistent with the 𝑓(𝑟) results. Additionally, the decrease rate for the outcrop sample was higher

379

than that of the drilled core sample. Since the outcrop sample had a higher accessible porosity, the

380

results indicated that the shale samples with a higher accessible porosity experienced a greater

381

compressive effect by hydrostatic pressurization. This phenomenon is similar to the results of the

Similar to the 𝑓(𝑟) of the shale core sample, the 𝑓(𝑟) also slightly decreased with

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

Page 22 of 31

382

two coal samples described in the previous section. In addition, weathering may reduce the shale

383

strength, and therefore, the outcrop sample with a significantly high quartz content (> 75 %)

384

experienced a greater degree of compression compared to the drilled core sample that was obtained

385

from deep in the underground.

386

3.3.2.2 Methane injection

387

Figure 9 shows the PSD and porosity evolution of the accessible pores under methane

388

injection for the two shale samples. The 𝑓(𝑟) slightly decreased with increasing methane pressure

389

for pore sizes greater than ~6 nm for the shale core sample (Figure 9a). The 𝑓(𝑟) peak at a pore

390

size of ~5 nm moved up first and then shifted left with increasing pressure, and another peak at a

391

pore size of ~3 nm moved down first and then shifted left with increasing pressure (insert figure

392

in Figure 9a). However, the accessible porosity first decreased at 20 bar and then remained constant

393

with increasing pressure up to a pressure of 68 bar, as shown in Figure 9c. This result suggests that

394

there may have been very little sorption-induced matrix swelling up to 40 bar and a negligible

395

swelling effect at pressures > 40 bar for the Marcellus drilled core sample. The sorption-induced

396

swelling for the shale was approximately one magnitude lower than that for the coal.18 This

397

phenomenon could be due to the low TOC content (~2 %) of the shale drilled core sample

398

compared to the San Juan coal (~64 %), even though the accessible porosity of the drilled core

399

sample was two times greater than that of the San Juan coal.

22 ACS Paragon Plus Environment

Page 23 of 31 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

400 401 402 403 404 405

Figure 9. Evolution of accessible pore size distribution under CD4 injection for (a) Marcellus drilled core shale and (b) Marcellus outcrop shale; Evolution of accessible porosity under CD4 injection for (c) Marcellus drilled core shale and (d) Marcellus outcrop shale.

406

increasing methane pressure when the pore size was greater than ~6 nm, as shown in Figure 9b.

407

The 𝑓(𝑟) peak at a pore size of ~4 nm shifted to a smaller pore size with increasing pressure.

408

Specifically, we found that a significant decrease in the 𝑓(𝑟) occurred in the pore size range

409

between 6 and 9 nm, where the 𝑓(𝑟) may have even vanished when the pressure was above 40 bar.

410

This result suggests that methane densification could occur in the nanopores within the pore size

411

range between 6 and 9 nm for the shale outcrop sample. The accessible porosity generally

412

decreased with increasing pressure, as shown in Figure 9d, which is different from the results of

413

the drilled core sample. We also found that the decrease rate of the accessible porosity as a function

However, for the Marcellus shale outcrop samples, the 𝑓(𝑟) obviously decreased with

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 31

414

of the gas pressure was higher for methane injection compared to argon injection at relatively

415

lower pressures for the outcrop sample. The latter indicates that there was not only hydrostatic

416

compression but also a combination of adsorption-induced matrix swelling and pore filling (in the

417

pore size range between 6 and 9 nm) during methane injection for the Marcellus shale outcrop

418

sample. These results are similar to those of the San Juan coal sample.

419

3.3.2.3 CO2 injection

420

Figure 10 shows the PSD and porosity evolution of the accessible pores under CO2

421

injection for the tested shale samples. As depicted in Figure 10a, the 𝑓(𝑟) obviously decreased

422

with increasing CO2 pressure for pore sizes greater than ~10 nm for the Marcellus shale drilled

423

core sample. However, large variations were observed for the 𝑓(𝑟) peaks at pore sizes of ~3 and

424

~5 nm with increasing pressure (inserted figure in Figure 10a). The accessible porosity decreased

425

with increasing CO2 pressure initially below 40 bar and slightly increased at a pressure of 68 bar,

426

as shown in Figure 10c. The result is different from the supercritical CO2 swelling of shale, where

427

the degree of swelling first increased and then decreased with increasing pressure.17 The decrease

428

rate was higher during CO2 injection compared to methane injection, indicating a higher degree of

429

CO2 sorption-induced matrix swelling compared to methane for the Marcellus shale drilled core

430

sample. Notably, CO2 became liquid at 68 bar at room temperature. Condensation of liquid CO2

431

could occur at pore sizes between ~3 and ~4 nm, as shown in the inserted figure in Figure 10a,

432

because the 𝑓(𝑟) peak shifted to pore sizes smaller than ~3 nm and the number of nanopores

433

decreased.

24 ACS Paragon Plus Environment

Page 25 of 31 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

434 435 436 437 438 439

Figure 10. Evolution of accessible pore size distribution under CO2 injection for (a) Marcellus drilled core shale and (b) Marcellus outcrop shale; Evolution of accessible porosity under CO2 injection for (c) Marcellus drilled core shale and (d) Marcellus outcrop shale.

440

pore sizes greater than ~10 nm, though there were variations in the 𝑓(𝑟) as a function of the

441

pressure at pore sizes < ~10 nm for the Marcellus shale outcrop sample (Figure 10b). As was

442

expected, the 𝑓(𝑟) showed a significant decrease when CO2 became liquid. Consistently, the

443

accessible porosity continuously decreased with increasing CO2 pressure, as shown in Figure 10d.

444

At 68 bar, the accessible porosity significantly decreased and became less than half of the

445

accessible porosity value at vacuum. The results indicate that the effect of CO2 sorption-induced

446

matrix swelling was more pronounced compared to methane for the outcrop sample. At 68 bar, the

The value of 𝑓(𝑟) only slightly decreased with increasing CO2 pressure below 40 bar for

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

Page 26 of 31

447

detectable pore size was greater than ~15 nm. Thus, the condensation of and pore filling with liquid

448

CO2 could occur when the pore size was smaller than ~15 nm for this typical shale outcrop.

449

Additionally, nearly half of the accessible pore volume in the pore size range between 3.8 and 82

450

nm could be filled with a combination of liquid and adsorbed CO2 for the shale outcrop sample.

451

In general, for the two shale samples, CO2 had a profound sorption-induced swelling effect

452

compared to methane adsorption, which is similar to the coal results.9

453

3.3.2.4 Stability

454

Figure 11 shows the PSD and porosity of the total pores before and after gas injection for

455

the tested shale samples. There were variations in the 𝑓(𝑟) of the total pores initially, after methane

456

and argon injection, and before CO2 injection for the Marcellus shale drilled core sample, as shown

457

in Figure 11a. Specifically, there may have been a slight increase in the 𝑓(𝑟) when the pore size

458

was greater than ~30 nm for the shale drilled core sample (Figure 11a). There were fluctuations in

459

the 𝑓(𝑟) in the pore size range between ~5 and ~30 nm. Additionally, the 𝑓(𝑟) peak at a pore size

460

of ~1 nm slightly shifted to a smaller size (inserted figure in Figure 11a). As shown in Figure 11c,

461

the total porosity slightly increased after methane and argon injection and before CO2 injection.

462

We observed that the increased porosity could occur for a pore size greater than ~10 nm (Figure

463

11a). Relatively larger variations in the 𝑓(𝑟) were observed when the pore size was smaller than

464

~20 nm as was found for the Marcellus shale outcrop sample (Figure 11b). The total porosity also

465

exhibited larger variations, which decreased after methane and argon injection and before CO2

466

injection, as shown in Figure 11d. The decreased total porosity may be associated with the

467

decreased 𝑓(𝑟) for certain pore sizes: between ~5 and 10 nm, at ~1.5 nm, and at ~0.5 nm.

26 ACS Paragon Plus Environment

Page 27 of 31 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

468 469 470 471 472 473

Figure 11. The stability of total pore size distribution during gas injection for (a) Marcellus drilled core shale and (b) Marcellus outcrop shale; The stability of total porosity before and after gas injection for (c) Marcellus drilled core shale and (d) Marcellus outcrop shale.

474

SANS measurements were conducted to quantify the PSD and porosity of the total and

475

accessible pores for San Juan coal, Hazleton coal, Marcellus shale drilled core and outcrop samples.

476

The in situ SANS measurements were conducted to investigate the evolution of the PSD and

477

porosity of the accessible pores under various pressurization and adsorption environments. Based

478

on the characterization of the accessible pore structure evolution, several conclusions can be drawn

479

below:

480

(1) The total porosity ranges between 0.25 and 5.8 % at the detected pore size range (< ~82 nm)

481

for the tested coal and shale samples. Less than 50 % of the total porosity is accessible to

4. Conclusions

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

Page 28 of 31

482

methane for the two coals compared to greater than 75 % of the fraction of accessible porosity

483

for the two shales, indicating that the organic matter pores tend to be disconnected in this

484

nanopore range.

485

(2) For the two coal samples with different ranks, the hydrostatic compression has a higher effect

486

on coal with a higher accessible porosity and with a lower rank. For the San Juan

487

subbituminous coal, the nanopores with pore sizes < ~20 nm are affected by a combination of

488

compression and sorption-induced matrix swelling during the methane injection. Methane

489

densification and pore filling could occur above 68 bar. CO2 has a more profound effect on

490

matrix swelling, gas densification and pore filling for pore sizes smaller than ~7 nm, which

491

could occur above 40 bar. For the Hazleton semi-anthracite coal, there is no obvious effect of

492

matrix swelling, but the fraction of accessible pores could slightly increase with increasing

493

pressure during methane and CO2 injection. Methane and argon penetration could permanently

494

create nanopores at pore sizes of ~1 nm for the San Juan coal but could damage the nanopores

495

with a pore size smaller than 10 nm for the Hazleton coal.

496

(3) For the two Marcellus shale samples, the drilled core sample has a lower effect of compression

497

compared to the outcrop sample with a higher accessible porosity. The effect of methane

498

sorption-induced matrix swelling is small below 20 bar and negligible above 20 bar for the

499

Marcellus shale drilled core. However, there is notable matrix swelling for the Marcellus shale

500

outcrop sample, where densification and pore filling could occur in the nanopores with a size

501

ranging between 6 and 9 nm. Similar to the coal samples, CO2 also has a greater effect on

502

sorption-induced matrix swelling for the two shale samples. Densification of and pore filling

503

with CO2 occurs in pores with a size < ~15 nm at 68 bar for the outcrop sample. In addition, it

28 ACS Paragon Plus Environment

Page 29 of 31 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

504

is found that the total porosity increases for Marcellus shale drilled core but decreases for the

505

Marcellus shale outcrop after methane and argon injection.

506

Acknowledgments

507

We want to thank Jitendra Bahadur, Lilin He and the late Yuri Melnichenko for SANS

508

measurement and data reduction. This research used resources at the High Flux Isotope Reactor, a

509

DOE Office of Science User Facility operated by the Oak Ridge National Laboratory. This

510

material is, in part, based upon work supported by the U.S. Department of Energy, Office of

511

Science, Office of Basic Energy Sciences, under contract DE-SC0006883. This paper describes

512

objective technical results and analysis. Any subjective views or opinions that might be expressed

513

in the paper do not necessarily represent the views of the U.S. Department of Energy or the United

514

States Government. Sandia National Laboratories is a multimission laboratory managed and

515

operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned

516

subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear

517

Security Administration under contract DE-NA-0003525. The authors declare no competing

518

financial interests.

519

References

520 521 522 523 524 525 526 527 528 529 530 531 532 533

(1) Outlook, A. E., Energy information administration. Department of Energy 2016. (2) Karacan, C. O., Swelling-induced volumetric strains internal to a stressed coal associated with CO2 sorption. International Journal of Coal Geology 2007, 72, (3-4), 209-220. (3) Mao, L.; Hao, N.; An, L.; Chiang, F.-p.; Liu, H., 3D mapping of carbon dioxide-induced strain in coal using digital volumetric speckle photography technique and X-ray computer tomography. International Journal of Coal Geology 2015, 147, 115-125. (4) Karacan, C. O., Heterogeneous sorption and swelling in a confined and stressed coal during CO2 injection. Energy & Fuels 2003, 17, (6), 1595-1608. (5) Toishi, A.; Yamazaki, Y.; Uchida, A.; Saito, Y.; Aoki, H.; Nomura, S.; Arima, T.; Kubota, Y.; Hayashizaki, H., Quantitative Evaluation for Relationship between Development of Pore Structure and Swelling of Coal in Carbonization of Single Coal Particle. Isij International 2013, 53, (10), 1739-1748. (6) Feng, Z.; Zhou, D.; Zhao, Y.; Cai, T., Study on microstructural changes of coal after methane adsorption. Journal of Natural Gas Science and Engineering 2016, 30, 28-37. 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

534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578

Page 30 of 31

(7) Brochard, L.; Vandamme, M.; Pelenq, R. J. M.; Fen-Chong, T., Adsorption-Induced Deformation of Microporous Materials: Coal Swelling Induced by CO2-CH4 Competitive Adsorption. Langmuir 2012, 28, (5), 2659-2670. (8) Zang, J.; Wang, K.; Zhao, Y., Evaluation of gas sorption-induced internal swelling in coal. Fuel 2015, 143, 165-172. (9) Day, S.; Fry, R.; Sakurovs, R., Swelling of moist coal in carbon dioxide and methane. International Journal of Coal Geology 2011, 86, (2-3), 197-203. (10) Staib, G.; Sakurovs, R.; Gray, E. M., Kinetics of coal swelling in gases: Influence of gas pressure, gas type and coal type. International Journal of Coal Geology 2014, 132, 117-122. (11) Fry, R.; Day, S.; Sakurovs, R., Moisture-induced swelling of coal. International Journal of Coal Preparation and Utilization 2009, 29, (6), 298-316. (12) Liu, S.; Wang, Y.; Harpalani, S., Anisotropy characteristics of coal shrinkage/swelling and its impact on coal permeability evolution with CO2 injection. Greenhouse Gases: Science and Technology 2016. (13) Battistutta, E.; van Hemert, P.; Lutynski, M.; Bruining, H.; Wolf, K.-H., Swelling and sorption experiments on methane, nitrogen and carbon dioxide on dry Selar Cornish coal. International Journal of Coal Geology 2010, 84, (1), 39-48. (14) Yang, K.; Lu, X.; Lin, Y.; Neimark, A. V., Effects of CO2 adsorption on coal deformation during geological sequestration. J Geophys Res-Sol Ea 2011, 116. (15) Yang, K.; Lu, X.; Lin, Y.; Neimark, A. V., Deformation of Coal Induced by Methane Adsorption at Geological Conditions. Energy & Fuels 2010, 24, 5955-5964. (16) Zhang, W.; Hu, H.; Li, X.; Fang, Z., Interplay of Montmorillonite-H2O-scCO(2) System between Mechanical Behavior and Adsorption: Molecular Dynamics. Journal of Physical Chemistry C 2015, 119, (38), 21959-21968. (17) Lu, Y.; Ao, X.; Tang, J.; Jia, Y.; Zhang, X.; Chen, Y., Swelling of shale in supercritical carbon dioxide. Journal of Natural Gas Science and Engineering 2016, 30, 268-275. (18) Chen, T.; Feng, X.-T.; Pan, Z., Experimental study of swelling of organic rich shale in methane. International Journal of Coal Geology 2015, 150, 64-73. (19) Wang, J. G.; Ju, Y.; Gao, F.; Peng, Y.; Gao, Y., Effect of CO2 sorption-induced anisotropic swelling on caprock sealing efficiency. Journal of Cleaner Production 2015, 103, 685-695. (20) Zhang, R.; Liu, S.; Wang, Y., Fractal evolution under in situ pressure and sorption conditions for coal and shale. Scientific Reports 2017, 7, (1), 8971. (21) Zhang, R.; Liu, S.; Bahadur, J.; Elsworth, D.; Melnichenko, Y.; He, L.; Wang, Y., Estimation and modeling of coal pore accessibility using small angle neutron scattering. Fuel 2015, 161, 323332. (22) Wignall, G. D.; Littrell, K. C.; Heller, W. T.; Melnichenko, Y. B.; Bailey, K. M.; Lynn, G. W.; Myles, D. A.; Urban, V. S.; Buchanan, M. V.; Selby, D. L.; Butler, P. D., The 40 m general purpose small-angle neutron scattering instrument at Oak Ridge National Laboratory. Journal of Applied Crystallography 2012, 45, 990-998. (23) He, L.; Melnichenko, Y. B.; Mastalerz, M.; Sakurovs, R.; Radlinski, A. P.; Blach, T., Pore accessibility by methane and carbon dioxide in coal as determined by neutron scattering. Energy & Fuels 2012, 26, (3), 1975-1983. (24) Radlinski, A. P.; Boreham, C. J.; Lindner, P.; Randl, O.; Wignall, G. D.; Hinde, A.; Hope, J. M., Small angle neutron scattering signature of oil generation in artificially and naturally matured hydrocarbon source rocks. Organic Geochemistry 2000, 31, (1), 1-14.

30 ACS Paragon Plus Environment

Page 31 of 31 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

579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614

Energy & Fuels

(25) Kline, S. R., Reduction and analysis of SANS and USANS data using IGOR Pro. Journal of Applied Crystallography 2006, 39, 895-900. (26) Wignall, G. D.; Bates, F. S., Absolute calibration of small-angle neutron scattering data. Journal of Applied Crystallography 1987, 20, 28-40. (27) Liu, S. M.; Harpalani, S., Compressibility of sorptive porous media: Part 2. Experimental study on coal. Aapg Bulletin 2014, 98, (9), 1773-1788. (28) Ilavsky, J.; Jemian, P. R., Irena: tool suite for modeling and analysis of small-angle scattering. Journal of Applied Crystallography 2009, 42, 347-353. (29) Radlinski, A. P., Small-angle neutron scattering and the microstructure of rocks. Rev Mineral Geochem 2006, 63, 363-397. (30) Gu, X.; Mildner, D. F. R.; Cole, D. R.; Rother, G.; Slingerland, R.; Brantley, S. L., Quantification of Organic Porosity and Water Accessibility in Marcellus Shale Using Neutron Scattering. Energy & Fuels 2016, 30, (6), 4438-4449. (31) Bahadur, J.; Medina, C. R.; He, L. L.; Melnichenko, Y. B.; Rupp, J. A.; Blach, T. P.; Mildner, D. F. R., Determination of closed porosity in rocks by small-angle neutron scattering. Journal of Applied Crystallography 2016, 49, 2021-2030. (32) Melnichenko, Y. B.; He, L.; Sakurovs, R.; Kholodenko, A. L.; Blach, T.; Mastalerz, M.; Radlinski, A. P.; Cheng, G.; Mildner, D. F. R., Accessibility of pores in coal to methane and carbon dioxide. Fuel 2012, 91, (1), 200-208. (33) 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, (2), 772-779. (34) Glatter, O.; Kratky, O., Small angle X-ray scattering. Academic Press: 1982. (35) Sakurovs, R.; Radlinski, A. P.; Melnichenko, Y. B.; Blach, T.; Cheng, G.; Lemmel, H.; Rauch, H., Stability of the Bituminous Coal Microstructure upon Exposure to High Pressures of Helium. Energy & Fuels 2009, 23, 5022-5026. (36) Zhang, R.; Liu, S., Experimental and theoretical characterization of methane and CO2 sorption hysteresis in coals based on Langmuir desorption. International Journal of Coal Geology 2017, 171, 49-60. (37) Zhou, D.; Feng, Z. C.; Zhao, D.; Zhao, Y. S.; Cai, T. T., Experimental study of mesostructural deformation of coal during methane adsorption-desorption cycles. Journal of Natural Gas Science and Engineering 2017, 42, 243-251. (38) Majewska, Z.; Majewski, S.; Zietek, J., Swelling of coal induced by cyclic sorption/desorption of gas: Experimental observations indicating changes in coal structure due to sorption of CO2 and CH4. International Journal of Coal Geology 2010, 83, (4), 475-483.

31 ACS Paragon Plus Environment