Numerical simulation of land subsidence induced by groundwater overexploitation in Su-Xi-Chang area, China

Su-Xi-Chang area is one of the typical regions in China which suffers from severe land subsidence. Various field monitoring records were integrated to study the characteristics and mechanisms of land subsidence in this region. The development of the land subsidence in this region shows a tight spatial and temporal correlation with the groundwater pumping. Based on the analysis of the field data, it is found that the deformation patterns of the hydrogeologic units are greatly related to the hydrogeologic properties and groundwater level variations. Some have an elastic behavior, others may have an elastic–plastic rheology. Hence, a 3D finite element numerical model considering the rheological properties of the soil was developed to simulate the groundwater level and land subsidence. Both hydraulic conductivity and specific storage were expected to vary with the porosity during the process of consolidation. Multiscale finite element method (MsFEM) was applied to solve the model during the period from 1996 to 2004. After calibrating the model with the observed groundwater level and subsidence data, the parameters of the multi-layers system were estimated. The calibrated model outputs fit reasonably well with the observed data. Consequently the model can be applied to predict groundwater level and land subsidence in future pumping scenarios. The model predictive results show that land subsidence rate can be controlled and even rebound may occur after the implementation of the groundwater exploitation prohibition.

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

[2]  S. Leake.  Interbed storage changes and compaction in models of regional groundwater flow , 1990 .

[3]  Jyh-Jong Liao,et al.  Characterization of land subsidence in the Choshui River alluvial fan, Taiwan , 2004 .

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

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

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

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

[8]  C. F. Lee,et al.  Characteristics and trends of land subsidence in Tanggu, Tianjin, China , 2002 .

[9]  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 .

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

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

[12]  Wallace E. Pratt,et al.  Local Subsidence of the Goose Creek Oil Field , 1926, The Journal of Geology.

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

[14]  M. Sneed,et al.  Aquifer-System compaction and land subsidence: Measurements, analyses, and simulations-the Holly Site, Edwards Air Force Base, Antelope Valley, California , 2000 .

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

[16]  Devin L. Galloway,et al.  Land subsidence in the United States , 1999 .

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

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