Irrigation
Akbar Khedri; Amir Saberinasr; Nasrollah Kalantari
Abstract
Introduction: The comprehension of the hydrogeological conditions of the aquifer and the determination of its hydraulic characteristics, such as hydraulic conductivity, transmissivity coefficient, and specific storage, are crucial for the management and preservation of groundwater resources. Various ...
Read More
Introduction: The comprehension of the hydrogeological conditions of the aquifer and the determination of its hydraulic characteristics, such as hydraulic conductivity, transmissivity coefficient, and specific storage, are crucial for the management and preservation of groundwater resources. Various conventional methods, including empirical formulas, laboratory techniques (constant and falling head), tracer tests, field tests (Lugeon, Lefranc, slug, flowmeter, and pumping tests), and groundwater inverse modeling, are employed to establish these characteristics, particularly hydraulic conductivity. Empirical formulas are limited to ideal conditions, and in laboratory methods, the sample must be kept undisturbed. Due to the impracticality of measuring large-scale effective factors, the hydraulic conductivity determined through laboratory methods is also the only representative of the hydraulic conductivity at the sampling point. Tracer studies encounter numerous constraints, such as time, cost, porosity determination, and tracer dispersion in multilayered aquifers. It is also difficult to determine the average hydrodynamic properties of the heterogeneous aquifer based on the data obtained from a specific section of the Lefranc and Slug tests. Consequently, pumping tests are commonly selected for hydraulic parameter estimation. Although costly and time-intensive, these tests provide more precise coefficients. Geophysical methods have been greatly developed during the last two decades and have shown a significant correlation with the hydraulic parameters of the aquifer derived from borehole pumping tests or direct laboratory measurements. This approach minimizes uncertainties in numerical model calibration, improves data coverage, and reduces the time and cost of regional hydrogeological investigations. The conventional approach, known as the electrical resistivity method, is still widely used in global and local research projects for evaluating aquifer hydraulic characteristics (Ige et al., 2018; Arétouyap et al., 2019; Youssef, 2020; Ullah et al., 2020; de Almeida et al., 2021; Lekone et al., 2023). Therefore, this study aims to use the integrated approach of the geophysical method and pumping test as a cost-effective and efficient alternative for estimating the hydraulic parameters of the alluvial aquifer in the northeast of Gachsaran city.
Material and Methods: The research area is an alluvial aquifer located 5 km to the northeast of Gachsaran, between coordinates 50-52 to 51-09 E longitude and 30-15 to 30-28 N latitude. Using 86 vertical electrical soundings, Archie's equations, and the IPI2win software, the hydraulic characteristics of the aquifer under investigation were estimated. Subsequently, these characteristics were then compared to the coefficients derived from the data of two pumping test wells, which were calculated using the Aquifer test software and obtained via the Cooper-Jacob and Neuman methods.
Results and discussion: The hydrodynamic coefficients of the aquifer were initially determined using the Cooper-Jacob method in this study. The hydraulic conductivity values for wells one and two are 4.9 m/day and 5.7 m/day, respectively. Correspondingly, the storage coefficient values for wells one and two are 0.015 and 0.021, respectively. Based on the Cooper-Jacob approach, it is deduced that if the storage coefficient values exceed 0.001, the aquifer is classified as unconfined. In this study, the storage coefficient values for both pumping wells suggest that the aquifer is unconfined. Since the vertical flow component and the delayed yield phenomenon should also be taken into account in unconfiend aquifers, the Neuman analytical model has been used in the studied aquifer. The values of specific yield (Sy) for pumping wells one and two, which are related to delayed yield, are 0.05 and 0.04, respectively. These values were calculated by analyzing the first segment of the curve derived from the Neuman logarithmic drawdown-time plot. The storage coefficient values for pumping wells one and two, extracted from the second section of the curve, are 0.015 and 0.021, respectively. Furthermore, the transmissivity value for well number 1 was 323 m2/day, while for well number 2, it was 655.5 m2/day. The vertical electrical sounding (VES) data were subsequently initially analyzed and interpreted using the IPI2win software and the equalization curve method (partial curve matching technique). The coefficients denoted as m and n, indicative of the degree of cementation of the sediments, were determined based on the sedimentary composition prevalent in the area. Archie's equations were employed to calculate the formation factor and porosity parameters. The aquifer exhibits a porosity range of approximately 0.15 in the eastern and southeastern parts (near the outlet of the plain) and around 0.41 in the centeral, northern, and northwestern sections of the area (next to the Asmari Formation). The specific yield (Sy) of the aquifer was calculated using the provided formula:
The minimum and maximum specific yield were estimated as 0.006 (in the eastern and southeastern regions) and 0.089 (in the western and northwestern regions of the plain), respectively, with an average value of 0.04. The transmissivity coefficients for the entire aquifer were then calculated based on the fitted relationship between hydraulic conductivity (K) and formation factor (F):
The range of transmissivity coefficients varies from a minimum of 63 m2/day (in the western and northwestern sections of the plain) to a maximum of 608.9 m2/day (in the eastern and southeastern areas). The average transmissivity coefficient is calculated as 323.7 m2/day.
To ensure the precision of the geoelectric method's coefficients, a comparative analysis was conducted with the hydrodynamic coefficients obtained from the two pumping test wells, as presented in the table below:
Conclusion: The evaluation and comparison of the hydrodynamic coefficients derived from the aforementioned methods indicate that the geoelectric method coefficients exhibit acceptable agreement with the pumping test coefficients. In other words, the analysis of the pumping test conducted using the Neuman technique in the unconfined aquifer revealed that well number two displayed a greater transmissivity coefficient, while well number one presented a higher specific yield. These findings are confirmed by the geoelectric approach. Consequently, such hybrid approaches, which include simultaneous analysis of geophysical methods (such as VES) and pumping tests will be a great alternative to multiple costly pumping tests for evaluating the hydrodynamic coefficients of an aquifer. Moreover, employing this hybrid technique enables the generation of dense hydrodynamic coefficients in an aquifer for use as inputs in the groundwater model.
R. Lalehzari; Saeid Boroomand Nasab; Hadi Moazed; A. Haghighi
Abstract
Introduction: Groundwater is the largest resource of water supplement and shortages of surface water supplies in drought conditions that requires an increase in groundwater discharge. Groundwater flow dependson the subsurface properties such as hydraulic gradient (water table gradient or head loss in ...
Read More
Introduction: Groundwater is the largest resource of water supplement and shortages of surface water supplies in drought conditions that requires an increase in groundwater discharge. Groundwater flow dependson the subsurface properties such as hydraulic gradient (water table gradient or head loss in artesian condition) and hydrodynamic coefficients. The flow treatment is analyzed with an accurate estimation of effective parameters in groundwater equation. This function is obtained using the continuous equation. Inlet and outlet flows of a cell are equal to storage amounts in the continuous equation. Analytical solution of this equation is complex, so numerical methods are developed including finite element and finite difference methods. For example, Feflow is a groundwater modeling tool that makesuse of finite element method (Reynolds and Marimuthu, 2007). Modflow as a finite difference three-dimensional model simulated underground flow under steady and unsteady conditions in anisotropic and non-homogeneous porous media. Modflow is designed to simulate aquifer systems in which saturated-flow conditions exist, Darcy’s Law applies, the density of groundwater is constant, and the principal directions of horizontal hydraulic conductivity or transmissivity do not vary within the system. In Modflow, an aquifer system is replaced by a discretized domain consisting of an array of nodes and the associated finite difference blocks. Groundwater modeling and water table prediction by this model have the acceptable results, because many different informations of water resource system are applied. Many people and organizations have contributed to the development of an effective groundwater monitoring system, as well as experimental and modeling studies (Lalehzari et al., 2013). The objective of this paper is investigation of hydraulic and physical conditions. So, a numerical model has to be developed by PMWIN software for Bagh-i Malek aquifer to calculate hydrodynamic coefficients and predict water table in the future.
Materials and Methods: Bagh-i Malek aquifer located in Khuzestan province is mainly recharged by inflow at the boundaries, precipitation, local rivers and return flows from domestic, industrial and agricultural sectors. The discharge from the aquifer is through water extraction from wells, springs, and qanats as well as groundwater outflow and evapotranspiration. In this study, conceptual model of Bagh-i Malek aquifer on the framework of finite difference numerical approach has been used in simulating groundwater flow treatment. Water table data of 8 piezometers was collected for the 10 year duration from 2002 to 2012. The study years are divided into 40 seasonal stress periods with daily time step. Hydraulic conductivity, specific yield and recharge were calibrated in these periods. Verification was made between the simulated and measured hydraulic heads in the next calibration year. To simulate the groundwater table elevation in this study area, the PMWIN model is used. Bagh-i Malek aquifer is considered as a single layered aquifer, and therefore only the horizontal hydraulic conductivity is estimated. Modflow was used to simulate both steady state and transient flow systems. In steady conditions it is assumed that the total of time simulation is a time period and it does not change inlet data in the modeling duration. In unsteady conditions,the duration of study is divided into some stress periods that data is changed in every period.
Results and Discussion: Estimation of hydraulic conductivity is the first step of calibration process at steady state conditions. The correct assignment of hydraulic conductivity has a main effect on other parameters accuracy. Hydraulic conductivity mapping indicated that the maximum values are in the Eastern North (6-7 m/day) of the aquifer. The twice calibrated parameter is specific yield in unsteady conditions. Specific yield mapping indicated that the values vary from 0.03 to 0.08 and are maximum in the Southern regions of the plain similar to hydraulic conductivity. The results confirm that the flow model has the tolerable simulation accuracy by variances of 3.1 and 3.84 in calibration and verification processes, respectively. The sensitivity of the flow model to decreasing the hydraulic conductivity is more than increasing it. 50 percentage declined into the hydraulic conductivity causes the increase of the variance from 3.1 of initial value to 44.
Conclusions: Mapping of calibrated hydraulic conductivity showed that the Eastern North of aquifer has the higher transmissivity and discharge capability in comparison to Southern parts. At last, the result show that the Bagh-i Malek aquifer model is sensitive to recharge, hydraulic conductivity and specific yield, respectively.
