基于不同地质统计方法的渗透系数场对污染物运移的影响

渗透系数场的空间变异性是影响污染物运移结果的决定因素,而地质统计方法是解决渗透系数空间变异性的主要技术手段。本文利用野外场地实测数据,采用普通克里格法和指示克里格法、顺序高斯模拟法和顺序指示模拟法四种地质统计方法,插值估测和模拟再现随机渗透系数场,进而对比研究四种渗透系数场对大尺度污染物运移的影响。研究结果表明,污染羽的质心位置（一阶矩）主要由渗透系数的平均值来决定;污染羽在空间上的展布范围（二阶矩）主要受渗透系数空间变异方差的影响;条件模拟克服了估计法的平滑效果,较好地再现真实曲线的波动性,渗透系数（ lnK）估计方差与污染羽空间二阶矩随着条件模拟次数的增加而减小,并且顺序指示模拟程度更加明显。     

Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.     

Sensitivity Analysis of a Combined Groundwater Flow and Solute Transport Model Using Local-Grid Refinement: A Case Study

Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.     

含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响

为探讨含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响,基于非均质含水层砂箱实验,分别用传统等效均质模型、克立金插值和水力层析刻画含水层渗透系数场,并探讨了先验信息对水力层析结果的影响.将不同方法估算的渗透系数场用以预测地下水流和溶质运移过程,以此判断不同方法估算结果的优劣,分析含水层非均质性对地下水流和溶质运移的影响.结果表明:与克立金插值法相比,水力层析法可以更好地刻画含水层非均质性,较准确地预测地下水流和溶质运移过程;钻孔岩心渗透系数样本值作为先验信息可以提高水力层析法估算结果的精度;传统等效均质模型无法准确预测地下水流和溶质运移过程.含水层非均质性的增强将导致溶质污染羽分布形态和运移路径的空间变异性增强,并且优势通道直接决定溶质的分布及运移路径.      

扬-泰-靖地区地下水系统水力联系与硫酸盐污染特征

本文对长江三角洲扬-泰-靖地区第四系松散层地下水中环境同位素(D、18O、34S)的分布特征进行了分析,旨在揭示大气降水、长江水、潜水及承压水之间的水力联系,辨别地下水中硫酸盐的来源及其污染状况。研究结果表明潜水含水层接受大气降水及长江水的补给,硫酸盐主要为农业污染来源或与海源硫酸盐的混合。承压含水层主要接受大气降水的补给,与潜水含水层及长江之间的水力联系较差,硫酸盐来源不同。在研究区顶部和沿江地段的浅层孔隙承压水中,硫酸盐来源于硫化物的氧化;在东部的深层孔隙承压水中,硫酸盐主要来源于硫酸盐岩的溶解或海源硫酸盐的滞留,基本未受到潜水或地表水中硫酸盐的污染。     

Field characterization and data integration to define the hydraulic heterogeneity of a shallow granular aquifer at a sub-watershed scale

Providing a sound basis for aquifer management or remediation requires that hydrogeological investigations carried out to understand groundwater flow and contaminant transport be based on representative data that capture the heterogeneous spatial distribution of aquifer hydraulic properties. This paper describes a general workflow allowing the characterization of the heterogeneity of the hydraulic properties of granular aquifers at an intermediate scale of a few km². The workflow involves characterization and data integration steps that were applied on a 12-km² study area encompassing a decommissioned landfill emitting a leachate plume and its main surface water receptors. The sediments composing the aquifer were deposited in a littoral–sublittoral environment and show evidence of small-scale transitional heterogeneities. Cone penetrometer tests (CPT) combined with soil moisture and electrical resistivity (SMR) measurements were thus used to identify and characterize spatial heterogeneities in hydraulic properties over the study area. Site-specific statistical relationships were needed to infer hydrofacies units and to estimate hydraulic properties from high-resolution CPT/SMR soundings distributed all over the study area. A learning machine approach was used due to the complex statistical relationships between colocated hydraulic and CPT/SMR data covering the full range of aquifer materials. Application of this workflow allowed the identification of hydrofacies units and the estimation of horizontal hydraulic conductivity, vertical hydraulic conductivity and porosity over the study area. The paper describes and discusses data acquisition and integration methodologies that can be adapted to different field situations, while making the aquifer characterization process more time-efficient and less labor-intensive.     

山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究

在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,?     

云应盆地东北部含水层结构特征及地下水转化模式

云应盆地东北部属鄂北贫水地区,赋存于古近系—第四系含水层中的地下水是当地生产、生活用水的主要来源,亟需查明含水层的结构、含水层间地下水的转化关系等基本条件,为研究区内合理开发利用地下水提供依据。本研究通过野外水文地质调查、水文地质钻探工作,将研究区划分为单层含水层与双层含水层结构两个亚区（6个小区）。并通过地下水水位动态长期监测,获取了区内不同含水层的水位动态变化特征,分析各含水层之间的水力联系,建立了区域地下水转化的概念模式,即:研究区地下水以接受山前降雨入渗及风化裂隙水侧向径流补给为主,主要以水平径流的形式经古近系孔隙-裂隙含水层及第四系孔隙承压含水层往澴水方向运移,而后进入第四系孔隙潜水含水层。地下水和地表水在不同季节补排模式不同,雨季地表水（澴水）补给地下水,旱季地下水向地表水（澴水）排泄。古近系孔隙-裂隙水与上覆第四系孔隙水联系密切互为补给,共同构成具有统一水力联系的垂向多层结构的含水系统。独特的含水层结构决定了区内地下水接受降水补给的条件较差,地下水可开采资源量总体较贫乏,建议重点利用区域地表水资源,适度开发地下水资源,推进农业节水灌溉工程,实现水资源可持续利用。     

The dependence of stream depletion by seasonal pumping on various hydraulic characteristics and engineering factors

Compensation pumping is used to alleviate deficiencies in streamflow discharge during dry seasons. Short-term groundwater pumping can use aquifer storage instead of catchment-zone water until the drawdown reaches the edge of the stream. The capacitance is a complex, dimensionless parameter of an aquifer system that defines the delayed effect on streamflow when there is groundwater pumping. This parameter is a function of aquifer hydraulic characteristics, pumping time, and distance between the well and stream edge; the latter can involve stream leakance and vertical leakance of an associated aquitard. Three typical hydraulic cases of combined water systems (major catchment-zone wells close to the stream and compensation pumping wells) were classified depending on their capacitance structure (i.e. the relationship between surface water and groundwater): (1) perfect hydraulic connection between the stream and aquifer; (2) imperfect hydraulic connection between the stream and aquifer; and (3) essentially imperfect hydraulic connection between the stream and the underlying confined aquifer. The impact of various hydraulic characteristics and engineering factors on stream depletion was examined by conceptual and numerical modeling. To predict the suitability and efficiency of a combined water system application, regression tests were undertaken on unit stream depletion and capacitance, and power dependencies were defined.     

A geoelectrical and hydrogeological study for the assessment of groundwater resources in Wadi Al Bih, UAE

Drilling information, historical water table levels, groundwater salinity records of the existing water wells in Wadi Al Bih area, United Arab Emirates, were stored in a geodatabase and used to characterize the geological and hydrogeological settings of this area. A 2D earth resistivity imaging survey was conducted for the first time in the Northern UAE to determine the potential of the Quaternary aquifer and its groundwater quality in the areas where there are no monitoring wells. The results of the chemical analyses of the collected groundwater samples together with the inversion results of the resistivity data were used to draw a total salinity map and determine the spatial variations in groundwater quality. The inversion results of the 2D earth resistivity imaging data indicated that the Quaternary aquifer in the study area is in a good connection with the underlying carbonate aquifer. It also indicated that the carbonate aquifer is of major regional and vertical extension and it contains the fresh water in this area. The data stored in the developed database were used to produce different types of geopotential maps.     

Using lithologic modeling techniques for aquifer characterization and groundwater flow modeling of the Sohag area, Egypt

Three dimensional lithologic modeling techniques have been used for detailed characterization and groundwater flow modeling of the Quaternary aquifer system of the Sohag area, Egypt. Well log data were used for building the lithologic model using RockWorks. A groundwater flow model, facilitated by MODFLOW 2000, was built using results of the lithologic model. The obtained lithologic model honored the well log data and revealed a complex sedimentary system, which is mainly composed of six lithologic categories: clay, clay and sand, fine sand, coarse sand, sand and gravel, and gravel. Inter-fingering and presence of lenses are the main characteristics of the sedimentary basin in the study area. A wide range of hydraulic conductivities is present, which vary spatially and control the groundwater flow. Heterogeneity of the aquifer system is spatially represented where different hydraulic conductivities are found in the different directions. Sandy layers tend to be connected. Hydraulic continuity is represented by inter-fingering and connection of sandy materials within the aquifer system.     

Mapping the hydraulic connection between a coalbed and adjacent aquifer: example of the coal-seam gas resource area,north Galilee Basin,Australia

Coal-seam gas production requires groundwater extraction from coal-bearing formations to reduce the hydraulic pressure and improve gas recovery. In layered sedimentary basins, the coalbeds are often separated from freshwater aquifers by low-permeability aquitards. However, hydraulic connection between the coalbed and aquifers is possible due to the heterogeneity in the aquitard such as the existence of conductive faults or sandy channel deposits. For coal-seam gas extraction operations, it is desirable to identify areas in a basin where the probability of hydraulic connection between the coalbed and aquifers is low in order to avoid unnecessary loss of groundwater from aquifers and gas production problems. A connection indicator, the groundwater age indictor (GAI), is proposed, to quantify the degree of hydraulic connection. The spatial distribution of GAI can indicate the optimum positions for gas/water extraction in the coalbed. Depressurizing the coalbed at locations with a low GAI would result in little or no interaction with the aquifer when compared to the other positions. The concept of GAI is validated on synthetic cases and is then applied to the north Galilee Basin, Australia, to assess the degree of hydraulic connection between the Aramac Coal Measure and the water-bearing formations in the Great Artesian Basin, which are separated by an aquitard, the Betts Creek Beds. It is found that the GAI is higher in the western part of the basin, indicating a higher risk to depressurization of the coalbed in this region due to the strong hydraulic connection between the coalbed and the overlying aquifer.     

Influence of Hydraulic Conductivity and Wellbore Design in the Fate and Transport of Nitrate in Multi-aquifer Systems

Nitrate concentrations in multi-aquifer systems are heavily affected by the presence of wellbores (active or abandoned) that are screened in several aquifers. The spatial variability of hydraulic conductivity in the confining layers has also an important impact on the concentrations. A synthetic three-dimensional flow and transport exercise was carried in a multi-aquifer system consisting of two aquifers separated by an aquitard in which 100 vertical wellbores had been drilled. To model the wellbores and the flow and transport connection between aquifers that they may induce, we assign a high vertical hydraulic conductivity and a low effective porosity to the cell blocks including the wells. With these parameters, a solute will travel quickly from one aquifer to the other without being stored in the well itself. The wellbores will act as preferential pathways, and the solute will move quickly between aquifers according to the hydrodynamic conditions. Not considering these preferential pathways could induce erroneous interpretations of the solute distribution in an aquifer. We also noted that when there are vertical wellbores that connect aquifers in a multi-aquifer system, low conductivity in the aquitard enhances the flow of solute through the wellbores. Time-varying pumping rates induce important fluctuations in nitrate concentrations; therefore, any estimate of the water quality of the aquifer will depend on the moment when the data has been recorded. Consequently, concentration maps obtained by interpolation of point samples are seldom a good indicator of the chemical status of groundwater bodies; alternatively, we recommend complementing the usual interpolated maps with numerical models to gain a true understanding of the spatial distribution of the solute concentration.     

孟加拉国巴拉普库利亚煤矿含水层水力联系研究

为进一步查明孟加拉国巴拉普库利亚煤矿水文地质特征及含水层间水力联系,为矿井水害防治提供理论依据,以19组水质化验数据为基础,结合含水层及隔水层空间展布特征、井田构造特征、水位历时曲线、水化学类型、氢氧同位素特征等,综合分析新近系UDT含水层、Ⅵ煤顶板含水层、Ⅵ煤含水层之间的水力联系。结果表明,井田北部LDT隔水层局部缺失,为UDT含水层水向含煤地层补给提供了条件;井田北翼煤层顶板含水层与UDT含水层水位变化规律密切相关,且水位相近,初步证明两者存在水力联系;各含水层水均为HCO₃-Ca·Na·Mg型,均为低矿化度水,进一步证明各含水层间存在水力循环;聚类分析结果表明各含水层水质存在一定关联度,推演含水层间水力联系程度;氢(δD)氧(δ18O)同位素特征点分布于全球大气降水线附近,表明大气降水是各含水层共同的补给水源。研究成果可以指导孟加拉国巴拉普库利亚煤矿Ⅵ煤层开采时水害防治方向。      

水力层析法与克立金法估算非均质含水层渗透系数场比较

含水层非均质性空间分布特征的识别,是准确模拟地下水流和污染物运移的前提.基于室内非均质含水层砂箱实验,分别利用水力层析法和克立金插值法刻画了非均质含水层渗透系数场.研究结果表明:(1) 水力层析法与克立金法相比,不仅可以更好地刻画非均质含水层渗透系数场,还可以更高精度地预测地下水流过程;(2) 水力层析抽水实验中,通过增加抽水实验组数可以有效地提高水力层析参数反演的精度,但是抽水实验组数增加到一定程度以后,再增加抽水组数不会显著提升参数反演的效果.后续需要进一步研究水力层析抽水实验合适的组数,进一步对抽水井进行优化布设.      

Numerical modeling of contaminant transport in a stratified heterogeneous aquifer with dipping anisotropy

The combined influence of dip angle and adsorption heterogeneity on solute transport mechanisms in heterogeneous media can be understood by performing simulations of steady-state flow and transient transport in a heterogeneous aquifer with dipping anisotropy. Reactive and non-reactive contaminant transport in various types of heterogeneous aquifer is studied by simulations. The hydraulic conductivity (K) of the heterogeneous aquifer is generated by HYDRO_GEN with a Gaussian correlation spectrum. By considering the heterogeneity of the adsorption distribution coefficient (K _d), a perfect negative correlation between lnK and lnK _d is obtained by using the spherical grains model. The generated K and K _d are used as input to groundwater flow and transport models to investigate the effects of dipping sedimentary heterogeneity on contaminant plume evolution. Simulation results showed that the magnitude of the dip angle strongly controls the plume evolution in the studied anisotropic and heterogeneous aquifer. The retarded average pore-water velocity (v/R) of the adsorption model significantly controls the horizontal spreading of the plume. The bottom plume is intensively retarded in the zones between the dipping lenses of lower hydraulic conductivity and the no-flow bottom boundary. The implications of these findings are very important for the management of contaminated heterogeneous aquifers.     

Flow and solute transport monitoring in the karst aquifer in SW Slovenia

The role of the unsaturated zone in the karst aquifer hydraulic behaviour was brought into focus in these studies of the catchment of the Hubelj spring (SW Slovenia). The variations of natural tracers in precipitation and in groundwater during a summer storm event made it possible to trace local flow and solute transport in the observed aquifer. The results produced data on the aquifer recharge, storage and discharge processes, as well as on mechanisms that affected them, which reflects a karst groundwater dynamics also at a regional scale. They point out the significance of effects of the fast preferential flow—epiflow that is the main factor controlling solute/contaminant transport towards the aquifer saturated zone. Numerous arguments indicate that the karst aquifer flow and solute transport mechanisms depend on the hydraulic behaviour of the epikarst zone.     

北京市丰台区东高地岩溶裂隙水与第四系地下水的水力联系

专家一般认为东高地地区寒武系岩溶裂隙水的补给来源以上游的佣向径流为主,与第四系地下水没有明显的水力联系;通过近来地下水位动态监测和岩溶裂隙水的同住素成分研究,我们认为寒武系岩溶裂隙水与第四系地下水之间有着明显的水力联系。     

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

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

Flow and transport model of a polder area in the Belgian coastal plain: example of data integration

Tracer data can serve to derive parameters related to groundwater systems or can be used in the calibration of transport and flow models. Temperature, ³H/³He (age dating) and hydrochemical data, together with head measurements, borehole logs and an aquifer test, were used to obtain insight into the flow and transport of a Quaternary aquifer system located in the Belgian polder area. Flow and transport model code SEAWAT acted as the integration medium for the different data. Each type of data has its own interpretation technique and adds components to the model. Additionally, different types of data aid in verification of the results. For instance, fluxes from temperature logs (used with the SEAWAT model), and water quality and age dating all provide information on flow to and velocities in the vicinity of drainage ditches. Different data also provide information on different scales. Temperature logs and aquifer tests act on a small scale, groundwater age is influenced by larger-scale flow, and water quality is determined by the general flow of the area. Integration of different kinds of geological, hydrological, geophysical and geochemical data is shown to be an important way forward in the efforts to model real-world cases.