Technical Notes - A Note on Transient Two-Phase Flow Calculations

Since the appearance of the paper, "Solution of the Equations of Un-steady State Two-Phase Flow in Oil Reservoirs," by W. J. West, W. W. Garvin, and J. W. Sheldon,' a two-fold investigation of this subject has been carried out. One objective of the investigation has been to deter-mine the feasibility of solving such problems 7on a medium-size com-puter such as the Datatron*, and the other objective has been to in-vestigate the application of such cal-culations to experimental and theo-retical petroleum reservoir research. In the first Datatron calculations, the fluid and rock properties published by West, et al, were used, together with the published equations describ-ing the system. Details of the formu-lation not given in the original paper are discussed in the Appendix to this note. Reference should be made to the subject paper for the complete equations and defintion of symbols. An unexpected result of this in-vestigation was the discovery that the linear solution published by West was in error. Thus, in addition to describing the Datatron solutions and to discussing certain numerical diffi-culties which will be encountered if one uses the published method of solution, the purpose of this note is to indicate the nature of this error. LINEAR FLOW Since the linear case requires a minimum amount of scaling, a fixed-decimal point Datatron program was written for the one-dimensional flow problem and an attempt was made to duplicate the solution described by West. In the case described, fluid was produced at a constant rate, Q, until such time as well pressure reached 0.04. Production was then continued at constant pressure. From the constants and curves given by West it was determined that the ini-tial constant production rate could be approximated by Q = 0.007. An ini-tial dimensionless time step ?t = 0.434 X 10 - "as used, and each suc-cessive time step was doubled until a value of = 0.444 was reached. This constant interval was then used for the remainder of the solution. In subsequent solutions, several varia-tions in the time schedule were em-ployed, including smaller time steps and slower rates of increase in the time steps. In all cases, almost identical results were obtained regardless of the time schedule employed. However, as described below, it was noted that the time schedule had some influence on the rate of convergence of the solutions. As a check on the accuracy of the solution, the cumulative production at each time step was calculated using the two methods described in the Appendix. Satisfactory agreement was observed with the differences in these two values of the order of two parts in 50,000. It should be noted that the mass balance check as described is of questionable value, particularly with regard to the well pressure and saturation. This is especially true in the radial solution where pressure and saturation values near the wellbore would make only a negligible contribution to the numer-ical integration. It is believed, how-ever, that such a comparison is ot value in determining the over-all accuracy of a solution. In comparing the Datatron solu-tion with that published by West it was discovered that in the later stages of depletion, the pressures near the well declined more rapidly in our solution than in the West solution, and that the limiting well pressure of 0.04 was reached at an earlier time than that originally reported. It thus became evident that it would be impossible to duplicate the production schedule described by West and a constant rate of production was maintained until the well pressure was equal to 0.0. A representative comparison of the results published by West with those obtained in this investigation is shown in Fig. 1, which is a plot of GOR as a function of cumulative recovery. These two curves should be in agreement until a cumulative recovery is reached which corresponds to a well pressure of 0.04 — for the Datatron solution, a recovery of approximately 5.6 per cent. Actually, a major disagreement is evident. Subsequent correspondence