Investigating the Effect of the Temperature and Pressure on Wettability

Aug 27, 2018 - The objective of the present work is to investigate the influence of the temperature and pressure on the wettability in crude oil–bri...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of South Dakota

Fossil Fuels

Investigating the Effect of Temperature and Pressure on Wettability in Crude oil-Brine-Rock Systems Yongchao Zhang, Jianhui Zeng, Juncheng Qiao, Xiao Feng, and Yuyang Dong Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.8b01404 • Publication Date (Web): 27 Aug 2018 Downloaded from http://pubs.acs.org on August 28, 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 38 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

Investigating the Effect of Temperature and Pressure on Wettability in Crude

2

oil-Brine-Rock Systems

3

Yongchao Zhang†, ‡, ⁋, Jianhui Zeng*, †, ‡, Juncheng Qiao†, ‡, Xiao Feng†, ‡, and Yuyang Dong†, ‡

4

†State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, PR China

5

‡ College of Geosciences, China University of Petroleum, Beijing 102249, PR China

6

⁋ Institute of Petroleum Engineering, Heriot-Watt University, Edinburgh EH14 4AS, UK

7 8 9

ABSTRACT Wettability is a key parameter that affects the petrophysical properties of reservoir

10

formations. The objective of the present work is to investigate the influence of temperature and

11

pressure on the wettability in crude oil-brine-rock (COBR) systems. By using a captive droplet

12

method, the contact angle results of seven minerals and two rock core samples over a range of

13

pressures and temperatures are reported. The data show that raising the pressure from 10 MPa to

14

70 MPa has no discernible effect on the contact angle regardless of mineral type. However, the

15

effect of temperature on the contact angle depends on the primary wettability types of the

16

minerals. For water-wet mineral surfaces, the temperature has a notable impact on the measured

17

contact angles, i.e., the contact angle decreases with increasing temperature; but for neutral-wet

18

or oil-wet samples, the influence of temperature on the measured contact angles is relatively

19

weak. To explain the behaviors of the wettability as a function of temperature, a mathematical

20

calculation is completed based on the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. The

21

calculation results show good consistency with the experimental measurements.

22

Keywords: Physical experiment; temperature; pressure; wettability; contact angle; DLVO theory

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

1

1. INTRODUCTION Wettability is defined as “the tendency of one fluid to adhere to a solid surface in the

2 3

presence of other immiscible fluids” 1, 2. The knowledge of wettability is very important for

4

many aspects in the petroleum industry, including estimation of oil reserves, prediction of

5

production performance, and implementation of enhanced oil recovery (EOR) methods 3, 4.

6

Under oil reservoir conditions, studies of wettability have focused on a crude oil-brine-rock

7

(COBR) system. Changes in the wettability of the rock surface affect the microscopic

8

distribution of formation liquid at the pore scale, which in turn affects the displacement process

9

and relative permeability characteristics 5-7. It has been clearly shown that the wettability

10

characteristics of reservoirs depend on several factors including oil composition, brine

11

composition, pH, temperature and pressure 8-10.

12

1.1 The effect of pressure and temperature on wettability Typically, it is assumed that the temperature may significantly affect the wettability

13 14

performance on rock surfaces, and the pressure dependence of the wettability is relatively weak 2,

15

11, 12

16

studies show different viewpoints. Some researchers 13-15 reported experimentally that the

17

wettability of a rock surface became more water-wet with increasing temperature. In contrast,

18

other researchers 16-20 assumed the wettability would change towards oil-wet rather than water-

19

wet with increasing temperature. In addition, some investigators believed that the changes in the

20

wettability with temperature were not simply oil-wet-towards or water-wet-towards. For example,

21

Wang and Gupta 12 and Rao 21 shared the same viewpoint that sandstones tended to become more

22

oil-wet with increasing temperature, while most carbonates tended to show water-wet behavior.

23

Escrochi et al. 22 experimentally showed that the wettability of Berea sandstone was altered to be

. However, regarding the specific variation in the wettability with temperature, existing

2

ACS Paragon Plus Environment

Page 2 of 38

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

strongly oil-wet within a temperature window of 150 to 400 °C but was subsequently restored its

2

original water-wetness with a further increase in temperature. In recent studies, Lu et al. 11 and

3

Alotaibi et al. 23 reported that the observed changes in the wettability with temperature depended

4

on several factors including the brine composition, salinity, and mineral composition of the

5

system studied.

6

In addition, controversy exists in the literature concerning the changes in the wettability

7

with pressure. Most investigators believe that wettability is not sensitive to pressure conditions,

8

and several studies have claimed a decreasing trend in the contact angles (a common method of

9

quantifying the wettability property) with increasing pressure 12, 18. This effect is also the reason

10 11

why the pressure condition is also considered in our experimental design. In general, the literature review reveals that the mechanisms of wettability changes with

12

temperature and pressure remain in question. Most previous studies have reported only the

13

observed trends without sufficient explanations provided. There have been no previous studies

14

addressing the roles of temperature and pressure on wettability in a systematic manner.

15

1.2 DLVO Theory

16

Several theories have been utilized for interpreting the wettability performance and its

17

controlling factors, which includes the molecular kinetic theory 24, 25, hydrodynamic theory 26,

18

and Derjaguin-Landau-Verwey-Overbeek (DLVO) theory 27-29. The DLVO theory was

19

developed to explain the interactions between colloidal particles and their aggregation behaviors.

20

Some researchers have applied this theory to study the phase interactions and surface energies in

21

a hydrocarbon-brine-rock system 30-34. Hirasaki 35 developed a calculation model based on the

22

DLVO theory and the augmented Laplace-Young equation to determine the stable and 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

1

metastable film-thickness profiles at three-phase contact regions. The interdependence of the

2

spreading, contact angle and the capillary pressure was studied in Hirasaki’s research. Schembre

3

et al. 36 experimentally studied the wettability alteration as a function of temperature. In their

4

work, the DLVO theory was developed for calculations of fines stability and the wettability of

5

silica surfaces. In Chaturvedi et al.’s research 37, the wettability of coal was investigated at scales

6

that ranges from the microscopic to the core. The microscopic wettability was evaluated based

7

on the DLVO theory. The estimates of the contact angles suggest a trend that goes through a

8

maximum at a pH approximately 4. The core-scale wetness obtained from the imbibition

9

experiments shows similar trends as those predicted at the pore scale. Sadeqi-Moqadam et al. 38

10

predicted the contact angle, the work of adhesion, and the stability of a wetting film on rock

11

surfaces using different calculation models based on the DLVO theory and the augmented

12

Laplace-Young equation. In Sadeqi-Moqadam et al.’s work, different factors and their

13

importance of improving the predicted results were thoroughly investigated. However, seen from

14

the published literature, two problems relating to the DLVO calculations of wettability in the

15

COBR system are still unsolved. (1) The obtained experimental results do not always favor the

16

calculation predictions 25. (2) An integrated framework of the combined experiments and

17

theoretical calculations has been rarely reported.

18

The aim of the present work is to investigate the effects of temperature and pressure on the

19

wettability in COBR systems. To achieve this goal, a set of experiments and a DLVO calculation

20

are conducted systematically. The paper is organized as follows. In Section 2, we present our

21

experimental materials and methods. Section 3 is devoted to a presentation of the experimental

22

results. In Section 4, an attempt is made using a mathematical model based on the DLVO theory

23

to explain the observed experimental results. Concluding remarks are presented in Section 5. 4

ACS Paragon Plus Environment

Page 4 of 38

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

2. MATERIALS AND METHODS

2

2.1 Materials

3

The oil sample was taken from the Chang 7 member of the Yanchang formation in the Ordos

4

basin, a typical tight-oil reservoir in China. Before the experiments, a filtration procedure was

5

made to ensure that there was no sand nor other solid deposits remaining in the oil samples. The

6

properties of the oil sample are listed in Table 1. Seven different types of minerals and two rock

7

core samples were used as the rock substrates in the experiments. The two rock core samples are

8

a sandstone core and a carbonate core, whose mineralogy compositions are given in Table 2. The

9

seven single minerals with high purity (more than 98 % based on X-ray diffraction analysis)

10

included calcite, dolomite, illite, chlorite, feldspar, smectite, and quartz. Synthetic brine was

11

prepared as the aqueous phase based on the geochemical analysis of the produced water from

12

Yanchang formation, whose property is shown in Table 3.

13

Table 1

14

Table 2

15

Table 3

16

2.2 Contact Angle Measurements

17

Fig. 1

18

Fig. 2

19

The experimental apparatus used in this work is DSA-100HP from KRÜSS GmbH,

20

Germany (Fig. 1). The DSA100HP combines a drop shape analyzer (DSA100) for contact angle

21

analysis with a high-quality measuring cell for high-pressure and high-temperature applications. 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

1

A schematic diagram of the apparatus is shown in Fig. 2. The experimental apparatus consists of

2

a rock sample holder, a pressure and temperature sealed cell, pressure and temperature control

3

system (injection pumps, oil storage cylinder, water tanks, thermostat, and regulators), and data

4

acquisition and processing system (camera, light source, and computer). In Fig. 2, injection

5

pump 1 is used for oil injection during the test; and injection pump 2 is used for injecting brine

6

and controlling the confining pressure of the cell. The temperature of the measuring cell is

7

controlled by a thermostat connected to the cell. The experimental parameters are orthogonally

8

designed to ensure that the two specified parameters (pressure and temperature) can be estimated

9

independent of each other. The designed pressure measurements increase from 10 MPa to 70

10 11

MPa, and the temperature measurements increase from 35 °C to 110 °C. Specific experimental procedures for the contact angle measurements under different

12

temperature and pressure conditions are given as follows.

13

(1) To avoid unwanted remaining impurities influencing the measurements, the high-pressure

14

and high-temperature cell and oil-flow lines (lines used for oil injection) are cleaned with toluene

15

first, then followed by methanol. The water-flow lines (lines used for brine injection) are flushed

16

with ethanol first, followed by synthetic brine. At least 10 pore volumes of liquid are used for the

17

cleaning or flushing procedures.

18

(2) The seven mineral samples and two rock core samples from hydrocarbons are cleaned using

19

Soxhlet extraction method 10. Afterward, the samples are heated in the oven at 120 °C for 72 h.

20

(3) The samples are cut into substrates with the diameter of Φ 2.5 cm × 0.2 cm. One side of each

21

substrate is polished with aluminum oxide powders with a particle size of 0.25 µm to minimize

6

ACS Paragon Plus Environment

Page 6 of 38

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

the effect of surface roughness on the measurements. Polishing is terminated with brine rinsing

2

and ultrasonic cleaning 39. The samples are heated in the oven at 120 °C for 72 h.

3

(4) The mineral substrates are submerged in the oil sample for aging. The aging time and

4

temperature are 7 days and 60 °C, respectively 40, 41. After aging, the substrates are cleaned with

5

synthetic brine and are submerged into the synthetic brine for 48 h.

6

(6) A substrate is placed on the sample holder. The cell is filled with synthetic brine using

7

injection pump 2. An oil droplet is deposited on the polished surface of the substrate by

8

controlling the injection needle vertically up or down. The volume of the oil droplet is 10 µL.

9

The temperature and pressure condition of the cell is set at a designed value (35 °C and 10 MPa

10

for the first measurement).

11

(7) The pressure is raised from 10 MPa to 70 MPa in steps by injecting brine into the cell. For

12

each pressure measurement, the system is left for at least three hours to ensure that the

13

mineral/brine/oil interfaces reach equilibrium. Three hours later, the images of static contact

14

angle are captured and recorded. The captured images are analyzed to obtain the contact angle

15

values using ADVANCE software by fitting the oil drop profile to the Young-Laplace equation.

16

(8) The substrate is held within the cell, and the measurement temperature is increased

17

incrementally. Step (7) is repeated until all the tests were finished.

18

(9) Steps (6) - (8) are repeated on the other samples.

19

2.3 Limitations of the measurements

20

For a better comparison with other works, the limitations of the experiments are addressed

21

before analyses. First, all the measured contact angles are static contact angles. For the studies

22

under reservoir conditions, the dynamic contact angle (advancing and receding contact angle) 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

1

may be more reliable and representative 42. However, limited by the experimental apparatus,

2

measurements of the dynamic contact angles under high-temperature (up to 110 °C) and high-

3

pressure (up to 70 MPa) conditions are very difficult. The surface roughness and surface

4

contamination, respectively, have a dramatic impact on the measured contact angle values. To

5

reduce the impact of the surface roughness on the results, different experimental polishing

6

procedures have been utilized for substrate preparation in literature. These procedures include

7

the use of carborundum powder 43, abrasive paper 44, alumina powder 39, etc. Drelich et al. 45

8

experimentally studied the effect that different polishing approaches produce on the contact

9

angle measurements. The necessity for the removal of alumina particles from the sample surface

10

after alumina powder polishing was also recognized in Drelich et al.’s research. In our work, the

11

samples were polished and cleaned according to the procedure of Drelich et al. before being used

12

for contact angle measurement. To eliminate the effect of some systematic errors on the results,

13

all the nine samples were polished and cleaned in the same manner. Furthermore, in the field of

14

surface chemistry, surface cleaning is important to reduce the impact of surface contamination on

15

wettability determination. Surface cleaning methods using different chemical reagents and

16

techniques can be found in the literature 46-48. These cleaning procedures guarantee highly

17

accurate and reproducible results for contact angle measurements. However, the main aim of the

18

contact angle measurements in our work is to study the wettability behaviors after oil aging as a

19

function of temperature and pressure using different rock samples. Thus, a simple surface

20

cleaning procedure of synthetic brine rinse and submersion was conducted after oil aging. This

21

simple procedure was probably sufficient to clean the surface properly for our research.

22

3. EXPERIMENTAL RESULTS

8

ACS Paragon Plus Environment

Page 8 of 38

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

Energy & Fuels

In this study, the dependence of the wettability on pressure and temperature is explained

2

according to the contact angle results from experiments. In the petroleum industry, contact

3

angles in a range of 0° to 75° are considered water-wet, angles of 75° to 105° are referred to as

4

neutral-wet and angles of 105° to 180° are considered oil-wet 9. Based on the measured results, it

5

can be calculated that the average contact angles of the samples, as shown in Fig. 3-i to 3-ix, are

6

34.75°, 117.11°, 105.57°, 87.42°, 57.70°, 56.86°, 42.80°, 84.78°, and 33.46°, respectively.

7

According to the petroleum industry standards, sandstone, illite, chlorite, feldspar and quartz can

8

be considered water-wet minerals; carbonate and calcite are oil-wet minerals; the remaining

9

dolomite, smectite are neutral-wet minerals.

10

3.1 The effect of pressure on the surface wettability

11

Fig. 3

12

Fig. 3 illustrates the effect of pressure on the measured contact angles. Figs. 3-i to 3-ix show

13

the contact angle results for sandstone core sample, carbonate core sample, calcite, dolomite,

14

illite, chlorite, feldspar, smectite, and quartz, respectively. From Fig. 3, contact angles of all the

15

samples show minor changes with the increase of pressure. The variations of contact angles

16

within a pressure window of 10 MPa – 70 MPa are 2.0°, 2.0°, 3.0°, 1.0°, 2.5°, 1.5°, 3.0°, 1.5°,

17

and 1° from Fig. 3-i to 3-ix, respectively. Compared with reported experimental results in other

18

works 4, 49, it can be concluded that the surrounding pressure has no discernible effect on the

19

measured contact angles, regardless of mineral type.

20

3.2 The effect of temperature on the surface wettability

21

Fig. 4

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

1

Fig. 4 shows the effect of temperature on the measured contact angles. The results show that

2

temperature has a noticeable impact on the measured contact angles but with different

3

characteristics for different samples. In Figs. 4-i, 4-vii and 4-ix, the contact angles of sandstone,

4

feldspar, and quartz show obvious decreasing trends with the increase of temperature; and the

5

variations of contact angles are 11°, 8°, and 9°, respectively, in a temperature range of 35 °C -

6

110 °C. In Figs. 4-ii, 4-iv, 4-vi, 4-viii, the measured contact angles of the samples (carbonate,

7

dolomite, chlorite, and smectite) remain fairly constant with increasing temperature. Even though

8

an overall decreasing trend with temperature can be found for carbonate, the variation in the

9

contact angle against the carbonate surface is much less than the others. Therefore, the contact

10

angle variation in carbonate is considered constant here. In Fig. 4-iii and Fig. 4-v, the measured

11

contact angles of calcite and illite, respectively, decrease with temperature first. A further

12

increase in temperature induces the contact angles to restore their original values. The contact

13

angle variations of calcite and illite are 5° and 4°, respectively.

14

In brief, the trends of the contact angles as a function of temperature are different for

15

different samples. However, when the wetting preferences are brought into consideration, a

16

general conclusion covering all the variational features can be reached. That is, for water-wet

17

minerals, the temperature has a notable impact on the measured contact angles, i.e., the contact

18

angle decreases with the increasing temperature; however, for neutral-wet or oil-wet samples, the

19

temperature dependence of the contact angles is relatively weak.

20

4. DLVO CALCULATIONS AND DISCUSSION

21 22

To explain the observed trends of the measured contact angles with temperature, a calculation model was built based on DLVO theory.

10

ACS Paragon Plus Environment

Page 10 of 38

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

4.1 Calculation method based on the DLVO theory

2

Fig. 5

3

In the COBR system, the wetting behavior of the rock surface depends on the stability of the

4

thin brine film between the rock surface and oil drop 50 (Fig. 5). This relationship can be

5

described by the augmented Laplace-Young equation 51-53 written as

6

Pc = Π + 2 γ H

,

(1)

is the capillary pressure, Π is disjoining pressure in the thin film, γ is the interfacial

7

where

8

tension, and H is the meniscus curvature. Representing the thin film as a Gibbs diving surface, a

9

relationship between the contact angle and the disjoining pressure can be deduced as 54

Pc

cos θ = 1 +

10

1

γ

Pc

∫ h.d Π ,

(2)

0

11

where h is the separation distance between two surfaces. The disjoining pressure of the brine

12

film 55 is described as the force that tends to disjoin or separate the two interfaces. A negative

13

value of the disjoining pressure means attraction between the contact interfaces. Based on Eq. 2,

14

the correlations between the contact angles and its controlling factors (e.g., temperature) can be

15

explained by calculating the disjoining pressure.

16

The DLVO theory assumes that the disjoining pressure per unit area between two contact

∏( h) ) can be well estimated by three additive contributions, as in Eq. 3 29. The

17

interfaces (

18

contributions are identified as the London van der Waals force ( W V D W ), the electrical double

19

layer force ( W EDL ), and the structural forces ( WS ).

11

ACS Paragon Plus Environment

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

Page 12 of 38

∏( h) = W ( h) +W ( h) +W ( h)

1

VDW

EDL

(3)

S

2

The structural forces are usually assumed to be constant with a change in the temperature and

3

pressure conditions 31. Therefore, the variation in the calculated disjoining pressure with

4

temperature can be explained by only considering the van der Waals (VDW) and electrical

5

double layer forces (EDL) forces.

6

4.1.1 Van der Waals (VDW) forces

7

VDW forces are usually present due to the effect of polarization between different particles.

8

A widely used approximation method to estimate the VDW forces is the Lifshitz continuum

9

method 56, which is given by WVDW = −

10

H . 12π h 2

(4)

11

In Eq. 4, H is the Hamaker constant 35, 57, which defines the strength of the VDW force. Typical

12

values of H range from 10-21 - 10-19 J. In our studied case, the Hamaker constant in the COBR

13

system is calculated through

H =

14

 ε −εw 3 k BT  o 4  εo + εw

 , 

is the Boltzmann constant, T is the surrounding temperature in Kelvin,

15

where

16

dielectric permittivity of crude oil,

17

dielectric permittivity of rock.

18

4.1.2 Electrical double layer (EDL) forces

kB

  εs − εw ⋅   εs + εw

εw

(5)

εo

is the

is the dielectric permittivity of water/brine, and ε s is the

12

ACS Paragon Plus Environment

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

The EDL forces are caused by the development of the charges between an interacting

1 2

surface and the liquid bulk media 31; these forces play an essential role on the wettability

3

determination. The EDL forces can be derived from the solution of the Poisson- Boltzmann

4

equation (PBE). In one dimension, the PBE is given in a form of d 2 Ψ −e = dh 2 ε wε

5

 − zi eψ  K BT

∑ z ρ exp  i

i

 , 

(6)

e is the electron charge, ε is the permittivity of the

6

where Ψ is the electrostatic potential,

7

vacuum, and

8

However, Eq. 6 is a non-linear secondary-order differential equation and cannot be easily

9

analyzed. To obtain a linear solution for the EDL forces and obtain the correlations between

zi

and

ρi

are the valency and charge density of different ions, respectively.

10

EDL forces and temperature, the Debye- Hückel (D-H) approximation is used here to simplify

11

the 1-D PBE. Notably, the D-H approximation is valid for the condition of

12

ψ < 2 5m v

13

59

14

58

z i eψ

K BT < 1

or

. Deviations from the conditional limitations can cause errors in the calculation 38, 54,

. Another essential condition for solving the PBE is the boundary condition. Usually, the

15

solution of the PBE can be obtained using two types of boundary conditions: (i) linear constant

16

potential- constant potential boundary condition (PP condition) and (ii) linear constant charge-

17

constant charge boundary condition (CC condition) 38. The PP condition is represented by those

18

surfaces in which charges can freely transfer through their surfaces. Therefore, their potentials

19

are sustained at a constant level. In contrast, the CC condition assumes a constant charge at their

20

interfaces, where the charge transfer is restricted. For a COBR system, the CC condition seems

21

to be closer to reality for the contact interfaces of the phases 60. However, the general solution for 13

ACS Paragon Plus Environment

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

Page 14 of 38

1

the PBE based on CC condition is very complicated and not convenient for further analysis. For

2

instance, the following equation proposed by Gregory 61 is probably the simplest expression for

3

the CC condition. This equation is reasonably accurate for the condition of a 1:1 electrolyte 27.

  1    ze ζ + ζ kT 2  2  ze (ζ − ζ ) kT  2 e− kh  (   1 2 1 2)   ∏ EDL ( h ) = ρ kT 2 1 +  ekh 2 − e−kh 2   − ze(ζ + ζ ) kT 2 − 2 ,         1 +  kh 12 2− kh 2     e −e 

4

5

6

where

ζ 1 and ζ 2

(7)

are the ζ -potentials of the different media interfaces.

In recent research, Mahani et al. 30 used a DLVO model by using the PP condition to

7

estimate the interaction potential at the crude oil and carbonate rock interfaces. A good

8

consistency between their predicted results and experimental results was obtained in their work.

9

Based on their experience, in our work, the solution of the PBE was calculated using the PP

10

condition assumption, which is believed to be simple and realistic enough for the study of the

11

correlation between wettability and temperature. In a COBR system, an expression for the EDL

12

potential based on the PP condition and D-H approximation is given as



13

EDL

(h) =

εε w k  2ζ sζ o − (ζ s2 + ζ 2 )e − kh  o

e −e kh

− kh

,

are the ζ -potentials of the rock/ brine and crude oil/ brine interfaces,

14

where

15

respectively. The inverse of k in Eq. 8 is known as the Debye length 27.

16

ζs

and

ζo

(8)

k −1 =

εεo KBT 2NAe2 I

,

14

ACS Paragon Plus Environment

(9)

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

∑z

1

In Eq. 9, I = 0.5

2

the ion of the type in brine,

3

Avogadro’s number.

4

4.2 Calculation Results

5

2 i

pi , which represents the ionic strength. The parameter zi is the valence of

pi is its concentration expressed in mol/L, and N A is the

The published ζ -potential data as a function of temperature are very sparse in previous

6

studies. Based on the literature reviewed, the ζ -potential data versus temperature for three

7

minerals (quartz, calcite, and carbonate rock) are carefully selected. Most of these ζ -potential

8

values used for the calculation are from Rodríguez and Araujo’s 62 and Al Mahrouqi et al.’s

9

studies 63. Scattered values are selected from other related literature 11, 19, 30, 36, 59, 64. All these ζ -

10

potential values are measured at a NaCl electrolyte concentration of 0.01 M at a pH of 7.

11

However, due to the limitation of the Zetasizer instrument, all the measuring temperatures of the

12

obtained ζ -potential data are restrictedly below the maximum of 70 °C. Fig. 6 gives the

13

selected ζ -potential values of quartz, calcite, carbonate, and crude oil samples obtained from

14

the literature. It is shown that the ζ -potential of quartz and calcite decrease with an increase in

15

temperature; the ζ -potential of the carbonate and oil sample shows an increasing trend with

16

increasing temperature. However, it should be noticed that the experimental instrument, material,

17

and procedures used in the literature for these selected data were not exactly the same. Several

18

contradictory measured values for the minerals (e.g., kaolinite mineral) in the literature were not

19

used for the DLVO calculation because of their effect on accuracy. Literature data 65, 66 also

20

report varying trend of zeta potential as a function of temperature for other minerals and rock

21

samples.

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

1

Fig. 6

2

Table 4

3

Based on Eqs. 4 – 5 and Eqs. 8 – 9, the disjoining pressure values for different samples are

4

calculated. Other essential parameters used in calculations are given in Table 4. Figs. 7 – 9 show

5

the calculated disjoining pressure isotherms for the quartz-brine-crude oil system, calcite-brine-

6

crude oil system, and carbonate-brine-crude oil system at three different temperatures. Using Eq.

7

9, it can be calculated that the separation distance for the three COBR systems is approximately

8

10 nm. The variation in the separation distance with temperature is negligible for the temperature

9

range studied. From Figs. 7 - 9, the calculated disjoining pressure values are negative for all

10

three cases, indicating the interfacial forces between the minerals and oil droplets are attractive.

11

With an increase in temperature, the disjoining pressure isotherm trends to the X-coordinate axis

12

(the separation distance coordinate axis) in Figs. 7 – 9, which means that the attraction between

13

the mineral and oil phases decreases with increasing temperature. Moreover, the variations in the

14

attraction decrease with the temperature at a separation distance of 10 nm are different for

15

various minerals. The variation for water-wet quartz is much larger than that of the other oil-wet

16

minerals (calcite and carbonate) when the temperature increased from 25 to 70 °C. This behavior

17

indicates that compared with the oil-wet calcite and carbonate samples, the wettability of quartz

18

is more likely to be altered toward water-wetness with an increase in temperature. The calculated

19

results in this section show good consistency with the observed experimental trends in Section 3.

20

Note that the contact angle measurements were conducted under temperature conditions of 35 –

21

110 °C. However, we cannot complete the calculations for temperature conditions above 70 °C

22

due to the absence of the ζ potential data. Therefore, the DLVO trends as a function of

23

temperature are partly extrapolated. 16

ACS Paragon Plus Environment

Page 16 of 38

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

Fig. 7

2

Fig. 8

3

Fig. 9

4

4.3 Discussions

5

Recently, Mahani et al.30 published a research study addressing the influence of brine

6

salinity and temperature on the wettability alteration in carbonates using core sample

7

experiments and DLVO calculations. A comparison between Mahani et al.’s research and our

8

presented work is first presented in the following discussion. As reported in Mahani et al.’s work,

9

increasing temperature will increase the measured contact angles of the carbonate samples.

10

These experimental results seem to contradict our presented results. Actually, the reported

11

contact angle results in Mahani et al.’s work were measured using sea water and 25-times diluted

12

sea water. It can also be found in the authors results that the ζ potential values of the formation

13

water as a function of temperature show obvious different behaviors from the measurements

14

using sea water. Furthermore, the main aim of Mahani et al.’s research was to evaluate the

15

influence of the brine salinity and temperature on the carbonate wettability alteration. Only, three

16

carbonate samples including one dolomite and two limestone specimens are used for the

17

wettability alteration experiment. The contact angle results of the three samples show similar

18

behaviors with temperature. However, in our work, seven mineral samples and two rock core

19

samples with various wettability values are used for studying the contact angle trend as a

20

function of temperature. It is experimentally reported that minerals with different wettability

21

show different alteration behaviors as a function of temperature. The DLVO calculation is used

22

for explaining the observed experimental results. 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

1

The presented work was completed on the basis of a relationship between the ζ -potential

2

and the contact angle. However, there are complexities and limitations in utilizing such a

3

relationship. The reason is given briefly as follows. The contact angle in the COBR system is

4

controlled by the thickness of the thin brine film (Fig. 5), which is affected by the charges at the

5

rock-brine and crude oil-brine interfaces of the thin film 67. The charges at the interfaces of the

6

thin brine film, however, are difficult to measure experimentally. Instead, the ζ potential data

7

was utilized for representing the behavior of the thin brine film. But, the problem is that the ζ

8

potential values do not necessarily represent the charges at the interfaces of the thin brine film 11.

9

This is the first limitation of this work. Another limitation of this research is that the ζ potential

10

measurements are not conducted using the same samples as the contact angle measurements.

11

Thus, the DLVO calculations can be explanations rather than validations for the observed contact

12

angle trends. In terms of accuracy, only three minerals including quartz, calcite, and carbonate

13

with ζ potential data obtained from the literature are used for the DLVO calculations. The

14

selected three samples seem representative since both water-wet mineral and oil-wet samples are

15

taken into consideration. In our research, a next step would include using experimentally

16

measured ζ potential data to estimate the contact angle behaviors of various rock samples with

17

temperature and then compare these results with the observed experimental results. Some

18

important aspects including surface roughness and oil aging time will be considered in the future

19

work.

20

5. CONCLUSION

21 22

Based on the experimental and DLVO calculation results, the following conclusions can be drawn: 18

ACS Paragon Plus Environment

Page 18 of 38

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

Energy & Fuels

1. Raising the pressure from 10 MPa to 70 MPa has no discernible influence on the mineral surface wettability regardless of mineral type in the COBR system.

2 3

2. In the COBR system, the variation in the mineral surface wettability as a function of

4

temperature depends on the wettability types. For water-wet mineral surfaces, the

5

temperature has a noticeable impact on the contact angles, i.e., the contact angle

6

decreases with increasing temperature; however, for neutral-wet or oil-wet samples, the

7

temperature dependence of the measured contact angle is relatively weak.

8

3. A DLVO calculation considering the van der Waals, electrostatic forces is successfully

9

utilized for explaining the temperature dependence of wettability in COBR system.

10

4. The DLVO calculations show that the disjoining pressure of water-wet mineral and oil-

11

wet minerals increases with increasing temperature, but the variation in the disjoining

12

pressure with temperature for water-wet minerals is much larger than that for oil-wet

13

minerals. These results agree with the observed experimental results.

14 15

AUTHOR INFORMATION

16

Corresponding Author

17

*E-mail: [email protected]

18

ORCID

19

Yongchao Zhang: 0000-0002-0288-9398

20

Jianhui Zeng: 0000-0003-1000-2205

21

Notes

22

The authors declare no competing financial interest. 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

1

ACKNOWLEDGMENTS

2

This work was financially supported by the National Natural Science Foundation of China

3

(Grant No. 41330319). I would like to express my sincere thanks to Heriot-Watt University for

4

their help in the preparation of the manuscript.

5 6

20

ACS Paragon Plus Environment

Page 20 of 38

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

Reference

2

1.

3

Fund of AIME: New York: 1971; Vol. 3.

4

2.

5

pressure, and crude oil composition on interfacial properties. SPE Reservoir Engineering 1986, 1

6

(04), 321-328.

7

3.

8

Shale Formations, SPE Western North American and Rocky Mountain Joint Meeting, Society of

9

Petroleum Engineers: 2014.

Craig, F. F., The reservoir engineering aspects of waterflooding. HL Doherty Memorial

Hjelmeland, O.; Larrondo, L., Experimental investigation of the effects of temperature,

Morsy, S.; Sheng, J. In Surfactant Preflood to Improve Waterflooding Performance in

10

4.

Najafi-Marghmaleki, A.; Barati-Harooni, A.; Soleymanzadeh, A.; Samadi, S. J.;

11

Roshani, B.; Yari, A., Experimental investigation of effect of temperature and pressure on

12

contact angle of four Iranian carbonate oil reservoirs. Journal of Petroleum Science and

13

Engineering 2016, 142, 77-84.

14

5.

15

Petroleum Science and Engineering 1998, 20 (3-4), 155-160.

16

6.

17

SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers: 2013.

18

7.

19

temperature, oil asphaltene content, and water salinity on wettability alteration, SPE/DGS Saudi

20

Arabia Section Technical Symposium and Exhibition, 2011; Society of Petroleum Engineers:

21

2011.

22

8.

23

oil/brine/rock interactions and oil recovery. Journal of Petroleum Science and Engineering 1999,

24

24 (2-4), 99-111.

Buckley, J.; Liu, Y., Some mechanisms of crude oil/brine/solid interactions. Journal of

Kathel, P.; Mohanty, K. K. In EOR in tight oil reservoirs through wettability alteration,

Al-Aulaqi, T.; Grattoni, C.; Fisher, Q.; Musina, Z.; Al-Hinai, S. In Effect of

Tang, G.; Morrow, N. R., Influence of brine composition and fines migration on crude

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

1

9.

Anderson, W., Wettability literature survey-part 2: Wettability measurement. Journal of

2

petroleum technology 1986, 38 (11), 1,246-1,262.

3

10.

4

on conventional and unconventional reservoir cores, Unconventional Resources Technology

5

Conference (URTEC), 2015.

6

11.

7

oil/brine/rock systems. Energy & Fuels 2017, 31 (5), 4989-4995.

8

12.

9

wettability using modified pendant drop method, SPE Annual Technical Conference and

Teklu, T. W.; Alameri, W.; Kazemi, H.; Graves, R. M. In Contact angle measurements

Lu, Y.; Najafabadi, N. F.; Firoozabadi, A., Effect of temperature on wettability of

Wang, W.; Gupta, A. In Investigation of the effect of temperature and pressure on

10

Exhibition, Society of Petroleum Engineers: 1995.

11

13.

12

irreducible water saturation and relative permeability of unconsolidated sands. Society of

13

Petroleum Engineers Journal 1970, 10 (02), 171-180.

14

14.

15

Penneabilities. Journal of Canadian Petroleum Technology 1974, 13 (04).

16

15.

17

waterflooding. SPE Reservoir Engineering 1997, 12 (04), 269-276.

18

16.

19

brine on oil recovery in sandstone reservoirs. Journal of Petroleum Science and Engineering

20

2014, 122, 384-395.

21

17.

22

evaluation of the effect of temperature on wettability of the grosmont carbonate reservoir in

Poston, S.; Ysrael, S.; Hossain, A.; Montgomery III, E., The effect of temperature on

McCaffery, F.; Bennion, D., The Effect OfWettability On Two-Phase Relative

Tang, G.; Morrow, N. R., Salinity, temperature, oil composition, and oil recovery by

Nasralla, R. A.; Nasr-El-Din, H. A., Impact of cation type and concentration in injected

Modaresghazani, J.; Moore, R. G.; Mehta, S. A.; Ursenbach, M. G., Experimental

22

ACS Paragon Plus Environment

Page 22 of 38

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

Alberta. In SPE Canada Heavy Oil Technical Conference, Society of Petroleum Engineers:

2

Calgary, Alberta, Canada, 2016.

3

18.

4

Interfacial Tension of Quartz-Water-Bitumen Systems, Canadian International Petroleum

5

Conference, Petroleum Society of Canada: 2009.

6

19.

7

alteration of carbonate reservoirs, SPE/DOE Symposium on Improved Oil Recovery, Society of

8

Petroleum Engineers: 2006.

9

20.

Rajayi, M.; Kantzas, A. In Effect of Temperature and Pressure on Contact Angle and

Hamouda; Anis, A.; Gomari, R.; Ali, K. In Influence of temperature on wettability

Hamouda, A. A.; Karoussi, O., Effect of temperature, wettability and relative

10

permeability on oil recovery from oil-wet chalk. Energies 2008, 1 (1), 19-34.

11

21.

12

Recovery Symposium, Society of Petroleum Engineers: 1996.

13

22.

14

at elevated temperatures: the consequences of asphaltene precipitation, SPE International

15

Symposium and Exhibition on Formation Damage Control, Society of Petroleum Engineers:

16

2008.

17

23.

18

salinity water in sandstone reservoirs. SPE Reservoir Evaluation & Engineering 2011, 14 (06),

19

713-725.

20

24.

21

Physicochemical and Engineering Aspects 1994, 93, 15-23.

Rao, D. N. In Wettability effects in thermal recovery operations, SPE/DOE Improved Oil

Escrochi, M.; Nabipour, M.; Ayatollahi, S. S.; Mehranbod, N. In Wettability alteration

Alotaibi, M. B.; Nasralla, R. A.; Nasr-El-Din, H. A., Wettability studies using low-

Hayes, R. A.; Ralston, J., The dynamics of wetting processes. Colloids and Surfaces A:

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

1

25.

De Ruijter, M.; Kölsch, P.; Voué, M.; De Coninck, J.; Rabe, J., Effect of temperature

2

on the dynamic contact angle. Colloids and Surfaces A: Physicochemical and Engineering

3

Aspects 1998, 144 (1-3), 235-243.

4

26.

5

molecular-kinetic? Journal of colloid and interface science 2005, 288 (1), 213-221.

6

27.

Israelachvili, J. N., Intermolecular and surface forces. Academic press: 2011.

7

28.

Trefalt, G.; Borkovec, M., Overview of DLVO theory. Laboratory of Colloid and Surface

8

Chemistry, University of Genebra 2014.

9

29.

Ranabothu, S. R.; Karnezis, C.; Dai, L. L., Dynamic wetting: hydrodynamic or

Derjaguin, B., Theory of the stability of strongly charged lyophobic sols and the adhesion

10

of strongly charged particles in solutions of electrolytes. Acta Physicochim. USSR 1941, 14, 633-

11

662.

12

30.

13

V., Insights into the impact of temperature on the wettability alteration by low salinity in

14

carbonate rocks. Energy & Fuels 2017, 31 (8), 7839-7853.

15

31.

16

crude oil-brine-rock systems, SPE/DOE Symposium on Improved Oil Recovery, Society of

17

Petroleum Engineers: 2004.

18

32.

19

pressure: Implications on geological CO2 storage. International Journal of Greenhouse Gas

20

Control 2017, 63, 26-36.

21

33.

22

surface force in low salinity water flooding. Journal of Molecular Liquids 2016, 221, 658-665.

Mahani, H.; Menezes, R.; Berg, S.; Fadili, A.; Nasralla, R.; Voskov, D.; Joekar-Niasar,

Busireddy, C.; Rao, D. N. In Application of DLVO theory to characterize spreading in

Zhou, Y.; Hatzignatiou, D. G.; Helland, J. O., On the estimation of CO2 capillary entry

Xie, Q.; Saeedi, A.; Pooryousefy, E.; Liu, Y., Extended DLVO-based estimates of

24

ACS Paragon Plus Environment

Page 24 of 38

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

34.

Tokunaga, T. K., DLVO-based estimates of adsorbed water film thicknesses in geologic

2

CO2 reservoirs. Langmuir 2012, 28 (21), 8001-8009.

3

35.

4

1991, 6 (02), 217-226.

5

36.

6

imbibition at elevated temperatures. Journal of Petroleum Science and Engineering 2006, 52 (1-

7

4), 131-148.

8

37.

9

characteristics of Powder River Basin coal. International Journal of Coal Geology 2009, 77 (1-

Hirasaki, G., Wettability: fundamentals and surface forces. SPE Formation Evaluation

Schembre, J.; Tang, G.; Kovscek, A., Wettability alteration and oil recovery by water

Chaturvedi, T.; Schembre, J.; Kovscek, A., Spontaneous imbibition and wettability

10

2), 34-42.

11

38.

12

behavior of oil/water/reservoir rock interfaces: The implication for improvement in wettability

13

prediction. Colloids and Surfaces A: Physicochemical and Engineering Aspects 2016, 490, 268-

14

282.

15

39.

16

determination by contact angle measurements: hvbB coal–water system with injection of

17

synthetic flue gas and CO2. Journal of colloid and interface science 2011, 364 (1), 237-247.

18

40.

19

International Symposium on Oilfield Chemistry, Society of Petroleum Engineers: 1991.

20

41.

21

fluid distribution in Middle Bakken Cores. Journal of Petroleum Science and Engineering 2018,

22

166, 1042-4052.

Sadeqi-Moqadam, M.; Riahi, S.; Bahramian, A., An investigation into the electrical

Kaveh, N. S.; Rudolph, E. S. J.; Wolf, K.-H. A.; Ashrafizadeh, S. N., Wettability

Jia, D.; Buckley, J.; Morrow, N. In Control of core wettability with crude oil, SPE

Karimi, S.; Saidian, M.; Kazemi, H., Experimental study of the effect of core aging on

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

1

42.

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

2

and temperature on dynamic contact angles and interfacial tension of supercritical

3

CO2/brine/quartz systems. International Journal of Greenhouse Gas Control 2014, 28, 147-155.

4

43.

5

Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases 1982, 78

6

(1), 61-73.

7

44.

8

coal surface and its effect on coal wettability. Journal of Colloid and Interface Science 2002, 247

9

(1), 107-116.

Lamb, R. N.; Furlong, D. N., Controlled wettability of quartz surfaces. Journal of the

Gosiewska, A.; Drelich, J.; Laskowski, J.; Pawlik, M., Mineral matter distribution on

10

45.

Drelich, J.; Laskowski, J.; Pawlik, M.; Veeramasuneni, S., Preparation of a coal surface

11

for contact angle measurements. Journal of adhesion science and technology 1997, 11 (11),

12

1399-1431.

13

46.

14

surfaces: Impact on water–CO2–quartz and glass contact angle measurements. International

15

Journal of Greenhouse Gas Control 2014, 22, 325-328.

16

47.

17

Oostrom, M., Correlation of oil–water and air–water contact angles of diverse silanized surfaces

18

and relationship to fluid interfacial tensions. Langmuir 2012, 28 (18), 7182-7188.

19

48.

20

CO2 at storage conditions. International Journal of Greenhouse Gas Control 2013, 12, 18-25.

21

49.

22

tension and contact angle of Khuff gas reservoir, Saudi Arabia, SPE/DGS Saudi Arabia Section

Iglauer, S.; Salamah, A.; Sarmadivaleh, M.; Liu, K.; Phan, C., Contamination of silica

Grate, J. W.; Dehoff, K. J.; Warner, M. G.; Pittman, J. W.; Wietsma, T. W.; Zhang, C.;

Farokhpoor, R.; Bjørkvik, B. J.; Lindeberg, E.; Torsæter, O., Wettability behaviour of

Okasha, T. M.; Al-Shiwaish, A.-J. A. In Effect of temperature and pressure on interfacial

26

ACS Paragon Plus Environment

Page 26 of 38

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

Technical Symposium and Exhibition, Al-Khobar, Saudi Arabia Society of Petroleum Engineers:

2

Al-Khobar, Saudi Arabia 2010.

3

50.

4

and silica. Journal of Colloid and Interface Science 1978, 65 (2), 296-306.

5

51.

6

imbibition based on Hagen–Poiseuille flow in tortuous capillaries with variably shaped apertures.

7

Langmuir 2014, 30 (18), 5142-5151.

8

52.

9

Journal of colloid and interface science 2018, 521, 226-231.

Aronson, M. P.; Petko, M.; Princen, H., On the stability of aqueous films between oil

Cai, J.; Perfect, E.; Cheng, C.-L.; Hu, X., Generalized modeling of spontaneous

Xiao, J.; Cai, J.; Xu, J., Saturated imbibition under the influence of gravity and geometry.

10

53.

Cai, J.; Yu, B.; Zou, M.; Luo, L., Fractal characterization of spontaneous co-current

11

imbibition in porous media. Energy & Fuels 2010, 24 (3), 1860-1867.

12

54.

13

solid surfaces. Journal of colloid and interface science 1996, 181 (2), 443-455.

14

55.

15

pressure of a thin water film being pressed by an oil droplet. Journal of colloid and interface

16

science 2010, 348 (1), 255-264.

17

56.

18

with applications to wetting phenomena. Advances in Colloid and Interface Science 1980, 14 (1),

19

3-41.

20

57.

21

1937, 4 (10), 1058-1072.

22

58.

23

electrical double layers. Journal of Colloid and Interface Science 1973, 44 (2), 347-355.

Basu, S.; Sharma, M. M., Measurement of critical disjoining pressure for dewetting of

Kumar, D.; Biswas, S. K., Contribution of different physical forces to the disjoining

Hough, D. B.; White, L. R., The calculation of Hamaker constants from Liftshitz theory

Hamaker, H., The London-van der Waals attraction between spherical particles. physica

Kar, G.; Chander, S.; Mika, T., The potential energy of interaction between dissimilar

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

1

59.

Buckley, J.; Takamura, K.; Morrow, N., Influence of electrical surface charges on the

2

wetting properties of crude oils. SPE Reservoir Engineering 1989, 4 (03), 332-340.

3

60.

4

Charging and aggregation properties of carboxyl latex particles: Experiments versus DLVO

5

theory. Langmuir 2000, 16 (6), 2566-2575.

6

61.

7

and Interface Science 1975, 51 (1), 44-51.

8

62.

9

reservoir minerals. Journal of colloid and interface science 2006, 300 (2), 788-794.

Behrens, S. H.; Christl, D. I.; Emmerzael, R.; Schurtenberger, P.; Borkovec, M.,

Gregory, J., Interaction of unequal double layers at constant charge. Journal of Colloid

Rodríguez, K.; Araujo, M., Temperature and pressure effects on zeta potential values of

10

63.

Al Mahrouqi, D.; Vinogradov, J.; Jackson, M. D., Temperature dependence of the zeta

11

potential in intact natural carbonates. Geophysical Research Letters 2016, 43 (22).

12

64.

13

of silicates. Colloids and Surfaces 1986, 21, 355-369.

14

65.

15

investigation of the effect of temperature, salinity and salt type on brine/mineral interfacial

16

properties. 2017, 59, 136-147.

17

66.

18

interfacial properties of low to high rank coal seams. Fuel 2017, 194, 211-221.

19

67.

20

salinity brine. Current Opinion in Colloid & Interface Science 2015, 20 (2), 105-114.

Ramachandran, R.; Somasundaran, P., Effect of temperature on the interfacial properties

Arif, M.; Jones, F.; Barifcani, A.; Iglauer, S. J. I. J. o. G. G. C., Electrochemical

Arif, M.; Jones, F.; Barifcani, A.; Iglauer, S., Influence of surface chemistry on

Myint, P. C.; Firoozabadi, A., Thin liquid films in improved oil recovery from low-

21 22

28

ACS Paragon Plus Environment

Page 28 of 38

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

Energy & Fuels

Table 1. Properties and compositions of the oil sample used in the experiments Properties/ compositions

Values

saturates content (wt %)

60.01

aromatics content (wt %)

22.97

resins content (wt %)

9.01

asphaltenes content (wt %)

2.84

density (g/ml)

0.95

viscosity at 25 °C (cP)

747.60

viscosity at 50 °C (cP)

79.54

viscosity at 70 °C (cP)

46.50

2 3

Table 2. Mineral compositions (in wt %) of the sandstone and carbonate core samples Mineral compositions

Sandstone

Carbonate

Quartz

48

trace

Feldspar

22

-

Calcite

6

80

Dolomite

3

18

Chlorite

5

-

Illite/smectite

10

-

Kaolinite

6

-

4 5 6

Table 3. Brine properties and compositions

(g/cm3) Salinity (mg/L) Density(

33634

1.051

Ionic concentration( (mg/L) )

pH

6.729

Na+

Ca2+

Mg2+

Cl-

10419

1750

318

21147

7 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

Page 30 of 38

1 2

Table 4. Parameters used for the disjoining pressure calculation Parameter

Unit

Value

KB

J/K

1.38×10-23

εw

F/m

ε

F/m

8.85×10-12

εo

F/m

2.00

εs

F/m

7.00

e

Coulomb

1.60×1019

NA

mol-1

6.02×1023

78.30 (298 K)

69.90 (323 K)

3 4 5

30

ACS Paragon Plus Environment

63.80 (343 K)

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

Fig. 1 Experimental apparatus DSA100HP from KRÜSS GmbH

3

4 5

Fig. 2 A schematic diagram of the wettability measurement apparatus

6 7 8 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

Page 32 of 38

i: sandstone (av.CA=34.75°)

ii: carbonate (av.CA=117.11°)

iii: calcite (av.CA=105.57°)

iv: dolomite (av.CA=87.42°)

v: illite (av.CA=57.70°)

vi: chlorite (av.CA=56.86°)

vii: feldspar (av.CA=42.80°)

viii: smectite (av.CA=84.78°)

ix: quartz (av.CA=33.46°)

Fig. 3 The effect of pressure on the contact angles. The samples in Fig. 3-i to 3-ix represent sandstone, carbonate, calcite, dolomite, illite, chlorite, feldspar, smectite, and quartz, respectively. The term “av.CA” refers to the average contact angle the sample 1

32

ACS Paragon Plus Environment

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

i: sandstone (av.CA=34.75°)

ii: carbonate(av.CA=117.11°)

iii: calcite(av.CA=105.57°)

iv: dolomite(av.CA=87.42°)

v: illite (av.CA=57.70°)

vi: chlorite (av.CA=56.86°)

vii: feldspar (av.CA=42.80°)

viii: smectite (av.CA=84.78°)

ix: quartz (av.CA=33.46°)

Fig. 4 The effect of temperature on the contact angles. The samples in Fig. 4-i to 4-ix, represent sandstone, carbonate, calcite, dolomite, illite, chlorite, feldspar, smectite, and quartz, respectively. The term “av.CA” refers to the average contact angle the sample 1

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

Fig. 5 Profile of the contact region in the COBR system 1 2 3 4 5 6 7 8

34

ACS Paragon Plus Environment

Page 34 of 38

Page 35 of 38

Temperature (°C)

Temperature (°C) 0

10

20

30

40

50

60

70

80

90

0

100

0

10

20

30

40

50

60

70

80

90

100

40 20

-20

0

ζ potiential (mV)

ζ potiential (mV)

-40 -60 -80 -100

-20 -40 -60 -80 -100 -120

-120

Calcite Estimated Calcite Oil Sample

-140

Quartz Estimated Quartz Oil Sample

-140 -160

-160

Quartz

Calcite

Temperature (°C) 0

10

20

30

40

50

60

70

80

90

100

0 -20 -40

ζ potiential (mV)

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

Energy & Fuels

-60 -80 -100 -120

Carbonate Estimated Carbonate Oil Sample

-140 -160

Carbonate Fig. 6 Estimated ζ potential values used in the calculations obtained from the literature. The asterisk symbols in the figures are the estimated data for our calculation; and the filled symbols depict the experimental data from the literature 1

35

ACS Paragon Plus Environment

Energy & Fuels

Separation Distance (h) (nm) 0

10

20

30

40

50

60

70

80

90

100

0.0E+00

The Disjoining Pressure (atm)

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

-2.0E-07

-4.0E-07

The attraction decreases with the increasing temperature

-6.0E-07

Quartz at 25 °C Quartz at 50 °C Quartz at 70 °C

-8.0E-07

-1.0E-06

-1.2E-06

Fig. 7 Calculated disjoining pressure results versus separation distance in quartz at different temperatures

36

ACS Paragon Plus Environment

Page 37 of 38

Separation Distance (h) (nm) 0

10

20

30

40

50

60

70

80

90

100

0.0E+00

The Disjoining Pressure (atm)

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

-2.0E-07

-4.0E-07

The attraction decreases with the increasing temperature

-6.0E-07

Calcite at 25 °C Calcite at 50 °C Calcite at 70 °C

-8.0E-07

-1.0E-06

-1.2E-06

Fig. 8 Calculated disjoining pressure results versus separation distance in calcite at different temperatures 1

37

ACS Paragon Plus Environment

Energy & Fuels

Separation Distance (h) (nm) 0

10

20

30

40

50

60

70

80

90

100

0.0E+00

The Disjoining Pressure (atm)

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

Page 38 of 38

-2.0E-07

-4.0E-07

The attraction decreases with the increasing temperature

-6.0E-07

Carbonate at 25 °C Carbonate at 50 °C Carbonate at 70 °C

-8.0E-07

-1.0E-06

-1.2E-06

Fig. 9 Calculated disjoining pressure results versus separation distance in carbonate at different temperatures 1

38

ACS Paragon Plus Environment