分层系统中无机污染物通过电动力屏障的二维运移：一种半解析方法

《Journal of Environmental Management》：Two-dimensional transport of inorganic contaminants through electro-kinetic barriers in a layered system: A semi-analytical method

　　

摘要指出，电动力屏障作为一种新型主动污染控制技术，在抑制无机污染物离子向邻近土壤和地下水迁移方面具有重要潜力。针对分层地下系统中污染源空间分布不均匀、传统一维模型难以描述二维羽流扩散与电动力屏障防控性能等问题，本研究建立了考虑非均匀污染源边界的二维半解析污染物迁移模型，模型耦合了电渗流、电动迁移以及平流—弥散—吸附过程，并基于拉普拉斯变换、有限余弦变换和传递矩阵法推导出拉普拉斯域的半解析解，再通过Talbot数值反演获得时域解。模型通过与已有的一维解析解、实验室土柱试验数据以及COMSOL数值模拟结果进行对比验证，进一步分析了模型维度、污染源空间分布及关键设计参数对三层系统中污染物迁移行为的影响。结果表明，现有一维模型往往高估污染物迁移距离和羽流轴线浓度，同时低估侧向污染范围；与固定浓度入口边界相比，非均匀分布污染源会延后突破时间，说明现有设计方法可能过于保守。基于上述结论，研究提出了电动力屏障厚度与电压梯度协同优化的设计策略，为工程应用中的设计与参数优化提供了实用工具。

引言部分首先指出，随着全球工业化进程加快，工业场地有毒化学品泄漏和尾矿坝渗滤液排放已造成严重的土壤与地下水污染，对生态系统和人体健康构成重大威胁。在传统工程实践中，垂直截水墙和复合防渗系统被广泛用于防止或减缓污染物迁移，但这些被动屏障通常受限于高施工成本、场地可达性差以及材料老化、不均匀沉降等长期渗漏风险。相比之下，电动力屏障作为一种主动控制技术，能够克服物理隔离限制，为污染控制工程设计提供新思路。

电动力系统中的物质传输过程主要包括电渗流、电动迁移、平流和扩散。当直流电场作用于埋设于土壤中的电极时，孔隙水、溶解离子和带电颗粒会在电场作用下定向迁移。该技术最早由Esrig用于低渗透黏土的脱水与固结，随后被拓展用于污染场地修复，被称为电动力修复。基于这一原理，电动力屏障能够阻止重金属离子向周围土壤继续迁移。相关研究还表明，在连续电流条件下施加较低电位梯度即可有效抑制镉、铜等重金属污染物的迁移，说明电动力屏障在矿区土壤与地下水污染防控中具有可行性和工程价值。

为定量刻画电动力系统中涉及多场耦合的复杂物理化学过程，已有研究建立了基于非平衡热力学的理论框架，描述水、电和离子在水力、电学和化学梯度共同作用下的耦合传输，并进一步引入土体体积压缩性以发展多场耦合有限元模型。虽然高精度数值模型能够处理复杂边界条件并给出细致预测，但其依赖离散网格、物理机制表达不够直观，限制了在方案初选和快速评估中的应用。解析解虽然做了几何简化，但能够清晰体现参数依赖关系，并用于敏感性分析和数值模拟验证，因此在工程设计中具有不可替代的作用。然而，已有电动力屏障研究仍较薄弱，尤其是多数模型采用半无限均匀介质和恒定浓度边界假设，难以反映实际工程中有限范围、非均匀分布的污染源特征，容易高估进入屏障体系的污染物通量，导致设计偏于保守。此外，现场常存在明显垂直分层，若将多层体系简化为均质单层介质，则难以准确捕捉层间界面的浓度连续与通量交换。已有的一维模型还忽略了横向水力弥散，进一步削弱了对非均匀污染源条件下污染物浓度场和迁移特性的描述能力。因此，建立能够同时反映非均匀污染源与分层地基地质条件的二维污染物迁移理论模型，对于阐明电动力屏障的运移机理和评估其防控性能具有重要意义。

在电极布置方面，研究指出，传统竖向电极虽施工简便，但阳极附近易产生开裂、界面阻力增加，系统效率受限；横向电极布置更有利于抑制裂隙扩展并提高电场均匀性。电动力土工合成材料的发展进一步增强了横向布置的工程可行性，其具备柔性好、铺设便捷、电化学稳定性强等优势，可避免传统金属电极腐蚀后释放金属离子造成的二次污染，因此更适用于长期稳定运行的大型分层地下工程。

本研究的三项主要目标包括：一是建立可同时考虑过渡层、电动力屏障层和原土层分层异质性以及非均匀污染源边界的二维半解析污染物迁移模型；二是通过与一维解析解、实验室土柱试验及数值模拟对比验证模型可靠性，并量化传统一维模型在非均匀污染羽流条件下的预测偏差；三是分析电压梯度、电动力屏障厚度及污染源空间分布对突破时间和稳态浓度的影响，进而提出电压梯度与屏障厚度的协同优化策略。

基本假设部分说明，研究模型采用二维笛卡尔坐标系，系统由过渡层、电动力屏障层和原土层组成，电极采用水平条状阴阳极阵列布置。研究重点考察污染物在多层介质中的二维迁移，并考虑水力平流、水力弥散、线性吸附、电渗流和电动迁移的耦合作用。由于文本中的图件和部分参数说明被截断，无法完整复现所有符号定义与边界条件细节，但核心物理框架和数学处理方式仍可明确把握。

模型验证部分采用三种情形进行校核：一维同层体系解析解对比、一维土柱试验数据对比以及二维数值模拟对比。文中虽未给出完整数值结果，但从描述可判断，验证目的旨在确认该半解析模型在分层介质和非均匀边界条件下能够再现主要运移趋势、突破行为及浓度分布规律。

参数分析部分强调，工业场地和尾矿库中低渗透土层或防渗屏障上的垂直水力梯度具有显著场地依赖性，可能与上覆积水、渗滤液下渗及地下水头差等因素有关。现场资料显示水力梯度通常在1到10之间，表明尾矿堆积条件下局部水力梯度可能偏高。分析重点包括模型维度差异、污染源空间分布特征以及电压梯度、屏障厚度等设计参数对污染物迁移的影响。

电压梯度与电动力屏障厚度集成设计部分指出，工程设计本质上是一个带约束的优化问题，既要满足预定防控指标，也要兼顾场地空间限制和长期运行成本。以往研究多集中于单一参数敏感性分析，忽略了设计变量之间的协同效应。本部分借助二维半解析模型的计算结果，探索电压梯度与屏障厚度的组合优化方法，以提高工程设计的经济性、合理性和实际可实施性。

局限性部分承认，所提出的二维解析模型建立在若干基本假设之上，虽便于数学推导和工程计算，但也带来一定限制。例如，将分层土壤简化为饱和、均匀、各向同性多孔介质，无法完全反映天然沉积层中的明显各向异性、非均质性以及浅层局部非饱和特征。这些限制意味着该模型更适合用于初步设计、方案筛选和参数趋势判断，而非替代针对复杂现场条件的高精度三维数值模拟。

结论部分总结，本研究建立了用于研究分层地下系统中多物理场耦合污染物迁移的二维半解析理论模型，模型考虑了空间非均匀污染源以及水力平流、水力弥散、线性吸附、电渗流和离子电动迁移的耦合作用。研究通过与一维解析解、实验室土柱试验和二维数值模拟的对比完成验证，并进一步讨论了维度差异、源项分布和设计参数对迁移行为的影响。研究结果表明，传统一维方法在复杂源项条件下存在系统性偏差，而采用二维视角和非均匀源项假设能够获得更贴近实际工程条件的预测结果。最终，研究提出了一种兼顾防控性能与工程经济性的电压梯度—屏障厚度协同优化思路，可为电动力屏障在污染场地阻隔与地下水保护工程中的应用提供理论支撑和设计参考。

总体来看，该文献的核心贡献在于突破了一维、均质、恒定源边界的传统简化，提出了适用于分层地下结构和非均匀污染源条件的二维半解析模型，并在工程优化层面给出了电压梯度与屏障厚度协同设计的方法，兼具理论深度和工程应用价值。