Modelling karst aquifer evolution in fractured, porous rocks

The removal of material in soluble rocks by physical and chemical dissolution is an important process enhancing the secondary porosity of soluble rocks. Depending on the history of the soluble rock, dissolution can occur either along fractures and bedding partings of the rock in the case of a telogenetic origin, or within the interconnected pore space in the case of eogenetic origin. In soluble rocks characterised by both fractures and pore space, dissolution in both flow compartments is possible. We investigate the dissolution of calcite both along fractures and within the pore space of a limestone rock by numerical modelling. The limestone rock is treated as fractured, porous aquifer, in which the hydraulic conductivity increases with time both for the fractures and the pore spaces. We show that enlargement of pore space by dissolution will accelerate the development of a classical fracture-dominated telogenetic karst aquifer, breakthrough occurs faster. In the case of a pore-controlled aquifer as in eogenetic rocks, enlargement of pores results in a front of enlarged pore spaces migrating into the karst aquifer, with more homogeneous enlargement around this dissolution front, and later breakthrough. (C) 2016 Elsevier B.V. All rights reserved.