Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

    填埋场好氧修复三维沉降变形特性及加速稳定化分析

    冯世进, 白真白, 郑奇腾

    冯世进, 白真白, 郑奇腾. 填埋场好氧修复三维沉降变形特性及加速稳定化分析[J]. 岩土工程学报, 2021, 43(11): 1976-1985. DOI: 10.11779/CJGE202111003
    引用本文: 冯世进, 白真白, 郑奇腾. 填埋场好氧修复三维沉降变形特性及加速稳定化分析[J]. 岩土工程学报, 2021, 43(11): 1976-1985. DOI: 10.11779/CJGE202111003
    FENG Shi-jin, BAI Zhen-bai, ZHENG Qi-teng. Three-dimensional settlement characteristics and accelerated stabilization of landfills under aerobic remediation[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(11): 1976-1985. DOI: 10.11779/CJGE202111003
    Citation: FENG Shi-jin, BAI Zhen-bai, ZHENG Qi-teng. Three-dimensional settlement characteristics and accelerated stabilization of landfills under aerobic remediation[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(11): 1976-1985. DOI: 10.11779/CJGE202111003

    填埋场好氧修复三维沉降变形特性及加速稳定化分析  English Version

    基金项目: 

    国家自然科学基金杰出青年基金项目 41725012

    国家自然科学基金重点基金项目 41931289

    上海市科委社发领域项目 20dz1203402

    详细信息
      作者简介:

      冯世进(1978— ),男,博士,教授,博士生导师,主要从事环境岩土工程、土动力学的教学与科研工作。E-mail:fsjgly@tongji.edu.cn

    • 中图分类号: TU43

    Three-dimensional settlement characteristics and accelerated stabilization of landfills under aerobic remediation

    • 摘要: 沉降变形是重要的填埋场稳定化指标,对于稳定化进程的评估具有重要意义,然而,由于氧浓度分布的不均匀,曝气作用下填埋场内生化降解和沉降变形高度耦合且极为复杂,已有厌氧填埋场沉降模型难以适用。建立了考虑厌氧-好氧生化降解、液气运移、多组分扩散和弹塑性-降解骨架变形的填埋场多场耦合三维模型,采用有限体积法和开源计算平台OpenFoam进行编程求解,揭示了垃圾填埋场好氧修复的沉降变形特性。结果表明,曝气易导致填埋场产生显著的不均匀沉降和垃圾体-曝气井相对位移,可高达约0.7 m(模型为10 m×10 m×高15 m),最大不均匀沉降出现在曝气中间阶段;依据填埋场土地高度利用的国家标准,好氧修复可降低88%的后期维护运营时间,建议采用填埋场90%降解度作为好氧修复的终止标准,以保证修复工后沉降速率满足标准。
      Abstract: The settlement is an important index for evaluating the stabilization of a landfill. However, the settlement of an aerated landfill is extremely complex and highly coupled with waste biodegradation due to non-homogeneous distribution of oxygen. Thus, the traditional models for settlement of anaerobic landfills are not applicable to aerobic ones. A three-dimensional multi-field couped model for landfills is established considering anaerobic-aerobic biodegradation, liquid-gas migration, multi-component diffusion and elastic-plastic-degradation skeleton deformation. The finite volume method and the open source computing platform OpenFoam are used to solve the model, and the settlement and deformation characteristics of an aerobic landfill are revealed. The results show that the aeration can easily induce significant uneven settlement and relative displacement of waste-aeration well, e.g., reaching up to 0.7 m for a 10 m×10 m×15 m model in this paper, and the largest uneven settlement occurs at the intermediate stage of aeration. According to the China's national standard of highly utilizing landfills, the aerobic remediation can reduce the post-maintenance time of landfills by 88%, and a termination condition of 90% degradation degree is suggested for aerobic remediation of landfills to satisfy the criteria in terms of post-aeration settlement rate.
    • 经典塑性力学中的下限定理由Drucker等[1]于1951提出,下限定理认为,对理想刚塑性材料,任意一个静力许可应力场都构成极限荷载的一个下限解。陈惠发[2]在其1969年出版的专著《极限分析与土体塑性》中,阐述了极限分析方法在土工问题中的应用。1970年,Lysmer[3]基于有限单元法的基本思想和Mohr-Coulomb屈服准则的线性化,建立了利用下限定理求解稳定性问题的线性规划模型。Sloan[4]于1988年基于最速下降主动集法,对地基承载力问题进行了求解,获得了稳定性问题的解答。

      下限定理为岩土工程中的稳定性问题提供了一种特别有用的分析方法,它避开了复杂的弹塑性分析过程,同样获得了问题的解,吸引了众多研究人员的兴趣[5-8],在边坡[9-11]、地下洞室[12-14]、基础[15-17]等各类岩土问题中得到了广泛应用。

      在工程实践中,通常首先需要知道岩土工程的稳定性,当岩土工程的稳定性不满足要求时,需要对其进行加固处理,此时还需要知道临界滑动面的形状与位置,因此,临界滑动面的确定对岩土工程有着重要意义。目前对于临界滑动面的研究多采用上限法,基于下限定理或下限模型的临界滑动面理论与数值方法目前尚未见到。

      本文在下限模型的基础上,提出了基于下限模型的临界滑动面理论,给出了基于下限模型的临界滑动面数值解法;并以黏性土体直立边坡的临界高度问题为例,验证了基于下限模型的临界滑动面数值解法的适用性。

      (1)应力场

      在研究对象所在区域内,应力场可表示为

      σx=σx(x,y,z),σy=σy(x,y,z),σz=σz(x,y,z),τyz=τyz(x,y,z),τzx=τzx(x,y,z),τxy=τxy(x,y,z)} (1)

      将应力场(σx,σy,σz,τyz,τzx,τxy)简记为σ

      (2)平衡方程

      静力平衡微分方程为

      σxx+τyxy+τzxz=Fx,τxyx+σyy+τzyz=Fy,τxzx+τyzy+σzz=Fz} (2)

      将静力平衡微分方程简记为

      E(σ)=F  (3)

      (3)应力边界条件

      在应力边界Sσ上,设边界法方向为(l,m,n),在每一点给定了分布的表面力ˉTx, ˉTy 与 ˉTz,应力边界条件可写为

      σxl+τyxm+τzxn=ˉTx,τxyl+σym+τzyn=ˉTy,τxzl+τyzm+σzn=ˉTz} (4)

      将应力边界条件简记为

      S(σ)=T  (5)

      (4)屈服准则

      屈服准则以应力形式表达,一般地可写为

      M(σ)0 (6)

      下限定理可表述为:任何一个静力许可应力场都构成极限荷载的一个下限解。

      静力许可应力场是指同时满足如下条件的应力场:①平衡方程;②应力边界条件;③屈服准则约束。

      与下限定理相应的数学模型简称为下限模型。下限模型可写为

      max (7)

      式中: \sigma 为应力场,自变量; Q(\sigma ) 为目标函数; E(\sigma ) = F 为平衡方程; S(\sigma ) = T 为应力边界条件; M(\sigma ) \leqslant 0 为屈服准则。

      下限模型是以应力场为变量,以下限荷载为目标函数,以平衡方程、应力边界条件和屈服准则为约束的最优化模型。

      下限模型是一个最优化模型,具有最优化模型的一般特性。

      满足最优化模型全部约束条件的解称为可行解。目标函数的最大值称为最优值;使得目标函数取得最优值的可行解,称为最优解。对某一确定的最优化模型,可行解可能存在,也可能不存在。如果可行解不存在,则最优值也就不存在;如果可行解存在,最优值可能存在(为有限值),也可能不存在(为无限大)。如果最优值存在,则最优解必然存在。最优值只有一个,最优解可能为有限个,也可能为无限多个。

      下限模型的每一个可行解,都对应着一个静力许可应力场。下限模型目标函数的最优值即为极限荷载。下限模型的每一个最优解,都对应着极限荷载作用下的一个静力许可应力场。为论述方便,称极限荷载作用下的静力许可应力场为极限许可应力场。也即下限模型的每一个最优解,都对应着一个极限许可应力场。

      对某一确定的下限模型,静力许可应力场可能存在,也可能不存在。如果静力许可应力场不存在,则下限模型所描述的岩土问题是不稳定的,此时极限荷载也就不存在;如果静力许可应力场存在,则下限模型所描述的岩土问题是稳定的,此时极限荷载可能存在(为有限值),也可能不存在(为无限大)。如果极限荷载存在,则极限许可应力场必然存在。极限荷载只有一个,极限许可应力场可能为有限个,也可能为无限多个。

      在某一给定的极限许可应力场中,对某一点的应力而言,该点应力要么处于屈服状态,要么处于非屈服状态,二者必居其一。如果某一点的应力在任一极限许可应力场中均处于屈服状态,则称该点为下限塑性点。全部的下限塑性点的集合称为下限塑性区。下限塑性区即为下限模型的临界滑动面。

      下限塑性点的直观物理意义就是在极限荷载作用下不得不屈服的点;下限塑性区的直观物理意义就是在极限荷载作用下不得不屈服的区域。对二维平面问题,临界滑动面可能是一条二维曲线,也可能是一个二维区域;对三维问题,临界滑动面可能是一个三维曲面,也可能是一个三维区域。由于屈服准则的外凸性,下限塑性区内的应力分布具有唯一性,非下限塑性区内的应力分布不具有唯一性。

      由于实际问题的复杂性,一般很难获得下限模型的解析解,工程中可考虑采用数值解法。以平面应变问题例,采用四边形单元进行网格划分,以全部的节点应力为自变量,并对目标函数、平衡方程、应力边界条件和屈服准则进行线性化处理,可得到一个线性化的下限模型(称之为下限线性规划模型),通过对下限线性规划模型的求解,可获得相应问题的极限荷载及临界滑动面,具体如下。

      对平面应变问题,采用四边形单元进行网格划分。设网格共有N个节点,M个单元。

      以全部节点的应力分量为自变量,记为{\sigma _n}。因网格共有N个节点,故{\sigma _n}为一个3N维列向量。

      采用线性化的目标函数,可一般地记为

      C{\sigma _n}。 (8)

      式中:C为目标向量。

      单元在体力与周边应力作用下,在x方向和y方向上应分别保持平衡。设单元周边应力为直线分布,单元边界上x方向的应力分布如图 1所示,单元边界上y方向的应力分布如图 2所示,则单元在x方向与y方向上的平衡方程可写为

      [{A^{\text{e}}}][{\sigma ^{\text{e}}}] = [{b^{\text{e}}}] 。 (9)
      图  1  单元边界上x方向的应力分布
      Figure  1.  Stresses along element boundary in x direction
      图  2  单元边界上y方向的应力分布
      Figure  2.  Stress along the element boundary in y direction

      式中:

      [{A}^{\text{e}}]=\left[\begin{array}{c}\begin{array}{cc}\begin{array}{cc}\begin{array}{ccc}{\eta }_{1}& 0& {\zeta }_{1}\end{array}& \begin{array}{ccc}{\eta }_{2}& 0& {\zeta }_{2}\end{array}\end{array}& \begin{array}{cc}\begin{array}{ccc}{\eta }_{3}& 0& {\zeta }_{3}\end{array}& \begin{array}{ccc}{\eta }_{4}& 0& {\zeta }_{4}\end{array}\end{array}\end{array}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{ccc}0& {\zeta }_{1}& {\eta }_{1}\end{array}& \begin{array}{ccc}0& {\zeta }_{2}& {\eta }_{2}\end{array}\end{array}& \begin{array}{cc}\begin{array}{ccc}0& {\zeta }_{3}& {\eta }_{3}\end{array}& \begin{array}{ccc}0& {\zeta }_{4}& {\eta }_{4}\end{array}\end{array}\end{array}\end{array}\right]\text{;} \\ \left[{\sigma }^{\text{e}}\right]={\left[\begin{array}{cc}\begin{array}{cc}\begin{array}{ccc}{\sigma }_{x1}& {\sigma }_{y1}& {\tau }_{xy1}\end{array}& \begin{array}{ccc}{\sigma }_{x2}& {\sigma }_{y2}& {\tau }_{xy2}\end{array}\end{array}& \begin{array}{cc}\begin{array}{ccc}{\sigma }_{x3}& {\sigma }_{y3}& {\tau }_{xy3}\end{array}& \begin{array}{ccc}{\sigma }_{x4}& {\sigma }_{y4}& {\tau }_{xy4}\end{array}\end{array}\end{array}\right]}^{\text{T}}\text{;} \\ [{b^{\text{e}}}] = {\left[ {\begin{array}{*{20}{c}} {A{g_x}}&{A{g_y}} \end{array}} \right]^{\text{T}}}\text{;}
      {\eta }_{1}=\frac{1}{2}({y}_{2}-{y}_{4})\text{;}{\eta }_{2}=\frac{1}{2}({y}_{3}-{y}_{1}) \text{;}\\ {\eta }_{3}=\frac{1}{2}({y}_{4}-{y}_{2})\text{;}{\eta }_{4}=\frac{1}{2}({y}_{1}-{y}_{3}) \text{;}
      {\zeta _1} = \frac{1}{2}({x_4} - {x_2});{\zeta _2} = \frac{1}{2}({x_1} - {x_3})\text{;}\\{\zeta _3} = \frac{1}{2}({x_2} - {x_4});{\zeta _4} = \frac{1}{2}({x_3} - {x_1})。

      式中:A为单元面积;{g_x}{g_y}分别为x方向与y方向的重力加速度;({x_1},{y_1})({x_4},{y_4})为单元4个节点的坐标;({\sigma _{x1}},{\sigma _{y1}},{\tau _{xy1}}) ({\sigma _{x4}},{\sigma _{y4}},{\tau _{xy4}}) 元为单元4个节点的应力分量。

      对全部单元平衡方程集成,可得到如下形式的总平衡方程:

      {A_1}{\sigma _n} = {b_1}。 (10)

      应力边界条件实质上为施加在应力上的等式约束。采用线性化的应力边界条件,可一般地记为

      {A_2}{\sigma _n} = {b_2}。 (11)

      对Mohr-Coulomb屈服准则,采用线性化近似,式(6)可写为[4]

      {A}_{k}{\sigma }_{x}+{B}_{k}{\sigma }_{y}+{C}_{k}{\tau }_{xy}\le D\text{ }。 (12)

      式中: {A}_{k}=\mathrm{cos}(2\text{π}k\text{/}p)+\mathrm{sin}φ \mathrm{cos}(\text{π/}p);{B}_{k}=\mathrm{sin}φ \mathrm{cos}(\text{π/}p)- \mathrm{cos}(2\text{π}k\text{/}p);{C}_{k}=2\mathrm{sin}(2\text{π}k\text{/}p);D=2c\mathrm{cos}φ \mathrm{cos}(\text{π/}p);k= 1,2,…,pp为内接多边形的边数;c为黏聚力;φ为内摩擦角。

      所有的节点应力均应满足屈服准则。将所有的线性化后的节点应力屈服准则约束集成,即可得到线性化的屈服准则约束,可一般地记为

      {A_3}{\sigma _n} \leqslant {b_3}。 (13)

      综合目标函数式(8)、平衡方程式(10)、应力边界条件式(11)及屈服准则式(13),下限线性规划模型可一般地写为

      \left.\begin{array}{ccc} \max & C^{\mathrm{T}} \sigma_n, \\ \text { st. } & A_1 \sigma_n=b_1 & , \\ & A_2 \sigma_n=b_2 & , \\ & A_3 \sigma_n \leqslant b_3 & 0 \end{array}\right\} (14)

      式中:{\sigma _n}为节点应力;C为目标向量;{A_1}{\sigma _n} = {b_1}为离散后的线性化平衡方程;{A_2}{\sigma _n} = {b_2}为离散后的线性化应力边界条件;{A_3}{\sigma _n} \leqslant {b_3}为离散后的线性化节点屈服准则约束。

      模型式(14)为带有等式与不等式约束的线性规划模型,其目标函数的最大值即为极限荷载。对下限线性规划模型的任一最优解,如果某一不等式约束均处于等式状态,称该约束为主动约束,称该约束所在的节点为下限塑性节点。全部下限塑性节点的集合构成下限模型的临界滑动面。据此,求解下限模型的临界滑动面问题就转化为求解线性规划模型的主动约束问题。关于线性规划模型最优值及主动约束的求解,可参见相关线性规划相关内容[18-19],本文不再做详细论述。

      以黏性土体直立边坡的临界高度问题为例,对本文所提出的临界滑动面理论与数值解法进行简单验证。

      设土体重度\gamma = 18 kN/m3,黏聚力c=20 kPa,内摩擦角\varphi = 15°。设土体竖直方向高度为h,水平方向宽度b=h;竖直方向与水平方向均划分为15个单元,如图 3所示,共计有256个节点,225个单元。其中左侧边界AD上节点的切应力与x向正应力为零,上侧边界AB上节点的切应力与y向正应力为零;采用Mohr-Coulomb屈服准则。

      图  3  黏性土体直立边坡的临界高度问题
      Figure  3.  Critical heights of vertical slope in cohesive soil

      目标向量给定为零,也即只要下限模型存在可行解,边坡就是稳定的,否则是不稳定的。根据本文中所述方法,建立黏性土体直立边坡临界高度问题的线性规划下限模型。

      数值求解时,先通过试算法求解直立边坡的临界高度。即先假定一个竖直高度,如果求解发现下限模型的可行解不存在,则减小假定的竖直高度;如果求解发现下限模型的可行解存在,则增加假定的竖直高度,直至高度范围小于给定的精度为止,这样就获得了黏性直立边坡的临界高度。在获得临界高度后,进一步求解下限模型的全部下限塑性节点。

      根据卡尔曼公式[2],其临界高度hcr1为5.79 m。采用本文所述方法,黏性土体直立边坡的临界高度hcr2为5.36 m,全部下限塑性节点(Plastic node)分布如图 3所示。卡尔曼公式假定临界滑动面为直线(Linear sliding)。根据本文提出的临界滑动面理论,全部塑性节点构成临界滑动面。由图 3中的塑性节点分布可以看出,其临界滑动面总体上呈曲线分布。

      在理论研究与工程实践中,下限定理通常用来求解极限荷载,临界滑动面通常通过上限法来求解。本文的研究表明,基于下限模型不但可以求解极限荷载,也可以求解临界滑动面。

      在建立起线性规划下限模型后,极限荷载的求解可借助成熟的商业软件进行。与临界滑动面求解相对应的是线性规划模型主动约束的求解,目前的商业软件并不直接提供主动约束求解这一功能,需要编程完成。主动约束求解的标准化与模块化应是下一步需要解决的问题。

      本文只是初步提出了基于下限模型求解临界滑动面理论和数值解法,更多的理论研究需要进一步的深入,更多的实践验证还有待开展。

    • 图  1   填埋场耦合模型求解示意图

      Figure  1.   Solution procedures of proposed coupled model for landfills

      图  2   填埋场好氧修复的几何示意图

      Figure  2.   Schematic of aerobic remediation of landfill

      图  3   曝气条件下填埋场的沉降变形空间分布

      Figure  3.   Spatial distribution of deformation in landfill under aeration

      图  4   好氧修复过程中垃圾体沉降变形网格示意图

      Figure  4.   Grid diagram of waste deformation under aeration

      图  5   不均匀沉降量随时间的变化规律

      Figure  5.   Variation of uneven settlement over time

      图  6   曝气井-垃圾相对位移随时间的变化

      Figure  6.   Variation of relative displacement of well-waste over time

      图  7   曝气压力对井周相对沉降量的影响

      Figure  7.   Effects of aeration pressure on relative settlement of well

      图  8   氧浓度和有机物降解度Dd的空间分布

      Figure  8.   Distribution of O2 concentration and degradation degree Dd

      图  9   两种降解模式下填埋场沉降速率

      Figure  9.   Settlement rates of landfill under two degradation modes

      图  10   好氧修复后地表沉降随时间的变化

      Figure  10.   Change of surface subsidence over time after aeration

      表  1   垃圾体好氧反应多场耦合模型的核心变量

      Table  1   Key variables of coupled model for anaerobic-aerobic landfill

      核心变量符号单位
      液相和气相压力pw, pgPa
      气相CH4,CO2质量分数ygCH4, ygCO2
      气相O2,N2质量分数ygO2, ygN2
      骨架位移增量dum
      下载: 导出CSV

      表  2   填埋场好氧修复的边界条件

      Table  2   Boundary conditions of aerobic remediaton of landfill

      边界PgPwyCH4 g,yCO2 gyO2 gyN2 gdu
      顶部0 Pa不透水零梯度零梯度零梯度应力0 Pa
      底部不透气自由排水零梯度零梯度零梯度固定
      曝气井6 kPa不透水0%26.6%73.4%固定
      回收井0 Pa不透水零梯度零梯度零梯度固定
      四周  对称边界  滑移
      下载: 导出CSV

      表  3   垃圾生化反应和甲烷氧化参数

      Table  3   Waste biodegradation and CH4 oxidation parameters

      参数单位
      好氧和厌氧降解常数km,A, km,Nd-10.1, 0.01
      好氧和厌氧生长率半饱和常数ks,A, ks,Nkg/m3100
      好氧反应氧气半饱和常数KO20.07
      初始好氧菌和厌氧菌浓度XA,0, XN,0kg/m30.15
      最大CH4消耗速率Vmaxkg/m3/s0.6×10-5
      甲烷氧化O2半饱和常数DO20.012
      甲烷氧化CH4半饱和常数DCH40.0066
      下载: 导出CSV
    • [1] 孔宪京, 孙秀丽. 城市固体废弃物沉降模型研究现状及其进展[J]. 大连理工大学学报, 2006, 46(4): 615-624. doi: 10.3321/j.issn:1000-8608.2006.04.029

      KONG Xian-jing, SUN Xiu-li. Current situation and advance of research on settlement model for municipal solid wastes[J]. Journal of Dalian University of Technology, 2006, 46(4): 615-624. (in Chinese) doi: 10.3321/j.issn:1000-8608.2006.04.029

      [2] 李秀兰, 刘东燕, 罗云菊. 山地城市垃圾填埋场渗滤液环境问题的探讨[J]. 地下空间与工程学报, 2006, 2(5): 855-858. doi: 10.3969/j.issn.1673-0836.2006.05.035

      LI Xiu-lan, LIU Dong-yan, LUO Yun-ju. Discussion on the environmental problems from leachate of solid waste in mountainous city[J]. Chinese Journal of Underground Space and Engineering, 2006, 2(5): 855-858. (in Chinese) doi: 10.3969/j.issn.1673-0836.2006.05.035

      [3]

      NIKOLAOU A, GIANNIS A, GIDARAKOS E. Comparative studies of aerobic and anaerobic treatment of MSW organic fraction in landfill bioreactors[J]. Environmental Technology, 2010, 31(12): 1381-1389. doi: 10.1080/09593331003743104

      [4]

      KO J H, MA Z Y, JIN X, et al. Effects of aeration frequency on leachate quality and waste in simulated hybrid bioreactor landfills[J]. Journal of the Air & Waste Management Association, 2016, 66(12): 1245-1256.

      [5]

      RAGA R, COSSU R. Landfill aeration in the framework of a reclamation project in Northern Italy[J]. Waste Management, 2014, 34(3): 683-691. doi: 10.1016/j.wasman.2013.12.011

      [6]

      HEYER K U, HUPE K, RITZKOWSKI M, et al. Pollutant release and pollutant reduction - Impact of the aeration of landfills[J]. Waste Management, 2005, 25(4): 353-359. doi: 10.1016/j.wasman.2005.02.007

      [7]

      HOSSAIN M S, GABR M A, BARLAZ M A. Relationship of compressibility parameters to municipal solid waste decomposition[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(12): 1151-1158. doi: 10.1061/(ASCE)1090-0241(2003)129:12(1151)

      [8]

      MARQUES A C M, FILZ G M, VILAR O M. Composite compressibility model for municipal solid waste[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(4): 372-378. doi: 10.1061/(ASCE)1090-0241(2003)129:4(372)

      [9]

      BEAVEN R, POWRIE W. Determination of the hydrogeological and geotechnical properties of refuse in relation to sustainable landfilling[C]//Proceedings of the Annual Madison Waste Conference, 1996, Madison.

      [10]

      BAREITHER C A, BENSON C H, EDIL T B. Compression behavior of municipal solid waste: immediate compression[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2012, 138(9): 1047-1062. doi: 10.1061/(ASCE)GT.1943-5606.0000672

      [11]

      BAREITHER C A, BENSON C H, EDIL T B. Compression of municipal solid waste in bioreactor landfills: mechanical creep and biocompression[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(7): 1007-1021. doi: 10.1061/(ASCE)GT.1943-5606.0000835

      [12]

      SIVAKUMAR BABU G L, REDDY K R, CHOUKSEY S K. Constitutive model for municipal solid waste incorporating mechanical creep and biodegradation-induced compression[J]. Waste Management, 2010, 30(1): 11-22. doi: 10.1016/j.wasman.2009.09.005

      [13]

      OLIVIER F, GOURC J P. Hydro-mechanical behavior of Municipal Solid Waste subject to leachate recirculation in a large-scale compression reactor cell[J]. Waste Management, 2007, 27(1): 44-58. doi: 10.1016/j.wasman.2006.01.025

      [14]

      MACHADO S L, CARVALHO M F, VILAR O M. Constitutive model for municipal solid waste[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2002, 128(11): 940-951. doi: 10.1061/(ASCE)1090-0241(2002)128:11(940)

      [15]

      MCDOUGALL J. A hydro-bio-mechanical model for settlement and other behaviour in landfilled waste[J]. Computers and Geotechnics, 2007, 34(4): 229-246. doi: 10.1016/j.compgeo.2007.02.004

      [16]

      SIVAKUMAR BABU G L, REDDY K R, CHOUSKEY S K, et al. Prediction of long-term municipal solid waste landfill settlement using constitutive model[J]. Practice Periodical of Hazardous, Toxic, and Radioactive Waste Management, 2010, 14(2): 139-150. doi: 10.1061/(ASCE)HZ.1944-8376.0000024

      [17] 柯瀚, 董鼎, 陈云敏, 等. 考虑剪缩性的城市固体废弃物非线性弹性模型[J]. 浙江大学学报(工学版), 2017, 51(11): 2158-2164. doi: 10.3785/j.issn.1008-973X.2017.11.009

      KE Han, DONG Ding, CHEN Yun-min, et al. Nonlinear elastic model for municipal solid waste considering dilatancy effect[J]. Journal of Zhejiang University (Engineering Science), 2017, 51(11): 2158-2164. (in Chinese) doi: 10.3785/j.issn.1008-973X.2017.11.009

      [18]

      KIM S Y, TOJO Y, MATSUTO T. Compartment model of aerobic and anaerobic biodegradation in a municipal solid waste landfill[J]. Waste Management & Research, 2007, 25(6): 524-537.

      [19]

      THEMELIS N J, KIM Y H. Material and energy balances in a large-scale aerobic bioconversion cell[J]. Waste Management & Research, 2002, 20(3): 234-242.

      [20]

      ABICHOU T, MAHIEU K, CHANTON J, et al. Scaling methane oxidation: From laboratory incubation experiments to landfill cover field conditions[J]. Waste Management, 2011, 31(5): 978-986. doi: 10.1016/j.wasman.2010.12.002

      [21]

      SPOKAS K, BOGNER J, CHANTON J. A process-based inventory model for landfill CH4 emissions inclusive of seasonal soil microclimate and CH4 oxidation[J]. Journal of Geophysical Research: Biogeosciences, 2011, 116(G4): G04017.

      [22] 郭亚丽, 赵由才, 徐迪民. 上海老港生活垃圾填埋场陈垃圾的基本特性研究[J]. 上海环境科学, 2002, 21(11): 669-671.

      GUO Ya-li, ZHAO You-cai, XU Di-ming. Study on basic characteristics of aged refuse at Shanghai Laogang landfill[J]. Shanghai Environmental Sciences, 2002, 21(11): 669-671. (in Chinese)

    图(10)  /  表(3)
    计量
    • 文章访问数:  281
    • HTML全文浏览量:  25
    • PDF下载量:  154
    • 被引次数: 0
    出版历程
    • 收稿日期:  2021-03-28
    • 网络出版日期:  2022-12-01
    • 刊出日期:  2021-10-31

    目录

    /

    返回文章
    返回