Numerical Simulation Of Fluid Flow In Porous/Fractured Media

INTRODUCTION Our growing concern for adequate and secure sources of energy and minerals has stimulated vigorous exploration for new sources, research toward a better understanding of geological processes, and development of new extraction technologies. The need for control, or at least prediction, of subsurface fluid flow is important for many of these technologies for example: primary, secondary and tertiary recovery of oil; ground water and waste management; in-situ fossil energy extraction (oil shale, coal, tar sands); and solution mining of uranium, copper and other minerals. These technologies, especially the last two, are characterized by highly complex systems. A partial list of the physical processes occurring would include flow in porous/fractured media, multi-phase and multi-component flow with heat transport, chemically active fluids and soils, tracers, diffusion and dispersion, fracturing and dissolution. A great deal of understanding of how these processes behave and interact can be obtained from models. MODELS Theoretical models are valuable because they: 1) provide a frame of reference for interpreting results of laboratory and field experiments; 2) once validated by experiments, allow a variety of geometries, injection/production strategies, etc., to be examined (at relatively low cost) for efficiency and stability; 3) can provide guidance to the design of experiments and field operations. Most models are based on the fundamental principles of mass, momentum and energy balance. But from this starting point, many paths can be taken. For example, there is the theory of Payatakes (1973) which concentrates on the microscale dynamics. In this approach, the rock or soil matrix is represented by a complex of characteristic channels (such as periodically constricted smooth tubes). Detail of flow within the characteristic channel is calculated very accurately. A difficulty with this model is that description of a representative channel can require several parameters which may not be easily measured. Also, it is not clear how the channels can be combined practically to model a large scale flow. Another type of model is the "global" one described by Bear (1972). Here, the continuum equations for conservation of mass, momentum and energy are averaged over a distribution of pore channels, resulting in a set of conservation equations in which the small scale structure of the medium is replaced by quantities such as porosity, permeability, dispersion coefficients and tortuosity. This approach allows computation of large scale flows. However, additional constitutive equations are needed which relate averaged quantities such as permeability to observables such as local saturation, particle size distributions, and others. An important difference between models is the way they handle the momentum equation. Payatakes' model solves the full equation. In others, it is replaced by a simpler relation such as Darcy's law (valid for slow flow rates) or by Forchheimer's equation which extends Darcy's law to higher flow rates. These simple relations can nevertheless match a great deal of experimental data (Dullien 1975). The permeability term which appears in Darcy's law and in Forchheimer's expression has been related to other quantities such as porosity and particle surface area. The Kozeny-Carmen equation is a well-known example, valid for some ranges of porosity and particle sizes and for some materials. Several other semi-theoretical, semi-empirical formulations have been devised, none of which are entirely satisfactory. One goal of researchers has been the ability to predict permeability of a porous material from basic measurable quantities such as grain size distributions without recourse to adjustable parameters. This effort has proceeded from consideration of distributions of idealized, non-intersecting channels, to intersecting channels, to consideration of both distributed pore radii and pore neck radii. This last approach has been used successfully to predict permeabilities (ranging over several orders of magnitude) in sandstones (Dullien 1975). Additionally, studies on explicit networks of channels (e.g., Fatt 1956) and percolation theory (e.g., Larson et. al. 1981) have been used to examine interconnectivity effects in porous media with the hope of eventually being able to predict permeability accurately. To be of more than theoretical interest, a model must be compared s controlled laboratory conditions where boundary and initial conditions and relevant material properties can be accurately determined. Also, extraneous forces can be eliminated so that the interactions between processes of interest can be clearly seen and flow models can be vigorously tested. In contrast to the well-defined environment and