宋鵬志 ，李翠霞 ，2

（1.東北農(nóng)業(yè)大學經(jīng)濟管理學院，哈爾濱150030；2.黑龍江省綠色食品科學研究院，哈爾濱，150028）

0 引 言

2001 年12 月11 日中國加入世界貿(mào)易組織（World Trade Organization，WTO）后，乳制品進口關稅大幅降低，乳制品進口量快速增長。2008 年-2017年，中國乳制品進口量從 38.7 萬 t 增加至 255.7 萬 t，年均復合增長率為21.5%。乳制品進口增加通過競爭效應對不同養(yǎng)殖主體產(chǎn)生差異化影響，對規(guī)模養(yǎng)殖主體有正向作用，對散戶有負向影響，帶動養(yǎng)殖業(yè)在更高水平上良性發(fā)展[1]，對國內(nèi)產(chǎn)業(yè)生產(chǎn)率增長有積極的正向影響。從統(tǒng)計數(shù)據(jù)分析，近年來，奶牛養(yǎng)殖散戶退出養(yǎng)殖市場的速度加快，規(guī)?；B(yǎng)殖顯著增加，規(guī)?；潭扔?007 年的46.64%提高到2018 年的61.4%，奶牛養(yǎng)殖越來越科學、精細化，規(guī)模化水平越來越高，質(zhì)量也得到保證，但是奶牛規(guī)?；B(yǎng)殖推進過快，規(guī)模養(yǎng)殖場的管理設備、防控設施未能同步改善，奶牛養(yǎng)殖主體面臨著糞污處理壓力，導致養(yǎng)殖環(huán)境效率降低[2]。因此，如何在提高產(chǎn)業(yè)率的同時，提升養(yǎng)殖環(huán)境效率、保證良好的生態(tài)環(huán)境是促進奶牛養(yǎng)殖業(yè)發(fā)展的關鍵[3-4]。本文在運用SBM-Undesirable 模型測定奶牛養(yǎng)殖環(huán)境效率的基礎上，Tobit 模型實證分析了乳制品進口對不同規(guī)模奶牛養(yǎng)殖環(huán)境效率的的異質(zhì)性，有助于更好的協(xié)調(diào)乳制品進口貿(mào)易與中國奶牛養(yǎng)殖環(huán)境之間健康快速的發(fā)展。

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

1.1 研究方法——SBM-Undesirable模型

Tone[5]在 2001 年提出 SBM 模型，SBM 模型作為DEA 模型中的非徑向、非角度模型，可以解決輸入松弛問題，在出現(xiàn)意外輸出的情況下解決效率評估問題，有效地避免徑向和角度選擇差異帶來的偏差和影響，更好地反映效率評估的實質(zhì)。因此，本文使用基于SBM 模型的非產(chǎn)出模型來衡量奶牛養(yǎng)殖環(huán)境效率。

本文以不同省份、不同規(guī)模的奶牛養(yǎng)殖場為決策單元，構建了奶牛養(yǎng)殖環(huán)境效益的可能邊界。假定奶牛養(yǎng)殖系統(tǒng)有n 個決策單元，且均有3 個投入產(chǎn)出向量：投入、期望產(chǎn)出、非期望產(chǎn)出，表示為x∈Rm、y

g∈Rs1及yb∈Rs2，定義矩陣X、Yg、Yb如下：

可變規(guī)模報酬下的生產(chǎn)可能性集P 表示為：

其中，λ∈Rn為權重向量。

根據(jù)Tone 等[6]定義，非期望產(chǎn)出的SBM-Undesirable模型為：

1.2 數(shù)據(jù)來源

本研究數(shù)據(jù)選取2008 年-2017 年《全國農(nóng)產(chǎn)品成本收益資料匯編》、《中國統(tǒng)計年鑒》、《第一次全國污染源普查——畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù)手冊》（2009年）[7]。其中奶牛養(yǎng)殖的投入和期望產(chǎn)出數(shù)據(jù)來自《全國農(nóng)產(chǎn)品成本收益資料匯編》（下文簡稱《匯編》）；奶牛養(yǎng)殖污染物排污系數(shù)來自《第一次全國污染源普查——畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù)手冊》（下文簡稱《手冊》）。

補充說明：第一，不同規(guī)模的奶牛養(yǎng)殖存在飼養(yǎng)方式、飼養(yǎng)周期等方面的差異，該差異可能會影響奶牛養(yǎng)殖的環(huán)境效率?！秴R編》中是基于散養(yǎng)、小規(guī)模、中規(guī)模和大規(guī)模奶牛養(yǎng)殖的統(tǒng)計數(shù)據(jù)，學者們通常對畜禽養(yǎng)殖污染按小規(guī)模、中規(guī)模和大規(guī)模進行分析[8-9]。由于奶牛散養(yǎng)模式的樣本相對較少，本文最終選取了小、中、大規(guī)模奶牛養(yǎng)殖為研究對象。第二，《手冊》中將畜禽養(yǎng)殖污染物排污系數(shù)按華北、東北、華東、中南、西南和西北6 個地區(qū)進行統(tǒng)計。因此，本文將不同規(guī)模的奶牛養(yǎng)殖樣本省份劃分為6 個區(qū)域：華北、東北、華東、中南、西南和西北。通過梳理，本文獲取的不同規(guī)模奶牛養(yǎng)殖樣本數(shù)依次為小規(guī)模16 個、中規(guī)模24個、大規(guī)模21個。

1.3 指標選擇

本文選取的奶牛養(yǎng)殖環(huán)境效率指標包括：投入指標（生產(chǎn)用工數(shù)量、精飼料數(shù)量和其他費用）、期望產(chǎn)出（主產(chǎn)品產(chǎn)量）和非期望產(chǎn)出（污染物）。

1.3.1 投入指標

生產(chǎn)用工數(shù)量（d）是指生產(chǎn)過程中生產(chǎn)者（包括其家庭成員）和雇傭工人直接勞動的天數(shù)；精飼料數(shù)量（kg）是奶牛養(yǎng)殖總投入的主要部分，是指實際耗用的各種精飼料的實物數(shù)量；其他費用（元）包含固定資產(chǎn)折舊、工具材料費、修理維護費。（其他費用以2008年為基期，利用《中國統(tǒng)計年鑒》的固定資產(chǎn)投資價格指數(shù)進行平減處理[10]）。

1.3.2 期望產(chǎn)出指標

奶牛養(yǎng)殖的主產(chǎn)品是牛奶，并且奶牛養(yǎng)殖收入來源主要來自主產(chǎn)品的產(chǎn)量，可以較為真實地反應奶牛養(yǎng)殖的產(chǎn)出狀態(tài)。

1.3.3 非期望產(chǎn)出指標

《手冊》中給出的畜禽養(yǎng)殖排污系數(shù)包括化學需氧量（CODCr）、總氮（TN）、總磷（TP）、銅（Cu）和鋅（Zn）等，其中化學需氧量的排放最大，在污染物總排放量中占90%以上，因此為了測算方便以及盡量避免多個非期望產(chǎn)出所帶來的結果偏差，本文確定化學需氧量的排放量（g/(頭·天)）為非期望產(chǎn)出。

1.4 實證結果與分析

1.4.1 不同規(guī)模奶牛養(yǎng)殖環(huán)境效率的差異分析

2008 年-2017 年，中國不同規(guī)模奶牛養(yǎng)殖場總體環(huán)境效率均在0.6 以下，如圖1 所示。其中，2008年-2017年，小規(guī)模、中規(guī)模和大規(guī)模奶牛養(yǎng)殖環(huán)境效率均值依次為0.473、0.486 和0.592，呈現(xiàn)出小、中、大規(guī)模逐步增長趨勢。隨著奶牛養(yǎng)殖規(guī)模的擴大，糞便排放物也持續(xù)增加，對環(huán)境造成越來越大的危害。如果沒有科學有效的處理，將嚴重影響奶牛養(yǎng)殖業(yè)的可持續(xù)發(fā)展。大規(guī)模養(yǎng)殖場更傾向于用先進的技術和設備處理糞污，因此環(huán)境效率比中、小規(guī)模奶牛養(yǎng)殖場略高[11]。

從變化趨勢來看，2008 年-2010 年，不同規(guī)模奶牛養(yǎng)殖環(huán)境效率均處于下降趨勢。2010 年-2012 年規(guī)模化初期，奶牛養(yǎng)殖環(huán)境效率處于較高水平，小規(guī)模養(yǎng)殖場投入和產(chǎn)出達到優(yōu)化所需時間相對較短，產(chǎn)奶量少、糞污處理等優(yōu)勢，導致其環(huán)境效率一直略高于中規(guī)模養(yǎng)殖場，中、大規(guī)模奶牛養(yǎng)殖場建設周期較長，效率達到最優(yōu)需要一定的時間。2011 年-2014年，不同規(guī)模奶牛養(yǎng)殖環(huán)境效率均值一直處于下降趨勢（環(huán)境效率值由0.555 下降至0.476），主要原因是中國奶牛養(yǎng)殖成本增加，政府發(fā)布了畜禽養(yǎng)殖污染處理標準，導致不合格的小、中規(guī)模養(yǎng)殖場逐漸退出養(yǎng)殖市場。為了確保市場供應，大規(guī)模養(yǎng)殖場增加牛奶產(chǎn)量，由于當時糞污治理能力有限，畜禽糞污排放量的增加在一定程度上降低了奶牛養(yǎng)殖的環(huán)境效率。2013 年-2014 年，生鮮乳市場價格下降，養(yǎng)殖場降低了養(yǎng)殖成本和糞污處理設施的積極性，小規(guī)模養(yǎng)殖場環(huán)境效率大幅下降[12]。自2015 年以來，由于中國政府近年來對生態(tài)和環(huán)境問題的重視以及加強環(huán)境污染治理，小規(guī)模和大規(guī)模奶牛養(yǎng)殖場的環(huán)境效率總體上一直波動上升提高，例如：2013 年，中央“一號文件”明確提出，要積極開展畜禽養(yǎng)殖污染防治工作；2014 年，《畜禽規(guī)模養(yǎng)殖污染防治條例》等；2016 年，中央政府成立“環(huán)境保護督察組”，以加強對各地區(qū)環(huán)境污染的監(jiān)督管理等。

圖1 2008年-2017年中國不同規(guī)模奶牛養(yǎng)殖場平均環(huán)境效率變化

1.4.2 省際差異分析

本文測算了2008 年-2017 年中國各省不同規(guī)模奶牛養(yǎng)殖場的環(huán)境效率，并求得不同規(guī)模奶牛養(yǎng)殖場環(huán)境效率的平均值，得到結果如表1 所示。根據(jù)測算，2008 年-2017 年，各省份奶牛養(yǎng)殖場平均環(huán)境效率值為0.513。其中，北京、天津、上海、黑龍江、甘肅、新疆均處在全國平均水平之上，這些省份屬于奶牛養(yǎng)殖環(huán)境效率較高地區(qū)，他們對糞污排放和環(huán)境治理等方面做的較好，能將有限資源有效轉(zhuǎn)化成“期望產(chǎn)出”，“非期望產(chǎn)出”相對較少。均低于全國平均水平以下的省份有，內(nèi)蒙古、陜西、遼寧、河南、江蘇、安徽、浙江、重慶、廣西、河北、廣東、青海，奶牛養(yǎng)殖環(huán)境效率有待提高。各省份小規(guī)模奶牛養(yǎng)殖場均處在中效率值階段，只有少量省份的環(huán)境效率值高于全國平均效率值；中規(guī)模奶牛養(yǎng)殖場，大部分省份處于中效率階段，只有浙江和廣西兩省處于低效率階段；大規(guī)模奶牛養(yǎng)殖場中處于高效率的省份最多，且大部分省份效率值均處于全國平均水平之上。

2 環(huán)境效率影響因素分析

2.1 模型選擇

為了探究哪些因素影響了中國奶牛養(yǎng)殖的環(huán)境效率，本文采用面板Tobit 模型進行驗證。由于SBM-Undesirable 模型測算的 2008 年-2017 年中國 28個省份奶牛養(yǎng)殖環(huán)境效率值為0～1，具有被切割和截斷的特點，本文采用隨機效應Tobit 模型能夠保證運算結果的無偏性，并消除誤差項的序列相關性。建立面板Tobit模型為：

其中下標 i 表示原奶主產(chǎn)省，t 為年份；Yit為 i 省 t年的奶牛養(yǎng)殖環(huán)境效率；IM 表示乳制品進口與國內(nèi)乳制品產(chǎn)量之比，反映乳制品進口貿(mào)易的水平；MP 表示原奶出售價格；LW 表示養(yǎng)殖勞工薪酬；PE 表示雇工比重；NEA 表示非環(huán)保意識；βi表示待定系數(shù)；uit表示隨機擾動項。

2.2 影響因素選擇

綜合學者們對環(huán)境效率影響因素的研究，主要涉及經(jīng)濟因素、外貿(mào)因素、人口因素、技術因素、產(chǎn)業(yè)結構因素、科技因素等[13-14]。本文根據(jù)研究實際情況選取乳制品進口占比、原奶市場價格、勞工薪酬、雇工比重和非環(huán)保意識，這5 個可能影響中國奶牛養(yǎng)殖環(huán)境效率的因素（如表2所示），采用Tobit回歸模型進行分析。

其中（1）乳制品進口是檢驗的關鍵變量，即乳制品進口量與國內(nèi)產(chǎn)量之比。這一指標可以更好的反映乳制品進口貿(mào)易水平及其對環(huán)境效率的影響。（2）原奶市場價格是由原奶在市場上的出售價格表示。原奶價格的上漲可能會使奶牛養(yǎng)殖業(yè)的收益增多，從而加大對環(huán)境治理技術的投入并提高治理水平。因此，原奶的出售價格可能會影響奶牛養(yǎng)殖的環(huán)境效率。（3）勞工工資是由奶牛養(yǎng)殖場付給勞工的薪酬表示。通常狀況下，勞工綜合能力水平越高，所得薪酬越高，奶牛養(yǎng)殖場中勞工水平越高相應對環(huán)境治理越重視，有助于提高奶牛養(yǎng)殖環(huán)境效率。（4）雇工比重是由每日雇工費用與人工總成本比值表示。由于雇傭工人從事奶牛養(yǎng)殖已有較長時間，因此他們通常具有一定的養(yǎng)殖經(jīng)驗，在一定程度上提高了環(huán)境效率。（5）非環(huán)保意識是由地區(qū)單位有效灌溉面積化肥施用量來表示。指標一方面可以反映出奶農(nóng)的環(huán)保意識，另一方面，奶牛產(chǎn)出的糞便可以代替化學肥料，有效減少化肥施用量，從而減少奶牛養(yǎng)殖的環(huán)境污染，非環(huán)保意識越高，奶牛養(yǎng)殖環(huán)境的效率越低。

表1 2008年-2017年中國各省份不同規(guī)模奶牛養(yǎng)殖場平均環(huán)境效率分布情況

表2 關鍵變量及其他變量描述性統(tǒng)計特征

2.3 回歸結果與分析

本文借助Stata16.0 軟件，對非比例變量進行對數(shù)化處理后，運用Hausman 檢驗判定采用隨機效應模型，相關分析表明，各影響因素之間的相關性不大。利用 Tobit 模型對 2008 年-2017 年中國 28 個省份小、中、大規(guī)模奶牛養(yǎng)殖場的影響因素進行回歸分析[16-17]，如表3所示。

表3 Tobit模型回歸結果

具體分析如下：（1）乳制品進口對小、中、大規(guī)模奶牛養(yǎng)殖場環(huán)境效率均為正向影響，且小規(guī)模、中規(guī)模至少通過了5%和1%水平的顯著檢驗。說明乳制品進口增加對規(guī)模奶牛養(yǎng)殖環(huán)境效率具有顯著正向影響，乳制品進口競爭對其均產(chǎn)生正向促進作用。小、中規(guī)模養(yǎng)殖場存在較大的進步空間，可以繼續(xù)推進規(guī)?；岣吣翀龉芾砟芰?，進一步提升養(yǎng)殖環(huán)境效率。與小、中規(guī)模養(yǎng)殖場相比，乳制品進口對大規(guī)模養(yǎng)殖場環(huán)境效率的正向影響十分微弱，說明在競爭中持續(xù)經(jīng)營的養(yǎng)殖場，盡管可以通過持續(xù)擴大養(yǎng)殖規(guī)模提高環(huán)境效率，但也并不是可以無限擴大，在養(yǎng)殖場達到適度規(guī)模后，繼續(xù)擴大養(yǎng)殖規(guī)模對環(huán)境效率的提升有限，即并不是養(yǎng)殖場規(guī)模越大，環(huán)境效率越高。（2）原奶出售價格和勞工薪資對小、中、大規(guī)模養(yǎng)殖場均具有負向影響，這說明提高原奶價格和養(yǎng)殖人員的工資并不能有效的提高養(yǎng)殖環(huán)境效率。一方面，國內(nèi)原奶價格提高使國內(nèi)乳企使用國外奶粉意愿更強，利潤空間更大，從而導致原奶銷量降低，奶牛養(yǎng)殖廠收益減少并忽視對環(huán)境的治理工作，且原奶價格全部通過顯著性檢驗；另一方面，勞工工資不能直接反映出養(yǎng)殖人員的飼養(yǎng)水平，只有中規(guī)模養(yǎng)殖場通過1%顯著性檢驗。（3）雇工比重在中、大規(guī)模奶牛養(yǎng)殖在中均顯著為正，且分別通過5%和10%顯著性檢驗，表明雇傭工人將提高奶牛養(yǎng)殖的環(huán)境效率；但是，雇傭工人的占比對小規(guī)模奶牛養(yǎng)殖產(chǎn)生負向影響且不顯著，小規(guī)模奶牛養(yǎng)殖可能主要是家庭自用工，雇傭工人較少。（4）非環(huán)保意識對小、中、大規(guī)模養(yǎng)殖場均具有負向的影響，且對小、大規(guī)模養(yǎng)殖場通過了1%水平的顯著性檢驗。說明化肥等非環(huán)保型生產(chǎn)材料投入越多，證明其環(huán)保意識越差，奶牛養(yǎng)殖場對奶牛養(yǎng)殖污染問題重視度越低，并且過度依賴化肥會降低有機肥的使用，不利于奶牛糞污的資源化利用和奶牛養(yǎng)殖環(huán)境效率的提高。

3 結論與建議

本文運用SBM-Undesirable 模型對2008 年-2017年中國28 個省份小規(guī)模、中規(guī)模和大規(guī)模奶牛養(yǎng)殖場環(huán)境效率進行測算與分析，并利用Tobit 模型探究乳制品進口增加對小規(guī)模、中規(guī)模和大規(guī)模養(yǎng)殖場環(huán)境效率的影響。

研究結果表明：（1）2008 年-2017 年，中國不同規(guī)模奶牛養(yǎng)殖場總體環(huán)境效率較低，均在0.6 以下，呈現(xiàn)出小規(guī)模、中規(guī)模和大規(guī)模逐步增長趨勢?？梢娭袊膛ｐB(yǎng)殖環(huán)境效率處于較低水平，目前仍存在重視奶牛養(yǎng)殖的經(jīng)濟效益最大化而忽視奶牛場糞污等污染物排放對環(huán)境影響的問題。根據(jù)奶牛養(yǎng)殖環(huán)境效率測算，隨著時間推移呈現(xiàn)出不同程度的增長，說明近幾年政府對奶牛養(yǎng)殖環(huán)保方面的政策制定與具體實施對環(huán)境效率提升有很大幫助，中國奶牛養(yǎng)殖環(huán)境效率的提升仍有很大潛力。（2）2008 年-2017 年，中國 28個省份奶牛養(yǎng)殖場平均環(huán)境效率值為0.513，發(fā)現(xiàn)不同省份之間均存在明顯差異。其中，北京、天津、上海、黑龍江、甘肅、新疆均處在全國平均水平之上，內(nèi)蒙古、陜西、遼寧、河南、江蘇、安徽、浙江、重慶、廣西、河北、廣東、青海，這些省份奶牛養(yǎng)殖環(huán)境效率需要加大改進力度。廣西省和安徽省環(huán)境效率較低，說明經(jīng)濟比較落后的省份忽視奶牛養(yǎng)殖環(huán)境治理問題同樣比較嚴重，國家和政府應針對不同省份養(yǎng)殖情況制定相關政策進行環(huán)保支持。（3）通過對奶牛養(yǎng)殖環(huán)境效率影響因素分析，發(fā)現(xiàn)乳制品進口對小、中、大規(guī)模奶牛養(yǎng)殖主體的環(huán)境效率具有顯著正向影響。乳制品進口整體推動了奶牛養(yǎng)殖業(yè)環(huán)境效率的提升，促進養(yǎng)殖業(yè)向更高水平發(fā)展，帶動產(chǎn)業(yè)整合。在進口競爭激勵的市場環(huán)境下，小規(guī)模養(yǎng)殖主體可通過擴大養(yǎng)殖規(guī)模提高環(huán)境效率水平，中規(guī)模養(yǎng)殖主體也有很大提升空間，而大規(guī)模養(yǎng)殖主體提升空間較為微弱。

因此，提出以下建議：首先應該科學客觀地審視中國乳制品進口貿(mào)易的現(xiàn)狀，制定適合中國乳制品市場的貿(mào)易開放政策，積極引導奶牛養(yǎng)殖業(yè)參與國際競爭，推進奶牛養(yǎng)殖業(yè)產(chǎn)業(yè)整合，提升養(yǎng)殖業(yè)總體環(huán)境效率，為奶牛和人類創(chuàng)造舒適、安全、無污染的居住和生活環(huán)境以及為乳業(yè)發(fā)展提供穩(wěn)定、健康的奶源。其次繼續(xù)促進奶牛養(yǎng)殖場規(guī)?；l(fā)展，加快小規(guī)模養(yǎng)殖主體的擴張，建立適度規(guī)模的奶牛養(yǎng)殖場，最大限度發(fā)揮規(guī)模效益，切忌“超大規(guī)模”帶來的負效應，盡快制定出適合中國現(xiàn)階段奶牛養(yǎng)殖場的環(huán)境治理方案和配套糞污處理設施設備，以減少污染物排放并提高環(huán)境效率。最后仍要加大對奶牛養(yǎng)殖場糞污等污染物處理相應政策法規(guī)的出臺，并根據(jù)不同省份的實際情況建立健全奶牛養(yǎng)殖環(huán)境治理補貼制度，促進種養(yǎng)結合、糞污等污染物資源化利用，進一步提高養(yǎng)殖場的環(huán)境效率，增強中國奶牛養(yǎng)殖業(yè)在國際競爭中的優(yōu)勢。