scholarly journals Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

Author(s):  
Peter Mora ◽  
Gabriele Morra ◽  
Dave A. Yuen ◽  
Ruben Juanes

AbstractWe present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ∈ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ∈ [ 180 ∘ , 0 ∘ ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

2021 ◽  
Author(s):  
Peter Mora ◽  
Gabriele Morra ◽  
Dave Yuen ◽  
Ruben Juanes

Abstract We present a suite of numerical simulations of two-phase flow through a 2D model of a porous medium using the Rothman-Keller Lattice Boltzmann Method to study the effect of viscous fingering on the recovery factor as a function of viscosity ratio and wetting angle. This suite involves simulations spanning wetting angles from non-wetting to perfectly wetting and viscosity ratios spanning from 0.01 through 100. Each simulation is initialized with a porous model that is fully saturated with a "blue" fluid, and a "red" fluid is then injected from the left. The simulation parameters are set such that the capillary number is 10, well above the threshold for viscous fingering, and with a Reynolds number of 0.2 which is well below the transition to turbulence and small enough such that inertial effects are negligible. Each simulation involves the "red" fluid being injected from the left at a constant rate such in accord with the specified capillary number and Reynolds number until the red fluid breaks through the right side of the model. As expected, the dominant effect is the viscosity ratio, with narrow tendrils (viscous fingering) occurring for small viscosity ratios with M ≪ 1, and an almost linear front occurring for viscosity ratios above unity. The wetting angle is found to have a more subtle and complicated role. For low wetting angles (highly wetting injected fluids), the finger morphology is more rounded whereas for high wetting angles, the fingers become narrow. The effect of wettability on saturation (recovery factor) is more complex than the expected increase in recovery factor as the wetting angle is decreased, with specific wetting angles at certain viscosity ratios that optimize yield. This complex phase space landscape with hills, valleys and ridges suggests the dynamics of flow has a complex relationship with the geometry of the medium and hydrodynamical parameters, and hence recovery factors. This kind of behavior potentially has immense significance to Enhanced Oil Recovery (EOR). For the case of low viscosity ratio, the flow after breakthrough is localized mainly through narrow fingers but these evolve and broaden and the saturation continues to increase albeit at a reduced rate. For this reason, the recovery factor continues to increase after breakthrough and approaches over 90% after 10 times the breakthrough time.


Author(s):  
Peter Mora ◽  
Gabriele Morra ◽  
Dave A. Yuen ◽  
Ruben Juanes

AbstractWe conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios $$M = \nu _\mathrm{r}/\nu _\mathrm{b} \in [0.001,100]$$ M = ν r / ν b ∈ [ 0.001 , 100 ] and wetting angles $$\theta _\mathrm{w} \in [ 0^\circ ,180^\circ ]$$ θ w ∈ [ 0 ∘ , 180 ∘ ] . As expected, at low-viscosity ratios $$M=\nu _\mathrm{r}/\nu _\mathrm{b} \ll 1$$ M = ν r / ν b ≪ 1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios $$M \gg 1$$ M ≫ 1 , we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.


1984 ◽  
Vol 24 (03) ◽  
pp. 325-327 ◽  
Author(s):  
L. Paterson ◽  
V. Hornof ◽  
G. Neale

Abstract This paper discusses the viscous fingering that occurs when water or a surfactant solution displaces oil in a porous medium. Such floods were visualized in an porous medium. Such floods were visualized in an oil-wet porous medium composed of fused plastic particles. The flow structure changed significantly within the range of capillary numbers between 10 -4 and 10 -3 . The addition of surfactant resulted in narrower fingers, which developed in a more dispersive fashion. Introduction In describing fluid/fluid displacements in porous media, a useful dimensionless quantity is the capillary number, (1) which corresponds to the ratio of viscous forces to capillary forces. Here, v is the specific fluid discharge or Darcy velocity, it is viscosity, and o is interfacial tension (IFT). It has been shown that the recovery of oil from an underground reservoir increases significantly if the capillary number can be increased beyond the range of 1 × 10 -4 to 2 × 10 -3 during water flooding (see Larson et al. 1 ). To this end, surfactants are used extensively in tertiary oil recovery operations with the objective of reducing IFT and consequently mobilizing the oil ganglia which otherwise would remain trapped. This paper is concerned with the viscous fingering that occurs when water displaces oil in a porous medium, and we present a brief consideration on the effects that surfactants have on fingering. Noting that Peters and Flock have visualized fingering within the range of capillary numbers between 1.6 × 10 -6 and 7.2 × 10 -4, we present here visualizations at capillary numbers of 7.7 × 10 5 and 1.0 × 10 -3. Both our visualizations and the experiments of Peters and Flock involve large viscosity ratios so that only the viscosity of the more viscous fluid is considered when determining the capillary number. In particular, it is observed that as the capillary number increases, ganglia or blobs of displacing fluid are created at the displacement front in correspondence with the capillary numbers at which trapped ganglia are mobilized. This creation of ganglia at capillary numbers above 10 -3 was noted briefly in a previous paper 3 in which heptane displacing glycerine previous paper 3 in which heptane displacing glycerine was described. A secondary objective of this work was to test the Chuoke et al. theory for predicting the wavelength of fingers, wavelength being the peak-to-peak distance between adjacent well-developed fingers. Experimental Procedure The apparatus for these studies was described in Ref. 3. Basically, it consists of a slab of consolidated plastic particles 1.34 × 0.79 × 0.0 1 8 ft [0.44 × 0.26 × 0.006 m] with particles 1.34 × 0.79 × 0.0 1 8 ft [0.44 × 0.26 × 0.006 m] with a porosity of 0.43 and a permeability of 7, 100 darcies. This high permeability, when compared with that of reservoir rocks, should not be important for this study since capillary numbers and residual saturations are independent of pore size. Water (viscosity 1 cp [1 mPa s]) was used to displace paraffin oil (viscosity 68 cp 168 mPa s] at 77F [25C]). The water was dyed with methylene blue (which acts as a mild surfactant). Without the dye, the oil/water IFT was 42 dyne/cm [42 mN/m]. The addition of dye lowered this value to 36 dyne/cm [36 mN/m] for the concentration of dye used. For the surfactant flood, a 1 % sodium alkyl aryl sulfonate solution was used, giving a surfactant-solution/paraffin-oil IFT of 3.0 dyne/cm [3.0 mN/m]. Water Displacing Oil To compare our experiments with previous investigations of fingering, the displacement of paraffin oil by water at an average specific fluid discharge of 1.34 × 10–4 ft/sec [4.1 × 10 -5 m/s], corresponding to a capillary number of 7.7 × 10 -5, was studied (Fig. 1). Chuoke et al .4 and later Peters and Flock 2 have presented a formula for predicting the wavelength of presented a formula for predicting the wavelength of finger, lambda m : (2) where k is permeability, C is a dimensionless parameter which Peters and Flock call the wettability number and suggest would have the value 25 for an oil-wet porous medium, and mu o and mu ware viscosities of the displaced oil and displacing water, respectively. It was observed that the plastic particles of the porous medium considered here were oil wet because of adsorption of oil. SPEJ P. 325


Author(s):  
Mehrdad Sepehri ◽  
Babak Moradi ◽  
Abolghasem Emamzadeh ◽  
Amir H. Mohammadi

Nowadays, nanotechnology has become a very attractive subject in Enhanced Oil Recovery (EOR) researches. In the current study, a carbonate system has been selected and first the effects of nanoparticles on the rock and fluid properties have been experimentally investigated and then the simulation and numerical modeling of the nanofluid injection for enhanced oil recovery process have been studied. After nanofluid treatment, experimental results have shown wettability alteration. A two-phase flow mathematical model and a numerical simulator considering wettability alteration have been developed. The numerical simulation results show that wettability alteration from oil-wet to water-wet due to presence of nanoparticles can lead to 8–10% increase in recovery factor in comparison with normal water flooding. Different sensitivity analyses and injection scenarios have been considered and assessed. Using numerical modeling, wettability alteration process and formation damage caused by entrainment and entrapment of nanoparticles in porous media have been proved. Finally, the net rate of nanoparticles’ loss in porous media has been investigated.


Author(s):  
Akshay C. Gunde ◽  
Sushanta K. Mitra ◽  
Tayfun Babadagli

Study of flow through porous media has been an area of major interest due to its application in diverse areas like Enhanced Oil Recovery. In order to gain a better understanding of the physical processes taking place inside a porous structure, a large number of attempts have been made to computationally simulate multiphase fluid flow at pore-scale. Recently, application of Lattice Boltzmann Method has gained popularity for this very purpose, considering its relative superiority in dealing with complex boundaries and multiphase flow. However, in order that such a numerical analysis is successful, a proper understanding of the geometry of the pore structure at the microscale is required. This paper uses a Micro-CT scan image of a Berea Sandstone core, which displays a two dimensional representation of pore network inside the scanned sample. The processed image has been imported and simulation of an immiscible two-phase flow has been carried out by using a Lattice Boltzmann program. The resident fluid (oil) has been displaced by the invading fluid (water) due to application of a pressure gradient. The pore surfaces have been treated as solid boundaries and bounce back scheme has been implemented on them to account for the no-slip condition. The ability of the code to import an arbitrary porous geometry and perform numerical analysis of fluid flow has been demonstrated.


2020 ◽  
Vol 8 ◽  
Author(s):  
Xuemei Wei ◽  
Wenchao Jiang ◽  
Yanyu Zhang ◽  
Zhao Wang ◽  
Xiaojun Li ◽  
...  

Clay minerals are usually regarded as an important factor affecting the results of low salinity water (LSW) flooding. However, experiments on clay minerals are mainly in qualitative stage, the mechanism of clay minerals has not been studied completely. In this paper, Zeta potential of four kinds of clay minerals (montmorillonite; chlorite; illite; kaolinite) in different brine was measured, microscopic models of these clay minerals were made to measured wetting angle in different brine, and montmorillonite and kaolinite were chosen to conduct microscopic displacement experiments through customized micro-glass etching models. From experiment results, the following conclusions can be get: 1). With the decrease of salinity of injected water, the negative zeta potential of clay minerals increases and the wetting angle decreases. 2). Clay minerals are more sensitive to monovalent Na+ than bivalent Ca2+. 3). The results of microscopic experiments show that LSW can effectively improve oil recovery, whether kaolinite or montmorillonite. The recovery of montmorillonite is better with a relatively high salinity of LSW and kaolinite is better with a relatively low salinity of LSW. The mechanism of LSW improves kaolinite recovery factor is the change of wettability while that of montmorillonite is the increase of water phase wettability. However, a lot of droplet-like residual oil cannot be displaced in the montmorillonite throat. In filed production, both kaolinite-rich and montmorillonite-rich reservoirs are suitable for LSW flooding to improve oil recovery. However, for kaolinite reservoirs, a lower salinity of injected water would produce a better result, while for montmorillonite reservoirs, residual oil droplets in the throat are noteworthy.


2021 ◽  
Author(s):  
Rumbidzai Nhunduru ◽  
Omid Shahrokhi ◽  
Krystian Wlodarczyk ◽  
Amir Jahanbakhsh ◽  
Susana Garcia ◽  
...  

<p>Immiscible fluid displacement and the trapping of residual oil and gas phases in the pore spaces of reservoir rocks is critical to geological operations such as carbon geo-sequestration and enhanced oil recovery. In carbon geo-sequestration, residual trapping is advantageous because it ensures long-term storage security of carbon dioxide (CO<sub>2</sub>). In contrast, residual trapping can pose significant challenges during waterflooding in oil recovery operations where large volumes of oil may remain trapped in the interstitial spaces of the porous reservoir rock and cannot be extracted, thereby reducing the efficiency of the recovery process. In such operations, residual trapping is strongly influenced by the inherent surface roughness of the solid rock matrix amongst many factors. Surface roughness occurs in natural reservoir rocks as a result of geological processes that physically, chemically or biologically convert sediments into sedimentary rock (known as diagenesis) and weathering.</p><p>The effects of surface roughness on immiscible two-phase flow are currently not well understood. Previous investigations into residual trapping in porous media have mainly focused on the influence of factors such as pore geometry, wettability, fluids interfacial tension, mobility ratio and injection scenarios. Although some of these studies acknowledge the potential effect of surface roughness, there is still a lack of quantitative characterization and understanding of the influence of surface roughness on immiscible two-phase displacements in porous media.</p><p>In this study, the impacts of surface roughness on immiscible two-phase displacement are quantified. Immiscible two-phase displacement of air by water was conducted in a custom laser-manufactured glass microfluidic chip (micromodel). The glass chip comprised a 2.5D micro-structure analogous to the pore network pattern (micro-structure) of a natural reservoir rock, Oolitic limestone. The pore network pattern consisted of cylindrical pillars 400 µm in diameter arranged in a rhombohedra type of packing, generated on to a glass substrate using an ultrafast, pulsed picosecond laser. Surface roughness is an innate characteristic of laser machined surfaces and as a result, small variations in depth of the porous micro-structure were observed (50 ± 8 µm). The average surface roughness (S<sub>a</sub>) of the laser-machined structure was measured to be 1.2 μm.</p><p>Experimental results for the rough micromodel exhibit high repeatability of fluid displacement patterns (preferential flow pathways) demonstrating that surface roughness has a strong influence on fluid invasion patterns and sweep efficiency and its effects must not be ignored. To ascertain the effects of surface roughness on the fluid displacement process, a direct numerical simulation (DNS) of the fluid displacement process was performed in OpenFoam using the Volume of Fluid (VOF) method assuming zero surface roughness. Comparing the experimental results with the numerical simulations, we show that surface roughness can significantly enhance residual trapping in porous media by up to 49.2%.</p><p> </p>


RSC Advances ◽  
2015 ◽  
Vol 5 (104) ◽  
pp. 85373-85382 ◽  
Author(s):  
Mingming Lv ◽  
Shuzhong Wang

The pore-scale behaviors of hot water displacement in a pore–throat microchannel were revealed by simulations for different wettability systems.


2017 ◽  
Vol 824 ◽  
pp. 550-573 ◽  
Author(s):  
Ioannis Zacharoudiou ◽  
Emily M. Chapman ◽  
Edo S. Boek ◽  
John P. Crawshaw

The aim of this work is to better understand fluid displacement mechanisms at the pore scale in relation to capillary-filling rules. Using specifically designed micro-models we investigate the role of pore body shape on fluid displacement during drainage and imbibition via quasi-static and spontaneous experiments at ambient conditions. The experimental results are directly compared to lattice Boltzmann (LB) simulations. The critical pore-filling pressures for the quasi-static experiments agree well with those predicted by the Young–Laplace equation and follow the expected filling events. However, the spontaneous imbibition experimental results differ from those predicted by the Young–Laplace equation; instead of entering the narrowest available downstream throat the wetting phase enters an adjacent throat first. Thus, pore geometry plays a vital role as it becomes the main deciding factor in the displacement pathways. Current pore network models used to predict displacement at the field scale may need to be revised as they currently use the filling rules proposed by Lenormandet al.(J. Fluid Mech., vol. 135, 1983, pp. 337–353). Energy balance arguments are particularly insightful in understanding the aspects affecting capillary-filling rules. Moreover, simulation results on spontaneous imbibition, in excellent agreement with theoretical predictions, reveal that the capillary number itself is not sufficient to characterise the two phase flow. The Ohnesorge number, which gives the relative importance of viscous forces over inertial and capillary forces, is required to fully describe the fluid flow, along with the viscosity ratio.


Sign in / Sign up

Export Citation Format

Share Document