中试尺度下可渗透反应墙位置优化模拟——以铬污染地下水场地为例

可渗透反应墙（PRB）是一种高效的地下水污染原位修复技术。不同水文地质条件下,污染场地墙体位置布设合理性影响其修复效果,而利用地下水数值模拟可实现墙体位置优化。文章以某Cr6+污染地下水场地为例,基于Visual Modflow建立了研究区平面二维稳定流数值模型,并通过模型检验。根据墙体的设计尺寸（长20 m×宽2 m×深12 m）及填充材料的渗透系数（80 m/d）,利用所建模型分别计算了4种布设方案（墙体尺寸大小和填充材料渗透系数相同,布设位置不同）下墙体的捕获区宽度、粒子滞留时间和通过墙体的Cr6+通量。结果表明:4种布设方案模拟的滞留时间和捕获区宽度取值差异性不大,变异系数小于2%;Cr6+通量差别较大,变异系数高达76.32%,主要由地下水中Cr6+浓度空间分布不均引起。对比分析4种方案的各评价指标,方案2求得的捕获区宽度为21.9 m,粒子滞留时间为4.1 d,Cr6+去除量可达127.7 mg/d,可作为最佳布设方案。本研究建立的地下水流数值模型符合场地实际情况,可有效评估PRB截获污染羽的范围和去除目标污染物的能力,为铬渣类污染场地PRB原位修复工程设计与实施提供技术支撑和参考依据。     

Applying numerical modeling for designing strategies of effective groundwater remediation

Selection of effective groundwater remediation scenarios is a complex issue that requires understanding of contaminants’ transport processes. The effectiveness of cleanup measures may be verified by fate and transport numerical modeling. The goal of this work was to present the usefulness of fate and transport modeling for planning, verification and fulfillment of effective groundwater remediation methods. Selection methodology was developed, which is based on results of numerical flow and transport modeling. A field site located in south-east Poland was selected as a case study, in which groundwater contamination of trichloroethene and tetrachloroethene was detected. The results indicated that “pump and treat” was the most effective among the studied remediation methods, followed by permeable reactive barrier and in situ chemical oxidation. Natural attenuation-based remediation was demonstrated to be the least suitable, as it requires the longest time to reach predefined remediation goals, principally due to low sorption capacity and unfavorable hydrogeochemical conditions for biodegradation. Fate and transport numerical modeling allowed simulating different remediation strategies, and thus the decision-making process was facilitated.     

基于多相流数值模拟的某石油污染场地地下水中VOCs自然衰减过程识别及能力评估

监测自然衰减(monitoring natural attenuation,MNA)技术是目前普遍认可的去除地下水中挥发性有机污染物(volatile organic compounds,VOCs)的技术。但受其修复周期长、监测费用昂贵等因素的影响,实地开展MNA技术修复污染场地具有一定的局限性。基于此,本研究运用多相流数值模拟手段识别了某石油污染场地内典型VOCs污染物(苯、甲苯、萘)在地下水中的自然衰减过程并评估了其自然衰减能力。结果表明:采用TMVOC所建立的多相流数值模拟模型能较好地预测和识别VOCs在地下水中的衰减规律;在研究区中,苯、甲苯和萘由于理化性质差异,在地下水中的污染羽分布特征不同,其自然衰减过程受挥发、吸附和生物降解作用的影响程度也不同;挥发和生物降解作用对VOCs自然衰减的影响程度均为苯>甲苯>萘,而吸附作用对VOCs自然衰减的影响程度为萘>甲苯>苯;在污染源被阻断的前提下,苯、甲苯和萘分别在泄漏发生7.0、6.5和6.0年后通过自然衰减达到理想去除效果。本文研究成果可以为水文地质条件类似的VOCs污染场地MNA修复方案的制定和修复效果评估提供理论支撑。     

Groundwater contaminant source identification based on iterative local update ensemble smoother

Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation, and is called groundwater contaminant source identifcation. This is a kind of typical groundwater inverse problem, and the solution is usually ill-posed. Especially considering the spatial variability of hydraulic conductivity field, the identification process is more challenging. In this paper, the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model (MT3DMS) and a data assimilation method (Iterative local update ensemble smoother, ILUES). In addition, Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction. In practical problems, the geostatistical method is usually used to characterize the hydraulic conductivity feld, and only the contaminant source information is inversely calculated in the identifcation process. In this study, the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field. The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site, and ILUES has good performance in solving high-dimensional parameter inversion problems.     

地下水污染场地风险管理与修复技术筛选

国际上对于地下水污染场地的控制与修复研究已经取得了许多成果,已有成功的修复实例。我国虽然起步晚,但非常重视地下水污染的防治,开展了全国范围的地下水污染调查,并进行了地下水污染的防治规划。地下水污染的控制与修复已经逐渐进入示范性研究阶段。面对地下水污染场地风险管理的不同方法,以及众多的污染修复技术,如何制定风险管理策略,如何在各种各样的修复技术中筛选合适的技术或技术组合,对于地下水污染场地的防治具有非常重要的意义。笔者分析了发达国家地下水污染风险管理策略,结合在地下水污染场地研究方面的经验,对一些主要的地下水污染修复技术进行了分析论述,提出了考虑污染物特征、场地水文地质条件的地下水污染修复技术的筛选过程和方法。     

Assessment of copper,cobalt and zinc contaminations in soils and plants of industrial area in Esfahan city (in Iran)

This paper discusses the artificial groundwater recharge effect of high-infiltration basins. For this purpose, the hydrogeological parameters of the study area are collected to construct a conceptualized physical model. The TOUGH2 numerical simulation software is then used to simulate the infiltration behavior of an artificial recharge into an underground aquifer. Four wells (MW-1, MW-2, MW-3, and MW-4) are observed at the field site, after which the groundwater levels are compared with the simulation results. It is found that good agreement exists between the observed and numerical data for MW-1 and MW-2. However, the observed groundwater level in MW-3 is higher than the simulated level. We also find that MW-3 is at the edge of the artificial recharge lake, and that the high groundwater level may well be the result of a portion of the infiltration load following the well border into the well screen. Conversely, the groundwater level in MW-4 is found to be lower than in the simulated well due to local permeability in the well location. Finally, the numerical results predict that the groundwater level will attain a steady state at approximately 47 h after the onset of infiltration.     

天然冷源对地下水源热泵的影响规律

东北某些地区长期存在着一定规模的地下冻土,在地下水源热泵的运行中可视为天然冷源。天然冷源对地下水源热泵的利用产生一定的影响。为了提高地下水源热泵运行效率,运用室内外实验和数值模拟方法探寻冷源对地下水源热泵影响规律。对加格达奇地区的冷源场地进行了室内热物性试验和现场岩土热响应试验。研究表明：该地区砂、砾的热物性参数（比热容、导热系数）随深度增加呈现离散现象,砂岩、花岗岩热物性参数变化不大;岩土体导热系数花岗岩最大,砾、砂最小,砂岩居中;比热容受深度影响波动较大,特别是不同深度的砾、砂层比热容差别较大。依据该场地的地质条件,构建数值模型,开展模拟研究,对多井抽灌系统下热突破现象以及影响范围进行模拟分析,模拟结果显示冷源对该场地地下水源热泵的影响范围在150 m以内,而且在小于75 m的范围内对地下水源热泵的影响最为显著,为未来该地区地下水源热泵的设计、选址提供理论基础和参考依据。     

中国地下水污染修复方法和技术应用展望

[研究目的]:中国地下水污染调查和修复日益受到科学界的重视,了解和掌握地下水污染修复方法和技术有助于对污染场地进行科学修复.[研究方法]:本文在系统分析国内外地下水污染修复案例的基础上,对中国地下水污染修复现场实施的技术方法进行总结.[研究结果]:结合中国区域经济发展特征和地下水污染调查评价成果认为,复杂的水文地质条件...     

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.     

基于水化学和同位素特征的新乡某地下水污染场地水文地质概念模型细化

构建合理的水文地质概念模型对于地下水数值模拟至关重要,概念模型如果确定不合理,整个数值模拟将毫无意义。通常水文地质调查可以大致确定水文地质概念模型,然而对于一些不易被发现的补给源,通过水文地质调查也无法确定。水化学和同位素特征对地下水的补给源及水流路径有很好的指示作用,可以根据这些信息对水文地质概念模型进行细化。文章通过对研究区地下水取样,分析地下水水化学和同位素特征,确定出研究区内不易识别的地下水补给源及顶层黏土-粉质黏土的透水性,最终细化了水文地质概念模型边界条件。这对于后期得到正确的地下水流和溶质运移模拟结果具有重要意义。     

Delineation of mafic intrusions near Bedford (Virginia,USA) using geological and geophysical methods

In the present work, geological and geophysical methods were used to delineate the locations of multiple mafic intrusions at the Claytor Nature Study Center (CNSC) near Bedford, VA. An investigation of the groundwater hydrology of CNSC was launched in 2007. As a first step in that project a preliminary geological survey revealed sparse evidence of a number of diabase intrusions in the area. While diabase intrusions are not particularly permeable features, contact metamorphism of the host rock could provide conduits for groundwater due to stress fractures and joints and high-temperature recrystallization of the rock matrix. Following the geological survey, geophysical surveys including seismic refraction, ground penetrating radar, and magnetic ground measurements were conducted to determine the location and extent of these multiple igneous intrusions. Seismic and radar surveys proved to be inconclusive, but the magnetic surveys showed strong magnetic anomalies. The magnetic data were obtained using a Geometrics G-856 proton precession magnetometer and were interpreted using the Mag2dc algorithm and SGeMS geostatistical software. The results show that the intrusions are dikes that cut across nearly perpendicular to the regional metamorphic structures and trend generally north–south with a dip of approximately 75°–90° to the west. These findings are consistent with one of the general directions of stresses associated with the North Atlantic seafloor spreading in late Triassic or early Jurassic period. Subsequent hydrologic testing and groundwater modeling confirm the role of the dikes in the overall hydrogeology of the site.     

Evaluation of remedial alternatives for a petroleum contaminated unconfined aquifer with fluctuating groundwater level

This study evaluated remedial alternatives for a petroleum-contaminated site where an unconfined aquifer composed of a sandy layer of about 3–3.5 m thickness is covered by alluvial deposits and reclaimed soil of about 1.5 m thickness. Precambrian gneiss, of low permeability, lies below the sandy layer. The shallow water table is about 3 m below the surface, but shows high fluctuations of up to 1.5 m in response to precipitation events. The unsaturated soil near the water table and the groundwater are highly contaminated with petroleum hydrocarbons, especially toluene, which have apparently leaked from storage tanks. Selection of the remedial alternatives required consideration of the relevant laws on soil and groundwater conservation in Korea, the results of risk analysis and the hydrogeological conditions. The contaminated area has been divided into zones in which different remediation goals are set based on risk analysis and the degree of natural attenuation. It is estimated that the clean-up goal can be achieved in two years by the combined use of a trench drain and well point pumping to collect the contaminated groundwater for treatment, and a dual air injection system for the contaminated soil.     

Pumping-rate control for contaminated groundwater with VOCs by using fuzzy inference model

In this paper, a new method is proposed for the remediation of contaminated groundwater. The method is based on fuzzy inference and risk evaluation. The effectiveness of the presented method is assessed in numerical simulations. The main results of the study are as follows: 1) A pump-and-treat control algorithm was proposed by using the fuzzy inference and the human health risk calculation model, 2) It was found that applying the proposed algorithm was likely to reduce the pumped quantity, and 3) It was found that the proposed model can be used to notify residents of the human health risk; the model is based on the ASTM RBCA model for residents. A new fuzzy control system for contaminated groundwater can be used as a useful model for characterising the effects of contaminants on human health and providing helpful information on the human risk assessment of the contaminated groundwater site.     

基于数值模拟的岩溶地下水源保护区划分技术研究

以山东省邹城市某岩溶水源地为例,研究数值模拟方法在岩溶地下水水源地保护区划分中的应用,以及水文地质参数对于水源地保护区划分结果的影响。基于GMS软件构建了研究区三维非稳定流地下水数值模型,采用质点追踪技术,通过反向追踪示踪粒子在运移100 d、1 000 d以及25 a后的位置及迁移轨迹,对该岩溶水源地进行保护区划分;改变数值模型中的主要水文地质参数,以变化后的保护区面积作为衡量参数敏感度的指标,并计算各参数的敏感度系数,研究水文地质参数对于保护区划分范围的影响。结果表明:当参数改变幅度在20%范围内时,一层垂向渗透系数(VK1)的敏感度系数最大可达到2.63×10-3,二层垂向渗透系数(VK2)的敏感度系数最大可达到3.64×10-3;垂向渗透系数的敏感度明显高于其他参数,说明垂向渗透系数对保护区划分范围的影响最大。因此,在应用数值模拟法对岩溶地下水源地进行保护区划分时,应尤其注意模型中各含水层垂向渗透系数取值的准确性和合理性。      

Small scale study of groundwater flow in a fractured carbonate-rock aquifer at the St-Eustache quarry, Québec, Canada

A small-scale hydrogeological study was conducted in a fractured carbonate-rock aquifer located in a quarry to relate groundwater flow to the fracture network. The field study in the St-Eustache quarry, which integrates structural surveys, well logging and hydraulic tests, showed that the most important features that affect groundwater flow in the sedimentary aquifer are high hydraulic conductivity horizontal bedding planes. Vertical fractures are abundant in the quarry and throughout the region, but they have a minor effect on groundwater flow. To have a significant impact on the flow regime and lead to vertical groundwater flow, the permeability of all vertical joints need to be enhanced compared to what was generally observed at the site. Such an increase in permeability could potentially occur where dissolution and fracturing is more intense or at stress release locations such as near the surface in the quarry.     

A topo-geochemical sequence study of groundwater in Asa drainage basin, Kwara State, Nigeria

 A strong geochemical gradient was observed in the thick overburden aquifer of the Asa drainage basin. Different types of groundwater occur at different (downslope) locations and groundwater table depths. The following sequence was noticed with increasing distance downslope or with increasing groundwater table depth: 1. Ca–Mg–HCO₃ water at about 390-m groundwater table elevations or upslope locations. 2. Ca–Mg–HCO₃–Cl water at middle-slope locations or groundwater table elevations of about 350 m above sea level; 3. Ca–Mg–SO₄–Cl water at downslope locations or groundwater table elevations of about 300 m above sea level. In this basin, changes in the type of water are expected at about every 40–50 m depth from the surface. Statistical analysis via the determination of the correlation coefficient (r) and regression analysis shows that about 80–99% of the variation in groundwater chemistry is accounted for by the topography, using the model presented in this paper. The rate of change in the sequence will depend on the permeability of the aquifer, which determines the rate of groundwater flow and the residence time, and the nature of recharge. Received: 4 February 1997 · Accepted: 22 July 1997     

Modeling of an MTBE plume at Pascoag, Rhode Island

A numerical groundwater flow and mass transport model was developed to predict the extent of impact from methyl tertiary butyl ether (MTBE) release on a down-gradient drinking water well field. An MTBE incident in Pascoag, Rhode Island, was used as a case study and the plume’s past and future development was simulated using scenario analysis. The numerical code used was GeoSys/Rockflow, which permits a coupled flow and transport simulation as one object, thus alleviating the need for simulating the MTBE fate with separate flow and transport codes. The numerical model was built on available hydrogeological and chemical data as well as on GIS information of the site. By comparing the simulated results with observed field data, it was found that the model could provide reliable results even when the simulated aquifer was simplified to a two-dimensional flow and transport domain. Finally, the calibrated model was used for exploring a location that may be suitable for a new well field. Despite the model limitations associated with uncertainties of data and simplifying assumptions, numerical modeling of this MTBE contaminated site proved a useful tool and provided guidance for future municipal well field operation strategies and aquifer remediation alternatives.     

Geologic uncertainty in a regulatory environment: An example from the potential Yucca Mountain nuclear waste repository site

Regulatory geologists are concerned with predicting the performance of sites proposed for waste disposal or for remediation of existing pollution problems. Geologic modeling of these sites requires large-scale expansion of knowledge obtained from very limited sampling. This expansion induces considerable uncertainty into the geologic models of rock properties that are required for modeling the predicted performance of the site.One method for assessing this uncertainty is through nonparametric geostatistical simulation. Simulation can produce a series of equiprobable models of a rock property of interest. Each model honors measured values at sampled locations, and each can be constructed to emulate both the univariate histogram and the spatial covariance structure of the measured data. Computing a performance model for a number of geologic simulations allows evaluation of the effects of geologic uncertainty. A site may be judged acceptable if the number of failures to meet a particular performance criterion produced by these computations is sufficiently low. A site that produces too many failures may be either unacceptable or simply inadequately described.The simulation approach to addressing geologic uncertainty is being applied to the potential high-level nuclear waste repository site at Yucca Mountain, Nevada, U.S.A. Preliminary geologic models of unsaturated permeability have been created that reproduce observed statistical properties reasonably well. A spread of unsaturated groundwater travel times has been computed that reflects the variability of those geologic models. Regions within the simulated models exhibiting the greatest variability among multiple runs are candidates for obtaining the greatest reduction in uncertainty through additional site characterization.     

Contamination,risk, and heterogeneity: on the effectiveness of aquifer remediation

The effectiveness of aquifer remediation is typically expressed in terms of a reduction in contaminant concentrations relative to a regulated maximum contaminant level (MCL), and is usually confirmed by sparse monitoring data and/or simple model calculations. Here, the effectiveness of remediation is re-examined from a more thorough risk-based perspective that goes beyond the traditional MCL concept. A methodology is employed to evaluate the health risk to individuals exposed to contaminated household water that is produced from groundwater. This approach explicitly accounts for differences in risk arising from variability in individual physiology and water use, the uncertainty in estimating chemical carcinogenesis for different individuals, and the uncertainties and variability in contaminant concentrations within groundwater as affected by transport through heterogeneous geologic media. A hypothetical contamination scenario is developed as a case study in a saturated, alluvial aquifer underlying an actual Superfund site. A baseline (unremediated) human exposure and health risk scenario, as induced by contaminated groundwater pumped from this site, is predicted and compared with a similar estimate based upon pump-and-treat exposure intervention. The predicted reduction in risk in the remediation scenario is not an equitable one—that is, it is not uniform to all individuals within a population and varies according to the level of uncertainty in prediction. The importance of understanding the detailed hydrogeologic connections that are established in the heterogeneous geologic regime between the contaminated source, municipal receptors, and remediation wells, and its relationship to this uncertainty is demonstrated. Using two alternative pumping rates, we develop cost-benefit curves based upon reduced exposure and risk to different individuals within the population, under the presence of uncertainty.     

Groundwater contaminant source identification based on QS-ILUES

When groundwater pollution occurs,to come up with an efficient remediation plan,it is particularly important to collect information of contaminant source(location and source strength)and hydraulic conductivity field of the site accurately and quickly.However,the information can not be obtained by direct observation,and can only be derived from limited measurement data.Data assimilation of observations such as head and concentration is often used to estimate parameters of contaminant source.As for hydraulic conductivity field,especially for complex non-Gaussian field,it can be directly estimated by geostatistics method based on limited hard data,while the accuracy is often not high.Better estimation of hydraulic conductivity can be achieved by solving inverse groundwater problem.Therefore,in this study,the multi-point geostatistics method Quick Sampling(QS)is proposed and introduced for the first time and combined with the iterative local updating ensemble smoother(ILUES)to develop a new data assimilation framework QS-ILUES.It helps to solve the contaminant source parameters and non-Gaussian hydraulic conductivity field simultaneously by assimilating hydraulic head and pollutant concentration data.While the pilot points are utilized to reduce the dimension of hydraulic conductivity field,the influence of pilot points’layout and the ensemble size of ILUES algorithm on the inverse simulation results are further explored.