Quantify the impact of irregular aquifer geometry on transmissivity estimates: a numerical pumping test study
(Work in progress)

The main objective of this study is whether aquifers with irregular thickness lead to errors in the transmissivity estimates when analyzing its pumping test data. In addition, quantify the magnitude of transmissivity estimations errors caused by this irregular geometry for different scenarios.

Skills and Tools Used

  • Numerical modeling.
  • Pumping test data analisys.
  • Aquifer characterization.
  • Python.
  • OpenGeoSys.
  • Gmsh.

Introduction to the Project

It is of great importance to have knowledge of how the hydraulics of an aquifer in which quality use will be projected, which ease and remediation of improvements to the most efficient areas. A pumping test is a commonly used technique for aquifer investigation. With this, it is possible to obtain data on the evolution of the water level drawdown curve over time and use them to characterize the hydraulic properties of the aquifer. Among the conventional methods, the most famous are Cooper-Jacob [1946] and the use of the analytical solution of Theis [1935] for the estimation of the transmissivity (T) and the calculation of the storativity (S). These conventional methods assume that the aquifers are homogeneous, that is, they have constant T and S. However, the vast majority of aquifers are heterogeneous.

To overcome the limitations of conventional methods in the analysis of heterogeneous aquifers, some techniques were used, such as the use of data from many observation wells during a pumping test in its stationary stage to create hydrographs and analyze the heterogeneity [Wen et al., 2010]. Methods were developed to characterize the geostatistical hydraulic parameters of a heterogeneous aquifer, for example by using an inverse geostatistical method in Hydraulic Tomography [Yeh and Liu, 2000] and by expanding the Theis and Thiem solution for heterogeneous aquifers with log-normal distribution [Zech et al., 2012; and Zech et al., 2016]. There are several theoretical studies investigating the characterization of the heterogeneity of an aquifer using numerical models [Wu et al., 2005] and sandbox models [Berg and Illman, 2011].

Investigating methods for characterizing aquifer heterogeneity Wu et al. [2005] used a heterogeneous T field and Berg and Illman [2011] a sandbox with layers of different hydraulic conductivity (K). In addition, analytical methods assume that the flow of water during a pumping test is laminar and horizontal. However, there are several aquifers of complex geometry [e. g., Gutentag et al., 1984; Viguier et al., 2018] that during a pumping test they can present vertical components of water flow and can cause erroneous estimates of the T of these aquifers.

With the hypothesis that aquifers of complex geometry can cause wrong estimates of T from the analysis of pumping test data, this project proposes to investigate possible aquifer geometries that result in a deviation between the estimated T values and the actual T values of an aquifer from the use of analytical solutions for heterogeneous aquifers since the variation of the thickness of an aquifer (complex geometry) characterizes it as a heterogeneous T aquifer. Therefore, the approach of this project will be through the use of numerical models for the creation of a field of heterogeneous thickness of aquifer of two dimensions (a horizontal and a vertical axis). It will be able to provoke water flows in the porous medium with vertical components during a pumping test. This could lead to a divergence of the estimated T values of the aquifer when compared to the actual T values from the use of analytical solutions. With this, it will be possible to measure the magnitude of the transmissivity estimation errors caused by the irregularity of the thickness of a confined aquifer.

References

Berg, S. J., and W. A. Illman, 2011. Capturing aquifer heterogeneity: Comparison of approaches through controlled sandbox experiments, Water Resour. Res., 47, W09514, doi:10.1029/2011WR010429.

Cooper, H. H., Jr., and C. E. Jacob, 1946. A generalized graphical method for evaluating formation constants and summarizing well-field history, Eos Trans. AGU, 27(4), 526–534.

Gutentag, E. G., F. J. Heimes, N. C. Krothe, R. R. Luckey, and J. B. Weeks. 1984. Geohydrology of the High Plains Aquifer in parts of Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming. U.S. Geological Survey Professional Paper 1400-B.

Theis, C. V. (1935), The relation between lowering the piezometric surface and the rate and duration of discharge of a well using ground water storage, Eos Trans. AGU, 16, 519–524.

Wen, J. C., C. M. Wu, T.-C. J. Yeh, and C. M. Tseng (2010), Estimation of effective aquifer hydraulic properties from an aquifer test with muti-well observations (Taiwan), Hydrogeol. J., 18, 1143–1155, doi:10.1007/ s10040-010-0577-1

Wu, C. M., T.-C. J. Yeh, J. Zhu, T. H. Lee, N. S. Hsu, C. H. Chen, and A. F. Sancho (2005), Traditional analysis of aquifer tests: Comparing apples to oranges?, Water Resour. Res., 41, W09402, doi:10.1029/2004WR003717.

Yeh, T. C. J., & Liu, S. (2000). Hydraulic tomography: Development of a new aquifer test method. Water Resources Research, 36(8), 2095-2105. https://doi.org/10.1029/2000WR900114

Zech, A., S. Muller, J. Mai, F. Heße, and S. Attinger. 2016. `Extending Theis' solution: Using transient pumping tests to estimate parameters of aquifer heterogeneity. Water Resources Research 52, no. 8: 6156–6170. https://doi.org/ 10.1002/2015WR018509

Zech, A., C.L. Schneider, and S. Attinger. 2012. The extended Thiem’s solution Including the impact of heterogeneity. Water Resources Research 48, no. 10: W10535. https://doi .org/10.1029/2012WR011852