廚余垃圾厭氧消化產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ偷男拚芯?/h1>
2011-12-20 09:10:56張代鈞重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心重慶400067重慶大學(xué)資源與環(huán)境科學(xué)學(xué)院重慶400030
中國環(huán)境科學(xué) 2011年5期
吳 云,張代鈞 (1.重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心,重慶 400067；.重慶大學(xué)資源與環(huán)境科學(xué)學(xué)院,重慶 400030)
廚余垃圾厭氧消化產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ偷男拚芯?/p>
吳 云1*,張代鈞2(1.重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心,重慶 400067；2.重慶大學(xué)資源與環(huán)境科學(xué)學(xué)院,重慶 400030)
分別通過單相和兩相厭氧消化實(shí)驗(yàn),對廚余垃圾產(chǎn)甲烷累積產(chǎn)量曲線進(jìn)行了擬合.結(jié)果表明,單相厭氧消化產(chǎn)甲烷累積產(chǎn)量曲線符合二次函數(shù)關(guān)系,相關(guān)系數(shù)(R2)為 0.9982;兩相厭氧消化產(chǎn)甲烷累積產(chǎn)量曲線符合線性函數(shù)關(guān)系,相關(guān)系數(shù)(R2)為 0.9972.以擬合函數(shù)為基礎(chǔ),在考慮系統(tǒng)負(fù)荷的前提下,建立了基準(zhǔn)狀態(tài)下的產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ?通過引入 pH抑制系數(shù)和氨氮抑制系數(shù)對建立的產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行了修正.對經(jīng)驗(yàn)修正模型預(yù)測結(jié)果進(jìn)行了實(shí)驗(yàn)驗(yàn)證及誤差分析.結(jié)果表明, 在檢驗(yàn)水平a=0.01條件下,單相厭氧第1組和第2組試驗(yàn)數(shù)據(jù)試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為486.22和268.63,均遠(yuǎn)大于臨界值(F0.99(1,58)=7.12);兩相厭氧第1組試驗(yàn)數(shù)據(jù)和第2組試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為421.69和372.56,也遠(yuǎn)高于其臨界值(F0.99(1,58)=7.12).說明該模型能很好地預(yù)測單相和兩相厭氧消化產(chǎn)甲烷累積產(chǎn)量,且計(jì)算過程簡單.
廚余垃圾；厭氧消化；產(chǎn)甲烷；經(jīng)驗(yàn)?zāi)Ｐ?/p>
對于以產(chǎn)甲烷為主要目的的廚余垃圾厭氧消化處理系統(tǒng)而言,甲烷產(chǎn)量決定著系統(tǒng)運(yùn)行的成敗.因此,長期以來,各國學(xué)者對厭氧消化產(chǎn)甲烷數(shù)學(xué)模型的研究從未停止[1-3].2002年,國際水協(xié)(IWA)正式公布了厭氧消化 1號數(shù)學(xué)模型(ADM1),該模型是集大成的研究成果,較為完整和準(zhǔn)確地描述了復(fù)雜厭氧消化系統(tǒng)的各個子過程,為厭氧工藝的建模和模擬提供了一個通用的平臺[4].ADM1模型推出后,一些關(guān)于模型補(bǔ)充和修正的研究工作陸續(xù)展開,如 Fedorovich等[5]提出了包含硫酸鹽還原的擴(kuò)展模型;Batstone[6]提出了利用氫氣的氧化描述硫酸鹽還原過程的簡便方法; Parker等[7]以ADM1為基礎(chǔ),提出了包含產(chǎn)氣及其釋放的擴(kuò)展模型,從而將產(chǎn)氣模型的機(jī)理性描述提高到一個新的水平.與此同時,一些新的建模思想和方法也在陸續(xù)出現(xiàn),如王凱軍等[8]采用系統(tǒng)動力學(xué)理論構(gòu)建了厭氧消化數(shù)學(xué)模型(SD模型),將變量間的各種復(fù)雜關(guān)系轉(zhuǎn)換為輸入輸出關(guān)系,較為簡便地實(shí)現(xiàn)了不同模擬方法之間的有效契合.王洪濤等[9]也以厭氧消化過程中碳素轉(zhuǎn)移機(jī)理為基本構(gòu)建思想,提出了基于碳素轉(zhuǎn)移機(jī)理的厭氧消化產(chǎn)氣模型.
但是,各種描述產(chǎn)甲烷過程的模型,都存在一個共同的問題,其計(jì)算工作十分繁雜,這給模型的實(shí)際應(yīng)用造成了不便.而在某些實(shí)際工作中,人們并不關(guān)心模型對厭氧消化過程細(xì)節(jié)的描述,只看重宏觀上的運(yùn)行結(jié)果,以便為可行性分析和投資決策提供相對準(zhǔn)確的依據(jù).通常,為了簡便快速地預(yù)測運(yùn)行結(jié)果,常常通過有機(jī)物總量或產(chǎn)甲烷潛能來粗略估算甲烷產(chǎn)量和產(chǎn)生速率,而這種估算的結(jié)果與實(shí)際運(yùn)行結(jié)果往往存在很大差異,導(dǎo)致投資決策出現(xiàn)偏差.本研究針對這些問題,在黑箱理論思想的指導(dǎo)下,通過實(shí)驗(yàn)對單相和兩相厭氧消化產(chǎn)甲烷規(guī)律進(jìn)行考察和分析,采用函數(shù)擬合的方法,構(gòu)建介于準(zhǔn)確描述和粗略估計(jì)之間的簡化經(jīng)驗(yàn)?zāi)Ｐ?以達(dá)到簡便預(yù)測甲烷產(chǎn)量和產(chǎn)生速率的目的.
1 材料和方法
1.1 材料
廚余垃圾樣品采自中餐廳、西餐廳、火鍋店和學(xué)生食堂.樣品采回后,剔除其中的泡沫飯盒、玻璃、石塊和金屬等雜物,并按重慶市主城區(qū)廚余垃圾基礎(chǔ)數(shù)據(jù)調(diào)研分析報(bào)告中提供調(diào)查數(shù)據(jù)進(jìn)行配比,即比例為:中餐類/湯鍋火鍋類/西餐類為0.65/0.30/0.05,將樣品用破碎機(jī)破碎后混勻,采用4分法取出,放入4℃冰箱中保存待用.廚余垃圾的基本物料性質(zhì)見表 1.實(shí)驗(yàn)用原生污泥取自重慶市唐家橋污水廠厭氧消化池,污泥含水率為94.23%,MLSS為47.36g/L,VSS為32.48g/L.經(jīng)馴化培養(yǎng)后作為接種污泥.
表1 廚余垃圾性質(zhì)(%)Table 1 Characteristics of kitchen wastes(%)
1.2 試驗(yàn)方法
單相厭氧消化試驗(yàn)采用試驗(yàn)裝置(圖1),試驗(yàn)前將物料破碎為粒徑小于10mm的顆粒.根據(jù)文獻(xiàn)中關(guān)于接種污泥混合比的報(bào)道[10-12],污泥和原料的配比從1:10到1:1不等,當(dāng)接種污泥比例低時,產(chǎn)氣滯后期較長,實(shí)驗(yàn)進(jìn)行時間較長,接種污泥比例高時,產(chǎn)氣波動較大,誤差也較大.而根據(jù)課題組前期的實(shí)驗(yàn)結(jié)果[13],采用2:1的配比時,消化情況較好,產(chǎn)氣情況也較穩(wěn)定.因此,選擇2:1的比例,即將1.5kg廚余垃圾和0.75kg活性污泥填入反應(yīng)器中,然后加水至 3L左右.試驗(yàn)過程中調(diào)節(jié)pH為7左右,控制溫度為35℃.試驗(yàn)進(jìn)行8d后,開始每天投加0.3kg新鮮垃圾.新鮮垃圾與當(dāng)天收集的50mL滲濾液混合,由反應(yīng)器頂部加入,并加水保持反應(yīng)物體積在 3L左右.穩(wěn)定運(yùn)行后每日由底部出料口出料1次,出料量為0.2～0.4kg.整個試驗(yàn)共持續(xù) 60d.試驗(yàn)過程中采用排飽和氯化鈉溶液法收集氣體,并跟蹤測定甲烷產(chǎn)量.
兩相厭氧消化采用圖 2所示試驗(yàn)裝置,試驗(yàn)開始時,將 2.86kg經(jīng)破碎處理的廚余垃圾填入水解產(chǎn)酸相反應(yīng)器,并加水 5.4L.活性污泥與廚余垃圾的接種比例為 1:2[10-12].試驗(yàn)過程中調(diào)節(jié)系統(tǒng)pH值為7左右,并控制溫度為35℃.系統(tǒng)啟動運(yùn)行8d后,每天將0.5kg新鮮垃圾與UASB出水混合后添加到產(chǎn)酸相反應(yīng)器,水解后的出水進(jìn)入水解液平衡池,并根據(jù)產(chǎn)酸相情況確定每天的出料量.
UASB反應(yīng)器內(nèi)接種污泥量為2.6L,為了防止水解液 pH值過低影響甲烷相運(yùn)行,進(jìn)水前調(diào)節(jié)水解液pH7左右,并由計(jì)量泵定時將水解液送入U(xiǎn)ASB反應(yīng)器.為了維持UASB反應(yīng)器內(nèi)一定的上升流速和溫度,設(shè)置了外循環(huán)水池,通過外循環(huán)加熱將反應(yīng)器水溫控制在 33～37℃.進(jìn)水泵間歇運(yùn)行,每小時運(yùn)行10min,由定時器控制,進(jìn)水泵流量為10～12mL/min.外循環(huán)泵24h運(yùn)行,循環(huán)流量為 12～24mL/min.試驗(yàn)過程中采用排飽和氯化鈉溶液法收集氣體,并跟蹤測定甲烷產(chǎn)量.
圖1 單相厭氧消化試驗(yàn)裝置示意Fig.1 Sketch of experiment equipment of single phrase anaerobic process
圖2 兩相厭氧消化試驗(yàn)裝置示意Fig.2 Sketch of experiment equipment of two phrase anaerobic process
試驗(yàn)采用GASBOARD紅外氣體分析儀(武漢光電科技)來測定氣體成分及含量(體積百分?jǐn)?shù)).甲烷氣體含量測定結(jié)果確定后,將收集氣體體積乘以氣體分析儀測定得到的體積百分?jǐn)?shù),即可計(jì)算出甲烷產(chǎn)量;氨氮濃度的測定采用納氏試劑比色法進(jìn)行測定(上海精科的 721分光光度計(jì));pH值采用PHS-3C型pH計(jì)進(jìn)行測定.
2 結(jié)果與分析
由圖3可見,單相厭氧消化甲烷累積產(chǎn)量曲線趨勢與二次函數(shù)曲線非常相似,兩相厭氧消化甲烷累積產(chǎn)量曲線趨勢則基本與線性函數(shù)曲線相似.因此,將3次實(shí)驗(yàn)結(jié)果取平均值,并分別采用二次函數(shù)和線性函數(shù)對甲烷累積產(chǎn)量的平均值曲線進(jìn)行擬合.
圖3 單相及兩相厭氧消化甲烷累積產(chǎn)量平行實(shí)驗(yàn)Fig.3 Cumulative CH4 curve of single and two phrase digestion repeated experiment
從圖3和圖4中可以看到,二次函數(shù)和線性函數(shù)的擬合效果很好,R2分別達(dá)到了 0.9982和0.9972.經(jīng)分析后認(rèn)為,在適宜的溫度(35℃)、pH6～7和有機(jī)負(fù)荷條件下,單相和兩相厭氧消化甲烷累積產(chǎn)量曲線的這種變化趨勢與其各自的反應(yīng)器運(yùn)行特性是基本吻合的.因?yàn)?對于單相厭氧消化而言,系統(tǒng)的環(huán)境條件會隨著有機(jī)物水解、產(chǎn)酸、產(chǎn)乙酸和產(chǎn)甲烷過程發(fā)生變化,由于單相厭氧未將產(chǎn)甲烷過程分置進(jìn)行,因此系統(tǒng)的這些過程直接影響著甲烷菌產(chǎn)甲烷速率的高低,且產(chǎn)甲烷速率整體上表現(xiàn)為前期低、后期高的趨勢,這就導(dǎo)致其甲烷累積曲線表現(xiàn)為二次曲線形態(tài).而兩相厭氧消化系統(tǒng)將產(chǎn)甲烷過程分置到甲烷相中進(jìn)行,甲烷菌所受系統(tǒng)環(huán)境條件的影響較小,基本能夠保持穩(wěn)定產(chǎn)氣狀態(tài),其甲烷累積曲線也就基本表現(xiàn)為線性曲線趨勢[13].
圖4 單相及兩相厭氧消化甲烷累積產(chǎn)量擬合曲線Fig.4 Cumulative CH4 AV curve and fitting curve of single and two phrase digestion
3 簡化模型的建立與修正
3.1 基準(zhǔn)模型的建立
根據(jù)上述實(shí)驗(yàn)條件,由于物料供給充足;對系統(tǒng)溫度、pH值、氨氮濃度等環(huán)境條件進(jìn)行了嚴(yán)格監(jiān)測和控制,使其保持在適宜的消化狀態(tài).因此,在建立基準(zhǔn)模型前,可作如下假設(shè):①假設(shè)厭氧消化過程中基質(zhì)濃度充足,厭氧菌群生長代謝不受基質(zhì)濃度限制;②假設(shè)厭氧消化進(jìn)行過程中,系統(tǒng)溫度、pH值、氨氮濃度等環(huán)境因素處于適宜狀態(tài),不對厭氧消化過程造成抑制影響;③假設(shè)系統(tǒng)未超負(fù)荷運(yùn)行時,累積產(chǎn)氣量與厭氧消化的有機(jī)物負(fù)荷量呈正比關(guān)系.
從圖3中可看到,不同試驗(yàn)條件下的甲烷累積曲線與簡單函數(shù)曲線相類似,圖3a的甲烷累積曲線與二次函數(shù)曲線線類似,圖3b的甲烷累積曲線則與線性函數(shù)曲線類似,擬合結(jié)果表明,實(shí)驗(yàn)結(jié)果與擬合曲線相關(guān)系數(shù)達(dá)到0.95以上,證實(shí)擬合曲線方程能夠很好地預(yù)測甲烷累積產(chǎn)量.因此,可將甲烷累積產(chǎn)量擬合曲線方程作為簡化的甲烷累積產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ头匠?經(jīng)驗(yàn)?zāi)Ｐ头匠痰谋磉_(dá)式如式(1)和式(2):
單相累積甲烷產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ头匠?
式中: t為厭氧消化時間; G為甲烷累積產(chǎn)量,L.
再將上述方程對時間求導(dǎo)即可得到簡化的產(chǎn)氣速率模型方程.對式(1)和式(2)求導(dǎo)得:
單相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
兩相累積甲烷產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ头匠?
需要說明的是,從式(4)可看出,在本實(shí)驗(yàn)條件下,建立的兩相厭氧消化的產(chǎn)氣速率方程接近常數(shù),這與劉國濤[14]和Sun等[15]的觀察結(jié)果基本一致.
根據(jù)假設(shè),在未超負(fù)荷情況下,甲烷產(chǎn)量與有機(jī)物負(fù)荷量呈正比.因此,在考慮系統(tǒng)有機(jī)負(fù)荷的情況下,可引入負(fù)荷比A/A0,A為實(shí)際有機(jī)負(fù)荷,A0為基準(zhǔn)有機(jī)負(fù)荷,此處的基準(zhǔn)有機(jī)負(fù)荷即為前述實(shí)驗(yàn)中的系統(tǒng)每日有機(jī)負(fù)荷量.則式(3)和式(4)可變?yōu)?
單相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
兩相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
兩相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
因此根據(jù)模型假設(shè)和實(shí)驗(yàn)條件,可以將式(5)和式(6)稱為基準(zhǔn)狀態(tài)下的產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ?
3.2 基準(zhǔn)模型的修正
由于式(5)和式(6)是在假設(shè)條件和特定實(shí)驗(yàn)條件下建立的基準(zhǔn)狀態(tài)產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?而在實(shí)際情況下,還應(yīng)考慮系統(tǒng) pH值和氨氮抑制性對產(chǎn)氣速率的影響,因此有必要對式(5)和式(6)進(jìn)行修正.根據(jù)ADM1模型中提出的關(guān)于pH值抑制系數(shù)和氨氮抑制系數(shù)的相關(guān)理論[4,16],pH值抑制系數(shù)和氨氮抑制系數(shù)可分別表達(dá)如下:
IpH=1,pH=7(中性條件產(chǎn)甲烷系統(tǒng)不受酸堿抑制)
式中:pH代表系統(tǒng)實(shí)際pH值;pHUL和pHLL是生物種群受到50%抑制時的pH值上限和下限,一般分別取7和6[4]; S1是氨氮濃度;K1是氨氮抑制參數(shù).因此,在實(shí)際情況下,可引入上述抑制系數(shù)對基準(zhǔn)經(jīng)驗(yàn)?zāi)Ｐ褪?5)和式(6)進(jìn)行修正.修正后的經(jīng)驗(yàn)?zāi)Ｐ捅磉_(dá)式如下:
單相產(chǎn)氣速率經(jīng)驗(yàn)修正模型:
兩相產(chǎn)氣速率經(jīng)驗(yàn)修正模型:
式中: IpH為 pH抑制系數(shù);INH3為氨氮抑制系數(shù),上述經(jīng)驗(yàn)?zāi)Ｐ偷膯挝粸長/d. 根據(jù)建立基準(zhǔn)經(jīng)驗(yàn)產(chǎn)氣速率模型的試驗(yàn)條件,可令單相厭氧的A0=0.3kg/d,兩相厭氧的A0=0.5kg/d.
4 經(jīng)驗(yàn)?zāi)Ｐ偷尿?yàn)證
4.1 試驗(yàn)曲線與預(yù)測曲線的比較
單相厭氧消化驗(yàn)證試驗(yàn)所用物料性質(zhì)和試驗(yàn)條件與1.1和1.2所述實(shí)驗(yàn)條件基本一致,但將有機(jī)負(fù)荷分別控制為 0.25,0.36kg,并在試驗(yàn)過程中控制系統(tǒng)pH值約為5.5左右(改變系統(tǒng)pH條件);因此當(dāng)系統(tǒng)pH值為5.5時,其系統(tǒng)氨氮平均濃度為 0.826g/L左右,另根據(jù)文獻(xiàn)報(bào)道,K1可取45.063[4],因此,經(jīng)抑制系數(shù)公式計(jì)算,IpH和INH3分別為0.389和0.982.
兩相厭氧消化驗(yàn)證試驗(yàn)所用物料性質(zhì)和試驗(yàn)條件與1.1和1.2所述實(shí)驗(yàn)條件基本一致,僅將有機(jī)負(fù)荷分別變?yōu)?0.3,0.75kg.在上述實(shí)驗(yàn)條件下,系統(tǒng) pH值為 7,甲烷相氨氮平均濃度為1.446g/L,經(jīng)抑制系數(shù)公式計(jì)算,得IpH和INH3分別為1和0.967.模型驗(yàn)證結(jié)果如圖7和圖8所示.
從圖5中可以看出,甲烷累積產(chǎn)量模型預(yù)測曲線與試驗(yàn)曲線基本吻合,說明簡化的產(chǎn)氣模型能夠用于預(yù)測單相和兩相厭氧消化甲烷產(chǎn)氣量.
圖5 單相及兩相厭氧消化甲烷累積產(chǎn)量和預(yù)測曲線Fig.5 cumulative CH4 Curve and predict curve of single and two phrase digestion
4.2 誤差分析
將不同試驗(yàn)條件下得到的單相厭氧消化和兩相厭氧消化甲烷累積產(chǎn)量試驗(yàn)數(shù)據(jù)與模型預(yù)測值進(jìn)行方差檢驗(yàn).檢驗(yàn)結(jié)果表明,單相厭氧第 1組和第 2組試驗(yàn)數(shù)據(jù)試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為486.22和268.63,在檢驗(yàn)水平α=0.01下均遠(yuǎn)大于臨界值[F0.99(1,58)=7.12];兩相厭氧第1組試驗(yàn)數(shù)據(jù)和第2組試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為421.69和372.56,也遠(yuǎn)高于其臨界值[F0.99(1,58)=7.12].說明簡化的厭氧消化產(chǎn)氣修正經(jīng)驗(yàn)?zāi)Ｐ偷幕貧w是顯著的,因此該模型能夠較好地預(yù)測厭氧消化甲烷累積產(chǎn)量.
5 結(jié)論
5.1 在一定的實(shí)驗(yàn)條件下,單相厭氧消化甲烷累積產(chǎn)量曲線與二次函數(shù)曲線相似,兩相厭氧消化甲烷累積產(chǎn)量曲線與線性函數(shù)曲線相似.擬合結(jié)果表明,二次函數(shù)和線性函數(shù)曲線能夠很好地對單相和兩相厭氧消化甲烷累積產(chǎn)量曲線進(jìn)行擬合,R2分別達(dá)到了0.9982和0.9972.
5.2 在模型假設(shè)基礎(chǔ)上,建立了基準(zhǔn)狀態(tài)下的甲烷累積產(chǎn)量簡化經(jīng)驗(yàn)?zāi)Ｐ秃彤a(chǎn)甲烷速率簡化經(jīng)驗(yàn)?zāi)Ｐ?通過引入 pH值抑制系數(shù)和氨氮抑制系數(shù)對模型進(jìn)行修正,建立了修正后的甲烷累積產(chǎn)量簡化經(jīng)驗(yàn)?zāi)Ｐ秃彤a(chǎn)甲烷速率簡化經(jīng)驗(yàn)?zāi)Ｐ?
5.3 模型驗(yàn)證結(jié)果表明,簡化的甲烷累積產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ皖A(yù)測曲線與試驗(yàn)曲線基本吻合,簡化經(jīng)驗(yàn)?zāi)Ｐ湍軌蛴糜陬A(yù)測單相和兩相厭氧消化甲烷累積產(chǎn)量;誤差分析結(jié)果表明,對簡化經(jīng)驗(yàn)?zāi)Ｐ颓€的方差檢驗(yàn)結(jié)果是顯著的,同樣也說明,簡化經(jīng)驗(yàn)?zāi)Ｐ湍軌蛴糜陬A(yù)測甲烷累積產(chǎn)量.
[1] Christensen T H.Landfill of waste:biogas [M]. London: E and FNSpon, 1996,56-74.
[2] Schumacher M M.Landfill methane recovery [M]. USA: NewJersey,1983,35-46.
[3] El-fadel M, Findikakis A N. A numerical model for methane production in managed sanitary landfills [J]. Waste Management and Research, 1989,7(1):31-42.
[4] 國際水協(xié)厭氧消化工藝數(shù)學(xué)模型課題組.厭氧消化數(shù)學(xué)模型[M]. 張亞雷,周雪飛譯.上海:同濟(jì)大學(xué)出版社, 2004:15-30.
[5] Fedorovich V, Lens P, KalyuzhnyiI S, et al. Extension of anaerobic digestion model No.1 with processes of sulfate reduction [J]. Applied Biochemistry and Biotechnology, 2003, 109(1-3):33-45.
[6] Batstone D J. Mathematical modelling of anaerobic reactors treating domestic wastewater: rational criteria for model use[J]. Reviews in Environmental Science and Bio-Technology, 2006,5(1):57-71.
[7] Parker W J, WU G H. Modifying ADM1 to include formation and emission of odourants [J]. Water Science and Technology, 2006,54(4):111-117.
[8] 閻 中,王凱軍.厭氧消化過程系統(tǒng)動力學(xué)模型構(gòu)建方法研究[J]. 中國沼氣, 2009,27(2):3-8.
[9] 楊渤京,王洪濤,陸文靜,等.厭氧條件下纖維素降解產(chǎn)氣的研究[J]. 中國環(huán)境科學(xué), 2008,28(3):255-259.
[10] 任南琪,王愛杰.厭氧生物技術(shù)原理與應(yīng)用 [M]. 北京:化學(xué)工業(yè)出版社, 2004.
[11] 趙杰紅,張 波,蔡偉民.溫度對廚余垃圾兩相厭氧消化中水解和酸化過程的影響 [J]. 環(huán)境科學(xué), 2006,27(8):1682-1686.
[12] 李 東,孫永明,袁振宏等.原料比例和pH值對廚余垃圾和廢紙聯(lián)合厭氧消化的影響 [J]. 過程工程學(xué)報(bào), 2009,9(1):53-57.
[13] 吳 云.廚余垃圾厭氧消化影響因素及動力學(xué)研究 [D]. 重慶:重慶大學(xué)資源及環(huán)境科學(xué)學(xué)院, 2009.
[14] 劉國濤.城鎮(zhèn)有機(jī)垃圾水解及其水解液厭氧消化工藝與性能研究 [D]. 重慶:重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院, 2007
[15] Sun K H, Sang H K, Hang S S. UASB treatment of wastewater with VFA and alcohol generated during hydrogen fermentation of food waste [J]. Process Biochemistry, 2005,40(8):2897-2905.
[16] Angelidaki I,Ellegaard L,Ahring B K. Amathematical model for dynamic simulation of anaerobic digestion of complex substrates: focusing on ammonia inhibition [J]. Biotech Bioeng, 1993,43: 159-166.
Revision of methane generation rate empirical model in anaerobic digestion of kitchen wastes.
WU Yun1*, ZHANG Dai-jun2(1.Engineering Research Centre for Waste Oil Recovery Technology and Equipment, Ministry of Education, Chongqing Technology and Business University, Chongqing 400067, China；2.Resource and Environment Science College, Chongqing University, Chongqing 400030,China). China Environmental Science, 2011,31(5)：789～794
Fitted for cumulative CH4curve by single and two phrase anaerobic digestion of kitchen wastes experiment, respectively. Cumulative CH4curve of single phrase digestion complied with quadratic function relation(R2=0.9982) and cumulative CH4curve of two phrase anaerobic digestion complied with linear function relation(R2=0.9972); On the premise of considering system load, the fiducial empirical model of methane generation rate was formulated and established, furthermore, the fiducial empirical model was revised by adding coefficients of pH inhibition and ammonia nitrogen inhibition; This model was verified by experimental data. The statistic of test of single phrase digestion were 486.22 and 268.63,in condition that level of test was 0.01; the statistic of test of two phrase digestion were 421.69 and 372.56,in condition that level of test was 0.01. These statistics all were higher than critical value [F0.99(1,58)=7.12]. The empirical model could predict the production of cumulative CH4in the single and two phrase anaerobic digestion of kitchen wastes. Moreover, the calculation process of empirical model was simple.
kitchen wastes；anaerobic digestion；methane generation；empirical model
X705
A
1000-6923(2011)05-0789-06
2010-09-18
國家“十一五”科技攻關(guān)項(xiàng)目(2003BA808A17-02-03)
* 責(zé)任作者, 副研究員, wooyun123@163.com
吳 云(1973-),男,重慶人,副研究員,博士,主要從事固體廢物資源化利用及油污染治理方面的研究.發(fā)表論文20余篇.
猜你喜歡
濕垃圾與病死豬混合厭氧消化產(chǎn)氣性能研究環(huán)境衛(wèi)生工程(2021年5期)2021-11-20 05:45:24 Meso-mechanical model of concrete under a penetration loadDefence Technology(2019年6期)2020-01-07 09:11:04 基于PI+重復(fù)控制的單相逆變器研究通信電源技術(shù)(2016年4期)2016-04-04 02:57:12 煤層氣排采產(chǎn)氣通道適度攜煤粉理論中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版）(2015年2期)2015-11-10 06:07:43 零級動力學(xué)模型產(chǎn)甲烷量估算值與實(shí)測值的比較研究河北建筑工程學(xué)院學(xué)報(bào)(2015年2期)2015-04-29 12:23:52 一種新型斬波AC/DC/AC變換的單相DVR電測與儀表(2015年11期)2015-04-09 11:46:34 普通一級模型產(chǎn)甲烷量估算河北建筑工程學(xué)院學(xué)報(bào)(2015年4期)2015-03-24 09:48:11 丙烯酸鹽及對甲苯磺酸鹽對乙酸、丙酸產(chǎn)甲烷活性的影響環(huán)境工程技術(shù)學(xué)報(bào)(2014年5期)2014-07-19 06:44:22 采用干擾觀測器PI控制的單相SPWM逆變電源電測與儀表(2014年11期)2014-04-04 09:21:36 產(chǎn)氣剖面測試在蘇里格氣田東區(qū)的應(yīng)用石油化工應(yīng)用(2014年8期)2014-03-11 17:40:02
吳 云,張代鈞 (1.重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心,重慶 400067；.重慶大學(xué)資源與環(huán)境科學(xué)學(xué)院,重慶 400030)
廚余垃圾厭氧消化產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ偷男拚芯?/p>
吳 云1*,張代鈞2(1.重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心,重慶 400067；2.重慶大學(xué)資源與環(huán)境科學(xué)學(xué)院,重慶 400030)
分別通過單相和兩相厭氧消化實(shí)驗(yàn),對廚余垃圾產(chǎn)甲烷累積產(chǎn)量曲線進(jìn)行了擬合.結(jié)果表明,單相厭氧消化產(chǎn)甲烷累積產(chǎn)量曲線符合二次函數(shù)關(guān)系,相關(guān)系數(shù)(R2)為 0.9982;兩相厭氧消化產(chǎn)甲烷累積產(chǎn)量曲線符合線性函數(shù)關(guān)系,相關(guān)系數(shù)(R2)為 0.9972.以擬合函數(shù)為基礎(chǔ),在考慮系統(tǒng)負(fù)荷的前提下,建立了基準(zhǔn)狀態(tài)下的產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ?通過引入 pH抑制系數(shù)和氨氮抑制系數(shù)對建立的產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行了修正.對經(jīng)驗(yàn)修正模型預(yù)測結(jié)果進(jìn)行了實(shí)驗(yàn)驗(yàn)證及誤差分析.結(jié)果表明, 在檢驗(yàn)水平a=0.01條件下,單相厭氧第1組和第2組試驗(yàn)數(shù)據(jù)試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為486.22和268.63,均遠(yuǎn)大于臨界值(F0.99(1,58)=7.12);兩相厭氧第1組試驗(yàn)數(shù)據(jù)和第2組試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為421.69和372.56,也遠(yuǎn)高于其臨界值(F0.99(1,58)=7.12).說明該模型能很好地預(yù)測單相和兩相厭氧消化產(chǎn)甲烷累積產(chǎn)量,且計(jì)算過程簡單.
廚余垃圾；厭氧消化；產(chǎn)甲烷；經(jīng)驗(yàn)?zāi)Ｐ?/p>
對于以產(chǎn)甲烷為主要目的的廚余垃圾厭氧消化處理系統(tǒng)而言,甲烷產(chǎn)量決定著系統(tǒng)運(yùn)行的成敗.因此,長期以來,各國學(xué)者對厭氧消化產(chǎn)甲烷數(shù)學(xué)模型的研究從未停止[1-3].2002年,國際水協(xié)(IWA)正式公布了厭氧消化 1號數(shù)學(xué)模型(ADM1),該模型是集大成的研究成果,較為完整和準(zhǔn)確地描述了復(fù)雜厭氧消化系統(tǒng)的各個子過程,為厭氧工藝的建模和模擬提供了一個通用的平臺[4].ADM1模型推出后,一些關(guān)于模型補(bǔ)充和修正的研究工作陸續(xù)展開,如 Fedorovich等[5]提出了包含硫酸鹽還原的擴(kuò)展模型;Batstone[6]提出了利用氫氣的氧化描述硫酸鹽還原過程的簡便方法; Parker等[7]以ADM1為基礎(chǔ),提出了包含產(chǎn)氣及其釋放的擴(kuò)展模型,從而將產(chǎn)氣模型的機(jī)理性描述提高到一個新的水平.與此同時,一些新的建模思想和方法也在陸續(xù)出現(xiàn),如王凱軍等[8]采用系統(tǒng)動力學(xué)理論構(gòu)建了厭氧消化數(shù)學(xué)模型(SD模型),將變量間的各種復(fù)雜關(guān)系轉(zhuǎn)換為輸入輸出關(guān)系,較為簡便地實(shí)現(xiàn)了不同模擬方法之間的有效契合.王洪濤等[9]也以厭氧消化過程中碳素轉(zhuǎn)移機(jī)理為基本構(gòu)建思想,提出了基于碳素轉(zhuǎn)移機(jī)理的厭氧消化產(chǎn)氣模型.
但是,各種描述產(chǎn)甲烷過程的模型,都存在一個共同的問題,其計(jì)算工作十分繁雜,這給模型的實(shí)際應(yīng)用造成了不便.而在某些實(shí)際工作中,人們并不關(guān)心模型對厭氧消化過程細(xì)節(jié)的描述,只看重宏觀上的運(yùn)行結(jié)果,以便為可行性分析和投資決策提供相對準(zhǔn)確的依據(jù).通常,為了簡便快速地預(yù)測運(yùn)行結(jié)果,常常通過有機(jī)物總量或產(chǎn)甲烷潛能來粗略估算甲烷產(chǎn)量和產(chǎn)生速率,而這種估算的結(jié)果與實(shí)際運(yùn)行結(jié)果往往存在很大差異,導(dǎo)致投資決策出現(xiàn)偏差.本研究針對這些問題,在黑箱理論思想的指導(dǎo)下,通過實(shí)驗(yàn)對單相和兩相厭氧消化產(chǎn)甲烷規(guī)律進(jìn)行考察和分析,采用函數(shù)擬合的方法,構(gòu)建介于準(zhǔn)確描述和粗略估計(jì)之間的簡化經(jīng)驗(yàn)?zāi)Ｐ?以達(dá)到簡便預(yù)測甲烷產(chǎn)量和產(chǎn)生速率的目的.
1 材料和方法
1.1 材料
廚余垃圾樣品采自中餐廳、西餐廳、火鍋店和學(xué)生食堂.樣品采回后,剔除其中的泡沫飯盒、玻璃、石塊和金屬等雜物,并按重慶市主城區(qū)廚余垃圾基礎(chǔ)數(shù)據(jù)調(diào)研分析報(bào)告中提供調(diào)查數(shù)據(jù)進(jìn)行配比,即比例為:中餐類/湯鍋火鍋類/西餐類為0.65/0.30/0.05,將樣品用破碎機(jī)破碎后混勻,采用4分法取出,放入4℃冰箱中保存待用.廚余垃圾的基本物料性質(zhì)見表 1.實(shí)驗(yàn)用原生污泥取自重慶市唐家橋污水廠厭氧消化池,污泥含水率為94.23%,MLSS為47.36g/L,VSS為32.48g/L.經(jīng)馴化培養(yǎng)后作為接種污泥.
表1 廚余垃圾性質(zhì)(%)Table 1 Characteristics of kitchen wastes(%)
1.2 試驗(yàn)方法
單相厭氧消化試驗(yàn)采用試驗(yàn)裝置(圖1),試驗(yàn)前將物料破碎為粒徑小于10mm的顆粒.根據(jù)文獻(xiàn)中關(guān)于接種污泥混合比的報(bào)道[10-12],污泥和原料的配比從1:10到1:1不等,當(dāng)接種污泥比例低時,產(chǎn)氣滯后期較長,實(shí)驗(yàn)進(jìn)行時間較長,接種污泥比例高時,產(chǎn)氣波動較大,誤差也較大.而根據(jù)課題組前期的實(shí)驗(yàn)結(jié)果[13],采用2:1的配比時,消化情況較好,產(chǎn)氣情況也較穩(wěn)定.因此,選擇2:1的比例,即將1.5kg廚余垃圾和0.75kg活性污泥填入反應(yīng)器中,然后加水至 3L左右.試驗(yàn)過程中調(diào)節(jié)pH為7左右,控制溫度為35℃.試驗(yàn)進(jìn)行8d后,開始每天投加0.3kg新鮮垃圾.新鮮垃圾與當(dāng)天收集的50mL滲濾液混合,由反應(yīng)器頂部加入,并加水保持反應(yīng)物體積在 3L左右.穩(wěn)定運(yùn)行后每日由底部出料口出料1次,出料量為0.2～0.4kg.整個試驗(yàn)共持續(xù) 60d.試驗(yàn)過程中采用排飽和氯化鈉溶液法收集氣體,并跟蹤測定甲烷產(chǎn)量.
兩相厭氧消化采用圖 2所示試驗(yàn)裝置,試驗(yàn)開始時,將 2.86kg經(jīng)破碎處理的廚余垃圾填入水解產(chǎn)酸相反應(yīng)器,并加水 5.4L.活性污泥與廚余垃圾的接種比例為 1:2[10-12].試驗(yàn)過程中調(diào)節(jié)系統(tǒng)pH值為7左右,并控制溫度為35℃.系統(tǒng)啟動運(yùn)行8d后,每天將0.5kg新鮮垃圾與UASB出水混合后添加到產(chǎn)酸相反應(yīng)器,水解后的出水進(jìn)入水解液平衡池,并根據(jù)產(chǎn)酸相情況確定每天的出料量.
UASB反應(yīng)器內(nèi)接種污泥量為2.6L,為了防止水解液 pH值過低影響甲烷相運(yùn)行,進(jìn)水前調(diào)節(jié)水解液pH7左右,并由計(jì)量泵定時將水解液送入U(xiǎn)ASB反應(yīng)器.為了維持UASB反應(yīng)器內(nèi)一定的上升流速和溫度,設(shè)置了外循環(huán)水池,通過外循環(huán)加熱將反應(yīng)器水溫控制在 33～37℃.進(jìn)水泵間歇運(yùn)行,每小時運(yùn)行10min,由定時器控制,進(jìn)水泵流量為10～12mL/min.外循環(huán)泵24h運(yùn)行,循環(huán)流量為 12～24mL/min.試驗(yàn)過程中采用排飽和氯化鈉溶液法收集氣體,并跟蹤測定甲烷產(chǎn)量.
圖1 單相厭氧消化試驗(yàn)裝置示意Fig.1 Sketch of experiment equipment of single phrase anaerobic process
圖2 兩相厭氧消化試驗(yàn)裝置示意Fig.2 Sketch of experiment equipment of two phrase anaerobic process
試驗(yàn)采用GASBOARD紅外氣體分析儀(武漢光電科技)來測定氣體成分及含量(體積百分?jǐn)?shù)).甲烷氣體含量測定結(jié)果確定后,將收集氣體體積乘以氣體分析儀測定得到的體積百分?jǐn)?shù),即可計(jì)算出甲烷產(chǎn)量;氨氮濃度的測定采用納氏試劑比色法進(jìn)行測定(上海精科的 721分光光度計(jì));pH值采用PHS-3C型pH計(jì)進(jìn)行測定.
2 結(jié)果與分析
由圖3可見,單相厭氧消化甲烷累積產(chǎn)量曲線趨勢與二次函數(shù)曲線非常相似,兩相厭氧消化甲烷累積產(chǎn)量曲線趨勢則基本與線性函數(shù)曲線相似.因此,將3次實(shí)驗(yàn)結(jié)果取平均值,并分別采用二次函數(shù)和線性函數(shù)對甲烷累積產(chǎn)量的平均值曲線進(jìn)行擬合.
圖3 單相及兩相厭氧消化甲烷累積產(chǎn)量平行實(shí)驗(yàn)Fig.3 Cumulative CH4 curve of single and two phrase digestion repeated experiment
從圖3和圖4中可以看到,二次函數(shù)和線性函數(shù)的擬合效果很好,R2分別達(dá)到了 0.9982和0.9972.經(jīng)分析后認(rèn)為,在適宜的溫度(35℃)、pH6～7和有機(jī)負(fù)荷條件下,單相和兩相厭氧消化甲烷累積產(chǎn)量曲線的這種變化趨勢與其各自的反應(yīng)器運(yùn)行特性是基本吻合的.因?yàn)?對于單相厭氧消化而言,系統(tǒng)的環(huán)境條件會隨著有機(jī)物水解、產(chǎn)酸、產(chǎn)乙酸和產(chǎn)甲烷過程發(fā)生變化,由于單相厭氧未將產(chǎn)甲烷過程分置進(jìn)行,因此系統(tǒng)的這些過程直接影響著甲烷菌產(chǎn)甲烷速率的高低,且產(chǎn)甲烷速率整體上表現(xiàn)為前期低、后期高的趨勢,這就導(dǎo)致其甲烷累積曲線表現(xiàn)為二次曲線形態(tài).而兩相厭氧消化系統(tǒng)將產(chǎn)甲烷過程分置到甲烷相中進(jìn)行,甲烷菌所受系統(tǒng)環(huán)境條件的影響較小,基本能夠保持穩(wěn)定產(chǎn)氣狀態(tài),其甲烷累積曲線也就基本表現(xiàn)為線性曲線趨勢[13].
圖4 單相及兩相厭氧消化甲烷累積產(chǎn)量擬合曲線Fig.4 Cumulative CH4 AV curve and fitting curve of single and two phrase digestion
3 簡化模型的建立與修正
3.1 基準(zhǔn)模型的建立
根據(jù)上述實(shí)驗(yàn)條件,由于物料供給充足;對系統(tǒng)溫度、pH值、氨氮濃度等環(huán)境條件進(jìn)行了嚴(yán)格監(jiān)測和控制,使其保持在適宜的消化狀態(tài).因此,在建立基準(zhǔn)模型前,可作如下假設(shè):①假設(shè)厭氧消化過程中基質(zhì)濃度充足,厭氧菌群生長代謝不受基質(zhì)濃度限制;②假設(shè)厭氧消化進(jìn)行過程中,系統(tǒng)溫度、pH值、氨氮濃度等環(huán)境因素處于適宜狀態(tài),不對厭氧消化過程造成抑制影響;③假設(shè)系統(tǒng)未超負(fù)荷運(yùn)行時,累積產(chǎn)氣量與厭氧消化的有機(jī)物負(fù)荷量呈正比關(guān)系.
從圖3中可看到,不同試驗(yàn)條件下的甲烷累積曲線與簡單函數(shù)曲線相類似,圖3a的甲烷累積曲線與二次函數(shù)曲線線類似,圖3b的甲烷累積曲線則與線性函數(shù)曲線類似,擬合結(jié)果表明,實(shí)驗(yàn)結(jié)果與擬合曲線相關(guān)系數(shù)達(dá)到0.95以上,證實(shí)擬合曲線方程能夠很好地預(yù)測甲烷累積產(chǎn)量.因此,可將甲烷累積產(chǎn)量擬合曲線方程作為簡化的甲烷累積產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ头匠?經(jīng)驗(yàn)?zāi)Ｐ头匠痰谋磉_(dá)式如式(1)和式(2):
單相累積甲烷產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ头匠?
式中: t為厭氧消化時間; G為甲烷累積產(chǎn)量,L.
再將上述方程對時間求導(dǎo)即可得到簡化的產(chǎn)氣速率模型方程.對式(1)和式(2)求導(dǎo)得:
單相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
兩相累積甲烷產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ头匠?
需要說明的是,從式(4)可看出,在本實(shí)驗(yàn)條件下,建立的兩相厭氧消化的產(chǎn)氣速率方程接近常數(shù),這與劉國濤[14]和Sun等[15]的觀察結(jié)果基本一致.
根據(jù)假設(shè),在未超負(fù)荷情況下,甲烷產(chǎn)量與有機(jī)物負(fù)荷量呈正比.因此,在考慮系統(tǒng)有機(jī)負(fù)荷的情況下,可引入負(fù)荷比A/A0,A為實(shí)際有機(jī)負(fù)荷,A0為基準(zhǔn)有機(jī)負(fù)荷,此處的基準(zhǔn)有機(jī)負(fù)荷即為前述實(shí)驗(yàn)中的系統(tǒng)每日有機(jī)負(fù)荷量.則式(3)和式(4)可變?yōu)?
單相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
兩相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
兩相產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?
因此根據(jù)模型假設(shè)和實(shí)驗(yàn)條件,可以將式(5)和式(6)稱為基準(zhǔn)狀態(tài)下的產(chǎn)甲烷速率經(jīng)驗(yàn)?zāi)Ｐ?
3.2 基準(zhǔn)模型的修正
由于式(5)和式(6)是在假設(shè)條件和特定實(shí)驗(yàn)條件下建立的基準(zhǔn)狀態(tài)產(chǎn)氣速率經(jīng)驗(yàn)?zāi)Ｐ?而在實(shí)際情況下,還應(yīng)考慮系統(tǒng) pH值和氨氮抑制性對產(chǎn)氣速率的影響,因此有必要對式(5)和式(6)進(jìn)行修正.根據(jù)ADM1模型中提出的關(guān)于pH值抑制系數(shù)和氨氮抑制系數(shù)的相關(guān)理論[4,16],pH值抑制系數(shù)和氨氮抑制系數(shù)可分別表達(dá)如下:
IpH=1,pH=7(中性條件產(chǎn)甲烷系統(tǒng)不受酸堿抑制)
式中:pH代表系統(tǒng)實(shí)際pH值;pHUL和pHLL是生物種群受到50%抑制時的pH值上限和下限,一般分別取7和6[4]; S1是氨氮濃度;K1是氨氮抑制參數(shù).因此,在實(shí)際情況下,可引入上述抑制系數(shù)對基準(zhǔn)經(jīng)驗(yàn)?zāi)Ｐ褪?5)和式(6)進(jìn)行修正.修正后的經(jīng)驗(yàn)?zāi)Ｐ捅磉_(dá)式如下:
單相產(chǎn)氣速率經(jīng)驗(yàn)修正模型:
兩相產(chǎn)氣速率經(jīng)驗(yàn)修正模型:
式中: IpH為 pH抑制系數(shù);INH3為氨氮抑制系數(shù),上述經(jīng)驗(yàn)?zāi)Ｐ偷膯挝粸長/d. 根據(jù)建立基準(zhǔn)經(jīng)驗(yàn)產(chǎn)氣速率模型的試驗(yàn)條件,可令單相厭氧的A0=0.3kg/d,兩相厭氧的A0=0.5kg/d.
4 經(jīng)驗(yàn)?zāi)Ｐ偷尿?yàn)證
4.1 試驗(yàn)曲線與預(yù)測曲線的比較
單相厭氧消化驗(yàn)證試驗(yàn)所用物料性質(zhì)和試驗(yàn)條件與1.1和1.2所述實(shí)驗(yàn)條件基本一致,但將有機(jī)負(fù)荷分別控制為 0.25,0.36kg,并在試驗(yàn)過程中控制系統(tǒng)pH值約為5.5左右(改變系統(tǒng)pH條件);因此當(dāng)系統(tǒng)pH值為5.5時,其系統(tǒng)氨氮平均濃度為 0.826g/L左右,另根據(jù)文獻(xiàn)報(bào)道,K1可取45.063[4],因此,經(jīng)抑制系數(shù)公式計(jì)算,IpH和INH3分別為0.389和0.982.
兩相厭氧消化驗(yàn)證試驗(yàn)所用物料性質(zhì)和試驗(yàn)條件與1.1和1.2所述實(shí)驗(yàn)條件基本一致,僅將有機(jī)負(fù)荷分別變?yōu)?0.3,0.75kg.在上述實(shí)驗(yàn)條件下,系統(tǒng) pH值為 7,甲烷相氨氮平均濃度為1.446g/L,經(jīng)抑制系數(shù)公式計(jì)算,得IpH和INH3分別為1和0.967.模型驗(yàn)證結(jié)果如圖7和圖8所示.
從圖5中可以看出,甲烷累積產(chǎn)量模型預(yù)測曲線與試驗(yàn)曲線基本吻合,說明簡化的產(chǎn)氣模型能夠用于預(yù)測單相和兩相厭氧消化甲烷產(chǎn)氣量.
圖5 單相及兩相厭氧消化甲烷累積產(chǎn)量和預(yù)測曲線Fig.5 cumulative CH4 Curve and predict curve of single and two phrase digestion
4.2 誤差分析
將不同試驗(yàn)條件下得到的單相厭氧消化和兩相厭氧消化甲烷累積產(chǎn)量試驗(yàn)數(shù)據(jù)與模型預(yù)測值進(jìn)行方差檢驗(yàn).檢驗(yàn)結(jié)果表明,單相厭氧第 1組和第 2組試驗(yàn)數(shù)據(jù)試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為486.22和268.63,在檢驗(yàn)水平α=0.01下均遠(yuǎn)大于臨界值[F0.99(1,58)=7.12];兩相厭氧第1組試驗(yàn)數(shù)據(jù)和第2組試驗(yàn)數(shù)據(jù)與模型預(yù)測值的顯著性檢驗(yàn)統(tǒng)計(jì)量F分別為421.69和372.56,也遠(yuǎn)高于其臨界值[F0.99(1,58)=7.12].說明簡化的厭氧消化產(chǎn)氣修正經(jīng)驗(yàn)?zāi)Ｐ偷幕貧w是顯著的,因此該模型能夠較好地預(yù)測厭氧消化甲烷累積產(chǎn)量.
5 結(jié)論
5.1 在一定的實(shí)驗(yàn)條件下,單相厭氧消化甲烷累積產(chǎn)量曲線與二次函數(shù)曲線相似,兩相厭氧消化甲烷累積產(chǎn)量曲線與線性函數(shù)曲線相似.擬合結(jié)果表明,二次函數(shù)和線性函數(shù)曲線能夠很好地對單相和兩相厭氧消化甲烷累積產(chǎn)量曲線進(jìn)行擬合,R2分別達(dá)到了0.9982和0.9972.
5.2 在模型假設(shè)基礎(chǔ)上,建立了基準(zhǔn)狀態(tài)下的甲烷累積產(chǎn)量簡化經(jīng)驗(yàn)?zāi)Ｐ秃彤a(chǎn)甲烷速率簡化經(jīng)驗(yàn)?zāi)Ｐ?通過引入 pH值抑制系數(shù)和氨氮抑制系數(shù)對模型進(jìn)行修正,建立了修正后的甲烷累積產(chǎn)量簡化經(jīng)驗(yàn)?zāi)Ｐ秃彤a(chǎn)甲烷速率簡化經(jīng)驗(yàn)?zāi)Ｐ?
5.3 模型驗(yàn)證結(jié)果表明,簡化的甲烷累積產(chǎn)量經(jīng)驗(yàn)?zāi)Ｐ皖A(yù)測曲線與試驗(yàn)曲線基本吻合,簡化經(jīng)驗(yàn)?zāi)Ｐ湍軌蛴糜陬A(yù)測單相和兩相厭氧消化甲烷累積產(chǎn)量;誤差分析結(jié)果表明,對簡化經(jīng)驗(yàn)?zāi)Ｐ颓€的方差檢驗(yàn)結(jié)果是顯著的,同樣也說明,簡化經(jīng)驗(yàn)?zāi)Ｐ湍軌蛴糜陬A(yù)測甲烷累積產(chǎn)量.
[1] Christensen T H.Landfill of waste:biogas [M]. London: E and FNSpon, 1996,56-74.
[2] Schumacher M M.Landfill methane recovery [M]. USA: NewJersey,1983,35-46.
[3] El-fadel M, Findikakis A N. A numerical model for methane production in managed sanitary landfills [J]. Waste Management and Research, 1989,7(1):31-42.
[4] 國際水協(xié)厭氧消化工藝數(shù)學(xué)模型課題組.厭氧消化數(shù)學(xué)模型[M]. 張亞雷,周雪飛譯.上海:同濟(jì)大學(xué)出版社, 2004:15-30.
[5] Fedorovich V, Lens P, KalyuzhnyiI S, et al. Extension of anaerobic digestion model No.1 with processes of sulfate reduction [J]. Applied Biochemistry and Biotechnology, 2003, 109(1-3):33-45.
[6] Batstone D J. Mathematical modelling of anaerobic reactors treating domestic wastewater: rational criteria for model use[J]. Reviews in Environmental Science and Bio-Technology, 2006,5(1):57-71.
[7] Parker W J, WU G H. Modifying ADM1 to include formation and emission of odourants [J]. Water Science and Technology, 2006,54(4):111-117.
[8] 閻 中,王凱軍.厭氧消化過程系統(tǒng)動力學(xué)模型構(gòu)建方法研究[J]. 中國沼氣, 2009,27(2):3-8.
[9] 楊渤京,王洪濤,陸文靜,等.厭氧條件下纖維素降解產(chǎn)氣的研究[J]. 中國環(huán)境科學(xué), 2008,28(3):255-259.
[10] 任南琪,王愛杰.厭氧生物技術(shù)原理與應(yīng)用 [M]. 北京:化學(xué)工業(yè)出版社, 2004.
[11] 趙杰紅,張 波,蔡偉民.溫度對廚余垃圾兩相厭氧消化中水解和酸化過程的影響 [J]. 環(huán)境科學(xué), 2006,27(8):1682-1686.
[12] 李 東,孫永明,袁振宏等.原料比例和pH值對廚余垃圾和廢紙聯(lián)合厭氧消化的影響 [J]. 過程工程學(xué)報(bào), 2009,9(1):53-57.
[13] 吳 云.廚余垃圾厭氧消化影響因素及動力學(xué)研究 [D]. 重慶:重慶大學(xué)資源及環(huán)境科學(xué)學(xué)院, 2009.
[14] 劉國濤.城鎮(zhèn)有機(jī)垃圾水解及其水解液厭氧消化工藝與性能研究 [D]. 重慶:重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院, 2007
[15] Sun K H, Sang H K, Hang S S. UASB treatment of wastewater with VFA and alcohol generated during hydrogen fermentation of food waste [J]. Process Biochemistry, 2005,40(8):2897-2905.
[16] Angelidaki I,Ellegaard L,Ahring B K. Amathematical model for dynamic simulation of anaerobic digestion of complex substrates: focusing on ammonia inhibition [J]. Biotech Bioeng, 1993,43: 159-166.
Revision of methane generation rate empirical model in anaerobic digestion of kitchen wastes.
WU Yun1*, ZHANG Dai-jun2(1.Engineering Research Centre for Waste Oil Recovery Technology and Equipment, Ministry of Education, Chongqing Technology and Business University, Chongqing 400067, China；2.Resource and Environment Science College, Chongqing University, Chongqing 400030,China). China Environmental Science, 2011,31(5)：789～794
Fitted for cumulative CH4curve by single and two phrase anaerobic digestion of kitchen wastes experiment, respectively. Cumulative CH4curve of single phrase digestion complied with quadratic function relation(R2=0.9982) and cumulative CH4curve of two phrase anaerobic digestion complied with linear function relation(R2=0.9972); On the premise of considering system load, the fiducial empirical model of methane generation rate was formulated and established, furthermore, the fiducial empirical model was revised by adding coefficients of pH inhibition and ammonia nitrogen inhibition; This model was verified by experimental data. The statistic of test of single phrase digestion were 486.22 and 268.63,in condition that level of test was 0.01; the statistic of test of two phrase digestion were 421.69 and 372.56,in condition that level of test was 0.01. These statistics all were higher than critical value [F0.99(1,58)=7.12]. The empirical model could predict the production of cumulative CH4in the single and two phrase anaerobic digestion of kitchen wastes. Moreover, the calculation process of empirical model was simple.
kitchen wastes；anaerobic digestion；methane generation；empirical model
X705
A
1000-6923(2011)05-0789-06
2010-09-18
國家“十一五”科技攻關(guān)項(xiàng)目(2003BA808A17-02-03)
* 責(zé)任作者, 副研究員, wooyun123@163.com
吳 云(1973-),男,重慶人,副研究員,博士,主要從事固體廢物資源化利用及油污染治理方面的研究.發(fā)表論文20余篇.
猜你喜歡
環(huán)境衛(wèi)生工程(2021年5期)2021-11-20 05:45:24
Defence Technology(2019年6期)2020-01-07 09:11:04
通信電源技術(shù)(2016年4期)2016-04-04 02:57:12
中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版）(2015年2期)2015-11-10 06:07:43
河北建筑工程學(xué)院學(xué)報(bào)(2015年2期)2015-04-29 12:23:52
電測與儀表(2015年11期)2015-04-09 11:46:34
河北建筑工程學(xué)院學(xué)報(bào)(2015年4期)2015-03-24 09:48:11
環(huán)境工程技術(shù)學(xué)報(bào)(2014年5期)2014-07-19 06:44:22
電測與儀表(2014年11期)2014-04-04 09:21:36
石油化工應(yīng)用(2014年8期)2014-03-11 17:40:02
