据统计,至“十二五”期末我国湿污泥量将突破4600万t,而污泥厌氧消化技术以其低能耗、高产出的经济优势成为污泥资源化利用的主要技术之一.除厌氧消化产甲烷以外,污泥产挥发性脂肪酸(fatty acids)也是实现污泥资源化的有效途径,近年来,越来越多的学者开始关注污泥厌氧发酵产挥发性脂肪酸.目前,有关污泥厌氧发酵产酸的研究主要集中在通过改进装置构型、产酸微生物生态、优化控制运行条件,如控温、pH等条件因素来提高产酸效率.已有研究表明,通过调节发酵污泥底物的C/N比可增加发酵产酸量并调控其产酸类型,然而,目前研究人员对污泥厌氧发酵产酸过程中不同C/N比与关键酶酶活及有机酸产酸量间的关系并不清楚.仅有为数不多的研究,如优化C/N比条件作为酒精发酵的实验模型研究,而对于数学模型则没有报道.数学模型法作为现代科学研究的重要手段,它有助于描述和理解生物处理系统的反应过程,可为工程设计提供理论上的指导;还有助于工艺的优化和控制,从而更好地指导实际生产运行.
多元线性回归是一种理想的描述多个因素(factor)之间关系的数学方法,能较好地确定被解释变量和解释变量之间的关系,在很多领域得到了应用.因此,本研究通过设置不同C/N比条件来调控污泥厌氧发酵产酸,在Matlab7.0平台上建立多元线性回归函数模型,拟合C/N比、关键酶酶活和产酸类型之间的关系,以期为今后污泥发酵产酸条件调控研究和工程(Engineering)放大提供参考.
2 材料与方法
2.1 实验材料
2.1.1 污泥与种泥
原始污泥取自无锡市太湖新城污水处理厂,发酵底物是经过热碱预处理的污泥上清液.污泥采集后置于阴凉处,风干10 d,采用机械打碎仪粉碎,再过30目筛,密封置于-15 ℃冰柜中保存.
接种污泥来源于无锡某柠檬酸厂上流式厌氧污泥反应器中的厌氧颗粒污泥.在100 ℃下煮沸2 h以杀死产甲烷菌,然后导入有效容积为2 L的UASB中进行驯化,活化种泥中的产酸菌,驯化温度为35 ℃.每日监测驯化种泥的pH值,待种泥驯化后出水pH值降低至4.0左右,稳定3~5 d后,认为种泥驯化成功.原始（Original)污泥和污泥预处理(chǔ lǐ)液及接种种泥的性质见表 1.
表 1 原始污泥、污泥预处理(chǔ lǐ)液和接种污泥的性质
2.1.2 厌氧发酵
调节热-碱预处理后离心液pH为10.0,取500 mL离心液置于1000 mL的厌氧反应瓶中,分别加入不同量葡萄糖,以使得底物混合液的初始C/N比为12、56、69、156.接入驯化后的种泥,种泥接种量为10%.污泥发酵前充氮气10 min以去除氧气,然后迅速密封置于转速为120 r?min-1和温度℃的摇床厌氧发酵.在发酵期间,每12 h调节pH为初始的10.0.每24 h取样1次,用针管吸出部分发酵液,取样完成后调节pH并充氮气保持厌氧.
2.1.3 实验药品和材料
本实验采用的药品包括4-甲基戊酸、磷酸(化学式： H3PO4)溶液、NaO
H、HCL等;主要仪器包括pH计、气相色谱仪、马弗炉、凯氏定氮仪,离心机等.
2.2 试验方法 2.2.1 污泥初始指标测定
污泥预处理前后的总固体、溶解性固体及污泥挥发性固体的测定采用重量法,具体操作详见GB11901-89和《水和废水监测分析方法》.pH值测定参照国标法.
污泥中的总脂类物质采用Bligh-Dyer方法提取后,在 80 ℃下干燥直至溶剂完全挥发后,采用重量法测定.总氮采用凯氏定氮法测定.总蛋白含量通过凯氏氮减去氨氮后再乘以 6.25计算得到.氨氮采用纳氏试剂比色法测定,污泥中的总碳水化合物采用甲醛离心法提取后,再用苯酚-硫酸法测定.用 Liquid TOC 分析(Analyse)仪测定总有机碳,详见《水质总有机碳的测定燃烧氧化-非分散红外吸收法》.污泥中的总磷(P)含量用钼酸铵分光光度法测定,详见《水质总磷的测定钼酸铵分光光度法》.
2.2.2 关键酶活测定
乙酸激酶的活性采用文献的方法提取并测定.磷酸转移乙酰酶的提取方法同乙酸激酶,活性测定参照文献方法.丁酸激酶微生物细胞的破壁方法和提取方法同乙酸激酶,活性测定采用文献方法.磷酸转移丁酰酶微生物细胞的破壁和提取方法同乙酸激酶,活性测定采用文献方法. 甲基丙二酰CoA变位酶活性测定采用文献方法.
2.2.3 挥发性短链脂肪酸的测定
采用GC法检测挥发性短链脂肪(fat)酸的质量浓度,样品处理及色谱条件等参见文献.为方便不同条件下产酸效率的比较,将测得的VFAs浓度折算成 COD值,换算方法参见文献.
3 结果与讨论
3.1 有机酸浓度的变化
分别设定底物初始C/N比值为12、56、156,进行厌氧发酵.发酵过程中时,体系中有机酸的产量分别如图 1所示.由图 1可以看出,底物初始C/N比不同,厌氧发酵产生的末端酸化产物也不同.C/N比为12时,产量最大的为乙酸,在第5 d达到9.45 kg?m-3;其次是丙酸,第5 d时可以达到3.55 kg?m-3;最低是丁酸,第5 d时产量约为2.35 kg?m-3.当C/N比为56时,产量最大的为丙酸,在第5 d可达到10.36 kg?m-3;其次是乙酸,可以达到7.79 kg?m-3;最低是丁酸,产量约为2.79 kg?m-3.当C/N比为156时,产量最大的为丁酸,在第5 d达到13.59 kg?m-3;其次是乙酸,可以达到5.89 kg?m-3;最低是丁酸,产量约为4.72 kg?m-3.
图 1 C/N比对发酵产酸的影响
图 2 底物发酵产酸的代谢途径
3.2 多项式关系的数学模型的建立
在发酵过程中,通过设定不同C/N比条件,测定不同C/N比下关键（解释:比喻事物的重要组成部分)酶的酶活和产酸量,则三者可以建立函数关系.短链脂肪酸的生成途径如图 2所示,乙酸生成的关键酶分别为乙酸激酶和磷酸转移乙酰酶,丁酸合成的关键酶有丁酸激酶和磷酸转移丁酰酶;丙酮酸转化为乙酸过程中,关键酶为甲基丙二酰CoA变位酶.由于底物C/N比的改变会导致产酸微生物体内的酶活性改变,从而改变微生物不同代谢途径的代谢通量,并最终导致各种短链脂肪酸的生成受到影响,由此产生了不同的产酸类型.因此,C/N比是自变量,而关键酶活性和产酸类型是因变量.
根据以上理论分析,为了建立C/N比、关键酶活、不同酸产量之间的函数关系模型,采用多元非线性回归的方法,建立了二元二次多项式模型、二元三次多项式模型、二元三次多项式模型.
式中,Z为因变量,表示有机酸产量,X和Y为自变量,分别表示C/N比和关键酶活性,bi表示函数中的常数.根据实测数据对以上3种多项式模型进行拟合优度的检验,结果如表 2所示.
表 2 模型拟合优度分析
表 2中,
A、
B、C分别表示3种关系式模型,乙酸-A
K、乙酸-PTA等分别表示各种产物酸与其对应的关键酶.根据优度拟合理论,R2大于 0.9的较好,因此,在多项式A乙酸-A
K、B乙酸-AK和C乙酸-AK的R2检验中,C乙酸-AK效果最好,为0.9327.同样的,C乙酸-PT
A、C丙酸-MC
M、C丁酸-B
K、C丁酸-PTB的R2分别为0.9348、0.9494、0.9880和0.9771,均具有最高的R2.模型拟合优度检验同时要求模型残差平方和越小越好,因此,在多项式A乙酸-A
K、B乙酸-AK和C乙酸-AK的残差平方和检验中,C乙酸-AK的效果最好,为0.0292.同样的,C乙酸-PT
A、C丙酸-MC
M、C丁酸-B
K、C丁酸-PTB的残差平方和分别为0.0076、1.8556、0.0038和3.1362,均具有最小的残差平方和.综合以上分析,可以判定二元三次多项式C模型拟合度最好.依据二元三次多项式C标准结构,结合实验测得的不同C/N比条件下的关键酶活和各种有机酸产量,建立了可视化的三者之间的曲面模型.
3.3 C/N比-关键酶-乙酸曲面模型的建立
根据上述多项式C的形式,利用Matlab软件进行拟合,得到多项式当中相应bi的值,进而得出C/N比-关键（解释:比喻事物的重要组成部分)酶AK-乙酸之间的定量关系表达式和C/N比-关键酶PTA-乙酸之间的定量关系表达式.
根据定量关系式和绘制得到图 3的可视化曲面模型.从图 3的曲面模型可以看出：在初始C/N比为10~50时,酶活性较低,乙酸激酶和磷酸转乙酰酶平均酶活性分别为1.06 U?mg-1和0.67 U?mg-1,而乙酸浓度从2.8 kg?m-3上升到7.8 kg?m-3.在初始C/N比为50~150时,AK和PTA酶活性较高,平均值分别为2.79 U?mg-1和1.08 U?mg-1,而乙酸浓度为2.6 kg?m-3.从以上数据可知,酶活性水平和乙酸产量不一致,说明在此C/N比条件下,乙酸主要不是通过丙酮酸途径合成,而可能是通过丙酸、丁酸的转化形成.其他学者的研究都表明,C/N比通过直接或间接影响产能过程及NADH/NAD+的氧(Oxygen)化还原偶联过程,促使不同发酵产酸类型的形成.在低C/N比条件下,乙酸的产生主要是通过氨基酸之间的Stickland反应形成,本文的曲面模型很好地解释了这一结论.
图 3 C/N比、产乙酸关键（解释:比喻事物的重要组成部分)酶和不同酸产量之间的曲面模型所示.由C/N比-关键酶MCM-丙酸曲面模型可以看出,当 C/N比为 55~70时,甲基丙二酰CoA变位酶活性都有大幅度增加,同时丙酸的产量也开始上升,这与任南琪等的研究一致,表明丙酸型发酵中,丙酸产生于糖酵解（glycolysis）丙酮酸途径.根据已有的理论研究,通常情况(Condition)下,丙酸型发酵代谢途径有利于NADH+H+的氧化,而丁酸型发酵缺乏对 NADH 的再生能力.所以当 C/N 值处于56~69时,微生物的细胞合成速率较小,丙酸型发酵比丁酸型发酵有更高的稳定性,产酸结果就形成丙酸型发酵,这也解释了本文曲面模型当中丙酮酸的产量和关键酶活性大幅提升的原因.
3.5 C/N比-关键酶BK-丁酸曲面关系的数学模型的建立
C/N比-关键酶BK-丁酸之间的定量关系如式所示,C/N比-关键酶PTB-丁酸之间的定量关系如式所示.
在初始C/N比为150及更高时,丁酸激酶和磷酸转移丁酰酶的活性均有大幅度增加.以前的研究（research)表明,丁酸的产生主要由糖酵解丙酮（tóng）酸形成.随着C/N值升高到156,碳作为微生物细胞的重要组成元素,促使了细胞合成代谢速率的提高,有机质分解代谢过程中产生的部分NADH可以被厌氧微生物迅速地用于细胞合成而得以再生,所以,呈现较稳定的丁酸型发酵类型.
已有研究表明,通过调控污泥预处理液的初始 C/N 值可实现乙酸、丙酸和丁酸不同厌氧发酵类型之间的转变.在初始C/N比为10~50时,此时发酵是乙酸型发酵.而当 C/N 值为 55~70时,此时厌氧发酵类型则转变为丙酸型.在初始C/N比为150及更高时,此时,丙酸型发酵类型则转变为丁酸型.本文图 3~5建立的曲面模型结果表明,所建立的模型能很好地反映这一趋势,说明模型能较好地模拟C/N比调控下的污泥发酵产酸实验结果.具体参见污水宝商城资料或
图 4 C/N比、丙酸关键酶-甲基丙二酰CoA变位酶和不同酸产量之间的曲面模型
图 5 C/N比、丁酸关键酶活和不同酸产量之间的曲面模型采用多元非线性回归的方法建立了不同C/N比、关键酶活和酸产量之间的函数关系模型.经拟合优度、显著性和残差平方和等指标的比较,发现二元三次多项式模型的拟合优度系数和显著性高、残差平方和较小,表明该模型能较好地拟合三因素之间的定量关系.
2)调节污泥水解液的C/N比分别为12、56、156,将测定的关键酶酶活和乙酸、丙酸、丁酸的结果代入建立的多项式模型,应用Matlab 7.0 软件得到了相应的定量关系式和可视化的曲面模型,曲面模型的结果很好地验证了不同C/N比条件下的产酸代谢途径.
