Effect of Pore Size Distribution on Dissociation Temperature

Mar 23, 2018 - Oil & Gas Research Center, Korea Institute of Geoscience and Mineral ... Advanced in Situ Combustion Enhanced Oil Recovery Processes...
1 downloads 0 Views 5MB Size
Subscriber access provided by Queen Mary, University of London

Article

Effect of pore size distribution on dissociation temperature depression and phase boundary shift of gas hydrate in various fine-grained sediments Taehyung Park, Joo Yong Lee, and Tae-Hyuk Kwon Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.8b00074 • Publication Date (Web): 23 Mar 2018 Downloaded from http://pubs.acs.org on March 24, 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.

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 46 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

Effect of pore size distribution on dissociation temperature

2

depression and phase boundary shift of gas hydrate in various

3

fine-grained sediments

4

5

6 7 8 9

Taehyung Park1, Joo Yong Lee2 and Tae-Hyuk Kwon1* 1

Department of Civil and Environmental Engineering, Korea Advanced Institute of Science and

Technology (KAIST), Daejeon, 34141, Korea 2

Oil & Gas Research Center, Korea Institute of Geoscience and Mineral Resources (KIGAM),

Daejeon, 34132, Korea.

10 11

Corresponding Author

12

* Tae-Hyuk Kwon, Ph.D, Corresponding Author

13

Assistant Professor, Department of Civil and Environmental Engineering,

14

Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Korea

15

Tel: +82-42-350-3628

16

Email: [email protected]

Fax: +82-42-869-3610

17 18 19 20 21 22

Number of words: Number of tables: Number of figures:

4633 (without References, Figures or Tables) 1 9

23

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

Abstract

25

Capillarity in small, confined pores has a pronounced effect on depression of dissociation

26

temperature of gas hydrates, known as the Gibbs-Thomson effect. However, this effect remains

27

poorly understood in natural fine-grained sediments with wide pore size distributions. This study

28

investigated the effect of pore size distributions of fine-grained sediments on the dissociation

29

temperature of gas hydrate. Gas hydrate was synthesized under partially water-saturated

30

conditions in nano-sized silica gels and in various natural fine-grained sediment samples,

31

including sand, silt, diatoms, a diatom-sand mixture, and clayey sediment. The synthesized

32

hydrate samples were thermally dissociated under isochoric conditions, while the melting

33

temperature depression and the shifted phase boundaries were monitored. We observed the

34

dissociation temperature depression of approximately 0.1–0.3°C in silt, 0.2–0.4°C in the diatom

35

sample, and 1.2–1.5°C in clayey silt, while no temperature depression was observed in sand. In a

36

particular case of diatom-sand mixture, the dual porosity condition with the submicron-scale

37

internal pores of diatoms and the macro-pores of sands, rendered dual phase boundaries, one

38

with a ~0.4°C temperature depression and one with no depression, respectively. Despite the wide

39

ranges of pore size, gas hydrates were preferentially formed in smaller pores; which comprise

40

less than 40% of the cumulative pore volumes. This was because the initial water loci

41

exacerbated the Gibbs-Thomson effect in partially water-saturated conditions. Our results

42

provide clear experimental evidence on and novel insights into the effect of pore size

43

distributions of fine-grained sediments on the dissociation behavior and phase boundaries of gas

44

hydrates, both in the presence of free gas, or in water-limiting conditions that exhibit a

45

considerable Gibbs-Thomson effect.

46

2

ACS Paragon Plus Environment

Page 2 of 46

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

47

Introduction

48

Gas hydrates are ice-like crystalline structures comprised of hydrogen-bonded water cages

49

containing gas molecules.1 Natural gas hydrates are considered to be potential energy resources

50

due to the vast amount accumulated on continental margins and in permafrost regions.2-4 At the

51

same time, they act as a carbon sink in the global carbon cycle,5 and as potential triggers and/or

52

accelerators of geological hazards.6-9 Furthermore, the natural occurrence of carbon dioxide

53

(CO2) hydrate in CO2 reservoirs in the deep ocean10 has aroused interest in a gas hydrate

54

production method involving a CH4-CO2 replacement; gas hydrates are a possible solution to the

55

problem of geologic CO2 storage.11-14 In all cases, gas hydrate dissociation reactions are

56

unavoidably involved, understanding of the dissociation behavior of gas hydrates in natural

57

sediments is thus essential. Accordingly, there have been extensive studies to understand the

58

chemo-thermo-mechanical changes during the dissociation of gas hydrates in natural sediments,

59

such as the phase changes,15, 16 pore water chemistry changes,17 and heat and mass transfer.18-20

60

The majority of natural gas hydrates are reportedly hosted in fine-grained sediments, in locations

61

such as Blake Ridge, the eastern coast of India (NGHP), Cascadia Margin, offshore Peru,

62

Okushiri Ridge, Ulleung Basin (UB), Orca Basin, Sea of Okhotsk, Atwater Valley in Gulf of

63

Mexico, and Hydrate Ridge.21-24 Such fine-grained sedimentary environments affect the

64

formation and dissociation of gas hydrates in various ways, through the influence of negatively

65

charged clay surfaces,25, 26 the release of adsorbed cations into pore water,27, 28 and the capillarity

66

generated by small pore sizes.25, 29-34 In particular, gas hydrates that are confined in small pores

67

are influenced by the capillary effect, which shifts the thermodynamic phase boundary to lower

68

temperatures (or higher pressures). This phenomenon is called melting point depression (or

69

dissociation temperature depression), often referred to as the Gibbs-Thomson effect.16, 29-33, 35-38 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

70

Water molecules in confined small pores have a lower activity than unconfined molecules; their

71

lowered activity requires greater thermodynamic driving forces to form and stabilize clathrate

72

gas hydrates. Therefore, gas hydrates in small pores are stable under higher pressures and lower

73

temperatures than those found in large pores. In turn, they dissociate at lower temperatures,

74

having the phase boundary shifted. As a result, it has been reported that submicron sediment pore

75

sizes cause the uplift of the base of the methane hydrate stability zone (MHSZ), and decrease its

76

thickness.39, 40 The prediction of the temperature depression of the hydrate phase boundary in

77

clay-rich sediments is critical to understanding the formation and dissociation of gas hydrates in

78

those sediments, as their small pore size affects the phase boundary, location of the base of the

79

MHSZ, hydrate morphology, and properties of hydrate-bearing sediments.

80

The dissociation temperature depression and phase boundary shift of gas hydrates in small pores

81

have been studied by using artificially synthesized porous silica gels, due to their well-controlled

82

and narrow range of pore diameters.29-32, 35-37 Natural fine-grained sediments can have a wide

83

pore size distribution over more than one order of magnitude. This is most apparent when the

84

sediment is composed of various particles, including sand, silt, and clay-sized particles, or where

85

coarse-grained sediments contain diatoms, which have their own sub-micron scale internal

86

pores.21,

87

offshore Japan, showed pore size distributions ranging from 10 nm to 1 µm and ~10 nm to 10

88

µm, respectively.40, 41 However, a limited number of studies have been conducted using natural

89

fine-grained sediments that have a wide pore size distribution; the dissociation behavior of gas

90

hydrates in such sediments remains poorly examined.

22

For instance, the sediments cored from the Kumano Basin and Nankai Trough,

4

ACS Paragon Plus Environment

Page 4 of 46

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

91

In this study, we investigated the Gibbs-Thomson effect and the dissociation temperature

92

depression of CO2 hydrate in fine-grained sediments with pore sizes ranging over two orders of

93

magnitude. CO2 hydrate was synthesized under water-limiting conditions in various sediment

94

samples, including porous silica gels, sand, silt, diatoms, a diatom-sand mixture, and clayey

95

sediment from Ulleung Basin, offshore Korea. Then, the hydrates in these sediments were

96

thermally dissociated while the pressure and temperature traces were monitored. The dissociation

97

temperature depression was identified and correlated to the pore size distribution of the host

98

sediments using the Gibbs-Thomson model. The implication of the observations and results are

99

discussed with respect to the occurrence of hydrate in the presence of free gas, the base of

100

MHSZ, and the dissociation behavior during hydrate production.

101

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

102

2

Materials and Methods

103

2.1

Description of the sediment samples used

104

In this study, two artificial porous samples and five natural sediment samples were used. We first

105

prepared two porous silica gels (SG; Sigma-Aldrich, St. Louis, MO, USA) with a narrow range

106

of pore sizes and mean pore diameters of 6 nm and 15 nm, hereafter denoted as 6SG and 15SG,

107

respectively. These control samples were chosen to validate the Gibbs-Thomson model and the

108

corresponding input parameters to determine the melting temperature shift. Five natural

109

sediments were prepared: fine quartz sand (Ottawa F110, U.S. Silica, Frederick, MD, USA),

110

crushed silica silt (Sil-Co-Sil, U.S. Silica, Frederick, MD, USA.), diatoms (Celite 545, Sigma-

111

Aldrich, St. Louis, MO, USA), a mixture of diatom and sand at a 1:9 mass ratio, and a fine-

112

grained sediment core (UBGH1-10B-2b) retrieved from hydrate deposits in the Ulleung Basin,

113

East Sea, offshore Korea, hereafter referred to as sand, silt, diatoms, sand-diatom mixture, and

114

UB sediment, respectively. The basic index and physical properties of the prepared samples,

115

including mean particle diameter (D50), specific surface area (Sa), and porosity, are listed in

116

Table 1. More details of the sand and the silt samples can be found in Kwon et al.,38 and details

117

of the UB sediment are available in Kim et al.42

118

The particle size distributions of the sediment samples were analyzed using a laser light

119

scattering analyzer (Saturn DigiSizer 5200, Micromeritics, Norcross, GA, USA), as shown in

120

Figure 1. The specific surface areas were measured using the nitrogen (N2) adsorption dry

121

method (TriStar II 3020, Micromeritics, Norcross, GA, USA) and the methylene blue adsorption

122

wet method43. Pore size distributions of the samples were estimated using mercury intrusion

123

porosimetry (MIP; MicroActive AutoPore V 9600, Micromeritics, Norcross, GA, USA) for

6

ACS Paragon Plus Environment

Page 6 of 46

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

124

meso-to-macro pores larger than 10 nm, and the Barrett-Joyner-Halenda (BJH) N2

125

adsorption/desorption method for nano-sized pores smaller than 20 nm, as shown in Figure 2.

126

The pore sizes of 6SG and 15SG were obtained from the BJH method (Figure 2a), and the pore

127

size distribution of the natural sediments (silt, diatoms, and UB sediment) were estimated by

128

integrating both the MIP and BJH results to identify the broad range of the pore size

129

distribution.44 The pore size distribution of the sand and sand-diatom mixture was estimated by

130

MIP only, owing to a minimal presence, or lack of, of nano-sized pores in those samples (Figure

131

2b).

132

Samples 6SG and 15SG had mean particle diameters (D50) of approximately 250 µm and 560

133

µm, and mean pore diameters of 6.5 nm and 16.2 nm, respectively. The sand (D50 = 120 µm)

134

represents coarse-grained sediments or sandy sediments frequently found in natural hydrate

135

deposits. The silt (D50 = 28 µm) was chosen to represent silty sediments in formations containing

136

hydrate with a low specific surface area and low plasticity (e.g., Nankai Trough, Japan and

137

Shenhu, South China Sea).45,

138

plasticity, D50 = 4.8 µm, a clay fraction greater than 20%, and a high specific surface area of ~90

139

m2/g (Table 1). The pore sizes of the silt and UB sediment differ by approximately one order of

140

magnitude, with the majority of pores falling into the ranges of 0.1–10 µm and 0.01–1 µm,

141

respectively. The diatom samples (D50 = 66.7 µm) had pores smaller than the sand, and had a low

142

specific surface area, similar to the silt. However, the diatoms showed a pore size distribution of

143

0.1–10 µm, in which two distinctive peaks were observed; one in the range of 100–400 nm,

144

which were of their internal pores, and the other in the range of 1–10 µm, from their particle size.

145

The diatom sample was selected in this study because many hydrate-containing formations in

146

deep oceans are reported to contain diatomaceous sediment.21, 47

46

The UB sediment represents clay-rich sediments with high

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

The morphology and microstructures of the sediment grains were imaged using a scanning

148

electron microscopy (SEM; Magellan 400, FEI Company, Hillsboro, OR, USA). Figure 3 shows

149

some representative images of the samples. The presence of nano-sized pores in 6SG, 15SG, and

150

the diatom sample can be seen in Figures 3a, 3b, and 3h. The SEM images of the UB sediment

151

confirmed that the cored UB sample contained clay minerals (Figure 3f) and diatoms (Figure 3g)

152

with nano-scale intrapores (~10 to 100 nm).

153 154

2.2

Experimental setup

155

In this study, CO2 (research-graded with 99.9% purity, Sam-O Gas Co., Daejeon, Korea) was

156

used to form gas hydrates. CO2 was selected as a suitable analogue for CH4 hydrate because of

157

the similarities between CO2 hydrate and CH4 hydrate, such as their crystal structures and

158

physicochemical characteristics. Both CO2 hydrate and CH4 hydrate form structure I hydrates,

159

and the size of CO2 and CH4 gas molecule are similar (4.36 and 5.12 Å, respectively).1 A

160

cylindrical polycarbonate high pressure column with a height of 80 mm, an outer diameter of 19

161

mm, an inner diameter of 6.4 mm, and an inner volume of 3.18 cm3 was fabricated to synthesize

162

and dissociate CO2 hydrate in the porous samples, as shown in Figure 4. One resistance

163

temperature detector (RTD; PT100; Hankook Electric Heater, Daejeon, Korea) was inserted into

164

the center of the polycarbonate column to monitor the temperature of the samples. In addition,

165

two pressure transducers (PX302; Omega Engineering, INC., Norwalk, CT, USA) were installed

166

to monitor the pressure inside, while the inlet (CO2 gas injection part) and outlet (CO2 gas

167

purging part) pressures were monitored to confirm that they were in a similar range. The pressure

168

inside the column, the CO2 gas pressure in particular, was regulated by a pressurized CO2

169

cylinder. The temperature was controlled by submerging the column in a temperature-controlled 8

ACS Paragon Plus Environment

Page 8 of 46

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

170

bath (RW3-2025; Jeio Tech, Daejeon, Korea). The temperature and pressure inside the column

171

were logged with a data acquisition unit (34970A; Keysight Technologies, Santa Rosa, CA,

172

USA) over the course of the hydrate formation and dissociation experiments.

173

174

2.3

Experimental procedures

175

Pre-wetted sediment sample preparation. All samples were pre-wetted with deionized water

176

(DIW) at a pre-determined water content to achieve partially water-saturated conditions. Then,

177

the pre-wetted samples were hand-tamped in the column to be approximately 8 cm high. The

178

porosities of the packed samples are listed in Table 1.

179

CO2 hydrate formation. The column was flushed with pure CO2 gas several times to remove the

180

residual air inside. CO2 gas was then injected and pressurized to 3 MPa at 8°C. Thereafter, the

181

temperature of the column was lowered to within the hydrate thermodynamic stability zone. The

182

temperature was lowered to −2°C for the 6SG and 15SG samples, and to 1°C for the other

183

samples. Upon confirmation of hydrate nucleation, exhibited by an exothermal temperature rise,

184

the pressure and temperature were maintained for more than 24 h to allow sufficient time for

185

CO2 hydrate formation and growth.

186

Stepwise thermal dissociation of CO2 hydrate. After the completion of the hydrate formation

187

process, the temperature of the column was increased in steps of 1°C every 5 h, while the

188

constant volume condition was maintained with no mass flux (i.e., isochoric heating). When the

189

temperature exceeded the dissociation temperature (or equilibrium temperature), CO2 hydrates

190

started to melt and release free CO2 gas, thus increasing the pressure. Such step-wise heating

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

191

continued until no distinct pressure increase was observed, which indicated the complete

192

dissociation of CO2 hydrate, referred to as the first heating cycle. The temperature and pressure

193

were logged every 10 s over the course of the experiments. As a result, the pressure and

194

temperature traces obtained during hydrate dissociation indicate the thermodynamic phase

195

equilibrium of the CO2 hydrate in the sample. Again, the reaction column was cooled to form

196

CO2 hydrate (i.e., cooling cycle). To confirm the consistency and repeatability of the experiment

197

results, the aforementioned, step-wise thermal dissociation procedure was repeated, referred to as

198

the second heating cycle. Figure 5 shows the typical temperature and pressure traces obtained

199

from the first heating, through cooling, and to the second heating cycle.

200

10

ACS Paragon Plus Environment

Page 10 of 46

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

201

3

Results and Analyses

202

3.1

Results of porous silica glass: Validation of Gibbs-Thomson model

203

Natural, gas hydrate-bearing clayey sediments are expected to have nano-scale pores that affect

204

the formation and dissociation behavior of gas hydrates. The capillarity exerted by small pores

205

hinders water molecules from associating with other water molecules, hence requiring lower

206

temperatures, or higher pressures, for hydrate formation. This can be represented by the reduced

207

water activity, which is a measure of the energy state (i.e., chemical potential) of the water

208

molecules. Therefore, to identify the effects of nano-scale pore size on gas hydrate dissociation,

209

CO2 hydrate stability tests with nano-sized porous silica gels (6SG, 15SG) were performed to

210

validate the Gibbs-Thomson model in phase boundary calculations.

211

Figure 6a shows the pressure-temperature (PT) traces of 6SG and 15SG, which are shifted

212

toward lower temperatures compared to the bulk phase boundary (when the gas hydrates are

213

unconfined). This confirms the depression of the CO2 hydrate dissociation temperature by small

214

pores of the silica gels, as corroborated by many previous studies.29, 31, 33 Temperature shifts

215

(∆Tdep) of ~5°C and ~2°C were observed for the 6SG and 15SG samples, respectively (Figure

216

6a). The PT traces of the first and second heating cycles were found to be consistent, as both

217

traces were in line with each other and relatively parallel to the bulk phase boundary (Figure 6b).

218

The dissociation temperature depression (∆Tdep) can be modeled using the Gibbs-Thomson

219

equation, as follows.31, 33, 38:

220

2  γ m cos θ ∆Tdep = −  hw h d  ρ h 0 L f

  ⋅ Tbulk 

for cylindrical shaped gas hydrates

11

ACS Paragon Plus Environment

(1)

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

221

and

222

4  γ m cos θ ∆Tdep = −  hw h d  ρ h 0 L f

  ⋅ Tbulk 

for spherical shaped gas hydrates,

Page 12 of 46

(2)

223

where γ hw is the interfacial tension between CO2 hydrate and water, mh is the molar mass of the

224

CO2 hydrate, ρh0 is the density of CO2 hydrate, L f is the latent heat of dissociation of CO2

225

hydrate, and Tbulk is the equilibrium temperature of CO2 hydrate in a bulk condition with no

226

confinement. Herein, due to the pore morphologies of the samples, the morphologies of gas

227

hydrates are assumed to be cylindrical crystals in artificial porous silica gels, and spherical

228

crystals in natural sediment samples.

229

The phase boundaries for various pore sizes, calculated by the Gibbs-Thomson equation, were

230

superimposed on the experimentally obtained phase boundaries of CO2 hydrate for 6SG and

231

15SG in Figure 6b. Herein, the properties of CO2 hydrate were assumed to be equal to those

232

gathered from previous literature,16, 33 that is: γ hw = 0.030 N/m, mh = 174 g/mol, cos θ = 1, ρh0

233

= 1065 kg/m3, and L f = 65.2 kJ/mol. It appears that the PT traces during hydrate dissociation in

234

6SG and 15SG fell within the phase boundaries of the 6–7 nm and 14–15 nm pores, respectively

235

(Figure 6b). This agrees well with the BJH porosimetry results of 6SG and 15SG in Figure 2a.

236

Thus, these results for 6SG and 15SG validate our method for identifying the phase boundary of

237

CO2 hydrate in confined pores, and confirm the use of the Gibbs-Thomson equation and the

238

physical parameters to predict the thermodynamic phase boundary of CO2 hydrate in nano-scale

239

pores.

240

12

ACS Paragon Plus Environment

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

241

3.2

Results of sand and silt: Coarse- versus Fine-grained sediments

242

Figure 7a shows the PT traces of the sand and silt samples, which highlights the dissociation

243

behavior of coarse-grained and fine-grained sediment, respectively. As expected, the PT traces of

244

CO2 hydrate formed in the sand sample follow the bulk phase boundary with no melting

245

temperature depression due to the sufficiently large pore sizes of sand, which are larger than 10

246

µm (Figure 2b). This result is consistent with previous works by Kwon et al.16, 38 and Jang and

247

Santamarina,48 which have shown that the melting temperature depression becomes observable

248

when the pore sizes are less than 1 µm.

249

On the other hand, CO2 hydrate formed in the silt sample was observed to dissociate at a lower

250

temperature than the bulk phase boundary, with a temperature depression (∆Tdep) of ~0.2–0.4°C.

251

The PT traces were placed between the predicted phase boundaries of 200 nm and 600 nm pores

252

(Figure 7b). The morphology of hydrate in silt was assumed to be spherical rather than

253

cylindrical because the majority of pores in fine-grained sediments are open-ended and are

254

mutually connected (Equation 2). This implies that gas hydrate was formed in the pores of ~200–

255

600 nm diameter and melted at the lower dissociation temperature. From the MIP analysis

256

(Figure 2b), the pore size of the silt sample was estimated to range from 100 nm to 20 µm.

257

Accordingly, it appears that CO2 hydrate was preferentially formed in small pores, likely due to

258

the capillarity caused by the water-limiting conditions (or excess gas conditions). Furthermore, it

259

was observed that the pressure began to increase near the phase boundary of the 200 nm pores

260

and gradually converged onto the phase boundary of the 600 nm pores, confirming that the

261

hydrates in smaller pores dissociate earlier as the temperature increases (see 1st heating of silt in

262

Figure 7b).

13

ACS Paragon Plus Environment

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

263

3.3

Results of Ulleung Basin sediment, diatoms, and sand-diatom mixture: Implications

264

to natural fine-grained sediment

265

Figure 8 shows the PT traces of the UB sediment, diatoms, and sand-diatom mixture. CO2

266

hydrate formed in the UB sediment dissociated at a lower temperature than the bulk phase

267

boundary, showing a pronounced dissociation temperature depression (∆Tdep) of ~1.3–1.6°C.

268

Assuming a spherical hydrate morphology, it was found that most of the CO2 hydrate

269

dissociation occurred between the phase boundaries for 45 nm and 55 nm pores. This indicates

270

that gas hydrate was formed in the pores of approximately 45–55 nm for the UB sediment, which

271

exhibited the melting temperature depression when dissociated.

272

It was also found that CO2 hydrate formed in the diatom sample and in the sand-diatom mixture

273

sample underwent melting temperature depression due to the capillarity caused by small pores

274

(Figures 8b and 8d). The temperature depressions (∆Tdep) were ~0.2–0.4°C in the diatoms and

275

0.4–0.5°C in the sand-diatom mixture, implying that the sizes of hydrate-containing pores are

276

~150–600 nm. It was interesting to note that for the sand-diatom mixture, there was a noticeable

277

pressure increase following the bulk phase boundary (as denoted in Figure 8d). This shows that

278

CO2 hydrate was also formed in the pores of sand that were unaffected by capillarity due to their

279

large size, though the majority of the hydrate was formed in the small internal pores offered by

280

diatoms.

281

14

ACS Paragon Plus Environment

Page 14 of 46

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

282

4

Discussion

283

4.1.

Analysis on sizes of hydrate-containing pores

284

The PT traces obtained during hydrate dissociation facilitated the estimation of the diameter of

285

hydrate-containing pores, as depicted in Figures 6-8. The normalized temperature depression

286

(∆Tdep/T) over the course of CO2 hydrate dissociation was plotted against the corresponding pore

287

diameter based on the Gibbs-Thomson equation. Then, the cumulative pore size distribution of

288

each sample, obtained from MIP and/or N2 adsorption methods, was superimposed, as shown in

289

Figure 9. This allowed us to examine the range of pore sizes where CO2 hydrate was formed in

290

each sample.

291

As shown in Figure 9a and 9b, CO2 hydrate formed in the dominant pores of the silica gel

292

samples. While the pore diameter of the 6SG sample mostly ranged from 4–20 nm, the hydrate-

293

containing pores were ~6–7 nm in diameter. These pores correspond to approximately 38–58%

294

of the cumulative pore size distribution, and comprise the middle of the pore size range (Figure

295

9a). It was found that hydrate was formed in these medium-sized pores in the 6SG sample.

296

Likewise, while the pores of 15SG sample ranged from 10–25 nm, most of the hydrates was

297

formed in the pores of ~14–16 nm, which correspond to the 50–75% range in the cumulative

298

pore size distribution. The experimental results of both silica gels indicated that the sediments

299

with a narrow pore size range are expected to preferentially contain gas hydrates in the dominant

300

pores.

301

The silt sample had a relatively wide pore size range (100 nm to 100 µm); however, the hydrate

302

was formed in pores of ~200–600 nm. These pores comprise just 2–7% of the cumulative

303

distribution, as shown in Figure 9c. This corroborates the fact that the excess gas condition (or

304

water-limiting condition) forces water to reside in small pores due to capillarity prior to hydrate

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

305

nucleation. Hence, this condition leads to nucleation of CO2 hydrate in place of water in small

306

pores, despite the fact that confined small pores reduce the water activity and inhibit gas hydrate

307

nucleation. Similarly, the hydrate-containing pores in the UB sediment and diatom samples were

308

also found to be fairly small, with diameters of 45–55 nm and 150–600 nm, respectively (Figures

309

9d and 9e). These pores correspond to the lower range (15–18%) of pore sizes in the UB

310

sediment sample and 3–40% in the diatom sample, confirming the preferential formation of gas

311

hydrates in small pores.

312

In contrast, it appears that the hydrate-containing pores in the sand-diatom mixture sample have

313

a wide size distribution, from 150 nm to tens of microns, although most of the hydrate was

314

formed in small pores less than 1 µm. As can be seen in Figure 2b, the pore size distribution of

315

the sand-diatom mixture sample is presumed to exhibit two dominant pore size regimes; sub-

316

micron pores that are present in diatomaceous grains, and macro-pores formed by sand-diatom

317

skeletal structures. We presume that some hydrate crystals were formed in the macro-pores

318

larger than 1 µm, which is supported by the dissociation of the PT trace following the bulk phase

319

boundary without the Gibbs-Thomson effect, as observed in Figure 8d. It was concluded that in

320

fine-grained sediments, with a wide pore size distribution and under water-limiting condition,

321

gas hydrates preferentially form in small pores over large pores, due to the initial loci of water.

322

In turn, the thermodynamic phase boundaries of such hydrates are heavily affected by the

323

capillarity-induced Gibbs-Thomson effect.

324

16

ACS Paragon Plus Environment

Page 16 of 46

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

325

4.2. Implications for hydrate systems with free gas

326

It has been generally accepted that the Gibbs-Thomson effect prevails in submicron-scale pores

327

and causes the preferential formation of gas hydrates in sandy layers over clay-rich fine-grained

328

sediments.25, 40 This is valid in many hydrate deposits because two-phase equilibrium conditions

329

with no free gas are mostly established above the base of the MHSZ, or above the bottom

330

simulating reflector (BSR). In such hydrate regions with advective, but dissolved, methane

331

fluxes (i.e., Blake Ridge), methane hydrate is expected to form and fill the large pores of sandy

332

sediments, or grow as grain-displacing hydrates in clay-rich sediments.49 However, some field

333

evidence has also suggested the presence of free gas within the regional MHSZ.50-54 A free

334

methane gas phase has often been found in the vicinity of gas vents, where a water-limiting

335

condition is expected, such as Hydrate Ridges offshore Oregon,51, 54, 55 Shenhu area, South China

336

Sea,56 Keathley Canyon Block 151 in northern Gulf of Mexico,57 Krishna‐Godavari Basin

337

offshore India,58 and Mackenzie Delta.59 When the host sediments in these environments are rich

338

in diatoms, clay, silt or ash, hydrate nucleation can take place in small pores due to the capillarity

339

induced by the water-limiting condition. As we observed in our tests, hydrate dissociation is,

340

therefore, presumed to exhibit a considerable Gibbs-Thomson effect with a maximum

341

dissociation temperature depression of ~0.5°C in silt- and diatom-rich sediments, and ~1.5°C in

342

clay-rich sediments.

343

The presence of three-phase equilibrium (vapor, liquid, hydrate) in hydrate deposits has been

344

suggested in some deep sediments near the BSR. The BSR is considered to be an interface

345

between overlying hydrate-bearing sediment and underlying gassy sediment.60 Previous studies

346

have pointed out that the free gas saturation near the BSR varied between 1% and 2% in well 955

347

at Blake Ridge51 and at the IODP Site C0002 in Kumano Basin, offshore Japan.40 In these 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

348

conditions, gas hydrates can be disseminated and the Gibbs-Thomson effect can be pronounced.

349

In fact, the frequent mismatch between the depths of the BSR and BMHSZ can be explained by

350

the presence of clay-dominated sediment and the resultant Gibbs-Thomson effect.25, 40 This is

351

also consistent with our observation of the depression of the dissociation temperature in the clay-

352

rich UB sediment sample.

353

Furthermore, the physical properties of hydrate-bearing fine-grained sediments are undergoing

354

continued investigation, and in many of the lab experiments, hydrate-bearing sediments can be

355

synthesized from water-limiting conditions.42, 50, 61, 62 It is worth noting that the gas hydrates

356

formed in fine-grained sediments can exhibit a considerable dissociation temperature depression

357

due to the Gibbs-Thomson effect, as we observed in this study.

358

18

ACS Paragon Plus Environment

Page 18 of 46

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

359

5

Conclusions

360

This study investigated the Gibbs-Thomson effect on the dissociation temperature depression

361

and the phase boundary of gas hydrates in fine-grained sediments, which have wide pore size

362

distributions of over more than one or two orders of magnitude. In fine-grained sediments with

363

free gas and limited water, gas hydrates preferentially form in smaller pores, less than 1 µm,

364

despite the thermodynamically unfavorable conditions due to the capillarity-induced Gibbs-

365

Thomson effect. Meanwhile, the experimental result of the sand-diatom mixture indicates that

366

even coarse-grained sediments with some diatoms can exhibit the Gibbs-Thomson effect, due to

367

the presence of submicron-scale internal pores of diatoms. It appears that dual porosity, formed

368

by the submicron diatom pores and macro-pores of the sand, is likely to produce dual phase

369

boundaries, one with a dissociation temperature depression and one without. Our experimental

370

results indicate the critical role of the pore size distribution range of fine-grained sediments in

371

governing the dissociation behavior and phase boundary of gas hydrates, in the presence of free

372

gas or in water-limiting conditions.

373 374

Acknowledgments

375

This research was supported by a grant (17CTAP-C129729-01) from Technology Advancement

376

Research Program (TARP) funded by Ministry of Land, Infrastructure and Transport of Korean

377

government, by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) and

378

the Ministry of Trade, Industry and Energy (MOTIE) of the Republic of Korea

379

(20152520100760), and by the MOTIE through the Project “Gas Hydrate Exploration and

380

Production Study” under the management of the Gas Hydrate Research and Development

381

Organization (GHDO) of Korea and the Korea Institute of Geoscience and Mineral Resources

382

(KIGAM). The data presented in this paper can be found in the text, figures or in the table, and 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

383

the data can be requested to the author by email ([email protected]).

388 389

References

390

(1) Sloan, E. D. Gas hydrates: review of physical/chemical properties. Energy Fuels 1998, 12,

391

191-196.

392

(2) Makogon, Y.; Holditch, S.; Makogon, T. Natural gas-hydrates—A potential energy source for

393

the 21st Century. J. Petrol. Sci. Eng. 2007, 56, 14-31.

394

(3) Milkov, A. V. Global estimates of hydrate-bound gas in marine sediments: how much is really

395

out there? Earth-Sci. Rev. 2004, 66, 183-197.

396

(4) Kvenvolden, K. A. Methane hydrate—a major reservoir of carbon in the shallow geosphere?

397

Chem. Geol. 1988, 71, 41-51.

398

(5) Dickens, G. R. Rethinking the global carbon cycle with a large, dynamic and microbially

399

mediated gas hydrate capacitor. Earth Planet. Sci. Lett. 2003, 213, 169-183.

400

(6) Fischer, D.; Mogollón, J. M.; Strasser, M.; Pape, T.; Bohrmann, G.; Fekete, N.; Spiess, V.;

401

Kasten, S. Subduction zone earthquake as potential trigger of submarine hydrocarbon seepage.

402

Nature Geosci. 2013, 6, 647-651.

403

(7) Hornbach, M. J.; Saffer, D. M.; Holbrook, W. S. Critically pressured free-gas reservoirs

404

below gas-hydrate provinces. Nature 2004, 427, 142-144.

405

(8) Phrampus, B. J.; Hornbach, M. J. Recent changes to the Gulf Stream causing widespread gas

406

hydrate destabilization. Nature 2012, 490, 527-530.

407

(9) Rothwell, R.; Thomson, J.; Kahler, G. Low-sea-level emplacement of a very large Late

408

Pleistocene'megaturbidite'in the western Mediterranean Sea. Nature 1998, 392, 377.

409

(10) Inagaki, F.; Kuypers, M. M.; Tsunogai, U.; Ishibashi, J.-i.; Nakamura, K.-i.; Treude, T.;

410

Ohkubo, S.; Nakaseama, M.; Gena, K.; Chiba, H. Microbial community in a sediment-hosted

411

CO2 lake of the southern Okinawa Trough hydrothermal system. Proc. Natl. Acad. Sci. 2006,

412

103, 14164-14169.

413

(11) Lee, S.; Lee, Y.; Lee, J.; Lee, H.; Seo, Y. Experimental verification of methane–carbon

414

dioxide replacement in natural gas hydrates using a differential scanning calorimeter. Environ.

415

Sci. Technol. 2013, 47, 13184-13190.

416

(12) House, K. Z.; Schrag, D. P.; Harvey, C. F.; Lackner, K. S. Permanent carbon dioxide storage 20

ACS Paragon Plus Environment

Page 20 of 46

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

Energy & Fuels

417

in deep-sea sediments. Proc. Natl. Acad. Sci. 2006, 103, 12291-12295.

418

(13) Ohgaki, K.; Takano, K.; Sangawa, H.; Matsubara, T.; Nakano, S. Methane Exploitation by

419

Carbon Dioxide from Gas Hydrates—Phase Equilibria for CO2-CH4 Mixed Hydrate System—.

420

J. Chem. Eng. Jpn. 1996, 29, 478-483.

421

(14) Lee, H.; Seo, Y.; Seo, Y. T.; Moudrakovski, I. L.; Ripmeester, J. A. Recovering methane

422

from solid methane hydrate with carbon dioxide. Angew. Chem. Int. Edit. 2003, 42, 5048-5051.

423

(15) Kimoto, S.; Oka, F.; Fushita, T.; Fujiwaki, M. A chemo-thermo-mechanically coupled

424

numerical simulation of the subsurface ground deformations due to methane hydrate

425

dissociation. Comput. Geotech. 2007, 34, 216-228.

426

(16) Kwon, T. H.; Cho, G. C.; Santamarina, J. C. Gas hydrate dissociation in sediments:

427

Pressure[temperature evolution. Geochem, Geophys, Geosys. 2008, 9, Q08019.

428

(17) Ussler III, W.; Paull, C. Effects of ion exclusion and isotopic fractionation on pore water

429

geochemistry during gas hydrate formation and decomposition. Geo-Mar. Lett. 1995, 15, 37-44.

430

(18) Wang, Y.; Feng, J.-C.; Li, X.-S.; Zhang, Y.; Li, G. Analytic modeling and large-scale

431

experimental study of mass and heat transfer during hydrate dissociation in sediment with

432

different dissociation methods. Energy 2015, 90, 1931-1948.

433

(19) Wang, Y.; Feng, J.-C.; Li, X.-S.; Zhang, Y. Experimental and modeling analyses of scaling

434

criteria for methane hydrate dissociation in sediment by depressurization. Appl. Energy 2016,

435

181, 299-309.

436

(20) Wang, Y.; Feng, J.-C.; Li, X.-S.; Zhang, Y.; Chen, Z.-Y. Fluid flow mechanisms and heat

437

transfer characteristics of gas recovery from gas-saturated and water-saturated hydrate reservoirs.

438

Int. J. Heat Mass Transfer 2018, 118, 1115-1127.

439

(21) Dai, S.; Santamarina, J. C.; Waite, W. F.; Kneafsey, T. J. Hydrate morphology: Physical

440

properties of sands with patchy hydrate saturation. J. Geophys. Res. Solid Earth 2012, 117,

441

B11205.

442

(22) Waite, W. F.; Santamarina, J. C.; Cortes, D. D.; Dugan, B.; Espinoza, D.; Germaine, J.; Jang,

443

J.; Jung, J.; Kneafsey, T. J.; Shin, H. Physical properties of hydrate[bearing sediments. Rev.

444

Geophys. 2009, 47.

445

(23) Boswell, R. Is gas hydrate energy within reach? Science 2009, 325, 957-958.

446

(24) Boswell, R.; Collett, T. S. Current perspectives on gas hydrate resources. Energy. Environ.

447

Sci. 2011, 4, 1206-1215. 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

448

(25) Clennell, M. B.; Hovland, M.; Booth, J. S.; Henry, P.; Winters, W. J. Formation of natural

449

gas hydrates in marine sediments: 1. Conceptual model of gas hydrate growth conditioned by

450

host sediment properties. J. Geophys. Res. Solid Earth 1999, 104, 22985-23003.

451

(26) Park, T.; Kwon, T. Effect of electric field on gas hydrate nucleation kinetics: Evidence for

452

the enhanced kinetics of hydrate nucleation by negatively charged clay surfaces. Environ. Sci.

453

Technol. 2018. 52, 3267-3274.

454

(27) Dickens, G. R.; Quinby[Hunt, M. S. Methane hydrate stability in pore water: a simple

455

theoretical approach for geophysical applications. J. Geophys. Res. Solid Earth 1997, 102, 773-

456

783.

457

(28) Park, T.; Kyung, D.; Lee, W., Effect of organic matter on CO2 hydrate phase equilibrium in

458

phyllosilicate suspensions. Environ. Sci. Technol. 2014, 48, 6597-6603.

459

(29) Handa, Y. P.; Stupin, D. Y. Thermodynamic properties and dissociation characteristics of

460

methane and propane hydrates in 70-. ANG.-radius silica gel pores. J. Phys. Chem. 1992, 96,

461

8599-8603.

462

(30) Uchida, T.; Ebinuma, T.; Ishizaki, T. Dissociation condition measurements of methane

463

hydrate in confined small pores of porous glass. J. Phys. Chem. B. 1999, 103, 3659-3662.

464

(31) Uchida, T.; Ebinuma, T.; Takeya, S.; Nagao, J.; Narita, H. Effects of pore sizes on

465

dissociation temperatures and pressures of methane, carbon dioxide, and propane hydrates in

466

porous media. J. Phys. Chem. B. 2002, 106, 820-826.

467

(32) Uchida, T.; Takeya, S.; Chuvilin, E. M.; Ohmura, R.; Nagao, J.; Yakushev, V. S.; Istomin, V.

468

A.; Minagawa, H.; Ebinuma, T.; Narita, H. Decomposition of methane hydrates in sand,

469

sandstone, clays, and glass beads. J. Geophys. Res. Solid Earth 2004, 109, (B5).

470

(33) Anderson, R.; Llamedo, M.; Tohidi, B.; Burgass, R. W. Experimental measurement of

471

methane and carbon dioxide clathrate hydrate equilibria in mesoporous silica. J. Phys. Chem. B.

472

2003, 107, 3507-3514.

473

(34) Henry, P.; Thomas, M.; Clennell, M. B. Formation of natural gas hydrates in marine

474

sediments: 2. Thermodynamic calculations of stability conditions in porous sediments. J.

475

Geophys. Res. Solid Earth 1999, 104, 23005-23022.

476

(35) Smith, D. H.; Wilder, J. W.; Seshadri, K. Thermodynamics of carbon dioxide hydrate

477

formation in media with broad pore-size distributions. Environ. Sci. Technol. 2002, 36, 5192-

478

5198. 22

ACS Paragon Plus Environment

Page 22 of 46

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

479

(36) Seo, Y.; Lee, H.; Uchida, T. Methane and carbon dioxide hydrate phase behavior in small

480

porous silica gels: three-phase equilibrium determination and thermodynamic modeling.

481

Langmuir 2002, 18, 9164-9170.

482

(37) Li, X.-S.; Zhang, Y.; Li, G.; Chen, Z.-Y.; Yan, K.-F.; Li, Q.-P. Gas hydrate equilibrium

483

dissociation conditions in porous media using two thermodynamic approaches. J. Chem.

484

Thermodyn. 2008, 40, 1464-1474.

485

(38) Kwon, T.-H.; Kim, H.-S.; Cho, G.-C. Dissociation behavior of CO2 hydrate in sediments

486

during isochoric heating. Environ. Sci. Technol. 2008, 42, 8571-8577.

487

(39) Liu, X.; Flemings, P. B. Capillary effects on hydrate stability in marine sediments. J.

488

Geophys. Res. Solid Earth 2011, 116, (B7).

489

(40) Daigle, H.; Dugan, B. Pore size controls on the base of the methane hydrate stability zone in

490

the Kumano Basin, offshore Japan. Geophys. Res. Lett. 2014, 41, 8021-8028.

491

(41) Daigle, H.; Dugan, B. In Data report: Permeability, consolidation, stress state, and pore

492

system characteristics of sediments from Sites C0011, C0012, and C0018 of the Nankai Trough,

493

Proc. Integr. Ocean Drill. Program. 2011, 333 (2).

494

(42) Kim, H. S.; Cho, G. C.; Kwon, T. H. Effect of CO2 hydrate formation on seismic wave

495

velocities of fine[grained sediments. Geochem. Geophy. Geosy. 2013, 14, 1787-1799.

496

(43) Santamarina, J.; Klein, K.; Wang, Y.; Prencke, E. Specific surface: determination and

497

relevance. Can. Geotech. J. 2002, 39, 233-241.

498

(44) Daigle, H.; Johnson, A. Combining mercury intrusion and nuclear magnetic resonance

499

measurements using percolation theory. Transp. Porous. Med. 2016, 111, 669-679.

500

(45) Uchida, T.; Lu, H.; Tomaru, H. Subsurface occurrence of natural gas hydrate in the Nankai

501

Trough area: Implication for gas hydrate concentration. Resour. Geol. 2004, 54, 35-44.

502

(46) Wang, X.; Hutchinson, D. R.; Wu, S.; Yang, S.; Guo, Y. Elevated gas hydrate saturation

503

within silt and silty clay sediments in the Shenhu area, South China Sea. J. Geophys. Res. Solid

504

Earth 2011, 116, (B5).

505

(47) Kwon, T.-H.; Lee, K.-R.; Cho, G.-C.; Lee, J. Y. Geotechnical properties of deep oceanic

506

sediments recovered from the hydrate occurrence regions in the Ulleung Basin, East Sea,

507

offshore Korea. Mar. Petrol. Geol. 2011, 28, 1870-1883.

508

(48) Jang, J.; Santamarina, J. C. Hydrate bearing clayey sediments: Formation and gas

509

production concepts. Mar. Petrol. Geol. 2016, 77, 235-246. 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

510

(49) Daigle, H.; Dugan, B. Origin and evolution of fracture[hosted methane hydrate deposits. J.

511

Geophys. Res. Solid Earth 2010, 115, (B11).

512

(50) Gei, D.; Carcione, J. M. Acoustic properties of sediments saturated with gas hydrate, free

513

gas and water. Geophys. Prospect. 2003, 51, 141-158.

514

(51) Ecker, C.; Dvorkin, J.; Nur, A. M. Estimating the amount of gas hydrate and free gas from

515

marine seismic data. Geophysics 2000, 65, 565-573.

516

(52) Wood, W. T.; Ruppel, C. Seismic and thermal investigations of the Blake Ridge gas hydrate

517

area: A synthesis, Proc. Ocean. Drill. Program Sci. Results 2000, 253-264.

518

(53) Ginsburg, G.; Soloviev, V. Methane migration within the submarine gas-hydrate stability

519

zone under deep-water conditions. Mar. Geol. 1997, 137, 49-57.

520

(54) Trehu, A. M.; Long, P. E.; Torres, M.; Bohrmann, G.; Rack, F.; Collett, T.; Goldberg, D.;

521

Milkov, A.; Riedel, M.; Schultheiss, P. Three-dimensional distribution of gas hydrate beneath

522

southern Hydrate Ridge: constraints from ODP Leg 204. Earth Planet. Sci. Lett. 2004, 222, 845-

523

862.

524

(55) Milkov, A. V.; Dickens, G. R.; Claypool, G. E.; Lee, Y.-J.; Borowski, W. S.; Torres, M. E.;

525

Xu, W.; Tomaru, H.; Tréhu, A. M.; Schultheiss, P. Co-existence of gas hydrate, free gas, and

526

brine within the regional gas hydrate stability zone at Hydrate Ridge (Oregon margin): evidence

527

from prolonged degassing of a pressurized core. Earth Planet. Sci. Lett. 2004, 222, 829-843.

528

(56) Guan, J.; Liang, D.; Wu, N.; Fan, S. The methane hydrate formation and the resource

529

estimate resulting from free gas migration in seeping seafloor hydrate stability zone. J. Asian.

530

Earth. Sci. 2009, 36, 277-288.

531

(57) Cook, A. E.; Goldberg, D.; Kleinberg, R. L. Fracture-controlled gas hydrate systems in the

532

northern Gulf of Mexico. Mar. Petrol. Geol. 2008, 25, 932-941.

533

(58) Lee, M.; Collett, T. Gas hydrate saturations estimated from fractured reservoir at site

534

NGHP[01[10, Krishna[Godavari basin, India. J. Geophys. Res. Solid Earth 2009, 114, (B7).

535

(59) Wright, J.; Nixon, F.; Dallimore, S.; Henninges, J.; Côté, M. Thermal conductivity of

536

sediments within the gas-hydrate-bearing interval at the JAPEX/JNOC/GSC et al. Mallik 5L-38

537

gas hydrate production research well. Bull. Geol. Surv. Can. 2005, 585, 129.

538

(60) Haacke, R. R.; Westbrook, G. K.; Hyndman, R. D. Gas hydrate, fluid flow and free gas:

539

Formation of the bottom-simulating reflector. Earth Planet. Sci. Lett. 2007, 261, 407-420.

540

(61) Kneafsey, T. J.; Tomutsa, L.; Moridis, G. J.; Seol, Y.; Freifeld, B. M.; Taylor, C. E.; Gupta, 24

ACS Paragon Plus Environment

Page 24 of 46

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

541

A. Methane hydrate formation and dissociation in a partially saturated core-scale sand sample. J.

542

Petrol. Sci. Eng. 2007, 56, 108-126.

543

(62) Sun, D.; Englezos, P. Storage of CO2 in a partially water saturated porous medium at gas

544

hydrate formation conditions. Int. J. Greenh. Gas. Con. 2014, 1-8.

545 546

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

547

List of Table

548 549

Table 1. Physical and index properties of sediment samples used in this study

26

ACS Paragon Plus Environment

Page 26 of 46

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

550 551

Energy & Fuels

Table 1. Physical and index properties of sediment samples used in this study Sample 6 nm silica gel (6SG)

Properties

15 nm silica gel (15SG) Sand

Silt

UB sediment

Diatom

Sand-diatom mixture (90%sand+10%diatom) 552 553 554 555 556 557

D50 = 247 µm Sa = 447 m2/g (dry), 1028 m2/g (wet) φ = 0.72 D50 = 555 µm Sa = 310 m2/g (dry), 872 m2/g (wet) φ = 0.70 D50 = 120 µm Sa = 0.02 m2/g φ = 0.46 D50 = 28 µm Sa = 0.85 m2/g (dry), 2.8 m2/g (wet) φ = 0.48 D50 = 5 µm Sa = 32 m2/g (dry), 90 m2/g (wet) φ = 0.62; PL = 46; LL = 62 D50 = 67 µm Sa = 1.8 m2/g (dry), 3.7 m2/g (wet) φ = 0.69 D50 = 123 µm Sa = 19 m2/g (dry) φ = 0.47

Note: D50 is the mean particle diameter; Sa is the specific surface area; φ is the porosity; PL is the plastic limit; and LL is the liquid limit. The specific surface area was measured by the gas adsorption method (or dry method) and/or the methylene blue adsorption method (or wet method; Santamarina et al., 2002).

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

558

List of Figures

559

Figure 1. Grain size distributions of the sediments used in the experiments.

560 561

Figure 2. Pore size distributions of (a) the artificial porous silica gels and (b) the natural sediments.

562 563

Figure 3. SEM images of artificial porous silica gels with 6 nm pores (a), 15 nm pores, sand (c), silt (d), UB sediment (e,f,g), and diatom (h).

564

Figure 4. A schematic drawing of the experimental setup.

565 566

Figure 5. The temperature and pressure traces of Sample SG15 during CO2 hydrate formation and dissociation.

567 568 569

Figure 6. (a) PT traces of CO2 hydrate formed in the 6SG and 15SG samples, and (b) the shifts of their phase boundaries during dissociation. The superimposed dotted lines were predicted by the Gibbs-Thomson equation.

570 571 572

Figure 7. (a) PT traces of CO2 hydrate formed in the sand and silt samples, and (b) the shifts of their phase boundaries during dissociation. The superimposed dotted lines were predicted by the Gibbs-Thomson equation.

573 574 575 576 577 578 579 580 581 582 583 584

Figure 8. (a) PT traces of CO2 hydrate formed in the UB sediment sample, (b) PT traces of CO2 hydrate formed in the diatom and sand-diatom mixture samples, (c) the shifts of the phase boundaries of CO2 hydrate in the UB sediment sample, and (d) the shifts of the phase boundaries of CO2 hydrate in the diatom and sand-diatom mixture samples. The superimposed dotted lines were predicted by the Gibbs-Thomson equation. Figure 9. Sizes of hydrate-containing pores. The normalized dissociation temperature depression (∆Tdep/T) and the cumulative pore volume distribution versus the pore diameter: (a) 6SG, (b) 15SG, (c) silt, (d) UB sediment, (e) diatom, and (f) sand-diatom mixture. The solid lines were calculated from the Gibbs-Thomson model, and the dashed lines were redrawn from the results in Figure 2. The data points indicate the experiment data, and the shaded areas represent the sizes of pores containing CO2 hydrate.

585 586

28

ACS Paragon Plus Environment

Page 28 of 46

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

587 588 589

Figure 1. Grain size distributions of the sediments used in the experiments.

29

ACS Paragon Plus Environment

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

590 591 592

Figure 2. Pore size distributions of (a) the artificial porous silica gels and (b) the natural sediments.

593 594

30

ACS Paragon Plus Environment

Page 30 of 46

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

595

596 597 598

Figure 3. SEM images of artificial porous silica gels with 6 nm pores (a), 15 nm pores, sand (c), silt (d), UB sediment (e,f,g), and diatom (h).

31

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

599

600

Figure 4. A schematic drawing of the experimental setup.

601

32

ACS Paragon Plus Environment

Page 32 of 46

Page 33 of 46 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

602

603 604 605

Figure 5. The temperature and pressure traces of Sample SG15 during CO2 hydrate formation and dissociation.

33

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

606

607 608 609 610

Figure 6. (a) PT traces of CO2 hydrate formed in the 6SG and 15SG samples, and (b) the shifts of their phase boundaries during dissociation. The superimposed dotted lines were predicted by the Gibbs-Thomson equation.

34

ACS Paragon Plus Environment

Page 34 of 46

Page 35 of 46 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

611 612 613 614 615

Figure 7. (a) PT traces of CO2 hydrate formed in the sand and silt samples, and (b) the shifts of their phase boundaries during dissociation. The superimposed dotted lines were predicted by the Gibbs-Thomson equation.

35

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

616

617 618 619 620 621 622

Figure 8. (a) PT traces of CO2 hydrate formed in the UB sediment sample, (b) PT traces of CO2 hydrate formed in the diatom and sand-diatom mixture samples, (c) the shifts of the phase boundaries of CO2 hydrate in the UB sediment sample, and (d) the shifts of the phase boundaries of CO2 hydrate in the diatom and sand-diatom mixture samples. The superimposed dotted lines were predicted by the Gibbs-Thomson equation.

36

ACS Paragon Plus Environment

Page 36 of 46

Page 37 of 46 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

623 624 625 626 627 628 629 630

Energy & Fuels

Figure 9. Sizes of hydrate-containing pores. The normalized dissociation temperature depression (∆Tdep/T) and the cumulative pore volume distribution versus the pore diameter: (a) 6SG, (b) 15SG, (c) silt, (d) UB sediment, (e) diatom, and (f) sand-diatom mixture. The solid lines were calculated from the Gibbs-Thomson model, and the dashed lines were redrawn from the results in Figure 2. The data points indicate the experiment data, and the shaded areas represent the sizes of pores containing CO2 hydrate.

37

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

Figure 1. Grain size distributions of the sediments used in the experiments. 321x246mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 38 of 46

Page 39 of 46 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

Figure 2. Pore size distributions of (a) the artificial porous silica gels and (b) the natural sediments. 643x246mm (300 x 300 DPI)

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

Figure 3. SEM images of artificial porous silica gels with 6 nm pores (a), 15 nm pores, sand (c), silt (d), UB sediment (e,f,g), and diatom (h). 162x237mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 46

Page 41 of 46 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

Figure 4. A schematic drawing of the experimental setup. 225x145mm (300 x 300 DPI)

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

Figure 5. The temperature and pressure traces of Sample SG15 during CO2 hydrate formation and dissociation. 500x208mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 42 of 46

Page 43 of 46 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

Figure 6. (a) PT traces of CO2 hydrate formed in the 6SG and 15SG samples, and (b) the shifts of their phase boundaries during dissociation. The superimposed dotted lines were predicted by the Gibbs-Thomson equation. 643x246mm (300 x 300 DPI)

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

Figure 7. (a) PT traces of CO2 hydrate formed in the sand and silt samples, and (b) the shifts of their phase boundaries during dissociation. The superimposed dotted lines were predicted by the Gibbs-Thomson equation. 643x246mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 44 of 46

Page 45 of 46 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

Figure 8. (a) PT traces of CO2 hydrate formed in the UB sediment sample, (b) PT traces of CO2 hydrate formed in the diatom and sand-diatom mixture samples, (c) the shifts of the phase boundaries of CO2 hydrate in the UB sediment sample, and (d) the shifts of the phase boundaries of CO2 hydrate in the diatom and sand-diatom mixture samples. The superimposed dotted lines were predicted by the Gibbs-Thomson equation. 379x290mm (300 x 300 DPI)

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

Figure 9. Sizes of hydrate-containing pores. The normalized dissociation temperature depression (∆Tdep/T) and the cumulative pore volume distribution versus the pore diameter: (a) 6SG, (b) 15SG, (c) silt, (d) UB sediment, (e) diatom, and (f) sand-diatom mixture. The solid lines were calculated from the GibbsThompson model, and the dashed lines were redrawn from the results in Figure 2. The data points indicate the experiment data, and the shaded areas represent the sizes of pores containing CO2 hydrate. 209x240mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 46 of 46