基于参数敏感性分析的轻质非水相液体在非均质包气带中迁移模拟研究

利用TOUGH2程序构建室内实验尺度多相流数值模型,基于扰动分析法,对模型中的介质孔隙度、绝对渗透率、残余水饱和度(S_wr)、残余非水相(NAPL)饱和度(S_nr)、进气值参数(α)进行敏感性分析,在此基础上设计中砂夹粗砂层和夹细砂层两种非均质条件情景,探究轻质非水相液体(LNAPL)以不同速率泄漏以及地下水位波动时的迁移过程。结果表明：(1)模型参数敏感度呈绝对渗透率>S_nr>α>S_wr>孔隙度的规律,S_nr、α、S_wr与LNAPL稳定饱和度呈正线性相关,绝对渗透率与其呈负线性相关。(2)主要由于介质间绝对渗透率的差异,中砂介质中的细砂夹层阻滞、粗砂夹层促进LNAPL迁移,泄漏速率的增大致使细砂夹层阻滞作用更为显著。(3)地下水位波动扩大LNAPL横、纵向污染范围,降低夹层内LNAPL饱和度,水位上升,细砂夹层下界面、粗砂夹层上界面阻滞LNAPL抬升,主要由于介质绝对渗透率的差异,两种非均质情景对地下水位波动的响应特征不同。     

颗粒层中非水相液体污染物冲洗排出速度的数值模拟简

在进行污染土壤的冲洗修复时，模型模拟是极具优势的预测净化效果的方法。本研究分别采用5种粒径的球形玻璃微珠模拟实际土壤，实验考察了玻璃微珠颗粒层中水和非水相液体的排出速度，给出了颗粒层中液含率分布函数，以及最大和最小颗粒层孔隙直径的定义，建立了一个描述颗粒层中水和非水相液体排出速度的多孔介质孔隙毛细管束模型。使用构建的模型模拟水和异辛烷在颗粒层中的排出速度，模拟结果与实验数据有很好的一致性，模型预测结果的平均相对误差小于6.4%，特别是对粒径较小颗粒层，预测的准确性更高。模型可作为进一步研究在非水相液体饱和土壤水冲洗修复时，水和非水相液体排出速度和修复效果的基础。     

不同粒径高浓度粉末活性炭组合UF膜工艺特征和过滤效果

采用一体式高浓度粉末活性炭-超滤(HCPAC-UF)组合工艺,处理高氨氮微污染水源水。选取3种不同粒径(50、3和0.2μm)的粉末活性炭考察PAC粒径变化对组合工艺处理效果的影响。结果表明,PAC粒径改变分别对UV254表征的腐殖质类有机物和DOC表征的溶解性有机物的去除影响较小;PAC粒径显著影响系统中的微生物的生长环境,PAC粒径越小,系统中微生物的生长受抑制越明显,同时,PAC的粒径越小对减缓HCPAC-UF工艺的膜污染越显著。研究结果对于优化HCPAC-UF工艺参数具有一定参考价值。     

流化床结晶软化高硬度地下水及其晶体生长动力学

采用流化床结晶软化的方法对北京市城区某地下水进行软化处理,研究考察了软化药剂投加量、诱晶颗粒粒径、诱晶颗粒填料高度、入水流速、碳酸钙过饱和度等参数对水质硬度去除的影响。研究发现,当药剂投加量为220 mg·L~(-1)时,出水钙硬度去除率可达90%以上,总硬度去除率达到60%,出水总硬度可达300 mg·L~(-1)左右。使用石英砂颗粒作为诱晶材料,粒径为0.21~0.43 mm。当诱晶材料填料高度为40 cm,入水流速控制在40~70 m·h-1时,均可得到较好的出水水质。优化反应参数,可控制处理后的出水总硬度在300 mg·L~(-1)左右,即去除235 mg·L~(-1)的暂时硬度,煮沸的水不再浑浊,无水垢生成,并且出水pH可控制在8.5以下。此外,探究了多相流化过程中诱晶颗粒的线性增长速率以及床层增长速率与过饱和度、颗粒粒径、入水流速的关系,通过计算得到晶体线性增长速率相关表达式。     

多孔介质中轻非水相液体挥发过程的电阻率成像法监测

在含有细砂透镜体的三维非均质砂槽中进行了轻非水相液体(LNAPLs)的挥发实验,并应用电阻率成像法进行同步动态监测,尝试通过电阻率图像的变化反映砂槽中LNAPLs的挥发过程。结果表明,砂槽注油后电阻率明显上升,LNAPLs污染中心区域的电阻率大于边缘区域,由上往下电阻率减小,其图像近似半球形。注油结束后,污染物逐渐挥发,污染区域LNAPLs浓度不断减小,表现为电阻率逐渐降低。监测过程表明,应用电阻率成像法监测LNAPLs挥发过程是有效的。     

重金属污染土壤旋流洗脱设备的设计及性能评估

针对高浓度重金属污染土壤,尤其是污染负荷较高的黏土土壤,传统的物化方法难以实现其高效的洗脱。利用旋流场中土壤颗粒高速自转/公转,实现土壤颗粒污染物的强化快速脱附。土壤旋流洗脱实验分为旋流器的分离性能和单一/复合污染物的脱附性能2部分。在土壤-水体系下,旋流器的最优操作条件为:进口流量0.7 m~3·h~(-1),分流比0.12,固液比为1:20。对于Pb污染物,底流脱附效率均能达到近85%,溢流也能够达到70%,对于Cu污染物,底流和溢流脱附均能达到90%左右。对于Cr(Ⅵ)污染物,底流脱附最高能达到60%左右,但溢流洗脱效率极低。复合污染能够在单次通过后脱除Pb、Cu、Cr(Ⅵ)等多种重金属污染,且洗脱效果与单一污染洗脱时基本一致。对实际的场地修复具有指导意义。     

废弃铬盐厂土壤中铬的赋存特征及异位淋洗修复可行性研究

土壤理化性质、污染物的分布及赋存特征是土壤异位淋洗技术是否可行的关键依据。以西北某废弃铬盐厂铬污染场地土壤作为研究对象,测定了土壤基本理化参数及各粒径污染物的浓度,探究总铬和六价铬在不同粒径土壤中的分布特征,并通过铬的赋存形态分析确定清水作为淋洗剂,开展异位清水淋洗实验。结果表明:供试土壤粗颗粒物料(粒径分别为:9.50 mm、2.00~9.50 mm、0.841~2.00 mm及0.25~0.841 mm)占土壤总质量的比例达到93.36%,总铬和六价铬污染物富集于土壤细颗粒(粒径0.25 mm);同时,土壤各级粒径颗粒中水溶态和弱酸提取态铬所占比例均大于50%;以清水做淋洗剂时,总铬和六价铬的去除率分别高于60%和80%,且淋洗后土壤各级粒径中六价铬含量均低于30 mg·kg~(-1),达到修复标准要求。     

不同地表雨水径流冲刷特性分析

选择北京市某城区典型地表(校园、自行车道、商业街和高架桥)的雨水径流作为研究对象,系统分析了不同地表雨水径流的冲刷特性。分析结果表明:不同地表初期冲刷现象都比较明显,控制不同地表20%的初期雨水可以削减污染物总量的40%～60%,高架桥初期雨水径流中COD、SS和TP浓度最高,校园初期雨水径流中TN浓度最高。不同地表雨水径流的水质优劣顺序为:校园>商业街>自行车道>高架桥。校园和自行车道在整个降雨过程中D50和D90粒径大小相近,分别为30μm和50μm左右,且变化幅度都较小。高架桥和商业街在雨水径流初期D50和D90粒径有明显减小,但后期均出现较大幅度的波动,其粒径明显大于校园和自行车道。上述研究成果对于城市不同地表雨水径流污染控制和场地开发具有重要意义。     

工艺条件对磷回收过程中鸟粪石沉淀颗粒粒径的影响

在鸟粪石沉淀法回收废水中磷的过程中,鸟粪石颗粒的大小将直接影响其沉淀的速率,进而影响鸟粪石的沉淀效果和磷的回收率。本文采用激光粒度分析仪测定鸟粪石的平均粒径,详细考察在小型连续搅拌反应-沉淀磷回收装置中不同的工艺条件下鸟粪石颗粒粒径的变化规律,并结合Stokes公式计算鸟粪石颗粒在废水出口处的沉降速率,为沉淀池的设计提供参考依据。结果表明:鸟粪石的平均粒径在12～25μm之间,沉降速率在5.46×10-5～2.37×10-4m/s之间。随着反应室水力停留时间的延长,鸟粪石颗粒的粒径逐渐增大,当停留时间超过18 min时,颗粒的粒径基本不变;随着沉淀室水力停留时间的延长,鸟粪石颗粒的平均粒径缓慢增大,当停留时间超过70 min后颗粒粒径的变化不大;鸟粪石颗粒的平均粒径在一定程度上受废水中磷初始浓度变化的影响,在磷初始浓度为62～128 mg/L时颗粒的粒径变化不大,当磷浓度为496 mg/L时粒径有较大增加,此时鸟粪石颗粒的沉降速率也大幅度增加;鸟粪石颗粒的平均粒径受pH值的影响不大;随氮磷摩尔比的增大,鸟粪石颗粒的平均粒径略有增加;随镁磷摩尔比的增大,鸟粪石颗粒的平均粒径逐渐减小,沉降速率则有明显的下降。     

纳米铁原位注入技术对六价铬污染地下水的修复

纳米零价铁原位注入是六价铬污染地下水的有效修复技术之一。为验证其场地修复效果,在北京某电镀厂搬迁后的遗留六价铬污染场地,现场制备纳米铁并通过原位注入对场地内六价铬污染地下水进行原位修复现场中试研究。选取6 m×6 m实验场地,在地下水中六价铬污染浓度最高为2 mg·L~(-1)的中心点设置注射井并在四周设置4口监测井。实验结果表明,原位注入纳米铁药剂具有良好修复效果,修复后该实验场地范围内地下水中六价铬浓度均低于地下水质量Ⅳ类标准0.1 mg·L~(-1),注射井中六价铬还原率达到99%。通过对注射井及监测井地下水中pH、溶解氧、氧化还原电位等检测指标进行跟踪监测和相关性分析,发现氧化还原电位与污染物浓度变化相关性最强,可以作为污染物变化的表征因子。该中试研究结果对于六价铬污染地下水的原位修复具有重要的参考价值。     

Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone

A quantitative two-dimensional laboratory experiment was conducted to investigate the immiscible flow of a light non-aqueous phase liquid (LNAPL) in the vadose zone. An image analysis technique was used to determine the two-dimensional saturation distribution of LNAPL, water and air during LNAPL infiltration and redistribution. Vertical water saturation variations were also continuously monitored with miniature resistivity probes. LNAPL and water pressures were measured using hydrophobic and hydrophilic tensiometers. This study is limited to homogeneous geological conditions, but the unique experimental methods developed will be used to examine more complex systems. The pressure measurements and the quantification of the saturation distribution of all the fluids in the entire flow domain under transient conditions provide quantitative data essential for testing the predictive capability of numerical models. The data are used to examine the adequacy of the constitutive pressure-saturation relations that are used in multiphase flow models. The results indicate that refinement of these commonly used hydraulic relations is needed for accurate model prediction. It is noted in particular that, in three-fluid phase systems, models should account for the existence of a residual NAPL saturation occurring after NAPL drainage. This is of notable importance because residual NAPL can act as a non negligible persistent source of contamination.     

Evaluating LNAPL contamination using GPR signal attenuation analysis and dielectric property measurements: practical implications for hydrological studies

Groundwater and sub-surface contamination by Light Non-Aqueous Phase Liquids (LNAPLs) is one of the industrial world's most pressing environmental issues and a thorough understanding of the hydrological, physical and bio-chemical properties of the sub-surface is key to determining the spatial and temporal development of any particular contamination event. Non-invasive geophysical techniques (such as electrical resistivity, electromagnetic conductivity, Ground-Penetrating Radar, etc.) have proved to be successful sub-surface investigation and characterisation tools with Ground-Penetrating Radar (GPR) being particularly popular. Recent studies have shown that the spatial/temporal variation in GPR signal attenuation can provide important information on the electrical properties of the sub-surface materials that, in turn, can be used to assess the physical and hydrological nature of the pore fluids and associated contaminants. Unfortunately, a high percentage of current LNAPL-related GPR studies focus on contaminant mapping only, with little emphasis being placed on characterising the hydrological properties (e.g., determining contaminant saturation index, etc.). By comparing laboratory-based, dielectric measurements of LNAPL contaminated materials with the GPR signal attenuation observed in both contaminated and 'clean' areas of an LNAPL contaminated site, new insights have been gained into the nature of contaminant distribution/saturation and the likely signal attenuation mechanisms. The results show that, despite some practical limitations of the analysis technique, meaningful hydrological interpretations can be obtained on the contaminant properties, saturation index and bio-degradation processes. A generalised attenuation/saturation model has been developed that describes the physical and attenuation enhancement characteristics of the contaminated areas and reveals that the most significant attenuation is related to smeared zone surrounding the seasonally changing water table interface. It is envisaged that the model will provide a basis for the interpretation of GPR data from analogous LNAPL contaminated sites and provide investigators with an appreciation of the merits and limitations of GPR-based, attenuation analysis techniques for hydrological applications.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

Multispectral image analysis method to determine dynamic fluid saturation distribution in two-dimensional three-fluid phase flow laboratory experiments

The need for measuring dynamic fluid saturation distribution in multi-dimensional three-fluid phase flow experiments is hampered by lack of appropriate techniques to monitor full field transient flow phenomena. There is no conventional technique able to measure dynamic three-fluid phase saturation at several array points of the flow field at the same time. A multispectral image analysis technique was developed to determine dynamic NAPL, water and air saturation distribution in two-dimensional three-fluid phase laboratory experiments. Using a digital near-infrared camera, images of sand samples with various degrees of NAPL, water and air saturation were taken, under constant lighting conditions and within three narrow spectral bands of the visible and near-infrared spectrum. It was shown that the optical density defined for the reflected luminous intensity was a linear function of the NAPL and the water saturation for each spectral band and for any two and three-fluid phase systems. This allowed the definition of dimensionless lump reflection coefficients for the NAPL and the water phase within each spectral band. Consequently, at any given time, two images taken within two different spectral bands provided two linear equations which could be solved for the water and the NAPL saturation. The method was applied to two-dimensional three-phase flow experiments, which were conducted to investigate the migration and the distribution of LNAPL in the vadose zone. The method was used to obtain continuous, quantitative and dynamic full field mapping of the NAPL saturation as well as the variation of the water and the air saturation during NAPL flow. The method provides a non-destructive and non-intrusive tool for studying multiphase flow for which rapid changes in fluid saturation in the entire flow domain is difficult to measure using conventional techniques.     

Large-scale effects on resistivity index of porous media

The estimation of humidity in the unsaturated zone of soils and NAPL saturation in contaminated aquifers may be based on the interpretation of electrical resistivity index logs. In the present work, concepts of the theory of the two-phase flow in pore networks are employed to interpret the form of the equilibrium and dynamic resistivity index curves of large porous samples. A resistivity cell is constructed to measure the capillary and electrical properties of large samples of unconsolidated porous media. The drainage capillary pressure and resistivity index curves of a sand column are measured by using the micropore membrane (porous plate) method, where a 0.5% wt/vol NaCl aqueous solution is displaced by n-dodecane. The dynamic resistivity index curves are measured by using the continuous injection technique for various orientations of the sand column. Based on concepts of the two-phase flow theory, concerning the dominant displacement growth pattern in a pore network and arising from the cooperative effects of capillary, buoyancy, and viscous forces, approximate relationships are developed for the resistivity index and saturation exponent as functions of the water saturation. The saturation exponent decreases as the displacement advances and the fluid distribution across the sand column tends to be homogenized after oil breakthrough. Both the resistivity index and saturation exponent increase as the displacement pattern tends to become compact and stable. In the destabilized flow pattern, as the Bond number decreases, the resistivity index may increase respectably within a narrow range of values of the Bond number. This happens when the thickness of the unstable capillary finger exceeds the lateral dimension of the porous sample and becomes a fractal percolation cluster. The saturation exponent becomes almost constant and independent of water saturation only over the destabilized displacement pattern at high values of the Bond number.     

地下水中轻质有机污染物(LNAPL)透镜体研究

在二维砂槽模型中模拟了轻质油在均匀多孔介质地下水非饱和区中的运移过程。模拟结果表明,地下水毛细区是轻质油污染的重点区,除了ＬＮＡＰＬ的残留以外,进入地下水饱和的ＬＮＡＰＬ终将被地下水顶托回到毛细区中,毛细区以上的约大多民将在重力作用下进入毛细区中,试验中观察到达稳定状态时ＬＮＡＰＬ透镜体的上边缘略微高出毛细区。利用多孔介质毛细管模型,建立了利用界面张力、接触角、介质特征孔隙直径等物理量估算不同位置     

非饱和土壤轻油污染多相流研究进展

轻油泄漏到地下后 ,在非饱和土壤中水、气、油三相相互作用、相互影响 ,共同构成多相流体系统。在简要分析轻油污染问题的产生、危害以及轻油泄漏在地表以下运移方式的基础上 ,介绍并总结了近些年来国内外土壤轻油污染有关的多相流研究方面的进展 ,包括国外的模型研究成果、国内外试验研究方法和研究成果 ,并指出目前研究中存在的问题和需要进一步探讨的方向     

Infiltration and redistribution of LNAPL into unsaturated layered porous media

Enhanced understanding of light non-aqueous phase liquid (LNAPL) infiltration into heterogeneous porous media is important for the effective design of remediation strategies. We used a 2-D experimental facility that allows for visual observation of LNAPL contours in order to study LNAPL redistribution in a layered porous medium. The layers are situated in the unsaturated zone near the watertable and they are inclined to be able to observe the effect of discontinuities in capillary forces and relative permeabilities. Two experiments were performed. The first experiment consisted of LNAPL infiltration into a fine sand matrix with a coarse sand layer, and the second experiment consisted of a coarse sand matrix and a fine sand layer. The numerical multi-phase flow model STOMP was validated with regard to the experimental results. This model is able to adequately reproduce the experimental LNAPL contours. Numerical sensitivity analysis was also performed. The capillarity contrast between sands was found to be the main controlling factor determining the final LNAPL distribution.     

A field trial to assess the performance of CO2-supersaturated water injection for residual volatile LNAPL recovery

A pilot scale field trial was conducted to evaluate the recovery of volatile, light non-aqueous phase liquids (LNAPLs) using a novel remediation method termed supersaturated water injection (SWI). SWI uses a patented technology to efficiently dissolve high concentrations of CO₂ into water at elevated pressures. This water is injected into the subsurface resulting in the nucleation of CO₂ bubbles at and away from the injection point. The nucleating bubbles coalesce, rise and volatilize residual LNAPL ganglia. In this study, an LNAPL composed of 103 kg of volatile pentane and hexane, and 30 kg of non-volatile Soltrol was emplaced below the water table at residual saturation. The SWI technology removed 78% of the pentane and 50% of the less volatile hexane. Contaminant mass was still being removed when the system was shut down for practical reasons. The mass removed is comparable to that expected for air sparging but a much smaller volume of gas was injected using the SWI system.     

滤料粒径对给水曝气生物滤池运行稳定性的影响

利用给水曝气生物滤池(UBAF)处理刘屋洲水源水,考察了不同陶粒粒径对氨氮、高锰酸盐指数的去除效果以及水头损失变化的影响。结果表明,陶粒3～5 mm的UBAF对氨氮、高锰酸盐指数的去除效果较高,但是其水头损失较大,且反冲洗前后水头损失变化偏大,不利于UBAF稳定运行。陶粒6～10 mm的UBAF运行状况良好,氨氮的去除率为73.8%,高锰酸盐指数的去除率为18.6%,滤池24 h过滤水头损失<0.4 m,反冲洗前后过滤水头损失变化量<5 cm。其出水氨氮和高锰酸盐指数均达到《生活饮用水卫生标准》(GB5749-2006),滤池的水头损失小,反冲洗前后也变化甚微。进一步的实验表明陶粒破碎会导致滤料粒径级配变化,使得水头损失增加。