基于MAROS平台的地下水污染监测网优化设计

以东北地区某石油烃污染场地为例,在充分掌握研究区内地质、水文地质条件以及相关污染监测数据的基础上,以最优化监测决策支持系统——MAROS为平台,对研究区内现有地下水污染监测网进行优化设计。在阐述数理统计方法、Delaunay方法以及修正的CES方法原理的基础上。分别对区内单井浓度变化趋势和监测网系统污染物质量的空间矩特征、监测井的冗余性以及单井监测频率的有效性进行了分析:基于MAROS平台优化结果,并结合现有地下水污染监测网中取样点的冗余性、空间分布以及监测频率的有效性的分析,提出了该研究区的地下水污染监测网的优化设计方案并分析了其相关影响因素及不确定性因素的影响。研究结果显示,在现有地下水污染监测网中可以去除5口监测井,并根据每口井中污染物浓度的变化趋势、空间分布及其冗余性调整为一年四次、一年一次、两年一次或不监测,进一步提高区内地下水污染监测网的监测效率、降低监测费用,并对地下水污染防治和修复工作具有重要的指导意义。     

地下水污染监测网的设计研究进展

地下水污染监测网的设计包括取样点在空间上的采样位置和时间上的取样频率这两方面的确定，其目的是为了准确刻画污染羽在含 水层中随时间的变化状况。概要地回顾了近20年来地下水污染监测网设计的研究成果。分别介绍了统计方法、模拟方法和模拟-优化模型 等监测网设计方法的研究进展。监测网设计方法的选择取决于最终的监测目的以及可供利用的基础资料。指出参数的不确定性是影响污染监 测设计结果的最重要因素。如何将地下水污染监测网设计的理论研究真正与实际应用相结合是今后的主要研究方向。     

介绍一种优化地下水污染监测网的软件--MAROS

MAROS是首个基于Windows平台利用MicrosoftAccess和MicrosoftExcel开发出来的用于地下水污染监测网优化设计的决策支持系统软件。该软件通过统计分析监测区域有关数据在时空上的相关分析结果,来确定监测网在时间上的取样频率和空间上的分布密度,从而保证以最省的监测费用获得地下水污染物随时间和空间变化的最全面信息。本文简单地介绍了MAROS的设计原理和模块组成,并对该软件的优缺点进行了评述。     

参数空间变异性下地下水污染监测网多目标优化机制研究

基于野外实际含水层参数存在空间变异性的客观事实,研发概率Pareto遗传算法(Probabilistic Pareto genetic algorithm,PPGA),用于求解考虑含水层参数空间变异性下地下水污染监测网多目标优化设计问题。PPGA在ε-改进非劣支配遗传算法(epsilon-dominance non-dominated sorted genetic algorithm II,ε-NSGAII)的基础上通过添加概率择优排序和概率拥挤度技术,寻求考虑参数空间变异条件下地下水污染监测网模拟—优化耦合模型的Pareto最优解。将优化结果与蒙特卡洛(Monte Carlo,MC)模拟分析结果进行对比,验证优化结果的可靠性。算例求解结果表明:在求解考虑参数空间变异性条件下地下水污染监测网多目标优化设计问题时,PPGA优化所得Pareto最优解变异性小,可靠性高,可为决策者提供一系列稳定可靠的监测方案。     

北京市朝阳区沙子营垃圾污染物在潜水含水层中迁移规律的现场模拟实验

选择北京平原区水文地质、环境地质等方面都比较典型的沙子营垃圾堆放场,建立了其水文地质模型。在充分收集资料、分析得出部分计算所需参数后,采用现场弥散实验、勘查取样测试等方法,求得了该含水层的弥散系数等参数;用二维非稳定流溶质运移方程对污染物在此含水层中的迁移扩散规律、速度和污染范围等进行了模拟计算;采用现场钻探、取样测试分析等方法,评价了该含水层的实际污染状况。实验模拟计算和现场调查结果表明:污染物在潜水含水层中的运移规律遵循二维非稳定流场中的溶质运移方程,污染物在潜水含水层中的运移速度约为86.25m/a,迁移扩散主要发生在地下水流向上,侧向扩散宽度极小,是地下水流向上的1/17。     

地下水污染监测网多目标优化设计模型及进化求解

采用模拟-优化方法建立了一个用于地下水污染监测网设计的多目标优化模型,该模型包括最小化监测费用、污染物质量评估误差、污染羽一阶矩评估误差和二阶矩评估误差等4个目标函数,以充分揭示减少地下水污染监测费用与提高污染监测精度之间的权衡关系.将改进小生境Pareto遗传算法与地下水流模拟程序和污染物运移模拟程序相耦合用于求解地下水污染监测网多目标设计模型.算例研究表明,采用进化算法求解监测网的多目标模型,能真实地反映各个目标函数间的权衡关系,并且不用考虑传统方法中惩罚因子的影响.与单目标优化模型相比,多目标优化模型可在较短的时间内得到优化问题的一系列Pareto权衡解,以利于相应条件下决策者选择最为经济有效的地下水污染监测方案.     

1:5万冲洪积扇含水层水质调查方法探讨

以豫北山前冲洪积扇含水层为研究对象,探讨1:5万冲洪积扇含水层水质调查方法。该方法由工作量预部署、野外调查和取样点布设3个阶段组成。工作量预部署阶段,以前期调查资料为基础,编制工作区地形地貌、潜在污染源、地下水流场和水化学类型演变图,并赋值叠加形成不同调查点密度分区的野外工作图,宏观控制调查工作量; 野外调查阶段,重视pH值、溶解氧、电导率和氧化还原电位的测定与质量控制,根据调查区块环境状况调整调查点密度,布点兼顾地下水点与非地下水点及不同深度水点间的比例; 取样点布设阶段,以地质环境对含水层进行分区,以不同深度水井地下水功能对含水层进行分层。利用调查点现场水化学参数,识别并圈定水质变异区,分别采用判断法和均布法布设水质变异区和水质正常演化区取样点。该方法可为类似地区开展高精度地下水水质或污染调查提供参考。     

基于水力走时反演刻画裂隙含水层非均质性

获取含水层水力参数空间非均质分布信息是研究地下水渗流、地下水污染物运移等诸多地下水问题的重要基础。然而,受常规勘察技术所限,含水层的非均质性难以直接刻画,尤其在裂隙介质含水层中,水力参数分布的非均质性更加突出,进一步增加了刻画的难度。针对该问题,本研究首先通过德国哥廷根大学北校区内的Neutra试验场地进行的64次抽水试验,获取了该场地不同深度观测点的水头响应曲线,然后使用解析方法对64组数据进行分析和参数评估,同时采用水力走时的方法对井间水力参数分布进行反演计算,最后将得到的结果分别与经典解析解参数估计和热示踪试验结果进行对比验证。结果表明,解析解参数估计的结果虽然能够在一定程度上展现含水层的垂向非均质性,但是无法刻画井间含水层水力参数的非均质分布;与热示踪试验结果的对比验证了基于水力走时反演的水力层析法在刻画裂隙介质水力参数的空间非均质分布的可靠性。     

地理信息系统在地下水监测网设计方面的应用

地下水中污染物检测的有监测设置可以通过分析候选样品位置和含水带对污染影响之间的空间关系来设计，该地带的实例位于污染源范围，可能污染迁移带和供水井区。地理信息系统最适应于完成主要的地下水监测网设计的任务，例如计算确定接近候选位置对高污染敏感区距离变量值和通过型利用这些变量确定相关监测值。通过一研究实例，本文勾划出ＧＩＳ基本测定地下水质监测网设计利用的轮廓，其结果认为ＧＩＳ分析空间参考数据性能可促进建     

某石化污染场地含水层自然降解BTEX能力评估

地下水污染问题是国家关注的重要环境问题之一,监测与评估含水层自然降解污染物能力是防治地下水污染的基本手段,也是国外地下水中修复技术研究热点课题。以华北平原某石油化工类场地为案例,通过调查场地水文地质条件、土壤及地下水污染现状,监测场地地下水中苯系物(BTEX)浓度及相关化学参数变化,运用微生物水文地球化学方法和水文地质方法,估算了该场地含水层自然降解苯系物量。这一研究成果为评估我国石油类场地地下水污染的自然修复能力提供了实证和基础数据。     

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.     

渗透系数的空间变异性对污染物运移的影响研究

随机水文地质学方法,较传统的确定性方法而言,是解决非均质含水层中水流和溶质运移问题的一种更为合理的手段。据以往研究,假设渗透系数场遵循对数正态分布,利用直接傅立叶变换方法来生成渗透系数随机场。应用基于随机理论的蒙特卡罗方法,来研究渗透系数的空间变异性对污染物运移结果的影响。实例研究表明,污染物在含水层中运移过程中污染羽的展布范围(二阶矩)随着渗透系数空间变异方差的增大而扩大,而污染羽在空间上的质心位置(一阶矩)基本不受方差的影响,仅取决于渗透系数随机场的均值大小。另外还分析了污染羽在各点的浓度变化方差和变异系数分别随渗透系数变异方差的变化状况。     

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

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

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.     

Application of electrical imaging to leachate plume evolution studies under in-situ and model conditions

Electrical imaging (resistivity tomography) is increasingly used as a geophysical exploration technique in contaminated land research. The present work demonstrates the efficiency of electrical imaging in monitoring pollution plume evolution under both in-situ and model conditions. The in-situ campaign was done at an unlined landfill site situated in the city of Patras, Hellas. A partially saturated zone of alluvial fan deposits underlies the site, which retards leachate percolation to the underlying water table. Electrical imaging provided details of the internal structure of the waste tip, and confirmed the presence of a leachate plume beneath the base of the landfill. This field data component provided the constraints for the design of a generic model of contaminant infiltration into partially saturated sand. The aim was to study leachate plume evolution in a laboratory environment. For this purpose, a miniaturised resistivity tomography technique, developed at Cardiff University, was used to image the resistivity distribution before during and after contaminant infiltration. Comparison of resulting two-dimensional tomography with observed plume geometry at the end of the test showed the miniaturised electrical imaging technique to be highly effective. Experiment showed that contaminant evolution taking place in the model was gravity-driven, with capillary water in the vadose zone being displaced by the denser contaminant solution.     

Boundary element analysis of contaminant transport in fractured porous media

A two-dimensional boundary integral method to analyse the flow of contaminant in fractured media having a two- or three-dimensional orthogonal fracture network is presented. The method assumes that the fractures provide the paths of least resistance for transport of contaminants while the matrix, because of its low permeability, acts as ‘storage blocks’ into which the contaminant diffuses. Laplace transform is used to eliminate the time variable in the governing equation in order to facilitate the formulation of a boundary integral equation in the Laplace transform space. Conventional boundary element techniques are applied to solve for the contaminant concentrations at specified locations in the spatial domain. The concentration in the time domain is then obtained by using an efficient inversion technique developed by Talbot. The method is able to analyse the behaviour of waste repositories which have diminishing concentration due to the mass transport of the contaminant into the surrounding fractured media.     

An assessment of miniaturised electrical imaging equipment to monitor pollution plume evolution in scaled centrifuge modelling

This paper demonstrates the capabilities and limitations of a miniaturised electrical imaging technique (resistivity tomography) developed at Cardiff University to image contaminant plumes in scaled centrifuge models of the vadose zone. For this purpose a generic model of contaminant infiltration into unsaturated sand was designed. The imaging technique produces two-dimensional contoured plots of the resistivity distribution before and during contaminant infiltration experiments. During the experiments, dyed NaCl solution was released into the model and the change in resistivity associated with the contaminant plume evolution was imaged as a function of time and g-level. Capillary pressure was monitored constantly by matrix potential probes (tensiometers) in order to investigate the effect of capillary forces on plume evolution. Tests at 1g (static conditions) and 10g are described in this paper. Comparison of resulting two-dimensional tomography with observed plume geometry at the end of the 1g test showed this imaging technique to be highly effective.

Contaminant plume evolution in the unsaturated sand model was observed to be mainly gravity-driven, with plume migration and geometry being strongly affected by a tenfold increase in gravity in the centrifuge experiment. It is concluded that miniaturised electrical imaging can be a useful tool for monitoring pollution plume evolution during centrifuge tests, but when plume evolution is rapid, the time taken to interrogate each array restricts the effectiveness of the technique in monitoring changes in plume geometry. However, in such cases, resistivity tomography does provide valuable information on residual levels of contaminant fluid retained within the soil after passage of the plume.     

浅层地下水氯代烃污染天然衰减速率的估算

天然衰减恢复技术是恢复和控制浅层地下水氯代烃污染的技术之一,如何简便获取可靠的氯代烃衰减速率常数是该技术应用的一个关键。趋势线分析方法是一种简便有效的方法,在污染羽状体稳定的条件下,通过地下水流向上至少3口监测井的资料,能够比较准确地估算出污染物的天然衰减速率常数和生物降解速率常数。某氯代烃污染典型区的应用实例研究表明,该区四氯乙烯(PCE)的天然衰减速率常数和生物降解速率常数分别为0·000925d-1和0·000537d-1,证实该区浅层地下水中的PCE存在天然生物降解,但降解速率比较缓慢。忽略吸附作用的天然衰减容量计算所得出的天然衰减速率常数明显小于实际结果,说明尽管典型区包气带及含水层介质的有机碳含量很少,但它们对PCE的吸附作用不容忽视。