Foam injection is a promising enhanced-oil-recovery (EOR) technology that significantly improves the sweep efficiency of gas injection. Simulation of foam/oil displacement in reservoirs is an expensive process for conventional simulation because of the strongly nonlinear physics, such as multiphase flow and transport with oil/foam interactions. In this work, an operator-based linearization (OBL) approach, combined with the representation of foam by an implicit-texture (IT) model with two flow regimes, is extended for the simulation of the foam EOR process. The OBL approach improves the efficiency of the highly nonlinear foam-simulation problem by transforming the discretized nonlinear conservation equations into a quasilinear form using state-dependent operators. The state-dependent operators are approximated by discrete representation on a uniform mesh in parameter space. The numerical-simulation results are validated by using three-phase fractional-flow theory for foam/oil flow. Starting with an initial guess depending on the fitting of steady-state experimental data with oil, the OBL foam model is regressed to experimental observations using a gradient-optimization technique. A series of numerical validation studies is performed to investigate the accuracy of the proposed approach. The numerical model shows good agreement with analytical solutions at different conditions and with different foam parameters. With finer grids, the resolution of the simulation is better, but at the cost of more expensive computations. The foam-quality scan is accurately fitted to steady-state experimental data, except in the low-quality regime. In this regime, the used IT foam model cannot capture the upward-tilting pressure gradient (or apparent viscosity) contours. 1D and 3D simulation results clearly demonstrate two stages of foam propagation from inlet to outlet, as seen in the computed-tomography (CT) coreflood experiments: weak foam displaces most of the oil, followed by a propagation of stronger foam at lower oil saturation. OBL is a direct method to reduce nonlinearity in complex physical problems, which can significantly improve computational performance. Taking its accuracy and efficiency into account, the data-driven OBL-based approach could serve as a platform for efficient numerical upscaling to field-scale applications.

A new capillary number (N _ca ) definition is proposed for 2-D etched micromodels. We derive the new definition from a force balance on a nonwetting ganglion trapped by capillarity. It incorporates the impact of pore microstructure on mobilization. The geometrical factors introduced can be estimated directly from image analysis of the pore network etched in the micromodel, without conducting flow experiments. The improved fit of the new N _ca to published data supports its validity. The new definition yields a consistent trend in the capillary-desaturation curve. The conventional N _ca definitions proposed for porous rock give a large scatter in the capillary-desaturation curve for data in micromodels. This is due to the different type of flow in micromodels, as 2-D networks, relative to 3-D geological porous media. In particular, permeability is dominated by channel depth in micromodels with shallow depth of etching, and generally, there is no simultaneous multiphase flow under capillary-dominated conditions. Applying the conventional definitions to results in micromodels may lead to misleading conclusions for fluid transport in geological formations.

Foam flow in porous media without oil shows two regimes depending on foam quality (gas fractional flow). Complexity and limited data on foam/oil interactions in porous media greatly restrict understanding of foam in contact with oil. Distinguishing which regimes are affected by oil is key to modeling the effect of oil on foam. We report steady-state corefloods to investigate the effect of oil on foam through its effect on the two flow regimes. We fit the parameters of a widely used local-equilibrium (LE) foam model to data for concurrent foam/oil flow. This research provides a practical approach and initial data for simulating foam enhanced oil recovery (EOR) in the presence of oil. To ensure steady state, oil is coinjected with foam at a fixed ratio of oil (U _o ) to water (U _w ) superficial velocities in a Bentheimer Sandstone core. Model oils used here consist of a composition of hexadecane, which is benign to foam stability, and oleic acid (OA), which can destroy foam. Varying the concentration of OA in the model oil allows one to examine the effect of oil composition on steady-state foam flow. Experimental results show that oil affects both high- and low-quality regimes, with the high-quality regime being more sensitive to oil. In particular, oil increases the limiting water saturation (S _w ) in the high-quality regime and also reduces gas-mobility reduction in the low-quality regime. Unevenly spaced !p contours in the high-quality regime suggest either strongly shear-thinning behavior or an increasingly destabilizing effect of oil. In some cases, the pressure gradient (!p) in the low-quality regime decreases with increasing U _w at fixed gas superficial velocity (U _g ), either with or without oil. This might reflect either an effect of oil, if oil is present, or easier flow of bubbles under wetter conditions. Increasing the OA concentration extends the high-quality regime to lower foam qualities, indicating more difficulty in stabilizing foam. Thus, oil composition plays as significant a role as oil saturation (S _o ). A model fit assuming a fixed S _w and including shear thinning in the low-quality regime does not represent the two regimes when the oil effect is strong enough. In such cases, fitting S _w to each !p contour and excluding shear thinning in the low-quality regime yield a better match to these data. The dependency of S _w on S _o is not yet clear because of the absence of oil-saturation data in this study. Furthermore, none of the current foam-simulation models captures the upward-tilting !p contours in the low-quality regime.

Foam flow in porous media without oil shows two regimes, depending on foam quality (gas fractional flow). Complexity and limited data on foam-oil interactions in porous media greatly restrict understanding of foam in contact with oil. Distinguishing which regimes are affected by oil is key to modelling the effect of oil on foam. We report steady-state corefloods to investigate the effect of oil on foam through its effect on the two flow regimes. We fit parameters of the widely used STARS foam model to data for foam-oil concurrent flow. This research provides a practical approach and initial data for simulating foam EOR in the presence of oil. To ensure steady state, oil is co-injected with foam at a fixed ratio of oil (Uo) to water (Uw) superficial velocities in a Bentheimer sandstone core. Model oils used here consist of two components: hexadecane, which is benign to foam stability, and oleic acid, which can destroy foam. Varying the concentration of oleic acid in the model oil allows one to examine the effect of oil composition on steady-state foam flow. Experimental results show that oil impacts both high- and low-quality regimes, with the high-quality regime more vulnerable to oil. In particular, oil increases the limiting water saturation (Sw^∗) in the high-quality regime and also lessens gas mobility reduction in the low-quality regime. The high-quality regime is strongly shear-thinning in the presence of oil. Pressure gradient (p) in the low-quality regime, in some cases, decreases with increasing Uw at fixed gas superficial velocity (Ug), either with or without oil. This may reflect either an effect of oil, if oil is present, or easier flow of bubbles under wetter conditions. Increasing oleic acid concentration extends the high-quality regime to lower foam qualities, indicating more difficulty in stabilizing foam. Thus oil composition plays as significant a role as oil saturation. A model fit assuming a fixed Sw∗ and including shear-thinning in the low-quality regime doesn't represent each regime when the oil effect is strong enough. In such cases, fitting Sw∗ to each p contour and excluding shear-thinning in the low-quality regime yields a better match to data. The dependency of Sw∗ on oil saturation is not yet clear owing to absence of oil-saturation data in this study. Furthermore, none of the current foam simulation models can capture the upward-tilting p contours in the low-quality regime.

Foam flow in porous media without oil shows two regimes, depending on foam quality (gas fractional flow). Complexity and limited data on foam-oil interactions in porous media greatly restrict understanding of foam in contact with oil. Distinguishing which regimes are affected by oil is key to modelling the effect of oil on foam. We report steady-state corefloods to investigate the effect of oil on foam through its effect on the two flow regimes. We fit parameters of the widely used STARS foam model to data for foam-oil concurrent flow. This research provides a practical approach and initial data for simulating foam EOR in the presence of oil. To ensure steady state, oil is co-injected with foam at a fixed ratio of oil (Uo) to water (Uw) superficial velocities in a Bentheimer sandstone core. Model oils used here consist of two components: hexadecane, which is benign to foam stability, and oleic acid, which can destroy foam. Varying the concentration of oleic acid in the model oil allows one to examine the effect of oil composition on steady-state foam flow. Experimental results show that oil impacts both high- and low-quality regimes, with the high-quality regime more vulnerable to oil. In particular, oil increases the limiting water saturation (Sw^∗) in the high-quality regime and also lessens gas mobility reduction in the low-quality regime. The high-quality regime is strongly shear-thinning in the presence of oil. Pressure gradient (p) in the low-quality regime, in some cases, decreases with increasing Uw at fixed gas superficial velocity (Ug), either with or without oil. This may reflect either an effect of oil, if oil is present, or easier flow of bubbles under wetter conditions. Increasing oleic acid concentration extends the high-quality regime to lower foam qualities, indicating more difficulty in stabilizing foam. Thus oil composition plays as significant a role as oil saturation. A model fit assuming a fixed Sw∗ and including shear-thinning in the low-quality regime doesn't represent each regime when the oil effect is strong enough. In such cases, fitting Sw∗ to each p contour and excluding shear-thinning in the low-quality regime yields a better match to data. The dependency of Sw∗ on oil saturation is not yet clear owing to absence of oil-saturation data in this study. Furthermore, none of the current foam simulation models can capture the upward-tilting p contours in the low-quality regime.