运用三维变密度潮汐效应模型确定滨海含水系统的海底等效边界——山东烟台夹河中下游地区为例

确定咸淡水界面的位置是滨海地区海水入侵研究的主要任务之一。对于天然条件承压含水层而言，含水层顶板向海底延伸的距离直接影响了咸淡水界面的位置，它可以通过承压含水层中地下水的潮汐效应信息来确定。考虑到咸淡水之间密度的差异，建立了山东省夹河中下游地区滨海含水系统地下水三维变密度潮汐效应模型。通过反复对比潮汐效应观测中的地下水水头波动与模型计算出的水头波动，确定了滨海承压含水系统的海底边界。同时，也初步估计出海区与近海陆区含水层的水文地质参数。     

利用观测孔中地下水位潮汐效应计算天津滨海新区含水层参数

在沿海地区,尤其是围海造陆工程形成的陆域地区地下水水位受潮汐影响较大,使传统水文地质试验求取含水层参数存在较大误差。因此通过合理概化地下水在潮汐作用下运动规律,建立数学模型,推导解析公式求取沿海含水层参数具有重要意义。分析天津滨海新区两处观测孔地下水位及潮汐波动特征,在滞后时间不明显的情况下,利用观测孔水位变幅数据计算了含水层水头扩散系数,并根据承压含水层储水系数经验值进一步获得含水层渗透系数。通过两个观测孔分别计算,对比计算结果互相验证发现,该方法取得了令人满意的结果。利用地下水潮汐效应计算含水层参数可以广泛应用于沿海地区水文地质工作中。     

用地下水的越流动态确定水文地质参数的初步研究

本文是在吉林西部白城地区进行地下水资源评价时,为了取得更多的点参数,根据该区水文地质条件和地下水的动态特征,经初步研究而提出的。 到目前为止,研究地下水流向集水构筑物的运动规律的一些主要理论公式都是以地下水流向井的水平运动为主要研究内容。并借助于抽水所产生的水位下降建立理论公式,确定水文地质参数。 本文试图讨论在由两个含水层(潜水含水层和下伏有不透水层的第一个半承压含水层)及其间的     

多层含水层抽水试验分层止水方法介绍

前言沙岭子电厂水源地进行水文地质勘测中,由于水源地地下水由潜水、承压水及承压自流水所组成的多层含水层类型。鉴于这样一个多含水层的结构特征,为取得各层含水层的水文地质参数及水位、水温、水质等资料,就必须解决井内、外分层止水进行抽水试验以及井内水位降深观测装置和观测孔的观测管分层安装、分层止水等一系列多工艺问题。     

永定东中石灰石矿区抽水试验

分析抽水试验观测记录数据,计算确定含水层水文地质参数,判断地下水补给来源、影响范围、含水层富水性和矿区地下水运动变化规律,为矿山生产建设科学发展形成技术依据要求。应用永定东中石灰石矿区抽水工作实例,为正确分析单孔稳定流抽水试验资料并建立合适的Q-S曲线方程以预测未来水平矿坑涌水量提供借鉴。     

承压水潮汐荷载效应研究进展

地下水潮汐效应是滨海地区水文地质研究的一个重要方面,潮汐荷载效应是承压水潮汐效应的动力学机制之一。介绍了国内外潮汐荷载效应的研究进程,系统地总结了解析法和数值法在潮汐荷载效应研究中的应用及取得的成果。解析法是潮汐荷载效应研究的主要方法,该方法多通过荷载系数来刻画潮汐荷载效应,能够较清晰、准确地揭示复杂含水层结构、多种影响因素、复杂边界条件下的地下水潮汐效应规律;数值法主要通过三种方法(现有地下水水流数值模型改进法、荷载系数法和半数值法)来实现对潮汐荷载效应的研究。在现有潮汐荷载效应研究中,上述两种方法都很少考虑应力和孔隙渗流的耦合。程序化、高效化和可视化使得数值法仍应作为滨海区复杂水文地质条件下地下水动态的主要研究方法。为此,通用的地下水数值模拟代码(MODFLOW)需要改进,使之能够实现对潮汐荷载效应的刻画。另外,承压水作为滨海区(如:大陆岛)重要的淡水资源,人为开采是地下水动态的主要影响因素之一,滨海区地下水动态研究应综合考虑人为因素和潮汐荷载效应的共同作用。     

利用滨海地下水系统对双频潮汐的响应识别含水层压力传导系数

含水层的压力传导系数是地下水溶质运移及扩散的重要参数,通常压力传导系数由抽水试验确定,即通过抽水给地下水系统一个扰动,然后监测地下水水位的响应,进而计算压力传导系数。对滨海地区含水层来说,潮汐波动就是天然的抽(注)水试验。潮汐由不同频率的分潮构成,前两项主频潮汐基本可代表潮汐主要特征。利用地下水水位对潮汐的响应,识别含水层压力传导系数。在推导出地下水对双频潮汐响应的解析解的基础上,提出了潮汐信号衰减法识别压力传导系数,并做了敏感性分析。通过对数值解叠加高斯噪音信号进行参数估算,证明该方法估算压力传导系数具有较高的可靠性和实用性。     

北京市朝阳区沙子营垃圾污染物在潜水含水层中迁移规律的现场模拟实验

选择北京平原区水文地质、环境地质等方面都比较典型的沙子营垃圾堆放场,建立了其水文地质模型。在充分收集资料、分析得出部分计算所需参数后,采用现场弥散实验、勘查取样测试等方法,求得了该含水层的弥散系数等参数;用二维非稳定流溶质运移方程对污染物在此含水层中的迁移扩散规律、速度和污染范围等进行了模拟计算;采用现场钻探、取样测试分析等方法,评价了该含水层的实际污染状况。实验模拟计算和现场调查结果表明:污染物在潜水含水层中的运移规律遵循二维非稳定流场中的溶质运移方程,污染物在潜水含水层中的运移速度约为86.25m/a,迁移扩散主要发生在地下水流向上,侧向扩散宽度极小,是地下水流向上的1/17。     

抽水试验确定承压含水层参数方法探讨

地下水资源评价与地下水可开采量计算,需要对地下含水层组参数进行分析确定。本文结合鲁北地区两个地下水水源地采用抽水试验确定承压含水层水文地质参数的实践,探讨定流量(单孔或多孔)抽水试验确定含水层参数的可行性,并对定降深抽水试验确定水文地质参数方法进行了探索。     

均质潜水含水层抽、灌水实验测求参数的对比

本文设计了在室内均质潜水含水层抽、灌地下水的实验，对实验准备工作和实验过程进行了规划，使用观测的实验数据分别对水文地质参数进行了对比计算和规律分析。最后对实施人工回灌的具体工作进行了归纳总结。     

A data-driven approach for analyzing dynamics of tide–aquifer interaction in coastal aquifer systems

Analysis of tidal effects on aquifer systems plays an important role in coastal aquifer management owing to various hydrological, engineering and environmental problems in coastal areas. Using the real-world data of unconfined and confined aquifers, a data-driven approach is presented in this study for the analysis of tide–aquifer interaction in coastal aquifers. Six analytical tide–aquifer interaction models were selected which take into account the effects of vertical beach, sloping beach, tidal loading, aquifer leakage, outlet capping, and combined leakage and outlet capping on tide-induced groundwater fluctuations. The tide–aquifer interaction datasets were obtained from the Konan groundwater basin (unconfined aquifer) of Japan and the Dridrate groundwater basin (confined aquifer) of Morocco. The analysis of the results obtained by the sloping beach model revealed that for a given beach slope, the amplitude of groundwater level increases with an increase in aquifer diffusivity and a decrease in aquifer thickness. However, no significant effect of beach slope was observed in this study at unconfined sites for all the datasets. The influence of tidal loading was found to be considerably less for all the three confined sites. Further, the analysis of the results of the leakage model indicated that with an increase in leakage into the aquifer, the amplitude of groundwater level as well as the phase shift (time lag) decreases. Of all the confined and unconfined datasets, only two confined sites were found to be affected by outlet capping. Overall, it is concluded that the coastal beach bordering the Konan basin is not significantly sloping, the contribution of tidal loading to tide-induced groundwater fluctuations in the Dridrate aquifer is not appreciable, and that the aquifer leakage and outlet capping do not exist at the unconfined sites under investigation.     

Groundwater mixing in a sand-island freshwater lens: density-dependent flow and stratigraphic controls

This paper focuses on a small back-barrier sand-island on the southeast coast of Queensland. The freshwater lens in the study area exhibits anomalously high short-range salinity gradients at shallow depths, which cannot be explained using a standard seawater intrusion model. The island groundwater system consists of two aquifers: a semiconfined aquifer hosting saline to hypersaline groundwater and an overlying unconfined freshwater aquifer. The deeper aquifer is semiconfined within an incised paleovalley, and groundwater flow is restricted to an east – west direction. Tidal response observations show that the tidal signal propagates far more rapidly and is of much higher magnitude in the semiconfined aquifer than the unconfined aquifer. The tidal wave-pulse amplitude is also subject to greater attenuation in the unconfined aquifer. A conceptual hydrogeological model illustrates how upwelling of hypersaline groundwater, induced by density-dependent flow and tidal pumping, has contaminated the shallow groundwater resource. Salinisation at shallow depths is restricted to an area proximal to the paleovalley aquifer. The spatial distribution of lithological heterogeneity is an initial limiting control on the movement of the upwelling saline plume. The extent of shallow groundwater contamination is also limited by the presence of a baroclinic field, resulting from lateral variations in fluid density. Hydrochemical signatures have been used to support the model hypothesis and link the salinisation of fresh groundwater with the semiconfined aquifer as opposed to the surrounding estuarine surface water. The geometry and thickness of the freshwater lens are further controlled by the presence of the largely impermeable bedrock paleosurface between 9 and 12 m depth. The combination of hypersaline groundwater and hydraulically restrictive lithology at shallow depths has produced excessive thinning of the freshwater lens, demonstrating that the application of a model such as the Dupuit – Ghyben – Herzberg relationship would grossly overestimate the available groundwater resource.     

Aquifer parameter estimation using tide-induced water-table fluctuations in the Biscayne Aquifer,Miami-Dade County,Florida (USA)

The Biscayne Aquifer (Florida, USA) is a coastal, shallow, unconfined, and heterogeneous aquifer with high water tables, composed of less-permeable sand to highly permeable karstic limestone. These properties make the Biscayne Aquifer one of the world’s most productive groundwater resources. The aquifer’s high yield and non-Darcian flow cause challenges for estimating aquifer parameters, which are essential for understanding groundwater processes and managing and protecting the groundwater resources. Water-table fluctuations in the Biscayne Aquifer are associated with astronomical tidal forces and gate operations on canal water-control structures. Analysis of observed groundwater level fluctuations can provide an understanding of the connectivity between the aquifer, Biscayne Bay, and the water level in the canals. Further, groundwater level fluctuations can be used for aquifer parameter estimation. In this research, observed ocean water levels measured at tidal stations and groundwater levels are fitted to Jacob’s analytical solution, where the amplitude of the groundwater head fluctuation decreases exponentially, and the time lag increases with distance from the shore. Observed groundwater levels were obtained from monitoring wells along the Miami-Dade shore and the barrier island of Miami Beach. Results indicate that Jacob’s solution is effective for aquifer parameter estimation in Miami Beach, where monitoring wells are closer to the shore. Estimated hydraulic conductivity appears to increase by four orders of magnitude to approximately 1 m s^–1 as the distance from shore increases. Constructing monitoring wells closer to the shore in Miami-Dade County and elsewhere would permit improved aquifer parameter estimation and support enhanced groundwater modeling efforts.

     

基于土柱排水实验探讨潜水含水层给水度的影响因素

给水度是重要的水文地质参数之一,为探讨潜水含水层给水度的影响因素,设计均质粗砂、均质细砂、上细下粗、上粗下细4种土柱的排水实验,按照一定间距分段降低地下水位,每段降深设计不同排水时间,求取不同水位埋深的给水度.当地下水位埋深小于支持毛细水高度,给水度受地下水位埋深影响显著,其关系可以用非线性函数来表示.排水时间越长,给水度越大,当土柱分段排水时间超过1 h,给水度值稳定.0.6~0.9 mm粒径的粗砂给水度大于0.2~0.4 mm的细砂;当地下水位在土层分界面下20 cm时,上细下粗土柱给水度显著增大,上粗下细土柱给水度显著减小.给水度的大小与地下水位埋深、排水时间、岩土颗粒和潜水面附近及之上非均质结构等因素有关.      

吉林省白城市地下水数值模拟及其变化趋势

为了给吉林省白城市地下水资源合理开发利用提供科学依据,运用地下水模拟技术对白城市地下水位变化趋势进行预测。首先在分析研究区水文地质条件的基础上建立水文地质概念模型,将含水系统概化为潜水含水层、越流层和承压含水层,并确定了系统的边界条件。运用地下水模型软件(GMS)对研究区地下水流进行了模拟和预测,预测结果表明在扩大开采的情况下,潜水含水层水位总体下降约1~2 m;随着开采量的增加和时间的增长,含水层在洮儿河扇形地、洮南、研究区东北部、镇赉、大安南部等地相继出现疏干现象,最大降幅达5~6 m。承压含水层在大安南部、镇赉北部等地水位下降较大,但是流场整体降幅不大,总体降幅在0.5 m以内。最后提出了地下水合理开发利用建议。     

基于地下水流数值模型的改进DRASTIC方法

地下水脆弱性评价作为地下水资源保护和地下水开发利用规划的一个重要工具,被广泛的应用于实际工作中。尝试利用地下水数值模型为改进的DRASTIC方法提供数据支持,并以北京市平原区为例探讨地下水脆弱性评价方法。评价结果与传统方法在高值区和低值区具有很好的对应性,而基于模型的方法在地下水水位计算、含水层介质和水力传导系数确定上较传统方法更具优势,如地下水位的计算上较传统方法更为客观地体现含水介质对地下水运动的影响,且能够方便地获得模拟期内任意时间的流场数据;经由模型调试后的含水层参数数据,较传统方法更为准确。评价结果分区之间的变化较传统方法更为平滑,更符合水文地质条件渐变的特性。     

Global curve-fitting for determining the hydrogeological parameters of leaky confined aquifers by transient flow pumping test

Proper management of groundwater resources requires an accurate evaluation of the parameters (hydraulic properties) that control the movement and storage of groundwater. Hydrogeological parameters are the basis of groundwater evaluation, modeling, and management and so on. A global curve-fitting method incorporating pumping test data and water table recovery data was introduced in the present study. The principal and procedures of the method were elucidated in detail. The drawdown and recovery data from two sets of transient flow pumping test conducted in no. 2 water source site of Shizuishan city were used to verify the calculation accuracy of the proposed method. The hydrogeological parameters were also estimated with traditional type curve-fitting method on the basis of formula derived by Hantush and Jacob. The hydrogeological parameters calculated by the two methods were compared and the results show that the parameters obtained by the global curve-fitting method are a little bigger than but very close to those obtained by the traditional type curve-fitting method. The proposed method which possesses three major advantages is feasible and reliable in aquifer parameter identification. A comparative study on various methods for parameter identification is required and expected in future study.     

Determination of the parameter pattern and values for a one-dimensional multi-zone unconfined aquifer

In an aquifer, heterogeneity plays an important role in governing groundwater flow. Hence, aquifer characterization should involve both the pattern and values of the hydrogeological parameters. A new analytical solution describing the one-dimensional groundwater flow in a multi-zone unconfined aquifer is presented, and a methodology developed from the analytical solution and a heuristic approach for determining the pattern and values of the aquifer parameters are proposed. The analytical solution demonstrates that the hydraulic head varies spatially and is influenced by aquifer heterogeneity. Simulated annealing, a heuristic approach, is incorporated with the solution to simultaneously identify the pattern and values of the hydraulic conductivity for a horizontal multi-zone unconfined aquifer. This approach may be used to give an approximate result for a two-dimensional problem by dividing the model area into a number of transects along the transverse direction, identifying the parameter values along the longitudinal direction for each transect, and then smoothing the identified results.     

利用岸边水头动态确定含水层在临海方向上的边界

本文建立了潮汐效率与越流共同作用条件下承压水头动态的数学模型.利用复变函数的基本原理,考虑到地下水流动的连续性,推导出了相应的承压水头动态公式.以此动态公式为基础,提出了直接利用承压水头的潮汐动态来确定延伸到海底的含水层长度(边界条件)的计算方法.最后结合广西北海市的地下水头长观资料,对上述方法进行了应用和分析研究.