中國(guó)不同規(guī)模奶牛養(yǎng)殖環(huán)境效率及其影響因素分析

王善高，許 昭

（南京農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，江蘇南京 210095）

伴隨著需求增長(zhǎng)，我國(guó)奶業(yè)發(fā)展迅速，奶牛養(yǎng)殖規(guī)?；?、專(zhuān)業(yè)化程度不斷提升。統(tǒng)計(jì)數(shù)據(jù)顯示，2007—2015年，我國(guó)奶牛養(yǎng)殖規(guī)?；潭扔?6.6%增長(zhǎng)至70.5%[1]。規(guī)?；B(yǎng)殖在節(jié)約成本、增加產(chǎn)量、保障供給等方面均具有比較優(yōu)勢(shì)，但其帶來(lái)的環(huán)境污染問(wèn)題也日益凸顯[2]。奶牛是僅次于生豬、肉牛的第三大產(chǎn)糞畜種，其排放的糞污富含化學(xué)需氧量、氮、磷等有機(jī)物質(zhì)，而這些有機(jī)物質(zhì)的流失不僅會(huì)造成水體富營(yíng)養(yǎng)化[3]，還會(huì)污染空氣和土壤，最終導(dǎo)致奶牛養(yǎng)殖的不可持續(xù)。環(huán)境效率反映了在產(chǎn)出和投入不變的情況下，污染物可以減少的程度，測(cè)度環(huán)境效率不僅可以了解畜禽養(yǎng)殖污染現(xiàn)狀，還可以測(cè)算出污染物的潛在削減量。在畜禽養(yǎng)殖污染日益嚴(yán)重的背景下，科學(xué)地測(cè)度出我國(guó)奶牛養(yǎng)殖的環(huán)境效率并計(jì)算出奶牛養(yǎng)殖污染物的潛在削減量顯得尤為重要。

關(guān)于畜禽養(yǎng)殖環(huán)境效率的測(cè)算，國(guó)內(nèi)外學(xué)者均進(jìn)行了大量研究。國(guó)外學(xué)者研究的相對(duì)較早，如Reinhard等[4]、Lansink等[5]和Yang[6]分別測(cè)算了荷蘭奶牛、荷蘭生豬、中國(guó)臺(tái)灣生豬的環(huán)境效率。隨著我國(guó)畜禽養(yǎng)殖環(huán)境污染問(wèn)題日益突出，國(guó)內(nèi)學(xué)者也開(kāi)始關(guān)注畜禽養(yǎng)殖環(huán)境效率，并針對(duì)生豬、蛋雞、奶牛等畜禽進(jìn)行了研究[7-9]。然而，梳理國(guó)內(nèi)外文獻(xiàn)發(fā)現(xiàn)，現(xiàn)有研究存在一些不足：第一，大多數(shù)研究?jī)H僅停留在測(cè)算環(huán)境效率表面，鮮有學(xué)者基于環(huán)境效率值做進(jìn)一步研究。第二，污染物的界定并不統(tǒng)一，使得測(cè)算出的環(huán)境效率難以真實(shí)地反映畜禽養(yǎng)殖環(huán)境現(xiàn)狀。第三，研究方法上大多采用DEA方法和SBM模型，這2種方法都是基于線性規(guī)劃方法求解效率，會(huì)導(dǎo)致隨機(jī)誤差和數(shù)據(jù)隨機(jī)擾動(dòng)的影響被忽視[10]。此外，DEA方法還增加了很強(qiáng)的凸性假設(shè)約束，這與實(shí)際生產(chǎn)過(guò)程不相符。

針對(duì)以上不足，本文在測(cè)算我國(guó)不同規(guī)模奶牛養(yǎng)殖環(huán)境效率時(shí)，擬做出以下改進(jìn)：第一，在測(cè)算出環(huán)境效率后，進(jìn)一步分析污染物的潛在削減量以及環(huán)境效率的提升路徑。第二，基于第一次全國(guó)污染源普查公報(bào)，將污染物的范圍拓展界定為化學(xué)需氧量、全氮、全磷、銅和鋅。第三，借鑒Reinhard等[4]的思路，利用SFA模型測(cè)算奶牛養(yǎng)殖的環(huán)境效率。

1 研究方法與數(shù)據(jù)說(shuō)明

1.1 研究方法 Reinhard等[4]提出環(huán)境效率的測(cè)算思路：假設(shè)合意產(chǎn)出和常規(guī)投入要素不變，通過(guò)計(jì)算污染物最小可能數(shù)量與當(dāng)前實(shí)際污染物數(shù)量的比值來(lái)衡量環(huán)境效率。如圖1所示，假設(shè)YR是等產(chǎn)量曲面，曲面上包含若干個(gè)等產(chǎn)量點(diǎn)，其中，A和B是2個(gè)等產(chǎn)量點(diǎn)。當(dāng)前，假設(shè)在A點(diǎn)進(jìn)行生產(chǎn)，此時(shí)需要XR個(gè)常規(guī)要素投入（如：勞動(dòng)、資本等），同時(shí)也會(huì)產(chǎn)生WR個(gè)污染物。根據(jù)Reinhard等[4]測(cè)算環(huán)境效率的思路，需要保持合意產(chǎn)出和常規(guī)投入要素不變，尋找出污染物的最小可能數(shù)量。因此，可以選擇在B點(diǎn)進(jìn)行生產(chǎn)，此時(shí)合意產(chǎn)出和常規(guī)投入要素不變，而污染物數(shù)量減少至WF。那么，環(huán)境效率可表示為：EE=OWF/OWR。

圖1 環(huán)境效率測(cè)算示意圖

在測(cè)算奶牛養(yǎng)殖環(huán)境效率時(shí)，本文將SFA模型設(shè)定為超越對(duì)數(shù)生產(chǎn)函數(shù)形式。一方面，超越對(duì)數(shù)生產(chǎn)函數(shù)形式相對(duì)靈活，能夠有效地降低模型設(shè)定誤差；另一方面，超越對(duì)數(shù)生產(chǎn)函數(shù)是變彈性生產(chǎn)函數(shù)，能夠較好地研究各投入要素的相互影響。模型構(gòu)造如下：

其中，Yit表示奶牛養(yǎng)殖的產(chǎn)出；Xit表示奶牛養(yǎng)殖的投入要素，諸如勞動(dòng)、資本等；Wit表示奶牛養(yǎng)殖排泄物中的污染物質(zhì)；T表示時(shí)間趨勢(shì)項(xiàng)，這里用來(lái)考察奶牛養(yǎng)殖的技術(shù)進(jìn)步；β為一系列待估計(jì)參數(shù)；vit為一般誤差項(xiàng)，表示隨機(jī)因素對(duì)產(chǎn)出的影響，服從獨(dú)立于uit的正態(tài)分布N(0，)；uit為技術(shù)非效率項(xiàng)，服從獨(dú)立于vit的半正態(tài)分布N+(0，)。

根據(jù)環(huán)境效率的測(cè)算思路，需要保持合意產(chǎn)出和常規(guī)投入要素不變，用污染物可能達(dá)到的最小數(shù)量WF來(lái)替代當(dāng)前污染物的實(shí)際數(shù)量W。此外，需要說(shuō)明的是，環(huán)境效率的測(cè)算是在不存在技術(shù)效率損失的假設(shè)下實(shí)現(xiàn)的，即uit=0。公式（1）可以改寫(xiě)為：

由上文知道，環(huán)境效率可以表示為EE=OWF/OWR，其對(duì)數(shù)形式為lnEE= lnWF-lnW。利用公式（2）減去公式（1），并整理成lnWF-lnW的形式，得：

由公式（3）可知，環(huán)境效率值實(shí)際上就是關(guān)于（lnWF-lnW）的一元二次方程的解，因此通過(guò)數(shù)學(xué)推導(dǎo)即可求出奶牛養(yǎng)殖的環(huán)境效率值。考慮到環(huán)境效率值通常在0～1，故需舍棄小于0或者大于1的解。奶牛養(yǎng)殖環(huán)境效率可表示為：

1.2 數(shù)據(jù)來(lái)源 本文數(shù)據(jù)來(lái)源于2004—2016年《全國(guó)農(nóng)產(chǎn)品成本收益匯編》（簡(jiǎn)稱(chēng)“農(nóng)本資料”）、《中國(guó)統(tǒng)計(jì)年鑒》以及2009年公布的《第一次全國(guó)污染源普查公報(bào)》（簡(jiǎn)稱(chēng)“公報(bào)”）。其中，奶牛養(yǎng)殖的投入、產(chǎn)出數(shù)據(jù)來(lái)自于《農(nóng)本資料》；相關(guān)價(jià)格指數(shù)數(shù)據(jù)來(lái)自《中國(guó)統(tǒng)計(jì)年鑒》；污染物的排放系數(shù)來(lái)自于《公報(bào)》。

需要說(shuō)明的是，《農(nóng)本資料》由國(guó)家發(fā)改委組織調(diào)查和匯編，是目前反映我國(guó)農(nóng)業(yè)生產(chǎn)成本收益的最權(quán)威資料。關(guān)于奶牛養(yǎng)殖的調(diào)查，《農(nóng)本資料》按散養(yǎng)、小規(guī)模、中規(guī)模和大規(guī)模進(jìn)行了統(tǒng)計(jì)，考慮到奶牛散養(yǎng)模式的樣本相對(duì)較少，因此本文針對(duì)規(guī)模養(yǎng)殖進(jìn)行研究，選取了小規(guī)模、中規(guī)模和大規(guī)模奶牛養(yǎng)殖作為研究對(duì)象。第二，《公報(bào)》由環(huán)境保護(hù)部、國(guó)家統(tǒng)計(jì)局以及農(nóng)業(yè)部聯(lián)合調(diào)查，是我國(guó)第一次全國(guó)污染源普查數(shù)據(jù)，具有一定的全面性和科學(xué)性。第三，根據(jù)《公報(bào)》數(shù)據(jù)，本文將不同規(guī)模奶牛養(yǎng)殖的樣本省份劃分為華北、東北、華東、中南、西南和西北六大地區(qū)進(jìn)行對(duì)比研究。

1.3 變量說(shuō)明 本文選取的變量包括合意產(chǎn)出、要素投入和污染物3類(lèi)，具體說(shuō)明如下：

1）合意產(chǎn)出（Y，kg）：一般而言，奶牛養(yǎng)殖的收益主要來(lái)自原料奶，因此本文用主產(chǎn)品產(chǎn)量來(lái)表示合意產(chǎn)出。

2）要素投入：①勞動(dòng)力投入（X1，工日）：本文采用用工數(shù)量表示勞動(dòng)力投入，包括家庭用工和農(nóng)業(yè)雇工兩部分。②飼料投入（X2，元）：飼料投入為精飼料費(fèi)用和青粗飼料費(fèi)之和?？紤]到價(jià)格波動(dòng)的影響，以2004年為基期用飼料價(jià)格指數(shù)對(duì)飼料投入進(jìn)行平減。③固定資產(chǎn)投入（X3，元）：通常來(lái)說(shuō)，畜禽的規(guī)模養(yǎng)殖需要進(jìn)行大量投資，因此本文將固定資產(chǎn)作為一個(gè)投入指標(biāo)。固定資產(chǎn)投入為固定資產(chǎn)折舊、修理維護(hù)費(fèi)、工具材料費(fèi)和飼料加工費(fèi)之和。同樣地，采用固定資產(chǎn)投資價(jià)格指數(shù)對(duì)固定資產(chǎn)投入進(jìn)行平減。④其他投入費(fèi)用（X4，元）：其他投入費(fèi)用主要包括水費(fèi)、燃料動(dòng)力費(fèi)和醫(yī)療防疫費(fèi)等，用畜禽生產(chǎn)價(jià)格指數(shù)進(jìn)行平減。

3）污染物（X5，kg）：污染物主要來(lái)自奶牛排泄物，本文采用排污系數(shù)，而非產(chǎn)污系數(shù)。排污系數(shù)即在典型的正常生產(chǎn)和管理?xiàng)l件下，單個(gè)畜禽產(chǎn)生的原始污染物經(jīng)處理設(shè)施處理或利用后，或未經(jīng)處理利用而直接排放到環(huán)境中的污染物量。產(chǎn)污系數(shù)即在典型的正常生產(chǎn)和管理?xiàng)l件下，一定時(shí)間內(nèi)（以天為單位）單個(gè)畜禽所產(chǎn)生的原始污染物量?？紤]到雖然有些地區(qū)奶牛養(yǎng)殖產(chǎn)生的污染物多，但經(jīng)過(guò)科學(xué)處理，污染物可以降到最低，本文依然認(rèn)為這是生態(tài)環(huán)保型養(yǎng)殖，因此采用排污系數(shù)更合理。需要說(shuō)明的是，《農(nóng)本資料》將奶牛規(guī)模養(yǎng)殖劃分為了小規(guī)模、中規(guī)模和大規(guī)模3類(lèi)，而《公報(bào)》中將奶牛規(guī)模養(yǎng)殖劃分為了養(yǎng)殖專(zhuān)業(yè)戶、規(guī)?；B(yǎng)殖場(chǎng)以及養(yǎng)殖小區(qū)3類(lèi)。借鑒朱寧等[8]、李翠霞等[9]做法，將小規(guī)模與養(yǎng)殖專(zhuān)業(yè)戶、中規(guī)模與規(guī)?；B(yǎng)殖場(chǎng)、大規(guī)模與養(yǎng)殖小區(qū)分別對(duì)應(yīng)起來(lái)進(jìn)行研究。此外，為了方便實(shí)證分析，本文將污染物加總折算為一個(gè)變量。最后，本文污染物的計(jì)算公式如下：

其中，i為污染物種類(lèi)，本文主要指化學(xué)需氧量、全氮、全磷、銅和鋅。此外，這5種污染物的排污系數(shù)i均來(lái)自《公報(bào)》；平均飼養(yǎng)天數(shù)來(lái)自《農(nóng)本資料》。

2 實(shí)證研究

2.1 模型檢驗(yàn)和估計(jì) 利用單邊似然比檢驗(yàn)（LR）分別對(duì)變差率γ以及函數(shù)設(shè)定形式進(jìn)行檢驗(yàn)，檢驗(yàn)結(jié)果見(jiàn)表1。首先，變差率γ檢驗(yàn)考察的是uit是否為0。如果uit=0，說(shuō)明方程不存在非效率項(xiàng)，即沒(méi)有必要采用SFA模型進(jìn)行分析。結(jié)果顯示，模型在1%的水平上通過(guò)檢驗(yàn)，說(shuō)明無(wú)論小規(guī)模、中規(guī)模還是大規(guī)模奶牛養(yǎng)殖，均存在養(yǎng)殖無(wú)效率情況，因此采用SFA模型進(jìn)行分析是有必要的。其次，函數(shù)設(shè)定形式檢驗(yàn)考察的是函數(shù)形式設(shè)定為超越對(duì)數(shù)形式還是設(shè)定為Cobb-Douglas形式。結(jié)果顯示，模型均在1%的水平上通過(guò)檢驗(yàn)，說(shuō)明無(wú)論小規(guī)模、中規(guī)模還是大規(guī)模奶牛養(yǎng)殖，將生產(chǎn)函數(shù)形式設(shè)定為超越對(duì)數(shù)形式相對(duì)更合理。利用Stata12軟件對(duì)小規(guī)模、中規(guī)模以及大規(guī)模奶牛養(yǎng)殖的隨機(jī)前沿生產(chǎn)函數(shù)進(jìn)行了估計(jì)，相關(guān)估計(jì)系數(shù)見(jiàn)表2。

2.2 不同規(guī)模奶牛養(yǎng)殖的環(huán)境效率 利用表2的參數(shù)估計(jì)值，根據(jù)公式（4）計(jì)算不同規(guī)模奶牛養(yǎng)殖的環(huán)境效率，結(jié)果見(jiàn)表3。首先，小規(guī)模、中規(guī)模、大規(guī)模奶牛養(yǎng)殖的環(huán)境效率值依次為0.467、0.503和0.603，呈現(xiàn)出小規(guī)模、中規(guī)模、大規(guī)模依次遞增的趨勢(shì)。一種可能的解釋是隨著奶牛養(yǎng)殖規(guī)模增大，糞污排放會(huì)越來(lái)越多，如果不科學(xué)地處理則會(huì)危及奶牛養(yǎng)殖的可持續(xù)發(fā)展[11]，因此大規(guī)模養(yǎng)殖場(chǎng)可能更傾向于購(gòu)買(mǎi)先進(jìn)的設(shè)備、采用先進(jìn)的技術(shù)處理糞污，從而更能實(shí)現(xiàn)清潔生產(chǎn)，所以大規(guī)模奶牛養(yǎng)殖場(chǎng)的環(huán)境效率會(huì)高于中、小規(guī)模奶牛養(yǎng)殖的環(huán)境效率。其次，從環(huán)境效率的變化趨勢(shì)來(lái)看，無(wú)論是小規(guī)模、中規(guī)模還是大規(guī)模，隨著時(shí)間推移，奶牛養(yǎng)殖的環(huán)境效率均呈現(xiàn)出不同程度的上升趨勢(shì)，且這種趨勢(shì)在中、小規(guī)模養(yǎng)殖場(chǎng)中體現(xiàn)得更明顯，說(shuō)明我國(guó)奶牛養(yǎng)殖的環(huán)境效率在逐步改善。這可能與我國(guó)政府近年來(lái)高度重視生態(tài)環(huán)境問(wèn)題并加強(qiáng)對(duì)環(huán)境污染監(jiān)管有關(guān)，早在2001年我國(guó)政府就對(duì)畜禽養(yǎng)殖污染物的排放標(biāo)準(zhǔn)以及排放技術(shù)做了相關(guān)規(guī)定。2013年中央“一號(hào)文件”又明確提出要開(kāi)展畜禽養(yǎng)殖污染防治工作。2016年中央政府成立了“環(huán)保督查組”，加強(qiáng)了對(duì)各地區(qū)環(huán)境污染的監(jiān)督和管理。

表1 模型設(shè)定形式檢驗(yàn)

表2 不同規(guī)模奶牛養(yǎng)殖隨機(jī)前沿生產(chǎn)函數(shù)的估計(jì)

表3 不同規(guī)模奶牛養(yǎng)殖的環(huán)境效率值

有意思的是，雖然不同規(guī)模奶牛養(yǎng)殖環(huán)境效率均出現(xiàn)上升趨勢(shì)，但相對(duì)于大規(guī)模奶牛養(yǎng)殖，中、小規(guī)模奶牛養(yǎng)殖環(huán)境效率提高的更加明顯?？赡茉蚴桥c大規(guī)模奶牛養(yǎng)殖場(chǎng)相比，中、小規(guī)模養(yǎng)殖場(chǎng)的糞污處理設(shè)備可能會(huì)相對(duì)落后，因此環(huán)境效率會(huì)處在相對(duì)較低水平，但近年來(lái)，政府加強(qiáng)了對(duì)畜禽養(yǎng)殖環(huán)境的治理，中、小規(guī)模養(yǎng)殖場(chǎng)可能會(huì)更新?lián)Q代糞污處理設(shè)備，因此環(huán)境效率提升相對(duì)明顯。相反，為了維持運(yùn)行的穩(wěn)定性以及養(yǎng)殖場(chǎng)的可持續(xù)發(fā)展，大規(guī)模養(yǎng)殖場(chǎng)可能一開(kāi)始就采用先進(jìn)的糞污處理設(shè)備，環(huán)境效率一直處在較高水平，而大型糞污處理設(shè)備的研發(fā)和創(chuàng)新需要一定過(guò)程，因此大規(guī)模養(yǎng)殖場(chǎng)的環(huán)境效率提升相對(duì)較慢。

此外，自2012年以來(lái)中規(guī)模養(yǎng)殖的環(huán)境效率出現(xiàn)了略微下滑，崔姹等[1]的研究中也出現(xiàn)了類(lèi)似現(xiàn)象。可能原因是2012年以來(lái)我國(guó)奶牛養(yǎng)殖成本增加、養(yǎng)殖效益下降，加上政府提高了畜禽養(yǎng)殖污染的治理標(biāo)準(zhǔn)，這變相淘汰了部分環(huán)境效率低下的小規(guī)模養(yǎng)殖場(chǎng)，使得小規(guī)模奶牛養(yǎng)殖效率出現(xiàn)了提升。散戶和小規(guī)模養(yǎng)殖場(chǎng)相繼退出，為了保證市場(chǎng)供給，中規(guī)模和大規(guī)模養(yǎng)殖場(chǎng)可能會(huì)提高牛奶產(chǎn)量，而產(chǎn)出的增加會(huì)帶來(lái)畜禽糞污的增加，但大規(guī)模養(yǎng)殖場(chǎng)具有較強(qiáng)的糞污處理能力，因此不會(huì)帶來(lái)環(huán)境效率下降，而中規(guī)模養(yǎng)殖場(chǎng)糞污處理能力相對(duì)有限，因此糞污增加在一定程度上會(huì)拉低奶牛養(yǎng)殖環(huán)境效率，所以環(huán)境效率出現(xiàn)了略微下降。

2.3 不同規(guī)模奶牛養(yǎng)殖污染物的潛在削減量 根據(jù)環(huán)境效率的定義，維持合意產(chǎn)出和常規(guī)投入要素不變，污染物能夠達(dá)到的最小可能數(shù)量與當(dāng)前實(shí)際污染物數(shù)量的比值即為環(huán)境效率，這意味著可以通過(guò)環(huán)境效率值測(cè)算出不同規(guī)模奶牛養(yǎng)殖污染物排放數(shù)量的潛在削減量，這為奶牛養(yǎng)殖污染物的處理提供了一定的理論參考依據(jù)，對(duì)實(shí)現(xiàn)環(huán)境友好型奶牛養(yǎng)殖具有重要意義。利用公式“污染物潛在削減量=污染物當(dāng)前實(shí)際排放量×（1-當(dāng)前環(huán)境效率值）”計(jì)算不同規(guī)模奶牛養(yǎng)殖污染物的潛在削減量，結(jié)果見(jiàn)表4。需要說(shuō)明的是，公式中“污染物當(dāng)前實(shí)際排放量”表示的是不同養(yǎng)殖規(guī)模內(nèi)部各省份2004—2016年奶牛養(yǎng)殖排放污染物的均值；而“當(dāng)前環(huán)境效率值”表示的是不同養(yǎng)殖規(guī)模內(nèi)部各省份2004—2016年奶牛養(yǎng)殖環(huán)境效率的均值，因此本文計(jì)算的污染物潛在削減量是2004—2016年的平均水平?？傮w來(lái)看，無(wú)論小規(guī)模、中規(guī)模還是大規(guī)模奶牛養(yǎng)殖，一頭牛每年產(chǎn)生的污染物還可以進(jìn)一步削減。平均而言，小規(guī)模、中規(guī)模和大規(guī)模奶牛養(yǎng)殖場(chǎng)養(yǎng)殖的一頭奶牛一年產(chǎn)生的污染物還可以依次減少819.70、773.39、538.92 kg，這說(shuō)明通過(guò)相關(guān)技術(shù)手段，我國(guó)奶牛養(yǎng)殖產(chǎn)生的污染物還有進(jìn)一步削減的可能。

分規(guī)模來(lái)看。首先，就環(huán)境效率值而言，無(wú)論小規(guī)模、中規(guī)模還是大規(guī)模奶牛養(yǎng)殖，其內(nèi)部省份奶牛養(yǎng)殖環(huán)境效率均有嚴(yán)重的兩級(jí)分化趨勢(shì)。以小規(guī)模為例，寧夏、新疆等西北地區(qū)奶牛養(yǎng)殖的環(huán)境效率相對(duì)較高，效率值在0.80左右，而河南、廣西、湖南等中南地區(qū)奶牛養(yǎng)殖的環(huán)境效率相對(duì)較低，僅在0.35左右，呈現(xiàn)出明顯的地域差異。其次，就污染物潛在削減量而言，一個(gè)有趣的現(xiàn)象是環(huán)境效率較低的地區(qū)往往具有較高的潛在削減量。如：小規(guī)模養(yǎng)殖地區(qū)中的安徽省，其環(huán)境效率僅有0.35，而其奶牛養(yǎng)殖污染物排放達(dá)到了2 050.00 kg，潛在削減量高達(dá)1 324.30 kg。同樣地，中規(guī)模奶牛養(yǎng)殖地區(qū)中的廣西省以及大規(guī)模奶牛養(yǎng)殖地區(qū)中的廣東省均具有類(lèi)似現(xiàn)象。值得注意的是，在不同養(yǎng)殖規(guī)模內(nèi)部，奶牛養(yǎng)殖污染物的潛在削減量存在明顯的區(qū)域差異，這可能與各地區(qū)糞污處理方式有關(guān)，如中、西部地區(qū)普遍采用干清糞工藝，糞污大多轉(zhuǎn)化為肥料還田，而東部地區(qū)普遍采用水沖清糞工藝，糞污大多用來(lái)發(fā)酵成生物能源，糞污處理方式的差異會(huì)影響到奶牛養(yǎng)殖環(huán)境效率，進(jìn)而影響到奶牛養(yǎng)殖污染物的潛在削減量。

表4 不同規(guī)模奶牛養(yǎng)殖污染物的潛在削減量

3 不同規(guī)模奶牛養(yǎng)殖環(huán)境效率的提升路徑

一般而言，影響環(huán)境效率的因素有很多，諸如養(yǎng)殖方式、技術(shù)投入、環(huán)境治理等。因此，本文擬從以下幾個(gè)方面考察影響因素：①養(yǎng)殖方式，其最直觀的表現(xiàn)是生產(chǎn)要素投入狀況。本文用勞動(dòng)力投入、飼料投入、固定資產(chǎn)投入等來(lái)表示。②技術(shù)投入，本文用奶牛養(yǎng)殖技術(shù)服務(wù)費(fèi)來(lái)表示。③雇工比重：由于雇工長(zhǎng)期從事奶牛養(yǎng)殖工作，通常具有一定的養(yǎng)殖經(jīng)驗(yàn)，雇工的使用在一定程度上會(huì)提升環(huán)境效率，本文用雇工數(shù)量占總投工數(shù)量的比重來(lái)表示。④環(huán)境治理：一直以來(lái)，我國(guó)采用“誰(shuí)污染、誰(shuí)治理、誰(shuí)破壞、誰(shuí)恢復(fù)”的原則治理環(huán)境污染，因此環(huán)境治理費(fèi)用通常由地方政府支出。環(huán)境治理費(fèi)用越多，擁有高環(huán)境效率的可能性越大。本文采用地方財(cái)政環(huán)境保護(hù)支出來(lái)表征環(huán)境治理。⑤地區(qū)虛擬變量：考慮到不同地區(qū)的養(yǎng)殖方式差異，本文在方程中也加入了地區(qū)虛擬變量。

基于以上分析，本文擬構(gòu)造如下的奶牛養(yǎng)殖環(huán)境效率的影響因素方程：

其中，EEit表示奶牛養(yǎng)殖環(huán)境效率；Xit表示影響奶牛養(yǎng)殖環(huán)境效率的諸多因素，如：養(yǎng)殖方式、技術(shù)投入、環(huán)境治理等；D表示地區(qū)虛擬變量；εit表示隨機(jī)誤差項(xiàng)；α0、α1、α2表示一系列待估計(jì)參數(shù)?？紤]到奶牛養(yǎng)殖環(huán)境效率值在0～1，如果直接用線性方程回歸，則會(huì)使得估計(jì)結(jié)果有偏，相反采用Tobit模型估計(jì)不僅無(wú)偏而且有效[12]，因此本文擬采用Tobit模型估計(jì)方程式（6）。

不同規(guī)模奶牛養(yǎng)殖環(huán)境效率影響因素的估計(jì)結(jié)果見(jiàn)表5。首先，就養(yǎng)殖方式而言，無(wú)論小規(guī)模、中規(guī)模還是大規(guī)模奶牛養(yǎng)殖，本文選取的勞動(dòng)力投入、飼料投入以及固定資產(chǎn)投入均對(duì)環(huán)境效率具有顯著影響。其中，總投工及固定資產(chǎn)投入的估計(jì)系數(shù)均在1%水平上顯著為正，說(shuō)明投工量增加及固定資產(chǎn)投資增加均會(huì)顯著提高奶牛養(yǎng)殖環(huán)境效率。此外，飼料費(fèi)用的估計(jì)系數(shù)在小規(guī)模及中規(guī)模奶牛養(yǎng)殖中顯著為正，而在大規(guī)模養(yǎng)殖中顯著為負(fù)，可能原因在于大規(guī)模奶牛養(yǎng)殖場(chǎng)具有產(chǎn)量?jī)?yōu)勢(shì)，如果單純比產(chǎn)量，中、小規(guī)模奶牛養(yǎng)殖場(chǎng)難以在競(jìng)爭(zhēng)中取勝，因此中、小規(guī)模養(yǎng)殖場(chǎng)可能會(huì)轉(zhuǎn)向追求產(chǎn)品品質(zhì)，合理使用飼料，而大規(guī)模養(yǎng)殖往往會(huì)追求產(chǎn)品產(chǎn)量，可能會(huì)存在過(guò)量使用飼料的現(xiàn)象，從而增加了奶牛排泄，降低了環(huán)境效率。其次，技術(shù)服務(wù)費(fèi)的估計(jì)系數(shù)僅在大規(guī)模養(yǎng)殖中顯著為正，而在中、小規(guī)模養(yǎng)殖中不顯著，這可能是由于中、小規(guī)模奶牛養(yǎng)殖的規(guī)模有限，難以很好體現(xiàn)技術(shù)服務(wù)的規(guī)模效應(yīng)。然后，雇工比重的估計(jì)系數(shù)在中規(guī)模和大規(guī)模奶牛養(yǎng)殖中均顯著為正，說(shuō)明雇工的使用會(huì)提升奶牛養(yǎng)殖環(huán)境效率；但雇工比重的估計(jì)系數(shù)在小規(guī)模奶牛養(yǎng)殖中不顯著，說(shuō)明雇工比重對(duì)小規(guī)模奶牛養(yǎng)殖的環(huán)境效率沒(méi)有影響，這可能是由于小規(guī)模養(yǎng)殖大多采用家庭自用工，而較少使用雇工。最后，地方財(cái)政環(huán)境保護(hù)支出的估計(jì)系數(shù)在小規(guī)模和大規(guī)模奶牛養(yǎng)殖中顯著為正，說(shuō)明地方政府環(huán)境保護(hù)支出的增加會(huì)提升奶牛養(yǎng)殖環(huán)境效率。

4 結(jié)果與討論

本研究利用2004—2016年奶牛養(yǎng)殖的省級(jí)宏觀數(shù)據(jù)，基于SFA模型測(cè)算了我國(guó)不同規(guī)模奶牛養(yǎng)殖的環(huán)境效率，并分析了奶牛養(yǎng)殖產(chǎn)生污染物的潛在削減量，在此基礎(chǔ)上進(jìn)一步考察了奶牛養(yǎng)殖環(huán)境效率的提升路徑。主要結(jié)論：①我國(guó)小規(guī)模、中規(guī)模、大規(guī)模奶牛養(yǎng)殖的環(huán)境效率值依次為0.467、0.503和0.603，呈現(xiàn)出小規(guī)模、中規(guī)模、大規(guī)模依次遞增的趨勢(shì)，且隨著時(shí)間推移，奶牛養(yǎng)殖的環(huán)境效率均呈現(xiàn)出不同程度的上升趨勢(shì)，說(shuō)明我國(guó)奶牛養(yǎng)殖的環(huán)境效率在逐步改善。②我國(guó)奶牛養(yǎng)殖產(chǎn)生的污染物還有進(jìn)一步削減的可能。平均而言，小規(guī)模、中規(guī)模和大規(guī)模奶牛養(yǎng)殖場(chǎng)養(yǎng)殖的一頭奶牛一年產(chǎn)生的污染物還可以依次減少819.70、773.39和538.92 kg。③養(yǎng)殖方式、技術(shù)投入、雇工比重以及環(huán)境治理等因素均會(huì)顯著影響奶牛養(yǎng)殖環(huán)境效率。

對(duì)此，本文提出以下建議：第一，建議采用科學(xué)的養(yǎng)殖方式飼養(yǎng)奶牛。合理配置勞動(dòng)力及飼料等投入要素，避免出現(xiàn)過(guò)量使用現(xiàn)象。第二，考慮到固定資產(chǎn)投入對(duì)環(huán)境效率具有正向作用，建議進(jìn)一步加強(qiáng)固定資產(chǎn)投資，引進(jìn)先進(jìn)的環(huán)保設(shè)施，提升污染物的處理能力。第三，由于雇工對(duì)環(huán)境效率具有正向作用，建議施行專(zhuān)業(yè)化養(yǎng)殖，雇傭高素質(zhì)技術(shù)人員，提升管理運(yùn)營(yíng)能力。第四，建議地方政府加強(qiáng)對(duì)環(huán)境污染的監(jiān)管，有效落實(shí)環(huán)境治理資金的發(fā)放與使用，保證環(huán)境治理資金落到實(shí)處，促進(jìn)我國(guó)奶牛養(yǎng)殖向環(huán)保型養(yǎng)殖方式轉(zhuǎn)變。

表5 不同規(guī)模奶牛養(yǎng)殖環(huán)境效率的影響因素

需要說(shuō)明的是，現(xiàn)實(shí)中影響奶牛養(yǎng)殖環(huán)境效率的因素太多，諸如飼料種類(lèi)、養(yǎng)殖方式、養(yǎng)殖技術(shù)、糞污處理設(shè)施、政府環(huán)境監(jiān)管力度等等，而且其中有些因素并不能通過(guò)技術(shù)外溢、人力資本投資等渠道消除，所以有些地區(qū)奶牛養(yǎng)殖的環(huán)境效率值并不能達(dá)到1，因此本文計(jì)算出的奶牛養(yǎng)殖污染物的潛在削減量是各地區(qū)實(shí)際奶牛養(yǎng)殖中可以達(dá)到的上限。實(shí)際中會(huì)存在有些地區(qū)污染物削減的多、有些地區(qū)污染物削減的少的情況。即便如此，本研究依然能夠?yàn)槲覈?guó)不同規(guī)模奶牛養(yǎng)殖污染物的削減提供理論參考。