含水层层状非均质对地下水流系统的影响

区域尺度上含水层非均质具有复杂的结构性和随机性,难以准确刻画,造成非均质对区域地下水流系统的影响机制研究不够深入。本文以鄂尔多斯盆地白垩系地下水流系统为研究实例,选择典型剖面,采用剖面二维随机数值模拟方法,通过对比不同非均质刻画方法下地下水流场的变化,探讨含水层层状非均质对地下水流系统的影响机制。结果显示,均质条件下模型各向异性（含水层水平和垂向渗透系数比值Kh/Kv）取值为1000时,地下水流场与实际条件较为接近;非均质条件下,渗透系数方差取值0.91,水平相关长度取值5000 m,Kh/Kv取值150时,接近实际条件。研究表明,在大尺度地下水流模拟研究中,采用水平相关长度、渗透系数方差和各向异性值三个变量生成的随机场能很好地刻画含水层的层状非均质特征及其对水流系统的影响控制作用。由于含水层不同尺度层状非均质的叠加效应,采用均质各向异性介质等效概化含水层层状非均质性会造成等效各向异性值偏大失真的效应。     

Effect of sand lenses on groundwater flow and contaminant migration

The effect of homogeneous sand lenses on the groundwater flow in an otherwise homogeneous clayey deposit is examined by performing Monte Carlo simulations using a finite element flow model. In the simulations, the locations of the sand lenses are assumed to be mutually independent. The paper examines the effect on the flow field in a clayey deposit of (a) different percentages of sand lenses, (b) different hydraulic conductivities of the sand lenses, (c) different average sand lens sizes, (d) non-uniformity of the sizes of the sand lenses, and (e) localization of the sand lenses. The effect of these non-uniform flow fields on contaminant migration is then examined using a finite element contaminant transport model. For the range of cases considered it is shown that: the volume of sand lenses present has a greater influence than the shape, size, location and hydraulic conductivity of the sand lenses; simplified calculations performed using the geometric and harmonic means of hydraulic conductivity bracket the behaviour evident from more complex analyses; and the maximum impact on an aquifer separated from a waste disposal facility by a deposit containing sand lenses can be modelled to sufficient accuracy, using quite simple flow and contaminant transport models.     

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.     

The effect of aqueous metal-chlorophenolate complexation on contaminant transport in groundwater systems

Aqueous metals and chlorophenols are common co-contaminants of groundwater systems. However, the importance of aqueous metal-chlorophenolate complexation cannot be accurately assessed because the stability constants for environmentally important aqueous metal-phenolate and metal-chlorophenolate complexes have not been measured. In order to determine the role these complexes play in contaminant transport, this study applies a correlation technique to the limited data that do exist to estimate the stabilities of metal-chlorophenolate complexes of environmental interest. Speciation calculations that are based on these estimated stability constants indicate that aqueous metal-chlorophenolate complexation may significantly affect both aqueous metal and dissolved chlorophenolate species distributions. Therefore, aqueous metal-chlorophenolate complexation may affect the extent of adsorption of both metals and chlorophenolates onto mineral surfaces. In addition, aqueous complexation may significantly enhance dissolution of aquifer aluminosilicate minerals. This study suggests that aqueous complexation between metals and chlorophenolates can significantly affect the mobility of metal and phenolic contaminants.     

The effect of aquifer heterogeneity on natural attenuation rate of BTEX

Monitored natural attenuation can be a viable option for remediation of groundwater contamination by BTEX compounds. Under the field conditions, the rate of contaminant mass attenuation through natural processes, such as biodegradation, to a large extent affected by the groundwater flow regime, which is primarily controlled by the aquifer heterogeneity. Numerical simulation techniques were used to describe quantitatively the relationship between biodegradation rate of BTEX and aquifer heterogeneity. Different levels of aquifer heterogeneity were described by random hydraulic conductivity fields (K) having different statistical parameters, the coefficient of variation (CV) and the correlation length (h). The Turning Bands Algorithm was used to generate such K fields. Visual MODFLOW/RT3D was used to simulate the fate and transport of dissolved BTEX plume within heterogeneous aquifers. The multispecies reactive transport approach described BTEX degradation using multiple terminal electron-accepting processes. First-order biodegradation rate constants were calculated from simulated BTEX plumes in heterogeneous flow fields. The results showed that aquifer heterogeneity significantly affected biodegradation rate; it decreased with increasing CV when h was in the range of up to 12 m, whereas it increased with increasing CV when h was greater than about 12 m. For well characterized aquifers, this finding could be of great value in assessing the effectiveness of natural attenuation during feasibility studies at BTEX contaminated sites.     

Adaptive state estimation of groundwater contaminant boundary input flux in a 2-dimensional aquifer

In many circumstances involving heat and mass transfer issues,it is considered impractical to measure the input flux and the resulting state distribution in the domain.Therefore,the need to develop techniques to provide solutions for such problems and estimate the inverse mass flux becomes imperative.Adaptive state estimator(ASE)is increasingly becoming a popular inverse estimation technique which resolves inverse problems by incorporating the semi-Markovian concept into a Bayesian estimation technique,thereby developing an inverse input and state estimator consisting of a bank of parallel adaptively weighted Kalman filters.The ASE is particularly designed for a system that encompasses independent unknowns and/or random switching of input and measurement biases.The present study describes the scheme to estimate the groundwater input contaminant flux and its transient distribution in a conjectural two-dimensional aquifer by means of ASE,which in particular is because of its unique ability to efficiently handle the process noise giving an estimation of keeping the relative error range within 10%in 2-dimensional problems.Numerical simulation results show that the proposed estimator presents decent estimation performance for both smoothly and abruptly varying input flux scenarios.Results also show that ASE enjoys a better estimation performance than its competitor,Recursive Least Square Estimator(RLSE)due to its larger error tolerance in greater process noise regimes.ASE's inherent deficiency of being slower than the RLSE,resulting from the complexity of algorithm,was also noticed.The chosen input scenarios are tested to calculate the effect of input area and both estimators show improved results with an increase in input flux area especially as sensors are moved closer to the assumed input location.     

The influence of conceptual model of sedimentary formation hydraulic heterogeneity on contaminant transport simulation

Development of heterogeneity model of layered sandy-clay formation and impact of this model on transport is considered. The lithological data of more than 250 wells that captured 300 meters formation at the investigated area of 40 km² are used for model of heterogeneity construction. Two models of heterogeneity were developed with using these well data: TP/MC model based on 3D Markov chain simulation for four hydrofacies and 2D kriging interpolation of thicknesses of elementary lithological layers. Simulation of conservative transport by particle tracking algorithm shows that horizontal transport along layers is similar for both models. The main difference is in vertical transport cross formation bedding. The kriging interpolation model gives more conservative results than TP/MC model due to larger characteristic horizontal length of layers in the kriging model. As the result vertical effective hydraulic conductivity of formation is in two times larger and the first particle arriving time is in four times faster in TP/MC model.     

Adaptive state estimation of groundwater contaminant boundary input flux in a 2-dimensional aquifer

In many circumstances involving heat and mass transfer issues, it is considered impractical to measure the input flux and the resulting state distribution in the domain. Therefore, the need to develop techniques to provide solutions for such problems and estimate the inverse mass flux becomes imperative. Adaptive state estimator (ASE) is increasingly becoming a popular inverse estimation technique which resolves inverse problems by incorporating the semi-Markovian concept into a Bayesian estimation technique, thereby developing an inverse input and state estimator consisting of a bank of parallel adaptively weighted Kalman filters. The ASE is particularly designed for a system that encompasses independent unknowns and /or random switching of input and measurement biases. The present study describes the scheme to estimate the groundwater input contaminant flux and its transient distribution in a conjectural two-dimensional aquifer by means of ASE, which in particular is because of its unique ability to efficiently handle the process noise giving an estimation of keeping the relative error range within 10% in 2-dimensional problems. Numerical simulation results show that the proposed estimator presents decent estimation performance for both smoothly and abruptly varying input flux scenarios. Results also show that ASE enjoys a better estimation performance than its competitor, Recursive Least Square Estimator (RLSE) due to its larger error tolerance in greater process noise regimes. ASE’s inherent deficiency of being slower than the RLSE, resulting from the complexity of algorithm, was also noticed. The chosen input scenarios are tested to calculate the effect of input area and both estimators show improved results with an increase in input flux area especially as sensors are moved closer to the assumed input location.     

Analysis on contaminant migration through vertical barrier walls in a landfill in China

Vertical barrier walls are often constructed to prevent contamination of ground water and soils by landfill leachate. The leachate water levels in landfills in southern China are generally high. Contaminants in such landfills may migrate through the vertical barrier walls and give rise to environmental problems. Qizishan landfill in Suzhou, China was taken as an example to investigate contaminant migration through the vertical barrier walls. Advection, diffusive and adsorption processes were considered in the analysis. Influences of permeability and depth of the barrier wall on contaminant migration were analyzed. The results show that it will be 13.5 years before breakthrough at 0.1% of the source concentration and 20.5 years before breakthrough at 10% of the source concentration. By and large, the contaminant has not passed through the barrier wall at present, and the contaminated zone is mainly located in the sandy clay layer near the earth dam, which is validated by testing on sampled soils. Hydraulic conductivity and depth of the barrier wall are critical to contaminant migration. Special attention need to be paid on them when building such a barrier wall. If bottom of the barrier wall is keyed into the aquitard and the hydraulic conductivity reaches 10⁻⁹ m/s, the time before breakthrough will be long enough to allow stabilization of the landfill. Pollution of the surroundings will be avoided, and therefore the requirement for contaminant control will be attained.     

Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.     

北京西山岩溶水中硝酸盐来源及迁移特征

西山岩溶水作为北京市重要的供水水源,在区域水文地质调查过程中发现西山岩溶水补给区一些样品硝酸盐浓度偏高,然而其周围岩溶水中硝酸盐浓度并未见明显变化（浓度均值为18.81 mg/L）,具点状硝酸盐源特征。岩溶含水层中硝酸盐分布特征和影响因素具有重要的研究意义。文章利用岩溶水化学和同位素分析结果研究了岩溶水硝酸盐分布和来源,利用IsoSource软件计算各来源的贡献率。结果表明:军庄—永丰屯径流区和军庄—古城—玉泉山泉径流区岩溶水中硝酸盐主要来源于生活废水、土壤有机氮矿化和大气沉降,贡献率值分别为37.1%、36.3%、26.6%;潭柘寺—四季青岩溶水中硝酸盐在潭柘寺补给区有养殖场废水下渗影响,径流区中硝酸盐含量低并发生反硝化降解作用。永定河水源硝酸盐不是岩溶水硝酸盐的主要来源。研究成果对西山岩溶水中硝酸盐的控制具有重要参考价值。     

基桩对含水层渗流阻挡作用效应研究

超高层建筑的桩基础通常会打入到深度较大的承压含水砂土层中;在承压含水层的上方通常会进行超深基坑的开挖构建超高层建筑的地下室。为了保持超深基坑的基底稳定,需要在承压含水层中进行降水施工。由于深基坑存在大量的桩基础,基桩的存在会影响降水时地下水的渗流路径,继而影响土体的渗透性状。应用室内试验研究基桩对含水层渗流阻挡作用的机制。试验采用PVC管模拟基桩,将颗粒较均匀的砂土作为试验土体,通过室内渗流试验来研究基桩对地下水渗流的影响。试验结果表明,基桩对于地下水渗流有明显的阻挡效应：基桩的数量以及布置形式会对渗透效应产生较大的影响。考虑砂土体积置换率的有效介质理论可以用来分析矩形布置的基桩对渗透系数的影响,但对梅花形布置的基桩应用有效介质理论计算所得的渗透系数会有将近20%的偏差;其原因是梅花形布置的基桩,会使渗流路径增加。基于对渗流路径的分析,提出了考虑有梅花形布置的基桩存在时等效渗透系数计算的修正方法。     

Aquifer heterogeneity: Hydrogeological and hydrochemical properties of the Botany Sands aquifer and their impact on contaminant transport

Detailed geological mapping and drilling has shown that the contact between the Cambrian volcano‐sedimentary sequence at Rosebery and the Mt Read Volcanics is formed by a major thrust fault dipping east at 40° and having a displacement of at least 1.5 km. The sedimentary sequence is part of the Dundas Group, a Middle to Late Cambrian forearc‐like sequence which unconformably overlaps the volcanics south of Rosebery. The Rosebery Thrust Fault marks the eastern boundary of a zone of folding, faulting and disruption which affects the Dundas Group and the tectonically interfingered and underlying basaltic greywacke‐mudstone sequence of the Crimson Creek Formation. At least some of this deformation occurred prior to deposition of the Ordovician Limestone, as evidenced by marked angular discordances. The complex area can be interpreted as a Cambrian accretionary prism‐forearc‐arc sequence developed above an east‐dipping subduction zone.

The Henty Fault System, which cuts obliquely through the Mt Read belt and encloses a misfit wedge of sediments, pillow lavas, gabbros and ultramafic rocks, is interpreted as a remnant of an inter‐arc basin. The fault system separates a dacitic‐andesitic arc segment to the northwest from a more rhyolitic segment to the southeast. The latter is overlain by a younger arc sequence, the Tyndall Group, which may have been the source for the Dundas Group volcanic detritus.     

苏锡常地区主采层地下水流数值模拟

本次数值模拟是从苏锡常整个大区域考虑,建立统一的主采层地下水流模拟模型,通过对水文地质概念模型的重新认识及边界条件的合理概化,较大程度地减少了人为因素引起的误差.利用已校正过的模型对该地区禁采后地下水水位恢复情况进行了预测,结果显示该模型较适用于实际.     

哈密盆地多层含水系统地下水流数值模拟

对具有多层含水结构的哈密盆地地下水系统, 提出了一类考虑弱透水层弹性释水作用的地下水流耦合数学模型。利用实际资料的数值模拟结果表明, 与忽略弱透水层弹性释水作用的数学模型相比, 其精度提高、更符合实际的水文地质条件, 为定量评价多层越流系统地下水资源提供了一种实用可靠的数值方法。     

Groundwater flow rate and contaminant migration in fissure-karstic aquifer of Opole Triassic system southwest Poland

 This paper presents hydrogeological problems occurring during municipal water exploitation and mine dewatering. These activities result in groundwater quantity and quality changes in the fissure-karstic aquifer. Increase of nitrate concentration up to 12 mg NNO₃/l due to intensive fertilizer use, and high tritium concentration, show water system impact up to 100 m depth. Intensive water exploitation produces large cones of depression with over 40-m water-level depletion in the Opole region. Flow rates of major components and isotopes have been verified by chemical migration history. Some aspects of the protection policy of this type of aquifer are also discussed. Received: 7 March 1997 · Accepted: 17 November 1998     

含水层渗透系数空间变异性对地下水数值模拟的影响

地下水数值模拟是目前定量研究地下水水量和水质的重要手段。使用基于随机理论的MonteCarlo方法来进行地下水数值模拟。这种方法能较好地考虑水文地质参数的空间变异性。主要将MonteCarlo方法和确定性模型模拟方法的模拟结果在渗透系数场、水头场、速度场和浓度场等方面进行了比较。结果表明:在模拟三维非均质含水层中的溶质运移问题时,充分考虑了含水层渗透系数空间变异性的MonteCarlo法比确定性方法更为有效,模拟精度提高了很多,且对模拟误差及误差来源有合理的数学解释。     

Characterizing vertical contaminant distribution in a thick unconfined aquifer,Hanford site,Washington, USA

Drilling-intensive aquifer characterization techniques have been used to obtain depth-discrete water samples from a thick, hydrogeologically continuous unconfined aquifer system; groundwater results indicate that carbon tetrachloride contamination is widespread and extends deeper and at concentrations much higher than detected in monitoring networks at the water table. Carbon tetrachloride, a legacy waste, was used in the plutonium extraction process at the Hanford site in south-central Washington State. Vertical, depth-discrete groundwater samples were collected during well drilling throughout a 28-km² region to determine the concentration of carbon tetrachloride present as a dissolved phase in the aquifer. Results indicate that high concentrations of carbon tetrachloride, three orders of magnitude above the allowable regulatory limit, are present at depths greater than 25 m beneath the water table. In support of future efforts to remediate the carbon tetrachloride contamination, it is imperative to locate the remaining chemical inventory, determine the vertical as well as the lateral distribution of this contaminant and its physical form. Depth-discrete aquifer characterization throughout the uppermost-unconfined aquifer system is providing this information and improving the understanding of the contaminant distribution and the hydrogeologic framework through which it moves.

Postaudit evaluation of conceptual model uncertainty for a glacial aquifer groundwater flow and contaminant transport model

Numerical groundwater flow and contaminant transport modeling incorporating three alternative conceptual models was conducted in 2005 to assess remedial actions and predict contaminant concentrations in an unconfined glacial aquifer located in Milford, Michigan, USA. Three alternative conceptual models were constructed and independently calibrated to evaluate uncertainty in the geometry of an aquitard underlying the aquifer and the extent to which infiltration from two manmade surface water bodies influenced the groundwater flow field. Contaminant transport for benzene, cis-DCE, and MTBE was modeled for a 5-year period that included a 2-year history match from July 2003 to May 2005 and predictions for a 3-year period ending in July 2008. A postaudit of model performance indicates that predictions for pumping wells, which integrated the transport signal across multiple model layers, were reliable but unable to differentiate between alternative conceptual model responses. In contrast, predictions for individual monitoring wells with limited screened intervals were less consistent, but held promise for evaluating alternative hydrogeologic models. Results of this study suggest that model conceptualization can have important practical implications for the delineation of contaminant transport pathways using monitoring wells, but may exert less influence on integrated predictions for pumping wells screened over multiple numerical model layers.