一种分析餐厨垃圾厌氧消化产甲烷能力的方法与流程

本发明涉及一种分析餐厨垃圾厌氧消化产甲烷能力的方法，属于有机固体废弃物处理技术领域。

背景技术：

餐厨垃圾的主要成分是糖类、蛋白质类、脂肪类和纤维素类物质，因为地域、季节和饮食习惯的不同，餐厨垃圾中各组分的含量可能有很大差别，其产甲烷性能也存在差异。另外，各组分的产甲烷特性相对稳定，可以通过测定某一实际餐厨垃圾中各组分的含量来模拟和预测其产甲烷特性。尽管有些文献报道了关于餐厨垃圾及其组分的产气研究，但是鲜有关于各组分对餐厨垃圾产气的影响以及相关模型的研究报道。

为获得有机物的产甲烷潜力，传统方法是进行产甲烷潜力实验，但是该实验耗时较长，通常需要30-60天。而采用数学模型可以快速获得产甲烷潜力的结果，比如buswell公式，它根据底物的化学特性计算理论的产甲烷量，然而该类模型不能提供有关底物降解的动力学信息。

技术实现要素：

为了解决上述问题，本发明提供了一种评估餐厨垃圾厌氧消化产甲烷的动力学模型，包括(a)，(b)，(c)，(d)，(e)；

(a)，y(t)＝8.66+1.63×b1×p1(t)+1.01×b2×p2(t)+1.16×b3×p3(t)+1.3×b4×p4(t)；

(b)p1(t)＝496.94×(1-e^-0.147×t)；

(c)p2(t)＝422.06×(1-e^-0.122×t)；

(d)p3(t)＝431.25×(1-e^-0.178×t)；

(e)

p1(t)、p2(t)、p3(t)、p4(t)分别表示淀粉、纤维素、蛋白质和脂肪在t时刻的产气量；y(t)表示t时刻餐厨垃圾的产甲烷量；bi(i＝1,2,3,4)分别表示餐厨垃圾中淀粉、蛋白质、纤维素和脂肪的百分含量；e为自然对数；t为消化反应时间。

本发明的第二个目的是提供一种评估餐厨垃圾产甲烷能力的方法，所述方法是应用所述动力学模型，根据餐厨垃圾中淀粉、蛋白质、脂肪和纤维素的百分含量，代入所述动力学模型中，计算餐厨垃圾产甲烷量随时间的变化以及最大的甲烷产量。

在本发明的一种实施方式中，所述餐厨垃圾ts为20～30％，vs为20～25％，淀粉含量为20～40％，蛋白质含量为15～25％，脂肪含量为15～25％，纤维素含量为15～30％。

在本发明的一种实施方式中，所述餐厨垃圾ts为24.13％，vs为22.60％，淀粉含量为31.87g/gts，蛋白质含量为21.02g/gts，脂肪含量为17.56g/gts，纤维素含量为23.21g/gts。

在本发明的一种实施方式中，所述餐厨垃圾ts为25.95g/gts，vs为24.46g/gts，淀粉含量为35.61g/gts，蛋白质含量为18.72g/gts，脂肪含量为19.11g/gts，纤维素含量为20.82g/gts。

本发明的第三个目的是提供一种产甲烷的方法，所述方法是应用所述动力学模型，根据餐厨垃圾中淀粉、蛋白质、脂肪和纤维素的百分含量，代入所述动力学模型中，计算餐厨垃圾产甲烷量为最大甲烷产量80～100％所需的发酵时间，并在35～37℃下发酵相应时间。

在本发明的一种实施方式中，所述方法接种污泥进行厌氧发酵；所述污泥的性质为ts13～18％，vs11～16％。

在本发明的一种实施方式中，所述方法接种污泥进行厌氧发酵；所述污泥的性质为ts16.21％，vs14.22％。

在一种实施方式中，所述厌氧消化是在37℃下进行。

在一种实施方式中，所述接种物/底物的ts比为2.66，含固率为8.06％，初始ph为8.87。

在一种实施方式中，所述调节含固率是用水调节；调整ph是用naoh溶液和/或hcl溶液调节。

在一种实施方式中，所述发酵在产甲烷潜力测试系统(amptsⅱ)中进行。

在一种实施方式中，所述底物量为8gts，所述底物为米饭、肉和蔬菜。

本发明还提供所述动力学模型在降解餐厨垃圾方面的应用。

有益效果：本发明提供的评估餐厨垃圾厌氧消化产甲烷的动力学模型能够根据餐厨垃圾的淀粉、蛋白质、脂肪、纤维素含量评估其产甲烷的过程及最大产甲烷量，以动力学模型评估的餐厨垃圾甲烷产量理论值与实际产量的相对误差可达0.14％，评估结果准确。

附图说明

图1为淀粉、蛋白质、脂肪、纤维素等特征组分的产气及拟合情况；

图2为餐厨垃圾各组分反应速率情况；

图3为餐厨垃圾各组分含量及产气情况；

图4为模拟餐厨垃圾产气及拟合情况；

图5为实际餐厨垃圾产气及拟合情况。

具体实施方式

实验装置采用全自动甲烷潜力分析系统(ampts)，接种ts16.21％，vs14.22％的污泥进行厌氧发酵；接种物与底物的ts比为2.66，用水调节含固率为8±1％，调节初始ph为8。0～9.0；反应温度为37℃，气体体积由amptsv5.0软件统计。淀粉、蛋白质、脂类、纤维素含量分别采用gb5009.9-2016，gb50095-2010，gb5009.6-85,gbt5009.10-2003文件中的方法进行测定。

表1为不同组分底物的组分及组成情况，共分为两组：一组为特征物质实验，分别以米饭(r1)、豆腐(r2)、肥肉(r3)和青菜(r4)为底物进行厌氧消化反应；另一组以餐厨垃圾为特征物质，分别以两种来源的实际餐厨垃圾(r5、r6)和三种模拟餐厨垃圾(r7、r8、r9)为底物进行厌氧消化反应。其中，餐厨垃圾r5的淀粉含量为31.87g/gts，蛋白质含量为21.02g/gts，脂肪含量为17.56g/gts，纤维素含量为23.21g/gts；餐厨垃圾r6的淀粉含量为35.61g/gts，蛋白质含量为18.72g/gts，脂肪含量为19.11g/gts，纤维素含量为20.82g/gts。

表1不同组分底物的组分及组成情况

用修正gompertz模型或一级动力学模型对特征物质原料的累积产气情况的拟合结果。

p(t)＝p×(1-exp(-k×t))(1)

以上式子中p(t)-t时刻累积甲烷产量，ml/gts；p-最大甲烷产量，ml/gts；rm-最大产甲烷速率，ml/gts/d；e-自然对数，为2.718；λ-延滞期；t-消化反应时间；k-一级底物降解速率。

此外，对甲烷产量而言，用一级动力学常数描述反应速率。

实施例1

表2、表3分别为修正gompertz模型和一级动力学模型对不同特征物质原料的累积产气情况的拟合结果。

表2修正gompertz模型参数

表3一级动力学模型参数

由修正gompertz模型和一级动力学模型拟合，淀粉、纤维素、蛋白质和脂肪的累积产气情况可分别用如下式子表示：

p1(t)＝496.94×(1-e^-0.147×t)(3)；

p2(t)＝422.06×(1-e^-0.122×t)(4)；

p3(t)＝431.25×(1-e^-0.178×t)(5)；

式中p1(t)、p2(t)、p3(t)、p4(t)分别表示淀粉、纤维素、蛋白质和脂肪在t时刻的产气量。

图1为特征物质累积产气及对应模型的拟合情况，各特征物质用相应的模型拟合效果很好，所有的r²均在0.99以上。厌氧反应22天后，淀粉、纤维素、蛋白质和脂肪的实际累积产量分别是467.44ml/gts、383.91ml/gts、424.53ml/gts、334.57ml/gts，采用式3～6预测甲烷量分别为477.36ml/gts、392.60ml/gts、422.66ml/gts、338.52ml/gts，相对误差分别为2.12％、2.26％、0.44％、1.18％，预测最大的甲烷潜力分别为496.94ml/gts、422.06ml/gts、431.25ml/gts、410.28ml/gts，生物降解能力分别为94.06％、90.97％、98.44％、81.55％。

实施例2

图2为各组原料的水解速率k值和最终的累积甲烷产量。由图2可知，模拟餐厨垃圾的累积甲烷产量明显高于各特征物质的累积甲烷产量。r8的累积产甲烷量为535.69ml/gts，比r7和r9的累积产甲烷量分别提高5.21％和2.51％，同时，与r7和r9相比，r8具有更高的k值，为0.129。餐厨垃圾的累积产甲烷能力在500-540ml/gts范围内，而特征物质的累积产甲烷量在330-470ml/gts范围内，餐厨垃圾的产甲烷量明显高于特征物质的产甲烷量。由图2还可以发现，虽然餐厨垃圾的产甲烷量比特征物质高，但是，动力学常数并不都高于特征物质，有些甚至低于特征物质，模拟餐厨垃圾的动力学常数在0.1～0.13之间，低于米饭和豆腐的k值，而实际餐厨垃圾的k值为0.168和0.189。实际餐厨垃圾的产甲烷量和模拟餐厨垃圾的产甲烷量差别不大。

实施例3

将餐厨垃圾的组分分为碳水化合物、蛋白质和脂肪三类，随着碳水化合物的含量增加，餐厨垃圾的产甲烷量增加，但当碳水化合物的含量达到70％时，随着碳水化合物含量的增加，餐厨垃圾的产甲烷量反而降低。当脂肪的含量低于30％时，餐厨垃圾的产甲烷量与脂肪含量呈正比，而超过30％时，餐厨垃圾的产甲烷量与脂肪含量呈反比。碳水化合物的含量为50％-70％，蛋白质的含量为25％-50％，脂肪在20％-30％时，餐厨垃圾的产甲烷量较高。

实施例4

将3组自制餐厨垃圾t时刻的产甲烷数据与对应餐厨垃圾中各组分理论产甲烷量数据，共3×23＝69组数据进行主效应相关性分析得到，米饭、豆腐、肥肉、青菜的理论产甲烷量与餐厨垃圾的产甲烷量的pearson相关性值分别为0.925、0.922、0.846、0.902，说明它们与餐厨垃圾的产甲烷量呈强相关性。因此，使用方程式(7)进行多元线性回归方程拟合分析，得到相应参数如表4所示。

y(t)＝a0+a1×b1×p1(t)+a2×b2×p2(t)+a3×b3×p3(t)+a4×b4×p4(t)(7)

其中，y(t)-t时刻餐厨垃圾的产甲烷量；ai-多元线性回归系数；bi-餐厨垃圾中淀粉、蛋白质、纤维素和脂肪的百分含量；pi(t)-t时刻淀粉、蛋白质、纤维素和脂肪的产甲烷量；

表4多元线性回归系数及拟合参数

由表4可知，模型拟合参数r²为0.99≈1，p值为0.0000<0.01，说明根据餐厨垃圾中淀粉、蛋白质、纤维素和脂肪的含量能够很好的预测餐厨垃圾的产气情况，预测表达式如方程式(8)所示。

y(t)＝8.66+1.63×b1×p1(t)+1.01×b2×p2(t)+1.16×b3×p3(t)+1.3×b4×p4(t)(8)

图4为模拟餐厨垃圾和实际餐厨垃圾的累积产气及用方程式(8)的拟合情况。由图4可知，方程式(8)可以很好地描述模拟餐厨垃圾的产甲烷过程。

采用方程式(8)对实际餐厨垃圾的产气情况进行拟合，与实际产气情况进行比对，结果如图5所示，方程式(8)可以很好地描述实际餐厨垃圾的产甲烷过程，两种实际餐厨垃圾的拟合r²分别为0.950和0.951，22天后餐厨垃圾a和餐厨垃圾b的累积产甲烷量分别为527.47ml/gts和522.1ml/gts，方程(8)拟合值分别为528.22ml/gts和545.29ml/gts，相对误差分别为0.14％和4.44％。

虽然本发明已以较佳实施例公开如上，但其并非用以限定本发明，任何熟悉此技术的人，在不脱离本发明的精神和范围内，都可做各种的改动与修饰，因此本发明的保护范围应该以权利要求书所界定的为准。