Secondary Recovery - Mathematical Description of Detergent Flooding in Oil Reservoirs

Physically absurd, triple-valued saturations appear in the straight-forward solution of the Buckley-Leverell equations for the displacement of oil by water or gas. From an engineering viewpoint, the triple value causen no difficulty. It is well known how it may be compensnted in order to obtain physically meaningful, numerical results. From a scientific viewpoint, the question still arises: What did the triple value mean? This paper explains how and why the triple value arose in non-capillary Buckley-Leverett theory. The discussion should serve as a background for the understanding and use of the modern method of characteristics in di.~pkzc(>ment theory. INTRODUCTION Displacement theory was introduced to petroleum technologists in 1 941 by Buckley and Leverett.' Their first equation, which they wrote down without derivation, was equivalent to the following: ?s/?t= - qr/F ?f/?r (1) where S = saturation of displacing fluid (a fraction) r = time q,. = total volumetric velocity of both the displacing phase and the phase being displaced (usually, but not necessarily, assumed to be a constant) F = porosity (a fraction) f - fraction of the total volumetric velocity that is the volumetric velocity of the displacing phase (assumed to he a function of S only x = distance. Buckley and Leverett called Eq. I "a material balance equation". Actually it is derived from both a continuity equation, or material balance equation, and Darcy's law. Buckley and Leverett jumped from their Eq. 1 immediately to their Eq. 2, which was the equivalent of the following: (?x/?t)x = qr/F df/ds = qr/F f'(s) (2) The jump from Eq. 1 to Eq. 2 involves apparently simple mathematics, but in that apparently simple mathematics lies a subtle point that is part of the key to the meaning of the triple value to be discussed here. Eq. 2 may be readily interpreted as saying that a point of constant saturation (on a saturation-vs-dis-tance curve) moves with a constant velocity that is proportional to the total volumetric rate, inversely proportional to the porosity, and is otherwise a function of the saturation itself. So that if one knows the derivativc with respect to saturation of the fractional flow function, f(S), for each saturation, one knows the velocity of each point of the moving saturation-vs-distance curve. The function f(S) can he determined experimentally and its derivative f'(S) can be calculated. It turns out that the experimentally derived f'(S) is not a monotonic function of S, but instead it has a definite maximum and declines away at both high and low saturations. This in turn means, in accordance with Eq. 2, that on a moving saturation-vs-distance curve, ccrtain intermediate saturations will travel faster than the saturations either higher or lower. The result is indicated in Fig. 1. In Fig. 1, the abscissas represent distances along a column of porous medium, which column is assumed to be uniform in cross-section perpendicular to the ab-scissal direction. The ordinates represent fractions of pore space occupied by the displacing phase, which may be either water or gas. It is assumed that the saturation is uniform over all planes perpendicular to the abscissal direction. The saturation is a function only of time and one space variable. Time variation of the saturation is represented by change in shape of a curve such as curve AGJ, which represents the initial saturation-vs-