Reservoir Engineering–General - Flow in Heterogeneous Porous Media

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. 1. 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. 2. 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. 3. The permeability determined from a pressure build-up curve lor a heterogeneous reservoir gives a reasonable value for the effective permeability of the drainage area. 4. 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. 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. Extensive research, both theoretical and experimental, has provided reasonable descriptions of the physical processes that are associated with primary and secondaryrecovery mechanisms. Similar progress in the fie1d of numerical analysis and the evolution of truly high-speed digital computers have made it feasible to calculate the behavior of these mechanisms under almost any environmental conditions. It appears to be logical, then, to conclude that the mechanistic performance of any reservoir can be predicted with some confidence. Unfortunate1y, a necessary condition for the validity of this conclusion is that the formation be adequately described in a macroscopic sense. This condition cannot be satisfied ill general; therefore, a frustrating situation currently exists — it is possible to compute the effect of reservoir heterogeneity on behavior, but it is not possible to specify the heterogeneity itself. In many cases, the predicted behavior of a project is so completely dominated by irregularities in the formation's physical properties that the gratuitous assumptior, of a particular form for the variation can reduce the solution of the problem to a mere tautological exercise. While all porous media are microscopically heterogeneous, only macroscopic fluctuations need be considered since the fundamental concepts of reservoir mechanics are based on macroscopic quantities. For the simplest case in which only one physical property varies from point to point, the degree of heterogeneity that is apparent depends on three distinct factors, i.e., the nature of the property's variation (frequency distribution), its disposition (spatial distribution) and the inherent stability of the mechanism being studied. Consequently, any attempt to devise an unambiguous scale of heterogeneity must take these factors into account. Since this investigation is intended to be exploratory, it will be restricted to single-phase flow because of the simplification implied; both steady-state and unsteady-state behavior will be considered. Furthermore, to keep the number of parameters