Abstract
A mathematical model describing one-dimensional (1D), isothermal flow of a ternary, two-phase surfactant system in isotropic porous media is presented along with numerical solutions of special cases. These solutions exhibit oil recovery profiles similar to those observed in laboratory tests of oil displacement by surfactant systems in cores. The model includes the effects of surfactant transfer between aqueous and hydrocarbon phases and both reversible and irreversible surfactant adsorption by the porous medium. The effects of capillary pressure and diffusion are ignored, however. The model is based on relative permeability concepts and employs a family of relative permeability curves that incorporate the effects of surfactant concentration on interfacial tension (IFT), the viscosity of the phases, and the volumetric flow rate. A numerical procedure was developed that results in two finite difference equations that are accurate to second order in the timestep size and first order in the spacestep size and allows explicit calculation of phase saturations and surfactant concentrations as a function of space and time variables. Numerical dispersion (truncation error) present in the two equations tends to mimic the neglected present in the two equations tends to mimic the neglected effects of capillary pressure and diffusion. The effective diffusion constants associated with this effect are proportional to the spacestep size. proportional to the spacestep size.
Introduction
In a previous paper we presented a system of differential equations that can be used to model oil recovery by chemical flooding. The general system allows for an arbitrary number of components as well as an arbitrary number of phases in an isothermal system. For a binary, two-phase system, the equations reduced to those of the Buckley-Leverett theory under the usual assumptions of incompressibility and each phase containing only a single component, as well as in the more general case where both phases have significant concentrations of both components, but the phases are incompressible and the concentration in one phase is a very weak function of the pressure of the other phase at a given temperature. pressure of the other phase at a given temperature. For a ternary, two-phase system a set of three differential equations was obtained. These equations are applicable to chemical flooding with surfactant, polymer, etc. In this paper, we present a numerical solution to these equations paper, we present a numerical solution to these equations for I D flow in the absence of gravity. Our purpose is to develop a model that includes the physical phenomena influencing oil displacement by surfactant systems and bridges the gap between laboratory displacement tests and reservoir simulation. It also should be of value in defining experiments to elucidate the mechanisms involved in oil displacement by surfactant systems and ultimately reduce the number of experiments necessary to optimize a given surfactant system.