本研究对北京某生活垃圾干式厌氧发酵工程进行了为期1 a的研究。该工程采用法国Valorga仓筒型干式厌氧发酵技术，该技术的特点为：1)适用于高TS含量(可高达60%)和高容积有机负荷(OLR: organic loading rate)(最高可达15 kg·(m³·d)⁻¹，以VS计)的厌氧发酵，可采用间歇或连续的进料方式；2)发酵温度通常采用中温(35 ℃)或高温(55~58 ℃)；3)通过气动搅拌即沼气循环方式实现厌氧发酵物料的搅拌、混合^[1]。

该生活垃圾干式厌氧发酵工程主要包含预处理、干式厌氧发酵、沼气处理与利用、发酵残渣处理4个工艺段。预处理段将组成复杂的生活垃圾进行分类、挤压破碎、过筛后，与厌氧发酵回流物料混合(回流量为进料量体积的3~5倍)，泵入厌氧发酵罐。核心工艺采用干式厌氧发酵，预处理后的浆料泵入厌氧罐进行发酵，产生的沼气经加压注入发酵罐，对物料进行气动搅拌，一方面加快沼气的逸出和维持发酵罐内的温度均匀；另一方面推动物料往出料口方向缓慢流动。沼气处理与利用工艺段对生成的沼气进行过滤、冷凝脱水等处理，所得的沼气一部分回注厌氧罐用于搅拌，其余的经后续脱硫处理后，可进入热电联产系统。发酵残渣经过压榨脱水处理后可进行堆肥或外运处理。

所研究规模化干式厌氧发酵工程处理对象为生活垃圾，原料来自于北京市某垃圾转运站的生活垃圾。在本研究开展期间，北京尚未严格实施垃圾分类，运入厂区的生活垃圾有机组分偏低。经预处理后，进入厌氧罐发酵的物料中易生化有机组分平均含量约为62%，较难生化有机组分平均含量约为27%，而惰性物质包括塑料、砂砾玻璃、金属物质、木竹类和骨头贝类等平均含量约为11%。原料中大块惰性物质(如织物、陶瓷制品、树枝等)经过预处理后作为废料运出，但仍有部分直径较小的惰性物质(碎骨头、玻璃、砂砾等)进入厌氧发酵系统，沉积在厌氧罐内占据空间，并可能造成进、出料口及管道的堵塞，降低厌氧发酵效率。

工程采用中温干式厌氧发酵工艺，启动了2个并联干式厌氧罐。单个厌氧罐(直径16.5 m，高21.8 m)设计有效容积为3 735 m³，单罐运行有效容积约2 000 m³，2个厌氧罐共用1套预处理、进出料设备、搅拌系统和沼气处理系统，总产气量为2个厌氧罐的产气量之和。

正常工况下每天进料2次，进料时间依据进料量在1.5~2.5 h，进料的同时控制出料量，并通过控制进、出料的体积调节罐内液位高度。在厌氧发酵过程中，当pH≤7.5、挥发性脂肪酸(VFAs)≥8 000 mg·L⁻¹等生物不平衡指征出现时减少进料量或停止进料，直至系统pH、平均产气量等指标恢复。工程运行时在环境低温期通过蒸汽加热进料提高厌氧发酵运行温度，其他季节不进行加热操作。

在剔除大块木料、玻璃、金属等物质后，有机生活垃圾及厌氧发酵物的TS、VS采用重量法测定；总挥发性脂肪酸(total volatile fatty acids, TVFAs)采用氢氧化钠滴定法测定，氨氮采用纳氏试剂分光光度法测定；碱度采用酸碱指示剂滴定法^[18]。沼气产量采用E+H超声波流量计(GFM 200)测定。

本研究采用候平均气温划分季节，当候平均气温小于10 ℃为冬季，10~22 ℃为春季或秋季(从10 ℃升到22 ℃是春季，从22 ℃降到10 ℃是秋季)，大于22 ℃划为夏季^[19]。因此，基于北京2019年全年气温数据，本研究中的四季划分为：春季3月14日至5月20日、夏季5月21日至10月2日、秋季10月3日至11月12日、冬季11月13日至3月13日。本研究所述的温差是不同季节下该生活垃圾干式厌氧发酵运行温度与气温的温差，并重点分析温差与厌氧发酵产气效能(VS产气率、容积产气率等)、操作参数(垃圾处理量、OLR、pH等)的相互关系。VS产气率、容积产气率、容积有机负荷的计算公式如式(1)~式(3)所示。

式中：$ {{p}}_{\rm{m}} $为当日VS产气率，m³·(kg·d)⁻¹；P是日产沼气量，m³·d⁻¹；m为进料中VS质量，kg；$ {p}_{\rm{v}} $为容积产气率，m³·( m³·d)⁻¹；V是厌氧有效容积，m³。

不同季节下该生活垃圾干式厌氧发酵工程的主要运行参数如表1所示。该工程在冬、夏2季的运行时间均较长，分别为121和135 d，且冬季的进料量、进料VS/TS比值最高，但冬季的进料TS含量最低。该工程全年厌氧发酵运行温度在24~30 ℃，四季环境气温与厌氧发酵温度的平均温差从高到低依次为冬季(26.63 ℃)>春季(12.45 ℃)>秋季(11.92 ℃)>夏季(3.43 ℃)。

由图1(a)可知，不同于严格受控的实验室小试研究，该干式厌氧发酵工程采用连续运行方式，其进料量通常根据所需处理的垃圾量和运行参数调控的需求进行调整，其日均进料量为(40.38 ± 26) t·d⁻¹，波动较大(变异系数0.65)，进料TS在35%～55%，出料TS约为15%。月进料量(即处理量)波动较大，不同季节进料量差异较显著，冬季垃圾处理量最多，占全年处理量的40.88%(图1(b)、图1(c))。干式厌氧发酵工程实际运行期间，处理未进行垃圾分类的原料时，包含预处理废料和发酵残渣的垃圾残渣日外运量基本上高于每日进罐物料量，即日处理量。因此，为提高干式厌氧发酵工程的原料品质和产气效能，亟需做好生活垃圾源头分类。

干式厌氧发酵工程受进料量、物料特性、运行温度等参数综合影响，不同季节中日产气量波动幅度变化较大，如图2(a)所示，前90天即1-3月波动幅度为全年最小，而年日均产气量为(9 243±1 248) m³·d⁻¹，日变异系数为0.14，最大日产气量14 498 m³·d⁻¹出现在夏季，是最少日产气量(冬季，6 728 m³·d⁻¹)的2.15倍。月均日产气量在7 989 m³·d⁻¹ 与10 623 m³·d⁻¹之间波动(图2(b))。其中，4月的变异系数最大(0.16)，最高月均日产气量出现在6月，最低月均日产气量是1月。按季节分析，如图2(c)，产气量总体按春夏秋冬顺序呈先增加后减少再逐步上升的循环特征，产气量均值由高到低依次为：夏季>春季>冬季>秋季。其中，夏季和冬季的产气量波动较大可能受到进料量波动的直接影响，对应变异系数分别为0.14和0.13。该干式厌氧发酵工程全年累计产沼气量约340×10⁴ m³，沼气中甲烷平均含量约为65%。经拟合分析，累积产气量与时间符合线性函数模型，数学模型见式(4)。

式中：y_沼气为全年累积产沼气量，m³；t为时间，d。R²为0.999 7。

根据候平均气温划分，北京季节特征呈现为夏、冬季时间长，春、秋季短(表1)。调研年度内夏、冬季累积产气量占全年产气量的71.33%。其中，夏季贡献39.20%，冬季贡献32.13%。厌氧发酵温度是影响产气的关键参数，中温干式厌氧发酵温度易受到季节性环境气温变化的影响，尤其北京冬季低温条件不利于厌氧发酵反应。因此，针对北京的季节特征，为保障该干式厌氧发酵工程的稳定高效运行，需做好厌氧发酵温度的调控，以实现冬季低温期、春秋季环境气温升高和降低过程的平稳产气。

如图3(a)所示，平均容积产气率与VS产气率的季节性变化规律均为先升高后降低，但最高平均容积产气率出现在夏季((2.44±0.33) m³·(m³·d)⁻¹)，而最高VS产气率则出现在秋季，秋季的进料量最少，因此有可能在牺牲容积产气率的情况下进行了更充分的厌氧发酵产气过程。进料垃圾成分具有显著异质性是导致不同季节OLR均有较大波动(变异系数在0.58~0.89)的主要因素，如图3(b)，冬季容积有机负荷(OLR)最高((3.42±1.98) kg·(m³·d)⁻¹)，而后依次是春季、夏季和秋季。不同批次的进料VS含量不同，这更凸显了垃圾分类的重要性。一方面，OLR的提高增加了新鲜物料对厌氧微生物的供给，在一定范围内可增加发酵产气；另一方面，随着OLR的提高，大量物料的分解、水解引起挥发性脂肪酸(VFAs)、氨氮等可能在系统中累积，进而抑制微生物的活性，打破厌氧发酵系统的微生物平衡，导致产气量下降。因此，冬季容积产气率较低的原因可能是OLR升高导致产气受到抑制。

如表2所示，该工程全年平均VS产气率为(1.28±0.66) m³·kg⁻¹，其中秋季平均VS产气率最高，是冬季均值的2.16倍，高于ROCAMORAL^[1]所报道的生活垃圾干式厌氧发酵工程。此外，提高OLR增加了产气量，但同时也降低了单位有机质产气率，表明VS产气率与OLR呈极显著负相关关系(p<0.001)。

1)运行温度。如图4(a)、图4(b)所示，自1月起，气温从冬季开始逐步升温，到夏季最高，再逐步降温，顺序经历春夏秋季，再到冬季，形成向下抛物线形的气温变化规律。不同于严格受控的中温和高温厌氧发酵需要持续外加热源，该中温干式厌氧发酵工程只在环境低温时(冬季和春季部分时间)进行外加热以维持厌氧微生物的活性，故其运行温度易受环境温度波动的影响。结合图4(a)和表1的温度参数可知，该工程全年运行温度变化范围为21~34 ℃。其中，夏季运行温度最高(图4(c))，在30 ℃运行维持了71 d，秋季运行温度最低，而冬季、春季初始阶段(图4(a)中1月至4月，11月中旬至12月底)通过蒸汽加热维持了一定的厌氧发酵运行温度，调控后运行温度维持在25~30 ℃，实际运行中温差大小依次是：冬季(26.63 ℃)>春季(12.45 ℃)>秋季(11.92 ℃)>夏季(3.43 ℃)。

工程上中温干式厌氧发酵具有规模效应，且易受环境温度的影响，季节性热量耗散率差异大，调节或维持厌氧运行温度难度大。因此，结合运行温度变化，该工程并不严格表现为随温度升高VS产气率提高。运行温度最低的秋季表现出最高的VS产气率，其原因可能是秋季的垃圾处理量较少，容积有机负荷低，处理单位垃圾的微生物菌群数量更多，使物料分解地更完全。此外，夏季较佳的厌氧发酵状态也持续影响到秋季，且生物产气有滞后效应，其结果表现为秋季VS产气率提高。秋季时间最短，又属于厌氧运行温度变化的过渡期，秋季的温度调控对干式厌氧发酵的平稳运行有重要影响。

在冬季采用外加蒸汽的操作，干式厌氧发酵运行温度均在25 ℃以上，甚至在11月和12月期间大多维持在30 ℃以上。然而，冬季的VS产气率和容积产气率均低于其他季节。其可能原因是，由垃圾原料的初始温度低引起的厌氧罐内外物料温差大，影响发酵效率。城市生活垃圾的收运通常在全天气温最低点的凌晨进行，这导致冬季进料时垃圾原料温度在0 ℃左右甚至更低，虽通过蒸汽加热将物料表面温度提升至39 ℃，但由于进料量大、有机负荷高、温差大等因素造成进料与罐内物料的传热与传质不均，加之微生物对温度变化敏感，故导致有机物未能及时被微生物分解利用，产气量低。考虑到进料量占整个厌氧罐物料的比例低，只加热进料的点源加热方式不能有效提升厌氧发酵产气效能，可考虑优化加热模式和气动搅拌频率、增加回流比、减缓进料速度、降低有机负荷等方式，提升进料与罐内物料的传热传质均匀性，从而提高冬季的VS产气率和容积产气率。

2)主要参数变化特征。VFAs、碱度和pH值是干式厌氧发酵主要参数，其中TVFAs和碱度的季节性变化如图5(a)所示。按季节分析，其变化总体表现为TVFAs增加，碱度降低。夏季TVFAs较高，相应的碱度((15 010±2 520) mg·L⁻¹)均低于春、秋、冬季，而春季具有较低的TVFAs浓度和较高的碱度。全年干式厌氧发酵体系的TVFAs在2 500~10 900 mg·L⁻¹变化，属于高 TVFAs浓度的厌氧体系。挥发性脂肪酸(VFAs)是厌氧发酵过程重要的中间产物，作为产甲烷菌的底物被利用生成甲烷。由于干式厌氧发酵的含固率高，物料搅拌难度大，因混合不均使代谢生成的VFAs不能很好地通过搅拌向周围扩散进而被产甲烷菌及时利用，造成VFAs在局部累积，局部酸浓度增高抑制产甲烷菌活性从而减少了产气。夏季环境温度和运行温度均高于其他3个季节，温度升高加速了有机质的分解和产酸速率，可能形成了局部抑制，是夏季VS产气率低于秋季的原因之一。VFAs代谢平衡对干式厌氧发酵体系稳定运行有重要作用。因此，在工程运行过程中应通过监测、预判、调控及时发现并缓解发酵体系酸化情况。此外，碱度也是评价厌氧发酵体系稳定运行的关键参数。

厌氧发酵体系的酸碱度受复杂微生物过程和化学过程控制，如有机酸、无机酸、钙等参与酸碱平衡反应，而pH是酸碱平衡综合呈现的指标^[20]。由图5(b)可知，全年运行过程中pH值在7.5~8.0波动。其中，秋季的平均pH值(7.94±0.13)最高。这表明，各阶段都处于偏碱性环境，有利于避免酸化，维持发酵系统的稳定性。产甲烷菌对pH值的变化比产酸菌更为敏感，当pH值下降引起产甲烷菌活性降低，VFAs因代谢缓慢而累积，导致pH值进一步降低。因此，VFAs和pH值的变化是综合判断厌氧发酵系统是否稳定运行的重要指征。如系统出现pH值连续下降的趋势，TVFAs有所升高，同时，碱度降低，表明出现酸化的前兆，此时应及时调整进料量，降低有机负荷，以调控厌氧系统恢复到新的平衡。高含固或干式厌氧发酵系统中的氨氮浓度高，由于氨氮呈弱碱性，不少研究者日益认识到高浓度氨氮也贡献了碱度，在厌氧发酵pH调控中需考虑氨氮-碳酸盐-VFA的三元酸碱缓冲体系，讨论厌氧发酵系统的pH缓冲能力及相应的pH变化规律^[21-23]，这可为干式厌氧发酵的pH调控提供新思路。

3)pH缓冲能力及系统稳定性。厌氧发酵系统的pH缓冲能力一般是指系统对酸碱的抵抗能力，是厌氧发酵过程稳定性及生物不平衡自我调节能力的重要指征，可以指示或预警厌氧发酵系统的抑制状况。因此，厌氧发酵过程常用总挥发性脂肪酸(TVFAs)和碳酸氢盐碱度的比值作为厌氧系统pH缓冲能力大小的指标之一，其指示厌氧发酵系统的稳定性和平衡性，比值越小，系统缓冲能力越大。不同季节下TVFAs/碱度的比值变化如图6所示。该中温干式厌氧发酵在春季和冬季具有较大的体系缓冲能力，TVFAs/碱度的比值分别为(0.23±0.13)和(0.29±0.17)，夏季和秋季时体系的缓冲能力相当，但最小缓冲能力出现在秋季为1.18，全年TVFAs/碱度大于1.0的天数仅2 d，均出现在秋季。这进一步表明，秋季运行调控对维持厌氧发酵工程整体产气稳定性至关重要。彭绪亚^[3]等通过餐厨垃圾的小试研究指出，当TVFAs/碱度< 0.4，厌氧发酵系统具有足够的缓冲能力；而当TVFAs/碱度= 0.4~0.8，厌氧发酵系统具备有限的缓冲能力；当TVFAs/碱度> 0.8时，则认为缓冲能力极小，有酸累积的可能。GAO等^[24]采用VFA/碱度评价厨余垃圾厌氧发酵性能及其稳定性。SONG等^[25]在污泥单相厌氧中发现，中温厌氧比高温厌氧的缓冲能力大。按照缓冲能力理论参考值分析，综合pH值变化，当缓冲能力值为1.18时，系统已经出现酸累积。在该项目的实际操作过程中，根据指标提示减少了50%进料量，避免造成厌氧菌大量失活，系统又快速恢复了平衡，并且工程上厌氧反应罐具备规模化效应，局部酸累积可由搅拌混合缓解而不影响整个厌氧发酵系统的稳定，能包容更多局部生物不平衡态，可耐受更高的VFAs，一定程度上也更易恢复系统平衡。

由表3可知，该干式厌氧发酵的产气能力与运行温度显著相关，春、夏、冬3季的平均温差与VS产气率显著相关。春季厌氧运行温度与VS产气率显著负相关，但与平均温差、OLR、pH和碱度表现出显著正相关；夏季厌氧运行温度与日产气量、容积产气率以及TVFA显著负相关，而与反映系统稳定性的参数pH、碱度以及缓冲能力(-TVFAs/碱度)的大小显著正相关；秋季和冬季的厌氧运行温度与产气效能、运行稳定性各参数均显著相关，秋季的温差与缓冲能力显著正相关，冬季的温差与VS产气率、日处理量、OLR和TVFAs显著相关。温差与厌氧运行温度基本上呈正相关关系。

干式厌氧发酵的稳定高效运行主要从运行效能和稳定性这2方面综合平衡，体系稳定性是提高运行效能的基础，高运行效能可实现更高的环境与经济效益。该干式厌氧发酵工程通常采用垃圾日处理量、进料TS含量、加热进料这3个指标的单个或组合变化指导日常操作和调控，使物料在厌氧发酵体系下稳定高效产气。依据环境温度变化，候气温15 ℃以下对中温干式厌氧发酵产气的影响较大。这主要是因为，自秋季降温期至初春升温期阶段(10-12月和1-3月)，此时干式厌氧发酵特征表现为秋季降温期，受厌氧温度逐步降低和夏季产气延后的影响，产气各参数指标间无显著相关性关系，且系统缓冲能力相对较小；冬季低温时，微生物活性降低，分解、产气速率均降低，增加处理量难以增加产气，反而使VS产气率等产气效能降低；春季升温期，微生物活性逐渐恢复，厌氧温度升，与体系的OLR呈正相关，体系发酵效能提高。而夏季厌氧运行温度较高时，有机质被微生物快速分解生成小分子酸、多糖等再进一步被利用产生沼气，容积产气率提高。进料TS因不同季节垃圾组成成份的显著差异而有较大波动，高TS时物料难以传质、传热均匀，导致在温差较大的冬季环境低温时，厌氧发酵局部低温而降低了产气能力，或者需要消耗更多能源使厌氧发酵正常运行。而在北方秋季的高TS中，木质素、纤维素等含量较高时，物料较难降解，为厌氧系统提供的可直接利用的有机物含量就少，故产气效能受到影响或延后产气。温度调控则体现为结合环境温度和厌氧发酵温度的全年变化规律，按季节性温差与厌氧值运行温度的关系，提出干式厌氧工程调控策略。

根据上文的分析以及各参数间的影响关系，在春季，中温干式厌氧发酵表现出较大的缓冲能力，可以适当增加处理量和提高厌氧发酵运行温度(增大温差)，既提高微生物活性和日产气量，又确保厌氧发酵系统的稳定性。在夏季，厌氧发酵运行温度较高，可增加处理量以提高厌氧发酵效能。由表3可知，日处理量增加，可提高日产气量和容积产气率，而pH值升高，表明调控没有导致厌氧发酵酸化。在秋季，日处理量最少，VS产气率最高，但体系缓冲能力较小。因此，为提高厌氧发酵效能，一方面可通过增加垃圾处理量提高容积产气率和日产气量；另一方面可增大温差即提高厌氧发酵运行温度，达到减少TVFAs、增加VS产气率的目的。冬季环境温度最低，温差最大，厌氧发酵运行平均温度为(27.20±2.99) ℃，pH缓冲能力充足，产气率有进一步提高空间。因此，根据表3分析结果，厌氧发酵运行温度为影响冬季各项参数的重要因素，建议在冬季提高厌氧发酵运行温度(增大温差)或提高回流比例，可提高干式厌氧体系的产气效能以及体系稳定性。

综上所述，为实现该中温干式厌氧发酵工程的稳定高效运行，在冬季低温期通过蒸汽加热提高运行温度的基础上，建议依据气温升温期(春季初期)和降温期(秋季后期)进行调控，在现有运行条件下提高春季初期和秋季后期的厌氧发酵运行温度；冬季还可考虑通过回流提升运行温度；夏季可维持现有运行温度，并通过提高日处理量提升处理效能。

1)中温生活垃圾干式厌氧发酵工程的有机容积负荷波动大，呈现冬季最高、秋季最低的特征；产气效能表现为容积产气率变化趋势与平均厌氧发酵运行温度一致，不同季节平均容积产气率从高到低依次为夏季(2.44±0.33 m³·(m³·d)⁻¹)>春季(2.31±0.28 m³·(m³·d)⁻¹)>冬季(2.25±0.29 m³·(m³·d)⁻¹)>秋季(2.10±0.14 m³·(m³·d)⁻¹)；VS产气率秋季最高，冬季最低。

2) 中温生活垃圾干式厌氧发酵过程呈现较高的TVFAs浓度和碱度浓度，大多数的TVFAs/碱度比值<0.8，pH在7.5~8.0波动，具备较大的pH缓冲能力。

3)干预条件下，中温生活垃圾干式厌氧发酵工程的全年运行温度在21~34 ℃，平均厌氧发酵运行温度表现为夏季高于春季和冬季，秋季最低。相关性分析结果表明，非严格受控中温干式厌氧发酵的产气能力与运行温度显著相关(p<0.05)，春、夏、冬3季的温差与VS产气率显著相关(p<0.05)。