響應(yīng)面法優(yōu)化干清豬糞半干法連續(xù)厭氧發(fā)酵試驗參數(shù)

陳文成 袁存亮 張重 李靖 裴海林 邢向欣 孫京 趙國明

摘要為實現(xiàn)干清豬糞連續(xù)高效厭氧發(fā)酵，以干清豬糞為原料進行連續(xù)厭氧發(fā)酵試驗，利用實驗室自制連續(xù)厭氧發(fā)酵反應(yīng)器探究水力滯留時間（HRT）、攪拌間隔與攪拌轉(zhuǎn)速對干清豬糞連續(xù)厭氧發(fā)酵的影響，采用響應(yīng)面分析法進行參數(shù)優(yōu)化，得到最優(yōu)參數(shù)組合：HRT 19 d、攪拌間隔2 h、攪拌轉(zhuǎn)速30 r/min。經(jīng)過干清豬糞連續(xù)厭氧發(fā)酵試驗驗證，其平均日產(chǎn)氣量達23.48 L，VS降解率56.47%，平均容積產(chǎn)氣率2.35 L/（L·d），能夠?qū)崿F(xiàn)干清豬糞連續(xù)厭氧發(fā)酵的高效穩(wěn)定運行。因此，干清豬糞的無害化處理可采用半干法連續(xù)厭氧發(fā)酵運行，從而實現(xiàn)干清豬糞的資源化利用。

關(guān)鍵詞干清豬糞；半干法；連續(xù)厭氧發(fā)酵；參數(shù)優(yōu)化

中圖分類號X 713文獻標識碼A

文章編號0517-6611（2023）24-0196-06

doi：10.3969/j.issn.0517-6611.2023.24.044

Optimization of the Continuous Anaerobic Fermentation Experiment Parameters of Semidried Pig Manure by Using Response Surface Method

CHEN Wencheng， YUAN Cunliang， ZHANG Zhong et al

（Jilin Academy of Agricultural Machinery， Changchun， Jilin 130022）

AbstractIn order to realize the continuous and efficient anaerobic fermentation of dry pig manure， the continuous anaerobic fermentation test was carried out with dry pig manure as raw materials. The effects of hydraulic retention time （HRT）， stirring interval and stirring speed on the continuous anaerobic fermentation of dry pig manure were studied by using the selfmade continuous anaerobic fermentation reactor in the laboratory.The response surface analysis method was used to optimize the experiment parameters，and the optimal parameter combination was obtained as follows： HRT 19 d， stirring interval 2 h， and stirring speed 30 r/min.The average daily gas yield was up to 23.48 L after the continuous anaerobic fermentation test of dry pig manure.VS degradation rate was 56.47%， and the average volumetric gas production rate was 2.35 L/（L· d ）， which could achieve efficient and stable operation of dry pig manure continuous anaerobic fermentation.Therefore， the harmless treatment of dry pig manure could be operated by semidry continuous anaerobic fermentation， so as to realize the resource utilization of dry pig manure.

Key wordsDry pig manure;Semidried method;Continuous anaerobic fermentation;Parameter optimization

隨著我國人民生活水平的不斷提高，對肉蛋奶的需求量也逐年遞增，從而導(dǎo)致我國畜牧業(yè)的飛速發(fā)展［1-3］，但目前生豬仍是畜牧業(yè)養(yǎng)殖的重要飼養(yǎng)品種之一［4］。截至2019年我國生豬出欄量為54 419.3萬頭，而每頭生豬每日排放廢水約13.5 kg、糞便1.5 kg［5-6］，養(yǎng)殖場排除的糞污及廢水若不經(jīng)處理將直接影響生態(tài)環(huán)境，污染水土資源，因此豬糞的無害化處理至關(guān)重要［7-9］。鄭盼［2］在豬糞干式厭氧消化試驗中發(fā)現(xiàn)，豬糞干發(fā)酵產(chǎn)氣率較濕發(fā)酵產(chǎn)氣率增加22.1%，豬糞干式厭氧發(fā)酵工藝最佳參數(shù)為發(fā)酵溫度34 ℃，接種量40%；朱圣權(quán)［10］通過對豬糞厭氧干發(fā)酵工藝研究及其效益分析，優(yōu)化出適宜的厭氧發(fā)酵工藝參數(shù)，經(jīng)過20～22 d厭氧發(fā)酵，其原料日產(chǎn)甲烷量可達1 600 mL/（kg·d）。目前我國大規(guī)模生豬養(yǎng)殖企業(yè)及畜禽養(yǎng)殖舍大多采用干清糞方式進行糞污收集，其優(yōu)點是節(jié)水、污染小，但會導(dǎo)致糞便干物質(zhì)質(zhì)量濃度高［11-13］。采取半干法連續(xù)厭氧發(fā)酵技術(shù)來處理干清豬糞是一種很好的技術(shù)模式，可以使產(chǎn)生的糞污得到及時有效地處理，并且產(chǎn)生清潔能源。為達到良好的處理效果和較高的生產(chǎn)效率，需要考慮各參數(shù)的綜合控制，以使微生物生長繁殖達到最佳平衡狀態(tài)，因此，為實現(xiàn)厭氧發(fā)酵的高產(chǎn)氣率及合格的降解率，對連續(xù)厭氧發(fā)酵的水力滯留時間、攪拌參數(shù)進行研究與優(yōu)化，從而實現(xiàn)豬糞的無害化處理。

1材料與方法

1.1試驗材料試驗所用干清豬糞取自吉林省長春市九臺區(qū)小袁家窩堡養(yǎng)殖戶的干清豬糞；接種物取自吉林省農(nóng)業(yè)機械研究院廢棄物利用實驗室，是牛糞經(jīng)過充分厭氧反應(yīng)后直至不再產(chǎn)氣的活性污泥。干清豬糞及接種物的理化性質(zhì)如表1所示。

1.2試驗儀器試驗所用儀器具體見表2。

1.3試驗裝置

試驗裝置為實驗室自主研制的臥式三軸攪拌厭氧發(fā)酵反應(yīng)器（圖1），反應(yīng)器容積約10 L。整個厭氧發(fā)酵系統(tǒng)（圖2）包括電控系統(tǒng)、水箱循環(huán)系統(tǒng)、反應(yīng)器、脫硫脫水裝置、濕式氣體流量計、集氣裝置等。通過電控系統(tǒng)，控制電機帶動反應(yīng)器內(nèi)部的攪拌葉翅來實現(xiàn)對發(fā)酵物料的攪拌。通過調(diào)速器調(diào)節(jié)電機轉(zhuǎn)速，通過時控開關(guān)來控制攪拌電機的開啟與關(guān)閉，從而滿足不同的攪拌間隔與攪拌時長；通過水箱循環(huán)系統(tǒng)以及測溫傳感器和控制系統(tǒng)，可實現(xiàn)對發(fā)酵物料的溫度控制。此外，反應(yīng)器還設(shè)有進料口和出料口，可實現(xiàn)連續(xù)運行中進出料的要求。

1.4試驗設(shè)計

1.4.1水力滯留時間（HRT）對干清豬糞半干法連續(xù)厭氧發(fā)酵的影響試驗。以干清豬糞為發(fā)酵原料，有效反應(yīng)容積為反應(yīng)器容積的60%，HRT逐漸降低（有機負荷逐漸升高），通過HRT換算出進料量，再通過進料量的改變來研究HRT對干清豬糞半干法連續(xù)厭氧發(fā)酵日產(chǎn)氣量和VS降解率的影響。HRT設(shè)定25、20、15、10 d 4個水平，每天固定時間進出料各1次，先出料，后進料。進料體積的計算如公式（1）所示，出料體積的計算如公式（2）所示。該試驗采取中溫發(fā)酵（35 ℃），料液濃度為半干法厭氧發(fā)酵濃度（干清豬糞+接種物混合濃度），攪拌參數(shù)為上述試驗篩選出適宜的攪拌參數(shù)，反應(yīng)初始以50%接種率啟動，產(chǎn)氣平穩(wěn)后開始連續(xù)進出料，每天同一時間記錄每日產(chǎn)氣量并測定出料時的VS含量，換算出其VS降解率。

式中：V進為反應(yīng)器連續(xù)運行過程中的每日進料體積，單位mL；V出為反應(yīng)器連續(xù)運行過程中的每日出料體積，單位mL；V下為出料前反應(yīng)器內(nèi)的發(fā)酵料液下降的體積，單位mL；V有效容積為反應(yīng)器內(nèi)的發(fā)酵料液體積，單位mL；HRT為反應(yīng)器運行過程中的水力滯留時間，單位d。

1.4.2干清豬糞半干法連續(xù)厭氧發(fā)酵參數(shù)優(yōu)化試驗。

采用Design Expert 8.0.6的Box-Behnken中心組合設(shè)計原理，在中溫（35 ℃）、發(fā)酵料液濃度為半干法厭氧發(fā)酵濃度的條件下，設(shè)計水力滯留時間（HRT）、攪拌間隔和攪拌轉(zhuǎn)速3個因素3個水平，以每組試驗在產(chǎn)氣穩(wěn)定后的平均日產(chǎn)氣量和VS降解率作為試驗結(jié)果評價指標，進行響應(yīng)面分析試驗，確定干清豬糞半干法連續(xù)厭氧發(fā)酵最佳工藝條件。正交試驗因素與水平設(shè)計如表3所示，對優(yōu)化后的參數(shù)進行試驗驗證。

1.5測定項目與方法每天09：00使用濕式氣體流量計進行讀數(shù)，記錄每日產(chǎn)氣情況；TS、VS含量采用灼燒法進行測定。

1.6數(shù)據(jù)處理與分析試驗數(shù)據(jù)通過Excel 2019軟件進行統(tǒng)計，并利用Design Expert 12軟件對試驗數(shù)據(jù)進行進一步分析處理。

2結(jié)果與分析

2.1HRT對干清豬糞半干法連續(xù)厭氧發(fā)酵的影響

圖3、圖4為HRT對干清豬糞半干法連續(xù)厭氧發(fā)酵日產(chǎn)氣量及VS降解率的影響。整個發(fā)酵時間為25 d，其中0～7 d為啟動階段。從圖3可以看出，所有試驗組均能正常啟動，HRT=15 d時的日產(chǎn)氣量明顯高于HRT=25、20、10 d時，在連續(xù)運行穩(wěn)定后的15 d內(nèi)，HRT=25、20、15、10 d時平均日產(chǎn)氣量分別為19.16、22.73、25.47和2.05 L，說明當HRT=25、20、15 d時隨著HRT的縮短，平均日產(chǎn)氣量升高。當HRT=10 d時，厭氧發(fā)酵系統(tǒng)內(nèi)部有機負荷過高，營養(yǎng)物質(zhì)不能及時消耗，酸積累過多，厭氧發(fā)酵系統(tǒng)遭到破壞，不再是一個正常的發(fā)酵系統(tǒng)。從圖4可以看出，在連續(xù)運行穩(wěn)定后的15 d內(nèi)，HRT=25、20、15、10 d時的平均VS降解率為62.42%、58.04%、51.51%和9.18%，說明隨著HRT的縮短，VS降解率逐漸降低。這是因為當HRT=25、20 d時，發(fā)酵物滯留期相對較長，日進料量相對較少，反應(yīng)充分，VS降解率也較高；當HRT=15 d時，日進料量增多，營養(yǎng)物質(zhì)充足，水解酸化產(chǎn)物更多，此時日產(chǎn)氣量高于HRT=25、20 d時，但是此時VS降解率有所降低，說明此時已不能完全消化營養(yǎng)物質(zhì)；當HRT=10 d時，日產(chǎn)氣量與VS降解率急劇下降，這可能是隨著進料量的增多，有機負荷過高，引起整個發(fā)酵系統(tǒng)酸化，形成病態(tài)池。綜上所述，干清豬糞在半干法連續(xù)厭氧發(fā)酵試驗中，在一定范圍內(nèi)HRT越長，降解效果越好。若縮短HRT，日產(chǎn)氣量會有所提高，但降解效果有所下降，若HRT過短則會引起發(fā)酵系統(tǒng)酸化，形成病態(tài)池，因此在干清豬糞半干法連續(xù)厭氧發(fā)酵中HRT選取20 d較為合適。

2.2干清豬糞半干法連續(xù)厭氧發(fā)酵參數(shù)的優(yōu)化

干清豬糞半干法連續(xù)厭氧發(fā)酵在中溫35℃條件下進行，按表3進行參數(shù)的響應(yīng)面優(yōu)化試驗，并以平均日產(chǎn)氣量與VS降解率為試驗結(jié)果評價指標，結(jié)果見表4。

利用Design Expert 8.0.6軟件對試驗數(shù)據(jù)進行分析擬合，并進行差異顯著性分析。平均日產(chǎn)氣量與VS降解率的回歸方程見公式（3）和（4）。

Y=24.06+0.23A+0.21B-0.11C-0.22AB-0.20AC+0.01BC-0.79A2-0.35B2-0.23C2（3）

Z=53.03+3.82A+0.70B+0.60C+0.14AB+0.57AC+0.35BC+1.00A2-1.87B2-3.17C2（4）

響應(yīng)面方差分析結(jié)果見表5～6。

由表5～6可知，該模型P值＜0.000 1，說明回歸模型極顯著。失擬項P值大于0.05，影響不顯著，說明該方程擬合度良好。通過分析各因素的一次項、二次項和交互項的P值可知，一次項因素A、B、C的P值均小于0.05，影響顯著；二次項因素A2、B2、C2的P值均小于0.05，影響顯著；交互項AB的P值小于0.05，影響顯著；交互項BC的P值大于0.05，說明影響不顯著。通過對F值的分析發(fā)現(xiàn)，各影響因素的主次順序為A>B>C，即HRT>攪拌間隔>攪拌轉(zhuǎn)速。

（1）對平均日產(chǎn)氣量的影響分析。基于回歸模型對平均日產(chǎn)氣量的方差分析，為了進一步研究相關(guān)因素及其交互作用并且確定最優(yōu)點，利用Design Expert 8.0.6軟件對回歸模型進行響應(yīng)面分析。當攪拌轉(zhuǎn)速一定時（攪拌轉(zhuǎn)速=30 r/min），平均日產(chǎn)氣量隨著HRT和攪拌間隔變化的三維曲面和等高線圖如圖5（a）、（b）所示。從圖5（a）、（b）可以看出，隨著HRT和攪拌間隔的增加，平均日產(chǎn)氣量先增加后減小，而且HRT對日產(chǎn)氣量的變化比攪拌間隔更為敏感，因此等高線圖呈現(xiàn)HRT為短軸、攪拌間隔為長軸的橢圓，所以HRT比攪拌間隔對平均日產(chǎn)氣量的影響要更為顯著，此結(jié)論與表5中HRT的F值大于攪拌間隔F值的結(jié)論相一致，證實了HRT的重要程度大于攪拌間隔。

圖6（a）、（b）為當攪拌間隔一定時（攪拌間隔=2 h），平均日產(chǎn)氣量隨著HRT和攪拌轉(zhuǎn)速變化的三維曲面和等高線圖。從圖6（a）可以看出，隨著HRT和攪拌轉(zhuǎn)速的增加，平均日產(chǎn)氣量先增加后減小。由圖6（b）可知，HRT比攪拌轉(zhuǎn)速對平均日產(chǎn)氣量的影響更為顯著。此結(jié)論與表5中HRT的F值大于攪拌轉(zhuǎn)速F值的結(jié)果相一致。

圖7（a）、（b）為HRT=16 d時平均日產(chǎn)氣量隨攪拌間隔和攪拌轉(zhuǎn)速變化的三維曲面和等高線圖。從圖7（a）可以看出，隨著攪拌間隔和攪拌轉(zhuǎn)速的增大，平均日產(chǎn)氣量變化不明顯。由圖7（b）可知，等高線圖大致呈圓形，說明攪拌轉(zhuǎn)速與攪拌間隔對平均日產(chǎn)氣量的影響基本相當，與表5中攪拌轉(zhuǎn)速F值與攪拌間隔F值大致相同的結(jié)果一致。

（2）對VS降解率的影響分析。

基于回歸模型對VS降解率的方差分析，利用Design Expert 8.0.6軟件對回歸模型進行響應(yīng)面分析。圖8（a）、（b）為攪拌轉(zhuǎn)速一定時（攪拌轉(zhuǎn)速=30 r/min），VS降解率隨著HRT和攪拌間隔變化的三維曲面和等高線圖。由圖8（a）可知，VS降解率隨著HRT的增大而增大。從8（b）可以看出，當HRT處于合適范圍內(nèi)，攪拌間隔的大幅變化對VS降解率的影響不大，反之當攪拌間隔處于合適范圍時，HRT小幅變化會對VS降解率產(chǎn)生較大影響，所以HRT比攪拌間隔對VS降解率的影響要更為顯著，此結(jié)論與表6中HRT的F值大于攪拌間隔F值的結(jié)果相一致。

圖9（a）、（b）為攪拌間隔一定時（攪拌間隔=2 h），VS降解率隨著HRT和攪拌轉(zhuǎn)速變化的三維曲面和等高線圖。由圖9（a）可知，VS降解率隨著HRT的增大而增大。從圖9（b）可以看出，HRT對VS降解率的影響要比攪拌轉(zhuǎn)速對VS降解率的影響更顯著，此結(jié)論可從表6中HRT的F值大于攪拌轉(zhuǎn)速F值得到驗證。

圖10（a）、（b）為HRT一定時（HRT=16 d），VS降解率隨著攪拌間隔和攪拌轉(zhuǎn)速變化的三維曲面圖和等高線圖。由圖10（b）可知，等高線圖大致呈圓形，說明攪拌轉(zhuǎn)速與攪拌間隔對VS降解率的影響基本相當，從表6中攪拌轉(zhuǎn)速F值與攪拌間隔F值大致相同，可以得到驗證。

2.3最佳條件的確定及試驗驗證

在干清豬糞半干法連續(xù)厭氧發(fā)酵試驗中，平均日產(chǎn)氣量和VS降解率越高越好，因此利用Design Expert 8.0.6軟件的優(yōu)化模塊將試驗指標都調(diào)至最大值，根據(jù)相關(guān)回歸模型分析得到以下參數(shù)：HRT為18.59d、攪拌間隔為2.11h、攪拌轉(zhuǎn)速為29.51r/min。預(yù)測最佳平均日產(chǎn)氣量為23.89 L、VS降解率為56.09%。

在進行最佳試驗條件的驗證時，考慮到實際應(yīng)用的簡便性，將干清豬糞半干法連續(xù)厭氧消化條件調(diào)整為HRT19 d、攪拌間隔2 h、攪拌轉(zhuǎn)速30 r/min，結(jié)果表明干清豬糞連續(xù)厭氧發(fā)酵的平均日產(chǎn)氣量實際值為23.48 L，VS降解率56.47%，實測值與理論值基本一致。

3結(jié)論

通過研究干清豬糞半干法連續(xù)厭氧發(fā)酵中HRT對平均日產(chǎn)氣量和VS降解率的影響，并利用響應(yīng)面分析法優(yōu)化干清豬糞半干法連續(xù)厭氧發(fā)酵中的最優(yōu)工藝條件，得到以下結(jié)論：

固定發(fā)酵溫度、攪拌參數(shù)、發(fā)酵料液濃度等條件不變，研究當HRT為25、20、15、10 d時對干清豬糞半干法連續(xù)厭氧發(fā)酵的影響，發(fā)現(xiàn)在一定范圍內(nèi)HRT越長，降解效果越好；縮短HRT，日產(chǎn)氣量會有所提高，但降解效果有所下降；若HRT過短時則會引起發(fā)酵系統(tǒng)酸化，形成病態(tài)池，因此在干清豬糞半干法連續(xù)厭氧發(fā)酵中HRT宜控制在20 d。

通過采用Design Expert 8.0.6的Box-Behnken 中心組合設(shè)計原理，在發(fā)酵溫度、料液濃度等其他因素一致的條件下，研究HRT、攪拌間隔和攪拌轉(zhuǎn)速3個因素對干清豬糞半干法連續(xù)厭氧發(fā)酵平均日產(chǎn)氣量和VS降解率的影響，發(fā)現(xiàn)各因素對日產(chǎn)氣量和VS降解率的影響重要性排序均為HRT＞攪拌間隔＞攪拌轉(zhuǎn)速，并通過模型尋優(yōu)與驗證，得到最優(yōu)參數(shù)組合為HRT 19 d、攪拌間隔2 h、攪拌轉(zhuǎn)速30 r/min，實際平均日產(chǎn)氣量23.48 L，VS降解率56.47%，平均容積產(chǎn)氣率2.35 L/（L·d）。

參考文獻

［1］李丹妮.豬場糞污厭氧干發(fā)酵酸氨抑制規(guī)律研究［D］.北京：中國農(nóng)業(yè)科學(xué)院，2021.

［2］ 鄭盼.豬糞干式厭氧消化實驗研究［D］.昆明：云南師范大學(xué)，2019.

［3］ 孔德望，張克強，房芳，等.豬糞厭氧發(fā)酵消化液回流體系微生物群落結(jié)構(gòu)特征與產(chǎn)氣關(guān)系研究［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2018，37（3）：559-566.

［4］ 許彩云，靳紅梅，杜靜，等.水力停留時間對豬糞厭氧發(fā)酵殘留物中磺胺類抗生素分布的影響［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2016，35（11）：2187-2194.

［5］ 蔡連.以豬糞為原料的高溫厭氧發(fā)酵產(chǎn)沼氣工藝研究［D］.南昌：南昌大學(xué)，2014.

［6］ 王悅超，雷中方.馴化接種對高固體濃度豬糞厭氧發(fā)酵的影響［J］.復(fù)旦學(xué)報（自然科學(xué)版），2012，51（1）：118-124.

［7］ 龐小平，牛明芬，王賽月，等.北方玉米秸稈與豬糞混合厭氧發(fā)酵條件［J］.生態(tài)學(xué)雜志，2011，30（1）：126-130.

［8］ 劉榮厚，郝元元，武麗娟.溫度條件對豬糞厭氧發(fā)酵沼氣產(chǎn)氣特性的影響［J］.可再生能源，2006，24（5）：32-35.

［9］ 張靖雨，汪邦穩(wěn)，夏小林，等.農(nóng)村規(guī)?；笄蒺B(yǎng)殖污染生態(tài)綜合治理技術(shù)研究進展［J］.安徽農(nóng)業(yè)科學(xué)，2020，48（19）：9-14，29.

［10］ 朱圣權(quán).豬糞厭氧干發(fā)酵工藝研究及其效益分析［D］.武漢：華中農(nóng)業(yè)大學(xué)，2009.

［11］ 沈其林，單勝道，周健駒，等.豬糞發(fā)酵沼液成分測定與分析［J］.中國沼氣，2014，32（3）：83-86.

［12］ 馬媛，鄧功成，趙洪，等.溫度對豬糞發(fā)酵產(chǎn)沼氣的影響［J］.廣東農(nóng)業(yè)科學(xué)，2012，39（18）：190-192.

［13］ 譚新，方熱軍.豬糞對環(huán)境的污染及其處理與利用技術(shù)［J］.飼料工業(yè)，2006，27（21）：58-60.