西南岩溶山区某石化类场地地下水污染风险评价

本论文以西南岩溶山区某石化类场地为例,以污染场地风险评价导则作为指导,结合国内外地下水污染风险评价方法,在完成场地地下水污染风险识别以及源项分析的基础上,建立了符合场区实际水文地质条件的地下水水流及溶质运移数值模型,对不同情形下污染物泄露对场区浅层孔隙-溶蚀裂隙地下水的影响进行预测,结合预测结果建立了场区地下水污染风险...     

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

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

郑州地下水均衡试验场的改建工程——试验场监测资料的推广应用核心问题

郑州地下水均衡试验场的改建,是中国地质环境监测院主持的国家地下水监测工程中地下水均衡试验场的改建项目之一,由河南省地质环境监测院具体负责实施。针对场地尺度均衡试验资料推广至野外生产实践核心问题,采用确定型模型(水分运移方程)和随机模型(区域性包气带参数)相结合的研究方法。拟选黄河南冲洪积区24 km~2作为野外试验场地,通过场地上有限数量的采样,标定该场地包气带参数变异性特征,并借助于区域气象资料,结合地下水均衡试验场所建土壤水分运移数学模型成果,计算区域降雨入渗补给量与潜水蒸发,为地下水资源评价、地质环境评价及地下水资源可持续开发利用提供科学依据。同时也为大范围土壤面状污染评价工作奠定基础。     

基于FEFLOW的三维土壤-地下水耦合铬污染数值模拟研究

`土壤-地下水耦合数值模拟是定量刻画水流和溶质运移的主要手段。现有大范围场地尺度的研究受到数据采集难度及模拟计算量的限制,多是将土壤和地下水分成两个系统,这种方式不利于模型之间的计算反馈,易出现计算误差,因此将土壤和地下水作为整体系统研究具有重要意义。为精确刻画实际场地土壤-地下水系统中污染物迁移规律,揭示变饱和反应溶质迁移模型的参数敏感性,以某铬污染场地为研究对象,基于现场试验及前人研究所获数据,采用Galerkin有限元法建立三维土壤-地下水模型,定量描述六价铬在土壤-地下水中的迁移规律。在此基础上,通过改变补给条件,研究潜水面在土壤-地下水系统中的波动。并讨论阻滞系数和反应常数对溶质运移的影响。结果表明:在土壤中,污染物最大水平迁移距离为场地东南侧300 m;地下水中污染晕最大分布面积约为1.632 km2;垂向上土壤中的六价铬仅需15.6 h即可下渗至潜水面,第6天贯穿含水层。当潜水面随着补给量变化而波动时,地下水中六价铬会随水流进入土壤,影响土壤中污染分布。对溶质运移参数的讨论显示,当反应常数由0增大至10?6 s-1时,迁移出场区边界时地下水中污染物浓度约减少2000 mg/L,较难迁移至涟水河。基于FEFLOW的数值模型,能够解决各系统之间交互性差的问题,提供较为精确的模拟结果。     

锑污染土壤环境风险评价研究

综合《场地环境评价导则》与美国超级基金法案的风险评价导则中环境风险评价的方法和要求,以武汉某化工厂污染场地为例,探讨土壤中锑的环境风险评价方法;并从人体健康评价与生态风险评价出发,就经口摄入、皮肤接触、呼吸吸入锑污染物对成人和儿童进行健康风险定量计算,并评价锑污染物的生态风险。结果表明：该废弃场地内污染土壤虽然不属于危险废物,土壤中有机物和大部分重金属含量均处于相应标准限制以下,但部分区域锑含量偏高,最高质量分数达到1 262.7×10-6,66%的土壤样品锑质量分数大于推荐修复值26×10-6;污染土壤的非致癌风险非常大;污染土壤对生物具有一定的生态毒理学影响,场地土壤须进行修复后才能使用。     

甲苯在土壤地下水中迁移规律研究

石油泄漏导致甲苯等有机污染物污染土壤地下水,并且甲苯等石油类有机污染物对土壤地下水污染具有隐蔽性强、难以逆转等特点,污染治理难度大。因此,研究甲苯有机污染物在土壤地下水中运移规律是治理的关键。该文以某污染场地为研究区,以甲苯有机污染物为研究对象,考虑了淋滤、吸附的影响,采用模型计算方法,分析了甲苯在污染场地土壤地下水中迁移规律。结果表明：甲苯在土壤地下水中迁移过程中,淋滤作用对甲苯浓度变化影响显著;吸附影响下扩散作用及地下水位变化对甲苯迁移产生的影响均受到抑制。该文获得的甲苯在土壤地下水中迁移规律可为甲苯有机污染物污染场地修复提供依据。     

弥散试验在罗河铁矿变更建设项目地下水水质污染范围预测中的应用

罗河铁矿变更建设工程在项目建设、生产运行和服务期满后的各个过程中,可能造成地下水水质污染。野外弥散试验用于研究污染物在地下水中运移时其浓度的变化规律,并通过试验获得进行地下水质量定量评价的弥散参数,从而计算出特定时间内地下水水质污染的范围。     

污染场地土壤风险基准值构建与评价方法研究

基于对国内外土壤风险基准值和风险评价模型的调查，采用多介质暴露评价模型（MMSOILS），分析不同国家和地区土壤风险基准值的差异性和污染场地风险评价的关键要素。选择砷和苯作为评价的目标污染物，利用最不利的场地条件，进行风险值的量化评价，分析风险值随地下水使用点到污染场地距离的远近而发生的变化。利用线性回归，建立土壤风险基准值与健康风险值的量化关系，确立土壤浓度值14．74mg／kg和1．00mg／kg作为砷和苯的土壤风险基准值。     

西宁市贵德县河滨公园林场第四系含水层弥散试验研究

地下水污染问题日益严重,研究溶质运移的弥散理论开始应用于实际问题。建立地下水溶质运移模型,对地下水中污染物的运移及发展趋势进行准确预测,是对地下水进行保护、对地下水污染进行控制的基础。而弥散参数的确定则是地下水溶质运移模型建立的关键环节之一,直接影响着模型预测结果的精度和准确性。 对西宁市贵德县地下水污染的水质运移规律进行分析,在贵德县河滨公园林场采用径向收敛流水动力弥散理论方法进行了第四系含水层现场弥散试验,计算了试验场地潜水含水层的弥散度,获得纵向弥散度(a_L)为0.843~0.998 cm,横向弥散度(a_T)经验推断值为0.17~0.20 cm,为进一步建立该地区的地下水溶质运移模型、预测地下水污染的发展趋势和评价该地区地下水环境质量提供了数据参考。     

考虑污染物扩散风险的场地地下水污染多层次风险评估方法

地下水中污染物具有迁移性和扩散性,会对下游敏感受体造成威胁,目前场地地下水风险评估主要关注人体健康风险,还未能综合考虑地下水污染的整体风险,尤其是忽视了污染物迁移引起的对下游敏感受体的风险。本研究基于“源-径-汇”模型构建了考虑污染物扩散风险的场地地下水污染风险评估的指标体系与风险评估模式。在指标体系构建方面,重点考虑场地地下水的污染源、迁移路径和敏感受体3个方面。在风险评估模式方面,根据场地不同地下水污染状态开展3个层次的风险评估。基于假想的铬污染场地开展了案例分析,设置了地下水污染状态的4种情景,利用Wexler溶质运移模型计算了地下水污染羽的时空变化,并针对地下水污染的不同层次开展了风险评估。结果表明,在场地地下水污染羽未到达场地边界的2种情形中,场地地下水的风险评分分别为4.0,6.2,分别属于低风险与中风险。在场地地下水污染羽到达或超出场地边界的2种情形中,场地地下水的风险评分分别为7.0,8.8,分别属于中风险与高风险。综合而言,本研究构建的方法能够用来对场地地下水进行系统全面的评估和对比,能够根据风险结果对污染场地进行有效的分级管控,为污染场地的风险管控提供技术支撑。     

A possible approach for Tier 2 risk assessments of polluted sites: Framework,computer spreadsheet and application

Owing to the increasing attention placed on problems concerning site pollution that environmental geotechnics deals with, an implementation of Tier2 risk assessment (forward and backward mode) is proposed and developed in a computer spreadsheet. Consistently with a Tier 2 approach, contaminant migration is described by analytical solutions of transport models using site-specific parameters. The calculations are implemented with Microsoft Excel® while the user interfaces, which manage the various worksheets, were built with Visual Basic®. The spreadsheet was validated by comparing it with other available software that implement the same model for a given migration pathway. In the present version, the computer tool is consistent with Italian guidelines for Tiers 2 risk assessment. However, the tool can be easily adapted to comply with different regulations and recommendations. Some illustrative examples of applications are given in the paper: a case study of risk assessment for contaminated site is illustrated and a sensitivity analyses of transport factors to site-specific parameters is presented.     

某农药厂场地土壤地下水污染修复指导值探讨

参考加拿大Alberta SGRG层次性分析模式,选择其中Tier 1、Tier 2两层次进行分析,结合适合我国的参数计算出了某农药厂土壤地下水中的有机污染物1,2-二氯乙烷、1,2-二氯丙烷修复值,得到了通用的Tier 1修复指导值分别为:渗坑土壤中1,2-二氯乙烷6.2μg/kg,1,2-二氯丙烷67μg/kg;场地地下水中1,2-二氯乙烷、1,2-二氯丙烷均为5μg/L。在此基础上,结合实地调查获取场地参数,利用Tier 2评价对指导值进行修正,确定了场地最终修复目标值:渗坑土壤中1,2-二氯乙烷为11μg/kg,1,2-二氯丙烷为45μg/kg;地下水中1,2-二氯乙烷、1,2-二氯丙烷仍为5μg/L。各途径修复值的改变较符合层次性评价的理论,即Tier 2修复值较Tier 1宽松;但1,2-二氯丙烷的保护DUA途径修复值较Tier 1更加严格,这是由于场地特定条件造成的。     

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.     

Bayesian hierarchical models for soil CO<Subscript>2</Subscript> flux and leak detection at geologic sequestration sites

Proper characterizations of background soil CO₂ respiration rates are critical for interpreting CO₂ leakage monitoring results at geologic sequestration sites. In this paper, a method is developed for determining temperature-dependent critical values of soil CO₂ flux for preliminary leak detection inference. The method is illustrated using surface CO₂ flux measurements obtained from the AmeriFlux network fit with alternative models for the soil CO₂ flux versus soil temperature relationship. The models are fit first to determine pooled parameter estimates across the sites, then using a Bayesian hierarchical method to obtain both global and site-specific parameter estimates. Model comparisons are made using the deviance information criterion (DIC), which considers both goodness of fit and model complexity. The hierarchical models consistently outperform the corresponding pooled models, demonstrating the need for site-specific data and estimates when determining relationships for background soil respiration. A hierarchical model that relates the square root of the CO₂ flux to a quadratic function of soil temperature is found to provide the best fit for the AmeriFlux sites among the models tested. This model also yields effective prediction intervals, consistent with the upper envelope of the flux data across the modeled sites and temperature ranges. Calculation of upper prediction intervals using the proposed method can provide a basis for setting critical values in CO₂ leak detection monitoring at sequestration sites.     

A GIS raster technique to optimise contaminated soil removal

A geographic information system (GIS) raster technique has been developed and used interactively with remediation designers to evaluate the optimum extent of excavating soil contaminated by chlorinated solvents. The technique and the results of its application are presented. The site was a former chemical storage plant for acids and solvents. Two distinct solvent plumes were detected within the ground using a photo-ionisation detector. The solvents were found to be dissolved in the groundwater and migrating in the general direction of groundwater flow. A remediation strategy was proposed involving the localised excavation of contamination ‘hot spots’ followed by the implementation of a groundwater remediation system. A number of excavation options were discussed and the GIS raster technique was developed to evaluate these options in terms of contaminant removed and excavation cost.

The plumes were initially mapped using a triangular irregular network (TIN). These TIN models were rasterised to produce a regular grid of rectangular cells, each cell having a value relating to the concentration of contaminant at that spatial point. The proposed excavation zones were then overlaid on to the raster models as masks. The relationship between the value of contaminant concentration of cells within the mask (or excavation zone) and the total value of contaminant concentration of cells within the solvent plume was used to determine the efficiency of the excavation.

The excavation options were compared taking into account the percentage of the contaminant plume removed, the excavation area (soil volumes) and related costs. Once the GIS raster technique had been developed, it proved very quick to rerun the analysis for the other excavation zones. The optimum excavation zone, based upon cost and contaminant recovery, was found for the site. The technique helped by targeting the worst area of contamination and provided the client with a cost-benefit analysis of the different remediation options.     

Sensitivity of advective travel time of contaminants to correlated formation porosities

An approach is presented to quantify sensitivity of advective contaminant travel time to porosities of hydrogeologic units (HGUs) along the flowpaths when the porosities of different HGUs are correlated. The approach is an extension of the previous sensitivity analysis technique for independent input porosity cases. It is applicable in situations where a calibrated groundwater flow model exists, but a full contaminant transport model is not available. Three sensitivity indices are introduced based on the decomposition of covariance between the advective contaminant travel time and individual input porosities of HGUs. When the input HGU porosities are correlated, the three sensitivity indices quantify the total, intrinsic and correlated contributions from each individual HGU porosity, which should be considered in order to determine the relative importance to the uncertainty in advective contaminant travel time of the input HGU porosities contributing either independently or correlatedly. Two simple one-dimensional flow examples are presented to illustrate the applications of the approach to scenarios when each HGU has distinct porosity and situations of spatially variable porosity field. The approach is applicable to more complex multi-dimensional cases where advective contaminant travel time can be calculated based on simulated flow results from groundwater flow models.     

Estimating greenhouse gas emissions from open-cut coal mining: application to the Sydney Basin

A new site-specific (Tier 3) method has been developed to determine greenhouse gas emissions from open coal mining. The Tier 3 method presented here is based on extensive measurement of gas emissions from open-cut coal mines and the physics of gas desorption from coal. It was adopted by Australian National Greenhouse and Energy Reporting in 2009 and since 2012 formed the scientific basis for the Australian Government guidelines on calculating greenhouse gas emissions from open cut mines. The main strength of this method is its site-specific nature and accuracy, as well as its ability to be integrated with routine coal exploration programs. New concepts were produced for the model: a coal mine is regarded as a ‘gas reservoir,’ with coal seam gas being emitted from a ‘gas release zone’ that consists of sedimentary geological units (emission layers) above and below the base of the mine. The primary data required for the method are the in situ gas content and gas composition of the coal and carbonaceous rocks contained within the gas-release zone. These data are obtained through direct measurement of gas desorption from bore cores. To reduce gas drilling, a mine lease is compartmentalised into ‘gas zones’ of similar gas content and reservoir properties. The outputs of the method are emission density (the potential volume of gas emitted from mining site per unit area of the ground surface) and emission factor (the gas volume emitted per tonne of raw coal extracted). Owing to spatial variability and errors of measurement, the estimate of emissions is associated with uncertainty. A simple method of calculating uncertainty of emissions is presented in this work.     

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.     

Probability of plume interception using conditional simulation of hydraulic head and inverse modeling

A procedure to estimate the probability of intercepting a contaminant groundwater plume for monitoring network design has been developed and demonstrated. The objective of the procedure is to use all available information in a method that accounts for the heterogeneity of the aquifer and the paucity of data. The major components of the procedure are geostatistical conditional simulation and parameter estimation that are used sequentially to generate flow paths from a suspected contaminant source location to a designated monitoring transect. From the flow paths, a histogram is constructed that represents the spatial probability distribution of plume centerlines. With an independent estimate of the plume width, a relationship between the total cost and the probability of detecting a plume can be made. The method uses geostatistical information from hydraulic head measurements and is conditioned by the data and the physics of groundwater flow. This procedure was developed specifically for the design of monitoring systems at sites where very few, if any, hydraulic conductivity data are available.