Modeling and Analysis of Magnetic Nanoparticles Injection in Water-Oil Two-Phase Flow in Porous Media under Magnetic Field Effect

1College of Engineering, Effat University, Jeddah 21478, Saudi Arabia 2Computational Transport Phenomena Laboratory (CTPL), Division of Physical Sciences and Engineering (PSE), King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia 3Mathematics Department, Faculty of Science, Aswan University, Aswan 81528, Egypt 4Faculty of Engineering, University of Regina, Regina, SK, Canada


Introduction
Industry is now looking seriously into using nanotechnology as a viable tool to solving new challenges in several fields.In particular, there has been interest among oil and gas production companies to explore using nanotechnology in solving challenges related to unconventional oil and gas reservoirs, such as those found in tight and shale formations [1][2][3][4][5][6][7].Nanotechnology has been used in different areas of the oil and gas industry from exploration, drilling, production to reservoir monitoring, and refining.The conventional Enhanced Oil Recovery (EOR) methods have several problems from high cost to low oil recovery in addition to operations problems especially in thermal and chemical methods.The nature of nanoparticles results in some useful characteristics such as increased surface area, which at the nanoscale size, does matter when it comes to how molecules react to and bond with each other.So, for example, nanoparticles can be used in EOR, because they are small enough to pass through pore throats in typical reservoirs, and they can be retained by the rock.Ju and Fan [8] calibrated a model for nanoparticles transport in two-phase flow in porous media based on the formulation of the colloid model of fine particles transport in two-phase flow in porous media [9].El-Amin et al. [10][11][12][13] have presented modeling and simulations of nanoparticles transport associated with two-phase flow in porous media.On the other hand, experimental studies of using nanoparticles in EOR have been conducted by Suleimanov et al. [14] and Hendraningrat et al. [15].
One of the prospective applications of nanotechnology is nanoferrofluids, as the flow of such fluids can be controlled by the introduction of external magnetic field.This opens the way for various applications from directing flow in reservoir monitoring, diverting flow in acid jobs to control and boost of injection fluid advancement during pressure maintenance to increase oil recovery.The idea of using a strong external magnetic field with large magnetic susceptibility fluid is to mobilize ferrofluid through porous media.Oil recovery can be increased by using nanoparticles with electromagnetic properties (such as iron oxide, Fe 2 O 3 , and zinc oxide, ZnO) under waves generated from an electromagnetic source.Both direct and alternating magnetic fields are under investigation; in our case here we will focus on direct magnetic field.The magnetization of the particles and their attraction toward the magnet causes flow of the magnetic particles suspension.The movement of the magnetic nanoparticles under the magnetic field effect is independent of the orientation of the magnet.In the last few years, a number of publications have been considered nanoferrofluids in oil and gas recovery or environmental applications (e.g., [16][17][18][19][20][21][22][23]).
Borglin et al. [24] conducted experiments to measure the magnetic induction, which converted to magnetic field strength, at various distances in a direction aligned with the poles.McCaig and Clegg [25] presented equations which indicate that ferrofluid magnetization is variable at all locations far away from the magnet due to the decreasing magnetic field strength.The gradient of magnetic field strength varies also with distance from the magnet.Moreover, the force on the ferrofluid decreases with distance from a magnet.Therefore, we may conclude that the distance from the magnet is important and should be reconsidered with taking into consideration its location in particular in the core-scale as well as the direction of flow.On the other hand, it is known that, in the presence of an external magnetic field, the nanoparticles in ferrofluid become magnetized and are pulled toward a magnet.In this work, we attempt to discuss the effect of location of the magnet on the transport of the nanoparticles in porous media.
In the current work, we develop a mathematical model to describe the magnetic nanoparticles-water suspension imbibition into an initially oil saturated porous domain under magnetic field effect.The porous medium is considered initially saturated totally with oil except for a residual amount of the other phase.We consider countercurrent imbibition into a small-scale porous core.This countercurrent imbibition refers to the case in which all the porous medium domain boundaries have no flow except one side.Physical variables are investigated under the influences of magnetic field with two different locations of magnet, namely, right and left to the porous core.Numerical experiments for the two cases are performed and results are introduced in graphical representations.

Modeling and Mathematical Formulation
Consider suspension of magnetic nanoparticles injected in an isothermal incompressible water-oil two-phase flow under an external magnetic field.In the following, we describe the mathematical modeling of the problem under consideration.
In the following subsections, we firstly introduce the magnetic force and other magnetic modeling.The magnetic body force, which acts as a body force on the nanoparticles suspension per unit volume, appears in the extended Darcy's law as presented in the second subsection.Other magnetic parameters such as magnetization and magnetic field strength are also presented in the first subsection.Then, the governing flow equations such as momentum (extended Darcy's law) and mass conservation (saturation equation) are provided in the second subsection.The third subsection is devoted to the modeling of nanoparticles transport which is coupled with the flow equations via velocity and saturation.Finally, the initial and boundary conditions are presented in the last subsection.

Magnetic Force.
The magnetization (polarization) of the nanoparticles suspension interacts with the external permanent magnetic field to produce attractive forces on each particle.The external magnetic force acts as a body force on the nanoparticles suspension per unit volume which can be expressed as [16] where  0 is the magnetic permeability,  is the magnetization, and  is the magnetic field strength.The magnetization  is a function of , approximated by where the parameters  1 and  1 depend on the particular type of the ferromagnetic material.The values of the initial susceptibility and the saturation magnetization are controlled by the parameters  1 and  1 , respectively.The larger  1 leads to a larger initial susceptibility which corresponds to larger particles or agglomeration of particles.The range of  1 may be of order 10 4 -10 5 , while the order of  1 may be of order 10 −6 -10 −5 .The magnetic field strength in 1D may be written as follows [16]: where   is the residual magnetization (in this work,   = 1.19 [T]) and  is the distance between the poles of the magnet.In Figure 1,  mag ,   , and   / are plotted against .
Oldenburg et al. [16] reported that magnetic particles (ferrofluids) could have both miscible and immiscible behaviors with water.For example, an oil-based ferrofluid is immiscible with water while aqueous ferrofluids are miscible in water.The immiscible case should be described by twophase aqueous/nonaqueous mixtures regardless of external magnetic field strength, while the case of miscible in water should be described as single-phase aqueous mixtures.Under no external magnetic field, the two fluids are miscible.When a magnetic field is applied, it affects the ferrofluid part of the aquatic mixture in such a way that it behaves as immiscible.
That is, the ferrofluids motion is dictated by the magnetic field.When the magnetic field is relatively strong, its effect on the ferrofluids may be influential.Therefore, a two-phase flow model is suitable for the flow of ferrofluids under strong external magnetic fields, but when the magnet is at larger distances from the flow (thus segregation effects are unimportant), the fully miscible single-phase liquid model is appropriate.Both of these models assume that the mixture is not so dilute.Pure ferrofluid and dilute colloidal suspension have approximately a thousand particles per cubic meter (see [16,26]).So, one may conclude that the miscible and immiscible models may be applicable for ferrofluid fractions of 10% in water, and they are probably satisfactory in much more dilute mixtures [16].
In the current model, we treat the water-nanoparticles suspension as a miscible mixture while it is immiscible with the oil phase.Oldenburg et al. [16] assumed that magnetization increases linearly with the ferrofluid mass fraction ,  () =  ( = 1) . (4) In the case of the immiscible two-phase (water-oil) model, Oldenburg et al. [16] assumed that magnetization increases linearly with the ferrofluid saturation, For the immiscible conceptualization, density and viscosity are considered properties of the respective phases, and no mixing relations are required.Therefore, in the current model we consider immiscible two-phase (water-oil) model; however, the water-nanoparticles suspension is a miscible mixture.So, one may expresses the the magnetization of this system as, Following the work of Reeves and Cranwell [27] and Herbert et al. [28] on brine-water mixtures, Oldenburg et al. [16] have assumed that the volumes of pure water and ferrofluid are additive and define for the mixture density as where  , [kg⋅m −3 ] is the density of water component and The viscosity of the nanoparticles-water mixture is calculated by the following linear relationship [16,28]: where  , is the viscosity of pure water.

Flow Model.
Considering the external magnetic force, which acts as a body force as explained in the above subsection, the vertical 1D (i.e., -direction) velocities of oil phase and nanoparticles-water mixture phase may be written as The mass conservation equation is where  [-] is the porosity,   [kg⋅m −3 ] is the density of phase ,   [-] is the saturation of phase , and u  [m⋅s −1 ] is the velocity of the phase . stands for the nanoparticles-water suspension phase, and  stands for the oil phase. [m 2 ] is the permeability,   [-] is the relative permeability of phase ,   [Pa] is the pressure of phase ,  [m⋅s −2 ] is the gravitational acceleration, and   [kg m −1 s −1 ] is the viscosity of phase .
The fluid saturations for the two-phase flow of water and oil are related by In countercurrent imbibition, the sum of the velocities of the wetting and nonwetting phases is zero, that is, The capillary pressure is defined as Since, in this case of study, we have two phases, one magnetized and the other nonmagnetized, the magnetostrictive effect should not be neglected.The magnetized phase pressure will have additional pressure term with the conventional thermodynamic pressure.The resulting pressure is called composite pressure [21], which can be given as where   is the ferrofluid phase dynamic pressure,   is the fluid magnetic pressure,   is the magnetostrictive pressure, and   is the magnetic normal pressure, which is neglected [21] in this study.Also, one can write The specific volume is given as So, where We may approximate  1 / and  1 /, when  1, and  1, are equal to zero in the nonmagnetized phase, as Therefore, Similarly, therefore, Substituting from ( 9), ( 13), (20), and ( 21) into (12) and eliminating   /, we may have So, the water velocity becomes Therefore, the saturation equation for the water phase becomes where where   is the entry pressure for the imbibition.Moreover, the relative permeabilities are defined as where  is the normalized water phase saturation, which is given as and  0  =   ( = 1) and  0  =   ( = 0) are the endpoint relative permeability of the water and oil phase, respectively. 2 and  2 are positive numbers.  is the irreducible water saturation and   is the residual oil saturation.

Nanoparticles Transport Model.
Assuming that the nanoparticles exist only in the water phase and have only one size interval.The Brownian diffusion is considered for the nanoparticles and the volumes of pure water and the particles suspension are additive.The transport equation of the nanoparticles-water suspension in the water phase can be written as where  [m 3 ⋅m −3 ] is the volume concentration of nanoparticles in the water phase. 1 [m 3 ⋅m −3 ] is the volume of the nanoparticles in contact with the water phase available on the pore surfaces per unit bulk volume of the porous medium. 2 [m 3 ⋅m −3 ] is the volume of the nanoparticles entrapped in pore throats from the water phase per unit bulk volume of porous medium due to plugging and bridging. diff [m 2 ⋅s −1 ] is the molecular diffusion coefficient which can be calculated using the Stokes-Einstein equation, and  is the tortuosity of the flow.  is the Boltzmann constant;  [K] absolute temperature.The mechanical dispersion coefficient  disp [m 2 ⋅s −1 ] is a function of Darcy's velocity and may be given as where  , is the longitudinal dispersion coefficient.
The modified Gruesbeck-Collins model [30] for the surface deposition is used in this study [8,12].In the presence of a critical velocity of the surface deposition only particle retention occurs while above it retention and entrainment of the nanoparticles take place simultaneously, which can be modeled as Also, the rate of entrapment of the nanoparticles in the water phase is given by where   [m −1 ] is the rate coefficient for surface retention of the nanoparticles in the water phase.  [m −1 ] is the rate coefficient for entrainment of the nanoparticles.  is the critical velocity for the water phase.  [m −1 ] is the pore throat blocking constant.
The porosity variation due to nanoparticles deposition is given as [8,9] where  0 is the initial porosity.Therefore, permeability can be changed as [8] where  0 is the initial permeability and  = I,  0 =  0 I, where  [m 2 ] and  0 [m 2 ] are positive constants and I is a unit matrix.  is constant for fluid seepage allowed by the plugged pores.The flow efficiency factor expressing the fraction of unplugged pores available for flow is given by where   is the coefficient of flow efficiency for the nanoparticles.The exponent  has a value of the range from 2.5 to 3.5.The variation of relative permeabilities due to the nanoparticles reposition is given by El-Amin et al. [12] as follows: where   is the ratio of the phase  relative permeability due to nanoparticles adhering; that is, where  , is the relative permeabilities of water/oil phase when the surfaces per unit bulk volume of the porous media is completely occupied by the nanoparticles.  =  tot / sp such that is a specific area of the sand core and  [m 2 ] is the crosssectional area.Also, the total surface area in contact with fluids per unit bulk volume is defined as [8] and  is the diameter of the nanoparticles in a given interval size.

Initial and Boundary Conditions.
At the beginning of the flow the following initial conditions are considered: where ℎ is the rock depth and  0  is the initial water saturation.The boundary conditions are where  0 is the concentration of nanoparticles-water suspension on the inlet boundary.
In the following we consider two different locations of the magnet, namely, right and left side, with respect to the core and its boundary conditions (see Figures 2 and 3).On other words, the magnet location and flow directions are very important factors in this problem.If the magnet location is on the right side of the no-flow boundary of the core, the magnetic field will be assisting the flow, while the opposite is true when the magnet is placed on the left side of the inflow boundary of the rock (opposing flow).4-9. Figure 4 illustrates the nanoparticles-water saturation profiles against the rock length with and without magnetic field effect for various values of the time of imbibition.It is interesting to note from this figure the effect of placing the magnet on the right side of the rock as it increases the saturation of nanoparticleswater suspension on this side while it decreases the saturation on the left side of the rock near to the flow inlet.This can be explained through an important fact that is the magnetic force is proportional to magnetic field strength, so fluid that is closer to the magnet is more strongly magnetized and pulled more strongly toward the magnet.For instance, this is in contrast to gravity driven flow in which gravitational force is independent of position and fluid is pulled uniformly under gravity.

Magnet at the Right Side of the Core. Firstly, we present the first case in Figures
The nanoparticles concentration is plotted against the rock length in Figure 5, with and without magnetic field effect for various values of the imbibition time.From this figure, it can be seen that after long time of imbibition,  the nanoparticles concentration decreases slightly under the effect of the magnetic field.Figure 6 shows the concentration of deposited nanoparticles on the pore wall as plotted against the core length with and without magnetic field effect for various values imbibition times.Also, Figure 7 shows the concentration of deposited nanoparticles on the pore throat which is plotted against the core length with and without magnetic field effect for various values imbibition times.As expected the deposited nanoparticles have an opposite behavior of the nanoparticles in the water; that is, when the nanoparticles concentration in the water decreases, the deposited nanoparticles concentration increases.
In Figure 8, the permeability variation is plotted against the distance, with and without the magnetic field effect for various values of imbibition time when the magnet is located on the right side of the core.It is clear from this figure that the reduction in permeability is related to the quantity of the precipitation of nanoparticles on the wall and the imbibition time.Also, in Figure 9, the porosity variation is plotted against the distance with and without the magnetic field effect for various values of imbibition time when the magnet is located on the right side of the core.It can be seen that the change in porosity is related to the quantity of nanoparticles plugging the throat of the pore and imbibition time.

Magnet at the Left
Side of the Core.Now, we consider the second case, when the magnet is to the left of the rock (which opposes the flow), in Figures 10-13.Saturations  profiles are plotted in Figure 10, against the core distance with and without magnetic field effect for various values of imbibition time when the magnet is placed on left side of the core.An interesting behavior can be seen in this figure.Unlike the case of the right magnet, the left magnet resists the flow and reduces water invasion.Moreover, Figure 11 shows the profiles of nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time, in the case of left location of the magnet.It is clear from this figure that the concentration of nanoparticles in water increases especially after long imbibition time.Figure 12 illustrates deposited nanoparticles concentration on the pore wall against the core length, with and without magnetic field effect for various values of imbibition time when the magnet is located on the left side of the core.This figure shows that the magnetic field reduces significantly the deposited nanoparticles on the pore surface.Finally, Figure 13 illustrates the nanoparticles concentration plugging pore throat against core length with and without magnetic field effect for various values of imbibition time when the magnet is located on the left side of the core.One may notice that the magnetic field reduces the nanoparticles plugging pore throat.
Figure 14 illustrates the permeability profiles against the distance with various values of imbibition time, with and without magnetic field effect when the magnet is on the left side.Similarly, the permeability decreases as a consequence of the precipitation of nanoparticles.The permeability has a minor reduction compared to the initial permeability.Figure 15 shows the porosity variation against the distance with various values of imbibition time, with and without    magnetic field effect when the magnet is located on the left side.It is interesting to note reduction in the porosity which is a consequence of the precipitation of nanoparticles on the pore walls.The porosity variation is negligibly small.

Conclusions
This paper was devoted to study the magnetic field effects on the the transport of magnetic nanoparticles injected into a two-phase water-oil system in porous media.Firstly, we developed the mathematical model which governs the transport of magnetic nanoparticles in a two-phase, wateroil system in porous media under the effect of an external magnetic field.The countercurrent imbibition in a smallscale porous medium core is considered as an example.Both Brownian diffusion and mechanical dispersion are taken into consideration, with the assumption that only one size interval of the nanoparticles exists.Variation of porosity, permeability, and relative permeabilities due to particles deposition are also considered.Immiscible mixture relationships have been used to determine the magnetization properties, density, and viscosity of this ferrofluid.In such flow field, the pressure includes additional terms, namely, the ferrofluid dynamic pressure, the fluid magnetic pressure, and the magnetostrictive pressure.Two different locations of the magnetic field with respect to the core inflow and no-flow boundaries are considered.Effects of the magnetic field on the physical variables such as saturation, nanoparticles concentrations, porosity, and permeability have been examined.The important result from this simulation is that the magnetic source location has a significant influence on the physical variables of the model.Based on the flow direction and the location of the magnet, the magnetic field can assist or oppose the flow of this two-phase system.Such observed behavior can be used for different applications depending on whether the goal is to aid or delay the injected fluid advancement.According to this investigation, if the magnet is placed next to the no-flow boundary of the studied countercurrent numerical experiment, the saturation of nanoparticleswater suspension increases.Moreover, the concentration of nanoparticles is observed to decrease slightly, which has been attributed to the slight increase of the deposition of nanoparticles.On the other hand, when the magnet is placed next to the inflow/outflow boundary, the magnet resists the flow of the ferrofluid suspension and decreases water invasion capacity, thereby.Furthermore, the nanoparticles concentration seems to increase under the effect of the magnetic field, and the deposited nanoparticles concentration decreases.Both porosity and permeability are reduced due to the nanoparticles adherence to the walls of the porous media.

Figure 4 :
Figure 4: Saturations profiles against the core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.

Figure 5 :
Figure 5: Nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.

3 •m − 3 )Figure 6 : 5 T 3 •m − 3 )Figure 7 :
Figure 6: Pore wall deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.

Figure 8 :
Figure 8: Permeability variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.

Figure 9 :Figure 10 :
Figure 9: Porosity variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the right side of the core.

Figure 11 :
Figure 11: Nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.

3 •m − 3 )Figure 12 : 3 •m − 3 )Figure 13 :
Figure12: Pore wall deposited nanoparticles concentration against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.

Figure 14 :Figure 15 :
Figure14: Permeability variation against core length with and without magnetic field effect for various values of imbibition time: magnet is placed on the left side of the core.
=   /  and   =   /  are mobility ratios of water and oil phases, respectively.=   +   is the total mobility.=   /  and   =   /  are the flow fraction of water and oil phases, respectively.Δ =   −   .The capillary pressure is a function of the normalized saturation which can be given as[29]