Document Type : Research Article
Authors
1 Director of Kusar Dam, Kohgiluyeh and Boyer-Ahmad Regional Water Authority
2 Faculty of Earth Sciences, Shahid Chamran University of Ahvaz, Ahvaz, Iran
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 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.
Keywords
Main Subjects
)
Send comment about this article