一种测试垃圾双孔隙度渗流参数的装置及方法

1.本发明涉及测量领域，具体涉及一种测试垃圾双孔隙度渗流参数的装置及方法。


背景技术：

2.我国城市及周边分布着2.8万余座生活垃圾填埋场，存量垃圾达60多亿吨。垃圾生化降解产生大量渗滤液和填埋气容易诱发填埋体滑坡、地下水土及大气污染等环境灾害，这些灾害的产生源头与液气在垃圾孔隙中的运移过程密切相关，因此准确认识生活垃圾的水气渗流特性是提高填埋场灾害防控水平的基础。
3.垃圾渗流定义为水(即渗滤液)和气(填埋气)在垃圾孔隙中的流动。和一般均质岩土体不同，生活垃圾组分多样，复杂的颗粒形状及不良的颗粒级配引起多尺寸孔隙并存。垃圾渗流影响因素包括垃圾的颗粒级配、孔隙比、含水率、水或者气本身的粘滞度等。渗滤液在垃圾堆体内的运移具有明显的不均质性，大孔隙区域为优势流，小孔隙区域为基质流，可以建立双孔隙度模型或单孔隙度模型分析垃圾堆体水气运移过程。传统单孔隙度渗流理论仅适用于均质岩土体，无法准确描述垃圾填埋场水气优势流行为。现有研究已经表明双孔隙度模型的精确度显著高于单孔隙度模型，例如ke等(analytical solution of leachate flow to vertical wells in municipal solid waste landfills using a dual
‑
porosity model.[j]engineering geology,2018,239,27
‑
40.)对上海某填埋场的竖井降水试验进行数值模拟，距离4.4m处的实测水位降深为1.66m，单孔隙度模型模拟计算水位降深(4.23m)远高于实测值，而双孔隙度模型模拟计算值(1.43m)与实测值较为接近(如图1所示)。
[0004]
因此需要建立双孔隙度模型来获取渗流参数，双孔隙度模型将垃圾多孔介质分为大孔隙优势流区域及小孔隙基质流区域并赋予不同渗流参数，渗滤液分别在大孔隙区域和小孔隙区域运移，同时大孔隙区域和小孔隙区域之间由于存在水头差发生水分交换。相比于单孔隙度模型，双孔隙度模型虽能更加准确描述流体运移过程，但其计算所需渗流参数较多。单孔隙度模型在进行计算时需要的渗流参数为6个，而双孔隙度模型需要的渗流参数为14个。对于所述双孔隙度模型中的渗流参数的获得一般有物理试验和数值模拟两种方式。数值模拟是对现场或室内的垃圾水气运移试验结果进行反演，通过数值模拟这14个参数里的10余个使得计算结果与测试值相吻合。这也是当前学者们获得双孔隙度渗流参数最常用的方法，如ke等(analytical solution of leachate flow to vertical wells in municipal solid waste landfills using a dual
‑
porosity model.[j]engineering geology,2018,239,27
‑
40.)和liu等(modeling the oxygen transport process under preferential flow effect in landfill.[j]environmental science and pollution research,2018,25,18559
‑
18569.)进行的竖井降水及注气试验无法区分出大孔隙和小孔隙区域的渗流过程，导致在试验过程中是按照单孔隙模型(即把垃圾堆体视作一个整体)来测试数据，在完成试验以后通过数值模拟的方式把测试结果拆分成双孔隙度模型，分别得到大孔隙区域渗流参数和小孔隙区域渗流参数。如ke等通过数值模拟竖井抽水试验确定其
中3个关键参数：大孔隙区域体积占比w
f
、大孔隙区域饱和渗透系数k
swf
、小孔隙区域饱和渗透系数k
swm
，而剩余的模型参数取值也只是进行简单假定或借鉴已有结果，这种数值模拟方式难以保证渗流参数取值的真实性及准确性，也显著降低了双孔隙度渗流模型的可靠性。
[0005]
物理试验相比于数值模拟具有更好的唯一性和可靠性，但是还没有通过物理试验直接确定双孔隙度模型中的渗流参数的测试装置及方法。目前文献中报道的物理试验无法区分出大孔隙和小孔隙区域的渗流过程，试验的直接测试结果仅仅只是大、小孔隙各自渗流参数的混合值，既不能代表大孔隙区域渗流参数，也不能代表小孔隙区域渗流参数。现有技术难以通过物理方式测试获取垃圾双孔隙度渗流参数取值，极大地限制了双孔隙度模型用于垃圾填埋场渗流分析的可行性及实际使用价值。


技术实现要素：

[0006]
为克服现有技术的不足，本发明所要解决的技术问题是提供一种测试垃圾双孔隙度渗流参数的装置及方法，巧妙地利用大、小孔隙区域持水能力的差异性，分阶段获得大、小孔隙区域渗流参数，分别根据自主设计的饱和状态下定水头渗流试验确定垃圾整体渗流参数、根据重力下自由排水试验确定大孔隙区域渗流参数，根据不同气压下分级排水试验确定小孔隙区域渗流参数，通过物理测试的方式直接确定了双孔隙度模型中13个渗流参数的取值，仅大孔隙与小孔隙区域的界面交换系数k
a
通过数值反演得到，减少了需要数值模拟的渗流参数数量，提高了双孔隙度渗流模型的精度和可靠度。
[0007]
本发明所采用的技术方案如下：一种测试垃圾双孔隙度渗流参数的方法，所述的方法包括以下步骤：
[0008]
(1
‑
1)在垃圾双孔隙度渗流参数测试装置中，按指定孔隙比分层均匀填充垃圾，获得指定孔隙比的垃圾柱，直接测量所述垃圾柱的体积v
waste
；
[0009]
(1
‑
2)进行饱和状态下定水头渗流试验，从底部通入气体排出所述垃圾柱孔隙中的空气，然后在所述测试装置顶部施加一定范围的负压；从所述测试装置底部对所述垃圾柱进行缓慢注水，当整个垃圾柱浸没在水中后即视为饱和状态，停止注水，直接测量得到饱和阶段累积注水量v
total
；在所述垃圾柱一直处于饱和状态下，在不同水头下(所述水头均高于所述垃圾柱顶部)，直接测量所述测试装置底部的出流速度v
total
以及相应水力梯度i
total
；作为优选，所述一定范围负压为20～50kpa；
[0010]
(1
‑
3)在重力下进行自由排水试验，实时直接测量自由排水阶段初始时刻(饱和状态下)所述测试装置底部的出水流速v
f0
及其渗流梯度i
f0
，自由排水阶段后续时刻(非饱和状态下)测得的所述出水流速v
ft
、渗流梯度i
ft
以及相应的大孔隙吸力，大孔隙体积含水量θ
wf
，以及自由排水阶段累积出流量v
f
；
[0011]
(1
‑
4)进行不同气压下分级排水试验，对所述的垃圾柱分级增加气压，直接测量不同气压下所述测试装置底部的累计出流量v
m
、出流速度v
m
，渗流梯度i
m
，相应的小孔隙吸力以及小孔隙体积含水量θ
wm
。
[0012]
作为优选，所述的孔隙比为所述垃圾柱中孔隙体积与其固体颗粒体积之比，反映了所述垃圾柱的密实程度，孔隙比可以通过现场垃圾实际测得；不同学者根据各自研究的实际情况对于垃圾的大小孔隙的划分界限不一样，一般以30
‑
145μm为界，其中直径在此界限及以上的为大孔隙，直径在此界限以下为小孔隙。
[0013]
作为优选，所述水头为计算点距基准面的高度，在本发明中指的是水位面距离所述垃圾柱底部的高度。所述水头的计算方法为：大孔隙水头h
f
＝大孔隙吸力的负数/每米水头的压力，小孔隙水头h
m
＝小孔隙吸力的负数/每米水头的压力；每米水头的压力10kpa/m。
[0014]
作为优选，利用步骤(1
‑
2)的测量结果，获得垃圾整体渗流参数，方法如下：
[0015]
(3
‑
1)通过式(1)获得垃圾整体饱和体积含水量θ
ws
；
[0016][0017]
(3
‑
2)根据测试装置底部的出流速度v
total
和相应水力梯度i
total
的测量结果，通过式(2)获得整体饱和渗透系数k
sw
；
[0018][0019]
作为优选，利用步骤(1
‑
2)、(1
‑
3)的测量结果，获得垃圾大孔隙区域渗流参数取值，方法如下：
[0020]
(4
‑
1)根据研究(fellner等，2010，modeling of leachate generation from msw landfills by a 2
‑
dimensional 2
‑
domain approach.waste management,30,2084
‑
2095)表明：在饱和状态下认定所述大孔隙全部被水占据，因此定义大孔隙区域饱和体积含水量θ
wsf
＝1；同时在自由排水阶段认定大孔隙中的水可以全部渗流排出，因此定义大孔隙区域残余体积含水量θ
wrf
＝0；
[0021]
(4
‑
2)在自由排水阶段初始时刻，可认为大孔隙区域仍处于饱和状态，根据自由排水阶段初始时刻测得的所述出水流速v
f0
及其渗流梯度i
f0
，通过式(3)计算得到大孔隙区域饱和渗透系数k
swf
；
[0022][0023]
(4
‑
3)通过式(4)计算得到大孔隙区域体积占比w
f
；
[0024][0025]
(4
‑
4)根据自由排水阶段获得的所述大孔隙体积含水量θ
wf
和水头h
f
，通过式(5)，采用最小二乘法(matlab软件)拟合得到大孔隙区域进气值的倒数α
f
及决定持水曲线斜率的孔隙尺寸分布参数n
f
；
[0026][0027]
式中：θ
wf
为大孔隙区域体积含水量，θ
wrf
为大孔隙区域残余体积含水量，θ
wsf
为大孔隙区域饱和体积含水量，s
ewf
为大孔隙区域有效饱和度，h
f
为大孔隙区域的水头(单位：m)，α
f
为大孔隙区域进气值的倒数(单位：m
‑1)，n
f
为大孔隙区域决定持水曲线斜率的孔隙尺寸分布参数，m
f
＝1
‑
1/n
f
；
[0028]
(4
‑
5)在自由排水阶段后续时刻，大孔隙区域处于非饱和状态，根据自由排水阶段后续时刻测得的所述出水流速v
ft
及其渗流梯度i
ft
，计算得到大孔隙区域非饱和渗透系数
k
wf
＝v
ft
/i
ft
，通过式(6)，采用最小二乘法(matlab软件)拟合得到大孔隙区域连通性参数l
f
；
[0029][0030]
式中：k
wf
为大孔隙区域非饱和渗透系数(单位：m/s)；l
f
为大孔隙区域连通性参数；
[0031]
作为优选，利用步骤(1
‑
4)的测量结果，获得垃圾小孔隙区域渗流参数，方法如下：
[0032]
(5
‑
1)通过式(7)计算得到小孔隙区域饱和渗透系数k
swm
；
[0033][0034]
(5
‑
2)通过式(8)计算得到小孔隙区域饱和体积含水量θ
wsm
；
[0035][0036]
(5
‑
3)在最后一次施加气压排水后所述垃圾柱处于残余含水量状态，对所述垃圾柱烘干，直接测量得到小孔隙区域残余体积含水量θ
wrm
取值；
[0037]
(5
‑
4)根据分级排水阶段获得的所述小孔隙体积含水量θ
wm
和水头h
m
，通过式(9)，采用最小二乘法(matlab软件)拟合得到小孔隙区域进气值的倒数α
m
及决定持水曲线斜率的孔隙尺寸分布参数n
m
；
[0038][0039]
式中：θ
wm
为小孔隙区域体积含水量，θ
wrm
为小孔隙区域残余体积含水量，θ
wsm
为小孔隙区域饱和体积含水量，s
ewm
为小孔隙区域有效饱和度，h
m
为小孔隙区域的水头(单位：m)，α
m
为小孔隙区域进气值的倒数(单位：m
‑1)，n
m
为小孔隙区域决定持水曲线斜率的孔隙尺寸分布参数，m
m
＝1
‑
1/n
m
；
[0040]
(5
‑
5)根据分级排水阶段测得的所述出水流速v
m
及其渗流梯度i
m
，计算得到小孔隙区域非饱和渗透系数k
wm
＝v
m
/i
m
，通过式(10)，采用最小二乘法(matlab软件)拟合得到小孔隙区域连通性参数l
m
；
[0041][0042]
式中：k
wm
为小孔隙区域非饱和渗透系数(单位：m/s)；l
m
为小孔隙区域连通性参数；
[0043]
作为优选，在确定整体渗流参数、大孔隙区域渗流参数和小孔隙区域渗流参数取值后，基于式(5)、(6)、(9)、(10)、(11)
‑
(13)所述的双孔隙度渗流模型，采用comsol multiphysics软件模拟底部出流量，反演得到大孔隙区域与小孔隙区域界面交换系数k
a
；
[0044]
大孔隙区域和小孔隙区域的渗流控制关系式为：
[0045]
大孔隙：
[0046]
小孔隙：
[0047]
交换项：γ
w
＝k
a
(h
f
‑
h
m
)
ꢀꢀ
(13)
[0048]
式中：t代表时间(单位：s)，z代表高度(单位：m)，г
w
为大孔隙区域与小孔隙区域质量交换项(单位：s
‑1)，k
a
为大小孔隙区域界面交换系数(单位：m
‑1s
‑1)。
[0049]
本发明还公开了一种垃圾双孔隙度渗流参数测试装置，所述测试装置包括装填结构、均匀布水结构、可调节供水供气结构、测量结构，其中：
[0050]
所述装填结构，用于装填垃圾柱的容纳物，相对密封，为所述可调节供水供气结构以及所述测量结构提供接口；
[0051]
所述均匀布水结构，用于给所述垃圾柱均匀供水，使得所述垃圾柱各部分含水量趋于一致；
[0052]
所述可调节供水供气结构，用于给所述垃圾柱注水或注气，来控制所述垃圾柱的含水量和顶部气压；
[0053]
所述测量结构：在饱和状态下定水头渗流试验中，测量饱和阶段累积注水量v
total
；在所述垃圾柱一直处于饱和状态下，在不同水头下(所述水头均高于所述垃圾柱顶部)，测量所述测试装置底部的出流速度v
total
以及相应的水力梯度i
total
；在重力下自由排水试验中，直接测量自由排水阶段初始时刻(饱和状态下)所述测试装置底部的出水流速v
f0
及其渗流梯度i
f0
，自由排水阶段后续时刻(非饱和状态下)测得的所述出水流速v
ft
、渗流梯度i
ft
以及相应的大孔隙吸力，大孔隙体积含水量θ
wf
，以及自由排水阶段累积出流量v
f
；在不同气压下分级排水试验中，测量不同气压下累计出流量v
m
、出流速度v
m
，渗流梯度i
m
，相应的小孔隙吸力以及小孔隙体积含水量θ
wm
；
[0054]
作为优选，所述装填结构包括：容器、气腔、密封盖；所述容器用于装填垃圾柱；所述气腔用于注气时缓冲气压以及注水时调节水头；所述密封盖，用于盖住所述容器，使得所述装填装置内的含水量和气压保持一定值，不受外界影响。
[0055]
作为优选，所述均匀布水结构包括：多孔板和碎石层；所述碎石层用于防止渗滤液流动时带动垃圾中细颗粒堵塞多孔板，并使得注水注气时流体能均匀分布在所述垃圾柱表面；所述垃圾柱顶部和底部碎石层铺设于所述多孔板与所述垃圾柱之间。
[0056]
作为优选，所述可调节供水供气结构包括：供水器、水阀门，注气器、调压阀、流量计和压力计；其中所述供水器和水阀门用于调节供水量；所述注气器、调压阀、流量计和压力计用于调节注气量；所述可调节供水供气结构通过管道与所述装填结构连接。
[0057]
作为优选，所述测量结构包括：时域反射仪(tdr)、张力计、储水箱和水量计；所述tdr用于监测垃圾柱的含水率；所述张力计用于监测垃圾柱吸力；所述储水箱用于储存从垃圾柱中渗流出来的水；所述水量计用于测量储水箱中的水量；在所述装填结构侧壁不同高度打孔，所述tdr、张力计和所述孔通过密封圈进行连接；所述储水箱与所述装填结构底部通过出水管路连接。
[0058]
作为优选，在所述垃圾柱顶部和底部均铺满碎石层及多孔板。
[0059]
作为优选，所述碎石层厚度为500～1000mm之间，所述多孔板的开孔大小为8～30mm之间。
[0060]
本发明的有益效果是：
[0061]
1)本发明巧妙地利用大、小孔隙区域持水能力的差异性，分阶段获得大、小孔隙区域渗流参数，分别根据自主设计的饱和状态下定水头渗流试验确定垃圾整体渗流参数、根
据重力下自由排水阶段试验确定大孔隙区域渗流参数，根据不同气压作用下分级排水试验确定小孔隙区域渗流参数，通过物理测试的方式直接确定了双孔隙度模型中13个渗流参数的取值，仅大孔隙与小孔隙区域的界面交换系数k
a
通过数值反演得到，减少了需要数值模拟的渗流参数数量，提高了双孔隙度模型的精度和可靠度。而目前文献中报道的物理试验无法区分出大孔隙和小孔隙区域的渗流过程，试验的直接测试结果仅仅只是大、小孔隙各自渗流参数的混合值，既不能代表大孔隙区域渗流参数，也不能代表小孔隙区域渗流参数，极大地限制了垃圾双孔隙度模型用于垃圾填埋场渗流分析的可行性及实际使用价值。本发明克服了垃圾双孔隙度渗流参数无法直接通过物理测试方式确定的技术瓶颈，提升了双孔隙度模型用于填埋场渗流分析的可行性及实际使用价值；
[0062]
2)由于现场不同深度填埋垃圾的应力水平和降解程度相差较大，导致垃圾的优先流程度在不同深度存在差异；本发明通过测试装置内填充不同孔隙比及降解龄期的垃圾柱，用以研究应力压缩和降解作用对垃圾双孔隙度渗流参数的影响，为准确分析不同深度的渗流过程提供依据；
[0063]
3)测试装置操作简便，同一垃圾柱在整个测试过程无需拆卸，避免了拆卸过程对垃圾柱孔隙结构产生的扰动。
附图说明
[0064]
图1是分别用单孔隙模型和双孔隙模型对竖井降水进行模拟的结果图。
[0065]
图2是本发明所述垃圾双孔隙度渗流参数测试装置的结构示意图；
[0066]
图3是本发明所述测试垃圾双孔隙度渗流参数的方法的流程示意图；
[0067]
图中：1、容器，2、垃圾柱，3、气腔，4、多孔板，5、碎石层，6、注气器，7、调压阀，8、流量计，9、压力计，10、阀门，11、tdr传感器，12、张力计，13、密封圈，14、密封盖，15、供水器，16、储水箱，17、水量计。
具体实施方式
[0068]
下面结合附图和具体实施例对本发明作进一步详细说明。
[0069]
如图2所述，本发明公开了一种垃圾双孔隙度渗流参数测试装置，所述测试装置包括装填结构、均匀布水结构、可调节供水供气结构、测量结构，其中：
[0070]
所述装填结构，用于装填垃圾柱的容纳物，相对密封，为所述可调节供水供气结构以及所述测量结构提供接口；
[0071]
所述均匀布水结构，用于给所述垃圾柱均匀供水，使得所述垃圾柱各部分含水量趋于一致；
[0072]
所述可调节供水供气结构，用于给所述垃圾柱注水或注气，来控制所述垃圾柱的含水量和顶部气压；
[0073]
所述测量结构：在饱和状态下定水头渗流试验中，测量饱和阶段累积注水量v
total
；在所述垃圾柱一直处于饱和状态下，在不同水头下(所述水头均高于所述垃圾柱顶部)，测量所述测试装置底部的出流速度v
total
以及相应的水力梯度i
total
；在重力下自由排水试验中，实时直接测量自由排水阶段初始时刻(饱和状态下)所述测试装置底部的出水流速v
f0
及其渗流梯度i
f0
，自由排水阶段后续时刻(非饱和状态下)测得的所述出水流速v
ft
、渗流梯度
i
ft
以及相应的大孔隙吸力，大孔隙体积含水量θ
wf
，以及自由排水阶段累积出流量v
f
；在不同气压下分级排水试验中，测量不同气压下累计出流量v
m
、出流速度v
m
，渗流梯度i
m
，相应的小孔隙吸力以及小孔隙体积含水量θ
wm
。
[0074]
作为优选，所述装填结构包括：容器1、气腔3、密封盖14；所述容器用于装填垃圾柱；所述气腔用于注气时缓冲气压以及注水时调节水头用；所述密封盖，用于盖住所述容器，使得所述装填装置内的含水量和气压保持一定值，不受外界影响。
[0075]
作为优选，所述均匀布水结构包括：多孔板4和碎石层5；所述碎石层用于防止渗滤液流动时带动垃圾中细颗粒堵塞多孔板，并使得注水注气时流体能均匀分布在所述垃圾柱表面。所述垃圾柱顶部和底部碎石层铺设于所述多孔板与所述垃圾柱之间。
[0076]
作为优选，在所述垃圾柱顶部和底部均铺满碎石层及多孔板。
[0077]
作为优选，所述碎石层厚度为500～1000mm之间，所述多孔板的开孔大小为8～30mm之间，所述顶部和底部碎石层铺设于所述多孔板与所述垃圾柱之间。
[0078]
作为优选，所述可调节供水供气结构包括供水器15、水阀门10，注气器6、调压阀7、流量计8和压力计9；其中所述供水器和水阀门用于调节供水量；所述注气器、调压阀、流量计和压力计用于调节注气量；所述可调节供水供气结构通过管道与所述装填结构连接。
[0079]
作为优选，所述测量结构包括时域反射仪(tdr)11、张力计12、储水箱16和水量计17；所述tdr用于监测垃圾柱的含水率；所述张力计用于监测垃圾柱吸力；所述储水箱用于储存从垃圾柱中渗流出来的水；所述水量计用于测量储水箱中的水量；在所述装填结构侧壁不同高度打孔，所述tdr、张力计和所述孔通过密封圈进行连接；所述储水箱与所述装填结构底部通过出水管路连接。
[0080]
具体地，本发明提供的垃圾双孔隙度渗流参数测试装置可以是如图2所示的结构。考虑到垃圾颗粒尺寸较大，容器1的内径取为400mm，高度为1200mm，壁厚8mm；垃圾柱2需进行多级出流试验，即上边界在不同气相压力条件下测试渗滤液出流情况，因而容器上部为气腔3，气腔高度设为200mm，注气器6中压缩气体经调压阀7、流量计8及压力计9后进入气腔，调压阀7通过控制进气量的多少对气腔压力进行调节；容器1下部为分层填充垃圾柱2，高度为1000mm，初始孔隙比e控制在2.5；垃圾柱2顶部及底部为厚度25mm的不锈钢多孔板4，开孔大小为20mm，沿径向均匀排列3圈，共37个孔，开孔率为9.3％；多孔板与垃圾柱之间碎石层5厚度为300mm，其有两个作用：一、防止渗滤液流动时带动垃圾中细颗粒堵塞多孔板4；二、使得注水注气时流体能均匀分布在垃圾柱2表面；试验过程中主要对垃圾柱的含水量、吸力及底部出水流量进行实时连续监测，其中含水量和吸力通过埋设在不同高度的tdr传感器11和张力计12进行测量，共布设有row a～c三排传感器，其距垃圾柱2底部的距离分别为700,500及300mm；tdr传感器11和张力计12是在垃圾分层填充完成后，在侧壁钻孔进行埋设，埋设后采用密封圈13密封；容器底部出流连接储水箱16，累积出流量通过水量计17进行测量；上部供水器15对垃圾柱2起供水作用，通过阀门10进行控制，可进行饱和及注水等过程。
[0081]
如图3所示，本发明还公开了一种测试垃圾双孔隙度渗流参数的方法，所述的方法包括以下步骤：
[0082]
(1
‑
1)准备阶段，在垃圾双孔隙度渗流参数测试装置中，按指定孔隙比分层均匀填充垃圾2，获得指定孔隙比的垃圾柱，填充过程中在相应高度埋设tdr传感器11和张力计12，
采用注气方式检查所述装置气密性，直接测量所述垃圾柱的体积v
waste
；
[0083]
(1
‑
2)进行饱和状态下定水头渗流试验：从所述测试装置底部以5l/min的速率通入co2气体，用以排出所述垃圾柱孔隙中的空气；然后在所述测试装置顶部施加20kpa负压，在所述测试装置底部注入去离子水进行饱和，当整个所述垃圾柱2浸没在水中之后即视为饱和状态，停止注水，同时缓慢释放负压，通过流量计测量饱和阶段累积注水量v
total
；在所述垃圾柱一直处于饱和状态下，通过供水器15向垃圾柱2上部注水，改变供水器15的高度(均高于所述垃圾柱顶部))，直接测量所述测试装置底部的出流速度v
total
以及相应水力梯度i
total
。
[0084]
(1
‑
3)在重力下进行自由排水实验：停止注水，首先将底部出水管管口高度抬升至与顶部多孔板4平齐，用以排出气腔3中多余水；然后将出水管口下放至垃圾柱2底部，开始进行重力下自由排水，即上部气腔3中不施加气压，待排水稳定后停止试验；这里的排水达到稳定不仅指出水流量不再增加，还包括tdr传感器11和张力计12读数不再变化(下同)；自由排水试验可视为大孔隙中水分出流的过程，测量自由排水阶段初始时刻(饱和状态下)所述测试装置底部的出水流速v
f0
及其渗流梯度i
f0
，测量自由排水阶段后续时刻(非饱和状态下)底部出水流速v
ft
、渗流梯度i
ft
以及相应的大孔隙吸力，大孔隙体积含水量θ
wf
，以及自由排水阶段累积出流量v
f
。
[0085]
(1
‑
4)进行不同气压下分级排水试验：对所述的垃圾柱分级增加气压，首先向气腔3中通气至气压为10kpa，底部出水口保持与底板平齐进行排水；待排水达到稳定后，再向气腔3中增大气压至30kpa，底部进行排水；依次类推，共进行五级排水过程，各级上部气压分别为10,30,50,70,90kpa；分级排水试验实质是保持底部边界条件不变，通过改变顶部气压边界条件，使垃圾柱2中小孔隙区域水分逐级排出，直接测量不同气压下所述测试装置底部的累计出流量v
m
、出流速度v
m
，渗流梯度i
m
，相应的小孔隙吸力以及小孔隙体积含水量θ
wm
。
[0086]
作为优选，所述水头为计算点距基准面的高度，在本发明中指的是水位面距离所述垃圾柱底部的高度。所述水头的计算方法为：大孔隙水头h
f
＝大孔隙吸力的负数/每米水头的压力，小孔隙水头h
m
＝小孔隙吸力的负数/每米水头的压力；每米水头的压力10kpa/m。
[0087]
作为优选，利用步骤(1
‑
2)的测量结果，获得垃圾整体渗流参数，方法如下：
[0088]
(3
‑
1)通过式(1)获得垃圾整体饱和体积含水量θ
ws
；
[0089][0090]
(3
‑
2)根据测试装置底部的出流速度v
total
和相应水力梯度i
total
的测量结果，通过式(2)获得整体饱和渗透系数k
sw
；
[0091][0092]
作为优选，利用步骤(1
‑
2)、(1
‑
3)的测量结果，获得垃圾大孔隙区域渗流参数，方法如下：
[0093]
(4
‑
1)现有的研究(fellner等，2010，modeling of leachate generation from msw landfills by a 2
‑
dimensional 2
‑
domain approach.waste management,30,2084
‑
2095)表明：在饱和状态下认定所述大孔隙全部被水占据，因此定义大孔隙区域饱和体积含水量θ
wsf
＝1；同时在自由排水阶段认定大孔隙中的水可以全部渗流排出，因此定义大孔隙
区域残余体积含水量θ
wrf
＝0；
[0094]
(4
‑
2)在自由排水阶段初始时刻，可认为大孔隙区域仍处于饱和状态，根据自由排水阶段初始时刻测得的所述出水流速v
f0
及其渗流梯度i
f0
，通过式(3)计算得到大孔隙区域饱和渗透系数k
swf
；
[0095][0096]
(4
‑
3)通过式(4)计算得到大孔隙区域体积占比w
f
；
[0097][0098]
(4
‑
4)根据自由排水阶段获得的所述大孔隙体积含水量θ
wf
和水头h
f
，通过式(5)，采用最小二乘法(matlab软件)拟合得到大孔隙区域进气值的倒数α
f
及决定持水曲线斜率的孔隙尺寸分布参数n
f
；
[0099][0100]
式中：θ
wf
为大孔隙区域体积含水量，θ
wrf
为大孔隙区域残余体积含水量，θ
wsf
为大孔隙区域饱和体积含水量，s
ewf
为大孔隙区域有效饱和度，h
f
为大孔隙区域的水头(单位：m)，α
f
为大孔隙区域进气值的倒数(单位：m
‑1)，n
f
为大孔隙区域决定持水曲线斜率的孔隙尺寸分布参数，m
f
＝1
‑
1/n
f
；
[0101]
(4
‑
5)在自由排水阶段后续时刻，大孔隙区域处于非饱和状态，根据自由排水阶段后续时刻测得的所述出水流速v
ft
及其渗流梯度i
ft
，计算得到大孔隙区域非饱和渗透系数k
wf
＝v
ft
/i
ft
，通过式(6)，采用最小二乘法(matlab软件)拟合得到大孔隙区域连通性参数l
f
；
[0102][0103]
式中：k
wf
为大孔隙区域非饱和渗透系数(单位：m/s)；l
f
为大孔隙区域连通性参数；作为优选，利用步骤(1
‑
4)的测量结果，获得垃圾小孔隙区域渗流参数，方法如下：
[0104]
(5
‑
1)通过式(7)计算得到小孔隙区域饱和渗透系数k
swm
；
[0105][0106]
(5
‑
2)通过式(8)计算得到小孔隙区域饱和体积含水量θ
wsm
；
[0107][0108]
(5
‑
3)在最后一次施加气压排水后垃圾柱处于残余含水量状态，对所述垃圾柱烘干，直接测量得到小孔隙区域残余体积含水量θ
wrm
取值；
[0109]
(5
‑
4)根据分级排水阶段测得的所述小孔隙体积含水量θ
wm
和水头h
m
，通过式(9)，采用最小二乘法(matlab软件)拟合得到小孔隙区域进气值的倒数α
m
及决定持水曲线斜率的孔隙尺寸分布参数n
m
；
[0110][0111]
式中：θ
wm
为小孔隙区域体积含水量，θ
wrm
为小孔隙区域残余体积含水量，θ
wsm
为小孔隙区域饱和体积含水量，s
ewm
为小孔隙区域有效饱和度，h
m
为小孔隙区域的水头(单位：m)，α
m
为小孔隙区域进气值的倒数(单位：m
‑1)，n
m
为小孔隙区域决定持水曲线斜率的孔隙尺寸分布参数，m
m
＝1
‑
1/n
m
；
[0112]
(5
‑
5)根据分级排水阶段测得的所述出水流速v
m
及其渗流梯度i
m
，计算得到小孔隙区域非饱和渗透系数k
wm
＝v
m
/i
m
，通过式(10)，采用最小二乘法(matlab软件)拟合得到小孔隙区域连通性参数l
m
；
[0113][0114]
式中：k
wm
为小孔隙区域非饱和渗透系数(单位：m/s)；l
m
为小孔隙区域连通性参数；
[0115]
作为优选，在确定整体渗流参数、大孔隙区域渗流参数和小孔隙区域渗流参数取值后，基于式(5)、(6)、(9)、(10)、(11)
‑
(13)所述的双孔隙度渗流模型，采用comsol multiphysics软件模拟底部出流量，反演得到大孔隙区域与小孔隙区域界面交换系数k
a
；
[0116]
大孔隙区域和小孔隙区域的渗流控制关系式为：
[0117]
大孔隙：
[0118]
小孔隙：
[0119]
交换项：γ
w
＝k
a
(h
f
‑
h
m
)
ꢀꢀ
(13)
[0120]
式中：t代表时间(单位：s)，z代表高度(单位：m)，г
w
为大孔隙区域与小孔隙区域质量交换项(单位：s
‑1)，k
a
为大小孔隙区域界面交换系数(单位：m
‑1s
‑1)。
[0121]
根据上述图2所示垃圾双孔隙度渗流参数测试装置以及图3所述测试方法，巧妙地利用大、小孔隙区域持水能力的差异性，分阶段测试了我国典型高厨余含量生活垃圾的大、小孔隙区域渗流参数，得到以下14个双孔隙度渗流参数，具体数值见表1，可见，根据本实施例的结果可知，本发明通过物理测试的方式直接确定了双孔隙度模型中13个渗流参数的取值，仅大孔隙与小孔隙区域的界面交换系数k
a
通过数值反演得到，本发明采用的技术方案能精准确定垃圾双孔隙渗流参数取值，从而提升双孔隙模型的实际使用价值。
[0122]
[0123][0124]
表1垃圾双孔隙度渗流参数取值
[0125]
以上所述的仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的范围之内，所做的任何修改、等效替换等，均应包含在本发明的保护范围之内。