IC-AnMBR处理村镇厨余渗滤液废水稳定运行调控机制

王拓; 郁达伟; 方一; 陈梅雪; 魏源送

doi:10.12030/j.cjee.202001108

IC-AnMBR处理村镇厨余渗滤液废水稳定运行调控机制

1.
中国科学院生态环境研究中心，水污染控制实验室，北京 100085
2.
中国科学院大学，北京 100049
3.
加州大学戴维斯分校，戴维斯 95616，美国

作者简介: 王拓(1995—)，男，硕士研究生。研究方向：膜生物反应器。E-mail：wangtuo17@mails.ucas.edu.cn

通讯作者: 郁达伟(1982—)，男，博士，助理研究员。研究方向：膜生物反应器等。E-mail：dwyu@rcees.ac.cn;

基金项目:
国家重点研发计划(2018YFD1100600)
中图分类号: X505

Regulation mechanism for steady operation of IC-AnMBR treating high strengh organic rural leachate from food waste

1.
Laboratory of Water Pollution Control Technology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
2.
University of Chinese Academy of Sciences, Beijing 100049, China
3.
University of California-Davis, Davis 95616, The United States

Corresponding author: YU Dawei, dwyu@rcees.ac.cn ;

摘要: 村镇厨余垃圾渗滤液等高浓度有机废水的高效处理是提升村镇环境卫生水平的一个重要方面。为满足村镇厨余垃圾渗滤液低能耗有机物排放达标的处理需求，构建了内循环厌氧膜生物反应器 (internal circulation anaerobic membrane bioreactor，IC-AnMBR)，并用来处理厨余渗滤率废水，重点分析了反应器的COD去除性能和调控机制；根据pH、VFAs/碱度、容积产气率、膜通量和出水有机污染物组分等指标，考察了COD在水解酸化、产甲烷和膜截留过程中的转移转化特征。结果表明：通过耦合膜擦洗曝气和沼气曝气循环，将VFAs/碱度和容积产气率分别从1.5和0.1优化到0.02和1.0；优化了COD稳定达标性能和去除负荷，将COD去除率和负荷从59%和0.3 kg·(m³·d)⁻¹分别提高到了97.7%和1.8 kg·(m³·d)⁻¹；采用沼气循环曝气擦洗陶瓷膜，控制了滤饼层积累，并将膜通量从0.6 L·(m²·h)^-1提高到2.1 L·(m²·h)⁻¹。IC-AnMBR短流程工艺能够实现村镇厨余垃圾渗滤液的稳定处理。
- 厨余垃圾渗滤液 /
- 内循环厌氧膜生物反应器(IC-AnMBR) /
- COD去除率 /
- 膜通量
Abstract: Efficient handling high strength organic wastewater such as food waste leachate is an important aspect for improving rural environmental sanitation level. An internal circulation anaerobic membrane bioreactor (IC-AnMBR) was constructed to treat the wastewater from rural food waste leachate and meet its organic discharge requirements with low energy consumption. The COD removal performance and regulation mechanism of the reactor were mainly analyzed. According to the indexes of pH, VFAs/ alkalinity, volumetric gas production rate, membrane flux and effluent organic pollutants, the migration and transformation characteristics of COD during hydrolysis and acidification, methane production and membrane separation were also investigated. The results show that by coupling membrane scouring aeration and biogas aeration cycle together, VFAs/ alkalinity and volumetric gas production rates were optimized from 1.5 and 0.1 to 0.02 and 1.0, respectively. The performance of stably meeting the discharging standard for effluent COD and COD removal loading were also optimized, the COD removal rate and loading increased from 59% and 0.3 kg·(m³·d)⁻¹ to 97.7% and 1.8 kg·(m³·d)⁻¹, respectively, achieving the short-cut compliance of IC-AnMBR for the food waste leachate. Adopting biogas circulation to scour the ceramic membrane mitigated the accumulation of the filter cake layer and increased the membrane flux from 0.6 L·(m²·h)⁻¹ to 2.1 L·(m²·h)⁻¹. So IC-AnMBR short process can achieve the stable treatment of rural food waste leachate.
- food waste leachate /
- internal circulation anaerobic membrane bioreactor (IC-AnMBR) /
- COD removal rate /
- membrane flux

随着产业升级和城市扩张，大量工厂搬迁或废弃后遗留的场地存在土壤污染问题，需要进行土壤修复才能再次开发^[1]。电阻加热技术具有对环境扰动小，受土壤异质性影响小，处理深度大等优点，尤其适合修复含有挥发性、半挥发性有机污染物的污染场地^[2-5]。但在实际修复工程中，ERH技术的工程参数设计，例如电极间距、电场强度等，都会显著地影响场地电阻加热的实际效果^[6]，导致修复周期和成本控制的不确定性。修复场地的污染物分布、地下水流场和土壤特性等往往差异较大，但工程师只能根据已有工程经验和有限的取样勘探结果，进行原位加热工程的参数设计。若能够采用建模的方法，对场地条件下的加热过程进行预测，将有利于减少设计的盲目性，帮助缩短工程周期和控制修复成本。

目前，对于ERH技术的数值模型已经有了一定的探索和应用。HIEBERT等^[7-8]开发了用于模拟单相电阻加热过程的二维有限差分模型，并研究了不同的横卧电极设置方式对非均质含油地层的加热效果的影响。CARRIGAN等^[9]将改进的欧姆加热模型与非等温多孔流动和传输模型进行了耦合，研究了电极阵列的电相位如何影响电阻加热的均匀性。MCGEE等^[10]进一步简化了模拟多相电阻加热的欧姆方程，并模拟了电阻加热从非均质油砂中回收沥青的过程。KROL等^[11]考虑温度对密度、粘度、扩散系数的影响，建立了二维有限差分模型，模拟了电阻加热到50 ℃的情况下对地下水流动的影响，发现地下水流动方向和流速发生显著变化。许丹芸等^[12]使用有限元方法模拟了电阻加热土壤过程。

尽管关于ERH技术的数值模型研究已有一定的开展，但一方面，以往的模型对电阻加热土壤过程中的水分蒸发缺乏关注和进一步的验证；另一方面，大部分模型是针对实际场地的验证评估，对如何运用模型指导ERH工艺参数的选取探讨不足。本研究使用COMSOL多物理场耦合软件，基于有限元计算方法开展原位电阻加热温度场模拟研究，建立了考虑土壤水分蒸发的模拟原位电阻加热温度场的数值模型。通过对比土柱装置小试实验和数值模拟的结果，验证了数值模型的准确性，并利用数值模型分析了场地尺度下电场强度、电极间距和地下水流动对电阻加热温度场的影响。本研究结果有助于预测修复周期和优化电极井布设，从而达到节约能源和降低修复成本的目的。

1. 数值模型与实验方法

1.1 控制方程及条件假设

本研究利用COMSOL的传热模块、电流模块、PDE模块以及电流和传热耦合的电磁热模块，构建土壤电阻加热模型。为简化模型概念，在模拟时做如下假设：1)忽略土壤和水在加热过程中密度、热容随温度的变化；2)将土体视为均质且各向同性的多孔介质，土壤初始温度均匀一致；3)忽略土体在加热过程中的热变形；4)忽略电极和土壤之间的接触电阻；5)由于实验土柱较短，需要考虑水分在低于沸点时的挥发^[13]，但在模拟大尺度的场地加热中，忽略水分的挥发；6)忽略水蒸气对热量传递的影响。

电阻加热土壤过程中能量的控制方程如式(1)所示。

$stringUtils.convertMath($ !{formula.content})

$(1)$

式中：ρ_eff为单元的有效密度，kg·m⁻³；C_eff为单元的有效热容，J·(kg·K)⁻¹；T为温度，K；t为时间，s；λ为导热系数，W·(m·K)⁻¹；σ为土壤电导率，S·m⁻¹；E为电势梯度，V·m⁻¹；m_LG为水的气化速率，kg·(m³·s)⁻¹；ΔH_vap为水的潜热，J·kg⁻¹；ρ_f为流体密度，kg·m⁻³；C_f为流体热容，J·(kg·K)⁻¹；u_f为流体流速，m·s⁻¹。

土壤单元的有效密度ρ_eff和有效热容C_eff由土壤中固，液，气三相的体积分数决定，如式(2)~式(3)所示。

$stringUtils.convertMath($ !{formula.content})

$(2)$

$stringUtils.convertMath($ !{formula.content})

$(3)$

式中：θ表示各相的体积分数；ρ为各相的密度，kg·m⁻³；C为各相的热容，J·(kg·K)⁻¹；下标S，L，G表示固，液，气三相。

土壤含水量的控制方程如式(5)所示。

$stringUtils.convertMath($ !{formula.content})

$(4)$

$stringUtils.convertMath($ !{formula.content})

$(5)$

式中：D_L为导水系数，m²·s⁻¹；α为比例常数，m²·s⁻¹；θ_L^*为残余饱和度。

液态水变为气态水的情况可分为2种，一种是低于水的沸点时的挥发，一种是到达沸点时的沸腾，用式(6)可以得到水的气化速率。需要注意，只有水的饱和蒸气压(p*)大于等于外部气压(p_G)且含水量大于0时，沸腾才会发生。可以通过安托因方程^[14](式(7))计算不同温度下水的饱和蒸气压，进而判断温度是否到达沸点。

$stringUtils.convertMath($ !{formula.content})

$(6)$

$stringUtils.convertMath($ !{formula.content})

$(7)$

式中：A、B、C为经验常数；m_vap为挥发速率，kg·(m³·s)⁻¹；k_vap为蒸发速率常数，s⁻¹。

土壤单元的热导率会随着温度和含水量的变化发生极大的变化，在此使用TARNAWSKI等^[15]推导出的经验公式，如式(8)~式(9)所示。

$stringUtils.convertMath($ !{formula.content})

$(8)$

$stringUtils.convertMath($ !{formula.content})

$(9)$

式中:下标sat和dry分别表示饱和和干燥状态的土壤；a~g为经验常数；S_w为水饱和度。

电阻加热一般使用低频率电压(50~60 Hz)，产生的电磁波长远大于系统的物理尺寸，位移电流可以忽略，因此可以假设电阻加热产生的电场为准静态电场，可以通过将欧姆定律代入电流连续性方程中来求解电势分布，电流连续性方程如公式(10)所示。

$stringUtils.convertMath($ !{formula.content})

$(10)$

土壤电导率则使用Archie定律进行计算^[16-17]，考虑温度对电导率的影响^[18]，如式(11)所示。

$stringUtils.convertMath($ !{formula.content})

$(11)$

式中: ψ为电势，V；φ为孔隙度；m、n和β为经验常数(m为胶结系数，n为饱和度系数，β为温度系数)；σ_L为土壤溶液电导率，S·m⁻¹。

模拟单相交流电加热时，其电势分布与直流电基本一致，可以设置2个电极分别为接地和施加的电压。但模拟三相交流电时，电极电势与直流电则完全不同，三相交流电的电极电势可以视为由虚部和实部组成^[19]，如式(12)所示。

$stringUtils.convertMath($ !{formula.content})

$(12)$

式中：ω为角频率，rad·s⁻¹；t为时间，s；θ为相角，°(三相交流电分别为0 °、120 °和240 °)；j为虚部；E₀为正弦交流电振幅的绝对值，一般为电压的2^1/2倍，V。

1.2 几何模型和网格划分

电阻加热土柱装置模型如图1(a)所示，装置高16 cm，内径4.5 cm。装置的罐体用不导电的聚四氟乙烯制成，内部装填细砂。装置2边为对称分布的用螺丝与螺母固定的不锈钢片电极，螺丝作为导电连接，电线连接到螺丝上以引入电压，热电偶从装置底部插入，用于监测土柱中心点的温度。使用Solidworks对电阻加热土柱装置进行建模，得到电阻加热土柱装置的几何模型，再导入到COMSOL中用于模拟，几何模型的网格划分均采用三角单元网格，单元大小选择细化。

图 1 实验装置与场地模型

Figure 1. (a) Electrical resistance heating soil column device and (b) site geometric model

下载: 全尺寸图片幻灯片

为了模拟实际场地中电阻加热的过程，构建了图1(b)所示的三相电阻加热场地模型。场地模型为长宽15 m、高5 m的长方体，在场地中部按等边三角形放置3根半径20 cm、高5 m的电极，电极间距为6 m，并在3个电极构成的等边三角形中心点设置了温度监测点。

1.3 边界条件和参数设置

电阻加热土柱装置数值模拟的边界条件为：1)流动边界条件，所有边界均设置为0通量边界条件，这是因为已在方程中添加了描述水蒸发的汇项；2)温度边界条件，所有边界均设置为自然对流热通量边界；3)电势边界条件，2个电极分别设置为接地和电势，其他边界设置为电绝缘边界。土壤的初始温度、含水量、水的电导率和电极电压根据实测得到的初始值设置。对于各项参数的取值如表1所示。

表 1 数值模拟参数

Table 1. Numerical simulation parameters

模拟参数	取值	模拟参数	取值
液体密度，ρ_L	1 000 kg·m⁻³	温度系数，β	0.02
固体密度，ρ_S	2 650 kg·m⁻³	比例常数，α	5×10⁻⁶ m²·s⁻¹
气体密度，ρ_G	1.9 kg·m⁻³	经验常数，A	8.07131
液体热容，C_L	4 200 J·(kg·K)⁻¹	经验常数，B	1 730.63
固体热容，C_S	920 J·(kg·K)⁻¹	经验常数，C	233.426
气体热容，C_G	1 000 J·(kg·K)⁻¹¹	胶结系数，m	1.44
湿导热系数，λ_sat	1.87 W·(m·K)⁻¹	饱和度系数，n	2
干导热系数，λ_dry	0.23 W·(m·K)⁻¹	水的潜热，ΔH_vap	2 257.2 kJ·kg⁻¹
孔隙度，φ	0.5	蒸发速率系数，k_vap	1×10⁻⁶ s⁻¹

| Show Table

DownLoad: CSV

电阻加热场地模拟的边界条件为：1)流动边界条件，模拟地下水位上涨时，底部边界设置为通量边界条件，其他边界设置为0通量边界条件；2)温度边界条件，地下水流入的边界设置为流入边界，流入温度10 ℃，其他边界设置为热绝缘边界；3)电势边界条件，3个电极分别设置为三相电势中的1相，其他边界设置为电绝缘边界。各项参数的设置与土柱装置实验相同。

1.4 实验条件

为验证模型在不同的土壤含水量和土壤溶液电导率时的准确性，使用电阻加热土柱装置进行了验证实验。由于装置较小，电场强度较大，故使用电导率较小的细砂充当模拟土壤。分别向300 g细砂中加入40、50、60、70 mL的去离子水，质量分数为0.3%的NaCl溶液和质量分数为0.5%的NaCl溶液，充分搅拌，并装填进装置中。为了保持砂的孔隙度一致，在装填时保证装填高度为15 cm。加入去离子水的砂在100 V的条件下进行电阻加热，加入NaCl溶液的砂在50 V的条件下进行电阻加热。由于模拟只探究电阻加热的升温过程，所以实验只进行到温度到达100 ℃就停止，未加热到100 ℃时则在加热2 h后停止。

为模拟实际工程环境下土壤内热量传递规律，用实际场地模型探究了地下水流动、电场强度和电极间距的变化对土壤温度变化的影响。模拟的工况如表2所示，其他参数与表1一致，土壤溶液电导率设为0.1 dS·m⁻¹。

表 2 数值模拟工况

Table 2. Numerical simulation conditions

工况	电场强度/ (V·m⁻¹)	初始水饱和度	地下水流速/ (m·s¹)	电极间距/m	地下水涨速/ (m·d¹)
A-1	30	0.6	0	6	0
A-2	60	0.6	0	6	0
A-3	90	0.6	0	6	0
B-1	90/60	0.6	0	3	0
B-2	90/60	0.6	0	6	0
B-3	90/60	0.6	0	9	0
C-1	90	1	0.1	6	0
C-2	90	1	0.2	6	0
C-3	90	1	0.3	6	0
D-1	90	0.6	0	6	0.05
D-2	90	0.6	0	6	0.1
D-3	90	0.6	0	6	0.2

| Show Table

DownLoad: CSV

2. 结果与讨论

2.1 电阻加热土柱装置实验及模型验证

图2给出了不同含水率和电导率的情况下，土柱中心热电偶监测的温度变化。可以看出，随着含水量和土壤溶液电导率的增加，中心点加热到100 ℃所需要的时间不断减小。而当含水量较小或者土壤溶液电导率较小时，升温速率较小，甚至出现温度平台的情况。升温速率较小是由于此时土壤电导率较小^[20]，电流产生的焦耳热较小。升温速率减小则是因为，土柱中挥发掉的水份随着加热时间的增长逐渐变多^[21]，进一步降低了土壤电导率，使电流产生的焦耳热进一步减少。此外，由图2(a)可以看到，加入不同体积的去离子水后的细砂依然可以被加热，而去离子水的电导率几乎可以忽略不计。这说明，细砂中的离子溶解进入了去离子水中，提高了去离子水的电导率。

图 2 土柱加入不同体积不同溶液中心点温度实测值与模拟值

Figure 2. Measured and simulated values of the center point temperature of the soil column, added with different volumes of different solutions

下载: 全尺寸图片幻灯片

为验证模型的可靠性，将不同含水量和土壤溶液电导率的电阻加热实测值和模拟值进行了对比。利用均方误差MSE和平均相对误差MRE评价模拟值和实测值的差异(式(12)~式(13))。

$stringUtils.convertMath($ !{formula.content})

$(13)$

$stringUtils.convertMath($ !{formula.content})

$(14)$

式中：n为实测数据个数；M_i、S_i分别为第i个实测和模拟得到的数据。

图2(a)中给出的模拟值趋势线与实测值差异较小，实测值和模拟值均方误差为0.05~0.66，平均相对误差为0.42%~1.97%。较小的均方误差和平均相对误差表明，实测数据与模拟数据之间的偏差较小，模型具有较好的准确性。加入不同体积的0.3%NaCl(图2(b))和0.5%NaCl(图2(c))溶液，土壤溶液电导率随着NaCl溶液的体积和质量分数的增加而增大。对比实测和模拟预测结果，均方误差为0.38~12.29，平均相对误差为1.15%~5.32%。这表明，在较宽的土壤电导率范围内(加入70 mL的0.5%NaCl溶液时，土壤的电导率为1.1 dS·m⁻¹)，实测数据与模拟数据之间的偏差也较小，模型用于预测不同电导率的土壤加热过程是可靠的。误差产生的原因可能是，对土壤原有的电导率，实验过程中空气热对流导致的热量散失以及水分挥发速率的估值存在一定的偏差。

2.2 场地条件下电场强度和电极间距对电阻加热的影响

电场强度是影响电阻加热升温速率的重要因素，由电极电压和电极间距共同决定，为了考察场地尺度下的温度场以及各条件对温度场的影响，建立了图1(b)的场地模型。图3和图4分别展示了电极间距6 m时，30、60、90 V·m⁻¹电场强度下加热70 d后，温度场的横截面以及位于3个电极中心点(见图1(b))温度的变化。从图3可以看出，电场强度越大，土壤升温速率越快，电极附近的升温更快，温度更高。这与MCGEE等^[10]和HAN等^[22]的研究结果是一致的。电流产生的焦耳热随电场强度的增大而增大，电场强度越大，土壤升温速率越快；电极附近的电流密度最高，产生的焦耳热最多，所以土壤升温速率更快。

图 3 不同电场强度下加热70 d模拟场地横截面温度场

Figure 3. Cross-sectional temperature field of the site heated for 70 days under different electric field intensities

下载: 全尺寸图片幻灯片

图 4 电场强度对模拟场地中心点温度变化的影响

Figure 4. Effect of electric field intensity on central point temperature of simulated site

下载: 全尺寸图片幻灯片

从图4可以看到，场地在加热90 d后，30、60和90 V·m⁻¹电场强度获得的中心温度分别为25.5、49.1和100.0 ℃，中心点的平均升温速率分别为0.079，0.420和1.210 ℃·d⁻¹。由式(1)可知，焦耳热与电场强度的平方呈正比，电场强度为30、60、90 V·m⁻¹时中心点升温速率之比应为1∶4∶9，小于模拟得到的比值，即1∶5.3∶15.5。这说明，中心点的温度可能是电流焦耳热和外部热传导叠加共同决定。此外，电场强度为90 V·m⁻¹时，当中心点温度达到97 ℃后，升温速率明显放缓。这是因为，此时的电极温度已经达到水的沸点，电极土壤水分蒸发带走了大量热量，导致土壤热导率和土壤电导率下降，从而使升温速率下降。

电极间距的设置决定了电极井的数量，会极大地影响修复场地的成本和热传导过程。为探究电极间距这一单因素对电阻加热过程的影响，在保持电场强度90 V·m⁻¹不变的条件下，模拟电极间距3、6和9 m时中心温度的变化，结果如图5所示。从图5(a)可以看出,在电场强度为90 V·m⁻¹，电极间距为3、6和9 m时，中心点的温度达到100 ℃的时间分别为84、62和75 d，达到100 ℃的能耗分别为10 418、23 375和51 311 kW·h。综上可知，6 m为最适宜的电极间距，此时中心点的升温速率最快，且相较于电极间距为3 m时场地修复需要布设的电极井数量更少，相较于电极间距为9 m时需要的能耗更少。从图5(a)还可以看出，随着电极间距的增大，中心点升温速率先增大后减小。这是因为，为了保证电场强度一定，电压随电极间距增大而增大，导致电极处的升温速率也随之增大，这一定程度上弥补了间距增大导致的中心点热传导距离增大的不足。但是，9 m间距下的电压增加，造成电极附近温度过早达到水的沸点，导致电极周围土壤水分过早蒸发，土壤热导率和电导率下降，升温速率略为下降。

图 5 90 V·m⁻¹电场强度下电极间距变化和电极补水对模拟场地中心点温度的影响

Figure 5. Effect of electrode spacing and replenishment of water near the electrode on central point temperature of simulated site under the electric field intensity of 90 V·m⁻¹

下载: 全尺寸图片幻灯片

为了考察电极附近土壤水分蒸发对中心点升温速率的影响，模拟了电场强度90 V·m⁻¹、电极间距为9 m时，对电极附近补水的情况下，中心点温度的变化。从图5(b)可以看到，在加热40 d后，补水的情况下中心点温度明显高于不补水的情况，补水的情况下中心点的温度达到100 ℃的时间为68 d，比不补水的情况下早8 d。此结果表明，电极周围土壤水分过早蒸发是造成升温速率下降的原因之一，补水可以使土壤含水率增大，增加土壤的电导率和热导率，提高了升温速率。这与葛松等^[3]的研究结果一致。监测电极的电流变化，并及时对电极附近进行补水，对于更快加热到目标温度至关重要。

2.3 场地条件下地下水流动对电阻加热的影响

使用原位电阻加热的场地一般修复深度都较大，有必要考虑地下水流动和地下水位变化对温度场的影响。为进行预测，分别设定了0、0.1、0.2、0.3 m·d⁻¹的地下水流速下加热30 d和0、0.05、0.10、0.20 m·d⁻¹的地下水上涨速度下加热30 d的条件进行模拟，流入的地下水出的温度假定为10 ℃，温度场的变化结果如图6和图7所示。从图6可以看到，随着地下水流速的增加，场地左侧的温度逐渐降低到10 ℃，温度场右移的趋势越明显，中心点的最高温度从100.0 ℃降低到38.5 ℃。这说明，较高的地下水流速会将场地中的热量带向下游，对场地升温速率造成不利影响，与MUNHOLLAND等^[23]在二维沙箱中得到的结果一致。实际工程中应尽量减小地下水流速，必要时可增加侧向的物理屏障来形成止水帷幕，或通过蒸汽注入等手段来增加流入地下水的温度。

图 6 地下水流速对场地横截面温度场的影响

Figure 6. Influence of groundwater velocity on temperature field of site

下载: 全尺寸图片幻灯片

图 7 地下水上涨速度对场地纵截面温度场的影响

Figure 7. Influence of groundwater rising speed on temperature field of site

下载: 全尺寸图片幻灯片

从图7可以看到，随着地下水位上涨速度提升至0.1 m·d⁻¹，场地下部和两侧的温度逐渐降低到10 ℃，但由不饱和区域变为饱和区域的土壤升温速率上升，中心点最高温度从51 ℃升高到100 ℃。随着地下水位上升速度达到0.2 m·d⁻¹，相较于0.1 m·d⁻¹时，温度场明显向上移动，中心点最高温度从100 ℃降低到54 ℃。这是因为，地下水位上涨可以起到一定的补水作用，提高不饱和区域土壤的电导率，从而提高土壤升温速率，但地下水上涨速度过快还是会导致流失的热量多于通过补水增加的焦耳热，使得场地温度从底部开始下降，不利于场地的修复。应当指出，由于实际场地地层的复杂性，地下水的流速和流向并不是均匀的。因此，对于非均质地层，可以建立相应的多模块的空间耦合模型，通过改变各地层模块的物理参数，从而实现数值模拟目标。

3. 结论

1)电阻加热温度场数值模型与土柱装置实验数据吻合度较好，模拟值和实测值均方误差为0.05~12.29，平均相对误差为0.42%~5.32%，数值模型具有较好的预测准确性。

2)场地模型研究发现，升温速率与电场强度成正比，电极处升温速率最快；电场强度为90 V·m⁻¹时6 m为最适宜的电极间距，此时中心点升温速率最快，能耗相对较低，需要的电极井较少；模型考虑水分蒸发的情况下，电极周围土壤水分会更早蒸发，导致土壤升温速率下降，此时对电极附近补水可以显著增加土壤升温速率，故在实际工程中应监测电流和功率，并及时对电极附近进行补水。

3)场地地下水流动会带走热量，不利于场地的热修复。但对于不饱和场地，地下水位上涨速度小于等于0.1 m·d⁻¹时可以起到补水作用，提高场地升温速率。故对于地下水流速过高的场地应采取水力平衡或设置止水帷幕等工程措施以缓解热量流失。

图 1 反应器结构

Figure 1. Structure of reactor

下载: 全尺寸图片幻灯片

图 2 运行时间对COD去除率与容积负荷的影响

Figure 2. Effect of operation time on COD removal rate and volumetric loading rate

下载: 全尺寸图片幻灯片

图 3 反应器运行期间的VFAs/碱度和容积产气率的变化

Figure 3. Change of VFAs/ alkalinity and biogas yield during operation time

下载: 全尺寸图片幻灯片

图 4 反应器运行期间酸碱缓冲体系的变化

Figure 4. Change of acid-base buffer system during operation time

下载: 全尺寸图片幻灯片

图 5 反应器运行期间膜通量与跨膜压差的变化

Figure 5. Change of TMP and flux during operation period

下载: 全尺寸图片幻灯片

图 6 厌氧污泥粒径分布

Figure 6. Size distribution of anaerobic sludge

下载: 全尺寸图片幻灯片

图 7 膜滤饼层与甲烷八叠球菌电镜照片

Figure 7. SEM images of filter cake layer and Methanosarcina thermophile

下载: 全尺寸图片幻灯片

表 1 废水水质

Table 1. Water quality of wastewater

阶段	运行时间/d	pH	COD/(mg·L⁻¹)	NH₃-N/(mg·L⁻¹)	HRT/d
Ⅰ	1~35	3.56±0.17	24 322±9 845	77.1±36.5	30
Ⅱ	35~108	3.53±0.77	18 998±3 479	47.9±35.8	45
Ⅲ	108~210	3.47±0.57	34 512±5 602	113.26±73.5	20

下载: 导出CSV

表 2 出水水质

Table 2. Water quality of effluent

阶段	运行时间/d	HRT/d	容积负荷/(kg·(m³·d)⁻¹)	COD去除率/%	COD/(mg·L⁻¹)	pH	TIC/(mg·L⁻¹)	NH₃−N/(mg·L⁻¹)	容积产气率/(L·(L·d)⁻¹)
Ⅰ	1~35	30	0.61±0.40	59.1±26.4	7 440±2 121	8.23±0.39	2 555.8±522.0	2 936.5±448.0	—
Ⅱ	35~108	45	0.30±0.14	84.1±9.5	3 064±2 168	8.17±0.14	1 456.1±358.6	1 531.0±451.6	0.111±0.100
Ⅲ	108~210	20	1.17±0.51	97.8±0.70	675±283	7.71±0.37	705.1±78.9	885.2±120.8	0.350±0.266

下载: 导出CSV

[1]	国家统计局. 国家数据[EB/OL]. [2020-01-05]. http://data.stats.gov.cn/.
[2]	李雪兵. 农业废弃物与厨余垃圾共发酵工艺的优化研究[D]. 西安: 西安建筑科技大学, 2018.
[3]	刘银秀, 董越勇, 边武英, 等. 浙江省农村沼气利用典型技术的表征和演进[J]. 浙江农业科学, 2019, 60(12): 2295-2299.
[4]	WANG K, YIN J, SHEN D S, et al. Anaerobic digestion of food waste for volatile fatty acids (VFAs) production with different types of inoculum: Effect of pH[J]. Bioresource Technology, 2014, 161: 395-401. doi: 10.1016/j.biortech.2014.03.088
[5]	BRAGUGLIA C M, GALLIPOLI A, GIANICO A, et al. Anaerobic bioconversion of food waste into energy: A critical review[J]. Bioresource Technology, 2018, 248(A): 37-56.
[6]	LI Y, CHEN Y G, WU J. Enhancement of methane production in anaerobic digestion process: A review[J]. Applied Energy, 2019, 240: 120-137. doi: 10.1016/j.apenergy.2019.01.243
[7]	KONDUSAMY D, KALAMDHAD A S. Pre-treatment and anaerobic digestion of food waste for high rate methane production: A review[J]. Journal of Environmental Chemical Engineering, 2014, 2(3): 1821-1830. doi: 10.1016/j.jece.2014.07.024
[8]	ULLAH K I, HAFIZ D O M, HASHIM H, et al. Biogas as a renewable energy fuel: A review of biogas upgrading, utilisation and storage[J]. Energy Conversion and Management, 2017, 150: 277-294. doi: 10.1016/j.enconman.2017.08.035
[9]	CHENG H, HIRO Y, HOJO T, et al. Upgrading methane fermentation of food waste by using a hollow fiber type anaerobic membrane bioreactor[J]. Bioresource Technology, 2018, 267: 386-394. doi: 10.1016/j.biortech.2018.07.045
[10]	迟文涛, 赵雪娜, 江翰, 等. 厌氧反应器的发展历程与应用现状[J]. 城市管理与科技, 2004, 6(1): 31-33. doi: 10.3969/j.issn.1008-2271.2004.01.011
[11]	国家环境保护总局. 水和废水监测分析方法[M]. 4版. 北京: 中国环境科学出版社, 2002.
[12]	GRIMBERG S J, HILDERBRANDT D, KINNUNEN M, et al. Anaerobic digestion of food waste through the operation of a mesophilic two-phase pilot scale digester-assessment of variable loadings on system performance[J]. Bioresource Technology, 2015, 178: 226-239. doi: 10.1016/j.biortech.2014.09.001
[13]	LUO G, LI J, LI Y, et al. Performance, kinetics behaviors and microbial community of internal circulation anaerobic reactor treating wastewater with high organic loading rate: Role of external hydraulic circulation[J]. Bioresource Technology, 2016, 222: 470-477. doi: 10.1016/j.biortech.2016.10.023
[14]	WANG Y Y, ZHANG Y L, WANG J B, et al. Effects of volatile fatty acid concentrations on methane yield and methanogenic bacteria[J]. Biomass and Bioenergy, 2009, 33(5): 848-853. doi: 10.1016/j.biombioe.2009.01.007
[15]	LI L, HE Q M, WEI Y M, et al. Early warning indicators for monitoring the process failure of anaerobic digestion system of food waste[J]. Bioresource Technology, 2014, 171: 491-504. doi: 10.1016/j.biortech.2014.08.089
[16]	WANG X M, BAI X, LI Z F, et al. Evaluation of artificial neural network models for online monitoring of alkalinity in anaerobic co-digestion system[J]. Biochemical Engineering Journal, 2018, 140: 85-92. doi: 10.1016/j.bej.2018.09.010
[17]	LANE N, MARTIN W F. The origin of membrane bioenergetics[J]. Cell, 2012, 151(7): 1406-1416. doi: 10.1016/j.cell.2012.11.050
[18]	YU D W, MENG X S, LIU J B, et al. Formation and characteristics of a ternary pH buffer system for in-situ biogas upgrading in two-phase anaerobic membrane bioreactor treating starch wastewater[J]. Bioresource Technology, 2018, 269: 57-66. doi: 10.1016/j.biortech.2018.08.072
[19]	YENIGVN O, DEMIREL B. Ammonia inhibition in anaerobic digestion: A review[J]. Process Biochemistry, 2013, 48(5/6): 901-911.
[20]	LEE D Y, XU K Q, KOBAYASHI T, et al. Effect of organic loading rate on continuous hydrogen production from food waste in submerged anaerobic membrane bioreactor[J]. International Journal of Hydrogen Energy, 2014, 39(30): 16863-16871. doi: 10.1016/j.ijhydene.2014.08.022
[21]	MONTES J A, LEIVAS R, MARTINEZ-PRIETO D, et al. Biogas production from the liquid waste of distilled gin production: Optimization of UASB reactor performance with increasing organic loading rate for co-digestion with swine wastewater[J]. Bioresource Technology, 2019, 274: 43-47. doi: 10.1016/j.biortech.2018.11.060
[22]	苑宏英, 李琦, 杨玉萍, 等. pH对蛋白类餐厨垃圾发酵产酸的影响[J]. 环境工程学报, 2018, 12(10): 2929-2934. doi: 10.12030/j.cjee.201804102
[23]	张玉静, 蒋建国, 王佳明. pH值对餐厨垃圾厌氧发酵产挥发性脂肪酸的影响[J]. 中国环境科学, 2013, 33(4): 680-684. doi: 10.3969/j.issn.1000-6923.2013.04.015
[24]	XIAO K, LIANG S, WANG X M, et al. Current state and challenges of full-scale membrane bioreactor applications: A critical review[J]. Bioresource Technology, 2019, 271: 473-481. doi: 10.1016/j.biortech.2018.09.061
[25]	JEONG Y, HERMANOWICZ S W, PARK C. Treatment of food waste recycling wastewater using anaerobic ceramic membrane bioreactor for biogas production in mainstream treatment process of domestic wastewater[J]. Water Research, 2017, 123: 86-95. doi: 10.1016/j.watres.2017.06.049

点击查看大图

图( 7) 表( 2)

计量

文章访问数: 5586
HTML全文浏览数: 5586
PDF下载数: 96
施引文献: 0

1. 数值模型与实验方法
1.1 控制方程及条件假设
1.2 几何模型和网格划分
1.3 边界条件和参数设置
1.4 实验条件
2. 结果与讨论
2.1 电阻加热土柱装置实验及模型验证
2.2 场地条件下电场强度和电极间距对电阻加热的影响
2.3 场地条件下地下水流动对电阻加热的影响
3. 结论

全文HTML

村镇承载着我国约40%的人口和农副食品的初级加工^[1]，厨余垃圾对村镇环境影响较大^[2]。由于厨余垃圾易腐败酸化，因此，必须进行及时处理，才能减少厨余垃圾腐败酸化带来的环境影响，这同时也限制了村收集-镇转运-县处理模式的正常运转。基于此，依托农村现有沼气工程^[3]，就地厌氧处理利用，探索可行有效的技术路线，具有重要现实意义。

厨余渗滤液中易降解COD含量较高，极易导致厨余废水的酸化，从而对厌氧工艺的稳定运行带来较大风险^[4]，并会进一步影响厌氧微生物的产甲烷和产乙酸过程，导致甲烷产量降低乃至反应器酸化崩溃。已有研究^[5]表明，加碱调pH、重新接种和混合厌氧发酵是常用的控制酸化方法，这些方法在村镇应用时存在药剂投加量大、运行维护复杂等缺点，造成运行成本上升15%~37%^[6-8]，因此，制约了上述方法在村镇处理厨余垃圾渗滤液过程中的推广应用。在酸化崩溃的早期，厌氧体系内CO₂大量溢出，不仅使氢营养型产甲烷菌底物浓度下降，限制氢营养型产甲烷途径，还会导致系统中碱度出现损失^[9]，影响系统运行稳定性。因此，通过沼气循环曝气，促进沼气中CO₂在厌氧体系中的再溶解，强化CO₂气液传质，提升厌氧系统中CO₂浓度，从而促进氢营养型产甲烷反应，可在酸化早期控制过度酸化现象。内循环(IC)厌氧反应器具有二次气液循环系统，通过循环管路增强反应器内气液混合，无需额外加压沼气曝气^[10]，有利于实现村镇厨余渗滤液废水的就地处理和利用；通过该系统与外置膜系统组合，克服IC反应器出水水质不稳定缺点，有望一步满足工程排放需求，扩大厌氧工艺在处理类似高浓度废水中的应用范围。

本研究针对我国村镇厨余渗滤液处理需求，构建了IC-AnMBR内循环厌氧膜生物反应器，通过反应器构型的设计和过程控制的优化，缓解由VFAs累积造成的酸化抑制，以短流程工艺实现厨余渗滤液的稳定处理，并考察沼气循环曝气对膜污染的缓解效果，以期为AnMBR工艺的高效运行提供参考。

3. 结论

1)本研究所构建的IC-AnMBR经合理调控，能够实现对村镇厨余等高浓度有机废水的有效处理，COD去除率、出水COD分别达到97.7%、(675±283) mg·L⁻¹。

2)通过氨氮与TIC对废水酸化中VFAs的缓冲作用，可减少进水冲击引起的pH剧烈变化，缓解反应器酸化崩溃风险，能够提升反应器负荷，实现反应器稳定运行与高效产气；反应器容积负荷最高能够达到2.1 kg·(m³·d)⁻¹、最佳容积产气率能够到达1.014 L·(L·d)⁻¹。

3)沼气曝气能够强化CO₂气液传质，增加系统缓冲能力；抑制膜表面滤饼层的形成，能够使膜通量稳定维持在2.19 L·(m²·h)⁻¹，实现对反应器的高负荷运行及对COD的高效去除。

参考文献 (25)

IC-AnMBR处理村镇厨余渗滤液废水稳定运行调控机制

作者简介: 王拓(1995—)，男，硕士研究生。研究方向：膜生物反应器。E-mail：wangtuo17@mails.ucas.edu.cn

通讯作者: 郁达伟(1982—)，男，博士，助理研究员。研究方向：膜生物反应器等。E-mail：dwyu@rcees.ac.cn;