Miscellaneous - Relaxation Methods Applied to Oilfield Research

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.' THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy's Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. 'Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If' .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A