王善高，許 昭

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

2 實證研究

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

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

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

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

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

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

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

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

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

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

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

一般而言，影響環(huán)境效率的因素有很多，諸如養(yǎng)殖方式、技術(shù)投入、環(huán)境治理等。因此，本文擬從以下幾個方面考察影響因素：①養(yǎng)殖方式，其最直觀的表現(xiàn)是生產(chǎn)要素投入狀況。本文用勞動力投入、飼料投入、固定資產(chǎn)投入等來表示。②技術(shù)投入，本文用奶牛養(yǎng)殖技術(shù)服務(wù)費來表示。③雇工比重：由于雇工長期從事奶牛養(yǎng)殖工作，通常具有一定的養(yǎng)殖經(jīng)驗，雇工的使用在一定程度上會提升環(huán)境效率，本文用雇工數(shù)量占總投工數(shù)量的比重來表示。④環(huán)境治理：一直以來，我國采用“誰污染、誰治理、誰破壞、誰恢復”的原則治理環(huán)境污染，因此環(huán)境治理費用通常由地方政府支出。環(huán)境治理費用越多，擁有高環(huán)境效率的可能性越大。本文采用地方財政環(huán)境保護支出來表征環(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表示隨機誤差項；α0、α1、α2表示一系列待估計參數(shù)?？紤]到奶牛養(yǎng)殖環(huán)境效率值在0～1，如果直接用線性方程回歸，則會使得估計結(jié)果有偏，相反采用Tobit模型估計不僅無偏而且有效[12]，因此本文擬采用Tobit模型估計方程式（6）。

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

4 結(jié)果與討論

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

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

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

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