Regional land subsidence simulation in Su-Xi-Chang area and Shanghai City, China

Su-Xi-Chang area and Shanghai City, located in the south of Yangtze Delta, China, has subsided due to groundwater overpumping. Because of the regional scale of the groundwater exploitation, cone of depression and land subsidence at present, Su-Xi-Chang area and Shanghai City are treated as a single area for land subsidence study to avoid the uncertainty of boundary condition due to the regionalism. The characteristics of aquifer system compaction are complex because of the difference in the types, compositions and structures of the soils that the hydrostratigraphic units are composed of, and in the histories of groundwater level change the hydrostratigraphic units have experienced. 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, elasto-plasticity, and visco-elasto-plasticity, at different sites of the cone of depression or in different periods, corresponding constitutive laws have been adopted. This avoids the shortcomings of the previous research that the same constitutive law was adopted in all the hydrostratigraphic units during the entire time period. A coupled flow and subsidence model, which includes a three-dimensional flow model with variable coefficients and a one-dimensional (vertical) subsidence model, is built according to the complicated hydrological condition in the region. The simulation model is calibrated using observed data, which include compression of individual strata from groups of extensometers and groundwater levels from observation wells from 1995 to 2002. The model reproduced that the primary subsidence layer in Shanghai shifts from the shallow aquitard to the fourth confined aquifer because of the groundwater yield variations and the change of exploitation aquifers. However the third aquitard was the primary subsidence layer in Su-Xi-Chang area and the compaction deformation of the sandy aquifers was remarkable. The simulation results could provide some reasonable advice about groundwater exploitation in the future.

[1]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[2]  Shujun Ye,et al.  Characteristics of aquifer system deformation in the Southern Yangtse Delta, China , 2007 .

[3]  Shujun Ye,et al.  Application of the multiscale finite element method to flow in heterogeneous porous media , 2004 .

[4]  J. Pacheco,et al.  Delimitation of ground failure zones due to land subsidence using gravity data and finite element modeling in the Querétaro valley, México , 2006 .

[5]  Giuseppe Gambolati,et al.  Mathematical simulation of the subsidence of Venice: 2. Results , 1974 .

[6]  T. R. Shearer.  A numerical model to calculate land subsidence, applied at Hangu in China , 1998 .

[7]  Donald C. Helm,et al.  One‐dimensional simulation of aquifer system compaction near Pixley, California: 2. Stress‐Dependent Parameters , 1976 .

[8]  Nonlinear Modeling of Groundwater Flow and Total Subsidence of the Mexico City Aquifer-Aquitard System , 1991 .

[9]  Xiaoqing Shi,et al.  Characterization of land subsidence induced by groundwater withdrawals in Su-Xi-Chang area, China , 2007 .

[10]  J. Jiao,et al.  Land subsidence caused by groundwater exploitation in Suzhou City, China , 2003 .

[11]  Thomas Y. Hou,et al.  Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients , 1999, Math. Comput..

[12]  Donald C. Helm,et al.  One‐dimensional simulation of aquifer system compaction near Pixley, California: 1. Constant parameters , 1975 .

[13]  Herman Bouwer,et al.  Land Subsidence and Cracking Due to Ground‐Water Depletiona , 1977 .

[14]  Viscoelastic aquifer model applied to subsidence due to pumping , 1977 .

[15]  S. P. Neuman,et al.  Adaptive explicit‐implicit quasi three‐dimensional finite element model of flow and subsidence in multiaquifer systems , 1982 .

[16]  Giuseppe Gambolati,et al.  Mathematical Simulation of the Subsidence of Ravenna , 1991 .

[17]  E. Frind,et al.  Hydraulic response of highly compressible aquitards during consolidation , 1991 .

[18]  Miguel A. Mariño,et al.  Prediction of optimal safe ground water yield and land subsidence in the Los Banos-Kettleman City area, California, using a calibrated numerical simulation model , 2001 .

[19]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[20]  Xiaoqing Shi,et al.  Characterization of regional land subsidence in Yangtze Delta, China: the example of Su-Xi-Chang area and the city of Shanghai , 2008 .

[21]  Giuseppe Gambolati,et al.  Mathematical simulation of the subsidence of Venice: 1 , 1973 .

[22]  J. F. Poland,et al.  LAND SUBSIDENCE DUE TO WITHDRAWAL OF FLUIDS , 1969 .