Theory and high-resolution finite element analysis of 2-D and 3-D effective permeabilities in strongly heterogeneous porous media

S. P. Neuman, S. Orr, O. Levin, E. Paleologos

Research output: Contribution to conferencePaper

23 Scopus citations

Abstract

Starting from the premise that Darcy's law applies with a random hydraulic conductivity field K(x) defined on a support ω, we consider steady state flow in a larger domain Ω driven by statistically independent random source and boundary functions. Writing K(x) as the sum of a slowly varying unbiased estimate κ(x) and a zero mean estimation error K′(x), the conditional ensemble moments 〈h(x)〉κ≡〈h(x)|κ(x)〉 and 〈q(x)〉κ≡〈q(x)|κ(x)〉 become conditionally unbiased predictors of the hydraulic head h(x) and Darcy flux q(x), respectively; 〈q(x)〉κ satisfies a standard continuity equation driven by ensemble mean source and boundary functions; and 〈q(x)〉κ = -κ(x)▽〈h(x)〉κ+rκ(x) where rκ(x) is a residual flux. We present an exact formal expression for 〈r(x)〉κ which demonstrates that 〈q(x)〉κ is nonlocal (depends on head gradients at points other than x) with a well-defined kernel, and non-Darcian (there is no effective or equivalent conductivity valid for arbitrary directions of conditional mean flow), except in special cases. In some of these cases the effective conductivity is κ(x), in some it is a symmetric or a nonsymmetric second-rank tensor, and in yet other situations it exists only as a set of directional scalars but not as a tensor. We describe a weak approximation for 〈r(x)〉κ which improves with the quality of the estimate κ(x), and an ancillary approximation for the effective hydraulic conductivity tensor in an infinite domain. We then use high-resolution finite element Monte Carlo simulation to verify aspects of these approximations under uniform 3-D and radial 2-D flows in strongly heterogenous, statistically homogeneous and isotropic porous media.

Original languageEnglish (US)
Pages117-136
Number of pages20
StatePublished - Jan 1 1992
EventProceedings of the 9th International Conference on Computational Methods in Water Resources - Denver, CO, USA
Duration: Jun 1 1992Jun 1 1992

Other

OtherProceedings of the 9th International Conference on Computational Methods in Water Resources
CityDenver, CO, USA
Period6/1/926/1/92

    Fingerprint

ASJC Scopus subject areas

  • Engineering(all)

Cite this

Neuman, S. P., Orr, S., Levin, O., & Paleologos, E. (1992). Theory and high-resolution finite element analysis of 2-D and 3-D effective permeabilities in strongly heterogeneous porous media. 117-136. Paper presented at Proceedings of the 9th International Conference on Computational Methods in Water Resources, Denver, CO, USA, .