Numerical modeling of land subsidence induced by groundwater extraction considering unsaturated effects and using element-free Galerkin (EFG) method

Document Type : Article

Authors

1 Geotechnical Engineering, Civil Engineering Department, Sharif University of Technology

2 Civil Engineering Department,, Khajeh Nasir Toosi University of Technology

Abstract

This study aims to provide a coupled flow-deformation (hydromechanical) model for simulating land subsidence associated with groundwater extraction in aquifers. For this simulation, we have adopted the element-free Galerkin (EFG) method and considered the unsaturated effects in the aquifers based on the aquifer's hydrologic and geotechnical characteristics. This model gives us a better understanding of the aquifer's geological characteristics, enabling us to forecast changes in the hydraulic head and land subsidence.


To ensure the credibility of our model and validate the code, we modeled unsaturated hydromechanical benchmark problems. Then, using the EFG method as a numerical tool, we modeled an isotropic aquifer to investigate the effects of groundwater pumping on land subsidence and hydraulic changes in the aquifer. To ascertain the reliability of the modeling, we compared EFG results with finite element method (FEM) models. The comparison results were satisfying, and the land subsidence and hydraulic head profiles demonstrate that the EFG method is capable of land subsidence simulation caused by water pumping. Furthermore, our findings highlight the nonlinear correlation between groundwater extraction and the subsequent decrease in hydraulic head and land subsidence augmentation. Finally, we performed a parametric study to better understand the effect of various characteristics of aquifers and observe the effect of the aquifer's parameters, such as hydraulic conductivity, elastic modulus, and Poisson's ratio. We investigated the effect of each parameter on land subsidence increase and hydraulic head decline. The results show that elastic modulus and Poisson's ratio have the most significant effect on land subsidence. Although hydraulic conductivity controls the hydraulic decrement and land subsidence increase time, it slightly affects the final hydraulic head and land subsidence at the steady-state stage. These results highlight the importance of in-situ measurement of elastic modulus and Poisson's ratio parameters with acceptable accuracy for groundwater extraction projects, as these parameters play a significant role in the feasibility studies.

Keywords

Main Subjects