Simulation of land subsidence using finite element method: Rafsanjan plain case study

Land subsidence caused by the excessive use of groundwater resources has caused serious damage to Rafsanjan area. In this study, using finite element method, a 2D plain strain simulation of land subsidence has been conducted. A linear elastic constitutive law has been used for the simulation of the soil material skeleton. Actual water level during the analysis period has been modeled via specifying nodal water pressure at piezometer wells situation. The solution procedure consists of two parts. First, an initial static analysis is carried out in order to find initial steady-state solution for the pore pressure and stress distribution. Then, the above solution is used as initial condition for dynamic computation of consolidation equations during pumpage period. Tectonic effect has been considered as a rigid body motion. Numerical results showed that if the rate of pumpage remains constant in the future, settlement due to water withdrawal near Rafsanjan city will reach up to 110 cm by the year 2022.