Large-scale properties for two-phase flow in random porous media

Abstract Hydrocarbon contaminants in the subsurface are important sources of pollution in aquifers. Hence, mathematical models of these flows have become key tools in environmental studies. In this paper we are interested in the flow in the saturated zone. Simulation of two-phase flow in large, complex heterogeneous domains often requires an unacceptably large number of computational grid blocks. Despite recent progress in computational methods and tools, we must call upon special techniques in order to use larger grid blocks while compensating for intracell variations in rock properties and fluid saturations. The use of pseudo-functions is one way of increasing grid dimensions to a more tractable level with a minimal loss of simulation representativeness. This change of scale problem has also been treated theoretically by different scaling-up techniques, such as large-scale averaging . This method calculates the transport equations and the effective properties at a given scale by an averaging process over the equations corresponding to a lower scale. This procedure leads to a closure problem which is very complex in the general case. Previously a first solution of the large-scale averaging problem was proposed in the quasi-static case corresponding to local capillary equilibrium . In the general case of a heterogeneous medium with a complex geometry, this boundary value problem (closure problem) can be solved by numerical methods. For this purpose, after having chosen a grid-block description of our system in accordance with the description used in reservoir simulators, we have implemented a three-dimensional (3-D) numerical resolution of the closure problem. The most important variations in the rock properties are associated with the permeability. We have therefore generated porous media with a random permeability distribution using different methods. Other multiphase properties are chosen to depend on the permeability. The properties of the closure problem in the case of randomly generated porous media are investigated. In particular, we give the conditions under which the general form of the local capillary pressure and relative permeability curves are recovered at the large scale. Particular properties related to each generation method are stated. The equivalent properties are calculated using averages over the results of many realizations of a given medium. The influence of the size of the averaging surface for a given correlation length as a function of the variance of the permeability is studied. We have therefore established some general rules for the calculation of the large-scale properties of random porous media.

[1]  John H. Cushman,et al.  Dynamics of fluids in hierarchical porous media , 1990 .

[2]  A. Mantoglou,et al.  The Turning Bands Method for simulation of random fields using line generation by a spectral method , 1982 .

[3]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[4]  Ann Muggeridge,et al.  GENERATION OF EFFECTIVE RELATIVE PERMEABILITIES FROM DETAILED SIMULATION OF FLOW IN HETEROGENEOUS POROUS MEDIA , 1991 .

[5]  Alain Bourgeat,et al.  Homogenized behavior of two-phase flows in naturally fractured reservoirs with uniform fractures distribution , 1984 .

[6]  Michel Quintard,et al.  Two-phase flow in heterogeneous porous media: The method of large-scale averaging , 1988 .

[7]  W. A. Bruce,et al.  Evaluation Of Capillary Character In Petroleum Reservoir Rock , 1949 .

[8]  G. Dagan.  Time‐dependent macrodispersion for solute transport in anisotropic heterogeneous aquifers , 1988 .

[9]  Peter K. Kitanidis,et al.  Determination of the effective hydraulic conductivity for heterogeneous porous media using a numerical spectral approach: 1. Method , 1992 .

[10]  Kenneth Stuart Sorbie,et al.  Permeability tensors for sedimentary structures , 1994 .

[11]  L. Dake.  Fundamentals of Reservoir Engineering , 1983 .

[12]  Gedeon Dagan,et al.  Models of groundwater flow in statistically homogeneous porous formations , 1979 .

[13]  A. K. Khurana,et al.  Pseudofunctions for Water Coning in a Three-Dimensional Reservoir Simulator , 1977 .

[14]  J. R. Kyte,et al.  New Pseudo Functions To Control Numerical Dispersion , 1975 .

[15]  Ahmadi Azita,et al.  Large-scale properties for flow through a stratified medium: Various approaches , 1993 .

[16]  G. Matheron.  Éléments pour une théorie des milieux poreux , 1967 .

[17]  W. Rose.  Measuring transport coefficients necessary for the description of coupled two-phase flow of immiscible fluids in porous media , 1988 .

[18]  L. Schwartz.  Théorie des distributions , 1966 .

[19]  C. C. Mattax,et al.  The Modeling of a Three-Dimensional Reservoir with a Two-Dimensional Reservoir Simulator-The Use of Dynamic Pseudo Functions , 1973 .

[20]  T. Kortekaas.  Water/Oil Displacement Characteristics in Crossbedded Reservoir Zones , 1985 .

[21]  L. Durlofsky.  Numerical calculation of equivalent grid block permeability tensors for heterogeneous porous media , 1991 .

[22]  Visualisation des propriétés capillaires des roches réservoir , 1982 .

[23]  John C. Martin.  Partial Integration of Equations of Multiphase Flow , 1968 .

[24]  J. R. Macmillan,et al.  Stochastic analysis of spatial variability in subsurface flows: 2. Evaluation and application , 1978 .

[25]  H. L. Stone,et al.  Rigorous Black Oil Pseudo Functions , 1991 .

[26]  Jacob Bear,et al.  Effective and relative permeabilities of anisotropie porous media , 1987 .

[27]  Benoît N∄tinger,et al.  The effective permeability of a heterogeneous porous medium , 1994 .

[28]  Charles A. Kossack,et al.  Scaling up heterogeneities with pseudofunctions , 1990 .

[29]  Larry W. Lake,et al.  Reservoir Characterization II , 1991 .

[30]  T. Hewett,et al.  A review of current trends in petroleum reservoir description and assessment of the impacts on oil recovery , 1992 .

[31]  A. E. Sáez,et al.  The effective homogeneous behavior of heterogeneous porous media , 1989 .

[32]  P. King.  The use of renormalization for calculating effective permeability , 1989 .

[33]  J. D. Huppler.  Numerical Investigation of the Effects of Core Heterogeneities on Waterflood Relative Permeabilities , 1970 .

[34]  G. Dagan.  Flow and transport in porous formations , 1989 .

[35]  A. S. Emanuel,et al.  Pseudo-relative permeability for well modeling , 1974 .

[36]  Two-phase flow in heterogeneous porous media II: Numerical experiments for flow perpendicular to a stratified system , 1990 .

[37]  H. S. Price,et al.  Flow in Heterogeneous Porous Media , 1961 .

[38]  R. Ababou,et al.  Three-dimensional flow in random porous media , 1988 .

[39]  Rachid Ababou,et al.  Random Porous Media Flow on Large 3-D Grids: Numerics, Performance, and Application to Homogenization , 1996 .

[40]  S. Whitaker.  Flow in porous media I: A theoretical derivation of Darcy's law , 1986 .

[41]  J. Auriault.  Nonsaturated deformable porous media: Quasistatics , 1987 .

[42]  C. L. Hearn,et al.  Simulation of stratified waterflooding by pseudo relative permeability curves , 1971 .

[43]  Aristotelis Mantoglou,et al.  Three-dimensional unsaturated flow in heterogeneous systems and implications on groundwater contamination: A stochastic approach , 1989 .

[44]  Michel Quintard,et al.  Two-phase flow in heterogeneous porous media I: The influence of large spatial and temporal gradients , 1990 .

[45]  S. Kirkpatrick.  Percolation and Conduction , 1973 .

[46]  Dominique Salin,et al.  Large-scale percolation theory of drainage , 1993 .

[47]  K. H. Coats,et al.  Simulation of three-dimensional, two-phase flow in oil and gas reservoirs , 1967 .

[48]  W. Rose.  Data interpretation problems to be expected in the study of coupled fluid flow in porous media , 1989 .

[49]  G. Chilingar,et al.  Oil and gas production from carbonate rocks , 1972 .