作者： 郭楷文 1 ， 张坤勇 1,2 ， 艾英钵 1,2 ， 李 威 3 ，周 江 3
单位：
1．河海大学岩土工程科学研究所， 南京210024
２．岩土力学与堤坝工程教育部重点实验室（河海大学）， 南京210024
３．南京水利科学研究院岩土工程研究所， 南京210029
４．中铁四局集团有限公司市政工程分公司， 合肥230023
摘要: 国内垃圾填埋场大多已处于长期服役状态，依据实时监测数据动态评估垃圾堆体稳定性，计算当日垃圾合理进场量与堆填位置，对填埋场日常安全运营具有重要的指导意义。 以国内某长期运营垃圾填埋场为工程背景，采用极限平衡方法，计算分析了渗沥液存在形式、液面高度、垃圾堆填高度及位置对垃圾堆体边坡稳定性的影响，得到警戒水位高程、极限垃圾堆填高度和日均垃圾进场量。 结果表明：渗沥液液面高程、新进场垃圾堆填高度双因素与安全系数间呈线性关系。 应合理布置新进场垃圾堆填位置，不得过于靠近坡底，并及时推平新进场垃圾以控制堆填高度。 为保证填埋场的长期稳定性，应及时检查填埋场的导排系统，确保其畅通。
城市垃圾指城市区域内固体废弃物的混合体，因城市垃圾的含水率较高，堆填后其中的液体会形成渗沥液（包含水、废气和有机废液）。随着我国快速城市化，城市生活垃圾产出量急剧增加，填埋法是当前处理生活垃圾的主要方法 [1] 。近年来，我国生活垃圾的处置方式从填埋向焚烧转变，但是由于垃圾焚烧设施处于建设阶段，大部分生活垃圾仍须进行填埋处理 [2] 。垃圾填埋达到一定高度后，会覆上一定厚度的黏土作为覆盖层而后继续填埋垃圾。由于中间覆盖层的存在，渗沥液会以滞水的形式存在于垃圾体内，渗沥液产生的压力最大为垃圾的层高。部分填埋场采用防渗膜作为中间覆盖，当防渗膜破裂，各层垃圾间的渗沥液可能相互连通，形成整体渗流，渗沥液产生的压力比有覆盖层间隔时大很多。渗沥液的存在形式不同，对垃圾堆体稳定性影响也不同， 容易引发边坡失稳问题 点 [3-4] 危害周围环境 [5] 。
Koerner等调查了1０个失稳填埋场，发现渗沥液液面升高是填埋场边坡失稳的主要原因 [6] 。陈云敏等通过计算发现填埋体中渗沥液饱和度越大、渗沥液液面越高，填埋体安全系数就越小 [7] 。冯世进利用有限元法分析了填埋体坡度、高度及渗沥液液面高程变化对填埋体安全系数的影响 [8] 。叶建忠等将填埋场所涉及的稳定问题分为五类并进行了分析 [9] 。沈磊在前人工作的基础上总结了三种渗沥液液面形成条件，并就滞液液面高程对边坡稳定性的影响进行了详尽分析 [10] 。邱贞波等研究了不同垃圾堆填高度和不同渗沥液液面条件下垃圾坝的稳定性，发现垃圾堆体液面高程对垃圾坝的稳定性影响较大 [11] 。黄志亮等建立了扩容后填埋场堆体边坡的计算模型，分析评价了三种不同工况边坡抗滑稳定性的影响 [12] 。杨逾等对液面高程、垃圾土的摩擦角、垃圾土的黏聚力进行敏感性分析，结果表明液面高程变化对安全系数的敏感度相对较大 [13] 。
某填埋场于1994年建成并投入使用，设计总库容约462 万 m 3 。经过十几年的垃圾堆填和城区的不断扩大，一期工程已不能满足需要，于2007、2014年两次在原工程顶部进行扩建，设计总库容达到900 万 m 3 ，如图１所示。
根据现场资料，选取典型断面如图 ２ 所示，典型断面走势较陡，附近存在隧道、污水池、水位管。
通过现场查勘确定最后一层黏土覆盖层厚度为43m，渗沥液液面距地表12.7m，垃圾体内渗沥液液面高程由定期对液面管内的液面测量所得。勘察报告中对不同年代垃圾土进行了分层分组试验，建议以各年份垃圾参数的平均值作为计算参数，密度为0.82g/c m 3 ，不排水剪切指标为19.33kpa、17.54°，试验结果与文献 [14] 结果也较为符合。新进场未及时推平的垃圾荷载、用于填埋垃圾的行车道板荷载、机械荷载等是坡顶荷载的主要组成部分。新进场垃圾堆填范围近似一长方形，长、宽约为38，23ｍ。生活垃圾填埋场使用平面为8.0ｍ×1.5ｍ、每块约３ｔ 的钢板路基箱铺在垃圾上作为行车道板，提高垃圾承载力，以此解决车辆进入填埋场的问题。实际用于堆放垃圾的范围长、宽约为30，15ｍ，荷载分布如图3所示，H代表新进场垃圾堆填高度。
极限平衡理论是经典的边坡稳定性分析方法 [15] ，以摩尔－库仑准则为基础，将潜在滑移面内的土体划分为若干土条，通过土条的平衡条件建立静力平衡方程，从而求解安全系数。若土坡沿某一滑面滑动，则在该滑面上的任一点的土体达到极限平衡状态，如式（１）所示：
极限 平 衡 分 析 常 用 的 方 法 有 瑞 典 条 分 法、Bishop法、Janbu法、Spencer 法、Morgenstern-Price法等，这些方法不同之处主要在于对平衡条件、条间力以及条间力关系的假定不同。目前关于各方法之间的对比分析研究已经相当充分 [16-19] ，计算采用较严密的Morgenstern-Price 法 [20] ，该方法考虑了土条间法向力和切向力，假定两相邻土条法向条间力与切向条间力之间存在函数关系，该方法同时满足整体力和力矩的平衡条件，考虑坡内静水和渗流对边坡稳定影响，适用于任意形状的滑动面。
安全系数通常有两种定义方法：一种是通过加大外力以达到极限平衡状态，具有超载系数的性质；另一种是通过降低材料的强度达到极限平衡状态，具有强度储备系数的性质。 Bishop 定义了如下的安全系数具有强度储备的性质，即安全系数为整个滑动面上的抗剪强度与达到平衡所须的剪应力之比：
式中：τｆ 为滑面上的抗剪强度；τ 为达到平衡所须的剪应力。
采用在Geo-Studio商业软件中的SLOPE/W模块下选择Morgenstern-Price法作为计算方法，条间力分布函数选择半正弦函数。 建立几何模型后，定义垃圾土参数并设置水位条件和滑移面进出口范围，对垃圾填埋场典型断面不同工况下的稳定性展开计算，使用滑动面自动搜索功能搜寻最危险滑动面（ 边坡稳定性系数最小的滑动面）。 根据CJJ176-2012《生活垃圾卫生填埋场岩土工程技术规范》规定，该工程安全等级为一级，边坡稳定安全系数应大于1.35。
2.1.1 计算工况
对典型断面考虑不同工况下的稳定性，各工况如表1所示。
工况1，不考虑渗沥液的影响，相当于堆填体内无地下水位影响情形；工况2，采用实测渗沥液液面高程，考虑2012年前垃圾填埋过程中的黏土覆盖层，土层间考虑黏性土覆盖层的滞液作用；工况3，考虑在极端不利条件下，防渗膜及黏土覆盖层滞液作用未发挥，每层垃圾土间渗沥液贯通，在填埋体内形成整体稳定渗流的工况。
2.1.2 计算结果
分别开展了圆弧滑动和非圆弧滑动计算分析，限于篇幅，只展示典型堆坡非圆弧滑动稳定性计算结果，如图４所示。
整理不同渗沥液赋存形式的计算结果如图5所示，典型断面在工况1、2下，堆填体内无论是非圆弧滑动破坏或圆弧滑动破坏，安全系数均大于1.35，满足现行标准的要求；对比工况1、2，堆填体内从无液到滞液，安全系数降低，可见堆填体内渗沥液的存在对填埋体的稳定性有重要影响；工况３在两种滑动方式下，计算安全系数均小于1.35，不满足 CJJ176-2012 的要求，说明填埋体极有可能出现失稳状况； 对比工况2与3，渗沥液的状态从滞液到形成整体渗流，安全系数再次降低，可见渗沥液存在的状态对堆填体稳定性有重要影响。
随着渗沥液液面高程的提高，孔压增大，滑动面上的垃圾强度指标降低，造成垃圾体稳定性的降低，故有必要考虑渗沥液液面高程变化对堆体稳定性的影响。 因此须考虑典型断面在工况2的基础上液面高程最大变化2ｍ 范围内的稳定性。 给出拟合曲线和警戒液面高程计算值随液面高程变化规律如图6所示。
堆填体内无论是非圆弧滑动破坏还是圆弧滑动破坏，随着渗沥液液面高程的上升，堆填体安全系数减小，稳定性降低。 根据最危险断面的计算结果，建议在运营中保证渗沥液液面高程不能比现有液面高程高2.4ｍ。
垃圾堆体的稳定性除受渗沥液影响之外，堆体上部的荷载对堆体稳定性也有重要影响。 堆体上部荷载主要包括新进场垃圾荷载、机械荷载、用于行车的道板荷载等。 考虑新进场垃圾的不断积累，随着堆填高度的增加，堆体稳定性将发生变化。对典型断面，考虑在现在断面形状及渗沥液液面高程的情况基础上，按前所述荷载填加方法及添加范围，分别计算垃圾堆填高度达5，10，15，20ｍ时堆体稳定性，并分析随堆填高度增加，堆体稳定性变化情况。 限于篇幅，只展示非圆弧滑动下垃圾堆填高度为5ｍ 的稳定性计算结果，如图 ７所示。 拟合曲线和警戒堆填高度计算值如图 ８所示。
从图 ８ 可见：堆填体内无论是非圆弧滑动破坏或圆弧滑动破坏，断面的安全系数随着垃圾堆填高度的增加而减小，稳定性降低，计算得到垃圾堆填高度警戒值为8ｍ。 建议堆高根据现场实际情况宜控制在4ｍ，最高不超过6ｍ，按照现控制标准计算（控制高度6ｍ，堆填面积15ｍ×30ｍ，垃圾密度0.7t／ m 3 ），每日进场垃圾量应控制在1900ｔ 以下。
前文给出了单因素影响下警戒液面高程值及警戒垃圾堆填高度值。当填埋场遭遇强降雨等不利因素引起渗沥液液面高程短时间内显著变化时，堆填过程中，考虑堆填垃圾高度变化的同时，须考虑渗沥液液面高程变化的影响。此外，通过分析垃圾堆填高度及渗沥液液面高程共同作用对堆体稳定性的影响，可以确定某一堆填高度下允许液面高程及某一液面高程下允许垃圾堆填高度。故下面对断面1进行双因素共同作用下堆体稳定性影响分析。
对典型断面分别计算垃圾堆填高度为 5，10，15ｍ 时，渗沥液液面高程分别为-2，-1，1，2ｍ 条件下堆体的安全系数。 计算结果如表 ２ 所示。
可见，安全系数与渗沥液液面高程及垃圾堆填高度间呈线性关系，双因素共同作用下，考虑安全系数与两者之间为线性关系，对表2计算结果进行线性回归分析：
非 圆弧滑 动下 ， 回归方程为 ： Ｆ ｓ ＝ -0.0339H-0.0968h+1.6027， 判定系数 Ｒ 2 ＝0.96 。
圆 弧滑动下 ， 回归方程为 ： Ｆ ｓ=- 0.028 1H-0.0997h＋1.7395 ， 判定系数 Ｒ 2 ＝ 0.98。
依据回归方程，给出任意垃圾堆填高度和液面高程就可以得到相应的安全系数，见图 9，以此评估施工的合理性。
新进场垃圾堆填位置对堆体稳定性亦有影响。 通过将垃圾荷载向坡体后方移动以反映堆填位置改变，分别计算荷载向后移动-6，-4，-2，0.5，10，15，20，25m情况下堆体的稳定性，计算结果如图10、11所示。
图10反映了垃圾堆填高度为 ５，１０，１５ ｍ 时，不同渗沥液液面高程下堆填位置后移对安全系数的影 响。 从中可看出：相同堆填高度位置下，安全系数随渗沥液液面高程的增大而减小。相同渗沥液液面高程下，堆填位置后移10m范围内安全系数略有下降，从后移10m开始，随着后移距离的增加，安全系数逐渐增大。
图11分别对应渗沥液液面高程为-1，1，1ｍ时，不同垃圾堆填高度下堆填位置后移对安全系数的影响。可见，同一垃圾堆填高度下，在对应的渗沥液液面高程下，堆填位置向前移动4ｍ 范围内安全系数有所降低，当前移距离大于4ｍ 时降低明显，故堆填范围最多前移 ４ ｍ，且当堆填高度过高或在对应的渗沥液液面高程上升明显时严禁堆填位置前移，应保持堆填位置不变或适当后移。向后移10ｍ范围内安全系数略有下降，从后移10ｍ 开始安全系数逐渐增大，不同堆填高度间安全系数差值减小，当后移距离达25ｍ 时，不同垃圾堆填高度下堆体安全系数相同， 说明此时上部荷载对坡体最小安全系数无影响，可不考虑上部堆载对坡体稳定性的影响。
根据以上计算分析，可能造成该填埋场安全问题主要原因为：１） 垃圾堆填速度快，填埋作业面有限，对下部垃圾堆体稳定性产生不利影响；２） 新填埋区渗沥液液面较高，由于垃圾本身密度小，渗沥液的浮力及渗透力对垃圾堆体稳定性产生不利影响；３）新填埋区垃圾堆填方式与老垃圾堆填方式不同，原来垃圾堆填过程中，每堆填 ８ ～ １０ ｍ 高就会设置一中间覆盖黏土层，可以隔开上、下两层垃圾中的渗沥液；而现在采用行车道板填埋垃圾，利用防渗膜作为中间覆盖，致使上、下层垃圾中的渗沥液容易联通，形成连续渗流，对垃圾堆体稳定性不利。
根据对不同因素影响下堆体稳定性分析结果，针对堆填现场可能出现的安全隐患，提出以下处置措施建议：
１）通过现场监测和理论分析，抽排竖井迫降水位是最直接、有效的应急抢险措施。 填埋过程中应定期对堆体内渗沥液液面高程变化情况进行监测，当液面接近警戒位时，应进行应急抢险：ａ．及时采取应急措施降低液面高程，如采用小口径抽排竖井快速迫降渗沥液液面高程。ｂ．停止垃圾进场，将坡顶堆填垃圾适当清除。
２）堆填过程中，新进场垃圾应及时推平，堆填高度不宜过高，最大填高控制在 ６ ｍ 以下，日均垃圾进场量控制在1900ｔ。 雨季垃圾堆填过程中，进场垃圾含水率很高，抗剪强度比老垃圾低且新填垃圾无凝聚力，呈现散粒体材料特性，为防止降雨导致局部堆体渗沥液液面升高，出现局部表层溜滑，应密切监测液面高程数据。 当液面比现状液面高程高出1.0ｍ 时，应将堆填高度控制在4ｍ 以下，对应日均垃圾进场量控制在1300ｔ 以下；当液面比现状液面高程高出1.8ｍ 时，应将堆填高度控制在2ｍ 以下，对应日均垃圾进场量控制在700ｔ以下。
３）堆填过程中，垃圾堆放严格按照设计断面施工，临时堆放不能超过设计堆填界限 ４ ｍ，且当堆填高度过高或水位上升明显时，需要严格控制垃圾堆填范围严格按照设计堆填位置作业，确保堆填位置不变或适当后移。
对某填埋场垃圾堆体稳定性分别考虑渗沥液赋存形式和高程、新进场垃圾堆填高度及位置，开展了极限平衡稳定分析，根据稳定分析结果，给出以下结论和建议：
１）垃圾体内的渗沥液对填埋场安全性有较大影响。 当渗沥液导排系统有效工作时，垃圾体内渗沥液及时排除，垃圾体稳定性能得到有效控制；坡体内渗沥液以分层滞液形式存在时的堆体稳定性，高于渗沥液在坡体内贯通后形成整体渗流时的稳定性；对于极端情况，原垃圾土层间黏性土及防渗膜的滞液作用受各种不利因素影响失效，垃圾土层内渗沥液贯通，原以滞液形式存在的渗沥液在坡体内形成整体渗流，垃圾体稳定性显著降低，可能导致边坡失稳。
２）渗沥液液面与坡体稳定性之间存在明显规律。 在大量计算基础上给出了不同断面坡体安全系数与渗沥液液面高程之间关系，确定了不同断面对应的警戒液面高程；根据最危险断面计算结果，建议在运营中控制渗沥液液面高程不得比现有液面高程高2.4ｍ。
３）堆体安全系数与垃圾堆填高度间存在良好的线性关系。 根据两者间关系确定不同断面对应的警戒垃圾堆填高度，最终计算警戒值为8ｍ；建议堆高根据现场实际情况宜控制在4ｍ，最高不超过6ｍ，每日进场垃圾量应控制在1900ｔ 以下。
４）考虑渗沥液液面高程和垃圾堆填高度共同作用下堆体的安全稳定性变化规律。 对三者关系进行回归分析，根据回归式和安全系数，可以确定某一堆填高度下允许渗沥液液面高程或某一渗沥液液面高程下允许垃圾堆填高度。
５）分析了堆填位置对堆体稳定性的影响。 堆填位置前移堆体稳定性降低，堆填位置后移一定距离后，堆填垃圾对堆体稳定性减小直到可忽略堆载对堆体稳定性的影响；根据计算结果给出在现有渗沥液液面高程条件下，堆填前移位置允许变化范围不宜超过现有界限4ｍ。
参考文献
［１］ 管仁秋． 城市固体废弃物填埋体边坡稳定分析及工程控制措施［Ｄ］．杭州：浙江大学，2010．
［２］ 刘彦辉． 城市生活垃圾焚烧处理污染防治方法探讨［Ｊ］． 中国环保产业，2015（5）：47-51．
［３］ 詹良通， 罗小勇， 管仁秋， 等． 某垃圾填埋场污泥坑外涌及其引发下游堆体失稳机理［ Ｊ］． 岩土工程学报，2013，35（ ７）：1189－1196．
［４］ ＪＡＨＡＮＦＡＲ Ｍ Ａ． Ｌａｎｄｆｉｌｌ ｓｌｏｐｅ ｓｔａｂｉｌｉｔｙ ｒｉｓｋ ａｓｓｅｓｓｍｅｎｔ ［ Ｄ］．Ｇｕｅｌｐｈ：Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｇｕｅｌｐｈ，2014．
［５］ 白丽荣， 龚航远， 徐敏等． 太原市某垃圾填埋场渗滤液及周边环境重金属污染及健康风险评估 ［ Ｊ］． 生态毒理学报，２０２１， １６（４）：３１３－３２２．
［６］ ＫＯＥＲＮＥＲ Ｒ Ｍ， ＳＯＯＮＧ Ｔ Ｙ． Ｌｅａｃｈａｔｅ ｉｎ ｌａｎｄｆｉｌｌｓ：ｔｈｅ ｓｔａｂｉｌｉｔｙ ｉｓｓｕｅｓ［Ｊ］． Ｇｅｏｔｅｘｔｉｌｅｓ ａｎｄ Ｇｅｏｍｅｍｂｒａｎｅｓ， ２０００， １８（ ５）：２９３ －309
［７］ 陈云敏， 王立忠， 胡亚元， 等． 城市固体垃圾填埋场边坡稳定分析［Ｊ］． 土木工程学报， ２０００（３）：９２－９７．
［８］ 冯世进．城市固体废弃物静动力强度特性及填埋场的稳定性分析［Ｄ］．杭州：浙江大学，２００５．
［９］ 叶建忠， 周健． 垃圾填埋场设计中的稳定分析问题［ Ｊ］． 工业建筑， ２００５，３５（增刊 １）：７３５－７３８，７４１．
［１０］ 沈磊． 城市固体废弃物填埋场渗滤液水位及边坡稳定分析［Ｄ］． 杭州：浙江大学， ２０１１．
［１１］ 邱贞波， 杲加俊，潘倩，等． 某生活垃圾填埋场垃圾坝稳定和渗流分析［Ｊ］．水利规划与设计， ２０２１（９）：１２６－１３１．
［１２］黄志亮，郭晓静，蔡小平，等．含污泥坑生活垃圾填埋场竖向扩容工程稳定沉分析［Ｊ］．环境工程学报， ２０２２，１６（２）：６３８－６４４．
［１３］ 杨逾， 何胤实， 陈锋， 等． 渗沥液与载荷作用下陈旧型垃圾填埋场坝体稳定性［Ｊ］． 辽宁工程技术大学学报（自然科学版），２０２２， ４１（１）：７－１５．
［１４］ 刘晓立． 降雨渗流作用下垃圾填埋场边坡稳定分析［Ｄ］． 天津：天津大学， ２００６．
［１５］ 刘波， 黄卫． 基于极限平衡理论的边坡稳定性分析方法对比研究［Ｊ］． 工业建筑， ２０１４， ４４（增刊 １）：７１５－７２１，７１４．
［１６］ 陈祖煜． 土质边坡稳定分析：原理·方法·程序［ Ｍ］． 北京：中国水利水电出版社， ２００３．
［１７］ 郑颖人， 杨明成． 边坡稳定安全系数求解格式的分类统一［Ｊ］．岩石力学与工程学报， ２００４（１６）：２８３６－２８４１．
［１８］ 王肇慧， 肖盛燮， 刘文方． 边坡稳定性计算方法的对比分析［Ｊ］． 重庆交通学院学报， ２００５（６）：９９－１０３．
［１９］ 朱大勇， 李焯芬， 黄茂松， 等． 对 ３ 种著名边坡稳定性计算方法的改进［Ｊ］． 岩石力学与工程学报， ２００５（２）：１８３－１９４．
［２０］ ＭＯＲＧＥＮＳＴＥＲＮ Ｎ Ｕ， ＰＲＩＣＥ Ｖ Ｅ． Ｔｈｅ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｓｔａｂｉｌｉｔｙ ｏｆ
ｇｅｎｅｒａｌ ｓｌｉｐ ｓｕｒｆａｃｅｓ［Ｊ］． Ｇｅｏｔｅｃｈｎｉｑｕｅ， １９６５， １５（１）：７９－９３．
