1. Introduction

Geomaterials

2161-7538

Scientific Research Publishing

10.4236/gm.2014.43008

GM-47299

Articles

EARTH & ENVIRONMENTAL SCIENCES

Effect of Fluid Flow Characteristics on Migration of Nano-Particles in Porous Media

Saeed

Sourani

¹^*Mohammad

Afkhami

²^*Yousef

Kazemzadeh

³^*Hooman

Fallah

⁴^*

Department of Petroleum Engineering, Shiraz University, Fars, Iran

Department of Petroleum Engineering, Islamic Azad University of Lamerd, Fars, Iran

Department of Petroleum Engineering, Islamic Azad University of Marvdasht, Fars, Iran

Department of Petroleum Engineering, Islamic Azad University, Firoozabad Branch, Firoozabad, Fars, Iran

* E-mail:Soorani.saeed@yahoo.com(SS);Mohammad.afkhami1@gmail.com(MA);yusefkazemzade@yahoo.com(YK);Hoomanfallah2@gmail.com(HF);

25062014

040373812 May 201427 May 2014 20 June 2014

2014

This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

Inter-particle and molecular forces are of most important forces affecting fines migration. Investigating the impact of these types of forces requires information about electrical properties of constitutive component and Nano-sized solids. Surface potential is a well-known parameter, which can be used to measure aforementioned forces. The force balance among electrical, gravity, drag and buoyance forces tracing on a particle can be estimated using fluid properties and physical properties of reservoir such as surface potential, pore size distribution and fine size distribution. The task is to set these forces at conditions that fine migration does not take place. This paper investigates the impact of several parameters that could influence pertinent forces. Effect of ionic strength of flowing fluid is taken into account to evaluate Debye length that determines double layer forces. The impact of injection rate on a parameter named erosion number is also studied. The results of this study show that how introducing salts and injection rate can affect stability of fines on their locations.

Fines Migration Reservoir Damage Inter-Particle Forces Electrical Double Layer Injection Rate

1. Introduction

Water injection in hydrocarbon reservoirs for EOR purposes generally leads to formation damage. Several mechanisms of this phenomenon have been explained by different authors [1] -[3] . Formation damage could be defined as formation permeability reduction. Formation damage reduces production and injectivity of formation. Figure 1 [4] shows a case in which injectivity of an Iranian oil field reduced as a result of fine migration.

Several parameters play role in the process of formation damage. These parameters can be classified in two categories: properties of porous media and the characteristics of flowing fluids. Pore size distribution, pore continuity, concentration and size distribution of fines are the parameters of the former category. Injection rate, salinity and pH of water (as a displacing phase) are the parameters of the latter one. Investigating the parameters of the first group is out of scope of this study. Therefore, ionic strength (concentration and ionic valence) and injection rate (as a function of velocity of fluid) are studied as two important parameters determining behavior of fine particles in presence of a liquid phase.

1.1. Inter Particle Forces

Some particles in reservoir porous media are very small so that their static and dynamic conditions can be controlled by surface phenomena. These types of forces include London-Van der Waals attraction, electrostatic repulsion, Hydration, Hydrophobic and several other forces. In this study, only two major inter-particle forces; i.e. London-Van der Waals and electrostatic double layer are focused on. The force balance of these inter-particle forces could be shown as bellow:

(1)

where, subscripts of LVA and DLR stand for Attractive London-Van der Waals and Repulsive Double Layer forces, respectively. And:

where x, γ, T, , pH and r are the distance from surface of solid, ionic strength, temperature, permittivity, pH of solution, and radius of particles, respectively. Also [5] :

where, n and z are number and valency of i^th ion, respectively. is dielectric constant and is permittivity of vacuum.

1.1.1. Electrostatic Forces

One of the surface forces controlling behavior of Nano-sized materials in porous media results by ion dispersion in a manner that two layers are formed in the liquid phase that surrounds the solid; one interior molecular layer with a charge opposite of solid surface charge, with an exterior layer made by distribution of ions with a charge opposite to the first layer. The second layer fades away with increasing distance from the solid surface [5] . These two layers are schematically shown in Figure 2. It should be noted that the whole solution is neutral; however, this type of ion distribution makes the liquid to have electrical potential in points close to the solid surface.

Instantly beyond the surface of a charged particle, opposite ions are attached and move as the particle moves (stern layer). Then, ions with the same sign as surface of particle slip over the first layer (slipping or diffuse layer) and their concentration decrease with distance from solid surface.

Figure 1

Water injection history in Iranian Siri oil field

Figure 2

Schematic of a double layer in vicinity of a negatively charged surface and ion distribution of negative and positive ions around the surface. The induced potential decreases vs. distance from the face of charged plate

Formation of these two layers leads to definition of double layer and double layer force. This force can be estimated using surface potential. Surface potential could be estimated by a parameter named zeta potential that is an electrical potential in a certain distance from surface of solid. Zeta potential is estimated by experimental methods. Here, the surface charge value was assumed to be in a normal range. Particle is considered to be made of the same matter that makes reservoir rock. Thus, the potential of particle surface is equal to pore wall potential because they are both surrounded by a same liquid.

Two assumptions are made in this study to determine electrostatic forces in porous media include:

• The force between a particle and pore wall considered as a force between sphere and plate. Potential of sphere and plate are constant. This force can be measured by Equation (2) [6] . Figure 3 shows a schematic indication of ion distribution that causes the related forces.

• All pores are wide enough that double layers induced by two close walls do not overlap; i.e. zeta potential of solution is generated due to one surface charge.

(2)

where, r_p is particle radius, k_B is Boltzmann’s constant, e is electron potential, and refer to surface potential of plate and sphere, and κ is the reciprocal Debye length.

According to the fact that surface charge depends on ionic strength and pH, therefore electrical forces and properties of double layer should be functions of these two parameters as well. Subsequently, the possibility of releasing fines from a wall would be strongly affected by ionic strength and pH.

1.1.2. London-Van der Waals Force

This force is due to coupling of electron clouds around adjacent atoms. The previous assumption of considering fines as spheres and pore wall as plates is also made here. The force between a sphere and a plate can be calculated using Equation (3) [7] .

(3)

where, A is the Hamaker constant; which is a function of permittivity of media, temperature, volume and molecule types. Values of A is tabulated in several texts [8] . Hamaker constant used in this study is an average for reservoir rock minerals and estimated about 5 × 10⁻²⁰ j.

1.2. Migration of Fines in Porous Media

Fluid flow through the pores makes several forces that could impact the movement of fines in the media. As shown in Figure 4 these forces are: 1) electrical forces, F_e; 2) drag force; F_d, 3) lifting (buoyance) force, F_l, and 4) gravity, F_g [9] .

Erosion number can be defined as Equation (4):

(4)

where l_n and l_d are normal and drag levers as shown in Figure 4. When the erosion number reaches unity, the particle is in an unstable condition and increasing the numerator of Equation (4) leads to releasing the particle. When this number is lower than unity the fine grips to the surface. These forces should act in a way that particle do not release from pore surface. i.e. (F_DLR + F_l)l_n + (F_d)l_d < (F_VDA + F_g)l_n.

E_DLR and E_VDA can be obtained using Equations (2) and (3), respectively. Other forces are determined as follows.

1.2.1. Drag Force [10]

For low amounts of Reynolds number (i.e. Re < 1) the drag force predominated by frictional force. Then:

(5)

where, μ is water viscosity, ω is proportionality factor (in the range of 10 - 60), H and h_c are the hight of pore and deposited particles on the pore wall, respectively.h_c is asseumed to be equal to only one layer of fine particle (2r_p). U and φ are the velocity of fluid and porosity, respectively.

Figure 3

Ion distribution that causes the Van der Waals Atraction and double layer repulsion forces between suspended particle and pore wall

Figure 4

An indication of forces acting on a Nano-sized particle as a result of fluid flow through porous media: Electrical force (F_e), Lifting force (F_l), Drag force (F_d), Gravity force (F_g) and normal (l_n) and drag (l_d) levers (Figure from [8] )

1.2.2. Lifting Force [11]

(6)

where ρ_w is water density.

1.2.3. Gravity Force

(7)

where g is gravitational constant and ρ_p is density of the particle.

2. Composition of Displacing Phase

As can be seen from Equations (2) and (3), electrostatic double layer force is a function of ionic strength and fluid composition because of its dependency on Debye length, κ, and permittivity,. Dependency of Van der Waals force on fluid composition is described by the Hamaker constant, which is a proportionality factor for forces caused by the transient induced dipoles associated with interatomic bonds [12] . Calculations made in this study are based on the data summarized in Table 1.

3. Results and Discussion3.1. Debye Length

Debye Length is a function of ion concentrations and their valences. It changes with alteration of solution composition as is discussed in the following sections.

3.2. Ionic Strength Effect

Figure 5 shows the impact of valance and concentration on Debye length caused a charged solid. As can be seen, increasing valence of ions increases the Debye length (i.e. decreasing diffuse layer distance from surface of solid). This implies that density of ions close to the solid surface increases when valence increases.

According to the results of Figure 5 that shows increasing ionic strength increases Debye length, Figure 6 shows that by increasing salinity double layer force between the particle and plate decreases. This value is a function of valence of ions. It should be noted that, at very low concentrations the forces of different valences show different trend compared to the general trend (that the higher ionic strengths yield higher double layer repulsive). Another important point is that according to Equation (2) at lower concentrations (0.0037 for 1:1, 0.0044 for 2:2, 0.0023 for 2:1, 0.0053 for 3:1, and 0.0023 for 3:2 solutions) the double layer force acts as an attractive force.

Table 1

. Values of Parameters used in this study

Parameter	Value	Parameter	Value
Surface Potential of Wall, ψ₀₁ [V]	0.100	Bulk Density of Solid, ρ_s [kg/m³]	2200
Surface Potential of Fine, ψ₀₂ [V]	0.100	Porosity, [percent] Φ	20
Hamaker Constant, A [j]	10⁻²⁰	Displacing Phase Velocity, V [m/s]	2 × 10⁻³
Boltzmann Constant, k_B [j/k]	1.38 × 10⁻²³	Particle Radius, r_p [m]	4 × 10⁻⁷
Permittivity of Media, ε	80	Channel (pore) Height [m]	1 × 10⁻⁶
Electron Charge, e [C]	1.602 × 10⁻¹⁹	Salinity [molar]	0 - 2
Temperature, T [k]	268.15	Ionic Valence	1, 2, 3
Water Viscosity, μ_w [cp]	1	Distance of Particle and Wall [m]	10⁻¹⁰
Water Density, ρ_w [kg/m³]	1000.1	-	-

Figure 5

Effect of ionic Strength as a function of concentration and ionic valence on Debye length expressed as [1/m]

Figure 6

Effect of ionic strength on double layer forces; at very low concentrations the behavior is in contrary with the general trend

3.3. Distance between Particle and Pore Wall

Finally, effect of distance between particle and plate is shown in Figure 7. This factor impacts entrapping of a moving particle to adsorb on pore wall and also releasing a particle according to its position on the wall.

3.4. pH Effect

H⁺ and OH⁻ are two potential determining components [5] , which can be predominant to define potential of solid surface. According to some of wettability studies, charge of some materials is not a linear function of pH. In experimental results of some researchers [6] [13] show that by increasing pH, contact angle (as a function of zeta potential) reaches to a certain point and starts to alter its trend.

3.5. Properties That Affect Fines Migration

As explained by Equations (2) to (8), separation of a Nano-sized particle from the solid face of pore wall as a result of fluid flow (because of oil/gas production or injection processes) is predictable using Equation (4). In this work, changes of ionic strength (valence and concentration), fluid velocity (flow rate), and pore wall―fine spacing are studied.

Figure 7

Decreasing double layer force by increasing distance from the charged surface

Separation of a particle from a pore wall and its migration happens when erosion number approaches to unity (where the repulsive forces overcome the attractive forces). Double layer force is determined by Debye length, which is a function of ionic strength. Thus, changing salinity would affect erosion number and consequently im- pacts fines migration. Figure 8 shows effect of ionic strength on erosion number; constant parameters are those in Table 1.

Figure 9(a) to Figure 9(f) show the effect of velocity (injection rate) on erosion number in different ionic strengths. During injection processes, an important factor is to optimize injection rate to minimize formation damage caused by fines migration.

In these figures, the velocities under the interception line of erosion number vs. velocity curve in are under the critical values and do not cause fines migration. Velocity is not a crucial parameter under the conditions in Figure 9(a) to Figure 9(f), since the curve is either bellow (i.e. migration does not happen, even under high injection rate conditions; consider limited range of fluid injection velocity in pore scale) or above (i.e. migration happens; presence of liquid causes separation of fine from pore wall, even if the liquid phase is under static condition). Figure 10 shows another case that in a critical velocity erosion number becomes more than 1. This shows that in some conditions the slightest changes of injection rate (velocity) can cause fines migration and reservoir damages.

4. Conclusions

Inter-molecular forces have important roles in determining behaviors of material in different fields. The results of this study explain circumstances of fine separation from a solid surface that is defined as fine migration in porous media. This phenomenon could cause serious damage in hydrocarbon reservoirs due to blocking pores and permeability reduction.

There are several parameters that effect migration of fines. These parameters include reservoir formation characteristics and injected agent properties. Some of these affecting parameters such as ionic strength (ion valences and their concentrations), velocity of injected fluid, particle and pore wall spacing and their influences on discussed phenomena were studied in this paper. The impact of each parameter on the electrical properties of media and the resultant forces that could cause releasing of fines from a wall is summarized as follow:

• Debye length increases with ionic strength of solution. That is, the ions are concentrated around the charged face and decrease the distance between diffuse layer and solid face.

• Increasing Debye length reduces the double layer repulsion; it means that at higher ionic strengths, repulsion is lower. This implies that in unconsolidated reservoir with high concentration of fines that inclined to migrate, high salinity water should be injected to avoid reservoir damage.

• By balancing involved forces and definition of a quantity named erosion number that is a ratio of forces tend

Figure 8

Erosion number change as a function of ionic strength. The values intersecting with ε = 1 imply critical concentration

Figure 9

Erosion number changes as a function of velocity for different ionic strengths

Figure 10

Velocity importance under a specific condition

to separate a particle to the forces try to settle it, this is predictable that whether a particle will release from a surface of solid or not. The particle stands until the erosion number reaches to unity; i.e. separating forces overcome to the settling forces.

• Dependency of the introduced forces to different parameters causes the erosion number to be different for different reservoir conditions and different displacing fluid properties. As shown in the plots, though increasing fine size causes increase in double layer repulsive force, however, according to presence of other forces, fines migration decreases.

• Injection rate is one of the most important parameters in EOR processes. The salinity of the displacing phase has also a key role in final recovery. Thus, to eliminate reservoir damages a study of injection rate and suitable salinity should be performed. The general trend is that by increasing the ionic strength the maximum allowable velocity increases.

References1

BEDRIKOVETSKY, P., FURTADO, C.A., DE SOUZA, A.L.S. AND SIQUEIRA, F.D. (2007) INTERNAL EROSION IN ROCKS DURING PRODUCED AND SEAWATER INJECTION. SPE 107513, SPE EUROPEC/EAGE ANNUAL CONFERENCE AND EXHIBITION, LONDON, 11-14 JUNE 2007.

KWAN, M.Y.M., CULLEN, M.P., JAMIESON, P.R. AND FORTIER, R.A. (1988) A STUDY OF FINES MIGRATION RELATED PERMEABILITY DAMAGE IN EXTRACTED COLD LAKE HEAVY OIL CORES. 39TH ANNUAL TECHNICAL MEETING OF THE PETROLEUM SOCIETY OF CIM, CALGARY, 12-16 JUNE 1988.

VALDYA

R.N.

, FOGLER

H.S.

,et al. (1992)VALDYA, R.N. AND FOGLER, H.S. FINES MIGRATION AND FORMATION DAMAGE: INFLUENCE OF PH AND ION EXCHANGE SPE PRODUCTION ENGINEERING 7, 325-330.

MOGHADASI, J., JAMIALAHMADI, M., MÜLLER-STEINHAGEN, H., SHARIF, A., IZADPANAH, M.R., MOTAEI, E. AND BARATI, R. (2002) FORMATION DAMAGE IN IRANIAN OIL FIELDS. SPE 73781.

ADAMSON, A.W. AND GAST, A.P. (1997) PHYSICAL CHEMISTRY OF SURFACES, CHAPTER 5. A WILEY-INTERSCIENCE PUBLICATION, 170.

TAKASHI, S. (2009) WATER IMBIBITION, ELECTRICAL SURFACE FORCES, AND WETTABILITY OF LOW PERMEABILITY FRACTURED POROUS MEDIA. PH.D. THESIS, STANFORD UNIVERSITY, STANFORD.

ROSS, S. AND MORRISON, I.D. (1988) COLLOIDAL SYSTEMS AND INTERFACES. JOHN WILEY PUBLICATION, NEW YORK.

FRENCH

R.H.

, MUÈLLEJANS

, JONES

D.J.

, DUSCHER

, CANNON

R.M.

, RUÈHLE

,et al. (1998)DISPERSION FORCES AND HAMAKER CONSTANTS FOR INTERGRANULAR FILMS IN SILICON NITRIDE FROM SPATIALLY RESOLVED-VALENCE ELECTRON ENERGY LOSS SPECTRUM IMAGING ACTA MATERIALIA 46, 2271-2287.HTTP://DX.DOI.ORG/10.1016/S1359-6454(98)80008-6

ZEINIJAHROMI, A., PHUONG, T.N. AND BEDRIKOVETSKY, P. (2011) TAKING ADVANTAGE OF FINES-MIGRATION-INDUCED FORMATION DAMAGE FOR IMPROVED WATERFLOODING (RESERVOIR SIMULATION USING POLYMER FLOOD OPTION). SPE 144009, SPE EUROPEAN FORMATION DAMAGE CONFERENCE, NOORDWIJK, 7-10 JUNE 2011.

FARAJZADEH, R. (2004) PRODUCED WATER RE-INJECTION (PWRI) AN EXPERIMENTAL INVESTIGATION INTO INTERNAL FILTRATION AND EXTERNAL CAKE BUILD UP. M.SC. THESIS, DELFT UNIVERSITY.

KHILAR

K.C.

, VAIDYA

R.N.

, FOGLER

H.S.

,et al. (1990)COLLOIDALLY-INDUCED FINES RELEASE IN POROUS MEDIA JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING 4, 213-221.HTTP://DX.DOI.ORG/10.1016/0920-4105(90)90011-Q

HOEIL

, S.

BARTH

, T.

BLOKHUS

, A.M.

SKAUGE

, A.

,et al. (2001)THE EFFECT OF CRUDE OIL ACID FRACTIONS ON WETTABILITY AS STUDIED BY INTERFACIAL TENSION AND CONTACT ANGLES JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING 30, 91-103.HTTP://DX.DOI.ORG/10.1016/S0920-4105(01)00106-1