接種比對豬糞中溫條件產(chǎn)甲烷潛力的影響

馬文林，王安吉

（1.北京建筑大學，北京應對氣候變化研究和人才培養(yǎng)基地，北京102616；2.北京建筑大學，環(huán)境與能源工程學院，北京102616）

為減緩氣候變化，《巴黎協(xié)定》確定了到2100 年末將全球平均氣溫升幅控制在不超過工業(yè)化前水平2 ℃的升溫目標，并努力限制在1.5 ℃之內(nèi)[1]。甲烷對于全球升溫的影響作用，在所有溫室氣體種類中排位第二，僅次于二氧化碳。聯(lián)合國糧農(nóng)組織統(tǒng)計表明，2020 年我國因畜禽糞便管理引起的甲烷排放占我國農(nóng)業(yè)活動甲烷排放總量的8.42%[2]。因此，畜禽糞污管理導致的甲烷排放量對全球升溫有不可忽視的潛在負面影響力，需要采取適當?shù)墓芾聿呗訹3]。我國是全球生豬飼養(yǎng)量最大的國家，Li 等[4]的研究得出一頭豬在其畜牧生命周期內(nèi)因糞便管理產(chǎn)生1.34 kg 甲烷排放量，準確評價我國豬糞管理的甲烷排放量，既利于明確我國在全球溫室氣體排放管理中的責任擔當，也利于摸清家底，更有針對性地制定豬糞管理甲烷減排措施，為我國農(nóng)業(yè)領(lǐng)域乃至全國實現(xiàn)碳中和提供科技支持。

IPCC《國家溫室氣體清單指南》給出了畜禽糞便管理甲烷排放量的計算方法，即糞便管理甲烷排放量等于活動數(shù)據(jù)[以糞便中揮發(fā)性固體（VS）質(zhì)量計的糞便量] 乘以排放因子（每單位質(zhì)量糞便VS 的甲烷排放質(zhì)量）。最大產(chǎn)甲烷能力（B0）是估算甲烷排放因子的關(guān)鍵參數(shù)[5]，IPCC《國家溫室氣體清單指南》針對全球各大洲給出了豬糞便管理B0缺省值（以VS 計），亞洲地區(qū)為0.29 m3·kg-1，使用亞洲水平的缺省值評估我國豬糞管理甲烷排放情況，估算結(jié)果會有較高的不確定性。

產(chǎn) 甲 烷 潛 力（Biochemical methane potential，BMP）可用于評估有機物的產(chǎn)甲烷活性[6-7]，它是直接獲取B0值的常用實驗方法。但BMP 測試具有條件復雜和周期長等特點[8-9]，且至今無統(tǒng)一的測試方法，這在很大程度上影響了BMP 值測試結(jié)果的準確性和可比性。劉翌晨等[10]研究得出，不同生長期豬糞的BMP（以VS計）有顯著不同，妊娠母豬（230 mL·g-1）顯著低于保育豬（273 mL·g-1）與生長豬（272 mL·g-1）；徐文倩等[11]于37 ℃中溫條件下得出豬糞的BMP（以VS計）為0.38 m3·kg-1；楊紫怡[12]以污水處理廠污泥為接種物在37 ℃中溫條件下得出豬糞的BMP 為427 mL·g-1；Wang 等[13]研究了中溫條件下預熱處理對豬糞BMP的影響；張丹志等[14]分析了豬糞厭氧發(fā)酵有機酸代謝與微生物的相關(guān)性；王思淇等[15]研究了豬糞厭氧發(fā)酵過程中揮發(fā)性脂肪酸對微生物群落的影響。以上研究表明，不同學者采用不同方法或不同來源接種液對不同豬糞進行BMP 和相關(guān)指標的分析，得出的測試結(jié)果存在較大差異，這將影響到B0值和甲烷排放量估算結(jié)果的準確性，有可能加大糞便管理甲烷排放量和減排項目減排量估算的誤差，導致誤判糞便管理甲烷減排潛力。因此，為使不同研究者得出的BMP 具有可比性，取得符合我國實際排放水平的豬糞管理B0值，有必要盡早制定規(guī)范統(tǒng)一的BMP 測試方法。

根據(jù)BMP 測試方法可知，測試溫度、接種液微生物活性、底物與接種液比例、甲烷氣體計量方法等多方面因素都可能影響B(tài)MP 測試結(jié)果。本研究選取有機物厭氧反應最普遍使用的中溫條件（35 ℃）、以微晶纖維素為參照物表征接種液微生物活性、利用能夠?qū)淄闅怏w體積自動計量并進行標準狀態(tài)換算的全自動產(chǎn)甲烷潛力儀，研究不同底物與接種液比例對豬糞BMP 的影響，以期為提出豬糞BMP 統(tǒng)一測試方法提供科學依據(jù)。該方法的提出，對降低豬糞管理甲烷排放量核算結(jié)果不確定性至關(guān)重要。同時，本研究還對不同接種比反應系統(tǒng)的產(chǎn)甲烷過程中消化液的理化指標進行連續(xù)監(jiān)測，并對反應結(jié)束后的微生物群落進行分析，以探索接種比對豬糞產(chǎn)甲烷過程中的理化指標變化規(guī)律和微生物菌群特征的影響，為開展相關(guān)機理分析提供參考借鑒。

1 材料與方法

1.1 樣品來源與處理

本實驗所用豬糞樣品采集于北京市東部郊區(qū)某規(guī)?；B(yǎng)殖場的母豬舍。該養(yǎng)殖場總存欄量1 055頭，包括仔豬380 頭、育肥生長豬570 頭和母豬105頭，豬舍采用干清糞方式。從母豬舍漏縫地板采集鮮母豬糞樣品，樣品于實驗室內(nèi)35 ℃下進行風干處理，除雜后研磨過40目篩備用。

本實驗以豬糞為底物，以北京市北部某規(guī)模養(yǎng)殖場沼氣工程（2 個容積為50 m3的沼氣發(fā)酵罐）排放的沼液為接種液。在實驗室內(nèi)對沼液進行除雜處理，并于35 ℃環(huán)境下進行為期1個月的培養(yǎng)后備用。

在實驗開始前分別用燃燒法測定接種液與豬糞的總固體（TS）和灰分含量（以質(zhì)量計），計算得到其VS值，樣品TS、VS與灰分值見表1。

表1 接種液與豬糞的TS、VS與灰分值（%）Table 1 TS，VS and ash values of inoculum and pig manure（%）

1.2 實驗裝置

本研究利用系統(tǒng)Ⅰ（AMPTS-Ⅱ，碧普）和系統(tǒng)Ⅱ（MultiTalent 203，碧臣）兩套全自動甲烷潛力測試裝置進行豬糞中溫厭氧發(fā)酵實驗，如圖1 所示，系統(tǒng)Ⅰ和系統(tǒng)Ⅱ的反應瓶反應體積分別為0.4 L和1.6 L。各系統(tǒng)均由厭氧消化單元、二氧化碳固定單元與氣體體積測定單元3 部分構(gòu)成。樣品消化單元由玻璃發(fā)酵瓶、恒溫水浴鍋與自動攪拌裝置組成。二氧化碳固定單元由若干個容積為100 mL 的玻璃瓶組成，實驗開始前向每個玻璃瓶裝入80 mL 3 mol·L-1的NaOH溶液以吸收厭氧消化單元排出的二氧化碳、硫化氫等酸性氣體。為保證不同環(huán)境條件下所產(chǎn)生甲烷的體積計量結(jié)果具有可比性，儀器將測得的產(chǎn)氣量數(shù)據(jù)自動換算為標準狀況（0 ℃、1個標準大氣壓）下的體積。

圖1 BMP裝置示意圖Figure 1 Schematic diagram of the BMP apparatus

1.3 實驗設計

系統(tǒng)Ⅰ與系統(tǒng)Ⅱ采用相同的接種液和豬糞樣品，各實驗組初始總VS 水平（豬糞+接種液）均設為3%，反應溫度以恒溫水浴控制在（35.0±0.5）℃。根據(jù)實驗設定的各反應瓶的反應體積、物料VS 含量、初始VS濃度以及物料接種比，計算各反應瓶的物料投入水平，接種比以VS 計，計算結(jié)果見表2。在各反應瓶添加相應物料后，用純水補充至相應體積。本實驗以豬糞為反應底物，沼液為接種液，將豬糞與接種液的比例分別設置為1∶2、1∶1 與2∶1，即接種比（S/I）分別為0.5、1.0和2.0。

表2 系統(tǒng)Ⅰ和系統(tǒng)Ⅱ各實驗組物料投入情況（g）Table 2 Material input for each experimental group in system Ⅰand system Ⅱ（g）

系統(tǒng)Ⅰ用于測試豬糞的BMP值，共15個反應瓶，各反應瓶反應體積設為0.4 L，分為空白組（CK 組，僅含接種液）、纖維素對照組（CG組，纖維素和接種液的比例為1∶2）和處理組（M1 組、M2 組和M3 組，S/I 分別為0.5、1.0 和2.0）3 類5 個實驗組，每組設3 個平行實驗。實驗期間，自動攪拌系統(tǒng)以120 r·min-1的速度間歇進行攪拌工作。BMP 測試裝置每日自動記錄各反應瓶的甲烷產(chǎn)生量，當甲烷日產(chǎn)量小于累積產(chǎn)氣量的1%時系統(tǒng)自動停止日產(chǎn)氣量記錄，同時終止實驗。

系統(tǒng)Ⅱ用于連續(xù)監(jiān)測豬糞的中溫厭氧消化過程，只設空白組和3 個處理組共4 個實驗組，不設對照組和平行實驗，各反應瓶反應體積設為1.6 L。自動攪拌系統(tǒng)以120 r·min-1的速度間歇進行攪拌工作。在反應過程中定期從反應瓶中取發(fā)酵混合液樣品，將其離心并過0.45 μm 濾膜后，對濾液進行pH、氨氮、總氮與葡萄糖濃度檢測，各項理化指標的檢測方法見表3。

表3 系統(tǒng)Ⅱ中各實驗組相關(guān)指標測試方法Table 3 Test methods of relevant indexes for each experimental group in System Ⅱ

在實驗結(jié)束后從系統(tǒng)Ⅱ各反應瓶中取出一定體積污泥進行16S rRNA 測試，檢測厭氧消化系統(tǒng)的細菌與產(chǎn)甲烷古菌的群落結(jié)構(gòu)，擴增區(qū)域為V4區(qū)，所用測序平臺為Novaseq-PE250，引物序列（5′→3′）如下：細菌引物為515F（GTGCCAGCMGCCGCGGTAA）和806R（GGACTACHVGGGTWTCT AAT）；產(chǎn)甲烷古菌 引 物 為1106F（TTWAGTCAGGCAACGAGC）和1378R（TGTGCAAGGAGCAGG GAC）。

1.4 厭氧消化動力學分析

將所得產(chǎn)甲烷數(shù)據(jù)利用修正的Gompertz 模型[16]與Logistic 模型[17]進行動力學擬合。以均方根誤差（RMSE）和決定系數(shù)（R2）作為模型的評估指標。RMSE 常用于評價模型精度，該值越小模型精度越高[18]，R2用于評價模型擬合優(yōu)度，R2越高模型擬合度越好[19]。在進行動力學擬合前，將所得豬糞的日甲烷產(chǎn)量與累積甲烷產(chǎn)量數(shù)據(jù)根據(jù)各組物料投入水平換算成均以VS 計的日產(chǎn)甲烷速率（R，mL·g-1·d-1）與累積產(chǎn)甲烷量（M，mL·g-1）。修正的Gompertz模型與Logistic模型如公式（1）與公式（2）所示。

式中：M（t）為t時的累計產(chǎn)甲烷量（以VS 計），mL·g-1；A為最大產(chǎn)甲烷潛力（以VS計），mL·g-1；Rmax為最大產(chǎn)甲烷速率（以VS 計），mL·g-1·d-1；e 為自然常數(shù)；λ為延滯期，d；t為發(fā)酵時間，d。

2 結(jié)果與分析

2.1 BMP測試

各實驗組日甲烷產(chǎn)量與累積甲烷產(chǎn)量變化如圖2所示。由圖2A可知，CK組第2天和第3天無甲烷產(chǎn)出，第4～6 天有少量甲烷產(chǎn)生，于第6 天自動結(jié)束產(chǎn)氣記錄，表明接種液中有機物匱乏，基本不再產(chǎn)甲烷。3 個豬糞處理組在整個產(chǎn)甲烷過程中有相似的產(chǎn)甲烷變化規(guī)律。M1 組、M2 組、M3 組3 個處理組產(chǎn)甲烷的第一個峰值出現(xiàn)在第2 天，第3 天產(chǎn)氣量出現(xiàn)小幅下降，第4 天開始回升，第5 天出現(xiàn)第二個產(chǎn)氣峰值。該峰值是整個產(chǎn)甲烷過程中日產(chǎn)氣量最高值，各處理組最大日甲烷產(chǎn)量從大到小依次排序為M3 組（410.05 mL）＞M2 組（326.90 mL）＞M1 組（208.40 mL），表明豬糞接種比例越高，日產(chǎn)甲烷最大峰值越高，最大日甲烷產(chǎn)量與豬糞接種比例呈顯著正相關(guān)關(guān)系（P＜0.05）。第6 天3 個處理組日產(chǎn)氣量驟降，下降幅度強于之前第一個峰值后出現(xiàn)的降幅。第7 天產(chǎn)氣量再次回升，但回升能力顯著低于第二次峰值時的情況，M1 組和M2 組于第7 天出現(xiàn)第3 個產(chǎn)氣峰值，M3組的第3 個峰值出現(xiàn)在第8 天，第三個峰值的產(chǎn)氣量由高到低的順序與甲烷最大產(chǎn)量順序相同。在第三個產(chǎn)氣峰值出現(xiàn)后各處理組的日產(chǎn)氣量均呈現(xiàn)穩(wěn)定下降趨勢，M1、M2 和M3 組分別于第15 天、第16 天和第18天結(jié)束產(chǎn)氣。

圖2 不同處理組日甲烷產(chǎn)量和累積甲烷產(chǎn)量Figure 2 Daily methane and cumulative methane yields in different treatment groups

由圖2B 可以看出，接種比例越高，總的產(chǎn)氣時間越長、累積產(chǎn)氣量越大。M1、M2 和M3 組累積甲烷產(chǎn)量分別為1 068.5、1 558.2 mL和2 054.6 mL，豬糞累積甲烷產(chǎn)量與其接種比例顯著正相關(guān)（P＜0.05）。

圖3 為豬糞處理組與微晶纖維素組的BMP 測試結(jié)果。由圖3 可知，M1、M2、M3 與CG 組產(chǎn)甲烷潛力（以VS 計）分別為（252.05±5.12）、（250.07±16.57）、（249.65±13.69）mL·g-1與（302.03±21.73）mL·g-1，不同接種比例處理組的BMP 無顯著差異（P＞0.05）。各豬糞處理組實驗所得豬糞BMP 平均值為（250.59±1.28）mL·g-1。

圖3 對照組和各處理組產(chǎn)甲烷潛力測試結(jié)果Figure 3 Results of BMP tests for control group and treatment groups

2.2 pH的變化

各實驗組厭氧發(fā)酵反應過程的pH變化情況見圖4。pH是影響厭氧發(fā)酵過程的重要條件因素，對微生物繁殖代謝、酶活性、物料溶解等有重要影響。由圖4可看出，不同實驗組組間的發(fā)酵上清液pH變化具有較大差異，且無明顯規(guī)律性。但各組在厭氧消化結(jié)束后發(fā)酵上清液的pH均較初始狀態(tài)有所提升。

圖4 不同實驗組的pHFigure 4 pH in different experimental groups

2.3 氮的變化

各實驗組厭氧發(fā)酵反應過程的氨氮與總氮變化情況見圖5。從圖5A 可以看出，M1、M2 和M3 這3 組的氨氮濃度在厭氧消化前4 d 都呈上升趨勢。M1 組第4 天出現(xiàn)一個峰值（772.50 mg·L-1），且在第6 天迅速下降，降幅為8.16%。M2 組與M3 組第6 天的氨氮濃度也較第4 天有所下降，降幅分別為2.08%與3.81%。在反應結(jié)束后，M1、M2 和M3 組的發(fā)酵上清液中的氨氮濃度較初始值分別增加28.37%、43.25%與60.85%。

圖5 不同實驗組的氮濃度Figure 5 Nitrogen concentration in different experimental groups

對于總氮（圖5B），各組BMP 測試完成后的發(fā)酵上清液中的總氮濃度均較初始狀態(tài)有所提高。M1、M2 和M3 組厭氧反應結(jié)束后發(fā)酵上清液中的總氮濃度較初始值分別增加了27.23%、37.15%與80.49%。上清液中氨氮與總氮濃度較初始值增長幅度與接種比例呈顯著正相關(guān)（P＜0.05）。

2.4 葡萄糖的變化

圖6 為不同實驗組葡萄糖濃度的變化趨勢。從圖6 可以看出，各處理組發(fā)酵上清液中的葡萄糖濃度在厭氧消化期間呈下降-上升-再下降的波動變化，且各處理組的葡萄糖濃度均于第6 天出現(xiàn)了明顯的峰值。

將圖6與圖2進行比較觀察發(fā)現(xiàn)，第2天與第4天葡萄糖濃度均較前一次檢測結(jié)果有所降低，且當日甲烷產(chǎn)量均較前一天有所提升，第4 天各豬糞處理組的發(fā)酵上清液的葡萄糖濃度處于最低值，各處理組日產(chǎn)甲烷量均較第3 天有所提高且均出現(xiàn)了日產(chǎn)甲烷峰值。當反應進行至第6 天時，各處理組均出現(xiàn)了葡萄糖濃度峰值，而日產(chǎn)甲烷量均呈現(xiàn)產(chǎn)甲烷低谷。在各組反應結(jié)束前1～2 d，上清液中葡萄糖濃度均有所升高。由此表明，產(chǎn)甲烷過程與葡萄糖消耗有關(guān)，日甲烷產(chǎn)量高時，當日葡萄糖濃度與前日相比較就處于下降情況，反之則升高。

2.5 細菌與產(chǎn)甲烷古菌的變化

參與厭氧發(fā)酵反應的微生物群系包括細菌和產(chǎn)甲烷古菌。因此，反應結(jié)束后，檢測了系統(tǒng)Ⅱ各實驗組細菌與產(chǎn)甲烷古菌的相對豐度，相對豐度位于前20 的菌種情況如圖7與圖8所示。

圖7 不同實驗組細菌相對豐度Figure 7 Relative abundance of bacterial in different experimental groups

圖8 不同實驗組產(chǎn)甲烷古菌相對豐度Figure 8 Relative abundance of methanogenic archaea in different experimental groups

由圖7 可以看出，基于門水平，CK 組前3 種優(yōu)勢菌群為Bacteroidota（擬桿菌門）、Proteobacteria（變形菌門）與Firmicutes（厚壁菌門）。各豬糞處理組中豐度居前3 位的菌種均為厚壁菌門、擬桿菌門與Cloacimonadota，變形菌門與擬桿菌門的相對豐度均較CK組明顯下降，而厚壁菌門與Cloacimonadota 的相對豐度明顯上升。其中變形菌門在M1、M2 和M3 組的相對豐度由CK 的21.32%分別下降至4.95%、4.37%和3.09%，下降程度與接種比例呈正相關(guān)。對于豬糞處理組，在M2 組檢測到最高的擬桿菌門相對豐度與最低的厚壁菌門相對豐度。Thermotogota 和Synergistota的相對豐度均與接種比例呈負相關(guān)性。

基于屬水平，CK 組豐度居前3 的菌種為Bacteroidetes_VC2.1_Bac22、Pseudomonas（假單胞菌屬）與Thiopseudomonas（硫代單胞菌屬），豬糞處理組的Bacteroidetes_VC2.1_Bac22、Pseudomonas與Thiopseudomonas的相對豐度均較CK 組明顯下降，而DMER64、Clostridium_sensu_stricto_1與Candidatus_Cloacimonas的相對豐度較CK 組明顯上升，其中DMER64的相對豐度與接種比呈正相關(guān)性。

從圖8 可知，從門水平看，無論是否添加豬糞，Halobacterota 均是產(chǎn)甲烷古菌中的絕對優(yōu)勢菌。M1、M2 和M3 組中的Halobacterota 較CK 組的相對豐度有所提升。從屬水平看，Methanosaeta在各組中均為最優(yōu)勢菌種。相較CK 組，豬糞處理組的Methanosaeta的相對豐度有明顯提升，Methanobacterium與Methanospirillum的相對豐度明顯下降。Methanobacterium在M1、M2 和M3 組的相對豐度分別為CK 組的71.22%、74.42%與84.79%，相對豐度下降幅度與接種比例呈負相關(guān)。

2.6 厭氧消化動力學分析

圖9 為中溫條件下采用兩種模型對不同接種比豬糞的累積產(chǎn)甲烷量的擬合結(jié)果。表4 為中溫條件不同接種比下兩種模型對累積產(chǎn)甲烷量的擬合參數(shù)。

由圖9 與表4 可知，在中溫條件下修正的Gompertz 模型預測的各接種比的最大產(chǎn)甲烷潛力分別為247.40 mL·g-1（S/I=0.5）、244.10 mL·g-1（S/I=1.0）和246.39 mL·g-1（S/I=2.0），其R2保持在0.991～0.997 之間，R2均大于0.99，擬合程度良好。對于Logistic 模型，預測的各接種比的最大產(chǎn)甲烷潛力分別為244.19 mL·g-1（S/I=0.5）、240.19 mL·g-1（S/I=1.0）和241.74 mL·g-1（S/I=2.0），其R2在0.981～0.991 之間，R2大于0.9，擬合程度良好。修正的Gompertz 模型的RMSE 數(shù) 值 范 圍 在4.67～6.97 之 間，Logistic 模 型 的RMSE 數(shù)值范圍在7.70～10.12 之間。通過對比R2與RMSE 可知，修正的Gompertz 模型較Logistic 模型能更好地對中溫條件下不同接種比例豬糞的產(chǎn)甲烷潛力進行擬合。

3 討論

3.1 實驗條件對測試結(jié)果的影響

本研究結(jié)果表明，針對同一物料，采用相同的測試溫度和初始VS 水平時，接種比對BMP 測試結(jié)果有影響，但差異不顯著。而Santos 等[20]以市政污水處理廠污泥為接種液，設定0.65 與1 兩種水平的接種比，測得葡萄牙中部地區(qū)養(yǎng)殖場豬糞的BMP 值分別為0.329 m3·kg-1與0.568 m3·kg-1，不同接種比下的測試結(jié)果有較大差異。兩項研究實驗方案的差異主要是Santos等[20]的實驗沒有采取固定的初始VS濃度。

微晶纖維素BMP 測試結(jié)果用來表征接種液的活性程度，因此，該實驗組被作為BMP 測試的對照組（CG 組）。本研究測得CG 組的BMP 值為（302.03±21.73）mL·g-1。Zhu 等[21]于37 ℃條件下測得的微晶纖維素BMP 值為183.32 mL·g-1，同時其依據(jù)Buswell 公式計算出微晶纖維素的理論BMP 值為390.21 mL·g-1，本研究測得的微晶纖維素BMP 值達到理論產(chǎn)值的77%，說明本次實驗所用的接種液活性良好。接種液活性對于BMP 測試結(jié)果有重要影響，規(guī)范接種液制備方法（包括接種液來源和接種液培養(yǎng)馴化方法），可使得不同研究對各種有機物料測得的BMP 值具有可比性。也有其他研究者進行了豬糞BMP 值的測定，例如：Cao 等[22]將食品廢棄物中溫消化器污泥在37 ℃條件下培養(yǎng)3 d 后作為接種液，以美國俄亥俄州養(yǎng)殖場豬糞為測試對象，在（37.0±0.5）℃中溫條件下測得其BMP 值為0.476 m3·kg-1；Wang 等[23]于35 ℃、6%初始TS濃度和7∶3的豬糞與接種液比例下得出豬糞的BMP 為239.00 mL·g-1；Xu 等[24]利用葡萄糖與醋酸鈉對污水處理廠污泥培養(yǎng)20 d 后作為接種液，在35 ℃、4%初始TS 濃度、底物與接種液的比例為1.14（以VS計）時得出豬糞的BMP為221.3 mL·g-1。

根據(jù)上述研究結(jié)果可以發(fā)現(xiàn)，不同研究者的研究結(jié)果間存在較大差異。歸納起來，豬糞來源、接種液來源和培養(yǎng)馴化方法、反應溫度、實驗初始VS 或TS水平是造成差異的主要原因。因此，為使得不同測試者對豬糞的BMP 測試結(jié)果具有可比性，應對接種液制備方法、實驗初始VS 水平和反應溫度這3 方面關(guān)鍵測試條件進行統(tǒng)一規(guī)定。

3.2 氨氮變化量與日甲烷產(chǎn)量變化量的相關(guān)關(guān)系

高濃度氨氮對厭氧消化過程有抑制作用[25]。通過分析本實驗厭氧消化液中氨氮含量與日甲烷產(chǎn)量的關(guān)系，發(fā)現(xiàn)氨氮濃度日變化量與日甲烷產(chǎn)量變化量之間存在一定的相互關(guān)系，因此分別計算了各處理組第2天較第0天、第4天較第2天、第6天較第4天、第9天較第6 天、第13 天較第11 天與第15 天較第13 天的日氨氮濃度變化量與甲烷日產(chǎn)量變化，結(jié)果見圖10。

圖10 不同處理組相關(guān)天數(shù)的氨氮濃度與日甲烷產(chǎn)量變化Figure 10 Variations of ammonia-N concentration and daily methane production on relevant days in different experimental groups

有機物通過厭氧消化反應轉(zhuǎn)化為甲烷，同時也伴隨氨氮的產(chǎn)生。從圖2 曲線可得出，在豬糞厭氧消化前9 d，厭氧反應體系經(jīng)歷了3 個產(chǎn)甲烷高峰。從圖10A 可看出，在整個反應過程中，氨氮總體上呈現(xiàn)上升趨勢，但期間存在往復上升與下降的波動情況。圖10 同時表示厭氧反應體系中氨氮含量日變化量與甲烷日產(chǎn)生量變化情況。從圖10 可以看出，厭氧消化早期（0～7 d），由于豬糞中可降解有機物含量多，厭氧消化反應快，氨氮濃度變化和日產(chǎn)甲烷量變化都大，且氨氮變化量與日產(chǎn)甲烷量變化情況呈負相關(guān)關(guān)系。例如，第0～6天，M1、M2和M3處理組的氨氮濃度變化量依次降低，對應各組的甲烷日產(chǎn)生量增幅則逐一加大；第4～6 天，M1 組氨氮濃度下降最大，對應的日產(chǎn)甲烷增長量最低；在厭氧消化中期（6～9 d），隨厭氧消化反應的進行，豬糞中可降解有機物減少，反應體系的日產(chǎn)氣量及其日變化幅度都相應減少，但同時間段氨氮濃度增加趨勢未減緩，對比各組間氨氮濃度變化和日產(chǎn)甲烷量變化之間的關(guān)系可知，二者呈正相關(guān)關(guān)系；厭氧消化末期（9 d以后）與中期不同，氨氮濃度變化率明顯下降，其變化量與產(chǎn)甲烷變化量之間依然呈正相關(guān)關(guān)系。

3.3 日產(chǎn)甲烷量與葡萄糖的相關(guān)性

本研究中不同處理組均在第2 天出現(xiàn)首個日產(chǎn)甲烷峰值，然后產(chǎn)氣量下降，這是因為豬糞中所含的易降解有機物經(jīng)厭氧水解反應轉(zhuǎn)化生成的產(chǎn)甲烷底物被消耗殆盡。隨反應時間延長，豬糞中的大分子有機物進一步被轉(zhuǎn)化為產(chǎn)甲烷菌群可以利用的反應底物，因此在第5天出現(xiàn)了第二個日產(chǎn)甲烷峰值。

乙酸是產(chǎn)甲烷的直接底物，而葡萄糖發(fā)酵是產(chǎn)生乙酸的途徑之一，葡萄糖經(jīng)乙酸路徑產(chǎn)甲烷的過程如下：

C6H12O6+2H2O→2CH3COOH+2CO2+4H2

4H2+CO2→2H2O+CH4

CH3COOH→CO2+CH4

由該方程組可知，發(fā)酵液中葡萄糖濃度的增加，有利于反應向產(chǎn)甲烷過程偏移。觀察圖2 中第5～7天日產(chǎn)甲烷量和圖6 中同期葡萄糖濃度變化可知，第6 天日甲烷產(chǎn)量較第5 天顯著降低的原因或為反應器內(nèi)大分子有機物分解生成的葡萄糖濃度不足所致。這也為第7 天各處理組日甲烷產(chǎn)量均較第6 天出現(xiàn)明顯增長提供了解釋基礎(chǔ)，第7 天發(fā)酵液中葡萄糖濃度增加促使乙酸濃度上升，促進了產(chǎn)甲烷菌的產(chǎn)甲烷作用，因此出現(xiàn)了日甲烷產(chǎn)量較第6 天有所增長的現(xiàn)象。

產(chǎn)甲烷菌利用底物產(chǎn)甲烷具有一定的濃度閾值，當?shù)孜餄舛鹊陀谝欢ㄖ禃r，不發(fā)生產(chǎn)甲烷過程[26-29]。此外，互營降解有機物對于代謝產(chǎn)物的轉(zhuǎn)移具有一定依賴性，代謝產(chǎn)物和甲烷產(chǎn)生速度與互營菌和產(chǎn)甲烷菌間的距離呈正相關(guān)性，有機質(zhì)的類型和降解產(chǎn)物也會對產(chǎn)甲烷過程進行調(diào)節(jié)控制[30]。在本研究中，當反應結(jié)束時各處理組葡萄糖濃度從大到小依次為M1[（43.92±0.11）mg·L-1] ＞M2[（38.96±0.11）mg·L-1] ＞M3[（27.54±0.11）mg·L-1] 。關(guān)于上清液中葡萄糖濃度與日產(chǎn)甲烷量的定量關(guān)系有待進一步探索。

4 結(jié)論

（1）中溫條件下，初始揮發(fā)性固體濃度為3%時，母豬糞產(chǎn)甲烷潛力（BMP）為（250.59±1.28）mL·g-1，不同接種比例豬糞的產(chǎn)甲烷潛力無顯著差異（P＞0.05）。從優(yōu)化物料用量與時間成本考慮，豬糞BMP 測試的豬糞與接種液比例（以VS計）取1∶2為宜。

（2）豬糞累積甲烷產(chǎn)量和最大日甲烷產(chǎn)量均與豬糞接種比例呈顯著正相關(guān)性（P＜0.05）。

（3）豬糞厭氧消化過程中，厚壁菌門、擬桿菌門與Cloacimonadota為三大優(yōu)勢細菌，Halobacterota在古菌中占絕對優(yōu)勢。接種比對細菌和古菌中不同菌種的豐度均有影響。

（4）修正的Gompertz 模型較Logistic 模型更適合中溫條件下豬糞的厭氧消化動力學分析。