化肥投入變化對糧食生產(chǎn)效率的影響

孫小鈞 伍國勇 任秀

摘要：【目的】探究化肥投入變化對貴州省糧食生產(chǎn)效率的影響，為貴州省實現(xiàn)農(nóng)業(yè)綠色化、現(xiàn)代化發(fā)展和保障糧食安全提供理論參考?！痉椒ā繕?gòu)建隨機前沿生產(chǎn)函數(shù)模型，基于2001、2006、2011和2016年貴州省糧食總產(chǎn)量及其生產(chǎn)要素投入的面板數(shù)據(jù)，測算貴州省糧食生產(chǎn)效率，運用空間自相關(guān)法分析其空間分布態(tài)勢和集聚態(tài)勢，利用重力模型分析化肥投入與糧食生產(chǎn)效率的空間耦合關(guān)系，并以Tobit模型實證分析化肥投入變化對糧食生產(chǎn)效率的影響?！窘Y(jié)果】2001—2016年期間，貴州省糧食生產(chǎn)效率均值為0.8355，效率總體上雖呈下降趨勢，但提升空間較大;糧食生產(chǎn)效率空間地理分布呈西北向東南方向遞減趨勢，空間集聚態(tài)勢呈離散分布;化肥投入對糧食生產(chǎn)效率的正向影響顯著，但增產(chǎn)的邊際效用總體較低，且即將進入“內(nèi)卷化”階段?！窘ㄗh】充分考慮貴州省實際情況，采取測土配方施肥、加大農(nóng)業(yè)生產(chǎn)科研投入、區(qū)域良種和技術(shù)研發(fā)有償合作、糧食分區(qū)生產(chǎn)政策，既能防止化肥施用“內(nèi)卷化”，又能提高糧食生產(chǎn)效率。

關(guān)鍵詞： 化肥投入;糧食生產(chǎn)效率;空間自相關(guān);重心遷移軌跡

中圖分類號： S147;F326.11? ? ? ? ? ? ? ? ? 文獻標(biāo)志碼： A 文章編號：2095-1191（2019）08-1869-09

The impact of chemical fertilizer input changes on grain production efficiency—An empirical analysis based on panel data of counties in Guizhou Province

SUN Xiao-jun1， WU Guo-yong1，2*， REN Xiu1

（1School of Economics， Guizhou University， Guiyang? 550025， China; 2China Karst Regional Rural Revitalization Research Institute， Guizhou University， Guiyang? 550025， China）

Abstract：【Objective】The aim of the study was to explore the impact of chemical fertilizer application on grain production efficiency in Guizhou Province， and provide guidance for Guizhou Province to achieve greening of agriculture， modernization and food security. 【Method】This paper constructed a stochastic frontier production function model， based on the four-stage panel data of total grain output and production factors input in Guizhou Province in 2001， 2006， 2011 and 2016， measured its grain production efficiency， and analyzed its spatial distribution situation and agglomeration situation by spatial autocorrelation method. The spatial coupling relationship between fertilizer input and its grain production efficiency was analyzed by gravity model. The Tobit model was used to empirically analyze the impact of fertilizer input change on grain production efficiency. 【Result】In 2001-2016， the average grain production efficiency in Guizhou Province was 0.8355， although the overall efficiency showed a downward trend， but the room for improvement was large. In 2001-2016， the spatial distribution of grain production efficiency showed a decreasing trend from northwest to southeast， and the spatial agglomeration situation was discrete. The positive impact of fertilizer input on grain production efficiency was significant， but the marginal utility of increasing production was generally low， and it was about to enter the involution stage. 【Suggestion】Taking full account of the actual situation in Guizhou Province， adopting soil testing and formula fertilization， increasing investment in scientific research in agricultural production， compensation for regional improved varieties and technology research and development， and sub-regional production policies of grain can prevent the involution of chemical fertilizer application and improve food production efficiency.

Key words： fertilizer input change; grain production efficiency; spatial autocorrelation; gravity shift trajectory

0 引言

【研究意義】化肥作為糧食生產(chǎn)投入要素，對糧食增產(chǎn)發(fā)揮了一定的促進作用。但過多的化肥投入會導(dǎo)致土壤有機元素結(jié)構(gòu)不平衡，使土壤肥力下降，出現(xiàn)化肥對糧食增產(chǎn)的邊際產(chǎn)出遞減現(xiàn)象（陸文聰和劉聰，2017），此現(xiàn)象又被稱為化肥施用“內(nèi)卷化”。同時，過量的化肥投入會使土壤酸化，造成生態(tài)環(huán)境污染（周亮等，2014）。近年來，貴州省糧食產(chǎn)量與農(nóng)用化肥施用量總體上呈上漲趨勢，且農(nóng)用化肥施用量增長速度快于糧食產(chǎn)量增加速度。貴州省獨特的山地地形使其耕地呈現(xiàn)小塊化、分散化特征，使得機械替代化肥潛能降低;尤其是喀斯特地貌的耕地土壤有機質(zhì)流失嚴(yán)重，土壤難以存儲肥力，且過量投入化肥對其糧食生產(chǎn)效率帶來負向影響。因此，通過探討化肥投入變化對其糧食生產(chǎn)效率的影響，對貴州省糧食生產(chǎn)持續(xù)實施科學(xué)施肥指導(dǎo)意見提供依據(jù)，為實現(xiàn)其農(nóng)業(yè)綠色化提供方向指引;同時對促進貴州省糧食生產(chǎn)要素投入結(jié)構(gòu)升級、提升糧食增產(chǎn)潛能保障糧食生產(chǎn)安全及實現(xiàn)鄉(xiāng)村振興具有重要意義。【前人研究進展】目前，有關(guān)糧食生產(chǎn)效率的研究較豐富，多數(shù)學(xué)者運用相關(guān)模型測算糧食生產(chǎn)效率及其影響因素。范麗霞（2016）運用非徑向、非角度的方向性距離函數(shù)（SBM）模型和曼奎斯特（ML）指數(shù)法探討技術(shù)效率、技術(shù)進步與糧食生產(chǎn)效率的內(nèi)在邏輯關(guān)聯(lián)，指出技術(shù)進步顯著推動了我國糧食全要素生產(chǎn)率的提高，而技術(shù)效率未產(chǎn)生實質(zhì)性的影響;彭代彥和文樂（2016）運用超越對數(shù)隨機前沿生產(chǎn)函數(shù)（STF）模型探討勞動力結(jié)構(gòu)對糧食生產(chǎn)效率的影響，認為老齡化抑制了我國北方地區(qū)糧食生產(chǎn)效率，而女性化推動了北方地區(qū)糧食生產(chǎn)效率的提升;唐建和Jose Vila（2016）、楊皓天等（2016）通過隨機前沿生產(chǎn)函數(shù)（SFA）模型分析指出，加強職業(yè)農(nóng)民培育、大力發(fā)展農(nóng)業(yè)新技術(shù)是提高糧食生產(chǎn)效率的重要抓手;高鳴等（2017）基于數(shù)據(jù)包絡(luò)分析（DEA）方法，運用樣本選擇模型探討了種糧補貼對糧食生產(chǎn)效率的影響，指出農(nóng)業(yè)補貼有利于減少糧食生產(chǎn)效率損失;曾雅婷等（2018b）運用Tobit模型探討了農(nóng)地流轉(zhuǎn)對糧食生產(chǎn)效率的影響，發(fā)現(xiàn)農(nóng)地流轉(zhuǎn)推動了糧食生產(chǎn)效率的提高;卓樂和曾福生（2018）采用系統(tǒng)廣義矩估計法（SYS-GMM）驗證了農(nóng)村基礎(chǔ)設(shè)施對糧食全要素生產(chǎn)效率的影響，認為水利基礎(chǔ)設(shè)施與交通基礎(chǔ)設(shè)施對糧食生產(chǎn)效率有正向推動作用。關(guān)于化肥對糧食生產(chǎn)影響的相關(guān)研究結(jié)果表明，化肥投入的增加對糧食生產(chǎn)起正向影響，以王祖力和肖海峰（2008）、王美兔和田明華（2016）、衛(wèi)龍寶等（2017）為主要代表，認為農(nóng)用化肥投入的增加有利于提高糧食生產(chǎn)能力，從而增加糧食產(chǎn)出。然而，化肥施用量的大幅增加會對糧食生產(chǎn)起負向影響，李想等（2017）進一步分析發(fā)現(xiàn)，化肥施用量有閾值限制，當(dāng)化肥施用量低于該閾值時，糧食產(chǎn)量將會保持較快增長;反之，糧食產(chǎn)量將不再增長，甚至出現(xiàn)減產(chǎn)。當(dāng)前，我國農(nóng)業(yè)生產(chǎn)中化肥施用量的年均增速明顯快于糧食產(chǎn)量增速，使得化肥實際施用量已極大地超出了科學(xué)施肥量（史常亮和朱俊峰，2016），導(dǎo)致糧食生產(chǎn)中化肥的正增產(chǎn)效應(yīng)正在逐漸降低，從而出現(xiàn)單位糧食播種面積內(nèi)化肥投入?yún)s帶來糧食實際產(chǎn)量不斷減少的現(xiàn)象（趙志堅和胡小娟，2013;張利庠等，2018）。因此，未來合理管控化肥投入是優(yōu)化糧食生產(chǎn)要素投入結(jié)構(gòu)的重要發(fā)力點。【本研究切入點】至今，鮮有文獻從省域?qū)用嫣接懟释度雽Z食生產(chǎn)效率的影響，相關(guān)文獻也極少引入空間地理維度探討糧食生產(chǎn)效率的空間集聚分布態(tài)勢?！緮M解決的關(guān)鍵問題】基于2001、2006、2011和2016年4個年份的貴州省糧食總產(chǎn)量與要素投入的面板數(shù)據(jù)，測算貴州省糧食生產(chǎn)效率，分析糧食生產(chǎn)效率的空間分布態(tài)勢和集聚態(tài)勢，探討化肥投入變化對貴州省糧食生產(chǎn)效率的影響，為優(yōu)化糧食生產(chǎn)要素投入結(jié)構(gòu)和提高糧食生產(chǎn)效率提供指導(dǎo)，同時為貴州省實現(xiàn)農(nóng)業(yè)綠色化、現(xiàn)代化和保障其糧食安全提供理論參考。

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

1. 1 隨機前沿生產(chǎn)函數(shù)（SFA）模型

結(jié)合貴州省實際情況，構(gòu)建測度糧食生產(chǎn)效率的SFA如下：

ln（Yit）=γ0+γ1ln（Lit）+γ2ln（Mit）+γ3ln（Pit）+γ4ln（Fit）+

γ5ln（Iit）+γ6ln（Ait）+νit-μit? ? ? ? ?（1）

式中，i為縣（區(qū)），t為年份，Y為糧食產(chǎn)量;L為勞動力;M為農(nóng)機總動力;P為農(nóng)藥施用量;F為化肥施用量;I為收入水平，采用農(nóng)戶家庭人均純收入替代;A為糧食播種面積;γ0表示截距項，γi（i=1，2，3，…，6）為待估系數(shù);ν為擾動項，用以測度系統(tǒng)非效率，并服從正態(tài)分布;μ為技術(shù)效率損失項，服從半正態(tài)分布。上述各變量描述性統(tǒng)計詳見表1。

1. 2 空間自相關(guān)法、重力模型和Tobit模型

利用全局和局部Moran’s I指數(shù)測度貴州省糧食生產(chǎn)效率的時空演化趨勢和空間集聚態(tài)勢，計算公式如下：

I=[i=1nj≠inωij（xi-x）（xj-x）S2i=1nj≠inωij]? （2）

Ii=[（xi-x）S2][j=inωij（xj-x）]? （3）

式中，i、j為研究對象;ω為空間權(quán)重;x為糧食生產(chǎn)效率值;`x為糧食生產(chǎn)效率均值;S2為糧食生產(chǎn)效率值方差。其中，I∈[-1，1]。當(dāng)I>0時，說明貴州省糧食生產(chǎn)效率為正空間相關(guān)性，此時以高—高（High-High）、低—低（Low-Low）集聚態(tài)勢分布為主;當(dāng)I<0時，說明貴州省糧食生產(chǎn)效率為負空間相關(guān)性，此時以低—高（Low-High）、高—低（High-Low）集聚態(tài)勢分布為主;當(dāng)I=0時，說明貴州省糧食生產(chǎn)效率并不存在空間相關(guān)性，其空間集聚態(tài)勢呈隨機分布。

運用重力模型（周亮等，2014）分析貴州省糧食生產(chǎn)效率與化肥投入變化的空間耦合關(guān)系及作用機理，計算公式如下：

X=[i=1nDixii=1nDi]，Y=[i=1nDiyii=1nDi]? （4）

式中，（xi，yi）表示次級地區(qū)i的中心經(jīng)緯坐標(biāo)，Di表示次級地區(qū)i的觀測值，（X，Y）為觀測值的重心經(jīng)緯坐標(biāo)。

另外，參考楊勇等（2017）的方法構(gòu)建Tobit回歸模型，分析化肥投入變化對貴州省糧食生產(chǎn)效率的影響。

1. 3 數(shù)據(jù)來源

本研究數(shù)據(jù)均來源于《中國統(tǒng)計年鑒》《中國農(nóng)村統(tǒng)計年鑒》《貴州統(tǒng)計年鑒》《新中國六十年統(tǒng)計資料匯編》中相應(yīng)年份數(shù)據(jù)。在數(shù)據(jù)收集過程中，由于部分縣（區(qū)）要素投入數(shù)據(jù)空缺，因此采用各縣（區(qū)）糧食播種面積占農(nóng)作物播種面積的比值與要素投入的乘積獲得相應(yīng)數(shù)據(jù)。

1. 4 統(tǒng)計分析

運用Frontier 4.1測算貴州省糧食生產(chǎn)效率;運用Arcgis 10.2繪制貴州省糧食生產(chǎn)效率空間分布、集聚態(tài)勢圖及重心遷移軌跡;以GeoDa分析貴州省糧食生產(chǎn)效率空間相關(guān)性，并采用Stata v15.0分析化肥投入變化對貴州省糧食生產(chǎn)效率存在的影響。

2 結(jié)果與分析

2. 1 要素投入對貴州省糧食生產(chǎn)的影響

如表2所示，除農(nóng)藥施用量對貴州省糧食生產(chǎn)呈負向影響外，其余要素均對貴州省糧食生產(chǎn)有正向影響。其中，糧食播種面積對貴州省糧食生產(chǎn)的正向影響最大，表明因地制宜提高糧食播種面積比例，發(fā)展適度規(guī)模經(jīng)營，有利于充分釋放糧食增產(chǎn)潛能。農(nóng)業(yè)機械總動力影響其次，表明增加農(nóng)機投入有利于提高糧食生產(chǎn)效率。從事農(nóng)業(yè)生產(chǎn)的勞動力對糧食生產(chǎn)的促進作用低于農(nóng)業(yè)機械總動力，其主要原因是糧食生產(chǎn)中勞動力投入趨于飽和，邊際產(chǎn)出下降;同時，政府對農(nóng)機的大力推廣，替代了一部分農(nóng)業(yè)勞動力。農(nóng)藥施用量對貴州省糧食生產(chǎn)有負向影響，表明在糧食生產(chǎn)中存在農(nóng)藥施用“內(nèi)卷化”現(xiàn)象。因此，實施農(nóng)藥減量投入政策既能防止過量施用農(nóng)藥對糧食生產(chǎn)帶來的負向影響，又能減少對生態(tài)環(huán)境的污染。化肥施用量對貴州省糧食生產(chǎn)有正向影響，但促進作用很弱，說明化肥投入也趨于飽和，導(dǎo)致化肥增產(chǎn)的邊際效應(yīng)降低。通過實施化肥減量投入政策，防止過量施用化肥帶來的負向影響，能促進糧食生產(chǎn)。農(nóng)戶家庭收入水平對貴州省糧食生產(chǎn)的正向促進作用也很弱，究其原因在于，糧食作物預(yù)期收益較低，出于理性決策，農(nóng)戶在保持基本口糧安全的前提下，會把更多的收入投入到經(jīng)濟作物生產(chǎn)中，使得糧食生產(chǎn)增長乏力。

2. 2 貴州省糧食生產(chǎn)效率時空演化趨勢

基于SFA估計結(jié)果，進一步描繪2001、2006、2011和2016年貴州省糧食生產(chǎn)效率空間分布態(tài)勢?？傮w而言，4個年份貴州省糧食生產(chǎn)效率平均值分別為0.8398、0.8370、0.8341和0.8312，其提升空間較大。圖1顯示，貴州省糧食生產(chǎn)效率的空間地理分布總體上呈由西北向東南方向遞減的趨勢，且空間變化趨勢平穩(wěn)。從局部來看，糧食生產(chǎn)效率在0.8000以上的高值區(qū)主要分布在貴州省東北部、北部、西北部、西部及西南部，約占全省縣域總數(shù)的49.38%，具體包括銅仁市（松桃縣、銅仁市、印江縣、沿河縣、德江縣和思南縣）、遵義市、黔南州（甕安縣、福泉市、獨山縣、惠水縣）、畢節(jié)市、六盤水市（水城縣、盤縣）、黔西南州（興義縣、興仁縣、安龍縣）、安順市（普定縣）及黔東南州（天柱縣、黎平縣、從江縣）;糧食生產(chǎn)效率在0.7100～0.8000的較高值區(qū)分布在貴州省的東北部和南部地區(qū)，約占全省縣域總數(shù)的25.93%;糧食生產(chǎn)效率值在0.6100～0.7000的低值區(qū)分布在貴州省的中部、中東部、中南部及南部，約占全省縣域總數(shù)的22.22%;糧食生產(chǎn)效率在0.6000以下的較低值區(qū)分布在貴州省的東北部、西部及南部，約占全省縣域總數(shù)的3.70%，具體包括銅仁市的萬山特區(qū)、黔西南州的望謨縣以及六盤水市的鐘山區(qū)。上述結(jié)果表明，貴州省糧食生產(chǎn)效率空間分布特征與其地形地貌分布特征總體一致，反映出地形特征極大制約了糧食生產(chǎn)效率的提高。

2. 3 貴州省糧食生產(chǎn)效率空間自相關(guān)集聚特征

運用空間自相關(guān)法分析糧食生產(chǎn)效率空間集聚特征，通過全局Moran’s I指數(shù)分析得出，2001、2006、2011和2016年貴州省Morain’s I指數(shù)均值為 -0.0091，表明貴州省糧食生產(chǎn)效率呈離散分布特征，空間自相關(guān)為負。而根據(jù)局部Moran’s I指數(shù)分析發(fā)現(xiàn)，2001、2006、2011和2016年貴州省糧食生產(chǎn)效率主要呈High-Low、Low-High集聚態(tài)勢分布（圖2）。隨著時間推移，貴州省糧食生產(chǎn)效率空間變化總體平穩(wěn)。

2. 4 化肥投入變化對貴州省糧食生產(chǎn)效率的影響分析

2. 4. 1 化肥投入與糧食生產(chǎn)效率的重心演化分析 將化肥投入指標(biāo)分解為化肥施用量、化肥施用增長率、單位化肥施用量、單位化肥施用增長率和化肥施用結(jié)構(gòu)比（氮、磷、鉀肥配施比例）。上述指標(biāo)描述及統(tǒng)計值詳見表3。

運用重力模型測算化肥投入與糧食生產(chǎn)效率的重心遷移軌跡，通過分析二者間的空間耦合關(guān)系，有利于從理論上進一步揭示其內(nèi)在邏輯關(guān)聯(lián)。結(jié)果表明，貴州省糧食生產(chǎn)效率重心由東南向西北方向遷移（圖3-a）;化肥施用量重心遷移軌跡在2001—2006年由東北向西南方向移動，2007—2016年由西南向東北方向移動（圖3-b）;化肥施用量增長率重心遷移軌跡在2001—2006年向西北方向移動，而2006年之后其軌跡由東北向西南方向移動（圖3-c）;單位化肥用量重心遷移軌跡在2001—2006年由東北向西南方向移動，而2006年之后由西南向東北方向移動（圖3-d）;單位化肥用量增長率重心遷移軌跡在2001—2006年向西北方向移動，而2006年之后其軌跡由東北向西南方向移動（圖3-e）;化肥施用結(jié)構(gòu)比重心遷移軌跡在2001—2006年向西北方向移動，2007—2011年由東北向西南移動，而2012—2016年由西南向東北方向移動（圖3-f）。

綜上所述，在2001—2016年化肥施用量增長率、單位化肥用量增長率與貴州省糧食生產(chǎn)效率存在相關(guān)性;在2001—2006年化肥施用量、單位化肥用量與糧食生產(chǎn)效率存在相關(guān)性;在2007—2011年化肥施用結(jié)構(gòu)比與糧食生產(chǎn)效率存在相關(guān)性。

2. 4. 2 化肥投入變化對糧食生產(chǎn)效率影響的實證分析 基于Tobit模型，選取化肥投入變量指標(biāo)分析化肥投入變化對其糧食生產(chǎn)效率的影響，模型設(shè)計如下：

GPEi=α0+λ1CFi+λ2GRFi+λ3CFPi+λ4GRPFi+

λ5SRFi+εi? （5）

式中，i為縣（區(qū)）;GRP為糧食生產(chǎn)效率;CF為化肥施用總量;GRF為化肥施用增長率;CFP為單位面積化肥施用總量;GRPF為單位化肥施用增長率;SRF為化肥施用結(jié)構(gòu)比;α為常數(shù)項;λ為待估系數(shù);ε為擾動項。

為反映各變量隨時間變化對糧食生產(chǎn)效率的影響，選擇2001、2006、2011和2016年4期截面數(shù)據(jù)進行Tobit回歸分析，模型Ⅰ～Ⅳ分別表示上述4期截面數(shù)據(jù)回歸分析結(jié)果（表4）。Tobit回歸分析結(jié)果表明：（1）化肥施用總量（CF）在模型Ⅰ～Ⅳ中對糧食生產(chǎn)效率的正向影響均在1%水平上顯著，但化肥增產(chǎn)的邊際效用總體上較低。（2）化肥施用增長率（GRF）對糧食生產(chǎn)效率的影響總體在5%水平上顯著（除模型Ⅲ外）。其中，在模型Ⅰ和模型Ⅱ中對糧食生產(chǎn)效率有負向影響，在模型Ⅲ和模型Ⅳ中對糧食生產(chǎn)效率有正向影響。其主要原因是2001—2006年貴州省經(jīng)濟社會發(fā)展水平較低，農(nóng)業(yè)生產(chǎn)資金投入不足，使得農(nóng)業(yè)生產(chǎn)機械化程度較低，一定程度上促使農(nóng)戶增加化肥投入以保障糧食產(chǎn)量。但在2006年之后，隨著全省經(jīng)濟社會的快速發(fā)展，農(nóng)業(yè)生產(chǎn)獲得大量資金投入，農(nóng)業(yè)機械和良種研發(fā)并推廣，一些傳統(tǒng)投入要素開始被農(nóng)機和良種替代，在一定程度上也推動了化肥減量投入。（3）單位化肥施用量（CFP）和化肥施用結(jié)構(gòu)比（SRF）對糧食生產(chǎn)效率的負向影響差異顯著，但隨時間推移影響總體上降低。一方面，貴州省糧食生產(chǎn)中存在不合理施肥行為，導(dǎo)致化肥施用存在過量現(xiàn)象;另一方面，隨著農(nóng)機和良種等要素投入替代與測土配方施肥政策的實施，減少了不合理施肥行為帶來的負向影響。（4）單位化肥施用增長率（GRPF）對糧食生產(chǎn)效率的影響差異顯著，但其影響逐漸降低。其中，在模型Ⅰ和模型Ⅱ中正向影響逐漸降低，表明單位化肥投入使得其正增產(chǎn)邊際效應(yīng)降低，此時存在施肥過量行為;在模型Ⅲ和模型Ⅳ中負向影響逐漸降低，則表明隨著農(nóng)機和良種等要素投入替代與測土配方施肥政策的實施，減少了不合理施肥行為帶來的負向影響。

3 討論

本研究構(gòu)建了SFA模型，并基于2001—2016年中4個年份貴州省糧食總產(chǎn)量及其生產(chǎn)要素投入的面板數(shù)據(jù)測算貴州省糧食生產(chǎn)效率，結(jié)果表明，2001—2016年貴州省糧食生產(chǎn)效率均值為0.8355，雖然貴州省糧食生產(chǎn)效率總體上呈下降趨勢，但提升空間較大。該結(jié)論與閔銳和李谷成（2013）、田旭和王善高（2016）、曾雅婷等（2018a）在分析我國糧食生產(chǎn)效率過程中測算出的貴州省糧食生產(chǎn)效率基本一致。此外，SFA模型估計結(jié)果進一步證實了唐建和Jose Vila（2016）關(guān)于糧食播種面積的變化顯著影響糧食生產(chǎn)效率的結(jié)論。

本研究運用空間自相關(guān)法分析貴州省糧食生產(chǎn)效率的空間分布態(tài)勢和集聚態(tài)勢，并采用重力模型分析化肥投入與其糧食生產(chǎn)效率的空間耦合關(guān)系，結(jié)果顯示，貴州省糧食生產(chǎn)效率空間地理分布呈由西北向東南方向遞減趨勢，空間集聚態(tài)勢呈離散分布;在研究時段內(nèi)化肥投入與其糧食生產(chǎn)效率保持相關(guān)性。同時，本研究運用Tobit模型實證分析化肥投入變化對其糧食生產(chǎn)效率的影響，結(jié)果表明，貴州省糧食生產(chǎn)中化肥的正增產(chǎn)效應(yīng)總體上較低，即將進入“內(nèi)卷化”階段，與趙志堅和胡小娟（2013）、李想等（2017）、張利庠等（2018）的研究結(jié)果基本一致，但與王祖力和肖海峰（2008）、王美兔和田明華（2016）的研究結(jié)論相反，也與化肥施用量與糧食生產(chǎn)效率通常呈現(xiàn)的倒“U”形關(guān)系有所不同?？傮w來看，主要原因是研究對象選取的不同，現(xiàn)有文獻大多基于全國層面數(shù)據(jù)進行分析，忽視了各省份間的具體差異，而本研究以省域為研究對象，選取縣域?qū)用鏀?shù)據(jù)進行分析。具體來看，其主要原因有：（1）相對于中東部發(fā)達地區(qū)，貴州省經(jīng)濟發(fā)展水平滯后，農(nóng)戶人均收入較低，難以購買農(nóng)用機械，而化肥相對便宜。同時，耕地呈現(xiàn)小塊化、分散化特征，使得糧食生產(chǎn)難以實現(xiàn)機械化、規(guī)模化經(jīng)營，極大限制了機械替代化肥的潛能釋放。農(nóng)戶為保證糧食正常生產(chǎn)，會增加化肥施用，但尚未極大超出最優(yōu)施肥量，使得化肥對本省糧食生產(chǎn)仍起正向促進作用。（2）鑒于研究時段內(nèi)貴州省糧食生產(chǎn)概況，本研究重點考察化肥投入與其糧食生產(chǎn)效率的線性關(guān)系，而兩者間的非線性關(guān)系未納入研究范疇。因此，本研究結(jié)論與通常呈現(xiàn)的倒“U”形關(guān)系不同。這也是本研究存在的不足，后續(xù)討論中可適當(dāng)引入平方項、立方項等考察化肥投入與糧食生產(chǎn)效率間的非線性關(guān)系。

此外，由于缺乏長期穩(wěn)定的調(diào)研數(shù)據(jù)，尚未引入農(nóng)戶化肥施用次數(shù)變量，在后續(xù)研究中可根據(jù)實際情況進行補充;本研究僅探討了化肥對貴州省糧食生產(chǎn)效率的影響，而糧食生產(chǎn)效率受勞動力資源稟賦、農(nóng)藥施用量、農(nóng)業(yè)基礎(chǔ)設(shè)施、氣候變化、種糧補貼、農(nóng)地流轉(zhuǎn)、城鎮(zhèn)化及耕地質(zhì)量等多種因素共同制約，因此，在后續(xù)研究中可適當(dāng)引入上述變量構(gòu)建一個綜合分析框架進行探討。

4 建議

4. 1 堅定實施測土配方施肥措施

繼續(xù)建立和完善適合貴州省的農(nóng)業(yè)土壤環(huán)境檢測體系，定期對土壤進行質(zhì)量評估，尤其加強對黔東北、黔北、黔西北和黔西等高糧食生產(chǎn)效率地區(qū)的土壤環(huán)境監(jiān)測與質(zhì)量評估，制訂一套有針對性的農(nóng)用化肥施用方案。同時，加大貴州省實施測土配方施肥補貼政策力度，積極鼓勵農(nóng)戶使用農(nóng)業(yè)有機肥逐步替代農(nóng)用化肥，以減緩化肥施用強度，防止化肥施用“內(nèi)卷化”。

4. 2 加大農(nóng)業(yè)生產(chǎn)科研投入

為促進貴州省糧食生產(chǎn)效率的提高，應(yīng)大力實施科技興農(nóng)、藏技于糧戰(zhàn)略。第一，繼續(xù)孵化和培育本土農(nóng)機企業(yè)，在財政上給予此類企業(yè)較大的扶持力度，以增強貴州省農(nóng)機企業(yè)“自主創(chuàng)新、自主研發(fā)、自主生產(chǎn)”一體化能力。第二，根據(jù)貴州省地形地貌特征，研究開發(fā)適合各類地形農(nóng)業(yè)生產(chǎn)的小型多功能農(nóng)機，構(gòu)建完善的農(nóng)機售后服務(wù)體系，進而提高農(nóng)業(yè)生產(chǎn)的農(nóng)機配套綜合能力。第三，依托省內(nèi)高校、農(nóng)業(yè)科學(xué)院等核心機構(gòu)開展優(yōu)質(zhì)玉米和冬馬鈴薯培育，并加大對超級稻的研究和開發(fā)投入。

4. 3 實施區(qū)域良種和技術(shù)研發(fā)有償合作

建立橫向區(qū)域農(nóng)業(yè)生產(chǎn)合作補償機制，科學(xué)引導(dǎo)黔中、黔東和黔南等地區(qū)的低糧食生產(chǎn)效率市縣與黔東北、黔北、黔西北和黔西等地區(qū)的高糧食生產(chǎn)效率市縣開展農(nóng)業(yè)生產(chǎn)技術(shù)有償合作，適當(dāng)引入優(yōu)質(zhì)、高產(chǎn)的糧食作物新品種與先進的農(nóng)作物種苗培育、作物栽種技術(shù)和經(jīng)營管理理念。推動貴州各地區(qū)糧食生產(chǎn)投入要素合理流動，實現(xiàn)糧食生產(chǎn)資源的科學(xué)配置，保進糧食生產(chǎn)低效率地區(qū)糧食生產(chǎn)效率的提高，共同推動貴州省糧食生產(chǎn)效率的提升，從而確保貴州省糧食生產(chǎn)安全。

4. 4 實施糧食分區(qū)生產(chǎn)

穩(wěn)定黔東北、黔北、黔西北和黔西等地區(qū)的高糧食生產(chǎn)效率市縣的糧食生產(chǎn)，加強該地區(qū)高標(biāo)準(zhǔn)農(nóng)田建設(shè)和現(xiàn)代化農(nóng)業(yè)生產(chǎn)體系建設(shè)，在全省范圍內(nèi)形成示范效應(yīng)。重點加強黔東北、黔南等地區(qū)的較高糧食生產(chǎn)效率市縣與黔中、黔中東和黔中南等地區(qū)的低糧食生產(chǎn)效率市縣的農(nóng)業(yè)生產(chǎn)投入，有步驟地引導(dǎo)該地區(qū)農(nóng)地流轉(zhuǎn)，適當(dāng)擴大糧食播種比例，實現(xiàn)機械化和規(guī)?；a(chǎn)。扶持黔東北、黔西、黔南等地區(qū)的較低糧食生產(chǎn)效率市縣的糧食生產(chǎn)，保障人均基本口糧安全。同時，對地處城郊邊緣的糧食生產(chǎn)區(qū)可適度引導(dǎo)其三產(chǎn)融合發(fā)展，推動產(chǎn)業(yè)結(jié)構(gòu)升級，發(fā)展現(xiàn)代農(nóng)業(yè)。

參考文獻：

范麗霞. 2016. 技術(shù)效率、技術(shù)進步與糧食生產(chǎn)率增長[J]. 經(jīng)濟經(jīng)緯，33（3）：31-36. [Fan L X. 2016. A study on growth and decomposition of China’s grain TFP growth[J]. Economic Survey，33（3）： 31-36.]

高鳴，宋洪遠，Michael Carter . 2017. 補貼減少了糧食生產(chǎn)效率損失嗎？——基于動態(tài)資產(chǎn)貧困理論的分析[J]. 管理世界，（9）：85-100. [Gao M，Song H Y，Carter M. 2017. Does subsidies reduce the loss of food production efficiency？ —Based on the analysis of dynamic asset poverty theory[J]. Management World，（9）： 85-100.]

李想，戴維，高紅菊，徐文平，魏小紅. 2017. 基于BP神經(jīng)網(wǎng)絡(luò)的糧食產(chǎn)量與化肥用量相關(guān)性研究[J]. 農(nóng)業(yè)機械學(xué)報，48（S）：186-192. [Li X，Dai W，Gao H J，Xu W P，Wei X H. 2017. Correlation between grain yield and fertilizer use based on BP neural network[J]. Transactions of The Chinese Society of Agricultural Machinery，48（S）： 186-192.]

陸文聰，劉聰. 2017. 化肥污染對糧食作物生產(chǎn)的環(huán)境懲罰效應(yīng)[J]. 中國環(huán)境科學(xué)，37（5）：1988-1994. [Lu W C，Liu C. 2017. The “environmental punishment” effect of ferti-lizer pollution in grain crop production[J]. China Environmental Science，37（5）： 1988-1994.]

閔銳，李谷成. 2013. “兩型”視角下我國糧食生產(chǎn)技術(shù)效率的空間分異[J]. 經(jīng)濟地理，33（3）：144-149. [Min R，Li G C. 2013. Grain’s production efficiency and its spatial distribution in China from the “Two-oriented（Resource and Environment）” perspective[J]. Economic Geography，33（3）： 144-149.]

彭代彥，文樂. 2016. 農(nóng)村勞動力老齡化、女性化降低了糧食生產(chǎn)效率嗎——基于隨機前沿的南北方比較分析[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟，（2）：32-44. [Peng D Y，Wen L. 2016. Is the rural labor force aging and feminization reducing the efficiency of grain production？—A comparative analysis of the north and south based on stochastic frontiers[J]. Journal of Agrotechnical Economics，（2）： 32-44.]

史常亮，朱俊峰. 2016. 我國糧食生產(chǎn)中化肥投入的經(jīng)濟評價和分析[J]. 干旱區(qū)資源與環(huán)境，30（9）：57-63. [Shi C L，Zhu J F. 2016. Economic evaluation and analysis of fertilizer input in China’s grain production[J]. Journal of Arid Land Resources and Environment，30（9）： 57-63.]

唐建，Jose Vila. 2016. 糧食生產(chǎn)技術(shù)效率及影響因素研究——來自1990—2013年中國31個省份面板數(shù)據(jù)[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟，（9）：72-83. [Tang J，Vila J. 2016. Study on technical efficiency and influencing factors of grain production—From panel data of 31 provinces in China from 1990 to 2013[J]. Journal of Agrotechnical Economics，（9）： 72-83.]

田旭，王善高. 2016. 中國糧食生產(chǎn)環(huán)境效率及其影響因素分析[J]. 資源科學(xué)，38（11）：2106-2116. [Tian X，Wang S G. 2016. Environmental efficiency and its determinants regarding China’s grain production[J]. Resources Science，38（11）： 2106-2116.]

王美兔，田明華. 2016. 基于VECM模型的中國農(nóng)用化肥施用強度與糧食貿(mào)易等因素的相互關(guān)系研究[J]. 生態(tài)經(jīng)濟，32（11）：98-102. [Wang M T，Tian M H. 2016. The relationship between the intensity of chemical fertilizer and cereal trade and other factors based on VECM model[J]. Ecological Economy，32（11）： 98-102.]

王祖力，肖海峰. 2008. 化肥施用對糧食產(chǎn)量增長的作用分析[J]. 農(nóng)業(yè)經(jīng)濟問題，（8）：65-68. [Wang Z L，Xiao H F. 2008. Analysis on the role of fertilization on grain yield growth[J]. Issues in Agricultural Economy，（8）： 65-68.]

衛(wèi)龍寶，張艷紅，高敘文. 2017. 我國農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食安全的影響——基于面板數(shù)據(jù)的實證分析[J]. 經(jīng)濟問題探索，（2）：160-167. [Wei L B，Zhang Y H，Gao X W. 2017. The impact of agricultural labor transfer on food security in China—An empirical analysis based on panel data[J]. Inquiry into Economic Issues，（2）： 160-167.]

楊皓天，劉秀梅，句芳. 2016. 糧食生產(chǎn)效率的隨機前沿函數(shù)分析——基于內(nèi)蒙古微觀農(nóng)戶層面1312戶調(diào)研數(shù)據(jù)[J]. 干旱區(qū)資源與環(huán)境，30（12）：82-88. [Yang H T，Liu X M，Ju F. 2016. Stochastic frontier function analysis on the grain production efficiency—Based on the research data of 1312 farmers in Inner Mongolia[J]. Journal of Arid Land Resources and Environment，30（12）： 82-88.]

楊勇，鄧祥征，李志慧，吳鋒，李小云. 2017. 2000—2015年華北平原土地利用變化對糧食生產(chǎn)效率的影響[J]. 地理研究，36（11）：2171-2183. [Yang Y，Deng X Z，Li Z H，Wu F，Li X Y. 2017. Impact of land use change on grain production efficiency in North China Plain during 2000-2015[J]. Geographical Research，36（11）： 2171-2183.]

曾雅婷，李賓，呂亞榮. 2018a. 中國糧食生產(chǎn)技術(shù)效率區(qū)域差異及其影響因素——基于超越對數(shù)形式隨機前沿生產(chǎn)函數(shù)的測度[J]. 湖南農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），19（6）：13-21. [Zeng Y T，Li B，Lü Y R. 2018a. Regional di-fferences of grain production technical efficiency and its influencing factors in China： A method of stochastic frontier analysis on trans-logarithmic production function[J]. Journal of Hunan Agricultural University （Social Scien-ces） ，19（6）： 13-21.]

曾雅婷，呂亞榮，劉文勇. 2018b. 農(nóng)地流轉(zhuǎn)提升了糧食生產(chǎn)技術(shù)效率嗎——來自農(nóng)戶的視角[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟，（3）：41-55. [Zeng Y T，Lü Y R，Liu W Y. 2018b. Does farmland transfer promote the technical efficiency of grains production—From the farmer’s perspective[J]. Journal of Agrotechnical Economics，（3）： 41-55.]

張利庠，彭輝，靳興初. 2018. 不同階段化肥施用量對我國糧食產(chǎn)量的影響分析——基于1952—2006年30個省份的面板數(shù)據(jù)[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟，（4）：85-94. [Zhang L X，Peng H，Jin X C. 2018. Analysis on the impact of agricultural fertilizer application on China’s grain output at different stages-based on panel data for 30 provinces in 1952-2006[J]. Journal of Agrotechnical Economics，（4）： 85-94.]

趙志堅，胡小娟. 2013. 我國糧食生產(chǎn)中化肥投入的環(huán)境成本研究[J]. 湖南大學(xué)學(xué)報（社會科學(xué)版），27（6）：52-56. [Zhao Z J，Hu X J. 2013. Research on environmental cost of fertilizer inputs in Chinese grain production[J]. Journal of Hunan University（Social Sciences），27（6）： 52-56.]

周亮，徐建剛，蔡北溟，蔣金亮，孫東琪. 2014. 淮河流域糧食生產(chǎn)與化肥消費時空變化及對水環(huán)境影響[J]. 自然資源學(xué)報，29（6）：1053-1064. [Zhou L，Xu J G，Cai B M，Jiang J L，Sun D Q. 2014. The Spatio-temporal changes of grain production and fertilizer consumption and its impact on water environment in the Huaihe River Basin[J]. Journal of Natural Resources，29（6）： 1053-1064.]

卓樂，曾福生. 2018. 農(nóng)村基礎(chǔ)設(shè)施對糧食全要素生產(chǎn)率的影響[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟，（11）：92-101. [Zhuo L，Zeng F S. 2018. Research on the impact of rural infrastructure on total factor productivity of grain[J]. Journal of Agrotechnical Economics，（11）： 92-101.]

（責(zé)任編輯 鄧慧靈）