非均质介质中重非水相污染物运移受泄漏速率影响数值分析

重非水相污染物(DNAPL)在地下介质中运移和分布受多种因素控制,包括DNAPL本身的物理化学性质,土的性质,泄漏条件等等。由于介质的非均质性,使得多相流运移行为更为复杂。基于地下水随机理论构建渗透率随机场,采用蒙特卡罗方法探讨泄漏速率对非均质饱和介质中DNAPL运移的影响。数值结果表明,在泄漏总量一定的情况下,泄漏速率越低,介质非均质性对DNAPL运移的影响程度越高。反之,DNAPL的渗漏速率越高,小尺度地层的非均质性影响越低。由于DNAPL运移过程中在垂直方向受重力的影响,污染羽在空间上的质心位置(一阶矩)以及展布范围(二阶矩)在垂直方向上的变异程度要高于水平方向。     

多孔介质界面对重非水相液体迁移过程影响的图像法研究

重非水相液体(dense nonaqueous phase liquid,DNAPL)污染土壤和地下水的问题已引起广泛关注,研究其在不同粒径多孔介质及其界面的运移特征形态是确定污染区域、修复治理土壤和地下水环境的前提。文章通过室内试验研究多孔介质界面对DNAPL运移与分布特性的影响。首先在二维砂槽上进行DNAPL污染物的入渗试验,试验过程中用数码相机拍照,将DNAPL扩散过程以图像的形式记录下来;然后用AutoCAD对图片进行处理,绘制出DNAPL迁移过程的锋面变化图。结果表明:DNAPL入渗过程中,迁移主要受到重力作用与毛细作用的控制,毛细作用力随着介质粒径的增大逐层减小,重力作用逐渐起主导作用使污染物入渗速度逐层增大;介质结构影响DNAPL的迁移形态,介质粒径逐层增大,DNAPL污染物的渗流面与指进扩散宽度逐层减小,扩散方式由面状变为指状;在不同粒径介质界面介质结构发生突变时,DNAPL迁移锋面线曲率也相应变大,此时DNAPL的迁移呈现“凸”型特征,另外,不同的界面横向扩散的滞留宽度不同,随着介质粒径的增大,界面的横向扩散宽度相对变短。     

非均质孔隙介质中两相流的光透法应用研究

地下水中非水相液体(Non-aqueous Phase Liquid,NAPL)的流动及其曝气修复技术是典型的两相流问题。基于实际地下水含水介质的普遍非均质性,本文应用光透法对非均质孔隙介质中两相流进行了定量试验研究,设计了两组砂箱实验,研究气体和重非水相液体(DNAPL)在非均质孔隙介质中的迁移规律,应用了水/气两相和水/NAPL两相饱和度计算模型。实验结果表明:气体主要由不规则通道向上运动,遇到低渗透性透镜体时在其下方堆积,并开始横向运动,绕过透镜体后继续向上运动,最终在砂箱顶部形成连续气体分布,注气速度越大,气体运移范围越宽;DNAPL在自身重力作用下克服毛管压力向下迁移至低渗透性透镜体,DNAPL无法克服该介质的毛管压力,停止垂向入渗,并在其表面堆积,开始横向运移,绕过透镜体后继续向下运动,最终在砂箱底部形成连续DNAPL污染池。均质介质中建立的计算流体饱和度的水/气模型及水/NAPL模型与实验结果较吻合,可用于非均质多孔介质中水/气相和水/NAPL相饱和度的计算。     

基于集合卡尔曼滤波的多相流模型参数估计——以室内二维砂箱中重质非水相污染物入渗为例

重质非水相有机污染物(DNAPL)泄漏到地下后,其运移与分布特征受渗透率非均质性影响显著。为刻画DNAPL污染源区结构特征,需进行参数估计以描述水文地质参数的非均质性。本研究构建了基于集合卡尔曼滤波方法(EnKF)与多相流运移模型的同化方案,通过融合DNAPL饱和度观测数据推估非均质介质渗透率空间分布。通过二维砂箱实际与理想算例,验证了同化方法的推估效果,并探讨了不同因素对同化的影响。研究结果表明:基于EnKF方法同化饱和度观测资料可有效地推估非均质渗透率场;参数推估精度随观测时空密度的增大而提高;观测点位置分布对同化效果有所影响,布置在污染集中区域的观测数据对于参数估计具有较高的数据价值。     

渗透率空间变异性对重非水相流体运移的影响

中针对多孔介质中水、气和非水相流体的多相流运移特点,建立了由包气带至饱和带的 NAPLs- 水 - 气三相流模 型。用 Karhunen-Loeve 展开方法构建 4 个空间变异性不同的渗透率随机场。在此基础上,分别模拟 4 种情景重非水相流体 (DNAPLs)的泄漏运移过程。通过分析 DNAPLs 饱和度的空间分布以及其空间矩随时间的变化来探讨渗透率空间变异性对 DNAPLs 运移的影响。结果显示,随着渗透率空间变异性的增强,DNAPLs 饱和度的空间变异性增大。渗透率场中的优势通 道对 DNAPLs 运移有着显著影响,直接决定 DNAPLs 的运移路径和饱和度分布。渗透率场的各向异性降低了 DNAPLs 在垂向 上的运移速度,影响污染羽的形态。     

多孔介质中非水相流体运移的数值模拟

针对多孔介质中水、气和非水相流体(NAPLs)的多相流动特点,建立了非水相流体(NAPLs)污染物迁移模型,分析了非水相流体在土壤非饱和区和地下水系统中的运移规律。通过有限元数值解对轻非水相流体和重非水相流体在土壤系统中的迁移过程进行模拟,得到了污染物的时空分布特征和污染范围。计算结果表明,数值模拟方法能够合理地描述非水相流体的运移过程和污染特征。土体渗透性和污染物残余饱和度是其重要影响因素。     

盐度对多孔介质中DNAPL运移和分布的影响

选用四氯乙烯（PCE）作为典型DNAPL污染物,以NaCl作为地下水中溶解盐代表,研究盐度对DNAPL在饱和多孔介质中运移和分布的影响。通过批次实验测定NaCl水溶液/石英砂/PCE三相体系下的接触角和界面张力,结果表明,PCE在石英砂表面的接触角随着水中NaCl浓度的增大而减小,而PCE和NaCl水溶液的界面张力随着NaCl浓度的增大而增大,尤其当氯化钠浓度较高时（>0.1 mol/L）,影响程度更为显著。在此基础上,采用透射光法监测不同介质情景下DNAPL在二维砂箱中的运移和分布,定量测定DNAPL在介质中的饱和度。实验结果表明,地下水盐度的增加将促进DNAPL的垂向入渗,减少被截留在运移路径上的DNAPL量,使得DNAPL运移路径及累积形成的池状DNAPL（pool）向水流方向偏移。在均质多孔介质和含有透镜体的非均质多孔介质中,随着盐度的增加,DNAPL在横向和垂向上的展布均呈现出增加趋势,导致污染源区变大,且介质中以离散状存在的DNAPL量明显增加。     

水动力和介质非均质性影响下DNAPLs运移数值模拟

选定氯苯为典型DNAPL建立理想的地表污染泄漏多相流运移模型,利用TOUGH系列软件TMVOC模块进行数值模拟,探讨地下水动力、介质非均质性及两种因素协同作用时对DNAPLs运移和分布的影响。研究结果表明:(1)水力梯度增大了DNAPL水平和垂向运移速率;当水平向水力梯度过大时,反而导致垂向速率减小。(2)透镜体主要发挥阻滞作用,绕过透镜体后,下游侧污染指入渗深度明显大于上游侧;水力梯度的存在进一步增强介质非均质性的影响,产生更多集聚和绕流,形成不规则的污染指和污染池交替组合形态。     

基于转移概率随机模拟的DNAPL运移二维和三维数值模拟对比研究

非均质介质的空间维度变化对重非水相流体（DNAPL） 的运移具有重要的影响。在充分考虑地质体的空间连续 性、不对称性以及各向异性等特征的基础上,采用基于马尔可夫链的转移概率（transition probability） 模型来构建非均质 随机场。该文通过TMVOC-MP软件来模拟DNAPL在非均质介质中的运移规律,探讨非均质随机场的水平空间连续性、 空间维度变化以及侧向运移过程对DNAPL运移的影响。结果表明,介质的水平空间连续性越好,DNAPL在水平方向的 迁移范围越大,在垂向的迁移范围越小;相比于三维模型,二维模型中DNAPL在水平方向的展布更大、在透镜体上的蓄 积量更多,在实际应用中以二维模型代替三维模型会加大模拟结果与实际污染情况之间的误差;侧向运移过程削弱了单 个平面的非均质性对DNAPL运移的控制,当存在侧向运移时,DNAPL绕过透镜体所运移的距离以及在透镜体上的蓄积 量会相应减小。     

结合蒸汽和空气注入修复多孔介质中DNAPL污染物的多目标多相流模拟优化

基于数值方法模拟对比分析仅注入蒸汽与同时注入蒸汽和空气两种情形下原位曝气方法修复多孔介质中DNAPL污染物(以TCE为例)的修复过程。结果显示蒸汽的注入速率和焓值对蒸汽锋面在含水层中的运移速率影响显著;同时注入蒸汽和空气能够有效抑制只注入蒸汽时出现的DNAPL再冷凝现象,从而提高修复效率。基于模拟-优化方法对同时注入蒸汽和空气修复的理想范例建立了多目标多相流优化模型,采用一种自适应通用多目标优化算法AMALGAM求解Pareto锋面。优化结果表明基于混合优化框架的AMALGAM优化算法能有效综合利用多种独立进化算法的优点,自适应地根据独立进化算法的优劣生成子代,保证子代种群最优,收敛速度较NSGA算法更快。     

重非水相液体与水相界面面积的二维实验研究

多孔介质中重非水相液体(Dense Non Aqueous Phase Liquid,DNAPL)与水相的界面面积是影响DNAPL在地下水中溶解速率的关键参数。本文通过二维均质砂箱实验,运用界面分配示踪法和体积分配示踪法分别测定DNAPL与水相的界面面积(单位介质体积内的面积,Anw)和DNAPL饱和度(Sn),并研究两者之间的相关关系。使用C18H29NaO3S(SDBS)作为界面分配示踪剂,Alcohol(2,4 Dimethyl 3 Pentanol)作为体积分配示踪剂,CaBr2作为保守示踪剂。利用Fe2+催化的Na2S2O8对DNAPL进行原位氧化修复,改变Anw和Sn。结果表明,通过定量(用孔隙体积数PV值衡量)加入Na2S2O8,Anw由206 cm2逐渐降为37 cm2,下降速率约为14 cm2/PV;Sn由1.34%逐渐降为0.33%,下降速率约为0.1%/PV;Anw和Sn之间具有良好的线性关系(Anw=146×Sn,R2=0.987)。此关系可用于估算多孔介质中DNAPL与水相之间的界面面积。对实验影响因素的分析表明,孔隙度的增大(变化率为20.7%)使示踪剂穿透曲线出现肩现象和拖尾现象,导致保守示踪剂的保留时间增加5.6%,界面面积减小6.8%。     

二维孔隙介质中重非水相液体饱和度测定与分析

在二维均质饱和孔隙介质中利用光透法与高密度电阻率法定量监测重非水相液体(DNAPL)的运移过程与饱和度。 将光透法与高密度电阻率法相结合,进行高密度电阻率法对 DNAPL 在饱和多孔介质中渗流过程的定量监测。在二维砂箱 中进行 DNAPL 的入渗试验,应用光透法与高密度电阻率法进行动态监测,分别利用 CCD 相机、LCR 数字电桥采集数据。 光透法的数据处理采用“水-DNAPL”两相中 DNAPL 饱和度的计算公式。而高密度电阻率法的数据处理是将获得的电阻值 数据转化为电阻率值,然后利用 Archie 公式获得重非水相液体的饱和度空间分布。Archie 公式中参数β值本文采用与光透法 相结合的新方法来获取,克服了传统方法获取β值时存在的困难。根据不同时刻的饱和度空间分布可以估算注入砂箱内的 DNAPL 总量,并将其与实测的入渗量进行对比。结果表明,与光透法结合后的高密度电阻率成像法能定量监测饱和孔隙介 质中 DNAPL 的入渗过程,且估算的 DNAPL 入渗量与实测值比较吻合,提出的基于光透法计算β值的新方法是有效且可行的。     

多孔介质中水气交替注入微观渗流模拟

针对水气交替注入(water-alternating-gas,WAG)过程中,油气水三相渗流的微观机理认识不足和油气水三相流体在多孔介质中分布规律认识不准确等问题,基于三维孔隙网络模型,应用孔隙级模拟方法,从微观角度模拟了不同润湿性多孔介质中的WAG驱替过程.结果表明:连通性较好的多孔介质中,原油主要在前两轮的WAG循环中被驱替出来;在前两轮WAG驱替之后,流体饱和度和分布规律达到比较稳定的状态,但在完全水湿模型中油相仍然在多孔介质中流动.得出的WAG驱替过程中各相流体饱和度的变化规律、各相流体分布规律和驱替类型,较好地阐述和解释了多孔介质中的微观驱替机理.     

多孔介质中毛细压力、饱和度和相对渗透率的确定方法

目前石油溢出或者地下储油罐泄漏等原因引起的土壤和地下水非水相流体（NAPLs）污染问题越来越引起人们的关注，由NAPLs、水和气所组成的两相或三相系统中的多相流问题亦是当前的研究热点。其中毛细压力（h），饱和度（S）和相对渗透率（k）是多孔介质多相流研究中的三个重要参数，在多相流室内试验研究中是主要的物理监测量，而且三者之间基本关系式的确定是多相流模拟时进行流动控制方程求解的前提条件。本文从室内试验和模型关系两个方面综述了土壤中NAPLs、水和气所组成的多相流系统中毛细压力、饱和度和相对渗透率以及它们之间相关关系的确定方法。     

裂隙介质中流速及倾斜度对四氯乙烯的运移影响研究

文章选用四氯乙烯（PCE）为代表性重非水相液体（DNAPL）,采用透射光法监测PCE在单裂隙介质中的运移过程和形态分布,探究流速及裂隙面倾斜度对PCE运移分布行为的影响。结果表明,单裂隙介质中PCE呈团块状脱离注样针,并在运移过程中呈近椭圆状。地下水流速越大,初始面积Fs越大,近椭圆的PCE离心率e越小;裂隙面倾斜度越大,初始面积Fs越小,近椭圆的PCE离心率e越小。不同流速或倾斜度条件下,入渗稳定后的PCE团块面积与离心率均明显下降。地下水流速及裂隙面倾斜度的增大对于PCE纵向运移具有促进作用,使PCE污染物前端锋面运移速率增大,在相同时间内更快运移到裂隙底部。PCE入渗停止后,流速增大促使污染源区PCE饱和度降低,使饱和度峰值随水流速度增大而减小;倾斜度增大抑制污染源区PCE饱和度降低,使饱和度峰值随倾斜度增大而增大。     

层状非均质多孔介质中LNAPLs污染物运移与分布的锋面扩展模型

通过建立二维物理模型,研究非饱和层状非均质多孔介质中轻非水相液体LNAPLs(Light Non-a-queous Phase Liquids)的入渗机制与变化特征,建立LNAPLs运移与分布的锋面扩展模型,探讨入渗阶段油流锋面扩展速率的变化规律。结果发现:模型中砂土的渗透系数、孔隙度及油的相对渗透系数对LNAPLs运移的锋面扩展速率影响较大,而油的饱和度与压头变化对其影响较小。当LNAPLs由中砂进入细砂与细砂倾斜透镜体时,砂的渗透系数、孔隙度及油的相对渗透系数减小,这种介质结构面的突变改变了LNAPLs锋面扩展速率,其在中砂中锋面扩展速率快而在细砂与细砂倾斜透镜体中锋面扩展速率慢。此外,当LNAPLs进入干湿界面后,由于毛细作用的增强与含水量的加大,其横向比垂向运移的平均锋面扩展速率大。     

非等温条件下有机污染物在粘土衬垫中的扩散分析

填埋场中垃圾降解作用产生的热量在黏土衬垫中的运移会形成温度场,温度场会对污染物的运移产生影响。文章在多孔介质热力学理论和污染物运移理论相结合的基础上,主要考虑温度场对污染物扩散的影响,建立了非等温条件下有机污染物在黏土衬垫中的一维扩散模型,并采用分离变量法得到了模型的解析解。基于文中模型,对比了非等温条件和等温条件下有机污染物在黏土衬垫中的运移。非等温条件的黏土衬垫底部浓度比等温条件下的浓度大,非等温模型比等温模型更偏于安全。解析解可为实验数据的拟合及分析、衬垫的设计提供参考。进一步的参数分析表明:增大阻滞因子Rd和黏土厚度L可以减小黏土衬垫底部由温度梯度引起的有机污染物浓度。     

非等温条件下有机污染物在黏土衬垫中的扩散分析

填埋场中垃圾降解作用产生的热量在黏土衬垫中的运移会形成温度场,温度场会对污染物的运移产生影响。文章在多孔介质热力学理论和污染物运移理论相结合的基础上,主要考虑温度场对污染物扩散的影响,建立了非等温条件下有机污染物在黏土衬垫中的一维扩散模型,并采用分离变量法得到了模型的解析解。基于文中模型,对比了非等温条件和等温条件下有机污染物在黏土衬垫中的运移。非等温条件的黏土衬垫底部浓度比等温条件下的浓度大,非等温模型比等温模型更偏于安全。解析解可为实验数据的拟合及分析、衬垫的设计提供参考。进一步的参数分析表明:增大阻滞因子Rd和黏土厚度L可以减小黏土衬垫底部由温度梯度引起的有机污染物浓度。     

层状非均质结构包气带入渗过程单相流与两相流数值模拟对比研究

包气带水分入渗过程受多种因素的影响。定量研究层状非均质岩性结构和入渗速率对其影响,有助于解决根据不同条件选择单相流模型或水气二相流模型模拟包气带水分入渗过程的问题。结合填埋场等场地地层条件及污废水入渗特征,分别建立了“上细下粗”和“上粗下细”包气带层状非均质岩性结构水分入渗单相流和水气二相流模型,探讨不同层状非均质岩性结构条件下模型的适用性。在“上粗下细”岩性结构模型基础上,进一步探究入渗速率对水气两相运移结果的影响。基于论文模型研究表明:（1）在包气带岩性结构为“上细下粗”的条件下,气相的影响基本可以忽略,可直接采用单相流模型对包气带水分运移进行模拟;在“上粗下细”岩性结构和本次模型设定的底部压力保持不变及污废水泄漏前场地未接受降水入渗补给等条件下,当包气带上下层介质渗透率比值大于16倍时,气相会对水相运移产生明显影响,且下层介质渗透率越小、上下层介质渗透率比值越大,单相流与两相流的运移结果差别越大,需要采用水气二相流模型模拟包气带水分运移。（2）在包气带“上粗下细”岩性结构条件下,入渗速率越大,气相对水流入渗的阻滞作用越明显,此时包气带水分运移模拟应采用水气二相流模型。     

LNAPL在二维含水介质中运移规律的研究

为研究LNAPL(轻质非水相液体)污染物在均匀含水介质中的运移规律,本文建立微元模型推导出LNAPL在均匀含水介质中的对流--弥散方程。针对不同的介质模型,对二维对流--弥散方程运用差分法进行数值模拟,求其数值解。分析LNAPL污染物的运移特征,得出在均匀介质和非均匀层状介质中运移规律:LNAPL污染物在含水介质中的运移规律遵循对流--弥散方程,介质的弥散系数是影响LNAPL污染物运移的主要因素。