Reservoir Engineering - General - The Use of the Method of Characteristics in Determining Boundar...

For many years the problem of increasing ultimate recovery of oil from a reservoir has been a subject of interest to the oil industry. At present, a standard secondary recovery technique is to flood the reservoir with water once natural forces of the system are depleted. The permeability characteristics of the rock, together with interfacial forces which exist between the two fluids and the porous material. allow only part of the oil in the rock formation to be recovered by the water-flooding process. Laboratory experiments have shown that when detergents are added to the flooding water, the flow properties of oil in a porous rock arc improved. This improvement results in greater oil recovery. Even though higher recoveries are obtained using detergents, the extra cost associated with the method may make it economically unsatisfactory. To establish economic feasibility requires experimental laboratory measurement of fluid and rock properties together with a calculational procedure to predict field results. Pilot tests in the field may also be required. This paper presents a mathematical procedure for predicting detergent flood behavior. A number of reasonable simplifying assumptions are made to obtain a tractable system of equations. The resulting method may be regarded as a generalization of Buckley-Lev-crett' theory to solve simultaneously for water saturation and detergent concentration profiles. Consequently this type of calculation can be used in much the same manner as Buckley-Leverett theory is used to evaluate a water flood. The mathematical method employed is the so-called method of characteristics for solving hyperbolic partial differential equations. This application of the method is instructive because several other problems in reservoir analysis may be analyzed in a similar manner. FLUID FLOW EQUATIONS The flow behavior of oil and water without detergent through porous reservoir rock is usually expressed by Darcy's law and a continuity equation for each phase. If gravitational forces are neglected, vrg = kr (s)µr v p rg.........(1) vrg = kr (s)µr v p rg.........(1) ? pa va = F ?/?t (sd) where v is the volumetric flux, k represents permeability, S is saturation, and P the pressure. Fluid viscosity is denoted by 11, p is the density, 4 is porosity, and I the time. The subscripts w and o refer to water and oil. The capillary forces in the system lead to a difference in pressure between phases. This may be expressed by the relation: Pm -P0 =Pc(S).........(5) The capillary pressure? P,, is determined experimentally as a function of saturation. DETERGENT BALANCE EQUATION Addition of detergent to the water phase affects the flow properties of the system. Capillary pressure and relative permeability become functions of detergent concentration, c, as well as fluid saturation. That is: P--P (S,c) and k,,,.- - -- k,,,. At low oil saturation values, the presence of detergent reduces capillary pressure and reduces the relative permeability ratio. Detergents are surface-active and adsorb on rock surfaces. The quantity adsorbed, a, will depend on the concentration c of detergent solution and may also depend on the saturation; thus, a = a(S,c). In this analysis we have assumed that the detergents used are water soluble but not oil soluble, and their movement is therefore restricted to the water phase. The dispersion of detergent in water is neglected by comparison with transfer by flow within the system. Detergent movement is determined by the water flux and adsorption: ? cv m = -F ?/?t (cs + a) We will assume that the quantity of detergent in solution is too small to significantly change water density or viscosity. When the permeability and capillary pres-