CO2 and CH4 Wettabilities of Organic-Rich Shale - Energy & Fuels

Nov 22, 2017 - School of Petroleum Engineering, China University of Petroleum (East China), No. 66, Changjiang West Road, Qingdao 266580, ... Experime...
0 downloads 6 Views 1MB Size
Subscriber access provided by READING UNIV

Article

CO2 and CH4 wettabilities of organic-rich shale Bin Pan, Yajun Li, Hongqian Wang, Franca Jones, and Stefan Iglauer Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b01147 • Publication Date (Web): 22 Nov 2017 Downloaded from http://pubs.acs.org on December 2, 2017

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 39

180 160 140 298 K, Advancing, CO₂

120

Contact angle [⁰]

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

298 K, Receding, CO₂ 298 K, Advancing, CH₄

100

298 K, Receding, CH₄ 323 K, Advancing, CO₂ 323 K, Advancing, CO₂

80

323 K, Advancing, CH₄ 323 K, Receding, CH₄

60

343 K, Advancing, CO₂ 343 K, Receding, CO₂

40

343 K, Advancing, CH₄ 343 K, Receding, CH₄

20 0

5

10

15 Pressure [MPa]

20

ACS Paragon Plus Environment

25

30

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

CO2 and CH4 wettabilities of organic-rich shale

1

2

Page 2 of 39

Bin Pana,b,c, Yajun Lia*, Hongqian Wanga, Franca Jonesd, Stefan Iglauerc

3 4

a

5 6

China University of Petroleum (East China), School of Petroleum Engineering, No. 66,Changjiang West Road, Qingdao, China , 266580

b Subsurface Fluidics and

7

Porous Media Laboratory, University of Calgary, Department of Chemical and Petroleum

Engineering, AB, Canada, T2N 1N4

8

c

Curtin University, Department of Petroleum Engineering, 26 Dick Perry Avenue, 6151 Kensington, Australia

9

d

Curtin University, Department of Chemistry, Nanochemistry Research Institute, PO Box U1987, Perth WA, 6845

10

Abstract

11

CO2 and CH4 wettabilities of organic-rich shale are important physico-chemical

12

parameters which significantly influence CO2 sequestration and CH4 production.

13

However, there is a serious lack of understanding of these aspects as the data available

14

is scarce. Thus, we evaluated organic-rich shale CO2- and CH4-wettabilities (i.e.

15

brine/shale/gas systems) through advancing and receding brine contact angle

16

measurements as a function of pressure, temperature, salinity and ion type (as these can

17

vary significantly in underground formations). The results indicated that the brine

18

contact angles for both CO2/CH4-brine-shale systems increased with pressure and

19

salinity, but decreased with temperature. However, these effects were much less

20

significant for CH4. Furthermore, the brine contact angles for CO2-brine-shale system

21

reached 180° (i.e. the shale was completely wetted by CO2) when the pressure reached 1

ACS Paragon Plus Environment

Page 3 of 39

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

22

30 MPa at 343 K and 9 MPa at 298 K. The brine contact angles for the analogue CH4

23

systems was much lower, 50-90°, only indicating weakly water-wet to intermediate-

24

wet conditions. Finally, the brine contact angles for CO2-brine-shale system were also

25

larger for divalent ions (Ca2+, Mg2+) than for monovalent ions (Na+, K+), while ion type

26

had no significant influence on CH4 wettability. However, a similar CO2/CH4 density

27

resulted in a similar wettability. Consequently CH4 could not be used as a proxy for

28

predicting CO2 storage capacities unless they have similar densities.

29 30

Key words: CO2 wettability, CH4 wettability, organic-rich shale.

31

32

1. INTRODUCTION

33

34

Global warming and the greenhouse effect are primarily caused by anthropogenic

35

CO2 emissions. To reduce these emissions, carbon geo-sequestration (CGS) has been

36

identified as an important solution.1,2 In CGS, CO2 is captured from large-point-source

37

emitters, purified, compressed, and injected into geological formations at depths of or

38

greater than 800 m, including saline aquifers,3-5 (depleted) oil and gas reservoirs,4,6-8 or 2

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 39

39

un-mineable coal seams,4,9 where CO2 is stored in a supercritical (sc) state.1 However,

40

the scCO2 density is, approximately, only 50%-80% of formation brine density, which

41

leads to CO2 upwards migration.10 Four main sequestration mechanisms in CGS are: (a)

42

structural trapping,11,12 (b) residual trapping,13-15 (c) dissolution trapping,5,16,17 and (d)

43

mineral trapping.18,19 Another way is to store CO2 in clay-rich formations by CO2

44

adsorption trapping.20-22 In CCS, the advancing brine angle is related to residual

45

trapping, while the receding brine angle is related to structural trapping.23

46

In addition, due to the improvement of hydraulic fracturing and horizontal drilling

47

technologies, shale gas has been recognized as a promising energy resource, which

48

can be exploited commercially.24-26 Shale gas primarily contains methane, which

49

generally exists in organic-rich pores of low permeability shale formations.27,28 Such

50

shale gas recovery can potentially be improved by CO2 injection, while

51

simultaneously CO2 is sequestered (adsorption trapping).22,29-36 In such a scenario,

52

CO2- and CH4 rock wettabilities play very important roles as they determine fluid

53

spreading and fluid dynamics in the reservoir.12,13,30,37-39

54

However, organic content in shale can vary dramatically,40-42 and thus the

55

coverage of shale surfaces with organics is expected to vary substantially. Moreover,

56

only a few research papers reported CO2 wettability of shale4,22,37,43-45, and available

57

data is scarce.

58 59

Furthermore, some scholars used CH4 as a proxy for predicting structural CO2 storage capacities, assuming that CO2 and CH4 wettabilities are similar.46-48 This

3

ACS Paragon Plus Environment

Page 5 of 39

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

60

assumption, however, has not yet been systematically verified, and it is clear that CO2

61

and CH4 have fundamentally different physical and molecular properties.23 Moreover,

62

as the best of our knowledge, the brine contact angles of CH4-brine-shale system have

63

not been reported, although they are of key importance in shale gas production, see

64

above.37-39

65

Thus, we measured the advancing and receding brine contact angles of CO2/CH4-

66

brine-shale systems as a function of pressure, temperature, salinity and ion type in

67

order to aid shale gas production and CO2 geo-storage capacity predictions.

68

69

2. EXPERIMENTAL SECTION

70

71

2.1. Materials.

72

CO2 (purity = 99.7 mol%), CH4 (purity = 99.7 mol%) and de-gassed brines (0 wt%,

73

1 wt%, and 5 wt% NaCl, CaCl2 and MgCl2 dissolved in deionized water; stirred for 2

74

hr and vaccumed for 12 hr at ambient temperature) were used in all experiments; This

75

brine composition was chosen to simulate Shengli shale formation water (which mainly

76

contains Na+, Ca2+ and Mg2+ cations,49,50) and the salinities range with formation

77

depth.51

78

The shale sample examined was acquired from well FY-1 in the Shengli shale play

79

formation, Jiyang depression, China. Geological and geo-chemical properties of the

4

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 39

80

shale sample are listed in Table 1. The shale sample was thoroughly characterized by

81

quantitative X-ray diffraction (XRD), total organic content (TOC), scanning electron

82

microscopy-energy dispersive X-ray spectroscopy (SEM-EDS) and surface roughness

83

tests (Tables 1 and 2 and Figure 1).

84

85

Table 1. Geological and geo-chemical properties of the shale sample examined. Lithology

laminar-mixed limestone and mudstone

Depth (m)

3178

TOC (mg/kg)

30000

Vitrinite reflectance (%)

0.76

Matrix density (g/cm3)

2.73

Porositya (%)

8.21

Average pore throat radiusb (nm) Average Permeabilityc (mD)

5 10-6

Chemical compositiond (wt%) Calcite

49

Quartz

19

Illite

16

Ankerite

10

Pyrite

3 5

ACS Paragon Plus Environment

Page 7 of 39

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

Albite a

Porosity measured by NMR,52 using MacroMR12-110H-1 at ambient conditions.

b

c

3

Pore throat radius measured by gas adsorption and mercury intrusion.53

Permeability measured by gas pulse method.54

d

Chemical composition measured with a Bruker AXS XRD instrument.

86

87

Table 2. EDS analysis of the shale sample. Element

Atomic Concentration (%)

Oxygen

53.9

Silicon

13.2

Carbon

28.1

Aluminum

4.0

Magnesium

0.8

88 89

90

(b)

(a)

6

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 39

91

92 93

(c)

94 95

(d)

96 97

(e) 7

ACS Paragon Plus Environment

Page 9 of 39

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

98

Figure 1. SEM and atomic force microscopy images of the shale sample used in the

99

experiments. (a-c) SEM images; (d) 3D topography; (e) deflection signal, (d) and (e):

100

different heights are shown in different colours (black is 0 nm, white is 4000 nm).

101

102

The shale sample consisted mainly of calcite (49 wt%), quartz (19 wt%) and illite

103

(16 wt%) and a few other minor minerals, Table 1. The SEM images, Figure 1a-c,

104

showed that organic matter filled some of the pores; microfractures and inorganic pores

105

were also identified. The EDS analysis in Table 2 showed that the total atomic carbon

106

concentration on the surface was 28.1 mol%. Surface topography of the sample was

107

measured via atomic force microscopy (instrument model AFM DSE 95-200), Figure

108

1d,e, and the sample was relatively rough when compared to single cystals (e.g. a

109

smooth alpha-quartz had a RMS surface roughness of 40 nm)55 and a root-mean-square

110

(RMS) surface roughness of 490 nm was measured.

111

112

2.2. Contact angle measurement.

113

For the contact angle experiments, a cuboid sample was cut with a high speed

114

diamond blade (to 17 mm x 12 mm x 3 mm dimensions). The sample was then cleaned

115

with acetone and exposed to air plasma for 3 min to remove any organic surface

116

contaminants.56 Subsequently, the advancing brine contact angle (𝜃𝑎 ) and the receding

117

brine contact angle (𝜃𝑟 ) were measured by a Krüss DSA 100 instrument, using the tilted

118

plate method.57 The cleaned substrate was placed into the pressure cell at a pre-set 8

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 39

119

temperature (298 K, 323 K or 343 K) and the system was flooded with CO2 or CH4 for

120

10 min at ambient pressure. Subsequently the gas pressure was increased to prescribed

121

values (0.1 MPa, 5 MPa, 10 MPa, 20 MPa and 30 MPa). After pressure stabilization

122

was achieved, a droplet (average volume of ~6 ±1 µL) of brine was dispensed onto the

123

surface of the shale substrate through a needle. 𝜃𝑎 was measured just before the

124

leading edge of the droplet started to move and 𝜃𝑟 was measured simultaneously at

125

the trailing edge of the droplet. A high resolution video camera recorded the dispensing

126

process as a movie, and 𝜃𝑎 and 𝜃𝑟 were measured on the images extracted from the

127

movies. To obtain reliable and reproducible results, all experiments were conducted

128

with vacuumed brines and repeated 3 times (All the experimental data were listed in

129

SI-Table 1 Supporting Information). Brines were not equilibrated with CO2 or CH4, as

130

it was shown earlier that this has only an insignificant influence on θ during the θ

131

measurement.55,58 The standard deviation of these measurements was ± 3° based on

132

replicate measurements.

133

134

2.3. Zeta potential measurement.

135

For the zeta potential experiments, a Zetasizer Nano instrument (Malvern, UK) was

136

used to measure the zeta potentials of a powder suspension prepared from the examined

137

shale sample. The shale suspension was prepared as follows: the shale rock sample was

138

crushed to very fine particles of less than 10 micrometer size. 0.2 wt% shale powder

139

was then added to NaCl, CaCl2 and MgCl2 brines, respectively (0.1 g shale powder into

9

ACS Paragon Plus Environment

Page 11 of 39

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

140

50 mL brine). Each dispersion was mixed with a magnetic stirrer for 15 min at a speed

141

of 4×100 rpm at 323 K. The solution was subsequently equilibrated for 12 hr. and

142

transferred to a measurement cell. The cell temperature was set to 323 K (at ambient

143

pressure), and the zeta potential was measured at least ten times for each solution and

144

the arithmetic average was taken. The standard deviation of these measurements was

145

±1 mV based on replicate measurements. The pH values ranged from 6.72-8.81

146

(measured by RHA-SHKY pH Meter Model 8601), Table 3.

147 148 149

Table 3. The pH values of different solutions at 323 K and ambient pressure. (Note that 0.2 wt% shale powder was added to each solution)

150

wt% NaCl

pH

wt% CaCl2

pH

wt% MgCl2

pH

1

8.81

1

7.4

1

8.06

2

8.55

2.5

6.99

2.5

7.89

5

8.19

5

6.72

5

7.62

151

152

3. RESULTS AND DISCUSSION

153

154

3.1. Influence of pressure.

155

As pressure and temperature are expected to vary with storage depth,59 both, the

156

advancing (𝜃𝑎 ) and the receding (𝜃𝑟 ) brine contact angles of the CO2/CH4-brine-shale 10

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

157

158

Page 12 of 39

systems were measured at various pressures and temperatures, Figure 2. For the CO2-brine-shale system, the brine contact angles of 𝜃𝑎 and 𝜃𝑟 at

159

ambient conditions (T = 298 K and 0.1 MPa) were measured as 73° and 72°,

160

respectively, which indicated intermediate-wet conditions. Note all the relationships

161

between contact angles and wet state in this work were based on literature standard.12

162

However, at 20 MPa and 343 K, 𝜃𝑎 and 𝜃𝑟 were measured as 130° and 120°,

163

respectively, indicating weakly CO2-wet conditions. At even higher pressure, 30MPa,

164

all brine contact angles reached 180°, i.e. a completely CO2-wet state. In this case, the

165

brine droplet did not stay on the shale substrate and it rolled away immediately once it

166

dropped onto the tilted surface of the shale substrate.

167

On the contrary, for the CH4-brine-shale system, 𝜃𝑎 and 𝜃𝑟 at ambient

168

conditions (T = 298 K and 0.1 MPa) were 54° and 46°, respectively, which indicated

169

strongly to weakly water-wet conditions. 𝜃𝑎 and 𝜃𝑟 at 20 MPa and 343 K were 74°

170

and 70°, respectively (an increase by 42° and 38° when the pressure increased from

171

0.1 MPa to 20 MPa), and this indicated intermediate-wet conditions.

172

Although, the brine contact angles of both, CO2 and CH4, increased with

173

pressure (due to lower density differences between fluids,23,60,61 and stronger

174

intermolecular interactions between CO2/CH4 and the rock surface,23,61-64), consistent

175

with literature data for CO2 on quartz,55,56,61,65,66 mica,65,67 shale,4,44,68 sandstone,44

176

coal,69 and for CH4 on coal,70,71 the CH4–shale system was far more water-wet than

177

the analogue CO2 system, consistent with data measured on coal.70,71 To further

11

ACS Paragon Plus Environment

Page 13 of 39

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

178

explore and quantify the reason why the brine contact angles for CO2 and CH4 at the

179

same thermo-physcial conditions were different, we plotted advancing and receding

180

brine contact angles versus the CO2 and CH4 densities, respectively, Figure 3. When

181

the CO2 and CH4 densities were close, the brine contact angles matched. Specifically,

182

when CO2 density was 38 kg/m3, its advancing and receding brine contact angles were

183

70° and 60°, respectively; when CH4 density was 35 kg/m3, its advancing and

184

receding brine contact angle were 65° and 59°. When CO2 and CH4 density reached

185

133 and 135 kg/m3, the CO2- and CH4- advancing brine contact angles were 90° and

186

91°, respectively. These results clearly showed that when the fluid densities were

187

similar, their wettabilities were also similar, consistent with measurements for N2, Ar,

188

SF6 and He on quartz.66 Extrapolations from CH4 behavior to CO2 behavior thus need

189

to be avoided unless they have similar densities.

190

12

ACS Paragon Plus Environment

Energy & Fuels

180 160 140 298 K, Advancing, CO₂

Contact angle [⁰]

120

298 K, Receding, CO₂ 298 K, Advancing, CH₄

100

298 K, Receding, CH₄ 323 K, Advancing, CO₂ 323 K, Advancing, CO₂

80

323 K, Advancing, CH₄ 323 K, Receding, CH₄

60

343 K, Advancing, CO₂ 343 K, Receding, CO₂

40

343 K, Advancing, CH₄ 343 K, Receding, CH₄

20 0

5

10

15 Pressure [MPa]

191

20

25

30

192

Figure 2. Effect of pressure and temperature on the brine contact angles of the

193

CO2/CH4-5 wt% NaCl brine-shale system. 180 160 140

Contact angle [⁰]

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 39

120 100

Advancing, CO₂

80

Receding, CO₂

60

Advancing, CH₄

40

Receding, CH₄

20 0 0

100

200

300

400

500

Density

600

700

800

900

1000

[kg/m3]

194 195

Figure 3. Effect of fluid density on the brine contact angles of the CO2/CH4-5 wt%

196

NaCl brine-shale system (CO2 and CH4 density were taken from literature data.72,73) 13

ACS Paragon Plus Environment

Page 15 of 39

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

197

198

3.2. Influence of temperature.

199

𝜃𝑎 and 𝜃𝑟 decreased with temperature, Figure 2. This is consistent with

200

molecular dynamics simulations for CO2 on quartz,62 and some experimental data for

201

CO2 on quartz,61,66 and mica,67,74 where 𝜃𝑎 and 𝜃𝑟 decreased with temperature.

202

Specifically, for the CO2-brine-shale system, when temperature increased from 298 K

203

to 343 K at 10 MPa pressure, 𝜃𝑎 decreased by 60° from 180° to 120° and 𝜃𝑟

204

decreased by 80° from 180° to 100°, thus changing from completely CO2-wet to weakly

205

CO2-wet. For the CH4-brine-shale system, when temperature increased from 298 K to

206

343 K at 10 MPa pressure, 𝜃𝑎 decreased by 12° from 79° to 67°, and 𝜃𝑟 decreased

207

by 12° from 65° to 53°, thus the shale substrate became more water-wet. Note that the

208

temperature had a stronger influence on CO2 wettability than on CH4 wettability.

209

However, other researchers found an opposite trend for the effect of temperature on

210

CO2 wettability,37,69,75,76 or no temperature influence.77 A recent study reported that

211

temperature influences contact angles through density and Van der Waals forces.31 For

212

CO2-wet surfaces, Van der Waals forces (solid-CO2) decrease with increasing

213

temperature (causing CO2 to detach more from the surface, which renders the surface

214

more water-wet); however, hydrogen bonds are increasingly broken on water-wet

215

surfaces with increasing temperature (thus reducing the water-wetness of the surface).

216

Therefore, here in our work, a CO2-wet shale surface became more water-wet with

217

temperature increase. Furthermore, we did not find any previous report on the effect of

218

temperature on CH4 wettability. 14

ACS Paragon Plus Environment

Energy & Fuels

219

220

3.3. Influence of salinity and salt type.

221

The effect of salinity on the brine contact angles was investigated at 298 K and 10

222

MPa, Figure 4. Salinities tested ranged from (0 to 5) wt%.

180 Advancing, CO₂

160

Receding, CO₂ Advancing, CH₄

140

Contact angle [°]

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 39

Receding, CH₄

120 100 80 60 40 20 0

1

2

3

4

5

Salinity [% wt]

223 224

Figure 4. Effect of salinity on the advancing and receding brine contact angles of the

225

CO2 and CH4-NaCl brine-shale systems at 298 K and 10 MPa.

226

227

For the CO2-brine-shale system, 𝜃𝑎 increased by 90° and 𝜃𝑟 by 97° when NaCl

228

salinity increased from 0 to 5 wt% at 298 K and 10 MPa. For the CH4-brine-shale

229

system, 𝜃𝑎 increased by 15° and 𝜃𝑟 by 8° when NaCl salinity increased from 0 to 5

230

wt% at the same thermodynamic condition, consistent with literature data for CO2 on

231

quartz,61 mica,67,74,78 shale37 and other rock minerals.65,75,79 Thus salinity had a larger

232

influence on CO2 wettability than on CH4 wettability. Furthermore the brine contact 15

ACS Paragon Plus Environment

Page 17 of 39

233

angles for divalent cations were higher than those for monovalent cations (Figure 5),

234

following the order MgCl2 > CaCl2 > NaCl, consistent with literature data for CO2 on

235

quartz,61 mica,78 coal,80 and shale.37 However, these ion types only had an insignificant

236

influence (2°) on CH4 wettability (which is within experimental error). Note that

237

literature data for the effect of salinity and salt type on CH4 wettability is scarce. 110 100 90 80

Contact angle [°]

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

70 60

343 K, Advancing, CO₂

50

343 K, Receding, CO₂

40

298 K, Advancing, CH₄

30

298 K, Receding, CH₄

20 10 0 Na⁺

Ca²⁺

Mg²⁺

Ion type

238 239

Figure 5. Effect of salt type on the advancing and receding brine contact angles at 10

240

MPa and 1 wt% salt concentration.

241

16

ACS Paragon Plus Environment

Energy & Fuels

10

Zeta Potential [mV]

5 0 -5 -10

NaCl CaCl₂

-15

MgCl₂ -20 0

1

2

3

4

5

Concentration of different brines [wt%]

242

(a)

243

244

10 5

Zeta Potential [mV]

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 39

0 -5

NaCl

-10

CaCl₂ -15

MgCl₂

-20 0

0.5

1

1.5

Ionic strength [mol/kg]

245

(b)

246

247

Figure 6. Shale zeta potentials as a function of (a) salt concentration and (b) ionic

248

strength of the brine (in which the shale was dispersed) at 323 K. The error bars show

249

the standard error of the measurements. 17

ACS Paragon Plus Environment

Page 19 of 39

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

250

251

The effect of salinity and ion type on the brine contact angles can be explained

252

with the diffuse double layer theory with specific adsorption of ions included.37,81,82 It

253

is clear that at higher pH values more protons are released from the surface OH

254

groups, which results in a more negative surface charge (= more negative zeta

255

potential). However, we cannot measure zeta potentials at high pressures (due to

256

instrument limitations), thus we measured and reported zeta potentials at ambient

257

pressure (Figure 6). Note that efforts were made to correlate such (ambient condition)

258

zeta potential measurements to real reservoir conditions.80 Trendwise, the contact

259

angles followed the same trends with temperature and salinity at ambient and high

260

pressure, while the zeta potentials also followed the same trend at ambient pressure;

261

thus it was assumed that zeta potentials probably also follow a similar trend at high

262

pressures.

263

Furthermore, with the addition of NaCl, CaCl2 or MgCl2 into the shale suspension,

264

the zeta potential increased, Figure 6. Mechanistically, this increase in zeta potential is

265

caused by adsorption of divalent ions onto the shale surface, which leads to a more

266

positive charge.83,84 However, sodium ions are not so strongly adsorbed, so here the 85,86

1 n 2 As the ionic strength ( I   ci zi , where I is the ionic 2 i 1

267

effect is less significant.

268

strength, c is the concentration and z is the ion charge) increased, the surface charge

269

was more and more screened which resulted in more positive zeta potentials. This

270

charge screening effect is stronger for Mg2+ and Ca2+ than for Na+ (due to their higher 18

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 20 of 39

271

ionic strength).23 At very high ionic strengths charge reversal occurred for Mg2+ and

272

Ca2+, presumably due to specific adsorption of the divalent ions on the clay or the

273

calcitic shale surface.83,84 The ranking Mg2+ > Ca2+ based on the de-wetting ability and

274

zeta potential is probably caused by the stronger adsorption for Mg2+ than Ca2+ onto the

275

shale surface, consistent with other experimental results (where the zeta potential for

276

Mg2+ solutions was more positive than that for Ca2+ solutions).87 In addition, one study

277

investigated the effect of individual ionic concentration on zeta potential in calcite

278

solutions and found that Ca2+ can only change zeta potentials from negative to positive

279

values at high concentrations, while small additions of Mg2+ into the calcite solution

280

can readily achieve a positive zeta potential,88 consistent with our results. Mg2+ has a

281

stronger effect on calcite surface charge due to ionic adsorption88; and recall that the

282

calcite content in our shale sample was high (49 wt%). However, some other study

283

reported that Ca2+ had stronger adsorption ability than Mg2+ on shale,26 and note that

284

the primary composition in its shale sample was quartz.26 The discrepancy of divalent

285

ion adsorption capacity here probably was due to the composition of the shale sample.

286

Furthermore, this charge reversal has been previously associated with higher contact

287

angles.86,89,90 Moreover, once the zeta potential approaches zero, the surface polarity

288

decreases and the affinity of the shale surface for water decreases so that the brine

289

contact angle increases.23,89,91-93

290

291

3.4. Influence of total organic carbon content.

19

ACS Paragon Plus Environment

Page 21 of 39

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

292

Energy & Fuels

We furthermore investigated which influence the total organic carbon content (TOC)

293

has on the brine contact angle of brine-CO2-shale, Table 4. However, only a few

294

studies reported the CO2-wettability of shale, 22,37,43,44 thus data is scarce.

295

296

Table 4. The effect of TOC on shale wettability. Chemical

Mineral

TOC

NaCl

T

P

𝜃𝑟

𝜃𝑎 Reference

Composition

wt%

Calcite

mg/kg

MPa

°

°

49

0.1

60

70

Quartz

19

10

107

125

Illite

16

20

120

130

Ankerite

10

30

180

180

Other

6

Illite

33

0.1

20

30

Quartz

31

10

21

24

Literature data37,44

30000

810

wt%

5

5.85

K

343

This work

343

Analcite

15

15

25

40

Other

21

20

40

48

297

298

Clearly the high TOC shale was much more CO2-wet. Precisely, 𝜃𝑎 increased by

299

18° from 30° to 48°, and 𝜃𝑟 increased by 10° from 20° to 30°, when the pressure

300

increased from 0.1 MPa to 20 MPa for the low TOC shale; while 𝜃𝑎 increased by 110°

301

from 70° to 180°, and 𝜃𝑟 increased by 120° from 60° to 180°, when pressure increased 20

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 39

302

from 0.1 MPa to 30 MPa for the high TOC shale. It is thus clear that higher TOC and

303

thus a more hydrophobic surface leads to higher CO2-wettability, consistent with data

304

measured for CO2 on other hydrophobic surfaces, i.e. alkylated mica,74 quartz,94 silica,95

305

coal,69,96 and minerals aged in crude oil.97,98 Furthermore, one recent study reported the

306

effect of TOC on shale wettability (TOC ranged from 0.16 wt% to 23.4 wt%), and a

307

consistent conclusion with ours was reached.22 However, our result was inconsistent

308

with another recent report where CO2 wettability of organic-rich shale had no

309

relationship with organic matter content, which may be attributed to the remaining

310

mineral composition or the spatial distribution of the organics in the shale.43

311

312

4. CONCLUSION.

313 314

Organic-rich shale formations are an important fuel resource and they are also

315

potential CO2 sinks.34,35,46,99,100 In this context, CO2 and CH4 wettabilities are important

316

parameters as they determine the security and capacity of CO2 sequestration and shale

317

production respectively.12,13,30,37-39 However, their wettabilities in organic-rich shale

318

formation are poorly understood. We thus measured them as a function of pressure,

319

temperature, salinity, and ion type in order to aid CO2 geo-storage and CH4 production

320

schemes. Following conclusions were reached:

321

322

The brine contact angles of both, CO2 and CH4, increased with pressure, consistent

323

with literature data,4,44,55,56,61,65-71 while the CH4–shale system was far more water-wet

21

ACS Paragon Plus Environment

Page 23 of 39

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

324

than the analogue CO2 system, consistent with data measured on coal,70,71 However, a

325

similar CO2/CH4 density resulted in a similar wettability, consistent with literature

326

data.66 Thus, extrapolations from CH4 behavior to CO2 behavior need to be avoided

327

unless they have similar densities. The brine contact angles for CO2 and CH4 decreased

328

with temperature, consistent with literature data,61,62,66,67,74 while the temperature had a

329

stronger influence on CO2 wettability than on CH4 wettability. Higher salinity resulted

330

in higher brine contact angles, consistent with literature data,37,61,65,67,74,75,78,79 while the

331

brine contact angles for divalent cations were higher than those for monovalent cations,

332

consistent with literature data.37,61,78,80 However, the ion type only had an insignificant

333

influence (2°) on CH4 wettability. The effect of salinity and ion type on the brine contact

334

angles can be explained with the diffuse double layer theory with specific adsorption

335

of ions included.37,81,82 Higher TOC led to higher CO2 wettability, consistent with one

336

recent report,22 and other literature data.69,74,94-98

337 338

Especially noteworthy is that we reported CH4 wettability of shale here for the first time, which is of key importance to shale gas production.

339

340

 SUPPORTING INFORMATION

341

342

SI-Table 1. Advancing and receding brine contact angles (°) of CO2/CH4-brine -

343

shale for all cases measured in this work. T/k

P/MPa

Advancing for CO2

Receding for CO2 22

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

298

Page 24 of 39

0.1

76

74

77

73

72

74

5

90

90

89

88

87

86

10

180

180

180

180

180

180

20

180

180

180

180

180

180

30

180

180

180

180

180

180

0.1

73

72

73

71

70

68

5

88

87

88

85

84

84

10

130

132

133

110

111

108

20

180

180

180

180

180

180

30

180

180

180

180

180

180

0.1

70

71

70

60

60

58

5

80

82

79

70

71

71

10

125

125

125

107

108

109

20

130

130

130

120

118

121

30

180

180

180

180

180

180

T/k

P/MPa

Advancing for CH4

298

0.1

54

55

56

46

46

45

5

65

65

64

59

60

60

10

79

82

81

65

65

65

15

92

93

95

79

82

81

20

96

94

96

84

85

86

0.1

52

50

51

42

43

44

5

63

64

63

50

51

50

10

67

65

67

58

60

57

15

81

83

80

75

74

75

20

91

90

92

78

80

81

0.1

50

50

49

40

39

38

5

63

62

60

51

48

49

10

67

64

64

53

52

53

15

67

68

69

60

61

62

20

74

75

77

70

71

70

323

343

323

343

Receding for CH4

344 345 346



AUTHOR INFORMATION

347

348

Corresponding Author

349

*E-mail: [email protected]. 23

ACS Paragon Plus Environment

Page 25 of 39

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

350

ORCID

351

Bin Pan:0000-0003-3482-5660

352

Yajun Li: 0000-0001-6693-017X

353

Franca Jones: 0000-0002-8461-8291

354

Stefan Iglauer: 0000-0002-8080-1590

355

Notes

356

The authors declare no competing financial interest.

357

358



ACKNOWLEDGMENTS

359

360

This project was supported by the National Science and Technology Major Project

361

(2016ZX05023-001, 2017ZX05049-006), 973 Program (2014CB239103) and the

362

Chinese Scholarship Council.

363

364



ABBREVIATIONS

365

366

CO2

carbon dioxide

367

N2

nitrogen

24

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

368

Ar

argon

369

SF6

sulfur hexafluoride

370

He

helium

371

NaCl

sodium chloride

372

CaCl2

calcium chloride

373

MgCl2 magnesium chloride

Page 26 of 39

374

375



REFERENCES

376 377

(1) Intergovernmental Panel on Climate Change, 2005.

378

(2) Benson, S.M.; Cole, D.R. CO2 sequestration in deep sedimentary formations.

379 380 381 382 383 384 385

Elements. 2008, 4(5), 325-331. (3) Chadwick, R.A.; Eiken, O. Offshore CO2 storage: Sleipner natural gas field beneath the North Sea. Geological Storage of Carbon Dioxide. 2013, 227-253. (4) Kaveh, N.S.; Barnhoorn, A.; Wolf, K.K. Wettability evaluation of silty shale caprocks for CO2 storage. Int. J. Greenhouse Gas Control. 2016, 49, 425–435. (5) Iglauer, S. Dissolution trapping of carbon dioxide in reservoir formation brine— A carbon storage mechanism. 2011, Mass Transfer, InTech.

386

(6) Gao, R.S.; Sun, A.Y.; Nicot, J.P. Identification of a representative dataset for long-

387

term monitoring at the Weyburn CO2-injection enhanced oil recovery site,

388

Saskatchewan, Canada. Int. J. Greenhouse Gas Control. 2016, 54, 454-465. 25

ACS Paragon Plus Environment

Page 27 of 39

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

389

(7) Benson, S.; Cook, P. Underground Geological Storage. IPCC Special Report on

390

Carbon Dioxide Capture and Storage; Metz, B.,Davidson, O., de Coninck, H., Loos,

391

M., Meyer, L., Eds.; Cambridge University Press: New York, 2005; Chapter 5, pp

392

196−276.

393

(8) Li, Z.W.; Dong, M.Z.; Li, S.L.; Huang, S. CO2 sequestration in depleted oil and

394

gas reservoirs—caprock characterization and storage capacity. Energy Convers.

395

Manage. 2006, 47(11-12), 1372–1382.

396

(9) Zhang, Y.H.; Lebedev, M.; Sarmadivaleh, M.; Barifcani, A.; Iglauer, S. Swelling-

397

induced changes in coal microstructure due to supercritical CO2 injection. Geophys.

398

Res. Lett. 2016, 43(17), 9077-9083.

399

(10) Intergovernmental Panel on Climate Change (IPCC). IPCC Special Report on

400

Carbon Dioxide Capture and Storage; Metz, B., Davidson, O., de Coninck, H., Loos,

401

M., Meyer, L., Eds.; Cambridge University Press: New York, 2005; pp 442.

402 403

(11) Hesse, M. A.; Orr, F.M. Jr; Tchelepl, H.A. Gravity currents with residual trapping. J. Fluid Mech. 2008, 611, 35–60.

404

(12) Iglauer, S.; Pentlend, C.H.; Busch, A. CO2 wettability of seal and reservoir rocks

405

and the implications for carbon geo-sequestration. Water Resour. Res. 2014, 51(1), 729–

406

774.

407

(13) Pentland, C.H.; El-Maghraby, R.; Iglauer, S.; Blunt, M.J. Measurements of the

408

capillary trapping of supercritical carbon dioxide in Berea Sandstone. Geophys. Res.

409

Lett. 2011, 38(6).

410

(14) Juanes, R.; Spiteri, E.J.; Orr F.M. Jr.; Blunt, M.J. Impact of relative permeability 26

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

411

Page 28 of 39

hysteresis on geological CO2 storage. Water Resour. Res. 2006, 42(12), 395-397.

412

(15) Krevor, S.; Blunt, M.J.; Benson, S.M.; Pentland, C.H.; Reynolds, C.; Al-

413

Menhali, A.; Niu, B. Capillary trapping for geologic carbon dioxide storage – From

414

pore scale physics to field scale implications. Int. J. Greenhouse Gas Control. 2015, 52,

415

221-237.

416

(16) Al-Khdheeawi, E.A.; Vialle, S.; Barifcani, A.; Sarmadivaleh, M.; Iglauer, S.

417

Influence of CO2-wettability on CO2 migration and trapping capacity in deep saline

418

aquifers. Greenhouse Gases: Science and Technology. 2016, DOI: 10.1002/ghg.1648.

419

(17) Al-Khdheeawi, E.A.; Vialle, S.; Barifcani, A.; Sarmadivaleh, M; Iglauer, S.

420

Impact of reservoir wettability and heterogeneity on CO2-plume migration and trapping

421

capacity. Int. J. Greenhouse Gas Control. 2017, 58, 142-158.

422 423 424 425

(18) Oelkers, E.H.; Gislason S.R.; Matter, J. Mineral Carbonation of CO2. Elements. 2008, 4(5), 333-337. (19) Gaus, I. Role and impact of CO2/rock interactions during CO2 storage in sedimentary rocks. Int. J. Greenhouse Gas Control. 2010, 4, 73–89.

426

(20) Busch, A.; Alles, S.; Gensterblum, Y.; Prinz, D.; Dewhurst, D.N.; Raven, M.D.;

427

Stanjek, H.; Krooss, B.M. Carbon dioxide storage potential of shales. Int. J. Greenhouse

428

Gas Control. 2008, 2, 297–308.

429

(21) Giesting, P.; Guggenheim, S; Busch, A. Interaction of carbon dioxide with Na-

430

exchanged montmorillonite at pressures to 640 bars: Implications for CO2 sequestration.

431

Int. J. Greenhouse Gas Control. 2012, 8, 73-81.

432

(22) Arif, M.; Lebedev, M.; Barifcani, A.; Iglauer, S. Influence of shale-total organic 27

ACS Paragon Plus Environment

Page 29 of 39

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

433 434 435

Energy & Fuels

content on CO2 geo-storage potential. Geophys. Res. Lett. 2017, 44, 8769–8775. (23) Iglauer, S. CO2-water-rock wettability: variability, influencing factors and implications for CO2 geo-storage. Acc. Chem. Res. 2017, 50 (5), 1134–1142.

436

(24) Kelvin B. Gregory, K.B.; Vidic, R.D.; Dzombak, D.A.; Water management

437

challenges associated with the production of shale gas by hydraulic fracturing. Elements.

438

2011, 7(3), 181-186.

439 440

(25) EIA (2011) Annual Energy Outlook 2011. United States Department of Energy, Energy Information Administration DOE/EIA-0383: www.eia.doe.gov/forecasts/aeo/

441

(26) Roshan, H.; Ehsani, S.; Marjo, C.E.; Andersen, M.S.; Acworth, R.I. Mechanisms

442

of water adsorption into partially saturated fractured shales: An experimental study.

443

Fuel. 2015, 159, 628-637.

444

(27) Curtis, J.B. Fractured shale-gas systems. AAPG Bull. 2002, 86(11), 1921-1938.

445

(28) Louk, K.; Ripepi, N.; Luxbacher, K.; Gilliland, E.; Tang, X.; Keles, C.; Schlosser,

446

C.; Diminick, E.; Keim, S.; Amante, J.; Michael, K. Monitoring CO2 storage and

447

enhanced gas recovery in unconventional shale reservoirs: Results from the Morgan

448

County, Tennessee injection test. J NAT GAS SCI ENG. 2017, 45, 11-25.

449

(29) Eshkalak, M.O.; Al-Shalabi, E.W.; Sanaei, A.; Aybar, U.; Sepehrnoori, K.

450

Enhanced gas recovery by CO2 sequestration versus re-fracturing treatment in

451

unconventional shale gas reservoirs. 2014, DOI: 10.2118/172083-MS.

452

(30) Xu, R.N.; Zeng, K.C.; Zhang, C.W.; Jiang, P.X. Assessing the feasibility and CO2

453

storage capacity of CO2 enhanced shale gas recovery using triple-porosity reservoir

454

model, Appl. Therm. Eng. 2017, 115 1306–1314. 28

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 30 of 39

455

(31) Al-Yaseri, A.Z.; Roshan, H.; Zhang, Y.; Rahman, T.; Lebedev, M.; Barifcani, A.;

456

Iglauer, S. Effect of the temperature on CO2/Brine/Dolomite Wettability: Hydrophilic

457

versus Hydrophobic Surfaces. Energy Fuels. 2017, 31, 6329−6333.

458

(32) Cui, X.J.; Bustin, R.M.; Chikatamaria, L. Adsorption-induced coal swelling and

459

stress: Implications for methane production and acid gas sequestration into coal seams

460

J. Geophys. Res., 2007, 112, DOI: 10.1029/2004JB003482.

461

(33) Perera, M.S.A.; Ranjith, P.G.; Viete, D.R.; Choi, S.K. The effects of injection

462

and production well arrangement on carbon dioxide sequestration in deep, unmineable

463

coal seams: a numerical study. International Journal of Coal Preparation and

464

Utilization. 2012, 32(5), 211-224.

465

(34) Schepers, K.C.; Nuttall, B.C.; Oudinot, A.Y.; Gonzalez, R.J. Reservoir modeling

466

and simulation of the Devonian gas shale of eastern Kentucky for Enhanced gas

467

recovery and CO2 storage. 2009, DOI: org/10.2118/126620-MS.

468

(35) Fathi, E.; Akkutlu, I.Y.; Multi-component gas transport and adsorption effects

469

during CO2 injection and enhanced shale gas recovery. Int. J. Coal Geol. 2014, 123, 52-

470

61.

471

(36) Oldenburg, C.M.; Pruess, K.; Benson, S.M. Process modeling of CO2 injection

472

into natural gas reservoirs for carbon sequestration and enhanced gas recovery. Energy

473

Fuels. 2001, 15(2), 293-298.

474 475 476

(37) Roshan, H.; Al-Yaseri, A.Z.; Sarmadivaleh, M.; Iglauer, S. On wettability of shale rocks. J. Colloid Interface Sci. 2016, 475, 104-111. (38) Liang, L.X.; Luo, D.X.; Liu, X.J.; Xiong, J. Experimental study on the wettability 29

ACS Paragon Plus Environment

Page 31 of 39

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

477

and adsorption characteristics of Longmaxi Formation shale in the Sichuan Basin,

478

China. J NAT GAS SCI ENG. 2016, 33, 1107–111.

479

(39) Durucan, S.; Ahsan, M.; Syed, A.; Shi, J.Q.; Korre, A. Two phase relative

480

permeability of gas and water in coal for enhanced coalbed methane recovery and CO2

481

storage. Energy Procedia. 2013, 37(37), 6730-6737.

482

(40) Meyers, P.A.; Bernasconi, S.M.; Forster, A. Origins and accumulation of organic

483

matter in expanded Albian to Santonian black shale sequences on the Demerara Rise,

484

South American margin. Org. Geochem. 2006, 37(12), 1816–1830.

485

(41) Heller, R.; Zoback, M. Adsorption of methane and carbon dioxide on gas shale

486

and pure mineral samples. Journal of Unconventional Oil and Gas Resources. 2013, 8,

487

14–24.

488

(42) Li, A.; Ding, W.L.; He, J.H.; Dai, P.; Yin, S.; Xie, F. Investigation of pore

489

structure and fractal characteristics of organic-rich shale reservoirs: A case study of

490

Lower Cambrian Qiongzhusi formation in Malong block of eastern Yunnan Province,

491

South China. Mar. Pet. Geol. 2016, 70, 46–57.

492

(43) Guiltinan E.J.; Cardenas, M.B.; Bennett, P.C.; Zhang, T.W.; Espinoza, D.N. The

493

effect of organic matter and thermal maturity on the wettability of supercritical CO2 on

494

organic shales. Int. J. Greenhouse Gas Control. 2017, 65, 15-22.

495

(44) Iglauer, S.; Al-Yaseri, A.Z.; Rezaee, R.; Lebedev, M. CO2 wettability of

496

caprocks: Implications for structural storage capacity and containment security.

497

Geophys. Res. Lett. 2015, 42, 9279–9284.

498

(45) Chaudhary, K.; Guiltinan, E.J.; Cardenas, M.B.; Maisano, J. A.; Ketcham, R.A.; 30

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 32 of 39

499

Bennett, P.C. Wettability measurement under high P-T conditions using X-ray imaging

500

with application to the brine-supercritical CO2 system. Geochem. Geophys.Geosyst.

501

2015, 16, 2858–2864.

502

(46) Zhiyuan, T.; Clarens, A. Estimating the carbon sequestration capacity of shale

503

formations using methane production rates. Environ. Sci. Technol. 2013, 47,

504

11318−11325.

505

(47) White, C.M.; Smith, D.H.; Jones, K.L.; Goodman, A.L.; Jikich, S.A.; Lacount,

506

R.B.; Dubose, S.B.; Ozdemir, E.; Morsi, B.I.; Schroeder, K.T. Sequestration of carbon

507

dioxide in coal with enhanced coalbed methane recovery—a review. Energy Fuels.

508

2005, 19(3), 659–724.

509

(48) Bachu, S.; Bonijoly, D.; Bradshaw, J.; Burruss, R.; Holloway, S.; Christensen,

510

N.P.; Mathiassen, O.M. CO2 storage capacity estimation: Methodology and gaps. Int. J.

511

Greenhouse Gas Control. 2007, 1, 430–443.

512

(49) Gong, H.J.; Zhang, H.; Xu, L.; Li, K.N.; Yu, L.; San, Q.; Li, Y.J.; Dong, M.Z. The

513

synergistic effect of branched-preformed particle gel and hydrolyzed polyacrylamide

514

on further-enhanced oil recovery after polymer flooding. Energy Fuels. 2017, 31,

515

7904−7910.

516

(50) Gong, H.J.; Zhang, H.; Xu, L.; Li, K.N.; Yu, L.; Li, Y.J.; Dong, M.Z. Further

517

enhanced oil recovery by branched-preformed particle gel/HPAM/surfactant mixed

518

solutions after polymer flooding in parallel-sandpack models. RSC Adv. 2017, 7,

519

39564–39575.

520

(51) Schmidt, G.W. Interstitial water composition and geochemistry of deep gulf coast 31

ACS Paragon Plus Environment

Page 33 of 39

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

521

Energy & Fuels

shales and sandstones. AAPG Bulletin. 1973, 57(2), 321-337.

522

(52) Lucas, A.J.; Peyron, M.; Pierens, G.K.; Carpenter, A.T.; Hall, L.D.; Stewart, R.C.;

523

Potter, G.F.; Phelps, D.W. Rigorous evaluation of the experimental requirements for

524

performing quantitative porosity measurements by NMR. 1993, SPE-27420-MS.

525

(53) Clarkson, CR.; Solano, N.; Bustin, R.M.; Bustin, A.M.M.; Chalmers, G.R.L.; He,

526

L.; Melnichenko, Y.B.; Radlinski, A.P.; Blach, T.P. Pore structure characterization of

527

North American shale gas reservoirs using USANS/SANS, gas adsorption, and mercury

528

intrusion. Fuel. 2013, 103(1), 606-616.

529 530

(54) John, S.C. A technique for faster pulse-decay permeability measurements in tight rocks. SPE formation evaluation. 1997, 12, 19-26.

531

(55) Sarmadivaleh, M.; Al-Yaseri, A.; Iglauer, S. Influence of temperature and

532

pressure on quartz-water-CO2 contact angle and CO2-water interfacial tension. J.

533

Colloid Interface Sci. 2015,441, 59-64.

534

(56) Iglauer, S.; Salamah, A.; Sarmadivaleh, M.; Liu, K.Y.; Phan C. Contamination of

535

silica surfaces: impact on water–CO2–quartz and glass contact angle measurements. Int.

536

J. Greenhouse Gas Control. 2014, 22, 325-328.

537 538

(57) Lander, L.M.; Siewierski, L.M.; Brittan, W.J.; Vogler, E.A. A systematic comparison of contact angle methods Langmuir. 1993, 9(8), 2237-2239.

539

(58) Al-Yaseri, A.Z.; Sarmadivaleh, M.; Saeedi, A.; Lebedev, M.; Barifcanic, A;

540

Iglauer, S. N2+CO2+NaCl brine interfacial tensions and contact angles on quartz at CO2

541

storage site conditions in the Gippsland basin, V ictoria/ Australia. J. Pet. Sci. Technol.

542

2015, 129, 58–62. 32

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 34 of 39

543

(59) Dake; L. P. Fundamentals of Reservoir Engineering. 1978, Elsevier, Amsterdam.

544

(60) Rahman, T.; Lebedev, M.; Barifcani, A.; Iglauer, S. Residual trapping of

545

supercritical CO2 in oil-wet sandstone. J. Colloid Interface Sci. 2016, 469, 63−68.

546

(61) Al-Yaseri, A. Z.; Lebedev, M.; Barifcani, A.; Iglauer, S. Receding and advancing

547

(CO2+ brine+ quartz) contact angles as a function of pressure, temperature, surface

548

roughness, salt type and salinity. J.Chem. Thermodyn. 2016, 93, 416−423.

549

(62) Iglauer, S.; Mathew, M.S.; Bresme, F. Molecular dynamics computations of

550

brine–CO2 interfacial tensions and brine–CO2–quartz contact angles and their effects

551

on structural and residual trapping mechanisms in carbon geo-sequestration. J. Colloid

552

Interface Sci. 2012, 386(1), 405-14.

553

(63) Javanbakht, G.; Sedghi, M.; Welch, W.; Goual, L. Molecular dynamics

554

simulations of CO2/water/quartz interfacial properties: Impact of CO2 dissolution in

555

water. Langmuir. 2015, 31, 5812–5819.

556

(64) Arif, M.; Barifcani, A.; Iglauer, S. Solid/CO2 and solid/water interfacial tensions

557

as a function of pressure, temperature, salinity and mineral type: Implications for CO2-

558

wettability and CO2 geo-storage. Int. J. Greenhouse Gas Control. 2016, 53, 263-273.

559 560

(65) Chiquet, P.; Broseta, D.; Thibeau, S. Wettability alteration of caprock minerals by carbon dioxide. Geofluids. 2007, 7, 112–122.

561

(66) Al-Yaseri, A.Z.; Roshan, H.; Lebedev, M.; Barifcani, A.; Iglauer, S.

562

Dependence of quartz wettability on fluid density. Geophys. Res. Lett. 2016, 43, 3771-

563

3776.

564

(67) Arif, M.; Al-Yaseri, A.Z.; Barifcani, A.; Lebedev, M.; Iglauer, S. Impact of 33

ACS Paragon Plus Environment

Page 35 of 39

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

565

pressure and temperature on CO2–brine–mica contact angles and CO2–brine interfacial

566

tension: Implications for carbon geo-sequestration. J. Colloid Interface Sci. 2016, 462,

567

208–215.

568

(68) Kaveh, N.S.; Barnhoon, A; Schoemaker, FC; Wolf, K.H. Study of shale

569

wettability for CO2 storage. 77th EAGE Conference& Exhibition. 2015, DOI:

570

10.3997/2214-4609.201413370.

571

(69) Arif, M.; Barifcani, A.; Lebedev, M.; Iglauer, S. CO2-wettability of low to high

572

rank coal seams: Implications for carbon sequestration and enhanced methane recovery.

573

Fuel. 2016, 181, 680–689.

574 575 576 577

(70) Saghafi, A.; Javanmard, H.; Pinetown, K. Study of coal gas wettability for CO 2 storage and CH4 recovery. Geofluids. 2014, 14, 310–325. (71) Saghafi, A.; Pinetown, K.; Javanmard, H. Gas wettability of coal and implications for gas desorption and drainage. Coal operators’ conference. 2014.

578

(72) Georgiadis, A.; Maitland, G.; Trusler, J.P.M.; Bismarck, A.; Interfacial tension

579

measurements of the (H2O+ CO2) system at elevated pressures and temperatures.

580

Euroeconomica. 2010, 94(6), 7-14.

581

(73) Setzmann, U; Wagner, W. A new equation of state and tables of thermodynamic

582

properties for methane covering the range from the melting line to 625 K at pressures

583

up to 1000 MPa. J. Phys. Chem. Ref. Data. 1991, 20(20), 1061-1151.

584

(74) Arif, M.; Barifcani, A.; Lebedev, M.; Iglauer, S. Structural trapping capacity of

585

oil-wet caprock as a function of pressure, temperature and salinity. Int. J. Greenhouse

586

Gas Control. 2016, 50, 112-120. 34

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

587 588

Page 36 of 39

(75) Farokhpoor, R.; Lindeberg, E. Wettability behaviour of CO2 at storage conditions. Int. J. Greenhouse Gas Control. 2013, 12, 18-25.

589

(76) Saraji, S.; Goual, L.; Piri, M.; Plancher, H. Wettability of supercritical carbon

590

dioxide/water/quartz systems: Simultaneous measurement of contact angle and

591

interfacial tension at reservoir conditions. Langmuir. 2013, 29(23), 6856–6866.

592

(77) Saraji, S.; Piri, M.; Goual, L. The effects of SO2 contamination, brine salinity,

593

pressure, and temperature on dynamic contact angles and interfacial tension of

594

supercritical CO2/brine/quartz systems. Int. J. Greenhouse Gas Control. 2014, 28, 147–

595

155.

596

(78) Arif, M.; Jones, F.; Barificani, A.; Iglauer, S. Electrochemical investigation of the

597

effect of temperature, salinity and salt type on brine/mineral interfacial properties. Int.

598

J. Greenhouse Gas Control. 2017, 59, 136-147.

599 600 601 602

(79) Broseta, D.; Tonnet, N.; Shah, V. Are rocks still water-wet in the presence of dense CO2 or H2S? Geofluids. 2012, 12(4), 280–94. (80) Arif, M.; Jones, F.; Barifcani, A.; Iglauer, S. Influence of surface chemistry on interfacial properties of low to high rank coal seams. Fuel. 2017, 194, 211–221.

603

(81) Hartley, P.G.; Grieser F.; Mulvaney P.; Stevens G.W. Surface forces and

604

deformation at the oil-water interface probed using AFM force measurement. Langmuir.

605

1999, 15(21), 7282-7289.

606 607 608

(82) Kang, K.H. How electrostatic fields change contact angle in electrowetting. Langmuir. 2002, 18(26), 10318-10322. (83) Vane, L.M.; Zang, G.M. Effect of aqueous phase properties on clay partical zeta 35

ACS Paragon Plus Environment

Page 37 of 39

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

609

potential and electro-osmotic permeability: Implications for electro-kinetic soil

610

remediation Processes. J. Hazard. Mater. 1997, 55, l-22.

611

(84) Alroudhan, A.; Vinogradov, J.; Jackson M.D. Zeta potential of intact natural

612

limestone: Impact of potential-determining ions Ca, Mg and SO4. Colloids Surf., A.

613

2016, 493, 83–98.

614

(85) Roshan, H.; Andersen, M.S.; Rutlidge, H., Marjo, C.E.; Acworth, R.I.

615

Investigation of the kinetics of water uptake into partially saturated shales. Water

616

Resour. Res.2016, 52(4), 2420-2438.

617

(86) Mugele, F.; Siretanu, I.; Kumar, N.; Wang, L.; Maestro, A.; Duits, M.; Van den

618

Ende, D.; Collins, I. Insights from ion adsorption and contact angle alteration at mineral

619

surfaces for low-salinity water flooding. 2016 SPE-169143-PA. DOI: 10.2118/169143-

620

PA.

621 622

(87) Saraji, S.; Goual, L.; Piri, M. Dynamic adsorption of asphaltenes on quartz and calcite packs in the presence of brine films, Colloids Surf., A. 2013, 434, 260–267.

623

(88) Kasha, A.; Al-Hashim, A.; Abdallah, W; Taherian, R.; Sauerer, B. Effect of Ca2+,

624

Mg2+ and SO42− ions on the zeta potential of calcite and dolomite particles aged with

625

stearic acid, Colloids Surf., A, 2015, 482, 290-299.

626

(89) Ebeling, D.; Van den Ende, D.; Mugele F. Electrostatic interaction forces in

627

aqueous salt solutions of variable concentration and valency. Nanotechnology, 2011,

628

22(30), 305706.

629 630

(90) Cavale, A.; Bera, B.; Van den Ende, D.; Mugele F. Analytical model for the electrowetting properties of oil-water-solid systems. Phys. Rev. E. 2016, 93, 042606. 36

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 38 of 39

631

(91) Kheirandish, S.; Jabbari, E. Effect of surface polarity on wettability and friction

632

coefficient of silicone rubber/poly (acrylic acid) hydrogel composite. Colloid. Polym.

633

Sci. 2006. 284 (12), 1411-1417.

634 635 636 637

(92) Binks, B.P.; Clint, J.H. Solid wettability from surface energy components: relevance to pickering emulsions. Langmuir. 2002, 18(4), 1270-1273. (93) Vroman, L. Effect of adsorbed proteins on the wettability of hydrophilic and hydrophobic solids. Nature. 1962,196, 476-477.

638

(94) Espinoza, D.N.; Santamarina, J.C. Water-CO2-mineral systems: Interfacial

639

tension, contact angle, and diffusion: Implications to CO2 geological storage. Water

640

Resour. Res. 2010, 46 (46), 58-72.

641 642

(95) Dickson, J.L.; Gupta, G.; Horozov, T.S.; Binks, B.P.; Johnston, K.P. Wetting phenomena at the CO2/water/glass interface. Langmuir. 2006, 22(5), 2161-2170.

643

(96) Arif, M.; Barifcani, A.; Zubair, T.; Lebedev, M; Iglauer, S. CO2 wettability of

644

shales and coals as a function of pressure, temperature and rank: Implications for CO2

645

sequestration and enhanced methane recovery. 2016 PAPG/SPE Pakistan Section ATCE.

646

DOI: 10.2118/185308-MS.

647

(97) Al-Menhali, A.S.; Menke, H.P.; Blunt, M.J.; Krevor, S.C. Pore scale observations

648

of trapped CO2 in mixed-wet carbonate rock: Applications to storage in oil fields.

649

Environ. Sci. Technol. 2016, 50, 10282–10290.

650

(98) Al-Menhali, A.S. and Krevor, S.C. Capillary trapping of CO2 in oil Reservoirs:

651

Observations in a mixed wet carbonate rock. Environ. Sci. Technol. 2016, 50,

652

2727−2734. 37

ACS Paragon Plus Environment

Page 39 of 39

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

653

(99) Yu, W.; Al-Shalabi, E.W.; Sepehrnoori, K. A sensitivity study of potential CO2

654

injection for enhanced gas recovery in Barnett shale reservoirs. 2014, DOI:

655

org/10.2118/169012-MS.

656

(100) Liu, F.Y.; Ellett, K.; Xiao, Y.T.; Rupp, J.A. Assessing the feasibility of CO2

657

storage in the New Albany Shale (Devonian–Mississippian) with potential enhanced

658

gas recovery using reservoir simulation. Int. J. Greenhouse Gas Control. 2013, 17, 111-

659

126.

38

ACS Paragon Plus Environment