Experimental-numerical simulation of soluble formations in reservoirs
Introduction
Water storage in reservoirs usually increases the water residence time that can lead to negative impacts on water quality (Kadivar, 2007; Noori et al., 2018; Noori et al., 2019a; Aradpour et al., 2021). Reservoir water salinization (RWS) emerges when soluble geological formations outcrop within the reservoir area, a relevant threat in many regions of Iran (Naderkhanloo et al., 2017) that is poorly understood. A recent example of this problem is the salinization of the Gotvand dam reservoir located in the southwest of Iran due to the presence of halite in the reservoir location, which has caused serious concerns about the operation of this dam in the future (Naderkhanloo et al., 2017; Jalali et al., 2019).
The Parsian dam, as a hydropower dam, is currently under the study phase on the Fahlian river that flows through gypsum-type formation, called Gachsaran geological formation (GGF) in Iran. Limestone compounds have been expanded in some parts of the reservoir and their main component is gypsum (CaSO4.2H2O). Considering the high solubility of CaSO4.2H2O, there are serious concerns regarding the RWS when the Parsian dam reservoir (PDR) is impounded (Baghdardokht and Heidari, 2005). The problem arises when the GGF dissolution rate is not fully understood. On the other hand, it is not recommended to use the suggested dissolution rates for different gypsum-type geological formations due to the high range of uncertainty reported around the world (Raines and Dewers, 1997; Klimchouk and Aksem, 2002; Baghdardokht and Heidari, 2005; Aljubouri and Al-Kawaz, 2007; Mbogoro et al., 2011; Valor et al., 2011; Lebedev, 2015; Domínguez-Villar et al., 2017; Feng et al., 2017; Hong et al., 2018; Tang et al., 2018). Therefore, attempts to understand the rock-water interaction and the GGF dissolution rates in the PDR should be made through laboratory tests or field measurements. Saline inflow, mainly from the Hamzeh-beigi reach (a tributary of the Fahlian river), can also contribute to the RWS once the PDR is impounded. Therefore, the question raised here is how great the effect of GGF dissolution outcropped in the reservoir would be on the RWS. The next question centers around the impact of Hamze-beigi saline inflow and its influence on the RWS. Since there is a hydroelectric power plant in Parsian dam, its activation in the operation phase makes constant oscillation at the higher levels of the reservoir which leads to increasing the vertical velocity at the top layers. Therefore, it is important to address the effect of water oscillation on RWS in top layers, where the water is released for hydropower generation. Answering these questions can help decision-makers to choose the optimum operation policy and consider the possible change in dam location to avoid the severe RWS in practice.
To answer the questions raised above, the salinity distribution should be simulated across the reservoir under different initial and boundary conditions using mathematical models. In this study, we explored the GGF dissolution rates in the PDR through laboratory tests using a physical model, under different conditions of flow velocity that exactly take place in different depths of Parsian reservoir. Then, we connected the physical model's results of the specified dissolution rates to CE-Qual-W2 (W2) model to further explore the RWS under different scenarios of reservoir operation. Actually, the definition of dissolution rate in the W2 is not possible simply, and it needs to improve the model's ability using a creative technique. Here, we fed the GGF dissolution as a pollution source through virtual tributary inflows with low flows and high concentrations of total dissolved solids (TDS) to the W2 model. It should be noted that we selected the W2 model, a well-known reservoir hydrodynamic and water quality model (Cole and Wells, 2015), as it has been successfully used for water quality simulation in different reservoirs around the world (Kim and Kim, 2006; Ma et al., 2015; Noori et al., 2015; Zhang et al., 2015; Ji et al., 2017; Al-Zubaidi and Wells 2017; Sadeghian et al., 2018; Kim et al., 2019; Lindenschmidt et al., 2019; Mi et al., 2019, 2020; Al-Zubaidi and Wells, 2020; Morales-Marin et al., 2020). This work aims to contribute to future recommendation on dam construction to account for saline inflows in regions with soluble geological formations outcropped in the lake area to limit RWS.
Section snippets
Study area
The Parsian dam is planned to be constructed on the Fahlian river in the Fars province, south of Iran. The main purpose of the dam is hydropower generation as well as some secondary goals, including supplying water for agricultural, industrial and domestic purposes. The geographic location of the dam is 51º 40´ E and 29º 45´ N (Fig. 1 and Movie S1). According to a 43-year historical data, the average annual precipitation, air temperature, and potential evaporation in the study area,
TDS equilibrium
After pilot preparation, GSg and GSmg samples were placed in the pilot and water circulation was started. TDS was recorded every six hours from the test onset until the TDS variation was very low (i.e., equilibrium state). Results showed that 6–8 days of water circulation was enough to reach an equilibrium state. Due to the number of experiments (32 tests) and the two available pilots, conducting the tests lasted more than four months. The TDS changes in the experimental tanks under different
Conclusion
The Parsian dam is under the study phase in the south of Iran. Given Ghachsaran geological formation in the reservoir area and a saline tributary inflow to the dam, water salinization is a serious concern for reservoir impoundment. Due to the importance of this issue, a physical pilot was conducted and regular tests were performed to estimate the Gachsaran formation dissolution rate. In the next step, a hydrodynamic and water quality mathematical model (CE-QUAL-W2) was applied to simulate the
Data availability
The raw data used in this study are available via Data Archive of the Iran Water Resources Management Company http://wrs.wrm.ir/amar/login.asp.
CRediT authorship contribution statement
Negar Tavoosi: Data curation, Formal analysis, Software, Writing – original draft, Writing – review & editing. Farhad Hooshyaripor: Conceptualization, Methodology, Formal analysis, Writing – original draft, Writing – review & editing. Roohollah Noori: Conceptualization, Methodology, Supervision, Writing – original draft, Writing – review & editing. Ashkan Farokhnia: Methodology, Validation, Software, Visualization, Writing – review & editing. Mohsen Maghrebi: Formal analysis, Visualization,
Declaration of Competing Interest
None.
References (46)
- et al.
In situ nanoscale observations of gypsum dissolution by digital holographic microscopy
Chem. Geol.
(2017) - et al.
Salinization of reservoirs in regions with exposed evaporites. The unique case of upper Gotvand Dam, Iran
Water Res.
(2019) - et al.
Surface controlled dissolution rates of gypsum in aqueous solutions exhibit nonlinear dissolution kinetics
Geochim. Cosmochim. Acta
(2001) - et al.
Impacts of water level rise on algal bloom prevention in the tributary of three gorges reservoir, China
Ecol. Eng.
(2017) - et al.
Modeling density currents in a typical tributary of the three gorges reservoir, China
Ecol. Model.
(2015) - et al.
Ensemble warming projections in Germany's largest drinking water reservoir and potential adaptation strategies
Sci. Total Environ.
(2020) - et al.
Temporal and depth variation of water quality due to thermal stratification in Karkheh Reservoir, Iran
Journal of Hydrology: Regional Studies
(2018) - et al.
A reduced-order based CE-QUAL-W2 model for simulation of nitrate concentration in dam reservoirs
J. Hydrol.
(2015) - et al.
Mixed transport/reaction control of gypsum dissolution kinetics in aqueous solutions and initiation of gypsum karst
Chem. Geol.
(1997) - et al.
An empirical rate law for gypsum powder dissolution
Chem. Geol.
(2018)
Integration of a benthic sediment diagenesis module into the two dimensional hydrodynamic and water quality model – CE-QUAL-W2
Ecol. Model.
Dissolution rate of gypsum under different environments
Iraqi J. Earth Sci.
Analytical and field verification of a 3D hydrodynamic and water quality numerical scheme based on the 2D formulation in CE-QUAL-W2
J. Hydraul. Res.
3D numerical temperature model development and calibration for lakes and reservoirs: a case study
World Environ. Water Resour. Congr.
Water Quality for Agriculture
Metal contamination assessment in water column and surface sediments of a warm monomictic man-made lake: Sabalan Dam Reservoir, Iran
Hydrology Research
Alarming carcinogenic and non-carcinogenic risk of heavy metals in Sabalan dam reservoir, Northwest of Iran
Environ. Pollut. Bioavailab.
Dissolution study of karstic formation of Kangir Rservoir Dam
CE-QUAL-W2: A two-dimensional, Laterally averaged, Hydrodynamic and Water Quality Model. Version 3.72 User Manual
CE-QUAL-W2: A two-dimensional, Laterally averaged, Hydrodynamic and Water Quality Model.” Version 4.0
The role of gypsum and/or dolomite dissolution in tufa precipitation: lessons from the hydrochemistry of a carbonate–sulphate karst system
Earth Surf. Process. Landf.
An experimental study simulating the dissolution of gypsum rock
Energy Explor. Exploit.
Cited by (5)
A water quality database for global lakes
2024, Resources, Conservation and RecyclingApplication of Optimization Techniques for Searching Optimal Reservoir Rule Curves: A Review
2023, Water (Switzerland)