破损土工膜对垃圾填埋场有多大影响
城市固废：日益严峻的全球挑战
在全球城市化进程不断加快的大背景下，城市固体废物（Municipal Solid Waste, MSW）的产生量也在持续攀升，这已经成为了一个亟待解决的全球性问题。根据相关统计数据显示，全球每年产生的固废数量惊人，超过 20 亿吨 。而我国在其中占比约 16%，年产量高达 3.2 亿吨。这些城市固体废物涵盖的范围十分广泛，像日常生活中产生的生活垃圾、厨房产生的厨余垃圾、建筑施工过程中出现的建筑废弃物，以及废旧电池、过期药品等有害垃圾，都包含在其中。
这些种类繁多的固废在自然环境中，受到降雨、温度变化以及微生物活动等多种因素的综合作用，会产生大量的渗滤液。渗滤液中主要含有挥发性有机物、酚类化合物、脂肪酸等有机污染物 ，成分非常复杂。目前，全球范围内处理固废的方式多种多样，但仍有大约 33% 的废物采用简易填埋的方式处理，而非正规填埋的比例约占 25%。
简易填埋这种方式操作简单，通常就是开挖凹坑，然后将垃圾直接填埋进去。然而，这种方式存在着极大的弊端，由于缺乏有效的防渗措施，在降雨等情况下，渗滤液会轻而易举地渗入土体，进而对地下水造成污染。据调查，许多采用简易填埋的区域，周边地下水的水质都受到了不同程度的污染，水中的有害物质含量超标，对当地居民的用水安全构成了严重威胁。
在非正规填埋场中，虽然通常会采取一定的措施，在下伏土层表面铺设由高密度聚乙烯（HighDensity Polyethylene, HDPE）制成的土工膜（Geomembrane, GM），将其作为底部防渗屏障。但是在实际的施工以及后续运行过程中，土工膜却面临着诸多挑战。因为土工膜的厚度相对较薄，在施工时可能会因操作不当受到机械损伤；在长期的使用过程中，还可能会受到化学物质的侵蚀发生化学降解，或者因温度的剧烈变化而破损。一旦土工膜出现破损，渗滤液就会如同找到了 “突破口”，通过破损部位源源不断地渗透至下伏土层和地下水中，导致下伏土层的土质恶化，影响植被生长，地下水的水质也会急剧下降，原本清澈的地下水变得浑浊不堪，充满异味，从而大大削弱了土工膜原本应有的阻隔性能，使得填埋场的污染风险大幅增加。
土工膜：垃圾填埋场的关键防线
在垃圾填埋场的众多设施中，土工膜无疑扮演着至关重要的角色，堪称垃圾填埋场的关键防线 。土工膜通常由聚乙烯、聚氯乙烯等聚合物通过挤出或压延工艺制成，是一种柔性薄膜，其厚度一般在零点几毫米至数毫米之间 。别看它厚度不大，却有着极低的渗透系数，能够有效阻止液体和气体的迁移。在垃圾填埋场中，土工膜就像是一个巨大的 “保护膜”，将垃圾与周边环境隔离开来，特别是在防止渗滤液下渗这一关键环节上，发挥着不可替代的作用。渗滤液中富含高浓度的污染物，一旦让它们渗入土壤和地下水，将会对周边生态环境和人类健康造成难以估量的威胁 。而土工膜凭借其致密的结构，成功地将这些污染物封闭在指定区域内，大大降低了污染扩散的风险。
在实际应用中，土工膜通常与其他材料组合使用，形成更为稳固的多层防渗结构 。比如在典型的防渗系统中，就包含基础层、土工膜层、保护层和排水层 。土工膜作为中间层，承担着主防渗的重任 。其应用形式丰富多样，包括单层衬垫和复合衬垫 。单层衬垫一般适用于对防渗要求相对较低的场合，而复合衬垫则是通过将土工膜与土工布等材料巧妙结合，大大增强了系统的耐久性和抗穿刺能力 。在垃圾场库区，土工膜除了铺设在底部和边坡作为防渗衬垫系统的主防渗层外，还常用于覆盖层，防止雨水进入垃圾体，从而减少渗滤液的产生量 。此外，在垃圾填埋场的临时覆盖和最终封场环节，土工膜也发挥着重要作用，助力实现垃圾场的全生命周期管理 。
然而，在非正规填埋场中，土工膜却面临着诸多严峻的挑战，其中最突出的问题就是容易破损 。由于土工膜厚度较薄，在施工过程中，可能会因为施工人员操作不慎，被一些尖锐的物体如石块、钢筋等划破；也可能因为大型机械设备的碾压、碰撞而受损 。在后续漫长的运行过程中，土工膜还会受到化学物质的侵蚀 。垃圾在分解过程中会产生各种复杂的化学物质，这些物质可能与土工膜发生化学反应，导致土工膜的材料性能下降，出现老化、脆化等现象，从而容易产生破损 。温度的剧烈变化也是一个不可忽视的因素 。在一些地区，昼夜温差大，夏季高温炎热，冬季寒冷低温，土工膜在这样的温度环境下反复热胀冷缩，其内部结构会逐渐受到破坏，进而出现裂缝、孔洞等破损情况 。一旦土工膜破损，原本被阻隔的渗滤液就会找到 “突破口”，通过这些破损部位渗透至下伏土层和地下水中，这不仅会削弱土工膜的阻隔性能，还会引发一系列严重的环境污染问题，如土壤污染、地下水污染等 。所以，研究污染物通过破损土工膜的运移规律就显得尤为重要，它对于准确评估填埋场的环境风险，采取有效的污染防治措施具有关键意义 。
前人探索：污染物运移的研究之路
随着土工膜在垃圾填埋场中广泛应用，针对污染物通过破损土工膜在下伏饱和土层运移的研究也不断深入 。Abbas 构建了污染物通过由破损土工膜和压实黏土层（Compacted Clay Liner, CCL）组成的复合衬垫的运移模型 。在推导污染物运移方程时，他充分考虑了狄利克雷和诺伊曼边界条件 。为了提高计算效率，他采用了数值方法，使得复杂的运移过程能够通过计算机模拟进行分析 。在后续研究中，Abbas 等进一步通过耦合稳态流动与反应扩散−平流方程，对破损区域的三维效应展开分析 。研究结果表明，土工膜的褶皱程度以及质量状况对污染物的渗漏率有着显著影响，同时也会改变膨润土防水毯（Geosynthetic Clay Liner, GCL）中的污染物浓度 。比如在一些实际案例中，土工膜褶皱较多的区域，污染物渗漏率明显高于褶皱较少的区域 。
Rowe 等在研究中发现，当破损土工膜与压实黏土衬垫接触时，一个关键的影响因素浮出水面，那就是破损处的渗漏主要由界面传导率决定，而非压实黏土衬垫的渗透系数 。这一发现为后续研究指明了新的方向，让人们更加关注界面传导率在污染物运移过程中的作用 。Xie 等针对二氯甲烷（Dichloromethane, DCM）在由破损土工膜和压实黏土衬垫组成的复合衬垫系统中的运移规律开展研究 。他们利用拉普拉斯变换法，成功求得有机污染物运移的解析解 。研究指出，渗滤液水头与土工膜的施工质量是影响污染物在复合衬垫底部达到地下水最大允许浓度所需时间的关键因素 。如果渗滤液水头较高，土工膜施工质量又存在问题，那么污染物达到最大允许浓度的时间就会大大缩短，对地下水的污染风险也会随之增加 。Xie 等的进一步研究表明，土工膜破损导致的平流效应及土工膜皱褶长度对污染物通量的影响更加显著 。在一些模拟实验中，当土工膜出现破损且皱褶长度较长时，污染物通量明显增大 。
此外，Sun 等提出了半解析−数值耦合模型 。在求解过程中，他们采用离散化方法和积分变换相结合的方式得到半解析解 。通过这个模型，研究了有机污染物通过具有条带状破损的土工膜后在下伏饱和土层中的二维运移规律 。随后，Li 等建立了考虑土工膜破损沿水平方向任意分布的二维污染物运移模型 。为了求解这个复杂模型的半解析解，他们采用了拉普拉斯变换、傅里叶变换和边界变换方法 。研究表明，土工膜破损可显著加速污染物运移，并对污染物浓度的空间分布产生重要影响 。在实际填埋场中，土工膜破损的位置和程度不同，污染物浓度的空间分布也会呈现出明显差异 。
然而，这些研究均是基于污染物通过破损土工膜在下伏饱和土层中的运移过程展开的 。但部分非正规填埋场建在干旱或半干旱地区，像我国西部地区 。在这些地区，地下水位埋置较深，使得填埋场的下伏土层通常为非饱和土层（Unsaturated Soil Layer, USL） 。相较于饱和土层，非饱和土层的水分分布和吸附特性更加复杂 。饱和土层中水分分布相对均匀，而在非饱和土层中，水分会受到多种因素的影响，在不同深度呈现出不同的分布状态 。非饱和土层的吸附特性也与饱和土层不同，这使得污染物通过破损土工膜后在下伏非饱和土层中的运移规律更具不确定性 。所以，研究污染物在非饱和土层中的运移规律，对于全面了解填埋场污染物的扩散情况，保护这些地区的土壤和地下水环境，具有十分重要的意义 。
非饱和土层：复杂的污染运移舞台
污染物在非饱和土层中的传质与运移，就像是一场被诸多因素共同操控的复杂 “演出”，受到含水率、温度变化、扩散系数和吸附特性等关键因素的协同调控 。这些因素各自发挥作用，共同影响着污染物的扩散行为和传输机制 。比如在含水率方面，当土壤含水率较高时，水分就像是为污染物搭建了更多的 “运输通道”，使得污染物能够更快速地在土层中扩散 。而当含水率较低时，污染物的扩散就会受到阻碍，仿佛陷入了 “泥沼” 。
在温度变化的影响上，温度升高时，分子运动变得更加活跃，污染物的扩散速度也会随之加快 。在一些炎热的地区，夏季高温时段，污染物在非饱和土层中的扩散范围明显比其他季节更广 。扩散系数则直接决定了污染物在非饱和土层中的扩散能力 。不同的土壤类型，其扩散系数也有所不同，这就导致污染物在不同土壤中的扩散速度存在差异 。吸附特性同样至关重要，土壤对污染物的吸附作用，就像是给污染物加上了一个 “束缚” 。当土壤对污染物的吸附能力较强时，污染物就会更多地被吸附在土壤颗粒表面，从而减缓其在土层中的运移速度 。
许多学者也针对这些因素开展了深入研究 。Fityus 等采用稳态非饱和水分分布模型，并结合有限层方法对质量传递方程进行简化，深入分析了污染物在饱和与非饱和土层中的运移差异 。研究表明，扩散系数对体积含水率的依赖性显著影响了污染物在垃圾填埋场下非饱和土层中的运移行为 。在一些垃圾填埋场周边的非饱和土层中，由于扩散系数与体积含水率之间的复杂关系，污染物的运移路径和浓度分布呈现出独特的规律 。
Diaw 等采用非连续有限元法与有限差分法分别处理平流和扩散项，并构建了结合非线性和对流−扩散方程的污染物传输模型，以揭示多孔非饱和土中复杂流场对污染物运移的影响 。通过这个模型，他们发现复杂的流场会使污染物的运移变得更加复杂，可能导致污染物在土层中出现不规则的扩散现象 。Lin 等建立了污染物通过完整土工膜在下伏非饱和土层中运移的解析模型，采用罗宾型边界条件模拟污染物运输，并通过与有限差分解及现有模型的对比验证了解析解的有效性 。此外，研究还分析了阻滞因子对污染物运移的影响，并利用敏感性分析评估了模型参数对污染物分布的调控作用 。通过这些研究，我们对污染物在非饱和土层中的运移机制有了更深入的认识 。
Chen 等研究了污染物通过完整土工膜后在下伏非饱和土层中的运移特性 。该研究假设线性水压分布，采用指数型土层保水曲线与线性非饱和扩散系数模型，将非线性控制方程转化为线性方程，并利用解析解方法进行分析 。研究结果表明，土层类型、水压比（即水压力与静水压力之比）条件及地下水深度对污染物扩散具有重要影响 。在不同类型的土层中，污染物的扩散速度和范围有很大差异；水压比的变化也会导致污染物在非饱和土层中的扩散路径发生改变 。
然而，上述这些研究虽然取得了一定成果，但大多存在一个共同的局限性，那就是主要关注污染物在非饱和土层中的扩散行为，并探讨其运移的关键影响因素 。然而，大多数研究均基于完整土工膜的假设，未考虑土工膜破损导致的对流作用对污染物运移的影响 。在实际的非正规填埋场中，土工膜破损是较为常见的现象 。一旦土工膜破损，渗滤液就会形成对流，这种对流作用会极大地改变污染物的运移路径和速度 。如果忽略了这一关键因素，就可能导致对衬垫系统长期阻隔性能的评估出现偏差，无法准确预测污染物的扩散范围和对周边环境的影响程度 。所以，深入研究土工膜破损情况下污染物在非饱和土层中的运移规律，对于准确评估填埋场的环境风险，制定有效的污染防治措施具有至关重要的意义 。
新模型登场：解析污染运移规律
为了深入探究非正规填埋场渗滤液中的有机污染物，在通过破损土工膜后，在下伏非饱和土层中的运移规律，本文构建了一个十分关键的一维分析模型。这个模型就像是一把 “精密的手术刀”，能够细致地剖析污染物在这一复杂环境中的运动轨迹和变化情况 。
在推导污染物运移行为解析解的过程中，本文采用了一系列精妙的数学方法 。边界转换方法就像是一个 “巧妙的化妆师”，能够将复杂的边界条件进行巧妙的转化，使其更易于后续的分析和计算 。位移定理则如同一个 “神奇的搬运工”，在数学变换中发挥着重要的作用，帮助我们更好地理解和处理污染物的运移过程 。拉普拉斯变换及逆变换方法更是关键，它们就像是一对 “魔法钥匙”，拉普拉斯变换能够将时间域的函数转换到复频域，从而简化复杂的运算；而逆变换则能将复频域的结果再转换回时间域，让我们得到有实际意义的解 。通过这些方法的协同作用，我们成功地推导出了污染物运移行为的解析解 。
为了验证所提出解析解的可靠性，我们进行了一项非常严谨的工作 。首先，将非饱和土层模型退化到饱和土层模型，这就像是将一个复杂的系统逐步简化，以便于我们利用已有的知识和经验进行分析 。然后，利用典型算例进行污染物浓度分布的计算 。在计算过程中，我们对各种参数进行了精确的设定和调整，以确保计算结果的准确性 。最后，将所得结果与现有研究数据进行细致的对比 。通过对比发现，我们的计算结果与现有研究数据具有高度的一致性，这充分验证了所提出解析解的可靠性，就像是给我们的研究成果打上了一个 “可靠” 的标签 。
在揭示污染物运移规律的过程中，我们系统地分析了 3 种关键因素对非饱和土层中污染物浓度分布的影响 。土工膜的褶皱率是一个重要因素，它就像是土工膜的 “皱纹” 。当褶皱率较高时，土工膜的表面积增加，这会使得污染物的接触面积增大，从而影响污染物的运移速度和浓度分布 。在一些实际案例中，我们可以观察到，褶皱较多的土工膜区域，污染物的浓度分布更加不均匀 。土工膜的破损孔洞数也不容忽视，它就像是土工膜上的 “漏洞” 。破损孔洞数越多，意味着污染物有更多的通道可以渗透到下伏土层中，这会显著加快污染物的运移速度，导致污染物浓度在土层中的分布发生明显变化 。在模拟实验中，当破损孔洞数增加时，污染物在短时间内就能扩散到更大的范围 。土体的非饱和指数反映了土体的干湿程度和孔隙结构等特性，它就像是土体的 “健康指标” 。非饱和指数的变化会影响土体的吸附性能和水分传输能力，进而对污染物的运移产生重要影响 。在不同非饱和指数的土体中，污染物的吸附和扩散行为会有很大差异，导致浓度分布也截然不同 。通过对这些因素的深入分析，我们更加全面地揭示了污染物在非饱和土层中的运移规律，为后续的研究和实际应用提供了坚实的理论基础 。
未来展望：守护地下水源
通过本文构建的一维分析模型以及对污染物运移规律的深入研究，我们在理解非正规填埋场渗滤液中有机污染物，通过破损土工膜在下伏非饱和土层的运移机制方面，取得了显著的进展 。研究成果明确了土工膜褶皱率、破损孔洞数以及土体非饱和指数等因素，对污染物浓度分布的重要影响 。这些认识对于准确评估填埋场衬垫系统的长期阻隔性能，具有至关重要的意义 。
在实际应用中，我们能够依据这些研究成果，更加科学地设计填埋场的防渗系统 。通过合理控制土工膜的质量和铺设工艺，减少褶皱和破损的出现，从而有效降低污染物的渗漏风险 。同时，对于已有的填埋场，也可以根据研究结论，对其防渗性能进行更准确的评估，及时发现潜在的污染隐患，并采取相应的修复措施 。这对于保护地下水环境，保障周边居民的用水安全，具有不可估量的价值 。
展望未来，我们的研究仍有广阔的拓展空间 。在模型方面，虽然本文构建的一维分析模型取得了一定的成果，但实际填埋场的情况更为复杂，未来可以考虑建立更复杂的多维模型 。将土工膜的破损形态、分布特征以及非饱和土层的空间变异性等因素，全面纳入模型中，从而更加真实地模拟污染物的运移过程 。在研究内容上，除了关注有机污染物，还可以进一步拓展到其他类型的污染物，如重金属、微生物等 。深入研究它们在非饱和土层中的运移规律，以及与有机污染物之间的相互作用，以实现对填埋场污染问题的全面认识 。
随着科技的不断进步，新的监测技术和分析方法也在不断涌现 。我们可以积极引入这些新技术，如高分辨率的地球物理探测技术、先进的传感器技术等，对填埋场的渗漏情况和污染物运移进行实时、精准的监测 。利用机器学习、人工智能等数据分析方法，对大量的监测数据进行深入挖掘，从而更准确地预测污染物的扩散趋势，为填埋场的环境管理和污染防治提供更有力的支持 。通过不断深入的研究和技术创新，我们有信心为保护地下水资源，创造更加清洁、安全的环境做出更大的贡献 。
