STOCHASTIC FINITE-ELEMENT ANALYSIS FOR MULTIPHASE FLOW IN HETEROGENEOUS POROUS-MEDIA

Authors
Citation
R. Ghanem et S. Dham, STOCHASTIC FINITE-ELEMENT ANALYSIS FOR MULTIPHASE FLOW IN HETEROGENEOUS POROUS-MEDIA, Transport in porous media, 32(3), 1998, pp. 239-262
Citations number
26
Categorie Soggetti
Engineering, Chemical
Journal title
ISSN journal
01693913
Volume
32
Issue
3
Year of publication
1998
Pages
239 - 262
Database
ISI
SICI code
0169-3913(1998)32:3<239:SFAFMF>2.0.ZU;2-B
Abstract
This study is concerned with developing a two-dimensional multiphase m odel that simulates the movement of NAPL in heterogeneous aquifers. He terogeneity is dealt with in a probabilistic sense by modeling the int rinsic permeability of the porous medium as a stochastic process. The deterministic finite element method is used to spatially discretize th e multiphase flow equations. The intrinsic permeability is represented in the model via its Karhunen-Loeve expansion. This is a computationa lly expedient representation of stochastic processes by means of a dis crete set of random variables. Further the nodal unknowns, water phase saturations and water phase pressures, are represented by their stoch astic spectral expansions. This representation involves an orthogonal basis in the space of random variables. The basis consists of orthogon al polynomial chaoses of consecutive orders. The relative permeabiliti es of water and oil phases, and the capillary pressure are expanded in the same manner, as well. For these variables, the set of determinist ic coefficients multiplying the basis in their expansions is evaluated based on constitutive relationships expressing the relative permeabil ities and the capillary pressure as functions of the water phase satur ations. The implementation of the various expansions into the multipha se flow equations results in the formulation of discretized stochastic differential equations that can be solved for the deterministic coeff icients appearing in the expansions representing the unknowns. This me thod allows the computation of the probability distribution functions of the unknowns for any point in the spatial domain of the problem at any instant in time. The spectral formulation of the stochastic finite element method used herein has received wide acceptance as a comprehe nsive framework for problems involving random media. This paper provid es the application of this formalism to the problem of two-phase how i n a random porous medium.