土木建筑与环境工程  2013, Vol. 35 Issue (5): 44-49   PDF    
城市生活垃圾填埋场气压分布一维稳态分析模型
马小飞1,2, 李育超1, 詹良通1, 陈云敏1    
1. 浙江大学 软弱土与环境土工教育部重点实验室, 杭州 310058;
2. 杭州市建筑设计研究院有限公司, 杭州 310001
收稿日期：2013-02-20
基金项目：国家自然科学基金（50878194、51010008、51009121）；教育部留学回国人员启动基金（教外司留2012-1707）；苏州市科技支撑计划项目（520324）
作者简介：李育超(1978-), 男, 博士, 主要从事环境岩土工程研究, (E-mail)yuchao_li@hotmail.com。
摘要：掌握填埋场气压分布是填埋气灾害控制和资源化利用的基础。提出了分层垃圾填埋体气压分布一维稳态分析模型及求解方法。该模型可分析含有给定抽气压力或给定抽气流量水平导气层的填埋体气压分布。采用该模型探讨了垃圾分层特征、封顶覆盖层下和填埋体内的高渗透性水平导气层、填埋体底部渗沥液导排系统兼作填埋气导排通道对填埋体气压分布的影响规律。通过参数分析得到以下结论：填埋体垃圾不考虑产气速率和固有渗透系数随垃圾埋深增加而减小（产气速率和固有渗透系数按埋深平均）会高估填埋体内气压；封顶覆盖层下和填埋体内设置水平导气层可有效降低填埋体内气压；填埋体底部填埋气导排通道对填埋气气压较大的深部垃圾降压效果较明显。分析结果表明，该模型可指导填埋气收集系统设计，如水平导气层位置和间距布置及抽气功率选择等。
关键词：填埋场    分层垃圾    气压    填埋气收集系统    导气层    
One-dimensional Steady-state Model for Gas Pressure Distribution in Municipal Solid Waste Landfills
Ma Xiaofei1,2, Li Yuchao1, Zhan Liangtong1, Chen Yunmin1    
1. MOE Key Laboratory of Soft Soils and Geoenvironmental Engineering, Zhejiang University, Hangzhou 310058, P. R. China;
2. Hangzhou Archtectural Design & Research Co, LTD, Hangzhou 310001, P. R. China
Received: 2013-02-20
Abstract: Gas pressure distribution in landfills is of great importance to control hazards and improve recovery of landfill gas. An one-dimensional steady-state model and its solving approach are presented for gas pressure distribution in layered waste landfills. The model can be used to investigate the gas pressure distribution of the landfill with horizontal gas collection layers having specified pressure or flux. The influences of layered properties, horizontal permeable layer under final cover and horizontal gas collection layers and bottom leachate drainage systems dually serving as gas collection layers on gas pressure distribution are investigated. The results of parameter analyses using the presented model show that the assumption that homogenizations of gas generation rate and gas permeability through a landfill considerably overestimate gas pressure in landfill; horizontal permeable layer under final cover and horizontal gas collection layers can effectively reduce gas pressure in landfills; the bottom leachate drainage systems dually serving as gas collection layers can significantly reduce the gas pressure in the lower portion of landfills. The presented model provides a useful tool in the designing of landfill gas collection systems, such as, the positions and spacing of horizontal gas collection layers and pump section.
Key Words: landfill    layered waste    gas pressure    landfill gas collection system    permeable layer    

城市生活垃圾中含有大量的有机质，填埋后发生生化降解，产生大量主要成分为温室气体甲烷和二氧化碳的填埋气(其中甲烷体积约占一半)。填埋气无序排放，会加剧全球气候变暖，同时造成巨大的能源浪费(甲烷具有资源化利用价值)。垃圾填埋体中填埋气导排运移受阻造成的气体累积和气压上升，易引发填埋体失稳和物理爆炸；另一方面，填埋气与空气混合极易引发化学爆炸[1-3]。现代卫生填埋场应当科学设计填埋气导排和收集系统，控制填埋体气压，防止填埋气相关灾害，有效收集填埋气以资源化利用。

掌握填埋场内气压分布规律是合理设计填埋气导排和收集系统设计的基础。Findikakis等[4]和El-Fadel等[5-6]提出了单相一维对流扩散运移模型，仅考虑填埋气从填埋体顶部排出情况，分析了填埋体中气压和气体组分随时间的变化规律。Townsend等[7]建立了均质垃圾填埋体气压分布模型，分析了垃圾产气速率、填埋厚度和固有渗透系数对气压分布的影响，并针对填埋场顶部导气层和底部渗沥液导排系统兼作填埋气收集的情况分析了填埋气收集对填埋体气压分布的影响。上述模型均假设填埋体为均质，然而实际填埋场垃圾因填埋作业顺序的原因，工程特性在深度方向具有明显的分层特征。填埋垃圾的产气速率和固有渗透系数随填埋深度增加而减小[8-9]，上述已有模型无法考虑填埋垃圾产气速率和固有渗透系数的分层特征对气压分布的影响。魏海云等[8]提出了分层填埋体气压的有限差分分析模型，该模型在填埋体沿深度方向空间离散求解。国内学者对于竖井抽气条件下的填埋场气压分布研究较多[10-11]，水平导气层对气压分布的影响规律研究还较少。

提出可考虑填埋体分层特征和水平导气层的气压一维稳态分析模型，采用线性表达式描述填埋体边界条件和层间连续条件，提出给定气压和给定填埋气流量两种分析水平导气层的方法，并给出了相应的求解方法，通过4个例子分析了垃圾产气速率和固有渗透系数分层特征、封顶覆盖层下和填埋体内的高渗透性水平导气层、填埋体底部渗沥液导排系统兼作填埋气导排通道对填埋体气压分布的影响规律，指导填埋气收集系统设计。

1 分析模型

考虑实际填埋场垃圾生化降解时间随埋深增大的规律，针对填埋垃圾的分层特征提出以下气压分析模型。将填埋体分为n层，各层厚度为hi，产气速率和固有渗透系数分别为ai和Ki，如图 1所示，填埋体总厚度为H，z为垃圾埋深，第i层和第i+1层层间界面埋深为zi。各层可为具有不同填埋龄期的垃圾层、水平导气层、封顶覆盖层或中间覆盖层。

图 1 分层垃圾填埋体模型图

填埋体中某一层材料(垃圾层、导气层或覆盖层)的填埋气流动可根据多孔介质流体力学理论[12]采用一维气体流动控制方程表示：

$ - \frac{{\partial \left[ {{K_i}{\rho _i}\left( {\frac{{\partial {p_i}}}{{\partial {z_i}}} - {\rho _i}g} \right)} \right]}}{{\mu \partial {z_i}}} + \frac{{\partial \left( {{n_{\rm{g}}}{\rho _i}} \right)}}{{\partial t}} = {a_i} $ (1)

式中：Ki为固有渗透系数，m2；ρi为填埋气密度，kg/m3；pi为气体压力，Pa；μ为填埋气粘滞系数，1.37×10-5 kg/(m·s)；ng为垃圾含气率，无量纲；ai为产气速率，即单位体积垃圾在单位时间内的产填埋气质量，kg/(m3·s)，其中导气层和覆盖层ai为零。

根据理想气体状态方程，填埋气密度与气压存在如下关系：

$ {\rho _i} = \frac{\omega }{{RT}}{p_i} $ (2)

式中：ω为填埋气摩尔质量，0.03 kg/mol；R为气体常数，8.314 3 kg·m2/(s2·mol·K)；T为气体开氏温度, K，取298 K。

根据学者们研究与分析[7, 13-16]，填埋垃圾生化降解过程一般长达数十年，其生化降解速率较慢，填埋垃圾的产气速率ai和固有渗透系数Ki随填埋时间变化缓慢，可用稳态方程分析一定时期内填埋体气压分布，忽略式(1)中的重力项(填埋气密度ρi较小)，则式(1)可简化为：

$ \frac{{{\partial ^2}p_i^2}}{{\partial z_i^2}} = - \frac{{2\mu RT}}{{\omega {K_i}}}{a_i} $ (3)

令${A_i} = \frac{{\mu RT}}{{\omega {K_i}}}{a_i}$，并对式(3)两次积分可得：

$ p_i^2 = - {A_i}z_i^2 + {B_i}{z_i} + {C_i} $ (4)

式中：Bi和Ci为待定积分常数，由填埋体各层边界条件确定。

层与层之间满足以下连续条件(以第i层和第i+1层为例)：1) 气压连续，即pi=pi+1；2) 气体流量连续，即qi=qi+1，其中流量qi可表示为

$ {q_i} = {v_i}{\rho _i} = - \frac{{{K_i}}}{\mu }\frac{{\partial {p_i}}}{{\partial {z_i}}}\frac{{{p_i}\omega }}{{RT}} = - \frac{{{K_i}\omega }}{{2\mu RT}}\frac{{\partial p_i^2}}{{\partial {z_i}}} $ (5)

式中vi为气体流速。将式(4)和(5)分别代入层间边界条件1)和2)可得以下方程：

$ - {A_i}z_i^2 + {B_i}{z_i} + {C_i} = - {A_{i + 1}}z_i^2 + {B_{i + 1}}{z_i} + {C_{i + 1}} $ (6)
$ {K_i}\left( { - 2{A_i}{z_i} + {B_i}} \right) = {K_{i + 1}}\left( { - 2{A_{i + 1}}{z_i} + {B_{i + 1}}} \right) $ (7)

填埋体的顶部边界条件一般可设为给定气压(即大气压)。底部边界条件通常填埋气流量为零，即底部无填埋气外漏情况，渗沥液水位雍高的中国填埋场多可采用该底部边界条件；另外，根据Townsend等[7, 17]和Fredenslund等[18]，其他国家低水位填埋场，如渗滤液导排层水头控制在30 cm以内[19]，底部渗沥液导排系统兼作填埋气导排通道，如在导排层上部铺设导气管单独导排填埋场深部填埋气，底部边界条件可设为给定气压(如直接与大气相连作为填埋气导排通道)或给定填埋气流量(与抽气措施相连作为填埋气导排通道)。上述边界条件可表示为：

$ {C_1} = p_0^2\left( {{p_1} = {p_0}\;在\;z = {z_1};其中\;{p_0}\;为大气压} \right) $ (8)
$ \begin{array}{l} - \frac{{{K_n}\omega }}{{2RT}}\left( { - 2{A_n}{z_n} + {B_n}} \right) = {q_b}\left( {{q_n} = {q_{\rm{b}}}\;在\;z = H;} \right.\\ \left. {其中\;{q_{\rm{b}}}\;为给定流量} \right) \end{array} $ (9a)
$ \begin{array}{l} - {A_n}z_n^2 + {B_n}{z_n} + {C_n} = p_{\rm{b}}^2\left( {{p_n} = {p_{\rm{b}}}\;在\;z = H;} \right.\\ \left. {其中\;{p_{\rm{b}}}\;为给定气压} \right) \end{array} $ (9b)

上述2n-2个层间边界条件方程以及顶部和底部边界条件方程，共有2n个线性方程，n层共有2n个待定参数(即Bi和Ci各n个)，因此该问题解析可求解。当底部边界为给定气压时，上述边界方程可用如下线性方程组表示：

$ \left[ {\begin{array}{*{20}{c}} {{K_1}}&0&{ - {K_2}}&0&0&0& \cdots &0&0&0&0& \cdots &0&0&0\\ {{z_1}}&1&{ - {z_1}}&1&0&0&{}&0&0&0&0&{}&0&0&0\\ 0&0&{{K_2}}&0&{ - {K_3}}&0& \cdots &0&0&0&0& \cdots &0&0&0\\ 0&0&{{z_2}}&1&{ - {z_2}}&{ - 1}&{}&0&0&0&0&{}&0&0&0\\ \cdots &{}& \cdots &{}& \cdots &{}& \ddots&\cdots &{}& \cdots &{}& \ddots &{}& \cdots &{}\\ 0&0&0&0&0&0& \cdots &{{K_i}}&0&{ - {K_{i + 1}}}&0& \cdots &0&0&0\\ 0&0&0&0&0&0& \cdots &{{z_i}}&1&{ - {z_i}}&{ - 1}& \cdots &0&0&0\\ \cdots &{}& \cdots &{}& \cdots &{}& \ddots&\cdots &{}& \cdots &{}& \ddots &{}& \cdots &{}\\ 0&0&0&0&0&0&{}&0&0&0&0&{}&{{K_{n - 1}}}&0&{ - {K_n}}\\ 0&0&0&0&0&0& \cdots &0&0&0&0& \cdots &{{z_{n - 1}}}&1&{ - {z_{n - 1}}}\\ 0&1&0&0&0&0&{}&0&0&0&0&{}&0&0&0\\ 0&0&0&0&0&0& \cdots &0&0&0&0& \cdots &0&0&{{z_n}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{B_1}}\\ {{C_1}}\\ {{B_2}}\\ {{C_2}}\\ \cdots \\ {{B_i}}\\ {{C_i}}\\ \cdots \\ {{B_{n - 1}}}\\ {{C_{n - 1}}}\\ {{B_n}}\\ {{C_n}} \end{array}} \right] \\ = \left[ {\begin{array}{*{20}{c}} {2{z_1}\left( {{K_1}{A_1} - {K_2}{A_2}} \right)}\\ {\left( {{A_1} - {A_2}} \right) \cdot z_1^2}\\ {2{z_2}\left( {{K_2}{A_2} - {K_3}{A_3}} \right)}\\ {\left( {{A_2} - {A_3}} \right) \cdot z_2^2}\\ \cdots \\ {2{z_i}\left( {{K_i}{A_i} - {K_{i + 1}}{A_{i + 1}}} \right)}\\ {\left( {{A_i} - {A_{i + 1}}} \right) \cdot z_i^2}\\ \cdots \\ {2{z_{n - 1}}\left( {{K_{n - 1}}{A_{n - 1}} - {K_n}{A_n}} \right)}\\ {\left( {{A_{n - 1}} - {A_n}} \right) \cdot z_{n - 1}^2}\\ {p_0^2}\\ {{A_n}z_n^2 + p_b^2} \end{array}} \right] $ (10)

式(10)可通过线性变换求解所有未知量Bi、Ci，当填埋层分层数n较小时可手算求解。当底部边界条件为给定流量时，可按下面方法直接求得Bi、Ci。由式(9a)直接求得Bn：

$ {B_n} = 2{A_n}{z_n} - \frac{{2RT}}{{{K_n}\omega }}{q_{\rm{b}}} $ (11)

再通过式(7)可依次求得Bn-1、Bn-2、…、B1：

$ {B_i} = \frac{{{K_{i + 1}}\left( { - 2{A_{i + 1}}{z_i} + {B_{i + 1}}} \right)}}{{{K_i}}} + 2{A_i}{z_i} $ (12)

求得全部Bi后，由顶部边界条件(即式(8))可得C1，进而通过式(6)可依次求得C2、C3、…、Cn：

$ {C_{i + 1}} = \left( {{A_{i + 1}} - {A_i}} \right)z_i^2 + \left( {{B_i} - {B_{i + 1}}} \right){z_i} + {C_i} $ (13)

求解所有未知量Bi、Ci后，通过式(4)得到各层的气压分布。

若填埋体中含有水平导气层(假设为第i层)，假定导气层中气压随深度无变化，即有Bi=0。对于水平导气层抽气速率为qs的情况，原方程组中第i层和第i-1层、第i层和第i+1层的4个层间边界条件方程可替换为水平导气层与上下层间的边界条件方程(水平导气层不产气，即有Ai=0)：

$ {B_i} = 0\left( {第\;i\;层气压不随深度变化} \right) $ (14)
$ \begin{array}{l} - {A_{i - 1}}z_{i - 1}^2 + {B_{i - 1}}{z_{i - 1}} + {C_{i - 1}} = {C_i}\left( {{p_{i - 1}} = {p_i}\;在} \right.\\ \left. {z = {z_{i - 1}}} \right) \end{array} $ (15)
$ \begin{array}{l} - \frac{{{K_{i - 1}}\omega }}{{2\mu RT}}\left( { - 2{A_{i - 1}}{z_{i - 1}} + {B_{i - 1}}} \right) + \frac{{{K_{i - 1}}\omega }}{{2\mu RT}}\left( { - 2{A_{i + 1}}{z_i}} \right.\\ \left. { + {B_{i + 1}}} \right) = {q_{\rm{s}}}\left( {{q_{i - 1}} - {q_{i + 1}} = {q_{\rm{s}}}} \right) \end{array} $ (16)
$ \begin{array}{l} - {A_{i + 1}}z_i^2 + {B_{i + 1}}{z_i} + {C_{i + 1}} = {C_i}\left( {{p_{i + 1}} = {p_i}\;在\;z} \right.\\ \left. { = {z_i}} \right) \end{array} $ (17)

若水平导气层为已知气压ps(气压定值模型)，同理可采用以下4个方程替换4个层间边界条件方程，即：

$ {B_i} = 0\left( {第\;i\;层气压不随深度变化} \right) $ (18)
$ \begin{array}{l} - {A_{i - 1}}z_{i - 1}^2 + {B_{i - 1}}{z_{i - 1}} + {C_{i - 1}} = p_{\rm{s}}^2\left( {{p_{i - 1}} = {p_{\rm{s}}}\;在} \right.\\ \left. {z = {z_{i - 1}}} \right) \end{array} $ (19)
$ {C_i} = p_{\rm{s}}^2\left( {第\;i\;层气压\;{p_i} = {p_{\rm{s}}}} \right) $ (20)
$ \begin{array}{l} - {A_{i + 1}}z_i^2 + {B_{i + 1}}{z_i} + {C_{i + 1}} = p_{\rm{s}}^2\left( {{p_{i + 1}} = {p_{\rm{s}}}\;在} \right.\\ \left. {z = {z_i}} \right) \end{array} $ (21)

求解得到填埋体气压分布后，可通过下式求得水平导气层抽气速率，指导水平导气层抽气泵选择：

$ \begin{array}{l} {q_{CL}} = - \frac{{{K_{i - 1}}\omega }}{{2\mu RT}}\left( { - 2{A_{i - 1}}{z_{i - 1}} + {B_{i - 1}}} \right) + \frac{{{K_{i + 1}}\omega }}{{2\mu RT}}\left( - \right.\\ \left. {2{A_{i + 1}}{z_i} + {B_{i + 1}}} \right) \end{array} $ (22)
2 模型运用

填埋场气压现场监测主要针对垃圾表层，鲜有随深度的分层气压监测资料。采用Townsend等[7]的单层均质填埋体算例及魏海云等[8]的分层填埋体算例验证了本文提出模型，并分析了垃圾填埋龄期、低渗透性的中间覆盖层因素对填埋气气压分布的影响。本文重点分析垃圾分层特征、封顶覆盖层下和填埋体内的高渗透性水平导气层、填埋体底部渗沥液导排系统兼作填埋气导排通道对填埋场气压分布的影响。填埋场堆体高度为30 m，填埋速率为3 m/a，分析时将填埋体分5层，每层厚度为6 m，各层垃圾龄期取其平均值。填埋垃圾的产气速率a随填埋龄期的变化规律根据《生活垃圾填埋场填埋气体收集处理及利用工程技术规范》(CJJ 133-2009)中的公式确定，即

$ a\left( t \right) = M{L_0}k{{\rm{e}}^{ - kt}} $ (23)

式中：t为填埋垃圾龄期，a；a(t)为填埋垃圾在填埋龄期t的产气速率，m3/a；M为填埋垃圾的重量, t；L0为单位重量垃圾的最大产气量, m3/t；k为填埋垃圾产气速率常数, 1/a。文中L0取150 m3/t，k取湿润、中等湿润气候范围内值0.12，填埋垃圾密度取800 kg/m3。

深度小于18 m的填埋垃圾固有渗透系数根据Jain等[9]的现场实验数据取值，深度超过18 m的填埋垃圾固有渗透系数根据魏海云[20]实验数据所得固有渗透系数随深度变小趋势变缓的规律取值。填埋场封场时各垃圾层产气速率a和固有渗透系数K具体取值如表 1所示：

表 1 各垃圾层产气速率和固有渗透系数取值

模型边界条件采用顶部边界条件为给定气压(大气压)，底部边界条件为给定零流量。下面分析除说明外均采用上述封场时垃圾分层产气速率、固有渗透系数取值和边界条件。

2.1 垃圾分层特征对气压分布的影响

首先，采用本文模型分析垃圾分层与否对封场时填埋体内气压计算结果的影响(图 2)。若产气速率和固有渗透系数均分层取值，气压梯度在层间处具有明显变化，这是由于层间固有渗透系数不同产生的。若固有渗透系数分层取值而产气速率取各层平均值，气压梯度仍在层间处呈现明显变化，但由于产气速率取值按埋深平均后，深部垃圾产气速率高于其实际值，加上深部垃圾的固有渗透系数较低，因此填埋体底部气压明显高于产气速率和固有渗透系数均分层取值的分析结果。若产气速率分层取值而固有渗透系数取各层平均值，气压梯度无上述层间突变现象，但由于固有渗透系数取值按埋深平均后，浅部垃圾固有渗透系数低于其实际值，而浅部垃圾因填埋龄期较短而产气速率较高，因此填埋体浅部气压明显高于产气速率和固有渗透系数均分层取值的分析结果。若产气速率和固有渗透系数均取各层平均值，填埋体气压整体明显高于产气速率和固有渗透系数均分层取值分析结果，例如深度12 m处高出88%。可见，产气速率和固有渗透系数按埋深平均取值计算会显著高估填埋体内气压。本文得到的无覆盖条件下填埋体气压值与文献报道的无覆盖条件下气压计算值相近[7-8]。

图 2 垃圾分层特征对气压分布的影响

2.2 土质封顶覆盖层下设置水平导气层对气压分布的影响

与垃圾相比，土质覆盖层固有渗透系数明显较小，假设填埋体表面设有一厚度为1 m[13, 21-23]、固有渗透系数为1.0×10-14m2[21]的土质封顶覆盖层。由图 3可见，设置土质封顶覆盖层后封场填埋体表面气压增大为10.7 kPa，该部位气压过大可影响土质封顶覆盖层稳定性。如在土质封顶覆盖层设置一厚度为60 cm、固有渗透系数为3.0×10-9m2的高渗透性水平导气层[16]。气压的减小值与水平导气层抽气速率q基本成正比例关系。当水平导气层抽气速率为9.0×10-6 kg/(m2·s)时，填埋体表面气压下降至0.9 kPa。可见，水平导气层可有效减小填埋体表面气压。填埋场工程设计中可根据土质封顶覆盖层下气压控制要求，根据本文模型分析结果指导合理选择抽气泵。

图 3 封顶覆盖层下设置水平导气层对气压分布的影响

2.3 填埋体内水平导气层对气压分布的影响

在填埋体内设置水平导气层是高效收集填埋气的一种常用手段。假定在填埋体埋深12和24 m处各设有一厚度为60 cm的水平导气层(导气层上下两层垃圾厚度各减小30 cm)，其固有渗透系数为3.0×10-9m2，并假设12和24 m处导气层抽气速率分别为q1和q2。由图 4可见，当q2=0 kg/(m2·s)时，q1从2.0×10-6 kg/(m2·s)增大到4.0×10-6 kg/(m2·s)，12 m埋深导气层以下垃圾的填埋气气压在相同埋深处均降低了0.24 kPa，即导气层以下气压分布曲线形态不随导气层抽气速率变化而变化。这是由于稳态条件下导气层下部气压分布形状由该部分的产气速率和固有渗透系数决定。当抽气速率为2.0×10-6 kg/(m2·s)的导气层从埋深12 m处下降至埋深24 m处时，底部最大气压降低幅度从13%增大到49%，但12 m以上垃圾气压分布没有变化。这是因为在稳态条件下填埋体底部至12 m处的垃圾产气速率a不变，且该部分填埋气出流量恒定，故对于12 m以上产气速率a不变的垃圾气压分布不变。由上述分析可见，填埋体内设置水平导气层可有效降低填埋体内的气压；相同抽气速率条件下，填埋体内水平导气层埋深越深降压效果越好；本文模型可指导水平导气层位置和间距布置及合理选择抽气功率。

图 4 填埋体内设置水平导气层对气压分布的影响

2.4 填埋体底部填埋气导排通道对气压分布的影响

国外低水位填埋场底部渗沥液导排系统兼作填埋气导排通道[7, 17-18]。下面分析填埋体底部填埋气导排通道集气速率q对气压分布的影响。由图 5可见，填埋气气压减小值随集气速率q的增大而增大。集气速率为2.0×10-6kg/(m2·s)时，底部填埋气气压由无气体导排时的2.05 kPa下降至0.43 kPa，下降幅度为79%；而埋深12 m处填埋气气压下降幅度仅为28%。其最大气压则是由无气体导排时的2.05 kPa下降至0.93 kPa，降幅为55%，且其位置由底部上升至埋深21 m处。由上述分析可见，填埋体底部填埋气导排通道可有效减小填埋体内填埋气气压，特别是接近导排通道的深部垃圾填埋气气压。

图 5 填埋体底部渗沥液导排系统兼作填埋气导排通道对气压分布的影响

3 结论

提出了分层垃圾填埋体气压分布一维稳态分析模型，并提出顶部边界为给定气压、底部边界为给定气压或流量工况的求解方法，分析了垃圾分层特征、封顶覆盖层下和填埋体内高渗透性的水平导气层、填埋体底部渗沥液导排系统兼作填埋气导排通道对填埋场气压分布的影响，并得到如下结论：填埋体均质化模型(产气速率和固有渗透系数按埋深平均)会显著高估填埋体内气压，填埋体气压分析宜采用分层模型；封顶覆盖层下和填埋体内设置水平导气层可有效减小填埋体气压，在相同抽气速率条件下水平导气层埋深越深降压效果越好；填埋体底部填埋气导排通道对填埋气气压较大的深部垃圾降压效果较明显；本文模型可指导填埋气导气系统布置及合理选择抽气功率。

参考文献
[1] El F M, Findikakis A N, Leckie J O. Environmental impacts of solid waste landfilling[J]. Journal of Environmental Management, 1997, 50: 1–25.
[2] Manfredi S, Christensen T H. Environmental assessment of solid waste land?lling technologies by means of LCA-modeling[J]. Waste Management, 2009, 29: 32–43. DOI:10.1016/j.wasman.2008.02.021
[3] Chen Y M, Xu X B, Zhan L T. Analysis of solid-liquid-gas interactions in landfilled municipal solid waste by a bio-hydro-mechanical coupled model[J]. Science China Technological Sciences, 2012, 55: 81–89. DOI:10.1007/s11431-011-4667-7
[4] Findikakis A N, Leckie J O. Numerical simulation of gas flow in sanitary landfills[J]. ASCE Journal of Environmental Engineering, 1979, 115: 927–945.
[5] El F M, Findikakis A N, Leckie J O. Numerical modeling of generation and transport of gas and heat in landfills I. model formulation[J]. Waste Management & Research, 1996, 14: 483–504.
[6] El F M, Findikakis A N, Leckie J O. Numerical modeling of generation and transport of gas and heat in landfills Ⅱ. model application[[J]. Waste Management & Research, 1996, 14: 537–551.
[7] Townsend T G, Wise W R, Jain P. One-dimensional gas flow model for horizontal gas collection systems at municipal solid waste landfills[J]. Journal of Environmental Engineering, 2005, 131(12): 1716–1723. DOI:10.1061/(ASCE)0733-9372(2005)131:12(1716)
[8] 魏海云, 詹良通, 陈云敏, 等. 分层填埋垃圾体中气体一维稳态运移规律[J]. 岩土工程学报, 2009, 31(11): 1665–1671.
Wei H Y, Zhan L T, Chen Y M, et al. One-dimensional gas migration in multi-layered landfills of municipal solid wastes[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(11): 1665–1671. DOI:10.3321/j.issn:1000-4548.2009.11.004 (in Chinese)
[9] Jain P, Powell J, Townsend T G, Reinhart D R. Air permeability of waste in a municipal solid waste landll[J]. Journal of Environmental Engineering, 2005, 131(11): 1565–1573. DOI:10.1061/(ASCE)0733-9372(2005)131:11(1565)
[10] 彭绪亚, 余毅, 刘国涛. 垃圾填埋场竖井抽气条件下的填埋气压力分布[J]. 重庆大学学报, 2003, 26(1): 92–95.
Peng X Y, Yu Y, Liu G T. Pressure distribution of LFG in a landfill with vertical extraction system[J]. Journal of Chongqing University, 2003, 26(1): 92–95. DOI:10.11835/j.issn.1000-582X.2003.01.024 (in Chinese)
[11] 刘磊, 薛强, 梁冰, 等. 垃圾填埋场封场后气体产出及释放规律研究[J]. 环境卫生工程, 2010, 18(1): 19–21.
Liu L, Xue Q, Liang B, et al. Landfill gas generation and emission in closed landfill sites[J]. Environmental Sanitation Engineering, 2010, 18(1): 19–21. (in Chinese)
[12] BEAR J. 多孔介质流体动力学[M]. 李竞生, 陈崇希. 译. 北京: 中国建筑工业出版社, 1983.
[13] Arigala S G, Tsotsis T T, Webster I A, et al. Gas generation, transport, and extraction in landfills[J]. Journal of Environmental Engineering, 1995, 121(1): 33–44. DOI:10.1061/(ASCE)0733-9372(1995)121:1(33)
[14] Chen Y C, Wu C H, Hu H Y. Numerical simulation of gas emission in a sanitary landfill equipped with a passive venting system[J]. Journal of Environmental Science and Health, 2000, A35(9): 1735–1747.
[15] Tinet A J, Oxarango L. Stationary gas flow to a vertical extraction well in MSW landfill considering the effect of mechanical settlement on hydraulic properties[J]. Chemical Engineering Science, 2010, 65: 6229–6237. DOI:10.1016/j.ces.2010.09.007
[16] Jung Y, Imhoff P T, Augenstein D C, et al. Influence of high-permeability layers for enhancing landfill gas capture and reducing fugitive methane emissions from landfills[J]. Journal of Environmental Engineering, 2009, 135(3): 138–146. DOI:10.1061/(ASCE)0733-9372(2009)135:3(138)
[17] Townsend T G, Miller W L. Landfill gas extraction from leachate collection systems[J]. Journal of Solid Waste Technology and Management, 1997, 24(3): 131–136.
[18] Fredenslund A M, Scheutz C, Kjeldsen P. Tracer method to measure landfill gas emissions from leachate collection systems[J]. Waste Management, 2010, 30(11): 2146–2152. DOI:10.1016/j.wasman.2010.03.013
[19] USEPA, 1985. Minimum technology guidance on double liner systems for landfills and surface impoundment-design, construction and operation[S]. EPA/530/SW-85-012, U.S. Environmental Protection Agency, Cincinnati, OH. https://nepis.epa.gov/Exe/ZyPURL.cgi?Dockey=2000E0AU.TXT
[20] 魏海云. 城市生活垃圾填埋场气体运移规律研究[D]. 杭州: 浙江大学, 2007: 156-157.
[21] Thomas H R, Ferguson W J. A fully coupled heat and mass transfer model incorporating contaminant gas transfer in an unsaturated porous medium[J]. Computers and Geotechnics, 1999, 24: 65–87. DOI:10.1016/S0266-352X(98)00030-5
[22] Sanchez R, Tsotsis T T, Sahimi M. Computer simulation of gas generation and transport in landfills.Ⅳ:Modeling of liquid-gas flow[J]. Chemical Engineering Science, 2010, 65: 1212–1226. DOI:10.1016/j.ces.2009.09.076
[23] Yu L, Batlle F, Carrera J, et al. Gas flow to a vertical gas extraction well in deformable MSW landfills[J]. Journal of Hazardous Materials, 2009, 168: 1404–1416. DOI:10.1016/j.jhazmat.2009.03.045
X