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

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

用地下水潮汐效应确定潜水含水层水文地质参数

目前对潜水含水层地下水潮汐效应和水文地质参数求解方法的研究相对较少.通过对福建古雷半岛滨海潜水含水层地下水潮汐效应和海水潮汐动态的观测, 运用Fourier频谱分析方法确定了研究区海水潮汐波动方程(波动特征参数), 并以此作为地下水的边界条件, 推导了潜水含水层地下水潮汐效应的波动方程, 利用最小二乘法以地下水水位波动观测值为目标函数对潜水含水层的渗透系数与重力给水度的比值进行了反演识别, 为类似地区水文地质参数的确定提供了借鉴, 也为该区后续地下水中溶质迁移规律的研究奠定了基础.      

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

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

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

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

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

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

地下水易污性评价方法——DRASTIC指标体系

本文较详细地介绍了目前欧美国家在地下水易污性评价中所广泛采用的ＤＲＡＳＴＩＣ指标体系方法,对其中所含的地下水埋深,含水层的净补给,岩性,土壤类型,渗流区介质,水力传导系数７个参数的评分以及在评分过程中所应注意的问题进行了具体的阐述,最后简要介绍了ＤＲＡＳＴＩＣ易污性指标体系法在大连沿海地区地下水易污性评价中应用的情况,实际应用表明该方法可适用于我国广大地区的地下水易污性评价工作。     

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

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

江苏姜堰孔隙承压地下水资源评价三维数值模拟

通过建立符合该地区承压地下水赋存条件和渗流特征的地下水三维数值模型,并在模型识别、验证的基础上,根据各承压含水层的水位控制要求,调整现有地下水的开采布局,以乡镇为单位,确定出各含水层的地下水可开采资源量及最优开采布局,并预测2011—2020年底逐年地下水水位动态变化情况。     

山东烟台夹河中、下游地区海水入侵三维水质数值模拟研究

建立了三维变密度对流弥散水质数学模型来研究山东省烟台夹河中、下游地区咸淡水界面的运移规律。以四面体为基本离散单元 ,推导出三维海水入侵变密度水质模型求解的数值方法 ,其中水流方程求解时运用了迦辽金有限单元法。溶质运移方程求解时运用了欧拉拉格朗日混合方法 ,将对流项与弥散项分离 ,用传统迦辽金有限元方法求解弥散项 ;采用自适应MOC MMOC法求解对流项 ,以消除人工过量和数值弥散。根据地下水的潮汐效应观测信息 ,确定了含水系统的海底延伸边界 ;利用该地区地下水水头及水质长观资料识别了模型的水文地质参数 ,探讨了夹河地区海水入侵的原因 :认为夹河下游地区滨海地带地下水过量开采是造成烟台地区海水入侵的主要原因。此外 ,海水随潮定期地倒灌进入夹河 ,通过局部岩性天窗侵入淡水含水层加剧了沿夹河河床两侧地下水的咸化。同时还预测了几种情况下地下水的水质演化趋势 ,为防止和减轻夹河地区海水入侵提供合理、科学的依据。     

大陆岛地下水动力学特征—以湛江东海岛为例

东海岛是一具有独特水文地质条件的大陆岛,浅层含水层与大陆以浅海湾相隔,中深层承压水含水层与大陆地下水系统紧密相连。为了深入地认识大陆岛地下水水动力学特征,以湛江东海岛为例,阐述了其水文地质条件,并分析了东海岛浅、中、深层地下水的流场和动态特征。分析结果表明,东海岛为一个典型且独特的大陆岛,岛内和大陆的部分浅层含水层由湛江湾相隔,岛内中、深层含水层和大陆中、深层含水层通过湛江湾相连,且具有统一的水位分布,并保持着密切的水力联系,岛内中、深层地下水由南向北径流补给湛江市区的降落漏斗中心;滨海及海水区域浅层含水层及其下伏的粘土层构成了防止海水入侵中、深层地下水的保护层;浅层地下水流场基本保持天然状态,水位动态特征主要为入渗径流型,水位变化与降雨量相关;中、深层地下水流场以人工流场为主,地下水由南向北径流,水位动态类型主要为开采动态型,水位变化主要受到开采量变化的影响;在近海岸地区,地下水动态表现为潮汐效应型,在潮汐作用下,地下水位动态具有周期性。     

Groundwater response to dual tidal fluctuations in a peninsula or an elongated island

Groundwater flow in a peninsula or an elongated island is influenced by tidal fluctuations on both sides of the peninsula or island, which is named as dual tidal fluctuations. In this study, semianalytical solutions of transient groundwater flow in response to dual tidal fluctuations in an aquifer–aquitard system were presented for cases with and without the aquitard storage. These solutions were first derived using the Laplace transform and subsequently computed by the Fourier series numerical inverse Laplace transform. The derived solutions were found to agree very well with the results of numerical simulations by MODFLOW. The solution ignoring the aquitard storage approached the quasi‐steady state solution quickly when the mean sea level initial condition was used. The solutions with and without the aquitard storage were nearly the same at the early time and were separated from each other during the intermediate time, and the difference of solutions became constant at late time for small aquifer/aquitard storativity ratio and large tidal frequency. The propagation bias, which is the departure from the theoretical ratio of tidal attenuation to tidal lag, was enhanced not only with increase of the dimensionless specific leakage (aquitard/aquifer hydraulic conductance ratio) but also with decrease of the aquifer/aquitard storativity ratio and with the increase of the dimensionless tidal frequency. The solution with the aquitard storage was more sensitive to these three parameters. The newly developed solutions were capable of handling realistic initial conditions that might be approximated by piecewise linear functions. Copyright © 2012 John Wiley & Sons, Ltd.     

含水介质对地下水污染曝气修复的影响

为进一步研究地下水曝气原位修复技术的影响因素,采用室内实验的方法系统地研究了含水介质对原位曝气技术修复地下水污染的影响。研究结果表明:大量水相中的柴油在曝气的开始阶段被去除,随着时间的延长,污染物去除率不断增加,但增加的幅度逐渐减小。污染物去除率与时间关系曲线符合对数曲线规律,相关方程为y=alnx+b,R2=0.802 4～0.907 1,相关性较好。含水介质的渗透系数对地下水污染修复的影响较大,渗透系数与污染物的去除率基本呈正相关关系,渗透系数越大,污染物的去除率越大;含水介质密度与污染物去除率基本呈负相关关系,含水介质密度越大,污染物的去除率越小。     

Groundwater sustainability assessment in coastal aquifers

The present work investigates the response of shallow, coastal unconfined aquifers to anticipated overdraft conditions and climate change effect using numerical simulation. The groundwater flow model MODFLOW and variable density groundwater model SEAWAT are used for this investigation. The transmissivity and specific yield estimated from the existing database range from 10 to 810 m ²/day and 0.08% to 10.92% respectively. After successful calibration with Nash–Sutcliffe efficiency greater than 0.80, the values of horizontal hydraulic conductivity and specific yield of the unconfined aquifer were set in the range 1.85–61.90 m/day and 0.006–0.24 respectively. After validating the model, it is applied for forecasting the aquifer’s response to anticipated future scenarios of groundwater draft, recharge rate and sea level rise. The findings of the study illustrate that saltwater intrusion is intensified in the area adjoining the tidal rivers, rather than that due to the sea alone. Of all the scenarios simulated, the immense negative impact on groundwater quality emerges due to overdraft conditions and reduced recharge with the areal extent of seawater intrusion exceeding about 67% (TDS >1 kg/m ³). The study also arrives at the conclusion that, regional sea level rise of 1 mm/year has no impact on the groundwater dynamics of the aquifer.     

Large-scale hydraulic conductivity distribution in an unconfined carbonate aquifer using the ocean tidal propagation

The hydraulic conductivity of an unconfined carbonate aquifer at the uplifted atoll of Minami-Daito, Japan, was evaluated by a combination of cross-spectral analysis, analytical solution, and density-dependent groundwater modeling based on observed groundwater levels in 15 wells and at sea level. The island area was divided into 10 subregions based on island morphology and on inland propagation of ocean tides. The hydraulic conductivity was obtained for each subregion using analytical solutions based on phase lags of M₂ constituents of ocean tides at each well by assuming two aquifer thicknesses (300 and 1,800 m) and two effective porosities (0.1 and 0.3). The density-dependent groundwater model evaluated the hydraulic conductivity of the subregions by reproducing observed groundwater levels. The hydraulic conductivity in the subregions was estimated as 3.46?×?10^?3 to 6.35?×?10^?2 m/s for aquifer thickness of 300 m and effective porosity of 0.1, and as 1.73?×?10^?3 to 3.17?×?10^?2 m/s for aquifer thickness of 1,800 m and the effective porosity of 0.3. It was higher in southern and northern areas, and higher in interior lowland than in the western and eastern areas. Fissures and dolomite distributions on the island control differences of the omnidirectional ocean tidal propagation and cause these differences in hydraulic conductivity. The method used for this study may also be applicable to other small islands that have few or no data for hydraulic conductivity.

     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.     

Control of sea-water intrusion by salt-water pumping: Coast of Oman

A shallow alluvial coastal aquifer in the Batinah area of Oman, with sea-water intrusion that extends several kilometres inland, has been studied experimentally, analytically and numerically. The water table is proved to have a trough caused by intensive pumping from a fresh groundwater zone and evaporation from the saline phreatic surface. Resistivity traverses perpendicular to the shoreline indicated no fresh groundwater recharge into the sea. Using an analytical Dupuit-Forchheimer model, developed for the plain part of the catchment, explicit expressions for the water table, sharp interface location and stored volume of fresh water are obtained. It is shown that by the pumping of salt water from the intruded part of the aquifer, this intrusion can be mitigated. Different catchment sizes, intensities of fresh groundwater pumping, evaporation rates, water densities, sea level, incident fresh water level in the mountains and hydraulic conductivity are considered. SUTRA code is applied to a hypothetical case of a leaky aquifer with line sinks modeling fresh water withdrawal and evaporation. The numerical code also shows that pumping of saline water can pull the dispersion zone back to the shoreline.