董法秀，萬大娟，王開心，蘭鵬鵬

(湖南師范大學 資源與環(huán)境科學學院，湖南 長沙 410081)

隨著我國經(jīng)濟社會的快速發(fā)展和人民消費水平的不斷提高，國民對畜產(chǎn)品的需求量呈不斷增長趨勢，從而促進了畜禽養(yǎng)殖業(yè)快速發(fā)展[1]。畜禽養(yǎng)殖業(yè)的發(fā)展導(dǎo)致大量污染物產(chǎn)生，現(xiàn)已成為生態(tài)環(huán)境的重要污染源[2]。環(huán)境承載力是衡量經(jīng)濟與環(huán)境協(xié)調(diào)發(fā)展程度的指標[3]。目前，國內(nèi)有很多關(guān)于畜禽養(yǎng)殖環(huán)境承載力的研究，如胡雪飆[2]、王延吉[4]、許翼[5]等分別對重慶市、延邊地區(qū)、沈陽市的畜禽養(yǎng)殖環(huán)境承載力進行了相關(guān)研究。本文以長沙市畜禽養(yǎng)殖為研究對象，利用灰色系統(tǒng)模型對長沙市畜禽養(yǎng)殖環(huán)境承載力進行了模型評價和分析預(yù)測，以期為該地區(qū)畜禽養(yǎng)殖業(yè)的污染治理及規(guī)劃布局提供科學依據(jù)。

1 材料與方法

1.1 指標體系的構(gòu)建

本文在各項指標的選擇上以實用、簡明、反映系統(tǒng)本質(zhì)、易于量化為原則[6]。結(jié)合長沙市現(xiàn)狀，指標體系主要由牧業(yè)總產(chǎn)值、地表水資源量、耕地面積、COD產(chǎn)生量、地下水資源量、TN產(chǎn)生量、BOD產(chǎn)生量、TP產(chǎn)生量共8項指標構(gòu)成。

畜禽養(yǎng)殖環(huán)境承載力指標可劃分為兩大類：一是社會、人口、經(jīng)濟活動方面的發(fā)展類指標，主要包括牧業(yè)總產(chǎn)值、BOD產(chǎn)生量、TN產(chǎn)生量、COD產(chǎn)生量和TP產(chǎn)生量，描述社會經(jīng)濟發(fā)展的規(guī)模大小；二是與畜禽養(yǎng)殖系統(tǒng)環(huán)境有關(guān)的變量，即限制類指標，用來表征畜禽養(yǎng)殖系統(tǒng)的質(zhì)量、狀態(tài)、發(fā)展方向等，主要包括地表水資源量、耕地面積、地下水資源量[7]。

1.2 模型的選取

1.2.1 系統(tǒng)分析模型 根據(jù)國內(nèi)外環(huán)境承載力評價模型研究狀況，以及長沙市畜禽養(yǎng)殖發(fā)展現(xiàn)狀和原始數(shù)據(jù)的可獲取性，本文選擇系統(tǒng)分析模型作為長沙市畜禽養(yǎng)殖環(huán)境承載力評價模型。

在畜禽養(yǎng)殖環(huán)境承載力的研究中，針對n個指標給出了不同年份的畜禽養(yǎng)殖環(huán)境承載力分量，假設(shè)j個年份環(huán)境承載力為bj(j=1,2,3,…,m)，如果j個年份環(huán)境承載力又由n個具體指標所確定的分量組成，即：Bij=(B1j,B2j,B3j,…,Bnj)，Bij代表第i個指標第j個年份下具體指標值，對限制類指標進行規(guī)格化處理后，即：bij=(b1j,b2j,b3j,…,bnj)[4]。具體計算公式如下：

限制類指標：

(1)

發(fā)展類指標：

(2)

這樣，第j年環(huán)境承載力的大小可以用歸一化的矢量模型來表示，即：

(3)

1.2.2 灰色數(shù)列預(yù)測模型 畜禽養(yǎng)殖業(yè)的發(fā)展具有持續(xù)發(fā)展性，本文運用灰色數(shù)列預(yù)測模型建立預(yù)測方程，并利用該方程對長沙市未來各年份的畜禽養(yǎng)殖環(huán)境承載力進行預(yù)測分析。

灰色數(shù)列預(yù)測模型是通過利用數(shù)據(jù)之間的內(nèi)在關(guān)系，對其中一個指標的發(fā)展變化情況做出預(yù)測，預(yù)測結(jié)果表示該指標在將來發(fā)展變化的各個時間段的具體數(shù)值。數(shù)列預(yù)測是在GM(1,1)模型的基礎(chǔ)上進行累加生成數(shù)列[8]。

1.3 數(shù)據(jù)的收集與處理

1.3.1 基礎(chǔ)數(shù)據(jù) 2010～2016年長沙市畜禽養(yǎng)殖基礎(chǔ)數(shù)據(jù)和耕地數(shù)據(jù)來自2011～2017年的《長沙市統(tǒng)計年鑒》(表1)。我國豬、家禽的平均飼養(yǎng)周期都在一年以內(nèi)，故以年出欄量為飼養(yǎng)量；牛、羊的飼養(yǎng)周期都在一年以上，故以年存欄量為飼養(yǎng)量。

1.3.2 數(shù)據(jù)處理 長沙市2010～2016年地表水資源量和地下水資源量數(shù)據(jù)來自《長沙市水資源公報》；根據(jù)《畜禽養(yǎng)殖業(yè)污染治理工程技術(shù)規(guī)范》(HJ 497─2009)及相關(guān)參考文獻中各類畜禽排泄系數(shù)，結(jié)合長沙市畜禽的實際生長周期，對文獻中的排泄系數(shù)加以修正，計算得出各污染物產(chǎn)生量[9-10]；2010～2016年長沙市畜禽糞尿中污染物產(chǎn)生量、地表水資源量、地下水資源量及牧業(yè)總產(chǎn)值的適宜值和警戒值分別取自2010～2016年的平均值和最小值；其中耕地面積取自長沙市規(guī)劃2020年全市耕地保有量26.90萬hm2；BOD、TN、COD、TP污染物產(chǎn)生量適宜值和警戒值依據(jù)：適宜值取15生豬當量(hm2)、警戒值取25生豬當量(hm2)計算而得(表2)。

表2 長沙市畜禽養(yǎng)殖評價指標體系的基本數(shù)據(jù)

2 結(jié)果與分析

2.1 長沙市畜禽養(yǎng)殖現(xiàn)狀

依據(jù)長沙市畜禽養(yǎng)殖統(tǒng)計數(shù)據(jù)分析，長沙市生豬養(yǎng)殖量最大，2013年生豬出欄量達到835.68萬頭；其次是家禽，2015年家禽存籠量達到6005.75萬羽。從空間分布情況來看，寧鄉(xiāng)縣畜禽養(yǎng)殖數(shù)量最多，瀏陽市次之，芙蓉區(qū)、岳麓區(qū)、天心區(qū)、雨花區(qū)及開福區(qū)已禁止規(guī)模化畜禽養(yǎng)殖。

2.2 長沙市畜禽養(yǎng)殖環(huán)境承載力現(xiàn)狀

根據(jù)長沙市畜禽養(yǎng)殖數(shù)量、地表水資源量、畜禽養(yǎng)殖排泄物等基礎(chǔ)數(shù)據(jù)，應(yīng)用系統(tǒng)分析方法計算該地區(qū)2010～2016年畜禽養(yǎng)殖環(huán)境承載力R(無量綱)(表3)。長沙市環(huán)境承載力R適宜值為0.372，警戒值為0.271，環(huán)境承載力R值越大，表示畜禽養(yǎng)殖環(huán)境承載力越強。當R≥0.372時，表明該年份環(huán)境承載能力較強，當前的畜禽養(yǎng)殖狀況較好；當0.271≤R<0.372時，表明該年份畜禽養(yǎng)殖環(huán)境承載能力較弱，在當前的畜禽養(yǎng)殖狀況下還可進行持續(xù)發(fā)展，但發(fā)展的同時要注意畜禽養(yǎng)殖產(chǎn)生的環(huán)境問題；當R<0.271時，表明該年份畜禽養(yǎng)殖環(huán)境承載力很弱，在合理控制畜禽養(yǎng)殖業(yè)發(fā)展的同時要加強環(huán)境保護。2010～2016年長沙市畜禽養(yǎng)殖環(huán)境承載力依次為0.327、0.291、0.325、0.290、0.320、0.312、0.344。由圖1可以看出，長沙市畜禽養(yǎng)殖環(huán)境承載力值處于警戒值與適宜值之間，說明目前長沙市畜禽養(yǎng)殖環(huán)境承載力較弱。

表3 長沙市不同年份畜禽養(yǎng)殖環(huán)境承載力分量及綜合值(無量綱)

圖1 長沙市畜禽養(yǎng)殖環(huán)境承載力的變化趨勢

2.3 長沙市畜禽養(yǎng)殖環(huán)境承載力預(yù)測分析

2.3.1 自然狀態(tài)下長沙市畜禽養(yǎng)殖環(huán)境承載力預(yù)測分析 利用GM(1，1)灰色數(shù)列預(yù)測模型對自然狀態(tài)下的長沙市畜禽養(yǎng)殖環(huán)境承載力進行預(yù)測。

按照灰色預(yù)測模型方法，以X(0)={0.327，0.291，0.325，0.290，0.320，0.312，0.344}為原始數(shù)據(jù)數(shù)列進行計算，最后得出時間響應(yīng)函數(shù)：

X(1)(K+1)=12.3741e0.0236K-12.0471

(4)

將數(shù)據(jù)代入式(4)中，計算得到長沙市畜禽養(yǎng)殖環(huán)境承載力的實際值與預(yù)測值(表4)。通過檢驗，生成數(shù)列的相對誤差最大值為1.902%，最小值為-0.013%；還原數(shù)列的相對誤差最大值為6.905%，最小值為0.881%。因此，可認為該地區(qū)的預(yù)測值較真實，預(yù)測方法是合理有效的。

表4 生成數(shù)列的還原數(shù)列誤差檢驗(自然狀態(tài))

通過時間響應(yīng)函數(shù)方程預(yù)測長沙市2017～2025年9年畜禽養(yǎng)殖環(huán)境承載力(表5)，結(jié)果表明：長沙市未來9年畜禽養(yǎng)殖環(huán)境承載力值R呈上升趨勢，到2021年畜禽養(yǎng)殖環(huán)境承載力值為0.374，大于適宜值。這說明在自然狀態(tài)下未來長沙市畜禽養(yǎng)殖業(yè)發(fā)展良好，但在發(fā)展的同時，應(yīng)盡可能做到對畜禽糞便的減量化及再利用，減少對環(huán)境的影響。

2.3.2 基于畜禽養(yǎng)殖污染物處理下的長沙市畜禽養(yǎng)殖環(huán)境承載力預(yù)測分析 從上述研究結(jié)果中發(fā)現(xiàn)，雖然長沙市畜禽養(yǎng)殖環(huán)境承載力呈上升趨勢發(fā)展，但到2021年長沙市的畜禽養(yǎng)殖環(huán)境承載力才超過適宜值。如果將畜禽養(yǎng)殖污染物采取一定的處理方法，將其轉(zhuǎn)換為肥料、能源等，那將會很大程度地減小對環(huán)境造成的壓力。

根據(jù)長沙市畜禽養(yǎng)殖數(shù)量、耕地面積、地表水資源量等自然資源情況，若將70%的畜禽養(yǎng)殖污染物轉(zhuǎn)變成化肥、能源等產(chǎn)物，采用上述相同方法計算出長沙市畜禽養(yǎng)殖環(huán)境承載力值R(表6)，2010～2016年間該地區(qū)畜禽養(yǎng)殖環(huán)境承載力適宜值為0.248，警戒值為0.214，R值均超過適宜值。

表5 長沙市畜禽養(yǎng)殖環(huán)境承載力綜合預(yù)測

運用灰色預(yù)測模型方法計算得到長沙市畜禽養(yǎng)殖環(huán)境承載力預(yù)測值(表7)，通過檢驗，生成數(shù)列的相對誤差最大值為1.724%；還原數(shù)列最大值為5.3256%。因此，可認為該地區(qū)的預(yù)測值較真實，預(yù)測方法是合理有效的。通過預(yù)測結(jié)果可知，若將長沙市70%的畜禽養(yǎng)殖污染物再利用，2017～2025年長沙市畜禽養(yǎng)殖環(huán)境承載力均超過適宜值，說明該地區(qū)環(huán)境承載能力較強，未來畜禽養(yǎng)殖業(yè)發(fā)展前景樂觀。

3 結(jié)論

本研究通過選取長沙市地表水資源量、牧業(yè)總產(chǎn)值、COD(化學需氧量)產(chǎn)生量、TN(總氮)產(chǎn)生量等8項指標構(gòu)建成長沙市畜禽養(yǎng)殖環(huán)境承載力評價指標體系，建立系統(tǒng)分析模型。數(shù)據(jù)結(jié)果顯示2010～2016年長沙市畜禽養(yǎng)殖環(huán)境承載力處于警戒值與適宜值之間，說明該地區(qū)畜禽養(yǎng)殖環(huán)境承載能力相對較弱。

表6 長沙市不同年份畜禽養(yǎng)殖環(huán)境承載力分量及綜合值(無量綱，處理狀態(tài))

表7 長沙市畜禽養(yǎng)殖環(huán)境承載力綜合預(yù)測(處理狀態(tài))

通過利用灰色數(shù)列預(yù)測模型對長沙市畜禽養(yǎng)殖環(huán)境承載力的預(yù)測模擬，可以看出在自然狀態(tài)下2017～2025年長沙市畜禽養(yǎng)殖環(huán)境承載力值呈上升發(fā)展趨勢，到2021年畜禽養(yǎng)殖環(huán)境承載力值大于適宜值0.372，說明長沙市未來9年畜牧業(yè)發(fā)展情況良好，控制在環(huán)境承載力范圍之內(nèi)。

通過模型預(yù)測結(jié)果可知，若長沙市70%的畜禽養(yǎng)殖污染物再利用，轉(zhuǎn)換成有機肥、能源等，未來9年畜禽養(yǎng)殖環(huán)境承載力均超過適宜值，長沙市畜禽養(yǎng)殖環(huán)境承載力將更強，畜禽養(yǎng)殖業(yè)良性發(fā)展更持久。