Simulation of Groundwater Flow in Fractured-Karst Aquifer with a Coupled Model in Maling Reservoir, China

A coupled model has been developed to simulate groundwater flow in fractured karst systems according to the complex geological and karst hydrogeological conditions of the dam site, where a 3D mathematical model based on Boussinesq equation was used to describe the movement of groundwater flow in fractured medium, and a 1D conduit model for karst medium. The model was solved with the continuous hydraulic heads at the common boundaries. The hydraulic conductivities of karst medium were determined by geometrical parameters and flux of pipes. Furthermore, the permeability parameters for fractured medium were calibrated by the measured and calculated groundwater levels. The calibrated model was employed to predict the variation of groundwater flow field and leakage from the karst pipes and underground powerhouse during the reservoir operation. The simulated results showed that the groundwater level of the powerhouse had decreased by about 2–5 m. The water level of conveyance pipeline had risen by 10–20 m, and the water level on both banks had risen by 15–25 m. The leakage of karst conduits for impervious failure was larger than that for normal seepage control. In addition, the leakage of the powerhouse was estimated to be about 1000–3000 m3/d, and the seepage control of karst pipes had little influence on the leakage of underground powerhouse.