范學(xué)科，鄭愛泉

(楊凌職業(yè)技術(shù)學(xué)院，陜西 楊凌 712100)

近些年，溫室效應(yīng)所帶來(lái)的環(huán)境問題已備受人們所關(guān)注。溫室氣體的大量排放是造成溫室效應(yīng)和全球變暖的主要原因。為緩解溫室效應(yīng)和全球變暖趨勢(shì)，各個(gè)國(guó)家和地區(qū)均出臺(tái)了一系列溫室氣體減排措施，各個(gè)領(lǐng)域?qū)＜覍W(xué)者也加強(qiáng)了減少溫室氣體排放的相關(guān)研究。清潔能源作為減少溫室氣體排放的重要手段，現(xiàn)已被廣泛研究與使用，其可以明顯地保護(hù)和改善生態(tài)環(huán)境[1]。沼氣作為農(nóng)村主要的清潔能源之一，主要成分甲烷其具有產(chǎn)熱高、燃燒效率和利用率高，并且清潔衛(wèi)生的特點(diǎn)[2,3]。此外，沼氣可以充分利用秸稈、畜禽糞便等農(nóng)業(yè)廢棄物，進(jìn)一步減少了環(huán)境污染[1]。我國(guó)非常重視戶用沼氣池的建設(shè)，到2016年底，戶用沼氣池的規(guī)模已經(jīng)達(dá)到4 193.3萬(wàn)戶[11]。沼氣的大力推廣和應(yīng)用可以有效減少因化石燃料燃燒產(chǎn)生的CO2排放和因畜禽糞便露天堆放產(chǎn)生的CH4排放，對(duì)減少農(nóng)村生活溫室氣體排放及環(huán)境的改善有重要的意義[5]，特別對(duì)于我國(guó)西北部經(jīng)濟(jì)欠發(fā)達(dá)地區(qū)更具有重大意義。此外，沼氣的大規(guī)模使用，還可以有效減少森林砍伐，對(duì)恢復(fù)生態(tài)環(huán)境有重要作用[6～8]。筆者研究以陜西省為研究對(duì)象，通過(guò)文獻(xiàn)調(diào)研和資料收集相結(jié)合的方法，研究未來(lái)陜西省沼氣產(chǎn)量的趨勢(shì)，根據(jù)沼氣減排因子，進(jìn)而估算整個(gè)陜西省溫室氣體減排量。研究將為陜西省戶用沼氣推廣應(yīng)用，促進(jìn)農(nóng)村農(nóng)業(yè)廢棄物資源化，推進(jìn)農(nóng)村能源結(jié)構(gòu)調(diào)整，改善農(nóng)村生產(chǎn)生活條件，促進(jìn)資源節(jié)約型、環(huán)境友好型社會(huì)建設(shè)提供數(shù)據(jù)和理論支撐。

1 材料與方法

1.1 陜西省概況

根據(jù)最新的全國(guó)農(nóng)村沼氣發(fā)展“十三五”規(guī)化，綜合考慮各地區(qū)畜禽糞便、農(nóng)作物秸稈等資源量，肥料化、飼料化、原料化、基料化等競(jìng)爭(zhēng)性利用途徑，以及地域分異規(guī)律、沼氣發(fā)展基礎(chǔ)、經(jīng)濟(jì)水平、清潔能源需求等因素，《規(guī)劃》將全國(guó)31個(gè)省(直轄市、自治區(qū))劃分為三類地區(qū)：I類地區(qū)為資源量豐富地區(qū)，II類地區(qū)為資源量中等地區(qū)，III類地區(qū)為資源量一般地區(qū))。陜西省屬于II類地區(qū)，并且屬于“鐮刀彎”地區(qū)，是玉米結(jié)構(gòu)調(diào)整的重點(diǎn)地區(qū)，也是草食動(dòng)物養(yǎng)殖優(yōu)勢(shì)區(qū)，適宜發(fā)展以規(guī)?；託鉃榧~帶的循環(huán)農(nóng)業(yè)模式，適度發(fā)展生物天然氣工程和中小型沼氣工程[9]。

1.2 回歸模型

為了準(zhǔn)確研究沼氣使用對(duì)溫室氣體排放的影響，研究采用文獻(xiàn)調(diào)研和資料收集相結(jié)合的方法獲取有關(guān)數(shù)據(jù)：期刊論文，中國(guó)統(tǒng)計(jì)年鑒[10]、中國(guó)農(nóng)業(yè)年鑒[11]等來(lái)源搜集基礎(chǔ)數(shù)據(jù)。通過(guò)基礎(chǔ)數(shù)據(jù)的綜合分析，運(yùn)用多變量線性回歸模型和逐步回歸模型對(duì)陜西省沼氣產(chǎn)量進(jìn)行預(yù)測(cè)，進(jìn)而根據(jù)未來(lái)沼氣產(chǎn)量估算出陜西省“養(yǎng)殖-沼氣”生態(tài)模式溫室氣體減排潛力。

1.2.1 變量的選取 年沼氣產(chǎn)量是衡量和評(píng)價(jià)農(nóng)村戶用沼氣發(fā)展水平的關(guān)鍵指標(biāo)。因此，將陜西省當(dāng)年沼氣產(chǎn)量作為因變量，將搜集到的人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄、沼氣池?cái)?shù)量分別作為自變量x1～x7。

1.2.2 模型適配度檢驗(yàn) 在回歸分析中，如果有兩個(gè)或兩個(gè)以上的自變量，就稱為多元回歸。事實(shí)上，一種現(xiàn)象常常是與多個(gè)因素相聯(lián)系的，由多個(gè)自變量的最優(yōu)組合共同來(lái)預(yù)測(cè)或估計(jì)因變量，比只用一個(gè)自變量進(jìn)行預(yù)測(cè)或估計(jì)更有效，更符合實(shí)際。

逐步回歸是一種線性回歸模型自變量選擇方法，其基本思想是將變量一個(gè)一個(gè)引入，引入的條件使其偏回歸平方和經(jīng)驗(yàn)是顯著的。同時(shí)，每引入一個(gè)新變量后，對(duì)已入選回歸模型的老變量逐個(gè)進(jìn)行檢驗(yàn)，將經(jīng)檢驗(yàn)認(rèn)為不顯著的變量刪除，以保證所得自變量子集中每一個(gè)變量都是顯著的。此過(guò)程經(jīng)過(guò)若干步直到不能再引入新變量為止。這時(shí)回歸模型中所有變量對(duì)因變量都是顯著的。

基于沼氣產(chǎn)量預(yù)測(cè)受多個(gè)因子的影響，筆者研究選取多元線性回歸模型和逐步回歸模型進(jìn)行沼氣產(chǎn)量預(yù)測(cè)。從圖1可知，多元線性回歸和逐步回歸的標(biāo)準(zhǔn)化殘差的絕對(duì)值相差不大，且均接近1，說(shuō)明二者適配性一致，均可用來(lái)預(yù)測(cè)沼氣產(chǎn)量。

圖1 回歸標(biāo)準(zhǔn)化殘差直方圖

從表1可以看出，多元線性回歸模型的復(fù)相關(guān)系數(shù)為0.992，判定系數(shù)為0.984，校正后R2為0.972，Durbin-Watson統(tǒng)計(jì)量數(shù)值為2.096，F(xiàn)測(cè)驗(yàn)值86.20(p<0.001)。逐步回歸模型的復(fù)相關(guān)系數(shù)為0.986，判定系數(shù)為0.972，校正后R2為0.969，達(dá)沃檢驗(yàn)統(tǒng)計(jì)量數(shù)值為1.594，F(xiàn)測(cè)驗(yàn)值264.10(p<0.001)。從模型適配參數(shù)和殘差分布看，兩個(gè)模型適配度均較優(yōu)，可以用來(lái)后續(xù)預(yù)測(cè)分析。

表1 回歸模型擬合適配度

1.3 減排潛力估算

根據(jù)李長(zhǎng)江等[11]關(guān)于陜西省農(nóng)村CDM項(xiàng)目溫室氣體減排效應(yīng)研究發(fā)現(xiàn)，沼氣項(xiàng)目戶均溫室氣體減排量在1.77～2.46 tCO2e。根據(jù)該研究調(diào)研結(jié)果，沼氣項(xiàng)目戶均沼氣產(chǎn)量為257.84～543.0 m3?；诖?，可估算得到每立方米沼氣可以減排3.25～9.54 kg CO2e溫室氣體。因此，整個(gè)陜西省“養(yǎng)殖-沼氣”生態(tài)模式溫室氣體減排潛力按照以下公式估算：

GS=Yi×F

式中：Gs代表當(dāng)年溫室氣體總的減排量，以二氧化碳當(dāng)量計(jì)，t·a-1CO2e；Yi代表當(dāng)年沼氣預(yù)測(cè)產(chǎn)量，單位m3；代表每產(chǎn)生1 m3沼氣所帶來(lái)的減排效應(yīng)，單位t·m3CO2e。

1.4 經(jīng)濟(jì)效益計(jì)算

沼氣的使用能夠在很大程度上減少原煤的消耗。國(guó)際上為了使用的方便，統(tǒng)一標(biāo)準(zhǔn)，在進(jìn)行能源數(shù)量、質(zhì)量的比較時(shí)，將煤炭、石油、天然氣等都按一定的比例統(tǒng)一換算成標(biāo)準(zhǔn)煤來(lái)表示?；诖耍P者研究按照燃料轉(zhuǎn)換系數(shù)，將產(chǎn)生的沼氣轉(zhuǎn)換為原煤量，按照當(dāng)年原煤售價(jià)，估算沼氣利用節(jié)約的原煤花銷。具體估算公式如下：

BS=Yi×Fi×Fj×S

式中：BS代表經(jīng)濟(jì)效益，元；Yi代表當(dāng)年沼氣預(yù)測(cè)產(chǎn)量，單位m3；Fi代表沼氣和標(biāo)準(zhǔn)煤轉(zhuǎn)化系數(shù)，筆者研究采用0.714[15]；Fj代表標(biāo)準(zhǔn)煤和原煤轉(zhuǎn)化系數(shù)，研究采用1.400[15]；S代表當(dāng)年原煤價(jià)格,研究以1 360 t·元-1。

1.5 數(shù)據(jù)處理

采用Microsoft Excel 365和SPSS軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，采用Origin 2016進(jìn)行數(shù)據(jù)可視化。

2 結(jié)果

2.1 陜西省沼氣產(chǎn)量及社會(huì)經(jīng)濟(jì)發(fā)展

研究主要參考中國(guó)農(nóng)業(yè)統(tǒng)計(jì)年鑒和中國(guó)統(tǒng)計(jì)年鑒得到2000-2017年的陜西省沼氣產(chǎn)量及社會(huì)經(jīng)濟(jì)發(fā)展相關(guān)數(shù)據(jù)，具體如表2所示。從表2可以看出，自2000年以來(lái)，陜西省沼氣年實(shí)際產(chǎn)量和沼氣池?cái)?shù)量逐年提升，近年趨緩。豬、牛、羊、家禽也表現(xiàn)出相同趨勢(shì)。整個(gè)陜西省人均生產(chǎn)總值和人均農(nóng)林牧漁總產(chǎn)值也逐年提升。根據(jù)計(jì)算，2017年沼氣產(chǎn)量是2000年的7.29倍，沼氣池?cái)?shù)量也由14.2萬(wàn)口增加到133萬(wàn)口，增幅836.62%。人均生產(chǎn)總值增加了1 052.70%，人均農(nóng)林牧漁總產(chǎn)值增加了527.34%，豬出欄量增加了44.21%，牛出欄量減少了-2.97%，羊出欄量增加了23.38%，家禽出欄量增加了13.76%。

表2 2000-2017陜西省沼氣產(chǎn)量及社會(huì)經(jīng)濟(jì)發(fā)展相關(guān)數(shù)據(jù)

2.2 回歸預(yù)測(cè)方程

通過(guò)多元線性回歸，將當(dāng)年沼氣產(chǎn)量作為因變量，影響沼氣產(chǎn)量的相關(guān)因素作為自變量，根據(jù)擬合系數(shù)建立的多元線性回歸方程為：

y1=-1 560.384+0.08x1-5.526x2+9.164x3-

153.026x4+31.018x5-0.750x6+388.030x7

式中：x1代表人均生產(chǎn)總值(元)；x2代表人均農(nóng)林牧漁總產(chǎn)值(元)；x3代表當(dāng)年豬出欄量(萬(wàn)頭)；x4代表當(dāng)年牛出欄量(萬(wàn)頭)；x5代表當(dāng)年羊出欄量(萬(wàn)頭)；x6代表當(dāng)年家禽出欄量(萬(wàn)只)；x7代表當(dāng)年沼氣池?cái)?shù)量(萬(wàn)口)。

通過(guò)逐步回歸，在擬合過(guò)程中模型自動(dòng)剔除了當(dāng)年人均生產(chǎn)總值(元，x1)，當(dāng)年豬出欄量(萬(wàn)頭，x3)，當(dāng)年牛出欄量(萬(wàn)頭，x4)，當(dāng)年羊出欄量(萬(wàn)頭，x5)，當(dāng)年家禽出欄量(萬(wàn)只，x6)將當(dāng)年沼氣產(chǎn)量作為因變量，根據(jù)擬合系數(shù)建立的逐步回歸方程為：

y2=3 575.111+456.681x7-5.084x2

式中：x2代表人均農(nóng)林牧漁總產(chǎn)值(元)；x7代表當(dāng)年沼氣池?cái)?shù)量(萬(wàn)口)。

2.3 回歸模型沼氣產(chǎn)量預(yù)測(cè)分析

為進(jìn)一步檢驗(yàn)?zāi)Ｐ蛿M合的優(yōu)劣，筆者研究通過(guò)2000-2017年沼氣預(yù)測(cè)產(chǎn)量和實(shí)際產(chǎn)量的對(duì)比發(fā)現(xiàn)，多元線性回歸模型和逐步回歸模型與實(shí)際沼氣產(chǎn)量的差別較小，有較高的預(yù)測(cè)優(yōu)度(圖2)。

圖2 沼氣產(chǎn)量模型預(yù)測(cè)情況

2.4 未來(lái)沼氣產(chǎn)量及自變量預(yù)測(cè)值

通過(guò)分析發(fā)現(xiàn)，人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄和沼氣池?cái)?shù)量在發(fā)展趨勢(shì)上均較穩(wěn)定，且有較強(qiáng)的變化規(guī)律，可以通過(guò)自變量回歸擬合進(jìn)行趨勢(shì)分析?；诖?，筆者研究通過(guò)自變量回歸分析，對(duì)當(dāng)年人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄、沼氣池?cái)?shù)量等自變量分別使用線性函數(shù)、S型函數(shù)等方法，分別進(jìn)行擬合求值，得到2020、2030、2040和2050年的自變量值，具體見表3。根據(jù)擬合結(jié)果，與2017年相比，2020年人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄、沼氣池?cái)?shù)量分別增長(zhǎng)7.13%、13.48%、16.34%、4.12%、9.05%、-1.29%和36.58%；2030年人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄、沼氣池?cái)?shù)量分別增長(zhǎng)61.14%、66.55%、42.18%、19.69%、39.66%、-20.30%和102.39%；2040年人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄、沼氣池?cái)?shù)量分別增長(zhǎng)115.12%、119.57%、68.03%、35.26%、70.27%、7.78%和168.20%；2050年人均生產(chǎn)總值、人均農(nóng)林牧漁總產(chǎn)值、豬出欄量、牛出欄、羊出欄、家禽出欄、沼氣池?cái)?shù)量分別增長(zhǎng)169.07%、169.99%、93.87%、50.84%、100.89%、-4.86%和234.02%。

表3 社會(huì)經(jīng)濟(jì)因子預(yù)測(cè)

在此基礎(chǔ)上，運(yùn)用多元線性回歸模型和逐步回歸模型對(duì)未來(lái)沼氣產(chǎn)量進(jìn)行了預(yù)測(cè)分析，未來(lái)沼氣預(yù)測(cè)產(chǎn)量見圖 3。多元線性回歸模型結(jié)果表明，與2017年相比，2020年沼氣產(chǎn)量達(dá)4.00×108m3, 增幅66.97%；2030年沼氣產(chǎn)量達(dá)5.99×108m3, 增幅149.69%；2040年沼氣產(chǎn)量達(dá)7.79×108m3, 增幅224.72%；2050年沼氣產(chǎn)量達(dá)9.86×108m3, 增幅311.28%。逐步回歸模型結(jié)果表明，與2017年相比，2020年沼氣產(chǎn)量達(dá)4.02×108m3, 增幅71.17%；2030年沼氣產(chǎn)量達(dá)5.85×108m3, 增幅149.10%；2040年沼氣產(chǎn)量達(dá)7.68×108m3, 增幅227.13%；2050年沼氣產(chǎn)量達(dá)9.62×108m3, 增幅309.67%。

圖3 未來(lái)沼氣產(chǎn)量預(yù)測(cè)數(shù)值

2.5 減排潛力和經(jīng)濟(jì)效益估算

根據(jù)回歸模型預(yù)測(cè)結(jié)果，隨著戶用沼氣在整個(gè)陜西省的推廣應(yīng)用，能夠顯著減少溫室氣體減排量，并且?guī)?lái)可觀的經(jīng)濟(jì)效益(表4)。通過(guò)預(yù)測(cè)沼氣產(chǎn)量，估算陜西省戶用沼氣總減排量為1.30×106～9.41×106tCO2e。并且戶用沼氣使用后，節(jié)省了煤炭的使用，降低了農(nóng)戶的經(jīng)濟(jì)支出，為農(nóng)戶獲得經(jīng)濟(jì)收入。因此，戶用沼氣使用后，陜西農(nóng)戶可獲得的經(jīng)濟(jì)效益為煤炭節(jié)省費(fèi)用。通過(guò)原煤轉(zhuǎn)換發(fā)現(xiàn)陜西省每年減少煤炭消耗4.0×105～9.9×105t，采用該數(shù)據(jù)結(jié)合當(dāng)年當(dāng)?shù)孛繃嵜禾? 360元(來(lái)自中國(guó)碳排放交易網(wǎng))，估算陜西省戶用沼氣使用后所獲得的經(jīng)濟(jì)效益為54 419.6×104～134 045.1×104元。

表4 關(guān)中地區(qū)戶用沼氣減排潛力及經(jīng)濟(jì)效益估算

3 討論與結(jié)論

陜西省“養(yǎng)殖-沼氣”生態(tài)模式可以實(shí)現(xiàn)種植業(yè)和養(yǎng)殖業(yè)的綠色生態(tài)發(fā)展，是減少溫室氣體排放，減少農(nóng)村人口生活成本的有效措施，同時(shí)也是實(shí)現(xiàn)“美麗鄉(xiāng)村”的重要途徑[6,7,8,17,16]。筆者研究依據(jù)2000-2017 年陜西農(nóng)村沼氣發(fā)展的相關(guān)數(shù)據(jù), 通過(guò)模型擬合預(yù)測(cè)陜西省2020，2030，2040和2050年沼氣產(chǎn)量，并通過(guò)溫室氣體減排量估算表明，陜西省大力發(fā)展農(nóng)村戶用沼氣，是實(shí)現(xiàn)農(nóng)村廢棄資源再利用，改善鄉(xiāng)村環(huán)境的重要途徑。研究發(fā)現(xiàn)，隨著人均農(nóng)林牧漁總產(chǎn)值和豬、牛、羊、家禽出欄數(shù)量的持續(xù)增加，沼氣池?cái)?shù)量和沼氣年產(chǎn)量也同步增加。主要是由于農(nóng)村沼氣的原料主要是畜禽糞便和秸稈[15]。隨著農(nóng)業(yè)秸稈的增多，刺激畜牧業(yè)快速發(fā)展，畜禽糞便是優(yōu)良的沼氣發(fā)酵原料，進(jìn)而促進(jìn)了沼氣產(chǎn)業(yè)的發(fā)展。此外，陜西區(qū)域特點(diǎn)十分明顯，有利于農(nóng)村戶用沼氣的發(fā)展[18]。關(guān)中地區(qū)玉米的大范圍種植和大型養(yǎng)殖廠為戶用沼氣發(fā)展提供了較為豐富的發(fā)酵原料；陜南目前在大力發(fā)展蠶和漁業(yè)，可以將沼氣池殘?jiān)?jīng)過(guò)無(wú)害化處理后支持蠶、漁的發(fā)展，解決了沼氣殘?jiān)幚韱栴}；陜北通過(guò)冬季沼氣池保溫技術(shù)也解決了發(fā)展戶用沼氣的瓶頸。

前人研究表明，農(nóng)村戶用沼氣能夠減少溫室氣體排放主要是由于改變了農(nóng)戶能源結(jié)構(gòu)，減少了煤炭、薪柴及秸稈的使用，可以保護(hù)林地、減少農(nóng)民對(duì)林地的砍伐，并且畜禽糞便也有室外堆放管理方式改變?yōu)樽鳛檎託庠线M(jìn)行集中處理，進(jìn)而減少農(nóng)戶溫室氣體的排放[17,18]。此外，農(nóng)戶傳統(tǒng)能源在燃燒時(shí)會(huì)釋放大量PM 2.5，能影響人和牲畜的健康及農(nóng)村環(huán)境的保護(hù)，沼氣充分燃燒的產(chǎn)物為CO2和水，能夠有效抑制PM2.5的產(chǎn)生[21,22]。筆者研究發(fā)現(xiàn)，戶用能源結(jié)構(gòu)改變減少了溫室氣體的排放，陜西省每年溫室氣體減排分別為1.30×106～9.18×106tCO2e；沼氣使用后，因節(jié)約原煤而產(chǎn)生的經(jīng)濟(jì)效益為5.44×108～1.34×109元。目前關(guān)于沼氣、沼渣、沼液的研究表明，秸稈和糞便經(jīng)厭氧發(fā)酵后的殘液、殘?jiān)且环N無(wú)公害的農(nóng)家肥，可以有效殺死部分害蟲，并能夠提供部分作物所需養(yǎng)分，從而顯著減少農(nóng)業(yè)投入，減少溫室氣體排放，并對(duì)農(nóng)村農(nóng)產(chǎn)品的綠色生產(chǎn)也有重要意義[20]。從能源結(jié)構(gòu)轉(zhuǎn)變來(lái)看，戶用沼氣的大量推廣，可以有效減少原煤的使用，同時(shí)能夠減少因原煤開采附帶的污染和環(huán)境破壞，且能夠減少購(gòu)買原煤的生活支出等[13]。因此，未來(lái)需要進(jìn)一步完善沼氣使用產(chǎn)生減排的來(lái)源及相關(guān)的估算方法。以上結(jié)果也表明戶用沼氣應(yīng)用在農(nóng)村中的重要意義。

針對(duì)目前陜西省沼氣發(fā)展現(xiàn)狀，本研究對(duì)該區(qū)域加快沼氣發(fā)展提出以下對(duì)策建議：①加速大型沼氣池建設(shè)。目前全省農(nóng)業(yè)合作社蓬勃興起，沼氣發(fā)展可以農(nóng)業(yè)合作社為龍頭，建設(shè)規(guī)模較大沼氣池，可以提升產(chǎn)氣量和廢棄資源再利用效率；②全省沼氣發(fā)展應(yīng)該分區(qū)域規(guī)化。沼氣發(fā)展受到多因素影響，可以依據(jù)氣候和現(xiàn)有資源條件，對(duì)最適宜區(qū)域重點(diǎn)發(fā)展，資源薄弱區(qū)域強(qiáng)度減緩，使效益最大化；③不斷改進(jìn)沼氣技術(shù)。目前陜西省沼氣發(fā)展主要以畜禽糞便發(fā)酵為主，后期應(yīng)該加大以秸稈發(fā)酵為主的技術(shù)的開發(fā)利用，使作物殘?bào)w有效利用。