河北科技大学环境科学与工程学院//河北省污染防治生物技术重点实验室，河北
石家庄 050018
摘要：为揭示垃圾渗滤液中污染物在包气带中的迁移转化规律，通过静态吸附、静态降解、动态土柱和数学模型预测等方法进行模拟实验。结果表明，包气带对污染物的吸附过程是线性的，即S=KdC，吸附系数Kd＝0.0976；降解曲线符合一级动力学方程，即C=C0e-λt，降解系数λ=0.0324 d-1；弥散过程符合对流－弥散迁移转化模型，弥散系数D＝0.00435 m2·d-1，由此确定了污染物迁移数学模型。通过动态土柱实验验证了模型的可靠性，并利用模型对垃圾渗滤液中有机污染物(COD)的时空分布进行预测。
关键词：垃圾渗滤液；模拟实验；包气带；迁移转化；预测
中图分类号：X131.3
文献标识码：A
文章编号：1672-2175（2006）05-0960-06
垃圾填埋渗滤液是指垃圾在堆积和填埋过程中由于发酵和雨水的淋浴、冲刷以及地表水和地下水的浸泡而滤出的污水。渗滤液中各类污染物经过包气带时，与土壤介质发生各种物理化学及生物作用，即发生复杂的迁移转化。渗滤液对包气带的影响不仅仅存在于表层，而且贯穿于数十米深并地下水造成严重污染。因此研究地下水污染问题，不仅应当分别研究污染物在表土层、下包气带土层以及含水层各单元中的迁移、转化规律，而且还应将上述单元作为地下水污染系统整体来考虑，研究污染物在整体系统的迁移转化，达到防止地下水污染的目的[1,2]。目前国内外研究环境污染物的转化和归宿问题，多采用实验室模拟，建立数学模型后，再通过大量现场实测数据进行验证或修正的方法[3-4]。
随着我国经济发展和人民生活水平的提高，垃圾填埋场渗滤液已成为表层土壤、浅层地下水污染的重要因素之一。为揭示垃圾渗滤液中有机污染物在包气带中的迁移转化规律，本文以石家庄峡石沟垃圾填埋场为例，通过室内模拟实验研究垃圾渗滤液在包气带中迁移和作用机理；确定污染物在地下环境中迁移模型，拟合实验结果求出模型参数，并通过模型来预测污染物在包气带中的时空分布。为最终提出污染物迁移的控制措施和方案提供理论依据。
1
室内模拟实验
实验采集的水样取自石家庄峡石沟垃圾填埋场，具有较强的刺激性气味，浑浊，呈黄褐色。试验用土样采集于石家庄市郊区，距地面1～2.5 m的包气带。样品采回后，在20 ℃～25 ℃室温下，将其均匀铺于地面上自然风干，将其全部过20目筛。
1.1
垃圾渗滤液中NH4-N的变化规律实验
氨氮以游离氨或氨盐形式存在于水中，两者的组成比取决于水的pH值[5]。当pH值偏高时，游离氨的比例较高。反之，则氨盐的比例为高。在无氧环境中，水中存在的亚硝酸盐可受微生物作用，还原为氨。在有氧环境中，水中氨亦可转变为亚硝酸盐、甚或继续转变为硝酸盐。
1.1.1
实验原理和条件
氨气敏电极为复合电极，以pH玻璃电极为指示电极，银—氯化银电极为参比电极。
(1)仪器和电极的准备：按使用说明书进行调试仪器。
(2)标准曲线的绘制：吸取10.00 mL浓度为0.1、1.0、10、100、1000 mg·L-1的铵标准溶液于25 mL小烧杯中，浸入电极后加入1.0 mL氢氧化钠溶液，在搅拌下，读取稳定的电位值(1 min内变化不超过1 mV时，即可读数)。在半对数坐标线上绘制E-lgc的标准曲线。
(3)水样的测定：吸取10.00 mL水样，以下步骤与标准曲线绘制相同。由测得的电位值，在标准曲线上直接查得水样中的氨氮含量(mg·L-1)。
1.1.2
结果与分析
实验曲线，如图1（下页）所示。
1#瓶中氨氮含量变化曲线：
2#瓶中氨氮含量变化曲线：
由图1可直观的看出NH4-N的变化。最初采集的垃圾渗滤液中NH4-N含量较高，基本上全部超标。而将其放在培养瓶中，使之在厌氧的条件下观察其含量的变化，结果发现在初始阶段NH4-N含量有明显的变化，由785.60 mg·L-1降为323.3 mg·L-1，此后一段时间，其含量没有明显波动。这是因为在起始状态，由于有氧的参与，NH4-N发生了硝化和反硝化，同时由于土壤的吸附，使得NH4-N含量有很大程度的变化；之后，由于处于厌氧的条件下，对于NH4-N来说，硝化和反硝化都无法进行，仅存在土壤的吸附，所以其含量变化不明显。
Fig. 1
NH4-N contents at different time
|
1.2
垃圾渗滤液中COD迁移转化规律实验
Fig. 2
Static adsorption diagram
|
在本研究中模拟实际的土地处理系统下包气带土层，通过实验室条件下的静态弥散、静态吸附、静态降解确定COD在包气带的迁移转化模型[6]，对其迁移转化规律进行研究。
1.2.1
静态吸附实验
1 实验条件
(1)水土比为1∶3(质量比)，因为土样吸附量比较小，故选择比较大的土水比，土样质量10 g，水样是30 mL。
(2)水样浓度21.8～654 mg·L-1。
(3)摇床振荡频率30次·min-1，振幅2 cm，振荡时间8 h。
(4)土水分离条件：置于锥形瓶混合振荡8 h后，取其上清液进行离心分离5 min。再对上清液进行抽滤(微孔滤膜孔径65 μm)，滤液由于有一定浊度，故需将滤液取出10 mL进行蒸馏，取馏出液5 mL测定COD值。
1.2.2
结果与分析
实验曲线，如图2所示：
由图1可以看出，分配系数Kd描述的吸附，其吸附等温线是线性的，即S=KdC
（1），
其中S是土中COD吸附浓度（mg·d-1），C为水中COD值（mg·d-1）。
由计算机拟合曲线方程为：S=0.0976C，因此吸附分配系数Kd＝0.0976
1.3
静态降解实验
1.3.1
实验条件
静态降解实验是在黑暗厌氧条件下，测定土样对COD的降解能力[7]。取两个5 L的培养瓶，模拟包气带土层情况，进行静态降解实验。为消除土壤吸附和自身溶出COD的影响，两个培养瓶，1#瓶蒸馏水作空白，2#瓶为渗滤液。土样比为1∶2，取土样1500 g，渗滤液3000 mL，其COD值为1127.69 mg·L-1。培养开始后，每隔一定时间取30 mL左右的上清液，用微孔滤膜抽滤后，测定COD值，同时测定1#培养瓶中的COD值。
1.3.2
结果与分析
在数据处理时，将1#瓶中的COD值减掉。在实验的第23天1#瓶的COD已降为零，实验所得数据用图3表示如下：
降解曲线符合一级动力学方程，即
C=C0e-λt
（2）
式中 —一级动力学降解系数（1·d-1）
C0——COD初始浓度（mg·L-1）
t——时间（d）
则由计算机拟合曲线方程为：
C/C0=0.9148e-0.0324x
（3）
Fig. 3
COD static degradation test diagram
|
所以，实验条件下的静态降解系数λ为：
1.4
弥散实验
1.4.1
实验条件
（1）装柱：实验土样的容重为1.5 g·cm-3，柱高为31.2 cm，内径8.12 cm，重量为1594 g。底部装上一层粗沙，底层垫上四层纱布。
（2）洗柱：用蒸馏水进行淋滤，直至淋滤液无杂质离子。用EDTA进行检验。如果加入EDTA溶液变红，表明淋滤液还未洗净，需要继续淋滤。如果加入EDTA溶液不变红，表明淋滤液已经洗净。
（3）将浓度C0为0.05 mol·L-1的示踪剂NaCl从上部连续注入土柱，测定不同时刻，取样口(x=15.6 cm)处的示踪剂NaCl浓度，作出C-t曲线即穿透曲线。
Fig. 5
COD concentration diagram of dynamic soil column test
|
1.4.2
实验结果与讨论
弥散实验曲线，如图4所示：
Fig. 4
C-t penetration diagram of dispersion test
|
对于一维问题，在不考虑吸附和降解作用下，污染物在土壤中的运移符合对流－弥散迁移转化模型[8-9]，其方程可以写作：
（4）
式中：c——水中COD值(mg·L-1)
x——垂向运移距离(m)
D——弥散系数(m2·d-1)
v——垂直方向渗透速度(m·d-1)
t——时间（d）
弥散系数D也可以根据穿透曲线，用下述公式计算：
（5）
式中，v为水的实际运移速度（m·d-1）； 分别为c/c0等于0.16，0.5，0.84的时间。用此式求出的弥散系数为D＝0.00435 m2·d-1。
1.5
动态土柱实验
1.5.1
实验条件
在动水条件下进行垃圾渗滤液在包气带中的迁移转化模拟实验，即动态土柱实验[10]。测定弥散系数后的土柱，再经蒸馏水洗净后，用污水淋滤，实验水样用恒水头装置自土柱上端通入，由下端出口流出，隔一段时间接一次淋滤液，测定其COD。
1.5.2
实验结果与讨论
由于土柱经过彻底清洗，柱内微生物量很少，加之实验时间较短，曲线（图5）主要表现的是弥散和吸附的共同作用。
2
数学模型的确定
2.1
数学模型的确立和其解析解
采用何种模型进行数学模拟，主要取决于研究问题的目的，所要求的精度和已有或可能取得的基础资料。渗滤液中COD的变化是污水在包气带土层中经过各种物理、化学和生物反应对其迁移的综合影响形成的。考虑污染物质主要是沿垂向运移，所以其运移模型常按垂向一维问题处理，一般认为水在土层中运移符合推流模式，依据以上实验的机理分析，则污染物质在土层中垂直向下迁移的基本方程为[11-15]：
| 
