Open Access
Issue
Oil & Gas Science and Technology - Rev. IFP Energies nouvelles
Volume 73, 2018
Article Number 2
Number of page(s) 12
DOI https://doi.org/10.2516/ogst/2017038
Published online 12 January 2018

© S. Kashefi et al., published by IFP Energies nouvelles, 2018

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Introduction

Crude oil is composed of light to heavy hydrocarbons that based on the polarity and molecular structure is separated to four groups: Saturates, Aromatics, Resins, and Asphaltenes. The asphaltenes are the heaviest and polar cyclic components of crude oil that are insoluble in low molecular weight alkanes and soluble in light aromatics such as toluene [1]. One common problem during oil recovery is asphaltene instability due to changes of thermodynamic conditions. During oil production, the asphaltene flocs are formed slowly because of self-association and precipitated out of the oil phase. As a result, the asphaltene is deposited on the rock surfaces and effects on the rock properties such as permeability, porosity and wettability. It causes formation damage during oil production and transportation [2]. One of the effective methods for controlling asphaltene deposition is to use chemical additives such as surfactants [3,4], polymeric inhibitor [5] and adsorbents [6]. Up to now, several adsorbents were used for asphaltene adsorption from oil solution such as mineral surfaces (kaolin, calcite and dolomite) [6], clay and reservoir rocks [7,8], metal surfaces (stainless steel and gold) [9,10] and several kinds of metal and metal oxide nanoparticles (MgO, TiO2, NiO, CaO, Fe3O4, Al2O3, etc) [1118]. The reasons of interest on using nanoparticles in recent years for asphaltene adsorption are the special properties of nanoparticles such as their high surface area to volume ratio, functionalizable surface and their smaller size that cause easier transportation in porous media and enhance the performance of fluid flow and oil recovery [1416]. The amount of adsorbed asphaltene onto nanoparticles strongly depends on the type, surface chemistry of the nanoparticles and the kind of the interaction force between asphaltene and nanoparticles [11,12]. The asphaltene surface charge, Van der Waals force and acid-base interaction between surface of nanoparticles and polar asphaltene particles are different kinds of forces for asphaltene adsorption onto nanoparticles [11,16,19].

Most of the previous studies showed that the nanoparticles are good choices for adsorption and catalytic oxidation of asphaltene in crude oil upgrading process [11,13,19] while, limited researches have been carried out on using nanoparticles in porous media for asphaltene adsorption and consequently reduction of asphaltene deposition through oil flooding [17,18,20].

The objective of this study was to investigate the potential application of zeolite beta nanoparticles as a new adsorbent for asphaltene adsorption from crude oil solution which reduces the asphaltene deposition during oil injection in a porous media. At first the capacity of zeolite beta nanoparticles for equilibrium adsorption of asphaltene was determined by UV-Vis spectrophotometer. Then two dynamic experiments including co-injection of n-heptane (as an asphaltene precipitant) and crude oil into a silica sand pack were conducted in the absence and presence of the zeolite beta nanoparticles. The changes of permeability and porosity ratios in a porous media were experimentally determined with the use of zeolite beta. The process was also mathematically modeled by a developed asphaltene deposition model. To our best knowledge, this is the first report on the use of zeolite beta nanoparticles in porous media for investigating of asphaltene deposition during oil injection.

1 Experimental

1.1 Materials

A crude oil from an Iranian reservoir on the south west of Iran (API = 33.74) was used in this study. Table 1 shows the specifications of the crude oil, which has about 1 wt% asphaltene content. The asphaltene particles were extracted from the crude oil according to IP-143 standard analysis [21]. The acid and base numbers of asphaltene were found to be 2.25 and 12.46 by potentiometric titration method, respectively [22]. These values confirmed the basic nature of asphaltene particles. The n-heptane (Merck, Germany, >99%) as asphaltene precipitant and toluene (Merck, Germany, >99.8%) as solvent were used.

The zeolite beta nanoparticles (Si/Al = 25) was supplied from Zeolyst International Company in USA. The properties of the zeolite beta are reported in Table 2. Moreover, the liquid paraffin (by-product of refining crude oil of Iran, density = 0.88 g/cm3) was applied as a base fluid for dispersing nanoparticles.

Table 1

Specifications of the crude oil.

Table 2

Properties of the zeolite beta nanoparticles.

1.2 Estimation of Proper Amount of Nanoparticles for Asphaltene Adsorption

Several amounts of zeolite beta nanoparticles (1, 2.5, 5, 7.5, 10, 15 and 20 g/L) were added to the 1500 mg/L asphaltene-toluene solutions. The initial concentration of 1500 mg/L is the good choice for this experiment. Because it is lower than critical clustering concentration of asphaltene in toluene (≈2000 mg/L) [24]. Therefore, with selection of this concentration, the contribution of the asphaltene self-association to the adsorbed mass of the asphaltene was decreased [25]. Next, the prepared solutions were being shaken at 200 rpm until the equilibrium was achieved. Based on the performed preliminary experiments, 4 h shaking time was found to be sufficient to reach equilibrium and the amount of adsorbed asphaltene onto zeolite beta nanoparticles did not change with increasing contact time. After that, the nanoparticles containing the adsorbed asphaltene were separated by centrifuging at 5000 rpm for 30 min [14]. The light absorbance values of the residual solutions were measured by UV-Vis spectrophotometer (Perkin-Elmer Instruments Lambda 35) at the wavelength of 410 nm. The standard solutions with the specified concentrations of asphaltenes were used to create a calibration curve of UV-Vis absorbance and pure toluene was used as a blank [14,20].

To determine the proper amount of nanoparticles for asphaltene adsorption, the adsorption efficiency (E) was calculated by following equation [26]: E=(C0Ce)C0×100(1) where C0 and Ce are the initial and equilibrium concentrations of asphaltene in solutions (mg/L), respectively.

1.3 Batch Adsorption Experiments

The proper amount of zeolite beta nanoparticles was added to asphaltene-toluene solutions with concentration range between 100 and 3000 mg/L. The solutions were shaken at 200 rpm for 4 h. Then the nanoparticles containing adsorbed asphaltene were separated by centrifuging at 5000 rpm for 30 min. The amount of asphaltene in residual solutions was measured by UV-Vis spectrophotometer at the wavelength of 410 nm [14]. The adsorption of asphaltene onto nanoparticles was measured from the change in the concentration of asphaltene in toluene before and after adding the nanoparticles. The adsorbed amount of asphaltene onto the nanoparticles was calculated by [20]: qe=C0Cem(VAn)(2) where qe is the equilibrium adsorbed amount of asphaltene (mg/m2), V is the volume of solution (L), m is the dry mass of nano-particles (g), C0 is the initial concentration of asphaltene in the solution (mg/L) and Ce is the equilibrium concentration of asphaltene in the solution (mg/L), An is the surface area of nanoparticles (m2/g).

The Langmuir isotherm model was used to correlate the experimental equilibrium data of asphaltene adsorption onto zeolite beta nanoparticles [11,13,14,20,25]. The Langmuir model is presented as follows [27]: qe=qmaxkLCe1+kLCe(3) where qe is the equilibrium adsorbed amount of asphaltene (mg/m2), qmax is the maximum amount of adsorbed asphaltene for the monolayer adsorption on adsorbent surfaces (mg/m2), Ce is the equilibrium concentration of asphaltene in the solution (mg/L), and kL is the Langmuir equilibrium adsorption constant (L/mg).

1.4 Sand Pack Preparation

In this research a silica sand pack was prepared and used as porous media for investigating the asphaltene deposition process with and without the use of nanoparticles. At first, the silica sands (Tavansilice, Iran) were washed by deionized water to remove any impurities and placed in an oven at 333 K for 12 h to dry them. The silica sand particles had 297–595 micron diameter that were obtained by passing them through sieves with 30–50 mesh sizes. A stainless steel cylinder with 21 cm length and 3.5 cm inside diameter was used and filled with approximately 500 g of prepared silica sands by simultaneously packing and shaking. Figure 1 shows a schematic view of the experimental setup used in this work.

The setup consists of several sections such as: The silica sand pack; cylinders for transferring samples (water, n-heptane and crude oil with or without nanofluid); positive displacement pumps for allowing fluid displacement under constant flow rate; differential pressure transducer that measures the pressure differences between the inlet and outlet of the sand pack; Sampling vial to collect produced fluids; check valves and data acquisition system.

thumbnail Figure 1

Schematic view of fluid injection setup.

1.5 Nanofluid Preparation

For nanofluid preparation, the proper amount of zeolite beta nanoparticles was dispersed in liquid paraffin as a base fluid. The base fluid should be a clear liquid for better observation of nanoparticle stability. Nanoparticles dispersion procedure was conducted using a probe ultrasonic (400 W and 24 kHz) to break agglomerates of nanoparticles in the fluid and make stable suspension [28,29]. It was observed with naked eyes that the nanofluid was uniformly dispersed for 24 h and the complete sedimentation happened after three days. The size of zeolite beta nanoparticles in the liquid solution was about 50 nm, as measured by dynamic light scattering analyzer (DLS, Malvern Instruments, Ltd., Worcestershire, UK).

1.6 Dynamic Experiments

To analyze asphaltene deposition process during fluid injection into the sand pack and to investigate the effect of zeolite beta nanoparticles in this process, two dynamic experiments, DE-1 and DE-2, were conducted. The first dynamic experiment (DE-1) was co-injection of crude oil and n-heptane as an asphaltene precipitant into the sand pack. The second dynamic experiment (DE-2) was co-injection of n-heptane and crude oil in the presence of zeolite beta nanoparticles. The DE-1 experiment involved the following steps:

  • At first, the sand pack was vacuumed and saturated with water to determine the pore volume of porous media. How much water it has been adsorbed will be considered as pore volume (PV).

  • The initial porosity of the sand pack was calculated by dividing the PV to the total volume .

  • The absolute permeability was measured by water injection into the sand pack with constant flow rate. The pressure transducer was used to record the pressure differences (ΔP) between two ends of the sand pack. The absolute permeability was estimated by following Darcy's law [30] which was equal to 775 md.

k=(QΔP)(LμA)(4) where k is permeability (d), Q is fluid flow rate (cm3/s), ΔP is pressure differences (bar), μis fluid viscosity (cp), L is length of the sand pack (cm), A is cross-sectional area of the sand pack (cm2).
  • The crude oil with 1 wt% initial asphaltene content was injected into the sand pack by positive displacement pump at a constant flow rate for displacing the water. The oil injection was continued until no more water produced in outlet and the irreducible water saturation was found by subtracting the outlet volume of water from initial volume of saturated water in the sand pack (swi ≈25%).

  • The crude oil was injected into the sand pack at a constant flow rate for measuring the initial oil permeability (k0 = 469md).

  • The co-injection of crude oil and n-heptane was started inside the sand pack at same injection rates (50 cm3/h). The ΔP was measured by pressure transducer and shown by the data acquisition system at a specified time intervals. The fluid injection continued until the pressure difference was stabled.

  • By using recorded ΔP and measured mixture viscosity, the permeability ratios were continuously calculated by Darcy's law (Eq. (4)).

  • Produced fluids from outlet end of the sand pack during several PV of fluid injection were collected in samplers to measure asphaltene content.

  • Finally, the toluene was being injected continuously to remove any precipitated asphaltene formed due to n-heptane injection into the sand pack and to prepare the setup for second experiment.

The second dynamic experiment (DE-2) was carried out to study the effect of nanoparticles on the formation damage as following steps:

  • The porosity and absolute permeability were calculated by water injection as above procedure (k = 1194 md).

  • The crude oil/nanofluid solution was injected into the sand pack for finding the irreducible water saturation (swi = 32%).

  • The initial oil permeability was calculated by crude oil/nanofluid solution injection at constant flow rate and using Darcy's law (k0 = 294 md).

  • At the end, co-injection of n-heptane and crude oil/nanofluid solution into the sand pack was started at same injection rates (50 cm3/h) to record pressure differences. The permeability ratios after miscibility were continuously calculated with using recorded ΔP, flow rate of miscible fluid (100 cm3/h) and mixture viscosity by Darcy's law.

  • It should be noted that before the miscibility condition in both dynamic experiments, the fluid composition was unknown and a reliable value for the fluid viscosity could not be estimated [31,32]. The viscosities of different solutions including crude oil solution in DE-1, crude oil/nanofluid solution in DE-2 and mixture of n-heptane and oil solution after miscibility in both dynamic experiments were measured by Stabinger Anton paar- SVM 3000 viscometer at ambient condition. The measured viscosities are reported in Table 3. It can be observed that the viscosity of crude oil/nanofluid solution in DE-2 experiment is higher than the viscosity of crude oil without nanoparticles in DE-1 experiment.

2 Modeling of Asphaltene Deposition in the Sand Pack

In this study, theWang and Civan model [33] was modified with adding Langmuir adsorption term to model the asphaltene deposition in the case of using nanoparticles during miscible fluid injection into the sand pack. This new term represented the monolayer asphaltene adsorption onto nanoparticles. The proposed model included surface deposition, entrainment and pore throat plugging mechanisms to describe the asphaltene deposition. Several assumptions including one-dimensional flow, no gravity force and single liquid phase were considered for mathematical modeling.

The momentum balance (Eq. (4)) and mass balance equations in liquid phase for miscible fluid and asphaltene (Eqs. (5) and (6)) considered as follows [33]: t(ρlwfl)+x(ρlulwfl)=0(5) t(ρlwSAl+ρlwAl)+x(ρlulwSAl+ρlulwAl)=ρAEAnt(6) EAn=EANA(7) NA=109CnqLPρA(8) qLP=qmaxAnkLCEP1+kLCEP(9) CEP=106ρlCA(10) where EA is the volume fraction of the deposited asphaltene in the bulk volume of the sand pack, EAn is the volume fraction of the deposited asphaltene in the bulk volume after adsorption of asphaltene with nanoparticles, NA is the volume fraction of the adsorbed asphaltene with nanoparticles in the bulk volume of the porous media, Cn is the concentration of nanoparticles in crude oil, An is the surface area of nanoparticles, qLp is the adsorption amounts of asphaltene with zeolite beta nanoparticles in porous media that was calculated by Langmuir model, CEP is the mass fraction of the asphaltene content in the in the bulk volume of the sand pack, CA is the mass fraction of the asphaltene content in the liquid phase, qmax and kL are the Langmuir adsorption parameters that were obtained from batch adsorption experiments, wflis the mass fraction of the miscible fluid in the liquid phase, ρl and ρA are the liquid and asphaltene densities, respectively,  is the porosity, ulis the liquid velocity, wAland wSAl are the mass fractions of the dissolved and suspended asphaltene in the liquid phase, respectively.

The amount of dissolved asphaltene in the liquid phase was calculated by multiphase flash calculations and the Peng-Robinson equation of state was used for this purpose [34]. Subsequently, the asphaltene deposition model was used as follow [33]: EAt=αwSAlρlρAβEA(ννcr)+γulwSAlρlρA(11) where α, β and γ are the surface deposition rate, entrainment rate and pore throat plugging rate coefficients, respectively, and νcrare the interstitial and critical interstitial velocity of liquid phase, respectively. The value of β is described as follows: β=βi:υ>υcr(12) β=0:υ<υcr

The value of γ is set as follows: γ=γi(1+σEA):DptDptcr(13) γ=0:Dpt>Dptcr where σ is the snowball-effect deposition constant, Dpt and Dptcr are the pore throat diameter and critical pore throat diameter, respectively.

The porosity and permeability changes during asphaltene deposition process were calculated as follows:(14) k=k0(0)(15) where k is the permeability,  is the porosity, n is an adjustable constant and subscript 0 is refer to initial value.

Above equations were coupled and discretized by implicit finite difference method and solved by using an iteration calculation program. Numerical solution was carried out in a cartesian system with 100 grid blocks and 85 time steps with a Δt of 100 s. The genetic algorithm technique was used for optimization and determination of the model parameters (α, β, γi, νcr, Dptcr, σ, n). The square root of the summation of the differences between experimental and calculated permeability ratio and porosity ratio data was defined as an objective function: Fobj=w1i=1Ne((kk0)exp.i(kk0)Model,i)2+(1w1)i=1Ne((φφ0)exp.i(φφ0)Model,i)2(16) where Ne is the number of experimental data, and w1 is the weight of the permeability sentence in the objective function.

The initial and boundary conditions for numerical solution were considered as follows: EA=0;=0;k=k0:0<x<L,t=0(17) υ=QA:x=0,t>0(18) p=pambient::x= L,t>0(19) where Qis fluid flow rate (cm3/s), and A is cross-sectional area of the sand pack (cm2).

Table 3

Measured viscosities of the fluids in DE-1 and DE-2 experiments.

3 Results and discussion

3.1 Equilibrium Adsorption Experiment Results

For estimating the proper amount of nanoparticles for asphaltene adsorption from asphaltene-toluene solution, different concentrations of zeolite beta nanoparticles in solution were prepared and used. Figure 2 shows the adsorption efficiency versus the concentrations of zeolite beta nanoparticles. All adsorption experiments were replicated three times with the maximum uncertainty of ±5%.

It can be observed that with increasing the nanoparticle concentrations until 10 g/L, the adsorption efficiency was increased. This may be due to the fact that the higher amounts of zeolite beta nanoparticles provided more active surface for adsorption [35]. However, at the concentrations higher than 10 g/L nanoparticles in solution, the adsorption efficiency did not change too much. There are several works in the literature [11,14,19,26,36] concerning the asphaltene adsorption by using nanoparticles with the concentration of 10 (g/L). Therefore, the proper amount of zeolite beta nanoparticles was selected 10 g/L.

Figure 3 shows the adsorption amounts of asphaltene with the use of 10 g/L zeolite beta nanoparticles versus the equilibrium concentrations of asphaltene in the solutions. This figure demonstrates that the adsorption of asphaltene increased widely at lower asphaltene concentrations and reached to constant value by increasing the asphaltene concentrations. This characteristic confirmed the monolayer adsorption of asphaltene onto nanoparticles [25]. Therefore the Langmuir adsorption equation was used to model the equilibrium asphaltene adsorption data. As can be seen in Figure 3, the asphaltene adsorption data were described very well by Langmuir model (R2 > 0.99).

The Langmuir parameters (kL, qmax) were determined by linearizing. The kL value was found to be 0.0089 L/mg and the maximum equilibrium adsorbed asphaltene onto zeolite beta nanoparticles (qmax) was found to be 1.98 mg/m2. In comparison of the ability of zeolite beta nanoparticles in this study to another nanoparticles, the maximum amount of adsorbed asphaltene for the Langmuir adsorption on some nanoparticles (mg/m2) are reported in Table 4. It can be observed that the zeolite beta nanoparticles had a relatively good adsorption capacity for asphaltene adsorption from solutions. This result may be due to high surface area and acidic surface of zeolite beta nanoparticles [25].

thumbnail Figure 2

Adsorption efficiency results versus concentration of zeolite beta nanoparticles in solution.

thumbnail Figure 3

Adsorption amounts of asphaltene onto zeolite beta nanoparticles (the symbol is experimental data, and the solid line is the best fit of the Langmuir model).

Table 4

Maximum amount of adsorbed asphaltene by several nanoparticles.

3.2 Dynamic Experiment Results

3.2.1 Measurement of Pressure Differences

The aim of two dynamic experiments (DE-2 and DE-1) was to investigate the asphaltene deposition in a sand pack during the co-injection of n-heptane and crude oil solution with and without the use of zeolite beta nanoparticles, respectively. Figure 4 shows the recorded ΔP between two ends of the sand pack during fluid injection. It could be observed in both experiments, the ΔP decreased until near 1 PV of fluid injection. Co-injection of n-heptane and oil into the sand pack, which was previously saturated with the oil, caused a decrease of the fluid viscosity. Therefore, the ΔP values decreased until the miscibility occurred and a single phase was formed. An increase in ΔP values after 1 PV of fluid injection in the DE-1 indicated the formation damage. The n-heptane injection causes changing the equilibrium conditions and deposition of asphaltene in the sand pack. The increasing of ΔP values in DE-1 continued until 3 PV of fluid injection. However, in DE-2, the ΔP was relatively constant until 2 PV of fluid injection. After that slow increasing was seen in the ΔP values until 3.5 PV of fluid injection. Therefore, it can be resulted that the zeolite beta nanoparticles can delay the asphaltene deposition about 1 PV of fluid injection. It can be seen, ΔP values in DE-1 and DE-2 experiments were stabled at about 3 and 3.5 PV of fluid injection, respectively. Actually, the zeolite beta nanoparticles reduced the slope of increasing trend of ΔP due to asphaltene deposition in the sand pack.

It must be noted that the higher values of the ΔP in DE-2 with respect to DE-1 (Fig. 4) can be justified due to higher viscosity value of the crude oil/nanofluid solution and it is not related to more formation damage in the presence of zeolite.

thumbnail Figure 4

ΔP results between two ends of the sand pack.

3.2.2 Calculation of Permeability Ratio

During the dynamic experiments, permeability was calculated from the measured ΔP and mixture viscosity by Darcy's law. It can be seen in Figure 5 that the permeability ratio values decreased by increasing the PV of fluid injection. This confirms the asphaltene deposition occurring in n-heptane injection experiment. Before miscibility condition (<1 PV of fluid injection in this research), the permeability data could not be estimated because the viscosity of fluid was unknown.

Comparing the permeability curve in this research and literature [37] showed that in DE-1, the pore plugging was occurred immediately after miscibility (>1 PV of fluid injection). But in DE-2, most of the asphaltene adsorbed onto nanoparticles surfaces and the remaining asphaltene deposited onto the surface of the silica sands. By increasing the PV of fluid injection in the DE-2 (>2 PV of fluid injection), the asphaltene deposition were increased and the plugging of pores was gradually occurred in the sand pack according to limited capacity of nanoparticles for asphaltene adsorption. After 3 PV of fluid injection in DE-1 and after 3.5 PV of fluid injection in DE-2, the permeability ratios were not reduced. It is maybe due to that the deposited asphaltene decreased the cross-section of the flow in the sand pack. So the fluid velocity became higher than critical velocity. Therefore, the entrainment of deposited asphaltene in the sand pack occurred and consequently the permeability ratio was not decreased. As observed in Figure 5, the zeolite beta nanoparticles can increase the permeability ratio about 22% during n-heptane injection into the sand pack.

thumbnail Figure 5

Permeability ratio results during fluid injection experiment.

3.2.3 Determination of Outlet Fluid's Asphaltene Content and Porosity Ratio

The asphaltenes were extracted from the outlet fluid samples after specified PV of fluid injection based on the IP-143 standard analysis [21]. The asphaltene content (wt%) of the outlet fluid was calculated by dividing the weight of the extracted asphaltenes to the weight of the crude oil in the outlet sample. The asphaltene content of the outlet fluids after miscibility condition are showed in Figure 6.

It can be observed that the amounts of outlet asphaltene decreased during initial PV of fluid injection both in DE-1 and DE-2. But decreasing rate of outlet asphaltene in DE-1 was higher than in DE-2. The asphaltene content of outlet fluid without nanoparticles (DE-1) decreased from 1 wt% for initial concentration of asphaltene to about 0.54 wt%. But in the case of using nanoparticles (DE-2), the outlet asphaltene content decreased from 1 wt% to about 0.76 wt% resulting in an increase of the outlet asphaltene. In addition, it can be seen after 3 PV of fluid injection in DE-1 and after 3.5 PV of fluid injection in DE-2, the outlet asphaltene content was slightly increased because of occurring entrainment of deposited asphaltene. As observed in Figure 6, the zeolite beta nanoparticles increased the outlet fluid's asphaltene content about 40% comparing to DE-1 results. Because a part of asphaltene content adsorbed onto nanoparticles surfaces during fluid injection in the DE-2 experiment and it exited from the sand pack by flow. Therefore, the amount of deposited asphaltene (≈difference between asphaltene content in the outlet and inlet fluids) in the sand pack was decreased with the use of zeolite beta nanoparticles in the fluid injection experiment (Fig. 7).

In addition, the porosity ratio at specified pore volume of the fluid injection was calculated by dividing the differences between initial pore volume and deposited asphaltene volume to the bulk volume. It must be noticed that the volume of deposited asphaltene was calculated by dividing the mass differences of inlet and outlet asphaltene to the asphaltene density. The asphaltene density was equal to 1.0106 g/cm3. The porosity ratio results are shown in Figure 7. It can be observed the porosity ratio reduced to about 0.985 and 0.991 in DE-1 and DE-2, respectively. The less porosity reduction comparing to permeability reduction during fluid injection is indicated to asphaltene deposition in the local positions instead along the sand pack [38].

thumbnail Figure 6

Outlet fluid's asphaltene content (wt%).

thumbnail Figure 7

Porosity ratio and asphaltene deposition results during fluid injection experiments.

3.3 Dynamic Modeling Results

A model based on the adding adsorption term to the Wang and Civan deposition model [33] was presented for calculating of permeability and porosity ratios in the case of using nanoparticles in porous media (DE-2). The calculated porosity and permeability ratios by using proposed model and Wang and Civan model are shown in Figures 8 and 9, respectively.

The modeling results in Figures 8 and 9 confirmed that the nanoparticles can improve the porosity and permeability ratios during fluid injection test. It can be observed that the results of permeability and porosity ratios by the proposed model were in good agreement with experimental data. The average absolute deviation (AAD%) of the calculated permeability and porosity ratios from the experimental results with using both models were determined and the adjusted parameters of the models that obtained from genetic algorithm technique are reported in Table 5. It was found that the proposed model is capable of predicting the permeability ratios in DE-2 with AAD% of 1.07. Whereas permeability ratio prediction based on the Wang and Civan model is less accurate with AAD% of 2.25%. In addition, the proposed model can predict the porosity ratios in DE-2 with high accuracy (AAD% = 0.07). But the Wang and Civan model predicted porosity ratios with AAD% of 0.94. So, it can be concluded that the proposed model in study, can use as a suitable model for modeling the asphaltene deposition during fluid injection with the use of nanoparticles in the porous media.

thumbnail Figure 8

Comparison of the porosity ratios experimental data and predictions of the proposed model and Wang and Civan model.

thumbnail Figure 9

Comparison of the permeability ratios experimental data and predictions of the proposed model and Wang and Civan model.

Table 5

Adjusted parameters of asphaltene deposition models in DE-2 experiment.

Conclusions

At first, a series of experiments including equilibrium asphaltene adsorption onto zeolite beta nanoparticles was conducted by UV-Vis spectrophotometer to determine the adsorption capacity of nanoparticles. The adsorption experimental results showed that the zeolite beta can efficiently adsorb asphaltene from solution. It was found that, an adsorption isotherm good matched by the Langmuir model, suggesting a monolayer asphaltene adsorption onto zeolite. Special properties of zeolite beta nanoparticles such as high surface area to volume ratio and high functionalizable surface cause good adsorption capacity for the zeolite beta. The results of adsorption efficiency showed that the proper concentration of nanoparticles for asphaltene adsorption was equal to 10 g/L zeolite beta in the oil. However, in order to determine the economically optimum amount of nanoparticles in the industry, the economic aspect should be evaluated. In addition, two dynamic experiments including co-injection of crude oil and n-heptane into a silica sand pack were conducted with and without the use of zeolite beta nanoparticles. The purpose of dynamic experiments was to investigate the effect of the zeolite beta nanoparticles on reduction of formation damage during oil injection in porous media. The dynamic experiments showed that zeolite beta could delay the increasing of the ΔP and reduce the asphaltene deposition. It must be noted that the use of nanoparticles increased the fluid's viscosity and as a result, it increased the total values of the ΔP compared to the experiment without the use of zeolite beta. Moreover, the results showed that the zeolite beta nanoparticles increased the permeability and porosity ratios and reduced asphaltene deposition during co-injection of crude oil and n-heptane into the silica sand pack. The interactions between the asphaltene particles with their basic nature and acidic surfaces of zeolite beta nanoparticles caused well adsorption of asphaltene and therefore improvement of the formation damage. It can be seen from the modeling results that the proposed model based on the monolayer adsorption of asphaltene onto nanoparticles and three mechanisms of asphaltene deposition (surface deposition, entrainment and pore throat plugging) is more accurate in the prediction of permeability and porosity ratios (AAD% of 1.07 and 0.07) than the Wang and Civan model (AAD% of 2.25 and 0.94). This study confirms that the zeolite beta nanoparticles could be used for reducing asphaltene deposition in porous media and finally enhancing oil recovery.

Nomenclature

AAD: average absolute deviation

A: cross-sectional area of the sand pack

An: surface area of nanoparticles

CA: mass fraction of the asphaltene content in the liquid phase

C0: initial concentration of asphaltene in solutions

Ce: equilibrium concentration of asphaltene in solutions

Cn: concentration of nanoparticles in crude oil

CEP: mass fraction of the asphaltene content in the in the bulk volume of the sand pack

DE-1: dynamic experiment 1

DE-2: dynamic experiment 2

Dpt: pore throat diameter

Dptcr: critical pore throat diameter

E: adsorption efficiency

EA: volume fraction of the deposited asphaltene in the bulk volume of the sand pack

EAn: volume fraction of the deposited asphaltene in the bulk volume after adsorption of asphaltene on nanoparticles

Fobj: objective function

k: permeability

k0: initial permeability

kL: Langmuir equilibrium adsorption constant (L/mg)

L: length of the sand pack

MW: molecular weight

m: dry mass of nano-particles

n: adjustable parameter

Ne: number of experimental data

NA: volume fraction of the adsorbed asphaltene with nanoparticles in the bulk volume of the porous media

P: pressure

ΔP: pressure differences

PV: pore volume

Q: fluid flow rate

qe: equilibrium adsorbed amount of asphaltene

qmax: maximum amount of adsorbed asphaltene for the Langmuir adsorption on adsorbent surfaces

qLP: Langmuir Adsorption amounts of asphaltene with nanoparticles in porous media

R2: correlation coefficient

SG: specific gravity

swi: irreducible water saturation

t: time

ul: liquid velocity

V: volume of solution

w1: weight of the permeability sentence in objective function

wAl: mass fractions of the dissolved asphaltene in the liquid phase

wSAl: mass fractions of the suspended asphaltene in the liquid phase

wfl: mass fraction of the miscible fluid in the liquid phase

wt%: weight fraction percent

Greek symbols

α: surface deposition rate coefficient

β: entrainment rate coefficient

γ: pore throat plugging rate coefficient

δ: snowball-effect deposition constant

μ: fluid viscosity

ρ: density

ρl: liquid density

ρA: asphaltene density

φ: porosity

φ0: Initial porosity

ν: interstitial velocity of liquid phase

νcr: critical interstitial velocity of liquid phase

References

  • Jafari Behbahani T., Ghotbi C., Taghikhani V., Shahrabadi A. (2015) Experimental study and mathematical modeling of asphaltene deposition mechanism in core samples, Oil Gas Sci. Technol. − Rev. IFP 70, 6, 1051–1074. [CrossRef] [Google Scholar]
  • Buckley J.S. (1998) Wetting alteration of solid surfaces by crude oils and their asphaltenes, Oil Gas Sci. Technol. − Rev. IFP 53, 3, 303–312. [Google Scholar]
  • Karambeigi M.A., Kharrat R. (2014) An investigation of inhibitors performance on asphaltene precipitation due to associated gas injection, Petrol. Sci. Technol. 32, 10, 1213–1218. [CrossRef] [Google Scholar]
  • Kashefi S., Shahrabadi A., Jahangiri S., Lotfollahi M.N., Bagherzadeh, H. (2016) Investigation of the performance of several chemical additives on inhibition of asphaltene precipitation, Energy Sour. Part A Recovery Util. Environ. Eff. 38, 24, 3647–3652. [CrossRef] [Google Scholar]
  • Kashefi S., Shahrabadi A., Lotfollahi M.N., Varamesh, A. (2016) A new polymeric additive as asphaltene deposition inhibitor in CO2 core flooding, Korean J. Chem. Eng. 33, 11, 3273–3280. [CrossRef] [Google Scholar]
  • Marczewski A.W., Szymula M. (2002) Adsorption of asphaltenes from toluene on mineral surface, Colloids Surf. A: Physicochem. Eng. Asp.208, 1, 259–266. [CrossRef] [Google Scholar]
  • Pernyeszi T., Patzko A., Berkesi O., Dékány I. (1998) Asphaltene adsorption on clays and crude oil reservoir rocks, Colloids Surf. A: Physicochem. Eng. Asp. 137, 1, 373–384. [Google Scholar]
  • Dudášová D., Simon S., Hemmingsen P.V., Sjöblom, J. (2008) Study of asphaltenes adsorption onto different minerals and clays: Part 1. Experimental adsorption with UV depletion detection, Colloids Surf. A: Physicochem. Eng. Aspects 317, 1, 1–9. [Google Scholar]
  • Ekholm P., Blomberg E., Claesson P., Auflem I.H., Sjöblom J., Kornfeldt A. (2002) A quartz crystal microbalance study of the adsorption of asphaltenes and resins onto a hydrophilic surface, J. Colloid Interface Sci. 247, 2, 342–350. [CrossRef] [PubMed] [Google Scholar]
  • Abdallah W.A., Taylor S.D. (2007) Surface characterization of adsorbed asphaltene on a stainless steel surface, Nucl. Instrum. Meth. B 258, 1, 213–217. [CrossRef] [Google Scholar]
  • Nassar N.N., Hassan A., Pereira-Almao P. (2011) Metal oxide nanoparticles for asphaltene adsorption and oxidation, Energy Fuels 25, 3, 1017–1023. [CrossRef] [Google Scholar]
  • Nassar N.N., Hassan A., Pereira-Almao P. (2011) Effect of surface acidity and basicity of aluminas on asphaltene adsorption and oxidation, J. Colloid Interface Sci. 360, 1, 233–238. [CrossRef] [PubMed] [Google Scholar]
  • Nassar N.N., Hassan A., Carbognani L., Lopez-Linares F., Pereira-Almao P. (2012) Iron oxide nanoparticles for rapid adsorption and enhanced catalytic oxidation of thermally cracked asphaltenes, Fuel 95, 257–262. [CrossRef] [Google Scholar]
  • Nassar N.N. (2010) Asphaltene adsorption onto alumina nanoparticles: kinetics and thermodynamic studies, Energy Fuels 24, 8, 4116–4122. [CrossRef] [Google Scholar]
  • Hashemi R., Nassar N.N., Pereira-Almao P. (2012) Transport behavior of multimetallic ultradispersed nanoparticles in an oil-sands-packed bed column at a high temperature and pressure, Energy Fuels 26, 3, 1645–1655. [CrossRef] [Google Scholar]
  • Mirzayi B., Shayan N.N. (2014) Adsorption kinetics and catalytic oxidation of asphaltene on synthesized maghemite nanoparticles, J. Petrol. Sci. Eng. 121, 134–141. [CrossRef] [Google Scholar]
  • Kazemzadeh Y., Eshraghi S.E., Kazemi K., Sourani S., Mehrabi M., Ahmadi Y. (2015) Behavior of asphaltene adsorption onto the metal oxide nanoparticles surface and its effect on the heavy oil recovery, Ind. Eng. Chem. Res. 54, 1, 233–239. [Google Scholar]
  • Hashemi S.I., Fazelabdolabadi B., Moradi S., Rashidi A.M., Shahrabadi A., Bagherzadeh H. (2016) On the application of NiO nanoparticles to mitigate in situ asphaltene deposition in carbonate porous matrix, Appl. Nanosci. 6, 1, 71–81. [CrossRef] [Google Scholar]
  • Franco C., Patiño E., Benjumea P., Ruiz M.A., Cortés F.B. (2013) Kinetic and thermodynamic equilibrium of asphaltenes sorption onto nanoparticles of nickel oxide supported on nanoparticulated alumina, Fuel 105, 408–414. [CrossRef] [Google Scholar]
  • Franco C.A., Nassar N.N., Ruiz M.A., Pereira-Almao P., Cortés F.B. (2013) Nanoparticles for inhibition of asphaltenes damage: adsorption study and displacement test on porous media, Energy Fuels 27, 6, 2899–2907. [CrossRef] [Google Scholar]
  • ASTM D6560. (2005) Standard test method for determination of asphaltenes (heptane insolubles) in crude petroleum and petroleum products. Annual Book of ASTM Standards, American Society for Testing and Materials, Philadelphia, PA, 5(3). [Google Scholar]
  • Fan T., Buckley J.S. (2006) Acid number measurements revisited, SPE/DOE Symposium on Improved Oil Recovery, Tulsa, Oklahoma, 22–26 April, SPE Paper 99884-MS. [Google Scholar]
  • Saadatjoo N., Golshekan M., Shariati S., Kefayati H., Azizi P. (2013) Organic/inorganic MCM-41 magnetite nanocomposite as a solid acid catalyst for synthesis of benzo [α] xanthenone derivatives, J. Mol. Catal. A: Chem. 377, 173–179. [CrossRef] [Google Scholar]
  • Yudin I.K., Nikolaenko G.L, Gorodetskii E.E., Kosov V.I., Melikyan V.R., Markhashov E.L., Frot D., Briolant Y. (1998) Mechanisms of asphaltene aggregation in toluene-heptane mixtures, J. Petrol. Sci. Eng. 20, 3, 297–301. [CrossRef] [Google Scholar]
  • Hosseinpour N., Khodadadi A.A., Bahramian A., Mortazavi Y. (2013) Asphaltene adsorption onto acidic/basic metal oxide nanoparticles toward in situ upgrading of reservoir oils by nanotechnology, Langmuir29, 46, 14135–14146. [CrossRef] [PubMed] [Google Scholar]
  • Shayan N.N., Mirzayi B. (2015) Adsorption and removal of asphaltene using synthesized maghemite and hematite nanoparticles, Energy Fuels 29, 3, 1397–1406. [CrossRef] [Google Scholar]
  • Langmuir I. (1916) The constitution and fundamental properties of solids and liquids. Part I. Solids, J. Am. Chem. Soc. 38, 11, 2221–2295. [Google Scholar]
  • Mohammadi M., Akbari M., Fakhroueian Z., Bahramian A., Azin R., Arya S. (2011) Inhibition of asphaltene precipitation by TiO2, SiO2, and ZrO2 nanofluids, Energy Fuels 25, 7, 3150–3156. [CrossRef] [Google Scholar]
  • Haddad Z., Abid C., Oztop H.F., Mataoui A. (2014) A review on how the researchers prepare their nanofluids, Int. J. Therm. Sci. 76, 168–189. [CrossRef] [Google Scholar]
  • Jafari Behbahani T., Ghotbi C., Taghikhani V., Shahrabadi A. (2012) Investigation on asphaltene deposition mechanisms during CO2 flooding processes in porous media: a novel experimental study and a modified model based on multilayer theory for asphaltene adsorption, Energy Fuels 26, 8, 5080–5091. [CrossRef] [Google Scholar]
  • Papadimitriou N.I., Romanos G.E., Charalambopoulou G. Ch., Kainourgiakis M.E., Katsaros F.K., Stubos A.K. (2007) Experimental investigation of asphaltene deposition mechanism during oil flow in core samples, J. Petrol. Sci. Eng. 57, 3, 281–293. [CrossRef] [Google Scholar]
  • Bagherzadeh H., Ghazanfari M.H., Kharrat R., Rashtchian D. (2014) Experimental investigation and modeling of permeability impairment mechanisms due to asphaltene precipitation under CO2 injection conditions, Energy Sour. Part A Recovery Util. Environ. Eff. 36, 6, 591–604. [CrossRef] [Google Scholar]
  • Wang S., Civan F. (2001) Productivity decline of vertical and horizontal wells by asphaltene deposition in petroleum reservoirs, SPE International Symposium on Oilfield Chemistry, Houston, Texas, 13–16 Feb., SPE paper 64991-MS. [Google Scholar]
  • Kashefi S., Lotfollahi M.N., Shahrabadi A. (2017) An investigation of asphaltene deposition mechanisms during natural depletion process by a two phase modeling using genetic algorithm technique, J. Petrol. Sci. Technol. 7, 2, 12–20. [Google Scholar]
  • Selvin R., Roselin L.S., Khayyat S.A., Umar A. (2013) Hierarchical zeolite beta: An efficient and eco-Friendly nanocatalyst for the friedel-crafts acylation of toluene, J. Nanosci. Nanotechnol. 13, 6, 4415–4420. [CrossRef] [PubMed] [Google Scholar]
  • Nassar N.N., Hassan A., Luna G., Pereira-Almao P. (2013) Kinetics of the catalytic thermo-oxidation of asphaltenes at isothermal conditions on different metal oxide nanoparticle surfaces, Catal. Today207, 127–132. [CrossRef] [Google Scholar]
  • Roque C., Chauveteau G., Renard M., Thibault G., Bouteca M., Rochon J. (1995) Mechanisms of formation damage by retention of particles suspended in injection water, SPE European Formation Damage Conference, Hague, Netherlands, 15–16 May, SPE paper 30110-MS. [Google Scholar]
  • Bagheri M., Kharrat R., Ghotbi C. (2011) Experimental investigation of the asphaltene deposition process during different production schemes, Oil Gas Sci. Technol. − Rev. IFP 66, 3, 507–519. [Google Scholar]

All Tables

Table 1

Specifications of the crude oil.

Table 2

Properties of the zeolite beta nanoparticles.

Table 3

Measured viscosities of the fluids in DE-1 and DE-2 experiments.

Table 4

Maximum amount of adsorbed asphaltene by several nanoparticles.

Table 5

Adjusted parameters of asphaltene deposition models in DE-2 experiment.

All Figures

thumbnail Figure 1

Schematic view of fluid injection setup.

In the text
thumbnail Figure 2

Adsorption efficiency results versus concentration of zeolite beta nanoparticles in solution.

In the text
thumbnail Figure 3

Adsorption amounts of asphaltene onto zeolite beta nanoparticles (the symbol is experimental data, and the solid line is the best fit of the Langmuir model).

In the text
thumbnail Figure 4

ΔP results between two ends of the sand pack.

In the text
thumbnail Figure 5

Permeability ratio results during fluid injection experiment.

In the text
thumbnail Figure 6

Outlet fluid's asphaltene content (wt%).

In the text
thumbnail Figure 7

Porosity ratio and asphaltene deposition results during fluid injection experiments.

In the text
thumbnail Figure 8

Comparison of the porosity ratios experimental data and predictions of the proposed model and Wang and Civan model.

In the text
thumbnail Figure 9

Comparison of the permeability ratios experimental data and predictions of the proposed model and Wang and Civan model.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.