郭永奇，侯林岐

（1.南陽(yáng)理工學(xué)院，河南南陽(yáng) 473000；2.石河子大學(xué)經(jīng)濟(jì)與管理學(xué)院，新疆石河子 832003）

作為重要戰(zhàn)略物資，糧食安全問題事關(guān)國(guó)計(jì)民生,是國(guó)家安全、社會(huì)穩(wěn)定與經(jīng)濟(jì)發(fā)展的重要基礎(chǔ)和保障。作為我國(guó)糧食生產(chǎn)的重要基地和糧食生產(chǎn)的核心區(qū)域,糧食主產(chǎn)區(qū)在保證國(guó)家糧食生產(chǎn)安全和有效供給方面具有重要的戰(zhàn)略地位。2017 年，全國(guó)糧食總產(chǎn)量61 791 萬(wàn)t，其中13 個(gè)糧食主產(chǎn)區(qū)糧食產(chǎn)量為47 076 萬(wàn)t，占全國(guó)糧食總產(chǎn)量的76.2%，可見，糧食主產(chǎn)區(qū)是保障糧食有效供給和實(shí)現(xiàn)糧食安全的關(guān)鍵，其直接影響著國(guó)家的穩(wěn)定與發(fā)展。為實(shí)現(xiàn)糧食產(chǎn)量和農(nóng)業(yè)產(chǎn)值的不斷攀升，各種新技術(shù)和化學(xué)投入品在農(nóng)產(chǎn)品的研發(fā)、種植、加工和流通環(huán)節(jié)被廣泛使用。這在帶來(lái)我國(guó)糧食產(chǎn)量逐年提高的同時(shí)，也極大地消耗了我國(guó)原本稀缺的生產(chǎn)資源，并對(duì)農(nóng)業(yè)環(huán)境造成了長(zhǎng)期不可逆轉(zhuǎn)的嚴(yán)重污染和破壞［1］。隨著資源環(huán)境約束，在現(xiàn)有的單產(chǎn)水平、播種面積以及勞動(dòng)人口結(jié)構(gòu)等農(nóng)業(yè)資源條件下，多數(shù)糧食主產(chǎn)區(qū)已經(jīng)接近產(chǎn)量極限，如何加快轉(zhuǎn)變農(nóng)業(yè)發(fā)展方式，確保糧食等重要農(nóng)產(chǎn)品有效供給，實(shí)現(xiàn)綠色發(fā)展和資源永續(xù)利用，是必須破解的現(xiàn)實(shí)難題［2］。

綠色農(nóng)業(yè)作為一種農(nóng)業(yè)可持續(xù)發(fā)展的發(fā)展模式，既是農(nóng)業(yè)現(xiàn)代化的必然方向，也是實(shí)施鄉(xiāng)村振興戰(zhàn)略的重要措施，還是解決人民群眾日益增長(zhǎng)的綠色食品和優(yōu)質(zhì)生態(tài)產(chǎn)品需求之間矛盾的現(xiàn)實(shí)途徑。因此，發(fā)展綠色農(nóng)業(yè)就需要了解我國(guó)尤其是糧食主產(chǎn)區(qū)的綠色農(nóng)業(yè)發(fā)展現(xiàn)狀，其農(nóng)業(yè)綠色全要素生產(chǎn)率處于何種水平?其演變有何顯著特征?受哪些主要因素影響?探討這些問題，對(duì)于識(shí)別農(nóng)業(yè)綠色發(fā)展的動(dòng)力源泉和區(qū)域差異的原因，科學(xué)研判我國(guó)糧食主產(chǎn)區(qū)現(xiàn)代農(nóng)業(yè)發(fā)展水平，探尋中國(guó)農(nóng)業(yè)綠色協(xié)調(diào)發(fā)展的模式均具有非常重要的實(shí)踐和參考價(jià)值。

學(xué)界關(guān)于農(nóng)業(yè)綠色全要素生產(chǎn)率主要體現(xiàn)在測(cè)算指標(biāo)、測(cè)算方法和研究的空間尺度三個(gè)方面。在測(cè)算指標(biāo)上，主要體現(xiàn)在對(duì)環(huán)境污染因素的處理上，Zhiyang Shen 等［3］把環(huán)境污染因素視為一種非期望產(chǎn)出變量；Reinhard 等［4］將其視為一種投入變量。在測(cè)算方法上，學(xué)界多偏向采用隨機(jī)前沿分析法(SFA),數(shù)據(jù)包絡(luò)分析法(DEA)等方法測(cè)度農(nóng)業(yè)綠色全要素生產(chǎn)率，杜江等［5］采用GML 指數(shù)方法、梁俊等［6］基于DEA 模型和Luenberger 指數(shù)構(gòu)建農(nóng)業(yè)綠色全要素生產(chǎn)率指數(shù)。由于傳統(tǒng)DEA 模型的松弛問題，基于松弛變量的SBM(Slacks-Based Measure)方向性距離函數(shù)近年來(lái)被學(xué)者廣泛采用［7］。在研究的空間尺度上，國(guó)內(nèi)現(xiàn)有的研究成果大體為分為三類，一是薛建良等［8］、崔曉等［9］對(duì)全國(guó)的農(nóng)業(yè)綠色全要素生產(chǎn)率進(jìn)行研究;二是李文華等［10］、李欠男等［11］、李婷［12］對(duì)省級(jí)層面的農(nóng)業(yè)綠色全要素生產(chǎn)率進(jìn)行研究；三是針對(duì)特定區(qū)域性的研究，如徐晶晶［13］對(duì)沿海地區(qū)、劉芳等［14］對(duì)絲綢之路經(jīng)濟(jì)帶、吳傳清等［15］對(duì)長(zhǎng)江經(jīng)濟(jì)帶等區(qū)域的農(nóng)業(yè)綠色全要素生產(chǎn)率研究。學(xué)界對(duì)于糧食主產(chǎn)區(qū)農(nóng)業(yè)問題研究主要集中在糧食主產(chǎn)區(qū)的全要素生產(chǎn)率測(cè)度、利益補(bǔ)償機(jī)制、土地利用及生態(tài)環(huán)境效率測(cè)評(píng)等。關(guān)于糧食主產(chǎn)區(qū)的農(nóng)業(yè)綠色全要素的研究尚不多見。

基于此，本研究采用考慮農(nóng)業(yè)碳排放的全要素生產(chǎn)率模型，通過構(gòu)建13 個(gè)糧食主產(chǎn)區(qū)的面板數(shù)據(jù)模型，來(lái)對(duì)其農(nóng)業(yè)綠色全要素生產(chǎn)率及其變化規(guī)律進(jìn)行分析，進(jìn)而找出其影響因素，以期為我國(guó)農(nóng)業(yè)綠色發(fā)展提供政策決策的依據(jù)。

1 研究方法與數(shù)據(jù)來(lái)源

1.1 研究方法

準(zhǔn)確測(cè)算GTFP 是分析其影響因素的必要前提條件。DEA 方法因能解決多投入多產(chǎn)出的效率測(cè)量問題，更好的模擬實(shí)際生產(chǎn)中多投入多產(chǎn)出過程，近年來(lái)在GTFP 測(cè)度方面得到廣泛應(yīng)用，如李谷成等［16］、王軍等［17］運(yùn)用經(jīng)典非參數(shù)DEA 來(lái)對(duì)農(nóng)業(yè)GTFP 進(jìn)行測(cè)度。但傳統(tǒng)DEA 的Malmquist 指數(shù)無(wú)法測(cè)量含有非期望產(chǎn)出的GTFP，Chung 等［18］基于方向距離函數(shù)（DDF）提出了Malmquist-Luenberger(ML)指數(shù)，該方法能夠有效解決傳統(tǒng)DEA測(cè)量模型對(duì)于非期望產(chǎn)出的不得當(dāng)處理問題［19］。然而基于徑向或角度的DDF 模型仍然存在缺陷，一方面基于徑向距離函數(shù)只能保障期望產(chǎn)出與非期望產(chǎn)出變化比例相同，但當(dāng)投入或產(chǎn)出存在非零松弛時(shí)，基于徑向距離方法的測(cè)算結(jié)果會(huì)高估生產(chǎn)效率［20］；另一方面，基于角度距離函數(shù)需要對(duì)測(cè)度角度進(jìn)行選擇，無(wú)論是基于產(chǎn)出角度還是投入角度的測(cè)算都會(huì)因忽略另一角度而造成測(cè)算結(jié)果有偏［20］。因此，為了克服DDF 模型存在的缺陷，本文采用基于非徑向、非角度方向性距離函數(shù)的松弛效率模型（SBM）、ML指數(shù)來(lái)測(cè)量糧食主產(chǎn)區(qū)的農(nóng)業(yè)綠色全要素生產(chǎn)率。

本文將每個(gè)糧食主產(chǎn)區(qū)作為一個(gè)生產(chǎn)決策單元，假設(shè)每個(gè)產(chǎn)區(qū)的農(nóng)業(yè)產(chǎn)生產(chǎn)存在N種投入可生產(chǎn)出Q 種期望產(chǎn)出和L 種非期望產(chǎn)出參考Tone 等［21］研究方法，構(gòu)建包含期望產(chǎn)出和非期望產(chǎn)出的非徑向、非角度SBM 方向性模型函數(shù)為：

1.2 變量選擇及來(lái)源

(1)投入變量。根據(jù)現(xiàn)有文獻(xiàn)研究，本文選擇土地、勞動(dòng)力、機(jī)械、水資源、柴油、化肥、和農(nóng)膜等因素作為投入變量。1）土地投入：用農(nóng)作物播種面積（103 hm2）表示；2）勞動(dòng)力投入：用種植業(yè)勞動(dòng)力投入（萬(wàn)人）表示；3）機(jī)械投入：用種植業(yè)機(jī)械總動(dòng)力（104 KW）表示；4）水資源投入：用有效灌溉面積（103 hm2）表示；5）柴油投入：用農(nóng)用柴油使用量（104 t）表示；6）化肥投入：用農(nóng)用化肥施用折純量（104 t）表示；7）農(nóng)業(yè)投入：用農(nóng)藥使用量（t）表示；8）農(nóng)膜投入：用農(nóng)用塑料薄膜（t）表示。

（2）期望產(chǎn)出變量。農(nóng)業(yè)產(chǎn)出變量選擇我國(guó)13 個(gè)糧食主產(chǎn)區(qū)2000—2017 年的農(nóng)業(yè)總產(chǎn)值，采用各省份農(nóng)產(chǎn)品生產(chǎn)價(jià)格總指數(shù)將農(nóng)業(yè)總產(chǎn)值調(diào)整為以2000 年不變價(jià)格計(jì)的數(shù)據(jù)，以消除價(jià)格波動(dòng)因素帶來(lái)的影響。

（3）非期望產(chǎn)出變量。非期望產(chǎn)出是指在農(nóng)業(yè)生產(chǎn)中產(chǎn)生的污染排放，本文參考吳傳清等［15］研究成果，選擇農(nóng)業(yè)碳排放量(104 t)作為非期望產(chǎn)出變量。本文根據(jù)投入產(chǎn)出指標(biāo)數(shù)據(jù)，利用軟件MAXDEA6.0 對(duì)我國(guó)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率進(jìn)行測(cè)算，將糧食主產(chǎn)區(qū)13 個(gè)省份的農(nóng)業(yè)綠色全要素生產(chǎn)率幾何平均數(shù)作為當(dāng)年糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率。相關(guān)數(shù)據(jù)來(lái)源于2000—2017 年的《中國(guó)統(tǒng)計(jì)年鑒》《中國(guó)農(nóng)業(yè)統(tǒng)計(jì)年鑒》以及部分省市統(tǒng)計(jì)年鑒。

2 測(cè)度結(jié)果與分析

2.1 全國(guó)糧食主產(chǎn)區(qū)層面下農(nóng)業(yè)綠色全要素生產(chǎn)效率

從全國(guó)糧食主產(chǎn)區(qū)總體農(nóng)業(yè)GTFP 的數(shù)值和變化情況來(lái)看（圖1），我國(guó)糧食產(chǎn)出區(qū)的農(nóng)業(yè)GTFP年均值為1.029，除2001—2002 年、2002—2003 年和2013—2014 年的農(nóng)業(yè)GTFP 均值小于1，出現(xiàn)惡化情況之外，其他年份在碳排放約束下的農(nóng)業(yè)GTFP均值均大于1，說明在大多數(shù)年份中我國(guó)糧食主產(chǎn)區(qū)的農(nóng)業(yè)綠色全要素生產(chǎn)效率都存在一定程度的改善。從全國(guó)糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 的波動(dòng)趨勢(shì)來(lái)看，可以分為五個(gè)階段。具體來(lái)看：（1）第一個(gè)階段為2002—2006 年，我國(guó)糧食主產(chǎn)區(qū)的農(nóng)業(yè)GTFP 呈現(xiàn)上升態(tài)勢(shì)，其原因可能在該期間，國(guó)家發(fā)布環(huán)保專項(xiàng)通知，對(duì)違反排污企業(yè)進(jìn)行深入整治，從而降低污染排放，促使農(nóng)業(yè)GTFP 得到改善；（2）第二個(gè)階段為2006—2007 年，糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 達(dá)到“峰值”，究其原因，2005 年我國(guó)將環(huán)境規(guī)制正式納入“十一五”規(guī)劃，對(duì)單位GDP 碳排放進(jìn)行限制性約束，在各種污染物約束下，糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 在2006—2007 年中得到顯著提高；（3）2008—2009 年糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 出現(xiàn)下降，該現(xiàn)象出現(xiàn)的原因可能在于受到金融危機(jī)的影響，經(jīng)濟(jì)增長(zhǎng)受到衰退的威脅，需求不足導(dǎo)致GTFP 中的期望產(chǎn)出下降，使得糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 受到波動(dòng)；（4）2009—2010 年我國(guó)政府為應(yīng)對(duì)金融危機(jī)，實(shí)施一系列經(jīng)濟(jì)刺激計(jì)劃，使得中國(guó)經(jīng)濟(jì)實(shí)現(xiàn)快速增長(zhǎng)，GTFP 中期望產(chǎn)出GDP 增長(zhǎng)，糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP再次達(dá)到“峰值”；（5）2013—2015 年為第五個(gè)階段，隨著我國(guó)經(jīng)濟(jì)進(jìn)入新常態(tài)，伴隨著產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型和經(jīng)濟(jì)增速減緩，糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP再次落入“低谷”。

從糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 的分解項(xiàng)上來(lái)看，碳排放約束下的農(nóng)業(yè)GTFP 增長(zhǎng)主要依靠綠色技術(shù)進(jìn)步的推動(dòng)，除2002—2003 年綠色技術(shù)進(jìn)步的推動(dòng)作用出現(xiàn)惡化外，綠色技術(shù)進(jìn)步指數(shù)對(duì)于農(nóng)業(yè)GTFP的作用一直領(lǐng)先于綠色技術(shù)效率，而綠色技術(shù)效率在其增速落后于綠色技術(shù)進(jìn)步的情況下，出現(xiàn)了持續(xù)惡化現(xiàn)象。具體來(lái) 看：（1）綠色技術(shù)進(jìn)步指數(shù)年均值為1.03，除2002—2003 年外，技術(shù)進(jìn)步指數(shù)始終大于1，在經(jīng)過2000—2008 年上升趨勢(shì)后，在2008 年達(dá)到“瓶頸”狀態(tài)，隨后開始波動(dòng)下降；（2）綠色技術(shù)效率指數(shù)年均值為0.999，在大部分時(shí)期都處于小于1 的狀態(tài),2004—2010 年波動(dòng)上升，隨后進(jìn)入下滑。當(dāng)技術(shù)效率增長(zhǎng)速度長(zhǎng)期低于技術(shù)進(jìn)步指數(shù)時(shí)，資源配置會(huì)向技術(shù)進(jìn)步效率高的領(lǐng)域傾斜，造成資源配置失衡，導(dǎo)致生產(chǎn)效率會(huì)出現(xiàn)惡化，這也與唐家龍的研究觀點(diǎn)一致［22］。

圖1 2000—2017 年中國(guó)糧食主產(chǎn)區(qū)綠色全要素生產(chǎn)率指數(shù)

2.2 糧食主產(chǎn)區(qū)下農(nóng)業(yè)綠色全要素生產(chǎn)效率

從各個(gè)糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 情況上來(lái)看（圖2），位于東、中、西部不同位置的糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 水平差異明顯，各糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 發(fā)展驅(qū)動(dòng)因素存在差異。

圖2 2000—2017 年各糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP、EC、TC 對(duì)比分析

具體來(lái)看：（1）在區(qū)域差異上，東部地區(qū)農(nóng)業(yè)GTFP、技術(shù)進(jìn)步指數(shù)和技術(shù)效率水平均大于全國(guó)糧食主產(chǎn)區(qū)水平，2000—2017 年?yáng)|部的綠色全要素生產(chǎn)率平均增長(zhǎng)率為6.84%，遠(yuǎn)大于中部和西部，是農(nóng)業(yè)綠色經(jīng)濟(jì)發(fā)展的排頭兵；從增長(zhǎng)驅(qū)動(dòng)因素上來(lái)看，東部地區(qū)的農(nóng)業(yè)GTFP 增長(zhǎng)是農(nóng)業(yè)綠色技術(shù)進(jìn)步和技術(shù)效率“雙元”驅(qū)動(dòng)型的集約型發(fā)展，在各個(gè)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色經(jīng)濟(jì)發(fā)展中扮演著“引領(lǐng)者”和“技術(shù)革新者”的角色；中部地區(qū)農(nóng)業(yè)綠色經(jīng)濟(jì)發(fā)展主要依靠綠色技術(shù)進(jìn)步，是對(duì)東部地區(qū)農(nóng)業(yè)綠色經(jīng)濟(jì)發(fā)展的“追逐者”和“模仿者”，通過發(fā)揮后發(fā)優(yōu)勢(shì)，實(shí)現(xiàn)綠色技術(shù)進(jìn)步，從而促進(jìn)農(nóng)業(yè)GTFP效率的提高，而西部則處于技術(shù)進(jìn)步和技術(shù)效率雙項(xiàng)“塌陷”狀態(tài)，缺乏農(nóng)業(yè)綠色經(jīng)濟(jì)增長(zhǎng)動(dòng)力。（2）在省份差異上，河北（1.135）和黑龍江（1.081）兩省的農(nóng)業(yè)綠色全要素生產(chǎn)效率領(lǐng)先于其他糧食主產(chǎn)區(qū)，江西、吉林和內(nèi)蒙古三省農(nóng)業(yè)綠色全要素生產(chǎn)效率則處于惡化狀態(tài)；從增長(zhǎng)驅(qū)動(dòng)因素上來(lái)看，只有黑龍江、湖北、四川、江蘇、遼寧和山東6 個(gè)糧食主產(chǎn)區(qū)的農(nóng)業(yè)GTFP 增長(zhǎng)動(dòng)力來(lái)源于農(nóng)業(yè)綠色技術(shù)進(jìn)步和技術(shù)效率的“雙元”驅(qū)動(dòng)；河北、河南、湖南和安徽4 個(gè)糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 推動(dòng)主要來(lái)源于農(nóng)業(yè)綠色技術(shù)進(jìn)步，而農(nóng)業(yè)綠色技術(shù)效率則在不同程度上存在退步；糧食主產(chǎn)區(qū)江西的農(nóng)業(yè)GTFP驅(qū)動(dòng)主要來(lái)源于技術(shù)效率的提高，其技術(shù)進(jìn)步水平則存在一定倒退，吉林和內(nèi)蒙古2 個(gè)糧食主產(chǎn)區(qū)則處于綠色技術(shù)進(jìn)步和技術(shù)效率雙項(xiàng)低水平狀態(tài)，缺乏農(nóng)業(yè)GTFP 增長(zhǎng)動(dòng)力，說明各糧食主產(chǎn)區(qū)之間農(nóng)業(yè)技術(shù)創(chuàng)新差距不斷擴(kuò)大，技術(shù)創(chuàng)新主要集中在東部地區(qū)的糧食主產(chǎn)區(qū)，中、西部地區(qū)的糧食主產(chǎn)區(qū)面臨較大的資源環(huán)境協(xié)調(diào)發(fā)展壓力，農(nóng)業(yè)綠色經(jīng)濟(jì)與區(qū)域協(xié)調(diào)發(fā)展道路任重而道遠(yuǎn)。

3 糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率影響因素實(shí)證分析

3.1 變量選取

影響各糧食主產(chǎn)區(qū)的糧食全要素生產(chǎn)率的因素有很多，因農(nóng)業(yè)生產(chǎn)具有地域性特點(diǎn)，不同地方對(duì)其影響因素所選取的指標(biāo)也略有不同，對(duì)已有研究農(nóng)業(yè)綠色全要素影響因素進(jìn)行歸納和總結(jié)后，兼顧數(shù)據(jù)可得性，本文選取糧食作物受災(zāi)水平、灌溉設(shè)施水平、機(jī)械化水平、單位農(nóng)地化肥、農(nóng)藥、農(nóng)膜負(fù)荷指數(shù)、勞均經(jīng)營(yíng)規(guī)模、第二三產(chǎn)業(yè)發(fā)展水平、財(cái)政支出水平等7 個(gè)變量作為影響糧食主產(chǎn)區(qū)綠色全要素生產(chǎn)率的解釋變量。

（1）糧食作物受災(zāi)水平（FD/%）:該指標(biāo)用農(nóng)作物受災(zāi)水平是指各糧食主產(chǎn)區(qū)農(nóng)作物受災(zāi)面積占其農(nóng)作物總播種面積的比重來(lái)表示。該指標(biāo)值越大，表明受災(zāi)程度越嚴(yán)重，進(jìn)而對(duì)提高糧食生產(chǎn)效率和綠色全要素生產(chǎn)率越不利。

（2）灌溉設(shè)施水平（IF/103 hm2）:該指標(biāo)用各糧食主產(chǎn)區(qū)的有效灌溉面積來(lái)表示。其水平越高，越有利于優(yōu)化農(nóng)業(yè)資源配置，進(jìn)而對(duì)提高農(nóng)業(yè)生產(chǎn)率越有利，反之，會(huì)對(duì)農(nóng)業(yè)綠色全要素率不利。

（3）機(jī)械化水平（MI/104 KW）:該指標(biāo)用各糧食主產(chǎn)區(qū)的單位農(nóng)作物播種面積的農(nóng)用機(jī)械總動(dòng)力表示。其水平越高，越有利于提高農(nóng)業(yè)生產(chǎn)率。

（4）單位農(nóng)地化肥、農(nóng)藥、農(nóng)膜負(fù)荷指數(shù)（LI/%）：該指標(biāo)用單位農(nóng)作物面積所使用化肥、農(nóng)藥和農(nóng)膜的綜合計(jì)算值來(lái)表示。該指數(shù)越大，表明農(nóng)業(yè)生產(chǎn)對(duì)土地生態(tài)系統(tǒng)的壓力越大，對(duì)綠色農(nóng)業(yè)發(fā)展越不利。計(jì)算單位土地化肥、農(nóng)藥、農(nóng)膜負(fù)荷指數(shù)時(shí)，為消除化肥、農(nóng)藥、農(nóng)膜之間不同污染方式的差別以便于相互間的綜合計(jì)算，首先對(duì)各因子進(jìn)行標(biāo)準(zhǔn)化。采用下列方法計(jì)算：

單位土地化肥、農(nóng)藥、農(nóng)膜負(fù)荷指數(shù)= 0.5P1+0.3P2+0.2P3。

式（2）中，xi（i=1,2,3）為第i個(gè)評(píng)價(jià)指標(biāo)的指標(biāo)值，P(xi)為第i個(gè)指標(biāo)的單位農(nóng)作物面積農(nóng)業(yè)污染指數(shù)的標(biāo)準(zhǔn)值，0 ≤P(xi)≤1,xsi為評(píng)價(jià)指標(biāo)的全國(guó)的平均值。進(jìn)而以0.5,0.3,0.2 的權(quán)重進(jìn)行加權(quán)平均計(jì)算各因子的綜合指數(shù)，最終計(jì)算出各糧食主產(chǎn)區(qū)單位農(nóng)作物面積農(nóng)地化肥、農(nóng)藥、農(nóng)膜負(fù)荷指數(shù)。

（5）勞均經(jīng)營(yíng)規(guī)模(XS/人/hm2):該指標(biāo)用各糧食主產(chǎn)區(qū)種糧勞動(dòng)力數(shù)量與其糧食播種面積之比來(lái)表示。其中，種糧勞動(dòng)力數(shù)量用農(nóng)業(yè)從業(yè)人員乘以糧食作物面積與農(nóng)作物面積之比來(lái)表示［23］。該指標(biāo)越小表明勞動(dòng)生產(chǎn)效率越高，進(jìn)而全要素生產(chǎn)率也就越高。

（6）第二三產(chǎn)業(yè)發(fā)展水平(ST/%):該指標(biāo)用各個(gè)糧食主產(chǎn)區(qū)的第二三產(chǎn)業(yè)增加值占其國(guó)內(nèi)生產(chǎn)總值比重表示。從理論來(lái)看，第二、三產(chǎn)業(yè)的發(fā)展會(huì)對(duì)農(nóng)業(yè)生產(chǎn)要素的流出產(chǎn)生一定的虹吸效應(yīng)，這對(duì)農(nóng)業(yè)發(fā)展會(huì)有一定不利影響，但另一方面，其也會(huì)通過技術(shù)的外溢效應(yīng)促進(jìn)農(nóng)業(yè)生產(chǎn)效率的提高。

（7）財(cái)政支出水平(FE/%):該指標(biāo)用各糧食主產(chǎn)區(qū)的農(nóng)林水事務(wù)支出占其財(cái)政一般預(yù)算支出比重表示。財(cái)政支出水平越高，越有利于提高農(nóng)業(yè)生產(chǎn)的積極性，對(duì)農(nóng)業(yè)的良性快速發(fā)展的引導(dǎo)作用就越強(qiáng)。

本文基于上述情況，分析各要素對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率的影響，構(gòu)建面板模型如下：

式(3)中，GTFP 參考李谷成的研究方法，將以上一年為基準(zhǔn)的環(huán)比指數(shù)轉(zhuǎn)換為以2000 年為1 的累計(jì)變化指數(shù)，i代表糧食主產(chǎn)區(qū)各省份，t表示年份，μ表示各糧食主產(chǎn)區(qū)不可觀測(cè)的固體效應(yīng)，ε為隨機(jī)擾動(dòng)項(xiàng)。

3.2 實(shí)證結(jié)果及分析

經(jīng)過檢驗(yàn)發(fā)現(xiàn)，模型中的隨機(jī)擾動(dòng)項(xiàng)存在異方差和自相關(guān)問題，因此在穩(wěn)健標(biāo)準(zhǔn)誤回歸的基礎(chǔ)上，使用加權(quán)最小二乘法（WLS）和可行廣義最小二乘法（FGLS）對(duì)模型進(jìn)行穩(wěn)健回歸。此外，按照地理位置將糧食主產(chǎn)省份分為東中西三個(gè)區(qū)域進(jìn)行回歸，以分析影響因素對(duì)不同區(qū)域間糧食主產(chǎn)省份農(nóng)業(yè)綠色全要素生產(chǎn)效率影響差異?；貧w結(jié)果如表1 所示，各因素對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素影響估計(jì)系數(shù)顯著性與符號(hào)基本保持一致，回歸結(jié)果具有較好穩(wěn)健性。

通過表1 的回歸結(jié)果可以發(fā)現(xiàn)，（1）糧食作物受災(zāi)水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率有顯著負(fù)向影響。自然災(zāi)害頻發(fā)和防災(zāi)減災(zāi)技術(shù)落后是造成自然災(zāi)害對(duì)農(nóng)業(yè)綠色全要素生產(chǎn)率影響的重要原因。從區(qū)域來(lái)看，受災(zāi)率對(duì)不同地區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)效率的影響存在差異，具體表現(xiàn)為對(duì)東部地區(qū)的糧食主產(chǎn)省份影響并不顯著，對(duì)中西部省份農(nóng)業(yè)綠色生產(chǎn)效率影響較大，其原因在于，中西部農(nóng)業(yè)生產(chǎn)規(guī)模大，災(zāi)害防控機(jī)制尚不完善，自然災(zāi)害應(yīng)對(duì)能力較弱。

（2）灌溉設(shè)備水平對(duì)糧食主產(chǎn)區(qū)整體的農(nóng)業(yè)綠色全要素生產(chǎn)率呈現(xiàn)顯著的負(fù)向影響，其原因在于我國(guó)農(nóng)業(yè)水利設(shè)施水平較低，水平設(shè)施建設(shè)較為薄弱，導(dǎo)致供給水平低于需求水平，使得水利設(shè)備作為一種公共產(chǎn)品出現(xiàn)“過度使用”和“擁擠效應(yīng)”，阻礙了農(nóng)業(yè)綠色全要素生產(chǎn)效率的提升。從區(qū)域差異來(lái)看，灌溉設(shè)備水平對(duì)東部地區(qū)省份農(nóng)業(yè)綠色全要素生產(chǎn)效率為顯著正向影響，而對(duì)中西部地區(qū)表現(xiàn)為顯著負(fù)向影響，一方面說明我國(guó)中西部省份水利灌溉設(shè)備供給不足問題較為嚴(yán)重，另一方面也揭示了我國(guó)農(nóng)業(yè)灌溉設(shè)施建設(shè)在區(qū)域間不平衡問題，應(yīng)加快各農(nóng)業(yè)省份灌溉水利設(shè)施建設(shè)。

（3）機(jī)械化水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)效率影響顯著為正。推動(dòng)機(jī)械化水平提高是降低農(nóng)業(yè)生產(chǎn)成本，轉(zhuǎn)變農(nóng)業(yè)發(fā)展方式，促進(jìn)產(chǎn)業(yè)可持續(xù)發(fā)展的必然選擇，機(jī)械化水平的提高能夠顯著提高糧食主產(chǎn)區(qū)省份農(nóng)業(yè)綠色全要素生產(chǎn)效率。從區(qū)域差異來(lái)看，機(jī)械化對(duì)東部省份農(nóng)業(yè)綠色全要素生產(chǎn)效率的影響并不顯著，原因在于規(guī)模化是推動(dòng)機(jī)械化的基礎(chǔ)，東部省份地區(qū)地形復(fù)雜，丘陵山區(qū)較多，機(jī)械作業(yè)較為復(fù)雜，因此在推動(dòng)農(nóng)業(yè)機(jī)械化的過程中要因地制宜選擇機(jī)械化策略。

（4）農(nóng)資化肥負(fù)荷指數(shù)對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)效率呈現(xiàn)顯著負(fù)向影響。農(nóng)戶投入較高的化肥、農(nóng)藥和農(nóng)膜來(lái)提高農(nóng)作物產(chǎn)量，對(duì)提高農(nóng)戶收入發(fā)揮了顯著作用。但是化肥和農(nóng)藥會(huì)對(duì)土地造成化學(xué)污染，無(wú)法清理的殘膜改變土壤理化性狀，對(duì)農(nóng)作物生長(zhǎng)造成了負(fù)面影響，廢舊農(nóng)膜焚燒、遺棄地頭等不合理回收方式造成了農(nóng)田的白色污染，阻礙了農(nóng)業(yè)綠色全要素生產(chǎn)效率的發(fā)展。

（5）勞均經(jīng)營(yíng)規(guī)模對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)效率表現(xiàn)為顯著負(fù)向影響，表明單位勞動(dòng)力效率增加有利于提高農(nóng)業(yè)綠色全要素生產(chǎn)效率，人力資本的積累對(duì)提高糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率發(fā)揮著重要作用，因此，應(yīng)加大糧食主產(chǎn)省份的農(nóng)業(yè)教育培訓(xùn)，提高農(nóng)村醫(yī)療健康投資水平，加強(qiáng)農(nóng)民綜合素質(zhì)，提高糧食主產(chǎn)區(qū)農(nóng)業(yè)人力資本水平和單位勞動(dòng)生產(chǎn)效率。

（6）第二三產(chǎn)業(yè)發(fā)展水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)效率有顯著負(fù)向影響。糧食主產(chǎn)區(qū)的第二三產(chǎn)業(yè)發(fā)展帶來(lái)的虹吸效應(yīng)導(dǎo)致農(nóng)業(yè)生產(chǎn)要素大量流失，其技術(shù)溢出對(duì)農(nóng)業(yè)綠色全要素生產(chǎn)率發(fā)展的溢出效應(yīng)尚未顯現(xiàn)，第二三產(chǎn)業(yè)對(duì)農(nóng)業(yè)的反哺效應(yīng)并不顯著。從區(qū)域差異來(lái)看，第二三產(chǎn)業(yè)發(fā)展水平對(duì)中西部省份的農(nóng)業(yè)綠色全要素生產(chǎn)效率影響更為明顯，原因在于該地區(qū)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平較低，亟需通過促進(jìn)農(nóng)村一二三產(chǎn)融合等方式來(lái)提高當(dāng)?shù)剞r(nóng)業(yè)勞動(dòng)力水平，促進(jìn)農(nóng)業(yè)綠色全要素生產(chǎn)效率提高。

（7）財(cái)政支出水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率表現(xiàn)為顯著正向影響。農(nóng)業(yè)財(cái)政支出有利于提高農(nóng)村勞動(dòng)力生產(chǎn)積極性，促進(jìn)農(nóng)業(yè)技術(shù)創(chuàng)新和推廣，加快農(nóng)業(yè)水利基礎(chǔ)建設(shè)，完善農(nóng)村基礎(chǔ)設(shè)施，促進(jìn)農(nóng)業(yè)綠色全要素生產(chǎn)率提高。從區(qū)域差異來(lái)看，財(cái)政支出對(duì)中西部省份農(nóng)業(yè)綠色全要素生產(chǎn)效率提高效果最為顯著，而對(duì)東部影響并不顯著，說明農(nóng)業(yè)財(cái)政支出結(jié)構(gòu)有待優(yōu)化，需要探索建立完善的農(nóng)業(yè)補(bǔ)貼制度，實(shí)現(xiàn)農(nóng)業(yè)綠色全要素生產(chǎn)效率提高。

4 研究結(jié)論與政策建議

本文基于 SBM-DDF 方法的Luenberger 指數(shù)對(duì)中國(guó)糧食主產(chǎn)區(qū)2000—2017 年的農(nóng)業(yè)綠色全要素生產(chǎn)率進(jìn)行了測(cè)算，并進(jìn)一步探討其農(nóng)業(yè)全要素生產(chǎn)率收斂問題，具結(jié)論如下：

（1）從整體來(lái)看，全國(guó)糧食主產(chǎn)區(qū)總體農(nóng)業(yè)GTFP 年均值為1.029，除個(gè)別年份外，其在碳排放約束下的農(nóng)業(yè)GTFP 均值均大于1，說明在大多數(shù)年份中我國(guó)糧食主產(chǎn)區(qū)的農(nóng)業(yè)綠色全要素生產(chǎn)效率都存在一定程度的改善。

（2）各糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 水平差異明顯，其中河北、江蘇、山東、遼寧的農(nóng)業(yè)GTFP、技術(shù)進(jìn)步指數(shù)和技術(shù)效率水平均大于全國(guó)糧食主產(chǎn)區(qū)水平，處于第一方陣；安徽、河南、湖北、湖南、江西、黑龍江、吉林的農(nóng)業(yè)GTFP、技術(shù)進(jìn)步指數(shù)和技術(shù)效率水平處于第二方陣，是對(duì)第一方陣的“追逐者”和“模仿者”；四川和內(nèi)蒙古的農(nóng)業(yè)GTFP、技術(shù)進(jìn)步指數(shù)和技術(shù)效率處于第三方陣，缺乏農(nóng)業(yè)綠色經(jīng)濟(jì)增長(zhǎng)動(dòng)力。

（3）從糧食主產(chǎn)區(qū)農(nóng)業(yè)GTFP 的分解項(xiàng)上來(lái)看，碳排放約束下的農(nóng)業(yè)GTFP 增長(zhǎng)主要依靠綠色技術(shù)進(jìn)步的推動(dòng)，綠色技術(shù)進(jìn)步指數(shù)對(duì)于農(nóng)業(yè)GTFP 的作用一直領(lǐng)先于綠色技術(shù)效率，而綠色技術(shù)效率在其增速落后于綠色技術(shù)進(jìn)步的情況下，出現(xiàn)了持續(xù)惡化現(xiàn)象。

（4）自然條件、技術(shù)條件、政策條件和社會(huì)經(jīng)濟(jì)條件是糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率的重要影響因素。機(jī)械化水平、財(cái)政支出水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率表現(xiàn)為顯著正向影響；糧食作物受災(zāi)水平、灌溉設(shè)備水平、農(nóng)資化肥負(fù)荷指數(shù)、勞均經(jīng)營(yíng)規(guī)模、第二三產(chǎn)業(yè)發(fā)展水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)效率有顯著負(fù)向影響。糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素存在差異，東部地區(qū)糧食主產(chǎn)區(qū)灌溉設(shè)備水平、人力資本存量影響顯著，中西部糧食主產(chǎn)區(qū)受糧食作物受災(zāi)水平、機(jī)械化水平、第二三產(chǎn)業(yè)發(fā)展水平、財(cái)政支出影響顯著。

從上述結(jié)論可以看出，糧食主產(chǎn)區(qū)的農(nóng)業(yè)綠色發(fā)展的具體對(duì)策如下：

（1）加強(qiáng)農(nóng)業(yè)技術(shù)創(chuàng)新，要注重技術(shù)效率的提升。而農(nóng)業(yè)技術(shù)進(jìn)步既是農(nóng)業(yè)現(xiàn)代化的必然要求，也是落實(shí)鄉(xiāng)村振興戰(zhàn)略的重大舉措。為此，要積極推進(jìn)農(nóng)業(yè)生產(chǎn)機(jī)械、化肥等生產(chǎn)資料革新，降低農(nóng)業(yè)碳排放，走清潔高效的綠色農(nóng)業(yè)發(fā)展之路。

（2）推進(jìn)技術(shù)效率改善，實(shí)現(xiàn)農(nóng)業(yè)資源的優(yōu)化配置。目前13 個(gè)糧食主產(chǎn)區(qū)的農(nóng)業(yè)還處于過度依賴技術(shù)進(jìn)步對(duì)生產(chǎn)率提高的狀態(tài)，而且技術(shù)效率尚未發(fā)揮其最大的作用，故此要在推進(jìn)技術(shù)的同時(shí)，還要注重其技術(shù)效率，實(shí)現(xiàn)農(nóng)業(yè)資源的優(yōu)化配置，進(jìn)而擴(kuò)大全要素對(duì)綠色農(nóng)業(yè)發(fā)展的貢獻(xiàn)份額。

（3）優(yōu)化區(qū)域發(fā)展格局，縮小區(qū)域發(fā)展差距。從結(jié)論上可以看出，東部地區(qū)的農(nóng)業(yè)GTFP 增長(zhǎng)是農(nóng)業(yè)綠色技術(shù)進(jìn)步和技術(shù)效率“雙元”驅(qū)動(dòng)型的集約型發(fā)展，扮演著“引領(lǐng)者”和“技術(shù)革新者”的角色；中部地區(qū)農(nóng)業(yè)綠色經(jīng)濟(jì)發(fā)展主要依靠綠色技術(shù)進(jìn)步，是對(duì)東部地區(qū)農(nóng)業(yè)綠色經(jīng)濟(jì)發(fā)展的“追逐者”和“模仿者”，西部則處于技術(shù)進(jìn)步和技術(shù)效率雙項(xiàng)“塌陷”狀態(tài)，缺乏農(nóng)業(yè)綠色經(jīng)濟(jì)增長(zhǎng)動(dòng)力。因此，要要制定科學(xué)合理的區(qū)域發(fā)展政策，激發(fā)落后地區(qū)的發(fā)展?jié)摿?，進(jìn)而縮小區(qū)域發(fā)展的差距。

（4）擴(kuò)大財(cái)政支農(nóng)支出規(guī)模，優(yōu)化財(cái)政支農(nóng)的結(jié)構(gòu)。從影響因素看，財(cái)政支出水平對(duì)糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率表現(xiàn)為顯著正向影響。因此，提高財(cái)政支農(nóng)支出比重能夠有效提高糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色全要素生產(chǎn)率。糧食主產(chǎn)區(qū)應(yīng)在加大財(cái)政支農(nóng)規(guī)模的同時(shí)，優(yōu)化財(cái)政支農(nóng)結(jié)構(gòu)，并建立健全綠色農(nóng)業(yè)發(fā)展為導(dǎo)向的農(nóng)業(yè)補(bǔ)貼制度。