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

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

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

中圖分類號(hào)X 713文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)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

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

1材料與方法

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

1.2試驗(yàn)儀器試驗(yàn)所用儀器具體見表2。

1.3試驗(yàn)裝置

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

1.4試驗(yàn)設(shè)計(jì)

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

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

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

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

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

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

2結(jié)果與分析

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

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

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

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

利用Design Expert 8.0.6軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行分析擬合，并進(jìn)行差異顯著性分析。平均日產(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，說明回歸模型極顯著。失擬項(xiàng)P值大于0.05，影響不顯著，說明該方程擬合度良好。通過分析各因素的一次項(xiàng)、二次項(xiàng)和交互項(xiàng)的P值可知，一次項(xiàng)因素A、B、C的P值均小于0.05，影響顯著；二次項(xiàng)因素A2、B2、C2的P值均小于0.05，影響顯著；交互項(xiàng)AB的P值小于0.05，影響顯著；交互項(xiàng)BC的P值大于0.05，說明影響不顯著。通過對(duì)F值的分析發(fā)現(xiàn)，各影響因素的主次順序?yàn)锳>B>C，即HRT>攪拌間隔>攪拌轉(zhuǎn)速。

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

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

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

（2）對(duì)VS降解率的影響分析。

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

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

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

2.3最佳條件的確定及試驗(yàn)驗(yàn)證

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

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

3結(jié)論

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

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

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

參考文獻(xiàn)

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

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

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

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

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

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

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

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

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

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

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

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

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