Loading [MathJax]/jax/output/HTML-CSS/jax.js

ABR-MBR工艺反硝化除磷微生物群落特征分析

蒋志云, 韦佳敏, 缪新年, 潘家成, 刘文如, 沈耀良. ABR-MBR工艺反硝化除磷微生物群落特征分析[J]. 环境工程学报, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
引用本文: 蒋志云, 韦佳敏, 缪新年, 潘家成, 刘文如, 沈耀良. ABR-MBR工艺反硝化除磷微生物群落特征分析[J]. 环境工程学报, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
JIANG Zhiyun, WEI Jiamin, MIAO Xinnian, PAN Jiacheng, LIU Wenru, SHEN Yaoliang. Analysis of microbial community characteristics of denitrifying phosphorus removal in the ABR-MBR process[J]. Chinese Journal of Environmental Engineering, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
Citation: JIANG Zhiyun, WEI Jiamin, MIAO Xinnian, PAN Jiacheng, LIU Wenru, SHEN Yaoliang. Analysis of microbial community characteristics of denitrifying phosphorus removal in the ABR-MBR process[J]. Chinese Journal of Environmental Engineering, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

ABR-MBR工艺反硝化除磷微生物群落特征分析

  • 蒋志云 1
  • 韦佳敏 1
  • 缪新年 1
  • 潘家成 1
  • 刘文如 1,2,3
  • 沈耀良 1,2,3
1. 苏州科技大学环境科学与工程学院,苏州 215009; 2. 江苏省水处理技术与材料协同创新中心,苏州 215009; 3. 江苏省环境科学与工程重点实验室,苏州 215009

中图分类号: X703.1

发布日期:2019-06-26

DOI:10.12030/j.cjee.201810005

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

为了揭示ABR-MBR组合工艺中反硝化除磷微生物种群演替规律,采用Miseq高通量测序技术考察了该工艺在不同运行阶段除磷功能区的微生物群落结构。结果表明,硝化液回流比逐步从150%提升至300%可促进反硝化除磷菌大量富集,促进系统的启动和稳定运行;系统在运行过程中始终保持较高的微生物多样性;优势微生物种群均以变形菌门(Proteobacteria)和拟杆菌门(Bacteroidetes)为主,最大丰度分别为55.13%和7.76%,且变形菌门功能性微生物主要集中在γ-变形菌纲(Gamaproteobacteria);功能性除磷菌属主要为气单胞菌属(Aeromonas),假单胞菌属(Pseudomonas)和黄杆菌属(Flavobacterium);其中在逐步提升硝化液回流比过程中气单胞菌属(Aeromonas)被大量富集,其在γ-变形菌纲(Gamaproteobacteria)的相对丰度由5.30%上升至41.49%并在系统后续运行中维持主导地位。系统除磷效果与功能性除磷微生物相对丰度的变化密切相关。系统中微生物种群的多样性和功能微生物的结构稳定性为ABR-MBR工艺的稳定运行和高效处理提供了保证。

低耗、高效的反硝化除磷技术一直是废水生物脱氮除磷领域的研究热点之[

1,2]。然而,已有研[3]过多地关注反硝化除磷工艺效能而忽略了在微生物学方面的深入研究。尽管以往有研究者对反硝化除磷菌(denitrifying polyphosphate-accumulating organisms,DPAOs)进行了研[4,5],但受限于当时的微生物检测技术,对反硝化除磷菌的研究仍不够深入。高通量测序技术作为新型微生物种群鉴定技术,具有分析结果准确、高速、高灵敏度和高自动化等特点,可以准确的分析生物多样性,故广泛应用于环境微生物鉴定领[6,7]。因此采用高通量测序技术作为考察反硝化除磷微生物群落演替规律的手段不失为一个良好的技术选择,但同时微生物在系统中的种群特征的研究也有待进一步的深入开展。已有的研究表明,在反硝化除磷系统中的一些优势菌属,如气单胞菌属(Aeromonas)、假单胞菌属(Pseudomonas)、肠杆菌科等都具有反硝化除磷功[8],它们的存在为解释系统效能的变化提供了理论支撑。

厌氧折流板反应器(anaerobic baffled reactor, ABR)具有构造简单、能耗低、抗冲击负荷能力强、处理效率高以及可实现生物相分离等优[

9,10],膜生物反应器(membrane bioreactor, MBR)则具有利于世代时间长的硝化细菌的增殖,从而提高硝化效率的优[11]。因此,采用ABR-MBR研究反硝化除磷工艺在培养DPAOs以及研究不同功能区微生物的细菌群落特征方面具有显著优[12]

本研究采用Miseq高通量测序手段对ABR-MBR组合工艺中除磷功能区的细菌群落特征进行分析,以期从微生物多样性、种群丰度及优势菌群分布方面为系统效能的变化提供理论支撑,为反硝化除磷工艺应用于工程实践提供理论依据,为开发新型的脱氮除磷工艺提供技术支持。

1 材料与方法

1.1 实验装置

本研究采用的ABR-MBR如图1所示,该一体化反应器由7格室ABR和好氧MBR串联而成,材质均采用有机玻璃,总有效容积11.4 L,其中ABR和MBR的有效容积分别为7.6 L和3.8 L。ABR中的A1格室进行水解发酵产生优质碳源,A2~A4格室均为厌氧释磷区,A5~A7格室均为缺氧吸磷区,MBR为好氧区。

图1 ABR-MBR工艺实验装置

Fig. 1 Schematic diagram of ABR-MBR process

1.2 实验水质和接种污泥

实验用水为苏州某高校生活污水,采用葡萄糖、淀粉、蛋白胨及碳酸氢铵、磷酸二氢钾作为补充碳源、氨氮和磷,如表1所示。

表1 ABR-MBR组合工艺进水水质 ( ABR-MBR combination process )
Table 1 Influent wastewater quality in ( ABR-MBR combination process )
水质指标范围均值
COD/(mg·L-1) 322~440 379
NH+4-N/(mg·L-1) 47.3~65.3 59.5
NO-2-N/(mg·L-1) 0~0.2 0.1
NO-3-N/(mg·L-1) 0~0.5 0.2
TN/(mg·L-1) 50.1~67.8 60.6
TP/(mg·L-1) 6.25~8.6 7.6
pH 7.1~7.4 7.2

本研究污泥样品取自ABR除磷功能区。为保证反硝化除磷菌适宜生存的温度,ABR水浴维持温度(30±2) ℃。启动初期,A1及A4接种取自本实验室自培养的具有良好除碳性能的污泥,MLSS约为7 800 mg·L-1,SVI为85 mL·g-1。MBR及除磷功能区(A2、A3、A5~A7)接种自苏州市某污水处理厂A2/O工艺的二沉池,具有良好的硝化及除磷能力,MLSS分别为3 200 mg·L-1和4 020 mg·L-1,SVI分别为80 mL·g-1 和75 mL·g-1,污泥沉降性能良好。

1.3 实验方案和运行参数

吕亮[

13]发现提高硝化液回流比对ABR-MBR系统中反硝化除磷效能有促进作用,故本研究系 统的启动主要采用固定污泥回流比,逐步提高硝化液回流比的方法。控制进水有机负荷(COD)为 0.8 kg·(m3·d)-1,ABR水力停留时间(HRT)为12 h,MBR的HRT为6 h,控制反硝化除磷功能区(A2、A3、A5~A7)污泥龄(SRT)为25 d,污泥回流比R1和R2为80%,pH在7.1~7.4,保持MBR的DO在1~ 2 mg·L-1,硝化液回流比R3从150%逐步提升至300%并最终稳定在300%。

1.4 水质指标的检测和方法

反应器运行期间,定期采集水样经0.45 μm中性滤纸过滤后,按照国家环保部规定的标准方法《水质和废水监测分析方法[

14]对COD、NH+4-N、NO-2-N、NO-3-N、PO3-4-P、TN等指标进行测定分析,COD采用快速消解法;NH+4-N采用纳氏试剂光度法;NO-2-N采用N-(1-萘基)-乙二胺分光光度法;NO-3-N采用紫外分光光度法;PO3-4-P采用钼锑抗分光光度法;TN采用过硫酸钾氧化-紫外分光光度法;MLSS采用滤纸称重法。

1.5 Miseq高通量测序分析方法

取启动初期以及稳定运行期的ABR除磷功能区的污泥,污泥经过多点采样混合作为该功能区样品,启动初期至稳定运行期(10、25和60 d)的污泥样品对应编号为M1、M2和M3。采用16S RNA 基因 V3~V4区通用引物(338F/806R),PCR扩增测序区域为338F/806R。引物名称和引物序列分别是338F(ACTCCTACGGGAGGCAGCAG)和806R(GGACTACHVGGGTWTCTAAT)[

15]。采样PCR仪(ABI GeneAmp 9700型)对细菌16S rRNA基因进行PCR扩增。

2 结果与讨论

2.1 ABR-MBR运行状况分析

ABR-MBR对系统正磷酸盐的去除效果如图2所示。由图2可知,进水磷酸盐浓度在6.25~ 8.65 mg·L-1,进水总氮浓度在50.11~67.78 mg·L-1,逐步提升硝化液回流比从150%至300%,在系统运行的第1个月,系统磷酸盐去除效果较差,ABR-MBR平均出水正磷酸盐浓度在4.11 mg·L-1左右,并且在启动初期的时候磷酸盐几乎没有去除。由于启动初期系统具备一定的缺氧吸磷能力,说明接种污泥中含有较多的聚磷菌并且含有一定量的DPAOs;系统运行至25 d时,对正磷酸盐的去除率为50%左右,对总氮的去除率则在30%左右,说明系统中已经存在一定量的反硝化除磷菌并且此时的条件是适宜其生存以及增长的;直至系统运行第40天之后,硝化液回流比稳定在300%,此时,ABR-MBR中平均出水正磷酸盐浓度由4.11 mg·L-1降为0.6 mg·L-1左右,平均出水总氮浓度由30.8 mg·L-1降为17.4 mg·L-1左右,系统对正磷酸盐和总氮的去除率分别为91.8%和70.2%,并且之后的正磷酸盐去除率一直维持在90%以上,平均出水磷浓度为0.45 mg·L-1,总氮的去除率一直维持在70%以上,平均出水氮浓度为15.6 mg·L-1,系统具有稳定高效的反硝化除磷能力。因此,逐步提升硝化液回流比是ABR-MBR组合工艺实现反硝化除磷的主要原因。

图2 ABR-MBR组合工艺中氮磷浓度的变化

Fig. 2 Variations of PO3-4-P and nitrogen concentration in ABR-MBR combined process

2.2 Miseq高通量测序结果分析

2.2.1 微生物种群多样性分析

利用Miseq平台对ABR除磷功能区启动初期至稳定运行期(10、25和60 d)的3个样品进行高通量测序,分别获得39 808、33 309、37 197条优化序列,将优化序列在97%的相似性下聚类,分别获得OTU 651、690、706个(见表2),且3个样品的物种覆盖度(coverage)均大于99%,表明了此次高通量测序的结果能够代表样本中微生物的真实情况。

表2 ABR反应器除磷功能区微生物种群丰度和多样性分析
Table 2 Analysis of microbial population abundance and diversity in phosphorus removal functional zone of ABR
样品序列数量/个OTU数量/个ACE指数Chao指数Shannon指数Simpson指数覆盖度/%
M1 39 808 651 678 681 4.483 0.052 99.79
M2 33 309 690 744 739 3.619 0.188 99.62
M3 37 197 706 758 770 3.659 0.193 99.64

Chao、ACE、Shannon、Simpson指数表明系统中除磷功能区细菌群落、物种的丰富度和多样性,其中丰富度指数Chao和ACE可以估算群落中含OTU数目的指数,反映系统种群的丰富度,在生态学中常用来估计物种总数,值越大代表物种总数越[

16]。由表2可知,3个样品的Chao、ACE指数均为M3>M2>M1,结果表明,在启动阶段至稳定运行期细菌群落的丰富度总体上呈上升趋势。Shannon指数反映基于物种数量的群落种类多样性,指数越大表明群落的复杂程度越[17],Shannon指数为M1>M3>M2,随着工艺的运行,Shannon指数在启动期至稳定前期出现了明显的下降,而在稳定期有所上升,这可能是由于接种污泥中的微生物种类较为丰富,随着DPAOs的富集,竞争能力弱的菌群将逐渐减少,造成群落结构逐渐集中化,但随着功能性反硝化除磷菌逐渐占据主导地位,且功能性微生物的多样性同比增加。Simpson指数也是用来估算样本中微生物多样性的指[7],该指数中稀有物种所起的作用较小,而普遍物种所起的作用较大。Simpson指数值越大,说明优势微生物占总生物量的比例越[18],如表2所示,Simpson指数逐渐增大,说明随着硝化液回流比的逐步提升,系统中优势微生物占总生物量比重增大。同时反应器运行25 d后较启动初期Simpson指数大幅度增大,当系统近乎稳定运行之后,Simpson指数缓慢上升,说明当反应器稳定运行之后,系统除磷功能区的优势微生物逐渐凸显并且微生物群落越来越稳定。由于该工艺具备的优势可为反硝化除磷菌提供良好的生存条件,使反硝化除磷菌逐渐富集并淘汰掉弱势菌群,可有效的进行脱氮除磷,因而在稳定运行期,系统对磷的去除率高达90%以上,对总氮的去除率高达72%以上。因此,随着硝化液回流比的逐步提升,系统除磷功能区的微生物种群的丰度指数和多样性指数均较高,并且系统的微生物的种群结构也越来越稳定,保证了ABR-MBR组合工艺的稳定运行。

2.2.2 门水平优势微生物分析

污泥样品在门水平的微生物种群组成如图3所示,系统中除磷功能区的细菌群落结构多样性较为丰富,主要存在8~9个门,按照丰富度排列主要存在变形菌门(Proteobacteria)、绿弯菌门(Chloroflexi)、厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、绿菌门(Chlorobi)等,并以变形菌门,绿弯菌门,厚壁菌门和拟杆菌门的微生物为主,四者的丰度比例约占总细菌的82.38%~84.32%,并且3个样品中的优势菌群均来自变形菌门。如图3所示,变形菌门的丰度分别为55.13%、29.51%、52.23%;绿弯菌门的丰度分别为8.36%、26.62%、13.98%;厚壁菌门丰度分别为11.13%、20.50%、9.59%;拟杆菌门的丰度分别为7.76%、7.69%、7.35%。从中可知,随着硝化液回流比的逐步提升,变形菌门一直为系统中的优势菌群,且占有系统最大丰度比例。国内外研究表明,大多数具有脱氮除磷功能以及反硝化除磷功能的微生物均归属于变形菌[

19],部分归属于拟杆菌[20],从门水平上的丰度变化来看,厚壁菌门在系统中也占据一定的比例,其中厚壁菌门中的芽孢杆菌具有较强的抵抗外界有害因子的能[21]并且其具有促硝化作用,能够降解废水中的氨[18],因此对废水中的氮的去除以及维持系统稳定具有较大的贡[22]

图3 微生物种群分类(门)的群落组成相对丰度比例

Fig. 3 Relative abundance percentage of microbial community at phylum level

本研究主要探讨除磷功能区的微生物种群结构,由测序结果可知,变形菌门和拟杆菌门为系统的优势菌群,并且具有除磷功能的微生物主要集中于变形菌门,少量存在于拟杆菌门。

为进一步探索反硝化除磷微生物所在的变形菌门在随着硝化液回流比逐步提升条件下微生物的 分布特征,故对除磷功能区变形菌门在纲分类水平的分布特征分别进行分析,结果见表3。3个样品中β-变形菌纲(Bataproteobacteria)和γ-变形菌纲(Gamaproteobacteria)是变形菌门中占据比例较大的菌群,占细菌总量的比例约为81.57%~88.26%。另外已有研[

23,24,25]表明,β-变形菌纲和γ-变形菌纲中均存在部分反硝化除磷菌,有利于去除废水中的磷和其他污染物,在系统生物脱氮除磷方面起到重要的作用,而变形菌门中δ-变形菌纲(Deltaproteobacteria)和α-变形菌纲(Alphaproteobacteria)的丰度相对较低。

表3 除磷功能区变形菌门微生物的群落组成相对百分比

Table 3 Percentage of microbial community of

Proteobacteria at class level in phosphorus

removal functional zone of ABR %

菌纲M1M2M3
Bataproteobacteria 74.29 23.23 68.87
Gamaproteobacteria 7.28 65.03 13.79
Deltaproteobacteria 16.38 6.50 14.40
Alphaproteobacteria 1.10 4.87 1.96

2.2.3 属水平优势微生物种群分析

为进一步阐明工艺运行过程中细菌群落的变化情况,在属的水平上,选取系统中占有比例最多的具有脱氮除磷功能的γ-变形菌纲进行分析。由图4可知,随着系统的运行,稳定运行期样品中功能性脱氮除磷微生物主要属类为气单胞菌属(Aeromonas)、假单胞菌属(Pseudomonas)、肠杆菌科细菌属(Enterobacteriaceae)等。此外,在β-变形菌纲还存在一定数量的ThaueraDenitratisoma和红环菌科细菌(Rhodocyclaceae),在拟杆菌门还检测出了少量的黄杆菌属(Flavobacterium)和鞘脂杆菌纲细菌。已有研[

26,27,28]显示,气单胞菌属、假单胞菌属、红环菌科细菌、黄杆菌属等均为常见的反硝化除磷优势菌;熊付[29]通过对反硝化除磷污泥特性的研究发现系统中主导菌属为Thauera, 并且发现该菌属具有反硝化功能;此外,Denitratisoma也是一种具有反硝化功能的菌[30],它们的存在为系统同时反硝化脱氮除磷效果提供了保障。

图4 γ-变形菌纲(Gamaproteobacteria)的群落组成(属)相对丰度比例

Fig. 4 Relative abundance percentage of microbial community of the γ-Proteobacteria phylum (Gamaproteobacteria)

由结果可知,当系统运行稳定之后,主要功能性除磷微生物主要来自变形菌门下的γ-变形菌纲,其中相对丰度最大的菌属为气单胞菌属,并且在启动初期至系统稳定运行的过程中,气单胞菌属在 γ-变形菌纲的相对丰度由5.30%上升至41.49%,与此同时,系统的正磷酸盐去除率逐步稳定在90%以上。此外,从图4可以看出,由M1至M2,Candidatus competibacter的数量显著增多,此时虽然系统中气单胞菌属的相对丰度达到了41.2%,但系统磷酸盐去除率却只有51.8%,据有研究报道,Candidatus competibacter是聚糖[

31]的一种,它的存在会与PAO竞争碳源从而不利于磷酸盐的去除,这也是导致系统脱氮除磷效果一直得不到提高的原因。随着硝化液回流比的逐步提升,Candidatus competibacter的相对丰度由44.5%下降至8.0%,系统除磷效果也逐渐提升。

图5所示为ABR-MBR中反硝化除磷菌菌属总量的分布规律。从启动初期(10 d)至稳定运行(80 d)过程中,反硝化除磷菌属占总生物量的比例在逐渐增加,由最初的1.57%上升至10.22%。相对来说,气单胞菌属的增加幅度较大,而其余功能性菌属增加的幅度普遍较小并且以较低的丰度水平存在于系统之中,也就是说系统中功能性除磷微生物丰度越来越大,同时随着系统中优势微生物丰度的增大,系统的磷酸盐去除率和总氮去除率也逐步上升并分别稳定在90%和72%以上。总的来说,功能性除磷菌属的丰度呈现大幅度上升的趋势,为系统整体脱氮除磷效果的提升提供了保证。

图5 ABR除磷功能区主要反硝化除磷菌菌属的分布

Fig. 5 Main DPAOs distribution in the

NOTE:

phosphorous removal functional zone of ABR

傅以钢[

32]对除磷工艺中的聚磷菌种群结构进行了研究,结果表明变形杆菌在除磷系统中占据了比较重要的地位,这与本研究中反硝化除磷工艺生物群落的特征是一致的,其在缺氧/好氧条件下属水平上占优势的菌属为Bdellovibrio sp.、Dechloromonas sp.、Neisseria sp.、Thauera sp.等,在厌氧/好氧条件下的优势菌属为Dechloromonas sp.、Pseudomonas sp.等,其中Thauera sp.以及Pseudomonas sp.菌属也是本研究中的优势菌属,而Bdellovibrio sp.Neisseria sp.等菌属为其工艺条件下独有的优势菌属。因此,反硝化除磷工艺与聚磷菌除磷工艺微生物群落虽有相似之处但也存在着差异。吴昌永[33]的研究结果也表明随着进水水质及工艺运行工况的改变,系统中微生物群落特征也会发生改变,同时微生物种群结构呈现出一种动态的演替过程。可见,工艺条件与微生物群落结构之间存在着显著的映射关系。

3 结论

1) 在ABR-MBR组合工艺中,控制温度为(30±2) ℃,pH为7.1~7.4,保持MBR中的DO在1~2 mg·L-1,逐步提升硝化液回流比从150%至300%,可促进反硝化除磷菌大量富集,实现系统的启动以及稳定运行。在硝化液回流比稳定在300%时,系统总磷和总氮去除率分别维持在90%和72%以上。

2) ABR-MBR工艺中反硝化除磷微生物以变形菌门(Proteobacteria)和拟杆菌门(Bacteroidetes)为主,分别占29.51%~55.13%和7.35%~7.76%,其中变形菌门功能性微生物主要集中在γ-变形菌纲(Gamaproteobacteria),以气单胞菌属(Aeromonas)和假单胞菌属(Pseudomonas)为主要组成;而拟杆菌门中以黄杆菌属(Flavobacterium)为主导。其中气单胞菌属(Aeromonas)增幅最大,其在γ-变形菌纲(Gamaproteobacteria)的相对丰度由5.30%上升至41.49%。

3) 工艺运行过程中细菌群落多样性逐渐增加,系统内部菌群结构也趋于稳定,系统中功能性除磷微生物的多样性和稳定性为系统脱氮除磷效果提供了保证。

4) ABR-MBR组合工艺运行过程中的去除效果与功能性反硝化除磷菌的相对丰度密切相关,当系统中反硝化除磷菌总量由1.57%上升至10.22%时,系统的总磷去除率上升至90%以上。

参 考 文 献

1

LV X M, SHAO M F, LI C L, et al. A comparative study of the bacterial community in denitrifying and traditional enhanced biological phosphorus removal processes[J]. Microbes & Environments, 2014, 29(4): 261-268.

2

吴继滨, 沈耀良, 郭海娟. SBR中短程反硝化除磷菌的培养驯化研究[J]. 环境科学与技术, 2011, 34(10): 127-129.

3

韦佳敏, 蒋志云, 程诚. ABR-MBR反硝化除磷工艺的启动及稳定运行[J]. 环境科学, 2019, 40(2): 808-815.

4

KUBA T, SMOLDERS G, LOOSDRECHT M C M V. Biological phosphorus removal from wastewater by anaerobic-anoxic sequencing batch reactor[J]. Water Science & Technology, 1993, 27(5/6): 241-252.

5

KUBA T, LOOSDRECHT M C M V, BRANDSE F A, et al. Occurrence of denitrifying phosphorus removing bacteria in modified UCT-type wastewater treatment plants[J]. Water Research, 1997, 31(4): 777-786.

6

夏围围, 贾仲君. 高通量测序和DGGE分析土壤微生物群落的技术评价[J]. 微生物学报, 2014, 54(12): 1489-1499.

7

汪瑶琪, 张敏, 姜滢, . 厌氧氨氧化启动过程及微生物群落结构特征[J]. 环境科学, 2017, 38(12): 5184-5191.

8

李相昆. 反硝化除磷工艺与微生物学研究[D]. 哈尔滨: 哈尔滨工业大学, 2006.

9

蒋永荣, 胡明成, 李学军, . ABR处理硫酸盐有机废水的相分离特性研究[J]. 环境科学, 2010, 31(7): 1544-1553.

10

程朝阳, 赵诗惠, 吕亮, . 基于ABR-MBR组合工艺优化反硝化除磷性能的研究[J]. 环境科学, 2016, 37(11): 4282-4288.

11

陈谊, 孙宝盛, 张斌, . 不同MBR反应器中硝化菌群落结构的研究[J]. 中国环境科学, 2010, 30(1): 69-75.

12

吴鹏, 陆爽君, 徐乐中, . ABR-MBR一体化工艺节能降耗措施优化研究[J]. 环境科学, 2015, 36(8): 2934-2938.

13

吕亮, 尤雯, 张敏, . 硝化液回流比对ABR-MBR工艺反硝化除磷效能的影响[J]. 环境科学, 2018, 39(3): 1309-1315.

14

国家环境保护总局. 水和废水监测分析方法[M]. 4版. 北京: 中国环境科学出版社, 2002.

15

ZHANG X L, TIAN X Q, MA L Y, et al. Biodiversity of the symbiotic bacteria associated with toxic marine dinoflagellate Alexandrium tamarense[J]. Journal of Biosciences and Medicines, 2015, 3(6): 23-28.

16

SHU D T, HE Y L, YUE H, et al. Metagenomic insights into the effects of volatile fatty acids on microbial community structures and functional genes in organotrophic anammox process[J]. Bioresource Technology, 2015, 196: 621-633

17

陈重军, 张海芹, 汪瑶琪, . 基于高通量测序的ABR厌氧氨氧化反应器各隔室细菌群落特征分析[J]. 环境科学, 2016, 37(7): 2652-2658.

18

赵诗惠, 吕亮, 蒋志云, . ABR-MBR组合工艺短程硝化过程的微生物种群[J]. 中国环境科学, 2018, 38(2): 566-573.

19

刘晖, 孙彦富, 周康群, . PCR-DGGE法研究Sludge bio-membrane(SB)系统中反硝化聚磷菌的变化[J]. 中南大学学报, 2011, 42(4): 1167-1174.

20

杨浩, 张国珍, 杨晓妮, . 16S rRNA高通量测序研究集雨窖水中微生物群落结构及多样性[J]. 环境科学, 2017, 38(4): 1704-1716.

21

NICHOLSON W L, MUNAKATA N, HORNECK G, et al. Resistance of Bacillus endospores to extreme terrestrial and extraterrestrial environments[J]. Microbiology & Molecular Biology Review, 2000, 64(3): 548-572.

22

YANG C, ZHANG W, LIU R, et al. Phylogenetic diversity and metabolic potential of activated sludge microbial communities in full-scale wastewater treatment plants[J]. Environmental Science & Technology, 2011, 45(17): 7408-7415.

23

JURETSCHKO S, LOY A, LEHNER A, et al. The microbial community composition of a nitrifying-denitrifying activated sludge from an industrial sewage treatment plant analyzed by the full-cycle rRNA approach[J]. Systematic & Applied Microbiology, 2002, 25(1): 84-99.

24

刘君寒, 胡光荣, 李福利, . 厌氧消化系统微生物菌群的研究进展[J]. 工业水处理, 2011, 31(10): 10-14.

25

WONG M T, MINO T, SEVIOUR R J, et al. In situ identification and characterization of the microbial community structure of full-scale enhanced biological phosphorous removal plants in Japan[J]. Water Research, 2005, 39(13): 2901-2914.

26

罗宁, 罗固源, 吉方英, . 新型双泥生物反硝化除磷脱氮系统中微生物的组成[J]. 给水排水, 2003, 29(8): 33-35.

27

刘立, 汤兵, 黄绍松, . 反硝化聚磷菌快速富集、培养及其荧光原位杂交技术鉴别[J]. 环境科学, 2013, 34(7): 2869-2875.

28

焦中志, 李相昆, 张立成, . 反硝化聚磷菌菌种筛分与除磷特性分析[J]. 沈阳建筑大学学报, 2009, 25(3): 535-540.

29

熊付娟. 反硝化除磷污泥除磷脱氮特性及菌群结构研究[D]. 哈尔滨哈尔滨工业大学, 2012.

30

常玉梅, 杨琦, 郝春博, . 城市污水厂活性污泥强化自养反硝化菌研究[J]. 环境科学, 2011, 32(4): 1210-1216.

31

CECH J S, HARTMAN P. Competition between polyphosphate and polysaccharide accumulating bacteria in enhanced biological phosphate removal systems[J]. Water Research, 1993, 27(7): 1219-1225.

32

傅以钢, 戴睿, 刘洪波, . 除磷工艺中含氧条件对聚磷菌种群结构影响研究[J]. 环境科学, 2008, 29(2): 474-481.

33

吴昌永, 彭永臻, 王淑莹, . 强化反硝化除磷对A2O工艺微生物种群变化的影响[J]. 化工学报, 2010, 61(1): 186-191.

蒋志云, 韦佳敏, 缪新年, 潘家成, 刘文如, 沈耀良. ABR-MBR工艺反硝化除磷微生物群落特征分析[J]. 环境工程学报, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
引用本文: 蒋志云, 韦佳敏, 缪新年, 潘家成, 刘文如, 沈耀良. ABR-MBR工艺反硝化除磷微生物群落特征分析[J]. 环境工程学报, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
JIANG Zhiyun, WEI Jiamin, MIAO Xinnian, PAN Jiacheng, LIU Wenru, SHEN Yaoliang. Analysis of microbial community characteristics of denitrifying phosphorus removal in the ABR-MBR process[J]. Chinese Journal of Environmental Engineering, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005
Citation: JIANG Zhiyun, WEI Jiamin, MIAO Xinnian, PAN Jiacheng, LIU Wenru, SHEN Yaoliang. Analysis of microbial community characteristics of denitrifying phosphorus removal in the ABR-MBR process[J]. Chinese Journal of Environmental Engineering, 2019, 13(7): 1653-1661. doi: 10.12030/j.cjee.201810005

ABR-MBR工艺反硝化除磷微生物群落特征分析

    通讯作者: 沈耀良(1961—), 男, 博士, 教授。研究方向:水污染控制与理论。E-mail:ylshen@mail.usts.edu.cn
    作者简介: 蒋志云(1993—), 女, 硕士研究生。研究方向:水污染控制与理论。E-mail:1712391710@qq.com
  • 1. 苏州科技大学环境科学与工程学院, 苏州 215009
  • 2. 江苏省水处理技术与材料协同创新中心, 苏州 215009
  • 3. 江苏省环境科学与工程重点实验室, 苏州 215009
基金项目:
国家自然科学基金面上资助项目(51578353)

摘要: 为了揭示ABR-MBR组合工艺中反硝化除磷微生物种群演替规律, 采用Miseq高通量测序技术考察了该工艺在不同运行阶段除磷功能区的微生物群落结构。结果表明, 硝化液回流比逐步从150%提升至300%可促进反硝化除磷菌大量富集, 促进系统的启动和稳定运行;系统在运行过程中始终保持较高的微生物多样性;优势微生物种群均以变形菌门(Proteobacteria)和拟杆菌门(Bacteroidetes)为主, 最大丰度分别为55.13%和7.76%, 且变形菌门功能性微生物主要集中在γ-变形菌纲(Gamaproteobacteria);功能性除磷菌属主要为气单胞菌属(Aeromonas), 假单胞菌属(Pseudomonas)和黄杆菌属(Flavobacterium);其中在逐步提升硝化液回流比过程中气单胞菌属(Aeromonas)被大量富集, 其在γ-变形菌纲(Gamaproteobacteria)的相对丰度由5.30%上升至41.49%并在系统后续运行中维持主导地位。系统除磷效果与功能性除磷微生物相对丰度的变化密切相关。系统中微生物种群的多样性和功能微生物的结构稳定性为ABR-MBR工艺的稳定运行和高效处理提供了保证。

English Abstract

  • 低耗、高效的反硝化除磷技术一直是废水生物脱氮除磷领域的研究热点之一[1-2]。然而, 已有研究[3]过多地关注反硝化除磷工艺效能而忽略了在微生物学方面的深入研究。尽管以往有研究者对反硝化除磷菌(denitrifying polyphosphate-accumulating organisms, DPAOs)进行了研究[4-5], 但受限于当时的微生物检测技术, 对反硝化除磷菌的研究仍不够深入。高通量测序技术作为新型微生物种群鉴定技术, 具有分析结果准确、高速、高灵敏度和高自动化等特点, 可以准确的分析生物多样性, 故广泛应用于环境微生物鉴定领域[6-7]。因此采用高通量测序技术作为考察反硝化除磷微生物群落演替规律的手段不失为一个良好的技术选择, 但同时微生物在系统中的种群特征的研究也有待进一步的深入开展。已有的研究表明, 在反硝化除磷系统中的一些优势菌属, 如气单胞菌属(Aeromonas)、假单胞菌属(Pseudomonas)、肠杆菌科等都具有反硝化除磷功能[8], 它们的存在为解释系统效能的变化提供了理论支撑。

    厌氧折流板反应器(anaerobic baffled reactor, ABR)具有构造简单、能耗低、抗冲击负荷能力强、处理效率高以及可实现生物相分离等优点[9-10], 膜生物反应器(membrane bioreactor, MBR)则具有利于世代时间长的硝化细菌的增殖, 从而提高硝化效率的优点[11]。因此, 采用ABR-MBR研究反硝化除磷工艺在培养DPAOs以及研究不同功能区微生物的细菌群落特征方面具有显著优势[12]

    本研究采用Miseq高通量测序手段对ABR-MBR组合工艺中除磷功能区的细菌群落特征进行分析, 以期从微生物多样性、种群丰度及优势菌群分布方面为系统效能的变化提供理论支撑, 为反硝化除磷工艺应用于工程实践提供理论依据, 为开发新型的脱氮除磷工艺提供技术支持。

  • 本研究采用的ABR-MBR如图 1所示, 该一体化反应器由7格室ABR和好氧MBR串联而成, 材质均采用有机玻璃, 总有效容积11.4 L, 其中ABR和MBR的有效容积分别为7.6 L和3.8 L。ABR中的A1格室进行水解发酵产生优质碳源, A2~A4格室均为厌氧释磷区, A5~A7格室均为缺氧吸磷区, MBR为好氧区。

  • 实验用水为苏州某高校生活污水, 采用葡萄糖、淀粉、蛋白胨及碳酸氢铵、磷酸二氢钾作为补充碳源、氨氮和磷, 如表 1所示。

    本研究污泥样品取自ABR除磷功能区。为保证反硝化除磷菌适宜生存的温度, ABR水浴维持温度(30±2) ℃。启动初期, A1及A4接种取自本实验室自培养的具有良好除碳性能的污泥, MLSS约为7 800 mg·L-1, SVI为85 mL·g-1。MBR及除磷功能区(A2、A3、A5~A7)接种自苏州市某污水处理厂A2/O工艺的二沉池, 具有良好的硝化及除磷能力, MLSS分别为3 200 mg·L-1和4 020 mg·L-1, SVI分别为80 mL·g-1和75 mL·g-1, 污泥沉降性能良好。

  • 吕亮等[13]发现提高硝化液回流比对ABR-MBR系统中反硝化除磷效能有促进作用, 故本研究系统的启动主要采用固定污泥回流比, 逐步提高硝化液回流比的方法。控制进水有机负荷(COD)为0.8 kg·(m3·d)-1, ABR水力停留时间(HRT)为12 h, MBR的HRT为6 h, 控制反硝化除磷功能区(A2、A3、A5~A7)污泥龄(SRT)为25 d, 污泥回流比R1和R2为80%, pH在7.1~7.4, 保持MBR的DO在1~2 mg·L-1, 硝化液回流比R3从150%逐步提升至300%并最终稳定在300%。

  • 反应器运行期间, 定期采集水样经0.45 μm中性滤纸过滤后, 按照国家环保部规定的标准方法《水质和废水监测分析方法》[14]对COD、NH4+-N、NO2--N、NO3--N、PO43--P、TN等指标进行测定分析, COD采用快速消解法; NH4+-N采用纳氏试剂光度法; NO2--N采用N-(1-萘基)-乙二胺分光光度法; NO3--N采用紫外分光光度法; PO43--P采用钼锑抗分光光度法; TN采用过硫酸钾氧化-紫外分光光度法; MLSS采用滤纸称重法。

  • 取启动初期以及稳定运行期的ABR除磷功能区的污泥, 污泥经过多点采样混合作为该功能区样品, 启动初期至稳定运行期(10、25和60 d)的污泥样品对应编号为M1、M2和M3。采用16S RNA基因V3~V4区通用引物(338F/806R), PCR扩增测序区域为338F/806R。引物名称和引物序列分别是338F(ACTCCTACGGGAGGCAGCAG)和806R(GGACTACHVGGGTWTCTAAT)[15]。采样PCR仪(ABI GeneAmp 9700型)对细菌16S rRNA基因进行PCR扩增。

  • ABR-MBR对系统正磷酸盐的去除效果如图 2所示。由图 2可知, 进水磷酸盐浓度在6.25~ 8.65 mg·L-1, 进水总氮浓度在50.11~67.78 mg·L-1, 逐步提升硝化液回流比从150%至300%, 在系统运行的第1个月, 系统磷酸盐去除效果较差, ABR-MBR平均出水正磷酸盐浓度在4.11 mg·L-1左右, 并且在启动初期的时候磷酸盐几乎没有去除。由于启动初期系统具备一定的缺氧吸磷能力, 说明接种污泥中含有较多的聚磷菌并且含有一定量的DPAOs; 系统运行至25 d时, 对正磷酸盐的去除率为50%左右, 对总氮的去除率则在30%左右, 说明系统中已经存在一定量的反硝化除磷菌并且此时的条件是适宜其生存以及增长的; 直至系统运行第40天之后, 硝化液回流比稳定在300%, 此时, ABR-MBR中平均出水正磷酸盐浓度由4.11 mg·L-1降为0.6 mg·L-1左右, 平均出水总氮浓度由30.8 mg·L-1降为17.4 mg·L-1左右, 系统对正磷酸盐和总氮的去除率分别为91.8%和70.2%, 并且之后的正磷酸盐去除率一直维持在90%以上, 平均出水磷浓度为0.45 mg·L-1, 总氮的去除率一直维持在70%以上, 平均出水氮浓度为15.6 mg·L-1, 系统具有稳定高效的反硝化除磷能力。因此, 逐步提升硝化液回流比是ABR-MBR组合工艺实现反硝化除磷的主要原因。

  • 利用Miseq平台对ABR除磷功能区启动初期至稳定运行期(10、25和60 d)的3个样品进行高通量测序, 分别获得39 808、33 309、37 197条优化序列, 将优化序列在97%的相似性下聚类, 分别获得OTU 651、690、706个(见表 2), 且3个样品的物种覆盖度(coverage)均大于99%, 表明了此次高通量测序的结果能够代表样本中微生物的真实情况。

    Chao、ACE、Shannon、Simpson指数表明系统中除磷功能区细菌群落、物种的丰富度和多样性, 其中丰富度指数Chao和ACE可以估算群落中含OTU数目的指数, 反映系统种群的丰富度, 在生态学中常用来估计物种总数, 值越大代表物种总数越多[16]。由表 2可知, 3个样品的Chao、ACE指数均为M3 > M2 > M1, 结果表明, 在启动阶段至稳定运行期细菌群落的丰富度总体上呈上升趋势。Shannon指数反映基于物种数量的群落种类多样性, 指数越大表明群落的复杂程度越高[17], Shannon指数为M1 > M3 > M2, 随着工艺的运行, Shannon指数在启动期至稳定前期出现了明显的下降, 而在稳定期有所上升, 这可能是由于接种污泥中的微生物种类较为丰富, 随着DPAOs的富集, 竞争能力弱的菌群将逐渐减少, 造成群落结构逐渐集中化, 但随着功能性反硝化除磷菌逐渐占据主导地位, 且功能性微生物的多样性同比增加。Simpson指数也是用来估算样本中微生物多样性的指数[7], 该指数中稀有物种所起的作用较小, 而普遍物种所起的作用较大。Simpson指数值越大, 说明优势微生物占总生物量的比例越大[18], 如表 2所示, Simpson指数逐渐增大, 说明随着硝化液回流比的逐步提升, 系统中优势微生物占总生物量比重增大。同时反应器运行25 d后较启动初期Simpson指数大幅度增大, 当系统近乎稳定运行之后, Simpson指数缓慢上升, 说明当反应器稳定运行之后, 系统除磷功能区的优势微生物逐渐凸显并且微生物群落越来越稳定。由于该工艺具备的优势可为反硝化除磷菌提供良好的生存条件, 使反硝化除磷菌逐渐富集并淘汰掉弱势菌群, 可有效的进行脱氮除磷, 因而在稳定运行期, 系统对磷的去除率高达90%以上, 对总氮的去除率高达72%以上。因此, 随着硝化液回流比的逐步提升, 系统除磷功能区的微生物种群的丰度指数和多样性指数均较高, 并且系统的微生物的种群结构也越来越稳定, 保证了ABR-MBR组合工艺的稳定运行。

  • 污泥样品在门水平的微生物种群组成如图 3所示, 系统中除磷功能区的细菌群落结构多样性较为丰富, 主要存在8~9个门, 按照丰富度排列主要存在变形菌门(Proteobacteria)、绿弯菌门(Chloroflexi)、厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、绿菌门(Chlorobi)等, 并以变形菌门, 绿弯菌门, 厚壁菌门和拟杆菌门的微生物为主, 四者的丰度比例约占总细菌的82.38%~84.32%, 并且3个样品中的优势菌群均来自变形菌门。如图 3所示, 变形菌门的丰度分别为55.13%、29.51%、52.23%;绿弯菌门的丰度分别为8.36%、26.62%、13.98%;厚壁菌门丰度分别为11.13%、20.50%、9.59%;拟杆菌门的丰度分别为7.76%、7.69%、7.35%。从中可知, 随着硝化液回流比的逐步提升, 变形菌门一直为系统中的优势菌群, 且占有系统最大丰度比例。国内外研究表明, 大多数具有脱氮除磷功能以及反硝化除磷功能的微生物均归属于变形菌门[19], 部分归属于拟杆菌门[20], 从门水平上的丰度变化来看, 厚壁菌门在系统中也占据一定的比例, 其中厚壁菌门中的芽孢杆菌具有较强的抵抗外界有害因子的能力[21]并且其具有促硝化作用, 能够降解废水中的氨氮[18], 因此对废水中的氮的去除以及维持系统稳定具有较大的贡献[22]

    本研究主要探讨除磷功能区的微生物种群结构, 由测序结果可知, 变形菌门和拟杆菌门为系统的优势菌群, 并且具有除磷功能的微生物主要集中于变形菌门, 少量存在于拟杆菌门。

    为进一步探索反硝化除磷微生物所在的变形菌门在随着硝化液回流比逐步提升条件下微生物的分布特征, 故对除磷功能区变形菌门在纲分类水平的分布特征分别进行分析, 结果见表 3。3个样品中β-变形菌纲(Bataproteobacteria)和γ-变形菌纲(Gamaproteobacteria)是变形菌门中占据比例较大的菌群, 占细菌总量的比例约为81.57%~88.26%。另外已有研究[23-25]表明, β-变形菌纲和γ-变形菌纲中均存在部分反硝化除磷菌, 有利于去除废水中的磷和其他污染物, 在系统生物脱氮除磷方面起到重要的作用, 而变形菌门中δ-变形菌纲(Deltaproteobacteria)和α-变形菌纲(Alphaproteobacteria)的丰度相对较低。

  • 为进一步阐明工艺运行过程中细菌群落的变化情况, 在属的水平上, 选取系统中占有比例最多的具有脱氮除磷功能的γ-变形菌纲进行分析。由图 4可知, 随着系统的运行, 稳定运行期样品中功能性脱氮除磷微生物主要属类为气单胞菌属(Aeromonas)、假单胞菌属(Pseudomonas)、肠杆菌科细菌属(Enterobacteriaceae)等。此外, 在β-变形菌纲还存在一定数量的ThaueraDenitratisoma和红环菌科细菌(Rhodocyclaceae), 在拟杆菌门还检测出了少量的黄杆菌属(Flavobacterium)和鞘脂杆菌纲细菌。已有研究[26-28]显示, 气单胞菌属、假单胞菌属、红环菌科细菌、黄杆菌属等均为常见的反硝化除磷优势菌; 熊付娟[29]通过对反硝化除磷污泥特性的研究发现系统中主导菌属为Thauera, 并且发现该菌属具有反硝化功能; 此外, Denitratisoma也是一种具有反硝化功能的菌属[30], 它们的存在为系统同时反硝化脱氮除磷效果提供了保障。

    由结果可知, 当系统运行稳定之后, 主要功能性除磷微生物主要来自变形菌门下的γ-变形菌纲, 其中相对丰度最大的菌属为气单胞菌属, 并且在启动初期至系统稳定运行的过程中, 气单胞菌属在γ-变形菌纲的相对丰度由5.30%上升至41.49%, 与此同时, 系统的正磷酸盐去除率逐步稳定在90%以上。此外, 从图 4可以看出, 由M1至M2, Candidatus competibacter的数量显著增多, 此时虽然系统中气单胞菌属的相对丰度达到了41.2%, 但系统磷酸盐去除率却只有51.8%, 据有研究报道, Candidatus competibacter是聚糖菌[31]的一种, 它的存在会与PAO竞争碳源从而不利于磷酸盐的去除, 这也是导致系统脱氮除磷效果一直得不到提高的原因。随着硝化液回流比的逐步提升, Candidatus competibacter的相对丰度由44.5%下降至8.0%, 系统除磷效果也逐渐提升。

    图 5所示为ABR-MBR中反硝化除磷菌菌属总量的分布规律。从启动初期(10 d)至稳定运行(80 d)过程中, 反硝化除磷菌属占总生物量的比例在逐渐增加, 由最初的1.57%上升至10.22%。相对来说, 气单胞菌属的增加幅度较大, 而其余功能性菌属增加的幅度普遍较小并且以较低的丰度水平存在于系统之中, 也就是说系统中功能性除磷微生物丰度越来越大, 同时随着系统中优势微生物丰度的增大, 系统的磷酸盐去除率和总氮去除率也逐步上升并分别稳定在90%和72%以上。总的来说, 功能性除磷菌属的丰度呈现大幅度上升的趋势, 为系统整体脱氮除磷效果的提升提供了保证。

    傅以钢等[32]对除磷工艺中的聚磷菌种群结构进行了研究, 结果表明变形杆菌在除磷系统中占据了比较重要的地位, 这与本研究中反硝化除磷工艺生物群落的特征是一致的, 其在缺氧/好氧条件下属水平上占优势的菌属为Bdellovibrio sp.、Dechloromonas sp.、Neisseria sp.、Thauera sp.等, 在厌氧/好氧条件下的优势菌属为Dechloromonas sp.、Pseudomonas sp.等, 其中Thauera sp.以及Pseudomonas sp.菌属也是本研究中的优势菌属, 而Bdellovibrio sp.和Neisseria sp.等菌属为其工艺条件下独有的优势菌属。因此, 反硝化除磷工艺与聚磷菌除磷工艺微生物群落虽有相似之处但也存在着差异。吴昌永等[33]的研究结果也表明随着进水水质及工艺运行工况的改变, 系统中微生物群落特征也会发生改变, 同时微生物种群结构呈现出一种动态的演替过程。可见, 工艺条件与微生物群落结构之间存在着显著的映射关系。

  • 1) 在ABR-MBR组合工艺中, 控制温度为(30±2) ℃, pH为7.1~7.4, 保持MBR中的DO在1~2 mg·L-1, 逐步提升硝化液回流比从150%至300%, 可促进反硝化除磷菌大量富集, 实现系统的启动以及稳定运行。在硝化液回流比稳定在300%时, 系统总磷和总氮去除率分别维持在90%和72%以上。

    2) ABR-MBR工艺中反硝化除磷微生物以变形菌门(Proteobacteria)和拟杆菌门(Bacteroidetes)为主, 分别占29.51%~55.13%和7.35%~7.76%, 其中变形菌门功能性微生物主要集中在γ-变形菌纲(Gamaproteobacteria), 以气单胞菌属(Aeromonas)和假单胞菌属(Pseudomonas)为主要组成; 而拟杆菌门中以黄杆菌属(Flavobacterium)为主导。其中气单胞菌属(Aeromonas)增幅最大, 其在γ-变形菌纲(Gamaproteobacteria)的相对丰度由5.30%上升至41.49%。

    3) 工艺运行过程中细菌群落多样性逐渐增加, 系统内部菌群结构也趋于稳定, 系统中功能性除磷微生物的多样性和稳定性为系统脱氮除磷效果提供了保证。

    4) ABR-MBR组合工艺运行过程中的去除效果与功能性反硝化除磷菌的相对丰度密切相关, 当系统中反硝化除磷菌总量由1.57%上升至10.22%时, 系统的总磷去除率上升至90%以上。

参考文献 (33)

返回顶部

目录

/

返回文章
返回