Effects of Dissolved Oxygen on Water Imbibition in Gas Shales

Mar 12, 2018 - According to Figure 5, potassium and sodium ion concentrations increase during both oxic and degassed water imbibition experiments and ...
2 downloads 8 Views 2MB Size
Subscriber access provided by UNIVERSITY OF TOLEDO LIBRARIES

Fossil Fuels

Effects of Dissolved Oxygen on Water Imbibition in Gas Shales Mingxiang Xu, Mojtaba Binazadeh, Ashkan Zolfaghari, and Hassan Dehghanpour Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b03955 • Publication Date (Web): 12 Mar 2018 Downloaded from http://pubs.acs.org on March 26, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

1

Effects of Dissolved Oxygen on Water Imbibition in Gas Shales

2

Mingxiang Xu1, Mojtaba Binazadeh1,2, Ashkan Zolfaghari1, and Hassan Dehghanpour1

3

1

4

University of Alberta, Edmonton, Alberta, T6G 2W2, Canada.

5

2

6

Shiraz, Fars, Iran

7

Abstract

8

Understanding the water uptake of gas shales is critical for designing and optimizing hydraulic

9

fracturing operations during which a large volume of fracturing water containing dissolved oxygen is

10

injected into tight reservoirs. Recent studies show that the dissolved oxygen may promote oxidation

11

reactions which can affect salinity and pH value of flowback water; however, the effects of dissolved

12

oxygen and oxidation reactions on water imbibition into the shale matrix and on the concentration of

13

individual ions in flowback water are still poorly understood. In this study, we conduct water

14

imbibition experiments under degassed and oxic conditions, and measure the imbibed water volume

15

and concentrations of different ions in water. The results show that the initial rate and final amount of

16

water imbibition are higher under degassed conditions compared with that under oxic conditions.

17

These differences are mainly due to the enhanced dissolution of air in the shale pore network into the

18

imbibing water under degassed conditions and the consequent increase in relative permeability of

19

water. The results also suggest that oxidation of pyrite by dissolved oxygen produces sulfate and iron

20

ions, as well as iron-compound precipitations. Pyrite oxidation is also supported by the abundance of

21

pores in the vicinity of pyrite minerals observed in the SEM/EDS images.

22

1. Introduction

23

Recent advances in horizontal drilling and hydraulic fracturing have unlocked the vast unconventional

Department of Civil and Environmental Engineering, School of Mining and Petroleum Engineering,

Department of Chemical Engineering, School of Chemical and Petroleum Engineering, Shiraz University,

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

24

hydrocarbon resources such as shale and tight formations. Hydraulic fracturing operations include

25

injection of approximately 10-20 million liters of water (cross-linked gel, slickwater, energized and

26

hybrids) into the shale/tight formations per well to create fractures.1-3 Some wells in some thick shale

27

formations (e.g., Horn River Basin) consume up to 75 million liters of slick water for hydraulic

28

fracturing operations.4 Then, the wells are sometimes shut-in for a period of time, (i.e., soaking period)

29

to increase the early-time oil and gas production rate.5,6 During the hydrocarbon production stage,

30

only a fraction (usually less than 30% within 100 days of flowback) of the injected fracturing water

31

(i.e., flowback water) can be recovered.7-9

32

Laboratory imbibition experiments have been widely used to investigate the effects of shale-water

33

interactions. For example, water imbibition in shales may cause dissolution of precipitated salts in

34

shale and raise osmosis effect10,11, increase water saturation in the vicinity of water-rock interface

35

(water blockage) 12-14 , and reduce hydrocarbon relative permeability. 15-18 It is also reported that

36

shale-water interactions may reduce fracture conductivity and effective length 19 , 20 , induce

37

micro-fractures in rock matrix due to clay swelling21,22, and finally, alter shale strength and cause

38

failures around the wellbore.23-25

39

A common challenge for applying laboratory imbibition data to evaluate fluid-rock interactions during

40

fracturing and soaking periods is the differences between laboratory and field conditions:

41

Temperature. Previous studies indicate that increasing the temperature can accelerate the water

42

imbibition by increasing the wetting affinity of rock to water.26,27 Increasing the temperature increases

43

the activity of divalent ions such as Mg2+ and SO  to pair with the surface anions of carbonate

44

minerals, and enhances the affinity of rock toward water.28

45

Pressure. Usually rock samples are free to expand during laboratory water imbibition

46

experiments.29-31 However, under the field conditions, overburden pressure limits rock expansion.32

47

Ghanbari and Dehghanpour22 conducted water imbibition experiments on confined/unconfined shale

48

samples and found that confined samples have less induced microfractures and less water imbibition. 2

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

49

Water-Rock Interface Area. The water-rock interface area differs dramatically at laboratory and field

50

conditions. Previous studies show that increasing shale-water interface area per unit volume of shale

51

significantly increases the rate and final volume of water imbibition.33 Roshan et al.34 showed that

52

the water-rock interface area affects the exchangeable cations in flowback water. Bearinger35 showed

53

that the salt concentration in flowback water is related to the complexity of fracture network.

54

Low oxygen-content under reservoir conditions should also be considered in imbibition experiments.

55

Shale reservoirs usually have no (or limited) oxygen content (anoxic conditions) before fracturing

56

process. Anoxic reservoir conditions can be due to oxygen consumption by organic matter and lack of

57

oxygen supply in the reservoir.36 Oxygen is introduced by fracturing water into shale reservoirs,

58

creating oxic conditions. Shales usually contain chemically-reactive compounds such as organic

59

matter and pyrite, which may participate in oxidation reactions. 37 - 41 Recent studies show that

60

oxidation of pyrite may change the pore structure and increase the pore connectivity.42,43 Rowan et

61

al.44 observed a positive correlation between oxygen and chloride content in flowback water, which

62

suggests that the oxygen dissolved in injected water may affect the geochemical reactions. Zolfaghari

63

et al.33 found that the concentration of dissolved oxygen in flowback water may be lower than that in

64

deionized water. They concluded that the dissolved oxygen may cause oxidation reactions, which

65

could affect salt dissolution, salt precipitation and pH change in flowback water. Thus, understanding

66

the effects of dissolved oxygen on water imbibition and ion dissolution is important for more accurate

67

interpretation of flowback water chemistry.

68

This article extends the previous studies33 to 1) investigate the role of dissolved oxygen on water

69

imbibition rate/volume and 2) evaluate the effects of dissolved oxygen on concentration of major ions

70

(i.e. potassium and sodium ions) and redox-sensitive ions (i.e. iron and sulfate ions) in flowback water.

71

We perform comparative imbibition experiments using air-saturated de-ionized water (oxic water) and

72

degassed de-ionized water (degassed water). The rest of the paper is organized as follows: Section 2

73

describes the materials and methods. Section 3 shows the imbibition profiles and concentration

74

profiles of different ions for the oxic and degassed water. Section 4 describes the conclusions of this 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 28

75

study.

76

2. Materials and Methods

77

2.1 Shale Samples

78

Three pairs of shale samples from three different depths of the Evie Formation of the Horn River

79

Basin are labeled as EV-1, EV-2, and EV-3 and used for imbibition experiments. To eliminate the

80

effect of sample size, all samples are cut into 35×35×25 mm cubes with approximate mass of 100 g.

81

Table 1 provides the average depth, porosity, permeability, and mineralogy of the samples. The

82

porosity, permeability, and mineralogy are measured by helium porosimeter, pulse decay permeameter,

83

and x-ray diffraction (XRD) respectively using different samples from similar depth.

84

Table 1. Approximate depth, porosity, permeability and average rock composition for EV-1, EV-2,

85

and EV-3 samples. Sample

EV-1

EV-2

EV-3

Depth, m

2672.5

2681.2

2688.2

Porosity, %

6.4

5.0

4.6

Permeability, nD

575

384

372

Non-Clay Content (wt. %) Quartz

78.1

51.9

65.4

K-feldspar

5.2

4.9

3.4

Plagioclase

1.3

5.7

5.0

Calcite

3.1

13.3

8.4

Dolomite

0.7

2.6

1.1

Pyrite

1.8

3.0

3.0

Total Non-Clay

90.5

81.4

86.5

Clay Content (wt. %)

4

ACS Paragon Plus Environment

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Illite/Smectite

2.9

6.6

5.2

Illite+Mica

6.7

12.0

8.3

Total Clay

9.6

18.6

13.5

86

2.2 Imbibition Experiments

87

In order to investigate the effect of dissolved oxygen on water imbibition in the shale samples, for

88

each pair of samples, we immerse one sample in the air-saturated de-ionized (DI) water (oxic water)

89

and the other sample in the degassed DI water (degassed water). We consider the DI water which was

90

in contact with air (air-saturated) for 1 week as oxic water. The dissolved oxygen content of oxic water

91

is about 8.3 mg/L, which is measured by rugged dissolved oxygen (RDO) optical electrode (Thermo

92

ScientificTM OrionTM). To prepare the degassed water, we place 500 mL of DI water in a desiccator

93

(vaccum chamber) as shown in Figure 1. The desiccator is then sealed and connected to a vacuum of 1

94

kPa for 1 week. The dissolved oxygen content of degassed water is about 0.7 mg/L. Generally, the

95

pressure difference can potentially impact the imbibition results. However, the pressure difference

96

between oxic and degassed experiments is relatively low (100 KPa) compared with that between

97

laboratory and field conditions (usually more than 10,000 KPa). For example, the reservoir pressure in

98

the Horn River Basin can be as high as 53,000 KPa.45 In this study, compared with the significant

99

pressure difference in laboratory and subsurface, we assume that the relatively-small pressure

100

difference between the oxic and degassed experiments has a negligible impact on the rock pore

101

volume and imbibed mass. The imbibition experiments are conducted according to the following

102

steps:

103

1.

104 105

samples. 2.

106 107

All the samples are dried in the oven at 100 ℃ for 24 hours to remove the moisture in the

For each pair of samples, one sample is placed in an imbibition cell filled with 500 mL oxic water. The other sample is placed in an imbibition cell filled with 500 mL of degassed water.

3.

All the imbibition cells filled with degassed water are placed in a vacuum chamber connected to 5

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

108

a vacuum pump at 1 kPa to ensure degassed conditions. It is assumed that vacuum pressure

109

during degassed conditions has a negligible effect on the pore volume and imbibed mass. The

110

schematic of the vacuum set-up is shown in Figure 1.

111

4.

At each time step, the shale samples are taken out of the imbibition cell to measure the samples’

112

mass. We reported the normalized imbibed mass by dividing the mass gain by the initial mass of

113

the dried sample. Sartorius™ Entris™ Toploading Balance with the accuracy of ±0.01 g is used

114

to measure the samples’ mass. The dissolved oxygen concentrations of the oxic and degassed

115

water are measured using RDO optical electrode.

116

5.

2 mL of the oxic and degassed water are periodically collected for measuring iron concentrations

117

using inductively coupled plasma mass spectrometry (ICP-MS) and ion chromatography (IC).

118

During the measurements in steps 4 and 5, the shale and water samples are exposed to air. This may

119

cause oxygen adsorption on the shale and water samples. We maintained vacuum pressure at 1 kPa

120

during the experiments to reduce potential oxygen adsorption on the shale samples during the

121

measurements. The 2 mL collected degassed water is sealed in 2 mL air-tight vials to prevent oxygen

122

contamination.

123 124

Figure 1. The schematic of the vacuum set-up used for the imbibition tests under degassed conditions.

125

2.3 SEM-EDS Imaging

126

Scanning

electron

microscope

(SEM),

scanning

helium

6

ACS Paragon Plus Environment

ion

microscope

(SHIM),

and

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

127

energy-dispersive x-ray spectroscopy (EDS) are used to visualize the pores and obtain the elemental

128

maps in EV-2 sample. One end piece from EV-2 sample is prepared for SEM-EDS imaging. We first

129

cut the end piece into 1×1×0.5 cm cubes. The top surface of the sample is then mechanically polished

130

using 600-, 1000-, and 2000-grit polishing pads. The sample surface is further polished using

131

argon-ion milling to minimize the influence of roughness and artifacts of the sample surface on the

132

SEM images. All SEM/SHIM images are obtained using ~15 Kv beam energy.

133

3. Results and Discussion

134

In this section, we first compare the imbibition profiles of the oxic and degassed imbibition

135

experiments. Then, we discuss the observed difference in ion concentration profiles of the oxic and

136

degassed water samples.

137

3.1 Imbibition Profiles

138

Figure 2 shows the water imbibition profiles versus time. The imbibed water mass is normalized

139

through dividing it by the mass of the dry shale sample. The imbibition profiles increase over time and

140

reach to a plateau for both oxic and degassed experiments. Figure 2 shows that 1) the imbibition rate

141

and the final normalized mas of imbibed water are higher for degassed experiments compared with

142

those for oxic experiments, and 2) EV-1 samples show the lowest gap between oxic and degassed

143

water imbibition profiles while EV-3 samples show the largest gap. We explain these two observations

144

by investigating the dissolved oxygen in Section 3.2 and ion concentration profiles for both oxic and

145

degassed water samples in Section 3.3. The measured shale sample mass is also added in Table A in

146

the Appendix.

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

147

Figure 2. Mass of imbibed oxic and degassed water for (a) EV-1, (b) EV-2, and (c) EV-3 samples.

148

3.2 Dissolved Oxygen Profiles

149 150

Figure 3. The surface of EV-1 shale samples in a. oxic water and b. degassed water 2 hours after

151

starting the imbibition experiments. The blue areas in the figure are reference areas used for

152

determining the sample depth.

153

As can be seen in Figure 3, gas bubbles are observed on the shale samples’ surface when oxic water

154

imbibes into the samples 2 hours after starting the imbibition experiments. We did not observe gas

155

bubbles for degassed experiments which could be due to the rapid dissolution of gas existing in the

156

shale sample pore space into degassed water. The low pressure in vacuum chamber not only reduces

157

the concentration of the dissolved oxygen in degassed water, but also reduces the concentration of

158

other dissolved gases such as carbon dioxide and nitrogen.46

159

The oxic water imbibing into the shale sample displaces the gas out of the pore system. The gas 8

ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

160

bubbles expelled out by the imbibed oxic water were also observed in previous studies.34,47 The gas

161

bubbles in Figure 3(a) are probably the air initially present in the pores of the shale sample. In addition,

162

the gas bubbles observed in Figure 3(a) may also contain the gas produced by reaction of rock

163

minerals with water. For example, carbon dioxide may be produced by reaction between pyrite,

164

dolomite, water and oxygen48,49:

165

() +

 



() + ( ) +    → ()() +  +  + 2  + 2

166

(1)

167

The possibility of pyrite oxidation reaction will be discussed in section 3.3.

168

Under degassed conditions, the gas existing in the pore space dissolves into water, while under

169

degassed conditions, the gas cannot completely dissolve in the imbibing water. For the oxic water

170

imbibition, some of the gas in pore space may be trapped and reduce water saturation and relative

171

permeability.12 The trapped gas during oxic water imbibition may also increase the resistance against

172

the capillary force for water imbibition which can be referred to as relative permeability effect.50

173

According to Figure 2, the imbibition rate of degassed water is higher than that of oxic water. For

174

example, for the EV-1 sample, the degassed water imbibition rate for the initial 18.5 hours is about

175

0.07 g/hr which is higher than that for the oxic water imbibition rate (0.05 g/hr).

176

When the gravity effect is ignored, the imbibed mass is proportional to the square root of time:51  = !"

177 178 179

#$ %%&' ∅)* + ,* -

. × √1

(2)

Here,  is the mass of imbibed water, 2 and 3 are viscosity and density of the imbibed water respectively, is the water saturation of shale samples at time 1, 45 is the capillary pressure at the

180

saturation of , 678 is water relative permeability, 6 and ∅ are absolute permeability and porosity

181

of the shale samples, respectively.

182

Figure 4 shows the normalized mass of imbibed water versus square root of time. Compared with the

183

degassed tests, the lower slope for oxic water imbibition can be due to the potential trapped gas in the 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

184

pores of the rock. More specifically, if the gas traps in pore space, krw reduces due to the declined

185

water saturation. The potential decline of krw can be responsible for the lower imbibition rates during

186

oxic water imbibition compared with that during degassed water imbibition.

187

Figure 4. Mass of imbibed oxic and degassed water versus square root of time for (a) EV-1, (b) EV-2,

188

and (c) EV-3 samples. The trendlines show the initial imbibition rate.

189

During degassed water imbibition, dissolution of gas into the water may provide additional pore space

190

for water imbibition, which explains the difference in imbibed mass observed in Figure 4. On the other

191

hand, the gas trapped in pore space during oxic water imbibition may partially block the path for water

192

imbibition. In Figure 4c, the equilibrated imbibed mass of degassed water is up to 2 times higher than

193

that of oxic water in EV-3 samples which suggests that up to 50% of the water-accessible pore volume

194

is occupied by trapped gas during oxic water imbibition.

195

3.3 Ion Concentration Profiles

196

In this section, ion concentrations in oxic and degassed water during the imbibition process are

197

presented. We measure concentration of major ions such as potassium and sodium and redox sensitive

198

ions such as iron and sulfate to investigate the chemical reactions promoted by the presence of

199

dissolved oxygen in water.

10

ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

200

Energy & Fuels

3.3.1 Potassium and Sodium Ions

201

202

203 204

Figure 5. Potassium (a, c, e) and sodium (b, d, f) concentration of degassed/oxic water in contact with

205

EV-1 (a, b), EV-2 (c, d), and EV-3 (e, f) shale samples measured by ICP-MS.

206

According to Figure 5, potassium and sodium ion concentrations increase during both oxic and

207

degassed water imbibition experiments and no considerable difference is observed in the concentration 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

208

profile of these ions. The source of potassium and sodium can be 1) dissolution of precipitated salts

209

(such as NaCl and KCl) in shale samples; 2) dissolution of potassium or sodium-bearing components

210

in the shale samples such as halite52 and 3) leaching of exchangeable potassium and sodium ions from

211

interlayer of clays.53,54 The similarity of ion concentration profiles between oxic and degassed water

212

imbibition suggests that the presence of oxygen in the imbibition water has negligible impact on the

213

concentrations of produced Na+ and K+ ions.

214

To investigate the potential sources of produced ions during the imbibition experiments, we plot the

215

molar ratio of K+/Cl- and Na+/Cl- in Figure 6. Both K+/Cl- and Na+/Cl- molar ratios are more than one.

216

Leaching of exchangeable cations from the interlayer of clay minerals is a possible reason for the

217

excess produced Na+ and K+.55 Both K+/Cl- and Na+/Cl- molar ratios show an early increase for oxic

218

water. This can be due to the fact that chloride can only be oxidized (not reduced).56 Thus, in an oxic

219

environment, the leached K+ and Na+ from the interlayer of clay minerals may be responsible for the

220

early increase in K+/Cl- and Na+/Cl- molar ratios. Moreover, chloride ion may oxidize and decrease the

221

concentration of Cl- in the solution. Thus, at early times, both K+/Cl- and Na+/Cl- molar ratios increase

222

in the oxic environment. At later times, the potential dissolution of chloride-bearing minerals (such as

223

NaCl and KCl) can be responsible for the declining trends of K+/Cl- and Na+/Cl- molar ratios.

224

At early times, both K+/Cl- and Na+/Cl- molar ratios are higher for degassed water compared with that

225

for oxic water. On the other hand, the imbibed volume is higher for the degassed conditions compared

226

with that for the oxic conditions (Figure 2). Enhanced imbibition can potentially enhance the

227

accessibility to the clay-rich zones, facilitating the ion-exchange reactions. Thus, the probability of

228

leaching exchangeable cations (i.e., K+ and Na+) from the interlayer of clay minerals is higher in the

229

degassed conditions compared with that in the oxic conditions. Similarly, the possible dissolution of

230

chloride-bearing minerals reduces K+/Cl- and Na+/Cl- molar ratios at later times.

12

ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

231

232

233 234

Figure 6. The K+/Cl- (a, c, e) and Na+/Cl- (b, d, f) molar ratio of degassed/oxic water in contact with

235

EV-1 (a, b), EV-2 (c, d), and EV-3 (e, f) shale samples.

236

3.3.2 Iron and Sulfate Ions

237

Previous studies have shown that the concentration of redox sensitive elements (such as Fe and S) can

238

be influenced by the presence of oxygen.33,57 Figure 7 shows concentration profiles of iron and sulfate 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

239

ions in oxic and degassed imbibition experiments. The concentration of iron ion increases over time

240

during the oxic water imbibition while the concentration of iron is almost zero during the degassed

241

water imbibition. This observation indicates that pyrite oxidation reactions is inhibited in the absence

242

of oxygen. According to Table 1, pyrite is the only source of iron in the shale samples. It should be

243

noted that XRD may not detect other possible sources of iron such as goethite, hematite, or magnetite

244

in our shale samples due to their low concentration. The pyrite reactions with water and oxygen may

245

be a possible reason for the increasing trend of iron concentration profile during oxic water imbibition

246

in Figure 7. SEM images in section 3.4 shows pores in the size of 5 to 200 nm near pyrite minerals.

247

These pores provide an interface at which imbibed water, dissolved oxygen, and pyrite co-exist and

248

react.

249

250

14

ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

251 252

Figure 7. Iron (a, c, e) and sulfate (b, d, f) concentration of degassed and oxic imbibing water in

253

contact with EV-1 (a, b), EV-2 (c, d), and EV-3 (e, f) shale samples measured by ICP-MS and IC.

254

In pyrite oxidation reactions, the oxygen atoms from dissolved oxygen in water may react with the

255

sulfur atoms in pyrite molecules at the rock-water interface 58 and produce ferrous ions (Fe2+)

256

according to the following aqueous reactions:59-62 9

257

() + () +   →  + 2 + 2  

258

2  + 2 + () → 2  +   



(3)

(4)

259

The ferrous ions (Fe2+) can react with dissolved oxygen (Eq. 4) to produce ferric ions (Fe3+ in Eq. 4).

260

Pyrite oxidation rate generally increases with increasing pH and dissolved oxygen concentration.63-65

261

In the imbibition experiments, the pH range for degassed and oxic imbibing water is between 8.0 to

262

9.0 and 7.5 to 8.5, respectively. The higher pH of degassed water can be due to the removal of some

263

dissolved carbon dioxide in degassed water by vacuum.66 The overall pyrite oxidation reaction in

264

alkaline solution is as follows:67,68

265

() +

 



() + 4 → ()() +    + 2 

(5)

266

According to Eq. 5, oxidation of one mole of pyrite produces two moles of sulfate ion. It explains why

267

iron and sulfate concentrations in oxic water are higher than that in degassed water, as observed in

268

Figure 7. The rate of pyrite oxidation (Rpyrite) can be expressed by the rate of sulfate concentration 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

269 270

Page 16 of 28

difference:69 R ?@A =

BCDE7FGHI BG

=

B(J+KL*MNOPQ$ J+KL*M NRSOPQ$ )

(6)

BG

271

However, the difference in sulfate concentration of oxic and degassed imbibition is around 2 orders of

272

magnitude larger than the corresponding difference in Fe2+/Fe3+ concentration. The reacted pyrite will

273

produce either soluble Fe3+ (Eq. 4) or insoluble Fe3+ hydroxide precipitate (Eq. 5). Based on Pourbaix

274

diagram of iron, 70 Fe3+-bearing compounds start to precipitate at low pH values. For example,

275

Fe3+-hydroxide starts to precipitate at pH values of about 1.71 Since the pH range in our experiments

276

is 7.5-9.0, the released iron during pyrite oxidation may ultimately form insoluble Fe3+-bearing

277

compounds which cannot be identified by ICP-MS tests. According to the Pourbaix diagram of iron,70

278

at pH values of 7.5 to 9.0, the iron oxidation products could be Fe2+, Fe(OH)2, Fe2O3·nH2O, Fe3O4,

279

which are listed in Eqs. 3, 7, 8, 10, 11. Highly reactive Fe2+/Fe3+ ion may participate in production of

280

water insoluble iron hydroxide, goethite, hematite or magnetite according to Eqs. 7 to 11:54

281

 + 3  → ()() + 3 

(7)

282

 + 2 → ()()

(8)

283

 + 3 → () +  

(9)

284

2 () →  () +  

(10)

285

()() + 2 () →  () +  

(11)

286

It is worth mentioning that the imbibition experiments are conducted at room temperature,

287

atmospheric pressure and nearly neutral pH, which are different compared with the reservoir

288

conditions. For example, Horn River reservoir gas contains about 9-14% of CO2,72 and the reaction of

289

CO2 with water can potentially produce carbonic acid. It must be noted that although carbonic acid

290

production can reduce the pH, natural buffer systems such as dolomite and calcite may neutralize the

291

pH under reservoir conditions.33 Furthermore, the reservoir temperature in the Horn River Basin is in

292

the range of 135 to 175 ℃ and the reservoir pressure is between 20,000 to 53,000 KPa.45 The 16

ACS Paragon Plus Environment

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

293

different in laboratory and reservoir conditions may lead to difference in “speciation and complexation”

294

of ions, which can alter the pH.

295

3.4 SEM and EDS Images

296

In Figure 8, SEM images of a focus area and the corresponding elemental map of sulfur and iron show

297

that pyrite is widely spread in shale samples. High density of sulfur and iron elements in Figure 8.c

298

and 8.d suggests a pyrite-rich area. Figure 8.b shows the magnified SEM image of the area surrounded

299

by dashed line in Figure 8.a, where many pores with a throat size of 50 to 200 nm exist in the

300

pyrite-rich area. Figure 9 shows a scanning helium ion microscope (SHIM) image of the pyrite area.

301

Two pyrite framboids with the approximate size of 1 µm are shown in Figure 9 and nanopores are

302

observed between the pyrite particles. These pores provide an interface at which imbibed water,

303

dissolved oxygen, and pyrite chemically interact. The interaction between these three components may

304

result in iron-bearing precipitation on the surface of nanopores which may narrow or block the pores

305

and may ultimately lower the permeability of the shale sample.

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

306 307

Figure 8. (a) SEM image of a focused area on a EV-2 shale sample; (b) magnified SEM image

308

showing pores in pyrite; (c) sulfur elemental map and (d) iron elemental map obtained by EDS of the

309

EV-2 shale sample. Darker area in elemental maps represents higher element density.

310 311

Figure 9. SHIM image of the EV-2 shale sample. Nano pores are observed around pyrite framboid. 18

ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

312

Finally, pyrite dissolution and oxidation may have a dual effect: 1) pyrite dissolution may create more

313

pore space and may enhance water imbibition and 2) precipitation of iron-bearing compounds formed

314

as a result of pyrite oxidation in the pore space may reduce the flow conductivity for water imbibition.

315

Reduced imbibition rate and lower final imbibed water mass under the oxic conditions in this study

316

suggest that precipitation of iron-bearing compounds formed by pyrite oxidation in the pore space has

317

a more pronounced effect in comparison with the pore space creation by pyrite dissolution.

318

4. Conclusions

319

Imbibition experiments using degassed and oxic water were conducted on samples from Evie

320

Formation, a shale member of the Horn River Basin. Concentration of different ions in imbibing water

321

were measured to investigate the effect of dissolved oxygen on water imbibition and produced ions.

322

Pore network and minerals of the shale samples were investigated using SEM-EDS imaging. The key

323

results of this study are summarized as follows:

324



Degassed water can dissolve the gas in the shale pore space, contributing to a higher

325

imbibition rate and final imbibed mass, compared with the oxic water. This observation

326

indicates that the initially-degassed conditions in reservoir may accelerate the rate of water

327

imbibition. Less gas entrapment under degassed conditions promote the water imbibition into

328

shale matrix.

329



The oxygen dissolved in water can react with pyrite-bearing components in shale and produce

330

iron and sulfate ions. Interactions of degassed water with shale do not result in iron

331

production. This observation indicates that dissolved oxygen in fracturing fluid can cause

332

pyrite oxidation and dissolution.

333



Pores observed by SEM are in the vicinity of both pyrite framboids and pyrite crystals.

334

Although pyrite constitutes only a small portion of shale minerals, the oxygen-water-pyrite

335

interactions cannot be ignored when interpreting the chemistry of water produced during

336

flowback processes. 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

337

Acknowledgements

338

The authors are grateful to (1) Nexen Energy ULC, INPEX Gas British Columbia Ltd and Natural

339

Sciences and Engineering Research Council (NSERC) for supporting this project, (2) Doug Bearinger

340

for useful discussions, and (3) Todd Kinnee for his help in conducting the experiments.

341

References

( 1 ) Harper J. 2008. The Marcellus Shale – An Old ― New Gas Reservoir in Pennsylvania. Pennsylvania Geology, 28(1). (2) Soeder, D. J., Kappel, W. M. 2009. Water resources and natural gas production from the Marcellus Shale (pp. 1-6). Reston, Virginia: US Department of the Interior, US Geological Survey. (3) Alessi, D. S., Zolfaghari, A., Kletke, S., Gehman, J., Allen, D. M., Goss, G. G. 2017. Comparative analysis of hydraulic fracturing wastewater practices in unconventional shale development: Water sourcing, treatment and disposal practices. Canadian Water Resources Journal/Revue canadienne des ressources hydriques, 42(2), 105-121. (4) Cherry, J., Ben-Eli, M., Bharadwaj, L., Chalaturnyk, R., Dusseault, M. B., Goldstein, B., et al. 2014. Evironmental impacts of shale gas extraction in Canada--The expert panel on harnessing science and technology to understand the environmental impacts of shale gas extraction. Council of Canadian Academies. (5) Cho, Y., Ozkan, E., Apaydin, O. G. 2013. Pressure-dependent natural-fracture permeability in shale and its effect on shale-gas well production. SPE Reservoir Evaluation & Engineering, 16(02), 216-228. (6) Fakcharoenphol, P., Torcuk, M., Kazemi, H., Wu, Y. S. 2016. Effect of shut-in time on gas flow 20

ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

rate in hydraulic fractured shale reservoirs. Journal of Natural Gas Science and Engineering, 32, 109-121. (7) Asadi, M., Woodroof, R.A. and Himes, R.E., 2008. Comparative study of flowback analysis using polymer concentrations and fracturing-fluid tracer methods: a field study. SPE Production & Operations, 23(02), pp.147-157. (8) Xu, Y., Adefidipe, O. A., Dehghanpour, H. 2015. Estimating fracture volume using flowback data from the Horn River Basin: A material balance approach. Journal of Natural Gas Science and Engineering, 25, 253-270. (9) Zhou, Q., Dilmore, R., Kleit, A. and Wang, J.Y., 2016. Evaluating fracture-fluid flowback in Marcellus using data-mining technologies. SPE Production & Operations, 31(02), pp.133-146. (10) Neuzil, C. E. 2000. Osmotic generation of 'anomalous' fluid pressures in geological environments. Nature, 403(6766), 182. (11) Fakcharoenphol, P., Kurtoglu, B., Kazemi, H., Charoenwongsa, S., Wu, Y. S. 2015. The Effect of Chemical Osmosis on Oil and Gas Production from Fractured Shale Formations. Fluid Dynamics in Complex Fractured-Porous Systems, 210, 85. (12) Bennion, D. B., Thomas, F. B., Bietz, R. F., Bennion, D. W. 1996. Water and hydrocarbon phase trapping in porous media-diagnosis, prevention and treatment. Journal of Canadian Petroleum Technology, 35(10). (13) Sharma, M., Agrawal, S. 2013. Impact of liquid loading in hydraulic fractures on well productivity. In SPE hydraulic fracturing technology conference. Society of Petroleum Engineers. (14) Bertoncello, A., Wallace, J., Blyton, C., Honarpour, M., Kabir, C. S. 2014, February. Imbibition and water blockage in unconventional reservoirs: well management implications during flowback and early production. In SPE/EAGE European Unconventional Resources Conference and Exhibition.

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

(15) Cheng, Y. 2012. Impact of water dynamics in fractures on the performance of hydraulically fractured wells in gas-shale reservoirs. Journal of Canadian Petroleum Technology, 51(02), 143-151. (16) Yan, Q., Lemanski, C., Karpyn, Z. T., Ayala, L. F. 2015. Experimental investigation of shale gas production impairment due to fracturing fluid migration during shut-in time. Journal of Natural Gas Science and Engineering, 24, 99-105. (17) Chakraborty, N., Karpyn, Z. T. 2015. Gas Permeability Evolution with Soaking Time in Ultra Tight Shales. In SPE Annual Technical Conference and Exhibition. Society of Petroleum Engineers. (18) Zhou, Z., Abass, H., Li, X., Teklu, T. 2016. Experimental investigation of the effect of imbibition on shale permeability during hydraulic fracturing. Journal of Natural Gas Science and Engineering, 29, 413-430. (19) Paktinat, J., Pinkhouse, J. A., Johnson, N., Williams, C., Lash, G. G., Penny, G. S., Goff, D. A. 2006. Case Studies: Optimizing Hydraulic Fracturing Performance in Northeastern Fractured Shale Formations, SPE 104306-MS, SPE Production and Operations. SPE Production and Operations. (20) Bahrami, H., Rezaee, R., Clennell, B. 2012. Water blocking damage in hydraulically fractured tight sand gas reservoirs: An example from Perth Basin, Western Australia. Journal of Petroleum Science and Engineering, 88, 100-106. (21) Ji, L., Geehan, T. 2013. Shale failure around hydraulic fractures in water fracturing of shale gas. In SPE Unconventional Resources Conference Canada. Society of Petroleum Engineers. (22) Ghanbari, E., Dehghanpour, H., 2015. Impact of Rock Fabric on Water Imbibition and Salt Diffusion in Gas Shales. International Journal of Coal Geology 138: 55–67. (23) Chenevert, M. E. 1970. Shale alteration by water adsorption. Journal of petroleum technology, 22(09), 1-141. (24) van Oort, E., Hale, A. H., Mody, F. K., Roy, S. 1996. Transport in shales and the design of 22

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

improved water-based shale drilling fluids. SPE Drilling & Completion, 11(03), 137-146. (25) Lal, M. 1999. Shale stability: drilling fluid interaction and shale strength. In SPE Asia Pacific Oil and Gas Conference and Exhibition. Society of Petroleum Engineers. (26) Tang, G. Q., Morrow, N. R. 1997. Salinity, temperature, oil composition, and oil recovery by waterflooding. SPE Reservoir Engineering, 12(04), 269-276. (27) Høgnesen, E. J., Standnes, D. C., Austad, T. 2006. Experimental and numerical investigation of high temperature imbibition into preferential oil-wet chalk. Journal of Petroleum Science and Engineering, 53(1), 100-112. (28) Strand, S., Puntervold, T., Austad, T. 2008. Effect of temperature on enhanced oil recovery from mixed-wet chalk cores by spontaneous imbibition and forced displacement using seawater. Energy & Fuels, 22(5), 3222-3225. (29) Santos, H., Diek, A., Roegiers, J. C., Fontoura, S. 1996. Can shale swelling be (easily) controlled?. In ISRM International Symposium-EUROCK 96. International Society for Rock Mechanics. (30) Wang, D., Butler, R., Liu, H., Ahmed, S. 2011. Flow-rate behavior and imbibition in shale. SPE Reservoir Evaluation & Engineering, 14(04), 485-492. (31) Zhou, Z., Hoffman, B. T., Bearinger, D., Li, X. 2014. Experimental and numerical study on spontaneous imbibition of fracturing fluids in shale gas formation. In SPE/CSUR Unconventional Resources Conference–Canada. Society of Petroleum Engineers. (32) Chenevert, M. E. K., Osisanya, S. O. K. 1992. Shale swelling at elevated temperature and pressure. In The 33th US Symposium on Rock Mechanics (USRMS). American Rock Mechanics Association. (33) Zolfaghari, A., Dehghanpour, H., Noel, M., Bearinger, D. 2016. Laboratory and field analysis of 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

flowback water from gas shales. Journal of Unconventional Oil and Gas Resources, 14, 113-127. (34) Roshan, H., Ehsani, S., Marjo, C. E., Andersen, M. S., Acworth, R. I. 2015. Mechanisms of water adsorption into partially saturated fractured shales: An experimental study. Fuel, 159, 628-637. (35) Bearinger, D. 2013. Message in a Bottle. In Unconventional Resources Technology Conference (pp. 1455-1460). Society of Exploration Geophysicists, American Association of Petroleum Geologists, Society of Petroleum Engineers. (36) Curtis, J. B. 2002. Fractured shale-gas systems. AAPG bulletin, 86(11), 1921-1938. (37) Evangelou, V. P., Zhang, Y. L. 1995. A review: pyrite oxidation mechanisms and acid mine drainage prevention. Critical Reviews in Environmental Science and Technology, 25(2), 141-199. (38) Hutcheon, I. 1998. The potential role of pyrite oxidation in corrosion and reservoir souring. Journal of Canadian Petroleum Technology, 37(01). (39) Elie, M., Faure, P., Michels, R., Landais, P., Griffault, L. 2000. Natural and laboratory oxidation of low-organic-carbon-content sediments: comparison of chemical changes in hydrocarbons. Energy & Fuels, 14(4), 854-861. (40) Marynowski, L., Szełęg, E., Jędrysek, M. O., Simoneit, B. R. 2011. Effects of weathering on organic matter: Part II: Fossil wood weathering and implications for organic geochemical and petrographic studies. Organic Geochemistry, 42(9), 1076-1088. (41) Marriott, R. A., Pirzadeh, P., Marrugo-Hernandez, J. J., Raval, S. 2015. Hydrogen sulfide formation in oil and gas. Canadian Journal of Chemistry, 94(4), 406-413. (42) Jin, L., Mathur, R., Rother, G., Cole, D., Bazilevskaya, E., Williams, J., Brantley, S. 2013. Evolution of porosity and geochemistry in Marcellus Formation black shale during weathering. Chemical Geology, 356, 50-63. (43) Chen, Q., Kang, Y., You, L., Yang, P., Zhang, X., Cheng, Q. 2017. Change in composition and 24

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

pore structure of Longmaxi black shale during oxidative dissolution. International Journal of Coal Geology, 172, 95-111. (44) Rowan, E. L., Engle, M. A., Kraemer, T. F., Schroeder, K. T., Hammack, R. W., Doughten, M. W. 2015. Geochemical and isotopic evolution of water produced from Middle Devonian Marcellus shale gas wells, Appalachian basin, Pennsylvania. AAPG Bulletin, 99(2), 181-206. (45) B. C. Oil & Gas Commission. 2014. Horn River Basin unconventional shale gas play atlas. Report, BC Oil and Gas Commission. (46) Kolev, N. I. 2011. Solubility of O2, N2, H2 and CO2 in water. In Multiphase Flow Dynamics 4 (pp. 209-239). Springer Berlin Heidelberg. (47) Dehghanpour, H., Zubair, H. A., Chhabra, A., Ullah, A. 2012. Liquid intake of organic shales. Energy & Fuels, 26(9), 5750-5758. (48) Descourvières, C., Hartog, N., Patterson, B. M., Oldham, C., Prommer, H. 2010. Geochemical controls on sediment reactivity and buffering processes in a heterogeneous aquifer. Applied Geochemistry, 25(2), 261-275. (49) Chermak, J. A., Schreiber, M. E. 2014. Mineralogy and trace element geochemistry of gas shales in the United States: Environmental implications. International Journal of Coal Geology, 126, 32-44. (50) Jamin, J. M. 1860. Memoir on equilibrium and movement of liquids in porous substances. Compt. Rend, 50, 172-176. (51) Handy L. L. 1960. Determination of effective capillary pressures for porous media from imbibition data. Trans. AIME, 219, 75-80. (52) Blauch, M. E., Myers, R. R., Moore, T., Lipinski, B. A., Houston, N. A. 2009. Marcellus shale post-frac flowback waters-Where is all the salt coming from and what are the implications? In SPE Eastern Regional Meeting. Society of Petroleum Engineers. 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

(53) Hensen, E. J., Smit, B. 2002. Why clays swell. The Journal of Physical Chemistry B, 106(49), 12664-12667. (54) Essington, M. E. 2015. Soil and water chemistry: an integrative approach. CRC press. (55) Keller, W. D., da Costa, L. M. 1989. Comparative chemical compositions of aqueous extracts from representative clays. American Mineralogist, 74(9-10), 1142-1146. (56) Huang, L., Wojciechowski, G., Ortiz de Montellano, P. R. 2005. Prosthetic heme modification during halide ion oxidation. Demonstration of chloride oxidation by horseradish peroxidase. Journal of the American Chemical Society, 127(15), 5345-5353. (57) Xie, X., Ellis, A., Wang, Y., Xie, Z., Duan, M., Su, C. 2009. Geochemistry of redox-sensitive elements and sulfur isotopes in the high arsenic groundwater system of Datong Basin, China. Science of the total environment, 407(12), 3823-3835. (58) Biegler, T., Swift, D. A. 1979. Anodic behaviour of pyrite in acid solutions. Electrochimica Acta, 24(4), 415-420. (59) Singer, P. C., Stumm, W. 1968. Kinetics of the oxidation of ferrous iron. In Second Symposium Coal Mine Drainage Research (pp. 12-34). (60) Bailey, L. K., Peters, E. 1976. Decomposition of pyrite in acids by pressure leaching and anodization: the case for an electrochemical mechanism. Canadian Metallurgical Quarterly, 15(4), 333-344. (61) Davison, W., Seed, G. 1983. The kinetics of the oxidation of ferrous iron in synthetic and natural waters. Geochimica et Cosmochimica Acta, 47(1), 67-79. (62) McKibben, M. A. 1985. Kinetics of aqueous oxidation of pyrite by ferric iron, oxygen, and hydrogen-peroxide from pH 1-4 and 20-40oC (Doctoral dissertation). (63) Moses, C. O., Nordstrom, D. K., Herman, J. S., Mills, A. L. 1987. Aqueous pyrite oxidation by 26

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

dissolved oxygen and by ferric iron. Geochimica et Cosmochimica Acta, 51(6), 1561-1571. (64) Holmes, P. R., Crundwell, F. K. 2000. The kinetics of the oxidation of pyrite by ferric ions and dissolved oxygen: an electrochemical study. Geochimica et Cosmochimica Acta, 64(2), 263-274. (65) Manaka, M. 2009. Comparison of rates of pyrite oxidation by dissolved oxygen in aqueous solution and in compacted bentonite. Journal of mineralogical and petrological sciences, 104(2), 59-68. (66) Millero, F. J. 1995. Thermodynamics of the carbon dioxide system in the oceans. Geochimica et Cosmochimica Acta, 59(4), 661-677. (67) Nicholson, R. V., Gillham, R. W., Reardon, E. J. 1988. Pyrite oxidation in carbonate-buffered solution: 1. Experimental kinetics. Geochimica et Cosmochimica Acta, 52(5), 1077-1085. (68) Ciminelli, V. S. T., Osseo-Asare, K. 1995. Kinetics of pyrite oxidation in sodium hydroxide solutions. Metallurgical and Materials Transactions B, 26(4), 677-685. (69) Moses, C. O., Herman, J. S. 1991. Pyrite oxidation at circumneutral pH. Geochimica et Cosmochimica Acta, 55(2), 471-482. (70) Beverskog, B., Puigdomenech, I. 1996. Revised Pourbaix diagrams for iron at 25–300 C. Corrosion Science, 38(12), 2121-2135. (71) Taylor, K. C., Nasr-El-Din, H. A., Al-Alawi, M. J. 1999. Systematic study of iron control chemicals used during well stimulation. SPE Journal, 4(01), 19-24. (72) Reynolds, M. M., Munn, D. L. 2010. Development update for an emerging shale gas giant field-Horn River Basin, British Columbia, Canada. In SPE Unconventional Gas Conference. Society of Petroleum Engineers.

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 28

342

Appendix

343

Table A. Measured sample mass of oxic and degassed imbibition experiments.

Time Hours 0.00 1.30 2.58 9.83 18.50 27.50 44.72 72.83 94.83 119.00 147.42 172.33 194.50 218.58 242.67 265.67 361.67 457.67 505.67 553.67 625.67

Measured Sample Mass, g EV-1

EV-2

EV-3

Oxic

Degassed

Oxic

Degassed

Oxic

Degassed

91.90 92.06 92.15 92.29 92.37 92.39 92.44 92.49 92.55 92.55 92.60 92.60 92.61 92.63 92.65 92.65 92.69 92.71 92.72 92.73 92.73

96.87 97.21 97.33 97.46 97.57 97.62 97.66 97.69 97.74 97.76 97.76 97.79 97.81 97.81 97.82 97.84 97.84 97.84 97.85 97.86 97.86

98.85 99.43 99.60 99.82 99.93 100.00 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02 100.02

101.61 102.47 102.73 103.08 103.17 103.17 103.18 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20 103.20

112.14 112.33 112.36 112.39 112.43 112.45 112.48 112.48 112.52 112.54 112.54 112.54 112.55 112.55 112.56 112.56 112.56 112.54 112.56 112.56 112.56

108.40 109.00 109.11 109.16 109.22 109.25 109.26 109.28 109.30 109.33 109.33 109.34 109.35 109.35 109.35 109.35 109.34 109.36 109.34 109.35 109.33

344

28

ACS Paragon Plus Environment