趙麗平， 王雅鵬， 何　可

(華中農(nóng)業(yè)大學 a.經(jīng)濟管理學院； b.湖北農(nóng)村發(fā)展研究中心，湖北 武漢 430070)

?

我國糧食生產(chǎn)的環(huán)境技術效率測度

趙麗平a，b， 王雅鵬a，b， 何可a，b

(華中農(nóng)業(yè)大學 a.經(jīng)濟管理學院； b.湖北農(nóng)村發(fā)展研究中心，湖北 武漢 430070)

摘要：采用方向距離函數(shù)和Kernel密度估計方法，對1992—2012年中國29個省(市、區(qū))的糧食生產(chǎn)環(huán)境技術效率進行了測度和分析，以考察城鎮(zhèn)化和工業(yè)化加速時期中國各省份以及各糧食生產(chǎn)功能區(qū)糧食生產(chǎn)與資源環(huán)境的協(xié)調情況。研究結果表明：(1)全國及各糧食生產(chǎn)功能區(qū)糧食生產(chǎn)環(huán)境狀況總體上呈現(xiàn)出逐漸惡化現(xiàn)象，其中，糧食主產(chǎn)區(qū)優(yōu)于糧食主銷區(qū)，糧食主銷區(qū)優(yōu)于糧食平衡區(qū)；(2)中國各省份糧食生產(chǎn)環(huán)境技術效率存在明顯的非均衡性；(3)總體而言，我國糧食生產(chǎn)環(huán)境技術效率的地區(qū)差距在逐漸擴大，呈現(xiàn)出兩級分布, 但是在各個糧食生產(chǎn)功能區(qū)形勢有所不同。

關鍵詞：糧食生產(chǎn)； 環(huán)境技術效率； 環(huán)境污染；生產(chǎn)資料

一、引言

隨著工業(yè)化、城鎮(zhèn)化進程的加快，以土地、勞動力以及農(nóng)資為代表的農(nóng)業(yè)生產(chǎn)要素成本不斷上升，糧食種植的比較收益開始下降，大量青壯年農(nóng)村勞動力開始流入到效益高的產(chǎn)業(yè)部門，農(nóng)業(yè)部門的“人口紅利”逐漸消失[1]。另一方面，長期以來形成的農(nóng)業(yè)粗放式經(jīng)營導致了較為嚴重的農(nóng)業(yè)“立體污染”，與“兩型農(nóng)業(yè)”建設要求漸行漸遠[2]。由于農(nóng)藥、化肥的過量施用和工業(yè)生產(chǎn)導致的土壤污染，使得我國糧食每年減產(chǎn)約100億公斤[3]。此外，糧食作物秸稈焚燒現(xiàn)象屢禁不止，不僅造成了嚴重的環(huán)境污染，而且還導致土壤質量下降，進而引起了糧食減產(chǎn)，這嚴重違背了社會主義新農(nóng)村建設中的“發(fā)展生態(tài)農(nóng)業(yè)生產(chǎn)方式，改善農(nóng)村環(huán)境”的要求。為了科學評估新形勢下糧食生產(chǎn)投入要素的利用效果，本文將環(huán)境因素納入到糧食生產(chǎn)技術效率測算體系，以便正確認識全國及各個糧食生產(chǎn)功能區(qū)的資源利用情況及其差異，進一步對各個糧食生產(chǎn)功能區(qū)的糧食生產(chǎn)進行統(tǒng)籌協(xié)調。

環(huán)境技術效率是評價生產(chǎn)者(或決策單元)在給定投入和技術條件下，同時追求最大產(chǎn)出和最小污染的效率大小[4]。環(huán)境技術效率測算的重點在于如何將環(huán)境變量納入到技術效率測算模型中。關于環(huán)境變量的處理，目前主要有兩種方法：一種方法是直接把環(huán)境污染物作為投入要素來處理[5]，但在實際操作中很難分清楚哪些要素投入是用于生產(chǎn)期望產(chǎn)出，哪些要素投入是用于治理環(huán)境污染；另一種方法是把環(huán)境污染物作為非期望產(chǎn)出引進到DEA模型，很多學者對此進行了探索,有學者通過單調遞減函數(shù)對非期望產(chǎn)出進行轉換，將轉換后的數(shù)據(jù)在技術條件不變的情況下納入到產(chǎn)出函數(shù)中[6]，也有學者采取環(huán)境污染價值核算方法，對環(huán)境污染物的污染治理成本和環(huán)境退化成本進行核算，并納入到生產(chǎn)率測算范疇進行計算[7]，更多的學者則采用方向距離函數(shù)對環(huán)境技術效率進行測算[8-10]。其中，方向距離函數(shù)不僅可以考慮環(huán)境污染損失和松弛量對效率的影響，而且具有非角度性質，無需選擇測算角度，同時兼顧投入的減少和產(chǎn)出的增加，與利潤最大化假設一致，具有更強的適應性[11]。

在農(nóng)業(yè)生產(chǎn)領域，李谷成[11]等運用方向距離函數(shù)對農(nóng)業(yè)環(huán)境技術效率進行了測算，也有學者基于方向距離函數(shù)的Malmquist-Luenberger 生產(chǎn)率指數(shù)對我國13個油菜主產(chǎn)區(qū)的油菜環(huán)境全要素生產(chǎn)率進行了測算[12]，但是很少有學者對各糧食生產(chǎn)功能區(qū)(主產(chǎn)區(qū)、主銷區(qū)和平衡區(qū))的糧食生產(chǎn)環(huán)境技術效率進行分區(qū)域研究。

糧食安全在我國具有重要的經(jīng)濟地位和戰(zhàn)略地位，但隨著工業(yè)化和城鎮(zhèn)化的快速推進，粗放式生產(chǎn)經(jīng)營方式所產(chǎn)生的城市污染物和工業(yè)污染物給我國糧食生產(chǎn)所依賴的水土等自然環(huán)境造成了極大破化。根據(jù)環(huán)境保護部與國土資源部開展的首次全國土壤污染狀況調查顯示，我國耕地環(huán)境質量堪憂，近20%耕地面積污染超標，其中鎘、汞、砷等無機物質等超標點位數(shù)達82.8%，而且其中不少地區(qū)為糧食主產(chǎn)區(qū)，工礦業(yè)、不合理的農(nóng)業(yè)生產(chǎn)等人類活動是造成土壤污染或超標的主要原因。在糧食生態(tài)環(huán)境所面臨的嚴峻形勢下，各糧食生產(chǎn)功能區(qū)的糧食生產(chǎn)環(huán)境有何變化？各功能區(qū)的糧食生產(chǎn)和環(huán)境是否是協(xié)調發(fā)展的？各功能區(qū)的協(xié)調情況是否存在一定的差異？這種差異在不同的時間段是如何演變的？為此，本文運用方向距離函數(shù)對我國城鎮(zhèn)化和工業(yè)化加速時期的糧食生產(chǎn)環(huán)境技術效率進行測算，并對各糧食生產(chǎn)功能區(qū)的糧食生產(chǎn)環(huán)境技術效率進行時空差異分析、協(xié)調分析和動態(tài)演變分析，以了解我國各功能區(qū)的糧食生產(chǎn)環(huán)境變化情況，以及各功能區(qū)域之間和每個功能區(qū)內部各省份之間的糧食生產(chǎn)環(huán)境差異，從而對我國糧食生產(chǎn)和環(huán)境的協(xié)調情況進行全面了解，以促進我國糧食生產(chǎn)的可持續(xù)發(fā)展和綠色轉型。

二、 糧食生產(chǎn)環(huán)境技術效率的測度

(一)研究方法

糧食生產(chǎn)過程中，生產(chǎn)單位除了生產(chǎn)人們所希望獲得的糧食產(chǎn)量這一期望產(chǎn)出，還會生產(chǎn)出人們所不愿意獲得的氮、磷等非期望產(chǎn)出，傳統(tǒng)的DEA模型只考慮了期望產(chǎn)出，如果存在非期望產(chǎn)出，就無法準確估算出生產(chǎn)單位的效率。Chung 等[13]在測度瑞典紙漿廠生產(chǎn)率時構造了方向距離函數(shù)(DDF)，很好地解決了這個問題。設產(chǎn)出的方向向量為g=(gy,-gc)，那么基于產(chǎn)出的方向距離函數(shù)可以表達為：

(1)

式(1)可通過以下線性規(guī)劃方程來對每個決策單元的效率值進行求解。

(2)

式(2)中，k表示決策單元個數(shù)，t表示時期，當方向距離函數(shù)值為0時，技術效率值為1，表示該決策單元的效率值最高，位于前沿生產(chǎn)面上，否則就表示技術沒有效率。

(二) 數(shù)據(jù)來源及指標分析

1． 數(shù)據(jù)來源

本文的研究樣本為1992—2012年中國大陸29個省(市、區(qū))?？紤]到數(shù)據(jù)的可獲得性，以及方便處理和分析，樣本中未包括西藏，并將重慶市并入四川省。根據(jù)鐘甫寧等[14]和陸文聰?shù)萚15]的研究成果以及改革開放以來各省(市、區(qū))的糧食生產(chǎn)總量和消費量，將全國29個省(市、區(qū))劃分為糧食主產(chǎn)區(qū)、主銷區(qū)和平衡區(qū)3個糧食生產(chǎn)功能區(qū)，其中，糧食主產(chǎn)區(qū)為河北、內蒙古、遼寧、吉林、黑龍江、江蘇、安徽、江西、山東、河南、湖北、湖南、四川等13個省(市、區(qū))；糧食主銷區(qū)為北京、天津、上海、浙江、福建、廣東、海南等7個省((市、區(qū))；糧食平衡區(qū)為山西、廣西、貴州、云南、陜西、新疆、甘肅、青海、寧夏等9個省(市、區(qū))。根據(jù)本文研究目的，選取土地投入、勞動力投入、機械投入、水資源投入、化肥投入和農(nóng)藥投入作為投入變量；選取糧食總產(chǎn)量為合意產(chǎn)出變量，糧食生產(chǎn)污染排放物為非合意產(chǎn)出變量，以上各變量相關數(shù)據(jù)均來自《中國農(nóng)業(yè)年鑒》和《中國農(nóng)村統(tǒng)計年鑒》。

2. 指標選取

(1)糧食生產(chǎn)投入變量。其中，土地投入以各年度各省(市、區(qū))的糧食作物播種面積來表示；勞動力投入以各年度各省(市、區(qū))的第一產(chǎn)業(yè)從業(yè)人員來表示；機械投入以各年度各省(市、區(qū))的農(nóng)用機械總動力來表示；水資源投入因1998年以前農(nóng)業(yè)用水數(shù)據(jù)缺失，因此以各年度各省(市、區(qū))的農(nóng)用土地有效灌溉面積來替代，包括能進行正常灌溉的水田和水澆地面積；化肥投入以各年度各省(市、區(qū))農(nóng)用化肥折純量來表示；農(nóng)藥投入以各年度各省(市、區(qū))的的農(nóng)藥施用量來表示。為了保證與糧食產(chǎn)出口徑的一致性，本文借鑒馬文杰等[16]的做法，使用權重系數(shù)法將糧食生產(chǎn)的要素投入從廣義的農(nóng)業(yè)投入中剝離出來，其中，權重系數(shù)A=(農(nóng)業(yè)產(chǎn)值/農(nóng)林牧漁總產(chǎn)值)×(糧食播種面積/農(nóng)作物播種總面積)，權重系數(shù)B=糧食播種面積/農(nóng)作物播種總面積。種糧土地投入仍然用糧食播種面積來表示，種糧勞動力投入用第一產(chǎn)業(yè)從業(yè)人員與系數(shù)A相乘可以得到，種糧機械投入、水資源投入、化肥投入和農(nóng)藥投入用相應的農(nóng)業(yè)投入與系數(shù)B相乘即可得到。

(2)糧食生產(chǎn)產(chǎn)出變量。包括合意產(chǎn)出和非合意產(chǎn)出。合意產(chǎn)出以各年度各個糧食主產(chǎn)區(qū)的糧食總產(chǎn)量來表示；非合意產(chǎn)出用兩個指標來核算，即總氮(TN )和總磷(TP)。盡管賴斯蕓等、陳敏鵬等[17-18]將農(nóng)業(yè)污染物指標定為總氮(TN)、總磷(TP) 和化學需氧量(CODcr)，但是張桂英等[19]的研究表明農(nóng)田排水中有機物質對水質污染很小，全國第一次污染源普查也得出種植業(yè)COD排放量不足農(nóng)業(yè)源COD總排放量的5%[20]。因此，本文未將化學需氧量納入非合意產(chǎn)出指標。糧食種植污染物排放量參照賴斯蕓等、陳敏鵬等[17-18]采用的清單分析法，具體計算公式如下所示：

(2)

式中，E為糧食種植所產(chǎn)生的污染排放量，主要包括總氮(TN )和總磷(TP)；EUi為單元i指標的統(tǒng)計數(shù)；ρij為單元i污染物j的排放強度；ηi為相關資源的利用效率系數(shù)；Cij為單元i污染物j的排放系數(shù)。EUi數(shù)據(jù)來自歷年相關統(tǒng)計年鑒，ρij、ηi和Cij等各參數(shù)是通過廣泛文獻調研及綜合比較所得，重點參考了第一次全國污染源普查領導小組辦公室《污染源普查農(nóng)業(yè)源系數(shù)手冊》中的不同省份各參數(shù)取值。本文的糧食作物統(tǒng)計采納稻谷、小麥、玉米、豆類和薯類5類糧食作物的相關數(shù)據(jù)，糧食作物草谷比根據(jù)可再生能源戰(zhàn)略研究組[21]，以及相關學者的研究成果[22]綜合比較整理得出。

三、測度結果分析

本文運用MAXDEA PRO 6.3專業(yè)版軟件對1992—2012年我國29個省(市、區(qū))的糧食生產(chǎn)環(huán)境技術效率進行測度，具體分析結果如下：

(一) 糧食生產(chǎn)環(huán)境技術效率的時空差異分析

從糧食生產(chǎn)環(huán)境技術效率變化的時間趨勢來看，1992—2012年我國糧食生產(chǎn)環(huán)境技術效率整體上呈波動式下降趨勢(見圖1)，從1992年的0.940下降為2012年的0.915。從其發(fā)展趨勢來看，可以分為4個階段：1992—1998年為波動式下降期，糧食生產(chǎn)環(huán)境技術效率從0.917下降到0.864，1998年跌落到最低點；1998—2001年為持續(xù)上升期，由0.864上升到0.922；2001—2005年為持續(xù)下降期，從0.922下降為0.866；2005—2012年為波動式上升期，從0.866上升為0.888。

圖1　中國糧食生產(chǎn)環(huán)境技術效率總體變化趨勢

圖2　中國各糧食生產(chǎn)功能區(qū)內糧食生產(chǎn)環(huán)境技術效率變化趨勢

從各糧食生產(chǎn)功能區(qū)來看，1992—2012年糧食主產(chǎn)區(qū)的糧食生產(chǎn)環(huán)境技術效率要高于糧食主銷區(qū)和糧食平衡區(qū)，這說明相對而言，糧食主產(chǎn)區(qū)更加注重糧食生產(chǎn)與環(huán)境的和諧發(fā)展。就具體變化趨勢來看，糧食主產(chǎn)區(qū)和糧食主銷區(qū)呈波動式下降趨勢，但糧食主銷區(qū)波動幅度要大于糧食主產(chǎn)區(qū)，主銷區(qū)的糧食生產(chǎn)環(huán)境技術效率從1992年的0.920下降為2012年的0.806，21年下降了14.76%，糧食主銷區(qū)大多是一些經(jīng)濟發(fā)達的省份，這說明隨著市場化改革的大力推進，城鎮(zhèn)化的加速發(fā)展，糧食主銷區(qū)的經(jīng)濟雖然得到了大力發(fā)展，但是糧食生產(chǎn)環(huán)境卻在逐步惡化。糧食平衡區(qū)的糧食生產(chǎn)環(huán)境技術效率在1999年前呈現(xiàn)出波動式下降趨勢，從1999年后開始波動式上升，這說明1999年在四川、陜西、甘肅率先啟動并于2000年正式全面啟動的退耕還林(草)試點工作取得了巨大成效，促進了糧食平衡區(qū)糧食生產(chǎn)環(huán)境的改善。

(二)省際糧食生產(chǎn)環(huán)境技術效率變化的協(xié)調分析

以上是從平均意義角度對全國以及各糧食生產(chǎn)功能區(qū)糧食生產(chǎn)環(huán)境技術效率變化特征進行的分析，可能會掩蓋各個省(市)在考察期的具體表現(xiàn)及演變特征，為此本文進一步根據(jù)各省(市)糧食生產(chǎn)環(huán)境技術效率的變化特征，對各省(市)糧食生產(chǎn)與環(huán)境的協(xié)調情況進行分析。

聚類分析法是根據(jù)相似性原則對事物進行分組，是一種能夠對事物的內在特點和規(guī)律進行分析的探索性分析方法，是數(shù)據(jù)挖掘中常用的一種技術[23]。本文將1992—2012年分為4個階段，運用SPSS軟件采用系統(tǒng)聚類分析方法，以離差平方和最小為分類依據(jù)，歐氏距離為距離函數(shù)，對各階段各省(市)歷年的糧食生產(chǎn)環(huán)境技術效率進行聚類分析，將29個省(市)分為較協(xié)調區(qū)、較不協(xié)調區(qū)和不協(xié)調區(qū)3類，分析結果如表1所示。

表1　我國29個省(市)糧食生產(chǎn)與環(huán)境協(xié)調分類結果

總體來看，在第一階段(1992—1996)時，我國各省(市)糧食生產(chǎn)與環(huán)境的協(xié)調性最好，處于較協(xié)調區(qū)的省(市)有17個，這些省(市)的糧食生產(chǎn)環(huán)境技術效率值相對較高。處于較不協(xié)調區(qū)的省(市)有9個，處于不協(xié)調區(qū)的省(市)有3個，這3個省(市)的環(huán)境技術效率值最低。到第二階段(1997—2001)和第三階段(2002—2006)時，各省(市)糧食生產(chǎn)與環(huán)境的協(xié)調性開始變差，處于較協(xié)調區(qū)的省市減少為12個，處于較不協(xié)調區(qū)和不協(xié)調區(qū)的省(市)增加到17個。到第4階段(2007—2012)時，各省(市)糧食生產(chǎn)與環(huán)境的協(xié)調性開始有所好轉，處于較協(xié)調區(qū)的省(市)數(shù)量恢復第一階段水平，但處于不協(xié)調區(qū)的省(市)數(shù)量仍然多于第一階段的數(shù)量。

從表中，我們發(fā)現(xiàn)不少糧食主產(chǎn)區(qū)的糧食生產(chǎn)與環(huán)境協(xié)調狀況日趨惡化，如冀、吉、黑、蘇、鄂、湘、川等7個糧食主產(chǎn)區(qū)，說明我國糧食主產(chǎn)區(qū)面臨的生態(tài)環(huán)境壓力逐漸加大，亟需轉變糧食生產(chǎn)方式。

(三)各地區(qū)糧食生產(chǎn)環(huán)境技術效率差距的動態(tài)演變特征分析

從上述分析可以看出，在各個階段各省(市)區(qū)之間糧食生產(chǎn)環(huán)境技術效率存在較大的差距，這種差距是否會隨著時間的推移而發(fā)生縮小呢？為此本文應用Kernel密度估計對各地區(qū)糧食生產(chǎn)環(huán)境技術效率差距的動態(tài)演變特征進行具體分析。

核密度估計是由雷科特與希爾沃曼[24-25]經(jīng)過深入研究發(fā)展而得到的一種重要的非參數(shù)估計方法。該方法的原理是，通過連續(xù)的密度曲線來描述隨機變量的分布形態(tài)，從數(shù)據(jù)本身出發(fā)研究其分布特征，可以避免預設函數(shù)形式而導致測算值與實際值產(chǎn)生差異。設隨機變量X的密度函數(shù)為f(x)，則在點x0的概率密度可以由公式(3)來估計：

(3)

式(3)中，N是樣本數(shù)據(jù)觀測值的個數(shù)；h是帶寬，也可稱為平滑參數(shù)，用于設定分布函數(shù)fh(x)的整體平滑特性，在核密度估計中是至關重要的參數(shù)，帶寬選擇過大或過小都會產(chǎn)生比較大的誤差，h值越大，平滑度就越大，帶寬的選擇在很大程度上決定了密度曲線的形態(tài)[26]。實際應用中，樣本越多，寬帶就要求越小，但也不能太小，應滿足公式(4)：

(4)

1. 全國29個省(市、區(qū))糧食生產(chǎn)環(huán)境技術效率的Kernel密度估計

圖3對全國29個省(市、區(qū))糧食生產(chǎn)環(huán)境技術效率在考察期內的分布進行了描述。就全國而言，糧食生產(chǎn)環(huán)境技術效率在考察期內一直呈現(xiàn)出“雙峰”分布狀態(tài)，其中，第一個主峰位于0.78附近，第二個主峰位于0.98附近，而且密度函數(shù)中心逐步向左轉移，波峰寬度也開始逐漸變寬。這說明全國29個省(市、區(qū))的糧食生產(chǎn)環(huán)境技術效率分布在1992—2012年呈現(xiàn)出兩級分化趨勢，只有少數(shù)地區(qū)(如黑龍江、吉林、內蒙古等)的糧食生產(chǎn)環(huán)境技術效率較高，其他大部分地區(qū)則位于低糧食生產(chǎn)環(huán)境技術效率區(qū)間，而且效率值在逐步降低，各省(市)效率差距也在慢慢擴大。

糧食生產(chǎn)環(huán)境的提升需要以經(jīng)濟實力作為支撐，而各省份的經(jīng)濟實力懸殊，對農(nóng)業(yè)廢棄物資源化利用和農(nóng)業(yè)面源污染防控等方面的態(tài)度以及采取的措施有所差異，加之不同時期各個地區(qū)氣候、地形、災害等因素的不一致，從而導致全國各個省份糧食生產(chǎn)環(huán)境技術效率差異顯著。

2. 糧食主產(chǎn)區(qū)糧食生產(chǎn)環(huán)境技術效率的Kernel密度估計

圖4對糧食主產(chǎn)區(qū)糧食生產(chǎn)環(huán)境技術效率在考察期內分布的動態(tài)特征進行了描述。我國13個糧食主產(chǎn)區(qū)的糧食生產(chǎn)環(huán)境技術效率在考察期內除了2007年呈“單峰”，其他時間都呈“雙峰”分布，在整個考察期第一個主峰位于0.98附近，第二個主峰位于0.83附近，密度函數(shù)中心逐漸向左移動，波動寬度逐漸變寬，峰值也在逐漸變小。這說明糧食主產(chǎn)區(qū)各省(市)的糧食生產(chǎn)環(huán)境技術效率分布基本上也呈兩級分化狀態(tài)，效率值在逐步變小，效率差距有所變大，如黑龍江、吉林等省份位于第一個主峰附近，河北、河南和安徽等省份的效率值相對偏低，而且逐漸降低，位于第二個主峰附近，兩大區(qū)域差距逐漸擴大。

在糧食主產(chǎn)區(qū)中，河北、河南和安徽等省份糧食生產(chǎn)環(huán)境技術效率不僅相對偏低，而且還出現(xiàn)降低趨勢，與之相比，黑龍江、吉林等省份的糧食生產(chǎn)環(huán)境技術效率一直為1，處于生產(chǎn)前沿面，從而造成主產(chǎn)區(qū)內部各個省份糧食生產(chǎn)環(huán)境技術效率呈現(xiàn)出兩級分化狀態(tài)。

圖3　全國29個省(市、區(qū))糧食生產(chǎn)環(huán)境技術效率的演進

圖4　主產(chǎn)區(qū)糧食生產(chǎn)環(huán)境技術效率的演進

圖5　主銷區(qū)糧食生產(chǎn)環(huán)境技術效率的演進

圖6　平衡區(qū)糧食生產(chǎn)環(huán)境技術效率的演進

3. 糧食主銷區(qū)糧食生產(chǎn)環(huán)境技術效率的Kernel密度估計

圖5對糧食主銷區(qū)糧食生產(chǎn)環(huán)境技術效率在考察期內分布的動態(tài)特征進行了描述。我國7個糧食主銷區(qū)的糧食生產(chǎn)環(huán)境技術效率在1997年呈“單峰”分布，其他時間段都呈“雙峰”分布狀態(tài)，主峰值由1992年的0.97降為2012年的0.75，密度函數(shù)中心大幅向左轉移，波動寬度逐年變大。這說明糧食主銷區(qū)各省(市)的糧食生產(chǎn)環(huán)境技術效率在考察期內基本上處于兩級分布狀態(tài)，而且效率值逐年降低，大部分地區(qū)的效率值分布從0.97附近降低到0.75附近，效率差距逐步擴大。

糧食主銷區(qū)中的大部分省(市)經(jīng)濟較為發(fā)達，不同的經(jīng)濟發(fā)展模式導致了不同的產(chǎn)業(yè)結構，從而對糧食生產(chǎn)環(huán)境也產(chǎn)生了差異性影響。如天津作為典型的工業(yè)城市，糧食生產(chǎn)環(huán)境技術效率持續(xù)走低；廣東的糧食生產(chǎn)環(huán)境技術效率在考察期的頭幾年比較高，但是隨著廣東以低附加值出口為導向的經(jīng)濟模式的轉型，農(nóng)業(yè)省份向工業(yè)省份的跨越，糧食生產(chǎn)環(huán)境也開始出現(xiàn)惡化的趨勢；高新技術產(chǎn)業(yè)比較發(fā)達的浙江和上海，對糧食生產(chǎn)環(huán)境則給予了高度重視。這一系列因素造成了糧食主銷區(qū)糧食生產(chǎn)環(huán)境技術效率分布兩極化格局的形成。

4. 糧食平衡區(qū)糧食生產(chǎn)環(huán)境技術效率的Kernel密度估計

圖6對糧食平衡區(qū)糧食生產(chǎn)環(huán)境技術效率在考察期內分布的動態(tài)特征進行了描述。我國9個糧食平衡區(qū)的糧食生產(chǎn)環(huán)境技術效率在考察期呈“雙峰”分布，主峰值由0.78提升到1，密度函數(shù)中心大幅右移，波動寬度逐步拉大。這說明總體來說，平衡區(qū)糧食生產(chǎn)環(huán)境技術效率在考察期內在逐步提升，但呈現(xiàn)出兩級分化狀態(tài)。貴州、青海和新疆由于自然環(huán)境優(yōu)越，環(huán)境污染比較小，糧食生產(chǎn)環(huán)境技術效率一直比較高，隨著退耕還林工程的實施，甘肅、陜西和寧夏的生態(tài)環(huán)境也得到了改善，而廣西、云南和山西等省份由于自然環(huán)境、生產(chǎn)條件惡化等原因出現(xiàn)了糧食生產(chǎn)環(huán)境技術效率降低，從而造成了糧食平衡區(qū)糧食生產(chǎn)環(huán)境技術效率兩級分化局面的形成。

四、 結論與建議

(一)主要研究結論

全國及各糧食生產(chǎn)功能區(qū)糧食生產(chǎn)環(huán)境狀況總體上呈現(xiàn)出逐漸惡化趨勢，但糧食主產(chǎn)區(qū)優(yōu)于糧食主銷區(qū)，糧食主銷區(qū)優(yōu)于平衡區(qū)，這說明糧食主產(chǎn)區(qū)糧食生產(chǎn)與環(huán)境的和諧度最高，糧食主銷區(qū)和平衡區(qū)的糧食生產(chǎn)造成的環(huán)境污染更大，未能實現(xiàn)經(jīng)濟效益和環(huán)境效益的雙贏。

中國糧食生產(chǎn)環(huán)境技術效率在區(qū)域上存在明顯的非均衡性。根據(jù)系統(tǒng)聚類分析方法，29個省(市)區(qū)被劃分為糧食生產(chǎn)與環(huán)境“較協(xié)調區(qū)”“較不協(xié)調區(qū)”和“不協(xié)調區(qū)”。1992—1996年，“較協(xié)調區(qū)”和“較不協(xié)調區(qū)”主要以主產(chǎn)區(qū)和主銷區(qū)各省(市)為主，“不協(xié)調區(qū)”以少數(shù)平衡區(qū)省份為主；1996年后，位于“較不協(xié)調區(qū)”的省(市)變少，主產(chǎn)區(qū)和主銷區(qū)部分省(市)的糧食生產(chǎn)與環(huán)境協(xié)調狀況有所惡化，逐漸由“較協(xié)調區(qū)”向“不協(xié)調區(qū)”轉移。

總體而言，我國糧食生產(chǎn)環(huán)境技術效率的地區(qū)差距在逐漸擴大，呈現(xiàn)出兩級分布, 各功能區(qū)有所差異。就各功能區(qū)而言，糧食主產(chǎn)區(qū)各省(市)的差距逐步縮小，而糧食主銷區(qū)和平衡區(qū)的差距則在逐步擴大，各功能區(qū)分布基本上呈兩級分化趨勢。

(二) 建議

1. 糧食利益補償機制要進一步向糧食主產(chǎn)區(qū)傾斜，同時糧食主銷區(qū)也要加大對糧食生產(chǎn)環(huán)境的治理

首先，要進一步加大中央財政對糧食主產(chǎn)區(qū)財政的轉移支付力度。和糧食主銷區(qū)相比，糧食主產(chǎn)區(qū)經(jīng)濟比較落后，財政實力較弱，而隨著國家支農(nóng)項目的日益增多，地方財政壓力越來越大，因此中央財政要從基礎設施建設、耕地質量保護、防治自然災害等方面進一步加大對糧食主產(chǎn)區(qū)的財政資金投入。其次，加快糧食主銷區(qū)對糧食主產(chǎn)區(qū)利益補償機制的實施。糧食主銷區(qū)根據(jù)其調入的糧食量以及相應的經(jīng)濟實力，從財政收入中拿出一部分補貼糧食主產(chǎn)區(qū)，糧食主產(chǎn)區(qū)則根據(jù)其調出的糧食量接受相應的補償。此外，面對糧食主銷區(qū)的環(huán)境技術效率逐年降低，要加大對糧食生產(chǎn)環(huán)境的治理，促進該地區(qū)經(jīng)濟增長和環(huán)境的和諧發(fā)展。

2. 加快糧食生產(chǎn)技術創(chuàng)新，促進糧食生產(chǎn)的可持續(xù)發(fā)展

傳統(tǒng)的依靠化肥農(nóng)藥投入的粗放型糧食生產(chǎn)方式雖然為我國糧食產(chǎn)量的增長做出了巨大貢獻，但是也給我國的糧食生產(chǎn)環(huán)境造成了巨大破壞，為促進我國糧食生產(chǎn)的可持續(xù)發(fā)展，需要改革原有的糧食生產(chǎn)方式，加大對糧農(nóng)科學種糧知識的培訓，提高化肥農(nóng)藥的利用率，并加快有機肥、緩釋肥料、生態(tài)農(nóng)藥以及安全低殘留農(nóng)藥等低污染環(huán)保型生產(chǎn)資料的試點和推廣。針對目前秸稈焚燒現(xiàn)象嚴重而行政監(jiān)管治標不治本的問題，應加快秸稈還田以及秸稈能源化利用等農(nóng)業(yè)環(huán)保技術的示范和推廣，并對農(nóng)民的環(huán)保技術利用行為給予相應的優(yōu)惠和補貼，提高農(nóng)民利用農(nóng)業(yè)環(huán)保技術的積極性。

3. 對于糧食生態(tài)環(huán)境惡化比較嚴重的“不協(xié)調區(qū)”，要加快對重污染農(nóng)業(yè)水土資源的修復和治理

對于這些地區(qū)來說，生態(tài)環(huán)境惡化已經(jīng)成為制約其糧食生產(chǎn)可持續(xù)發(fā)展的一個瓶頸，干旱、土壤沙漠化和水土流失等問題不斷降低耕地的生產(chǎn)能力，當務之急是停止毀林開荒等破壞性行為，大力開展退耕還林，加快對鹽堿地、沙漠化土壤以及重污染水土資源的修復和治理。從短期來看，生態(tài)保護政策的實施可能會減少糧食播種面積，降低糧食產(chǎn)量，但是，從長遠來看，可以涵養(yǎng)水源、抑制土壤沙漠化、固定植被、提高生態(tài)環(huán)境的承載力，從而促進糧食生態(tài)環(huán)境系統(tǒng)的良性循環(huán)。

參考文獻：

[1]何可，張俊飚，豐軍輝. 自我雇傭型農(nóng)村婦女的農(nóng)業(yè)技術需求意愿及其影響因素分析——以農(nóng)業(yè)廢棄物基質化產(chǎn)業(yè)技術為例[J].中國農(nóng)村觀察，2014, (4)：84-94.

[2]何可, 張俊飚, 蔣磊. 生物質資源減碳化利用需求及影響機理實證研究——基于SEM模型分析方法和TAM理論分析框架[J].資源科學，2013, 35(8): 1635-1642.

[3]蔣立科, 鄒文政等. 糧食產(chǎn)量的九年連增給人留下的思考[J].今日中國論壇，2013, (Z1)：47-50.

[4]涂正革, 諶仁俊. 傳統(tǒng)方法測度的環(huán)境技術效率低估了環(huán)境治理效率?——來自基于網(wǎng)絡DEA的方向性環(huán)境距離函數(shù)方法分析中國工業(yè)省級面板數(shù)據(jù)的證據(jù)[J].經(jīng)濟評論, 2013, (5): 89-99.

[5]陳詩一. 能源消耗、二氧化碳排放與中國工業(yè)的可持續(xù)發(fā)展[J].經(jīng)濟研究,2009, 44(4): 41-55.

[6]SCHEEL H. Undesirable Outputs in Efficiency Valuations[J]. European Journal of Operational Research, 2001, 132(2): 400-410．

[7]薛建良, 李秉龍. 基于環(huán)境修正的中國農(nóng)業(yè)全要素生產(chǎn)率度量[J].中國人口·資源與環(huán)境,2011, 21(5): 113-118.

[8]涂正革. 環(huán)境、資源與工業(yè)增長的協(xié)調性[J].經(jīng)濟研究, 2008, (2): 93-105.

[9]王海寧. FDI、出口密集度與環(huán)境技術效率——基于產(chǎn)業(yè)數(shù)據(jù)的分位數(shù)回歸方法分析[J]. 財貿(mào)研究, 2011, (2): 24-32, 129.

[10]李偉娜, 金曉雨. 中國制造業(yè)的環(huán)境技術效率研究[J].中國科技論壇, 2011, (2):33-38.

[11]李谷成, 范麗霞，閔銳.資源、環(huán)境與農(nóng)業(yè)發(fā)展的協(xié)調性——基于環(huán)境規(guī)制的省級農(nóng)業(yè)環(huán)境效率排名[J].數(shù)量經(jīng)濟技術經(jīng)濟研究, 2011, (10): 21-36, 49.

[12]吳麗麗, 鄭炎成, 李谷成. 碳排放約束下我國油菜全要素生產(chǎn)率增長與分解——來自 13個主產(chǎn)區(qū)的實證[J]. 農(nóng)業(yè)現(xiàn)代化研究, 2013, (1): 77-80, 94.

[13]CHUNG Y H, FARE R, GROSSKOPF S. Productivity and Undesirable Outputs: A Directional Distance Function Approach [J]. Journal of Environmental Management, 1997, 51 (3): 229-240.

[14]鐘甫寧，劉順飛. 中國水稻生產(chǎn)布局[J].中國農(nóng)村經(jīng)濟，2007, (9): 39-44.

[15]陸文聰，黃祖輝. 中國糧食供求變化趨勢預測：基于區(qū)域化市場均衡模型[J].經(jīng)濟研究，2004, (8)：94-104.

[16]馬文杰. 中國糧食綜合生產(chǎn)能力研究[M]. 北京：科學出版社，2010：6-59.

[17]賴斯蕓，杜鵬飛，陳吉寧. 基于單元分析的非點源污染調查評估方法[J].清華大學學報(自然科學版)，2004, 44(9): 1184-1187.

[18]陳敏鵬，陳吉寧，賴斯蕓．中國農(nóng)業(yè)和農(nóng)村污染的清單分析與空間特征識別[J]. 中國環(huán)境科學，2007，26(6): 751-755.

[19]張桂英, 汪祖強. 蘇南太湖地區(qū)農(nóng)田水中COD的調查研究[J].農(nóng)村生態(tài)環(huán)境，1994，10(2): 41-44.

[20]王輝, 欒維新，康敏捷. 遼河流域社會經(jīng)濟活動的COD污染負荷[J].地理研究，2013，32(10)：1802-1813.

[21]中國可再生能源發(fā)展戰(zhàn)略研究項目組.中國可再生能源發(fā)展戰(zhàn)略研究叢書: 生物能源卷[M]. 北京: 中國電力出版社，2008.

[22]王亞靜，畢于運，高春雨.中國秸稈資源可收集利用量及其適宜性評價[J]. 中國農(nóng)業(yè)科學，2010， 43(9): 1852-1859.

[23]郝春旭，董戰(zhàn)峰，葛察忠，等.基于聚類分析法的省級環(huán)境績效動態(tài)評估與分析[J].生態(tài)經(jīng)濟，2015,31(1):154-157.

[24]PLACKETT R L. An Introduction to the Theory of Statistics [M]. Edinburgh: Oliver and Boyd,1971.

[25]SILVERMAN B W. Density Estimation for Statistics and Data Analysis [M]. London: Chapman and Hall, 1986.

[26]文洋. 收入分配對我國出口貿(mào)易的影響[J].世界經(jīng)濟研究，2011,(10): 33-39.

Measurement and Analysis of Environmental Technology Efficiency Based on Grain Production

ZHAO Li-pinga,b,WANG Ya-penga,b,HE Kea,b

(a.CollegeofEconomicsandManagement,HuazhongAgriculturalUniversity,Wuhan430070,China；b.HubeiRuralDevelopmentResearchCenter,Wuhan430070,China)

Abstract：Using the directional distance function and Kernel density estimation method, this paper calculated and analyzed the environmental technology efficiency of grain production of 29 provinces (cities/areas) in China from 1992 to 2012, so as to investigate the coordination between the grain production and resources environment during the acceleration period of urbanization and industrialization. The research results show that: (1) The grain production environment of every functional area presents a gradual deterioration in general; Among which, the main grain producing area is better than the main selling area, which in turn is better than the balanced area; (2) The environmental technical efficiency of grain production in China differentiates in each province; (3) Overall, the regional disparity of the environmental technical efficiency of grain production gradually expands, showing a two level distribution trend, but the situation is different in each functional area of grain production.

Key Words：grain production; environmental technical efficiency; environmental pollution; production goods

中圖分類號：F307.11

文獻標識碼：A

文章編號：1672-0202(2016)03-0028-10

作者簡介:趙麗平(1984—)，女，湖北荊門人，華中農(nóng)業(yè)大學經(jīng)濟管理學院博士研究生，主要研究方向為農(nóng)業(yè)資源保護開發(fā)、糧食安全。E-mail:ccnuzhao@qq.com

基金項目:國家社會科學基金重點項目(13ACL004)；國家自然科學基金青年項目(71303243)；國家自然科學基金青年項目(71203169)。

收稿日期:2016-01-24

DOI:10.7671/j.issn.1672-0202.2016.03.004