三维区域地面沉降数值模拟

由于地裂缝研究及地表变形监测技术（例如GPS,InSAR等）的快速发展,抽取地下水引起的地层水平位移受到关注。传统区域地面沉降模型虽然求解快速但不能模拟水平位移;比奥模型虽然能够模拟土体的三维变形,但模型求解的计算量较大,较少应用于大尺度的区域地面沉降数值模拟。为解决以上问题,推导了解耦三维地面沉降数学模型,模型推导过程显示：比奥模型假设法向总应力和不变,则可简化为解耦三维地面沉降模型;解耦三维地面沉降模型假设土体仅有垂向一维变形,则可简化为传统区域地面沉降模型。同时通过数值试验验证了解耦三维地面沉降模型可以作为比奥模型的替代模型和传统区域地面沉降模型的改进模型,用来模拟抽取地下水引起的三维区域地面沉降。     

Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015)

The North China Plain (NCP) has been suffering from groundwater storage (GWS) depletion and land subsidence for a long period. This paper collects data on GWS changes and land subsidence from in situ groundwater-level measurements, literature, and satellite observations to provide an overview of the evolution of the aquifer system during 1971–2015 with a focus on the sub-regional variations. It is found that the GWS showed a prolonged declining rate of ?17.8?±?0.1 mm/yr during 1971–2015, with a negative correlation to groundwater abstraction before year ~2000 and a positive correlation after ~2000. Statistical correlations between subsidence rate and the GWS anomaly (GWSA), groundwater abstraction, and annual precipitation show that the land subsidence in three sub-regions (Beijing, Tianjin, and Hebei) represents different temporal variations due to varying driver factors. Continuous drought caused intensive GWS depletion (?76.1?±?6.5 mm/yr) and land subsidence in Beijing during 1999–2012. Negative correlations between total groundwater abstraction and land subsidence exhibited after the 1980s indicate that it may be questionable to infer subsidence from regional abstraction data. Instead, the GWSA generally provides a reliable correlation with subsidence. This study highlights the spatio-temporal variabilities of GWS depletion and land subsidence in the NCP under natural and anthropogenic impacts, and the importance of GWS changes for understanding land subsidence development.     

Microwave D-InSAR technique for assessment of land subsidence in Kolkata city,India

An effective microwave Differential Interferometric Synthetic Aperture Radar (D-InSAR) technique was used to rapidly assess the potential land subsidence with high precision by exploiting the phase difference of two temporally separated SAR data in the region of Kolkata city, India. The objective of this study is to assess land subsidence using D-InSAR technique and to delineate the regions of land subsidence caused by over exploitation of groundwater by minimising the errors by applying topographic and atmospheric corrections. The study area forms a part of Indo-Gangetic plain. Three ENVISAT Advanced Synthetic Aperture Radar ?(ASAR) images of the years 2003, 2007 and 2010 were acquired to study the temporal evolution of land subsidence in the study area. The phase changes due to topography in the interferograms were removed by using Shuttle Radar Topography Mission (SRTM) degital elevation model data. Medium Spectral Resolution Imaging Spectrometer (MERIS) data were applied to remove the atmospheric noise in the interferogram. The deformation fringes were observed in the northern and central part of the study area where the land subsidence was 12 and 18 mm during the years of 2003–2007 and 2007–2010. The regional variation in the piezometric head compares well with the fringes of the interferogram. This confirms over extraction of groundwater is the main cause for land subsidence in this region. Hence, it is necessary to reduce groundwater pumping and to augment rainfall recharge in northern part of the study area.     

Sustainable development and utilization of groundwater resources considering land subsidence in Suzhou,China

Suzhou is located at the lower reaches of the Yangtze River in southeastern Jiangsu, China. It is part of the Su-Xi-Chang area including Suzhou, Wuxi and Changzhou. As one of the most developed areas in China, this region has suffered from severe land subsidence caused by extensive groundwater exploitation since 1980s. The land subsidence was controlled by prohibition of groundwater exploration in the past several years. However, the surface water pollution prompted a new task of how to sustainably utilize the groundwater resource, especially to satisfy the emergency demands of water supply. In this paper, we took Suzhou as a representative case to discuss how to develop groundwater resources while controlling the land subsidence. The relationship between the deformation and the groundwater level was analyzed, with focus on the deformation features after the period of groundwater exploitation ban. The results confirmed the conclusion by Shi et al. (2007, 2008a): even in the period of rising groundwater level, same units may manifest different deformation characteristics, such as elasticity, elasto-plasticity, and visco-elasto-plasticity, at different locations of the cone of depression. A land subsidence model that couples a 3-D groundwater model and a 1-D deformation model was developed to simulate the groundwater level and deformation. A high-resolution local grid (child model) for Suzhou was built based on the regional land subsidence model of Su-Xi-Chang area by Wu et al. (2009). The model was used for a number of predictive scenarios up to the year of 2012 to examine how to develop sustainable use of groundwater resources under the conditions of land subsidence control. Our results indicated that about 3.08 × 10⁷ m³/a groundwater could be provided as emergency and standby water source while meeting the land subsidence control target of 10 mm/a.     

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. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

上海市地下水位大幅抬升条件下土层变形特征分析

过去对地下水位持续下降条件的地面沉降研究较多,但对水位大幅持续抬升过程中的地面沉降研究较少。本文根据 上海大量地面沉降、水位观测和钻孔资料,系统分析上海市90年代末以来地下水位大幅抬升条件下各土层的变形特征。自 1998年以来,上海市通过大幅压缩开采量、回灌地下水等措施使第二、三、四和五承压含水层水位分别平均抬升2.1 m, 3.6 m, 12.4 m, 12.7 m。水位的抬升使上海市地面沉降平均速率由1998年的12.2 mm/a减小到2011年的1.83 mm/a,减少85%。 通过对27组分层标数据分析发现：现阶段主要压缩层位在第一、二软土层,年沉降速率为2~4 mm/a;而第二含水层以下土 层已经有少量回弹。在水位持续大幅抬升过程中,本文总结了两种变形特征：1） 变形和水位变化基本同步,残余变形量非 常小,变形可概化为线弹性变形,这种变形主要发生在第一、二、三和五承压含水层、第五和六弱透水层;2） 压缩速率逐 渐减小,无明显持续回弹趋势,有较大残余压缩量且存在变形滞后现象,变形可概化为弹塑性变形,这种变形主要发生在 地第二、三和四弱透水层。第四承压含水层变形较复杂,两种变形特征都有。其中较大残余变形量主要由塑性贮水率比弹 性贮水率大2个数量级引起;变形滞后主要由弱透水层中超孔隙水压力消散较慢引起。本文研究成果对于掌握水位抬升过 程中土层变形方式、发生发展机理、预测未来地面沉降及地下水科学管理和资源评价具有重要意义。     

高层建筑引发地面沉降模拟预测三维流固全耦合模型

为了准确预测由高层建筑引发土体应力场和渗流场变化而导致的地面沉降,以比奥固结理论为基础,结合土体非线性流变理论,将比奥固结理论中的本构关系拓展到黏弹塑性,并考虑了土体孔隙度、渗透系数及变形参数随有效应力的动态变化关系。以河北省沧州市为例,建立了沧州市高层建筑荷载、地下水渗流与土体变形三维流固全耦合数学模型。在对模型进行识别、验证的基础上,模拟预测了沧州市在地下水停采、仅存在高层建筑荷载的影响下,从2010年12月底到2025年12月底逐年的各含水层组地下水流场变化特征和地面沉降发展趋势。结果表明：沧州市由高层建筑荷载引发的最大地面沉降量为40.57 mm,最大地面沉降速率为2.7 mm/a,位于沧州市区。     

Land subsidence and earth fissures due to groundwater withdrawal in the Southern Yangtse Delta,China

Excessive extraction of groundwater has caused severe land subsidence and earth fissures in the Southern Yangtse Delta, China. Based on field data, the temporal and spatial distribution of land subsidence is investigated and the causes for earth fissures are analyzed. The areal distribution of the land subsidence is closely related to the cones of depression in the main exploited aquifers. The compaction of a hydrostratigraphic unit depends on its mechanical behavior, thickness, compressibility, and the piezometric level changing. The primary subsidence layers in Shanghai have been the first aquitard before 1990 and the third confined aquifer since then. But the second aquitard unit was the primary subsidence layer in Changzhou. Earth fissures, trending in several directions, occurred in the Husu tectonic zone. They were mainly caused by differential subsidence and horizontal displacement that resulted from tensile stress and shear stress in units. The majority of fissures in the study area are tensile.     

沧州市地面沉降灾害预测预警

为了准确模拟评价沧州市地下水开采对地面沉降的影响,为沧州市政规划和地下水资源管理提供决策依据,基于比奥固结理论,建立了地下水开采与地面沉降三维全耦合数学模型。在对模型进行识别和校正的基础上,模拟预测了在地下水现状开采情况下,从2010年12月31日至2025年12月31日逐年的地面沉降变化趋势,并根据地面沉降速率对地面沉降进行了地质灾害预警分区。结果表明,到2025年12月31日,沧州市累计最大地面沉降量为466.82 mm,最小地面沉降量为241.54 mm,大部分地区为四级预警区和五级预警区,仅肃宁县为三级预警区。     

地下水位变化模式下含水砂层变形特征及上海地面沉降特征分析

土体变形特征与其经历的应力状态有关。由于抽灌水位置和水量的变化,同一土层中不同时期的地下水位可以呈现不同的变化模式,土层表现出不同的变形特征。论文根据上海1400多个水位孔近40a的水位观测资料和各土层的变形资料,从土层变形角度将地下水位的变化方式划分为5种模式。分析了每种地下水位变化模式下土层的变形特征,并进一步分析了上海地面沉降在时间和空间上的特征。分析结果表明:地下水位的变化模式对上海土层的变形有显著影响。同一土层在不同的水位变化模式下可表现为弹性、弹塑性或粘弹塑性的变形特征;地面沉降与地下水开采量、地下水开采层次与主要沉降层具有密切的关系,开采地下水是上海地面沉降的主要原因;与现阶段含水层的水位变化模式相联系,第四承压含水层是上海最近几年来地面沉降的主要沉降层。     

Threat of land subsidence in and around Kolkata City and East Kolkata Wetlands,West Bengal,India

This paper attempts to estimate the possible rate of land subsidence of Kolkata City including Salt Lake City and the adjoining East Kolkata Wetlands located at the lower part of the deltaic alluvial plain of South Bengal basin. Demand of groundwater for drinking, agricultural and industrial purposes has increased due to rapid urbanization. The subsurface geology consists of Quaternary sediments comprising a succession of clay, silty clay and sand of various grades. Groundwater occurs mostly under confined condition except in those places where the top aquitard has been obliterated due to the scouring action of past channels. Currently, the piezometric head shows a falling trend and it may be accelerated due to further over-withdrawal of groundwater resulting in land subsidence. The estimated mean land subsidence rate is 13.53 mm/year and for 1 m drop in the piezometric head, the mean subsidence is 3.28 cm. The surface expression of the estimated land subsidence is however, cryptic because of a time lag between the settlement of the thick low-permeable aquitard at the top and its surface expression. Therefore, groundwater of the cities and wetland areas should be developed cautiously based on the groundwater potential to minimize the threat of land subsidence.     

Shearing along faults and stratigraphic joints controlled by land subsidence in the Valley of Queretaro,Mexico

Slip of nearly vertical faults or horizontal stratigraphic joints has provoked the shearing of at least 16 well casings in a period of over 10 years in the Valley of Queretaro aquifer, Mexico. Evidence integrated from field observations, remote surface-deformation monitoring, in-situ monitoring, stratigraphic correlation, and numerical modeling indicate that groundwater depletion and land subsidence induce shearing. Two main factors conditioning the stress distribution and the location of sheared well casings have been identified: (1) slip on fault planes, and (2) slip on stratigraphic joints. Additionally, the distribution of piezometric gradients may be a factor that enhances shearing. Slip on faults can be generated either by the compaction of sedimentary units (passive faulting) or by slip of blocks delimited by pre-existing faults (reactivation). Major piezometric-level declines and the distribution of hydraulic gradients can also be associated with slip at stratigraphic joints. Faults and hydraulic contrasts in the heterogeneous rock sequence, along with groundwater extraction, influence the distribution of the gradients and delimit the compartments of groundwater in the aquifer. Analogue modeling allowed assessment of the distribution of stress–strain and displacements associated with the increase of the vertical stress. Fault-bounded aquifers in grabens are common in the central part of Mexico and the results obtained can be applied to other subsiding, structurally controlled aquifer systems elsewhere.     

济宁市地下水与地面沉降三维有限元模拟

在分析济宁市水文地质条件的基础上 ,对引起地面沉降诸因素进行了分析 ,并建立了地面沉降量与地下水位变幅之间的相关关系。集中过量开采地下水是引起济宁市地面沉降的主要原因。在此基础上 ,建立了准三维地下水流模型和一维地面沉降模型。通过水力联系建立地下水与地面沉降耦合数值模型 ,运用有限元法对地下水渗流场和地面沉降量进行模拟 ,并对 2 0 0 0年和 2 0 10年地面沉降进行了预测     

我国地面沉降若干问题研究

摘要：地面沉降主要有（1）主要由过量开采地下水引起的和（2）在新构造活动构造沉降的基础上，由开采地下水、油气、地下热水形成的地面沉降两种类型。不仅不同的土层有不同的变形特点，在不同地点的同一层土也有不同的变形特点。同一地点同一层土在不同时期的变形特点也不尽相同。建立数学模型时要根据土层实际变形特点分别建立相应的模型，才能模拟出真实的地面沉降情况。解决耦合问题和垂向、水平向尺度不同带来的困难也是正确模拟必不可少的。     

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.     

Hydrogeological investigation of groundwater artificial recharge by treated wastewater in semi-arid regions: Korba aquifer (Cap-Bon Tunisia)

Korba aquifer is one of the most typical examples of overexploited coastal aquifer in the Mediterranean countries. In fact, from 1985, a considerable piezometric level drop, water salinization, and seawater intrusion were registered in the aquifer. In December 2008, Tunisian authorities initiated a general plan to groundwater management in order to augment groundwater resources, restore the piezometric levels, and improve water quality. The plan consists of artificial recharge of groundwater used treated wastewater through three infiltration basins. During the first 4 years (from December 2008 to December 2012), 1.41 Mm³ of treated wastewater was injected to the Korba aquifer. This study presents a hydrogeological assessment of groundwater evolution during the recharge processes. In this study, 32 piezometric and chemical surveys of 70 piezometers and observed wells are used to present hydrogeological investigation and water quality evolution of wastewater reuse through artificial recharge in Korba coastal aquifer. The piezometric evolution maps are used to specify the positive effect in groundwater level that exceeding 1.5 m in some regions. The interpretation of salinity evolution maps are used to indicate the improving of groundwater quality.     

The southeastern border of the Upper Rhine Graben: a 3D geological model and its importance for tectonics and groundwater flow

A 3D geological model of the area east of Basel on the southeastern border of the Upper Rhine Graben, consisting of 47 faults and six stratigraphic horizons relevant for groundwater flow, was developed using borehole data, geological maps, geological cross sections, and outcrop data. This model provides new insight into the discussions about the kinematics of the area between the southeastern border of the Upper Rhine Graben and the Tabular Jura east of Basel. A 3D analysis showed that both thin-skinned and thick-skinned tectonic elements occur in the modeled area and that the Anticline and a series of narrow graben structures developed simultaneously during an extensional stress-field varying from E–W to SSE–NNW, which lasted from the Middle Eocene to Late Oligocene. In a new approach the faults and horizons of the 3D geological model were transferred into discrete elements with distributed hydrogeological properties in order to simulate the 3D groundwater flow regime within the modeled aquifers. A three-layer approach with a horizontal regularly spaced grid combined with an irregular property distribution of transmissivity in depth permitted the piezometric head of the steady-state model to be automatically calibrated to corresponding measurements using more than 200 piezometers. Groundwater modeling results demonstrated that large-scale industrial pumping affected the groundwater flow field in the Upper Muschelkalk aquifer at distances of up to 2 km to the south. The results of this research will act as the basis for further model developments, including salt dissolution and solute transport in the area, and may ultimately help to provide predictions for widespread land subsidence risks.     

Comparison of water-level,extensometric, DInSAR and simulation data for quantification of subsidence in Murcia City (SE Spain)

Subsidence due to groundwater overexploitation has been recognized in the metropolitan area of Murcia (25 km²) in south-eastern Spain since the early 1990s. Previous published works have focused their attention on land subsidence that occurred during the drought period between 1995 and 2008. This work first analyzes the groundwater recovery that has occurred since 2008 and then determines the kind of associated ground deformation detected by the new extensometric data. Subsequently, subsidence time series are computed on 24 geotechnical boreholes scattered throughout the study area by means of a hydro-mechanical finite element code and a linear-elastic constitutive law. A spatio-temporal interpolation of the numerically modeled surface displacements is performed over the whole domain and compared with extensometers and DInSAR-derived displacement maps in two different periods: the drought period from 2004 to 2008, and the recovery period from 2008 to 2012. In spite of the limited information on the geomechanical parameters characterizing the modelled geological formations, the proposed approach is able to discriminate areas where the soils have an elastic behavior (small differences in the comparisons) or an elasto-plastic behavior (large differences in the comparisons). This zonation enhances the understanding of the subsidence phenomenon in Murcia City and could prevent, from a quantitatively point of view, future severe subsidence due to aquifer overexploitation.     

弹塑性变形条件下抽水引起的地面沉降三维数值模拟

基于部分耦合原理,采用TOUGH2和FLAC3D建立抽水引起的三维地面沉降弹塑性模型,模型中综合考虑土体的弹塑性变形特征、渗流-应力的双向耦合作用以及参数的非线性,探讨了持续抽水和脉冲抽水两种抽水过程中地面沉降发展演化过程。研究结果表明:（1）集中抽水停止后地面沉降会发生回弹,抽水中心沉降量不断减小。由于水平方向存在水力梯度,地下水继续向地下水位漏斗中心渗流从而导致沉降漏斗的范围仍继续扩大;（2）脉冲抽水导致土体的孔隙水压力、渗透系数以及沉降量均呈周期性波动变化,地面沉降会局部回弹,但总体仍随着抽水的持续,沉降量不断增加;（3）在抽水量相同前提下,对比持续抽水与脉冲抽水两种方式引发的塑形沉降量可知,抽水速率小、脉冲式多次开采导致的塑性沉降量较小,持续抽水的抽水速率越小、脉冲抽水间隔越短越有利于控制地面沉降。研究成果为地面沉降数值模拟提供了一种新方法,其中算例研究能为抽水条件下地面沉降的控制提供参考。