A New Approach to Forecasting Miscible WAG Performance at the Field Scale

1998 ◽  
Vol 1 (03) ◽  
pp. 192-200 ◽  
Author(s):  
R.M. Giordano ◽  
R.S. Redman ◽  
F. Bratvedt

Abstract Full-field EOR performance predictions are generally obtained from scale-up tools, since three-dimensional finite-difference simulations would be too CPU intensive. Existing scale-up techniques require the user to define pattern elements and then to derive performance curves to apply to each injector-producer pair in the elements. Accurate assignment of these elements is difficult because the actual shape and size of the swept volumes are sensitive to reservoir faulting, well rate changes, and regional flux. In reality, the actual sweep region is not an input parameter, but should be determined by the regional pressure field which changes as well rates vary and new wells are drilled. Thus, a major source of error in using existing scale-up tools is trying to define representative pattern elements. In the current paper, we describe a scale-up technique in which the user does not have to define pattern elements or injector-producer pairs. In the new technique, the pressure field is computed at each time step and then a front-tracking algorithm propagates water and miscible injectant throughout the reservoir. By using an analogy between oil mobilization and adsorption/desorption of tracers, the miscible-gas process is modeled. The parameters for the model are obtained by fine-scale, two-dimensional, compositional, finite-difference simulations in a vertical cross-section. In the new approach, the injected solvent is divided into an effective and an ineffective portion. This approach reduces a three-dimensional problem to a two-dimensional, areal one in which the declining displacement efficiency of the solvent, which is caused by vertical effects, is captured by decreasing the injected concentration of effective solvent with time. In this paper, we show how the new scale-up tool has been used to model the miscible WAG process in the Eastern Peripheral Wedge Zone of the Prudhoe Bay field. We show a comparison between field response and model predictions. Introduction Good reservoir management requires the prediction of reliable oil and gas rates. In general, the degree of difficulty in making these predictions depends on the displacement process. For example, good predictions of primary depletion or gravity drainage by gas-cap expansion can usually be obtained by coarsely gridded finite-difference simulations. However, processes where injected water or gas must be tracked from injector to producer typically require finely-gridded simulations. Accurate prediction of oil and gas rates frequently require finely gridded simulations which contain (1) rock-measured relative permeabilities and (2) a reservoir description that accurately predicts high-permeability zones (thieves) and low-permeability barriers (e.g., shale location, size, and continuity). At Prudhoe Bay, modeling of miscible gas processes generally requires vertical grid blocks of the order of one foot to match field-measured saturation profiles. At the present time, three-dimensional, compositional modeling of gas displacement processes that satisfy these two requirements require at least a week of CPU time on IBM-590 workstation for a single pattern. Thus, it is not currently practical to use finely gridded finite-difference simulators to model large sections of a field. Traditionally, three approaches have been used to address this problem -- pseudo relative permeabilities, tank models, and streamtubes. Pseudo relative permeabilities are generally successful only when the saturation history experienced in the coarse-grid simulation will always be similar to the fine-grid simulation. Tank models can be difficult to apply when the original pattern changes by infill drilling or well conversions, and streamtube models have had difficulty when the initial conditions are not homogeneous along each streamline. To address the above problems, a new approach was created that can reproduce the response and timing characteristics of the produced components, but also has the ability to propagate and track injected fronts. In addition, the model does not require user-supplied injector-producer allocation factors. We explain, below, our new front-tracking technique and how this new scale-up tool has been used to model the miscible WAG process in the Eastern Peripheral Wedge Zone of Prudhoe Bay. P. 329

2017 ◽  
Vol 7 (17) ◽  
pp. 3707-3714 ◽  
Author(s):  
Carlos A. Carrero ◽  
Samuel P. Burt ◽  
Fangying Huang ◽  
Juan M. Venegas ◽  
Alyssa M. Love ◽  
...  

Dispersing two-dimensional VOx species on β-SiC offers a new approach to scale up propane ODH.


Author(s):  
V. Vlasenko ◽  
A. Shiryaeva

New quasi-two-dimensional (2.5D) approach to description of three-dimensional (3D) flows in ducts is proposed. It generalizes quasi-one-dimensional (quasi-1D, 1.5D) theories. Calculations are performed in the (x; y) plane, but variable width of duct in the z direction is taken into account. Derivation of 2.5D approximation equations is given. Tests for verification of 2.5D calculations are proposed. Parametrical 2.5D calculations of flow with hydrogen combustion in an elliptical combustor of a high-speed aircraft, investigated within HEXAFLY-INT international project, are described. Optimal scheme of fuel injection is found and explained. For one regime, 2.5D and 3D calculations are compared. The new approach is recommended for use during preliminary design of combustion chambers.


1951 ◽  
Vol 3 (3) ◽  
pp. 193-210 ◽  
Author(s):  
M.J. Lighthill

SummaryThe general technique for rendering approximate solutions to physical problems uniformly valid is here applied to the simplest form of the problem of correcting the theory of thin wings near a rounded leading edge. The flow investigated is two-dimensional, irrotational and incompressible, and therefore the results do not materially add to our already extensive knowledge of this subject, but the method, which is here satisfactorily checked against this knowledge, shows promise of extension to three-dimensional, and compressible, flow problems.The conclusion, in the problem studied here, is that the velocity field obtained by a straightforward expansion in powers of the disturbances, up to and including either the first or the second power, with coefficients functions of co-ordinates such that the leading edge is at the origin and the aerofoil chord is one of the axes, may be rendered a valid first approximation near the leading edge, as well as a valid first or second approximation away from it, if the whole field is shifted downstream parallel to the chord for a distance of half the leading edge radius of curvature ρL. It follows that the fluid speed on the aerofoil surface, as given on such a straightforward second approximation as a function of distance x along the chord, similarly is rendered uniformly valid (see equation (52)) if the part singular like x-1 is subtracted and the remainder is multiplied by .


1973 ◽  
Vol 13 (03) ◽  
pp. 175-185 ◽  
Author(s):  
Hugh H. Jacks ◽  
Owen J.E. Smith ◽  
C.C. Mattax

Abstract Dynamic pseudo-relative permeabilities derived from cross-section models can be used to simulate three-dimensional flow accurately in a two-dimensional areal model of a reservoir Techniques are presented for deriving and using dynamic pseudos that are applicable over a wide range of rates and initial fluid saturations. Their validity is demonstrated by showing calculated results from comparable runs in a vertical cross-section model and in a one-dimensional areal model using the dynamic pseudo-relative permeabilities and vertical equilibrium (VE) pseudo-capillary pressures. Further substantiation is provided by showing the close agreement in calculated performance for a three-dimensional model and corresponding two-dimensional areal model representing a typical pattern on the flanks of a large reservoir. The areal pattern on the flanks of a large reservoir. The areal model gave comparable accuracy with a substantial savings in computing and manpower costs. Introduction Meaningful studies can be made for almost all reservoirs now that relatively efficient three-dimensional reservoir simulators are available. In many instances, however, less expensive two-dimensional areal (x-y) models can be used to solve the engineering problem adequately, provided the nonuniform distribution and flow of fluids in the implied third, or vertical, dimension of the areal model is properly described. This is accomplished through the use of special saturation-dependent functions that have been labeled pseudo-relative permeability (k ) and pseudo-capillary pressure permeability (k ) and pseudo-capillary pressure (P ) or, for simplicity "pseudo functions", to distinguish them from the conventional laboratory measured values that are used in their derivation. Two types of reservoir models have been suggested in the past to derive pseudo functions: the vertical equilibrium (VE) model of Coats et al., which is based on gravity-capillary equilibrium in the vertical direction; and the stratified model of Hearn, which assumes that viscous forces dominate vertical fluid distribution. Neither of these models accounts for the effects of large changes in flow rate that take place as a field is developed, approaches and place as a field is developed, approaches and maintains its peak rate, and then falls into decline. This paper presents an alternative method for developing pseudo functions that are applicable over a wide range of flow rates and over the complete range of initial fluid saturations. The functions may be both space and time dependent and, again for clarity and convenience in nomenclature, we have labeled them "dynamic pseudo functions". DESCRIPTION OF PSEUDO-RELATIVE PERMEABILITY FUNCTIONS PERMEABILITY FUNCTIONS Methods for developing pseudo functions have been presented in the literature. The distinction between our method and those used by others lies in the technique for deriving the vertical saturation distribution upon which the pseudo-relative permeabilities are based. In our approach, the permeabilities are based. In our approach, the vertical saturation distribution is developed through detailed simulation of the fluid displacement in a vertical cross-section (x-z) model of the reservoir. The simulation is run under conditions that are representative of those to be expected during the period to be covered in the areal model simulations. period to be covered in the areal model simulations. Results of the cross-section simulation are then processed to give depth-averaged fluid saturations processed to give depth-averaged fluid saturations (S ) and "dynamic" pseudo-relative permeability values (k ) for each column of blocks in the cross-section model at each output time. The above approach can result in a different set of dynamic pseudo functions for each column of blocks due to differences in initial saturation, rate of displacement, reservoir stratification, and location. However, differences between columns are frequently minor or they can be accounted for by correlation of the data. In this and several other reservoir studies, it was possible to reduce the complexity of the pseudo function sets through correlations with initial fluid saturations and fluid velocities. SPEJ P. 175


Author(s):  
Feng Liang ◽  
Quanyong Xu ◽  
Xudong Lan ◽  
Ming Zhou

The thermohydrodynamic analysis of oil film bearing is essential for high speed oil film bearing. The temperature field is coupled with the pressure field. The numerical model can be built or chosen according to the complexity of the objects and requirement of the accuracy. In this paper, four pressure-thermal (P-T) models are proposed, which are zero-dimensional temperature field coupled with Reynolds equation (0D P-T model), two-dimensional temperature field coupled with Reynolds equation (2D P-T model), two-dimensional temperature with third dimensional correction coupled with Dawson equation (2sD P-T model), three-dimensional temperature field coupled with Dawson equation (3D P-T model). The non-Newtonian and temperature-viscosity effects of the lubrication oil are considered in all the four models. Two types of cylindrical journal bearing, the bearing with/without axial grooves, are applied for the simulation. All the simulated cases are compared with the solutions of the CFX. The results show that the 0D P-T model fails to predict the behavior of high speed bearing; The 2D and 2sD P-T model have an acceptable accuracy to predict the performance of the bearing without grooves, but are not able to simulate the P-T field of the bearing with grooves because of the under-developed thermal boundary layer; The 2sD P-T model shows a great improvement when calculating the pressure field compared with the 2D P-T model; the 3D P-T model coincides well with the CFX at any condition. The comparison of these four models provides a reference to help designer choose a proper numerical model for a certain project.


1968 ◽  
Vol 8 (02) ◽  
pp. 185-194 ◽  
Author(s):  
J.E. Briggs ◽  
T.N. Dixon

Abstract A study was made of numerical techniques for solving the large sets of simultaneous equations that arise in the mathematical modeling of oil reservoir behavior. It was found that noniterative techniques, such as the Alternating Direction Implicit (ADI) method, as well as some other finite difference approximations, produce oscillatory or unsmooth results for large time steps. Estimates of time step sizes sufficient to avoid such behavior are given. A comparison was made of the Point Successive Over-Relaxation (PSOR), Two-Line Cyclic Chebyshev Semi-Iterative SOR (2LCC), and iterative ADI methods, with respect to speed of solution of a test problem. It was found that, when applicable, iterative ADI is fastest for problems involving many points, while 2LCC is preferable for smaller problems. Introduction With the advent of high speed, large memory, digital computers, there has been an increasing emphasis on the development of improved methods for simulating and predicting reservoir performance. Two-dimensional, three-phase reservoir models with various combinations of PVT effects, as well as gravity and capillary forces, are common throughout the industry. Such models are also available through consulting firms, to anyone desiring to use them. Three-dimensional models will probably be practical in only a few years. We conducted a study of some of the numerical methods used for solving the large sets of simultaneous equations that arise in such models. A typical set of equations for a reservoir model is shown below: ..............(1) .............(2) .............(3) ..............(4) where a = 5.615 cu ft/bbl. In addition to Eqs. 1 through 4, one also would have to specify the conditions at the boundaries of the reservoir or aquifer being studied. Equations such as these are normally approximated by finite difference techniques and solved numerically because of their complexity. In deciding how to solve such equations, a number of decisions must be made. It is not our intention to cover all facets of the problem, but rather to concentrate on one of the important aspects, such as solving Eq. 1. SPEJ P. 185ˆ


1967 ◽  
Vol 7 (04) ◽  
pp. 377-388 ◽  
Author(s):  
K.H. Coats ◽  
R.L. Nielsen ◽  
Mary H. Terhune ◽  
A.G. Weber

COATS, K.H., THE U. OF TEXAS, AUSTIN, TEX. NIELSEN, R.L., ESSO PRODUCTION RESEARCH CO., HOUSTON, TEX. MEMBERS AIME TERHUNE, MARY H., AMERICAN AIRLINES, TULSA, OKLA., WEBER, A.G., ESSO PRODUCTION RESEARCH CO., HOUSTON, TEX. MEMBER AIME Abstract Two computer-oriented techniques for simulating the three-dimensional flow behavior of two fluid phases in petroleum reservoirs were developed. Under the first technique the flow equations are solved to model three-dimensional flow in a reservoir. The second technique was developed for modeling flow in three-dimensional media that have sufficiently high permeability in the vertical direction so that vertical flow is not seriously restricted. Since this latter technique is a modified two-dimensional areal analysis, suitably structured three-dimensional reservoirs can be simulated at considerably lower computational expenses than is required using the three-dimensional analysis. A quantitative criterion is provided for determining when vertical communication is good enough to permit use of the modified two-dimensional areal analysis. The equations solved by both techniques treat both fluids as compressible, and, for gas-oil applications, provide for the evolution of dissolved gas. Accounted for are the effects of relative permeability, capillary pressure and gravity in addition to reservoir geometry and rock heterogeneity. Calculations are compared with laboratory waterflood data to indicate the validity of the analyses. Other results were calculated with both techniques which show the equivalence of the two solutions for reservoirs satisfying the vertical communication criterion. Introduction Obtaining the maximum profits from oil and gas reservoirs during all stages of depletion is the fundamental charge to the reservoir engineering profession. In recent years much quantitative assistance in evaluating field development programs has been goaded by computerized techniques for predicting reservoir flow behavior. Because of the spatially distributed and dynamic nature of producing operations, automatic optimization procedures, such as those now in use for planning refining operations, are not now available for planning reservoir development. However, present mathematical simulation techniques do furnish powerful means for making case studies to help in planning primary recovery operations and in selecting and timing supplemental recovery operations. A number of methods have been reported which simulate the flow of one, two or three fluid phases within porous media of one or two effective spatial dimensions. However, applying computer analyses to actual reservoirs have been limited mostly to two-dimensional areal or cross-sectional flow studies for two immiscible reservoir fluids. To obtain a three-dimensional picture of reservoir performance using such two-dimensional techniques, it has been necessary to interpret the calculations by combining somehow the results from essentially independent areal and cross-sectional studies. To the author's knowledge, the only other three-dimensional computational procedure, in addition to those presented here, was developed by Peaceman and Rachford to simulate the behavior of a laboratory waterflood. Two computational techniques which may be used to simulate three-dimensional flow of two fluid phases are described in this paper. The first method, called the "three-dimensional analysis", employs a fully three-dimensional mathematical model that treats simultaneously both the areal and cross-sectional aspects of reservoir flow. SPEJ P. 377ˆ


Author(s):  
Mirka Deza ◽  
Francine Battaglia ◽  
Theodore J. Heindel

Fluidized beds are used to gasify materials such as coal or biomass in the production of producer gas. Modeling these reactors using computational fluid dynamics is advantageous when performing parametric studies for design and scale-up. While two-dimensional simulations are easier to perform than three-dimensional simulations, they may not capture the proper physics. This paper compares two- and three-dimensional simulations with experiments for a reactor geometry with side port air injection. The side port is located within the bed region so that the injected air can help promote mixing. Of interest in this study is validating the hydrodynamics of fluidizing biomass. Two operating conditions of the fluidized bed are studied for superficial gas velocities of 1.5Umf and 3.0Umf, where Umf is the minimum fluidization velocity. The material used to represent biomass is ground walnut shell because it tends to fluidize uniformly and falls within the Geldart type B classification. The simulations are compared to experimental data of time-averaged local gas holdup values using X-ray computed tomography. Results indicate that for the conditions of this study, two-dimensional simulations overpredict the gas holdup trends when compared to the experiments. However, the three-dimensional simulations compare exceptionally well with the experiments, thus predicting the fluidization hydrodynamics, irrespective of flowrate or complexity due to the side air port. Furthermore, the study demonstrates the importance of using a three-dimensional model for bubbling fluidized beds with complex physics.


Sign in / Sign up

Export Citation Format

Share Document