黄河三角洲典型地段地面沉降机理研究及预测预警

随着社会经济的快速发展,黄河三角洲地区普遍发生了地面沉降,已经成为了影响区域社会经济可持续发展的重要因素之一。本文以黄河三角洲地面沉降典型地段—滨州地区为研究对象,以该区2005~2019年年间共计6期地面沉降监测数据为基础,对区内多年来地面沉降发展演变规律进行分析,对地面沉降成因机理进行研究,并根据滨州地下水系统的水文地质条件以及地面沉降机理,概化出滨州地下水系统水位地质概念模型及地面沉降模型。在此模型基础上推演建立滨州地下水开采与地面沉降三维数值模型。利用该模型对地面沉降与地下水位监测数据进行内插外推计算,综合计算结果分析和研究,对地面沉降和地下水动态进行预测预警,从而为黄河三角洲地区地面沉降防治提供依据。     

北京地区基于不同地面沉降阈值的地下水位控制分析

城市地面沉降是对城市规划建设、经济发展和人民生活构成威胁的地质灾害。为了探究地下水位与地面沉降的关系,本文对北京顺义地区天竺地面沉降监测站多年分层地面沉降及对应含水层组地下水位监测数据进行统计分析,建立了该地区基于累计沉降量与含水层组水位标高、水位变幅及水位波动的多元回归模型,并对所建立的回归方程进行了检验,研究分层地面沉降与地下水位变化的定量关系,并结合《全国地面沉降防治规划(2011—2020)》中北京市的控沉目标,设置该地区不同地面沉降速率控制阈值,计算得到各层位达到控制阈值时所对应的地下水位,为下一步合理调整地下水开采层位,开展地面沉降防控工作提供科学依据。     

以控制地面沉降为约束的地下水资源评价--以上海地区为例

以上海市区第V承压含水层(埋深约350m)以上的整个第四系松散沉积层为研究体,将上海地下水系统作为一个完整且不可分割的系统来研究.由于开采地下水资源引起地下水系统内部物质能量变化,因而产生一系列的环境地质问题(如地面沉降、地下水位下降、资源开采条件恶化等).根据研究区地下水系统特征和已有的研究程度,在技术方法和途径上,着重研究上海地区第四系含水层结构、水文地质工程地质条件,建立水文地质工程地质概念模型;在此基础上建立以地下水资源开采量最大、水位降幅及地面沉降较小、社会经济效益较好为目标的地下水拟三维渗流一维地面沉降的评价与管理模型,并为我国大中城市(特别是在诸如开采地下水引起的地面沉降问题的地区)的地下水资源评价与管理提供可借鉴的经验和方法.     

广州金沙洲岩溶地面塌陷灾害预警预报研究

受武广客运专线金沙洲隧道施工抽排大量地下水的影响, 2007—2012年广州金沙洲发生了大规模岩溶地面塌陷、地面沉降地质灾害,造成了重大经济损失。为防治地质灾害,采用综合地质调查、地下水位监测、地面沉降监测、物探、钻探和试验等手段勘查和研究,建立了完善的地下水位和地面沉降监测网络。在此基础上,根据地质环境条件、岩溶发育程度和地质灾害分布特征进行岩溶塌陷易发性分区,选择地下水位变化量和地面沉降量这两个与岩溶塌陷最直接、最敏感和最重要的参数,从预警预报参数选择、时间尺度选取、判据分析计算、预警预报方法、模型系统的建立等方面进行研究,选取临界日综合地下水位变化量和临界日综合地面沉降量作为预警预报判据,结合三维地质结构模型建立岩溶塌陷预警预报模型系统。成功地预警预报了两起岩溶地面塌陷,取得较好的预报效果,说明该岩溶地面塌陷预警预报模型系统实用、有效,为金沙洲防灾减灾提供地质科学依据,为当地的社会稳定和经济可持续发展提供技术支撑。     

苏锡常地区地面沉降地质结构三维可视化模型虚拟现实系统研究

从苏锡常地区地面沉降研究实际需要出发,介绍了利用虚拟现实技术构建地面沉降地质结构三维可视化模型虚拟现实系统的设计思路、开发模式、实现方法和功能特点。该系统真实表达了研究区三维地层结构特征,逼真再现了地面沉降模型所预测的地下水流场、地面沉降发展变化过程,并虚拟表现出地面沉降可能造成的后果,为苏锡常地区地面沉降机理及预测研究提供了一个全新的三维可视化平台。     

天津市地面沉降数值模拟研究

地面沉降是天津市的主要环境地质问题。造成天津市地面沉降的原因很多,但主要原因是大量开采地下水。本文重点研究开采地下水条件下的地面沉降数学模型及模拟,为控制地面沉降提供科学依据。通过对第四系的水文地质条件分析,建立第四系的水文地质概念模型和地层压缩概念模型,并建立相应的准三维水流数学模型和一维地层压缩数学模型,两者通过含水层水头和弱透水层粘性土中孔隙水头的内在联系耦合在一起。利用以往系列水位数据和沉降数据对模型进行标定,模型基本上反映了天津市地面沉降的实际情况。     

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

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

天津市地面沉降及地下水位监测自动化系统的设计与应用

地面沉降是天津市当前最为主要的环境地质问题之一,影响着各类城市基础设施的建设 和维护,不利于城市的快速发展.天津地面沉降研究结果业已表明地下水超采是导致沉降的主要原因,大范围、长时间监测地面沉降发展过程及地下水开采状况是进一步开展地面沉降机理分析、预测地面沉降趋势并采取合理防治措施地基础.本文将详细介绍天津市一项地面沉降和地下水位监测自动化系统的设计和应用情况,为我国地面沉降监测与防治工作摸索新方法.     

广州金沙洲岩溶区地下水位变化与地面塌陷及地面沉降关系探讨

金沙洲可溶性灰岩分布面积广,岩溶洞隙发育,洞隙及地下水的连通性强,上覆第四系松散土体中软土广泛分布,客观存在岩溶地面塌陷及地面沉降的地质环境条件。2007年4月起,受某高铁隧道施工抽排地下水的影响,金沙洲地下水出现异常波动,引发了地面塌陷及地面沉降。文章根据监测数据,经对比分析结果表明,区内岩溶地面塌陷及地面沉降受控于地下水位的变化,地下水位波动至基岩面附近时,是地面塌陷较活跃的时期,地面沉降与地下水位变化呈正向相关。文章进一步对地面沉降与地下水位变化关系的机理进行了探讨,认为目前地下水位尚未恢复正常的区域仍存在地面塌陷及地面沉降的隐患。     

台湾浊水溪冲积扇地下水位变化及地面沉降对高铁之影响

浊水溪冲积扇是台湾水资源最为丰富但也是地面沉降最严重地区。近年来该区域的地面沉降因有可能威胁高铁行车安全而备受关注。本文整合历年累积的地下水位及地面沉降等相关监测数据,验证了地下水位变化与含水层补给之机制,探讨了地下水周期性波降条件下土层压缩特性以及高铁路堤与桩基础工程结构的沉降行为。认为对于设置桩基础的线型高架结构而言,区域性地下水位波降不致增加桥墩间的差异沉降,但桩基础若承受邻近局部的额外载重,则可能伴随地下水位波降产生持续性的差异沉降,其长期效应将对线型交通结构物的平整度及安全性造成负面影响。     

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

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

地面沉降对含水层参数及给水能力的影响研究

地面沉降主要由深层地下水开采造成,含水层的压缩释水是深层地下水开采量主要构成来源,同时随其压缩变形孔隙比减小造成储水系数、渗透系数的减小,对弱透水层非弹性释水量、越流补给量造成影响。地面沉降过程中的水文地质参数非线性变化及对承压含水层系统的反馈作用也成为水文地质领域的前沿问题。为对地下水开采量、沉降量、地质参数变化以及给水能力变化之间的关系做一个较为定量定性的探究,以含水层压缩过程中的物理机制为依据,并基于沧州地区深层承压含水层的地下水位变化过程和水文地质参数,采用数学手段构建储水系数和渗透系数变化的一维非线性沉降模型。模拟结果显示随承压水头下降,储水系数最高可减小77%,含水层系统的给水能力和储水能力会随压缩变形减小50%甚至更多。研究成果为深入认识水文地质参数与应力变化相关关系、科学评价承压含水层地下水储水调节能力有重要借鉴意义。     

基于MODFLOW参数不确定性的地下水水流数值模拟方法

考虑到模型不确定性引起的地下水数值模拟不确定性对模拟过程的影响,在简要介绍含水层水文地质参数变异性研究进展和地质统计学的基础上,基于常用的确定性地下水流数值模拟软件MODFLOW开发了MODFLOW-Gslib软件,相较于传统的数值模拟方法,将地质统计学与数值模拟结合的方法能够模拟非均质含水层中的参数变异性问题。将MODFLOW-Gslib软件运用于模拟实例中,选择常见的不确定性因素进行模拟,并对其模拟产生的数据进行统计分析,结果表明,软件转化后的参数符合水文地质参数不确定性的相关特征;与原模拟结果进行对比,该软件能够更加真实地刻画含水层参数变异性特征。     

Modelling the cutoff behavior of underground structure in multi-aquifer-aquitard groundwater system

The quaternary deposit of Shanghai is composed of an alternated multi-aquifer-aquitard system (MAAS) consisting of a sequence of aquitards laid over aquifers one by one. In the MAAS, any drawdown of groundwater head in an aquifer may cause consolidation of the overburden aquitard. When underground structures penetrate those aquifers, groundwater seepage path changes and drawdown occurs at the side characterized by the lower hydraulic potential along the flow direction (hereafter refers as to the lower side). This drawdown may cause additional subsidence at the lower side and unbalanced load between the two sides of the underground structure. In order to evaluate the cutoff effect of an underground structure on groundwater seepage in a MAAS representative of the underground of the city of Shanghai, a numerical analysis based on a groundwater flow model has been carried out. The simulated results have shown that underground structures which cut off groundwater flow locally change both magnitude and direction of the flow velocity field. The induced changes in the groundwater field are highly sensitive to the penetration depth and width of the underground structure. Design recommendations for underground structures in aquifers belonging to a MAAS are also presented, which has not yet been considered in the engineering practice of Shanghai.     

Hydrogeological model of the Baltic Artesian Basin

The Baltic Artesian Basin (BAB) is a complex multi-layered hydrogeological system in the south-eastern Baltic covering about 480,000 km². The aim of this study is to develop a closed hydrogeological mathematical model for the BAB. Heterogeneous geological data from different sources were used to build the geometry of the model, i.e. geological maps and stratigraphic information from around 20,000 boreholes. The finite element method was used for the calculation of the steady-state three-dimensional (3D) flow of unconfined groundwater. The 24-layer model was divided into about 1,000,000 finite elements. A simple recharge model was applied to describe the rate of infiltration, and the discharge was set at the water-supply wells. Variable hydraulic conductivities were used for the upper (Quaternary) deposits, while constant hydraulic conductivity values were assumed for the deeper layers. The model was calibrated on the statistically weighted borehole water-level measurements, applying L-BFGS-B (automatic parameter optimization method) for the hydraulic conductivities of each layer. The principal flows inside the BAB and the integral flow parameters were analyzed. The modeling results suggest that deeper aquifers are characterized by strong southeast–northwest groundwater flow, which is altered by the local topography in the upper, active water-exchange aquifers.     

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.     

Numerical simulation of groundwater flow and contamination transport in Nahavand plain aquifer,west of Iran

Numerical simulation of groundwater flow used for the estimation of hydraulic and hydrologic parameters which is an important tool for the management of aquifers. This study presents the results of a mathematical model developed for the simulation of groundwater flow in Nahavand plain aquifer in the southwest Hamadan province. For this purpose Groundwater Modeling Software (GMS) was used which supports the MODFLOW-2000 code. After gathering required data such as the hydrological, hydrogeological and topography maps, a 3D hydrogeological model of plain was constructed with borehole and surface elevation data. Then MODFLOW was used for simulation of flow. After initial simulation of the flow, the model was calibrated in steady state with trial-and-error and parameter estimation methods the observed head of groundwater table monitoring data of 1997. Results of calibration show that error between observed head and computed head is in allowable range. Also results of computed head with model show that groundwater flow is in the direction of the dominate slope (southeast to northwest). Finally MODPATH code which simulates advective transport of particles was used for estimation of flow path and source of contaminants.     

Coastal aquifer assessment based on geological and geophysical survey,northwestern Crete,Greece

Groundwater preservation comprises a major problem in water policy. The comprehension of the groundwater/hydraulic systems can provide the means to approach this problem. Generally, drilling is expensive and time-consuming. On the other hand, new techniques have been applied during the last few decades that provide useful information on the depth and quality of aquifers. Among them, transient electromagnetic method (TEM) is an appealing method that provides fast results with minimum field crew and solves several hydrogeological problems. Many portable systems for single-site measurements are commercially available. The TEM-Fast 48HPC was used for acquiring 106 soundings in the northwestern Crete in Greece for defining the hydrogeological characteristics of the study area, since there were no available data from boreholes. Detailed geological, hydrolithological and tectonic survey was applied prior to the geophysical measurements. All the data were integrated to produce a secure and reliable hydrogeological model for the study area prior to any future hydrowell. Specifically, geometrical and hydraulic data of the study area groundwater were acquired. Two unconnected aquifers were detected and their possible contamination due to saltwater intrusion was analyzed and eliminated. Moreover, a location for borehole construction and groundwater pumping based on the potential of the aquifer system was proposed. Finally, the contribution of TEM (and electrical resistivity tomography) geophysical methods in studying complex coastal aquifers is shown by this work.     

Evaluation of alternative conceptual models for groundwater modelling

This study evaluates the alternative conceptual models for groundwater modelling. A true model was created with a synthetic alluvial fan-plain hydrogeological framework. Various alternative conceptual models were evaluated for groundwater flow simulations. The first alternative model is a single aquifer layer model; the second alternative model is a 3-layer aquifer model; and the third model is a 5-layer model consisting of 3 aquifers separated by 2 aquitards. All models could fit very well to the observations with optimized values of hydraulic conductivities. However, the single aquifer layer model can only compute water balance components with good accuracy. The 3-layer aquifer model can be used for water balance computation and groundwater head simulation with small errors. The 5-layer model is capable of simulating water budget, groundwater head distribution and travel times with high accuracy. Multi-model analysis found only the 3rd alternative model superior.