Techniques for studying the performance characteristics of heterogeneous reservoirs have been developed. The effect of permeability variation on both the steady-state and the transient flow of a single fluid has been investigated. Limited comparisons with field and laboratory data have been made. The physical models studied consist of random three-dimensional arrays of homogeneous porous blocks. The permeabilities of the individual elements are assigned according to a specific distribution function; uniform anisotropy is introduced by varying the relative dimensions of the blocks. A particular model is perturbed simply by re-arranging its elements at random. The behavior of a physical model is determined by digitally solving its numerical analogue. Based on computational experiments, subject to the restrictions implied by the assumptions that were made, the following general conclusions have been drawn.The most probable behavior of a heterogeneous system approaches that of a homogeneous system with a permeability equal to the geometric mean of the individual permeabilities.The effects of flow geometry and anisotropy on the most probable value of the effective permeability of a heterogeneous medium are finite but not significant.The permeability determined from a pressure build-up curve for a heterogeneous reservoir gives a reasonable value for the effective permeability of the drainage area.A qualitative measure of the degree of heterogeneity and its spatial configuration are obtained from a comparative study of core analysis and pressure build-up data. It has been indicated that these conclusions are predicated on the assumption that the core analysis and the pressure build-up data represent true reservoir characteristics. Common sources of error in these data have been discussed. Introduction One of the most significant problems of reservoir engineering is that of determining the nature and the disposition of the heterogeneities that inevitably occur in petroliferous formations. It is with this problem that this paper is concerned.

No abstract

A numerical {formulation o~high older accuracy, based on variational metipods, is proposed /or the so[ution 0/ multidimensional di//usion-convectionty/re equations. Accurate 10,

During the past decade, efforts in reservoir modeling have focused on the three areas of capability efficiency, and reliability. Capability means the ability to handle larger and more complex problems where complexity includes physical problems where complexity includes physical phenomena, such as gas percolation and variable phenomena, such as gas percolation and variable PVT properties, and severe heterogeneity due to PVT properties, and severe heterogeneity due to property variation or geometry, or both. Efficiency property variation or geometry, or both. Efficiency is increased by improving model formulations anti solution techniques to increase tolerable time-step size and reduce computer time per time step. Reliability refers to ease of use and minimum burden in selecting or experimenting with time-step size, solution technique options, iteration parameters, and closure tolerances. parameters, and closure tolerances. The single facet of a reservoir simulator that has the greatest combined influence in all three categories is the technique used to solve the large systems of equations arising from the numerical approximation of the nonlinear fluid flow equations. Available techniques include both direct solution and iterative methods such as ADIP, SOR, and SIP. Iterative methods are currently used almost to the exclusion of direct solution because of the significantly higher computer storage and time requirements of the latter. This paper describes some new ordering schemes for Gaussian elimination that reduce computing time and storage requirements by factors as large as 6 and 3, respectively, relative to more standard orderings. Computational work estimates are given for these methods, for the standard Gaussian ordering, and for several iterative methods. These work estimates are checked by comparisons of actual run times using different solution techniques. Numerical examples are given to illustrate the increased efficiency and reliability that can be achieved in many cases through use of the new direct solution methods. Introduction It is well known that the way we number or order the unknowns of a sparse system of linear algebraic equations can drastically affect the amount of computation and storage for a direct solution. However, until recently the best ordering scheme that appeared in the literature numbered the points of a three-dimensional grid first along the shortest direction - i.e., the dimension with the fewest number of grid points - then in the next shortest direction, and finally in the longest direction. This ordering, which we shall call the standard ordering for Gaussian elimination is still widely used even though it is substantially slower than many other orderings. Ogbuobiri et al. present a survey of the literature related to ordering schemes that exploit matrix sparsity. These schemes are grouped into the two classes of matrix-banding schemes and optimal or pseudo-optimal schemes. The latter schemes pseudo-optimal schemes. The latter schemes purport to yield generally greater efficiency. purport to yield generally greater efficiency. In a recent paper, Georges has shown that for five-point difference approximations on square n x n two-dimensional grids, the total work for certain orderings of the grid points is less than C1n3 and the storage is less than C2n2 log n, compared with n4 and n3, respectively, for the standard ordering. Moreover, George has shown that no ordering scheme can require less work than the order of n3. For the special case of n - 21 he shows that work W is less than 10n3 and the storage S is less than 8ln2 for symmetric matrices. For nonsymmetric matrices these results become W less than 20n3 and S less than 16ln2, respectively. In this paper we describe some specific orderings in the matrix-banding class. Analyses of work and storage requirements are given for these orderings as applied to the diffusivity-type pressure equation that arises in reservoir simulation problems. These work and storage requirements are compared with those of the standard Gaussian ordering and of some iterative methods. These comparisons are performed for problems ranging from simple performed for problems ranging from simple homogeneous squares to practical reservoir problems of typical heterogeneity and irregular problems of typical heterogeneity and irregular geometry. The work requirements of the orderings presented here are also compared experimentally with those of one of the leading pseudo-optimal schemes. SPEJ P. 295