Enhanced Oil Recovery Driven by Nanofilm Structural Disjoining

Apr 1, 2016 - At first glance, these results are opposite of what one expects in the conventional EOR, where oil recovery is known to increase progres...
12 downloads 0 Views 3MB Size
Download PDF
Article pubs.acs.org/EF

Enhanced Oil Recovery Driven by Nanofilm Structural Disjoining Pressure: Flooding Experiments and Microvisualization Hua Zhang,†,‡ T. S. Ramakrishnan,‡ Alex Nikolov,† and Darsh Wasan*,† †

Department of Chemical and Biological Engineering, Illinois Institute of Technology, Chicago, Illinois 60616, United States Schlumberger-Doll Research, 1 Hampshire Street, Cambridge, Massachusetts 02139, United States

‡

Energy Fuels 2016.30:2771-2779. Downloaded from pubs.acs.org by KAOHSIUNG MEDICAL UNIV on 10/19/18. For personal use only.

S Supporting Information *

ABSTRACT: Nanofluids composed of liquid suspensions of nanoparticles may soon permit accelerated recovery of hydrocarbons from oil and gas reservoirs. Here, we present a series of flooding experiments at different capillary numbers to quantify the performance of a polymeric nanofluid compared to brine using the sintered glass-beads. A high resolution X-ray microtomography (micro-CT) was used to visualize oil and brine distribution in a sintered bead pack before and after nanofluid flooding. The results of flooding experiments showed that an additional oil recovery of approximately 15% is possible with nanofluids compared to brine at low capillary numbers and is as effective as high capillary number brine flooding. Nanofluid induced additional oil recovery decreases as we increase the capillary number, and the total oil recovered shows a marginal decrease. At first glance, these results are opposite of what one expects in the conventional EOR, where oil recovery is known to increase progressively with increasing capillary number. These results cannot be explained based on mobilization theories due to the reduced capillarity. Our results however are consistent with the mechanism of wettability alteration caused by structural disjoining pressure leading to the formation of the wetting nanofluid film between oil and substrate. Unlike brine displacement, X-ray micro-CT images show that mobilization of oil from the porous medium by a nanofluid is fairly uniform even at low capillary numbers and large-scale clusters of oil are absent. Our findings in this paper are expected to promote the understanding of EOR by nanofluids.

1. INTRODUCTION Nanofluids, liquid suspensions of nanoparticles such as polymer latexes, globular proteins, and surfactant micelles, have recently become attractive agents for EOR from oil reservoirs.1−3 Researchers have observed positive results and an increase in ultimate oil recovery in laboratory experiments by injecting nanofluids.4−18 Ju et al.4−6 conducted experimental and theoretical studies on wettability alteration and permeability change caused by adsorption of lipophobic and hydrophilic polysilicon nanoparticles (LHP) on the surface of sandstone cores. From core displacement experiments, they observed LHP in the size range of 10−500 nm could improve oil recovery by about 9% compared to water. Based on their numerical simulation, an LHP concentration of 2−3% percent by volume was suggested for EOR, since porosity and permeability declined through retention of LHP within the medium at higher concentrations. Babadagli et al.7 studied capillary imbibition using different surfactants and polymer solutions for enhancing oil recovery. They conjectured that the reduction in the interfacial tension between the aqueous phase and oil was the cause of faster and improved oil recovery. Onyekonwu and Ogolo8 reported the performance of polysilicon nanoparticles for enhancing oil recovery. Three different types of polysilicon nanoparticles were used for the displacement experiments on water-wet rocks: lipophobichydrophilic (LHPN), hydrophobic-lipophilic (HLPN), and neutrally wet nanoparticles (NWPN). They found that NWPN and HLPN were good EOR agents in water-wet formations, while LHPN yield poor recovery factors indicating © 2016 American Chemical Society

that its use for EOR should be restricted to oil-wet formations. These investigators inferred that EOR was due to both the reduction of interfacial tension and wettability change. Hendraningrat et al.9−13 examined the ability of hydrophilic silica nanoparticles to improve oil recovery. Some parameters influencing EOR processes, such as the size and concentration of nanoparticles, initial core wettability, nanofluid injection rate, and temperature were investigated. They concluded that oil recovery improves with increasing temperature, decreasing nanoparticle size, and decreasing injection rate. Additionally, the highest oil recovery was obtained from intermediate-wet rocks. They showed that EOR using a silica based nanofluid resulted from the medium becoming water-wet from its oil-wet state. Moghaddam et al.14 compared eight different nanoparticles previously investigated in the literature. These were zirconium dioxide (ZrO2), calcium carbonate (CaCO3), titanium dioxide (TiO2), silicon dioxide (SiO2), magnesium oxide (MgO), aluminum oxide (Al2O3), cerium oxide (CeO2), and carbon nanotube (CNT). The ability to alter wettability of carbonate rocks was studied. Based on their impact on wettability alteration and stability in brine, CaCO3 and SiO2 were selected for spontaneous imbibition tests and core flooding experiments. The results of both imbibition and core flooding experiments demonstrated that SiO2 nanoparticles treatment improved oil Received: January 6, 2016 Revised: March 27, 2016 Published: April 1, 2016 2771

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels recovery more than the treatment by CaCO3 nanoparticles but over a longer duration. In summary, two classical mechanisms of EOR using nanofluids have been proposed, namely, the reduction of the interfacial tension between the aqueous phase and oil phase7,15 and rock wettability alteration.4−6,9−13,16,19 In some instances, both mechanisms are believed to be operational.8,17,18 Despite recently widely conducted research using a nanofluid for EOR, the underlying operating mechanism of recovery by a nanofluid is still not well understood. A new view of oil displacement from a solid substrate using a nanofluid was recently proposed by Wasan and Nikolov.20 In their paper, the nanoparticle structuring (layering) in the wedge film was advocated as the driving mechanism for the wettability alteration. They showed that nanoparticles (or surfactant micelles) form two-dimensional (2-D) layered structures in the confined three-phase (solid-oil-aqueous phase) contact region of the wedge film (Figure 1). The nanoparticle

Figure 2. Photomicrograph taken by reflected-light interferometry depicting the nanofluid structural disjoining pressure on oil displacement dynamics. The “outer line” is the macroscopic three-phase contact line, and the “inner line” depicts the advancing nanofluid film driven by structural disjoining pressure. Reprinted from ref 24. Copyright 2012 American Chemical Society.

additional 30% of crude oil was recovered using an IIT nanofluid at a reservoir temperature of 55 °C. High resolution X-ray microtomography (micro-CT) has been used extensively in the oil and gas industry for imaging, quantifying properties, and determining the distribution of fluids in porous rocks.28 The main advantage of X-ray microCT is the nondestructive nature of the technique that allows 3D monitoring of internal structural changes at resolutions in the range of a few microns.29,30 Extensive research on imaging the residual nonwetting phase using micro-CT scanning at ambient conditions where the porous medium was strongly water-wet or oil-wet has been done.31−38 Karpyn et al.34 experimentally investigated trapped oil clusters in a water-wet bead pack subject to different flow conditions using X-ray microtomography. They presented oil cluster shape and size distribution and found an average cluster size of approximately five pores. Our present paper is the continuation of the previous exploratory study27 and is related to the understanding enhanced oil recovery from a porous medium using an IIT nanofluid (patent applied): a polymeric nanofluid (polyethylene glycol dispersed in brine) at different capillary numbers. Since these polymeric nanofluids consist of particles in the range of about 10 nm (see below) and appear to be stable to salinity, they may be considered to be a nanofluid consisting of deformable particles. The primary purpose of our work with these fluids was to understand the removal of trapped oil ganglia (or clusters) formed by low capillary number brine flooding and comparing the results to brine injection at various capillary numbers, which mimic a conventional tertiary recovery process with reduced interfacial tension. Thus, the present study differs from the previous ones, since the objective here is to distinguish interfacial tension induced oil recovery from that due to wettability change by the structural disjoining pressure. The former reduces capillarity, whereas the latter may even enhance capillarity. Although we found that additional oil recovery was enabled by the presence of a nanofluid, the incremental recovery decreased with an increasing capillary number, while the total oil recovered showed a slight decrease. This nonintuitive result is understood by considerations of the nanofluid film formation time scale in relation to the advection time scale. We provide detailed images

Figure 1. Schematic presentation of mechanism of oil displacement driven by nanofluid structural disjoining pressure: nanofluid layering leads to film tension gradient and drives oil displacement. Reprinted with permission from ref 20. Copyright 2003 Nature Publishing Group.

structuring phenomenon gives rise to the structural disjoining pressure (a force normal to the interface) in the wedge film with a higher disjoining pressure near the tip of the wedge than in the bulk meniscus. As a result, the oil-nanofluid interface moves forward with a progressive cleaving of the oil-solid contact eventually detaching the oil drop. Kao et al.21 observed two distinct contact lines during the separation of crude oil droplets from a solid silica surface in the presence of a nanofluid (micellar solution) using a differential interference microscope: an outer one (between the oil, solid and water film) and an inner one (between the oil, solid, and mixed oil/water film) (see Figure 2). The spreading of a mixed oil/water film was later understood to be driven by the structural disjoining pressure gradient arising from the ordering of the micelles in the wedge-film region by Wasan and Nikolov.20 Wasan, Nikolov, and their co-workers22−26 further noted that the dynamics of the inner contact line depend on the combination of the nanoparticle formulation, contact angle, and the capillary pressure. A suitable combination of these factors accelerates the spreading of the nanofluid on the solid surface, thereby detaching the oily soil from the substrate. In the case of high salinity, the effective diameter (including electrical double layer) (and thereby the particle volume fraction) decreases, owing to the shrinkage of the electrical double layer around each particle. Therefore, the magnitude of the structural disjoining pressure diminishes, thus reducing the driving force for the detachment of the oil drop. The application of nanofluids was also explored for hydrocarbon recovery from oil saturated Berea sandstone.27 Compared to brine (containing 0.32 mol/L NaCl), an 2772

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels

by volume ratio) and sintered in an oven at about 830 °C. Figure 3 shows a picture of an assembled bead pack.

of residual oil from X-ray micro-CT, both with a nanofluid and brine, and contrast the results by relating them to wettability.

2. EXPERIMENTAL SECTION 2.1. Materials. 2.1.1. Oil. Oil used in this study is Cargille immersion liquid (Cargille Laboratories, NJ, USA). Properties of the Cargille immersion liquid are given in Table 1.

Table 1. Oil Propertiesa

a

property

value

density, g mL−1 viscosity, Pa s surface tension, mN m−1 refractive index

0.854 0.0181 29.4 1.474

Figure 3. Photo of a sintered bead pack with a mixture of 1 mm and 250−300 μm borosilicate glass beads.

25 °C, 1 atm pressure.

The overall porosity of the sintered bead pack was determined in an AccuPyc 1330 Pyconometer (Micromeritics Instrument Corporation) and double checked by an X-ray micro-CT scanned images and threshold processing using Otsu’s algorithm.41 The porosity measured by the pyconometer is 28.2 ± 1.2% compared to 28.6% obtained by thresholding micro-CT scan rendered image. A permeability of 54.73 Darcy for the bead pack was measured experimentally by draining water from a hydrostatic column through the pack. At the conclusion of each flooding experiment, the bead pack was cleaned by sequentially flushing the medium with isopropyl alcohol, toluene, and isopropyl alcohol to restore wettability. After drying the bead pack overnight, the same pack was used for the next injection sequence. Thus, the pore structure of the porous medium was unaltered in the sequence of experiments. The cleaning procedure was chosen to make sure that the wettability of glass beads was the same for each run. This was confirmed by measuring water contact angle of 27 ± 3° on clean glass beads in air. After exposure to oil, and following the above-mentioned flushing steps and air drying, the water contact angle (with respect to air) was measured to be 29 ± 3°, essentially restoring the contact angle to its original state. 2.2. Interfacial Tension Measurement. The classical method of drop-shape analysis was used to calculate the interfacial tension of oil/ brine, oil/polymeric nanofluid, and oil/SDS solution.27 The measured interfacial tensions of oil/brine, oil/polymeric nanofluid, and oil/SDS are 43 ± 2 mN m−1, 8.8 ± 0.3 mN m−1, and 8.7 ± 0.3 mN m−1 respectively. 2.3. Oil/Brine (Nanofluid)/Borosilicate Glass Plate ThreePhase Contact Angle Measurement. The precleaned borosilicate glass substrates were placed in a transparent plastic cuvette filled with oil (without dye). The glass substrates were cleaned sequentially with isopropyl alcohol, toluene, and isopropyl alcohol followed by air drying overnight. A sessile droplet of brine or a nanofluid (0.3 ± 0.05 mL, with methylene blue dye) was placed on the glass substrate immersed in the oil with a Hamilton microsyringe (B-D needle, 25G 5/8). After placing the droplet on the solid substrate, the image of the droplet on the glass substrate was captured with a digital camera (Canon A720 IS) until the droplet stopped spreading. The static three-phase contact angle for brine and a nanofluid is 65 ± 4° and 25 ± 3° respectively. The advancing and receding contact angle was measured by the tilting plate method. Basically, after the deposited droplet stopped spreading on the solid substrate, we slowly tilted the sample until the sessile drop began to move in the downhill direction. The downhill contact angle is the advancing angle, and the uphill angle is the receding contact angle. The contact angle was calculated by image analysis (Image Pro, version 6). The measured brine advancing and receding contact angles are 88 ± 2° and 42 ± 3°, and the nanofluid advancing and receding contact angles are 25 ± 3° and 22 ± 3°. All of the experiments were conducted at a room temperature of 22 ± 1 °C. 2.4. Flooding Experiments. 2.4.1. Experimental Setup. The aim of the experiments was to compare and contrast the residual oil reduction with and without a nanofluid as a function of capillary number (Ca). All of the flooding experiments were conducted at room

A small amount (100 ppm) of Oil Red EGN dye (Aldrich Chemical Co., Milwaukee, WI) was added to the oil sample to increase visual contrast for observational ease. 2.1.2. Brine Solution. Brine was prepared by dissolving sodium chloride (NaCl, Fisher Scientific, USA) in deionized (DI) water. The concentration of NaCl was 0.25 wt % of solution. Density of brine was 0.99 ± 0.01 g/cm3, and pH was 6.5 ± 0.5 at 22 °C and 1 atm pressure. The viscosity of brine is 1.022 mPa s at 20 °C. 2.1.3. SDS Solution. We prepared sodium dodecyl sulfate (SDS, VWR Scientific, USA) solution in brine to a concentration of 1.35 mM so that the SDS solution has the same interfacial tension with oil as that of a nanofluid. 2.1.4. Polymeric Nanofluid. Most reservoir environments are at high temperature, pressure, and salinity. Most nanofluids, such as silica nanoparticle dispersions, are unstable and agglomerate in such environments. In contrast, a polymeric nanofluid containing polyethylene glycol 8000 (Fisher Scientific, USA) dispersed in brine is relatively stable with respect to electrolytes or temperature. In order to enhance the effect of the structural disjoining pressure on the oil recovery process, the nanofluid composition was selected based on a multistep process.27,39 The size and polydispersity of the polymeric nanofluid in brine were characterized by the dynamic light scattering method (Malvern Instruments, UK). The average diameter was found to be 9.5 ± 0.5 nm, with a polydispersity of around 8−10% at 25 °C. The concentration of the nanofluid used in this study is 0.277 wt % or 10 v% calculated based on the molar concentration of the polymer.40 The density of the nanofluid was 1.00 ± 0.01 g mL−1 with a pH = 6.5 ± 0.5 at 22 °C at 1 atm pressure, and the viscosity of the nanofluid is 1.097 mPa s at 20 °C. Table 2 shows the size and polydispersity of the polymeric nanofluid at two different salinities.

Table 2. Nanofluid Size and Polydispersity Characterized by Dynamic Light Scattering at Different Salinities at 25 °C concn of NaCl (wt %)

size (nm)

polydispersity

0.025 0.25

10.5 ± 1.0 9.5 ± 0.5

0.12−0.15 0.08−0.1

The near constancy of average diameter and narrow dispersity even with varying salinity illustrates that the polymeric dispersion is stable. As seen in Table 2, the effective diameter of the polymer particle is nearly the same at two different salinities. The small change in size might be due to hydrogen bonding of water molecules to electron-rich oxygen atoms in the polymer chain causing the polymer to be slightly charged. The polymer coil is “relaxed (stretched)” in low salinity and “constrained” in high salinity. The apparent small reduction in size is conjectured to be driven by a reduced electrical double layer with ionic strength. 2.1.5. Sintered Bead Pack. A quartz tube was filled with borosilicate glass beads of two different diameters of 1 mm and 250−300 μm (1:1 2773

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels temperature (22 ± 1 °C). A schematic diagram of the experimental setup is presented in Figure 4.

the volume of oil recovered has to be equal to the volume of brine injected according to volume conservation. Oil production and the pressure drop across the bead pack were recorded as a function of pore volumes. Thus, secondary water-flood residual oil saturation was established. The three steps of evacuation, oil intrusion, and low capillary number water-flood were the same for each flooding experiment; the results of which are presented below. Following water-flood, in the first experiment, a nanofluid was injected at a Ca = 10−7 ranging from three to seven PV mimicking a tertiary recovery process. The process was however unusual because Ca was held to values typical of a secondary water-flood. Nanofluid injection was terminated when no discernible oil was evident in the outlet stream. Oil cut and pressure were recorded. As mentioned previously, at the conclusion of the flooding experiment at Ca = 10−7, the bead pack was cleaned by flushing the medium with isopropyl alcohol, toluene, and isopropyl alcohol, essentially restoring the contact angle to its original state. To delineate the effect of reduced capillarity from that of wettability (other than inevitable changes through dynamic contact angle), no surfactant was added in the brine for the high Ca tertiary brine displacement. The flow rate was used to set Ca. Through two different tertiary injections at the same Ca with and without the nanofluid, a distinction in the outcome would be indicative of the effect of the disjoining pressure on oil removal. Figure 5 presents a schematic sequence of the experimental procedure. The displacement efficiency was then evaluated and calculated based on

Figure 4. Schematic diagram of the experimental setup. The experimental setup consisted of syringe pumps, pressure transducers, the bead pack, stainless syringes, and 1/8″ (3.175 mm) PFA tubings. The syringe pumps (Harvard Apparatus) were operated at constant flow rate. The pump reproducibility is ±0.05%. Using a 20 mL stainless steel syringe, the pump’s minimum and the maximum flow rates are 52.86 nL/min and 54.89 mL/min, respectively. The inlet and outlet of the bead pack were connected to pressure transducers (accuracy class 0.1%, supplied by HBM) for measuring pressure drop across the bead pack during the flooding experiments. 2.4.2. Flooding Procedure. The capillary number is proportional to the ratio of viscous flow induced pressure drop to capillary pressure. By definition, in this paper, the capillary number is

Ca =

μU γ

⎡ ⎛ S ⎞⎤ ED = ⎢1 − ⎜ or 2 ⎟⎥ ⎢⎣ ⎝ Sor1 ⎠⎥⎦

(2)

where Sor1 represents the oil saturation after brine injection at Ca = 10−7, and Sor2 represents the residual oil saturation after tertiary brine or nanofluid injection, i.e., at high Ca. 2.5. Imaging Technique. 2.5.1. 3D Visualization Using X-ray Microtomography. A high-resolution X-ray microtomography (X-ray micro-CT) was used to determine oil and brine (nanofluid) distributions in the sintered bead pack before and after injecting a nanofluid. The micro-CT system consists of an X-ray generator, a detector, a translation system, and a computer system that controls motions, data acquisition, and reconstruction. Once the sample is brought to its position and the scanner is activated, the projections of the magnified object are captured by the image intensifier. When the source is active, the X-rays penetrate through the object and reach the image intensifier, which converts X-ray energy into a form of light that is captured by a digital camera. The digitized data are sent to the computer and turned into raw files that can be processed into images.43 The images are stored as a stack of virtual 2D slices (a total of 3142 slices). The sample is rotated 360° in the X-ray beam, while the detector is providing attenuation views to the data acquisition computer. After the scan is complete, the 2D slices are reconstructed to 3D volume using VG studio Max 2.2. The standard scanning parameters in our cases were as follows: four frames per projection, exposure time one second, angle of rotation 360°, X-ray tube voltage 160 kV, and current 65 μA. 2.5.2. Imaging Procedures. Initially, the dry bead pack was scanned to obtain the overall porosity. After the core reached oil and brine saturation, it was scanned again to get initial oil and brine distribution. The flooding sequence was then followed by consecutive cycles of brine and nanofluid injection. The core was scanned before and after

(1)

where μ is the fluid viscosity,U is the superficial velocity, and γ is the interfacial tension between the aqueous phase and oil. The contact angle is not included in this definition, since the primary purpose is to delineate the effects of wettability alteration as discussed further below. We changed the pump flow rate to change the superficial velocity to realize different capillary numbers. During oil recovery, the displacing fluid is synonymous with brine (w is used as a subscript), and the displaced fluid is identified as oil (o is the subscript). The enhanced recovery experiment operates on a tertiary recovery mode as opposed to secondary enhanced recovery, meaning that a low Ca flow always precedes a high Ca flood. A pre-evacuated bead pack is filled with brine in order to avoid trapped air bubbles. Oil is then injected at a sufficiently high rate so that the capillary number is about 10−4. The primary purpose of injecting oil at this high rate is to reduce the zone of end-effect.42 The effluent brine and oil are collected, and thus the saturation of oil and brine are known. At this point, the bead pack is filled predominantly with oil except for residual brine saturation, Swr, of 15−20%. For oil recovery experiments, we injected three to four pore volumes (PV) brine at Ca = 10−7 into the oil-filled bead pack to simulate secondary water-flooding. At the end of this step, no oil production was observed. The breakthrough point was observed at about 0.48−0.5 PV for the brine flooding (slightly change due to different runs). Since the fractional flow of oil is unity until this point, recovery could be computed from the injection rate, and is shown accordingly in the oil recovery figures. Before the breakthrough point,

Figure 5. Schematic sequence of the experimental procedure showing the flooding process. 2774

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels

Figure 6. (a) Oil recovery vs injected PV for brine and nanofluid secondary and tertiary injection at Ca = 10−7 and T = 22 ± 1 °C. (b) Pressure drop across the bead pack during brine and nanofluid injection at Ca = 10−7 and T = 22 ± 1 °C.

Figure 7. Oil recovery vs injected PV and EOR process with a nanofluid compared to brine at (a) Ca = 10−5 and (b) Ca = 10−4 at a temperature of 22 ± 1 °C. nanofluid injection to map the oil distribution in order to see the performance of a nanofluid compared to that of brine. A digital camera (Canon PowerShot A720 IS) equipped with a 58 mm macrolens was also used to take pictures from the side at the end of each PV flooding.

absent, although a slight decay is evident. The trapped oil ganglia inside the bead pack started to break into disconnected oil blobs and was continuously perturbed and displaced as we injected the nanofluid. This will be clearer in Section 3.2 where micro-CT images highlight differences between brine and nanofluid flooding. As mentioned previously, the second set of experiments comprises a sequence of low and high capillary numbers, with the latter again being both brine and a nanofluid. High capillary numbers were achieved by increasing the injection rate. The outcome of the second sets of experiments was expected to clearly delineate the difference between reduced capillarity without wettability change and increased oil recovery including wettability change. Figure 7 shows oil recovery during EOR at Ca = 10−5 and Ca = 10−4. Additional oil due to nanofluid injection is 6.9% at Ca = 10−5 and 2.4% at Ca = 10−4. This is a downward trend with respect to incremental oil production, but the interesting observation is that with respect to Ca, even the total oil produced showed a slight decrease to no change after tertiary injection of a nanofluid. In Table 3 we have summarized the results of the flooding experiments.

3. RESULTS AND DISCUSSION 3.1. Flooding Results for Different Flooding Scenarios. Figure 6(a) shows the oil recovery versus pore volume of the fluid injected for secondary (waterflood) and EOR processes. From Figure 6(a), we observe that around 65% oil was displaced by brine. Oil production ceased at about one PV, but in order to standardize the procedure, we continued injection until three PV brine was injected. After brine injection, we continued with nanofluid injection as the tertiary mode but at the same capillary number (Ca = 10−7). Interestingly, 14.6% of additional oil was recovered using the nanofluid compared to 0% using brine at Ca = 10−7. Figure 6(b) shows the pressure drop across the porous medium during injection at Ca = 10−7. It was observed that the average pressure drop was almost constant during the brine injection, and the value is close to zero since the permeability of the core is large and is approximately 55 D (Darcy). The pressure drop increased and is oscillatory when the nanofluid was injected to the bead pack. The viscosity of the polymeric nanofluid was 1.097 mPa s. Given that the viscosity of oil was about 18 times higher, one would expect a higher pressure when oil flow dominates within the porous medium. Our observations are counter to this. We find that the pressure drop is about 0.1 to 0.15 psi higher during nanofluid injection, and about 20 lamellae in the flow line are sufficient to induce this pressure drop. Oscillations in pressure drop are coincident with oil recovery (compare Figures 6a and 6b). Interestingly, the pressure drop of 0.1 psi is retained (the reason for this 0.1 psi is not clear-cut) even after oil lamellae in the tube are visually

Table 3. Oil Recovery at Secondary Process and EOR Using a Nanofluid and Brine vs Ca at a Temperature of 22 ± 1 °C secondary process

2775

oil recovery after 3 PV brine injection at Ca = 10−7, % 65.0 ± 3 63.2 ± 2 62.2 ± 2

tertiary process (EOR)

Ca

oil recovery after brine flooding, %

oil recovery after nanofluid flooding, %

additional oil compared with brine, %

10−7 10−5 10−4

0 9.1 12.3

14.6 ± 1.5 16.0 14.7

14.6 ± 1.5 6.9 2.4

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels

nanofluid can perform better when it is in contact with the solid substrate and the oil for a longer time, two sets of flooding experiments were conducted. In the first set, we flooded the pack at Ca = 10−4 until oil production ceased, subsequent to which we decreased the capillary number to Ca = 10−7. With respect to pore volumes injected, Figure 8 shows an additional 4.1% oil removal. The time scale of course is considerably larger than with the higher Ca.

As seen from Table 3, oil recovery due to nanofluid tertiary recovery decreases compared to brine as we increase the capillary number contrary to intuitive expectation of improved oil displacement. It is in the limit of low Ca that the highest additional recovery of 14.6% was obtained. This is to be tempered by the fact that the enhanced brine flood also recovers extra oil, and therefore the residual oil at the onset of nanofluid flooding is lower than in the lower Ca floods; but even as a percentage, oil recovered by a nanofluid with respect to oil left behind by brine, we see that enhancing Ca decreases efficacy. Similar results have been observed recently by other researchers9 showing that increasing silica nanofluid injection rate significantly decreases incremental oil recovery. In some cases, they observed that there was no incremental oil recovery. Their explanation is that increasing the injection rate might affect nanoparticle accumulation near the core inlet rather than flowing through the pore throat. They observed visually that nanoparticles’ “cake” at core inlet was more noticeable at a higher injection rate. The situation is different in our case. The bead pack is more porous than the core samples used in ref 9 (permeability is around 1000 times larger and porosity is around nearly twice); in addition, we do not observe any aggregates at the inlet. The reason why a nanofluid performs best at a low capillary number of 10−7 is that a nanofluid requires time for a film to form as a consequence of structural disjoining pressure to perform; for example, it takes around 9 h for one PV of nanofluid flooding at Ca = 10−7, while it takes only around 30 s for one PV of nanofluid flooding at Ca = 10 −4. From our previous experimental studies on dynamic spreading of nanofluids on solids24 and a recent mechanistic study of oil detachment from glass substrate using a nanofluid,27 the time scale for structural disjoining pressure mechanism to operate must be taken into account since the inner contact line moves much slower than that of the outer contact line, the latter being driven by capillarity. In our case, the velocity of the inner contact line as measured is 0.032 mm/min and is 0.32 mm/min for the outer contact line. The apparent decrease in efficacy of a nanofluid with increasing Ca is best understood by comparing advection and nanofilm-formation time scales, Ta and Tf, respectively. In order for the structural disjoining pressure to operate, the characteristic time for nanofluid advection inside the bead pack should be much larger than that for nanofluid film formation, i.e., Ta ≫ Tf. For a bead pack with length L, porosity ϕ, and superficial velocity U, the characteristic time for nanofluid advection inside Lϕ the bead pack is Ta = U . The characteristic time for nanofluid film formation is Tf =

DP , vin

Figure 8. Oil recovery vs injected PV and EOR process with nanofluid flooding at different capillary numbers.

Since the time scale is of primary importance, in the Ca = 10−7 experiments postnanofluid flooding, the bead pack was left in contact with the injected fluid, and the pump turned off (Ca = 0) for different periods (57 h and 10 days). This was followed by additional injection at Ca = 10−7. Figure 9 shows the resulting additional oil recovery of 3.2% (10 day soak-time) and 1.9% (57 h soak-time).

Figure 9. Oil recovery vs injected PV and EOR process with different nanofluid soak-times.

where Dp is bead diameter, and νin is

These two sets of experiments further verify our structural disjoining pressure mechanism. The longer the time the nanofluid was in contact with the solid substrate and the oil, the better was the performance. The nanofluid can perform better at a low capillary number and as effectively as brine flooding at a higher capillary number. Similar experiments with soaking followed by displacement with brine or SDS (the same IFT as polymeric nanofluid) resulted in no additional recovery. To exclude interfacial tension induced oil recovery from that due to nanofluid structural disjoining pressure, an additional experiment was conducted by flooding with 1.35 mM SDS solution (the same IFT as polymeric nanofluid) at Ca = 10−7 after brine injection at Ca = 10−7. Figure 10 shows that 14.6% of additional oil was recovered using a nanofluid compared to

inner contact line velocity. Table 4 summarizes the results. From Table 4, we see that only the lowest Ca = 10−7 has Ta ≫ Tf, which also corresponds to the highest EOR by a nanofluid compared with brine. In order to verify that a Table 4. Characteristic Time for Nanofluid Advection along the Bead Pack and Nanofluid Film Formation vs Ca Ca

Ta, min

Tf, min

10−7 10−5 10−4

363.70 3.64 0.36

31.30

2776

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels

gravitational effect during the flooding process and oil was more uniformly displaced along the porous medium. The density of the nanofluid is substantially the same as that of brine, but the interfacial tension with respect to oil is considerably smaller. Therefore, one would expect capillarity to dominate gravity for brine in comparison to a nanofluid. The results are however quite contrary, and it is useful to consider the Bond number. Since the only variable not considered thus far is the contact angle, we analyze its effect through Bond number

Bo =

ΔρgR2 2γ cos θA

(3)

during where Δρ is the density difference between brine and oil, g is the gravitational acceleration, R is the bead radius, γ is the interfacial tension between aqueous phase and oil, and θA is the advancing contact angle. The calculated Bo for brine is 0.119 compared to that of 0.022 for a nanofluid. Clearly gravity may not be neglected during brine advance but plays a negligible role during nanofluid injection. Even with an advancing contact angle of 88°, Ca is still sufficiently small for ensuring oil trapping. The receding contact angle for brine is 42 ± 3°, indicative of water wetness. Video Clip 1 (see the Supporting Information) depicts oil (red) and brine (light blue) distribution after 3 PV brine flooding at Ca = 10−7, and Video Clip 2 (see the Supporting Information) shows oil and brine distribution after 3 PV nanofluid flooding at Ca = 10−7. Figure 12 presents trapped oil ganglia (or blobs) before and after 3 PV nanofluid flooding. After 3 PV brine flooding, oil

Figure 10. Oil recovery vs injected PV and EOR process with a nanofluid compared to 1.35 mM SDS solution at Ca = 10−7 at a temperature of 22 ± 1 °C.

2% using SDS solution at Ca = 10−7. The additional oil 2% by SDS solution is purely due to the interfacial tension reduction. Since the capillary number was kept the same, the role of the interfacial tension would be through an alteration in contact angle (glass/oil/aqueous three-phase contact angle changes from 115° to 155°). Compared to oil recovery by SDS solution, the additional oil by a nanofluid is 12.6%, which is due to the structural disjoining pressure. 3.2. Oil and Brine Distribution Using X-ray Micro-CT. Figure 11(a) shows the side views from the digital camera. The

Figure 12. Oil blobs before and after 3 PV nanofluid flooding at Ca = 10−7. Figure 11. Visualization of oil and brine distribution inside the bead pack after brine and nanofluid flooding (Ca = 10−7). (a) Side view from the digital camera and (b) 3D view from X-ray micro-CT; red is oil and light blue is aqueous phase.

blobs are in large clusters and predominantly on the upper part of the porous medium. Nanofluid flooding breaks the large structures into small disconnected oil blobs of a few bead lengths.

4. CONCLUSIONS In this work, a series of flooding experiments on EOR using a nanofluid at different capillary numbers (Ca = 10−7, 10−5, 10−4) were performed and compared with brine after secondary flooding at a Ca = 10−7. Additional oil recovery of around 15% may be achieved by a nanofluid compared to brine at a Ca = 10−7. The additional oil recovery due to the nanofluid decreases as we increase the capillary number due to two reasons: (i) there is less oil to be recovered because of mobilization due to enhanced Ca brine flooding and (ii) insufficient time for the nanofluid film to advance through structural disjoining pressure and displace oil. Allowing the bead pack to be treated with a

porous medium is predominantly filled with oil when brine was displaced by oil. The oil saturation was between 80 and 85%. After 3 PV of brine flooding at Ca = 10−7, we see that the gravitational effect is quite noticeable, and oil was mostly displaced from the bottom and is retained toward the top. This apparently gravity induced entrapped oil was displaced (although incompletely) when injection was switched to a nanofluid at Ca = 10−7. Figure 11(b) shows 3D images taken from X-ray micro-CT. In the 3D view, it is easier to see that gravity played an important role during brine flooding, although at a Ca = 10−7 one would expect capillarity to dominate over gravity. The nanofluid caused an apparent reduction in the 2777

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels

(7) Babadagli, T. Dynamics of capillary imbibition when surfactant, polymer, and hot water are used as aqueous phase for oil recovery. J. Colloid Interface Sci. 2002, 246, 203−213. (8) Onyekonwu, M.; Ogolo, N. A. Investigating the use of nanoparticles in enhancing oil recovery. SPE 2010, 10.2118/140744MS. (9) Torsæter, O.; Li, S.; Hendraningrat, L. Effect of some parameters influencing enhanced oil recovery process using silica nanoparticles: An experimental investigation. SPE 2013, 10.2118/165955-MS. (10) Torsæter, O.; Li, S.; Hendraningrat, L. Enhancing oil recovery of low-permeability Berea sandstone through optimized nanofluids concentration. SPE 2013, 10.2118/165283-MS. (11) Torsæter, O.; Li, S.; Hendraningrat, L. A coreflood investigation of nanofluid enhanced oil recovery in low-medium permeability Berea sandstone. SPE 2013, 10.2118/164106-MS. (12) Hendraningrat, L.; Li, S.; Torsæter, O. A coreflood investigation of nanofluid enhanced oil recovery. J. Pet. Sci. Eng. 2013, 111, 128− 138. (13) Hendraningrat, L.; Torsæter, O. Effects of the initial rock wettability on silica-based nanofluid-enhanced oil recovery processes at reservoir temperatures. Energy Fuels 2014, 28, 6228−6241. (14) Moghaddam, R. N.; Bahramian, A.; Fakhroueian, Z.; Karimi, A.; Arya, S. Comparative study of using nanoparticles for enhanced oil recovery: Wettability alteration of carbonate rocks. Energy Fuels 2015, 29, 2111−2119. (15) Suleimanov, B. A.; Ismailov, F. S.; Veliyev, E. F. Nanofluid for enhanced oil recovery. J. Pet. Sci. Eng. 2011, 78, 431−437. (16) Karimi, A.; Fakhroueian, Z.; Bahramian, A.; Khiabani, N. P.; Darabad, J. B.; Azin, R.; Arya, S. Wettability alteration in carbonates using zirconium oxide nanofluids: EOR implications. Energy Fuels 2012, 26, 1028−1036. (17) Shahrabadi, A.; Bagherzadeh, H.; Roustaei, A.; Golghanddashti, H. Experimental investigation of HLP nanofluid potential to enhance oil recovery: A mechanistic approach. SPE 2012, 10.2118/156642-MS. (18) Aoudia, M.; Al-Maamari, R. S.; Nabipour, M.; Al-Bemani, A. S. Laboratory study of alkyl ether sulfonates for improved oil recovery in high-salinity carbonate reservoirs: A case study. Energy Fuels 2010, 24, 3655−3660. (19) Zhang, H.; Nikolov, A. D.; Wasan, D. T. Dewetting film dynamics inside a capillary using a micellar nanofluid. Langmuir 2014, 30, 9430−9435. (20) Wasan, D. T.; Nikolov, A. D. Spreading of nanofluids on solids. Nature 2003, 423, 156−159. (21) Kao, R. L.; Wasan, D. T.; Nikolov, A. D.; Edwards, D. A. Mechanisms of oil removal form a solid surface in the presence of anionic micellar solutions. Colloids Surf. 1988, 34, 389−398. (22) Chengara, A.; Nikolov, A. D.; Wasan, D. T.; Trokhymchuk, A.; Henderson, D. Spreading of nanofluids driven by the structural disjoining pressure gradient. J. Colloid Interface Sci. 2004, 280, 192− 201. (23) Kondiparty, K.; Nikolov, A. D.; Wu, S.; Wasan, D. T. Wetting and spreading of nanofluids on solid surfaces driven by the structural disjoining pressure: Statics analysis and experiments. Langmuir 2011, 27, 3324−3335. (24) Kondiparty, K.; Nikolov, A. D.; Wu, S.; Wasan, D. T.; Liu, K. Dynamic spreading of nanofluids on solids Part I: Experimental. Langmuir 2012, 28, 14618−14623. (25) Liu, K.; Kondiparty, K.; Nikolov, A. D.; Wasan, D. T. Dynamic spreading of nanofluids on solids Part II: Modeling. Langmuir 2012, 28, 16274−16284. (26) Wu, S.; Nikolov, A. D.; Wasan, D. T. Cleansing dynamics of oily soil using nanofluids. J. Colloid Interface Sci. 2013, 396, 293−306. (27) Zhang, H.; Nikolov, A.; Wasan, D. T. Enhanced oil recovery (EOR) using nanoparticle dispersions: Underlying mechanism and imbibition experiments. Energy Fuels 2014, 28, 3002−3009. (28) Dunsmuir, J. H.; Ferguson, S. R.; D’Amico, K. L.; Stokes, J. P. Xray microtomography: a new tool for the characterization of porous media. In Paper SPE 22860, SPE 1991, 10.2118/22860-MS.

nanofluid and then displacing oil confirms that the second hypothesis is true. Improved oil recovery due to a nanofluid is effective when Ta ≫ Tf. To separate the interfacial tension induced oil recovery from that due to wettability alteration by structural disjoining pressure, a flooding experiment with SDS solution (the same IFT as polymeric nanofluid) was conducted at Ca = 10−7 after brine injection at Ca = 10−7. Compared to oil recovery by SDS solution, the additional oil by a nanofluid is 12.6% and is due to the structural disjoining pressure. High-resolution X-ray microtomography 3D views showed that gravity played a more pronounced role during brine injection compared to nanofluid injection and was initially unexpected. It appears that this was primarily driven by reduced capillarity during brine advance. Nanonfluids are as effective as or better than high Ca brine flooding. Allowing a nanofluid sufficient contact time appears to be the most effective for total oil recovered. Our findings in this study are expected to promote the understanding of EOR by nanofluids.

■

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.energyfuels.6b00035. Video Clip 1 of oil (red) and brine (light blue) distribution after 3 PV brine flooding at Ca = 10−7 (AVI) Video Clip 2 of oil and brine distribution after 3 PV nanofluid flooding at Ca = 10−7 (AVI)

■

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This research is supported by Schlumberger-Doll Research Center. One of the authors, H.Z., wishes to thank Schlumberger for support through an internship. Assistance from Raji Shankar and Albert Perez, Jr. is much appreciated.

■

REFERENCES

(1) Kong, X.; Ohadi, M. M. Application of micro and nano technologies in the oil and gas industry-An overview of the recent progress. SPE 2010, 10.2118/138241-MS. (2) Wong, K. V.; Leon, O. Applications of nanofluids: Current and future. Adv. Mech. Eng. 2010, 2, 519659. (3) Ayatollahi, S.; Zerafat, M. M. Nanotechnology-assisted EOR techniques: New solutions to old challenges. SPE 2012, 10.2118/ 157094-MS. (4) Ju, B.; Fan, T.; Ma, M. Enhanced oil recovery by flooding with hydrophilic nanoparticles. China Particuol. 2006, 4, 41−46. (5) Ju, B.; Fan, T. Experimental study and mathematical model of nanoparticle transport in porous media. Powder Technol. 2009, 192, 195−202. (6) Ju, B.; Dai, S.; Luan, Z.; Zhu, T.; Su, X.; Qiu, X. Study of wettability and permeability change caused by adsorption of nanometer structured polysilicon on the surface of porous media. SPE 2002, 10.2118/77938-MS. 2778

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779

Article

Energy & Fuels (29) Landis, E.; Keane, D. X-ray microtomography. Mater. Charact. 2010, 61, 1305−1316. (30) Auzerais, F. M.; Dunsmuir, J.; Ferreol, B. B.; Martys, N.; Olson, J.; Ramakrishnan, T. S.; Rothman, D. H.; Schwartz, L. M. Transport in sandstone: a study based on three dimensional microtomography. Geophys. Res. Lett. 1996, 23, 705−708. (31) Wildenschild, D.; Sheppard, A. P. X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems. Adv. Water Resour. 2013, 51, 217−246. (32) Prodanović, M.; Lindquist, W. B.; Seright, R. S. 3D image-based characterization of fluid displacement in a Berea core. Adv. Water Resour. 2007, 30, 214−226. (33) Kumar, M.; Senden, T. J.; Sheppard, A. P. Visualizing and quantifying the residual phase distribution in core material. Petrophysics 2010, 51, 323−332. (34) Karpyn, Z.; Piri, M.; Singh, G. Experimental investigation of trapped oil clusters in a water-wet bead pack using X-ray microtomography. Water Resour. Res. 2010, 46, W04510. (35) Knackstedt, M.; Senden, T.; Carnerup, A.; Fogden, A. Improved characterization of EOR processes in 3D. Characterizing mineralogy, wettability and residual fluid phases at the pore scale. SPE 2011, 10.2118/145093-MS. (36) Iglauer, S.; Favretto, S.; Spinelli, G.; Schena, G.; Blunt, M. J. Xray tomography measurements of power-law cluster size distributions for the nonwetting phase in sandstones. Phys. Rev. E 2010, 82, 056315. (37) Landry, C. J.; Karpyn, Z. T.; Piri, M. Pore-scale analysis of trapped immiscible fluid structures and fluid interfacial areas in oil-wet and water-wet bead packs. Geofluids 2011, 11, 209−227. (38) Iglauer, S.; Fernø, M. A.; Shearing, P.; Blunt, M. J. Comparison of residual oil cluster size distribution, morphology and saturation in oil-wet and water-wet sandstone. J. Colloid Interface Sci. 2012, 375, 187−192. (39) Nikolov, A.; Zhang, H. The dynamics of capillary-driven twophase flow: The role of nanofluid structural forces. J. Colloid Interface Sci. 2015, 449, 92−101. (40) Lim, S.; Zhang, H.; Wu, P.; Nikolov, A.; Wasan, D. The Dynamic Spreading of Nanofluids on Solid Surfaces–Role of the Nanofilm Structural Disjoining Pressure. J. Colloid Interface Sci. 2016, 470, 22−30. (41) Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans. Sys., Man. Cyber. 1979, 9, 62−66. (42) Ramakrishnan, T. S.; Cappiello, A. A new technique to measure static and dynamic properties of a partially saturated porous medium. Chem. Eng. Sci. 1991, 46, 1157−1163. (43) Alvarado, F. E.; Grader, A. S.; Karacan, O.; Halleck, P. M. Visualization of three phases in porous media using microcomputed tomography. Petrophysics 2004, 48, 490−498.

2779

DOI: 10.1021/acs.energyfuels.6b00035 Energy Fuels 2016, 30, 2771−2779