Water Sorption and Distribution Characteristics in ... - ACS Publications

Oct 10, 2016 - Characteristics of sorption and distribution of water in nanoporous shale are topics of great interest .... Earth-Science Reviews 2018 ...
0 downloads 0 Views 1MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Water Sorption and Distribution Characteristics in Clay and Shale: Effect of Surface Force Jing Li, Xiangfang Li, Keliu Wu, Xiangzeng Wang, Juntai Shi, Liu Yang, Hong Zhang, Zheng Sun, Rui Wang, and Dong Feng Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.6b00927 • Publication Date (Web): 10 Oct 2016 Downloaded from http://pubs.acs.org on October 10, 2016

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

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

Page 1 of 26

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

Water Sorption and Distribution Characteristics in Clay and Shale: Effect of Surface Force

2

Jing Li a, Xiangfang Li a, Keliu Wu a,b*, Xiangzeng Wang c, Juntai Shi a, Liu Yang d,

3

Hong Zhang a, Zheng Sun a, Rui Wang a, Dong Feng a

4

a

5

102249, P.R. China.

6

b

Chemical and Petroleum Engineering, University of Calgary, Calgary, Alberta T2N1N4, Canada

7

c

Shaanxi Yanchang Petroleum (Group) Corp. Ltd., Xi’an 710075, P.R. China.

8

d

9

Beijing 102249, P.R. China

MOE Key Laboratory of Petroleum Engineering, China University of Petroleum (Beijing), Beijing

State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum ,

10

* Corresponding author: Keliu Wu

11

Tel: +1 4039663673; E-mail address: [email protected]

12

Abstract

13

Characteristics of sorption and distribution of water in nanoporous shale are topics of great

14

interest to evaluate unconventional reservoirs. Also, a study of surface force of water/solid

15

interaction at nanoscale is significant for understanding the storage of initial water and the fate of

16

residual treatment liquid in shale systems. In this work, thickness and stability of water film were

17

investigated by vapor sorption experiments on clay and shale samples. Meanwhile, an approach

18

based on surface forces (disjoining pressure), which resulted in the instability of adsorbed film

19

transition into condensed bulk liquid, was developed to describe molecule/pore-wall interactions.

20

Our experimental results directly demonstrated the occurrence of capillary condensation in

21

hydrophilic clay minerals, however, water would not entirely fill in shale nanopores even under high

22

moist conditions. This remarkable finding mainly due to the inaccessibility of water molecules to

23

micropores of hydrophobic organic matter. In addition, the water distribution characteristics are also

24

significantly influenced by pore scale. Under a moist condition with certain relative humidity (e.g.

25

RH=0.98), the water distributed in hydrophilic inorganic pores with different sizes was mainly

26

classified as: (i) capillary water in small pores (e.g. < 6~7 nm), and (ii) water film in large pores

27

(e.g. > 6~7 nm). In contrast, the surface repulsion prevents water condensing and likely results in a

28

monolayer water film sorption in hydrophobic organic pores (e.g. θ=100°). Therefore, in actual shale

29

system with initial moisture content, the inorganic microporosity totally blocked by water might be

30

unable for gas transport or storage, while the hydrophobic organic pores mainly provide effective

31

space for gas accumulation.

32

Keywords: Shale, Clay, Water sorption, Water distribution, Surface force. 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

1 2

Page 2 of 26

1 Introduction Shale gas, as a typical unconventional energy resource, has become an increasingly important [1]

3

part of the world’s gas resource

. The production of recoverable reserves from shale formation

4

mainly depends on gas storage ability and flow capacity

5

within shale gas reservoir under actual condition may be a significant factor that influences the

6

evaluation of gas deliverability. Present studies showed that the methane adsorption capacity of

7

moist shale samples would significantly reduce when compared with dry samples

8

the presence of pre-adsorbed water would obviously decrease the apparent permeability (or diffusion

9

coefficient) and increase the stress sensitivity of shale matrix

[2-5]

. However, the initial water saturation

[6~7]

. Meanwhile,

[8~10]

. Therefore, understanding how

10

water stores and distributes in shale system is essential to reserve estimations and production

11

predictions.

12

Gas shales are heterogeneous rocks with complex mineralogical compositions

[11]

. Meanwhile,

13

a significant difference of wettability exists between organic matters and inorganic minerals, which

14

directly results in a much more complicated water distribution characteristic inside porous shale.

15

Generally, the organic pores (e.g. kerogen pores) forming during hydrocarbon accumulation and

16

generation processes are hydrophobic and almost without water

17

Hu et al. (2014) [13] showed that water molecules could gather as clusters near functional groups in

18

kerogen pores, and experimental results by Prinz et al. (2005) [14] indicated that water could adsorb

19

on hydrophilic sites in organic pores of low-rank coals. Contrast with organic pores, water storage in

20

inorganic pores seems to be much clearer. Especially for the clay minerals, they are usually

21

hydrophilic and present a strong affinity to water. Korb et al. (2014)

22

samples collected form oil shale formation and found that the initial water mainly stored in inorganic

23

pores. Ruppert et al. (2013) [16] investigated the accessibility of pores to methane and water by

24

ultra-small-angle neutron scattering study (USANS) and indicated that inorganic pores less than 30

25

nm display a strong affinity to water. Based on the spontaneous imbibition experiments, liquid water

26

presented a strong capacity to imbibe into shale samples due to the effect of inorganic matter [17~19].

27

Meanwhile, the additional driving forces, such as electrostatic interactions (a part of disjoining

28

pressure) [20], chemical osmosis [17] and adsorption effect [18] might significantly contribute to the

29

strong intake of water or brine by shale during imbibition process. Based on water or vapor sorption

30

experiments, Yua. et al. (2012)

31

nanopores of shale inorganic minerals (such as quartz and chlorite) by field-emission scanning

32

electron microscope (FESEM) images, and clay-rich shale presented more moisture uptake than

[21]

[12]

. However, simulation results by

[15]

analyzed the original

directly showed the water condensation behavior within

ACS Paragon Plus Environment 2

Page 3 of 26

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

1 2

Energy & Fuels

organic-rich shale [22-23]. Generally, the spatial distribution of water in shale gas reservoir has been categorized as two

3

forms, namely (i) hydration water attached to clay surface area, also referred as clay-bound water; (ii)

4

capillary water trapped in pore network system, also referred as capillary-bound water [24~25]. Due to

5

the electric charge in clay surface, water film could strongly adsorb on the external surface of

6

particles or internal surface inside pores though hydrogen bond and electrostatic force. Based on

7

tight rock analysis (TRA) technology, the clay bound water content (CBW) in gas shales could reach

8

about 2.63~7.19% of the total sample volume

9

pores should be paid more attention. However, to the best of our knowledge, the present studies only

10

illustrated a qualitative description of water occurrence in shale system, but the quantitative

11

characteristics of water distribution inside both inorganic and organic microporosity with different

12

surface chemistry (e.g. wettability) are still poorly understood.

[26]

. Thus the water distribution inside these inorganic

13

In this work, firstly, the surface interactions which dominate the stability, thickness and

14

wettability of thin liquid film on substrate were introduced, and a mathematical model for

15

determining sorption characters of water inside different pore sizes was established (in section 2).

16

Subsequently, water sorption and distribution in pure clay and shale samples were investigated by

17

our experimental approach (in section 3), and the comparison between calculated results and

18

experimental data was conducted to verify our proposed model (in section 4). Finally, sorption and

19

distribution characteristics of water within clay and shale porous media have been further discussed

20

as complete wetting case and partial wetting case, respectively. Our present work tries to pave a path

21

for characterizing water distribution in both hydrophilic inorganic and hydrophobic organic porosity,

22

and further demonstrates a better insight on the impact of water on gas storage and transport in shale

23

systems.

24

2 Theory

25

2.1 Surface Forces

26

For thin wetting films, the additional forces caused by interactions between film and solid

27

surface are becoming significant [27]. This action of surface forces is always the primary reason

28

dominating the stability, thickness and wettability of adsorbed film on substrate. Considering the

29

additional forces, the pressure in the thin film could be different from that in bulk liquid, and this

30

difference can be described as disjoining pressure Π(h). Based on Derjaguin et al. (1987)’s theory [28],

31

the total disjoining pressure is characterized by the sum of molecular component Πm, electrical

32

component Πe, and structural component Πs. More details for each term can be found in other ACS Paragon Plus Environment 3

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 26

1

published works [27-29], and the relationship between film thickness h and disjoining pressure Π(h) is

2

given by:

Π ( h ) = Π m ( h) + Π e ( h ) + Π s ( h) h − εε (ζ − ζ ) A Π m ( h ) = H3 ; Π e (h ) = 0 1 2 2 ; Π s ( h) = ke λ h 8π h 2

3

(1)

4

Where, AH is the Hamaker constant, J. ε0 is the permittivity of vacuum, F/m; εr is the relative

5

permittivity of film, dimensionless; ζ1 and ζ2 are the electric potentials of two interfaces between

6

film, respectively, mV; k is the constant of structural force, N/m2; λ is the characteristic thickness of

7

hydration layer, nm.

8

In this work, the values of these parameters refer to other published works. For molecular

9

component, the value of Hamaker constant AH for describing water sorption on inorganic adsorbents

10

[30-32]

11

from 0.1×10-20 to 1×10−20 J. For electrical component, the value of potentials difference △ζ=ζ1-ζ2

12

mainly depends on the electric charge of solid surface and ion inside electrolyte liquid

13

an empirical value of △ζ=0~80 mV is used to described nano-scale water film adsorption on

14

substrate [34-36]. For structural component, the characteristic thickness of hydration layer λ for water

15

approximates to 1~2 nm [34-36], however the value of coefficient k, which is related to the wettability

16

of solid surface, is much more complicated, and it will be discussed later.

17

2.2 Contact angle

18 19

(e.g. sand stone, clay, silt and siliceous shale) and organic adsorbents [33] (e.g. graphite) ranges

[20, 33]

. Here,

Considering the surface force interactions on spreading and wetting dynamics, the contact angle θ of liquid drop on flat surface by Frumkin-Derjaguin approach is given as [36]: h − A H εε 0 (ζ 1 − ζ 2 ) λ + kλe ] cos θ = 1 + (1 / γ l g ) ∫ Π (h)dh = 1 + (1 / γ l g ) ⋅ [ 2 + 2h0 8π h0 h0 2



20

(2)

21

Where, γlg is liquid-gas surface tension, the value of 72.5 mN/m is used in this study; h0 is the

22

equilibrium thickness of wetting film. For the contact angles measured by macroscopic liquid drops,

23

the equilibrium condition can be simplified as Πm(h0)+Πe(h0)+Πs(h0)=0

24

case, this equilibrium thickness approach to infinity in condition of Π(h0)=0, and the contact angle

25

calculated by Eq.(2) is 0°. However, in the partial wetting case, there are several roots of equation

26

Π(h0)=0 and the smallest one is chosen as the equilibrium thickness due to thicker films are

27

metastable or unstable [37].

28 29

[36]

. In the complete wetting

Generally, the behavior of wetting film on substrate surface can be categorized as three states [28-29]

: stable, metastable, and unstable. The stable case, representing continuous positive disjoining ACS Paragon Plus Environment 4

Page 5 of 26

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

pressure, refers as complete wetting of solid surface, while the unstable state, representing

2

continuous negative disjoining pressure, indicates ideal complete non-wetting (however this ideal

3

surface completely non-wetted by water cannot be found in nature

4

characterized by both of positive and negative forces, refers as partial wetting.

[38]

). The metastable state,

5

Based on the research of Churaev et al. (1995a, b and c) [34-36], the changed structure of water

6

molecule near to the solid-liquid interface mainly depends on the wettability of substrate surface. As

7

shown in Fig.1, two different orientations of polar molecules on hydrophilic and hydrophobic

8

surface are illustrated. The dipole moment of polar molecules are oriented preferably vertical to

9

hydrophilic surface which is caused by an attractive structural force (k>0), whereas the parallel

10

orientation to hydrophobic surface is result of a repulsive structural force (k