RESUMEN
Immiscible viscous fingering in porous media occurs when a high viscosity fluid is displaced by an immiscible low viscosity fluid. This paper extends a recent development in the modelling of immiscible viscous fingering to directly simulate experimental floods where the viscosity of the aqueous displacing fluid was increased (by the addition of aqueous polymer) after a period of low viscosity water injection. This is referred to as tertiary polymer flooding, and the objective of this process is to increase the displacement of oil from the system. Experimental results from the literature showed the very surprising observation that the tertiary injection of a modest polymer viscosity could give astonishingly high incremental oil recoveries (IR) of ≥100% even for viscous oils of 7000 mPa.s. This work seeks to both explain and predict these results using recent modelling developments. For the 4 cases (µo/µw of 474 to 7000) simulated in this paper, finger patterns are in line with those observed using X-ray imaging of the sandstone slab floods. In particular, the formation of an oil bank on tertiary polymer injection is very well reproduced and the incremental oil response and water cut drops induced by the polymer are very well predicted. The simulations strongly support our earlier claim that this increase in incremental oil displacement cannot be explained solely by a viscous "extended Buckley-Leverett" (BL) linear displacement effect; referred to in the literature simply as "mobility control". This large response is the combination of this effect (BL) along with a viscous crossflow (VX) mechanism, with the latter VX effect being the major contributor to the recovery mechanism.
RESUMEN
Several experimental studies have shown significant improvement in heavy oil recovery with polymers displaying different types of rheology, and the effect of rheology has been shown to be important. These experimental studies have been designed to investigate why this is so by applying a constant flow rate and the same polymer effective viscosity at this injection rate. The types of rheology studied vary from Newtonian and shear thinning behavior to complex rheology involving shear thinning and thickening behavior. The core flood experiments show a significantly higher oil recovery with polyacrylamide (HPAM), which exhibits shear thinning/thickening behavior compared to biopolymers like Xanthan, which is purely shear thinning. Various reasons for these observed oil recovery results have been conjectured, but, to date, a clear explanation has not been conclusively established. In this paper, we have investigated the theoretical rationale for these results by using a dynamic pore scale network model (DPNM), which can model imbibition processes (water injection) in porous media and also polymer injection. In the DPNM, the polymer rheology can be shear thinning, shear thinning/thickening, or Newtonian (constant viscosity). Thus, the local effective viscosity in a pore within the DPNM depends on the local shear rate in that pore. The predicted results using this DPNM show that the polymer causes changes in the local flow velocity field, which, as might be expected, are different for different rheological models, and the changes in the velocity profile led to local diversion of flow. This, in turn, led to different oil recovery levels in imbibition. However, the critical result is that the DPNM modelling shows exactly the same trend as was observed in the experiments, viz. that the shear thinning/thickening polymer gave the highest oil recovery, followed by the Newtonian Case and the purely shear thinning polymer gave the lowest recover, but this latter case was still above the waterflood result. The DPNM simulations showed that the shear-thinning/thickening polymer show a stabilized frontal velocity and increased oil mobilization, as observed in the experiments. Simulations for the shear-thinning polymer show that, in high-rate bonds, the average viscosity is greatly reduced, and this causes enhanced water fingering compared to the Newtonian polymer case. No other a priori model of the two-phase fluid physics of imbibition, coupled with the polymer rheology, has achieved this degree of predictive explanation, of these experimental observations, to our knowledge.
RESUMEN
Using X-ray computed microtomography, we have visualized and quantified the in situ structure of a trapped nonwetting phase (oil) in a highly heterogeneous carbonate rock after injecting a wetting phase (brine) at low and high capillary numbers. We imaged the process of capillary desaturation in 3D and demonstrated its impacts on the trapped nonwetting phase cluster size distribution. We have identified a previously unidentified pore-scale event during capillary desaturation. This pore-scale event, described as droplet fragmentation of the nonwetting phase, occurs in larger pores. It increases volumetric production of the nonwetting phase after capillary trapping and enlarges the fluid-fluid interface, which can enhance mass transfer between the phases. Droplet fragmentation therefore has implications for a range of multiphase flow processes in natural and engineered porous media with complex heterogeneous pore spaces.
RESUMEN
The in situ rheology of polymeric solutions has been studied experimentally in etched silicon micromodels which are idealizations of porous media. The rectangular channels in these etched networks have dimensions typical of pore sizes in sandstone rocks. Pressure drop/flow rate relations have been measured for water and non-Newtonian hydrolyzed-polyacrylamide (HPAM) solutions in both individual straight rectangular capillaries and in networks of such capillaries. Results from these experiments have been analyzed using pore-scale network modeling incorporating the non-Newtonian fluid mechanics of a Carreau fluid. Quantitative agreement is seen between the experiments and the network calculations in the Newtonian and shear-thinning flow regions demonstrating that the 'shift factor,'alpha, can be calculated a priori. Shear-thickening behavior was observed at higher flow rates in the micromodel experiments as a result of elastic effects becoming important and this remains to be incorporated in the network model.
