土体自重固结压缩对地面沉降影响研究

以南通市为例,建立了土体自重固结压缩、地下水开采与地面沉降三维全耦合模型.在对模型进行识别、验证的基础上,模拟预测了南通市在地下水停止开采,仅在土体自重固结压缩影响下,自2010年12月31日-2025年12月31日各含水层组地下水流场变化和地面沉降的发展趋势.结果表明,在此期间,南通市由土体自重固结压缩引起的最大地面沉降量为2.42 mm,最大地面沉降速率为0.16 mm/a,土体白重固结压缩对南通市地面沉降的影响极为有限.     

Mechanics of land subsidence due to groundwater pumping

This paper presents the formulation of the basic mechanics governing the changes in stress states from groundwater pumping and comparisons among predicted land subsidence from this mechanics with existing analyses and field data. Land subsidence is a growing, global problem caused by petroleum and groundwater withdrawal, mining operations, natural settlement, hydro‐compaction, settlement of collapsible soils, settlement of organic soils and sinkholes. This paper is concerned with the land subsidence due to groundwater level decline by groundwater pumping. It is shown that the stress state consists of asymmetric stresses that are best simulated by a Cosserat rather than a Cauchy continuum. Land subsidence from groundwater level decline consists of vertical compression (consolidation), shear displacement and macro‐rotation. The latter occurs when conditions are favorable (e.g. at a vertical interface) for the micro‐rotation imposed by asymmetric stresses to become macro‐rotation. When the length of the cone of depression is beyond √2 times the thickness of the aquifer, simple shear on vertical planes with rotation is the predominant deformation mode. Otherwise, simple shear on horizontal planes is present. The predicted subsidence using the mechanics developed in this paper compares well with data from satellite‐borne interferometric synthetic aperture radar. Copyright © 2009 John Wiley & Sons, Ltd.     

地下水渗流与地面沉降耦合模拟

为了准确模拟由地下水开采导致渗流场和应力场发生变化而引起的地面沉降问题,根据Terzaghi有效应力原理,建立了地下水三维渗流与一维垂向固结的地下水渗流与地面沉降耦合数值模拟模型和以比奥固结理论为基础,并结合土体非线性流变理论,将土体本构关系推广到粘弹塑性,同时考虑土体力学参数及水力参数的动态变化关系的地下水渗流与地面沉降三维全耦合数值模拟模型.通过对比分析,结果表明:基于Terzaghi有效应力原理建立的地下水三维渗流与一维垂向固结地下水渗流与地面沉降耦合数值模拟模型模拟所得地面沉降与地下水位呈现出同步变化的趋势,并且当地下水位逐步回升至初始水位时,地面沉降也逐步回升到初始的零沉降状态.而以比奥固结为基础建立的地下水渗流与地面沉降三维全耦合数值模拟模型模拟所得的地面沉降变化趋势滞后于地下水位的变化趋势,并且当地下水位逐步回升至初始水位时,地面沉降虽也逐步得到回升,但回不到初始的零沉降状态,存在一个永久的残余沉降量.在土体参数变化方面,土体的孔隙度、渗透系数及泊松比均呈现先减小后增大的变化趋势,而弹性模量则呈现先增大后减小的变化趋势,与地面沉降的变化相对应.      

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

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

Excessive groundwater withdrawal and resultant land subsidence in the Su-Xi-Chang area,China

Excessive groundwater withdrawal has caused severe land subsidence in the Su-Xi-Chang (SXC) area, China. The restriction and prohibition on groundwater pumping have been carried out since the late 1990s. Based on the latest updated field data, the changing pattern of groundwater level and the distribution of land subsidence are analyzed. The distribution of land subsidence in SXC is closely related to that of the cone of depression in the second confined aquifer in time and space. But land subsidence is not in synchronization with the changing groundwater level. Both aquitards and aquifers compacted continuously in the early period of groundwater level rising and behaved as creep materials. A series of laboratory tests were conducted on aquifer sands, which indicated that the creep deformation under virgin compression is much greater than that under recompression and unloading, and that the creep of sands decreases rapidly with the cycles of repeating load. The test results reveal the mechanism of sand creep under the condition of long-term groundwater pumpage. As a consequence of the restriction and prohibition on groundwater pumping, groundwater level has obviously recovered in the vast majority of the SXC area, and land subsidence has slowed down and even a little rebound has occurred in some sites in Suzhou and Changzhou. If the pumpage is strictly limited continuously, the groundwater level will not decline below the historical lowest value but fluctuate within a certain range. In such a case, land subsidence in SXC will no longer develop obviously.     

Estimation of land subsidence caused by loss of smectite-interlayer water in shallow aquifer systems

Traditionally, land subsidence that results from groundwater over-pumping has often been described by the theory of consolidation. The mechanism of land subsidence due to the dehydration of clay minerals is not well addressed. A model of the “hydration state of smectite”, and a “solid solution model of smectite dehydration”, incorporating a thermodynamic solid solution model and laboratory results concerning clay-water systems of swelling pressure, hydration state and basal spacing in smectite interlayer, are employed to examine the effect of the release of water from the smectite interlayer on land subsidence in the coastal area of the Chou-Shui River alluvial fan and the Yun Lin offshore industrial infrastructure complex in Taiwan. The results indicate that 9.56–22.80% of the total cumulative land subsidence to a depth of 300 m is consistent with smectite dehydration following the over-pumping of groundwater. This dehydration-related land subsidence occurred to a depth of 0–60 m, with subsidence due to smectite dehydration accounting for 6.20–13.32% of the primary consolidation. Additionally, the total amount of subsidence resulting from both smectite dehydration and primary consolidation is consistent with the subsidence observed in the field. This study reveals that smectite dehydration appears to be important in assessing and predicting land subsidence in shallow aquifer systems.     

Investigation of land subsidence due to groundwater withdraw in Rafsanjan plain using GIS software

Nowadays, the purpose of predicting land subsidence is to manage the optimum usage of groundwater, which is considered according to irregular use of groundwater. Digging deep and semi-deep wells and continuous drought, mainly in wasteland and semi-wasteland zone in recent years causes the land subsidence in Rafsanjan plain. The Rafsanjan basin is located in the nearly central part of Iran in the Kerman province, with a general elevation between 1,400–1,500 m above sea level. In this research, first, the deep and semi-deep wells were investigated and groundwater table data were colleted. Second, these informations were analyzed and corrected. These data were used to create great bank of information data, to manage and program the geographic information system (GIS) software. Then by investigation of an existing land subsidence data, which were collected by GPS in August 1998 and April 1999, by the GIS software, the results show that discharging of groundwater is the main factor of the land subsidence in Rafsanjan zone. Therefore, the critical land subsidence zone of the Rafsanjan plain was determined, and precaution and recommendations are presented.     

不同方法求解疏排水引起的地面沉降对比研究

地下水开发利用和疏排水过程引起的地面沉降是一个渗流场与应力场耦合的过程。基于算例模型,分别采用目前常用的基于渗流-沉降分步计算的土力学经验公式、渗流-沉降部分耦合的GMS软件中SUB模型以及目前少用的基于渗流-沉降全耦合的COMSOL Multiphysics模型对疏排水引起的地面沉降量进行了计算。对比不同方法计算结果表明,COMSOL Multiphysics计算疏排水引起的地面沉降可行且计算结果较经验公式法、GMS中SUB模型的结果更为合理、更符合实际沉降特征。在此基础上,采用COMSOL Multiphysics和GMS中SUB模型对比研究了不同渗透系数对地面沉降计算的影响,结果进一步验证了COMSOL Multiphysics求解疏排水引起地面沉降的可靠性。     

Groundwater resources management under environmental constraints in Shiroishi of Saga plain, Japan

A sinking of the land surface due to the pumping of groundwater has long been recognized as an environmental issue in the Shiroishi plain of Saga, Japan. Land subsidence can have several negative economic and social implications such as changes in groundwater and surface water flow patterns, restrictions on pumping in land subsidence prone areas, localized flooding, failure of well casings as well as shearing of structures. To minimize such an environmental effect, groundwater management should be considered in this area. In this study, a new integrated numerical model that integrates a three-dimensional numerical groundwater flow model coupled with a one-dimensional soil consolidation model and a groundwater optimization model was developed to simulate groundwater movement, to predict ground settlement and to search for optimal safe yield of groundwater without violating physical, environmental and socio-economic constraints. It is found that groundwater levels in the aquifers greatly vary from season to season in response to the varying climatic and pumping conditions. Consequently, land subsidence has occurred rapidly throughout the area with the Shiroishi plain being the most prone. The predicted optimal safe yield of the pumping amount is about 5 million m³. The study also suggests that pumping with this optimal amount will minimize the rate of land subsidence over the entire area. An erratum to this article can be found at     

The Characteristics and Causes of Land Subsidence in Tanggu Based on the GPS Survey System and Numerical Simulation

Land subsidence is a severe hazard threatening Tanggu, a flat lowland area, and evidences of land subsidence can be seen throughout the city. A new reasonable GPS network was set up in this area from 2008 to 2010. The monitoring data show that land subsidence was serious and two main subsidence cones were obviously formed in the region. One emerged at Hujiayuan, with the maximum rate reaching 60 mm/a, and the influence region enlarged prominently from 2005 to 2010. The other one occurred at Kaifaqu, which became obvious only after 2005, and it showed a decreasing tendency with time. To analyze the causes of ground settlement, a correlation between groundwater withdrawal and land subsidence was firstly made. The results confirmed that over-exploitation of groundwater was the major cause for the severe settlement in Hujiayuan. Meanwhile, the subsidence of Kaifaqu was also related to groundwater withdrawal before 2005. However, the relationship became unconspicuous after 2005. To find the cause of this abnormity, a three-dimensional finite element numerical model, coupled with groundwater flow and subsidence, was built. The simulation results indicate that the subsidence induced by high-rise buildings is serious, but the affected range is limited and it also shows a decreasing trend with time, corresponding to the subsidence characteristics at Kaifaqu. Therefore, more attention should be paid to this hazard induced by engineering construction besides groundwater withdrawal, as more high-rise buildings are under construction in Tanggu.     

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.     

天津市地下水流-地面沉降耦合模型

天津市平原区地面沉降主要由地下水大量开采引起,影响范围广、危害大,已成为天津市主要的环境地质问题。分析了研究区的水文地质条件,结合地下水开发利用状况,将研究区概化为6个含水层组,地下水流考虑三维非稳定流,地面沉降选用一维固结压缩模型,运用地下水流模型Modflow 2005和地面沉降模拟模块 Sub,建立了天津市平原区地下水流-地面沉降数值耦合模型,模型面积为1.1×10⁴ km²,利用1998-2008年地下水位等值线、过程线、地面沉降过程线等资料对模型进行了识别。模拟期的地下水均衡分析表明,在多年开采条件下,越流补给、压缩释水、侧向边界流入分别占深层含水层补给量的41.84％、32.15％和24.17％。将调试后的模型应用于南水北调实施后地下水控采条件下的地面沉降趋势预测,显示出停采或减少地下水的开采,有利于减缓地面沉降下降速度,且表现出开采层位越往下,地面沉降恢复难度越大的变化趋势。     

基于地下水渗流方程的三维地面沉降模型

针对地下水抽取引起的大面积地面沉降的问题,提出了一种计算方法。该方法基于饱和与不饱和岩土介质中地下水渗流理论,计算出三维状态下的大面积地面沉降。提出的方法被结合到地下水渗流的三维有限元法分析中,计算时考虑了地层的可压缩性,并被用来分析软土的固结,将该法与太沙基理论及固结试验的结果进行比较,结果表明其误差小于2 %。该法还被用来分析单井抽取承压地下水引起的周围地基的沉降及承压水层由于挡水板的阻断作用而引起的下游侧的地基的沉降。     

基于随机介质理论的抽水地表变形时空耦合预测模型

抽取地下水引起的地表变形对地面建(构)筑物的正常使用和结构安全构成了严重威胁。深入研究抽水地表变形预测理论对于沉降灾害防治具有重要意义。由于土体模型选择和土工参数测定上的困难,基于比奥固结理论的数值求解方法在计算抽水地表变形上尚未取得理想效果。文章将随机介质理论与一维固结理论结合,建立了新的抽水地表变形时空预测模型。首先,利用一维固结微分方程,建立反映地面沉降时间效应的半经验计算模型;其次,在分析抽水地面沉降空间分布规律的基础上,利用随机介质理论研究抽水地面沉降空间分布特征;再次,综合考虑抽水地面沉降的时间效应和空间分布形态,建立抽水地面沉降的时空耦合预测模型及抽水地表倾斜、水平移动、水平变形、曲率的时空计算模型。利用这些模型计算地表变形共需5个计算参数,介绍了参数求解方法。最后,利用上述时空耦合计算模型预测某地单井抽水引起地表变形的时空规律。研究表明,所建立的抽水地表变形预测模型能准确地反映抽水地表变形的时空规律,能方便、快捷地预测地下水开采引起的地表变形。     

地层形变与地下水位分层监测数据的交叉小波分析

地面沉降通常由于地下水的超采而引发,其发生发展相对于地下水位的变化具有一定滞后性。如何获取准确的地面沉降滞后时间一直是地面沉降研究的重要课题。基于北京平各庄地面沉降监测站2008—2018年地面沉降和地下水位长时间序列的分层监测数据,采用Mann-Kendall趋势检验、连续小波变换、交叉小波变换等方法,分析了不同层位地层形变对地下水位动态的滞后特征。结果表明:中–深层承压水具有1 a左右的主震荡周期,潜水和浅部承压水在大部分时域无显著周期;深部严重沉降层的形变量具有1 a左右的主震荡周期,地下水位与形变量共振周期显著,地层由浅到深形变时滞分别为(16.58±8.91)、(7.16±7.09)和(9.66±6.62) d;浅部弱沉降层中,埋深在32~63 m地层形变量具有1 a左右的主震荡周期,与中层承压水存在显著共振周期,形变时滞为(32.02±9.67) d,其他地层形变量与地下水位无显著周期及相关性。研究成果为构建地面沉降精细化模型、提高地面沉降预测精度以及研究更有效的地面沉降防控措施提供了新的技术思路。      

Finite-difference model of land subsidence caused by cluster loads in Zhengzhou,China

Groundwater exploitation has been regarded as the main reason for land subsidence in China and thus receives considerable attention from the government and the academic community.Recently,building loads have been identified as another important factor of land subsidence,but researches in this sector have lagged.The effect of a single building load on land subsidence was neglected in many cases owing to the narrow scope and the limited depth of the additional stress in stratum.However,due to the superposition of stresses between buildings,the additional stress of cluster loads is greater than that of a single building load under the same condition,so that the land subsidence caused by cluster loads cannot be neglected.Taking Shamen village in the north of Zhengzhou,China,as an example,a finite-difference model based on the Biot consolidation theory to calculate the land subsidence caused by cluster loads was established in this paper.Cluster loads present the characteristics of large-area loads,and the land subsidence caused by cluster loads can have multiple primary consolidation processes due to the stress superposition of different buildings was shown by the simulation results.Pore water migration distances are longer when the cluster loads with high plot ratio are imposed,so that consolidation takes longer time.The higher the plot ratio is,the deeper the effective deformation is,and thus the greater the land subsidence is.A higher plot ratio also increases the contribution that the deeper stratigraphic layers make to land subsidence.Contrary to the calculated results of land subsidence caused by cluster loads and groundwater recession,the percentage of settlement caused by cluster loads in the total settlement was 49.43%and 55.06%at two simulated monitoring points,respectively.These data suggest that the cluster loads can be one of the main causes of land subsidence.     

基于DFOS的排灌水条件下土体变形响应模型试验研究

黏土层和砂土层交替变化的多层土体在强烈开采地下水作用下极易产生压密固结而引发地面沉降灾害。本文针对含水层释水引起地面沉降问题,研制了地面沉降试验装置,进行了排灌水条件下含水层系统的沉降及回弹试验。采用分布式光纤感测技术对土体内部应变分布及含水率变化进行耦合监测,并分析了各分层对水位变化的响应特征。结果表明:黏土层和砂土层均表现出了排水压缩和灌水回弹特点,黏土层变形较砂土层明显。各层变形与含水率变化具有良好的对应关系,表现为砂土层变形和含水率变化基本同步,而黏土层变形略微滞后于含水率变化。黏性土压缩曲线具有明显的分段特征,排水时当含水率低于液限后迅速减小,黏土层压缩速率明显加快;回灌时当含水率高于液限后,回弹速率明显加快。试验结果对研究地面沉降机理、评价地面沉降潜力及地下水利用具有重要意义。     

深基坑降水三维变参数非稳定渗流与地面沉降耦合模型

为了精确模拟预测松散沉积层中深基坑降水引起的地下水渗流场和地面沉降的变化特征,考虑土体孔隙度、渗透系数、储水率随地下水位下降发生的动态变化,建立了深基坑降水三维变参数非稳定渗流与太沙基一维固结理论的地面沉降耦合模型,并采用有限元数值分析方法对模型进行求解。以南京地铁三号线浦珠路站深基坑降水为例进行模拟计算。结果表明：采用15口坑内抽水井,抽水井过滤器埋深为22.0～37.0 m,基坑围护连续墙底部埋深至41.5 m为最优降水方案;不仅使基坑内地下水位满足开挖要求,又使基坑外地面沉降在控制范围内。经验证,所建立的模型合理,计算结果可靠,研究理论用于模拟预测此类地区深基坑降水引起的地下水流场变化具有较高的可信度。     

Characterization of land subsidence induced by groundwater withdrawals 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 tools of field monitoring were integrated to study the characteristics and mechanisms of land subsidence in this region. The occurrence and the development of the land subsidence in this region are strongly related to the groundwater pumping both in time and space. The main consolidation layers are the soft mud layers; however, the compressibility of the confined sandy layers should not be ignored. The second and third confined aquifers contributed more than 30% of total subsidence. Meanwhile, irrecoverable deformations were also observed in the sandy layers. Different sandy layers deform diversely under different stress conditions. Some have the elastic feature. But the soil strata, including both sandy layers and clayey layers, located in the center of the groundwater level depression cone exhibited obvious viscous mechanical behavior which caused the common lag phenomenon. The sand composition (mingled with small clay particles or interbeds) and sand rheology are the two main reasons for the lag phenomena in sandy layers. A series of laboratory tests for modeling the effective stress changes due to groundwater withdrawals, were conducted to investigate the mechanism of the lag phenomenon. Based on the test results, the relationship of stress–strain–time for saturated sands is obtained; and it could be expressed as power functions. The results also showed that the compression of the sandy layers was time dependent, and its deformation could be remarkable. When establishing land subsidence model, the deformation for the similar soil formation could be elastic, visco-elastic and even visco-elastic–plastic, because of the different groundwater level fluctuation experienced.     

Simulation of land subsidence using finite element method: Rafsanjan plain case study

Land subsidence caused by the excessive use of groundwater resources has caused serious damage to Rafsanjan area. In this study, using finite element method, a 2D plain strain simulation of land subsidence has been conducted. A linear elastic constitutive law has been used for the simulation of the soil material skeleton. Actual water level during the analysis period has been modeled via specifying nodal water pressure at piezometer wells situation. The solution procedure consists of two parts. First, an initial static analysis is carried out in order to find initial steady-state solution for the pore pressure and stress distribution. Then, the above solution is used as initial condition for dynamic computation of consolidation equations during pumpage period. Tectonic effect has been considered as a rigid body motion. Numerical results showed that if the rate of pumpage remains constant in the future, settlement due to water withdrawal near Rafsanjan city will reach up to 110 cm by the year 2022.