伍國勇 包瀅暉

摘 要：本文測度城鎮(zhèn)化對化肥施用強度的影響，為降低化肥施用強度和實現(xiàn)糧食生產(chǎn)綠色化、現(xiàn)代化提供指引。本文基于2001—2017年全國31省面板數(shù)據(jù)，采用熵值法對構(gòu)建的城鎮(zhèn)化綜合指標進行賦權(quán)。構(gòu)建面板VAR模型和面板Tobit模型分別對城鎮(zhèn)化與化肥施用強度進行脈沖響應與回歸分析。結(jié)果表明：（1）人口城鎮(zhèn)化對化肥施用強度影響顯著。其中，在全國層面、主產(chǎn)區(qū)、平衡區(qū)、中部和西部地區(qū)有正向影響，在主銷區(qū)和東部地區(qū)有負向影響。（2）土地城鎮(zhèn)化對化肥施用強度影響差異顯著。其中，在全國層面和西部地區(qū)影響在統(tǒng)計學意義上并不成立;在主銷區(qū)、平衡區(qū)和東部地區(qū)有正向影響;在主產(chǎn)區(qū)耕地面積和人均綠地面積對其有正向影響;在中部地區(qū)耕地面積和城市征用土地面積對其有負向影響。（3）經(jīng)濟城鎮(zhèn)化對化肥施用強度影響差異顯著。其中，在全國層面、主產(chǎn)區(qū)和中部地區(qū)人均地區(qū)生產(chǎn)總值與二三產(chǎn)業(yè)值比重對其有正向影響，人均非農(nóng)產(chǎn)值對其有負向影響;在主銷區(qū)、平衡區(qū)、東西部地區(qū)有正向影響?；趯嵶C分析結(jié)果，提出實現(xiàn)糧食施肥減量投入和糧食生產(chǎn)綠色化與農(nóng)業(yè)現(xiàn)代化的建議。

關(guān)鍵詞：化肥施用強度;城鎮(zhèn)化;糧食生產(chǎn);脈沖響應

中圖分類號：F326.11;S147.2

文獻標識碼：A

文章編號：1000-5099（2022）02-0049-12

化肥作為糧食生產(chǎn)要素投入的重要組成部分，對糧食增產(chǎn)發(fā)揮了至關(guān)重要的作用，但隨著化肥施用強度的增大[1]，使得糧食生產(chǎn)施肥過量[2-4]，對糧食生產(chǎn)和農(nóng)業(yè)生態(tài)環(huán)境造成負面影響。一方面，會破壞土壤有機元素平衡結(jié)構(gòu)，導致土壤存肥能力下降，使得化肥增產(chǎn)邊際效用下降[5]。另一方面，導致土壤酸化與水質(zhì)下降，嚴重污染了生態(tài)環(huán)境[6]。隨著城市化的快速推進，大量農(nóng)業(yè)勞動力非農(nóng)化成為趨勢，城市用地規(guī)模擴大引起的耕地面積減少，經(jīng)濟發(fā)展帶來的收入增加與工業(yè)反哺農(nóng)業(yè)等都影響了糧食生產(chǎn)的化肥投入。那么，城鎮(zhèn)化快速發(fā)展對化肥施用強度的具體影響是提高了還是降低了？在糧食功能區(qū)與地理分區(qū)的影響是否存在差異？因此，深入探討二者關(guān)系既回答了上述問題，也有利于優(yōu)化糧食生產(chǎn)要素投入結(jié)構(gòu)，為實現(xiàn)化肥減量投入，達到糧食生產(chǎn)綠色化與現(xiàn)代化相統(tǒng)一有重要意義。

當前，學者們對城鎮(zhèn)化與化肥施用強度的影響關(guān)系研究主要體現(xiàn)在如下層面：一是勞動力結(jié)構(gòu)對化肥施用的影響。顏姜慧和朱舜[7]、劉穎等[8]認為，城市群發(fā)展吸引了眾多農(nóng)村人口向城市聚集，實現(xiàn)了農(nóng)業(yè)勞動力兼業(yè)化與非農(nóng)化。進一步，劉華軍和劉傳明[9]、向晶和鐘甫寧[10]研究指出，農(nóng)村青壯年勞動力的配置非農(nóng)化使得農(nóng)村農(nóng)業(yè)勞動力老齡化速度加快，對農(nóng)業(yè)有序生產(chǎn)有消極影響。李光泗和朱麗莉[11]、史常亮等[12]認為，大量農(nóng)村人口外流，導致農(nóng)業(yè)勞動力對糧食生產(chǎn)的邊際產(chǎn)出效應開始上升，農(nóng)業(yè)勞動力不足使化肥、機械等要素開始替代勞動。胡浩和楊泳冰[13]、張艦等[14]、欒江[15]進一步研究指出，隨著城市群的發(fā)展，勞動力市場交易機制逐漸完善，引起農(nóng)業(yè)勞動力用工成本上升，且其增長速度遠遠高于化肥價格的增長速度，這使得農(nóng)戶愿意提高化肥施用強度以替代必要的勞動。張聰穎和霍學喜[16]則認為，農(nóng)村勞動力轉(zhuǎn)移導致的農(nóng)業(yè)生產(chǎn)勞動力不足會促使農(nóng)戶采取測土配方施肥技術(shù)，從而減少了農(nóng)業(yè)生產(chǎn)化肥施用量。不僅如此，張軍偉等[17]研究指出，糧食生產(chǎn)勞動力投入數(shù)量也會對化肥施用存在抑制效應。二是耕地面積變化對化肥施用的影響。潘丹[18]認為，城鎮(zhèn)化擴大了城市用地規(guī)模，導致耕地面積減少，對化肥消費有重要影響。諸培新等[19]、紀龍等[20]研究指出，農(nóng)地經(jīng)營規(guī)模越大、土地集中程度越高，越有利于促使農(nóng)戶減少化肥投入。而黎孔清和馬豆豆[21]則認為，地塊破碎程度越高，農(nóng)地經(jīng)營越不集中，越不利于減少化肥投入。三是經(jīng)濟發(fā)展對化肥施用的影響。城鎮(zhèn)化發(fā)展促進三產(chǎn)結(jié)構(gòu)調(diào)整、經(jīng)濟發(fā)展、技術(shù)進步和收入的增加，這些因素也會對化肥施用產(chǎn)生顯著影響。潘丹[18]研究指出，三產(chǎn)結(jié)構(gòu)的調(diào)整將提高化肥施用強度。張軍偉等[17]研究發(fā)現(xiàn)，地區(qū)經(jīng)濟發(fā)展水平也顯著促進了化肥投入的增加。而楊建輝[23]則認為，農(nóng)業(yè)經(jīng)濟發(fā)展水平與農(nóng)業(yè)化學品投入存在絕對脫鉤關(guān)系，即農(nóng)業(yè)經(jīng)濟發(fā)展對其沒有影響。這與郭相平等[24]研究結(jié)論相反，他認為隨著經(jīng)濟發(fā)展可以減少農(nóng)作物化肥施用量，進而減少農(nóng)作物水足跡。

綜上所述，雖然相關(guān)研究較為豐富，但城鎮(zhèn)化對化肥施用強度影響的文獻較少，且多數(shù)文獻的研究視角較為單一，也未指出城鎮(zhèn)化對化肥施用強度的影響存在區(qū)域差異。因此，本文構(gòu)建城鎮(zhèn)化水平（人口、土地、經(jīng)濟城鎮(zhèn)化）對糧食生產(chǎn)化肥施用強度影響的綜合分析框架，基于2001—2017年全國31省糧食生產(chǎn)的面板數(shù)據(jù)，采用面板VAR模型和Tobit模型對城鎮(zhèn)化與化肥施用強度進行脈沖響應分析和回歸分析。以期為化肥減量投入，實現(xiàn)糧食生產(chǎn)現(xiàn)代化、綠色化相統(tǒng)一和保障糧食安全提供參考。

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

（一）研究方法

1.熵值法。其基本原理是基于各觀測指標的差異度，合理確定其權(quán)重。賦權(quán)步驟如下：

（1）數(shù)據(jù)標準處理。若i個對象有j個觀測指標Zj，對其標準處理后有（假定標準化后指標值用Xij表示）：

Xij=Zij-minZijmaxZij-minZij（1）

（2）信息熵計算。對于i個對象有j個觀測指標Zj，則其信息熵為：

βij=Xij∑ni=1XijIj=-lnn-1∑ni=1βijlnβij（2）

式（2）中，I為信息熵。當且僅當β=0時，有βijlnβij收斂于0值。

（3）權(quán)重確定。假定ωj為各觀測指標權(quán)重，由式（2）可得：

ωj=1-Ij∑1-Ij（3）

2.面板VAR模型。相較于VAR模型，該模型引入面板結(jié)構(gòu)，注重個體間的差異，動態(tài)估計參數(shù)方程。模型如下：

Hit=Hi，t-1B1+Hi，t-2B2+...+Hi，t-k+1Bk-1+Hi，t-kBk+RitD+ui+εit（4）

Ehitεit=Eritεit=Euiεit=0（5）

式（4）中，i為省份，t為年份，H為n個內(nèi)生變量組成的行向量，R為m個外生變量組成的行向量，B、D為待估參數(shù)，ui為個體效應，εit為誤差項。對于面板VAR模型，采用動態(tài)GMM估計，使結(jié)果更為穩(wěn)健。

3.Tobit模型。該模型運用MLE估計的結(jié)果更穩(wěn)健，而傳統(tǒng)OLS估計忽視了因變量的約束條件，難以避免估計偏差。故采用如下Tobit模型：

yit=α0+x′itγ+ui+εit，若α0+x′itγ+ui+εit>00，若α0+x′itγ+ui+εit≤0 （6）

（二）數(shù)據(jù)來源

綜合考慮，截取2001—2017年作為研究時段。文中數(shù)據(jù)包括：《中國統(tǒng)計年鑒（2002—2018）》《中國人口和就業(yè)統(tǒng)計年鑒（2002—2018）》《中國農(nóng)村統(tǒng)計年鑒（2002—2018）》《新中國六十年統(tǒng)計資料匯編》，部分數(shù)據(jù)來源于《中國城市統(tǒng)計年鑒（2002—2017）》和歷年各省統(tǒng)計年鑒及統(tǒng)計公報等。文中缺失數(shù)據(jù)采用歷年增長率進行補充，為消除各變量統(tǒng)計單位差異，對數(shù)據(jù)標準化。

二、城鎮(zhèn)化綜合指標體系構(gòu)建

（一）指標選取

城鎮(zhèn)化正在經(jīng)歷前所未有的大浪潮，為經(jīng)濟高質(zhì)量發(fā)展注入了強勁動力，推動了三次產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和收入增加。但城鎮(zhèn)化發(fā)展吸引了大量農(nóng)村青壯年向城市轉(zhuǎn)移，實現(xiàn)農(nóng)業(yè)勞動力就業(yè)非農(nóng)化，影響了農(nóng)業(yè)勞動力生產(chǎn)要素投入;城市用地規(guī)模擴大，不可避免地會影響農(nóng)用可耕地面積。因此，選取人口城鎮(zhèn)化、土地城鎮(zhèn)化、經(jīng)濟城鎮(zhèn)化綜合衡量城鎮(zhèn)化水平。

1.人口城鎮(zhèn)化。城市經(jīng)濟的快速發(fā)展，使得農(nóng)業(yè)比較利益降低，出于理性決策，大量農(nóng)村青壯年向城市聚集，使得城市人口比重上升，客觀上實現(xiàn)了農(nóng)業(yè)勞動力非農(nóng)化配置，但影響了農(nóng)業(yè)生產(chǎn)要素投入結(jié)構(gòu)平衡。對于農(nóng)業(yè)現(xiàn)代化程度較低、經(jīng)濟落后或偏遠地區(qū)，其糧食生產(chǎn)對勞動力和化肥等傳統(tǒng)要素投入依賴程度大，勞動力非農(nóng)化使得糧食生產(chǎn)要素投入不足。一方面，為彌補勞動力投入不足帶來的不利影響，以保證基本糧食生產(chǎn)，農(nóng)戶會在一定程度上增施化肥。另一方面，農(nóng)戶會減少糧食生產(chǎn)或?qū)崿F(xiàn)耕地用途非農(nóng)化，也影響了化肥施用。

2.土地城鎮(zhèn)化。城市用地規(guī)模增加，會占用耕地，使得糧食生產(chǎn)耕地面積減少，進而減少化肥施用。同時，城市人口增加對糧食需求增大，在有限的耕地面積內(nèi)進行糧食生產(chǎn)面臨重要挑戰(zhàn)。因此，農(nóng)業(yè)現(xiàn)代化程度較高的發(fā)達地區(qū)會調(diào)整化肥等要素投入結(jié)構(gòu)以保證糧食產(chǎn)量，而對于農(nóng)業(yè)現(xiàn)代化程度較低的偏遠落后地區(qū)，由于難以優(yōu)化糧食生產(chǎn)要素投入結(jié)構(gòu)，且受糧食生產(chǎn)要素投入慣性影響，農(nóng)戶可能增加化肥投入，影響了化肥施用。

3.經(jīng)濟城鎮(zhèn)化。城鎮(zhèn)化發(fā)展優(yōu)化了產(chǎn)業(yè)結(jié)構(gòu)，提高了人均可支配收入水平。一方面，人均可支配收入增加，有利于優(yōu)化消費結(jié)構(gòu)，提升消費水平，對糧食生產(chǎn)要素消費結(jié)構(gòu)產(chǎn)生了影響。另一方面，產(chǎn)業(yè)結(jié)構(gòu)升級帶來的是二三產(chǎn)業(yè)產(chǎn)值比重增加，有利于實現(xiàn)工業(yè)反哺農(nóng)業(yè)、城市支持農(nóng)村的戰(zhàn)略目標，使傳統(tǒng)農(nóng)業(yè)向現(xiàn)代農(nóng)業(yè)轉(zhuǎn)變，客觀上也影響了糧食化肥施用。

（二）指標賦權(quán)

基于熵值法，我們對城鎮(zhèn)化各項指標進行賦權(quán)（表1）?？傮w而言，土地城鎮(zhèn)化（51.8%）>經(jīng)濟城鎮(zhèn)化（33.31%）>人口城鎮(zhèn)化（14.89%）。具體來看，城市征用土地面積占比最大（32.5%），人均非農(nóng)產(chǎn)值其次（17.66%），耕地面積和人均地區(qū)生產(chǎn)總值占比緊隨其后（15.5%），二三產(chǎn)業(yè)產(chǎn)值比占比最低（0.16%）。

三、城鎮(zhèn)化對化肥施用強度分析

（一）城鎮(zhèn)化對化肥施用強度的脈沖響應分析

基于赤池信息準則（Akaike information? criterion，AIC）、貝葉斯信息準則（Bayesian information criterion，BIC）、HQ準則（hannan-quinn information criterion，HQIC）等信息準則，綜合確定最優(yōu)滯后期數(shù)，并將研究時段劃分為17期，對面板VAR模型進行了脈沖響應分析（圖1）。

就全國而言，城鎮(zhèn)化對糧食生產(chǎn)化肥施用強度脈沖響應總體正向顯著。其中，在城鎮(zhèn)化開始沖擊后，化肥施用強度顯正響應;在第1期至第3期城鎮(zhèn)化對化肥施用強度的顯正響應水平下降;第3期開始變得不顯著;第3期之后又轉(zhuǎn)變?yōu)檎蝻@著響應;在第10期之后又開始緩慢下降，并逐漸趨于不顯著。這表明，在研究時段的前幾年城鎮(zhèn)化綜合水平提高了化肥施用強度，但在2001至2003年城鎮(zhèn)化綜合水平對化肥施用強度的增幅效應開始下降，并在2003年下降為0，之后緩慢增長，最后在2010年開始下降。

在主產(chǎn)區(qū)，當城鎮(zhèn)化綜合水平開始沖擊后，化肥施用強度并沒有響應，而是隨著時間推移開始顯正響應;在第2期顯正響應最大;在第3期化肥施用強度沒有響應城鎮(zhèn)化的沖擊，之后又顯正響應;在第10期開始下降，逐漸趨于不顯著。這表明，在初期城鎮(zhèn)化綜合水平的提高對化肥施用強度沒有影響，隨著時間推移，提高了化肥施用強度;在2002年城鎮(zhèn)化的增幅效應達到最大;在2003年其增幅效應下降為0，之后緩慢增長;最后在2010年又開始下降，并逐漸趨于0。

在主銷區(qū)，當城鎮(zhèn)化綜合水平開始沖擊，化肥施用強度響應水平不斷減弱;在第1期減弱為0，于第2期顯著負向沖擊;于第4期沖擊為0，隨后又開始正向沖擊;在第9期沖擊力度減弱，逐漸趨于0。這表明，在城鎮(zhèn)化沖擊初期，城鎮(zhèn)化綜合水平的提高促進了化肥施用強度的提高;在2001至2004年，化肥施用強度下降，并在2002年其減幅效應達到最大;在2004年之后開始轉(zhuǎn)為增幅效應，但在2009年其增幅效應又開始減弱，并逐漸趨于0。

在平衡區(qū)，當城鎮(zhèn)化綜合水平開始沖擊，化肥施用強度響應水平不斷減弱;在第2期不響應，于第3期顯著負向沖擊;于第4期沖擊為0，隨后正向沖擊;在第7期顯著正向沖擊，隨后減弱并逐漸趨于0。這表明，在城鎮(zhèn)化沖擊初期，城鎮(zhèn)化綜合水平的提高促進了化肥施用強度的提高;在2002至2004年，化肥施用強度下降，并在2003年其減幅效應達到最大;在2004年之后開始轉(zhuǎn)為增幅效應，但在2007年其增幅效應開始減弱，并逐漸趨于0。

此外，地理區(qū)位優(yōu)勢造就的經(jīng)濟發(fā)展水平對化肥施用強度的沖擊也會有差異。因此，基于上述三種信息準則，綜合確定最優(yōu)滯后期數(shù)，同樣將研究時段劃分為17期，對面板VAR模型繼續(xù)進行脈沖響應分析（圖2）。

就東部地區(qū)而言，當城鎮(zhèn)化綜合水平開始沖擊后，化肥施用強度響應水平有小幅增加，在第1期達到最大;隨后開始減弱，在第2至第3期化肥施用強度響應水平顯著為負;之后繼續(xù)上升，并在第5期出現(xiàn)第二個正向響應峰值;但隨后又開始減弱，逐漸趨于0。這表明，在沖擊初期，城鎮(zhèn)化水平提高能增大化肥施用強度，但在2002至2003年降低了化肥施用強度;在2005年又顯著增大了化肥施用強度;2005年之后化肥施用強度增幅效應開始減弱，并逐漸變得不顯著。

就中部地區(qū)而言，在第1至第17期內(nèi)，當城鎮(zhèn)化開始沖擊后，化肥施用強度基本沒有響應。這說明在研究時段內(nèi)，城鎮(zhèn)化水平的提高對化肥施用強度的影響不顯著或沒有影響。

就西部地區(qū)而言，當城鎮(zhèn)化開始沖擊后化肥施用強度的響應速度減弱，在第1至第2期響應速度基本為0，在第3期響應速度顯著為負，在第5期及之后其響應速度顯著為正。這說明在研究時段的前幾年隨著城鎮(zhèn)化水平不斷提高，推動了化肥施用強度的增加，但增幅效應有所下降;在2003年其對化肥施用強度的影響顯著為負;在2005年之后又轉(zhuǎn)變?yōu)轱@著正向影響。

（二）城鎮(zhèn)化對化肥施用強度的回歸分析

1.Tobit模型構(gòu)建

基于城鎮(zhèn)化綜合評價指標體系，構(gòu)建如下面板Tobit模型：

Pit=0.176 1P1it+0.149 9P2it+0.674 0P3it（7）

Lit=0.627 5L1it+0.298 1L2it+0.074 4L3it（8）

Eit=0.464 9E1it+0.530 3E2it+0.004 9E3it（9）

URit=0.148 9Pit+0.518 0Lit+0.333 1Eit（10）

Fit=β0+βitURit+ui+εit（11）

式中，i為省份;t為年份;P為人口城鎮(zhèn)化綜合指數(shù)，其中，P1為城鎮(zhèn)人口比重，P2為非農(nóng)就業(yè)人口比重，P3為城市人口密度;L為土地城鎮(zhèn)化綜合指數(shù)，其中，L1為城市征用土地面積，L2為耕地面積，L3為人均公共綠地面積;E為經(jīng)濟城鎮(zhèn)化綜合指數(shù)，其中，E1為人均地區(qū)生產(chǎn)總值，E2為人均非農(nóng)產(chǎn)值，E3為二三產(chǎn)業(yè)產(chǎn)值比重;UR為城鎮(zhèn)化綜合水平;F為化肥施用強度（用單位面積糧食生產(chǎn)化肥施用量表示）;β為常系數(shù);ui為個體效應;εit為誤差項。各變量描述性統(tǒng)計如下（表2）。

2.單位根檢驗與協(xié)整檢驗

為保證數(shù)據(jù)回歸結(jié)果真實，先采用HT與IPS檢驗方法進行面板單位根檢驗。由于部分變量在零階和一階處不能拒絕單位根存在的假設，故進行二階單位根測試。結(jié)果表明，各變量存在二階單整（表3）。

采用Pedroni協(xié)整檢驗以驗證城鎮(zhèn)化與化肥施用強度的長期均衡關(guān)系，在協(xié)整檢驗中必須保證各變量同階單整。同時，由于變量較多，將所有變量分為三組進行檢驗（第一組：F與P1、P2、P3、L1、L2;第二組：F與P2、P3、L1、L2、L3、E1;第三組F與P3、L1、L2、L3、E1、E2、E3），檢驗結(jié)果均強烈拒絕沒有協(xié)整關(guān)系的原假設，表明城鎮(zhèn)化與化肥施用強度存在長期均衡關(guān)系（表4）。

3.Tobit回歸分析

首先進行OLS回歸作為Tobit回歸結(jié)果的對照（表5、表7）。關(guān)于Tobit模型是否采取固定效應或隨機效應模型回歸，通過似然比LR檢驗（表6、表8），結(jié)果表明強烈拒絕固定效應假設，認為有個體效應存在，所以采取包含個體的隨機效應面板回歸。不同城鎮(zhèn)化水平在三大糧食產(chǎn)區(qū)對化肥施用強度也會有差異，結(jié)果如表6所示。

（1）全國層面。人口城鎮(zhèn)化對化肥施用強度有顯著正向影響，即，隨著城鎮(zhèn)人口比重上升會增加化肥施用強度。這說明城市化的快速發(fā)展吸引了大量農(nóng)業(yè)勞動力向城市轉(zhuǎn)移，農(nóng)村勞動力就業(yè)非農(nóng)化，導致糧食生產(chǎn)勞動力要素投入不足。對于全國農(nóng)業(yè)現(xiàn)代化水平總體偏低的局面，且受傳統(tǒng)農(nóng)業(yè)生產(chǎn)要素投入慣性影響，農(nóng)戶會通過增加化肥投入以盡可能減少由于勞動力非農(nóng)轉(zhuǎn)移對糧食生產(chǎn)帶來的不利影響。土地城鎮(zhèn)化對化肥施用強度沒有影響，說明城市用地規(guī)模增加對糧食生產(chǎn)的耕地投入影響不顯著，對化肥施用沒有產(chǎn)生影響。經(jīng)濟城鎮(zhèn)化對化肥施用強度影響顯著。其中，人均地區(qū)生產(chǎn)總值與二三產(chǎn)業(yè)產(chǎn)值比重對其有顯正影響，人均非農(nóng)產(chǎn)值對化肥施用強度有顯著負向影響。前者說明，經(jīng)濟發(fā)展帶來的人均收入增加在理性決策下可能被用于非糧領(lǐng)域，并沒有改善糧食生產(chǎn)要素投入結(jié)構(gòu);二三產(chǎn)業(yè)產(chǎn)值比重上升帶來的工業(yè)反哺農(nóng)業(yè)的激勵效果不明顯，使得農(nóng)戶仍采用經(jīng)驗性的施肥習慣，農(nóng)業(yè)現(xiàn)代化之路任重道遠;城市化快速發(fā)展使得城市人口增加對糧食需求增加，在有限耕地面積約束條件下，糧食生產(chǎn)任務重大。因此，農(nóng)戶會通過大量增加化肥、農(nóng)藥等要素投入刺激糧食生產(chǎn)，以保障糧食供需平衡。后者說明，農(nóng)業(yè)比較利益較低，且二三產(chǎn)業(yè)產(chǎn)值比重增加，農(nóng)戶出于理性決策，會實現(xiàn)勞動力非農(nóng)轉(zhuǎn)移或兼業(yè)化，從而減少或放棄一部分糧食生產(chǎn)，減緩了化肥施用強度。

（2）主產(chǎn)區(qū)層面。人口城鎮(zhèn)化對化肥施用強度有顯著正向影響。這說明在主產(chǎn)區(qū)勞動力要素配置對糧食生產(chǎn)至關(guān)重要，隨著大量農(nóng)業(yè)人口非農(nóng)轉(zhuǎn)移，影響了要素投入平衡結(jié)構(gòu)，農(nóng)戶為確保糧食生產(chǎn)會在一定程度上增施化肥彌補勞動力投入不足的影響。土地城鎮(zhèn)化對化肥施用強度影響顯著。其中，人均公共綠地面積對其有顯著正向影響，而耕地面積對其也有顯著正向影響。前者說明主產(chǎn)區(qū)作為全國糧食安全保障重心，人均公共綠地面積擴大帶來的耕地面積減少使其在短時間內(nèi)難以優(yōu)化要素投入結(jié)構(gòu)，為了保障糧食安全，客觀上會增施化肥促進生產(chǎn)。后者說明城鎮(zhèn)化水平的提高引起城市人口快速增加，人均公共綠地面積減少。為促進城鎮(zhèn)均衡發(fā)展，會擴大城市用地規(guī)模，導致耕地面積與糧食生產(chǎn)用地減少，在一定程度上減少了化肥施用。經(jīng)濟城鎮(zhèn)化對化肥施用強度影響顯著。其中，人均地區(qū)生產(chǎn)總值與全國相比其正向影響減弱。這說明隨著收入增加主產(chǎn)區(qū)農(nóng)民會改善糧食生產(chǎn)條件，增加產(chǎn)出效率更高的要素投入，客觀上替代了化肥、農(nóng)藥等部分傳統(tǒng)要素，使得化肥施用強度減弱。人均非農(nóng)產(chǎn)值負向影響減弱，這說明非農(nóng)領(lǐng)域的豐厚收益吸引了大批農(nóng)業(yè)勞動力。然而，主產(chǎn)區(qū)較豐厚的種糧補貼在一定程度上減緩勞動力的外流規(guī)模，對糧食生產(chǎn)要素投入減少有一定抑制作用，使得化肥減施速度放緩。二三產(chǎn)業(yè)產(chǎn)值比重正向影響顯著，與全國層面一樣，非農(nóng)產(chǎn)值比重越低，化肥施用強度越大。

（3）主銷區(qū)層面。人口城鎮(zhèn)化對化肥施用強度有較為顯著的負向影響，這表明農(nóng)業(yè)勞動力非農(nóng)轉(zhuǎn)移，會引起糧食生產(chǎn)勞動力投入不足，使得農(nóng)民減少糧食生產(chǎn)，降低了化肥施用強度。同時，主銷區(qū)大多處于經(jīng)濟發(fā)達地區(qū)，人均購買力較高。因此，農(nóng)民也會減少糧食生產(chǎn)，也降低了化肥施用強度。土地城鎮(zhèn)化對化肥施用強度影響較為顯著。其中，耕地面積對其有較為顯著的負向影響，人均公共綠地面積對化肥施用強度正向影響顯著。由于主銷區(qū)經(jīng)濟發(fā)展程度高，土地資源寶貴，耕地面積減少影響了基本糧食生產(chǎn)。隨著城市人口不斷增加，對糧食需求不斷擴大，促使農(nóng)民增施化肥增加糧食產(chǎn)量。經(jīng)濟城鎮(zhèn)化對化肥施用強度有顯著正向影響。與二三產(chǎn)業(yè)相比農(nóng)業(yè)比重下降，此時農(nóng)戶可能會增施化肥以維護本地糧食生產(chǎn)。

（4）平衡區(qū)層面。人口城鎮(zhèn)化對化肥施用強度正向影響顯著，如前所述，農(nóng)戶在一定程度上增施化肥保障糧食生產(chǎn)。在土地城鎮(zhèn)化中，耕地面積對化肥施用強度有較為顯著的負向影響。由于平衡區(qū)土地城鎮(zhèn)化減少了耕地面積，其糧食生產(chǎn)受阻。且受其農(nóng)業(yè)現(xiàn)代化程度低和農(nóng)業(yè)生產(chǎn)慣性的雙重影響，農(nóng)戶可能會增加化肥施用，以保障糧食正常生產(chǎn)。經(jīng)濟城鎮(zhèn)化對化肥施用強度正向影響顯著，如前所述，農(nóng)戶也會增加化肥投入。

地理區(qū)位優(yōu)勢也會影響化肥施用強度，綜合地區(qū)經(jīng)濟發(fā)展水平與糧食生產(chǎn)特征進行回歸，分析不同地理區(qū)位對化肥施用強度的影響，結(jié)果如表8所示。

（1）東部地區(qū)。人口城鎮(zhèn)化對化肥施用強度有顯著的負向影響?？赡茉颍阂皇菛|部地區(qū)靠近沿海，經(jīng)濟發(fā)展程度高，勞動力以非農(nóng)就業(yè)為主，使得農(nóng)業(yè)生產(chǎn)勞動力少，在一定程度上減少了糧食生產(chǎn)。二是人均收入高，購買力較強，也減少了糧食生產(chǎn)。三是東部地區(qū)大多處于主銷區(qū)，糧食以外購為主。因此，綜合原因促使東部地區(qū)人口城鎮(zhèn)化水平越高而化肥施用強度越低。土地城鎮(zhèn)化對化肥施用強度影響顯著。其中，耕地面積對其有較顯著的負向影響，表明城市規(guī)模的擴大，使得耕地減少，糧食生產(chǎn)面臨重要挑戰(zhàn)，此時農(nóng)戶會增施化肥以求短時間穩(wěn)定糧食生產(chǎn)。同樣，人均公共綠地面積的增加也使得化肥投入增加。經(jīng)濟城鎮(zhèn)化對化肥施用強度正向影響顯著，如前所述，農(nóng)戶也會增施化肥，以保證糧食基本生產(chǎn)。

（2）中部地區(qū)。人口城鎮(zhèn)化對化肥施用強度有正影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持），原因是中部地區(qū)大多位于糧食主產(chǎn)區(qū)，農(nóng)業(yè)勞動力資源寶貴，如前所述，其大量非農(nóng)轉(zhuǎn)移會促使農(nóng)戶使用化肥替代，客觀上增加了化肥施用強度。土地城鎮(zhèn)化對化肥施用強度有顯著影響。其中，城市征用土地面積對其有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持），耕地面積對其有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。一方面，城市用地規(guī)模擴大，導致耕地面積減少，使得糧食生產(chǎn)下降，減少了化肥施用。另一方面，城市人口大量增加對糧食需求增大，促使農(nóng)民增施化肥提高糧食產(chǎn)出效率。經(jīng)濟城鎮(zhèn)化對化肥施用強度影響顯著，其中人均地區(qū)生產(chǎn)總值和二三產(chǎn)業(yè)產(chǎn)值比重對其有正向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持），人均非農(nóng)產(chǎn)值對其有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。如前所述，前者引起糧食供需失衡和農(nóng)業(yè)產(chǎn)值比重下降，會推動化肥投入的增加。后者引起生產(chǎn)要素的非農(nóng)配置，從而減少化肥投入。

（3）西部地區(qū)。人口城鎮(zhèn)化對化肥施用強度正向影響顯著，原因是西部地區(qū)經(jīng)濟發(fā)達程度和農(nóng)業(yè)現(xiàn)代化程度低，農(nóng)業(yè)勞動力大量非農(nóng)化，使得該地區(qū)糧食生產(chǎn)勞動力投入不足，客觀上推動了農(nóng)戶利用化肥替代勞動力，引起化肥投入增大。土地城鎮(zhèn)化對化肥施用強度影響不顯著，表明城市用地規(guī)模增加沒有影響到糧食生產(chǎn)，對其要素投入影響不顯著。經(jīng)濟城鎮(zhèn)化對化肥施用強度有正影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持），與糧食產(chǎn)區(qū)和東中部地區(qū)一樣，工業(yè)反哺農(nóng)業(yè)力度有待提高，且資本的非農(nóng)配置，導致化肥施用增加。

四、結(jié)論與啟示

本文基于省際面板數(shù)據(jù)，構(gòu)建包含人口城鎮(zhèn)化、土地城鎮(zhèn)化、經(jīng)濟城鎮(zhèn)化的綜合分析框架，采用熵值法對城鎮(zhèn)化綜合指標賦權(quán)，構(gòu)建面板VAR模型和面板Tobit模型對城鎮(zhèn)化與化肥施用強度進行實證分析。結(jié)果表明：1.人口城鎮(zhèn)化對化肥施用強度影響顯著。就全國層面而言，人口城鎮(zhèn)化對化肥施用強度有正向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。就糧食功能區(qū)而言，在主產(chǎn)區(qū)和平衡區(qū)人口城鎮(zhèn)化對化肥施用強度有正向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持），在主銷區(qū)人口城鎮(zhèn)化對化肥施用強度有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。就地理區(qū)位而言，在中部和西部地區(qū)人口城鎮(zhèn)化對化肥施用強度正向影響，在東部地區(qū)人口城鎮(zhèn)化對化肥施用強度有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。2.土地城鎮(zhèn)化對化肥施用強度影響存在顯著差異。就全國層面而言，土地城鎮(zhèn)化對化肥施用強度影響未取得統(tǒng)計數(shù)據(jù)支持。就糧食功能區(qū)而言，土地城鎮(zhèn)化在主銷區(qū)和平衡區(qū)均對其有正向影響，但在主產(chǎn)區(qū)存在顯著差異（該結(jié)果取得統(tǒng)計數(shù)據(jù)支持）。具體表現(xiàn)為：耕地面積減少降低了化肥施用強度，人均公共綠地面積增加增大了化肥施用強度。就地理分區(qū)而言，在東部地區(qū)土地城鎮(zhèn)化對化肥施用強度正向影響顯著;在中部地區(qū)土地城鎮(zhèn)化對其影響存在顯著差異。具體表現(xiàn)為：城市征用土地面積引起化肥施用強度降低，耕地面積減少引起化肥施用強度增大;在西部地區(qū)土地城鎮(zhèn)化對化肥施用強度影響不顯著（該結(jié)果未取得統(tǒng)計數(shù)據(jù)的支持）。3.經(jīng)濟城鎮(zhèn)化對化肥施用強度影響差異顯著。就全國層面而言，經(jīng)濟城鎮(zhèn)化對化肥施用強度影響顯著，其中人均地區(qū)生產(chǎn)總值與二三產(chǎn)業(yè)產(chǎn)值比重對化肥施用強度有正向影響，人均非農(nóng)產(chǎn)值對化肥施用強度有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。就糧食功能區(qū)而言，在主產(chǎn)區(qū)經(jīng)濟城鎮(zhèn)化對化肥施用強度影響顯著（該結(jié)果取得統(tǒng)計數(shù)據(jù)的支持）。表現(xiàn)為：人均地區(qū)生產(chǎn)總值的正向影響減弱，人均非農(nóng)產(chǎn)值負向影響減弱;在主銷區(qū)與平衡區(qū)經(jīng)濟城鎮(zhèn)化均對化肥施用強度有正向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。就地理區(qū)位而言，在東西部地區(qū)經(jīng)濟城鎮(zhèn)化均對化肥施用強度正向影響顯著，但在中部地區(qū)存在差異。表現(xiàn)為：人均地區(qū)生產(chǎn)總值和二三產(chǎn)業(yè)產(chǎn)值比重對其有正向影響，人均非農(nóng)產(chǎn)值對其有負向影響（該結(jié)果得到統(tǒng)計數(shù)據(jù)支持）。

綜上所述，為了降低化肥施用強度，實現(xiàn)糧食生產(chǎn)綠色化和農(nóng)業(yè)現(xiàn)代化，可從以下幾個方面入手：一是合理引導農(nóng)業(yè)勞動力非農(nóng)轉(zhuǎn)移，對勞動力資源配置進行科學規(guī)劃，推動實現(xiàn)農(nóng)業(yè)勞動力就地就業(yè)，防止農(nóng)業(yè)勞動力過度轉(zhuǎn)移引起化肥替代勞動力現(xiàn)象的發(fā)生。二是適度推進城市化進程，在規(guī)劃城市用地應優(yōu)先保障糧食生產(chǎn)用地，防止過多占用田地導致糧食播種面積減少，引發(fā)糧食供需矛盾而使得化肥施用增加。三是加大力度實施工業(yè)反哺農(nóng)業(yè)，城市支持農(nóng)村的戰(zhàn)略。加大政府支農(nóng)財政力度，大力實施科技興農(nóng)、藏糧于地、藏糧于技的重大戰(zhàn)略，研發(fā)優(yōu)質(zhì)、高產(chǎn)、綠色的糧食作物新品種，推動傳統(tǒng)農(nóng)業(yè)向現(xiàn)代農(nóng)業(yè)轉(zhuǎn)變。

參考文獻：

[1]欒江，仇煥廣，井月，等.我國化肥施用量持續(xù)增長的原因分解及趨勢預測[J].自然資源學報，2013（11）：1869-1878.

[2]仇煥廣，欒昊，李瑾，等.風險規(guī)避對農(nóng)戶化肥過量施用行為的影響[J].中國農(nóng)村經(jīng)濟，2014（3）：85-96.

[3]紀月清，張惠，陸五一，等.差異化、信息不完全與農(nóng)戶化肥過量施用[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2016（2）：14-22.

[4]孔凡斌，郭巧苓，潘丹.中國糧食作物的過量施肥程度評價及時空分異[J].經(jīng)濟地理，2018（10）：201-210.

[5]陸文聰，劉聰.化肥污染對糧食作物生產(chǎn)的環(huán)境懲罰效應[J].中國環(huán)境科學，2017（5）：1988-1994.

[6]周亮，徐建剛，蔡北溟，等.淮河流域糧食生產(chǎn)與化肥消費時空變化及對水環(huán)境影響[J].自然資源學報，2014（6）：1053-1064.

[7]顏姜慧，朱舜.農(nóng)村人口轉(zhuǎn)移趨勢及空間指向研究[J].中國人口·資源與環(huán)境，2017（5）：146-152.

[8]劉穎，鄧偉，宋雪茜，等.基于綜合城鎮(zhèn)化視角的省際人口遷移格局空間分析[J].地理科學，2017（8）：1151-1158.

[9]劉華軍，劉傳明.城鎮(zhèn)化與農(nóng)村人口老齡化的雙向反饋效應：基于中國省際面板數(shù)據(jù)聯(lián)立方程組的經(jīng)驗估計[J].農(nóng)業(yè)經(jīng)濟問題，2016（1）：45-52.

[10]向晶，鐘甫寧.農(nóng)村人口轉(zhuǎn)移、工業(yè)化和城鎮(zhèn)化[J].農(nóng)業(yè)經(jīng)濟問題，2018（12）：51-56.

[11]李光泗，朱麗莉.農(nóng)村勞動力流動背景下我國糧食生產(chǎn)技術(shù)變動分析[J].中國科技論壇，2014（7）：143-148.

[12]史常亮，朱俊峰，欒江.我國小麥化肥投入效率及其影響因素分析：基于全國15個小麥主產(chǎn)省的實證[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2015（11）：69-78.

[13]胡浩，楊泳冰.要素替代視角下農(nóng)戶化肥施用研究：基于全國農(nóng)村固定觀察點農(nóng)戶數(shù)據(jù)[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2015（3）：84-91.

[14]張艦，亞伯拉罕·艾賓斯坦，瑪格麗特·麥克米倫，等.農(nóng)村勞動力轉(zhuǎn)移、化肥過度使用與環(huán)境污染[J].經(jīng)濟社會體制比較，2017（3）：149-160.

[15]欒江.農(nóng)業(yè)勞動力轉(zhuǎn)移與化肥施用存在要素替代關(guān)系嗎？：來自我國糧食主要種植省份的經(jīng)驗證據(jù)[J].西部論壇，2017（4）：12-21.

[16]張聰穎，霍學喜.勞動力轉(zhuǎn)移對農(nóng)戶測土配方施肥技術(shù)選擇的影響[J].華中農(nóng)業(yè)大學學報（社會科學版），2018（3）：65-72.

[17]張軍偉，張錦華，吳方衛(wèi).糧食生產(chǎn)中化肥投入的影響因素研究：基于Durbin模型的分析[J].經(jīng)濟地理，2018（11）：174-182.

[18]潘丹.中國化肥消費強度變化驅(qū)動效應時空差異與影響因素解析[J].經(jīng)濟地理，2014（3）：121-126.

[19]諸培新，蘇敏，顏杰.轉(zhuǎn)入農(nóng)地經(jīng)營規(guī)模及穩(wěn)定性對農(nóng)戶化肥投入的影響：以江蘇四縣（市）水稻生產(chǎn)為例[J].南京農(nóng)業(yè)大學學報（社會科學版），2017（4）：85-94.

[20]紀龍，徐春春，李鳳博，等.農(nóng)地經(jīng)營對水稻化肥減量投入的影響[J].資源科學，2018（12）：2401-2413.

[21]黎孔清，馬豆豆.生態(tài)脆弱區(qū)農(nóng)戶化肥減量投入行為及決策機制研究：以山西省4縣421戶農(nóng)戶為例[J].南京農(nóng)業(yè)大學學報（社會科學版），2018（5）：138-145.

[22]潘丹.中國化肥施用強度變動的因素分解分析[J].華南農(nóng)業(yè)大學學報（社會科學版），2014（2）：24-31.

[23]楊建輝.農(nóng)業(yè)化學投入與農(nóng)業(yè)經(jīng)濟增長脫鉤關(guān)系研究：基于華東6省1市數(shù)據(jù)[J].自然資源學報，2017（9）：1517-1527.

[24]郭相平，高爽，吳夢洋，等.中國農(nóng)作物水足跡時空分布與影響因素分析[J].農(nóng)業(yè)機械學報，2018（5）：295-302.

（責任編輯：楊 洋）

The Influence of Urbanization on the Intensity of Fertilizer Usage：

the Analysis Based on the Fact of Grain Production in 31

Provinces（Autonomous Ragions and Municipalities）

WU Guoyong1，2，BAO Yinghui3

（1.China Institute of Rural Revitalization in Karst Areas，Guizhou University，Guiyang，Guizhou，China，550025;2.Guizhou High-end Think Tank for Grassroots Social Governance，Guizhou University，Guiyang，Guizhou，China，550025;3.Subsidiary Unit of Groupe ESC Clermont-Graduate School of Management，Clermont-Ferrand，Pue-de-Dome，F(xiàn)rance，63000）

Abstract：

This paper aims to measure the impact of urbanization on fertilizer usage intensity so as to provide guidance for reducing fertilizer usage intensity to achieve green and modern grain production.Based on the panel data of 31 provinces in China from 2001 to 2017，this paper adopts the entropy method to weight the constructed comprehensive urbanization indicators，and constructs panel VAR model and panel Tobit model to analyze respectively the impulse response and regression of the relation between urbanization and fertilizer usage intensity.The results find that： firstly，population urbanization has a significant influence on fertilizer usage intensity. positive effects are across the whole country，main production areas，balance areas，central and western regions，while the negative effects can be seen in the main sales areas and eastern regions.Secondly，land urbanization has a significant impact on the intensity of fertilizer usage.The influence was not statistically tenable in the whole country and the western regions，while the main sales areas，balance areas and eastern regions are influenced in a positive way.Besides，the cultivated land area and per capita green land area in the main producing areas have a positive impact on the fertilizer usage intensity，while the cultivated land area，as well as the urban land requisition area in central China has negative influence on it.Thirdly，the influence of economic urbanization on fertilizer usage intensity differs significantly.The ratio of GDP per capita to the value of secondary and tertiary industries in the whole country，major producing areas and the central regions has a positive impact，while per capita non-agricultural value has a negative impact.It has positive influence in main sales area，balance areas，eastern and western areas.Based on the results of empirical analysis，some suggestions are put forward to realize the reduction of fertilizer input，the green grain production and agricultural modernization.

Key words：

fertilizer usage intensity; urbanization; food production; impulse response

3185500589256