本实验所用AnMBR装置如图1所示，反应器主体为有机玻璃材质的长方体形容器(长×宽×高=30 cm×20 cm×25 cm)，有效容积为11 L。反应器内设置聚偏氟乙烯(PVDF)材质的中空纤维膜(MBR-0.1F，海科滤膜)，膜组件构型采用U型，膜丝数量共80根，单根膜丝的内、外径分别为0.8 mm和1.8 mm，膜孔径为0.1 μm，膜组件总有效面积0.19 m²。反应器进出水采用蠕动泵控制(BT100-2J，兰格恒流泵有限公司)，出水泵与膜组件之间设有压力变送器(SIN-P300，中国杭州联测自动化技术有限公司)，用于监测跨膜压差的变化。反应器外部覆盖有加热带(鸿泰电热配件)，用于维持反应器内部的温度，内部温度由反应器顶部的温度计读取。反应器顶部设有阀门用于连接沼气袋和收集厌氧反应所产气体。反应器两侧均设有取样口用于采样分析。

AnMBR每个运行周期为10 min，其中出水8 min，松弛2 min。反应器设有2种不同的混合方式：从出水侧将泥水混合液泵回进水侧以进行循环混合；内部设置1个造浪泵(JVP-100AB)进行浪式脉冲混合。

反应器接种污泥为黑色厌氧颗粒污泥，取自升流式厌氧污泥床反应器。实验进水采用由葡萄糖、尿素、磷酸二氢钾等人工配制而成的模拟废水，具体水质成分如下：0.23 g·L^–1 尿素、0.11 g·L^–1 KH₂PO₄、1 g·L^–1 K₂HPO₄、0.05 g·L^–1 Na₂SO₄、0.1 g·L^–1 MgCl₂·6H₂O、0.05 g·L^–1 CaCl₂·2H₂O、10 mL·L^–1 微量元素、10 mL·L^–1 维生素(注：微量元素和维生素具体组成参考文献^[20])。缓冲物质为质量浓度为1 g·L^–1的碳酸氢钠。葡萄糖的用量根据实验进行调整。

控制AnMBR反应器温度为(35±1) ^oC，水力停留时间(HRT)为20 h，通过改变进水COD来改变反应器的负荷。反应器启动后共运行87 d，运行过程根据负荷和混合方式分为6个阶段，其中3个负荷水平(A、B、C)对应的进水COD分别为500、1 000、1 500 mg·L^–1，每个负荷下混合方式分为循环混合与浪式脉冲混合，各运行工况条件见表1。

每天取进出水水样测定COD值(哈希(Hach)试剂盒，2125815-CN)，取反应器内混合液样品用于DCH₄的浓度测定(见1.4)；沼气袋内收集气体的体积用注射器抽取测量，并用气相色谱(天美GC7900，中国)分析其组分及含量。各运行阶段结束时，取反应器内混合液用于生物量(MLSS和MLVSS)的测定(重量法)。

1)水中DCH₄的测定。水中的DCH₄采用顶空平衡法进行测定。在20 mL的密封小瓶中注入3.5 mol·L^–1的盐酸1 mL，用于抑制产甲烷菌的生理活性，从反应器内采集9 mL混合液样品注入小瓶中。将小瓶在25 ^oC的振荡器中振荡5 min，并放置3 h使CH₄达到气液平衡，平衡后取1 mL顶空气体注入气相色谱分析气体组成。反应器混合液中DCH₄的浓度根据亨利定律计算(式(1))^[21]所示。根据亨利定律推出式(2)。

式中：C_L为液体中DCH₄的质量浓度，mg·L^–1；C_G,eq是平衡后样品瓶顶空中CH₄的质量浓度，mg·L^–1；C_L,eq是样品瓶液体中与顶空平衡的CH₄质量浓度，mg·L^–1；V_G和V_L分别为样品瓶中液体和气体的体积，L。

式中：P_P,eq是CH₄在顶空中与液相平衡时的分压，Pa；${{M}}_{{\text{CH}}_{\text{4}}}$是CH₄的相对分子质量，16 g·mol^–1；H_cp为亨利常数，mol·(m³·Pa)^–1，在25 ℃时为1.4×10^–5 mol·(m³·Pa)^–1。

2) COD质量平衡计算。在各工况运行稳定后，根据式(3)^[18]进行COD的物料衡算，以分析有机物分布去向。

式中：C_in和C_eff分别为进出水COD，mg·L^–1；${\text{C}}_{{\text{DCH}}_{\text{4}}}$为DCH₄质量浓度(以COD计)，mg·L^–1；${\text{C}}_{{\text{GCH}}_{\text{4}}}$为以COD计气态CH₄质量浓度，mg·L^–1；C_other是其他所有包括生物质合成、微生物代谢产物(可溶性微生物产物(SMP)、胞外聚合物(EPS))生成和硫酸盐还原等用到的耗氧有机物(以COD计)，mg·L^–1。将产生的CH₄转化为COD(${\text{CH}}_{\text{4}}\text{+2}{\text{O}}_{\text{2}}\text{=}{\text{CO}}_{\text{2}}\text{+2}{\text{H}}_{\text{2}}\text{O}$)。值的注意的是，CO₂不会贡献COD，故不出现在COD质量平衡中^[22]。

3)CH₄传质系数的计算。AnMBR中CH₄由液相向气相的传质系数可由式(4)^[10]进行计算。

式中：K_La为由液相向气相的传质系数，h^–1；C_L为液体中DCH₄的质量浓度，mg·L^–1；Q_g为甲烷产气速率，L·h^–1；V是反应器的有效体积，L；P_CH4为AnMBR顶空中CH₄的分压，Pa；R是理想气体常数，8.31×10³ L·Pa·(mol·K)^–1；T是温度，K；Q_V (Q_g/V)为反应器的CH₄产气量，L·(L·h)^–1；C_eq是与气相中CH₄分压处于热力学平衡的DCH₄的质量浓度，mg·L^–1；K_H是CH₄的亨利常数，35 ^oC时为1.18×10^–7 mol·(L·Pa)^–1。

4)能量计算。CH₄的热值约为37.8 MJ·m^–3，按电能转换效率为35%计算CH₄转换成的电能(式(5))^[23]。反应器运行过程中混合所用设备的能量需求通过设备的功率和运行时间得出(式(6))。反应器运行过程中的净能量(E)由设备能源消耗与回收甲烷所得能源的差值计算(式(7))。

式中：${{W}}_{{\text{CH}}_{\text{4}}}$为CH₄转化的电能，kWh；${{V}}_{{\text{CH}}_{\text{4}}}$为CH₄的体积，m³

式中：W_消耗为设备消耗的能量，kWh；P为设备的功率，浪式脉冲设备功率为3 W，混合液循环设备功率为20 W；t为设备运行时间，h。

本研究中AnMBR共运行87 d，固定反应器的HRT为20 h，通过增加进水COD来提升负荷，在此期间研究了反应器不同混合形式对有机物去除性能的影响，各运行工况下进出水的COD值及去除率变化如图2所示。结果表明，AnMBR可实现对有机物的高效去除，COD平均去除率超过90%，提升进水负荷AnMBR对有机物仍可保持着良好的去除效果，具有一定的抗冲击负荷性。如图2所示，各负荷下AnMBR浪式脉冲混合时对有机物的去除效果均好于循环混合，在进水COD为500 mg·L^–1，负荷为0.6 kg·(m³·d)^–1时，循环混合模式下反应器出水的COD值为(48±4.9) mg·L^–1，对应COD去除率为(90.7±0.8)%；改用浪式脉冲混合模式后，COD去除率提升至(96.8±0.4)%。这是因为浪式脉冲混合的系统扰动强度更高，促进了微生物与底物的接触与利用，从而提高了有机物的去除效果^[24]。

在AnMBR的运行过程中，测定各运行工况下的CH₄产气量(Q_V)和DCH₄，结合亨利定律和式(4)计算出理论DCH₄和甲烷传质系数(K_La)。如图3所示，各阶段反应器液相中DCH₄均处于过饱和状态，DCH₄浓度为理论浓度的1.3~2.1倍，该水平位于已报道的AnMBR中DCH₄过饱和度范围之内^[25-26]。在相同负荷下，浪式脉冲混合时DCH₄的浓度和过饱和度均低于循环混合时(图3(c)~(d))。在负荷为0.6 kg·(m³·d)^–1下，循环混合时DCH₄为20.2 mg·L^–1，对应过饱和度为2.1；而浪式脉冲混合时DCH₄为12.6 mg·L^–1，对应过饱和度为1.3。上述结果是由于不同条件下CH₄的K_La不同导致的，K_La较大时，CH₄能更好地从液相扩散至气相，从而有利于降低DCH₄浓度及过饱和度。K_La通常与甲烷生成速率、混合效率、流体黏度和温度等因素有关^[10,27]。在本研究中，反应器内流体黏度与温度基本保持不变，而浪式脉冲混合时CH₄的K_La更大，这主要由2方面原因造成：一方面，浪式脉冲混合时体系的混匀度更高；另一方面，由上文可知，浪式脉冲混合时有机物的去除效果更好，对应CH₄的Qv更大。在相同负荷水平，YEO等采用沼气循环混合时DCH₄的过饱和度为2.2~2.5^[10]，对应K_La为(0.018±0.003)~(0.088±0.017) h^–1，由此可见，浪式脉冲混合在降低DCH₄浓度及过饱和度方面比液体循环或气体循环混合具有更好的效果。

如图3所示，随着进水负荷的提升，2种混合方式的DCH₄浓度均有升高。其中，循环混合时的DCH₄由20.2 mg·L^–1上升至23.4 mg·L^–1(图3(c))，对应DCH₄的过饱和度由2.1降至1.9；浪式脉冲混合时的DCH₄由12.6 mg·L^–1升至15.6 mg·L^–1(图3(d))，对应DCH₄的过饱和度由1.3降至1.2。可以看出，随着负荷的提升，DCH₄的过饱和度逐渐降低。这一结果与YEO等^[10]的研究结果一致。CH₄的Q_V随负荷的提升而增加(图3(a)~(b))，Qv与负荷呈现较好的正相关关系(R²=0.998)，WIJEKOON等^[28]的研究中也得出了相似的结论。因为CH₄的K_La正比于Q_V(式(4))，所以K_La随负荷提升也不断增大，且与负荷也呈现出较好的正相关关系(R²=0.997)。高负荷(1.8 kg·(m³·d)^–1)时CH₄的K_La ((0.57±0.03) h^–1)约为低负荷(0.6 kg·(m³·d)^–1)时K_La ((0.15±0.01) h^–1)的3~4倍。综上所述，CH₄的Qv和K_La均与负荷显著相关，随着负荷提升DCH₄不断升高，但DCH₄的变化程度(增加倍数)小于Qv，这是因为随负荷提升K_La也同时增大，导致更多的DCH₄扩散至气相，从而使DCH₄在过饱和度上呈现下降趋势。

在AnMBR的运行过程中，计算不同运行工况下CH₄在气/液相中的分布及在进水COD中的占比(见图4)。同负荷下，循环混合时气态CH₄占总CH₄及进水COD的比例均低于浪式脉冲混合时。如负荷为0.6 kg·(m³·d)^–1，循环混合时气态CH₄占总CH₄及进水COD的比例分别为76.0%和50.2%；浪式脉冲混合时气态CH₄占总CH₄及进水COD的比例分别为84.5%和57.0%。这是因为浪式脉冲混合(K_La更大)可促进CH₄从液相向气相的传质，这与上文的分析结果一致。

如图4所示，随着进水负荷的提升，DCH₄占总CH₄及进水COD的比例逐渐下降，而气态CH₄占总CH₄及进水COD的比例不断升高，表明AnMBR在高负荷下更容易对有机物中的能量进行回收。此外，AnMBR中可回收的气态CH₄占进水COD的比例(50.2%~60.0%)低于通常文献报道的CH₄含量70%~85%^[29-31]，这是由于本研究关注的处理对象为低浓度废水，所用的进水COD和负荷水平要低于上述报道，根据上文分析此时的K_La较低从而导致了较多的进水COD以DCH₄形式存在。

进水COD在AnMBR中的转化形式可能为气相CH₄、DCH₄、生物量增长和出水残留等，为了进一步阐明进水底物的电子流向，对各运行工况下AnMBR进行了COD质量衡算，结果见图5。由图5可知，总CH₄ (DCH₄+气态CH₄)占进水COD的主要部分，为61.0%~67.4%，DCH₄占比为4.3%~16.1%。以COD计的气态CH₄、DCH₄和出水COD加和后与进水COD的差值表示为C_other，在3种不同的负荷下的占比为24.4%~30.6%，这部分COD主要用于生物量合成、微生物代谢产物(SMP、EPS)的产生和微量的硫酸盐还原等，在其他AnMBR中也观察到了类似的COD转化^[17,32]。整个运行过程没有对AnMBR进行排泥，但反应器内生物量(MLVSS)变化较小，从6 020.4 mg·L^–1增加到6 820.7 mg·L^–1。生物量(MLVSS)与COD的换算关系^[33]约为1 g=1.4 g，计算后可知，用于生物增长的COD只占了C_other的3.1%~3.5%，以上结果表明大部分C_other可能以微生物代谢产物(SMP、EPS)的形式存在。微生物代谢产物主要为大分子的多糖和蛋白质^[34-35]，出水膜可能对这些化合物有很强的截留作用，进而导致了膜污染。关于不同负荷与混合方式对微生物代谢产物及膜污染的特征影响需要后续进一步的研究。

净能量平衡是评价系统应用时能量回收潜力的重要指标，通过能量需求与能量回收的差值来计算。本研究中，AnMBR运行过程保持HRT与温度恒定，各阶段的所需能量主要来自混合方式的电耗；各阶段的回收能量来自气相CH₄。各运行条件下的CH₄回收量和相关能量计算见图6和图7。

如图7所示，相比于循环混合，浪式脉冲混合时的CH₄回收量更大且总能量消耗更低。浪式脉冲混合的CH₄产量比循环混合条件下提高了6.5%~9.2%(图7)。由上文分析可知，这是由更多有机物转化和更高的气相CH₄导致。另一方面，由图6和图7可知，浪式脉冲混合所需的能耗只有循环混合的15%，但甲烷回收率提升了6.5%~9.2%，因此，计算分析后，浪式脉冲混合时AnMBR总能量消耗比循环混合时降低了85.9%~88.0%。随着负荷的提升，AnMBR总能量消耗不断降低，高负荷(1.8 kg·(m³·d)^–1)时总能耗比低负荷(0.6 kg·(m³·d)^–1时下降了约16.9%。由此可知，通过调整运行参数如缩减HRT等进一步提升负荷有利于实现AnMBR处理废水过程的能量平衡甚至能量盈余。ASLAM等^[36]在负荷为5.8 kg·(m³·d)^–1的条件下，在流态化厌氧陶瓷膜生物反应器中实现了0.099 kWh·m^–3的能量盈余，与本研究所得高负荷运行时AnMBR易实现能量平衡或盈余的结论一致。

1) HRT维持在20 h，进水负荷为0.6~1.8 kg·(m³·d)^–1的运行条件下，AnMBR可实现高效稳定的有机物去除效果(COD去除率>90%)，浪式脉冲混合时的有机物去除效果比循环混合时更好。

2) AnMBR在浪式脉冲混合时(K_La大)，CH₄能更好地从液相扩散到气相，其DCH₄的浓度和过饱和度及占进水COD的比例均比循环混合时低，且随负荷提升，这些指标进一步降低。

3) AnMBR进水大部分COD (61.0%~67.4%)转化为了CH₄，其中气相CH₄占比为50.2%~60.0%，剩余进水COD主要转化为微生物代谢产物。

4) AnMBR浪式脉冲混合的总能耗比循环混合时低85.9%~88.0%，提高负荷有利于实现AnMBR处理污水过程的能量平衡甚至能量盈余。