The state of land subsidence and prediction approaches due to groundwater withdrawal in China

This article gives a general introduction to land subsidence with the prediction approaches due to withdrawal of groundwater in three subsided/subsiding regions in China: the deltaic plain of Yangtse River (YRDP), North China Plain (NCP), and Fenwei Plain (FP). On YRDP, Shanghai is the typical subsided/subsiding city; on NCP Tianjin is the typical subsided/subsiding city, and on FP Taiyuan is the typical subsided/subsiding city. The subsided area with subsidence over 200 mm on YRDP is about 10,000 km² and the maximum subsided value reached 2.9 m at Shanghai; on NCP the subsided area reached 60,000 km² with the maximum subsidence of 3.9 m at Tianjing; on FP the subsided area is relatively smaller than that on the other two plains and is about 1,135 km² with maximum subsidence of 3.7 m at Taiyuan city. In order to protect the civil and industrial facilities, it is necessary to predict the future development of land subsidence based on present state. Many researchers proposed several approaches to predict the land subsidence due to groundwater withdrawal according to different geological conditions and groundwater withdrawal practice. This article classifies these approaches into five categories: (i) statistical methods; (ii) 1D numerical method; (iii) quasi-3D seepage model; (iv) 3D seepage model; (v) fully coupled 3D model. In China, the former four categories are presently employed in the prediction practice and their merits and demerits are discussed. According to the prediction practice, 3D seepage model is the best method presently.     

Estimation of Land Subsidence Based on Groundwater Flow Model

This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well.     

Identification of geological characteristics from construction parameters during shield tunnelling

Acta Geotechnica - This paper proposes a framework to identify geological characteristics (GC) based on borehole data and operational data during shield tunnelling using a fuzzy C-means...     

An efficient optimization method for identifying parameters of soft structured clay by an enhanced genetic algorithm and elastic–viscoplastic model

Soft structured clays usually exhibit complex behaviors, which can lead to difficulties in the determination of parameters and high testing costs. This paper aims to propose an efficient optimization method for identifying the parameters of advanced constitutive model for soft structured clays from only limited conventional triaxial tests. First, a new real-coded genetic algorithm (RCGA) is proposed by combining two new crossover and mutation operators for improving the performance of optimization. A newly developed elastic–viscoplastic model accounting for anisotropy, destructuration and creep features is enhanced with the cross-anisotropy of elasticity and is adopted for test simulations during optimization. Laboratory tests on soft Wenzhou marine clay are selected, with three of them being used as objectives for optimization and others for validation. The optimization process, using the new RCGA with a uniform sampling initialization method, is carried out to obtain the soil parameters. A classic genetic algorithm (NSGA-II)-based optimization is also conducted and compared to the RCGA for estimating the performance of the new RCGA. Finally, the optimal parameters are validated by comparing with other measurements and test simulations on the same clay. All comparisons demonstrate that a reliable solution can be obtained by the new RCGA optimization combined with the appropriate soil model, which is practically useful with a reduction in testing costs.     

Estimation of Land Subsidence Based on Groundwater Flow Model

This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well.     

Analysis of urbanisation-induced land subsidence in Shanghai

Since 1980, land subsidence has accelerated and groundwater levels have decreased in the centre of Shanghai, although the net withdrawn volume of groundwater has not increased. Theoretical analysis of the monitored data shows that the decrease in the groundwater level is the primary reason for the observed land subsidence. Meanwhile, the net withdrawn volume of groundwater in the urban centre of Shanghai has not increased during this period. Many underground structures have been constructed in the multi-aquifer-aquitard system of Shanghai since 1980. This paper discusses the factors related to the development of land subsidence during the process of urbanisation in Shanghai during the past 30?years. These factors include additional load during and after structure construction, the cut-off and/or partial cut-off effect of underground structures in aquifers, the decrease in the groundwater level due to leakage of underground structures and the reduction in recharge of groundwater from the surroundings.     

Evaluation of land subsidence by considering underground structures that penetrate the aquifers of Shanghai, China

Underground structures that penetrate into aquifers can cause groundwater-level drawdown and land subsidence. Numerical analyses, based on a three-dimensional (3-D) groundwater flow model incorporated with a 1-D consolidation model, have been conducted to assess the behaviour of seepage and effect on subsidence by considering underground structures in the multi-aquifer?Caquitard system of Shanghai, China. Two extreme scenarios were examined: (1) distributed underground structures, and (2) concentrated underground structures around the heavily urbanized area. In the first scenario, the aquifer with underground structures was substituted with another material that possessed a lower hydraulic conductivity, established using the effective-medium theory; when the ratio of the volume of the underground structures to that of selected aquifer layers??(1) low-pressure partially-confined aquifer (Aq02), (2) the first confined aquifer (AqI), and (3) the second confined aquifer (AqII)??increases by 10?%, subsidence increases by about 3, 3 and 32?%, respectively. In the second scenario, part of the aquifer material was directly replaced by the structure material (very low hydraulic conductivity). In this situation, when the ratio of the volume of the underground structure to the volume of aquifers Aq02, AqI or AqII increases by 10?%, subsidence increases by about 5, 8 or 20?%, respectively.     

Laboratory Studies on Property Changes in Surrounding Clays Due to Installation of Deep Mixing Columns

This paper has identified six major factors causing property changes in surrounding soils during and after installation of deep mixing columns: soil thixotropy, soil fracturing, cement penetration and diffusion, cementation, consolidation, and heating. Laboratory tests were performed to investigate the effects of soil thixotropy, soil fracturing, and cementation in a soft marine clay, Ariake clay. Laboratory tests were conducted to evaluate property changes in surrounding clays due to installation of deep mixing columns. Test results showed that an influential zone of property changes existed in surrounding clay ranging from the edge of the columns to the distance of about twice the radius of the columns. Within this influential zone, water content decreased as samples neared the columns, while pH values and electric conductivity increased. Test results also showed that undrained shear strengths of the surrounding clays decreased during mixing but regained after a 7-day curing period and continued increasing during 28 days in this study.