The effects of aquifer heterogeneity on the 3D numerical simulation of soil and groundwater contamination at a chlor-alkali site in China

The simulation of groundwater flow and solute transport at contaminated sites often neglects the important influence that aquifer heterogeneity can have on the sub-surface distribution of contaminants. In this paper, the method of transition probability for geological statistics (T-PROGS) included in the Groundwater Model System (GMS) was applied to a chlor-alkali-contaminated site that was sampled with 68 soil borings and 15 groundwater monitoring wells. A 3-D groundwater numerical model and solute transport model was developed that was constrained by soil and groundwater data from the site. The spatial distribution of chloroethylene concentrations was simulated for a number of times using the levels measured in the field as a baseline. The results of these simulations showed that shapes and distribution of contaminant plumes are irregular both vertically and horizontally. The solute-transport simulations indicated that much of the contamination will preferentially move in groundwater through silt and fine-sands whereas flow is largely blocked in clays. Consequently, fine sand and silts become the most seriously polluted zones at the site, whereas, areas underlain by clays are largely uncontaminated. Heterogeneous lithologies beneath a site increase the complexity of coupling simulations of soil and groundwater.     

纵向弥散作用与渗透介质非均质性定量关系的模拟研究

在地下水污染模拟预报中,弥散参数是很难确定的一个模型参数。因实验室小尺度弥散规律一般不能用于大尺度弥散过程,而野外示踪试验却耗资大、周期长,限制了其实用性。文中利用随机数值模拟手段、基于随机理论的蒙特卡罗方法及序贯高斯模拟技术来生成渗透系数随机场,并研究渗透系数对数场的方差、相关长度以及变异函数类型在不同尺度上对纵向弥散度的影响,进而建立纵向弥散度与随机分布渗透系数场的方差和相关长度的统计定量关系,并与Gelhar理论计算结果进行比较。数值模拟结果表明,经过一定迁移距离后纵向弥散度与随机分布渗透系数对数场的方差和相关长度具有良好的线性统计关系,与Gelhar理论公式表达的关系类型类似。但对于较大的方差,纵向弥散度模拟结果明显大于Gelhar理论计算值,而对于较大相关长度在迁移距离不很大时,纵向弥散度模拟结果明显小于Gelhar理论计算值。本研究可为野外大尺度地下水污染预报模型中水动力弥散参数的确定提供方法借鉴。     

The effects of a gas well blow out on groundwater chemistry

Water wells were sampled near North Madison, Ohio, following a gas well blow out that injected large amounts of CH₄ into near-surface groundwater Chemical analyses showed elevated levels of Fe⁺², Mn⁺², Ca⁺², sulfide, alkalinity, and pH, and low levels of dissolved oxygen, SO₄ ⁻², and NO₃ ⁻ in CH₄-affected wells compared to unaffected wells. Sulfate reduction is quantitatively the most important vehicle for CH₄ oxidation Equilibrium thermodynamic computer models were used to simulate groundwaters from the North Madison area Model results showed that CH₄ is oxidized to HCO₃ ⁻, SO₄ ⁻² is reduced, iron and manganese oxides are reduced and dissolved, and pH increases These simulations are in excellent agreement with trends observed in the field data A laboratory experiment was designed to simulate CH₄ ⁻ perturbed groundwater in the methane-perturbed system, sulfide increased significantly, providing direct evidence for methane oxidation by sulfate reduction Although suitable anaerobic methane-oxidizing bacteria have not been isolated from groundwater aquifers, the combination of field data, laboratory experiment, and computer simulation form a convincing argument that CH₄ perturbation of aquifers can and does affect groundwater chemistry     

地下水支撑生态系统研究综述

地下水支撑生态系统（GDEs）是指部分或全部依赖于地下水的生态系统,开展GDEs研究,不仅对于了解GDEs的生态水文特征和过程具有重要学术价值,而且对于保护地下水生态环境、促进地下水可持续利用具有重要的现实意义。通过GDEs的自然和人类活动影响因素、GDEs的类型划分与识别方法、GDEs的生态水文过程模拟等方面的综述,探讨了当前存在的问题,并对未来研究趋势进行了展望。综述认为,GDEs主要受气候变化和地下水利用的影响;GDEs的类型划分要方便野外调查和管理;GDEs的识别,要综合利用水文地质调查、"3S"技术和地下水动物取样分析的方法。建议今后GDEs生态水文过程的模拟要侧重于地表水和地下水的相互作用以及地下水对植被、溶质输移和热交换等方面的贡献大小的研究。     

ICT-Based Remediation with Knowledge Information Management for Contaminated Groundwater

ICT-based remediation with knowledge information management is presented for the pump-and-treat method of contaminated groundwater. The usefulness of ICT is discussed for monitoring contaminants and groundwater level, transferring data between the remediation well and the remote remediation center, and decision support analysis for controlling the remediation well. A prototype system was developed and applied to field measurement. The prototype system performed reliably for ~600 days. As a decision support analysis, a fuzzy inference model is discussed. The membership functions were determined based on simple reliability theory. The effectiveness of the proposed method was assessed by numerical simulations. The simulation results suggest that the proposed method is likely to reduce the pumped quantity compared to PID control or an engineer’s empirical knowledge. Analysis results are also shown for cancer risks from contaminants and ground settlement risks due to excess pumping up of groundwater.     

Simulation of groundwater flow and contaminant transport at a landfill site using models

Numerical and analytical modelling studies were conducted for the analysis of groundwater flow and contaminant transport at the Innisfil landfill site in the Town of Innisfil, County of Simcoe, in Ontario, Canada. Previously conducted field studies categorized the upper stratigraphy at the site into three units: upper sand unit, upper silt/clay unit and Intermediate Sand unit. Essentially horizontal groundwater movement in the two sand units and vertical downward flow in the silt/clay unit were reported by the field hydrogeologists. In the following, application of three computer models (FLOWPATH, USGS MOC and POLLUTE) for the simulation of the groundwater flow and contaminant transport processes at the Innisfil landfill site is described. The paper focuses on the calibration of groundwater flow and contaminant transport systems, and demonstrates how the insight gained during the contaminant transport calibration was used to improve the initial groundwater flow characterization of the hydrogeological system.     

基于随机理论的地下水环境风险评价

针对确定性模型难以描述含水层非均质空间分布的问题,提出基于随机理论的地下水环境风险评价方法。以矩形场地地下水污染风险评价为例,采用蒙特卡罗法生成大量渗透系数随机场,模拟含水层参数各种可能的非均质空间分布,在此基础上建立场地地下水流模型与溶质运移模型,分别计算污染物在地下水中的迁移转化情况。统计大量随机模拟中污染事故发生的频率,当模拟次数足够多时,污染频率收敛于污染概率,污染风险即通过污染概率体现出来。该方法将模型参数设为满足一定分布特征的随机变量,避免了确定性方法得出的武断的评价结果,可为工厂的选址、水源地的选址等工作提供科学指导。     

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

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

利用趋势化随机参数场的地下水流数值模拟优化方法

水文地质参数场的刻画是建立地下水流数值模拟模型的关键问题和难点问题。通常来讲,参数场合理性程度越高,模型拟合精度越高。本次研究将随机方法和参数空间分布表达进行结合,提出了趋势化随机参数场的构建方法。以渗透系数为研究对象,首先利用MCMC采样和样本数据特征确定水文地质参数的基本数据结构,进而根据样本空间分布特征对其进行趋势化处理,最终形成趋势化的渗透系数场。通过算例分析,利用趋势化处理后的渗透系数场能够大幅提高模拟精度,相比传统赋均值方法其误差可降至原来的1/3。在北京大兴跌隆起地区进行的实例应用说明,趋势化渗透系数场对提升岩性粒径较大(中砂以上)地区模拟精度效果显著,案例中粗砂区域渗透系数经趋势化处理后平均拟合误差由2.76 m下降至0.64 m;而对岩性以细砂及以下粒径为主的区域模拟精度提升并不明显。总体来说,该方法可为地下水流数值模型的优化提供借鉴,提升模型拟合精度,从而更加合理地刻画地下水流系统。     

常州市滆湖典型地块浅层地下水数值模拟

苏锡常地区是我国经济最发达的地区之一。多年来，该地区由于超量开采深层地下水已经引起了严重的地质灾害。在当前严禁开采深层地下水的形势下，对该地区的浅层地下水资源进行科学评价，查明该地区浅层地下水的开发潜力具有十分重要的意义。本文采用GMS软件对常州市滆湖典型地块浅层地下水进行了精细的数值模拟，经模型识别和校正之后，对典型地块的浅层地下水可开采量做出评价。本次对浅层地下水的数值模拟在苏锡常地区内尚属首次，模拟结果可为今后全区合理开发利用浅层微承压地下水提供依据。     

Mass balance simulation and its application to refining flow field in Binchang area,China

Groundwater flow fields in aquifers are often determined by water level data measured in monitoring wells. The flow field can be further refined by mass balance simulations, especially when groundwater level data is limited. The mass balance simulation is based on the principle of mass conservation and relies on water quality data in the same aquifer. The approach is applied to the Luohe aquifer in the Binchang area, China. The water-rock interactions and the hydrogeochemical evolution were studied along four typical flow paths. The study indicates that groundwater in the Luohe formation flows from the southern border to the interior of the Ordos Basin. The southern border, approximately 1,400 km², is a recharge zone, where the Luohe formation outcrops. The total dissolved solids of the groundwater in the southern boarder are less than 1 g/l, and the hydrochemistry type is HCO₃–Na. This new finding refines the flow field of the water-bearing formation, and an additional 1,400 km² is included in the water resource planning of the area.     

黄河源区水文模拟中地形和融雪影响

为分析地形和融雪对水文模拟影响,利用SWAT模型在黄河源区进行1960—1990年水文模拟。二者的影响分别通过对子流域划分高程带和运用融雪模块来体现。研究表明,单独考虑融雪影响或考虑融雪、地形的共同影响下,模拟结果均能达到合格;单独考虑地形的影响下能达到良好。相比融雪影响,地形对流域水量平衡的模拟起主要作用。划分高程带对子流域的气温影响较大,对降水影响小。子流域气温的降低使得蒸散发减小,从而增加产水量,其中地下水增加最多,其次是地表径流和壤中流。地形和融雪的影响使地下水的补给来源改变,调整地下水参数后,最终模拟结果达到优秀。地形和融雪在黄河源区水文模拟中的影响可以为其他高山地区水文研究提供参考。     

Detection of contaminant plumes released from landfills: numerical versus analytical solutions

Contaminant leaks released from landfills are a threat to groundwater quality. The groundwater monitoring systems installed in the vicinity of such facilities are vital. In this study the detection probability of a contaminant plume released from a landfill has been investigated by means of both a simulation and an analytical model for both homogeneous and heterogeneous aquifer conditions. Since the detection probability is a sensitive quantity, we first compare the two methods for homogeneous aquifer conditions to assess the errors that are encountered by performing simulations. The analysis shows that the simulation model yields the detection probabilities of a contaminant plume at a given monitoring well quite well in the homogeneous case. For heterogeneous aquifers we apply the approximated analytical model based on macro-dispersivities. Here we find that this model is insufficient in monitoring system design, since the obtained analytical values of the detection probabilities at a given monitoring well differ significantly from those computed by simulation.     

Assessment of alluvial aquifer heterogeneity and development of stochastic hydrofacies models for the Hat Yai Basin in Southern Thailand

In previous studies, the groundwater flow models formulated for the Hat Yai Basin were conventional and deterministic because the geologic heterogeneity of the alluvial aquifer system in the basin had not yet been assessed. This paper describes an effort to develop hydrofacies models, such that the spatial variability of the aquifer system can be represented in a systematic way. Variogram parameters that characterize the alluvial aquifer heterogeneity were determined. Based on these variogram parameters, an indicator-based geostatistical approach was used to develop hydrofacies models using sequential indicator simulation. The hydrofacies models indicate three distinct aquifer units, namely Hat Yai, Khu Tao, and Kho Hong aquifers, which is in good agreement with a conceptual model, and incorporates spatial variability as observed in field data from borehole logs. The hydrofacies models can be used in groundwater modeling and simulations.     

水流冲刷过程中的边坡临界滑动场及河岸崩塌问题研究

针对河岸崩塌问题分析和研究,在考虑江河水位升降引起坡外水压力变化及坡内非稳定渗流基础上,同时考虑水流冲刷引起的河床冲深及河岸后退,提出了水流冲刷过程中的边坡临界滑动场和适用于天然江河崩岸的数值模拟,并对水流冲刷过程中的崩岸问题进行了分析。通过对两类不同土质岸坡的崩岸数值模拟,分析了水流冲刷引起的河床冲深及河岸后退过程中坡体的稳定性变化,探讨了不同土质岸坡的崩岸类型及崩塌模式。结果表明,坡度较陡的黏性岸坡崩塌时趋近于平面破坏且通过坡脚;坡度较缓的粉土岸坡崩塌时沿曲面破坏,且在水位骤降过程中易发生局部崩塌。     

松散承压含水层水文地质参数分区及水流场数值模拟

为了揭示华北型煤田松散承压含水层水文地质参数及其对地下水流数值模拟的意义,以安徽淮北煤田宿南矿区祁东煤矿松散层承压第四含水层(简称“四含”)为研究示范,对多个影响因素综合分析,采用层次分析-模糊综合评价法,对研究区四含水文地质参数进行分区,合理确定各分区的水文地质参数,并用于采煤情景下地下水流动态数值模拟。数值模拟结果表明：祁东煤矿2008-2012年浅部煤层开采过程中四含地下水流场没有明显改变,但南北分区水头差逐年增大,地下水流向始终从南向北,而且在井田北部水力梯度逐渐减小,在井田南部水力梯度逐渐增大,南北区以倾角为60°~70°、断距为10~320 m的魏庙断层为分界线,基岩中大倾角与大落差断层的出现是引起上覆松散承压含水层水力梯度异常的根本原因。     

北京房山应急水源地地下水数值模拟及河流渗漏量预测

以北京房山岩溶水应急水源地为例,采用等效连续介质体概化地下水流系统,建立地下水流数值模型,评价地下水源汇项,分析地下水系统特征。运用MODFLOW中的河流（渠）子程序包,刻画地表河流和泉水与地下水的交换关系。经参数识别和验证,计算结果与实测数据拟合较好,能够比较准确的模拟水源地岩溶水系统。模拟和预测结果表明,水源地地区地下水与地表水联系密切,应急开采后,地下水系统由向河流排泄转变为河流渗漏补给地下水,可为应急水源地合理开发利用提供科学依据。     

综合物探方法在应急抗旱找水中的应用

通过结合区域地质资料,采用重力、高密度电法、激电测深及音频大地电磁测深等多种地球物理勘探技术,计算机数值模拟反演计算等工作手段,为河南省鹤壁市气象干旱缺水区开展应急抗旱找水。多种物探方法的运用,对于应急抗旱找水打井的井位合理的布置、快速定井,保证成井后的出水率和地下水的合理利用提供了重要的依据,为综合物探方法的应用和类似应急找水工作提供了一定的借鉴。     

区域地下水动态数值模拟

在系统综述区域地下水数值模拟的基础上,运用有限差分和有限体积法,结合济南泉域裂隙岩溶地下水含水介质系统,采用二维非稳定流数学模型进行了数值模拟,对泉域地下水动态变化进行了评价,为泉域地下水资源评价与管理提供了一种有效的分析手段,为区域地下水资源的科学管理提供了有意义的参考。     

新疆玛河安集海灌区给水度变化规律及试验研究

在对新疆玛纳斯河流域安集海灌区进行地表水地下水转换机理和联合调度研究及利用各类集总式、分布式水文模型进行地下水数值模拟时,给水度都是一个十分重要的输入参数。分析研究区的基本地质情况,选取比较有代表性的地块,利用给水度参数测定装置分别测定常值给水度、变值给水度,揭示该区域给水度随地下水埋深的变化规律。试验得出：在计算地下水降雨入渗补给量和灌溉入渗补给量时要用到给水度这一参数,特别是在地下水埋深不大的地区,对计算成果有比较明显的影响;在地表排水沟设计中,如果考虑了给水度随埋深而变的情况,排水沟的尺寸就会发生变化。