Approximate Pore-Level Modeling for Apparent Viscosity of Polymer-Enhanced Foam in Porous Media

SPE Journal ◽  
2008 ◽  
Vol 13 (01) ◽  
pp. 17-25 ◽  
Author(s):  
Chun Huh ◽  
William R. Rossen

Summary Foam is used in the oil industry in a variety of applications, and polymer is sometimes added to increase foam's stability and effectiveness. A variety of surfactant and polymer combinations have been employed to generate polymer-enhanced foam (PEF), typically anionic surfactants and anionic polymers, to reduce their adsorption in reservoir rock. While addition of polymer to bulk foam is known to increase its viscosity and apparent stability, polymer addition to foams for use in porous media has not been as effective. In this pore-level modeling study, we develop an apparent viscosity expression for PEF at fixed bubble size, as a preliminary step to interpret the available laboratory coreflood data. To derive the apparent viscosity, the pressure-drop calculation of Hirasaki and Lawson (1985) for gas bubbles in a circular tube is extended to include the effects of shear-thinning polymer in water, employing the Bretherton's asymptotic matching technique. For polymer rheology, the Ellis model is employed, which predicts a limiting Newtonian viscosity at the low-shear limit and the well-known power-law relation at high shear rates. While the pressure drop caused by foam can be characterized fully with only the capillary number for Newtonian liquid, the shear-thinning liquid requires one additional grouping of the Ellis-model parameters and bubble velocity. The model predicts that the apparent viscosity for PEF shows behavior more shear-thinning than that for polymer-free foam, because the polymer solution being displaced by gas bubbles in pores tends to experience a high shear rate. Foam apparent viscosity scales with gas velocity (Ug) with an exponent [-a/(a+2)], where a, the Ellis-model exponent, is greater than 1 for shear-thinning fluids. With a Newtonian fluid, for which a = 1, foam apparent viscosity is proportional to the (-1/3) power of Ug, as derived by Hirasaki and Lawson. A simplified capillary-bundle model study shows that the thin-film flow around a moving foam bubble is generally in the high-shear, power-law regime. Because the flow of polymer solution in narrower, water-filled tubes is also governed by shear-thinning rheology, it affects foam mobility as revealed by plot of pressure gradient as a function of water and gas superficial velocities. The relation between the rheology of the liquid phase and that of the foam is not simple, however. The apparent rheology of the foam depends on the rheology of the liquid, the trapping and mobilization of gas as a function of pressure gradient, and capillary pressure, which affects the apparent viscosity of the flowing gas even at fixed bubble size. Introduction When a gas such as CO2 or N2 is injected into a mature oil reservoir for improved oil recovery, its sweep efficiency is usually very poor because of gravity segregation, reservoir heterogeneity, and viscous fingering of gas, and foam is employed to improve sweep efficiency with better mobility control (Shi and Rossen 1998; Zeilinger et al. 1996). When oil is produced from a thin oil reservoir overlain with a gas zone, a rapid coning of gas can drastically reduce oil production rate, and foam is used to delay the gas coning (Aarra et al. 1997; Chukwueke et al. 1998; Dalland and Hanssen 1997; Thach et al. 1996). During a well stimulation operation with acid, a selective placement of acid into a low-permeability zone from which oil has not been swept is desired, which can be accomplished with use of foam (Cheng et al. 2002). For environmental remediation of subsurface soil using surfactant, foam is used to improve displacement of contaminant, such as DNAPL, from heterogeneous soil (Mamun et al. 2002).

SPE Journal ◽  
2019 ◽  
Vol 24 (06) ◽  
pp. 2793-2803 ◽  
Author(s):  
Arthur Uno Rognmo ◽  
Sunniva Brudvik Fredriksen ◽  
Zachary Paul Alcorn ◽  
Mohan Sharma ◽  
Tore Føyen ◽  
...  

Summary This paper presents an ongoing CO2–foam upscaling research project that aims to advance CO2–foam technology for accelerating and increasing oil recovery, while reducing operational costs and lessening the carbon footprint left during CO2 enhanced oil recovery (EOR). Laboratory CO2–foam behavior was upscaled to pilot scale in an onshore carbonate reservoir in Texas, USA. Important CO2–foam properties, such as local foam generation, bubble texture, apparent viscosity, and shear–thinning behavior with a nonionic surfactant, were evaluated using pore–to–core upscaling to develop accurate numerical tools for a field–pilot prediction of increased sweep efficiency and CO2 utilization. At pore–scale, high–pressure silicon–wafer micromodels showed in–situ foam generation and stable liquid films over time during no–flow conditions. Intrapore foam bubbles corroborated high apparent foam viscosities measured at core scale. CO2–foam apparent viscosity was measured at different rates (foam–rate scans) and different gas fractions (foam–quality scans) at core scale. The highest mobility reduction (foam apparent viscosity) was observed between 0.60 and 0.70 gas fractions. The maximum foam apparent viscosity was 44.3 (±0.5) mPa·s, 600 times higher than that of pure CO2, compared with the baseline viscosity (reference case, without surfactant), which was 1.7 (±0.6) mPa·s, measured at identical conditions. The CO2–foam showed shear–thinning behavior with approximately 50% reduction in apparent viscosity when the superficial velocity was increased from 1 to 8 ft/D. Strong foam was generated in EOR corefloods at a gas fraction of 0.70, resulting in an apparent viscosity of 39.1 mPa·s. Foam parameters derived from core–scale foam floods were used for numerical upscaling and field–pilot performance assessment.


Author(s):  
Ram P. Bharti ◽  
Dalton J. E. Harvie ◽  
Malcolm R. Davidson

Pressure drop and electroviscous effects in the axisymmetric, steady, fully developed, pressure-driven flow of incompressible power-law fluids through a cylindrical microchannel at low Reynolds number (Re = 0.01) have been investigated. The Poisson-Boltzmann equation (describing the electrical potential) and the momentum equations in conjunction with electrical force and power-law fluid rheology have been solved numerically using the finite difference method. The pipe wall is considered to have uniform surface charge density (S = 4) and the liquid is assumed to be a symmetric electrolyte solution. In particular, the influence of the dimensionless inverse Debye length (K = 2, 20) and power-law flow behaviour index (n = 0.2, 1, 1.8) on the EDL potential, ion concentrations and charge density profiles, induced electrical field strength, velocity and viscosity profiles and pressure drop have been studied. As expected, the local EDL potential, local charge density and electrical field strength increases with decreasing K and/or increasing S. The velocity profiles cross-over away from the charged pipe wall with increasing K and/or decreasing n. The maximum velocity at the center of the pipe increases with increasing n and/or increasing S and/or decreasing K. The shear-thinning fluid viscosity is strongly dependent on K and S, whereas the shear-thickening viscosity is very weakly dependent on K and S. For fixed K, as the fluid behaviour changes from Newtonian (n = 1) to shear-thinning (n < 1), the induced electrical field strength increases and maximum velocity reduces. On the other hand, the change in fluid behaviour from Newtonian (n = 1) to shear-thickening (n > 1) decreases the electrical field strength and increases the maximum velocity. The non-Newtonian effects on maximum velocity and pressure drop are stronger in shear-thinning fluids at small K and large S, the shear-thickening fluids show opposite influence. Electroviscous effects enhance with decreasing K and/or increasing S. The electroviscous effects show complex dependence on the non-Newtonian tendency of the fluids. The shear-thickening (n > 1) fluids and/or smaller K show stronger influence on the pressure drop and thus, enhance the electroviscous effects than that in shear-thinning (n < 1) fluids and/or large K where EDL is very thin.


2012 ◽  
Vol 170-173 ◽  
pp. 2601-2608
Author(s):  
Xue Li Xia ◽  
Hong Fu Qiang ◽  
Wang Guang

To evaluate the effect of a converging injector geometry, volumetric flow rate and gallant content on the pressure drop, the velocity and viscosity fields, the governing equations of the steady, incompressible, isothermal, laminar flow of a Power-Law, shear-thinning gel propellant in a converging injector were formulated, discretized and solved. A SIMPLEC numerical algorithm was applied for the solution of the flow field. The results indicate that the mean apparent viscosity decreases with increasing the volumetric flow rate and increasing the gallant content results in an increase in the viscosity. The results indicate also that the convergence angle can produce additional decrease in the mean apparent viscosity of the fluid. The mean apparent viscosity decreases significantly with increasing the convergence angle of the injector, and its value is limited by the Newtonian viscosity η∞. The effect of the convergence angle on the mean apparent viscosity is more significant than the effect of the volumetric flow rate and the gallant content on the mean apparent viscosity. Additional decreasing the viscosity results in increasing the pressure drop with increasing convergence angle. It is important to injector design that the viscosity decreasing and the pressure drop increasing are took into account together.


2018 ◽  
Vol 838 ◽  
pp. 573-605 ◽  
Author(s):  
Fatima-Ezzahra Moukhtari ◽  
Brice Lecampion

We use the Carreau rheological model which properly accounts for the shear-thinning behaviour between the low and high shear rate Newtonian limits to investigate the problem of a semi-infinite hydraulic fracture propagating at a constant velocity in an impermeable linearly elastic material. We show that the solution depends on four dimensionless parameters: a dimensionless toughness (function of the fracture velocity, confining stress, material and fluid parameters), a dimensionless transition shear stress (related to both fluid and material behaviour), the fluid shear-thinning index and the ratio between the high and low shear rate viscosities. We solve the complete problem numerically combining a Gauss–Chebyshev method for the discretization of the elasticity equation, the quasi-static fracture propagation condition and a finite difference scheme for the width-averaged lubrication flow. The solution exhibits a complex structure with up to four distinct asymptotic regions as one moves away from the fracture tip: a region governed by the classical linear elastic fracture mechanics behaviour near the tip, a high shear rate viscosity asymptotic and power-law asymptotic region in the intermediate field and a low shear rate viscosity asymptotic far away from the fracture tip. The occurrence and order of magnitude of the extent of these different viscous asymptotic regions are estimated analytically. Our results also quantify how shear thinning drastically reduces the size of the fluid lag compared to a Newtonian fluid. We also investigate simpler rheological models (power law, Ellis) and establish the small domain where they can properly reproduce the response obtained with the complete rheology.


2011 ◽  
Vol 45 (8) ◽  
pp. 883-893 ◽  
Author(s):  
Jun Feng Wang ◽  
Wook Ryol Hwang

This work is an extension of our previous work [Wang JF and Hwang WR. Permeability prediction of fibrous porous media in a bi-periodic domain. J Compos Mater 2007; 42: 909—929], in which a finite-element fictitious-domain mortar-element technique was developed to investigate the permeability of fibrous porous media in the bi-periodic domain, to non-Newtonian shear-thinning fluid. Considering the amount of shear-thinning, the pressure drop, the fiber microstructure, and the porosity as parameters, we investigate (i) the (normalized) mobility and its dependence on both the amount of shear-thinning and the given pressure drop; (ii) mechanisms leading to the main flow path in a highly shear-thinning fluid in randomly distributed fiber problems, and (iii) inter-tow and intra-tow non-Newtonian flow characteristics in a fiber bundle problem. The dependence of the mobility on shear-thinning has been found to appear completely opposite according to given pressure drop values.


2021 ◽  
Vol 3 ◽  
Author(s):  
Subhadeep Roy ◽  
Santanu Sinha ◽  
Alex Hansen

Immiscible two-phase flow of Newtonian fluids in porous media exhibits a power law relationship between flow rate and pressure drop when the pressure drop is such that the viscous forces compete with the capillary forces. When the pressure drop is large enough for the viscous forces to dominate, there is a crossover to a linear relation between flow rate and pressure drop. Different values for the exponent relating the flow rate and pressure drop in the regime where the two forces compete have been reported in different experimental and numerical studies. We investigate the power law and its exponent in immiscible steady-state two-phase flow for different pore size distributions. We measure the values of the exponent and the crossover pressure drop for different fluid saturations while changing the shape and the span of the distribution. We consider two approaches, analytical calculations using a capillary bundle model and numerical simulations using dynamic pore-network modeling. In case of the capillary bundle when the pores do not interact to each other, we find that the exponent is always equal to 3/2 irrespective of the distribution type. For the dynamical pore network model on the other hand, the exponent varies continuously within a range when changing the shape of the distribution whereas the width of the distribution controls the crossover point.


SPE Journal ◽  
2020 ◽  
Vol 25 (02) ◽  
pp. 883-894 ◽  
Author(s):  
Chen Qian ◽  
Ali Telmadarreie ◽  
Mingzhe Dong ◽  
Steven Bryant

Summary The major challenge in enhanced oil recovery (EOR) by gas injection is poor volumetric sweep efficiency, mainly due to the high gas mobility and reservoir heterogeneity. Injecting gas as a foam increases sweep efficiency, but maintaining foam stability within the reservoir remains a challenge. This research evaluates the synergistic interaction of one type of nanoparticle and a surfactant to increase foam stability, using the concentration ratio of the two components to tune the affinity of the nanoparticle for the gas/liquid interface. We test the capability of the synergistic two-component system to stabilize methane foam and compare it with foam stabilized by surfactants only. A key distinction is the foam stability upon contact with oil, and we explain the observations in static and dynamic conditions. Foam stability was measured in both static (bubble structure) and dynamic (flow through porous media) conditions. In the static test, foam is generated by the shaking method, and foam texture (bubble size and shape) and the decay of foam height with time are indicators of foam stability. To test static stability in the presence of oil, heavy oil is injected into the foam/liquid interface. In dynamic tests, foam is pregenerated before flowing at elevated pressures into sandpacks containing various oil saturations. Normalized pressure gradient and apparent viscosity are the indicators of foam stability and effectiveness for improving oil recovery. The extent to which nanoparticles are covered with surfactant governs the foam stability, in both static and dynamic conditions. Static foam is stable in the presence of oil only if the nanoparticles are partially covered by the surfactant. In the dynamic test, foam stabilized with only the surfactant collapses in the porous media when oil is present. Nanoparticles alone could not generate foam regardless of the presence of oil or salinity, but foam stabilized with nanoparticles partially covered by surfactant is stable in the presence of both residual and initial oil, and foam apparent viscosity could reach up to 400 cp at residual heavy oil condition. In both static and dynamic conditions, nanoparticles completely covered with a bilayer of surfactant do not stabilize foam in the presence of oil. Partially covered nanoparticle foam also demonstrated salt tolerance in both static and dynamic tests, and foam apparent viscosity can reach up to 200 cp with high salinity and residual heavy oil presented. Thus, at appropriate surface coverage, the combination of nanoparticles and surfactant is more effective than either stabilizer alone. The result shows that interaction of surfactant and nanoparticles is important in foam stability in the porous media with oil. In particular, this interaction is synergistic at certain coverage. This type of synergy can provide much more robust mobility control for EOR processes involving gas injection.


2021 ◽  
Author(s):  
Ali Telmadarreie ◽  
Christopher Johnsen ◽  
Steven L. Bryant

Abstract This study designs a novel complex fluid (foam/emulsion) using as main components gas, low-toxicity solvents (green solvents) which may promote oil mobilization, and synergistic foam stabilizers (i.e. nanoparticles and surfactants) to improve sweep efficiency. This nanoparticle-enabled green solvent foam (NGS-foam) avoids major greenhouse gas emissions from the thermal recovery process and improves the performance of conventional green solvent-based methods (non-thermal) by increasing the sweep efficiency, utilizing less solvent while producing more oil. Surfactants and nanoparticles were screened in static tests to generate foam in the presence of a water-soluble/oil-soluble solvent and heavy crude oil from a Canadian oil field (1600 cp). The liquid phase of NGS-foam contains surfactant, nanoparticle, and green solvent (GS) all dispersed in the water phase. Nitrogen was used as the gas phase. Fluid flow experiments in porous media with heterogeneous permeability structure mimicking natural environments were performed to demonstrate the dynamic stability of the NGS-foam for heavy oil recovery. The propagation of the pre-generated foam was monitored at 10 cm intervals over the length of porous media (40 cm). Apparent viscosity, pressure gradient, inline measurement of effluent density, and oil recovery were recorded/calculated to evaluate the NGS-foam performance. The outcomes of static experiments revealed that surfactant alone cannot stabilize the green solvent foam and the presence of carefully chosen nanoparticles is crucial to have stable foam in the presence of heavy oil. The results of NGS-foam flow in heterogeneous porous media demonstrated a step-change improvement in oil production such that more than 60% of residual heavy oil was recovered after initial waterflood. This value of residual oil recovery was significantly higher than other scenarios tested in this study (i.e. GS- water and gas co-injection, conventional foam without GS, GS-foam stabilized with surfactant only and GS-waterflood). The increased production occurred because NGS-foam remained stable in the flowing condition, improves the sweep efficiency and increases the contact area of the solvent with oil. The latter factor is significant: comparing to GS-waterflood, NGS-foam produces a unit volume of oil faster with less solvent and up to 80% less water. Consequently, the cost of solvent per barrel of incremental oil will be lower than for previously described solvent applications. In addition, due to its water solubility, the solvent can be readily recovered from the reservoir by post flush of water and thus re-used. The NGS-foam has several potential applications: recovery from post-CHOPS reservoirs (controlling mobility in wormholes and improving the sweep efficiency while reducing oil viscosity), fracturing fluid (high apparent viscosity to carry proppant and solvent to promote hydrocarbon recovery from matrix while minimizing water invasion), and thermal oil recovery (hot NGS-foam for efficient oil viscosity reduction and sweep efficiency improvement).


2011 ◽  
Vol 14 (9) ◽  
pp. 761-776 ◽  
Author(s):  
Hamid Emami Meybodi ◽  
Riyaz Kharrat ◽  
Benyamin Yadali Jamaloei

Sign in / Sign up

Export Citation Format

Share Document