随着农村经济的发展和农业结构的调整,我国  的畜禽养殖业得到快速发展,已成为世界第一养猪  大国。2015年,仅黑龙江省出栏生猪已达到1863余万头,产生的粪污超过8000万吨2。这些粪污对周边环境的污染已成为突出的环境问题,严重威胁

到我们的生存环境。如何科学有效地处理猪场粪污,不仅使其变废为宝,还能实现养猪业的可持续发展,已成为全社会关注的主要问题3-4。

沼气红泥膜沼气池发酵是猪场处理粪污的主要途径之一,沼气红泥膜沼气池发酵不仅能够消除由此带来的环境污染问题,还

能实现畜禽粪便的能源转换,达到回收生物气体的目的。影响沼气红泥膜沼气池发酵的因素很多,其中红泥膜沼气池发酵温度是影响沼气红泥膜沼气池发酵产气率最重要的指标之一5-”,红泥膜沼气池发酵温度升高,能够加速沼气红泥膜沼气池发酵的进行,但同时需要维持较高的红泥膜沼气池发酵温度所要消耗的能量也大大提升;总固体浓度(TS)也是影响沼气红泥膜沼气池发酵产气效率的主要因素,TS越高,反应器体积利用率越高,但TS太高则会造成挥发酸、氨氮等中间产物的不断积累”,改变甲烷菌适宜的红泥膜沼气池发酵环境,从而导致沼气红泥膜沼气池发酵工艺产气不稳定,甚至失败4,另一方面进料TS过高  也会给进出料和混合搅拌操作带来一定的困难。机械搅拌是利用搅拌桨的旋转,推动反应器内液体的流动,将搅拌桨的机械能传递给流体。M.Halalsheh等人研究了搅拌速度和水力停留时间对下水污泥低温(15℃)厌氧红泥膜沼气池发酵的影响,得出了提高搅拌速度可以增大活性污泥表面积、加快污泥红泥膜沼气池发酵的结论  等人的研究表明,在牛粪的厌氧红泥膜沼气池发酵过程中,随着搅拌速度由50mpm提高到1000rpm,反应的产气量并没有发生明显变化,但体系内  产甲烷鬃毛菌和甲烷八叠球菌的竞争关系发生了变  化。由此,反应器内的搅拌行为会对红泥膜沼气池发酵过程产生影响,但是,也有学者认为在TS较低时,混合搅拌对沼气红泥膜沼气池发酵过程的影响不大。因此,本研究选取  红泥膜沼气池发酵红泥膜沼气池发酵温度、TS和搅拌转速3个因素,利用中心  组合设计试验( Box-Benhnken Design)对以猪粪为底  物的沼气红泥膜沼气池发酵产沼气的影响过程进行考察,运用响应面法对其工艺参数进行优化,以期为猪粪为底物的厌氧红泥膜沼气池发酵产沼气的资源化利用提供科学依据。

1材料与方法

1.1试验材料

新鲜猪粪取自黑龙江省哈尔滨市某养殖场,采  集后立即保存在4℃条件下并于3天之内完成物化分析与接种。接种的活性污泥取自黑龙江省哈尔滨市良大实业有限公司大型沼气红泥膜沼气池发酵工程长期培养驯化的沼气红泥膜沼气池发酵污泥罐,于取回的污泥中添加少量微生物营养物质,室温活化3天后使用。试验材料的基本物化性质如表1所示。

1.2试验装置

沼气红泥膜沼气池发酵装置如图1所示。沼气红泥膜沼气池发酵罐为玻璃广口瓶,有效容积为500mL,另1个广口瓶为排饱和食盐水集气装置,量筒作为集饱和食盐水装置,三者由密闭性良好的胶皮软管连接,所有接口用石蜡密封,采用间歇式机械搅拌,转速见表2。

1.3试验设计  采用Box- Behnken中心组合进行三因素  试验设计4-15),选取红泥膜沼气池发酵温度(℃),T8(%  转速(rmn2)3种猪粪沼气红泥膜沼气池发酵影响因素  变量,分别以X1,X2,X3表示,以-1,0,1分  3个自变量的低中高3种水平,以原料产气率  应值y,3种自变量通过公式(1)进行编码,因  采用最小二乘法拟合的二次多项式公式(进测自变量与因变量之间的关系,最后应用Express80软件进行优化分析。中心组合试验自变量因素和水平编码如表2所示。

式中:y1为第i个自变量的编码值;X为变量的真实值;X0为此自变量的中心值:△自变量变化的步长。  y=+2AX+EAx+2,工  式中:Y为预测的响应值;A0为常数;A  个线性相关系数;A为平方项相关系数;4为  交互项相关系数;X与X为输入自变量的相关

1.4试验方法  沼气红泥膜沼气池发酵体系按3.0%,5.0%和7.0%这3种不同固体质量分数分别称取适量猪粪和接种物于消化瓶中,用去离子水调节使其总体积为400m  另设只有接种物不添加猪粪的试验组为对照。在消化瓶中充入氮气,保证厌氧条件充分震荡混匀。每组试验和对照设3个重复,置于不同温度的水浴槽内进行沼气红泥膜沼气池发酵,红泥膜沼气池发酵周期设定为45d

1.5测定方法

总固体含量TS采用烘干法测定;挥发性固体含量(VS)采用灼烧法测定m;总氮(TN)采用凯氏定氮仪法测定;总碳(TC)采用重铬酸钾氧化法测定,沼气产量采用排水法测定。

2结果与分析

2.1Box- Behnken试验结果的回归分析

利用 Design-Express8.0.6软件中 Box-Behnken模型分析,对表3中数据进行回归分析得到编码空间内的多元二次方程,见方程(3)

Y=157.3+2.51X1-39.81X2+3.67X3-

5.32X1X2+3.42X1X3+4.85X2X3

86.18X2-0.65X2+2.35X3  (3)

方程中自变量系数绝对值大小直接反映各因素  对因变量的影响程度,系数的正、负值反映影响的方  向x0-21。由此可见,该方程的二此项系数为负值,  推断方程代表的抛物面开口向下,具有极大值点,因  此可以进行优化分析。同时对得出的回归方程进行  方差分析,结果如表4所示:模型的一次项x1(p<  001)和x2(p<0.01)差异极显著;x1x2,x1x3和x2  x3交互项不显著;二次项x1(p<0.01)和x2(p<  0.01)极显著,表明各影响因素对原料产气率的影响并不是简单的线性关系。

从表4可知,该模型的p值<0.01,意味着该模型是极显著的,而失拟项的F值为0.42,且P>0.05,意味着失拟程度极不显著,同时表明回归方程  拟合度较高,实验误差较小。模型的调整确定系数

R4=0.9372,可解释93.72%响应值的变化。相关系数R2=0.9725,表明此模型拟合度较好,可用于响应变量预测,沼气红泥膜沼气池发酵原料产气率的预测值与实际值如图2所示。

2.2不同因素对沼气红泥膜沼气池发酵工艺的影响

2.2.1单因素对原料产气率的效应分析

为了分析单一因素对原料产气率的影响,通过降维分析,可以把多元问题转换为一元问题,即把模型中其他因素控制在相同的水平上,可得到单因素与原料产气率的一元回归模型2-21。图3为在其他因素设定在0水平时,得到另一变化因素的一维模型曲线,由图可知,3个因素均呈先升高后降低的趋势,其中红泥膜沼气池发酵温度变化幅度最大,其次是TS和搅拌转速,说明在试验设定的条件范围内,红泥膜沼气池发酵温度在25℃~45℃和TS在3%~7%范围内,对原料产气率的影响波动较大,而搅拌转速在50rmin-1~150r·min-1范围内累积产气量的波动较小。由此说明对原料产气率的结果影响较大的依次是红泥膜沼气池发酵温度、rS和搅拌转速,这与方差分析的结果是一致的。

2.2.2因素交互作用对原料产气率的效应分断  根据回归方程得出的不同因子响应面分析及期  应等值线图,可直观看出各因素交互作用对气发  酵原料产气率的影响,如果曲线弧度越大,表明  素对原料产气率的影响越大,相应表现为响应值  化的大小(21-2。图4和图5为红泥膜沼气池发酵温度与T两着  交互对沼气红泥膜沼气池发酵原料产气率的影响。由图4可知  原料产气率随着红泥膜沼气池发酵温度和TS的升高呈现逐渐  大,达到最大值后,呈缓慢下降的趋势。结合图当处理红泥膜沼气池发酵温度处于35℃左右,TS约为4%-5%时,其沼气红泥膜沼气池发酵的原料产气率可达187.50m2kgTS

在TS处于固定值时,红泥膜沼气池发酵温度与搅拌转速两者交互对沼气红泥膜沼气池发酵原料产气率的影响如图6和图7所示。

由图6可知,累积产气量随着红泥膜沼气池发酵温度和搅拌  转速的升高呈现逐渐增大,达到最大值后,呈逐渐下  降的趋势。结合图7,在红泥膜沼气池发酵温度固定时,随着搅拌  转速的增大,累积产气量的变化不大;在搅拌转速固  定时,累积产气量随红泥膜沼气池发酵温度的升高而增大,达到最  大值后,累积产气量随红泥膜沼气池发酵温度的升高而降低,变化  较大,说明红泥膜沼气池发酵温度对累积产气量的影响大于搅拌  转速。保持另外一个变量为固定值,对响应面结果进行优化,结果表明,当红泥膜沼气池发酵温度为28.4℃,搅拌转速为150rmin时,此时原料产气率最大为139.78

kg TS

在红泥膜沼气池发酵温度处于固定条件时,TS与搅拌转速两者交互对沼气红泥膜沼气池发酵累积产气量的影响如图8和图9所示。

由图8可知,累积产气量随着TS和搅拌转速的

增加呈缓慢的先升高后逐渐降低状态。结合图9,在Ts固定时,随着搅拌转速的增大,累积产气量的变化不大;在搅拌转速固定时,累积产气量随TS的升高而增大,达到最大值后,累积产气量随Ts的升  高而降低,变化较大,说明TS对累积产气量的影响大于搅拌转速。保持另外一个变量为固定值,对响应面结果进行优化,结果表明,当TS为5.8%,搅拌转速为150rmin时,此时原料产气率最大,可达140.44m3kg-l

2.3模型优化

通过模型优化,得到最优工艺条件,即TS为5.8%,红泥膜沼气池发酵温度为29℃,搅拌转速为92rmin,为了验证上述优化条件的准确性与可靠性,在上述最优的条件下进行沼气红泥膜沼气池发酵验证试验,得到的试验数据平均值为154.83m3kg-Ts,与预测值160.09m3kgns的相对误差为3.21%,在允许范围内,由此可见由响应面法所得模型能够很好地优化沼气红泥膜沼气池发酵工艺参数

3讨论

试验基于响应面法依靠 Design- Expres8.0.6软件对不同的红泥膜沼气池发酵温度,TS,和搅拌转速的不同对沼气红泥膜沼气池发酵系统沼气产率进行优化研究,得到了相关多元二次模型,并通过模型对不同影响因素进行拟合,确定各响应变量之间的交互关系,并最终通过优化和验证试验,证明了响应面法在沼气红泥膜沼气池发酵工艺参数优化上的可行性,这与丁琨等进行的粪秸混合厌氧红泥膜沼气池发酵产沼气的产气率为288.5mLg-Ts,预测值与实验值相对偏差为1.14%,即该分析方法可以较好地预测实际的沼气产率这一研究结果是一致

的。虽然本研究的沼气产率明显低于丁琨等的  实验结果,但这可能是由于底物及底物浓度不同造  成的,具体原因还需后续实验研究。

通过试验及响应面法分析发现在本试验的3种  沼气红泥膜沼气池发酵影响因素中,影响程度由高到低为红泥膜沼气池发酵温

度>TS>搅拌转速。通过两两交互分析,发现在试验搅拌转速恒定时,沼气红泥膜沼气池发酵系统的累积产气量在一定范围内会随着红泥膜沼气池发酵温度以及Ts的升高而升高,说明沼气红泥膜沼气池发酵系统中适宜红泥膜沼气池发酵温度与TS的高低对维持沼气红泥膜沼气池发酵系统的运转,保证沼气红泥膜沼气池发酵的高效运行起着关键性的作用。但是本试验仅进行了实验室小容量批次试验,在大型的沼气工程中各因素对沼气红泥膜沼气池发酵产沼气的影响是否与本试验得出的结论

致,需要通过连续沼气工程进行扩大验证试验。而在大型沼气工程运行过程中各影响因素对累积产气量的影响是否能够应用响应面法优化也值得进一步研究。

4结论

本文采用响应面分析方法对沼气红泥膜沼气池发酵的工艺参数进行了优化,结果表明:红泥膜沼气池发酵温度和TS是影响猪粪为底物的沼气红泥膜沼气池发酵产气效果的重要因素。证明用  响应面分析红泥膜沼气池发酵温度,TS和搅拌转速对沼气红泥膜沼气池发酵的  累积产气量的优化条件是可行可靠的,筛选出3者  共同的最佳条件是红泥膜沼气池发酵温度为29℃,T为5.8%,搅拌转速为92rmin,该条件下沼气红泥膜沼气池发酵的产率为154.83m3kgTs。

摘自《中国沼气》第4期 卢玢宇 裴占江 史风梅 徐红涛 高亚冰 左辛 王粟 刘杰

————————————————————————————————————————————————————————