Estimation of the Effective Permeability of Heterogeneous Porous Media by Using Percolation Concepts

In this paper we present new methods to estimate the effective permeability (k_eff) of heterogeneous porous media with a wide distribution of permeabilities and various underlying structures, using percolation concepts. We first set a threshold permeability (k_th) on the permeability density function (pdf) and use standard algorithms from percolation theory to check whether the high permeable grid blocks (i.e. those with permeability higher than k_th) with occupied fraction of “p” first forms a cluster connecting two opposite sides of the system in the direction of the flow (high permeability flow pathway). Then we estimate the effective permeability of the heterogeneous porous media in different ways: a power law (k_eff=k_th p^m), a weighted power average (k_eff=[p.k_th^m+(1-p).k_g^m ]^(1/m) with k_g the geometric average of the permeability distribution) and a characteristic shape factor multiplied by the permeability threshold value. We found that the characteristic parameters (i.e. the exponent “m”) can be inferred either from the statistics and properties of percolation sub-networks at the threshold point (i.e. high and low permeable regions corresponding to those permeabilities above and below the threshold permeability value) or by comparing the system properties with an uncorrelated random field having the same permeability distribution. These physically based approaches do not need fitting to the experimental data of effective permeability measurements to estimate the model parameter (i.e. exponent m) as is usually necessary in empirical methods. We examine the order of accuracy of these methods on different layers of 10th SPE model and found very good estimates as compared to the values determined from the commercial flow simulators

Multiscale Computations for Flow and Transport in Heterogeneous Media

Many problems of fundamental and practical importance have multiple scale solutions. The direct numerical solution of multiple scale problems is difficult to obtain even with modern supercomputers. The major difficulty of direct solutions is due to disparity of scales. From an engineering perspective, it is often sufficient to predict macroscopic properties of the multiple-scale systems, such as the effective conductivity, elastic moduli, permeability, and eddy diffusivity. Therefore, it is desirable to develop a method that captures the small scale effect on the large scales, but does not require resolving all the small scale features. The purpose of this lecture note is to review some recent advances in developing multiscale finite element (finite volume) methods for flow and transport in strongly heterogeneous porous media. Extra effort is made in developing a multiscale computational method that can be potentially used for practical multiscale for problems with a large range of nonseparable scales. Some recent theoretical and computational developments in designing global upscaling methods will be reviewed. The lectures can be roughly divided into four parts. In part 1, we review some homogenization theory for elliptic and hyperbolic equations. This homogenization theory provides a guideline for designing effective multiscale methods. In part 2, we review some recent developments of multiscale finite element (finite volume) methods. We also discuss the issue of upscaling one-phase, two-phase flows through heterogeneous porous media and the use of limited global information in multiscale finite element (volume) methods. In part 4, we will consider multiscale simulations of two-phase flow immiscible flows using a flow-based adaptive coordinate, and introduce a theoretical framework which enables us to perform global upscaling for heterogeneous media with long range connectivity