Numerical Simulation of Viscoelastoplastic Land Subsidence due to Groundwater Overdrafting in Shanghai, China

Abstract

Land subsidence caused by groundwater overpumping in Shanghai is becoming a serious geological hazard. Due to its important economic position, the field data, including compression of individual stratum from extensometer groups and the groundwater levels from observation wells, have been collected over the past 45 years. Considering the fact that different hydrostratigraphic units have different kinds of deformation and that an identical unit may also present different deformation characteristics, such as elasticity, elastoplasticity, and viscoelastoplasticity, at different sites of the cone of depression or in different periods, a nonlinear coupled regional land subsidence model is developed. The coupled model consists of a three-dimensional groundwater flow model and a one-dimensional vertical deformation model, both based on a viscoelastoplastic constitutive laws (called modified Merchant model), and then solved using a multiscale iterative finite-element method. The model is calibrated using 28,184 hy...