No abstract

The Xiluodu Dam is a concrete double-curvature arch dam with a crest elevation of 610 m and a height of 285.5 m. Since the impoundment of the Xiluodu reservoir, remarkable river valley contractions (RVCs) have been observed upstream and downstream of the reservoir, potentially threatening the safety of the dam. However, the cause of these RVCs remains unclear. Based on an analysis of hydrogeological conditions, the RVCs were determined a result of the expansion of the aquifer, within which the effective stress decreased due to an increase in the hydraulic head after reservoir impoundment. Referring to the hydrostatic seasonal time (HST) model, a groundwater hydrostatic seasonal (GHS) model is proposed for simulating and predicting the development of the RVCs. Unlike the HST model, the GHS model can provide information on aquifer hydraulic diffusivity. The calibration results illustrate that the GHS model can accurately fit the observed RVCs data. The calculation results revealed that the RVCs were mainly affected by the hydraulic head of the confined aquifer, and that seasonal effects gave rise to less than 10% of the total RVCs. Finally, the development of RVCs were predicted using the GHS model. The prediction results demonstrated that the RVCs of most monitoring lines in the Xiluodu reservoir would gradually approach a convergence condition after 6 February 2021. Until the deadline of the prediction on 1 May 2035, there is still one monitoring line that has not reached a convergence condition (whose RVCs are 157.6 mm, and where the RVC growth rate will decrease to 0.005 mm/d by that time). Considering the large amount of RVCs, we think the safety of the dam requires closer consideration.

The response laws of groundwater dynamics on the riverbank to river level variations are highly dependent on the river level fluctuation process. Analytical solutions are widely used to infer the groundwater flow behavior. In analytical calculations, the river level variation is usually generalized as instantaneous uplift or stepped, and then the analytical solution of the unsteady groundwater flow in the aquifer is derived. However, the river level generally presents a complex, non-linear, continuous change, which is different from the commonly used assumptions in groundwater theoretical calculations. In this article, we propose a piecewise-linear approximation to describe the river level fluctuation. Based on the conceptual model of the riverbank aquifer system, an analytical solution of unsteady groundwater flow in an unconfined aquifer under complex boundary conditions is derived. Taking the Xiluodu Hydropower Station as an example, firstly, the monitoring data of the river level during the period of non-impoundment in the study area are used to predict the groundwater dynamics with piecewise-linear and piecewise-constant step approximations, respectively, and the long-term observation data are used to verify the calculation accuracy for the different mathematical models mentioned above. During the reservoir impoundment period, the piecewise-linear approximation is applied to represent the reservoir water level variation, and to predict the groundwater dynamics of the reservoir bank.

Although specific storage (S S ) is a crucial parameter for characterizing the dynamic water-release process in an aquitard, the Hvorslev analysis method is only capable of estimating the hydraulic conductivity (K), rather than S S . The Hvorslev method is based on an assumption of steady-state flow and uses the falling-head test. Considering the elastic water-release effect in an aquitard and the actual groundwater flow, two mathematical unsteady-state groundwater flow models (i.e., the constant-head and variable-head models), describing the head distribution are presented for a constant-head injection test (CHIT) and variable-head injection test (VHIT) performed in a multi-layered aquifer-aquitard system. The analytical solutions for the models are derived by Fourier transforms to estimate K and S S . The obtained K and S S values from the CHIT and VHIT were generally consistent, although the VHIT is recommended in practice. For the VHIT, estimating K using the Hvorslev analysis method under steadystate groundwater flow conditions resulted in a small relative error compared to the obtained result from the variable-head model. It was demonstrated that a certain relation existed between the constant-head and variable-head methods when the drawdown rate was lower than 1%, and so the constant-head method was applied to the VHIT data to determine K and S S . The synthetic error in the parameters was smaller than 10%, which indicated that the methods could be used interchangeably under these circumstances, thus simplifying the hydraulic parameter estimation process.