固体废弃物堆填场覆盖层水气耦合运移的计算方法

1766 编辑：中冶有色技术网来源：华中科技大学

2022-11-08 15:50:42

权利要求

1.固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于，包括以下步骤：步骤1，建立固废堆填场单层覆盖层土体中水气运移的模型，并设置模型参数及边界条件；步骤2，计算覆盖层中水分和气体的分布情况；步骤3，评价覆盖层的防渗闭气性能。2.根据权利要求1所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述步骤1中的模型参数包括土体性质参数和气体性质参数。 3.根据权利要求1所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述边界条件包括覆盖层底部水分边界、覆盖层顶部水分边界、覆盖层底部气体边界和覆盖层顶部气体边界。 4.根据权利要求1所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述步骤2中，覆盖层水分分布计算包括 B(T)＝exp(αh 0+T/4)， Z＝αz，T＝αk st/(θ s-θ r)，L＝αl，2λ kcos(λ kL)+sin(λ kL)＝0；式中，q(Z，T)为覆盖层顶部入渗或蒸发速率，入渗时为正值，蒸发时为负值，单位为mm/day；α为土体减饱和系数，单位为m；z为距离覆盖层底部的坐标高度，单位为m；t为时间，单位为day；k s为土体饱和渗透系数，单位为m/s；θ s为土体饱和体积含水率；θ r为土体残余体积含水率；θ w为土体体积含水率；l为覆盖层厚度，单位为m；h 0为覆盖层底部孔隙水压力水头，单位为m；h i为覆盖层初始孔隙水压力水头，(m)；λ k为特征值。 5.根据权利要求1所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述步骤2中覆盖层气体浓度分布计算包括式中，C g为覆盖层中气体浓度，单位为mOl/m 3；Y为由组成的n×1阶矩阵；H为由H nr组成的n×r阶矩阵；T为由T n组成的n×1阶矩阵；G为由G n组成的n×1阶矩阵；θ g为土体体积含气率；D s为气体有效扩散系数，单位为m 2/s；τ为土壤孔隙对气体运移的表观弯曲系数；D g为气体扩散系数，单位为m 2/s；k＝k se αh为水分不饱和渗透系数，单位为m/s；k rg为相对气体渗透系数；μ g为气体动力粘滞系数，单位为Pa·s；k i为土体固有渗透率，单位为m 2；P g为气体压强，单位为Pa；S w为土体饱和度；φ为土体孔隙率；H g为亨利常数；v w为水分渗流速度，单位为m/s；v g为气体对流速度，单位为m/s；C 0为覆盖层底部气体浓度，单位为mol/m 3；C t为覆盖层顶部气体浓度，单位为m/s；F 0为覆盖层底部气体流入速率，单位为mol/m 2/day。 6.根据权利要求5所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述步骤2中覆盖层气体浓度分布计算还包括 k取正整数。 7.根据权利要求5所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述步骤2中覆盖层气体浓度分布计算还包括 k取正整数。 8.根据权利要求1所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述步骤3中评价参数包括覆盖层底部水分渗漏速度、覆盖层底部水分累计渗漏量、覆盖层顶部气体排放速率和覆盖层气体累计排放量。 9.根据权利要求8所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述覆盖层底部水分累计渗漏量计算公式为：式中，Q为覆盖层底部水分渗漏量，单位为mm；t 0为计算时间，单位为day；S f为覆盖层土体储水能力，单位为mm。 10.根据权利要求8所述的固体废弃物堆填场覆盖层水气耦合运移的计算方法，其特征在于：所述覆盖层顶部气体排放速率计算公式为：式中F t为覆盖层顶部气体排放速率，单位为mol/m 2/day。

说明书

固体废弃物堆填场覆盖层水气耦合运移的计算方法

技术领域

本发明涉及环境岩土工程技术领域，具体为一种固体废弃物堆填场覆盖层水气耦合运移的计算方法。

背景技术

填埋是固体废弃物的主要处置方式之一。近年来，随着我国固体废弃物产量的不断增加以及环境保护与治理意识的增强，固废堆填场的污染问题愈发受到重视。堆填场中固体废弃物由于成分复杂，经过一系列生物化学反应，会产生污染性的渗滤液和气体。这些污染物若排放到地下水和大气中，会对固废堆填场周围的生态环境造成严重破坏。覆盖层是堆填场重要的组成部分，就像堆填场的“人工皮肤”。覆盖层不仅隔绝人和动物与废弃物之间的直接接触，同时还需具有防渗和闭气功能，即有效阻止雨水入渗到固废堆体中从而减少渗滤液的产生；防止固废堆体产生的有害气体排放到大气中造成空气污染以及隔绝大气中的氧气入侵到固废堆体中。但是，由于覆盖层直接暴露在复杂多变的大气环境中，须经受长期干湿循环、冷暖和冻融循环的作用，与此同时，堆填场中不断产生气体导致堆填场内气压增大，不断考验着覆盖层的防渗闭气性能。因此，计算堆填场覆盖层中水气运移过程对于评价覆盖层的防渗闭气性能、指导覆盖层设计具有重要意义。

目前，本领域相关技术人员已经做了许多研究，但主要侧重于覆盖层的防渗性能研究，即只考虑水分在覆盖层中的运移，而不考虑填埋气的运移以及水分运移气体运移的耦合作用的影响，这在实际工程应用中会产生较大的误差。因此，本领域需要开发一种堆填场覆盖层水气耦合运移的计算方法。

发明内容

鉴于现有技术中所存在的问题，本发明公开了一种固体废弃物堆填场覆盖层水气耦合运移的计算方法，采用的技术方案是，包括以下步骤：

步骤1，建立固废堆填场单层覆盖层土体中水气运移的模型，并设置模型参数及边界条件；

步骤2，计算覆盖层中水分和气体的分布情况；

步骤3，评价覆盖层的防渗闭气性能。

作为本发明的一种优选技术方案，所述步骤1中的模型参数包括土体性质参数和气体性质参数。

作为本发明的一种优选技术方案，所述边界条件包括覆盖层底部水分边界、覆盖层顶部水分边界、覆盖层底部气体边界和覆盖层顶部气体边界。

优选的，针对水分运移，覆盖层底部水分边界条件为常孔隙水压头边界；覆盖层顶部水分边界条件在降雨期t p时间内为水分入渗边界，在t p时间后为水分蒸发边界。针对气体运移，覆盖层底部气体边界条件可设置为常浓度边界或常通量边界；覆盖层顶部气体边界条件为常浓度边界。

作为本发明的一种优选技术方案，所述步骤2中，覆盖层水分分布计算包括 B(T)＝exp(αh 0+T/4)， Z＝αz，T＝αk st/(θ s-θ r)，L＝αl，2λ kcos(λ kL)+sin(λ kL)＝0；式中，q(Z，T)为覆盖层顶部入渗或蒸发速率,入渗时为正值，蒸发时为负值，单位为mm/day；α为土体减饱和系数，单位为m -1；z为距离覆盖层底部的坐标高度，单位为m；t为时间，单位为day；k s为土体饱和渗透系数，单位为m/s；θ s为土体饱和体积含水率；θ r为土体残余体积含水率；θ w为土体体积含水率；l为覆盖层厚度，单位为m；h 0为覆盖层底部孔隙水压力水头，单位为m；h i为覆盖层初始孔隙水压力水头，单位为m；λ k为特征值。

作为本发明的一种优选技术方案，所述步骤2中覆盖层气体浓度分布计算包括

D s＝τD g，式中，C g为覆盖层中气体浓度，单位为mol/m 3；Y为由组成的n×1阶矩阵；H为由H nr组成的n×r阶矩阵；T为由T n组成的n×1阶矩阵；G为由G n组成的n×1阶矩阵；θ g为土体体积含气率；D s为气体有效扩散系数，单位为m 2/s；τ为土壤孔隙对气体运移的表观弯曲系数；D g为气体扩散系数，单位为m 2/s；k＝k se αh为水分不饱和渗透系数，单位为m/s；k rg为相对气体渗透系数；μ g为气体动力粘滞系数，单位为Pa·s；k i为土体固有渗透率，单位为m 2；P g为气体压强，单位为Pa；S w为土体饱和度；φ为土体孔隙率；H g为亨利常数；v w为水分渗流速度，单位为m/s；v g为气体对流速度，单位为m/s；C 0为覆盖层底部气体浓度，单位为mol/m 3；C t为覆盖层顶部气体浓度，单位为m/s；F 0为覆盖层底部气体流入速率，单位为mol/m 2/day。

作为本发明的一种优选技术方案，所述步骤2中覆盖层气体浓度分布计算还包括

k取正整数。

作为本发明的一种优选技术方案，所述步骤2中覆盖层气体浓度分布计算还包括

k取正整数。

作为本发明的一种优选技术方案，所述步骤3中评价参数包括覆盖层底部水分渗漏速度、覆盖层底部水分累计渗漏量、覆盖层顶部气体排放速率和覆盖层气体累计排放量，根据污染排放标准或要求，当覆盖层底部渗漏量和顶部气体排放超过标准要求时，该覆盖层被击穿，由此得到覆盖层的水分和气体击穿时间；若覆盖层底部渗漏量和顶部气体排放在达到稳定时仍小于排放要求，则表示该覆盖层满足设计要求，防渗闭气性能较好。

作为本发明的一种优选技术方案，所述覆盖层底部水分累计渗漏量计算公式为：式中，Q为覆盖层底部水分渗漏量，单位为mm；t 0为计算时间，单位为day；S f为覆盖层土体储水能力，单位为mm。

作为本发明的一种优选技术方案，所述覆盖层顶部气体排放速率计算公式为：式中F t为覆盖层顶部气体排放速率，单位为mol/m 2/day。

本发明的有益效果：本发明综合考虑了水分入渗、水分蒸发、气体对流、气体扩散、水分与气体运移之间的相互影响，可以模拟不同工况和边界条件下覆盖层中水气耦合运移的过程，计算降雨和蒸发等复杂气候变化对气体运移的影响，根据覆盖层底部累积渗漏量和顶部气体排放速率评估覆盖层的防渗闭气性能；能够准确计算覆盖层中水气耦合运移的过程，为实际堆填场覆盖层设计提供了有价值的依据和数据；解决了不考虑覆盖层中水分和气体运移之间相互耦合作用而导致计算结果误差较大的问题，且简单方便、易于实现。

附图说明

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中，类似的元件或部分一般由类似的附图标记标识。附图中，各元件或部分并不一定按照实际的比例绘制。

图1为本发明计算固废堆填场覆盖层水气耦合运移过程的流程示意图；

图2为本发明固废堆填场覆盖层计算模型图；

图3为本发明覆盖层底部水分累积渗漏量图；

图4为本发明覆盖层顶部气体排放速率图。

具体实施方式

实施例1

如图1至图4所示，本发明公开了一种固体废弃物堆填场覆盖层水气耦合运移的计算方法，采用的技术方案是，包括以下步骤：

步骤1为覆盖层中水分和气体耦合运移模型的建立，主要输入土体性质参数、水分运移相关参数、气体运移相关参数、边界条件；

步骤2为计算覆盖层中水气耦合运移过程，得到孔隙水压力水头分布和气体浓度分布；

步骤3为评价覆盖层防渗闭气性能，主要通过计算覆盖层底部渗漏量和覆盖层顶部气体排放，结合相关标准来评价覆盖层的防渗闭气性能。

步骤1，建立堆填场覆盖层土体中水气运移的模型，并设置模型的参数及边界条件；

具体参数取值和边界条件设置如下：覆盖层土体厚度l＝1m，土体饱和体积含水率θ s＝0.43，土体残余体积含水率θ r＝0.0956，土体饱和渗透系数k s＝5×10 -7m/s，土体减饱和系数α＝0.3m -1；土体固有渗透率k i＝0.2×10 -15m 2；气体动力粘滞系数μ g＝1.1×10 -5Pa·s。底部孔隙水压力水头h 0＝-2.5m；覆盖层初始孔隙水压力水头h i＝-2.5m；q(A，T)为覆盖层顶部入渗或蒸发速率(入渗时为40mm/day，蒸发时为-2mm/day，降雨持续时间t p＝5day)。甲烷亨利常数H g＝0.0316；气体扩散系数D g＝2.137×10 -5m 2/s；覆盖层底部气体浓度C 0＝10mol/m 3；覆盖层顶部气体浓度C t＝0mol/m 3；初始气体浓度为0mol/m 3。