Abstract

The regional movement of shallow groundwater in the fractured rock aquifer is examined through a conceptual-deterministic modeling approach. The computer program FRACNET represents the fracture zones as straight laminar flow conditions in connection to regional constant head boundaries within an impermeable rock matrix. Regional scale fracture zones are projected onto the horizontal plane, invoking the Dupuit Forchheimer assumption for flow. The steady state flow solution for the two dimensional case is achieved by requiring nodal flow balances using a Gauss-Seidel iteration. Computer experiments based on statistically generated fracture networks demonstrate the emergence of preferred flow paths due to connectivity of fractures to sources or sinks of water, even in networks of uniformly distributed fractures of constant length and aperture. The implication is that discrete flow, often associated with the local scale, may maintain itself even at a regional scale. The distribution of uniform areal recharge is computed using the Analytic Element Method, and then coupled to the network flow solver to complete the regional water balance. The areal recharge weakens the development of preferential flow pathways. The possible replacement of a discrete fracture network by an equivalent porous medium is also investigated. A Mohr's circle analysis is presented to characterize the tensor relationship between the discharge vector and the piezometric gradient vector, even at scales below the representative elementary volume (REV). A consistent permeability tensor is sought in order to establish the REV scale and justify replacement of the discrete fracture network by an equivalent porous medium. Finally, hydrological factors influencing the chemical dissolution and initiation of conduits in carbonate (karst) terrain are examined. Based on hydrological considerations, and given the appropriate geochemical and hydrogeological conditions, the preferred flow paths are expected to develop with time into caves.