INTRODUCTION

Asphaltenes are the heaviest molecular weight components that are found in crude oil, along with saturates, aromatic and resins, which are soluble in toluene and insoluble in n-heptane [13]. The presence of asphaltenes in oil poses serious problems during oil recovery which is a major concern in all the world’s oil fields [4, 5]. One of these popular oil fields that encounter asphaltene deposition is Hassi Messaoud which is located in southern Algeria [1]. Indeed, as a result of asphaltenes deposition and precipitation, many oil wells have been entirely shut down in this oil field [6]. Currently, the remediation of asphaltenes includes several methods such as mechanical scratch treatment, tube coating, hot fluid removal, and solvent cleaning [7]. However, these classical techniques remain inefficient in asphaltenes removal. Recently, nanotechnology using novel nanoparticles has been considered to be an emerging technique with potential applications in the oil industry [8, 9]. Metal oxide nanoparticles like iron oxides, for instance, are one of particular interest because they tend to associate with asphaltene molecular aggregates and can behave as efficient nanoadsorbents [10, 11]. A recent study on Hamaca asphaltenes adsorption showed that iron oxide had an uptake of 3.5 to 4 mg asphaltenes/m2 at a temperature near to room temperature using a high dosage of adsorbent (50 to 30 g/L) [10]. Impregnation of Fe2O3 on the surface of kaolinite, carbon nanotubes (CNT), montmorillonite and SiO2 strongly improved the adsorption capacity relative to the surfaces of unmodified sorbent bases [12, 13]. Nassar et al. for instance, used six different metal oxide nanoparticles for the adsorption and oxidation of asphaltenes, namely Fe3O4, Co3O4, TiO2, MgO, CaO, and NiO [14, 15]. The results showed that the adsorption of asphaltenes was affected by the types of metal oxides and the asphaltenes uptake followed the order CaO > Co3O4 > Fe3O4 > MgO > NiO > TiO. Additionally, iron oxide nanoparticles have shown also to be an excellent adsorbent and a catalyst for catalytic oxidation of thermally cracked asphaltenes [14]. In another study, magnetic iron oxide nanoparticles capped with sulfonated asphaltenes have been proposed as natural and inexpensive materials that could act as oil collectors for Arabian heavy crude oil [16]. In that study, three types of metal oxide nanoparticles were considered for asphaltenes oxidation, namely: NiO, Co3O4, and Fe3O4, and the results showed that the presence of nanoparticles decreases the activation energy of asphaltenes oxidation and increases the rate of reaction [17]. The impact of iron oxide nanoparticles on the precipitation of asphaltenes during CO2 injection was studied by Kazemzadeh et al. [18]. The synthesized maghemite (γ-Fe2O3) and hematite (α-Fe2O3) nanoparticles were used for the adsorption and elimination of asphaltenes. The reported results showed that the maximum adsorption capacity of maghemite and hematite was 108.1 and 45.8 mg/g, respectively. Maghemite and hematite could be used as nanoadsorbents of asphaltenes, where the maghemite is more effective towards their adsorption [19]. Moreover, Franco et al. studied the effect of resins on the adsorption behavior of n-C7 asphaltenes on silica and hematite nanoparticles [20]. The results indicated that resin had no significant influence on the adsorbed amount of asphaltenes, which showed that resin I has a solvent-like behavior, such as toluene, mainly at low concentrations (<3000 mg/L) [20].

However, despite these exertions attempting to eliminate or prohibit asphaltenes precipitation, these challenges remain open and unsolved, especially for Algerian asphaltenes. This because the asphaltenes precipitation is among the most important challenges in Algerian oil fields which significantly affect oil recovery. It should be noted here that Algeria is among the world giants in the hydrocarbons sector, especially in Hassi Messaoud oil field. This oil field is the largest oil deposits in Algeria discovered in 1956 (2500 km2) and contains about 71% of proved reserves in Algeria [21]. In fact, with this field; Algeria is classified as the third in Africa continent and the 16th amongst more than 55 countries in the world regarding oil resources [22]. It produces up to 65.3 million tonnes of oil and 12.2 billion barrels of proved oil reserves for 2018 [22]. Based on the literature survey it seems clear that no work has been conducted to treat the asphaltenes problems using adsorption with nanoparticles in Hassi Messaoud oil field. Thus, accordingly, at the level of the author’s knowledge, it is the first time that such a study focuses on the adsorption of the asphaltenes contained in the Algerian crude oil using such novel synthesized maghemite iron oxide (γ-Fe2O3) nanoparticles. Herein, in the present study, the adsorption of Algerian asphaltenes from Hassi Messaoud onto novel improved maghemite iron oxide (γ-Fe2O3) nanoparticles was investigated. The nanoparticles were prepared and routinely characterized using HRTEM, BET and XRD techniques. The nanoparticles used in this study are proposed as potential candidates for asphaltenes elimination from Algerian crude oils. This approach of asphaltenes uptake by nanoparticles shall extensively improve the oil recovery and minimize the problematic issues caused by asphaltenes aggregation and precipitation.

EXPERIMENTAL

Materials. An oil sample was obtained from Hassi Messaoud oil field and employed to extract the asphaltenes. n-heptane (99%) was used as a solvent for asphaltenes extraction from Algerian light oil and toluene HPLC grade was used as a solvent to solubilize the asphaltene sample. Both solvents were obtained from Sigma-Aldrich, Ontario, Canada and used as received. Table S1 in the Supplementary Material provides a full saturates, aromatics, resins, asphaltenes (SARA) analysis of the oil. For γ-Fe2O3 nanoparticles preparation, the following chemicals and reagents were used, ammonium iron(II) sulfate hexahydrate (NH4)2Fe(SO4)2·6H2O, ammonium hydroxide solution (26%) NH4OH, hydrogen peroxide (H2O2, 30 wt %), and nitric acid (HNO3). All of these materials were received from Sigma-Aldric, Ontario, Canada and used without further modifications.

Extraction of asphaltenes. Asphaltenes were extracted from the Algerian Hassi Messaoud oil field following the standard method ASTM D6560 (i.e., 1 : 40 oil/solvent ratio) [23]. The extraction process was carried out by adding n-heptane to the oil well sample with an oil/solvent ratio (O/S) of 1 : 40 (g/mL). The resulting mixture was heated to 90°C in a flask. After that, the mixture was sonicated for 45 min at room temperature for stabilization. The solution was then filtered using a 0.45 µm pore diameter filter paper [24]. The precipitated asphaltenes were washed with n-heptane at a ratio of 1 : 5 (g/mL) (O/S) in a Soxhlet apparatus until the filtrate became colorless. Finally, the filtered asphaltenes were dried overnight under vacuum at 80°C.

Preparation of maghemite iron oxide nanoparticles (γ-Fe2O3). The γ-Fe2O3 nanoparticles were in-housed prepared using the precipitation method at room temperature [25]. The preparation method was firstly started by dissolving 6 g (NH4)2Fe(SO4)2·6H2O in 100 mL of deionized water (DIW). After complete dissolution of the iron salt (~5 min), 5 mL of NH4OH solution (26%) was added dropwise at a rate of 1 cm3/s using a syringe pump. Subsequently, 15 drops of H2O2 (30 wt %) were added to the solution to partially oxidize part of the Fe2+ to Fe3+ and thus producing a precipitated black nanocrystalline of Fe3O4. After that, vacuum filtration and washing were performed using 300−400 mL DIW at room temperature. The produced black materialwas allowed to dry under vacuum overnight. The next day, the dried Fe3O4 nanomaterial was calcined in an oven at 250°С for 2 h. This step is important to release some hydroxyl ions (OH) from the surface and to oxidize the Fe2+ remaining in the magnetite to produce the desired maghemite. Then, the obtained sample was immersed in the acidic medium (HNO3) to remove any excess Fe3+ produced during the transformation of magnetite into maghemite. A suspension of 3 g/L of nano-γ-Fe2O3 in DIW was prepared, and the pH of this suspension was adjusted to 3 using 0.10 M HNO3 acid. The solution was kept for at least 5 min at pH of 3. Then, the temperature was increased up to 90°C and kept for 30 min under magnetic stirring (300–400 rpm). Finally, the sample was filtered, washed with DIW and dried to produce Fe2O3 adsorbent nanoparticles [26].

Characterization of maghemite nanoparticles (γ-Fe2O3). The crystalline domain structures and sizes of the as-prepared γ-Fe2O3 nanoparticles were determined using an Ultima III multi-purpose diffraction system (Rigaku Corp., The Woodlands, TX) with CuKα radiation operating at 40 kV and 44 DmA and scan range of 3°–90° 2θ using a 0.05° degree and a counting time of 0.2°/min. The average crystalline domain sizes were measured using the Debye–Scherrer equation [Eq. (1)] [19] implemented in JADE commercial software (supplied with the diffractometer) by calculating the full width at half height, full width at half-maximum (FWHM) peaks corresponding to the experimental profile to a pseudo-Voigt profile function.

$$D = {{K\lambda } \over {\beta \cos \theta }},$$
((1))

where K = 0.9 is Scherrer’s constant, λ = 0.154 nm is the wavelength of X-rays, θ is the diffraction angle in degrees, and β is the width of the peak at the half-maximal height of the peak.

The surface area and pore size distributions were measured using Brunauer–Emmett–Teller (BET) method by performing nitrogen physisorption at –196°C using a porosity analyzer and TriStar II 3020 specific surface area (Micrometrics Instrument Corporation, Norcross, GA). The sample was firstly degassed at 150°C under N2 flow overnight to remove any existing moisture. Particle size was calculated with the BET surface area according to Eq. (2):

$$d = {{6000} \over {(SA \cdot {\rho _{{\rm{mag}}}})}},$$
((2))

where d is the particle size in nm, SA is the experimentally measured surface area (m2/g), and ρmag is the density of maghemite (4.9 g/cm3).

The surface morphology of the nanoparticles was also analyzed by FEI Tecnai F20 FEG Transmission electron microscopy (TEM) using a 200 kV acceleration voltage. About 0.5 mg of nanoparticles were dispersed in 1 mL of pure ethanol, and then a drop of the solution was deposited on the sample holder. After evaporating the ethanol, the dried powder, on the grid support, is ready for imaging.

Adsorption experiments. Adsorption kinetic experiments. Adsorption kinetic experiments were carried out with a UV-Vis spectrophotometer (Nicolet Evolution 100, Thermo Instruments Canada, Inc., Mississauga, Ontario, Canada). Three initial concentrations of asphaltenes of 100, 500, and 1000 mg/L were prepared, and the temperature of the solution was 25°C. 10 mL of asphaltenes-toluene solution at a ratio of 1 : 10 (L/g) asphaltenes mass solution/nanoparticle was used. The sample was shaken using a Wrist Action (Burrel, Model 75-BB) shaker, and thus the change of absorbance for asphaltenes in the supernatant was recorded continuously up to 4 h. The kinetic curves were plotted as a function of time of 1–150 min, especially the amount adsorbed increased rapidly during the first 15 min and remained unchanged after 60 min of contact.

Equilibrium adsorption isotherms. The batch adsorption experiments were carried out at a temperature of 25°C by preparing different solutions of asphaltenes, dissolved in toluene. All adsorption experiments were conducted in 10 mL of asphaltenes-toluene solution at a ratio of 1 : 10 (L/g) asphaltenes mass solution/ nanoparticle. The tested vials were sealed firmly to prevent toluene evaporation. The adsorption was allowed to proceed by shaking the flasks for 24 h on a Wrist Action (Burrel, Model 75-BB) shaker, ensuring equilibrium was achieved. Then, the nanoparticles containing asphaltenes were separated by centrifugation at 5000 rpm for 30 min. The concentrations of asphaltenes in the solution were evaluated using a UV-Vis spectrophotometer (Nicolet Evolution 100, Thermo Instruments Canada, Inc., Mississauga, Ontario, Canada). The calibration curve was constructed at a wavelength of 298 nm. The selected experiments were repeated twice to confirm their reproducibility. The adsorption of asphaltenes onto iron oxide nanoparticles was determined from the change in the concentration of asphaltenes in the model oil solution, before and after mixing with the nanoparticles. The adsorbed amount of n-C7 asphaltenes (mg/m2) was calculated based on the mass balance analysis, as indicated in Eq. (3):

$$Q = {{({c_0} - {c_{\rm{e}}})V} \over {SA \cdot m}},$$
((3))

where c0 is the initial concentration of n-C7 asphaltenes in solution (mg/L), ce is the equilibrium concentration of n-C7 asphaltenes in the supernatant (mg/L), V is the volume of solution (L), SA is the BET surface area of the γ-Fe2O3 prepared nanoparticles (m2/g), and m is the mass of the nanoparticles (g).

RESULTS AND DISCUSSIONS

Characterization of γ-Fe2O3 nanoparticles. The crystalline structure, particle size and specific surface area of the prepared nanoparticles (γ-Fe2O3) were measured and summarized in Table 1. As shown in the table, after measuring the specific surface area using BET equation (120 m2/g) and the external surface using the t-plot method (114 m2/g), there is no significant difference observed between values, thus suggesting that the prepared iron oxide nanoparticles are non or macro-porous and maintaining a high external surface area [27]. Nitrogen adsorption-desorption isotherms for iron oxide nanoparticles, provided in the Supplementary Material (Fig. S1), follows the isotherms type II according to the IUPAC classification, for macropore solids [28]. To get more information about the accessible surface area, one the most important parameters that plays an important role on the reactivity of the materials for the adsorption of asphaltenes, the estimated pore diameters using the Barrett–Joynes–Halenda (BJH) method ranged from 20 Å to more than one hundred nanometers with a sharp peak at 82 Å, as shown in (Fig. S2) in the Supplementary Material. At the pore size distribution below 50 Å, the micro and mesoporosity can barely observed and their contribution is very small. The surface area microporous area (t-plot Micropore Area) was around 6 m2/g. On the other hand, a significant contribution to pore size distribution can be observed above 50 Å, indicating for the macro-porosity of the material. Accordingly, the as-prepared nanoparticle has great potential to accommodate large molecules such as asphaltenes. In addition, for purposes of comparison with the measured XRD particle size, the particle size (assuming spherical particles) was also calculated based on the surface area measurements, as previously explained in Eq. (2). Thus, the obtained value for d agrees well with the one obtained by XRD, indicating that the material may indeed be composed of spherical-type nanoparticles. Figure 1 shows the X-ray diffractogram of the γ-Fe2O3 nanoparticles. The identification of the pattern confirms that the material is maghemite, as reported by the preparation protocol. The structure was identified by comparing the XRD signals with those reported in the pdf map 01-076-4113 ICDD 2005 (International Center for Diffraction Data) database included in the JADE V.7.5.1 software. As seen from Fig. 1, the peaks are very broad indicating very small crystalline domain sizes for the used γ-Fe2O3. The presence of the tetragonal vacancy-ordered superstructure in these nanoparticles is clearly marked by the appearance of weak XRD signals a 2θ CuKα of 15°, 24°, and 26°.

Table 1. Textural properties and average particle size of γ-Fe2O3 nanoparticles
Fig. 1.
figure 1

X-ray diffraction pattern of synthesized γ-Fe2O3 nanoparticles (blue lines). The vertical lines (black color) are the reference COD database of Maghemite-Q-Fe2O3 (01-076-4113 ICDD 2005).

To gain a better understanding of the morphology and size of γ-Fe2O3 nanoparticles used in this study, HRTEM images are depicted in (Fig. 2), which confirms the presence of nanometric particles of γ-Fe2O3. Most of the nanoparticles are spherical with diameters range between 4 and 11 nm, while some are in the form of short clusters. Large aggregates of particles allowed observation of morphology types and measurement of individual sizes when autonomous nanoparticles were expelled by the electron beam when an attempt was made to analyze them.

Fig. 2.
figure 2

HRTEM images for the γ-Fe2O3 studied nanoparticles (the line mark in the images corresponds to 20 and 50 nm).

Adsorption of asphaltenes onto γ-Fe2O3 nanoparticles. Kinetics of adsorption. Figure 3 shows the change in the amount of asphaltenes adsorbed as a function of contact time. In all cases (100, 500, and 1000 ppm), the amount of asphaltenes adsorbed was rapidly increased in the first 15 min and remained unchanged after 60 min of contact. This indicates that steady-state was almost reached within 15 min. This rapid adsorption of asphaltenes is certainly linked to the high degree of dispersion of the nanoparticles as well as to the great availability of the external surface because the γ-Fe2O3 selected is a non-porous adsorbent [14, 15]. Therefore, it takes a short time to adsorb the molecules, because of the lack of intra-particle diffusion which limits the adsorption rate. The equilibrium results are in agreement with what has been reported by Dr. Nassar group for the adsorption of asphaltenes on iron oxide nanoparticles [14]. However, it is different than the other conventional porous adsorbents, in which very slow adsorption of asphaltenes occurs through the pore diffusion steps, which in turn requires a longer equilibrium time [29, 30]. It is noteworthy that the adsorption of asphaltenes onto surfaces of nanoparticles depends on the type and strength of interactions between the asphaltene species and the solid surface [9]. It was reported that there are a number of interparticle forces are responsible for successful interaction between nanoparticles and asphaltenes functionality groups such as carboxylic, pyrrolic, pyridininc, thiophenic, and sulfite groups [31, 32]. The major forces that could contribute to these interactions are van der Waals, electrostatic, charge transfer, hydrogen bonding, and steric interactions [32]. The polar interactions, for example, are acid-base reactions occurring between adsorbate and adsorbent. However, in the case of iron oxide as a type of amphoteric oxides demonstrate higher adsorption capacity as compared with acidic oxides and similar to that of basic oxides. Similar observations were reported elsewhere for adsorption of asphaltenes onto different minerals and clays [3336]. Mohammadi et al. [37], for instance, used different types of metal oxide nanoparticles, such as SiO2, TiO2, and ZrO2, for enhancing the stability of asphaltene nanoaggregates by forming hydrogen bonds at acidic conditions. The authors confirmed that TiO2 nanofluids, under strongly acidic conditions, might be acting as a dispersant, thereby enhancing the stability of the asphaltenes and leading to a higher precipitation onset point. Similarly, Hosseinpour et al. [3] reported the adsorption of asphaltene onto different acid/base metal oxides nanoparticles, such as NiO, Fe2O3, WO3, MgO, CaCO3, and ZrO2, and they found that the adsorption isotherms follow Langmuir-type behavior with an adsorption capacity of 1.23–3.67 mg/m2. In another study using palladium and nickel oxide incorporated into fumed silica nanoparticles, Franco et al. [38] reported a significant enhancement in the adsorption Colombian n-C7 asphaltenes.

Fig. 3.
figure 3

Kinetics of asphaltenes adsorption onto nanoparticles of γ-Fe2O3 obtained at three initial concentrations of (a) 100, (b) 500, and (c) 1000 ppm at 25°C.

Adsorption isotherm. Figure 4 shows the amount of asphaltenes adsorbed onto the nanoparticles as a function of the equilibrium concentration (ce) at 25°C. The figure shows that the adsorption of the studied asphaltenes increases at low concentration and begins to stabilize at a higher concentration, which suggests that γ-Fe2O3 nanoparticles have a good adsorption affinity for our asphaltenes. The adsorption isotherm profile (Fig. 4) is similar to a previous study on iron oxide nanoparticles for rapid adsorption of asphaltenes [14]. Nazila and Behruz confirmed that the maximum adsorption capacity of maghemite γ-Fe2O3 is greater than that of hematite (α-Fe2O3) [19].

Fig. 4.
figure 4

Adsorption isotherm of asphaltenes on nanoparticles of γ-Fe2O3. Dosage of nanoparticles: 10 g/L; contact time: 24 h; and temperature: 25°C. (1) Experimental data and (2) prediction from the SLE model [Eq. (4)].

The difference in adsorption capacity was attributed to the smaller size and larger surface area of maghemite compared to hematite [19]. Another important factor that may contribute to the adsorption capacity could be the acidity of the maghemite surface of γ-Fe2O3 [39]. To better understand this adsorption behavior, the experimental isotherm data were fitted to the SLE model (Solid−Liquid−Equilibrium) [40]. This model is based on the chemical theory for gas adsorption on solid at different load ranges, and thus it improves the understanding of asphaltenes aggregation and interactions onto the solid surface at equilibrium and provides a realistic representation of asphaltenes aggregation [41] as expressed by the following equation:

$$c = {{\psi H} \over {1 + K\psi }}\exp \left( {{\psi \over {{q_{\rm{m}}}A}}} \right),$$
((4))

where

$$\psi = {{ - 1 + \sqrt {1 + 4K\xi } } \over {2K}},\xi = \left( {{{{q_{\rm{m}}}q} \over {{q_{\rm{m}}} - q}}} \right)A,$$

q (mg/m2) and qm (mg/m2) are the amounts of asphaltenes adsorbed on the surface of nanoparticles and the maximum adsorption capacity of the nanoparticles, respectively. A (m2/mg) is the BET surface measured for the dried nanoparticles, and c (mg/L) is the equilibrium concentration of asphaltenes. The parameters of SLE model K (g/g) and H (mg/g) are the adsorption constants related to the degree of self-association of asphaltenes on the surface of the nanoparticles and Henry’s law constant related to asphaltene preference in adsorbed liquids, respectively [32]. The lower the value of H, the higher the affinity for the adsorbed phase. However, the degree of asphaltenes aggregation at the adsorption site is higher as the K value increases, thus affecting the degree of saturation of the adsorption site [35]. The SLE model parameters along with the correlation coefficient (R2) and the nonlinear analysis of the chi-squared (χ2) were used to assess the fit quality of the SLE model to the experimental data are listed in Table 2. All analyses were performed using Data Fit software (version 8.2.79, Oakdale Engineering, Oakdale, PA, USA). As shown in (Fig. 4), an excellent agreement was obtained between the SLE model and the experimental data, which is also indicated by the values of R2 (close to 1.0) and the low values of χ2 given in Table 2, which shows a strong tendency for asphaltene molecules to dimerize and aggregate. Consequently, it can be deduced that the adsorption of asphaltenes on the nanoparticles was multilayer adsorption. In other words, after the formation of the first layer, due to the high affinity of the asphaltene molecules for aggregation, the subsequent molecules are attached to the first layer and the multilayer adsorption takes place. The use of iron oxide nanoparticles for asphaltenes adsorption has been reported by many researcher groups. Nassar et al. [14], for instance, employed iron oxide nanoparticles, 43 mg/m2 surface area and ~25 nm particle size, for C7-asphaltenes adsorption extracted from Athabasca vacuum residue. It was found that the maximum adsorption capacity was 1.7 mg/m2 using the Langmuir isotherm model. The authors confirmed that the adsorption mechanism is dominated by polar interactions between adsorbate and adsorbent [15]. In another study using maghemite and hematite iron oxide for Iranian asphaltenes, the results showed that maximum adsorption capacity is 108.1 and 45.8 mg/g, respectively, using the SLE model [19]. Similarly, using Iranian asphaltenes, the maximum adsorption capacity of asphaltene onto Fe2O3, 44 nm particle size, was 3.52 mg/m2 followed Langmuir-type behavior. In comparison, the prepared iron oxide nanoparticles, in this study, provided a higher surface area with a higher affinity toward Algerian asphaltenes adsorption compared to those reported iron oxide nanoparticles. These results obtained can support the objective of this study to synthesize and improve maghemite iron oxide to absorb solid asphaltenes.

Table 2. SLE isotherm parameters for the asphaltenes adsorption onto nanoparticles of γ-Fe2O3 at a temperature of 25°C

CONCLUSIONS

A simple synthetic route was used to prepare γ-Fe2O3 nano-adsorbent and applied for the removal of asphaltenes in a toluene solution. The nanomaterial was routinely characterized, and the structure was confirmed. Unprecedented adsorption of asphaltenes was observed onto the γ-Fe2O3 nanoparticles. The iron oxide nanoparticles are promising nano-adsorbents that have a high affinity to remove the asphaltenes and the equilibrium was recorded after only 15 min. The experimental adsorption isotherm matched well and showed a good fit to the SLE model. This study proposes the use of those nanomaterials as potential adsorbents for asphaltenes to avoid their precipitation in pipelines or auxiliary equipment during the transportation of Algerian Hassi Messaoud oil.

ADDITIONAL INFORMATION

H. Abbas, ORCID: http://orcid.org/0000-0001-7169-2030

A.D. Manasrah, ORCID: http://orcid.org/0000-0002-2534-1739

A. Abidi Saad, ORCID: http://orcid.org/0000-0001-8189-9473

K.O. Sebakhy, ORCID: http://orcid.org/0000-0001-6620-0951

Y. Bouhadda, ORCID: http://orcid.org/0000-0003-0376-7453