Three-Phase Fluid Flow Including Gravitational, Viscous and Capillary Forces

1969 ◽  
Vol 9 (02) ◽  
pp. 255-269 ◽  
Author(s):  
M. Sheffield

Abstract This paper presents a technique for predicting the flow, of oil, gas and water through a petroleum reservoir. Gravitational, viscous and capillary forces are considered, and all fluids are considered to be slightly compressible. Some theoretical work concerning the fluid flow in one-, two- and three-space dimensions is given along with example performance predictions in one- and two-space dimensions. predictions in one- and two-space dimensions Introduction Since the introduction of high speed computing equipment one of the goals of reservoir engineering research has been to develop more accurate methods of describing fluid movement through underground reservoirs. Various mathematical methods have been developed or used by reservoir engineers to predict reservoir performance. The work reported in this paper extends previously published work on three-phase fluid flow (1) by including a rigorous treatment of capillary forces and (2) by showing certain theoretical mathematical results proving that these equations can be approximated by certain numerical techniques and that a unique solution exists. Discussion The method of predicting three-phase compressible fluid flow in a reservoir can be summarized briefly by the following steps. 1.The reservoir, or a section of a reservoir, is characterized by a series of mesh points with varying rock and fluid properties simulated at each mesh point. point. 2.Three partial differential equations are written to describe the movement at any point in the reservoir of each of the three compressible fluids. All forces influencing movement are considered in the equations. 3. At each mesh point, the partial differential equations are replaced by a system of analogous difference equations. 4. A numerical technique is used to solve the resulting system of difference equations. Capillary forces have been included in two-phase flow calculations. The literature, however, does not contain examples of prediction techniques for three-phase flow that include capillary forces. Where capillary Races are considered, each of the three partial differential equations previously discussed has a different dependent variable, namely pressure in one of the three fluid phases. Therefore, pressure in one of the three fluid phases. Therefore, three difference equations must be solved at each point in the reservoir. Where large systems of point in the reservoir. Where large systems of equations must be solved simultaneously, an engineer might question whether a unique solution to this system of equations actually exists and, if so, what numerical techniques may be used to obtain a good approximation to the solution. It is shown in the Appendix that a unique solution to the three-phase flow problem, as formulated, always exists. It is also shown that several methods may be used to obtain a good approximation to the solution. The partial differential equations and difference equations partial differential equations and difference equations used are shown in the Appendix. Matrix notation has been used in developing the mathematical results. Two sample problems were solved on a CDC 1604. They illustrate the type of problems that can be solved using a three-phase prediction technique. SAMPLE ONE-DIMENSIONAL RESERVOIR PERFORMANCE PREDICTION A hypothetical reservoir was studied to provide an example of a one-dimensional problem that can be solved. Of the several techniques available, the direct method A solution as shown in the Appendix was used. The reservoir section studied was a truncated, wedge-shaped section, 2,400 ft long, with a 6' dip. (A schematic is shown in Fig. 1.) This section was represented by 49 mesh points, uniformly spaced at 50-ft intervals. The upper end of the wedge was 2 ft wide, and the lower end was 6 ft wide. SPEJ P. 255

Author(s):  
Renfrey B. Potts

AbstractOrdinary difference equations (OΔE's), mostly of order two and three, are derived for the trigonometric, Jacobian elliptic, and hyperbolic functions. The results are used to derive partial difference equations (PΔE's) for simple solutions of the wave equation and three nonlinear evolutionary partial differential equations.


This chapter describes the pdepe command, which is used to solve spatially one-dimensional partial differential equations (PDEs). It begins with a description of the standard forms of PDEs and its initial and boundary conditions that the pdepe solver uses. It is shown how various PDEs and boundary conditions can be represented in standard forms. Applications to the mechanics are presented in the final part of the chapter. They illustrate how to solve: heat transfer PDE with temperature dependent material properties, startup velocities of the fluid flow in a pipe, Burger's PDE, and coupled FitzHugh-Nagumo PDE.


1970 ◽  
Vol 10 (02) ◽  
pp. 192-202 ◽  
Author(s):  
R.B. Lantz

Abstract In the past miscible displacement calculations have been approximated with two-phase reservoir simulators. Such calculations have neglected diffusional mixing between miscible components. In fact, no analog bas been proposed for rigorously treating miscible simulations with two-phase programs. This paper describes requirements that programs. This paper describes requirements that permit such a rigorous simulation. permit such a rigorous simulation. The sets of partial differential equations describing each of the displacement processes are shown to be exactly analogous if relative permeability and capillary pressure functions are permeability and capillary pressure functions are adjusted in a special manner. Application of the "miscible" analogy in a two-phase simulator, however, has several limitations, the most severe of which is the truncation error (numerical diffusion) typical of an immiscible formulation. Since this error is time-step and/or block-size dependent, numerical smearing can, in principle, be made as small as necessary. But this feature limits the practical applicability of the "miscible" analogy practical applicability of the "miscible" analogy to cases with rather large physical diffusion. The range of applicability and other limitations are outlined in the paper. Also, illustrative sample calculations are presented for linear, radial and layer-cake systems. Component densities and viscosities are varied in the linear model. Introduction In recent years, use of two- and three-phase reservoir simulators to calculate the performance of immiscible fluid displacement has become widespread. Reservoir simulators capable of calculating miscible displacement problems, however, have been limited to special use programs. The primary reason for this limitation has been the significant truncation error (numerical diffusion) typical of ordinary finite difference approximations to the miscible equations. The method of characteristics provided a means of making miscible displacement calculations without significant truncation error. A recently proposed second calculation technique, based on variational methods, also reduces numerical diffusion. Both of these calculational techniques can be used for immiscible calculations. Still, general miscible displacement applications such as gas cycling, enriched-gas injection, or tracer injection have historically required use of immiscible reservoir simulators for performance predictions. Larson et al. have reported an example of such use of a two-phase computer program. Displacement involving two components flowing within a single phase would appear to be analogous to a two-phase displacement. Yet, past miscible calculations using immiscible simulators made the two-phase saturation profile as near piston-like as possible and neglected component mixing due to possible and neglected component mixing due to diffusional processes. The capillary pressure function was chosen to minimize capillary flow. Also, in these miscible approximations, no provision had been made for viscosity variations provision had been made for viscosity variations with component concentration. Though mixing due to diffusional processes had been neglected, countercurrent diffusion due to component concentration differences in a miscible process should be essentially analogous to countercurrent capillary flow due to saturation differences in a two-phase system. This paper describes a method by which two- and three-phase reservoir simulators can be made to calculate miscible displacement rigorously. The only requirement of the method is that relative permeability and capillary pressure be special permeability and capillary pressure be special functions of saturation. With these properly chosen functions, the set of partial differential equations describing immiscible displacement becomes completely analogous to the partial differential equations describing miscible displacement. SPEJ P. 192


2000 ◽  
Vol 43 (3) ◽  
pp. 485-510 ◽  
Author(s):  
Derek W. Holtby

AbstractThe purpose of this work is to establish a priori C2, α estimates for mesh function solutions of nonlinear positive difference equations in fully nonlinear form on a uniform mesh, where the fully nonlinear finite-difference operator ℱh is concave in the second-order variables. The estimate is an analogue of the corresponding estimate for solutions of concave fully nonlinear elliptic partial differential equations. We deal here with the special case that the operator does not depend explicitly upon the independent variables. We do this by discretizing the approach of Evans for fully nonlinear elliptic partial differential equations using the discrete linear theory of Kuo and Trudinger. The result in this special case forms the basis for a more general result in part II. We also derive the discrete interpolation inequalities needed to obtain estimates for the interior C2, α semi-norm in terms of the C0 norm.


Sign in / Sign up

Export Citation Format

Share Document